Patterns of positive selection in six Mammalian genomes

Genome-wide scans for positively selected genes (PSGs) in mammals have provided insight into the dynamics of genome evolution, the genetic basis of differences between species, and the functions of individual genes. However, previous scans have been limited in power and accuracy owing to small numbers of available genomes. Here we present the most comprehensive examination of mammalian PSGs to date, using the six high-coverage genome assemblies now available for eutherian mammals. The increased phylogenetic depth of this dataset results in substantially improved statistical power, and permits several new lineage- and clade-specific tests to be applied. Of approximately 16,500 human genes with high-confidence orthologs in at least two other species, 400 genes showed significant evidence of positive selection (FDR<0.05), according to a standard likelihood ratio test. An additional 144 genes showed evidence of positive selection on particular lineages or clades. As in previous studies, the identified PSGs were enriched for roles in defense/immunity, chemosensory perception, and reproduction, but enrichments were also evident for more specific functions, such as complement-mediated immunity and taste perception. Several pathways were strongly enriched for PSGs, suggesting possible co-evolution of interacting genes. A novel Bayesian analysis of the possible "selection histories" of each gene indicated that most PSGs have switched multiple times between positive selection and nonselection, suggesting that positive selection is often episodic. A detailed analysis of Affymetrix exon array data indicated that PSGs are expressed at significantly lower levels, and in a more tissue-specific manner, than non-PSGs. Genes that are specifically expressed in the spleen, testes, liver, and breast are significantly enriched for PSGs, but no evidence was found for an enrichment for PSGs among brain-specific genes. This study provides additional evidence for widespread positive selection in mammalian evolution and new genome-wide insights into the functional implications of positive selection.     

Cytochromeb gene of marine mammals: Phylogeny and evolution

The DNA sequences of the mitochondrial cytochromeb gene of marine mammals (Cetacea, Pinnipedia, Sirenia) were compared with cytochromeb genes of terrestrial mammals including the semiaquatic hippopotamus. The comparison included 28 sequences, representing 22 families and 10 orders. The dugong (order Sirenia) sequence associated with that of the elephant, supporting the Tethytheria clade. The fin whale and dolphin (order Cetacea) sequences are more closely related to those of the artiodactyls, and the comparison suggests that the hippopotamus may be the extant artiodactyl species that is most closely related to the cetaceans. The seal sequence may be more closely related to those of artiodactyls, cetaceans, and perissodactyls than to tethytheres, rodents, lagomorphs, or primates. The cytochromeb proteins of mammals do not evolve at a uniform rate. Human and elephant cytochromeb amino acid sequences were found to evolve the most rapidly, while those of myomorph rodents evolved slowest. The cytochromeb of marine mammals evolves at an intermediate rate. The pattern of amino acid substitutions in marine mammals is similar to that of terrestrial mammals.     

Positive selection at the ASPM gene coincides with brain size enlargements in cetaceans

The enlargement of cetacean brain size represents an enigmatic event in mammalian evolution, yet its genetic basis remains poorly explored. One candidate gene associated with brain size evolution is the abnormal spindle-like microcephaly associated (ASPM), as mutations in this gene cause severe reductions in the cortical size of humans. Here, we investigated the ASPM gene in representative cetacean lineages and previously published sequences from other mammals to test whether the expansion of the cetacean brain matched adaptive ASPM evolution patterns. Our analyses yielded significant evidence of positive selection on the ASPM gene during cetacean evolution, especially for the Odontoceti and Delphinoidea lineages. These molecular patterns were associated with two major events of relative brain size enlargement in odontocetes and delphinoids. It is of particular interest to find that positive selection was restricted to cetaceans and primates, two distant lineages both characterized by a massive expansion of brain size. This result is suggestive of convergent molecular evolution, although no site-specific convergence at the amino acid level was found.     

Increased rate of hair keratin gene loss in the cetacean lineage

Background

