Adaptive molecular evolution in the opsin genes of rapidly speciating cichlid species

Cichlid fish inhabit a diverse range of environments that vary in the spectral content of light available for vision. These differences should result in adaptive selective pressure on the genes involved in visual sensitivity, the opsin genes. This study examines the evidence for differential adaptive molecular evolution in East African cichlid opsin genes due to gross differences in environmental light conditions. First, we characterize the selective regime experienced by cichlid opsin genes using a likelihood ratio test format, comparing likelihood models with different constraints on the relative rates of amino acid substitution, across sites. Second, we compare turbid and clear lineages to determine if there is evidence of differences in relative rates of substitution. Third, we present evidence of functional diversification and its relationship to the photic environment among cichlid opsin genes. We report statistical evidence of positive selection in all cichlid opsin genes, except short wavelength-sensitive 1 and short wavelength-sensitive 2b. In all genes predicted to be under positive selection, except short wavelength-sensitive 2a, we find differences in selective pressure between turbid and clear lineages. Potential spectral tuning sites are variable among all cichlid opsin genes; however, patterns of substitution consistent with photic environment-driven evolution of opsin genes are observed only for short wavelength-sensitive 1 opsin genes. This study identifies a number of promising candidate-tuning sites for future study by site-directed mutagenesis. This work also begins to demonstrate the molecular evolutionary dynamics of cichlid visual sensitivity and its relationship to the photic environment.     

Adaptive evolution in the photosensory domain of phytochrome A in early angiosperms

Flowering plant diversity now far exceeds the combined diversity of all other plant groups. Recently identified extant remnants of the earliest-diverging lines suggest that the first angiosperms may have lived in shady, disturbed, and moist understory habitats, and that the aquatic habit also arose early. This would have required the capacity to begin life in dimly lit environments. If so, evolution in light-sensing mechanisms may have been crucial to their success. The photoreceptor phytochrome A is unique among angiosperm phytochromes in its capacity to serve a transient role under conditions where an extremely high sensitivity is required. We present evidence of altered functional constraints between phytochrome A (PHYA) and its paralog, PHYC. Tests for selection suggest that an elevation in nonsynonymous rates resulted from an episode of selection along the branch leading to all angiosperm PHYA sequences. Most nucleotide sites (95%) are selectively constrained, and the ratio of nonsynonymous to synonymous substitutions on branches within the PHYA clade does not differ from the ratio on the branches in the PHYC clade. Thus, positive selection at a handful of sites, rather than relaxation of selective constraints, apparently has played a major role in the evolution of the photosensory domain of phytochrome A. The episode of selection occurred very early in the history of flowering plants, suggesting that innovation in phyA may have given the first angiosperms some adaptive advantage.     

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.     

趋化因子及其受体家族分子适应性进化分析

Branch-Site模型是检测基因序列中单个密码子位点是否具有选择作用的统计学方法。该模型能有效地检测基因在进化历程中是否受到选择作用, 并预测出那些在进化过程中对功能分化有重要贡献的、受正选择作用的密码子位点。趋化因子是一类控制免疫细胞定向迁移的细胞因子, 其功能行使由趋化因子受体介导。该文用Branch-Site模型分析趋化因子及其受体基因家族的分子适应性, 发现只有少数种类基因受到正选择作用, 如RANTES、CCR5等。并预测出一些可能受到正选择作用的位点, 蛋白3D分析显示, 它们均位于趋化因子和相应受体相互作用的结构区域。     

Adaptive protein evolution and regulatory divergence in Drosophila

Two recent studies demonstrated a positive correlation between divergence in gene expression and protein sequence in Drosophila. This correlation could be driven by positive selection or variation in functional constraint. To distinguish between these alternatives, we compared patterns of molecular evolution for 1,862 genes with two previously reported estimates of expression divergence in Drosophila. We found a slight negative trend (nonsignificant) between positive selection on protein sequence and divergence in expression levels between Drosophila melanogaster and Drosophila simulans. Conversely, shifts in expression patterns during Drosophila development showed a positive association with adaptive protein evolution, though as before the relationship was weak and not significant. Overall, we found no strong evidence for an increase in the incidence of positive selection on protein-coding regions in genes with divergent expression in Drosophila, suggesting that the previously reported positive association between protein and regulatory divergence primarily reflects variation in functional constraint.     

