Emergence of long‐term balanced polymorphism under cyclic selection of spatially variable magnitude

A fundamental question in evolutionary biology is what promotes genetic variation at nonneutral loci, a major precursor to adaptation in changing environments. In particular, balanced polymorphism under realistic evolutionary models of temporally varying environments in finite natural populations remains to be demonstrated. Here, we propose a novel mechanism of balancing selection under temporally varying fitnesses. Using forward‐in‐time computer simulations and mathematical analysis, we show that cyclic selection that spatially varies in magnitude, such as along an environmental gradient, can lead to elevated levels of nonneutral genetic polymorphism in finite populations. Balanced polymorphism is more likely with an increase in gene flow, magnitude and period of fitness oscillations, and spatial heterogeneity. This polymorphism‐promoting effect is robust to small systematic fitness differences between competing alleles or to random environmental perturbation. Furthermore, we demonstrate analytically that protected polymorphism arises as spatially heterogeneous cyclic fitness oscillations generate a type of storage effect that leads to negative frequency dependent selection. Our findings imply that spatially variable cyclic environments can promote elevated levels of nonneutral genetic variation in natural populations.     

Local adaptation at the range peripheries of Sitka spruce

High‐dispersal rates in heterogeneous environments and historical rapid range expansion can hamper local adaptation; however, we often see clinal variation in high‐dispersal tree species. To understand the mechanisms of the species’ distribution, we investigated local adaptation and adaptive plasticity in a range‐wide context in Sitka spruce, a wind‐pollinated tree species that has recently expanded its range after glaciations. Phenotypic traits were observed using growth chamber experiments that mimicked temperature and photoperiodic regimes from the limits of the species realized niche. Bud phenology exhibited parallel reaction norms among populations; however, putatively adaptive plasticity and strong divergent selection were seen in bud burst and bud set timing respectively. Natural selection appears to have favoured genotypes that maximize growth rate during available frost‐free periods in each environment. We conclude that Sitka spruce has developed local adaptation and adaptive plasticity throughout its range in response to current climatic conditions despite generally high pollen flow and recent range expansion.     

Ecological specialization in populations adapted to constant versus heterogeneous environments*

Populations vary in their degree of ecological specialization. An intuitive, but often untested, hypothesis is that populations evolving under greater environmental heterogeneity will evolve to be less specialized. How important is environmental heterogeneity in explaining among‐population variation in specialization? We assessed juvenile viability of 20 Drosophila melanogaster populations evolving under one of four regimes: (1) a salt‐enriched environment, (2) a cadmium‐enriched environment, (3) a temporally varying environment, and (4) a spatially varying environment. Juvenile viability was tested in both the original selective environments and a set of novel environments. In both the original and novel environments, populations from the constant cadmium regime had the lowest average viability and the highest variance in viability across environments but populations from the other three regimes were similar. Our results suggest that variation in specialization among these populations is most simply explained as a pleiotropic by‐product of adaptation to specific environments rather than resulting from a history of exposure to environmental heterogeneity.     

The degree of adaptive phenotypic plasticity is correlated with the spatial environmental heterogeneity experienced by island populations of Rana temporaria

Although theoretical models have identified environmental heterogeneity as a prerequisite for the evolution of adaptive plasticity, this relationship has not yet been demonstrated experimentally. Because of pool desiccation risk, adaptation of development rate is important for many amphibians. In a simulated pool-drying experiment, we compared the development time and phenotypic plasticity in development time of populations of the common frog Rana temporaria, originating from 14 neighbouring islands off the coast of northern Sweden. Drying regime of pools used by frogs for breeding differed within and among the islands. We found that the degree of phenotypic plasticity in development time was positively correlated with the spatial variation in the pool-drying regimes present on each island. In addition, local adaptation in development time to the mean drying rate of the pools on each island was found. Hence, our study demonstrates the connection between environmental heterogeneity and developmental plasticity at the island population level, and also highlights the importance of the interplay between local specialization and phenotypic plasticity depending on the local selection pressures.     

Testing the role of phenotypic plasticity for local adaptation: growth and development in time-constrained Rana temporaria populations

