The PI/GLO-like locus in orchids: duplication and purifying selection at synonymous sites within Orchidinae (Orchidaceae)

Positive selection and relaxation of purifying constraints after duplication events have driven the functional diversification of gene families involved in development. One example of this occurred within the plant MADS-box genes. The evolution of the orchid flower was driven by duplication events followed by sub- and neo-functionalization of class B DEF-like MADS-box genes, which are present at three to four copies in the orchid genome. In contrast, the orchid PI/GLO-like class B MADS-box genes have been reported thus far as single-copy loci, with the only exception of Habenaria radiata.We isolated a novel PI/GLO-like gene (OrcPI2) in Orchis italica, which is different than the previously characterized OrcPI locus. The presence of two functional paralogs of PI/GLO-like genes in orchids is detectable only within the tribe Orchidinae. Evolutionary analyses revealed an apparent relaxation of purifying selection acting on the two PI/GLO-like paralogs of the Orchidinae when compared to the single-copy PI/GLO-like genes found in other orchid species. Furthermore, by measuring dN/dS (ω) ratios, we show that a high percentage of sites between the two PI/GLO-like paralogs have different evolutionary pressures. Interestingly, the apparent relaxation of selective constraints on the two PI/GLO-like paralogs is due to strong purifying selection at synonymous sites rather than to a high value of nonsynonymous substitution rate. This peculiar evolutionary pattern might be related to molecular processes such as mRNA folding and/or translational efficiency control. These processes could potentially be involved in or predate the functional diversification of the two PI/GLO-like paralogs within Orchidinae.     

Positive selection, relaxation, and acceleration in the evolution of the human and chimp genome

For years evolutionary biologists have been interested in searching for the genetic bases underlying humanness. Recent efforts at a large or a complete genomic scale have been conducted to search for positively selected genes in human and in chimp. However, recently developed methods allowing for a more sensitive and controlled approach in the detection of positive selection can be employed. Here, using 13,198 genes, we have deduced the sets of genes involved in rate acceleration, positive selection, and relaxation of selective constraints in human, in chimp, and in their ancestral lineage since the divergence from murids. Significant deviations from the strict molecular clock were observed in 469 human and in 651 chimp genes. The more stringent branch-site test of positive selection detected 108 human and 577 chimp positively selected genes. An important proportion of the positively selected genes did not show a significant acceleration in rates, and similarly, many of the accelerated genes did not show significant signals of positive selection. Functional differentiation of genes under rate acceleration, positive selection, and relaxation was not statistically significant between human and chimp with the exception of terms related to G-protein coupled receptors and sensory perception. Both of these were over-represented under relaxation in human in relation to chimp. Comparing differences between derived and ancestral lineages, a more conspicuous change in trends seems to have favored positive selection in the human lineage. Since most of the positively selected genes are different under the same functional categories between these species, we suggest that the individual roles of the alternative positively selected genes may be an important factor underlying biological differences between these species.     

An Evolutionary Analysis of Antigen Processing and Presentation across Different Timescales Reveals Pervasive Selection

The antigenic repertoire presented by MHC molecules is generated by the antigen processing and presentation (APP) pathway. We analyzed the evolutionary history of 45 genes involved in APP at the inter- and intra-species level. Results showed that 11 genes evolved adaptively in mammals. Several positively selected sites involve positions of fundamental importance to the protein function (e.g. the TAP1 peptide-binding domains, the sugar binding interface of langerin, and the CD1D trafficking signal region). In CYBB, all selected sites cluster in two loops protruding into the endosomal lumen; analysis of missense mutations responsible for chronic granulomatous disease (CGD) showed the action of different selective forces on the very same gene region, as most CGD substitutions involve aminoacid positions that are conserved in all mammals. As for ERAP2, different computational methods indicated that positive selection has driven the recurrent appearance of protein-destabilizing variants during mammalian evolution. Application of a population-genetics phylogenetics approach showed that purifying selection represented a major force acting on some APP components (e.g. immunoproteasome subunits and chaperones) and allowed identification of positive selection events in the human lineage.We also investigated the evolutionary history of APP genes in human populations by developing a new approach that uses several different tests to identify the selection target, and that integrates low-coverage whole-genome sequencing data with Sanger sequencing. This analysis revealed that 9 APP genes underwent local adaptation in human populations. Most positive selection targets are located within noncoding regions with regulatory function in myeloid cells or act as expression quantitative trait loci. Conversely, balancing selection targeted nonsynonymous variants in TAP1 and CD207 (langerin). Finally, we suggest that selected variants in PSMB10 and CD207 contribute to human phenotypes. Thus, we used evolutionary information to generate experimentally-testable hypotheses and to provide a list of sites to prioritize in follow-up analyses.     