Hair represents an evolutionary innovation that appeared early on mammalian evolutionary history, and presumably contributed significantly to the rapid radiation of the group. An interesting event in hair evolution has been its secondary loss in some mammalian groups, such as cetaceans, whose hairless phenotype appears to be an adaptive response to better meet the environmental conditions. To determine whether different repertoire of keratin genes among mammals can potentially explain the phenotypic hair features of different lineages, we characterized the type I and II clusters of alpha keratins from eight mammalian species, including the hairless dolphin and minke whale representing the order Cetacea.

Results

We combined the available genomic information with phylogenetic analysis to conduct a comprehensive analysis of the evolutionary patterns of keratin gene clusters. We found that both type I and II gene clusters are fairly conserved among the terrestrial mammals included in this study, with lineage specific gene duplication and gene loss. Nevertheless, there is also evidence for an increased rate of pseudogenization in the cetacean lineage when compared to their terrestrial relatives, especially among the hair type keratins.

Conclusions

Here we present a comprehensive characterization of alpha-keratin genes among mammals and elucidate the mechanisms involved in the evolution of this gene family. We identified lineage-specific gene duplications and gene loss among the Laurasiatherian and Euarchontoglires species included in the study. Interestingly, cetaceans present an increased loss of hair-type keratin genes when compared to other terrestrial mammals. As suggested by the ‘less-is-more’ hypothesis, we do not rule out the possibility that the gene loss of hair-type keratin genes in these species might be associated to the hairless phenotype and could have been adaptive in response to new selective pressures imposed by the colonization of a new habitat. Our study provides support for the idea that pseudogenes are not simply ‘genomic fossils’ but instead have adaptive roles during the evolutionary process.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-869) contains supplementary material, which is available to authorized users.     

Adaptive molecular evolution of HINTW,a female-specific gene in birds

It is well established that many genes on the male-specific Y chromosome of organisms such as mammals are involved in male reproduction and may evolve rapidly because of positive selection on male reproductive traits. In contrast, very little is known about the function and evolution of W-linked genes restricted to the female genome of organisms with female heterogamety. For birds (males ZZ, females ZW), only one W-linked gene (HINTW) is sufficiently different from its Z-linked homolog to indicate a female-specific function. Here, we report that HINTW shows evidence of adaptive molecular evolution, implying strong positive selection for new functional properties in female birds. Moreover, because HINTW is expressed in the gonads of female birds just before sexual differentiation and is thus a candidate for sex determination, it suggests adaptive evolution related to female development. This provides the first example of Darwinian evolution of a gene restricted to the female genome of any organism. Given that HINTW exists in multiple copies on W, similar to some testis-specific genes amplified on mammalian Y, avian HINTW may thus potentially represent a female parallel to the organization and evolution of Y chromosome genes involved in male reproduction and development.     

A genomics approach to the detection of positive selection in cattle: adaptive evolution of the T-cell and natural killer cell-surface protein CD2

The detection of adaptive evolution at the molecular level is of interest not only as an insight into the process of evolution but also because of its functional implications for genes of interest. Here, we present the first genomics approach to detecting positive selection operating on the Bos taurus lineage, an important domestic species. This analysis led to the identification of the T-cell and natural killer (NK) cell receptor cluster of differentiation 2 (CD2) as having a strong signal of selection. Further detailed investigation of CD2 revealed that this gene was subject to positive selection during the evolution of a number of mammalian lineages. Moreover, we show that selection has operated primarily on the extracellular domain of CD2 and discuss the implications of this for an important regulator of the adaptive immune response.     

Discordant evolution of the adjacent antiretroviral genes TRIM22 and TRIM5 in mammals

