The evolution of phenotypic plasticity when environments fluctuate in time and space

Genome-wide mapping of regulatory variants for temperature- and salinity-adaptive genes reveals genetic basis of genotype-by-environment interaction in Crassostrea ariakensis

Regulatory variants in gene expression serve as bridges linking genetic variation and phenotypic plasticity. Environmental conditions typically influence the effects of regulatory variants on phenotypic plasticity; however, such genotype-by-environment interactions (G × E) are poorly understood. This study aimed to investigate the genetic basis of G × E in estuarine oyster (Crassostrea ariakensis), which is an important model animal for studying environmental adaption owing to its high plasticity and large intraspecific divergence. Genome-wide mapping of expression quantitative trait loci (eQTLs) for 23 environmental adaptive genes was performed for 256 estuarine oysters. We identified 1194 eQTL single nucleotide polymorphisms (eSNPs), including 433 cis-eSNPs in four genes and 722 trans-eSNPs in eight genes. The expression variation explanation of cis-eSNPs (9.95%) was significantly higher than that of trans-eSNPs (9.15%). We specifically showed cis- and trans-eSNPs with high linkage disequilibrium (LD) for Traf7, Slc6a5, Ggt, and Dap3. For example, we identified a cis-regulatory LD block containing 68 cis-eSNP and a trans-regulatory LD block, including 20 trans-eSNPs in Traf7. A high proportion (85%) of 40 vital eSNPs exhibited significant G × E effects. We identified crossing and nonparallel interactions of G × E, with the tag cis-eSNPs of Baat and Slc6a5 as representatives. Our results indicated that cis-eQTLs are highly conserved. This study provides insights into the understanding of adaptive evolutionary mechanisms and phenotypic response prediction to variable environments, as well as the genetic improvement for superior adaptive traits for genetic resource conservation and aquaculture.

Role of genetic architecture in phenotypic plasticity

Phenotypic plasticity, the ability of an organism to display different phenotypes across environments, is widespread in nature. Plasticity aids survival in novel environments. Herein, we review studies from yeast that allow us to start uncovering the genetic architecture of phenotypic plasticity. Genetic variants and their interactions impact the phenotype in different environments, and distinct environments modulate the impact of genetic variants and their interactions on the phenotype. Because of this, certain hidden genetic variation is expressed in specific genetic and environmental backgrounds. A better understanding of the genetic mechanisms of phenotypic plasticity will help to determine short- and long-term responses to selection and how wide variation in disease manifestation occurs in human populations.

Antagonistic Effects of Assortative Mating on the Evolution of Phenotypic Plasticity along Environmental Gradients

AbstractPrevious theory has shown that assortative mating for plastic traits can maintain genetic divergence across environmental gradients despite high gene flow. Yet these models did not examine how assortative mating affects the evolution of plasticity. We here describe patterns of genetic variation across elevation for plasticity in a trait under assortative mating, using multiple-year observations of budburst date in a common garden of sessile oaks. Despite high gene flow, we found significant spatial genetic divergence for the intercept, but not for the slope, of reaction norms to temperature. We then used individual-based simulations, where both the slope and the intercept of the reaction norm evolve, to examine how assortative mating affects the evolution of plasticity, varying the intensity and distance of gene flow. Our model predicts the evolution of either suboptimal plasticity (reaction norms with a slope shallower than optimal) or hyperplasticity (slopes steeper than optimal) in the presence of assortative mating when optimal plasticity would evolve under random mating. Furthermore, a cogradient pattern of genetic divergence for the intercept of the reaction norm (where plastic and genetic effects are in the same direction) always evolves in simulations with assortative mating, consistent with our observations in the studied oak populations.

Bet-hedging via dispersal aids the evolution of plastic responses to unreliable cues

