Identifying Signatures of Natural Selection in Tibetan and Andean Populations Using Dense Genome Scan Data

High-altitude hypoxia (reduced inspired oxygen tension due to decreased barometric pressure) exerts severe physiological stress on the human body. Two high-altitude regions where humans have lived for millennia are the Andean Altiplano and the Tibetan Plateau. Populations living in these regions exhibit unique circulatory, respiratory, and hematological adaptations to life at high altitude. Although these responses have been well characterized physiologically, their underlying genetic basis remains unknown. We performed a genome scan to identify genes showing evidence of adaptation to hypoxia. We looked across each chromosome to identify genomic regions with previously unknown function with respect to altitude phenotypes. In addition, groups of genes functioning in oxygen metabolism and sensing were examined to test the hypothesis that particular pathways have been involved in genetic adaptation to altitude. Applying four population genetic statistics commonly used for detecting signatures of natural selection, we identified selection-nominated candidate genes and gene regions in these two populations (Andeans and Tibetans) separately. The Tibetan and Andean patterns of genetic adaptation are largely distinct from one another, with both populations showing evidence of positive natural selection in different genes or gene regions. Interestingly, one gene previously known to be important in cellular oxygen sensing, EGLN1 (also known as PHD2), shows evidence of positive selection in both Tibetans and Andeans. However, the pattern of variation for this gene differs between the two populations. Our results indicate that several key HIF-regulatory and targeted genes are responsible for adaptation to high altitude in Andeans and Tibetans, and several different chromosomal regions are implicated in the putative response to selection. These data suggest a genetic role in high-altitude adaption and provide a basis for future genotype/phenotype association studies necessary to confirm the role of selection-nominated candidate genes and gene regions in adaptation to altitude.     

Adaptation genomics: the next generation

Understanding the genetics of how organisms adapt to changing environments is a fundamental topic in modern evolutionary ecology. The field is currently progressing rapidly because of advances in genomics technologies, especially DNA sequencing. The aim of this review is to first briefly summarise how next generation sequencing (NGS) has transformed our ability to identify the genes underpinning adaptation. We then demonstrate how the application of these genomic tools to ecological model species means that we can start addressing some of the questions that have puzzled ecological geneticists for decades such as: How many genes are involved in adaptation? What types of genetic variation are responsible for adaptation? Does adaptation utilise pre-existing genetic variation or does it require new mutations to arise following an environmental change?     

Why the evolution of resistance to anthropogenic toxins normally involves major gene changes: the limits to natural selection

Standard population genetic theory suggests that adaptation should normally be achieved by the spread of many genes each of small effect (polygenes), and that adaptation by major genes should be unusual. Such models depend on consideration of the rates of acquisition of adaptation. In practice, adaptation to pollutants and anthropogenic toxins has most frequently been achieved by the spread of major genes. A simple model is developed to explain this discrepancy, in which the determining factor is not the rate of spread, but the maximum response achievable under the two contrasting models of polygenic or major gene inheritance. In the short term, for a given mean and genetic variance, characters in which the additive genetic variance is produced by the segregation of many genes of small effect at intermediate gene frequencies are unable to produce as large a response to directional selection as characters in which the variance is caused by genes of large effect at low frequency. If the target for selection is a long way from the mean prior to selection (as it may well be for adaptation to novel anthropogenic stresses) then adaptation can only be achieved by species possessing major genes. The model is discussed with reference to the example of heavy metal tolerance in plants.     

Comparative genomics reveals insights into cyanobacterial evolution and habitat adaptation

