A new method to uncover signatures of divergent and stabilizing selection in quantitative traits

Error rates in QST-FST comparisons depend on genetic architecture and estimation procedures

Genetic and phenotypic variation among populations is one of the fundamental subjects of evolutionary genetics. One question that arises often in data on natural populations is whether differentiation among populations on a particular trait might be caused in part by natural selection. For the past several decades, researchers have used QST-FST approaches to compare the amount of trait differentiation among populations on one or more traits (measured by the statistic QST) with differentiation on genome-wide genetic variants (measured by FST). Theory says that under neutrality, FST and QST should be approximately equal in expectation, so QST values much larger than FST are consistent with local adaptation driving subpopulations' trait values apart, and QST values much smaller than FST are consistent with stabilizing selection on similar optima. At the same time, investigators have differed in their definitions of genome-wide FST (such as "ratio of averages" vs. "average of ratios" versions of FST) and in their definitions of the variance components in QST. Here, we show that these details matter. Different versions of FST and QST have different interpretations in terms of coalescence time, and comparing incompatible statistics can lead to elevated type I error rates, with some choices leading to type I error rates near one when the nominal rate is 5%. We conduct simulations under varying genetic architectures and forms of population structure and show how they affect the distribution of QST. When many loci influence the trait, our simulations support procedures grounded in a coalescent-based framework for neutral phenotypic differentiation.

Detecting diversifying selection for a trait from within and between-species genotypes and phenotypes

To quantify selection acting on a trait, methods have been developed using either within or between-species variation. However, methods using within-species variation do not integrate the changes at the macro-evolutionary scale. Conversely, current methods using between-species variation usually discard within-species variation, thus not accounting for processes at the micro-evolutionary scale. The main goal of this study is to define a neutrality index for a quantitative trait, by combining within- and between-species variation. This neutrality index integrates nucleotide polymorphism and divergence for normalizing trait variation. As such, it does not require estimation of population size nor of time of speciation for normalization. Our index can be used to seek deviation from the null model of neutral evolution, and test for diversifying selection. Applied to brain mass and body mass at the mammalian scale, we show that brain mass is under diversifying selection. Finally, we show that our test is not sensitive to the assumption that population sizes, mutation rates and generation time are constant across the phylogeny, and automatically adjust for it.

Error rates in Q ST -F ST comparisons depend on genetic architecture and estimation procedures

Genetic and phenotypic variation among populations is one of the fundamental subjects of evolutionary genetics. One question that arises often in data on natural populations is whether differentiation among populations on a particular trait might be caused in part by natural selection. For the past several decades, researchers have usedQ S T - F S T approaches to compare the amount of trait differentiation among populations on one or more traits (measured by the statisticQ S T ) with differentiation on genome-wide genetic variants (measured byF S T ). Theory says that under neutrality,F S T andQ S T should be approximately equal in expectation, soQ S T values much larger thanF S T are consistent with local adaptation driving subpopulations' trait values apart, andQ S T values much smaller thanF S T are consistent with stabilizing selection on similar optima. At the same time, investigators have differed in their definitions of genome-wideF S T (such as "ratio of averages" vs. "average of ratios" versions ofF S T ) and in their definitions of the variance components inQ S T . Here, we show that these details matter. Different versions ofF S T andQ S T have different interpretations in terms of coalescence time, and comparing incompatible statistics can lead to elevated type I error rates, with some choices leading to type I error rates near one when the nominal rate is 5%. We conduct simulations under varying genetic architectures and forms of population structure and show how they affect the distribution ofQ S T . When many loci influence the trait, our simulations support procedures grounded in a coalescent-based framework for neutral phenotytpic differentiation.

Ecological trade-offs drive phenotypic and genetic differentiation of Arabidopsis thaliana in Europe

Plant diversity is shaped by trade-offs between traits related to competitive ability, propagule dispersal, and stress resistance. However, we still lack a clear understanding of how these trade-offs influence species distribution and population dynamics. In Arabidopsis thaliana, recent genetic analyses revealed a group of cosmopolitan genotypes that successfully recolonized Europe from its center after the last glaciation, excluding older (relict) lineages from the distribution except for their north and south margins. Here, we tested the hypothesis that cosmopolitans expanded due to higher colonization ability, while relicts persisted at the margins due to higher tolerance to competition and/or stress. We compared the phenotypic and genetic differentiation between 71 European genotypes originating from the center, and the south and north margins. We showed that a trade-off between plant fecundity and seed mass shapes the differentiation of A. thaliana in Europe, suggesting that the success of the cosmopolitan groups could be explained by their high dispersal ability. However, at both north and south margins, we found evidence of selection for alleles conferring low dispersal but highly competitive and stress-resistance abilities. This study sheds light on the role of ecological trade-offs as evolutionary drivers of the distribution and dynamics of plant populations.

Resource Variation Within and Between Patches: Where Exploitation Competition, Local Adaptation, and Kin Selection Meet

AbstractIn patch- or habitat-structured populations, different processes can favor adaptive polymorphism at different scales. While spatial heterogeneity can generate spatially disruptive selection favoring variation between patches, local competition can lead to locally disruptive selection promoting variation within patches. So far, almost all theory has studied these two processes in isolation. Here, we use mathematical modeling to investigate how resource variation within and between habitats influences the evolution of variation in a consumer population where individuals compete in finite patches connected by dispersal. We find that locally and spatially disruptive selection typically act in concert, favoring polymorphism under a wider range of conditions than when in isolation. But when patches are small and dispersal between them is low, kin competition inhibits the emergence of polymorphism, especially when the latter is driven by local competition for resources. We further use our model to clarify what comparisons between trait and neutral genetic differentiation (Q ST / F ST comparisons) can tell about the nature of selection. Overall, our results help us understand the interaction between two major drivers of polymorphism: locally and spatially disruptive selection, and how this interaction is modulated by the unavoidable effects of kin selection under limited dispersal.

An allele-sharing, moment-based estimator of global, population-specific and population-pair FST under a general model of population structure

