Evidence of adaptive evolution of accessory gland proteins in closely related species of the Drosophila repleta group

Accessory gland proteins (Acps) are part of the seminal fluid of Drosophila species. These proteins have important reproductive functions, being responsible for the proper functioning of several steps of the fertilization process. Acps also contribute indirectly for the reproductive success of males by modulating female behavior. Evidence that Acps participate in sperm competition and sexual conflict includes findings that, on average, Acps have fast evolutionary rates, suggestive of adaptive evolution. This is especially true in species of the Drosophila repleta group. Nevertheless, only in a few occasions have robust statistical tests been used to determine whether observed evolutionary rates are in fact due to positive selection on amino acid substitutions between related species. Here we apply maximum likelihood tests for positive selection on 14 Acps of the D. repleta group. To increase statistical robustness, we use at least 8 sequences, all belonging to species of the Drosophila mulleri complex, for each gene analyzed. We found significant evidence of adaptive evolution for 10 of the tested genes. Among these, the ones with a conserved protein domain had positively selected sites within the functional region of the sequence. We also detected one instance of lineage-specific adaptive evolution in a clade formed by 2 sister species.     

Plastome organization and evolution of chloroplast genes in Cardamine species adapted to contrasting habitats

Background

Plastid genomes, also known as plastomes, are shaped by the selective forces acting on the fundamental cellular functions they code for and thus they are expected to preserve signatures of the adaptive path undertaken by different plant species during evolution. To identify molecular signatures of positive selection associated to adaptation to contrasting ecological niches, we sequenced with Solexa technology the plastomes of two congeneric Brassicaceae species with different habitat preference, Cardamine resedifolia and Cardamine impatiens.

Results

Following in-depth characterization of plastome organization, repeat patterns and gene space, the comparison of the newly sequenced plastomes between each other and with 15 fully sequenced Brassicaceae plastomes publically available in GenBank uncovered dynamic variation of the IR boundaries in the Cardamine lineage. We further detected signatures of positive selection in ten of the 75 protein-coding genes of the examined plastomes, identifying a range of chloroplast functions putatively involved in adaptive processes within the family. For instance, the three residues found to be under positive selection in RUBISCO could possibly be involved in the modulation of RUBISCO aggregation/activation and enzymatic specificty in Brassicaceae. In addition, our results points to differential evolutionary rates in Cardamine plastomes.

Conclusions

Overall our results support the existence of wider signatures of positive selection in the plastome of C. resedifolia, possibly as a consequence of adaptation to high altitude environments. We further provide a first characterization of the selective patterns shaping the Brassicaceae plastomes, which could help elucidate the driving forces underlying adaptation and evolution in this important plant family.

Electronic supplementary material

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

Comparing patterns of natural selection across species using selective signatures

Comparing gene expression profiles over many different conditions has led to insights that were not obvious from single experiments. In the same way, comparing patterns of natural selection across a set of ecologically distinct species may extend what can be learned from individual genome-wide surveys. Toward this end, we show how variation in protein evolutionary rates, after correcting for genome-wide effects such as mutation rate and demographic factors, can be used to estimate the level and types of natural selection acting on genes across different species. We identify unusually rapidly and slowly evolving genes, relative to empirically derived genome-wide and gene family-specific background rates for 744 core protein families in 30 γ-proteobacterial species. We describe the pattern of fast or slow evolution across species as the “selective signature” of a gene. Selective signatures represent a profile of selection across species that is predictive of gene function: pairs of genes with correlated selective signatures are more likely to share the same cellular function, and genes in the same pathway can evolve in concert. For example, glycolysis and phenylalanine metabolism genes evolve rapidly in Idiomarina loihiensis, mirroring an ecological shift in carbon source from sugars to amino acids. In a broader context, our results suggest that the genomic landscape is organized into functional modules even at the level of natural selection, and thus it may be easier than expected to understand the complex evolutionary pressures on a cell.     

Rapid evolution of coral proteins responsible for interaction with the environment

