The spatial scale of local adaptation in a stochastic environment

The distribution of phenotypes in space will be a compromise between adaptive plasticity and local adaptation increasing the fit of phenotypes to local conditions and gene flow reducing that fit. Theoretical models on the evolution of quantitative characters on spatially explicit landscapes have only considered scenarios where optimum trait values change as deterministic functions of space. Here, these models are extended to include stochastic spatially autocorrelated aspects to the environment, and consequently the optimal phenotype. Under these conditions, the regression of phenotype on the environmental variable becomes steeper as the spatial scale on which populations are sampled becomes larger. Under certain deterministic models – such as linear clines – the regression is constant. The way in which the regression changes with spatial scale is informative about the degree of phenotypic plasticity, the relative scale of effective gene flow and the environmental dependency of selection. Connections to temporal models are discussed.     

Host environment and local genetic adaptation determine phenotype in parasitic Rhinanthus angustifolius

Plants have a remarkable capacity to adapt to local environmental conditions, which can result in ecotypic differentiation. Patterns of differentiation can, however, also be influenced by the extensive phenotypic plasticity exhibited by many plant species. In this study, we evaluated the distinctness of two putative ecotypes of the parasitic herb Rhinanthus angustifolius. We compared population means of characters commonly used to distinguish between the putative ecotypes after growing individuals of R. angustifolius with a variety of host species in a common garden. Resulting data were also pooled over environments to study how phenotypic plasticity affects the distinctness of ecotypes and individual populations. Except for node number, most of the characters were plastic. The pattern of differentiation was consistent with the existence of two, or possibly three, habitat‐related groups of populations; however, we observed considerable overlap in flowering time and morphological characters after pooling data across host environments. The results show that the complex phenological and morphological variation in R. angustifolius is caused by a combination of genetically determined ecotypic differentiation and plastic responses to the host environment and other factors. © 2015 The Linnean Society of London, Botanical Journal of the Linnean Society, 2016, 180 , 89–103.     

Phenology and temperature‐dependent development of Ceutorhynchus assimilis,a potential biological control agent for Lepidium draba

Lepidium draba (Brassicaceae) is a major concern for agriculture and biodiversity in the western United States. As current control methods do not provide long‐term, sustainable solutions, research has been conducted to find biological control agents. Ceutorhynchus assimilis is one of the currently investigated candidates. Known as oligophagous in the literature, a specialist clade of this root‐galling weevil exists in southern Europe. This raised the question of its ability to survive in colder climates in the target range. We investigated the phenology of C. assimilis in the field in southern France (specialist clade) and Romania (generalist clade) and measured various temperature‐dependent parameters in the laboratory. In both ranges, weevils were univoltine. Oviposition in autumn started later in France compared to Romania, while mature larvae exited galls (to pupate in the soil) earlier the following year. On average, 25% and 32% of galls from France and Romania were completely below the soil surface, respectively, and this appeared to depend on soil substrate. Weevils transported from France to Romania were able to develop, but at a much lower rate than Romanian weevils. Mortality of overwintering larvae of both clades increased with decreasing temperature and exposure time. At ?5°C, lethal times Lt₅₀ and Lt₉₅ were 15 and 42 days for the specialist clade and 26 and 72 days for the generalist clade. A higher proportion of third instar larvae compared to first and second instar larvae survived. Pupation time at different temperatures did not differ between weevils from France or Romania. A climate match model (comparing winter temperatures) indicated that the specialist clade of C. assimilis from France has the potential to establish in some parts of the target range (e.g. Washington, Oregon, California). However, temperature extremes and winters without snow cover will likely limit its establishment unless rapid adaptive evolution takes place.     

The role of environment and core‐margin effects on range‐wide phenotypic variation in a montane grasshopper