Being able to properly quantify genetic differentiation is key to understanding the evolutionary potential of a species. One central parameter in this context is FST, the mean coancestry within populations relative to the mean coancestry between populations. Researchers have been estimating FST globally or between pairs of populations for a long time. More recently, it has been proposed to estimate population-specific FST values, and population-pair mean relative coancestry. Here, we review the several definitions and estimation methods of FST, and stress that they provide values relative to a reference population. We show the good statistical properties of an allele-sharing, method of moments based estimator of FST (global, population-specific and population-pair) under a very general model of population structure. We point to the limitation of existing likelihood and Bayesian estimators when the populations are not independent. Last, we show that recent attempts to estimate absolute, rather than relative, mean coancestry fail to do so.

A review on Q ST-F ST comparisons of seed plants: Insights for conservation

Increased access to genome-wide data provides new opportunities for plant conservation. However, information on neutral genetic diversity in a small number of marker loci can still be valuable because genomic data are not available to most rare plant species. In the hope of bridging the gap between conservation science and practice, we outline how conservation practitioners can more efficiently employ population genetic information in plant conservation. We first review the current knowledge about neutral genetic variation (NGV) and adaptive genetic variation (AGV) in seed plants, regarding both within-population and among-population components. We then introduce the estimates of among-population genetic differentiation in quantitative traits (Q ST) and neutral markers (F ST) to plant biology and summarize conservation applications derived from Q ST-F ST comparisons, particularly on how to capture most AGV and NGV on both in-situ and ex-situ programs. Based on a review of published studies, we found that, on average, two and four populations would be needed for woody perennials (n = 18) to capture 99% of NGV and AGV, respectively, whereas four populations would be needed in case of herbaceous perennials (n = 14). On average, Q ST is about 3.6, 1.5, and 1.1 times greater than F ST in woody plants, annuals, and herbaceous perennials, respectively. Hence, conservation and management policies or suggestions based solely on inference on F ST could be misleading, particularly in woody species. To maximize the preservation of the maximum levels of both AGV and NGV, we suggest using maximum Q ST rather than average Q ST. We recommend conservation managers and practitioners consider this when formulating further conservation and restoration plans for plant species, particularly woody species.

Genomic data and common garden experiments reveal climate-driven selection on ecophysiological traits in two Mediterranean oaks

Improving our knowledge of how past climate-driven selection has acted on present-day trait population divergence is essential to understand local adaptation processes and improve our predictions of evolutionary trajectories in the face of altered selection pressures resulting from climate change. In this study, we investigated signals of selection on traits related to drought tolerance and growth rates in two Mediterranean oak species (Quercus faginea and Q. lusitanica) with contrasting distribution ranges and climatic niches. We genotyped 182 individuals from 24 natural populations of the two species using restriction-site-associated DNA sequencing and conducted a thorough functional characterization in 1602 seedlings from 21 populations cultivated in common garden experiments under contrasting watering treatments. Our genomic data revealed that both Q. faginea and Q. lusitanica have very weak population genetic structure, probably as a result of high rates of pollen-mediated gene flow among populations and large effective population sizes. In contrast, common garden experiments showed evidence of climate-driven divergent selection among populations on traits related to leaf morphology, physiology and growth in both species. Overall, our study suggests that climate is an important selective factor for Mediterranean oaks and that ecophysiological traits have evolved in drought-prone environments even in a context of very high rates of gene flow among populations.

Neutral and adaptive genetic diversity in plants: An overview

Genetic diversity is a prerequisite for evolutionary change in all kinds of organisms. It is generally acknowledged that populations lacking genetic variation are unable to evolve in response to new environmental conditions (e.g., climate change) and thus may face an increased risk of extinction. Although the importance of incorporating genetic diversity into the design of conservation measures is now well understood, less attention has been paid to the distinction between neutral (NGV) and adaptive (AGV) genetic variation. In this review, we first focus on the utility of NGV by examining the ways to quantify it, reviewing applications of NGV to infer ecological and evolutionary processes, and by exploring its utility in designing conservation measures for plant populations and species. Against this background, we then summarize the ways to identify and estimate AGV and discuss its potential use in plant conservation. After comparing NGV and AGV and considering their pros and cons in a conservation context, we conclude that there is an urgent need for a better understanding of AGV and its role in climate change adaptation. To date, however, there are only a few AGV studies on non-model plant species aimed at deciphering the genetic and genomic basis of complex trait variation. Therefore, conservation researchers and practitioners should keep utilizing NGV to develop relevant strategies for rare and endangered plant species until more estimates of AGV are available.

Evidence of climate-driven selection on tree traits and trait plasticity across the climatic range of a riparian foundation species

Selection on quantitative traits by heterogeneous climatic conditions can lead to substantial trait variation across a species range. In the context of rapidly changing environments, however, it is equally important to understand selection on trait plasticity. To evaluate the role of selection in driving divergences in traits and their associated plasticities within a widespread species, we compared molecular and quantitative trait variation in Populus fremontii (Fremont cottonwood), a foundation riparian distributed throughout Arizona. Using SNP data and genotypes from 16 populations reciprocally planted in three common gardens, we first performed Q_ST -F_ST analyses to detect selection on traits and trait plasticity. We then explored the environmental drivers of selection using trait-climate and plasticity-climate regressions. Three major findings emerged: 1) There was significant genetic variation in traits expressed in each of the common gardens and in the phenotypic plasticity of traits across gardens, both of which were heritable. 2) Based on Q_ST -F_ST comparisons, there was evidence of selection in all traits measured; however, this result varied from no effect in one garden to highly significant in another, indicating that detection of past selection is environmentally dependent. We also found strong evidence of divergent selection on plasticity across environments for two traits. 3) Traits and/or their plasticity were often correlated with population source climate (R² up to 0.77 and 0.66, respectively). These results suggest that steep climate gradients across the Southwest have played a major role in shaping the evolution of divergent phenotypic responses in populations and genotypes now experiencing climate change.

Local adaptation contributes to gene expression divergence in maize

