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.     

Landscape genomics reveals genetic evidence of local adaptation in a widespread tree,the Chinese wingnut (Pterocarya stenoptera)

Understanding the genetic basis underpinning local adaptation is one of the fundamental issues in ecological and evolutionary biology. In this study, we investigated the genomic basis underlying local adaptation of the Chinese wingnut (Pterocarya stenoptera C. DC). Our population genomic analyses revealed nine spatial genetic clusters across the current distribution range of this species. Based on the assessment of genetic–environment association, we found that adaptive divergence of the P. stenoptera populations were mainly shaped by solar radiation during fruit development, temperature seasonality, annual temperature, precipitation, and air humidity. In particular, our genome-wide scanning identified a total of 801 candidate single nucleotide polymorphisms (SNPs) that are highly correlated with diverse environmental factors. Further functional annotation of the SNPs identified some candidate genes that are involved into temperature, water, and light adaptation. Taken together, our results suggest that natural selection during local adaptation has contributed to the success of survival to diverse heterogenous environmental conditions. Our study provides important insights into the fundamental knowledge of the genetic basis underlying the local adaptation of non-model species.     

Genomic Architecture of Rapid Parallel Adaptation to Fresh Water in a Wild Fish

Rapid adaptation to novel environments may drive changes in genomic regions through natural selection. However, the genetic architecture underlying these adaptive changes is still poorly understood. Using population genomic approaches, we investigated the genomic architecture that underlies rapid parallel adaptation of Coilia nasus to fresh water by comparing four freshwater-resident populations with their ancestral anadromous population. Linkage disequilibrium network analysis and population genetic analyses revealed two putative large chromosome inversions on LG6 and LG22, which were enriched for outlier loci and exhibited parallel association with freshwater adaptation. Drastic frequency shifts and elevated genetic differentiation were observed for the two chromosome inversions among populations, suggesting that both inversions would undergo divergent selection between anadromous and resident ecotypes. Enrichment analysis of genes within chromosome inversions showed significant enrichment of genes involved in metabolic process, immunoregulation, growth, maturation, osmoregulation, and so forth, which probably underlay differences in morphology, physiology and behavior between the anadromous and freshwater-resident forms. The availability of beneficial standing genetic variation, large optimum shift between marine and freshwater habitats, and high efficiency of selection with large population size could lead to the observed rapid parallel adaptive genomic change. We propose that chromosomal inversions might have played an important role during the evolution of rapid parallel ecological divergence in the face of environmental heterogeneity in C. nasus. Our study provides insights into the genomic basis of rapid adaptation of complex traits in novel habitats and highlights the importance of structural genomic variants in analyses of ecological adaptation.     

A Nonparametric Test Reveals Selection for Rapid Flowering in the Arabidopsis Genome

The detection of footprints of natural selection in genetic polymorphism data is fundamental to understanding the genetic basis of adaptation, and has important implications for human health. The standard approach has been to reject neutrality in favor of selection if the pattern of variation at a candidate locus was significantly different from the predictions of the standard neutral model. The problem is that the standard neutral model assumes more than just neutrality, and it is almost always possible to explain the data using an alternative neutral model with more complex demography. Today's wealth of genomic polymorphism data, however, makes it possible to dispense with models altogether by simply comparing the pattern observed at a candidate locus to the genomic pattern, and rejecting neutrality if the pattern is extreme. Here, we utilize this approach on a truly genomic scale, comparing a candidate locus to thousands of alleles throughout the Arabidopsis thaliana genome. We demonstrate that selection has acted to increase the frequency of early-flowering alleles at the vernalization requirement locus FRIGIDA. Selection seems to have occurred during the last several thousand years, possibly in response to the spread of agriculture. We introduce a novel test statistic based on haplotype sharing that embraces the problem of population structure, and so should be widely applicable.     

A Nonparametric Test Reveals Selection for Rapid Flowering in the Arabidopsis Genome