The integration of genetic information with ecological and phenotypic data constitutes an effective approach to gain insight into the mechanisms determining interpopulation variability and the evolutionary processes underlying local adaptation and incipient speciation. Here, we use the Pyrenean Morales grasshopper (Chorthippus saulcyi moralesi) as study system to (i) analyse the relative role of genetic drift and selection in range‐wide patterns of phenotypic differentiation and (ii) identify the potential selective agents (environment, elevation) responsible for variation. We also test the hypothesis that (iii) the development of dispersal‐related traits is associated with different parameters related to population persistence/turnover, including habitat suitability stability over the last 120 000 years, distance to the species distribution core and population genetic variability. Our results indicate that selection shaped phenotypic differentiation across all the studied morphological traits (body size, forewing length and shape). Subsequent analyses revealed that among‐population differentiation in forewing length was significantly explained by a temperature gradient, suggesting an adaptive response to thermoregulation or flight performance under contrasting temperature regimes. We found support for our hypothesis predicting a positive association between the distance to the species distribution core and the development of dispersal‐related morphology, which suggests an increased dispersal capability in populations located at range edges that, in turn, exhibit lower levels of genetic variability. Overall, our results indicate that range‐wide patterns of phenotypic variation are partially explained by adaptation in response to local environmental conditions and differences in habitat persistence between core and peripheral populations.     

Latitudinal clines in body size,but not in thermal tolerance or heat‐shock cognate 70 (HSC70), in the highly‐dispersing intertidal gastropod Littorina keenae (Gastropoda: Littorinidae)

Natural populations of widely‐distributed animals often exhibit clinal variation in phenotypic traits or in allele frequencies of a particular gene over their geographical range. A planktotrophic intertidal snail, Littorina keenae is broadly distributed along the north‐eastern Pacific coast through a large latitudinal range (24°50′N–43°18′N). We tested for latitudinal clines in two complex phenotypic traits – thermal tolerance and body size – and one single locus trait – heat shock cognate 70 (HSC70) – in L. keenae along almost its entire geographical range. We found only weak evidence for a latitudinal cline in the thermal tolerance and no evidence for a cline in allele frequencies at HSC70. However, as predicted by Bergmann's rule, we detected a strong latitudinal cline that accounted for 60% of the variance in body size (R² = 0.598; P < 0.001). In contrast, body size did not significantly affect thermal tolerance. HSC70 showed no genetic differentiation among the populations, supporting our previous mitochondrial gene‐based estimate of high gene flow during this snail's free‐swimming larval stage. Given that L. keenae experiences panmixia along its species range, the observed size cline may be partially or entirely caused by a phenotypically plastic response to local thermal environments rather than by genetic divergence in body size among populations in response to locally optimizing natural selection. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 100 , 494–505.     

The relative importance of plasticity versus genetic differentiation in explaining between population differences; a meta‐analysis

Both plasticity and genetic differentiation can contribute to phenotypic differences between populations. Using data on non‐fitness traits from reciprocal transplant studies, we show that approximately 60% of traits exhibit co‐gradient variation whereby genetic differences and plasticity‐induced differences between populations are the same sign. In these cases, plasticity is about twice as important as genetic differentiation in explaining phenotypic divergence. In contrast to fitness traits, the amount of genotype by environment interaction is small. Of the 40% of traits that exhibit counter‐gradient variation the majority seem to be hyperplastic whereby non‐native individuals express phenotypes that exceed those of native individuals. In about 20% of cases plasticity causes non‐native phenotypes to diverge from the native phenotype to a greater extent than if plasticity was absent, consistent with maladaptive plasticity. The degree to which genetic differentiation versus plasticity can explain phenotypic divergence varies a lot between species, but our proxies for motility and migration explain little of this variation.     

How mechanisms of habitat preference evolve and promote divergence with gene flow

Habitat preference may promote adaptive divergence and speciation, yet the conditions under which this is likely are insufficiently explored. We use individual‐based simulations to study the evolution and consequence of habitat preference during divergence with gene flow, considering four different underlying genetically based behavioural mechanisms: natal habitat imprinting, phenotype‐dependent, competition‐dependent and direct genetic habitat preference. We find that the evolution of habitat preference generally requires initially high dispersal, is facilitated by asymmetry in population sizes between habitats, and is hindered by an increasing number of underlying genetic loci. Moreover, the probability of habitat preference to emerge and promote divergence differs greatly among the underlying mechanisms. Natal habitat imprinting evolves most easily and can allow full divergence in parameter ranges where no divergence is possible in the absence of habitat preference. The reason is that imprinting represents a one‐allele mechanism of assortative mating linking dispersal behaviour very effectively to local selection. At the other extreme, direct genetic habitat preference, a two‐allele mechanism, evolves under restricted conditions only, and even then facilitates divergence weakly. Overall, our results indicate that habitat preference can be a strong reproductive barrier promoting divergence with gene flow, but that this is highly contingent on the underlying preference mechanism.     