Cyanobacteria are photosynthetic prokaryotes that inhabit diverse aquatic and terrestrial environments. However, the evolutionary mechanisms involved in the cyanobacterial habitat adaptation remain poorly understood. Here, based on phylogenetic and comparative genomic analyses of 650 cyanobacterial genomes, we investigated the genetic basis of cyanobacterial habitat adaptation (marine, freshwater, and terrestrial). We show: (1) the expansion of gene families is a common strategy whereby terrestrial cyanobacteria cope with fluctuating environments, whereas the genomes of many marine strains have undergone contraction to adapt to nutrient-poor conditions. (2) Hundreds of genes are strongly associated with specific habitats. Genes that are differentially abundant in genomes of marine, freshwater, and terrestrial cyanobacteria were found to be involved in light sensing and absorption, chemotaxis, nutrient transporters, responses to osmotic stress, etc., indicating the importance of these genes in the survival and adaptation of organisms in specific habitats. (3) A substantial fraction of genes that facilitate the adaptation of Cyanobacteria to specific habitats are contributed by horizontal gene transfer, and such genetic exchanges are more frequent in terrestrial cyanobacteria. Collectively, our results further our understandings of the adaptations of Cyanobacteria to different environments, highlighting the importance of ecological constraints imposed by the environment in shaping the evolution of Cyanobacteria.Subject terms: Phylogenetics, Microbial ecology     

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.     

Habitat variability correlates with duplicate content of Drosophila genomes

The factors limiting the habitat range of species are crucial in understanding their biodiversity and response to environmental change. Yet the genetic and genomic architectures that produce genetic variation to enable environmental adaptation have remained poorly understood. Here we show that the proportion of duplicated genes (P(D)) in the whole genomes of fully sequenced Drosophila species is significantly correlated with environmental variability within the habitats measured by the climatic envelope and habitat diversity. Furthermore, species with a low P(D) tend to lose the duplicated genes owing to their faster evolution. These results indicate that the rapid relaxation of functional constraints on duplicated genes resulted in a low P(D) for species with lower habitat diversity, and suggest that the maintenance of duplicated genes gives organisms an ecological advantage during evolution. We therefore propose that the P(D) in a genome is related to adaptation to environmental variation.     

Molecular Parallelism Underlies Convergent Highland Adaptation of Maize Landraces

Convergent phenotypic evolution provides some of the strongest evidence for adaptation. However, the extent to which recurrent phenotypic adaptation has arisen via parallelism at the molecular level remains unresolved, as does the evolutionary origin of alleles underlying such adaptation. Here, we investigate genetic mechanisms of convergent highland adaptation in maize landrace populations and evaluate the genetic sources of recurrently selected alleles. Population branch excess statistics reveal substantial evidence of parallel adaptation at the level of individual single-nucleotide polymorphism (SNPs), genes, and pathways in four independent highland maize populations. The majority of convergently selected SNPs originated via migration from a single population, most likely in the Mesoamerican highlands, while standing variation introduced by ancient gene flow was also a contributor. Polygenic adaptation analyses of quantitative traits reveal that alleles affecting flowering time are significantly associated with elevation, indicating the flowering time pathway was targeted by highland adaptation. In addition, repeatedly selected genes were significantly enriched in the flowering time pathway, indicating their significance in adapting to highland conditions. Overall, our study system represents a promising model to study convergent evolution in plants with potential applications to crop adaptation across environmental gradients.     

Molecular genetic mapping and plant evolutionary biology

With the advent of molecular genetic mapping, it is possible to study the genetic basis of natural heritable variation in new ways. Here, three potential uses of molecular genetic mapping in plant ecology and evolutionary biology are discussed; (1) accurate estimation of genetic parameters, (2) understanding speciation and/or adaptation, and (3) investigating whole genome organization. Basic methods for mapping genes and important mapping strategies are outlined. Recent studies are introduced to illustrate progress so far in applying the new methods in ecological and evolutionary research.     

Identification of candidate genes related to highland adaptation from multiple Chinese local chicken breeds by whole genome sequencing analysis

