Influence of density and salinity on larval development of salt‐adapted and salt‐naïve frog populations

Environmental change and habitat fragmentation will affect population densities for many species. For those species that have locally adapted to persist in changed or stressful habitats, it is uncertain how density dependence will affect adaptive responses. Anurans (frogs and toads) are typically freshwater organisms, but some coastal populations of green treefrogs (Hyla cinerea) have adapted to brackish, coastal wetlands. Tadpoles from coastal populations metamorphose sooner and demonstrate faster growth rates than inland populations when reared solitarily. Although saltwater exposure has adaptively reduced the duration of the larval period for coastal populations, increases in densities during larval development typically increase time to metamorphosis and reduce rates of growth and survival. We test how combined stressors of density and salinity affect larval development between salt‐adapted (“coastal”) and nonsalt‐adapted (“inland”) populations by measuring various developmental and metamorphic phenotypes. We found that increased tadpole density strongly affected coastal and inland tadpole populations similarly. In high‐density treatments, both coastal and inland populations had reduced growth rates, greater exponential decay of growth, a smaller size at metamorphosis, took longer to reach metamorphosis, and had lower survivorship at metamorphosis. Salinity only exaggerated the effects of density on the time to reach metamorphosis and exponential decay of growth. Location of origin affected length at metamorphosis, with coastal tadpoles metamorphosing slightly longer than inland tadpoles across densities and salinities. These findings confirm that density has a strong and central influence on larval development even across divergent populations and habitat types and may mitigate the expression (and therefore detection) of locally adapted phenotypes.     

Adaptive divergence vs. environmental plasticity: tracing local genetic adaptation of metamorphosis traits in salamanders

In order to assess the significance of local adaptation relative to environmental plasticity on the evolution of life history traits, we analysed the possible genetic basis of differences between pond- and stream-breeding fire salamanders (Salamandra salamandra) in Germany. These salamanders typically deposit their larvae in small streams, where they grow until they are sufficiently large to metamorphose. However, some populations in Western Germany use ponds as larval habitat. Because habitat quality of streams differs from that of ponds one expects life history differences in the pond animals, which may result either from a plastic response or through genetic differentiation (i.e. local adaptation). Using a phylogeographical analysis of mitochondrial D-loop sequences, we show that both stream and pond populations in Western Germany are derived from a single lineage that recolonized following the last glaciation. This finding suggests that pond breeding originated very recently. Our studies of habitat quality and metamorphic behaviour of larvae in natural ponds and streams disclosed that pond larvae experience a significantly reduced food supply and greater risk of drying than do stream larvae. Pond larvae metamorphose earlier at the cost of reduced mass. Common-environment experiments with pond and stream larvae show that metamorphic behaviour of pond larvae under limited-food conditions is determined genetically and is not simply a plastic response to the differing habitat conditions. These results show that phenotypic plasticity is less important than local adaptation in explaining differences in ecological diversification within this species and suggests the possibility of rapid evolution of genetic adaptations when new habitats are exploited.     

The ubiquity of phenotypic plasticity in plants: a synthesis

Adaptation to heterogeneous environments can occur via phenotypic plasticity, but how often this occurs is unknown. Reciprocal transplant studies provide a rich dataset to address this issue in plant populations because they allow for a determination of the prevalence of plastic versus canalized responses. From 31 reciprocal transplant studies, we quantified the frequency of five possible evolutionary patterns: (1) canalized response–no differentiation: no plasticity, the mean phenotypes of the populations are not different; (2) canalized response–population differentiation: no plasticity, the mean phenotypes of the populations are different; (3) perfect adaptive plasticity: plastic responses with similar reaction norms between populations; (4) adaptive plasticity: plastic responses with parallel, but not congruent reaction norms between populations; and (5) nonadaptive plasticity: plastic responses with differences in the slope of the reaction norms. The analysis included 362 records: 50.8% life‐history traits, 43.6% morphological traits, and 5.5% physiological traits. Across all traits, 52% of the trait records were not plastic, and either showed no difference in means across sites (17%) or differed among sites (83%). Among the 48% of trait records that showed some sort of plasticity, 49.4% showed perfect adaptive plasticity, 19.5% adaptive plasticity, and 31% nonadaptive plasticity. These results suggest that canalized responses are more common than adaptive plasticity as an evolutionary response to environmental heterogeneity.     

Gene expression plasticity evolves in response to colonization of freshwater lakes in threespine stickleback