Hotspots of Biased Nucleotide Substitutions in Human Genes

Genes that have experienced accelerated evolutionary rates on the human lineage during recent evolution are candidates for involvement in human-specific adaptations. To determine the forces that cause increased evolutionary rates in certain genes, we analyzed alignments of 10,238 human genes to their orthologues in chimpanzee and macaque. Using a likelihood ratio test, we identified protein-coding sequences with an accelerated rate of base substitutions along the human lineage. Exons evolving at a fast rate in humans have a significant tendency to contain clusters of AT-to-GC (weak-to-strong) biased substitutions. This pattern is also observed in noncoding sequence flanking rapidly evolving exons. Accelerated exons occur in regions with elevated male recombination rates and exhibit an excess of nonsynonymous substitutions relative to the genomic average. We next analyzed genes with significantly elevated ratios of nonsynonymous to synonymous rates of base substitution (d_N/d_S) along the human lineage, and those with an excess of amino acid replacement substitutions relative to human polymorphism. These genes also show evidence of clusters of weak-to-strong biased substitutions. These findings indicate that a recombination-associated process, such as biased gene conversion (BGC), is driving fixation of GC alleles in the human genome. This process can lead to accelerated evolution in coding sequences and excess amino acid replacement substitutions, thereby generating significant results for tests of positive selection.     

Molecular Evolution of the Odorant and Gustatory Receptor Genes in Lepidopteran Insects: Implications for Their Adaptation and Speciation

Lepidoptera (comprised of butterflies and moths) is one of the largest groups of insects, including more than 160,000 described species. Chemoreception plays important roles in the adaptation of these species to a wide range of niches, e.g., plant hosts, egg-laying sites, and mates. This study investigated the molecular evolution of the lepidopteran odorant (Or) and gustatory receptor (Gr) genes using recently identified genes from Bombyx mori, Danaus plexippus, Heliconius melpomene, Plutella xylostella, Heliothis virescens, Manduca sexta, Cydia pomonella, and Spodoptera littoralis. A limited number of cases of large lineage-specific gene expansion are observed (except in the P. xylostella lineage), possibly due to selection against tandem gene duplication. There has been strong purifying selection during the evolution of both lepidopteran odorant and gustatory genes, as shown by the low ω values estimated through CodeML analysis, ranging from 0.0093 to 0.3926. However, purifying selection has been relaxed on some amino acid sites in these receptors, leading to sequence divergence, which is a precursor of positive selection on these sequences. Signatures of positive selection were detected only in a few loci from the lineage-specific analysis. Estimation of gene gains and losses suggests that the common ancestor of the Lepidoptera had fewer Or genes compared to extant species and an even more reduced number of Gr genes, particularly within the bitter receptor clade. Multiple gene gains and a few gene losses occurred during the evolution of Lepidoptera. Gene family expansion may be associated with the adaptation of lepidopteran species to plant hosts, especially after angiosperm radiation. Phylogenetic analysis of the moth sex pheromone receptor genes suggested that chromosomal translocations have occurred several times. New sex pheromone receptors have arisen through tandem gene duplication. Positive selection was detected at some amino acid sites predicted to be in the extracellular and transmembrane regions of the newly duplicated genes, which might be associated with the evolution of the new pheromone receptors.     

Testing the neutral theory of molecular evolution using genomic data: a comparison of the human and bovine transcriptome

