The genetics of adaptation to discrete heterogeneous environments: frequent mutation or large‐effect alleles can allow range expansion

Range expansions are complex evolutionary and ecological processes. From an evolutionary standpoint, a populations' adaptive capacity can determine the success or failure of expansion. Using individual‐based simulations, we model range expansion over a two‐dimensional, approximately continuous landscape. We investigate the ability of populations to adapt across patchy environmental gradients and examine how the effect sizes of mutations influence the ability to adapt to novel environments during range expansion. We find that genetic architecture and landscape patchiness both have the ability to change the outcome of adaptation and expansion over the landscape. Adaptation to new environments succeeds via many mutations of small effect or few of large effect, but not via the intermediate between these cases. Higher genetic variance contributes to increased ability to adapt, but an alternative route of successful adaptation can proceed from low genetic variance scenarios with alleles of sufficiently large effect. Steeper environmental gradients can prevent adaptation and range expansion on both linear and patchy landscapes. When the landscape is partitioned into local patches with sharp changes in phenotypic optimum, the local magnitude of change between subsequent patches in the environment determines the success of adaptation to new patches during expansion.     

Role of phenotypic plasticity and population differentiation in adaptation to novel environmental conditions

Species can adapt to new environmental conditions either through individual phenotypic plasticity, intraspecific genetic differentiation in adaptive traits, or both. Wild emmer wheat, Triticum dicoccoides, an annual grass with major distribution in Eastern Mediterranean region, is predicted to experience in the near future, as a result of global climate change, conditions more arid than in any part of the current species distribution. To understand the role of the above two means of adaptation, and the effect of population range position, we analyzed reaction norms, extent of plasticity, and phenotypic selection across two experimental environments of high and low water availability in two core and two peripheral populations of this species. We studied 12 quantitative traits, but focused primarily on the onset of reproduction and maternal investment, which are traits that are closely related to fitness and presumably involved in local adaptation in the studied species. We hypothesized that the population showing superior performance under novel environmental conditions will either be genetically differentiated in quantitative traits or exhibit higher phenotypic plasticity than the less successful populations. We found the core population K to be the most plastic in all three trait categories (phenology, reproductive traits, and fitness) and most successful among populations studied, in both experimental environments; at the same time, the core K population was clearly genetically differentiated from the two edge populations. Our results suggest that (1) two means of successful adaptation to new environmental conditions, phenotypic plasticity and adaptive genetic differentiation, are not mutually exclusive ways of achieving high adaptive ability; and (2) colonists from some core populations can be more successful in establishing beyond the current species range than colonists from the range extreme periphery with conditions seemingly closest to those in the new environment.     

Influence of learning on range expansion and adaptation to novel habitats

Learning has been postulated to ‘drive’ evolution, but its influence on adaptive evolution in heterogeneous environments has not been formally examined. We used a spatially explicit individual‐based model to study the effect of learning on the expansion and adaptation of a species to a novel habitat. Fitness was mediated by a behavioural trait (resource preference), which in turn was determined by both the genotype and learning. Our findings indicate that learning substantially increases the range of parameters under which the species expands and adapts to the novel habitat, particularly if the two habitats are separated by a sharp ecotone (rather than a gradient). However, for a broad range of parameters, learning reduces the degree of genetically‐based local adaptation following the expansion and facilitates maintenance of genetic variation within local populations. Thus, in heterogeneous environments learning may facilitate evolutionary range expansions and maintenance of the potential of local populations to respond to subsequent environmental changes.     

Phenotypic variation of life‐history traits in native,invasive, and landrace populations of Brassica tournefortii