The detection of footprints of natural selection in genetic polymorphism data is fundamental to understanding the genetic basis of adaptation, and has important implications for human health. The standard approach has been to reject neutrality in favor of selection if the pattern of variation at a candidate locus was significantly different from the predictions of the standard neutral model. The problem is that the standard neutral model assumes more than just neutrality, and it is almost always possible to explain the data using an alternative neutral model with more complex demography. Today's wealth of genomic polymorphism data, however, makes it possible to dispense with models altogether by simply comparing the pattern observed at a candidate locus to the genomic pattern, and rejecting neutrality if the pattern is extreme. Here, we utilize this approach on a truly genomic scale, comparing a candidate locus to thousands of alleles throughout the Arabidopsis thaliana genome. We demonstrate that selection has acted to increase the frequency of early-flowering alleles at the vernalization requirement locus FRIGIDA. Selection seems to have occurred during the last several thousand years, possibly in response to the spread of agriculture. We introduce a novel test statistic based on haplotype sharing that embraces the problem of population structure, and so should be widely applicable.     

Partially repeatable genetic basis of benthic adaptation in threespine sticklebacks

The extent to which convergent adaptation to similar ecological niches occurs by a predictable genetic basis remains a fundamental question in biology. Threespine stickleback fish have undergone an adaptive radiation in which ancestral oceanic populations repeatedly colonized and adapted to freshwater habitats. In multiple lakes in British Columbia, two different freshwater ecotypes have evolved: a deep‐bodied benthic form adapted to forage near the lake substrate, and a narrow‐bodied limnetic form adapted to forage in open water. Here, we use genome‐wide linkage mapping in marine × benthic F2 genetic crosses to test the extent of shared genomic regions underlying benthic adaptation in three benthic populations. We identify at least 100 Quantitative Trait Loci (QTL) harboring genes influencing skeletal morphology. The majority of QTL (57%) are unique to one cross. However, four genomic regions affecting eight craniofacial and armor phenotypes are found in all three benthic populations. We find that QTL are clustered in the genome and overlapping QTL regions are enriched for genomic signatures of natural selection. These findings suggest that benthic adaptation has occurred via both parallel and nonparallel genetic changes.     

Genomics of local adaptation with gene flow

Gene flow is a fundamental evolutionary force in adaptation that is especially important to understand as humans are rapidly changing both the natural environment and natural levels of gene flow. Theory proposes a multifaceted role for gene flow in adaptation, but it focuses mainly on the disruptive effect that gene flow has on adaptation when selection is not strong enough to prevent the loss of locally adapted alleles. The role of gene flow in adaptation is now better understood due to the recent development of both genomic models of adaptive evolution and genomic techniques, which both point to the importance of genetic architecture in the origin and maintenance of adaptation with gene flow. In this review, we discuss three main topics on the genomics of adaptation with gene flow. First, we investigate selection on migration and gene flow. Second, we discuss the three potential sources of adaptive variation in relation to the role of gene flow in the origin of adaptation. Third, we explain how local adaptation is maintained despite gene flow: we provide a synthesis of recent genomic models of adaptation, discuss the genomic mechanisms and review empirical studies on the genomics of adaptation with gene flow. Despite predictions on the disruptive effect of gene flow in adaptation, an increasing number of studies show that gene flow can promote adaptation, that local adaptations can be maintained despite high gene flow, and that genetic architecture plays a fundamental role in the origin and maintenance of local adaptation with gene flow.     

Adaptation in the age of ecological genomics: insights from parallelism and convergence

Parallel phenotypic diversification in closely related species is a rigorous framework for testing the role of natural selection in evolution. Do parallel phenotypes always diversify by parallel genetic bases or does selection pave many alternative genomic routes to the same phenotypic ends? In this review, we show that the advent of next-generation sequencing technologies and the growing use of genomic approaches make it increasingly feasible to answer these fundamental questions using ecological and evolutionary 'non-model' populations of vertebrates in nature. While it is generally expected, and often observed, that closely related populations or species have parallel genetic bases to parallel phenotypes, exceptions are not rare and show that alternative genetic routes can result in similar phenotypes. Ultimately, this framework may illuminate the ecological conditions, evolutionary histories and genetic architectures that result in recurrent phenotypes and rapid adaptation.     