Phenotypic plasticity is predicted to facilitate individual survival and/or evolve in response to novel environments. Plasticity that facilitates survival should both permit colonization and act as a buffer against further evolution, with contemporary and derived forms predicted to be similarly plastic for a suite of traits. On the other hand, given the importance of plasticity in maintaining internal homeostasis, derived populations that encounter greater environmental heterogeneity should evolve greater plasticity. We tested the evolutionary significance of phenotypic plasticity in coastal British Columbian postglacial populations of threespine stickleback (Gasterosteus aculeatus) that evolved under greater seasonal extremes in temperature after invading freshwater lakes from the sea. Two ancestral (contemporary marine) and two derived (contemporary freshwater) populations of stickleback were raised near their thermal tolerance extremes, 7 and 22 °C. Gene expression plasticity was estimated for more than 14 000 genes. Over five thousand genes were similarly plastic in marine and freshwater stickleback, but freshwater populations exhibited significantly more genes with plastic expression than marine populations. Furthermore, several of the loci shown to exhibit gene expression plasticity have been previously implicated in the adaptive evolution of freshwater populations, including a gene involved in mitochondrial regulation (PPARAa). Collectively, these data provide molecular evidence that highlights the importance of plasticity in colonization and adaptation to new environments.     

Variable salinity tolerance in ascidian larvae is primarily a plastic response to the parental environment

Both phenotypic plasticity and local genetic adaptation may contribute to a species’ ability to inhabit different environmental conditions. While phenotypic plasticity is usually considered costly, local adaptation takes generations to respond to environmental change and may be constrained by strong gene flow. The majority of marine species have complex life-cycles with pelagic stages that might be expected to promote gene flow and plastic responses, and yet several notable examples of local adaptation have been found in species with broadcast larvae. In the ascidian, Ciona intestinalis (Linnaeus, 1767),—a common marine species with broadcast spawning and a short larval stage—previous studies have found marked differences in salinity tolerance of early life-history stages among populations from different salinity regimes. We used common-garden experiments to test whether observed differences in salinity tolerance could be explained by phenotypic plasticity. Adult ascidians from two low salinity populations [2–5 m depth, ~25 practical salinity units (PSU)], and two full salinity populations (25–27 m depth, ~31 PSU) were acclimated for 2–4 weeks at both 25 and 31 PSU. Gametes were fertilized at the acclimation salinities, and the newly formed embryos were transferred to 10 different salinities (21–39 PSU) and cultured to metamorphosis. Adult acclimation salinity had an overriding and significant effect on larval metamorphic success: tolerance norms for larvae almost fully matched the acclimation salinity of the parents, independent of parental origin (deep or shallow). However we also detected minor population differences that could be attributed to either local adaptation or persistent environmental effects. We conclude that differences in salinity tolerance of C. intestinalis larvae from different populations are driven primarily by transgenerational phenotypic plasticity, a strategy that seems particularly favourable for an organism living in coastal waters where salinity is less readily predicted than in the open oceans.     

Quantitative genetic analysis of larval life history traits in two alpine populations of Rana temporaria

We estimated genetic and maternal variance components of larval life history characters in alpine populations of Rana temporaria (the common frog) using a full-sib/half-sib breeding design. We studied trait plasticity by raising tadpoles at 14 or 20°C in the laboratory. Larval period and metamorphic mass were greater at 14°C. Larval period did not differ between populations, but high elevation metamorphs were larger than low elevation metamorphs. Significant additive variation for larval period was detected in the low altitude population. No significant additive variation was detected for mass at metamorphosis (MM), which instead displayed significant maternal effects. Plasticity in metamorphic mass of froglets was greater in the high altitude population. The plastic response of larval period to temperature did not differ between the populations. Evolution of metamorphic mass is likely constrained by lack of additive genetic variation. In contrast, significant heritability for larval period suggests this trait may evolve in response to environmental change. These results differ from other studies on R. temporaria, suggesting that populations of this broadly distributed species present substantial geographic variation in the genetic architecture and plasticity of tadpole life history traits.     

Costs of phenotypic plasticity