Adaptive plasticity is expected to evolve when informative cues predict environmental variation. However, plastic responses can be maladaptive even when those cues are informative, if prediction mistakes are shared across members of a generation. These fitness costs can constrain the evolution of plasticity when initial plastic mutants use of cues of only moderate reliability. Here, we model the barriers to the evolution of plasticity produced by these constraints and show that dispersal across a metapopulation can overcome them. Constraints are also lessened, though not eliminated, when plastic responses are free to evolve gradually and in concert with increased reliability. Each of these factors be viewed as a form of bet-hedging: by lessening correlations in the fates of relatives, dispersal acts as diversifying bet-hedging, while producing submaximal responses to a cue can be understood as a conservative bet-hedging strategy. While poor information may constrain the evolution of plasticity, the opportunity for bet-hedging may predict when that constraint can be overcome.

Population Genomic Time Series Data of a Natural Population Suggests Adaptive Tracking of Fluctuating Environmental Changes

Natural populations are constantly exposed to fluctuating environmental changes that negatively affect their fitness in unpredictable ways. While theoretical models show the possibility of counteracting these environmental changes through rapid evolutionary adaptations, there have been few empirical studies demonstrating such adaptive tracking in natural populations. Here, we analyzed environmental data, fitness-related phenotyping and genomic time-series data sampled over 3 years from a natural Chironomus riparius (Diptera, Insecta) population to address this question. We show that the population's environment varied significantly on the time scale of the sampling in many selectively relevant dimensions, independently of each other. Similarly, phenotypic fitness components evolved significantly on the same temporal scale (mean 0.32 Haldanes), likewise independent from each other. The allele frequencies of 367,446 SNPs across the genome showed evidence of positive selection. Using temporal correlation of spatially coherent allele frequency changes revealed 35,574 haplotypes with more than one selected SNP. The mean selection coefficient for these haplotypes was 0.30 (s.d. = 0.68). The frequency changes of these haplotypes clustered in 46 different temporal patterns, indicating concerted, independent evolution of many polygenic traits. Nine of these patterns were strongly correlated with measured environmental variables. Enrichment analysis of affected genes suggested the implication of a wide variety of biological processes. Thus, our results suggest overall that the natural population of C. riparius tracks environmental change through rapid polygenic adaptation in many independent dimensions. This is further evidence that natural selection is pervasive at the genomic level and that evolutionary and ecological time scales may not differ at all, at least in some organisms.

Experimental evolution of environmental tolerance, acclimation, and physiological plasticity in a randomly fluctuating environment

Environmental tolerance curves, representing absolute fitness against the environment, are an empirical assessment of the fundamental niche, and emerge from the phenotypic plasticity of underlying phenotypic traits. Dynamic plastic responses of these traits can lead to acclimation effects, whereby recent past environments impact current fitness. Theory predicts that higher levels of phenotypic plasticity should evolve in environments that fluctuate more predictably, but there have been few experimental tests of these predictions. Specifically, we still lack experimental evidence for the evolution of acclimation effects in response to environmental predictability. Here, we exposed 25 genetically diverse populations of the halotolerant microalgae Dunaliella salina to different constant salinities, or to randomly fluctuating salinities, for over 200 generations. The fluctuating treatments differed in their autocorrelation, which determines the similarity of subsequent values, and thus environmental predictability. We then measured acclimated tolerance surfaces, mapping population growth rate against past (acclimation) and current (assay) environments. We found that experimental mean and variance in salinity caused the evolution of niche position (optimal salinity) and breadth, with respect to not only current but also past (acclimation) salinity. We also detected weak but significant evidence for evolutionary changes in response to environmental predictability, with higher predictability leading notably to lower optimal salinities and stronger acclimation effect of past environment on current fitness. We further showed that these responses are related to the evolution of plasticity for intracellular glycerol, the major osmoregulatory mechanism in this species. However, the direction of plasticity evolution did not match simple theoretical predictions. Our results underline the need for a more explicit consideration of the dynamics of environmental tolerance and its underlying plastic traits to reach a better understanding of ecology and evolution in fluctuating environments.

Adaptive phenotypic plasticity is under stabilizing selection in Daphnia