BACKGROUND: Corals worldwide are in decline due to climate change effects (e.g., rising seawater temperatures), pollution, and exploitation. The ability of corals to cope with these stressors in the long run depends on the evolvability of the underlying genetic networks and proteins, which remain largely unknown. A genome-wide scan for positively selected genes between related coral species can help to narrow down the search space considerably. METHODOLOGY/PRINCIPAL FINDINGS: We screened a set of 2,604 putative orthologs from EST-based sequence datasets of the coral species Acropora millepora and Acropora palmata to determine the fraction and identity of proteins that may experience adaptive evolution. 7% of the orthologs show elevated rates of evolution. Taxonomically-restricted (i.e. lineage-specific) genes show a positive selection signature more frequently than genes that are found across many animal phyla. The class of proteins that displayed elevated evolutionary rates was significantly enriched for proteins involved in immunity and defense, reproduction, and sensory perception. We also found elevated rates of evolution in several other functional groups such as management of membrane vesicles, transmembrane transport of ions and organic molecules, cell adhesion, and oxidative stress response. Proteins in these processes might be related to the endosymbiotic relationship corals maintain with dinoflagellates in the genus Symbiodinium. CONCLUSION/RELEVANCE: This study provides a birds-eye view of the processes potentially underlying coral adaptation, which will serve as a foundation for future work to elucidate the rates, patterns, and mechanisms of corals' evolutionary response to global climate change.     

Divergent adaptation to Qinghai-Tibetan Plateau implicated from transciptome study of Gymnocypris dobula and Schizothorax nukiangensis

The Schizothoracine fishes are widely distributed in the Qinghai-Tibetan Plateau (QTP) area and its peripheral regions, which provide a prime example of adaptation in highland aquatic environments. Recent progresses have revealed various genetic adaptations of these fishes by comparing to distantly related lowerland species, however, comparative studies on closely-related species of different altitudes are still lacking. In this study, we sequenced and annotated a primitive Schizothoracine fish Schizothorax nukiangensis Tsao and a highly specialized one Gymnocypris dobula. We performed evolutionary analyses to investigate the candidate genes and signaling pathways involved QTP highland adaptation in both Schizothoracine fishes. Analysis of the 11,007 one-copy orthologs to the primitive cyprinid species, Danio rerio, revealed that both G. dobula and S. nukiangensis showed elevated evolutionary rates. A large number of genes related to hypoxia, including genes involved metabolic processes and cardiovascular system development, exhibited signatures of positive selection in both Schizothoracine fishes, but very few positively selected genes were found overlapping among these Schizothoracines. Our results indicated divergent genetic adaptation to highland environment for aquatic species living in QTP.     

Hearing Aid for Vertebrates via Multiple Episodic Adaptive Events on Prestin Genes

Auditory detection is essential for survival and reproduction of vertebrates, yet the genetic changes underlying the evolution and diversity of hearing are poorly documented. Recent discoveries concerning prestin, which is responsible for cochlear amplification by electromotility, provide an opportunity to redress this situation. We identify prestin genes from the genomes of 14 vertebrates, including three fishes, one amphibian, one lizard, one bird, and eight mammals. An evolutionary analysis of these sequences and 34 previously known prestin genes reveals for the first time that this hearing gene was under positive selection in the most recent common ancestor (MRCA) of tetrapods. This discovery might document the genetic basis of enhanced high sound sensibility in tetrapods. An investigation of the adaptive gain and evolution of electromotility, an important evolutionary innovation for the highest hearing ability of mammals, detects evidence for positive selections on the MRCA of mammals, therians, and placentals, respectively. It is suggested that electromotility determined by prestin might initially appear in the MRCA of mammals, and its functional improvements might occur in the MRCA of therian and placental mammals. Our patch clamp experiments further support this hypothesis, revealing the functional divergence of voltage-dependent nonlinear capacitance of prestin from platypus, opossum, and gerbil. Moreover, structure-based cdocking analyses detect positively selected amino acids in the MRCA of placental mammals that are key residues in sulfate anion transport. This study provides new insights into the adaptation and functional diversity of hearing sensitivity in vertebrates by evolutionary and functional analysis of the hearing gene prestin.     

核糖核酸酶基因超家族分子进化

核糖核酸酶基因(Ribonuclease A, RNASE A)超家族是进化生物学中研究新基因起源及新功能演变的重要模式系统之一。RNASE A超家族中的很多成员表现出基因复制的进化模式, 而且在适应性(正)选择的驱动下, 发生了功能分化。文章综述了RNASE A超家族成员在不同动物类群中进化模式的研究进展, 包括近年来越来越多在基因组水平上开展的相关研究, 显示该基因超家族可能具有比人们以往认识的更为复杂的基因进化模式。随着越来越多动物基因组数据的产生, 对更多动物代表类群进行RNASE A超家族研究, 将有望揭示新的进化机制和功能分化, 为系统认识动物适应进化的遗传机制奠定基础。     

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.     