Phenotypically plastic organisms display alternative phenotypes in different environments. It is widely appreciated that possessing alternative phenotypes can affect fitness. However, some investigators have suggested that simply carrying the ability to be plastic could also affect fitness. Evolutionary models suggest that high costs of plasticity could constrain the evolution of optimal phenotypes. However, costs (and limits) of plasticity are primarily hypothetical. Little empirical evidence exists to show that increased plasticity leads to reduced growth and development, leads to increased developmental instability, or limits the ability of organisms to produce more extreme phenotypes. I used half-sib families of larval wood frogs (Rana sylvatica) reared in outdoor mesocosms to examine how tadpoles altered behavioral, morphological, and life-historical traits in response to larval dragonfly predators (Anax longipes). The predators induced lower activity and the development of relatively large tails and small bodies in wood frogs. As a result, wood frogs experienced reduced growth and development. I then examined whether tadpole sibships with higher plasticity experienced fitness costs (above and beyond the costs of expressing a particular phenotype) and whether they were limited in producing extreme phenotypes. Fitness effects of plasticity were widespread. Depending on the trait examined and the environment experienced, increased plasticity had either positive effects, negative effects, or no effects on tadpole mass, development, and survivorship. I found no relationship between increased plasticity and greater developmental instability. There was also no evidence that sibships with increased plasticity produced less extreme phenotypes; the most extreme trait states were always produced by the most plastic genotypes. This work suggests that costs of plasticity may be pervasive in nature and may substantially impact the evolution of optimal phenotypes in organisms that live in heterogeneous environments.     

Variation in the degree and costs of adaptive phenotypic plasticity among Rana temporaria populations

Adaptive phenotypic plasticity in the form of capacity to accelerate development as a response to pond drying risk is known from many amphibian species. However, very little is known about factors that might constrain the evolution of this type of plasticity, and few studies have explored to what degree plasticity might be constrained by trade-offs dictated by adaptation to different environmental conditions. We compared the ability of southern and northern Scandinavian common frog (Rana temporaria) larvae originating from 10 different populations to accelerate their development in response to simulated pond drying risk and the resulting costs in metamorphic size in a factorial laboratory experiment. We found that (i) northern larvae developed faster than the southern larvae in all treatments, (ii) a capacity to accelerate the response was present in all five southern and all five northern populations tested, but that the magnitude of the response was much larger (and less variable) in the southern than in the northern populations, and that (iii) significant plasticity costs in metamorphic size were present in the southern populations, the plastic genotypes having smaller metamorphic size in the absence of desiccation risk, but no evidence for plasticity costs was found in the northern populations. We suggest that the weaker response to pond drying risk in the northern populations is due to stronger selection on large metamorphic size as compared with southern populations. In other words, seasonal time constraints that have selected the northern larvae to be fast growing and developing, may also constrain their innate ability for adaptive phenotypic plasticity.     

Adaptive phenotypic plasticity and genetics of larval life histories in two Rana temporaria populations

Phenotypic plasticity provides means for adapting to environmental unpredictability. In terms of accelerated development in the face of pond-drying risk, phenotypic plasticity has been demonstrated in many amphibian species, but two issues of evolutionary interest remain unexplored. First, the heritable basis of plastic responses is poorly established. Second, it is not known whether interpopulational differences in capacity to respond to pond-drying risk exist, although such differences, when matched with differences in desiccation risk would provide strong evidence for local adaptation. We investigated sources of within- and among-population variation in plastic responses to simulated pond-drying risk (three desiccation treatments) in two Rana temporaria populations originating from contrasting environments: (1) high desiccation risk with weak seasonal time constraint (southern population); and (2) low desiccation risk with severe seasonal time constraint (northern population). The larvae originating from the environment with high desiccation risk responded adaptively to the fast decreasing water treatment by accelerating their development and metamorphosing earlier, but this was not the case in the larvae originating from the environment with low desiccation risk. In both populations, metamorphic size was smaller in the high-desiccation-risk treatment, but the effect was larger in the southern population. Significant additive genetic variation in development rate was found in the northern and was nearly significant in the southern population, but there was no evidence for genetic variation in plasticity for development rates in either of the populations. No genetic variation for plasticity was found either in size at metamorphosis or growth rate. All metamorphic traits were heritable, and additive genetic variances were generally somewhat higher in the southern population, although significantly so in only one trait. Dominance variances were also significant in three of four traits, but the populations did not differ. Maternal effects in metamorphic traits were generally weak in both populations. Within-environment phenotypic correlations between larval period and metamorphic size were positive and genetic correlations negative in both populations. These results suggest that adaptive phenotypic plasticity is not a species-specific fixed trait, but evolution of interpopulational differences in plastic responses are possible, although heritability of plasticity appears to be low. The lack of adaptive response to desiccation risk in northern larvae is consistent with the interpretation that selection imposed by shorter growing season has favored rapid development in north (approximately 8% faster development in north as compared to south) or a minimum metamorphic size at the expense of phenotypic plasticity.     

