The evolution of different maternal investment strategies in two closely related desert vertebrates

We compared egg size phenotypes and tested several predictions from the optimal egg size (OES) and bet‐hedging theories in two North American desert‐dwelling sister tortoise taxa, Gopherus agassizii and G. morafkai, that inhabit different climate spaces: relatively unpredictable and more predictable climate spaces, respectively. Observed patterns in both species differed from the predictions of OES in several ways. Mean egg size increased with maternal body size in both species. Mean egg size was inversely related to clutch order in G. agassizii, a strategy more consistent with the within‐generation hypothesis arising out of bet‐hedging theory or a constraint in egg investment due to resource availability, and contrary to theories of density dependence, which posit that increasing hatchling competition from later season clutches should drive selection for larger eggs. We provide empirical evidence that one species, G. agassizii, employs a bet‐hedging strategy that is a combination of two different bet‐hedging hypotheses. Additionally, we found some evidence for G. morafkai employing a conservative bet‐hedging strategy. (e.g., lack of intra‐ and interclutch variation in egg size relative to body size). Our novel adaptive hypothesis suggests the possibility that natural selection favors smaller offspring in late‐season clutches because they experience a more benign environment or less energetically challenging environmental conditions (i.e., winter) than early clutch progeny, that emerge under harsher and more energetically challenging environmental conditions (i.e., summer). We also discuss alternative hypotheses of sexually antagonistic selection, which arise from the trade‐offs of son versus daughter production that might have different optima depending on clutch order and variation in temperature‐dependent sex determination (TSD) among clutches. Resolution of these hypotheses will require long‐term data on fitness of sons versus daughters as a function of incubation environment, data as yet unavailable for any species with TSD.     

Bet hedging in desert winter annual plants: optimal germination strategies in a variable environment

In bet hedging, organisms sacrifice short‐term success to reduce the long‐term variance in success. Delayed germination is the classic example of bet hedging, in which a fraction of seeds remain dormant as a hedge against the risk of complete reproductive failure. Here, we investigate the adaptive nature of delayed germination as a bet hedging strategy using long‐term demographic data on Sonoran Desert winter annual plants. Using stochastic population models, we estimate fitness as a function of delayed germination and identify evolutionarily stable strategies for 12 abundant species in the community. Results indicate that delayed germination meets the criteria as a bet hedging strategy for all species. Density‐dependent models, but not density‐independent ones, predicted optimal germination strategies that correspond remarkably well with observed patterns. By incorporating naturally occurring variation in seed and seedling dynamics, our results present a rigorous test of bet hedging theory within the relevant environmental context.     

Short‐term insurance versus long‐term bet‐hedging strategies as adaptations to variable environments

Understanding how organisms adapt to environmental variation is a key challenge of biology. Central to this are bet‐hedging strategies that maximize geometric mean fitness across generations, either by being conservative or diversifying phenotypes. Theoretical models have identified environmental variation across generations with multiplicative fitness effects as driving the evolution of bet‐hedging. However, behavioral ecology has revealed adaptive responses to additive fitness effects of environmental variation within lifetimes, either through insurance or risk‐sensitive strategies. Here, we explore whether the effects of adaptive insurance interact with the evolution of bet‐hedging by varying the position and skew of both arithmetic and geometric mean fitness functions. We find that insurance causes the optimal phenotype to shift from the peak to down the less steeply decreasing side of the fitness function, and that conservative bet‐hedging produces an additional shift on top of this, which decreases as adaptive phenotypic variation from diversifying bet‐hedging increases. When diversifying bet‐hedging is not an option, environmental canalization to reduce phenotypic variation is almost always favored, except where the tails of the fitness function are steeply convex and produce a novel risk‐sensitive increase in phenotypic variance akin to diversifying bet‐hedging. Importantly, using skewed fitness functions, we provide the first model that explicitly addresses how conservative and diversifying bet‐hedging strategies might coexist.     

Ethanol resistance in Drosophila melanogaster has increased in parallel cold‐adapted populations and shows a variable genetic architecture within and between populations