Population connectivity,dispersal, and swimming behavior in Daphnia

The water flea Daphnia has the capacity to respond rapidly to environmental stressors, to disperse over large geographical scales, and to preserve its genetic material by forming egg banks in the sediment. Spatial and temporal distributions of D. magna have been extensively studied over the last decades using behavioral or genetic tools, although the correlation between the two has rarely been the focus. In the present study, we therefore investigated the population genetic structure and behavioral response to a lethal threat, ultraviolet radiation (UVR), among individuals from two different water bodies. Our results show two genetic populations with moderate gene flow, highly correlated with geographical location and with inheritable traits through generations. However, despite the strong genetic differences between populations, we show homogeneous refuge demand between populations when exposed to the lethal threat solar UVR.     

Lack of local adaptation of invasive crofton weed (Ageratina adenophora) in different climatic areas of Yunnan Province,China

Aims Crofton weed, with a subtropical origin, has successfully invaded in diverse habitats that belong to different climate zones in Southwest China. We tested whether local adaptation plays an important role in the successful invasion of crofton weed in heterogeneous environments.Methods Five populations from different habitats with an altitude ranging from 678 to 2356 m were selected. Plant height, biomass, seed yield and seed germination capability of these populations were investigated in the field. Greenhouse and reciprocal transplant experiments with the five populations were conducted, and all the above characters were measured and compared among these populations.Important findings Plant height, biomass, seed yield and seed germination rate were each significantly different among the five populations in field. However, there was no difference among these populations in the greenhouse experiment. In the reciprocal transplant experiment, plants from the five populations responded similarly to different habitats in the field, indicating lack of local adaptation. Instead, phenotypic plasticity likely plays a key role in the invasion success of crofton weed in different habitats.     

Phenotypic plasticity in two marine snails: constraints superseding life history

In organisms encountering predictable environments, fixed development is expected, whereas in organisms that cannot predict their future environment, phenotypic plasticity would be optimal to increase local adaptation. To test this prediction we experimentally compared phenotypic plasticity in two rocky-shore snail species; Littorina saxatilis releasing miniature snails on the shore, and Littorina littorea releasing drifting larvae settling on various shores, expecting L. littorea to show more phenotypic plasticity than L. saxatilis. We compared magnitude and direction of vectors of phenotypic difference in juvenile shell traits after 3 months exposure to different stimuli simulating sheltered and crab-rich shores, or wave-exposed and crab-free shores. Both species showed similar direction and magnitude of vectors of phenotypic difference with minor differences only between ecotypes of the nondispersing species, indicating that plasticity is an evolving trait in L. saxatilis. The lack of a strong plastic response in L. littorea might be explained by limits rather than costs to plasticity.     

Plasticity in the timing of a major life‐history transition and resulting changes in the age structure of populations of the salamander Hynobius retardatus

Variation in age and size at life‐history transitions is a reflection of the diversifying influence of biotic or abiotic environmental change. Examples abound, but it is not well understood how such environmental changes influence the age structure of a population. I experimentally investigated the effects of water temperature and food type on age and body size at metamorphosis in larvae of the salamander Hynobius retardatus. In individuals grown at a cold temperature (15 °C) or given Chironomidae as prey, the time to metamorphosis was significantly prolonged, and body size at metamorphosis was significantly enlarged, compared with individuals grown at a warmer temperature (20 °C) or fed larvae. I also examined whether larval density (a possible indicator of cannibalism in natural habitats) generated variation in the age structure of natural populations in Hokkaido, Japan, where the climate is subarctic. Natural ponds in Hokkaido may contain larvae that have overwintered for 1 or 2 years, as well as larvae of the current year, and I found that the number of age classes was related to larval density. Although cool water temperatures prolong the larval period and induce later metamorphosis, in natural ponds diet‐based enhancement of development translated into a shorter larval duration and earlier metamorphosis. Geographic variation in the frequency of cannibalism resulted in population differences in metamorphic timing in H. retardatus larvae. It is important to understand how environmental effects are ultimately transduced through individual organisms into population‐level phenomena, with the population response arising as the summation of individual responses. Without a thorough comprehension of the mechanisms through which population and individual responses to environmental conditions are mediated, we cannot interpret the relationship between population‐level and individual‐level phenomena. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 102 , 100–114.     

Resilience in reef‐building corals: The ecological and evolutionary importance of the host response to thermal stress