Adaptive evolution of newly emerged micro-RNA genes in Drosophila

How often micro-RNA (miRNA) genes emerged and how fast theyevolved soon after their emergence are some of the central questionsin the evolution of miRNAs. Because most known miRNA genes areancient and highly conserved, these questions can be best answeredby identifying newly emerged miRNA genes. Among the 78 miRNAgenes in Drosophila reported before 2007, only 5 are confirmedto be newly emerged in the genus (although many more can befound in the newly reported data set; e.g., Ruby et al. 2007;Stark et al. 2007; Lu et al. 2008). These new miRNA genes haveundergone numerous changes, even in the normally invariant maturesequences. Four of them (the miR-310/311/312/313 cluster, denotedmiR-310s) were duplicated from other conserved miRNA genes.The fifth one (miR-303) appears to be a very young gene, originatingde novo from a non-miRNA sequence recently. We sequenced these5 miRNA genes and their neighboring regions from a worldwidecollection of Drosophila melanogaster lines. The levels of divergenceand polymorphism in these miRNA genes, vis-à-vis thoseof the neighboring DNA sequences, suggest that these 5 genesare evolving adaptively. Furthermore, the polymorphism patternof miR-310s in D. melanogaster is indicative of hitchhikingunder positive selection. Thus, a large number of adaptive changesover a long period of time may be essential for the evolutionof newly emerged miRNA genes.     

POTION: an end-to-end pipeline for positive Darwinian selection detection in genome-scale data through phylogenetic comparison of protein-coding genes

Background

Detection of genes evolving under positive Darwinian evolution in genome-scale data is nowadays a prevailing strategy in comparative genomics studies to identify genes potentially involved in adaptation processes. Despite the large number of studies aiming to detect and contextualize such gene sets, there is virtually no software available to perform this task in a general, automatic, large-scale and reliable manner. This certainly occurs due to the computational challenges involved in this task, such as the appropriate modeling of data under analysis, the computation time to perform several of the required steps when dealing with genome-scale data and the highly error-prone nature of the sequence and alignment data structures needed for genome-wide positive selection detection.

Results

We present POTION, an open source, modular and end-to-end software for genome-scale detection of positive Darwinian selection in groups of homologous coding sequences. Our software represents a key step towards genome-scale, automated detection of positive selection, from predicted coding sequences and their homology relationships to high-quality groups of positively selected genes. POTION reduces false positives through several sophisticated sequence and group filters based on numeric, phylogenetic, quality and conservation criteria to remove spurious data and through multiple hypothesis corrections, and considerably reduces computation time thanks to a parallelized design. Our software achieved a high classification performance when used to evaluate a curated dataset of Trypanosoma brucei paralogs previously surveyed for positive selection. When used to analyze predicted groups of homologous genes of 19 strains of Mycobacterium tuberculosis as a case study we demonstrated the filters implemented in POTION to remove sources of errors that commonly inflate errors in positive selection detection. A thorough literature review found no other software similar to POTION in terms of customization, scale and automation.

Conclusion

To the best of our knowledge, POTION is the first tool to allow users to construct and check hypotheses regarding the occurrence of site-based evidence of positive selection in non-curated, genome-scale data within a feasible time frame and with no human intervention after initial configuration. POTION is available at http://www.lmb.cnptia.embrapa.br/share/POTION/.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1765-0) contains supplementary material, which is available to authorized users.     

Adaptive evolution of metabolic pathways in Drosophila