Understanding the genetic properties of adaptive trait evolution is a fundamental crux of biological inquiry that links molecular processes to biological diversity. Important uncertainties persist regarding the genetic predictability of adaptive trait change, the role of standing variation, and whether adaptation tends to result in the fixation of favored variants. Here, we use the recurrent evolution of enhanced ethanol resistance in Drosophila melanogaster during this species’ worldwide expansion as a promising system to add to our understanding of the genetics of adaptation. We find that elevated ethanol resistance has evolved at least three times in different cooler regions of the species’ modern range—not only at high latitude but also in two African high‐altitude regions. Applying a bulk segregant mapping framework, we find that the genetic architecture of ethanol resistance evolution differs substantially not only between our three resistant populations, but also between two crosses involving the same European population. We then apply population genetic scans for local adaptation within our quantitative trait locus regions, and we find potential contributions of genes with annotated roles in spindle localization, membrane composition, sterol and alcohol metabolism, and other processes. We also apply simulation‐based analyses that confirm the variable genetic basis of ethanol resistance and hint at a moderately polygenic architecture. However, these simulations indicate that larger‐scale studies will be needed to more clearly quantify the genetic architecture of adaptive evolution and to firmly connect trait evolution to specific causative loci.     

Early‐life patterns of growth are linked to levels of phenotypic trait covariance and postfledging mortality across avian species

Life history studies have established that trade‐offs between growth and survival are common both within and among species. Identifying the factor(s) that mediate this trade‐off has proven difficult, however, especially at the among‐species level. In this study, we examined a series of potentially interrelated traits in a community of temperate‐zone passerine birds to help understand the putative causes and consequences of variation in early‐life growth among species. First, we examined whether nest predation risk (a proven driver of interspecific variation in growth and development rates) was correlated with species‐level patterns of incubation duration and nestling period length. We then assessed whether proxies for growth rate covaried with mean trait covariance strength (i.e., phenotypic correlations ( ^rp), which can be a marker of early‐life stress) among body mass, tarsus length, and wing length at fledging. Finally, we examined whether trait covariance strength at fledging was related to postfledging survival. We found that higher nest predation risk was correlated with faster skeletal growth and that our proxies for growth corresponded with increased trait covariance strength ( ^rp), which subsequently, correlated with higher mortality in the next life stage (postfledging period). These results provide an indication that extrinsic pressures (nest predation) impact rates of growth, and that there are costs of rapid growth across species, expressed as higher mean ^rp and elevated postfledging mortality. The link between higher levels of trait covariance at fledging and increased mortality is unclear, but increased trait covariance strength may reflect reduced phenotypic flexibility (i.e., phenotypic canalization), which may limit an organism''s capacity for coping with environmental or ecological variability.     

Neutral and adaptive drivers of genomic change in introduced brook trout (Salvelinus fontinalis) populations revealed by pooled sequencing

Understanding the drivers of successful species invasions is important for conserving native biodiversity and for mitigating the economic impacts of introduced species. However, whole‐genome resolution investigations of the underlying contributions of neutral and adaptive genetic variation in successful introductions are rare. Increased propagule pressure should result in greater neutral genetic variation, while environmental differences should elicit selective pressures on introduced populations, leading to adaptive differentiation. We investigated neutral and adaptive variation among nine introduced brook trout (Salvelinus fontinalis) populations using whole‐genome pooled sequencing. The populations inhabit isolated alpine lakes in western Canada and descend from a common source, with an average of ~19 (range of 7–41) generations since introduction. We found some evidence of bottlenecks without recovery, no strong evidence of purifying selection, and little support that varying propagule pressure or differences in local environments shaped observed neutral genetic variation differences. Putative adaptive loci analysis revealed nonconvergent patterns of adaptive differentiation among lakes with minimal putatively adaptive loci (0.001%–0.15%) that did not correspond with tested environmental variables. Our results suggest that (i) introduction success is not always strongly influenced by genetic load; (ii) observed differentiation among introduced populations can be idiosyncratic, population‐specific, or stochastic; and (iii) conservatively, in some introduced species, colonization barriers may be overcome by support through one aspect of propagule pressure or benign environmental conditions.     