Coral reefs are under extreme threat due to a number of stressors, but temperature increases due to changing climate are the most severe. Rising ocean temperatures coupled with local extremes lead to extensive bleaching, where the coral‐algal symbiosis breaks down and corals may die, compromising the structure and function of reefs. Although the symbiotic nature of the coral colony has historically been a focus of research on coral resilience, the host itself is a foundational component in the response to thermal stress. Fixed effects in the coral host set trait baselines through evolutionary processes, acting on many loci of small effect to create mosaics of thermal tolerance across latitudes and individual coral reefs. These genomic differences can be strongly heritable, producing wide variation among clones of different genotypes or families of a specific larval cross. Phenotypic plasticity is overlaid on these baselines and a growing body of knowledge demonstrates the potential for acclimatization of reef‐building corals through a variety of mechanisms that promote resilience and stress tolerance. The long‐term persistence of coral reefs will require many of these mechanisms to adjust to warmer temperatures within a generation, bridging the gap to reproductive events that allow recombination of standing diversity and adaptive change. Business‐as‐usual climate scenarios will probably lead to the loss of some coral populations or species in the future, so the interaction between intragenerational effects and evolutionary pressure is critical for the survival of reefs.     

Evolution of positive and negative density-dependent dispersal

Understanding the evolution of density-dependent dispersal strategies has been a major challenge for evolutionary ecologists. Some existing models suggest that selection should favour positive and others negative density-dependence in dispersal. Here, we develop a general model that shows how and why selection may shift from positive to negative density-dependence in response to key ecological factors, in particular the temporal stability of the environment. We find that in temporally stable environments, particularly with low dispersal costs and large group sizes, habitat heterogeneity selects for negative density-dependent dispersal, whereas in temporally variable environments, particularly with high dispersal costs and small group sizes, habitat heterogeneity selects for positive density-dependent dispersal. This shift reflects the changing balance between the greater competition for breeding opportunities in more productive patches, versus the greater long-term value of offspring that establish themselves there, the latter being very sensitive to the temporal stability of the environment. In general, dispersal of individuals out of low-density patches is much more sensitive to habitat heterogeneity than is dispersal out of high-density patches.     

Local selection modifies phenotypic divergence among Rana temporaria populations in the presence of gene flow

In ectotherms, variation in life history traits among populations is common and suggests local adaptation. However, geographic variation itself is not a proof for local adaptation, as genetic drift and gene flow may also shape patterns of quantitative variation. We studied local and regional variation in means and phenotypic plasticity of larval life history traits in the common frog Rana temporaria using six populations from central Sweden, breeding in either open‐canopy or partially closed‐canopy ponds. To separate local adaptation from genetic drift, we compared differentiation in quantitative genetic traits (Q_ST) obtained from a common garden experiment with differentiation in presumably neutral microsatellite markers (F_ST). We found that R. temporaria populations differ in means and plasticities of life history traits in different temperatures at local, and in F_ST at regional scale. Comparisons of differentiation in quantitative traits and in molecular markers suggested that natural selection was responsible for the divergence in growth and development rates as well as in temperature‐induced plasticity, indicating local adaptation. However, at low temperature, the role of genetic drift could not be separated from selection. Phenotypes were correlated with forest canopy closure, but not with geographical or genetic distance. These results indicate that local adaptation can evolve in the presence of ongoing gene flow among the populations, and that natural selection is strong in this system.     

Local adaptation in populations of a Brachionus group plicatilis cryptic species inhabiting three deep crater lakes in Central Mexico