Varying environments can result in different patterns of adaptive phenotypes. By performing a common greenhouse experiment, we identified phenotypic differentiation on phenology, leaf morphology, branch architecture, size, and reproduction, among native, invasive, and landrace ranges of Brassica tournefortii. We first compared trait means and fitness functions among ranges, then we analyzed how trait means and selection strength of populations respond to varying aridity. Most traits varied such that landrace > invasive > native. Excluding reproduction, which was positively selected, most trait PCs experienced nonlinear selection in the native range but frequently shifted to directional selection in invasive and/or landrace ranges. The absence of strong clines for trait means in landrace and invasive populations suggest that agricultural practices and novel environments in source locations affected adaptive potential. Selection strength on faster reproductive phenology (negative directional) and leaf margin trait (disruptive) PCs coincided with increasing moisture. In native populations, higher aridity was associated with more days to reproduction, but landrace and invasive populations show stable mean time to reproduction with increasing moisture. A stable adaptive trait can increase range expansion in the invasive range, but stability can be beneficial for future harvest of B. tournefortii seed crops in the face of climate change.     

Life history variation in an annual plant under two opposing environmental constraints along an aridity gradient

Environmental gradients represent an ideal framework for studying adaptive variation in the life history of plant species. However, on very steep gradients, largely contrasting conditions at the two gradient ends often limit the distribution of the same species across the whole range of environmental conditions. Here, we study phenotypic variation in a winter annual crucifer Biscutella didyma persisting along a steep gradient of increasing rainfall in Israel. In particular, we explored whether the life history at the arid end of the gradient indicates adaptations to drought and unpredictable conditions, while adaptations to the highly competitive environment prevail at the mesic Mediterranean end. We examined several morphological and reproductive traits in four natural populations and in populations cultivated in standard common environment. Plants from arid environments were faster in phenological development, more branched in architecture and tended to maximize reproduction, while the Mediterranean plants invested mainly in vertical vegetative growth. Differences between cultivation and field in diaspore production were very large for arid populations as opposed to Mediterranean ones, indicating a larger potential to increase reproduction under favorable conditions. Our overall findings indicate two strongly opposing selective forces at the two extremes of the aridity gradient, which result in contrasting strategies within the studied annual plant species.     

Behavioral constraints on local adaptation and counter‐gradient variation: Implications for climate change

Resource allocation to growth, reproduction, and body maintenance varies within species along latitudinal gradients. Two hypotheses explaining this variation are local adaptation and counter‐gradient variation. The local adaptation hypothesis proposes that populations are adapted to local environmental conditions and are therefore less adapted to environmental conditions at other locations. The counter‐gradient variation hypothesis proposes that one population out performs others across an environmental gradient because its source location has greater selective pressure than other locations. Our study had two goals. First, we tested the local adaptation and counter‐gradient variation hypotheses by measuring effects of environmental temperature on phenotypic expression of reproductive traits in the burying beetle, Nicrophorus orbicollis Say, from three populations along a latitudinal gradient in a common garden experimental design. Second, we compared patterns of variation to evaluate whether traits covary or whether local adaptation of traits precludes adaptive responses by others. Across a latitudinal range, N. orbicollis exhibits variation in initiating reproduction and brood sizes. Consistent with local adaptation: (a) beetles were less likely to initiate breeding at extreme temperatures, especially when that temperature represents their source range; (b) once beetles initiate reproduction, source populations produce relatively larger broods at temperatures consistent with their local environment. Consistent with counter‐gradient variation, lower latitude populations were more successful at producing offspring at lower temperatures. We found no evidence for adaptive variation in other adult or offspring performance traits. This suite of traits does not appear to coevolve along the latitudinal gradient. Rather, response to selection to breed within a narrow temperature range may preclude selection on other traits. Our study highlights that N. orbicollis uses temperature as an environmental cue to determine whether to initiate reproduction, providing insight into how behavior is modified to avoid costly reproductive attempts. Furthermore, our results suggest a temperature constraint that shapes reproductive behavior.     

Competition and evolution along environmental gradients: patterns, boundaries and sympatric divergence