The adaptive significance of enzyme variation has been of central interest in population genetics. Yet, how natural selection operates on enzymes in the larger context of biochemical pathways has not been broadly explored. A basic expectation is that natural selection on metabolic phenotypes will target enzymes that control metabolic flux, but how adaptive variation is distributed among enzymes in metabolic networks is poorly understood. Here, we use population genetic methods to identify enzymes responding to adaptive selection in the pathways of central metabolism in Drosophila melanogaster and Drosophila simulans. We report polymorphism and divergence data for 17 genes that encode enzymes of 5 metabolic pathways that converge at glucose-6-phosphate (G6P). Deviations from neutral expectations were observed at five loci. Of the 10 genes that encode the enzymes of glycolysis, only aldolase (Ald) deviated from neutrality. The other 4 genes that were inconsistent with neutral evolution (glucose-6-phosphate dehydrogenase [G6pd]), phosphoglucomutase [Pgm], trehalose-6-phosphate synthetase [Tps1], and glucose-6phosphatase [G6pase] encode G6P branch point enzymes that catalyze reactions at the entry point to the pentose-phosphate, glycogenic, trehalose synthesis, and gluconeogenic pathways. We reconcile these results with population genetics theory and existing arguments on metabolic regulation and propose that the incidence of adaptive selection in this system is related to the distribution of flux control. The data suggest that adaptive evolution of G6P branch point enzymes may have special significance in metabolic adaptation.     

Beyond directed evolution: Darwinian selection as a tool for synthetic biology

Synthetic biology is an engineering approach that seeks to design and construct new biological parts, devices and systems, as well as to re-design existing components. However, rationally designed synthetic circuits may not work as expected due to the context-dependence of biological parts. Darwinian selection, the main mechanism through which evolution works, is a major force in creating biodiversity and may be a powerful tool for synthetic biology. This article reviews selection-based techniques and proposes strict Darwinian selection as an alternative approach for the identification and characterization of parts. Additionally, a strategy for fine-tuning of relatively complex circuits by coupling them to a master standard circuit is discussed.     

Codon-based tests of positive selection,branch lengths,and the evolution of mammalian immune system genes

Using basic probability theory, we show that there is a substantial likelihood that even in the presence of strong purifying selection, there will be a number of codons in which the number of synonymous nucleotide substitutions per site (d (S)) exceeds the number of non-synonymous nucleotide substitutions per site (d (N)). In an empirical study, we examined the numbers of synonymous (b (S)) and non-synonymous substitutions (b (N)) along branches of the phylogenies of 69 single-copy orthologous genes from seven species of mammals. A pattern of b (N) > b (S) was most commonly seen in the shortest branches of the tree and was associated with a high coefficient of variation in both b (N) and b (S), suggesting that high stochastic error in b (N) and b (S) on short branches, rather than positive Darwinian selection, is the explanation of most cases where b (N) is greater than b (S) on a given branch. The branch-site method of Zhang et al. (Zhang, Nielsen, Yang, Mol Biol Evol, 22:2472-2479, 2005) identified 117 codons on 35 branches as "positively selected," but a majority of these codons lacked synonymous substitutions, while in the others, synonymous and non-synonymous differences per site occurred in approximately equal frequencies. Thus, it was impossible to rule out the hypothesis that chance variation in the pattern of mutation across sites, rather than positive selection, accounted for the observed pattern. Our results showed that b (N)/b (S) was consistently elevated in immune system genes, but neither the search for branches with b (N) > b (S) nor the branch-site method revealed this trend.     

Sexual selection and the adaptive evolution of mammalian ejaculate proteins

An elevated rate of substitution characterizes the molecular evolution of reproductive proteins from a wide range of taxa. Although the selective pressures explaining this rapid evolution are yet to be resolved, recent evidence implicates sexual selection as a potentially important explanatory factor. To investigate this hypothesis, we sought evidence of a high rate of adaptive gene evolution linked to postcopulatory sexual selection in muroid rodents, a model vertebrate group displaying a broad range of mating systems. Specifically, we sequenced 7 genes from diverse rodents that are expressed in the testes, prostate, or seminal vesicles, products of which have the potential to act in sperm competition. We inferred positive Darwinian selection in these genes by estimation of the ratio of nonsynonymous (d(N), amino acid changing) to synonymous (d(S), amino acid retaining) substitution rates (omega = d(N)/d(S)). Next, we tested whether variation in this ratio among lineages could be attributed to interspecific variation in mating systems, as inferred from the variation in these rodents' relative testis sizes (RTS). Four of the 7 genes examined (Prm1, Sva, Acrv1, and Svs2, but not Svp2, Msmb, or Spink3) exhibit unambiguous evidence of positive selection. One of these, the seminal vesicle-derived protein Svs2, also shows some evidence for a concentration of positive selection in those lineages in which sperm competition is common. However, this was not a general trend among all the rodent genes we examined. Using the same methods, we then reanalyzed previously published data on 2 primate genes, SEMG1 and SEMG2. Although SEMG2 also shows evidence of positive selection concentrated in lineages subject to high levels of sperm competition, no such trend was found for SEMG1. Overall, despite a high rate of positive selection being a feature of many ejaculate proteins, these results indicate that the action of sexual selection potentially responsible for elevated evolutionary rates may be difficult to detect on a gene-by-gene basis. Although the extreme diversity of reproductive phenotypes exhibited in nature attests to the power of sexual selection, the extent to which this force predominates in driving the rapid molecular evolution of reproductive genes therefore remains to be determined.     