Correlated evolution of phenotypic plasticity in metamorphic timing

Phenotypic plasticity has long been a focus of research, but the mechanisms of its evolution remain controversial. Many amphibian species exhibit a similar plastic response in metamorphic timing in response to multiple environmental factors; therefore, more than one environmental factor has likely influenced the evolution of plasticity. However, it is unclear whether the plastic responses to different factors have evolved independently. In this study, we examined the relationship between the plastic responses to two experimental factors (water level and food type) in larvae of the salamander Hynobius retardatus, using a cause-specific Cox proportional hazards model on the time to completion of metamorphosis. Larvae from ephemeral ponds metamorphosed earlier than those from permanent ponds when kept at a low water level or fed conspecific larvae instead of larval Chironomidae. This acceleration of metamorphosis depended only on the permanency of the larvae's pond of origin, but not on the conspecific larval density (an indicator of the frequency of cannibalism) in the ponds. The two plastic responses were significantly correlated, indicating that they may evolve correlatively. Once plasticity evolved as an adaptation to habitat desiccation, it might have relatively easily become a response to other ecological factors, such as food type via the pre-existing developmental pathway.     

Testing for beneficial reversal of dominance during salinity shifts in the invasive copepod Eurytemora affinis,and implications for the maintenance of genetic variation

Maintenance of genetic variation at loci under selection has profound implications for adaptation under environmental change. In temporally and spatially varying habitats, non‐neutral polymorphism could be maintained by heterozygote advantage across environments (marginal overdominance), which could be greatly increased by beneficial reversal of dominance across conditions. We tested for reversal of dominance and marginal overdominance in salinity tolerance in the saltwater‐to‐freshwater invading copepod Eurytemora affinis. We compared survival of F1 offspring generated by crossing saline and freshwater inbred lines (between‐salinity F1 crosses) relative to within‐salinity F1 crosses, across three salinities. We found evidence for both beneficial reversal of dominance and marginal overdominance in salinity tolerance. In support of reversal of dominance, survival of between‐salinity F1 crosses was not different from that of freshwater F1 crosses under freshwater conditions and saltwater F1 crosses under saltwater conditions. In support of marginal overdominance, between‐salinity F1 crosses exhibited significantly higher survival across salinities relative to both freshwater and saltwater F1 crosses. Our study provides a rare empirical example of complete beneficial reversal of dominance associated with environmental change. This mechanism might be crucial for maintaining genetic variation in salinity tolerance in E. affinis populations, allowing rapid adaptation to salinity changes during habitat invasions.     

Distribution of genetic variation underlying adult migration timing in steelhead of the Columbia River basin

Fish migrations are energetically costly, especially when moving between freshwater and saltwater, but are a viable strategy for Pacific salmon and trout (Oncorhynchus spp.) due to the advantageous resources available at various life stages. Anadromous steelhead (O. mykiss) migrate vast distances and exhibit variation for adult migration phenotypes that have a genetic basis at candidate genes known as greb1L and rock1. We examined the distribution of genetic variation at 13 candidate markers spanning greb1L, intergenic, and rock1 regions versus 226 neutral markers for 113 populations (n = 9,471) of steelhead from inland and coastal lineages in the Columbia River. Patterns of population structure with neutral markers reflected genetic similarity by geographic region as demonstrated in previous studies, but candidate markers clustered populations by genetic variation associated with adult migration timing. Mature alleles for late migration had the highest frequency overall in steelhead populations throughout the Columbia River, with only 9 of 113 populations that had a higher frequency of premature alleles for early migration. While a single haplotype block was evident for the coastal lineage, we identified multiple haplotype blocks for the inland lineage. The inland lineage had one haplotype block that corresponded to candidate markers within the greb1L gene and immediately upstream in the intergenic region, and the second block only contained candidate markers from the intergenic region. Haplotype frequencies had similar patterns of geographic distribution as single markers, but there were distinct differences in frequency between the two haplotype blocks for the inland lineage. This may represent multiple recombination events that differed between lineages where phenotypic differences exist between freshwater entry versus arrival timing as indicated by Micheletti et al. (2018a). Redundancy analyses were used to model environmental effects on allelic frequencies of candidate markers, and significant variables were migration distance, temperature, isothermality, and annual precipitation. This study improves our understanding of the spatial distribution of genetic variation underlying adult migration timing in steelhead as well as associated environmental factors and has direct conservation and management implications.     