The present study examined how competitive interactions and environmental conditions generate species boundaries and determine species distributions. A spatially explicit, quantitative genetic, two-species competition model was used to manipulate the strengths of competition, gene flow and local adaptation along environmental gradients. This allowed us to assess the long-term persistence of each species and whether the ranges they inhabited had boundaries in space or were unlimited. We found that a species boundary arises along less steep environmental gradients when the strength of stabilizing selection and diversifying selection are similar. We also found that a species boundary may arise along shallow environmental gradients if interspecific competition is more intense than intraspecific, which relaxes previous requirements for steep gradients for generating range limits. We determined an analytical form for the critical environmental gradient as a function of ecological and genetic parameters at which a species boundary is expected to arise by competition. Results suggest an alternative to resource competition as an explanation for phenotypic divergence between sympatric competitors. Competitors sharing a trait that is under stabilizing selection along an environmental gradient may segregate spatially and evolve in different regions, with phenotypic sympatric divergence reflecting the resulting clines. Along various types of environmental gradients, variation in stabilizing selection intensities could lead to contrasting patterns in the distribution of species. For stabilizing selection strengths in accord with field data estimates, this study predicts that the level of sympatric character divergence would be limited along environmental gradients.     

Architecture of coastal and desert Encelia farinosa (Asteraceae): consequences of plastic and heritable variation in leaf characters

The shrub Encelia farinosa (Asteraceae) exhibits geographic variation in aboveground architecture and leaf traits in parallel with environmental variation in temperature and moisture. Measurements of plants occurring across a natural gradient demonstrated that plants in desert populations produce smaller, more pubescent leaves and are more compact and branched than plants in more mesic coastal environments. This phenotypic variation is interpreted in part as adaptive genetic differentiation; small size and pubescence reduce leaf temperature and thus increase water-use efficiency but at the cost of lower photosynthetic rate, which results in slower growth and more compact growth form. We explored the basis of phenotypic variation by planting seed offspring from coastal and desert populations in common gardens in both environments. Phenotypic differences among populations persisted in both common gardens, suggesting a genetic basis for trait variation. Desert offspring outperformed coastal offspring in the desert garden, suggesting superior adaptation to hot, dry conditions. Herbivore damage was greater for all offspring in the coastal garden. Phenotypic characters also showed plastic responses; all offspring had smaller, more pubescent leaves and more compact growth form in the desert garden. Our results confirm that leaf size and pubescence are heritable characters associated with pronounced variation in plant architecture.     

Quantifying natural seasonal variation in mutation parameters with mutation accumulation lines

Mutations create novel genetic variants, but their contribution to variation in fitness and other phenotypes may depend on environmental conditions. Furthermore, natural environments may be highly heterogeneous. We assessed phenotypes associated with survival and reproductive success in over 30,000 plants representing 100 mutation accumulation lines of Arabidopsis thaliana across four temporal environments at a single field site. In each of the four assays, environmental variance was substantially larger than mutational variance. For some traits, whether mutational variance was significantly varied between seasons. The founder genotype had mean trait values near the mean of the distribution of the mutation accumulation lines in all field experiments. New mutations also contributed more phenotypic variation than would be predicted, given phenotypic and sequence‐level divergence among natural populations of A. thaliana. The combination of large environmental variance with a mean effect of mutation near zero suggests that mutations could contribute substantially to standing genetic variation.     

Multiple common garden experiments suggest lack of local adaptation in an invasive ornamental plant

Aims Adaptive evolution along geographic gradients of climatic conditions is suggested to facilitate the spread of invasive plant species, leading to clinal variation among populations in the introduced range. We investigated whether adaptation to climate is also involved in the invasive spread of an ornamental shrub, Buddleja davidii, across western and central Europe.Methods We combined a common garden experiment, replicated in three climatically different central European regions, with reciprocal transplantation to quantify genetic differentiation in growth and reproductive traits of 20 invasive B. davidii populations. Additionally, we compared compensatory regrowth among populations after clipping of stems to simulate mechanical damage.Important findings Our results do not provide evidence for clinal variation among invasive B. davidii populations: populations responded similarly to the different environments, and trait values were not correlated to climatic conditions or geographic coordinates of their home sites. Moreover, we did not detect differences in the compensatory ability of populations. We suppose that the invasive spread of B. davidii has been facilitated by phenotypic plasticity rather than by adaptation to climate and that continent-wide shuffling of cultivars due to horticultural trade may have limited local adaptation so far.     

Neutral and adaptive genomic signatures of rapid poleward range expansion