Rapid Evolution, Genetic Variations, and Functional Association of the Human Spermatogenesis-Related Gene NYD-SP12

NYD-SP12 is a recently identified spermatogenesis-related gene with a pivotal role in human testis development. In this study, we analyzed between-species divergence and within-species variation of NYD-SP12 in seven representative primate species, four worldwide human populations, and 124 human clinical subjects. Our results indicate that NYD-SP12 evolves rapidly in both the human and the chimpanzee lineages, which is likely caused by Darwinian positive selection and/or sexual selection. We observed significant interpopulation divergence among human populations, which might be due to the varied demographic histories. In the association analysis, we demonstrated significant frequency discrepancy of a synonymous sequence polymorphism among the clinical groups with different sperm traits. [Reviewing Editor: Dr. Manyuan Long]     

Natural selection in avian protein-coding genes expressed in brain

The evolution of birds from theropod dinosaurs took place approximately 150 million years ago, and was associated with a number of specific adaptations that are still evident among extant birds, including feathers, song and extravagant secondary sexual characteristics. Knowledge about the molecular evolutionary background to such adaptations is lacking. Here, we analyse the evolution of > 5000 protein-coding gene sequences expressed in zebra finch brain by comparison to orthologous sequences in chicken. Mean d _N/ d _S is 0.085 and genes with their maximal expression in the eye and central nervous system have the lowest mean d _N/ d _S value, while those expressed in digestive and reproductive tissues exhibit the highest. We find that fast-evolving genes (those which have higher than expected rate of nonsynonymous substitution, indicative of adaptive evolution) are enriched for biological functions such as fertilization, muscle contraction, defence response, response to stress, wounding and endogenous stimulus, and cell death. After alignment to mammalian orthologues, we identify a catalogue of 228 genes that show a significantly higher rate of protein evolution in the two bird lineages than in mammals. These accelerated bird genes, representing candidates for avian-specific adaptations, include genes implicated in vocal learning and other cognitive processes. Moreover, colouration genes evolve faster in birds than in mammals, which may have been driven by sexual selection for extravagant plumage characteristics.     

Adaptive evolution of a key gene affecting queen and worker traits in the honey bee, Apis mellifera

The vitellogenin egg yolk precursor protein represents a well-studied case of social pleiotropy in the model organism Apis mellifera. Vitellogenin is associated with fecundity in queens and plays a major role in controlling division of labour in workers, thereby affecting both individual and colony-level fitness. We studied the molecular evolution of vitellogenin and seven other genes sequenced in a large population panel of Apis mellifera and several closely related species to investigate the role of social pleiotropy on adaptive protein evolution. We found a significant excess of nonsynonymous fixed differences between A. mellifera, A. cerana and A. florea relative to synonymous sites indicating high rates of adaptive evolution at vitellogenin. Indeed, 88% of amino acid changes were fixed by selection in some portions of the gene. Further, vitellogenin exhibited hallmark signatures of selective sweeps in A. mellifera, including a significant skew in the allele frequency spectrum, extreme levels of genetic differentiation and linkage disequilibrium. Finally, replacement polymorphisms in vitellogenin were significantly enriched in parts of the protein involved in binding lipid, establishing a link between the gene's structure, function and effects on fitness. Our case study provides unequivocal evidence of historical and ongoing bouts of adaptive evolution acting on a key socially pleiotropic gene in the honey bee.     