Understanding the genetic mechanism of highland adaptation is of great importance for breeding improvement of Tibetan chickens (TBC). Some studies of TBC have identified some candidate genes and pathways from multiple subgroups, but the related genetic mechanisms remain largely unknown. Different genetic backgrounds and the independent genetic basis of highland adaptation make it difficult to identity the selective region of highland adaptation with all TBC samples. In this study, we conducted pre-analysis in a large-scale population to select a TBC subgroup with the purest and highest level of highland-specific lineage for the further analysis. Finally, the 37 samples from a TBC subgroup and 19 Lahsa White chickens were used to represent the highland group for further analysis with 80 samples from five Chinese local lowland breeds as controls. Population structure analysis revealed that highland adaptation significantly affected population stratification in Chinese local chicken breeds. Genome-wide selection signal analysis identified 201 candidate genes associated with highland adaptation of TBC, and these genes were significantly enriched in calcium signaling, vascular smooth muscle contraction and the cellular response to oxidative stress pathways. Additionally, we identified a narrow 1.76 kb region containing an overlapping region between HBZ and an active enhancer, and our identified region showed a highly significant signal. The highland group selected the haplotype with high activity to improve the oxygen-carrying capacity, thus being adapted to a hypoxic environment. We also found that STX2 was significantly selected in the highland group, thus potentially reducing the oxidative stress caused by hypoxia, and that STX2 exhibited the opposite effects on highland adaptation and reproductive traits. Our findings advance our understanding of extreme environment adaptation of highland chickens, and provide some variants and genes beneficial to TBC genetic breeding improvement.     

Adaptation to climate through flowering phenology: a case study in Medicago truncatula

Local climatic conditions likely constitute an important selective pressure on genes underlying important fitness‐related traits such as flowering time, and in many species, flowering phenology and climatic gradients strongly covary. To test whether climate shapes the genetic variation on flowering time genes and to identify candidate flowering genes involved in the adaptation to environmental heterogeneity, we used a large Medicago truncatula core collection to examine the association between nucleotide polymorphisms at 224 candidate genes and both climate variables and flowering phenotypes. Unlike genome‐wide studies, candidate gene approaches are expected to enrich for the number of meaningful trait associations because they specifically target genes that are known to affect the trait of interest. We found that flowering time mediates adaptation to climatic conditions mainly by variation at genes located upstream in the flowering pathways, close to the environmental stimuli. Variables related to the annual precipitation regime reflected selective constraints on flowering time genes better than the other variables tested (temperature, altitude, latitude or longitude). By comparing phenotype and climate associations, we identified 12 flowering genes as the most promising candidates responsible for phenological adaptation to climate. Four of these genes were located in the known flowering time QTL region on chromosome 7. However, climate and flowering associations also highlighted largely distinct gene sets, suggesting different genetic architectures for adaptation to climate and flowering onset.     

Adaptive and neutral genetic differentiation among Scottish and endangered Irish red grouse (<Emphasis Type="Italic">Lagopus lagopus scotica</Emphasis>)

Studying patterns of intra-specific genetic variation among populations allows for a better understanding of population structure and local adaptation. However, those patterns may differ according to the genetic markers applied, as neutral genetic markers reflect demographic processes and random genetic drift, whereas adaptive markers also carry the footprint of selection. In combination, neutral and adaptive genetic markers permit to assess the relative roles of drift and selection in shaping population structure. Among the best understood adaptive genetic loci are the genes of the major histocompatibility complex (MHC). We here study variation and differentiation at neutral SNP markers and MHC class II genes in red grouse (Lagopus lagopus scotica) from Ireland and Scotland. Irish red grouse populations are fragmented and drastically declining, but red grouse are abundant in Scotland. We find evidence for positive selection acting on the MHC genes and variation in MHC gene copy numbers among Irish individuals. Furthermore, there was significant population differentiation among red grouse from Ireland and Scotland at the neutral SNP markers (F_ST = 0.084) and the MHC-BLB genes (F_ST: BLB1 = 0.116, BLB2 = 0.090, BLB3 = 0.104). Differentiation at the MHC-BLB1 was significantly higher than at the neutral SNP markers, suggesting that selection plays an important role in shaping MHC variation, in addition to genetic drift. We speculate that the observed differentiation pattern might be due to local adaptation to different parasite regimes. These findings have strong conservation implications and we advise against the introduction of Scottish red grouse to supplement Irish populations.     

Host–Pathogen Coevolution: The Selective Advantage of Bacillus thuringiensis Virulence and Its Cry Toxin Genes