The adaptive nature of phenotypic plasticity is widely documented. However, little is known about the evolutionary forces that shape genetic variation of plasticity within populations. Whether genetic variation in plasticity is driven by stabilizing or diversifying selection and whether the strength of such forces remains constant through time, remain open questions. Here, we address this issue by assessing the evolutionary forces that shape genetic variation in antipredator developmental plasticity of Daphnia pulex. Antipredator plasticity in D. pulex is characterized by the growth of a pedestal and spikes in the dorsal head region upon exposure to predator cue. We characterized genetic variation in plasticity using a method that describes the entire dorsal shape amongst >100 D. pulex strains recently derived from the wild. We observed the strongest reduction in genetic variation in dorsal areas where plastic responses were greatest, consistent with stabilizing selection. We compared mutational variation (Vm) to standing variation (Vg) and found that Vg/Vm is lowest in areas of greatest plasticity, again consistent with stabilizing selection. Our results suggest that stabilizing selection operates directly on phenotypic plasticity in Daphnia and provide a rare glimpse into the evolution of fitness-related traits in natural populations.

The role of phenotypic plasticity in the establishment of range margins

It has been argued that adaptive phenotypic plasticity may facilitate range expansions over spatially and temporally variable environments. However, plasticity may induce fitness costs. This may hinder the evolution of plasticity. Earlier modelling studies examined the role of plasticity during range expansions of populations with fixed genetic variance. However, genetic variance evolves in natural populations. This may critically alter model outcomes. We ask: how does the capacity for plasticity in populations with evolving genetic variance alter range margins that populations without the capacity for plasticity are expected to attain? We answered this question using computer simulations and analytical approximations. We found a critical plasticity cost above which the capacity for plasticity has no impact on the expected range of the population. Below the critical cost, by contrast, plasticity facilitates range expansion, extending the range in comparison to that expected for populations without plasticity. We further found that populations may evolve plasticity to buffer temporal environmental fluctuations, but only when the plasticity cost is below the critical cost. Thus, the cost of plasticity is a key factor involved in range expansions of populations with the potential to express plastic response in the adaptive trait. This article is part of the theme issue 'Species' ranges in the face of changing environments (part I)'.

Directional selection and the evolution of breeding date in birds, revisited: Hard selection and the evolution of plasticity

The mismatch between when individuals breed and when we think they should breed has been a long‐standing problem in evolutionary ecology. Price et al. is a classic theory paper in this field and is mainly cited for its most obvious result: if individuals with high nutritional condition breed early, then the advantage of breeding early may be overestimated when information on nutritional condition is absent. Price at al.'s less obvious result is that individuals, on average, are expected to breed later than the optimum. Here, we provide an explanation of their non‐intuitive result in terms of hard selection, and go on to show that neither of their results are expected to hold if the relationship between breeding date and nutrition is allowed to evolve. By introducing the assumption that the advantage of breeding early is greater for individuals in high nutritional condition, we show that their most cited result can be salvaged. However, individuals, on average, are expected to breed earlier than the optimum, not later. More generally, we also show that the hard selection mechanisms that underpin these results have major implications for the evolution of plasticity: when environmental heterogeneity becomes too great, plasticity is selected against, prohibiting the evolution of generalists.

Spatial variation in high temperature-regulated gene expression predicts evolution of plasticity with climate change in the scarlet monkeyflower

A major way that organisms can adapt to changing environmental conditions is by evolving increased or decreased phenotypic plasticity. In the face of current global warming, more attention is being paid to the role of plasticity in maintaining fitness as abiotic conditions change over time. However, given that temporal data can be challenging to acquire, a major question is whether evolution in plasticity across space can predict adaptive plasticity across time. In growth chambers simulating two thermal regimes, we generated transcriptome data for western North American scarlet monkeyflowers (Mimulus cardinalis) collected from different latitudes and years (2010 and 2017) to test hypotheses about how plasticity in gene expression is responding to increases in temperature, and if this pattern is consistent across time and space. Supporting the genetic compensation hypothesis, individuals whose progenitors were collected from the warmer-origin northern 2017 descendant cohort showed lower thermal plasticity in gene expression than their cooler-origin northern 2010 ancestors. This was largely due to a change in response at the warmer (40°C) rather than cooler (20°C) treatment. A similar pattern of reduced plasticity, largely due to a change in response at 40°C, was also found for the cooler-origin northern versus the warmer-origin southern population from 2017. Our results demonstrate that reduced phenotypic plasticity can evolve with warming and that spatial and temporal changes in plasticity predict one another.