Restoration of species‐rich grasslands by transfer of local plant material and its impact on species diversity and genetic variation—Findings of a practical restoration project in southeastern Germany

Restoration of species‐rich grasslands is a key issue of conservation. The transfer of seed‐containing local plant material is a proven technique to restore species‐rich grassland, since it potentially allows to establish genetically variable and locally adapted populations. In our study, we tested how the transfer of local plant material affected the species diversity and composition of restored grasslands and the genetic variation of the typical grassland plant species Knautia arvensis and Plantago lanceolata.For our study, we selected fifteen study sites in southeastern Germany. We analyzed species diversity and composition and used molecular markers to investigate genetic variation within and among populations of the study species from grasslands that served as source sites for restoration and grasslands, which were restored by transfer of green hay and threshed local plant material.The results revealed no significant differences in species diversity and composition between grasslands at source and restoration sites. Levels of genetic variation within populations of the study species Knautia arvensis and Plantago lanceolata were comparable at source and restoration sites and genetic variation among populations at source and their corresponding restoration sites were only marginal different.Our study suggests that the transfer of local plant material is a restoration approach highly suited to preserve the composition of species‐rich grasslands and the natural genetic pattern of typical grassland plant species.     

Environmental gradients of selection for an alpine-obligate bird,the white-tailed ptarmigan (Lagopus leucura)

The warming climate will expose alpine species adapted to a highly seasonal, harsh environment to novel environmental conditions. A species can shift their distribution, acclimate, or adapt in response to a new climate. Alpine species have little suitable habitat to shift their distribution, and the limits of acclimation will likely be tested by climate change in the long-term. Adaptive genetic variation may provide the raw material for species to adapt to this changing environment. Here, we use a genomic approach to describe adaptive divergence in an alpine-obligate species, the white-tailed ptarmigan (Lagopus leucura), a species distributed from Alaska to New Mexico, across an environmentally variable geographic range. Previous work has identified genetic structure and morphological, behavioral, and physiological differences across the species’ range; however, those studies were unable to determine the degree to which adaptive divergence is correlated with local variation in environmental conditions. We used a genome-wide dataset generated from 95 white-tailed ptarmigan distributed throughout the species’ range and genotype–environment association analyses to identify the genetic signature and environmental drivers of local adaptation. We detected associations between multiple environmental gradients and candidate adaptive loci, suggesting ptarmigan populations may be locally adapted to the plant community composition, elevation, local climate, and to the seasonality of the environment. Overall, our results suggest there may be groups within the species’ range with genetic variation that could be essential for adapting to a changing climate and helpful in guiding conservation action.Subject terms: Ecological genetics, Evolutionary ecology     

Playing smart vs. playing safe: the joint expression of phenotypic plasticity and potential bet hedging across and within thermal environments

Adaptive phenotypic plasticity evolves when cues reliably predict fitness consequences of life‐history decisions, whereas bet hedging evolves when environments are unpredictable. These modes of response should be jointly expressed, because environmental variance is composed of both predictable and unpredictable components. However, little attention has been paid to the joint expression of plasticity and bet hedging. Here, I examine the simultaneous expression of plasticity in germination rate and two potential bet‐hedging traits – germination fraction and within‐season diversification in timing of germination – in seeds from multiple seed families of five geographically distant populations of Lobelia inflata (L.) subjected to a thermal gradient. Populations differ in germination plasticity to temperature, in total germination fraction and in the expression of potential diversification in the timing of germination. The observation of a negative partial correlation between the expression of plasticity and germination variance (potential diversification), and a positive correlation between plasticity and germination fraction is suggestive of a trade‐off between modes of response to environmental variance. If the observed correlations are indicative of those between adaptive plasticity and bet hedging, we expect an optimal balance to exist and differ among populations. I discuss the challenges involved in testing whether the balance between plasticity and bet hedging depends on the relative predictability of environmental variance.     