Phenotypic plasticity can be important for local adaptation, because it enables individuals to survive in a novel environment until genetic changes have been accumulated by genetic accommodation. By analysing the relationship between development rate and growth rate, it can be determined whether plasticity in life-history traits is caused by changed physiology or behaviour. We extended this to examine whether plasticity had been aiding local adaptation, by investigating whether the plastic response had been fixed in locally adapted populations. Tadpoles from island populations of Rana temporaria, locally adapted to different pool-drying regimes, were monitored in a common garden. Individual differences in development rate were caused by different foraging efficiency. However, developmental plasticity was physiologically mediated by trading off growth against development rate. Surprisingly, plasticity has not aided local adaptation to time-stressed environments, because local adaptation was not caused by genetic assimilation but on selection on the standing genetic variation in development time.     

Thermal adaptation and phenotypic plasticity in a warming world: Insights from common garden experiments on Alaskan sockeye salmon

An important unresolved question is how populations of coldwater‐dependent fishes will respond to rapidly warming water temperatures. For example, the culturally and economically important group, Pacific salmon (Oncorhynchus spp.), experience site‐specific thermal regimes during early development that could be disrupted by warming. To test for thermal local adaptation and heritable phenotypic plasticity in Pacific salmon embryos, we measured the developmental rate, survival, and body size at hatching in two populations of sockeye salmon (Oncorhynchus nerka) that overlap in timing of spawning but incubate in contrasting natural thermal regimes. Using a split half‐sibling design, we exposed embryos of 10 families from each of two populations to variable and constant thermal regimes. These represented both experienced temperatures by each population, and predicted temperatures under plausible future conditions based on a warming scenario from the downscaled global climate model (MIROC A1B scenario). We did not find evidence of thermal local adaptation during the embryonic stage for developmental rate or survival. Within treatments, populations hatched within 1 day of each other, on average, and among treatments, did not differ in survival in response to temperature. We did detect plasticity to temperature; embryos developed 2.5 times longer (189 days) in the coolest regime compared to the warmest regime (74 days). We also detected variation in developmental rates among families within and among temperature regimes, indicating heritable plasticity. Families exhibited a strong positive relationship between thermal variability and phenotypic variability in developmental rate but body length and mass at hatching were largely insensitive to temperature. Overall, our results indicated a lack of thermal local adaptation, but a presence of plasticity in populations experiencing contrasting conditions, as well as family‐specific heritable plasticity that could facilitate adaptive change.     

Limited potential for adaptation to climate change in a broadly distributed marine crustacean

The extent to which acclimation and genetic adaptation might buffer natural populations against climate change is largely unknown. Most models predicting biological responses to environmental change assume that species' climatic envelopes are homogeneous both in space and time. Although recent discussions have questioned this assumption, few empirical studies have characterized intraspecific patterns of genetic variation in traits directly related to environmental tolerance limits. We test the extent of such variation in the broadly distributed tidepool copepod Tigriopus californicus using laboratory rearing and selection experiments to quantify thermal tolerance and scope for adaptation in eight populations spanning more than 17° of latitude. Tigriopus californicus exhibit striking local adaptation to temperature, with less than 1 per cent of the total quantitative variance for thermal tolerance partitioned within populations. Moreover, heat-tolerant phenotypes observed in low-latitude populations cannot be achieved in high-latitude populations, either through acclimation or 10 generations of strong selection. Finally, in four populations there was no increase in thermal tolerance between generations 5 and 10 of selection, suggesting that standing variation had already been depleted. Thus, plasticity and adaptation appear to have limited capacity to buffer these isolated populations against further increases in temperature. Our results suggest that models assuming a uniform climatic envelope may greatly underestimate extinction risk in species with strong local adaptation.     

Adaptive plasticity and niche expansion in an invasive thistle

