Recent expansion and adaptive evolution of the carcinoembryonic antigen family in bats of the Yangochiroptera subgroup
© The Author(s). 2017
Received: 22 March 2017
Accepted: 1 September 2017
Published: 11 September 2017
Expansions of gene families are predictive for ongoing genetic adaptation to environmental cues. We describe such an expansion of the carcinoembryonic antigen (CEA) gene family in certain bat families. Members of the CEA family in humans and mice are exploited as cellular receptors by a number of pathogens, possibly due to their function in immunity and reproduction. The CEA family is composed of CEA-related cell adhesion molecules (CEACAMs) and secreted pregnancy-specific glycoproteins (PSGs). PSGs are almost exclusively expressed by trophoblast cells at the maternal-fetal interface. The reason why PSGs exist only in a minority of mammals is still unknown.
Analysis of the CEA gene family in bats revealed that in certain bat families, belonging to the subgroup Yangochiroptera but not the Yinpterochiroptera subgroup an expansion of the CEA gene family took place, resulting in approximately one hundred CEA family genes in some species of the Vespertilionidae. The majority of these genes encode secreted PSG-like proteins (further referred to as PSG). Remarkably, we found strong evidence that the ligand-binding domain (IgV-like domain) of PSG is under diversifying positive selection indicating that bat PSGs may interact with structurally highly variable ligands. Such ligands might represent bacterial or viral pathogen adhesins. We have identified two distinct clusters of PSGs in three Myotis species. The two PSG cluster differ in the amino acids under positive selection. One cluster was only expanded in members of the Vespertilionidae while the other was found to be expanded in addition in members of the Miniopteridae and Mormoopidae. Thus one round of PSG expansion may have occurred in an ancestry of all three families and a second only in Vespertilionidae. Although maternal ligands of PSGs may exist selective challenges by two distinct pathogens seem to be likely responsible for the expansion of PSGs in Vespertilionidae.
The rapid expansion of PSGs in certain bat species together with selection for diversification suggest that bat PSGs could be part of a pathogen defense system by serving as decoy receptors and/or regulators of feto-maternal interactions.
KeywordsCarcinoembryonic antigen gene family Pregnancy-specific glycoproteins Positive selection Chiroptera Immunoglobulin superfamily Bats
Gene families are predestinated to rapid genetic adaptation to environmental cues in vertebrates. Accelerated gene family expansion, therefore, may provide hints for environmental forces on vertebrate species. For example, in species with an extraordinary large number of γ/δ T-cells, like cattle, sheep and chicken, the CD163 family of immune receptors which are pivotal for the function of γ/δ T-cells is expanded . These species are not very closely related, indicating that the evolution of this gene family is not due to an expansion in a common ancestor but happened independently probably driven by pathogens . Such a species-specific expansion is seen in various gene families, most of them being involved in immunity and/or reproduction. However, the precise function of members of these gene families is often unknown.
This is also true for the carcinoembryonic antigen (CEA) gene family, which belongs to the immunoglobulin superfamily and represents one of the fastest evolving gene families in mammals . In mammals, the ancestral CEA gene family was composed of five genes, i.e. CEACAM1, CEACAM16, CEACAM18, CEACAM19 and CEACAM20. These genes can be identified in almost all mammalian species. The ancestral CEACAM1 was subject to multiple duplications which led to species-specific expansion of CEACAM1-related members of the CEA gene family. CEACAM1 is a transmembrane inhibitory receptor composed of one N-terminal immunoglobulin variable (IgV)-like (also called N domain) and three Ig constant (IgC)-like extracellular domains (also named A1, B, and A2 domains). The IgV-like domain is the primary ligand-binding domain, which was shown to interact with other CEACAMs and other cell surface receptors such as galectins, integrins and TIM-3 as well as with various pathogen adhesins [3, 4]. The cytoplasmic tail of CEACAM1 contains one to two immunoreceptor tyrosine-based inhibition motifs (ITIM). CEACAM1 is expressed by various cell types including, endothelial, epithelial and immune cells. In immune cells CEACAM1 is an important regulator of cell activation [5–7]. In primates and rodents the CEACAM1-related genes belong either to the CEA-related cell adhesion molecule (CEACAM) or the pregnancy-specific glycoprotein (PSG) subgroups. While several CEACAMs are receptors involved in immunity, PSGs are expressed nearly exclusively in trophoblast cells and most likely play a role in maternal-fetal communication . Surprisingly, PSGs do not exist in various mammals including most of the members of the superorder Laurasiatheria [9–11]. However, more recently we and others found that in bats, namely in Myotis lucifugus (M. lucifugus) and Myotis davidii (M. davidii) which also belong to the superorder Laurasiatheria, a considerable gene amplification in the CEA gene family occurred [11, 12]. However, information on the structure and expression of the CEA family members in bats are completely missing.
Bats belong to the order Chiroptera, which is the second largest order of mammals, only rodents contain more species. Traditionally the order Chiroptera was divided into the two suborders Megachiroptera (fruit-eating, non-echolocating bats) and Microchiroptera (insectivorous, echolocating bats). However, current molecular evidence rather favors the division in the new subgroups Yinpterochiroptera and Yangochiroptera which have diverged approximately 60 million years ago [13, 14]. Yinpterochiroptera contain in addition to the old world fruit bats (Pteropodidae), four families of echolocating insectivorous bats. Common to all bats is that they play an important role as reservoirs for viruses. Currently more than 100 viruses have been detected in bats some of them, like lyssa, corona and ebola viruses, are of extraordinary importance for human health . As a consequence it has been speculated that the immune system of bats has unique features making them tolerant to several virus infections . Indeed the continuous threat by various pathogens may have a strong influence on the evolution of immune proteins, including the CEACAM receptors of the CEA family. In addition, we have recently speculated that the expansion of the PSG subgroup of the CEA gene family requires the presence of a hemochorial placenta as it is found in primates and rodents. In contrast to endotheliochorial and epitheliochorial placentae the hemochorial placenta allows direct contact of fetal cells with maternal blood and immune cells . Bats have either an endotheliochorial or a hemochorial placenta, depending on the bat species. Therefore, the high number of CEACAMs in certain bat species raises the question whether these CEACAMs may represent PSGs.
In this report we show that the vast majority of CEACAMs in bats, which have an extended CEA family, are secreted glycoproteins. These glycoproteins are not expressed in a number of tissues in which usually non-PSG CEACAMs are expressed suggesting that the secreted CEACAMs in bats have a restricted expression pattern. Interestingly, the IgV-like domain which is responsible for the interaction with almost all extracellular ligands is under strong positive selection in bats. Selection for diversification points to rapidly evolving ligands, like viruses and other microorganisms or to a family of closely related receptors, like members of a protein family. We hypothesize that these PSG-like proteins (further referred to as PSGs) are expressed at the maternal-fetal interface and that they play a role either in counteracting infection or regulating maternal-fetal communication.
Phylogeny of bat orthologous CEACAM genes
Tremendous expansion of the CEA gene family in certain bat species of the Yangochiroptera subgroup
Next we determined the number of CEACAM1 paralogs in each bat species as described in “Materials and Methods”. Sequences without an open reading frame (ORF) were considered to be part of a pseudogene. Within Yinpterochiroptera maximally 11 N domain exon sequences were identified per species with a maximum of six N domain exons with an ORF within one species (Fig. 1b). Miniopterus natalensis (M. natalensis) a species of the Yangochiroptera suborder has also six N domain exon sequences with an ORF. In all other species of the Yangochiroptera group investigated, a tremendous expansion of CEA family member N domain exons was observed. In M. lucifugus 102 different N domain exons were found, nearly half of them contained an ORF (Fig. 1b). Roughly, a one to one ratio of N domain exons with internal stop codons and N domain exons with an ORF were also found in Myotis brandii (M. brandii), while in M. davidii, Eptesicus fuscus (E. fuscus) and Pteronotus parnellii (P. parnellii) N domain exon sequences with internal stop codons were less frequent (Fig. 1b). Importantly, every N domain exon sequence was separated from another N domain sequence by CEA family-related exon sequences encoding other domains, strongly indicating that in bats each CEA family gene contains only a single N domain exon. Thus the number of N domain exons may indicate the number of CEA gene family members in bats.
A balanced expansion of genes coding for inhibitory and activation CEACAM receptors frequently took place in bats
Genes encoding secreted CEA family members expanded in bat species of the Yangochiroptera suborder
PSG N domains are more closely related with each other than with N domains of membrane anchored CEACAMs in M. lucifugus
PSG N domains of Yangochiroptera bat species
No PSG-like CEACAMs are found in members of the Yinpterochiroptera suborder
The ligand-binding domains of PSGs exhibit positive selection
Sites under positive selection differ between PSG I, PSG II and CEACAM N domains
Bat PSGs have a restricted expression pattern
Myotis myotis CEACAM mRNAs identified by RNA sequencing
In most species of the superorder Laurasiatheria, the CEA gene family is relatively small . However, a few exceptions of this rule exist . Recently, we have identified an expanded CEA gene family in the horse the expansion of which is due to the amplification of genes coding for secreted PSG-like CEACAMs . In addition, there are reports indicating that the CEA gene family has been expanded in certain microbats , namely in Myotis lucifugus  and in Myotis davidii . Such a co-expansion of a gene family in otherwise distantly related species may point to similar selective pressures working on these species during evolution and thereby may provide clues to the function of the gene family. In order to get a comprehensive knowledge of the CEA gene family in bats we have analyzed the CEA gene family in 12 different bat species. According to the phylogenetic tree based on the sequences of orthologous CEACAMs six of these species belong to the Yangochiroptera and six to the Yinpterochiroptera suborders which is in accordance with phylogenetic trees previously reported by other authors [19, 23]. The maximal number of CEA family-related N domains containing an open reading frame in a single species of Yinpterochiroptera suborder was six. In contrast in Yangochiroptera in particular in the Myotis genus up to tenfold as many CEA family member N domains were found. According to the currently proposed phylogeny of Yangochiroptera Miniopteridae are more closely related to Vespertilionidae than to Mormoopidae [19, 24]. Therefore, it is surprising that we found more CEA gene family members in P. parnellii than in M. natalensis. The most plausible explanation is that the expansion of the CEA gene family occurred in both linages independently. Furthermore the difference in the size of the CEA gene family in Myotis and Eptesicus indicates that a second boost of CEA gene family expansion occurred between 20 and 25 Mya years ago in the Vespertilionidae [19, 24].
Multiple CEACAM genes were found to code for transmembrane proteins with signaling capacities through ITIMs and ITAMs in the cytoplasmic tails. Interestingly, amplification of both gene types coding for ITAM- and ITIM-containing CEACAMs occurred. In the CEA gene families of other mammals described previously, a preferential expansion of either genes coding for ITAM-containing proteins (dog) or for ITIM-containing proteins occurred (mouse, horse and, opossum) [9, 11]. Interestingly, a balanced expansion of transmembrane CEACAMs could be found in both Yangochiroptera and Yinpterochiroptera. This may be explained by the fact that multiple ITIM- and ITAM-containing signaling CEACAMs were already present in the last common ancestor of Yangochiroptera and Yinpterochiroptera bat suborders. In the pooled tissues including spleen, thymus, lymph node and intestine we found strong mRNA expression of these signaling CEACAMs, consistent with the view that these CEACAMs are important for immune function. We have previously reported that pairs of ITIM- and ITAM-containing CEACAMs with very similar ligand-binding domains exist in most mammals and even in amphibians [11, 25]. The most plausible explanation for the evolution of these paired receptors is that activation receptors evolved as a countermeasure to the use of the inhibitory receptor by pathogens as cellular receptors. One indication for such a mechanism is that the ligand-binding domain is most similar between ITIM-containing CEACAMs and an ITAM-containing CEACAM which is the case for CEACAMs of Myotis lucifugus (CEACAM1L1/CEACAM1L2 and CEACAM3L1) (Fig. 4).
However, the enormous expansion of the CEA gene family in certain bat species is due to the amplification of genes coding for a single IgV-like domain followed by an IgC-like domain. The IgC-like domains are of the A2 type (named according to the most similar IgC-like domain of CEACAM1) and are encoded by exons which have either stop codons (at the beginning of the exon) or mutated splice donor sites. The mutation of the splice donor site creates a stop codon which is followed in near proximity by a polyadenylation signal. The presence of a leader sequence indicates that these molecules are secreted. Indeed, this is besides expression by trophoblast cells (which could not be demonstrated directly due to lack of bat placental tissues), the most important classification criterion for PSGs. The closer relationship of the ligand-binding domains (IgV-like domain) between PSGs than to other CEACAMs in a given species represents an additional criterion. Indeed this is the case for bat PSGs. Remarkably, the structure of bat PSGs is very similar to the PSGs recently found in the horse  suggesting that both have a common ancestor. Indeed the phylogenetic relationship of bats within Laurasiatheria is still a matter of debate, however several lines of evidence point to a close relationship of bats and horses . For example, Zhang and colleagues used 2492 nuclear-encoded genes to perform maximum-likelihood and Bayesian phylogenomic analysis. Their results vigorously supported bats as a member of Pegasoferae (Chiroptera + Perissodactyla + Carnivora), with the bat lineage diverging from the Equus (horse) lineage ~88 million years ago . Similar findings were obtained on transcriptome level by Papenfuss and coworkers .
Interestingly, phylogenetic analysis indicated that two main groups of PSGs exist in bats. Both groups contain an almost equal number of PSGs in Vespertilionidae. In the more distantly related Pteronotus parnellii PSGs of subgroup II expanded preferentially (Fig. 5). Taken together our data imply that the PSG group II is the more primordial PSG group being present in all analyzed representatives of the Yangochiroptera. The ancestor of PSGs group I may have arisen in a common ancestor of the Miniopteridae and the Vespertilionidae.
What is the reason for the selective expansion of the PSGs in certain families of Yangochiroptera? Until now PSGs are only described in species having a hemochorial placenta . This is the case for all bat species with PSGs. However, hemochorial placentae are also common within the Yinpterochiroptera suborder; this may suggest that a second prerequisite is needed to lead to PSG expansion. This view is further supported by our previous observation that PSGs did not evolve in hedgehogs, which have also a hemochorial placenta . The most obvious prerequisite is that a primordial PSG is created by duplication of a CEA gene family member. This would be a random event with a limited frequency, only occurring in restricted number of mammals with a hemochorial placenta. A second possibility is that a hemochorial placenta is not sufficient to drive PSG evolution but additional specific features of the hemochorial placenta, like special blood flow conditions, invasion depth or immunological challenges are necessary. Indeed, it is well known that placentation of bats is extremely diverse and therefore even placentae with a hemochorial interface may differ considerably [27, 28].
On the other hand the very recent and massive expansion of PSGs makes maternal-fetal communication as the only driving force for PSG evolution in microbats questionable. In particular, positive selection point to an interaction with fast evolving ligands. Ligands fulfilling such requirements are for example pathogen receptors. Indeed several pathogens were described to bind to certain CEACAMs. In humans a variety of bacterial pathogens were identified that bind to various human CEACAMs [29–35]. Furthermore, mouse hepatitis virus, which belongs to the corona viridae group 2 uses CEACAM1 as a cellular receptor to infect susceptible hosts [36, 37]. Bats are known to be prominent reservoirs for corona viruses and therefore it is worthwhile to speculate that in microbats viruses exist or have existed that interact with bat CEACAMs. Indeed, recently a bat corona virus of group 2 was isolated from the common vampire bat Desmodus rotundus . Secreted proteins with some similarities to these CEACAMs may function as decoy receptors and thereby limit virus binding to their cellular receptor. Such an interpretation would be consistent with a rapid expansion and a positive selection of the decoy receptors. This hypothesis is even more exciting for secreted proteins at the maternal-fetal interface, which could be involved in the prevention of transplacental infection. We further speculate that such a mechanism of innate immunity may be especially beneficial for an order of mammals that live in large colonies with synchronized pregnancies and an extraordinary close contact to other individuals. The rapid expansion of PSGs in certain bat species together with selection for diversification suggest that bat PSGs could be part of a pathogen defense system by serving as decoy receptors and/or regulators of feto-maternal interactions.
PSGs are a subgroup of the CEA family. We and others have suggested that maternal-fetal interactions are the drivers of the expansion of PSGs in some mammalian species, including humans and rodents. Both higher primates and rodents have a hemochorial placenta type and the close contact of semi- allogeneic fetal cells with the maternal immune system seems to be responsible for the expansion of PSGs. However, in numerous species although having a hemochorial placenta no expansion of PSGs is observed, arguing against a sole reason of maternal fetal communication for the expansion of PSGs. Our analyses of the CEA gene family in bats suggest that the expansion of PSGs could also be pathogen-driven. Therefore, we favor the hypothesis that a hemochorial placenta is a prerequisite for the expansion of PSGs but additional conditions are needed, for example a continuous threat by pathogens, to initiate PSG expansion. The identification of bat PSGs opens now the possibility to further determine the tissue of bat PSG expression as well as the screening for pathogens that bind to PSGs. Future investigations are warrant to test if PSGs play a role in preventing trans-placental infections.
Data sets and nomenclature of genes
Sequence similarity searches were performed using the NCBI BLAST tools “blastn” http://blast.ncbi.nlm.nih.gov/Blast.cgi and Ensembl BLAST/BLAT search programs http://www.ensembl.org/Multi/Tools/Blast?db=core using default parameters. For identification of bat CEACAM exons, exon and cDNA sequences from known CEACAM and PSG genes were used to search “whole-genome shotgun contigs (wgs)” databases limited to organism “Chiroptera (taxid:9397)”. Hits were considered to be significant if the E-value was < e-10 and the query cover was >50%. Once a wgs contig was identified that contained CEACAM-related sequences we confirmed manually the presence of the complete exon by the number of nucleotides and identification of CEACAM-typical splice site sequences. Only sequences which were considered to be complete exons were used for further analyses. In a second step we used the identified exon sequences to search the database limited to this bat species in order to identify all existing paralogous CEACAM genes. In some species we performed several rounds of searches using sequences of distantly related CEACAMs in a given species. Once we had identified individual exons we predicted the gene structure according to known CEACAMs. The location of different exons on the same contig was a prerequisite for considering that these exons belong to the same gene. Gene predictions were further supported by the identification of “expressed sequence tags (est)” and or predictions in genome builds at NCBI and Ensemble, if available. Short exons, like exons coding for the cytoplasmic tail, were identified by alignments of downstream sequences of identified transmembrane exons with cytoplasmic exon sequences of human CEACAMs. Sequence alignments for exon identification was performed using clustalw (http://www.genome.jp/tools/clustalw/). The following wgs data sets were used: Myotis lucifugus AAPE02 (Genome Coverage (GC): 7×; Sequencing Technology (ST): Sanger); Myotis brandtii ANKR01 (GC: 120×; ST: Illumina HiSeq 2000); Myotis davidii ALWT01 (GC: 110×; ST: Illumina HighSeq 2000); Eptesicus fuscus ALEH01 (GC: 84×; ST: Illumina Hi-Seq); Pteropus vampyrus ABRP02 (GC: 188×; ST: Illumina); Pteropus alecto ALWS01 (GC: 110×; ST: Illumina HighSeq 2000); Pteronotus parnellii AWGZ01 (GC: 17×; ST: Illumina HiSeq); Rhinolophus ferrumequinum AWHA01 (GC: 17×; ST: Illumina HiSeq); Megaderma lyra AWHB01 (GC: 18×; ST: Illumina HiSeq); Eidolon helvum AWHC01 (GC: 18×; ST: Illumina HiSeq); Miniopterus natalensis LDJU01 (GC: 77×; ST: Illumina HiSeq); Rousettus aegyptiacus LOCP02 (GC: 169.2×; Illumina HiSeq; PacBio).
The CEA gene family in bats is not well annotated; therefore, we adopted the nomenclature according to the one previously used for the CEA gene family of other mammals . Gene names and corresponding sequences are summarized in (Additional file 1).
Phylogenetic analyses based on nucleotide and amino acid sequences were conducted using MEGA5 and MEGA6. Sequence alignments were performed using “Muscle”. The maximum likelihood (ML) method with bootstrap testing (500 replicates) was applied for the construction of phylogenetic trees. To determine the selective pressure on the maintenance of the nucleotide sequences, the number of nonsynonymous nucleotide substitution per nonsynonymous site (dN) and the number of synonymous substitutions per synonymous site (dS) were determined for N domain exons. The dN/dS ratios were calculated after manual editing of sequence gaps or insertions guided by the amino acid sequences for all branches of the resulting phylogenetic trees using the Datamonkey web interface. The mean dN/dS ratios were calculated using the single likelihood ancestor counting (SLAC) algorithm. The synonymous nonsynonymous analysis program (SNAP; http://www.hiv.lanl.gov/content/sequence/SNAP/SNAP.html) allowed the calculation of cumulative average synonymous and nonsynonymous substitutions along coding regions of N domain exons from paralogous and orthologous genes. For the identification of N domain-wide episodic selection we used a branch-site unrestricted statistical test for episodic diversification (BUSTED) approach . For the detection of individual sites under positive selection we used the mixed effects model of evolution software (MEME)  after screening for recombination using the genetic algorithm for recombination detection (GARD) software .
M. myotis individuals, which could not survive in nature because of injuries, were used for isolation of tissues. Samples from immune organs including spleen, thymus, intestine and lymph nodes were stored in RNAlater at −80 °C until RNA extraction. Total RNA was extracted and contaminating DNA was removed by DNase I treatment using the RNeasy Mini Kit (Qiagen, Germany). The RNA concentration was determined with the NanoDrop 2000/2000c spectrophotometer (Thermo Fisher Scientific, USA), and RNA integrity was tested by measurement of the 28S/18S rRNA ratio using bioanalyzer Agilent2100 (Agilent Technologies Inc., USA). RNA from spleen, thymus, intestine and lymph nodes were pooled in a mass ratio of 1:1:1:1 and used for de novo transcriptome sequencing by Illumina Hi-SeqTM2000. RNA processing, cDNA library construction, sequencing and data processing were performed in the Beijing Genomics Institute (BGI), Shenzhen, China. For de novo assembly, raw reads were first filtered to remove adaptor sequences and reads with more than 5% unknown bases (N) and more than 20% low quality bases (bases with quality value ≤10), and then clean data were assembled using the short reads assembling program Trinity into non-redundant unigenes . Next, all of the unigenes were annotated by the best hits out of BLASTX alignments against protein databases of non-redundant proteins (NR) (http://www.ncbi.nlm.nih.gov), Swiss-Prot protein (http://www.uniprot.org/uniprot/), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway (http://www.genome.jp/kegg) and Cluster of Orthologous Groups (COG) (http://www.ncbi.nlm.nih.gov/COG) (E value <0.00001). Those without any match in above databases were further aligned by blastn to nucleotide databases (NT) (E value <0.00001). With NR annotations, GO functional annotations and classifications were obtained using the Blast2GO program  and the WEGO software , respectively.
This study was supported by the BMWi (KF2875802UL2), GIZ (Contract no. 81170269; Project No. 13.1432.7---001.00) and DFG (HE 6249/4–1) (to R.K.) and FAO (to S.M.)
Availability of data and materials
Nucleotide sequences from the N domains of newly described Myotis lucifugus genes which can be used as gene identifiers for search in data bases are provided (Additional file 1). Phylogenetic trees, sequence data and alignments used to produce the results were deposited to TreeBASE (https://www.treebase.org/). They are available using the following link http://purl.org/phylo/treebase/phylows/study/TB2:S21362. The RNA-seq datasets supporting the conclusions of this article are available upon request and will be available from the NCBI Gene Expression Omnibus database.
RK conceived the study, carried out most of the data analysis and drafted the manuscript. MM, JH, SM, contributed to data analysis and manuscript writing, XH, BK, provided critical reagents and RNA sequencing data, SM, WZ performed data mining and contributed to manuscript writing. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Ethical approval for all capturing and sampling was granted by the competent authorities in Germany and Czech Republic. The Czech Academy of Sciences Ethics Committee reviewed and approved the animal use protocol No. 169/2011 in compliance with Law No. 312/2008 on Protection of Animals Against Cruelty adopted by the Parliament of the Czech Republic. The capture and sampling of a M. myotis specimen in the Moravian Karst in November 2012 was in compliance with Law No. 114/1992 on Nature and Landscape Protection, and was based on permit 01662/MK/2012S/00775/MK/2012 issued by the Nature Conservation Agency of the Czech Republic.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Herzig CT, Waters RW, Baldwin CL, Telfer JC. Evolution of the CD163 family and its relationship to the bovine gamma delta T cell co-receptor WC1. BMC Evol Biol. 2010;10:181.View ArticlePubMedPubMed CentralGoogle Scholar
- Chuong EB, Tong W, Hoekstra HE. Maternal-fetal conflict: rapidly evolving proteins in the rodent placenta. Mol Biol Evol. 2010;27(6):1221–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Tchoupa AK, Schuhmacher T, Hauck CR. Signaling by epithelial members of the CEACAM family - mucosal docking sites for pathogenic bacteria. Cell Commun Signal. 2014;12:27.View ArticlePubMedPubMed CentralGoogle Scholar
- Singer BB, Opp L, Heinrich A, Schreiber F, Binding-Liermann R, Berrocal-Almanza LC, et al. Soluble CEACAM8 interacts with CEACAM1 inhibiting TLR2-triggered immune responses. PLoS One. 2014;9(4):e94106.View ArticlePubMedPubMed CentralGoogle Scholar
- Kammerer R, Hahn S, Singer BB, Luo JS, von Kleist S. Biliary glycoprotein (CD66a), a cell adhesion molecule of the immunoglobulin superfamily, on human lymphocytes: structure, expression and involvement in T cell activation. Eur J Immunol. 1998;28(11):3664–74.View ArticlePubMedGoogle Scholar
- Kammerer R, Stober D, Singer BB, Obrink B, Reimann J. Carcinoembryonic antigen-related cell adhesion molecule 1 on murine dendritic cells is a potent regulator of T cell stimulation. J Immunol. 2001;166(11):6537–44.View ArticlePubMedGoogle Scholar
- Huang YH, Zhu C, Kondo Y, Anderson AC, Gandhi A, Russell A, et al. CEACAM1 Regulates TIM-3-mediated tolerance and exhaustion. Nature. 2015;517(7534):386–90.View ArticlePubMedGoogle Scholar
- Moore T, Dveksler GS. Pregnancy-specific glycoproteins: complex gene families regulating maternal-fetal interactions. The International journal of developmental biology. 2014;58(2–4):273–80.View ArticlePubMedGoogle Scholar
- Kammerer R, Popp T, Hartle S, Singer BB, Zimmermann W. Species-specific evolution of immune receptor tyrosine based activation motif-containing CEACAM1-related immune receptors in the dog. BMC Evol Biol. 2007;7:196.View ArticlePubMedPubMed CentralGoogle Scholar
- Kammerer R, Popp T, Singer BB, Schlender J, Zimmermann W. Identification of allelic variants of the bovine immune regulatory molecule CEACAM1 implies a pathogen-driven evolution. Gene. 2004;339:99–109.View ArticlePubMedGoogle Scholar
- Kammerer R, Zimmermann W. Coevolution of activating and inhibitory receptors within mammalian carcinoembryonic antigen families. BMC Biol. 2010;8:12.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang G, Cowled C, Shi Z, Huang Z, Bishop-Lilly KA, Fang X, et al. Comparative analysis of bat genomes provides insight into the evolution of flight and immunity. Science. 2013;339(6118):456–60.View ArticlePubMedGoogle Scholar
- Teeling EC. Hear, hear: the convergent evolution of echolocation in bats? Trends Ecol Evol. 2009;24(7):351–4.View ArticlePubMedGoogle Scholar
- Teeling EC, Springer MS, Madsen O, Bates P, O'Brien SJ, Murphy WJ. A molecular phylogeny for bats illuminates biogeography and the fossil record. Science. 2005;307(5709):580–4.View ArticlePubMedGoogle Scholar
- Baker ML, Schountz T, Wang LF. Antiviral immune responses of bats: a review. Zoonoses Public Health. 2013;60(1):104–16.View ArticlePubMedGoogle Scholar
- Schountz T. Immunology of bats and their viruses: challenges and opportunities. Viruses-Basel. 2014;6(12):4880–U4410.View ArticleGoogle Scholar
- Grosser O. Vergleichende Anatomie und Entwicklungsgeschichte der Eihäute und der Placenta mit besonderer Berücksichtigung des Menschen. Wien: Verlag Braumüller; 1909.Google Scholar
- Platt RN 2nd, Mangum SF, Ray DA. Pinpointing the vesper bat transposon revolution using the Miniopterus Natalensis genome. Mob DNA. 2016;7:12.View ArticlePubMedPubMed CentralGoogle Scholar
- Agnarsson I, Zambrana-Torrelio CM, Flores-Saldana NP, May-Collado LJ. A time-calibrated species-level phylogeny of bats (Chiroptera, Mammalia). PLoS currents. 2011;3:RRN1212.View ArticlePubMedPubMed CentralGoogle Scholar
- Kawai K, Nikaido M, Harada M, Matsumura S, Lin LK, Wu Y, et al. The status of the Japanese and east Asian bats of the genus Myotis (Vespertilionidae) based on mitochondrial sequences. Mol Phylogenet Evol. 2003;28(2):297–307.View ArticlePubMedGoogle Scholar
- Aleksic D, Blaschke L, Missbach S, Hanske J, Weiss W. Handler J, et al. Reproduction: Convergent evolution of pregnancy-specific glycoproteins in the human and horse; 2016.Google Scholar
- Pavlopoulou A, Scorilas A. A comprehensive phylogenetic and structural analysis of the carcinoembryonic antigen (CEA) gene family. Genome Biol Evol. 2014;6(6):1314–26.View ArticlePubMedPubMed CentralGoogle Scholar
- Lei M, Dong D. Phylogenomic analyses of bat subordinal relationships based on transcriptome data. Sci Rep. 2016;6:27726.View ArticlePubMedPubMed CentralGoogle Scholar
- Miller-Butterworth CM, Murphy WJ, O'Brien SJ, Jacobs DS, Springer MS, Teeling EC. A family matter: conclusive resolution of the taxonomic position of the long-fingered bats, miniopterus. Mol Biol Evol. 2007;24(7):1553–61.View ArticlePubMedGoogle Scholar
- Zimmermann W, Kammerer R. Coevolution of paired receptors in Xenopus carcinoembryonic antigen-related cell adhesion molecule families suggests appropriation as pathogen receptors. BMC Genomics. 2016;17(1):928.View ArticlePubMedPubMed CentralGoogle Scholar
- Papenfuss AT, Baker ML, Feng ZP, Tachedjian M, Crameri G, Cowled C, et al. The immune gene repertoire of an important viral reservoir, the Australian black flying fox. BMC Genomics. 2012;13:261.View ArticlePubMedPubMed CentralGoogle Scholar
- Carter AM, Mess A. Evolution of the placenta and associated reproductive characters in bats. J Exp Zool B Mol Dev Evol. 2008;310(5):428–49.View ArticlePubMedGoogle Scholar
- Gopalakrishna A, Karim KB. Fetal membranes and placentation in Chiroptera. J Reprod Fertil. 1979;56(1):417–29.View ArticlePubMedGoogle Scholar
- Muenzner P, Bachmann V, Zimmermann W, Hentschel J, Hauck CR. Human-restricted bacterial pathogens block shedding of epithelial cells by stimulating integrin activation. Science. 2010;329(5996):1197–201.View ArticlePubMedGoogle Scholar
- Slevogt H, Zabel S, Opitz B, Hocke A, Eitel J, N'Guessan PD, et al. CEACAM1 Inhibits toll-like receptor 2-triggered antibacterial responses of human pulmonary epithelial cells. Nat Immunol. 2008;9(11):1270–8.View ArticlePubMedGoogle Scholar
- Sintsova A, Wong H, MacDonald KS, Kaul R, Virji M, Gray-Owen SD. Selection for a CEACAM receptor-specific binding phenotype during Neisseria gonorrhoeae infection of the human genital tract. Infect Immun. 2015;83(4):1372–83.View ArticlePubMedPubMed CentralGoogle Scholar
- Javaheri A, Kruse T, Moonens K, Mejias-Luque R, Debraekeleer A, Asche CI, et al. Helicobacter pylori adhesin HopQ engages in a virulence-enhancing interaction with human CEACAMs. Nature microbiology. 2016;2:16189.View ArticlePubMedGoogle Scholar
- Koniger V, Holsten L, Harrison U, Busch B, Loell E, Zhao Q, et al. Helicobacter pylori exploits human CEACAMs via HopQ for adherence and translocation of CagA. Nature microbiology. 2016;2:16188.View ArticlePubMedGoogle Scholar
- Hill DJ, Toleman MA, Evans DJ, Villullas S, Van Alphen L, Virji M. The variable P5 proteins of typeable and non-typeable Haemophilus influenzae target human CEACAM1. Mol Microbiol. 2001;39(4):850–62.View ArticlePubMedGoogle Scholar
- Berger CN, Billker O, Meyer TF, Servin AL, Kansau I. Differential recognition of members of the carcinoembryonic antigen family by Afa/Dr adhesins of diffusely adhering Escherichia Coli (Afa/Dr DAEC). Mol Microbiol. 2004;52(4):963–83.View ArticlePubMedGoogle Scholar
- Dveksler GS, Pensiero MN, Cardellichio CB, Williams RK, Jiang GS, Holmes KV, et al. Cloning of the mouse hepatitis virus (MHV) receptor: expression in human and hamster cell lines confers susceptibility to MHV. J Virol. 1991;65(12):6881–91.PubMedPubMed CentralGoogle Scholar
- Taguchi F, Hirai-Yuki A. Mouse hepatitis virus receptor as a determinant of the mouse susceptibility to MHV infection. Front Microbiol. 2012;3:68.View ArticlePubMedPubMed CentralGoogle Scholar
- Brandao PE, Scheffer K, Villarreal LY, Achkar S, Oliveira Rde N, Fahl Wde O, et al. A coronavirus detected in the vampire bat Desmodus Rotundus. The Brazilian journal of infectious diseases : an official publication of the Brazilian Society of Infectious Diseases. 2008;12(6):466–8.View ArticleGoogle Scholar
- Murrell B, Weaver S, Smith MD, Wertheim JO, Murrell S, Aylward A, et al. Gene-wide identification of episodic selection. Mol Biol Evol. 2015;32(5):1365–71.View ArticlePubMedPubMed CentralGoogle Scholar
- Murrell B, Wertheim JO, Moola S, Weighill T, Scheffler K, Kosakovsky Pond SL. Detecting individual sites subject to episodic diversifying selection. PLoS Genet. 2012;8(7):e1002764.View ArticlePubMedPubMed CentralGoogle Scholar
- Kosakovsky Pond SL, Posada D, Gravenor MB, Woelk CH, Frost SD. GARD: a genetic algorithm for recombination detection. Bioinformatics. 2006;22(24):3096–8.View ArticlePubMedGoogle Scholar
- Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644–52.View ArticlePubMedPubMed CentralGoogle Scholar
- Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M. Blast2GO: A universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6.View ArticlePubMedGoogle Scholar
- Ye J, Fang L, Zheng H, Zhang Y, Chen J, Zhang Z, et al. WEGO: a web tool for plotting GO annotations. Nucleic acids research 2006, 34(Web Server issue): W293–297.Google Scholar