TRIM5alpha provides a cytoplasmic block to retroviral infection, and orthologs encoded by some primates are active against HIV. Here, we present an evolutionary comparison of the TRIM5 gene to its closest human paralogs: TRIM22, TRIM34, and TRIM6. We show that TRIM5 and TRIM22 have a dynamic history of gene expansion and loss during the evolution of mammals. The cow genome contains an expanded cluster of TRIM5 genes and no TRIM22 gene, while the dog genome encodes TRIM22 but has lost TRIM5. In contrast, TRIM6 and TRIM34 have been strictly preserved as single gene orthologs in human, dog, and cow. A more focused analysis of primates reveals that, while TRIM6 and TRIM34 have evolved under purifying selection, TRIM22 has evolved under positive selection as was previously observed for TRIM5. Based on TRIM22 sequences obtained from 27 primate genomes, we find that the positive selection of TRIM22 has occurred episodically for approximately 23 million years, perhaps reflecting the changing pathogenic landscape. However, we find that the evolutionary episodes of positive selection that have acted on TRIM5 and TRIM22 are mutually exclusive, with generally only one of these genes being positively selected in any given primate lineage. We interpret this to mean that the positive selection of one gene has constrained the adaptive flexibility of its neighbor, probably due to genetic linkage. Finally, we find a striking congruence in the positions of amino acid residues found to be under positive selection in both TRIM5alpha and TRIM22, which in both proteins fall predominantly in the beta2-beta3 surface loop of the B30.2 domain. Astonishingly, this same loop is under positive selection in the multiple cow TRIM5 genes as well, indicating that this small structural loop may be a viral recognition motif spanning a hundred million years of mammalian evolution.     

Use of a mammalian interspersed repetitive (MIR) element in the coding and processing sequences of mammalian genes.

The mammalian interspersed repetitive (MIR) element was amplified in mammals 130 million years ago. The MIR element is at least 260 bp in length and is found in approximately 105 copies in the mammalian genome. We analyzed copies of the MIR element in the DNA of various mammals to determine its relationship to the structure and function of genes, in an attempt to identify specific uses of the MIR element within the mammalian genome. We found that alternative splicing within the acetylcholine receptor gene in humans takes place within the MIR element and results in the incorporation of part of the MIR element into the coding sequence of this gene. Furthermore, the polyadenylation signal (AATAAA) at the 3' end of four different mammalian genes is derived from the MIR element. These uses of the MIR element suggest that other regulatory sequences found within the mammalian genome originated from ancient transposable elements, many of which may no longer be recognizable.     

Phocid seal leptin: tertiary structure and hydrophobic receptor binding site preservation during distinct leptin gene evolution

The cytokine hormone leptin is a key signalling molecule in many pathways that control physiological functions. Although leptin demonstrates structural conservation in mammals, there is evidence of positive selection in primates, lagomorphs and chiropterans. We previously reported that the leptin genes of the grey and harbour seals (phocids) have significantly diverged from other mammals. Therefore we further investigated the diversification of leptin in phocids, other marine mammals and terrestrial taxa by sequencing the leptin genes of representative species. Phylogenetic reconstruction revealed that leptin diversification was pronounced within the phocid seals with a high dN/dS ratio of 2.8, indicating positive selection. We found significant evidence of positive selection along the branch leading to the phocids, within the phocid clade, but not over the dataset as a whole. Structural predictions indicate that the individual residues under selection are away from the leptin receptor (LEPR) binding site. Predictions of the surface electrostatic potential indicate that phocid seal leptin is notably different to other mammalian leptins, including the otariids. Cloning the grey seal leptin binding domain of LEPR confirmed that this was structurally conserved. These data, viewed in toto, support a hypothesis that phocid leptin divergence is unlikely to have arisen by random mutation. Based upon these phylogenetic and structural assessments, and considering the comparative physiology and varying life histories among species, we postulate that the unique phocid diving behaviour has produced this selection pressure. The Phocidae includes some of the deepest diving species, yet have the least modified lung structure to cope with pressure and volume changes experienced at depth. Therefore, greater surfactant production is required to facilitate rapid lung re-inflation upon surfacing, while maintaining patent airways. We suggest that this additional surfactant requirement is met by the leptin pulmonary surfactant production pathway which normally appears only to function in the mammalian foetus.     

The voltage-gated potassium channel subfamily KQT member 4 (KCNQ4) displays parallel evolution in echolocating bats