Effective population size is positively correlated with levels of adaptive divergence among annual sunflowers

The role of adaptation in the divergence of lineages has long been a central question in evolutionary biology, and as multilocus sequence data sets have become available for a wide range of taxa, empirical estimates of levels of adaptive molecular evolution are increasingly common. Estimates vary widely among taxa, with high levels of adaptive evolution in Drosophila, bacteria, and viruses but very little evidence of widespread adaptive evolution in hominids. Although estimates in plants are more limited, some recent work has suggested that rates of adaptive evolution in a range of plant taxa are surprisingly low and that there is little association between adaptive evolution and effective population size in contrast to patterns seen in other taxa. Here, we analyze data from 35 loci for six sunflower species that vary dramatically in effective population size. We find that rates of adaptive evolution are positively correlated with effective population size in these species, with a significant fraction of amino acid substitutions driven by positive selection in the species with the largest effective population sizes but little or no evidence of adaptive evolution in species with smaller effective population sizes. Although other factors likely contribute as well, in sunflowers effective population size appears to be an important determinant of rates of adaptive evolution.     

Adaptive evolution of fertilization proteins within a genus: variation in ZP2 and ZP3 in deer mice (Peromyscus)

Rapid evolution of reproductive proteins has been documented in a wide variety of taxa. In internally fertilized species, knowledge about the evolutionary dynamics of these proteins between closely related taxa is primarily limited to accessory gland proteins in the semen of Drosophila. Investigation of additional taxa and functional classes of proteins is necessary in order to determine if there is a general pattern of adaptive evolution of reproductive proteins between recently diverged species. We performed an evolutionary analysis of 2 egg coat proteins, ZP2 and ZP3, in 15 species of deer mice (genus Peromyscus). Both of these proteins are involved in egg-sperm binding, a critical step in maintaining species-specific fertilization. Here, we show that Zp2 and Zp3 gene trees are not consistent with trees based on nonreproductive genes, Mc1r and Lcat, where species formed monophyletic clades. In fact, for both of the reproductive genes, intraspecific amino acid variation was extensive and alleles were sometimes shared across species. We document positive selection acting on ZP2 and ZP3 and identify specific amino acid sites that are likely targets of selection using both maximum likelihood approaches and patterns of parallel amino acid change. In ZP3, positively selected sites are clustered in and around the region implicated in sperm binding in Mus, suggesting changes may impact egg-sperm binding and fertilization potential. Finally, we identify lineages with significantly elevated rates of amino acid substitution using a Bayesian mapping approach. These findings demonstrate that the pattern of adaptive reproductive protein evolution found at higher taxonomic levels can be documented between closely related mammalian species, where reproductive isolation has evolved recently.     

A genome‐wide scan for selection signatures in Yorkshire and Landrace pigs based on sequencing data

Pigs have experienced dramatic selection due to domestication, which has led to many different phenotypes when compared to their wild counterparts, especially in the last several decades. Currently, genome‐wide scans in both cattle and humans showing positive selection footprints have been investigated. However, few studies have focused on porcine selection footprints, particularly on a genome‐wide scale. Surveying for selection footprints across porcine genomes can be quite valuable for revealing the genetic mechanisms of phenotypic diversity. Here, we employed a medium sequencing depth (5–20x/site per individual, on average) approach called genotyping by genome reducing and sequencing (GGRS) to detect genome‐wide selection signatures of two domestic pig breeds (Yorkshire and Landrace) that have been under intensive selection for traits of muscle development, growth and behavior. The relative extended haplotype homozygosity test, which identifies selection signatures by measuring the characteristics of haplotypes’ frequency distribution within a single population, was also applied to identify potential positively selected regions. As a result, signatures of positive selection were found in each breed. However, most selection signatures were population specific and related to genomic regions containing genes for biological categories including brain development, metabolism, growth and olfaction. Furthermore, the result of the gene set enrichment analysis indicated that selected regions of the two breeds presented a different over‐representation of genes in the Gene Ontology annotations and Kyoto Encyclopedia of Genes and Genomes pathways. Our results revealed a genome‐wide map of selection footprints in pigs and may help us better understand the mechanisms of selection in pig breeding.     

Genome‐wide adaptive evolution to underground stresses in subterranean mammals: Hypoxia adaption,immunity promotion,and sensory specialization