Counter-gradient variation of reproductive effort in a widely distributed temperate oak (Quercus petraea)

The genetic and phenotypic variability of life history traits determines the demographic attributes of tree populations and, thus, their responses to anthropogenic climate change. Growth- and survival-related traits have been widely studied in forest ecology, but little is known about the determinism of reproductive traits.Using an elevation gradient experiment in the Pyrenees we assessed the degree to which variations in reproductive effort along climatic gradients are environmentally or genetically driven, by comparing oak populations (Quercus petraea) growing under field and common garden conditions.In situ monitoring revealed a decline in reproductive effort with increasing elevation and decreasing temperature. In common garden conditions, significant genetic differentiation was observed between provenances for reproduction and growth: trees from cold environments (high elevations) grew more slowly, and produced larger acorns in larger numbers. Our observations show that genetic and phenotypic clines for reproductive traits have opposite signs (counter-gradient) along the environmental gradient as opposed to growth, for which genetic variation parallels phenotypic variation (co-gradient).The counter-gradient found here for reproductive effort reveals that genetic variation partly counteracts the phenotypic effect of temperature, moderating the change in reproductive effort according to temperature. We consider the possible contribution to this counter-gradient in reproductive effort as an evolutionary trade-off between reproduction and growth.

Individual reversible plasticity as a genotype-level bet-hedging strategy

Reversible plasticity in phenotypic traits allows organisms to cope with environmental variation within lifetimes, but costs of plasticity may limit just how well the phenotype matches the environmental optimum. An additional adaptive advantage of plasticity might be to reduce fitness variance, in other words: bet-hedging to maximize geometric (rather than simply arithmetic) mean fitness. Here, we model the evolution of plasticity in the form of reaction norm slopes, with increasing costs as the slope or degree of plasticity increases. We find that greater investment in plasticity (i.e. a steeper reaction norm slope) is favoured in scenarios promoting bet-hedging as a response to multiplicative fitness accumulation (i.e. coarser environmental grains and fewer time steps prior to reproduction), because plasticity lowers fitness variance across environmental conditions. In contrast, in scenarios with finer environmental grain and many time steps prior to reproduction, bet-hedging plays less of a role and individual-level optimization favours evolution of shallower reaction norm slopes. However, the opposite pattern holds if plasticity costs themselves result in increased fitness variation, as might be the case for production costs of plasticity that depend on how much change is made to the phenotype each time step. We discuss these contrasting predictions from this partitioning of adaptive plasticity into short-term individual benefits versus long-term genotypic (bet-hedging) benefits, and how this approach enhances our understanding of the evolution of optimum levels of plasticity in examples from thermal physiology to advances in avian lay dates.

A population-genomic approach for estimating selection on polygenic traits in heterogeneous environments