1. Salinity is a strong selective force for many aquatic organisms, affecting both ecological and evolutionary processes. Most of our knowledge on the effects of salinity on rotifers in the Brachionus plicatilis species complex is based mainly on populations from waterbodies that experience broad environmental changes both seasonally and annually. We tested the hypothesis that, despite the supposedly high potential for gene flow among rotifers inhabiting neighbouring environments, constant salinity has promoted local adaptation, genetic population divergence and even cryptic speciation in B. plicatilis complex populations from three deep maar lakes of distinct salinities [1.1, 6.5 and 9.0 g L^?1 total dissolved solids (TDS)] in Central Mexico. 2. To look for local adaptation, we performed common garden experiments to test the effect of different salinities on population density and intrinsic growth rate (r). Then, we evaluated the genetic divergence by sequencing the cytochrome c oxidase subunit I (COI) gene and performed reproductive trials to assess the potential gene flow among the three populations and with other closely related B. plicatilis complex species. 3. We confirmed that the rotifer populations have phenotypic plasticity in tolerance of salinity, but only rotifers from the least saline lake are adapted to low salinity. Among the populations, sequence divergence at COI was very low (just a single haplotype was found), suggesting a persistent founder effect from a relatively recent single colonisation event and a subsequent dispersal from one lake to the others, and a very restricted immigration rate. In the phylogenetic analysis, rotifers from this area of Mexico clustered in the same clade with the middle‐sized species Brachionus ibericus and B. sp. ‘Almenara’. Mexican rotifers showed successful recognition, copulation and formation of hybrids among them, but interpopulation breeding with the Spanish B. ibericus and B. sp. ‘Almenara’ was unsuccessful. 4. We conclude that the B. plicatilis complex populations from these three lakes belong to a new biological species not yet described (presently named B. sp. ‘Mexico’). To our knowledge, this is the first report of local adaptation of a natural B. plicatilis complex population living in freshwater conditions (1.1 g L^?1 TDS).     

Exposure to a common environment erodes inherited between‐population trophic morphology differences in Arctic charr

Arctic charr Salvelinus alpinus juveniles reared from eggs collected from four Scottish populations showed inherited variation in their expression of trophic morphology, measured as a suite of eight characteristics of the head and mouth, before their first exogenous feeding. This demonstrated a genetic component to trophic morphology expression. During a period of 5 months following first feeding, typified by rapid growth, the differential between groups exposed to a common rearing and feeding environment was reduced significantly. It was concluded that this was the result of common environmental exposure acting on phenotypic plasticity in trophic morphology.     

High-performing plastic clones best explain the spread of yellow monkeyflower from lowland to higher elevation areas in New Zealand

The relative contribution of adaptation and phenotypic plasticity can vary between core and edge populations, with implications for invasive success. We investigated the spread of the invasive yellow monkeyflower, Erythranthe gutatta in New Zealand, where it is spreading from lowland agricultural land into high-elevation conservation areas. We investigated the extent of phenotypic variation among clones from across the South Island, looked for adaptation and compared degrees of plasticity among lowland core versus montane range-edge populations. We grew 34 clones and measured their vegetative and floral traits in two common gardens, one in the core range at 9 m a.s.l. and one near the range-edge at 560 m a.s.l. Observed trait variation was explained by a combination of genotypic diversity (as identified through common gardens) and high phenotypic plasticity. We found a subtle signature of local adaptation to lowland habitats but all clones were plastic and able to survive and reproduce in both gardens. In the range-edge garden, above-ground biomass was on average almost double and stolon length almost half that of the same clone in the core garden. Clones from low-elevation sites showed higher plasticity on average than those from higher elevation sites. The highest performing clones in the core garden were also top performers in the range-edge garden. These results suggest some highly fit general-purpose genotypes, possibly pre-adapted to New Zealand montane conditions, best explains the spread of E. gutatta from lowland to higher elevation areas.     

Evidence of adaptive divergence in plasticity: density- and site-dependent selection on shade-avoidance responses in Impatiens capensis

We investigated the conditions under which plastic responses to density are adaptive in natural populations of Impatiens capensis and determined whether plasticity has evolved differently in different selective environments. Previous studies showed that a population that evolved in a sunny site exhibited greater plasticity in response to density than did a population that evolved in a woodland site. Using replicate inbred lines in a reciprocal transplant that included a density manipulation, we asked whether such population differentiation was consistent with the hypothesis of adaptive divergence. We hypothesized that plasticity would be more strongly favored in the sunny site than in the woodland site; consequently, we predicted that selection would be more strongly density dependent in the sunny site, favoring the phenotype that was expressed at each density. Selection on internode length and flowering date was consistent with the hypothesis of adaptive divergence in plasticity. Few costs or benefits of plasticity were detected independently from the expressed phenotype, so plasticity was selected primarily through selection on the phenotype. Correlations between phenotypes and their plasticity varied with the environment and would cause indirect selection on plasticity to be environment dependent. We showed that an appropriate plastic response even to a rare environment can greatly increase genotypic fitness when that environment is favorable. Selection on the measured characters contributed to local adaptation and fully accounted for fitness differences between populations in all treatments except the woodland site at natural density.