Phenotypic differentiation in size and fecundity between native and invasive populations of a species has been suggested as a causal driver of invasion in plants. Local adaptation to novel environmental conditions through a micro‐evolutionary response to natural selection may lead to phenotypic differentiation and fitness advantages in the invaded range. Local adaptation may occur along a stress tolerance trade‐off, favoring individuals that, in benign conditions, shift resource allocation from stress tolerance to increased vigor and fecundity and, therefore, invasiveness. Alternately, the typically disturbed invaded range may select for a plastic, generalist strategy, making phenotypic plasticity the main driver of invasion success. To distinguish between these hypotheses, we performed a field common garden and tested for genetically based phenotypic differentiation, resource allocation shifts in response to water limitation, and local adaptation to the environmental gradient which describes the source locations for native and invasive populations of diffuse knapweed (Centaurea diffusa). Plants were grown in an experimental field in France (naturalized range) under water addition and limitation conditions. After accounting for phenotypic variation arising from environmental differences among collection locations, we found evidence of genetic variation between the invasive and native populations for most morphological and life‐history traits under study. Invasive C. diffusa populations produced larger, later maturing, and therefore potentially fitter individuals than native populations. Evidence for local adaptation along a resource allocation trade‐off for water limitation tolerance is equivocal. However, native populations do show evidence of local adaptation to an environmental gradient, a relationship which is typically not observed in the invaded range. Broader analysis of the climatic niche inhabited by the species in both ranges suggests that the physiological tolerances of C. diffusa may have expanded in the invaded range. This observation could be due to selection for plastic, “general‐purpose” genotypes with broad environmental tolerances.     

An experimental evolutionary study on adaptation to temporally fluctuating pH in Escherichia coli

In this study, we use the bacterium Escherichia coli to examine evolutionary responses to environmental acidity fluctuating temporally among pH 5.3, 6.3, 7.0, and 7.8 (5,000-15 nM [H(+)]). Two experimental protocols of temporal variation were used. One group (six replicate lines) of populations evolved for 2,000 generations during exposure to a cycled regime fluctuating daily between pH 5.3 and 7.8. The other group (also in six replicate lines) evolved during exposure for 2,000 generations to a randomly shifting regime fluctuating stochastically each day among pH 5.3, 6.3, 7.0, and 7.8. Adaptation to these fluctuating acidity regimes was measured as a change in fitness relative to the common ancestor by direct competition experiments in both constant and fluctuating pH regimes. For comparisons with constant pH evolution, a group evolved at a constant pH of 5.3 and another group evolved at pH 7.8 were also tested. This study initiated the first long-term laboratory natural selection experiment on adaptation to variable acidity and addressed key questions concerning patterns of adaptation (trade-offs, specialists, generalists, plasticity, transitions, and acclimation) in temporally fluctuating environments.     

Spatial and temporal dynamics in a sexual selection mosaic

Selective regimes and phenotypic optima could either change smoothly and in a clinal fashion or be spatially organized in a more unpredictable mosaic pattern over the geographic landscape. When natural or sexual selection is driven by intra- or interspecific biotic interactions, fine-grained spatial variation in selective regimes could result in selection mosaics rather than clinal variation in selection. We investigated temporal variation and spatial organization in sexual selection on male body size along an ecological coastal-inland gradient of a polymorphic damselfly Ischnura elegans. Body size increased in a clinal fashion along this gradient: animals were smaller in size at the coast, but became larger in the inland areas. In contrast, the sexual selection regimes on male body size showed evidence of more fine-grained spatial organization with no evidence for a clinal pattern and low spatial autocorrelations between populations. These spatially fine-grained sexual selection regimes varied in sign and magnitude and were driven by a combination of the densities of heritable female color morphs and local female body sizes. We suggest that the spatial organization of the selective regimes can be interpreted as a sexual selection mosaic that is influenced by highly localized density- and frequency-dependent social interactions.     

A multilocus perspective on colonization accompanied by selection and gene flow