Bats are the only mammals that use highly developed laryngeal echolocation, a sensory mechanism based on the ability to emit laryngeal sounds and interpret the returning echoes to identify objects. Although this capability allows bats to orientate and hunt in complete darkness, endowing them with great survival advantages, the genetic bases underlying the evolution of bat echolocation are still largely unknown. Echolocation requires high-frequency hearing that in mammals is largely dependent on somatic electromotility of outer hair cells. Then, understanding the molecular evolution of outer hair cell genes might help to unravel the evolutionary history of echolocation. In this work, we analyzed the molecular evolution of two key outer hair cell genes: the voltage-gated potassium channel gene KCNQ4 and CHRNA10, the gene encoding the α10 nicotinic acetylcholine receptor subunit. We reconstructed the phylogeny of bats based on KCNQ4 and CHRNA10 protein and nucleotide sequences. A phylogenetic tree built using KCNQ4 amino acid sequences showed that two paraphyletic clades of laryngeal echolocating bats grouped together, with eight shared substitutions among particular lineages. In addition, our analyses indicated that two of these parallel substitutions, M388I and P406S, were probably fixed under positive selection and could have had a strong functional impact on KCNQ4. Moreover, our results indicated that KCNQ4 evolved under positive selection in the ancestral lineage leading to mammals, suggesting that this gene might have been important for the evolution of mammalian hearing. On the other hand, we found that CHRNA10, a gene that evolved adaptively in the mammalian lineage, was under strong purifying selection in bats. Thus, the CHRNA10 amino acid tree did not show echolocating bat monophyly and reproduced the bat species tree. These results suggest that only a subset of hearing genes could underlie the evolution of echolocation. The present work continues to delineate the genetic bases of echolocation and ultrasonic hearing in bats.     

Evolution of cd59 gene in mammals

The CD59-coding sequences were obtained from 5 mammals by PCR and BLAST, and combined with the available sequences in GenBank, the nucleotide substitution rates of mammalian cd59 were calcu- lated. Results of synonymous and nonsynonymous substitution rates revealed that cd59 experienced negative selection in mammals overall. Four sites experiencing positive selection were found by using "site-specific" model in PAML software. These sites were distributed on the molecular surface, of which 2 sites located in the key functional domain. Furthermore, "branch-site-specific" model detected 1 positive site in cd59a and cd59b lineages which underwent accelerated evolution caused by positive selection after gene duplication in mouse.     

A Mitochondrial Genome Sequence of the Tibetan Antelope （ Pantholops hodgsonii ）

To investigate genetic mechanisms of high altitude adaptations of native mammals on the Tibetan Plateau, we compared mitochondrial sequences of the endangered Pantholops hodgsonii with its lowland distant relatives Ovis ames and Capra hircus, as well as other mammals. The complete mitochondrial genome of P. hodgsonii （16,498 bp） revealed a similar gene order as of other mammals. Because of tandem duplications, the control region of P. hodgsonii mitochondrial genome is shorter than those of O. ames and C. hircus, but longer than those of Bos species. Phylogenetic analysis based on alignments of the entire cytochrome b genes suggested that P. hodgsonii is more closely related to O. ames and C. hircus, rather than to species of the Antilopinae subfamily. The estimated divergence time between P. hodgsonii and O. ames is about 2.25 million years ago. Eutther analysis on natural selection indicated that the COXI （cytochrome c oxidase subunit I） gene was under positive selection in P. hodgsonii and Bos grunniens. Considering the same climates and environments shared by these two mammalian species, we proposed that the mitochondrial COXI gene is probably relevant for these native mammals to adapt the high altitude environment unique to the Tibetan Plateau.     

Genomic V exons from whole genome shotgun data in reptiles