Interspecific trait variability and local soil conditions modulate grassland model community responses to climate

Medium‐to‐high elevation grasslands provide critical services in agriculture and ecosystem stabilization, through high biodiversity and providing food for wildlife. However, these ecosystems face elevated risks of disruption due to predicted soil and climate changes. Separating the effects of soil and climate, however, is difficult in situ, with previous experiments focusing largely on monocultures instead of natural grassland communities. We experimentally exposed model grassland communities, comprised of three species grown on either local or reference soil, to varied climatic environments along an elevational gradient in the European Alps, measuring the effects on species and community traits. Although species‐specific biomass varied across soil and climate, species'' proportional contributions to community‐level biomass production remained consistent. Where species experienced low survivorship, species‐level biomass production was maintained through increased productivity of surviving individuals; however, maximum species‐level biomass was obtained under high survivorship. Species responded directionally to climatic variation, spatially separating differentially by plant traits (including height, reproduction, biomass, survival, leaf dry weight, and leaf area) consistently across all climates. Local soil variation drove stochastic trait responses across all species, with high levels of interactions occurring between site and species. This soil variability obscured climate‐driven responses: we recorded no directional trait responses for soil‐corrected traits like observed for climate‐corrected traits. Our species‐based approach contributes to our understanding of grassland community stabilization and suggests that these communities show some stability under climatic variation.     

Lineage‐specific plastid degradation in subtribe Gentianinae (Gentianaceae)

The structure and sequence of plastid genomes is highly conserved across most land plants, except for a minority of lineages that show gene loss and genome degradation. Understanding the early stages of plastome degradation may provide crucial insights into the repeatability and predictability of genomic evolutionary trends. We investigated these trends in subtribe Gentianinae of the Gentianaceae, which encompasses ca. 450 species distributed around the world, particularly in alpine and subalpine environments. We sequenced, assembled, and annotated the plastomes of 41 species, representing all six genera in subtribe Gentianinae and all main sections of the species‐rich genus Gentiana L. We reconstructed the phylogeny, estimated divergence times, investigated the phylogenetic distribution of putative gene losses, and related these to substitution rate shifts and species’ habitats. We obtained a strongly supported topology consistent with earlier studies, with all six genera in Gentianinae recovered as monophyletic and all main sections of Gentiana having full support. While closely related species have very similar plastomes in terms of size and structure, independent gene losses, particularly of the ndh complex, have occurred in multiple clades across the phylogeny. Gene loss was usually associated with a shift in the boundaries of the small single‐copy and inverted repeat regions. Substitution rates were variable between clades, with evidence for both elevated and decelerated rate shifts. Independent lineage‐specific loss of ndh genes occurred at a wide range of times, from Eocene to Pliocene. Our study illustrates that diverse degradation patterns shape the evolution of the plastid in this species‐rich plant group.     

Bet‐hedging germination in annual plants: a sound empirical test of the theoretical foundations

Seed dormancy is thought to be a key mechanism allowing annual plants to spread extinction risk in unpredictably varying environments. Theory predicts increasing germination fractions with increasing probability of reproductive success but solid empirical evidence is scarce and often confounded with environmental factors. Here we provide an empirical test of bet‐hedging via delayed germination for three annual plant species along a ‘predictability gradient’ in Israel. We excluded confounding environmental and maternal effects by raising inbred seed families and germinating them under controlled conditions. Additionally, we germinated field‐collected seeds in three consecutive seasons to compare their germination with inbred families where maternal effects were removed. Risk of reproductive failure was quantified using demographic data from the field and from second‐generation inbred lines raised in a rainfall gradient in the greenhouse. Our findings were consistent with bet‐hedging theory in that germination fraction was negatively related to species‐ and site‐specific risk of reproductive failure. Both field and hand‐raised seeds of one species exhibited higher dormancy with increasing risk of reproductive failure across sites, and hand‐raised seeds of another species showed the same pattern. The third species exhibited a rather random pattern of germination between years and sites, corresponding to the lack of site‐specific risk of reproductive failure. Species‐specific patterns of dormancy and risk could be related to alternative risk‐spreading strategies such as high adult survival, but were also affected by phylogeny. We provide strong empirical evidence for seed dormancy being a mechanism to reduce the risk of reproductive failure in highly variable environments, but a larger number of rigorous experimental tests of bet hedging germination are needed. Specifically, the genetic basis of bet‐hedging must be shown in species with different life histories, for demonstrating that dormancy is adaptive and how it is modified by other risk‐spreading traits.     