The colonization of novel habitats involves complex interactions between founder events, selection, and ongoing migration, and can lead to diverse evolutionary outcomes from local extinction to adaptation to speciation. Although there have been several studies of the demography of colonization of remote habitats, less is known about the demographic consequences of colonization of novel habitats within a continuous species range. Populations of the Eastern Fence Lizard, Sceloporus undulatus, are continuously distributed across two dramatic transitions in substrate color in southern New Mexico and have undergone rapid adaptation following colonization of these novel environments. Blanched forms inhabit the gypsum sand dunes of White Sands and melanic forms are found on the black basalt rocks of the Carrizozo lava flow. Each of these habitats formed within the last 10,000 years, allowing comparison of genetic signatures of population history for two independent colonizations from the same source population. We present evidence on phenotypic variation in lizard color, environmental variation in substrate color, and sequence variation for mitochondrial DNA and 19 independent nuclear loci. To confirm the influence of natural selection and gene flow in this system, we show that phenotypic variation is best explained by environmental variation and that neutral genetic variation is related to distance between populations, not partitioned by habitat. The historical demography of colonization was inferred using an Approximate Bayesian Computation (ABC) framework that incorporates known geological information and allows for ongoing migration with the source population. The inferences differed somewhat between mtDNA and nuclear markers, but overall provided strong evidence of historical size reductions in both white sand and black lava populations at the time of colonization. Populations in both novel habitats appear to have undergone partial but incomplete recovery from the initial bottleneck. Both ABC analyses and measures of mtDNA sequence diversity also suggested that population reductions were more severe in the black lava compared to the white sands habitat. Differences observed between habitats may be explained by differences in colonization time, habitat geometry, and strength or response to natural selection for substrate matching. Finally, effective population size reductions in this system appear to be more dramatic when colonization is accompanied by a change in selection regime. Our analyses are consistent with a demographic cost of adaptation to novel environments and show that it is possible to infer aspects of the historical demography of local adaptation even in the presence of ongoing gene flow.     

Adaptive divergence in body size overrides the effects of plasticity across natural habitats in the brown trout

The evolution of life-history traits is characterized by trade-offs between different selection pressures, as well as plasticity across environmental conditions. Yet, studies on local adaptation are often performed under artificial conditions, leaving two issues unexplored: (i) how consistent are laboratory inferred local adaptations under natural conditions and (ii) how much phenotypic variation is attributed to phenotypic plasticity and to adaptive evolution, respectively, across environmental conditions? We reared fish from six locally adapted (domesticated and wild) populations of anadromous brown trout (Salmo trutta) in one semi-natural and three natural streams and recorded a key life-history trait (body size at the end of first growth season). We found that population-specific reaction norms were close to parallel across different streams and Q_ST was similar – and larger than F_ST – within all streams, indicating a consistency of local adaptation in body size across natural environments. The amount of variation explained by population origin exceeded the variation across stream environments, indicating that genetic effects derived from adaptive processes have a stronger effect on phenotypic variation than plasticity induced by environmental conditions. These results suggest that plasticity does not “swamp” the phenotypic variation, and that selection may thus be efficient in generating genetic change.     

Genomic signature of natural and anthropogenic stress in wild populations of the waterflea Daphnia magna: validation in space, time and experimental evolution

Natural populations are confronted with multiple selection pressures resulting in a mosaic of environmental stressors at the landscape level. Identifying the genetic underpinning of adaptation to these complex selection environments and assigning causes of natural selection within multidimensional selection regimes in the wild is challenging. The water flea Daphnia is a renowned ecological model system with its well-documented ecology, the possibility to analyse subfossil dormant egg banks and the short generation time allowing an experimental evolution approach. Capitalizing on the strengths of this model system, we here link candidate genome regions to three selection pressures, known to induce micro-evolutionary responses in Daphnia magna: fish predation, parasitism and land use. Using a genome scan approach in space, time and experimental evolution trials, we provide solid evidence of selection at the genome level under well-characterized environmental gradients in the wild and identify candidate genes linked to the three environmental stressors. Our study reveals differential selection at the genome level in Daphnia populations and provides evidence for repeatable patterns of local adaptation in a geographic mosaic of environmental stressors fuelled by standing genetic variation. Our results imply high evolutionary potential of local populations, which is relevant to understand the dynamics of trait changes in natural populations and their impact on community and ecosystem responses through eco-evolutionary feedbacks.     

Evidence for Fine-Scale Natal Homing Among Island Beach Spawning Sockeye Salmon, Oncorhynchus nerka