Life underground has provided remarkable examples of adaptive evolution in subterranean mammals; however, genome‐wide adaptive evolution to underground stresses still needs further research. There are approximately 250 species of subterranean mammals across three suborders and six families. These species not only inhabit hypoxic and dark burrows but also exhibit evolved adaptation to hypoxia, cancer resistance, and specialized sensory systems, making them an excellent model of evolution. The adaptive evolution of subterranean mammals has attracted great attention and needs further study. In the present study, phylogenetic analysis of 5,853 single‐copy orthologous gene families of five subterranean mammals (Nannospalax galili, Heterocephalus glaber, Fukomys damarensis, Condylura cristata, and Chrysochloris asiatica) showed that they formed fou distinct clusters. This result is consistent with the traditional systematics of these species. Furthermore, comparison of the high‐quality genomes of these five subterranean mammalian species led to the identification of the genomic signatures of adaptive evolution. Our results show that the five subterranean mammalian did not share positively selected genes but had similar functional enrichment categories, including hypoxia tolerance, immunity promotion, and sensory specialization, which adapted to the environment of underground stresses. Moreover, variations in soil hardness, climate, and lifestyles have resulted in different molecular mechanisms of adaptation to the hypoxic environment and different degrees of visual degradation. These results provide insights into the genome‐wide adaptive evolution to underground stresses in subterranean mammals, with special focus on the characteristics of hypoxia adaption, immunity promotion, and sensory specialization response to the life underground.     

Genomic convergence underlying high-altitude adaptation in alpine plants

Evolutionary convergence is one of the most striking examples of adaptation driven by natural selection.However, genomic evidence for convergent adaptation to extreme environments remains scarce.Here, we assembled reference genomes of two alpine plants, Saussurea obvallata(Asteraceae)and Rheum alexandrae(Polygonaceae), with 37,938 and 61,463 annotated protein-coding genes. By integrating an additional five alpine genomes,we elucidated genomic convergence underlying high-altitude adaptation in al...     

The genomic footprint of climate adaptation in Chironomus riparius

The gradual heterogeneity of climatic factors poses varying selection pressures across geographic distances that leave signatures of clinal variation in the genome. Separating signatures of clinal adaptation from signatures of other evolutionary forces, such as demographic processes, genetic drift and adaptation, to nonclinal conditions of the immediate local environment is a major challenge. Here, we examine climate adaptation in five natural populations of the harlequin fly Chironomus riparius sampled along a climatic gradient across Europe. Our study integrates experimental data, individual genome resequencing, Pool‐Seq data and population genetic modelling. Common‐garden experiments revealed significantly different population growth rates at test temperatures corresponding to the population origin along the climate gradient, suggesting thermal adaptation on the phenotypic level. Based on a population genomic analysis, we derived empirical estimates of historical demography and migration. We used an F_ST outlier approach to infer positive selection across the climate gradient, in combination with an environmental association analysis. In total, we identified 162 candidate genes as genomic basis of climate adaptation. Enriched functions among these candidate genes involved the apoptotic process and molecular response to heat, as well as functions identified in studies of climate adaptation in other insects. Our results show that local climate conditions impose strong selection pressures and lead to genomic adaptation despite strong gene flow. Moreover, these results imply that selection to different climatic conditions seems to converge on a functional level, at least between different insect species.     

Unusual signatures of highly adaptable R-loci in closely-related Arabidopsis species

Plant resistance genes (R-genes) evolve rapidly in response to changing environments. What are the most remarkable signatures of fast adaptive genes, besides the commonly revealed rapid divergence and high non-synonymous substitution rate? Here we investigated these changes in five R-loci following recent differentiation between Arabidopsis thaliana and Arabidopsis lyrata. Extreme differences in evolutionary rates were observed: e.g., an overall 5.46-9.83-fold different nucleotide diversity at two R-loci between species, ten-fold higher non-synonymous substitution rates within one species versus the other, significantly different Ka/Ks ratios between species for the same R-gene, and high interspecific divergence at one R-locus. Particularly, we observed an elevated level of trans-specific polymorphism at one R-locus and a differentially maintained presence/absence polymorphism at another. The high frequency of ancestral polymorphisms amongst R-genes suggests that the persistence of some functional variation is an important evolutionary mechanism shaping genetic variation in R-genes, while the variation of presence/absence polymorphisms provides a potential mechanism for malleable activation of adaptive resistance response pathways. The distinct patterns among R-genes suggest that the same R-gene ortholog can be quickly shaped by different evolutionary processes, e.g., purifying selection in one species but positive selection in a closely-related species.     