The probability of parallel evolution

Abstract How often will natural selection drive parallel evolution at the DNA sequence level? More precisely, what is the probability that selection will cause two populations that live in identical environments to substitute the same beneficial mutation? Here I show that, under fairly general conditions, the answer is simple: if a wild‐type sequence can mutate to n different beneficial mutations, replicate populations will on average fix the same mutation with probability P= 2/(n + 1). This probability, which is derived using extreme value theory, is independent of most biological details, including the length of the gene in question and the precise distribution of fitness effects among alleles. I conclude that the probability of parallel evolution under natural selection is nearly twice as large as that under neutrality.     

Mosaic evolution of ruminant stomach lysozyme genes

The genomes of ruminant artiodactyls, such as cow and sheep, have approximately 10 lysozyme genes, 4 of which are expressed in the stomach. Most of the duplications of the lysozyme genes occurred 40-50 million years ago, before the divergence of cow and sheep. Despite this, the coding regions of stomach lysozyme genes within a species (e.g., cow, sheep, or deer) are more similar to each other than to lysozyme genes in other ruminants. This observation suggests that the coding regions of the stomach lysozyme genes have evolved in a concerted fashion. Our previous characterization of 3 cow stomach lysozyme genes suggested that it was only the coding exons that had participated in concerted evolution. To determine whether the introns and flanking regions of ruminant stomach lysozyme genes are evolving in a concerted or a divergent fashion, we have isolated and characterized 2 sheep stomach lysozyme genes. Comparison of the sequences of the sheep and cow stomach lysozyme genes clearly shows that the introns and flanking regions have evolved, like the 3' untranslated region of the mRNAs, in a divergent manner. Thus, if the four coding exons are evolving by concerted evolution, then a mosaic pattern of concerted and divergent evolution is occurring in these genes. The independent concerted evolution of coding exons of the ruminant stomach lysozyme gene may have assisted in the accelerated adaptive evolution of the lysozyme to new function in the early ruminant.     

Adaptive Evolution of cry Genes in Bacillus thuringiensis： Implications for Their Specificity Determination

The cry gene family, produced during the late exponential phase of growth in Bacillus thuringiensis, is a large, still-growing family of homologous genes, in which each gene encodes a protein with strong specific activity against only one or a few insect species. Extensive studies are mostly focusing on the structural and functional relationships of Cry proteins, and have revealed several residues or domains that are important for the target recognition and receptor attachment. In this study, we have employed a maximum likelihood method to detect evidence of adaptive evolution in Cry proteins, and have identified 24 positively selected residues, which are all located in Domain Ⅱ or Ⅲ. Combined with known data from mutagenesis studies, the majority of these residues, at the molecular level, contribute much to the insect specificity determination. We postulate that the potential pressures driving the diversification of Cry proteins may be in an attempt to adapt for the ＂arm race＂ between δ-endotoxins and the targeted insects, or to enlarge their target spectra, hence result in the functional divergence. The sites identified to be under positive selection would provide targets for further structural and functional analyses on Cry proteins.     

蚂蚁入侵适应性进化及其机制

入侵种被传至一个新的区域后,在面对当地某些生态因子的胁迫时如何发生适应性变化,是入侵生物学研究的一个重要课题。这方面信息对评价物种的入侵性、制订防控策略具有重要指导意义。入侵性蚂蚁作为无脊椎陆生动物中入侵最成功、危害最严重的一类,其生物学、生态学已得到大量研究,许多报道表明它们在定殖和扩张过程中会发生适应性变化。基于30多年的文献资料,从蚁群社会结构、超级蚁群形成能力、个体大小、觅食行为、种间竞争、天敌防御、适应低温能力、生态位转变等多个角度进行了全面总结,并针对这些方面适应性变化的机制,从表型可塑性、适应性进化、自身遗传学基础、与共生微生物的互利关系等角度进行了剖析。提出了今后有待重点研究的几个问题。将有助于深化对蚂蚁成功入侵机制的认识,为推进我国红火蚁(Solenopsis invicta Buren)、小火蚁(Wasmannia auropunctata Roger)等入侵性蚂蚁的研究与治理提供参考。