Gene expression links genotypes to phenotypes, so identifying genes whose expression is shaped by selection will be important for understanding the traits and processes underlying local adaptation. However, detecting local adaptation for gene expression will require distinguishing between divergence due to selection and divergence due to genetic drift. Here, we adapt a QST−FST framework to detect local adaptation for transcriptome-wide gene expression levels in a population of diverse maize genotypes. We compare the number and types of selected genes across a wide range of maize populations and tissues, as well as selection on cold-response genes, drought-response genes, and coexpression clusters. We identify a number of genes whose expression levels are consistent with local adaptation and show that genes involved in stress response show enrichment for selection. Due to its history of intense selective breeding and domestication, maize evolution has long been of interest to researchers, and our study provides insight into the genes and processes important for in local adaptation of maize.

Imprints of Past Habitat Area Reduction on Extant Taxonomic, Functional, and Phylogenetic Composition

Past environmental changes have shaped the evolutionary and ecological diversity of extant organisms. Specifically, climatic fluctuations have made environmental conditions alternatively common or rare over time. Accordingly, most taxa have undergone restriction of their distribution to local refugia during habitat contraction, from which they could expand when suitable habitat became more common. Assessing how past restrictions in refugia have shaped species distributions and genetic diversity has motivated much research in evolutionary biology and biogeography. But there is still lack of clear synthesis on whether and how the taxonomic, functional and phylogenetic composition of extant multispecies assemblages retains the imprint of past restriction in refugia. We devised an original eco-evolutionary model to investigate the temporal dynamics of a regional species pool inhabiting a given habitat today, and which have experienced habitat reduction in the past. The model includes three components: (i) a demographic component driving stochastic changes in population sizes and extinctions due to habitat availability, (ii) a mutation and speciation component representing how divergent genotypes emerge and define new species over time, and (iii) a trait evolution component representing how trait values have changed across descendants over time. We used this model to simulate dynamics of multispecies assemblages that occupied a restricted refugia in the past and could expand their distribution subsequently. We characterized the past restriction in refugia in terms of two parameters representing the ending time of past refugia, and the extent of habitat restriction in the refugia. We characterized extant patterns of taxonomic, functional and phylogenetic diversity depending on these parameters. We found that extant relative abundances reflect the lasting influence of more recent refugia on demographic dynamics, while phylogenetic composition reflects the influence of more ancient habitat change. Extant functional diversity depends on the interplay between diversification dynamics and trait evolution, offering new options to jointly infer current trait adaptation and past trait evolution dynamics.

Correlational selection in the age of genomics

Ecologists and evolutionary biologists are well aware that natural and sexual selection do not operate on traits in isolation, but instead act on combinations of traits. This long-recognized and pervasive phenomenon is known as multivariate selection, or-in the particular case where it favours correlations between interacting traits-correlational selection. Despite broad acknowledgement of correlational selection, the relevant theory has often been overlooked in genomic research. Here, we discuss theory and empirical findings from ecological, quantitative genetic and genomic research, linking key insights from different fields. Correlational selection can operate on both discrete trait combinations and quantitative characters, with profound implications for genomic architecture, linkage, pleiotropy, evolvability, modularity, phenotypic integration and phenotypic plasticity. We synthesize current knowledge and discuss promising research approaches that will enable us to understand how correlational selection shapes genomic architecture, thereby linking quantitative genetic approaches with emerging genomic methods. We suggest that research on correlational selection has great potential to integrate multiple fields in evolutionary biology, including developmental and functional biology, ecology, quantitative genetics, phenotypic polymorphisms, hybrid zones and speciation processes.

Characterizing the effect of background selection on the polygenicity of brain-related traits

BACKGROUND

Genome-wide association studies (GWAS) have demonstrated that psychopathology phenotypes are affected by many risk alleles with small effect (polygenicity). It is unclear how ubiquitously evolutionary pressures influence the genetic architecture of these traits.

METHODS

We partitioned SNP heritability to assess the contribution of background (BGS) and positive selection, Neanderthal local ancestry, functional significance, and genotype networks in 75 brain-related traits (8411 ≤ N ≤ 1,131,181, mean N = 205,289). We applied binary annotations by dichotomizing each measure based on top 2%, 1%, and 0.5% of all scores genome-wide. Effect size distribution features were calculated using GENESIS. We tested the relationship between effect size distribution descriptive statistics and natural selection. In a subset of traits, we explore the inclusion of diagnostic heterogeneity (e.g., number of diagnostic combinations and total symptoms) in the tested relationship.

RESULTS

SNP-heritability was enriched (false discovery rate q < 0.05) for loci with elevated BGS (7 phenotypes) and in genic (34 phenotypes) and loss-of-function (LoF)-intolerant regions (67 phenotypes). These effects were strongest in GWAS of schizophrenia (1.90-fold BGS, 1.16-fold genic, and 1.92-fold LoF), educational attainment (1.86-fold BGS, 1.12-fold genic, and 1.79-fold LoF), and cognitive performance (2.29-fold BGS, 1.12-fold genic, and 1.79-fold LoF). BGS (top 2%) significantly predicted effect size variance for trait-associated loci (σ2 parameter) in 75 brain-related traits (β = 4.39 × 10-5, p = 1.43 × 10-5, model r2 = 0.548). Considering the number of DSM-5 diagnostic combinations per psychiatric disorder improved model fit (σ2 ~ BTop2% × Genic × diagnostic combinations; model r2 = 0.661).

CONCLUSIONS

Brain-related phenotypes with larger variance in risk locus effect sizes are associated with loci under BGS. We show exploratory results suggesting that diagnostic complexity may also contribute to the increased polygenicity of psychiatric disorders.

Two centuries of monarch butterfly collections reveal contrasting effects of range expansion and migration loss on wing traits

