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.     

Enhanced AFLP genome scans detect local adaptation in high-altitude populations of a small rodent (Microtus arvalis)

Adaptation to adverse environmental conditions such as high altitude requires physiological and/or morphological changes. Genome scans provide a means to identify the genetic basis of such adaptations without previous knowledge about the particular genetic variants or traits under selection. In this study, we scanned 3027 amplified fragment length polymorphisms (AFLP) in four populations of the common vole Microtus arvalis for loci associated with local adaptation and high altitude. We investigated voles from two populations at high elevation (~2000 m a.s.l.) representing the upper limit of the altitudinal distribution of the species and two geographically close low-altitude populations (<600 m a.s.l.). Statistical analysis incorporated a new Bayesian F(ST) outlier approach specifically developed for AFLP markers, which considers the intensity of AFLP bands instead of mere presence/absence and allows to derive population-based estimates of allele frequencies and F(IS) values. Computer simulations showed that this approach increases the statistical power of the detection of AFLP markers under selection almost to the power of single nucleotide polymorphism (SNP) data without compromising specificity. Our enhanced genome scan resulted in 20 prime candidate markers for positive selection, which show mostly extremely high allele frequency differences between the low- and high-altitude populations. The comparison of global- and pairwise-enhanced genome scans demonstrated further that very strong selective signatures may also be associated with single populations suggesting the importance of local adaptation in alpine populations of common voles.     

Grasshopper ontogeny in relation to time constraints: adaptive divergence and stasis

1. Life history theory generally predicts a trade-off between shortjuvenile development and large adult size, assuming invariant growth rates within species. This pivotal assumption has been explicitly tested in few organisms. 2. We studied ontogeny in 13 populations of Omocestus viridulus grasshoppers under common garden conditions. High-altitude populations, facing short growing seasons and thus seasonal time constraints, were found to grow at a similar rate to low altitude conspecifics. 3. Instead, high-altitude grasshoppers evolved faster development, and the correlated change in body size led to an altitudinal size cline mediating a trade-off with female fecundity. 4. An additional juvenile stage occurred in low- but not high-altitude females. This difference is probably due to the evolution of lowered critical size thresholds in high-altitude grasshoppers to accelerate development. 5. We found a strikingly lower growth rate in males than females that we interpret as the outcome of concurrent selection for protandry and small male size. 6. Within populations, large individuals developed faster than small individuals, suggesting within-population genetic variation in growth rates. 7. We provide evidence that different time constraints (seasonal, protandry selection) can lead to different evolutionary responses in intrinsic growth, and that correlations among ontogenetic traits within populations cannot generally be used to predict life history adaptation among populations. Moreover, our study illustrates that comparisons of ontogenetic patterns can shed light on the developmental basis underlying phenotypic evolution.     

Broad Concordance in the Spatial Distribution of Adaptive and Neutral Genetic Variation across an Elevational Gradient in Deer Mice

When species are continuously distributed across environmental gradients, the relative strength of selection and gene flow shape spatial patterns of genetic variation, potentially leading to variable levels of differentiation across loci. Determining whether adaptive genetic variation tends to be structured differently than neutral variation along environmental gradients is an open and important question in evolutionary genetics. We performed exome-wide population genomic analysis on deer mice sampled along an elevational gradient of nearly 4,000 m of vertical relief. Using a combination of selection scans, genotype−environment associations, and geographic cline analyses, we found that a large proportion of the exome has experienced a history of altitude-related selection. Elevational clines for nearly 30% of these putatively adaptive loci were shifted significantly up- or downslope of clines for loci that did not bear similar signatures of selection. Many of these selection targets can be plausibly linked to known phenotypic differences between highland and lowland deer mice, although the vast majority of these candidates have not been reported in other studies of highland taxa. Together, these results suggest new hypotheses about the genetic basis of physiological adaptation to high altitude, and the spatial distribution of adaptive genetic variation along environmental gradients.     

Widespread Signals of Convergent Adaptation to High Altitude in Asia and America

Living at high altitude is one of the most difficult challenges that humans had to cope with during their evolution. Whereas several genomic studies have revealed some of the genetic bases of adaptations in Tibetan, Andean, and Ethiopian populations, relatively little evidence of convergent evolution to altitude in different continents has accumulated. This lack of evidence can be due to truly different evolutionary responses, but it can also be due to the low power of former studies that have mainly focused on populations from a single geographical region or performed separate analyses on multiple pairs of populations to avoid problems linked to shared histories between some populations. We introduce here a hierarchical Bayesian method to detect local adaptation that can deal with complex demographic histories. Our method can identify selection occurring at different scales, as well as convergent adaptation in different regions. We apply our approach to the analysis of a large SNP data set from low- and high-altitude human populations from America and Asia. The simultaneous analysis of these two geographic areas allows us to identify several candidate genome regions for altitudinal selection, and we show that convergent evolution among continents has been quite common. In addition to identifying several genes and biological processes involved in high-altitude adaptation, we identify two specific biological pathways that could have evolved in both continents to counter toxic effects induced by hypoxia.     