Does local adaptation along a latitudinal cline shape plastic responses to combined thermal and nutritional stress?

Thermal and nutritional stress are commonly experienced by animals. This will become increasingly so with climate change. Whether populations can plastically respond to such changes will determine their survival. Plasticity can vary among populations depending on the extent of environmental heterogeneity. However, theory conflicts as to whether environmental heterogeneity should increase or decrease plasticity. Using three locally adapted populations of Drosophila melanogaster sampled from a latitudinal gradient, we investigated whether plastic responses to combinations of nutrition and temperature increase or decrease with latitude for four traits: egg-adult viability, egg-adult development time, and two body size traits. Employing nutritional geometry, we reared larvae on 25 diets varying in protein and carbohydrate content at two temperatures: 18 and 25°C. Plasticity varied among traits and across the three populations. Viability was highly canalized in all three populations. The tropical population showed the least plasticity for development time, the sub-tropical showed the highest plasticity for wing area, and the temperate population showed the highest plasticity for femur length. We found no evidence of latitudinal plasticity gradients in either direction. Our data highlight that differences in thermal variation and resource predictability experienced by populations along a latitudinal cline are not sufficient to predict their plasticity.     

Drought‐induced saltwater incursion leads to increased wetland nitrogen export

Coastal wetlands have the capacity to retain and denitrify large quantities of reactive nitrogen (N), making them important in attenuating increased anthropogenic N flux to coastal ecosystems. The ability of coastal wetlands to retain and transform N is being reduced by wetland losses resulting from land development. Nitrogen retention in coastal wetlands is further threatened by the increasing frequency and spatial extent of saltwater inundation in historically freshwater ecosystems, due to the combined effects of dredging, declining river discharge to coastal areas due to human water use, increased drought frequency, and accelerating sea‐level rise. Because saltwater incursion may affect N cycling through multiple mechanisms, the impacts of salinization on coastal freshwater wetland N retention and transformation are not well understood. Here, we show that repeated annual saltwater incursion during late summer droughts in the coastal plain of North Carolina changed N export from organic to inorganic forms and led to a doubling of annual NH₄⁺ export from a 440 hectare former agricultural field undergoing wetland restoration. Soil solution NH₄⁺ concentrations in two mature wetlands also increased with salinization, but the magnitude of increase was smaller than that in the former agricultural field. Long‐term saltwater exposure experiments with intact soil columns demonstrated that much of the increase in reactive N released could be explained by exchange of salt cations with sediment NH₄⁺. Using these findings together with the predicted flooding of 1661 km² of wetlands along the NC coast by 2100, we estimate that saltwater incursion into these coastal areas could release up to 18 077 Mg N, or approximately half the annual NH₄⁺ flux of the Mississippi River. Our results suggest that saltwater incursion into coastal freshwater wetlands globally could lead to increased N loading to sensitive coastal waters.     

Evolution of growth by genetic accommodation in Icelandic freshwater stickleback

Classical Darwinian adaptation to a change in environment can ensue when selection favours beneficial genetic variation. How plastic trait responses to new conditions affect this process depends on how plasticity reveals to selection the influence of genotype on phenotype. Genetic accommodation theory predicts that evolutionary rate may sharply increase when a new environment induces plastic responses and selects on sufficient genetic variation in those responses to produce an immediate evolutionary response, but natural examples are rare. In Iceland, marine threespine stickleback that have colonized freshwater habitats have evolved more rapid individual growth. Heritable variation in growth is greater for marine full-siblings reared at low versus high salinity, and genetic variation exists in plastic growth responses to low salinity. In fish from recently founded freshwater populations reared at low salinity, the plastic response was strongly correlated with growth. Plasticity and growth were not correlated in full-siblings reared at high salinity nor in marine fish at either salinity. In well-adapted lake populations, rapid growth evolved jointly with stronger plastic responses to low salinity and the persistence of strong plastic responses indicates that growth is not genetically assimilated. Thus, beneficial plastic growth responses to low salinity have both guided and evolved along with rapid growth as stickleback adapted to freshwater.     

Consequences of genetic erosion on fitness and phenotypic plasticity in European tree frog populations (Hyla arborea)