Reciprocal coevolution between host and pathogen is widely seen as a major driver of evolution and biological innovation. Yet, to date, the underlying genetic mechanisms and associated trait functions that are unique to rapid coevolutionary change are generally unknown. We here combined experimental evolution of the bacterial biocontrol agent Bacillus thuringiensis and its nematode host Caenorhabditis elegans with large-scale phenotyping, whole genome analysis, and functional genetics to demonstrate the selective benefit of pathogen virulence and the underlying toxin genes during the adaptation process. We show that: (i) high virulence was specifically favoured during pathogen–host coevolution rather than pathogen one-sided adaptation to a nonchanging host or to an environment without host; (ii) the pathogen genotype BT-679 with known nematocidal toxin genes and high virulence specifically swept to fixation in all of the independent replicate populations under coevolution but only some under one-sided adaptation; (iii) high virulence in the BT-679-dominated populations correlated with elevated copy numbers of the plasmid containing the nematocidal toxin genes; (iv) loss of virulence in a toxin-plasmid lacking BT-679 isolate was reconstituted by genetic reintroduction or external addition of the toxins. We conclude that sustained coevolution is distinct from unidirectional selection in shaping the pathogen''s genome and life history characteristics. To our knowledge, this study is the first to characterize the pathogen genes involved in coevolutionary adaptation in an animal host–pathogen interaction system.     

Symbiont-mediated adaptation by planthoppers and leafhoppers to resistant rice varieties

For over 50 years, host plant resistance has been the principal focus of public research to reduce planthopper and leafhopper damage to rice in Asia. Several resistance genes have been identified from native varieties and wild rice species, and some of these have been incorporated into high-yielding rice varieties through conventional breeding. However, adaptation by hoppers to resistant rice has been phenomenally rapid, and hopper populations with virulence against several resistance genes are now widespread. Directional genetic selection for virulent hoppers seems unlikely given the rapid pace of adaptation reported from field and laboratory studies. Among the alternative explanations for rapid hopper adaptation are changes (genetic, epigenetic, or community structure) in endosymbiont communities that become advantageous for planthoppers and leafhoppers that feed on resistant rice varieties. This review examines the nature of these symbiont communities and their functions in planthoppers and leafhoppers—focusing on their likely roles in mediating adaptation to plant resistance. Evidence from a small number of experimental studies suggests that bacterial and eukaryotic (including yeast-like) symbionts can determine or mediate hopper virulence on rice plants and that symbiont functions could change over successive generations of selection on both resistant and susceptible plants. The review highlights the potential complexity of rice hopper–symbiont interactions and calls for a more careful choice of research materials and methods to help reduce this complexity. Finally, the consequences of symbiont-mediated virulence adaptation for future rice breeding programs are discussed.     

RUNAWAY COEVOLUTION: ADAPTATION TO HERITABLE AND NONHERITABLE ENVIRONMENTS

Populations evolve in response to the external environment, whether abiotic (e.g., climate) or biotic (e.g., other conspecifics). We investigated how adaptation to biotic, heritable environments differs from adaptation to abiotic, nonheritable environments. We found that, for the same selection coefficients, the coadaptive process between genes and heritable environments is much faster than genetic adaptation to an abiotic nonheritable environment. The increased rate of adaptation results from the positive association generated by reciprocal selection between the heritable environment and the genes responding to it. These associations result in a runaway process of adaptive coevolution, even when the genes creating the heritable environment and genes responding to the heritable environment are unlinked. Although tightening the degree of linkage accelerates the coadaptive process, the acceleration caused by a comparable amount of inbreeding is greater, because inbreeding has a cumulative effect on reducing functional recombination over generations. Our results suggest that that adaptation to local abiotic environmental variation may result in the rapid diversification of populations and subsequent reproductive isolation not directly but rather via its effects on heritable environments and the genes responding to them.     

Genomic Analysis of Natural Selection and Phenotypic Variation in High-Altitude Mongolians