A Chromosome-level Genome Assembly of Wild Castor Provides New Insights into its Adaptive Evolution in Tropical Desert

Wild castor grows in the high-altitude tropical desert of the African Plateau,a region known for high ultraviolet radiation,strong light,and extremely dry condition.To investigate the potential genetic basis of adaptation to both highland and tropical deserts,we generated a chromosome-level genome sequence assembly of the wild castor accession WT05,with a genome size of 316 Mb,a scaffold N50 of 31.93 Mb,and a contig N50 of 8.96 Mb,respectively.Compared with cultivated castor and other Euphorbiac...     

Explorative genome scan to detect candidate loci for adaptation along a gradient of altitude in the common frog (Rana temporaria)

Today, with the rapid development of population genomics, the genetic basis of adaptation can be unraveled directly at the genome level, without any prerequisites about the selectively advantageous genes or traits. For nonmodel species, it is now possible to screen many markers randomly scattered across the genome and to distinguish between the neutral genetic background and outlier loci displaying an atypical behavior (e.g., a higher differentiation between populations). This study investigated the genetic frame of adaptation to a gradient of altitude in the common frog (Rana temporaria) by means of a genome scan based on 392 amplified fragment length polymorphism markers. Using two outlier detection methods never applied to dominant data so far, we sought for loci with a genetic differentiation diverging from neutral expectations when comparing populations from different altitudes. All the detected loci were sorted out according to their most probable cause for outlier behavior and classified as false positives, outliers due to local effects, or outliers associated with altitude. Altogether, eight good candidate loci were identified as potentially involved in adaptation to altitude because they were picked out in several independent interaltitude comparisons. This result illustrated the potential of genome-wide surveys to reveal selection signatures along selection gradients, where the association between environmental variables and fitness-related traits may be complex and/or cryptic. In this article, we also underlined the need for confirmation of the selection footprints for the outlier loci. Finally, we provided some preliminary insights into the genetic basis of adaptation along an altitudinal cline in the common frog.     

Introgression,admixture, and selection facilitate genetic adaptation to high-altitude environments in cattle

Domestication and subsequent selection of cattle to form breeds and biological types that can adapt to different environments partitioned ancestral genetic diversity into distinct modern lineages. Genome-wide selection particularly for adaptation to extreme environments left detectable signatures genome-wide. We used high-density genotype data for 42 cattle breeds and identified the influence of Bos grunniens and Bos javanicus on the formation of Chinese indicine breeds that led to their divergence from India-origin zebu. We also found evidence for introgression, admixture, and migration in most of the Chinese breeds. Selection signature analyses between high-altitude (≥1800 m) and low-altitude adapted breeds (<1500 m) revealed candidate genes (ACSS2, ALDOC, EPAS1, EGLN1, NUCB2) and pathways that are putatively involved in hypoxia adaptation. Immunohistochemical, real-time PCR and CRISPR/cas9 ACSS2-knockout analyses suggest that the up-regulation of ACSS2 expression in the liver promotes the metabolic adaptation of cells to hypoxia via the hypoxia-inducible factor pathway. High altitude adaptation involved the introgression of alleles from high-altitude adapted yaks into Chinese Bos taurus taurus prior to their formation into recognized breeds and followed by selection. In addition to selection, adaptation to high altitude environments has been facilitated by admixture and introgression with locally adapted cattle populations.     

Whole‐genome resequencing uncovers molecular signatures of natural and sexual selection in wild bighorn sheep

The identification of genes influencing fitness is central to our understanding of the genetic basis of adaptation and how it shapes phenotypic variation in wild populations. Here, we used whole‐genome resequencing of wild Rocky Mountain bighorn sheep (Ovis canadensis) to >50‐fold coverage to identify 2.8 million single nucleotide polymorphisms (SNPs) and genomic regions bearing signatures of directional selection (i.e. selective sweeps). A comparison of SNP diversity between the X chromosome and the autosomes indicated that bighorn males had a dramatically reduced long‐term effective population size compared to females. This probably reflects a long history of intense sexual selection mediated by male–male competition for mates. Selective sweep scans based on heterozygosity and nucleotide diversity revealed evidence for a selective sweep shared across multiple populations at RXFP2, a gene that strongly affects horn size in domestic ungulates. The massive horns carried by bighorn rams appear to have evolved in part via strong positive selection at RXFP2. We identified evidence for selection within individual populations at genes affecting early body growth and cellular response to hypoxia; however, these must be interpreted more cautiously as genetic drift is strong within local populations and may have caused false positives. These results represent a rare example of strong genomic signatures of selection identified at genes with known function in wild populations of a nonmodel species. Our results also showcase the value of reference genome assemblies from agricultural or model species for studies of the genomic basis of adaptation in closely related wild taxa.     