Strong selection can cause rapid evolutionary change, but temporal fluctuations in the form, direction and intensity of selection can limit net evolutionary change over longer time periods. Fluctuating selection could affect molecular diversity levels and the evolution of plasticity and ecological specialization. Nonetheless, this phenomenon remains understudied, in part because of analytical limitations and the general difficulty of detecting selection that does not occur in a consistent manner. Herein, I fill this analytical gap by presenting an approximate Bayesian computation (ABC) method to detect and quantify fluctuating selection on polygenic traits from population genomic time-series data. I propose a model for environment-dependent phenotypic selection. The evolutionary genetic consequences of selection are then modelled based on a genotype-phenotype map. Using simulations, I show that the proposed method generates accurate and precise estimates of selection when the generative model for the data is similar to the model assumed by the method. The performance of the method when applied to an evolve-and-resequence study of host adaptation in the cowpea seed beetle (Callosobruchus maculatus) was more idiosyncratic and depended on specific analytical choices. Despite some limitations, these results suggest the proposed method provides a powerful approach to connect the causes of (variable) selection to traits and genome-wide patterns of evolution. Documentation and open-source computer software (fsabc) implementing this method are available from github (https://github.com/zgompert/fsabc.git).

Differences in spatial versus temporal reaction norms for spring and autumn phenological events

To do the right thing at the right time, organisms need to glean cues from their environment. How they respond can then be described by reaction norms, i.e., by the relationship between the phenotype expressed (the phenology of an event) and the environment (the date when a given number of degree-days are achieved). We use information on 178 phenological events across the former Soviet Union. We found the timing of events to differ more between sites in spring and less in autumn. These patterns of local adaptation translate to a massive imprint on nature’s calendar: geographic variation in phenology is more pronounced in spring and less pronounced in autumn than if organisms were to respond equally everywhere.

For species to stay temporally tuned to their environment, they use cues such as the accumulation of degree-days. The relationships between the timing of a phenological event in a population and its environmental cue can be described by a population-level reaction norm. Variation in reaction norms along environmental gradients may either intensify the environmental effects on timing (cogradient variation) or attenuate the effects (countergradient variation). To resolve spatial and seasonal variation in species’ response, we use a unique dataset of 91 taxa and 178 phenological events observed across a network of 472 monitoring sites, spread across the nations of the former Soviet Union. We show that compared to local rates of advancement of phenological events with the advancement of temperature-related cues (i.e., variation within site over years), spatial variation in reaction norms tend to accentuate responses in spring (cogradient variation) and attenuate them in autumn (countergradient variation). As a result, among-population variation in the timing of events is greater in spring and less in autumn than if all populations followed the same reaction norm regardless of location. Despite such signs of local adaptation, overall phenotypic plasticity was not sufficient for phenological events to keep exact pace with their cues—the earlier the year, the more did the timing of the phenological event lag behind the timing of the cue. Overall, these patterns suggest that differences in the spatial versus temporal reaction norms will affect species’ response to climate change in opposite ways in spring and autumn.

Unveiling the roles of interspecific competition and local adaptation in phenotypic differentiation of parapatric frogs

Understanding how ecological processes affect phenotypic evolution has been and continues to be an important goal of ecology and evolutionary biology. Interspecific competition for resources can be a selective force driving phenotypic differentiation that reduces competition among sympatric species (character divergence), enabling closely-related species to coexist. However, although patterns of character divergence are well documented in both empirical and theoretical researches, how local adaptation to abiotic environment affects trait evolution in the face of interspecific competition is less known. Here, we investigate how patterns in morphological traits of 2 parapatric frog species, Feirana quadranus and F. taihangnica, vary among allopatric and sympatric regions using range-wide data derived from extensive field surveys. Feirana quadranus was overall larger than F. taihangnica in body size (i.e., snout–vent length [SVL]), and the difference between SVL of both species in sympatry was larger than that in allopatry. From allopatry to sympatry, the 2 species diverged in foot and hand traits, but converged in eye size and interorbital span, even when we controlled for the effects of geographic gradients. Sympatric divergence in SVL, hand and foot traits is likely acting as a case of evolutionary shift caused by interspecific competition. In contrast, sympatric convergence of eye-related traits may derive at least partly from adaptation to local environments. These results imply the relative roles of interspecific competition and local adaptation in shaping phenotypic diversification. Our findings illustrate how traits evolve in parapatric species pair due to sympatric divergent and convergent evolution. It thus provides insights into understanding underlying evolutionary processes of parapatric species, that is, competition and local adaptation.

Natural and human-induced environmental changes and their effects on adaptive potential of wild animal populations

A major challenge of evolutionary ecology over the next decades is to understand and predict the consequences of the current rapid and important environmental changes on wild populations. Extinction risk of species is linked to populations' evolutionary potential and to their ability to express adaptive phenotypic plasticity. There is thus a vital need to quantify how selective pressures, quantitative genetics parameters, and phenotypic plasticity, for multiple traits in wild animal populations, may vary with changes in the environment. Here I review our previous research that integrated ecological and evolutionary theories with molecular ecology, quantitative genetics, and long-term monitoring of individually marked wild animals. Our results showed that assessing evolutionary and plastic changes over time and space, using multi-trait approaches, under a realistic range of environmental conditions are crucial steps toward improving our understanding of the evolution and adaptation of natural populations. Our current and future work focusses on assessing the limits of adaptive potential by determining the factors constraining the evolvability of plasticity, those generating covariation among genetic variance and selection, as well as indirect genetic effects, which can affect population's capacity to adjust to environmental changes. In doing so, we aim to provide an improved assessment of the spatial and temporal scale of evolutionary processes in wild animal populations.

Natural Selection Shapes Variation in Genome-wide Recombination Rate in Drosophila pseudoobscura

While recombination is widely recognized to be a key modulator of numerous evolutionary phenomena, we have a poor understanding of how recombination rate itself varies and evolves within a species. Here, we performed a comprehensive study of recombination rate (rate of meiotic crossing over) in two natural populations of Drosophila pseudoobscura from Utah and Arizona, USA. We used an amplicon sequencing approach to obtain high-quality genotypes in approximately 8,000 individual backcrossed offspring (17 mapping populations with roughly 530 individuals each), for which we then quantified crossovers. Interestingly, variation in recombination rate within and between populations largely manifested as differences in genome-wide recombination rate rather than remodeling of the local recombination landscape. Comparing populations, we discovered individuals from the Utah population displayed on average 8% higher crossover rates than the Arizona population, a statistically significant difference. Using a Q_ST-F_ST analysis, we found that this difference in crossover rate was dramatically higher than expected under neutrality, indicating that this difference may have been driven by natural selection. Finally, using a combination of short- and long-read whole-genome sequencing, we found no significant association between crossover rate and structural variation at the 200-400 kb scale. Our results demonstrate that (1) there is abundant variation in genome-wide crossover rate in natural populations, (2) at the 200-400 kb scale, recombination rate appears to vary largely genome-wide, rather than in specific intervals, and (3) interpopulation differences in recombination rate may be the result of local adaptation.

Where is the optimum? Predicting the variation of selection along climatic gradients and the adaptive value of plasticity. A case study on tree phenology

Many theoretical models predict when genetic evolution and phenotypic plasticity allow adaptation to changing environmental conditions. These models generally assume stabilizing selection around some optimal phenotype. We however often ignore how optimal phenotypes change with the environment, which limit our understanding of the adaptive value of phenotypic plasticity. Here, we propose an approach based on our knowledge of the causal relationships between climate, adaptive traits, and fitness to further these questions. This approach relies on a sensitivity analysis of the process‐based model phenofit, which mathematically formalizes these causal relationships, to predict fitness landscapes and optimal budburst dates along elevation gradients in three major European tree species. Variation in the overall shape of the fitness landscape and resulting directional selection gradients were found to be mainly driven by temperature variation. The optimal budburst date was delayed with elevation, while the range of dates allowing high fitness narrowed and the maximal fitness at the optimum decreased. We also found that the plasticity of the budburst date should allow tracking the spatial variation in the optimal date, but with variable mismatch depending on the species, ranging from negligible mismatch in fir, moderate in beech, to large in oak. Phenotypic plasticity would therefore be more adaptive in fir and beech than in oak. In all species, we predicted stronger directional selection for earlier budburst date at higher elevation. The weak selection on budburst date in fir should result in the evolution of negligible genetic divergence, while beech and oak would evolve counter‐gradient variation, where genetic and environmental effects are in opposite directions. Our study suggests that theoretical models should consider how whole fitness landscapes change with the environment. The approach introduced here has the potential to be developed for other traits and species to explore how populations will adapt to climate change.