Migratory animals exhibit traits that allow them to exploit seasonally variable habitats. In environments where migration is no longer beneficial, such as oceanic islands, migration-association traits may be selected against or be under relaxed selection. Monarch butterflies are best known for their continent-scale migration in North America but have repeatedly become established as nonmigrants in the tropical Americas and on Atlantic and Pacific Islands. These replicated nonmigratory populations provide natural laboratories for understanding the rate of evolution of migration-associated traits. We measured >6,000 museum specimens of monarch butterflies collected from 1856 to the present as well as contemporary wild-caught monarchs from around the world. We determined 1) how wing morphology varies across the monarch's global range, 2) whether initial long-distance founders were particularly suited for migration, and 3) whether recently established nonmigrants show evidence for contemporary phenotypic evolution. We further reared >1,000 monarchs from six populations around the world under controlled conditions and measured migration-associated traits. Historical specimens show that 1) initial founders are well suited for long-distance movement and 2) loss of seasonal migration is associated with reductions in forewing size and elongation. Monarch butterflies raised in a common garden from four derived nonmigratory populations exhibit genetically based reductions in forewing size, consistent with a previous study. Our findings provide a compelling example of how migration-associated traits may be favored during the early stages of range expansion, and also the rate of reductions in those same traits upon loss of migration.

Adaptation to drought is coupled with slow growth, but independent from phenology in marginal silver fir (Abies alba Mill.) populations

Drought is one of the most important selection pressures for forest trees in the context of climate change. Yet, the different evolutionary mechanisms, and their environmental drivers, by which certain populations become more drought tolerant than others is still little understood. We studied adaptation to drought in 16 silver fir (Abies alba Mill.) populations from the French Mediterranean Alps by combining observations on seedlings from a greenhouse experiment (N = 8,199) and on adult tress in situ (N = 315). In the greenhouse, we followed half-sib families for four growing seasons for growth and phenology traits, and tested their water stress response in a "drought until death" experiment. Adult trees in the field were assessed for δ 13C, a proxy for water use efficiency, and genotyped at 357 SNP loci. SNP data was used to generate a null expectation for seedling trait divergence between populations in order to detect the signature of selection, and 31 environmental variables were used to identify the selective environment. We found that seedlings originating from populations with low soil water capacity grew more slowly, attained a smaller stature, and resisted water stress for a longer period of time in the greenhouse. Additionally, adult trees of these populations exhibited a higher water use efficiency as evidenced by their δ 13C. These results suggest a correlated evolution of the growth-drought tolerance trait complex. Population divergence in bud break phenology was adaptive only in the second growing season, and evolved independently from the growth-drought tolerance trait complex. Adaptive divergence in bud break phenology was principally driven by the inter- and intra-annual variation in temperature at the geographic origin of the population. Our results illustrate the different evolutionary strategies used by populations to cope with drought stress at the range limits across a highly heterogeneous landscape, and can be used to inform assisted migration programs.

Detecting selection with a genetic cross

Distinguishing which traits have evolved under natural selection, as opposed to neutral evolution, is a major goal of evolutionary biology. Several tests have been proposed to accomplish this, but these either rely on false assumptions or suffer from low power. Here, I introduce an approach to detecting selection that makes minimal assumptions and only requires phenotypic data from ∼10 individuals. The test compares the phenotypic difference between two populations to what would be expected by chance under neutral evolution, which can be estimated from the phenotypic distribution of an F2 cross between those populations. Simulations show that the test is robust to variation in the number of loci affecting the trait, the distribution of locus effect sizes, heritability, dominance, and epistasis. Comparing its performance to the QTL sign test-an existing test of selection that requires both genotype and phenotype data-the new test achieves comparable power with 50- to 100-fold fewer individuals (and no genotype data). Applying the test to empirical data spanning over a century shows strong directional selection in many crops, as well as on naturally selected traits such as head shape in Hawaiian Drosophila and skin color in humans. Applied to gene expression data, the test reveals that the strength of stabilizing selection acting on mRNA levels in a species is strongly associated with that species' effective population size. In sum, this test is applicable to phenotypic data from almost any genetic cross, allowing selection to be detected more easily and powerfully than previously possible.

The roles of climate, geography and natural selection as drivers of genetic and phenotypic differentiation in a widespread amphibian Hyla annectans (Anura: Hylidae)

The role of geological events and Pleistocene climatic fluctuations as drivers of current patterns of genetic variation in extant species has been a topic of continued interest among evolutionary biologists. Nevertheless, comprehensive studies of widely distributed species are still rare, especially from Asia. Using geographically extensive sampling of many individuals and a large number of nuclear single nucleotide polymorphisms (SNPs), we studied the phylogeography and historical demography of Hyla annectans populations in southern China. Thirty-five sampled populations were grouped into seven clearly defined genetic clusters that closely match phenotype-based subspecies classification. These lineages diverged 2.32-5.23 million years ago (Ma), a timing that closely aligns with the rapid and drastic uplifting of the Qinghai-Tibet Plateau and adjacent southwest China. Demographic analyses and species distribution models indicate that different populations of this species have responded differently to past climatic changes. In the Hengduan Mountains, most populations experienced a bottleneck, whereas the populations located outside of the Hengduan Mountains have gradually declined in size since the end of the last glaciation. In addition, the levels of phenotypic and genetic divergence were strongly correlated across major clades. These results highlight the combined effects of geological events and past climatic fluctuations, as well as natural selection, as drivers of contemporary patterns of genetic and phenotypic variation in a widely distributed anuran in Asia.

Potential adaptive divergence between subspecies and populations of snapdragon plants inferred from QST -FST comparisons

Phenotypic divergence among natural populations can be explained by natural selection or by neutral processes such as drift. Many examples in the literature compare putatively neutral (F_ST ) and quantitative genetic (Q_ST ) differentiation in multiple populations to assess their evolutionary signature and identify candidate traits involved with local adaptation. Investigating these signatures in closely related or recently diversified species has the potential to shed light on the divergence processes acting at the interspecific level. Here, we conducted this comparison in two subspecies of snapdragon plants (eight populations of Antirrhinum majus pseudomajus and five populations of A. m. striatum) in a common garden experiment. We also tested whether altitude was involved with population phenotypic divergence. Our results identified candidate phenological and morphological traits involved with local adaptation. Most of these traits were identified in one subspecies but not the other. Phenotypic divergence increased with altitude for a few biomass-related traits, but only in A. m. striatum. These traits therefore potentially reflect A. m. striatum adaptation to altitude. Our findings imply that adaptive processes potentially differ at the scale of A. majus subspecies.