Population genomics of natural and experimental populations of guppies (Poecilia reticulata)

Convergent evolution represents one of the best lines of evidence for adaptation, but few cases of phenotypic convergence are understood at the genetic level. Guppies inhabiting the Northern Mountain Range of Trinidad provide a classic example of phenotypic convergent evolution, where adaptation to low or high predation environments has been found for a variety of traits. A major advantage of this system is the possibility of long‐term experimental studies in nature, including transplantation from high to low predation sites. We used genome scans of guppies from three natural high and low predation populations and from two experimentally established populations and their sources to examine whether phenotypic convergent evolution leaves footprints at the genome level. We used population‐genetic modelling approaches to reconstruct the demographic history and migration among sampled populations. Naturally colonized low predation populations had signatures of increased effective population size since colonization, while introduction populations had signatures of decreased effective population size. Only a small number of regions across the genome had signatures of selection in all natural populations. However, the two experimental populations shared many genomic regions under apparent selection, more than expected by chance. This overlap coupled with a population decrease since introduction provides evidence for convergent selection occurring in the two introduced populations. The lack of genetic convergence in the natural populations suggests that convergent evolution is lacking in these populations or that the effects of selection become difficult to detect after a long‐time period.     

Rapid parallel evolution of standing variation in a single,complex, genomic region is associated with life history in steelhead/rainbow trout

Rapid adaptation to novel environments may drive changes in genomic regions through natural selection. Such changes may be population-specific or, alternatively, may involve parallel evolution of the same genomic region in multiple populations, if that region contains genes or co-adapted gene complexes affecting the selected trait(s). Both quantitative and population genetic approaches have identified associations between specific genomic regions and the anadromous (steelhead) and resident (rainbow trout) life-history strategies of Oncorhynchus mykiss. Here, we use genotype data from 95 single nucleotide polymorphisms and show that the distribution of variation in a large region of one chromosome, Omy5, is strongly associated with life-history differentiation in multiple above-barrier populations of rainbow trout and their anadromous steelhead ancestors. The associated loci are in strong linkage disequilibrium, suggesting the presence of a chromosomal inversion or other rearrangement limiting recombination. These results provide the first evidence of a common genomic basis for life-history variation in O. mykiss in a geographically diverse set of populations and extend our knowledge of the heritable basis of rapid adaptation of complex traits in novel habitats.     

Gene expression analysis for the identification of selection and local adaptation in fishes

In recent years, variation in gene expression has been recognized as an important component of environmental adaptation in multiple model species, including a few fish species. There is, however, still little known about the genetic basis of adaptation in gene expression resulting from variation in the aquatic environment (e.g. temperature, salinity and oxygen) and the physiological effect and costs of such differences in gene expression. This review presents and discusses progress and pitfalls of applying gene expression analyses to fishes and suggests simple frameworks to get started with gene expression analysis. It is emphasized that well-planned gene expression studies can serve as an important tool for the identification of selection in local populations of fishes, even for non-traditional model species where limited genomic information is available. Recent studies focusing on gene expression variation among natural fish populations are reviewed, highlighting the latest applications that combine genetic evidence from neutral markers and gene expression data.     

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.     

Identifying loci with breeding potential across temperate and tropical adaptation via EigenGWAS and EnvGWAS