A comparison in species distribution model performance of succulents using key species and subsets of environmental predictors

Identifying the environmental drivers of the global distribution of succulent plants using the Crassulacean acid metabolism pathway of photosynthesis has previously been investigated through ensemble‐modeling of species delimiting the realized niche of the natural succulent biome. An alternative approach, which may provide further insight into the fundamental niche of succulent plants in the absence of dispersal limitation, is to model the distribution of selected species that are globally widespread and have become naturalized far beyond their native habitats. This could be of interest, for example, in defining areas that may be suitable for cultivation of alternative crops resilient to future climate change. We therefore explored the performance of climate‐only species distribution models (SDMs) in predicting the drivers and distribution of two widespread CAM plants, Opuntia ficus‐indica and Euphorbia tirucalli. Using two different algorithms and five predictor sets, we created distribution models for these exemplar species and produced an updated map of global inter‐annual rainfall predictability. No single predictor set produced markedly more accurate models, with the basic bioclim‐only predictor set marginally out‐performing combinations with additional predictors. Minimum temperature of the coldest month was the single most important variable in determining spatial distribution, but additional predictors such as precipitation and inter‐annual precipitation variability were also important in explaining the differences in spatial predictions between SDMs. When compared against previous projections, an a posteriori approach correctly does not predict distributions in areas of ecophysiological tolerance yet known absence (e.g., due to biotic competition). An updated map of inter‐annual rainfall predictability has successfully identified regions known to be depauperate in succulent plants. High model performance metrics suggest that the majority of potentially suitable regions for these species are predicted by these models with a limited number of climate predictors, and there is no benefit in expanding model complexity and increasing the potential for overfitting.     

Experimental evolution of bet hedging in rotifer diapause traits as a response to environmental unpredictability

The adaptive response of organisms to unpredictable environments is increasingly recognized as a central topic in fundamental and applied evolutionary ecology. Selection due to environmental unpredictability can act on multiple traits of an organism's life cycle to reduce the impact of high environmental variance. The aim of this research was to study how unpredictability selects for diapause traits: 1) the timing of sex (a proxy of the timing of diapausing egg production), and 2) the diapausing egg hatching fraction (a proxy of diapause duration). We used an experimental evolution approach with the facultative sexual rotifer Brachionus plicatilis. Laboratory populations experiencing two contrasting regimes of environmental fluctuation (predictable versus unpredictable) evolved divergently over a short time span (< 77 days). The populations under the unpredictable regime showed an earlier initiation of sexual reproduction and a lower hatching fraction of diapausing eggs than populations under the predictable regime. These findings demonstrate empirically the existence of bet‐hedging strategies in B. plicatilis regarding both traits, consistent with theoretical predictions of bet‐hedging evolution under conditions of unpredictable environmental variance. Given that scenarios of increased environmental variability are expected to occur in the near future, a comprehensive understanding of the role of bet‐hedging strategies is necessary for predicting population responses to environmental change.     

Spatially heterogeneous stochasticity and the adaptive diversification of dormancy

Diversified bet‐hedging, a strategy that leads several individuals with the same genotype to express distinct phenotypes in a given generation, is now well established as a common evolutionary response to environmental stochasticity. Life‐history traits defined as diversified bet‐hedging (e.g. germination or diapause strategies) display marked differences between populations in spatial proximity. In order to find out whether such differences can be explained by local adaptations to spatially heterogeneous environmental stochasticity, we explored the evolution of bet‐hedging dormancy strategies in a metapopulation using a two‐patch model with patch differences in stochastic juvenile survival. We found that spatial differences in the level of environmental stochasticity, restricted dispersal, increased fragmentation and intermediate survival during dormancy all favour the adaptive diversification of bet‐hedging dormancy strategies. Density dependency also plays a major role in the diversification of dormancy strategies because: (i) it may interact locally with environmental stochasticity and amplify its effects; however, (ii) it can also generate chaotic population dynamics that may impede diversification. Our work proposes new hypotheses to explain the spatial patterns of bet‐hedging strategies that we hope will encourage new empirical studies of this topic.     