Many species are expanding their range polewards, and this has been associated with rapid phenotypic change. Yet, it is unclear to what extent this reflects rapid genetic adaptation or neutral processes associated with range expansion, or selection linked to the new thermal conditions encountered. To disentangle these alternatives, we studied the genomic signature of range expansion in the damselfly Coenagrion scitulum using 4950 newly developed genomic SNPs and linked this to the rapidly evolved phenotypic differences between core and (newly established) edge populations. Most edge populations were genetically clearly differentiated from the core populations and all were differentiated from each other indicating independent range expansion events. In addition, evidence for genetic drift in the edge populations, and strong evidence for adaptive genetic variation in association with the range expansion was detected. We identified one SNP under consistent selection in four of the five edge populations and showed that the allele increasing in frequency is associated with increased flight performance. This indicates collateral, non‐neutral evolutionary changes in independent edge populations driven by the range expansion process. We also detected a genomic signature of adaptation to the newly encountered thermal regimes, reflecting a pattern of countergradient variation. The latter signature was identified at a single SNP as well as in a set of covarying SNPs using a polygenic multilocus approach to detect selection. Overall, this study highlights how a strategic geographic sampling design and the integration of genomic, phenotypic and environmental data can identify and disentangle the neutral and adaptive processes that are simultaneously operating during range expansions.     

The evolution of conditional dispersal and reproductive isolation along environmental gradients

Dispersal modulates gene flow throughout a population's spatial range. Gene flow affects adaptation at local spatial scales, and consequently impacts the evolution of reproductive isolation. A recent theoretical investigation has demonstrated that local adaptation along an environmental gradient, facilitated by the evolution of limited dispersal, can lead to parapatric speciation even in the absence of assortative mating. This and other studies assumed unconditional dispersal, so individuals start dispersing without regard to local environmental conditions. However, many species disperse conditionally; their propensity to disperse is contingent upon environmental cues, such as the degree of local crowding or the availability of suitable mates. Here, we use an individual-based model in continuous space to investigate by numerical simulation the relationship between the evolution of threshold-based conditional dispersal and parapatric speciation driven by frequency-dependent competition along environmental gradients. We find that, as with unconditional dispersal, parapatric speciation occurs under a broad range of conditions when reproduction is asexual, and under a more restricted range of conditions when reproduction is sexual. In both the asexual and sexual cases, the evolution of conditional dispersal is strongly influenced by the slope of the environmental gradient: shallow environmental gradients result in low dispersal thresholds and high dispersal distances, while steep environmental gradients result in high dispersal thresholds and low dispersal distances. The latter, however, remain higher than under unconditional dispersal, thus undermining isolation by distance, and hindering speciation in sexual populations. Consequently, the speciation of sexual populations under conditional dispersal is triggered by a steeper gradient than under unconditional dispersal. Enhancing the disruptiveness of frequency-dependent selection, more box-shaped competition kernels dramatically lower the speciation-enabling slope of the environmental gradient.     

Experimental ecology of Dryas octopetala ecotypes

Summary Growth, survival and reproduction of adult plants from three reciprocally-transplanted populations of Dryas octopetala were followed over a seven year period, from 1979 to 1986. The two most different ecotypes, from snowbed and fellfield environments, were strongly selected against in extreme foreign environments via differential growth, survival, and reproduction. The more phenotypically-plastic ecotype from the snowbed experienced rapid, early mortality on the fellfield, prior to morphological adjustment, but no mortality in the last five years of the experiment. The less plastic fellfield ecotype showed no mortality for the first three years followed by a steady decline in numbers for the next four years. The ecotypic difference in plasticity may account for the contrasting mortality patterns. A population of plants with intermediate morphology was intermediate in fitness in extreme environments, but showed no superiority in its home site, suggesting that it is not a locally-adapted population. Natural selection maintains genetic differences between snowbed and fellfield populations in the face of gene flow. However, based on the response of intermediates, there appear to be limits to the degree to which selection can fine-tune adaptation along environmental gradients.     

Detecting adaptive evolution based on association with ecological gradients: Orientation matters!