Despite growing evidence of rapid evolution in protein coding genes, the contribution of positive selection to intra- and interspecific differences in protein coding regions of the genome is unclear. We attempted to see if genes coding for secreted proteins and genes with narrow expression, specifically those preferentially expressed in the mammary gland, have diverged at a faster rate between domestic cattle (Bos taurus) and humans (Homo sapiens) than other genes and whether positive selection is responsible. Using a large data set, we identified groups of genes based on secretion and expression patterns and compared them for the rate of nonsynonymous (dN) and synonymous (dS) substitutions per site and the number of radical (Dr) and conservative (Dc) amino acid substitutions. We found evidence of rapid evolution in genes with narrow expression, especially for those expressed in the liver and mammary gland and for genes coding for secreted proteins. We compared common human polymorphism data with human-cattle divergence and found that genes with high evolutionary rates in human-cattle divergence also had a large number of common human polymorphisms. This argues against positive selection causing rapid divergence in these groups of genes. In most cases dN/dS ratios were lower in human-cattle divergence than in common human polymorphism presumably due to differences in the effectiveness of purifying selection between long-term divergence and short-term polymorphism.     

A phylogenetic approach to test for evidence of parental conflict or gene duplications associated with protein-encoding imprinted orthologous genes in placental mammals

There are multiple theories on the evolution of genomic imprinting. We investigated whether the molecular evolution of true orthologs of known imprinted genes provides support for theories based on gene duplication or parental conflicts (where mediated by amino-acid changes). Our analysis of 34 orthologous genes demonstrates that the vast majority of mammalian imprinted genes have not undergone any subsequent significant gene duplication within placental species, suggesting that selection pressures against gene duplication events could be operating for imprinted loci. As antagonistic co-evolution between imprinted genes can regulate offspring growth, proteins mediating this interaction could be subject to rapid evolution via positive selection. Supporting this, we detect evidence of site specific positive selection for the imprinted genes OSBPL5 (and GNASXL), and detect lineage-specific positive selection for 14 imprinted genes where it is known that the gene is imprinted in a specific lineage, namely for: PLAGL1, IGF2, SLC22A18, OSBPL5, DCN, DLK1, RASGRF1, IGF2R, IMPACT, GRB10, NAPIL4, UBE3A, GATM and GABRG3. However, there is an overall lack of concordance between the known imprinting status of each gene (i.e. whether the gene is imprinted or biallelically expressed in a particular mammalian lineage) and positive selection. While only a small number of orthologs of imprinted loci display evidence of positive selection, we observe that the majority of orthologs of imprinted loci display high levels of micro-synteny conservation and have undergone very few cis- or trans-duplications in placental mammalian lineages.     

Tissue-Specific Evolution of Protein Coding Genes in Human and Mouse

Protein-coding genes evolve at different rates, and the influence of different parameters, from gene size to expression level, has been extensively studied. While in yeast gene expression level is the major causal factor of gene evolutionary rate, the situation is more complex in animals. Here we investigate these relations further, especially taking in account gene expression in different organs as well as indirect correlations between parameters. We used RNA-seq data from two large datasets, covering 22 mouse tissues and 27 human tissues. Over all tissues, evolutionary rate only correlates weakly with levels and breadth of expression. The strongest explanatory factors of purifying selection are GC content, expression in many developmental stages, and expression in brain tissues. While the main component of evolutionary rate is purifying selection, we also find tissue-specific patterns for sites under neutral evolution and for positive selection. We observe fast evolution of genes expressed in testis, but also in other tissues, notably liver, which are explained by weak purifying selection rather than by positive selection.     

The Roles of Gene Duplication,Gene Conversion and Positive Selection in Rodent Esp and Mup Pheromone Gene Families with Comparison to the Abp Family