Microgeographic adaptation and the effect of pollen flow on the adaptive potential of a temperate tree species

In species with long-distance dispersal capacities and inhabiting a large ecological niche, local selection and gene flow are expected to be major evolutionary forces affecting the genetic adaptation of natural populations. Yet, in species such as trees, evidence of microgeographic adaptation and the quantitative assessment of the impact of gene flow on adaptive genetic variation are still limited. Here, we used extensive genetic and phenotypic data from European beech seedlings collected along an elevation gradient, and grown in a common garden, to study the signature of selection on the divergence of eleven potentially adaptive traits, and to assess the role of gene flow in resupplying adaptive genetic variation. We found a significant signal of adaptive differentiation among plots separated by < 1 km, with selection acting on growth and phenological traits. Consistent with theoretical expectations, our results suggest that pollen dispersal contributes to increase genetic diversity for these locally differentiated traits. Our results thus highlight that local selection is an important evolutionary force in natural tree populations and suggest that management interventions to facilitate movement of gametes along short ecological gradients would boost genetic diversity of individual tree populations, and enhance their adaptive potential to rapidly changing environments.

Limited thermal plasticity and geographical divergence in the ovipositor of Drosophila suzukii

Phenotypic plasticity has been repeatedly suggested to facilitate adaptation to new environmental conditions, as in invasions. Here, we investigate this possibility by focusing on the worldwide invasion of Drosophila suzukii: an invasive species that has rapidly colonized all continents over the last decade. This species is characterized by a highly developed ovipositor, allowing females to lay eggs through the skin of ripe fruits. Using a novel approach based on the combined use of scanning electron microscopy and photogrammetry, we quantified the ovipositor size and three-dimensional shape, contrasting invasive and native populations raised at three different developmental temperatures. We found a small but significant effect of temperature and geographical origin on the ovipositor shape, showing the occurrence of both geographical differentiation and plasticity to temperature. The shape reaction norms are in turn strikingly similar among populations, suggesting very little difference in shape plasticity among invasive and native populations, and therefore rejecting the hypothesis of a particular role for the plasticity of the ovipositor in the invasion success. Overall, the ovipositor shape seems to be a fairly robust trait, indicative of stabilizing selection. The large performance spectrum rather than the flexibility of the ovipositor would thus contribute to the success of D. suzukii worldwide invasion.

Common gardens in teosintes reveal the establishment of a syndrome of adaptation to altitude

In plants, local adaptation across species range is frequent. Yet, much has to be discovered on its environmental drivers, the underlying functional traits and their molecular determinants. Genome scans are popular to uncover outlier loci potentially involved in the genetic architecture of local adaptation, however links between outliers and phenotypic variation are rarely addressed. Here we focused on adaptation of teosinte populations along two elevation gradients in Mexico that display continuous environmental changes at a short geographical scale. We used two common gardens, and phenotyped 18 traits in 1664 plants from 11 populations of annual teosintes. In parallel, we genotyped these plants for 38 microsatellite markers as well as for 171 outlier single nucleotide polymorphisms (SNPs) that displayed excess of allele differentiation between pairs of lowland and highland populations and/or correlation with environmental variables. Our results revealed that phenotypic differentiation at 10 out of the 18 traits was driven by local selection. Trait covariation along the elevation gradient indicated that adaptation to altitude results from the assembly of multiple co-adapted traits into a complex syndrome: as elevation increases, plants flower earlier, produce less tillers, display lower stomata density and carry larger, longer and heavier grains. The proportion of outlier SNPs associating with phenotypic variation, however, largely depended on whether we considered a neutral structure with 5 genetic groups (73.7%) or 11 populations (13.5%), indicating that population stratification greatly affected our results. Finally, chromosomal inversions were enriched for both SNPs whose allele frequencies shifted along elevation as well as phenotypically-associated SNPs. Altogether, our results are consistent with the establishment of an altitudinal syndrome promoted by local selective forces in teosinte populations in spite of detectable gene flow. Because elevation mimics climate change through space, SNPs that we found underlying phenotypic variation at adaptive traits may be relevant for future maize breeding.

Effects of marker type and filtering criteria on Q ST-F ST comparisons

Comparative studies of quantitative and neutral genetic differentiation (Q ST-F ST tests) provide means to detect adaptive population differentiation. However, Q ST-F ST tests can be overly liberal if the markers used deflate F ST below its expectation, or overly conservative if methodological biases lead to inflated F ST estimates. We investigated how marker type and filtering criteria for marker selection influence Q ST-F ST comparisons through their effects on F ST using simulations and empirical data on over 18 000 in silico genotyped microsatellites and 3.8 million single-locus polymorphism (SNP) loci from four populations of nine-spined sticklebacks (Pungitius pungitius). Empirical and simulated data revealed that F ST decreased with increasing marker variability, and was generally higher with SNPs than with microsatellites. The estimated baseline F ST levels were also sensitive to filtering criteria for SNPs: both minor alleles and linkage disequilibrium (LD) pruning influenced F ST estimation, as did marker ascertainment. However, in the case of stickleback data used here where Q ST is high, the choice of marker type, their genomic location, ascertainment and filtering made little difference to outcomes of Q ST-F ST tests. Nevertheless, we recommend that Q ST-F ST tests using microsatellites should discard the most variable loci, and those using SNPs should pay attention to marker ascertainment and properly account for LD before filtering SNPs. This may be especially important when level of quantitative trait differentiation is low and levels of neutral differentiation high.

Parasite-mediated selection in a natural metapopulation of Daphnia magna