Population genetic signatures of local adaptation are frequently investigated by identifying loci with allele frequencies that exhibit high correlation with ecological variables. One difficulty with this approach is that ecological associations might be confounded by geographic variation at selectively neutral loci. Here, we consider populations that underwent spatial expansion from their original range, and for which geographical variation of adaptive allele frequency coincides with habitat gradients. Using range expansion simulations, we asked whether our ability to detect genomic regions involved in adaptation could be impacted by the orientation of the ecological gradients. For three ecological association methods tested, we found, counter-intuitively, fewer false-positive associations when ecological gradients aligned along the main axis of expansion than when they aligned along any other direction. This result has important consequences for the analysis of genomic data under non-equilibrium population genetic models. Alignment of gradients with expansion axes is likely to be common in scenarios in which expanding species track their ecological niche during climate change while adapting to changing environments at their rear edge.     

THE GENETIC STRUCTURE OF LARGE-LAKE DAPHNIA POPULATIONS

Cyclical parthenogenesis allows study of the genetic and evolutionary characteristics of groups exhibiting both asexual and sexual reproduction. The cladoceran genus Daphnia contains species which vary with respect to the relative incidence of sexual reproduction; pond species tend to undergo sexual reproduction more regularly than species found in large lakes. Previous genetic studies have focused on pond populations, generating expectations about large-lake populations that have not been fully met by recent studies. The present study of the Palearctic species Daphnia galeata further examines the genetic structure of large-lake populations. Nine local populations, from lakes in northern Germany, are examined for genetic variation at seven enzyme loci. Populations exhibit similar allelic arrays and often similar allele frequencies at the five polymorphic loci; values of Nei's genetic distance (D) ranged from 0.002 to 0.239, with a mean of 0.084. F_ST values range from 0.012 to 0.257, and spatial autocorrelation coefficients range from -0.533 to 0.551, for the eight alleles analyzed. With few exceptions, within-population genotypic frequencies were in Hardy-Weinberg equilibrium. There was, however, significant heterogeneity in genotypic frequencies among populations. The number of coexisting clonal groups, as determined by three locus genotypes, is high within populations. Clonal groups are widely distributed among localities. The amount of genetic divergence observed among these large-lake populations is smaller than that previously observed among pond populations and suggests that different processes may be important in determining the genetic structure and subsequent phenotypic divergence of lake versus pond populations.     

The genetic architecture of adaptation under migration-selection balance

Many ecologically important traits have a complex genetic basis, with the potential for mutations at many different genes to shape the phenotype. Even so, studies of local adaptation in heterogeneous environments sometimes find that just a few quantitative trait loci (QTL) of large effect can explain a large percentage of observed differences between phenotypically divergent populations. As high levels of gene flow can swamp divergence at weakly selected alleles, migration-selection-drift balance may play an important role in shaping the genetic architecture of local adaptation. Here, we use analytical approximations and individual-based simulations to explore how genetic architecture evolves when two populations connected by migration experience stabilizing selection toward different optima. In contrast to the exponential distribution of allele effect sizes expected under adaptation without migration (Orr 1998), we find that adaptation with migration tends to result in concentrated genetic architectures with fewer, larger, and more tightly linked divergent alleles. Even if many small alleles contribute to adaptation at the outset, they tend to be replaced by a few large alleles under prolonged bouts of stabilizing selection with migration. All else being equal, we also find that stronger selection can maintain linked clusters of locally adapted alleles over much greater map distances than weaker selection. The common empirical finding of QTL of large effect is shown to be expected with migration in a heterogeneous landscape, and these QTL may often be composed of several tightly linked alleles of smaller effect.     

The adaptive value of epigenetic mutation: Limited in large but high in small peripheral populations