Reptiles and mammals diverged over 300 million years ago, creating two parallel evolutionary lineages amongst terrestrial vertebrates. In reptiles, two main evolutionary lines emerged: one gave rise to Squamata, while the other gave rise to Testudines, Crocodylia, and Aves. In this study, we determined the genomic variable (V) exons from whole genome shotgun sequencing (WGS) data in reptiles corresponding to the three main immunoglobulin (IG) loci and the four main T cell receptor (TR) loci. We show that Squamata lack the TRG and TRD genes, and snakes lack the IGKV genes. In representative species of Testudines and Crocodylia, the seven major IG and TR loci are maintained. As in mammals, genes of the IG loci can be grouped into well-defined IMGT clans through a multi-species phylogenetic analysis. We show that the reptilian IGHV and IGLV genes are distributed amongst the established mammalian clans, while their IGKV genes are found within a single clan, nearly exclusive from the mammalian sequences. The reptilian and mammalian TRAV genes cluster into six common evolutionary clades (since IMGT clans have not been defined for TR). In contrast, the reptilian TRBV genes cluster into three clades, which have few mammalian members. In this locus, the V exon sequences from mammals appear to have undergone different evolutionary diversification processes that occurred outside these shared reptilian clans. These sequences can be obtained in a freely available public repository (http://vgenerepertoire.org).     

Evolution of <Emphasis Type="Italic">cd59</Emphasis> gene in mammals

The CD59-coding sequences were obtained from 5 mammals by PCR and BLAST, and combined with the available sequences in GenBank, the nucleotide substitution rates of mammalian cd59 were calculated. Results of synonymous and nonsynonymous substitution rates revealed that cd59 experienced negative selection in mammals overall. Four sites experiencing positive selection were found by using “site-specific” model in PAML software. These sites were distributed on the molecular surface, of which 2 sites located in the key functional domain. Furthermore, “branch-site-specific” model detected 1 positive site in cd59a and cd59b lineages which underwent accelerated evolution caused by positive selection after gene duplication in mouse.     

Phylogenetic relationships of artiodactyls and cetaceans as deduced from the comparison of cytochrome b and 12S rRNA mitochondrial sequences

A data set of complete mitochondrial cytochrome b and 12S rDNA sequences is presented here for 17 representatives of Artiodactyla and Cetacea, together with potential outgroups (two Perissodactyla, two Carnivora, two Tethytheria, four Rodentia, and two Marsupialia). We include seven sequences not previously published from Hippopotamidae (Ancodonta) and Camelidae (Tylopoda), yielding a total of nearly 2.1 kb for both genes combined. Distance and parsimony analyses of each gene indicate that 11 clades are well supported, including the artiodactyl taxa Pecora, Ruminantia (with low 12S rRNA support), Tylopoda, Suina, and Ancodonta, as well as Cetacea, Perissodactyla, Carnivora, Tethytheria, Muridae, and Caviomorpha. Neither the cytochrome b nor the 12S rDNA genes resolve the relationships between these major clades. The combined analysis of the two genes suggests a monophyletic Cetacea +Artiodactyla clade (defined as "Cetartiodactyla"), whereas Perissodactyla, Carnivora, and Tethytheria fall outside this clade. Perissodactyla could represent the sister taxon of Cetartiodactyla, as deduced from resampling studies among outgroup lineages. Cetartiodactyla includes five major lineages: Ruminantia, Tylopoda, Suina, Ancodonta, and Cetacea, among which the phylogenetic relationships are not resolved. Thus, Suiformes do not appear to be monophyletic, justifying their split into the Suina and Ancodonta infraorders. An association between Cetacea and Hippopotamidae is supported by the cytochrome b gene but not by the 12S rRNA gene. Calculation of divergence dates suggests that the Cetartiodactyla could have diverged from other Ferungulata about 60 MYA.      

Molecular evolution of candidate sour taste receptor gene PKD1L3 in mammals

The PKD1L3 gene encodes an ion channel protein that can interact with the PKD2L1 protein to form a candidate sour taste receptor. In the present study, we have analyzed the evolutionary patterns of PKD1L3 genes from 10 mammalian species. The results showed that PKD1L3 genes have evolved under a dominant purifying selection force. However, for some branches and sites, PKD1L3 genes were detected to have been operated by positive selection. Moreover, some of these positive evolutionary sites are likely to participate in acid stimulus recognition. In rodents, PKD1L3 genes evolved more rapidly than other mammalian lineages. Combined with other functional research reports, our results suggest that rodents may not be the most appropriate model for functional research on the PKD1L3 gene.     

An ancient history of gene duplications, fusions and losses in the evolution of APOBEC3 mutators in mammals

ABSTRACT: BACKGROUND: The APOBEC3 (A3) genes play a key role in innate antiviral defense in mammals by introducing directed mutations in the DNA. The human genome encodes for seven A3 genes, with multiple splice alternatives. Different A3 proteins display different substrate specificity, but the very basic question on how discerning self from non-self still remains unresolved. Further, the expression of A3 activity/ies shapes the way both viral and host genomes evolve. RESULTS: We present here a detailed temporal analysis of the origin and expansion of the A3 repertoire in mammals. Our data support an evolutionary scenario where the genome of the mammalian ancestor encoded for at least one ancestral A3 gene, and where the genome of the ancestor of placental mammals (and possibly of the ancestor of all mammals) already encoded for an A3Z1-A3Z2-A3Z3 arrangement. Duplication events of the A3 genes have occurred independently in different lineages: humans, cats and horses. In all of them, gene duplication has resulted in changes in enzyme activity and/or substrate specificity, in a paradigmatic example of convergent adaptive evolution at the genomic level. Finally, our results show that evolutionary rates for the three A3Z1, A3Z2 and A3Z3 motifs have significantly decreased in the last 100 Mya. The analysis constitutes a textbook example of the evolution of a gene locus by duplication and sub/neofunctionalization in the context of virus-host arms race. CONCLUSIONS: Our results provide a time framework for identifying ancestral and derived genomic arrangements in the APOBEC loci, and to date the expansion of this gene family for different lineages through time, as a response to changes in viral/retroviral/retrotransposon pressure.     

Contrasting modes of evolution between vertebrate sweet/umami receptor genes and bitter receptor genes

Taste reception is fundamental to diet selection in many animals. The genetic basis underlying the evolution and diversity of taste reception, however, is not well understood. Recent discoveries of T1R sweet/umami receptor genes and T2R bitter receptor genes in humans and mice provided an opportunity to address this question. Here, we report the identification of 20 putatively functional T1R genes and 167 T2R genes from the genome sequences of nine vertebrates, including three fishes, one amphibian, one bird, and four mammals. Our comparative genomic analysis shows that orthologous T1R sequences are relatively conserved in evolution and that the T1R gene repertoire remains virtually constant in size across most vertebrates, except for the loss of the T1R2 sweet receptor gene in the sweet-insensitive chicken and the absence of all T1R genes in the tongueless western clawed frog. In contrast, orthologous T2R sequences are more variable, and the T2R repertoire diverges tremendously among species, from only three functional genes in the chicken to 49 in the frog. These evolutionary patterns suggest the relative constancy in the number and type of sweet and umami tastants encountered by various vertebrates or low binding specificities of T1Rs but a large variation in the number and type of bitter compounds detected by different species. Although the rate of gene duplication is much lower in T1Rs than in T2Rs, signals of positive selection are detected during the functional divergences of paralogous T1Rs, as was previously found among paralogous T2Rs. Thus, functional divergence and specialization of taste receptors generally occurred via adaptive evolution.     

中国海兽研究进展

中国是海兽大国,共有46种海兽在中国水域有分布,包括鲸类38种（Cetartiodactyla, Cetacea）,海牛目儒艮科（Sirenia, Dugongidae）1种,食肉目鳍足类（Carnivora, Pinnipedia）5种和水獭亚科（Lutrinae）2种。从20世纪20年代开始,中国科学家围绕中国水域海兽的生物学与保护,通过多学科手段开展了一系列研究,特别是在白暨豚(Lipotes vexillifer)、江豚(Neophocaena spp.)和中华白海豚(Sousa chinensis)的生态学与保护、鲸类次生性水生适应的演化历史与适应机制、种群遗传多样性等方面,取得了一系列重要的创新进展,不仅促进了对中国海兽的认识,也推动中国海兽的保护工作取得成效。本文从海兽多样性、生态与保护生物学、形态与解剖学、饲养与繁殖生物学、声行为学、遗传与演化生物学以及基因组学等7个方面,较为系统地综述了中国海兽的研究进展。同时指出,海兽次生性水生适应历史及其演化机制,将是今后值得特别关注和需要深入揭示的重要科学问题;并且通过多学科手段揭示海兽的濒危机制与特点以提出促进其资源保护与管理措施,也具有十分重要的意义。