Genomic insights into adaptation to high-altitude environments

Elucidating the molecular genetic basis of adaptive traits is a central goal of evolutionary genetics. The cold, hypoxic conditions of high-altitude habitats impose severe metabolic demands on endothermic vertebrates, and understanding how high-altitude endotherms cope with the combined effects of hypoxia and cold can provide important insights into the process of adaptive evolution. The physiological responses to high-altitude stress have been the subject of over a century of research, and recent advances in genomic technologies have opened up exciting opportunities to explore the molecular genetic basis of adaptive physiological traits. Here, we review recent literature on the use of genomic approaches to study adaptation to high-altitude hypoxia in terrestrial vertebrates, and explore opportunities provided by newly developed technologies to address unanswered questions in high-altitude adaptation at a genomic scale.     

EP300 contributes to high-altitude adaptation in Tibetans by regulating nitric oxide production

The genetic adaptation of Tibetans to high altitude hypoxia likely involves a group of genes in the hypoxic pathway,as suggested by earlier studies.To test the adaptive role of the previously reported candidate gene EP300 (histone acetyltransferase p300),we conducted resequencing of a 108.9 kb gene region of EP300 in 80 unrelated Tibetans.The allele-frequency and haplotype-based neutrality tests detected signals of positive Darwinian selection on EP300 in Tibetans,with a group of variants showing allelic divergence between Tibetans and lowland reference populations,including Han Chinese,Europeans,and Africans.Functional prediction suggested the involvement of multiple EP300 variants in gene expression regulation.More importantly,genetic association tests in 226 Tibetans indicated significant correlation of the adaptive EP300 variants with blood nitric oxide (NO) concentration.Collectively,we propose that EP300 harbors adaptive variants in Tibetans,which might contribute to high-altitude adaptation through regulating NO production.     

Altitudinal gradients,plant hybrid zones and evolutionary novelty

Altitudinal gradients are characterized by steep changes of the physical and biotic environment that present challenges to plant adaptation throughout large parts of the world. Hybrid zones may form where related species inhabit different neighbouring altitudes and can facilitate interspecific gene flow and potentially the breakdown of species barriers. Studies of such hybrid zones can reveal much about the genetic basis of adaptation to environmental differences stemming from changes in altitude and the maintenance of species divergence in the face of gene flow. Furthermore, owing to recombination and transgressive effects, such hybrid zones can be sources of evolutionary novelty. We document plant hybrid zones associated with altitudinal gradients and emphasize similarities and differences in their structure. We then focus on recent studies of a hybrid zone between two Senecio species that occur at high and low altitude on Mount Etna, Sicily, showing how adaptation to local environments and intrinsic selection against hybrids act to maintain it. Finally, we consider the potential of altitudinal hybrid zones for generating evolutionary novelty through adaptive introgression and hybrid speciation. Examples of homoploid hybrid species of Senecio and Pinus that originated from altitudinal hybrid zones are discussed.     

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.     

Adaptive introgression of the beta-globin cluster in two Andean waterfowl

Introgression of beneficial alleles has emerged as an important avenue for genetic adaptation in both plant and animal populations. In vertebrates, adaptation to hypoxic high-altitude environments involves the coordination of multiple molecular and cellular mechanisms, including selection on the hypoxia-inducible factor (HIF) pathway and the blood-O₂ transport protein hemoglobin (Hb). In two Andean duck species, a striking DNA sequence similarity reflecting identity by descent is present across the ~20 kb β-globin cluster including both embryonic (HBE) and adult (HBB) paralogs, though it was yet untested whether this is due to independent parallel evolution or adaptive introgression. In this study, we find that identical amino acid substitutions in the β-globin cluster that increase Hb-O₂ affinity have likely resulted from historical interbreeding between high-altitude populations of two different distantly-related species. We examined the direction of introgression and discovered that the species with a deeper mtDNA divergence that colonized high altitude earlier in history (Anas flavirostris) transferred adaptive genetic variation to the species with a shallower divergence (A. georgica) that likely colonized high altitude more recently possibly following a range shift into a novel environment. As a consequence, the species that received these β-globin variants through hybridization might have adapted to hypoxic conditions in the high-altitude environment more quickly through acquiring beneficial alleles from the standing, hybrid-origin variation, leading to faster evolution.Subject terms: Evolutionary genetics, Evolutionary biology     

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.     

鹿科动物基因组学研究进展