Alpine lichen diversity in an isolated sky island in the Colorado Plateau,USA—Insight from an integrative biodiversity inventory

Lichens are major components of high altitude/latitude ecosystems. However, accurately characterizing their biodiversity is challenging because these regions and habitats are often underexplored, there are numerous poorly known taxonomic groups, and morphological variation in extreme environments can yield conflicting interpretations. Using an iterative taxonomic approach based on over 800 specimens and incorporating both traditional morphology‐based identifications and information from the standard fungal DNA barcoding marker, we compiled a voucher‐based inventory of biodiversity of lichen‐forming fungi in a geographically limited and vulnerable alpine community in an isolated sky island in the Colorado Plateau, USA—the La Sal Mountains. We used the newly proposed Assemble Species by Automatic Partitioning (ASAP) approach to empirically delimit candidate species‐level lineages from family‐level multiple sequence alignments. Specimens comprising DNA‐based candidate species were evaluated using traditional taxonomically diagnostic phenotypic characters to identify specimens to integrative species hypotheses and link these, where possible, to currently described species. Despite the limited alpine habitat (ca. 3,250 ha), we document the most diverse alpine lichen community known to date from the southern Rocky Mountains, with up to 240 candidate species/species‐level lineages of lichen‐forming fungi. 139 species were inferred using integrative taxonomy, plus an additional 52 candidate species within 29 different putative species complexes. Over 68% of sequences could not be assigned to species‐level rank with statistical confidence, corroborating the limited utility of current sequence repositories for species‐level DNA barcoding of lichen‐forming fungi. By integrating vouchered specimens, DNA sequence data, and photographic documentation, we provide an important baseline of lichen‐forming fungal diversity for the limited alpine habitat in the Colorado Plateau. These data provide an important resource for subsequent research in the ecology and evolution of lichens alpine habitats, including DNA barcodes for most putative species/species‐level lineages occurring in the La Sal Mountains, and vouchered collections representing any potentially undescribed species that can be used for future taxonomic studies.     

Variable maternal control of facultative egg diapause in the bushcricket Ephippiger ephippiger

1. Diapause is a term used to describe an arrest in the development of insects as adults, eggs, or embryos, and allows survival during adverse environmental conditions. 2. The influence of maternal age on the proportion of eggs that entered an initial facultative diapause and on fecundity and egg mortality was investigated in the bushcricket Ephippiger ephippiger. 3. In the absence of variable abiotic cues such as temperature and photoperiod, most E. ephippiger females increase the proportion of their eggs that enter diapause as they age, however there are large differences among females, with a minority showing different trends. Fecundity and egg survival decline with female age and also differ among females. 4. The influence of age on the incidence of diapause is likely to be an adaptive reproductive strategy, as those eggs laid later in the season are unlikely to complete egg development before the onset of winter. Females are probably hedging their bets by diversifying the rate of development among their offspring and effectively spreading development over several seasons. The low predictability of future environmental conditions might explain the heterogeneity of female reproductive strategies.     

Demographic consequences of delayed germination in two annual grasses from two locations of contrasting aridity