Three proteinaceous pheromone families, the androgen-binding proteins (ABPs), the exocrine-gland secreting peptides (ESPs) and the major urinary proteins (MUPs) are encoded by large gene families in the genomes of Mus musculus and Rattus norvegicus. We studied the evolutionary histories of the Mup and Esp genes and compared them with what is known about the Abp genes. Apparently gene conversion has played little if any role in the expansion of the mouse Class A and Class B Mup genes and pseudogenes, and the rat Mups. By contrast, we found evidence of extensive gene conversion in many Esp genes although not in all of them. Our studies of selection identified at least two amino acid sites in β-sheets as having evolved under positive selection in the mouse Class A and Class B MUPs and in rat MUPs. We show that selection may have acted on the ESPs by determining K_a/K_s for Exon 3 sequences with and without the converted sequence segment. While it appears that purifying selection acted on the ESP signal peptides, the secreted portions of the ESPs probably have undergone much more rapid evolution. When the inner gene converted fragment sequences were removed, eleven Esp paralogs were present in two or more pairs with K_a/K_s >1.0 and thus we propose that positive selection is detectable by this means in at least some mouse Esp paralogs. We compare and contrast the evolutionary histories of all three mouse pheromone gene families in light of their proposed functions in mouse communication.     

Molecular evolution of <Emphasis Type="Italic">PKD2</Emphasis> gene family in mammals

PKD2 gene encodes a critical cation channel protein that plays important roles in various developmental processes and is usually evolutionarily conserved. In the present study, we analyzed the evolutionary patterns of PKD2 and its homologous genes (PKD2L1, PKD2L2) from nine mammalian species. In this study, we demonstrated the orthologs of PKD2 gene family evolved under a dominant purifying selection force. Our results in combination with the reported evidences from functional researches suggested the entire PKD2 gene family are conserved and perform essential biological roles during mammalian evolution. In rodents, PKD2 gene family members appeared to have evolved more rapidly than other mammalian lineages, probably resulting from relaxation of purifying selection. However, positive selection imposed on synonymous sites also potentially contributed to this case. For the paralogs, our results implied that PKD2L2 genes evolved under a weaker purifying selection constraint than PKD2 and PKD2L1 genes. Interestingly, some loop regions of transmembrane domain of PKD2L2 exhibited higher P _N/P _S ratios than expected, suggesting these regions are more functional divergent in organisms and worthy of special attention. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. Chun Ye, Huan Sun have contributed equally to this work.     

Characterization of the definitive classical calpain family of vertebrates using phylogenetic,evolutionary and expression analyses

The calpains are a superfamily of proteases with extensive relevance to human health and welfare. Vast research attention is given to the vertebrate ‘classical’ subfamily, making it surprising that the evolutionary origins, distribution and relationships of these genes is poorly characterized. Consequently, there exists uncertainty about the conservation of gene family structure, function and expression that has been principally defined from work with mammals. Here, more than 200 vertebrate classical calpains were incorporated in phylogenetic analyses spanning an unprecedented range of taxa, including jawless and cartilaginous fish. We demonstrate that the common vertebrate ancestor had at least six classical calpains, including a single gene that gave rise to CAPN11, 1, 2 and 8 in the early jawed fish lineage, plus CAPN3, 9, 12, 13 and a novel calpain gene, hereafter named CAPN17. We reveal that while all vertebrate classical calpains have been subject to persistent purifying selection during evolution, the degree and nature of selective pressure has often been lineage-dependent. The tissue expression of the complete classic calpain family was assessed in representative teleost fish, amphibians, reptiles and mammals. This highlighted systematic divergence in expression across vertebrate taxa, with most classic calpain genes from fish and amphibians having more extensive tissue distribution than in amniotes. Our data suggest that classical calpain functions have frequently diverged during vertebrate evolution and challenge the ongoing value of the established system of classifying calpains by expression.     

Whole-gene positive selection, elevated synonymous substitution rates, duplication, and indel evolution of the chloroplast clpP1 gene

Background

Synonymous DNA substitution rates in the plant chloroplast genome are generally relatively slow and lineage dependent. Non-synonymous rates are usually even slower due to purifying selection acting on the genes. Positive selection is expected to speed up non-synonymous substitution rates, whereas synonymous rates are expected to be unaffected. Until recently, positive selection has seldom been observed in chloroplast genes, and large-scale structural rearrangements leading to gene duplications are hitherto supposed to be rare.

Methodology/Principle Findings