The detrimental effects of genetic erosion on small isolated populations are widely recognized contrary to their interactions with environmental changes. The ability of genotypes to plastically respond to variability is probably essential for the persistence of these populations. Genetic erosion impact may be exacerbated if inbreeding affects plastic responses or if their maintenance were at higher phenotypic costs. To understand the interplay 'genetic erosion-fitness-phenotypic plasticity', we experimentally compared, in different environments, the larval performances and plastic responses to predation of European tree frogs (Hyla arborea) from isolated and connected populations. Tadpoles from isolated populations were less performant, but the traits affected were environmental dependant. Heterosis observed in crosses between isolated populations allowed attributing their low fitness to inbreeding. Phenotypic plasticity can be maintained in the face of genetic erosion as inducible defences in response to predator were identical in all populations. However, the higher survival and developmental costs for isolated populations in harsh conditions may lead to an additional fitness loss for isolated populations.     

Variation in gene expression along a salinity gradient in wild populations of the euryhaline black-chinned tilapia Sarotherodon melanotheron

This study evaluated variation in expression of 11 genes within and among six wild populations of the black-chinned tilapia Sarotherodon melanotheron distributed along a salinity gradient from 0 to 100. Previous laboratory studies had shown that expression of these genes was sensitive to water salinity; the current study confirmed that a number of them also varied in expression in wild populations along the salinity gradient. Principal component analysis (PCA) first distinguished two, not mutually exclusive, sets of genes: trade-off genes that were highly expressed at one or other extreme of the salinity gradient and stress genes that were up-regulated at the two salinity extremes (i.e. a U-shaped expression pattern). The PCA clearly partitioned the populations into three groups based on their gene expression patterns and their position along the salinity gradient: a freshwater (GL; 0) population, four brackish and seawater (GB, HB, SM, SF; ranging from 20 to 50) populations and a hypersaline (SK, 100) population. Individual variation in gene expression was significantly greater within the populations at the extreme compared to intermediate salinities. These results reveal phenotypically plastic regulation of gene expression in S. melanotheron, and greater osmoregulatory and plasticity costs at extreme salinities, where fitness-related traits are known to be altered.     

Experimentally induced host‐shift changes life‐history strategy in a seed beetle

Expansion of the host range in phytophagous insects depends on their ability to form an association with a novel plant through changes in host‐related traits. Phenotypic plasticity has important effects on initial survival of individuals faced with a new plant, as well as on the courses of evolutionary change during long‐term adaptation to novel conditions. Using experimental populations of the seed beetle that evolved on ancestral (common bean) or novel (chickpea) host and applying reciprocal transplant at both larval and adult stage on the alternative host plant, we studied the relationship between the initial (plastic) phases of host‐shift and the subsequent stages of evolutionary divergence in life‐history strategies between populations exposed to the host‐shift process. After 48 generations, populations became well adapted to chickpea by evolving the life‐history strategy with prolonged larval development, increased body mass, earlier reproduction, shorter lifespan and decreased plasticity of all traits compared with ancestral conditions. In chickpea‐adapted beetles, negative fitness consequences of low plasticity of pre‐adult development (revealed as severe decrease in egg‐to‐adult viability on beans) exhibited mismatch with positive effects of low plasticity (i.e. low host sensitivity) in oviposition and fecundity. In contrast, beetles adapted to the ancestral host showed high plasticity of developmental process, which enabled high larval survival on chickpea, whereas elevated plasticity in adult behaviour (i.e. high host sensitivity) resulted in delayed reproduction and decreased fecundity on chickpea. The analysis of population growth parameters revealed significant fluctuation during successive phases of the host‐shift process in A. obtectus.     

Plasticity of plant defense and its evolutionary implications in wild populations of Boechera stricta

Phenotypic plasticity is thought to impact evolutionary trajectories by shifting trait values in a direction that is either favored by natural selection (“adaptive” plasticity) or disfavored (“nonadaptive” plasticity). However, it is unclear how commonly each of these types of plasticity occurs in natural populations. To answer this question, we measured glucosinolate defensive chemistry and reproductive fitness in over 1500 individuals of the wild perennial mustard Boechera stricta, planted in four common gardens across central Idaho, United States. Glucosinolate profiles—including total glucosinolate concentration as well as the relative abundances and overall diversity of different compounds—were strongly plastic both among habitats and within habitats. Patterns of glucosinolate plasticity varied greatly among genotypes. Plasticity among sites was predicted to affect fitness in 27.1% of cases; more often than expected by chance, glucosinolate plasticity increased rather than decreased relative fitness. In contrast, we found no evidence for within‐habitat selection on glucosinolate reaction norm slopes (i.e., plasticity along a continuous environmental gradient). Together, our results indicate that glucosinolate plasticity may improve the ability of B. stricta populations to persist after migration to new habitats.