Understanding the genomic basis of adaptation in maize is important for gene discovery and the improvement of breeding germplasm, but much remains a mystery in spite of significant population genetics and archaeological research. Identifying the signals underpinning adaptation are challenging as adaptation often coincided with genetic drift, and the base genomic diversity of the species in massive. In this study, tGBS technology was used to genotype 1,143 diverse maize accessions including landraces collected from 20 countries and elite breeding lines of tropical lowland, highland, subtropical/midaltitude and temperate ecological zones. Based on 355,442 high‐quality single nucleotide polymorphisms, 13 genomic regions were detected as being under selection using the bottom‐up searching strategy, EigenGWAS. Of the 13 selection regions, 10 were first reported, two were associated with environmental parameters via EnvGWAS, and 146 genes were enriched. Combining large‐scale genomic and ecological data in this diverse maize panel, our study supports a polygenic adaptation model of maize and offers a framework to enhance our understanding of both the mechanistic basis and the evolutionary consequences of maize domestication and adaptation. The regions identified here are promising candidates for further, targeted exploration to identify beneficial alleles and haplotypes for deployment in maize breeding.     

Signatures of selection in the three‐spined stickleback along a small‐scale brackish water – freshwater transition zone

Local adaptation is often obvious when gene flow is impeded, such as observed at large spatial scales and across strong ecological contrasts. However, it becomes less certain at small scales such as between adjacent populations or across weak ecological contrasts, when gene flow is strong. While studies on genomic adaptation tend to focus on the former, less is known about the genomic targets of natural selection in the latter situation. In this study, we investigate genomic adaptation in populations of the three‐spined stickleback Gasterosteus aculeatus L. across a small‐scale ecological transition with salinities ranging from brackish to fresh. Adaptation to salinity has been repeatedly demonstrated in this species. A genome scan based on 87 microsatellite markers revealed only few signatures of selection, likely owing to the constraints that homogenizing gene flow puts on adaptive divergence. However, the detected loci appear repeatedly as targets of selection in similar studies of genomic adaptation in the three‐spined stickleback. We conclude that the signature of genomic selection in the face of strong gene flow is weak, yet detectable. We argue that the range of studies of genomic divergence should be extended to include more systems characterized by limited geographical and ecological isolation, which is often a realistic setting in nature.     

Evolution of polygenic traits under global vs local adaptation

Observations about the number, frequency, effect size, and genomic distribution of alleles associated with complex traits must be interpreted in light of evolutionary process. These characteristics, which constitute a trait’s genetic architecture, can dramatically affect evolutionary outcomes in applications from agriculture to medicine, and can provide a window into how evolution works. Here, I review theoretical predictions about the evolution of genetic architecture under spatially homogeneous, global adaptation as compared with spatially heterogeneous, local adaptation. Due to the tension between divergent selection and migration, local adaptation can favor “concentrated” genetic architectures that are enriched for alleles of larger effect, clustered in a smaller number of genomic regions, relative to expectations under global adaptation. However, the evolution of such architectures may be limited by many factors, including the genotypic redundancy of the trait, mutation rate, and temporal variability of environment. I review the circumstances in which predictions differ for global vs local adaptation and discuss where progress can be made in testing hypotheses using data from natural populations and lab experiments. As the field of comparative population genomics expands in scope, differences in architecture among traits and species will provide insights into how evolution works, and such differences must be interpreted in light of which kind of selection has been operating.     

A genomic selection component analysis characterizes migration‐selection balance

The genetic differentiation of populations in response to local selection pressures has long been studied by evolutionary biologists, but key details about the process remain obscure. How rapidly can local adaptation evolve, how extensive is the process across the genome, and how strong are the opposing forces of natural selection and gene flow? Here, we combine direct measurement of survival and reproduction with whole‐genome genotyping of a plant species (Mimulus guttatus) that has recently invaded a novel habitat (the Quarry population). We renovate the classic selection component method to accommodate genomic data and observe selection at SNPs throughout the genome. SNPs showing viability selection in Quarry exhibit elevated divergence from neighboring populations relative to neutral SNPs. We also find that nonsignificant SNPs exhibit a subtle, but still significant, change in allele frequency toward neighboring populations, a predicted effect of gene flow. Given that the Quarry population is most probably only 30–40 generations old, the alleles conferring local advantage are almost certainly older than the population itself. Thus, local adaptation owes to the recruitment of standing genetic variation.     

Genome‐wide signature of local adaptation linked to variable CpG methylation in oak populations