Salmonid fishes aggregate for breeding at spatially defined, suitable habitats. These aggregations may evolve into discrete populations when precise natal homing leads to reproductive isolation, and local regimes of selection lead to adaptation. Population structure is often defined by persistent differences in selectively neutral genetic markers and in mean values of morphological and life-history traits between locations. This approach is limited by the spatial scale at which traits diverge; low levels of reproductively successful straying, combined with similar selective pressures on life-history traits resulting from similar habitat features and environmental conditions, can significantly reduce the power of these discriminatory methods. We compared data on three life-history traits and polymorphism of DNA microsatellites for evidence of population subdivision among sockeye salmon spawning on spatially discrete but physically similar beaches on islands in Iliamna Lake, Alaska. We found small but significant differences in average body length, body depth and age composition between sites as well as significant interactions between site and year. These interactions, reflecting random variation in growth or recruitment among sites, are a powerful tool for discriminating populations with similar mean trait values. These results suggest fine-scale homing to natal sites, but the microsatellite data revealed no evidence of restricted gene flow among sites. There seems to be enough straying among the populations to prevent differentiation at neutral traits but enough homing for them to be functionally distinct.     

GENE FLOW AND SELECTION ON PHENOTYPIC PLASTICITY IN AN ISLAND SYSTEM OF RANA TEMPORARIA

Gene flow is often considered to be one of the main factors that constrains local adaptation in a heterogeneous environment. However, gene flow may also lead to the evolution of phenotypic plasticity. We investigated the effect of gene flow on local adaptation and phenotypic plasticity in development time in island populations of the common frog Rana temporaria which breed in pools that differ in drying regimes. This was done by investigating associations between traits (measured in a common garden experiment) and selective factors (pool drying regimes and gene flow from other populations inhabiting different environments) by regression analyses and by comparing pairwise F_ST values (obtained from microsatellite analyses) with pairwise Q_ST values. We found that the degree of phenotypic plasticity was positively correlated with gene flow from other populations inhabiting different environments (among‐island environmental heterogeneity), as well as with local environmental heterogeneity within each population. Furthermore, local adaptation, manifested in the correlation between development time and the degree of pool drying on the islands, appears to have been caused by divergent selection pressures. The local adaptation in development time and phenotypic plasticity is quite remarkable, because the populations are young (less than 300 generations) and substantial gene flow is present among islands.     

Modifying effects of phenotypic plasticity on interactions among natural selection, adaptation and gene flow

Divergent natural selection, adaptive divergence and gene flow may interact in a number of ways. Recent studies have focused on the balance between selection and gene flow in natural populations, and empirical work has shown that gene flow can constrain adaptive divergence, and that divergent selection can constrain gene flow. A caveat is that phenotypic diversification may be under the direct influence of environmental factors (i.e. it may be due to phenotypic plasticity), in addition to partial genetic influence. In this case, phenotypic divergence may occur between populations despite high gene flow that imposes a constraint on genetic divergence. Plasticity may dampen the effects of natural selection by allowing individuals to rapidly adapt phenotypically to new conditions, thus slowing adaptive genetic divergence. On the other hand, plasticity may promote future adaptive divergence by allowing populations to persist in novel environments. Plasticity may promote gene flow between selective regimes by allowing dispersers to adapt to alternate conditions, or high gene flow may result in the selection for increased plasticity. Here I expand frameworks for understanding relationships among selection, adaptation and gene flow to include the effects of phenotypic plasticity in natural populations, and highlight its importance in evolutionary diversification.     

Comparing the consequences of natural selection,adaptive phenotypic plasticity,and matching habitat choice for phenotype–environment matching,population genetic structure,and reproductive isolation in meta‐populations

Organisms commonly experience significant spatiotemporal variation in their environments. In response to such heterogeneity, different mechanisms may act that enhance ecological performance locally. However, depending on the nature of the mechanism involved, the consequences for populations may differ greatly. Building on a previous model that investigated the conditions under which different adaptive mechanisms (co)evolve, this study compares the ecological and evolutionary population consequences of three very different responses to environmental heterogeneity: matching habitat choice (directed gene flow), adaptive plasticity (associated with random gene flow), and divergent natural selection. Using individual‐based simulations, we show that matching habitat choice can have a greater adaptive potential than plasticity or natural selection: it allows for local adaptation while protecting genetic polymorphism despite global mating or strong environmental changes. Our simulations further reveal that increasing environmental fluctuations and unpredictability generally favor the emergence of specialist genotypes but that matching habitat choice is better at preventing local maladaptation by individuals. This confirms that matching habitat choice can speed up the genetic divergence among populations, cause indirect assortative mating via spatial clustering, and hence even facilitate sympatric speciation. This study highlights the potential importance of directed dispersal in local adaptation and speciation, stresses the difficulty of deriving its operation from nonexperimental observational data alone, and helps define a set of ecological conditions which should favor its emergence and subsequent detection in nature.     