Deedu (DU) Mongolians, who migrated from the Mongolian steppes to the Qinghai-Tibetan Plateau approximately 500 years ago, are challenged by environmental conditions similar to native Tibetan highlanders. Identification of adaptive genetic factors in this population could provide insight into coordinated physiological responses to this environment. Here we examine genomic and phenotypic variation in this unique population and present the first complete analysis of a Mongolian whole-genome sequence. High-density SNP array data demonstrate that DU Mongolians share genetic ancestry with other Mongolian as well as Tibetan populations, specifically in genomic regions related with adaptation to high altitude. Several selection candidate genes identified in DU Mongolians are shared with other Asian groups (e.g., EDAR), neighboring Tibetan populations (including high-altitude candidates EPAS1, PKLR, and CYP2E1), as well as genes previously hypothesized to be associated with metabolic adaptation (e.g., PPARG). Hemoglobin concentration, a trait associated with high-altitude adaptation in Tibetans, is at an intermediate level in DU Mongolians compared to Tibetans and Han Chinese at comparable altitude. Whole-genome sequence from a DU Mongolian (Tianjiao1) shows that about 2% of the genomic variants, including more than 300 protein-coding changes, are specific to this individual. Our analyses of DU Mongolians and the first Mongolian genome provide valuable insight into genetic adaptation to extreme environments.     

Reciprocal adaptation of rice and Xanthomonas oryzae pv. oryzae: cross-species 2D GWAS reveals the underlying genetics

A 1D/2D genome-wide association study strategy was adopted to investigate the genetic systems underlying the reciprocal adaptation of rice (Oryza sativa) and its bacterial pathogen, Xanthomonas oryzae pv. oryzae (Xoo) using the whole-genome sequencing and large-scale phenotyping data of 701 rice accessions and 23 diverse Xoo strains. Forty-seven Xoo virulence-related genes and 318 rice quantitative resistance genes (QR-genes) mainly located in 41 genomic regions, and genome-wide interactions between the detected virulence-related genes and QR genes were identified, including well-known resistance genes/virulence genes plus many previously uncharacterized ones. The relationship between rice and Xoo was characterized by strong differentiation among Xoo races corresponding to the subspecific differentiation of rice, by strong shifts toward increased resistance/virulence of rice/Xoo populations and by rich genetic diversity at the detected rice QR-genes and Xoo virulence genes, and by genome-wide interactions between many rice QR-genes and Xoo virulence genes in a multiple-to-multiple manner, presumably resulting either from direct protein–protein interactions or from genetic epistasis. The observed complex genetic interaction system between rice and Xoo likely exists in other crop–pathogen systems that would maintain high levels of diversity at their QR-loci/virulence-loci, resulting in dynamic coevolutionary consequences during their reciprocal adaptation.

A complex system of genetic interactions leads to reciprocal adaptation between rice and its bacterial pathogen, Xanthomonas oryzae pv. oryzae.     

Impact of genetic architecture on the relative rates of X versus autosomal adaptive substitution

Molecular evolutionary theory predicts that the ratio of autosomal to X-linked adaptive substitution (K(A)/K(x)) is primarily determined by the average dominance coefficient of beneficial mutations. Although this theory has profoundly influenced analysis and interpretation of comparative genomic data, its predictions are based upon two unverified assumptions about the genetic basis of adaptation. The theory assumes that 1) the rate of adaptively driven molecular evolution is limited by the availability of beneficial mutations, and 2) the scaling of evolutionary parameters between the X and the autosomes (e.g., the beneficial mutation rate, and the fitness effect distribution of beneficial alleles, per X-linked versus autosomal locus) is constant across molecular evolutionary timescales. Here, we show that the genetic architecture underlying bouts of adaptive substitution can influence both assumptions, and consequently, the theoretical relationship between K(A)/K(x) and mean dominance. Quantitative predictions of prior theory apply when 1) many genomically dispersed genes potentially contribute beneficial substitutions during individual steps of adaptive walks, and 2) the population beneficial mutation rate, summed across the set of potentially contributing genes, is sufficiently small to ensure that adaptive substitutions are drawn from new mutations rather than standing genetic variation. Current research into the genetic basis of adaptation suggests that both assumptions are plausibly violated. We find that the qualitative positive relationship between mean dominance and K(A)/K(x) is relatively robust to the specific conditions underlying adaptive substitution, yet the quantitative relationship between dominance and K(A)/K(x) is quite flexible and context dependent. This flexibility may partially account for the puzzlingly variable X versus autosome substitution patterns reported in the empirical evolutionary genomics literature. The new theory unites the previously separate analysis of adaptation using new mutations versus standing genetic variation and makes several useful predictions about the interaction between genetic architecture, evolutionary genetic constraints, and effective population size in determining the ratio of adaptive substitution between autosomal and X-linked genes.