We found high substitution rates in the exons of the plastid clpP1 gene in Oenothera (the Evening Primrose family) and three separate lineages in the tribe Sileneae (Caryophyllaceae, the Carnation family). Introns have been lost in some of the lineages, but where present, the intron sequences have substitution rates similar to those found in other introns of their genomes. The elevated substitution rates of clpP1 are associated with statistically significant whole-gene positive selection in three branches of the phylogeny. In two of the lineages we found multiple copies of the gene. Neighboring genes present in the duplicated fragments do not show signs of elevated substitution rates or positive selection. Although non-synonymous substitutions account for most of the increase in substitution rates, synonymous rates are also markedly elevated in some lineages. Whereas plant clpP1 genes experiencing negative (purifying) selection are characterized by having very conserved lengths, genes under positive selection often have large insertions of more or less repetitive amino acid sequence motifs.

Conclusions/Significance

We found positive selection of the clpP1 gene in various plant lineages to correlated with repeated duplication of the clpP1 gene and surrounding regions, repetitive amino acid sequences, and increase in synonymous substitution rates. The present study sheds light on the controversial issue of whether negative or positive selection is to be expected after gene duplications by providing evidence for the latter alternative. The observed increase in synonymous substitution rates in some of the lineages indicates that the detection of positive selection may be obscured under such circumstances. Future studies are required to explore the functional significance of the large inserted repeated amino acid motifs, as well as the possibility that synonymous substitution rates may be affected by positive selection.     

Tightly Constrained Genome Reduction and Relaxation of Purifying Selection during Secondary Plastid Endosymbiosis

Endosymbiosis, the establishment of a former free-living prokaryotic or eukaryotic cell as an organelle inside a host cell, can dramatically alter the genomic architecture of the endosymbiont. Plastids or chloroplasts, the light-harvesting organelle of photosynthetic eukaryotes, are excellent models to study this phenomenon because plastid origin has occurred multiple times in evolution. Here, we investigate the genomic signature of molecular processes acting through secondary plastid endosymbiosis—the origination of a new plastid from a free-living eukaryotic alga. We used phylogenetic comparative methods to study gene loss and changes in selective regimes on plastid genomes, focusing on green algae that have given rise to three independent lineages with secondary plastids (euglenophytes, chlorarachniophytes, and Lepidodinium). Our results show an overall increase in gene loss associated with secondary endosymbiosis, but this loss is tightly constrained by the retention of genes essential for plastid function. The data show that secondary plastids have experienced temporary relaxation of purifying selection during secondary endosymbiosis. However, this process is tightly constrained, with selection relaxed only relative to the background in primary plastids. Purifying selection remains strong in absolute terms even during the endosymbiosis events. Selection intensity rebounds to pre-endosymbiosis levels following endosymbiosis events, demonstrating the changes in selection efficiency during different origin phases of secondary plastids. Independent endosymbiosis events in the euglenophytes, chlorarachniophytes, and Lepidodinium differ in their degree of relaxation of selection, highlighting the different evolutionary contexts of these events. This study reveals the selection–drift interplay during secondary endosymbiosis and evolutionary parallels during organellogenesis.     

FGF gene family characterization provides insights into its adaptive evolution in Carnivora

Fibroblast growth factors (FGFs) encoded by the FGF gene family can regulate development and physiology in animals. However, their evolutionary characteristics in Carnivora are largely unknown. In this study, we identified 660 sequences of three types of FGF genes from 30 unannotated genomes of Carnivora animals (before 7th May 2020), and the FGF genes from 52 Carnivora species were analyzed through the method of comparative genomics. Phylogenetic and selective pressure analyses were carried out based on the FGF genes of these 52 Carnivora species. The phylogenetic analysis results demonstrated that the FGF gene family was divided into 10 subfamilies and that FGF5 formed one clade rather than belonging to the subfamilies of FGF4 and FGF6. The evolutionary analysis results showed that the FGF genes were prominently subjected to purifying selection and were highly conserved in the process of Carnivora evolution. We also carried out phylogenetic comparative analyses, which indicated that the habitat was one of the factors that shaped the evolution of Carnivora FGF genes. The FGF1 and FGF6 genes were positively selected in the Carnivora animals, and positive selection signals were detected for the FGF19 gene in semiaquatic Carnivora animals. In summary, we clarified the phylogenetic and evolutionary characteristics of Carnivora FGF genes and provided valuable data for future studies on evolutionary characterization of Carnivora animals.