鹿科动物是世界上最丰富的大型哺乳动物之一,在极北地区、热带地区和高海拔地区都有分布.中国占世界鹿科动物40％以上,是鹿科动物进化的主战场.鹿科动物除了具有反刍动物常见的独特表型特征外,更是进化出周期性再生新器官——鹿茸角.鹿科动物是研究生态学、行为学和进化生物学非常有价值的动物模型,特别是在研究哺乳动物器官再生方面具有...     

Selection Signatures in Worldwide Sheep Populations

The diversity of populations in domestic species offers great opportunities to study genome response to selection. The recently published Sheep HapMap dataset is a great example of characterization of the world wide genetic diversity in sheep. In this study, we re-analyzed the Sheep HapMap dataset to identify selection signatures in worldwide sheep populations. Compared to previous analyses, we made use of statistical methods that (i) take account of the hierarchical structure of sheep populations, (ii) make use of linkage disequilibrium information and (iii) focus specifically on either recent or older selection signatures. We show that this allows pinpointing several new selection signatures in the sheep genome and distinguishing those related to modern breeding objectives and to earlier post-domestication constraints. The newly identified regions, together with the ones previously identified, reveal the extensive genome response to selection on morphology, color and adaptation to new environments.     

Genomewide scan for adaptive differentiation along altitudinal gradient in the Andrew's toad Bufo andrewsi

Recent studies of humans, dogs and rodents have started to discover the genetic underpinnings of high altitude adaptations, yet amphibians have received little attention in this respect. To identify possible signatures of adaptation to altitude, we performed a genome scan of 15 557 single nucleotide polymorphisms (SNPs) obtained with restriction site‐associated DNA sequencing of pooled samples from 11 populations of Andrew's toad (Bufo andrewsi) from the edge of the Tibetan Plateau, spanning an altitudinal gradient from 1690 to 2768 m.a.s.l. We discovered significant geographic differentiation among all sites, with an average F_ST = 0.023 across all SNPs. Apart from clear patterns of isolation by distance, we discovered numerous outlier SNPs showing strong associations with variation in altitude (1394 SNPs), average annual temperature (1859 SNPs) or both (1051 SNPs). Levels and patterns of genetic differentiation in these SNPs were consistent with the hypothesis that they have been subject to directional selection and reflect adaptation to altitudinal variation among the study sites. Genes with footprints of selection were significantly enriched in binding and metabolic processes. Several genes potentially related to high altitude adaptation were identified, although the identity and functional significance of most genomic targets of selection remain unknown. In general, the results provide genomic support for results of earlier common garden and low coverage genetic studies that have uncovered substantial adaptive differentiation along altitudinal and latitudinal gradients in amphibians.     

Genetic adaptation to high altitude in the Ethiopian highlands

Background

Genomic analysis of high-altitude populations residing in the Andes and Tibet has revealed several candidate loci for involvement in high-altitude adaptation, a subset of which have also been shown to be associated with hemoglobin levels, including EPAS1, EGLN1, and PPARA, which play a role in the HIF-1 pathway. Here, we have extended this work to high- and low-altitude populations living in Ethiopia, for which we have measured hemoglobin levels. We genotyped the Illumina 1M SNP array and employed several genome-wide scans for selection and targeted association with hemoglobin levels to identify genes that play a role in adaptation to high altitude.

Results

We have identified a set of candidate genes for positive selection in our high-altitude population sample, demonstrated significantly different hemoglobin levels between high- and low-altitude Ethiopians and have identified a subset of candidate genes for selection, several of which also show suggestive associations with hemoglobin levels.

Conclusions

We highlight several candidate genes for involvement in high-altitude adaptation in Ethiopia, including CBARA1, VAV3, ARNT2 and THRB. Although most of these genes have not been identified in previous studies of high-altitude Tibetan or Andean population samples, two of these genes (THRB and ARNT2) play a role in the HIF-1 pathway, a pathway implicated in previous work reported in Tibetan and Andean studies. These combined results suggest that adaptation to high altitude arose independently due to convergent evolution in high-altitude Amhara populations in Ethiopia.     

Adaptive responses for seed and leaf phenology in natural populations of sessile oak along an altitudinal gradient

We assessed the adaptive potential of seed and leaf phenology in 10 natural populations of sessile oak (Quercus petraea) sampled along two altitudinal transects using common garden experiments. Population differentiation for both phenological traits was observed with high-altitude populations germinating and flushing later than low altitude ones. However, high genetic variation and heritability values were also maintained within populations, despite slightly decreasing for dates of leaf unfolding with increasing altitude. We suggest that biotic and abiotic fluctuating selection pressures within populations and high gene flow are the main mechanisms maintaining high genetic variation for these fitness related traits. Moreover, changes in selection intensity and/or selection pressures along the altitudinal gradient can explain the reduction in genetic variation observed for leaf phenology. We anticipate that the maintenance of high genetic variation will be a valuable resource for future adaptation of sessile oak populations undergoing an upslope shift caused by climate change.