Crawl-out behaviour in response to predation cues in an aquatic gastropod: insights from artificial selection

Local adaptation to predation often occurs in populations experiencing stable predator regimes. Under such conditions, prey species may respond by fine-tuning their behavioural defences towards a local optimum, although it is often difficult to ascertain whether such local adaptation is due to selection on fixed traits, developmental plasticity that is dependent on relatively long term exposure to environmental cues or phenotypic plasticity that can respond rapidly to a changing environment. Here we investigate whether anti-predator behaviour in two populations of the freshwater gastropod Lymnaea stagnalis responded to artificial selection. Previous work had shown that populations of this species showed a higher level of innate avoidance behaviour (crawling above the water line) in the presence of predatory fish compared with sites lacking this predation threat. By selectively breeding from high and low response selection lines, we demonstrated that this crawl-out behaviour responds rapidly to artificial selection: high response selection lines showed a significant increase and low response selection lines a significant decrease in avoidance compared with non-selected control lines. This suggests that the crawl out response in this species is heritable, and that there is potential for a response to selection in natural populations, which may produce the divergence in the plasticity of crawl out behaviour found between gastropod populations experiencing high and low predation intensity.     

Adaptation to a seasonally varying environment: a strong latitudinal cline in reproductive diapause combined with high gene flow in Drosophila montana

Adaptation to seasonal changes in the northern hemisphere includes an ability to predict the forthcoming cold season from gradual changes in environmental cues early enough to prepare for the harsh winter conditions. The magnitude and speed of changes in these cues vary between the latitudes, which induces strong selection pressures for local adaptation.We studied adaptation to seasonal changes in Drosophila montana, a northern maltfly, by defining the photoperiodic conditions leading to adult reproductive diapause along a latitudinal cline in Finland and by measuring genetic differentiation and the amount of gene flow between the sampling sites with microsatellites. Our data revealed a clear correlation between the latitude and the critical day length (CDL), in which half of the females of different cline populations enter photoperiodic reproductive diapause. There was no sign of limited gene flow between the cline populations, even though these populations showed isolation by distance. Our results show that local adaptation may occur even in the presence of high gene flow, when selection for locally adaptive life-history traits is strong. A wide range of variation in the CDLs of the fly strains within and between the cline populations may be partly due to gene flow and partly due to the opposing selection pressures for fly reproduction and overwinter survival. This variation in the timing of diapause will enhance populations' survival over the years that differ in the severity of the winter and in the length of the warm period and may also help them respond to long-term changes in environmental conditions.     

Ontogenetic changes in genetic variances of age-dependent plasticity along a latitudinal gradient

The expression of phenotypic plasticity may differ among life stages of the same organism. Age-dependent plasticity can be important for adaptation to heterogeneous environments, but this has only recently been recognized. Whether age-dependent plasticity is a common outcome of local adaptation and whether populations harbor genetic variation in this respect remains largely unknown. To answer these questions, we estimated levels of additive genetic variation in age-dependent plasticity in six species of damselflies sampled from 18 populations along a latitudinal gradient spanning 3600 km. We reared full sib larvae at three temperatures and estimated genetic variances in the height and slope of thermal reaction norms of body size at three points in time during ontogeny using random regression. Our data show that most populations harbor genetic variation in growth rate (reaction norm height) in all ontogenetic stages, but only some populations and ontogenetic stages were found to harbor genetic variation in thermal plasticity (reaction norm slope). Genetic variances in reaction norm height differed among species, while genetic variances in reaction norm slope differed among populations. The slope of the ontogenetic trend in genetic variances of both reaction norm height and slope increased with latitude. We propose that differences in genetic variances reflect temporal and spatial variation in the strength and direction of natural selection on growth trajectories and age-dependent plasticity. Selection on age-dependent plasticity may depend on the interaction between temperature seasonality and time constraints associated with variation in life history traits such as generation length.