The fate of populations during range expansions, invasions and environmental changes is largely influenced by their ability to adapt to peripheral habitats. Recent models demonstrate that stable epigenetic modifications of gene expression that occur more frequently than genetic mutations can both help and hinder adaptation in panmictic populations. However, these models do not consider interactions between epimutations and evolutionary forces in peripheral populations. Here, we use mainland–island mathematical models and simulations to explore how the faster rate of epigenetic mutation compared to genetic mutations interacts with migration, selection and genetic drift to affect adaptation in peripheral populations. Our model focuses on cases where epigenetic marks are stably inherited. In a large peripheral population, where the effect of genetic drift is negligible, our analyses suggest that epimutations with random fitness impacts that occur at rates as high as 10^–3 increase local adaptation when migration is strong enough to overwhelm divergent selection. When migration is weak relative to selection and epimutations with random fitness impacts decrease adaptation, we find epigenetic modifications must be highly adaptively biased to enhance adaptation. Finally, in small peripheral populations, where genetic drift is strong, epimutations contribute to adaptation under a wider range of evolutionary conditions. Overall, our results suggest that epimutations can change outcomes of adaptation in peripheral populations, which has implications for understanding conservation and range expansions and contractions, especially of small populations.     

Trait differentiation and adaptation of plants along elevation gradients

Studies of genetic adaptation in plant populations along elevation gradients in mountains have a long history, but there has until now been neither a synthesis of how frequently plant populations exhibit adaptation to elevation nor an evaluation of how consistent underlying trait differences across species are. We reviewed studies of adaptation along elevation gradients (i) from a meta‐analysis of phenotypic differentiation of three traits (height, biomass and phenology) from plants growing in 70 common garden experiments; (ii) by testing elevation adaptation using three fitness proxies (survival, reproductive output and biomass) from 14 reciprocal transplant experiments; (iii) by qualitatively assessing information at the molecular level, from 10 genomewide surveys and candidate gene approaches. We found that plants originating from high elevations were generally shorter and produced less biomass, but phenology did not vary consistently. We found significant evidence for elevation adaptation in terms of survival and biomass, but not for reproductive output. Variation in phenotypic and fitness responses to elevation across species was not related to life history traits or to environmental conditions. Molecular studies, which have focussed mainly on loci related to plant physiology and phenology, also provide evidence for adaptation along elevation gradients. Together, these studies indicate that genetically based trait differentiation and adaptation to elevation are widespread in plants. We conclude that a better understanding of the mechanisms underlying adaptation, not only to elevation but also to environmental change, will require more studies combining the ecological and molecular approaches.     

Nonrandom larval dispersal can steepen marine clines

Sharp and stable clinal variation is enigmatic when found in species with high gene flow. Classical population genetic models treat gene flow as a random homogenizing force countering local adaptation across habitat discontinuities. Under this view, dispersal over large spatial scales will lower the effectiveness of adaptation by natural selection at finer spatial scales. Thus, random gene flow will create a shallow phenotypic cline across an ecotone in response to a steep selection gradient. In sedentary marine species that disperse primarily as larvae, nonrandom dispersal patterns are expected due to coastal hydrodynamics. Surprisingly sharp phenotypic and genotypic clines have been documented in marine species with high gene flow. We are interested in the extent to which nonrandom dispersal could accentuate such clines. We model a linear species range in which populations have stable and uniform densities along a selection gradient; in contrast to random dispersal, convergent advection of larvae can amplify phenotypic differentiation if coupled with a semipermeable dispersal barrier in the convergence zone. The migration load caused by directional dispersal pushes the phenotypic mean away from the local trait optimum in downstream populations, that is, near the convergence zone. A dispersal barrier is possible as a result of colliding currents if the water and larvae are mostly displaced offshore, away from suitable settlement habitat. Disjunctions in a quantitative trait were enlarged in the convergence zone by faster current flows or a more complete dispersal barrier. With advection of larvae per generation one-third as far as the average dispersal distance by diffusion, convergence on a dispersal barrier with 40% permeability generated a trait disjunction across the convergence zone of two phenotypic standard deviations. Without directional dispersal, similar clines also developed across a habitat gap, where population density was low, or across dispersal barriers with less than 1% permeability. These findings suggest that the types of hydrographic phenomena often associated with marine transition zones can strongly affect the balance between gene flow and selection and generate surprisingly steep clines given the large-scale gene flow expected from larvae.