Parasite-mediated selection varying across time and space in metapopulations is expected to result in host local adaptation and the maintenance of genetic diversity in disease-related traits. However, nonadaptive processes like migration and extinction-(re)colonization dynamics might interfere with adaptive evolution. Understanding how adaptive and nonadaptive processes interact to shape genetic variability in life-history and disease-related traits can provide important insights into their evolution in subdivided populations. Here we investigate signatures of spatially fluctuating, parasite-mediated selection in a natural metapopulation of Daphnia magna. Host genotypes from infected and uninfected populations were genotyped at microsatellite markers, and phenotyped for life-history and disease traits in common garden experiments. Combining phenotypic and genotypic data a Q_ST -F_ST -like analysis was conducted to test for signatures of parasite mediated selection. We observed high variation within and among populations for phenotypic traits, but neither an indication of host local adaptation nor a cost of resistance. Infected populations have a higher gene diversity (Hs) than uninfected populations and Hs is strongly positively correlated with fitness. These results suggest a strong parasite effect on reducing population level inbreeding. We discuss how stochastic processes related to frequent extinction-(re)colonization dynamics as well as host and parasite migration impede the evolution of resistance in the infected populations. We suggest that the genetic and phenotypic patterns of variation are a product of dynamic changes in the host gene pool caused by the interaction of colonization bottlenecks, inbreeding, immigration, hybrid vigor, rare host genotype advantage and parasitism. Our study highlights the effect of the parasite in ameliorating the negative fitness consequences caused by the high drift load in this metapopulation.

Thinking About the Evolution of Complex Traits in the Era of Genome-Wide Association Studies

Many traits of interest are highly heritable and genetically complex, meaning that much of the variation they exhibit arises from differences at numerous loci in the genome. Complex traits and their evolution have been studied for more than a century, but only in the last decade have genome-wide association studies (GWASs) in humans begun to reveal their genetic basis. Here, we bring these threads of research together to ask how findings from GWASs can further our understanding of the processes that give rise to heritable variation in complex traits and of the genetic basis of complex trait evolution in response to changing selection pressures (i.e., of polygenic adaptation). Conversely, we ask how evolutionary thinking helps us to interpret findings from GWASs and informs related efforts of practical importance.

Population genomics perspectives on convergent adaptation

Convergent adaptation is the independent evolution of similar traits conferring a fitness advantage in two or more lineages. Cases of convergent adaptation inform our ideas about the ecological and molecular basis of adaptation. In judging the degree to which putative cases of convergent adaptation provide an independent replication of the process of adaptation, it is necessary to establish the degree to which the evolutionary change is unexpected under null models and to show that selection has repeatedly, independently driven these changes. Here, we discuss the issues that arise from these questions particularly for closely related populations, where gene flow and standing variation add additional layers of complexity. We outline a conceptual framework to guide intuition as to the extent to which evolutionary change represents the independent gain of information owing to selection and show that this is a measure of how surprised we should be by convergence. Additionally, we summarize the ways population and quantitative genetics and genomics may help us address questions related to convergent adaptation, as well as open new questions and avenues of research. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.

Adaptation to local climate in multi-trait space: evidence from silver fir (Abies alba Mill.) populations across a heterogeneous environment

Heterogeneous environments, such as mountainous landscapes, create spatially varying selection pressure that potentially affects several traits simultaneously across different life stages, yet little is known about the general patterns and drivers of adaptation in such complex settings. We studied silver fir (Abies alba Mill.) populations across Switzerland and characterized its mountainous landscape using downscaled historical climate data. We sampled 387 trees from 19 populations and genotyped them at 374 single-nucleotide polymorphisms (SNPs) to estimate their demographic distances. Seedling morphology, growth and phenology traits were recorded in a common garden, and a proxy for water use efficiency was estimated for adult trees. We tested whether populations have more strongly diverged at quantitative traits than expected based on genetic drift alone in a multi-trait framework, and identified potential environmental drivers of selection. We found two main responses to selection: (i) populations from warmer and more thermally stable locations have evolved towards a taller stature, and (ii) the growth timing of populations evolved towards two extreme strategies, 'start early and grow slowly' or 'start late and grow fast', driven by precipitation seasonality. Populations following the 'start early and grow slowly' strategy had higher water use efficiency and came from inner Alpine valleys characterized by pronounced summer droughts. Our results suggest that contrasting adaptive life-history strategies exist in silver fir across different life stages (seedling to adult), and that some of the characterized populations may provide suitable seed sources for tree growth under future climatic conditions.

The Effects of Demography and Genetics on the Neutral Distribution of Quantitative Traits

Neutral models for quantitative trait evolution are useful for identifying phenotypes under selection. These models often assume normally distributed phenotypes. This assumption may be violated when a trait is affected by relatively few variants or when the effects of those variants arise from skewed or heavy tailed distributions. Molecular phenotypes such as gene expression levels may have these properties. To accommodate deviations from normality, models making fewer assumptions about the underlying genetics and patterns of variation are needed. Here, we develop a general neutral model for quantitative trait variation using a coalescent approach. This model allows interpretation of trait distributions in terms of familiar population genetic parameters because it is based on the coalescent. We show how the normal distribution resulting from the infinitesimal limit, where the number of loci grows large as the effect size per mutation becomes small, depends only on expected pairwise coalescent times. We then demonstrate how deviations from normality depend on demography through the distribution of coalescence times as well as through genetic parameters. In particular, population growth events exacerbate deviations while bottlenecks reduce them. We demonstrate the practical applications of this model by showing how to sample from the neutral distribution of [Formula: see text], the ratio of the variance between subpopulations to that in the overall population. We further show it is likely impossible to distinguish sparsity from skewed or heavy tailed mutational effects using only sampled trait values. The model analyzed here greatly expands the parameter space for neutral trait models.

Detecting Adaptive Differentiation in Structured Populations with Genomic Data and Common Gardens