Delayed seed germination is considered to be a bet-hedging strategy, but experimental evidence of its adaptive role as an inherited trait is still lacking. In each of two co-occuring annual grass species, populations of Mediterranean and desert origin were studied during three consecutive years for population demography and seed germination in the reciprocally introduced experimental soil seed banks. The two environments strikingly differed in productivity (annual rainfall) and predictability (variation in amount and timing of annual rainfall). The two species exhibited highly similar pattern of seed size and dormancy across the two environments. In both species, a higher proportion of dormant seeds was observed at the desert location and for the seeds of desert origin, consistent with bet-hedging buffering against unpredictability of rainfall and high probability of drought in this environment. In addition, in both species seed mass was significantly less in plants of desert origin than in plants of Mediterranean origin. The two environments differed in demographic consequences of temporal variation in precipitation. In the Mediterranean population, even in the year of least precipitation, adults grew to maturity and seeds were produced. These seeds served to maintain population size. In contrast, in the desert population, in the year of least rainfall no seedlings survived to maturity and the soil seed bank was the only source of population persistence. Altogether, the results concur with predicted by adaptive bet hedging importance of delayed germination under marginal precipitation.     

Population differences in Chinook salmon (Oncorhynchus tshawytscha) DNA methylation: Genetic drift and environmental factors

Local adaptation and phenotypic differences among populations have been reported in many species, though most studies focus on either neutral or adaptive genetic differentiation. With the discovery of DNA methylation, questions have arisen about its contribution to individual variation in and among natural populations. Previous studies have identified differences in methylation among populations of organisms, although most to date have been in plants and model animal species. Here we obtained eyed eggs from eight populations of Chinook salmon (Oncorhynchus tshawytscha) and assayed DNA methylation at 23 genes involved in development, immune function, stress response, and metabolism using a gene‐targeted PCR‐based assay for next‐generation sequencing. Evidence for population differences in methylation was found at eight out of 23 gene loci after controlling for developmental timing in each individual. However, we found no correlation between freshwater environmental parameters and methylation variation among populations at those eight genes. A weak correlation was identified between pairwise DNA methylation dissimilarity among populations and pairwise F _ST based on 15 microsatellite loci, indicating weak effects of genetic drift or geographic distance on methylation. The weak correlation was primarily driven by two genes, GTIIBS and Nkef. However, single‐gene Mantel tests comparing methylation and pairwise F _ST were not significant after Bonferroni correction. Thus, population differences in DNA methylation are more likely related to unmeasured oceanic environmental conditions, local adaptation, and/or genetic drift. DNA methylation is an additional mechanism that contributes to among population variation, with potential influences on organism phenotype, adaptive potential, and population resilience.     

Bet hedging in a warming ocean: predictability of maternal environment shapes offspring size variation in marine sticklebacks

Bet hedging at reproduction is expected to evolve when mothers are exposed to unpredictable cues for future environmental conditions, whereas transgenerational plasticity (TGP) should be favoured when cues reliably predict the environment offspring will experience. Since climate predictions forecast an increase in both temperature and climate variability, both TGP and bet hedging are likely to become important strategies to mediate climate change effects. Here, the potential to produce variably sized offspring in both warming and unpredictable environments was tested by investigating whether stickleback (Gasterosteus aculeatus) mothers adjusted mean offspring size and within‐clutch variation in offspring size in response to experimental manipulation of maternal thermal environment and predictability (alternating between ambient and elevated water temperatures). Reproductive output traits of F1 females were influenced by both temperature and environmental predictability. Mothers that developed at ambient temperature (17 °C) produced larger, but fewer eggs than mothers that developed at elevated temperature (21 °C), implying selection for different‐sized offspring in different environments. Mothers in unpredictable environments had smaller mean egg sizes and tended to have greater within‐female egg size variability, especially at 21 °C, suggesting that mothers may have dynamically modified the variance in offspring size to spread the risk of incorrectly predicting future environmental conditions. Both TGP and diversification influenced F2 offspring body size. F2 offspring reared at 21 °C had larger mean body sizes if their mother developed at 21 °C, but this TGP benefit was not present for offspring of 17 °C mothers reared at 17 °C, indicating that maternal TGP will be highly relevant for ocean warming scenarios in this system. Offspring of variable environment mothers were smaller but more variable in size than offspring from constant environment mothers, particularly at 21 °C. In summary, stickleback mothers may have used both TGP and diversified bet‐hedging strategies to cope with the dual stress of ocean warming and environmental uncertainty.