It has long been known that adaptive evolution can occur through genetic mutations in DNA sequence, but it is unclear whether adaptive evolution can occur through analogous epigenetic mechanisms, such as through DNA methylation. If epigenetic variation contributes directly to evolution, species under threat of disease, invasive competition, climate change or other stresses would have greater stores of variation from which to draw. We looked for evidence of natural selection acting on variably methylated DNA sites using population genomic analysis across three climatologically distinct populations of valley oaks. We found patterns of genetic and epigenetic differentiations that indicate local adaptation is operating on large portions of the oak genome. While CHG methyl polymorphisms are not playing a significant role and would make poor targets for natural selection, our findings suggest that CpG methyl polymorphisms as a whole are involved in local adaptation, either directly or through linkage to regions under selection.     

Habitat-specific natural selection at a flowering-time QTL is a main driver of local adaptation in two wild barley populations

Understanding the genetic basis of local adaptation requires insight in the fitness effects of individual loci under natural field conditions. While rapid progress is made in the search for genes that control differences between plant populations, it is typically unknown whether the genes under study are in fact key targets of habitat-specific natural selection. Using a quantitative trait loci (QTL) approach, we show that a QTL associated with flowering-time variation between two locally adapted wild barley populations is an important determinant of fitness in one, but not in the other population's native habitat. The QTL mapped to the same position as a habitat-specific QTL for field fitness that affected plant reproductive output in only one of the parental habitats, indicating that the genomic region is under differential selection between the native habitats. Consistent with the QTL results, phenotypic selection of flowering time differed between the two environments, whereas other traits (growth rate and seed weight) were under selection but experienced no habitat-specific differential selection. This implies the flowering-time QTL as a driver of adaptive population divergence. Our results from phenotypic selection and QTL analysis are consistent with local adaptation without genetic trade-offs in performance across environments, i.e. without alleles or traits having opposing fitness effects in contrasting environments.     

A genome scan to detect candidate regions influenced by local natural selection in human populations

As human populations dispersed throughout the world, they were subjected to new selective forces, which must have led to local adaptation via natural selection and hence altered patterns of genetic variation. Yet, there are very few examples known in which such local selection has clearly influenced human genetic variation. A potential approach for detecting local selection is to screen random loci across the genome; those loci that exhibit unusually large genetic distances between human populations are then potential markers of genomic regions under local selection. We investigated this approach by genotyping 332 short tandem repeat (STR) loci in Africans and Europeans and calculating the genetic differentiation for each locus. Patterns of genetic diversity at these loci were consistent with greater variation in Africa and with local selection operating on populations as they moved out of Africa. For 11 loci exhibiting the largest genetic differences, we genotyped an additional STR locus located nearby; the genetic distances for these nearby loci were significantly larger than average. These genomic regions therefore reproducibly exhibit larger genetic distances between populations than the "average" genomic region, consistent with local selection. Our results demonstrate that genome scans are a promising means of identifying candidate regions that have been subjected to local selection.     

CONVERGENCE AND DIVERGENCE DURING THE ADAPTATION TO SIMILAR ENVIRONMENTS BY AN AUSTRALIAN GROUNDSEL

Adaptation to replicate environments is often achieved through similar phenotypic solutions. Whether selection also produces convergent genomic changes in these situations remains largely unknown. The variable groundsel, Senecio lautus, is an excellent system to investigate the genetic underpinnings of convergent evolution, because morphologically similar forms of these plants have adapted to the same environments along the coast of Australia. We compared range‐wide patterns of genomic divergence in natural populations of this plant and searched for regions putatively affected by natural selection. Our results indicate that environmental adaptation followed complex genetic trajectories, affecting multiple loci, implying both the parallel recruitment of the same alleles and the divergence of completely different genomic regions across geography. An analysis of the biological functions of candidate genes suggests that adaptation to coastal environments may have occurred through the recruitment of different genes participating in similar processes. The relatively low genetic convergence that characterizes the parallel evolution of S. lautus forms suggests that evolution is more constrained at higher levels of biological organization.