Adaptation in quantitative traits often occurs through subtle shifts in allele frequencies at many loci-a process called polygenic adaptation. While a number of methods have been developed to detect polygenic adaptation in human populations, we lack clear strategies for doing so in many other systems. In particular, there is an opportunity to develop new methods that leverage datasets with genomic data and common garden trait measurements to systematically detect the quantitative traits important for adaptation. Here, we develop methods that do just this, using principal components of the relatedness matrix to detect excess divergence consistent with polygenic adaptation, and using a conditional test to control for confounding effects due to population structure. We apply these methods to inbred maize lines from the United States Department of Agriculture germplasm pool and maize landraces from Europe. Ultimately, these methods can be applied to additional domesticated and wild species to give us a broader picture of the specific traits that contribute to adaptation and the overall importance of polygenic adaptation in shaping quantitative trait variation.

Contrasting patterns of divergence at the regulatory and sequence level in European Daphnia galeata natural populations

Understanding the genetic basis of local adaptation has long been a focus of evolutionary biology. Recently, there has been increased interest in deciphering the evolutionary role of Daphnia's plasticity and the molecular mechanisms of local adaptation. Using transcriptome data, we assessed the differences in gene expression profiles and sequences in four European Daphnia galeata populations. In total, ~33% of 32,903 transcripts were differentially expressed between populations. Among 10,280 differentially expressed transcripts, 5,209 transcripts deviated from neutral expectations and their population-specific expression pattern is likely the result of local adaptation processes. Furthermore, a SNP analysis allowed inferring population structure and distribution of genetic variation. The population divergence at the sequence level was comparatively higher than the gene expression level by several orders of magnitude consistent with strong founder effects and lack of gene flow between populations. Using sequence homology, the candidate transcripts were annotated using a comparative genomics approach. Additionally, we also performed a weighted gene co-expression analysis to identify population-specific regulatory patterns of transcripts in D. galeata. Thus, we identified candidate transcriptomic regions for local adaptation in this key species of aquatic ecosystems in the absence of any laboratory-induced stressor.

A multispecies coalescent model for quantitative traits

We present a multispecies coalescent model for quantitative traits that allows for evolutionary inferences at micro- and macroevolutionary scales. A major advantage of this model is its ability to incorporate genealogical discordance underlying a quantitative trait. We show that discordance causes a decrease in the expected trait covariance between more closely related species relative to more distantly related species. If unaccounted for, this outcome can lead to an overestimation of a trait's evolutionary rate, to a decrease in its phylogenetic signal, and to errors when examining shifts in mean trait values. The number of loci controlling a quantitative trait appears to be irrelevant to all trends reported, and discordance also affected discrete, threshold traits. Our model and analyses point to the conditions under which different methods should fare better or worse, in addition to indicating current and future approaches that can mitigate the effects of discordance.

Determining the evolutionary forces shaping G × E

Contents Summary 31 I. Introduction 31 II. The maintenance of genetic variation for plasticity 32 III. Why is there environmental variation for genetic effects? 33 IV. Conclusions 35 Acknowledgements 35 References 35 SUMMARY: Phenotypic plasticity is common in nature, yet we lack a comprehensive understanding of the evolutionary forces that shape genetic variation for plasticity. This endeavor is especially important because variation for plasticity will result in genotype-by-environment interactions (G × E), a crucial component of variation in quantitative traits. Here, I review our understanding of the evolutionary forces shaping G × E, focusing specifically on: what evolutionary forces maintain variation for plasticity; and what forces maintain different genetic architectures across environments. My specific goal is to show that genomic data can be leveraged to explain the maintenance of G × E by contrasting patterns of genetic variation for plasticity with neutral expectations.

Geographic clines in wing morphology relate to colonization history in New World but not Old World populations of yellow dung flies

Geographic clines offer insights about putative targets and agents of natural selection as well as tempo and mode of adaptation. However, demographic processes can lead to clines that are indistinguishable from adaptive divergence. Using the widespread yellow dung fly Scathophaga stercoraria (Diptera: Scathophagidae), we examine quantitative genetic differentiation (QST ) of wing shape across North America, Europe, and Japan, and compare this differentiation with that of ten microsatellites (FST ). Morphometric analyses of 28 populations reared at three temperatures revealed significant thermal plasticity, sexual dimorphism, and geographic differentiation in wing shape. In North America morphological differentiation followed the decline in microsatellite variability along the presumed route of recent colonization from the southeast to the northwest. Across Europe, where S. stercoraria presumably existed for much longer time and where no molecular pattern of isolation by distance was evident, clinal variation was less pronounced despite significant morphological differentiation (QST >FST ). Shape vector comparisons further indicate that thermal plasticity (hot-to-cold) does not mirror patterns of latitudinal divergence (south-to-north), as might have been expected under a scenario with temperature as the major agent of selection. Our findings illustrate the importance of detailed phylogeographic information when interpreting geographic clines of dispersal traits in an adaptive evolutionary framework.

Detecting Polygenic Adaptation in Admixture Graphs

An open question in human evolution is the importance of polygenic adaptation: adaptive changes in the mean of a multifactorial trait due to shifts in allele frequencies across many loci. In recent years, several methods have been developed to detect polygenic adaptation using loci identified in genome-wide association studies (GWAS). Though powerful, these methods suffer from limited interpretability: they can detect which sets of populations have evidence for polygenic adaptation, but are unable to reveal where in the history of multiple populations these processes occurred. To address this, we created a method to detect polygenic adaptation in an admixture graph, which is a representation of the historical divergences and admixture events relating different populations through time. We developed a Markov chain Monte Carlo (MCMC) algorithm to infer branch-specific parameters reflecting the strength of selection in each branch of a graph. Additionally, we developed a set of summary statistics that are fast to compute and can indicate which branches are most likely to have experienced polygenic adaptation. We show via simulations that this method-which we call PolyGraph-has good power to detect polygenic adaptation, and applied it to human population genomic data from around the world. We also provide evidence that variants associated with several traits, including height, educational attainment, and self-reported unibrow, have been influenced by polygenic adaptation in different populations during human evolution.

Diversifying selection drives parallel evolution of gill raker number and body size along the speciation continuum of European whitefish