Duplication and diversification of trehalase confers evolutionary advantages on lepidopteran insects

Gene duplication provides a major source of new genes for evolutionary novelty and ecological adaptation. However, the maintenance of duplicated genes and their relevance to adaptive evolution has long been debated. Insect trehalase (Treh) plays key roles in energy metabolism, growth, and stress recovery. Here, we show that the duplication of Treh in Lepidoptera (butterflies and moths) is linked with their adaptation to various environmental stresses. Generally, two Treh genes are present in insects: Treh1 and Treh2. We report three distinct forms of Treh in lepidopteran insects, where Treh1 was duplicated into two gene clusters (Treh1a and Treh1b). These gene clusters differ in gene expression patterns, enzymatic properties, and subcellular localizations, suggesting that the enzymes probably underwent sub‐ and/or neofunctionalization in the lepidopteran insects. Interestingly, selective pressure analysis provided significant evidence of positive selection on duplicate Treh1b gene in lepidopteran insect lineages. Most positively selected sites were located in the alpha‐helical region, and several sites were close to the trehalose binding and catalytic sites. Subcellular adaptation of duplicate Treh1b driven by positive selection appears to have occurred as a result of selected changes in specific sequences, allowing for rapid reprogramming of duplicated Treh during evolution. Our results suggest that gene duplication of Treh and subsequent functional diversification could increase the survival rate of lepidopteran insects through various regulations of intracellular trehalose levels, facilitating their adaptation to diverse habitats. This study provides evidence regarding the mechanism by which gene family expansion can contribute to species adaptation through gene duplication and subsequent functional diversification.     

The Visayan Warty Pig (Sus cebifrons) Genome Provides Insight Into Chromosome Evolution and Sensory Adaptation in Pigs

It is largely unknown how mammalian genomes evolve under rapid speciation and environmental adaptation. An excellent model for understanding fast evolution is provided by the genus Sus, which diverged relatively recently and lacks postzygotic isolation. Here, we present a high-quality reference genome of the Visayan warty pig, which is specialized to a tropical island environment. Comparing the genome sequences and chromatin contact maps of the Visayan warty pig (Sus cebifrons) and domestic pig (Sus scrofa), we characterized the dynamics of chromosomal structure evolution during Sus speciation, revealing the similar chromosome conformation as the potential biological mechanism of frequent postdivergence hybridization among Suidae. We further investigated the different signatures of adaptive selection and domestication in Visayan warty pig and domestic pig with specific emphasize on the evolution of olfactory and gustatory genes, elucidating higher olfactory diversity in Visayan warty pig and positive and relaxed evolution of bitter and fat taste receptors, respectively, in domestic pig. Our comprehensive evolutionary and comparative genome analyses provide insight into the dynamics of genomes and how these change over relative short evolutionary times, as well as how these genomic differences encode for differences in the phenotypes.     

Allelic imbalance in Drosophila hybrid heads: exons, isoforms, and evolution

Unraveling how regulatory divergence contributes to species differences and adaptation requires identifying functional variants from among millions of genetic differences. Analysis of allelic imbalance (AI) reveals functional genetic differences in cis regulation and has demonstrated differences in cis regulation within and between species. Regulatory mechanisms are often highly conserved, yet differences between species in gene expression are extensive. What evolutionary forces explain widespread divergence in cis regulation? AI was assessed in Drosophila melanogaster-Drosophila simulans hybrid female heads using RNA-seq technology. Mapping bias was virtually eliminated by using genotype-specific references. Allele representation in DNA sequencing was used as a prior in a novel Bayesian model for the estimation of AI in RNA. Cis regulatory divergence was common in the organs and tissues of the head with 41% of genes analyzed showing significant AI. Using existing population genomic data, the relationship between AI and patterns of sequence evolution was examined. Evidence of positive selection was found in 30% of cis regulatory divergent genes. Genes involved in defense, RNAi/RISC complex genes, and those that are sex regulated are enriched among adaptively evolving cis regulatory divergent genes. For genes in these groups, adaptive evolution may play a role in regulatory divergence between species. However, there is no evidence that adaptive evolution drives most of the cis regulatory divergence that is observed. The majority of genes showed patterns consistent with stabilizing selection and neutral evolutionary processes.