Adaptive radiation is the evolution of ecological and phenotypical diversity. It arises via ecological opportunity that promotes the exploration of underutilized or novel niches mediating specialization and reproductive isolation. The assumed precondition for rapid local adaptation is diversifying natural selection, but random genetic drift could also be a major driver of this process. We used 27 populations of European whitefish (Coregonus lavaretus) from nine lakes distributed in three neighboring subarctic watercourses in northern Fennoscandia as a model to test the importance of random drift versus diversifying natural selection for parallel evolution of adaptive phenotypic traits. We contrasted variation for two key adaptive phenotypic traits correlated with resource utilization of polymorphic fish; the number of gill rakers and the total length of fish, with the posterior distribution of neutral genetic differentiation from 13 microsatellite loci, to test whether the observed phenotypic divergence could be achieved by random genetic drift alone. Our results show that both traits have been under diversifying selection and that the evolution of these morphs has been driven by isolation through habitat adaptations. We conclude that diversifying selection acting on gill raker number and body size has played a significant role in the ongoing adaptive radiation of European whitefish morphs in this region.

Selection on the morphology-physiology-performance nexus: Lessons from freshwater stickleback morphs

Conspecifics inhabiting divergent environments frequently differ in morphology, physiology, and performance, but the interrelationships amongst traits and with Darwinian fitness remains poorly understood. We investigated population differentiation in morphology, metabolic rate, and swimming performance in three-spined sticklebacks (Gasterosteus aculeatus L.), contrasting a marine/ancestral population with two distinct freshwater morphotypes derived from it: the "typical" low-plated morph, and a unique "small-plated" morph. We test the hypothesis that similar to plate loss in other freshwater populations, reduction in lateral plate size also evolved in response to selection. Additionally, we test how morphology, physiology, and performance have evolved in concert as a response to differences in selection between marine and freshwater environments. We raised pure-bred second-generation fish originating from three populations and quantified their lateral plate coverage, burst- and critical swimming speeds, as well as standard and active metabolic rates. Using a multivariate QST-FST framework, we detected signals of directional selection on metabolic physiology and lateral plate coverage, notably demonstrating that selection is responsible for the reduction in lateral plate coverage in a small-plated stickleback population. We also uncovered signals of multivariate selection amongst all bivariate trait combinations except the two metrics of swimming performance. Divergence between the freshwater and marine populations exceeded neutral expectation in morphology and in most physiological and performance traits, indicating that adaptation to freshwater habitats has occurred, but through different combinations of traits in different populations. These results highlight both the complex interplay between morphology, physiology and performance in local adaptation, and a framework for their investigation.

Evolutionary dynamics of the leaf phenological cycle in an oak metapopulation along an elevation gradient

It has been predicted that environmental changes will radically alter the selective pressures on phenological traits. Long-lived species, such as trees, will be particularly affected, as they may need to undergo major adaptive change over only one or a few generations. The traits describing the annual life cycle of trees are generally highly evolvable, but nothing is known about the strength of their genetic correlations. Tight correlations can impose strong evolutionary constraints, potentially hampering the adaptation of multivariate phenological phenotypes. In this study, we investigated the evolutionary, genetic and environmental components of the timing of leaf unfolding and senescence within an oak metapopulation along an elevation gradient. Population divergence, estimated from in situ and common-garden data, was compared to expectations under neutral evolution, based on microsatellite markers. This approach made it possible (1) to evaluate the influence of genetic correlation on multivariate local adaptation to elevation and (2) to identify traits probably exposed to past selective pressures due to the colder climate at high elevation. The genetic correlation was positive but very weak, indicating that genetic constraints did not shape the local adaptation pattern for leaf phenology. Both spring and fall (leaf unfolding and senescence, respectively) phenology timings were involved in local adaptation, but leaf unfolding was probably the trait most exposed to climate change-induced selection. Our data indicated that genetic variation makes a much smaller contribution to adaptation than the considerable plastic variation displayed by a tree during its lifetime. The evolutionary potential of leaf phenology is, therefore, probably not the most critical aspect for short-term population survival in a changing climate.

Population genetic differentiation of height and body mass index across Europe

Across-nation differences in the mean values for complex traits are common, but the reasons for these differences are unknown. Here we find that many independent loci contribute to population genetic differences in height and body mass index (BMI) in 9,416 individuals across 14 European countries. Using discovery data on over 250,000 individuals and unbiased effect size estimates from 17,500 sibling pairs, we estimate that 24% (95% credible interval (CI) = 9%, 41%) and 8% (95% CI = 4%, 16%) of the captured additive genetic variance for height and BMI, respectively, reflect population genetic differences. Population genetic divergence differed significantly from that in a null model (height, P < 3.94 × 10(-8); BMI, P < 5.95 × 10(-4)), and we find an among-population genetic correlation for tall and slender individuals (r = -0.80, 95% CI = -0.95, -0.60), consistent with correlated selection for both phenotypes. Observed differences in height among populations reflected the predicted genetic means (r = 0.51; P < 0.001), but environmental differences across Europe masked genetic differentiation for BMI (P < 0.58).

Common garden experiments in the genomic era: new perspectives and opportunities

The study of local adaptation is rendered difficult by many evolutionary confounding phenomena (for example, genetic drift and demographic history). When complex traits are involved in local adaptation, phenomena such as phenotypic plasticity further hamper evolutionary biologists to study the complex relationships between phenotype, genotype and environment. In this perspective paper, we suggest that the common garden experiment, specifically designed to deal with phenotypic plasticity, has a clear role to play in the study of local adaptation, even (if not specifically) in the genomic era. After a quick review of some high-throughput genotyping protocols relevant in the context of a common garden, we explore how to improve common garden analyses with dense marker panel data and recent statistical methods. We then show how combining approaches from population genomics and genome-wide association studies with the settings of a common garden can yield to a very efficient, thorough and integrative study of local adaptation. Especially, evidence from genomic (for example, genome scan) and phenotypic origins constitute independent insights into the possibility of local adaptation scenarios, and genome-wide association studies in the context of a common garden experiment allow to decipher the genetic bases of adaptive traits.