An ecologist's guide for studying DNA methylation variation in wild vertebrates

Variable performance of widely used bisulfite sequencing methods and read mapping software for DNA methylation

DNA methylation (DNAm) is the most commonly studied marker in ecological epigenetics, yet the performance of popular library preparation strategies and bioinformatic tools is seldom assessed in genetically variable natural populations. We profiled DNAm using reduced representation bisulfite sequencing (RRBS) and whole genome bisulfite sequencing (WGBS) of technical and biological replicates from threespine stickleback ( Gasterosteus aculeatus ) liver tissue. We then compared how the most commonly used methylation caller (Bismark) performed relative to two alternative pipelines (BWA mem or BWA meth read mappers analyzed with MethylDackel). BWA meth provided 50% and 45% higher mapping efficiency than BWA mem and Bismark, respectively. Despite differences in mapping efficiency, BWA meth and Bismark produced similar methylation profiles, while BWA mem systematically discarded unmethylated cytosines. Depth filters had large impacts on CpG sites recovered across multiple individuals, particularly with WGBS data. Notably, the prevalence of CpG sites with intermediate methylation levels is greatly reduced in RRBS, which may have important consequences for functional interpretations. We conclude by discussing how library construction and bisulfite sequence alignment software can influence the abundance and reliability of data available for downstream analysis. Our analyses suggest that researchers studying genetically variable populations will benefit from deeply sequencing a few initial individuals to identify the amount of genomic coverage necessary for mean methylation estimates to plateau, a value that may differ by species and population. We additionally advocate for paired end sequencing on RRBS libraries to filter SNPs that may bias methylation metrics, which is counter to conventional wisdom.

The interplay between epigenomic and transcriptomic variation during ecotype divergence in stickleback

BACKGROUND

Populations colonizing contrasting environments are likely to undergo adaptive divergence and evolve ecotypes with locally adapted phenotypes. While diverse molecular mechanisms underlying ecotype divergence have been identified, less is known about their interplay and degree of divergence.

RESULTS

Here we integrated epigenomic and transcriptomic data to explore the interactions among gene expression, alternative splicing, DNA methylation, and microRNA expression to gauge the extent to which patterns of divergence at the four molecular levels are aligned in a case of postglacial divergence between marine and freshwater ecotypes of nine-spined sticklebacks (Pungitius pungitius). Despite significant genome-wide associations between epigenomic and transcriptomic variation, we found largely non-parallel patterns of ecotype divergence across epigenomic and transcriptomic levels, with predominantly nonoverlapping (ranging from 43.40 to 87.98%) sets of differentially expressed, spliced and methylated genes, and candidate genes targeted by differentially expressed miRNA between the ecotypes. Furthermore, we found significant variation in the extent of ecotype divergence across different molecular mechanisms, with differential methylation and differential splicing showing the highest and lowest extent of divergence between ecotypes, respectively. Finally, we found a significant enrichment of genes associated with ecotype divergence in differential methylation.

CONCLUSIONS

Our results suggest a nuanced relationship between epigenomic and transcriptomic processes, with alignment at the genome-wide level masking relatively independent effects of different molecular mechanisms on ecotype divergence at the gene level.

A Palearctic view of a bat fungal disease

The fungal infection causing white-nose disease in hibernating bats in North America has resulted in dramatic population declines of affected species, since the introduction of the causative agent Pseudogymnoascus destructans. The fungus is native to the Palearctic, where it also infects several bat species, yet rarely causes severe pathology or the death of the host. Pseudogymnoascus destructans infects bats during hibernation by invading and digesting the skin tissue, resulting in the disruption of torpor patterns and consequent emaciation. Relations among pathogen, host, and environment are complex, and individuals, populations, and species respond to the fungal pathogen in different ways. For example, the Nearctic Myotis lucifugus responds to infection by mounting a robust immune response, leading to immunopathology often contributing to mortality. In contrast, the Palearctic M. myotis shows no significant immunological response to infection. This lack of a strong response, resulting from the long coevolution between the hosts and the pathogen in the pathogen's native range, likely contributes to survival in tolerant species. After more than 15 years since the initial introduction of the fungus to North America, some of the affected populations are showing signs of recovery, suggesting that the fungus, hosts, or both are undergoing processes that may eventually lead to coexistence. The suggested or implemented management methods of the disease in North America have encompassed, for example, the use of probiotics and fungicides, vaccinations, and modifying the environmental conditions of the hibernation sites to limit the growth of the pathogen, intensity of infection, or the hosts' responses to it. Based on current knowledge from Eurasia, policy makers and conservation managers should refrain from disrupting the ongoing evolutionary processes and adopt a holistic approach to managing the epizootic.

Epigenetic variation in light of population genetic practice

The evolutionary impact of epigenetic variation depends on its transgenerational stability and source - whether genetically determined, environmentally induced, or due to spontaneous, genotype-independent mutations. Here, we evaluate current approaches for investigating an independent role of epigenetics in evolution, pinpointing methodological challenges. We further identify opportunities arising from integrating epigenetic data with population genetic analyses in natural populations. Efforts to advance data quality, study design, and statistical treatment are encouraged to consolidate our understanding of the source of heritable epigenetic variation, quantify its autonomous potential for evolution, and enrich population genetic analyses with an additional layer of information.

Benchmarking sample pooling for epigenomics of natural populations

DNA methylation (DNAm) is a mechanism for rapid acclimation to environmental conditions. In natural systems, small effect sizes relative to noise necessitates large sampling efforts to detect differences. Large numbers of individually sequenced libraries are costly. Pooling DNA prior to library preparation may be an efficient way to reduce costs and increase sample size, yet there are to date no recommendations in ecological epigenetics research. We test whether pooled and individual libraries yield comparable DNAm signals in a natural system exposed to different pollution levels by generating whole-epigenome data from two invasive molluscs (Corbicula fluminea, Dreissena polymorpha) collected from polluted and unpolluted localities (Italy). DNA of the same individuals were used for pooled and individual epigenomic libraries and sequenced with equivalent resources per individual. We found that pooling effectively captures similar genome-wide and global methylation signals as individual libraries, highlighting that pooled libraries are representative of the global population signal. However, pooled libraries yielded orders of magnitude more data than individual libraries, which was a consequence of higher coverage. We would therefore recommend aiming for a high initial coverage of individual libraries (15×) in future studies. Consequently, we detected many more differentially methylated regions (DMRs) with the pooled libraries and a significantly lower statistical power for regions from individual libraries. Computationally pooled data from the individual libraries produced fewer DMRs and the overlap with wet-lab pooled DMRs was relatively low. We discuss possible causes for discrepancies, list benefits and drawbacks of pooling, and provide recommendations for future epigenomic studies.

Changes in global methylation patterns of Mytilus galloprovincialis exposed to microplastics

Microplastics (MPs) disturb the normal activity of aquatic organisms at different levels, causing physiological stress and altering feeding, growth, and reproduction. Alterations of epigenetic patterns due to exposure to MPs have scarcely been studied in invertebrates. In this study, Mytilus galloprovincialis mussels (N = 61) were intermittently exposed to different concentrations of pure polystyrene microbeads for three weeks. The concentrations used in this research were similar to those currently found in certain polluted environments (E1), as well as higher doses to which mussels could be further exposed (E2 and E3). After exposure period, the global methylation patterns were investigated using Amplified Fragment Length Polymorphism (AFLPs). Significantly lower methylation was found in exposed groups compared to the control group. The level of hypomethylation increased with the concentration of microbeads. Similar results were found from field samples inhabiting two sites differentially MPs-polluted. The implications of this discovery were analysed and discussed, noting the already known effects of MPs on metabolism and cell division. Further studies on this and other sentinel organisms are recommended to understand the response of the aquatic species to the currently increasing MPs pollution.

Corticosterone exposure is associated with long-term changes in DNA methylation, physiology and breeding decisions in a wild bird

When facing challenges, vertebrates activate a hormonal stress response that can dramatically alter behaviour and physiology. Although this response can be costly, conceptual models suggest that it can also recalibrate the stress response system, priming more effective responses to future challenges. Little is known about whether this process occurs in wild animals, particularly in adulthood, and if so, how information about prior experience with stressors is encoded. One potential mechanism is hormonally mediated changes in DNA methylation. We simulated the spikes in corticosterone that accompany a stress response using non-invasive dosing in tree swallows (Tachycineta bicolor) and monitored the phenotypic effects 1 year later. In a subset of individuals, we characterized DNA methylation using reduced representation bisulfite sequencing shortly after treatment and a year later. The year after treatment, experimental females had stronger negative feedback and initiated breeding earlier-traits that are associated with stress resilience and reproductive performance in our population-and higher baseline corticosterone. We also found that natural variation in corticosterone predicted patterns of DNA methylation. Finally, corticosterone treatment influenced methylation on short (1-2 weeks) and long (1 year) time scales; however, these changes did not have clear links to functional regulation of the stress response. Taken together, our results are consistent with corticosterone-induced priming of future stress resilience and support DNA methylation as a potential mechanism, but more work is needed to demonstrate functional consequences. Uncovering the mechanisms linking experience with the response to future challenges has implications for understanding the drivers of stress resilience.

Non-concordant epigenetic and transcriptional responses to acute thermal stress in western mosquitofish (Gambusia affinis)

Climate change is intensifying the frequency and severity of extreme temperatures. Understanding the molecular mechanisms underlying the ability to cope with acute thermal stress is key for predicting species' responses to extreme temperature events. While many studies have focused on the individual roles of gene expression, post-transcriptional processes and epigenetic modifications in response to acute thermal stress, the relative contribution of these molecular mechanisms remains unclear. The wide range of thermal limits of western mosquitofish (Gambusia affinis) provides an opportunity to explore this interplay. Here, we quantified changes in gene expression, alternative splicing, DNA methylation and microRNA (miRNA) expression in muscle tissue dissected from mosquitofish immediately after reaching high (CTmax) or low thermal limit (CTmin). Although the numbers of genes showing expression and splicing changes in response to acute temperature stress were small, we found a possibly larger and non-redundant role of splicing compared to gene expression, with more genes being differentially spliced (DSGs) than differentially expressed (DEGs), and little overlap between DSGs and DEGs. We also identified a small proportion of CpGs showing significant methylation change (i.e. differentially methylated cytosines, DMCs) in fish at thermal limits; however, there was no overlap between DEGs and genes annotated with DMCs in both CTmax and CTmin experiments. The weak interplay between epigenetic modifications and gene expression was further supported by our discoveries of no differentially expressed miRNAs. These findings provide novel insights into the relative role of different molecular mechanisms underlying immediate responses to extreme temperatures and demonstrate non-concordant responses of epigenetic and transcriptional mechanisms to acute temperature stress.

Early-life environmental effects on birds: epigenetics and microbiome as mechanisms underlying long-lasting phenotypic changes

Although the long-lasting effects of variation in early-life environment have been well documented across organisms, the underlying causal mechanisms are only recently starting to be unraveled. Yet understanding the underlying mechanisms of long-lasting effects can help us predict how organisms will respond to changing environments. Birds offer a great system in which to study developmental plasticity and its underlying mechanisms owing to the production of large external eggs and variation in developmental trajectories, combined with a long tradition of applied, physiological, ecological and evolutionary research. Epigenetic changes (such as DNA methylation) have been suggested to be a key mechanism mediating long-lasting effects of the early-life environment across taxa. More recently, changes in the early-life gut microbiome have been identified as another potential mediator of developmental plasticity. As a first step in understanding whether these mechanisms contribute to developmental plasticity in birds, this Review summarizes how changes in early-life environment (both prenatal and postnatal) influence epigenetic markers and the gut microbiome. The literature shows how both early-life biotic (such as resources and social environment) and abiotic (thermal environment and various anthropogenic stressors) factors modify epigenetic markers and the gut microbiome in birds, yet data concerning many other environmental factors are limited. The causal links of these modifications to lasting phenotypic changes are still scarce, but changes in the hypothalamic-pituitary-adrenal axis have been identified as one putative pathway. This Review identifies several knowledge gaps, including data on the long-term effects, stability of the molecular changes, and lack of diversity in the systems studied, and provides directions for future research.

Spatial covariation between genetic and epigenetic diversity in wild plant and animal populations: a meta-analysis

Epigenetic variation may be crucial in understanding the structure of wild populations, thereby aiding in their management and conservation. However, the relationship between epigenetic and genetic variation remains poorly understood, especially in wild populations. To address this, we conducted a meta-analysis of studies that examined the genetic and epigenetic structures of wild plant and animal populations. We aimed to determine whether epigenetic variation is spatially independent of genetic variation in the wild and to highlight the conditions under which epigenetic variation might be informative. We show a significant positive correlation between genetic and epigenetic pairwise differentiation, indicating that in wild populations, epigenetic diversity is closely linked to genetic differentiation. The correlation was weaker for population pairs that were weakly differentiated genetically, suggesting that in such cases, epigenetic marks might be independent of genetic marks. Additionally, we found that global levels of genetic and epigenetic differentiation were similar across plant and animal populations, except when populations were weakly differentiated genetically. In such cases, epigenetic differentiation was either higher or lower than genetic differentiation. Our results suggest that epigenetic information is particularly relevant in populations that have recently diverged genetically or are connected by gene flow. Future studies should consider the genetic structure of populations when inferring the role of epigenetic diversity in local adaptation in wild populations. Furthermore, there is a need to identify the factors that sustain the links between genetic and epigenetic diversity to improve our understanding of the interplay between these two forms of variation in wild populations.

Potential Role of DNA Methylation as a Driver of Plastic Responses to the Environment Across Cells, Organisms, and Populations

There is great interest in exploring epigenetic modifications as drivers of adaptive organismal responses to environmental change. Extending this hypothesis to populations, epigenetically driven plasticity could influence phenotypic changes across environments. The canonical model posits that epigenetic modifications alter gene regulation and subsequently impact phenotypes. We first discuss origins of epigenetic variation in nature, which may arise from genetic variation, spontaneous epimutations, epigenetic drift, or variation in epigenetic capacitors. We then review and synthesize literature addressing three facets of the aforementioned model: (i) causal effects of epigenetic modifications on phenotypic plasticity at the organismal level, (ii) divergence of epigenetic patterns in natural populations distributed across environmental gradients, and (iii) the relationship between environmentally induced epigenetic changes and gene expression at the molecular level. We focus on DNA methylation, the most extensively studied epigenetic modification. We find support for environmentally associated epigenetic structure in populations and selection on stable epigenetic variants, and that inhibition of epigenetic enzymes frequently bears causal effects on plasticity. However, there are pervasive confounding issues in the literature. Effects of chromatin-modifying enzymes on phenotype may be independent of epigenetic marks, alternatively resulting from functions and protein interactions extrinsic of epigenetics. Associations between environmentally induced changes in DNA methylation and expression are strong in plants and mammals but notably absent in invertebrates and nonmammalian vertebrates. Given these challenges, we describe emerging approaches to better investigate how epigenetic modifications affect gene regulation, phenotypic plasticity, and divergence among populations.

Genomic Approaches Are Improving Taxonomic Representation in Genetic Studies of Speciation

Until recently, our understanding of the genetics of speciation was limited to a narrow group of model species with a specific set of characteristics that made genetic analysis feasible. Rapidly advancing genomic technologies are eliminating many of the distinctions between laboratory and natural systems. In light of these genomic developments, we review the history of speciation genetics, advances that have been gleaned from model and non-model organisms, the current state of the field, and prospects for broadening the diversity of taxa included in future studies. Responses to a survey of speciation scientists across the world reveal the ongoing division between the types of questions that are addressed in model and non-model organisms. To bridge this gap, we suggest integrating genetic studies from model systems that can be reared in the laboratory or greenhouse with genomic studies in related non-models where extensive ecological knowledge exists.

The Evolutionary Complexities of DNA Methylation in Animals: From Plasticity to Genetic Evolution

The importance of DNA methylation in plastic responses to environmental change and evolutionary dynamics is increasingly recognized. Here, we provide a Perspective piece on the diverse roles of DNA methylation on broad evolutionary timescales, including (i) short-term transient acclimation, (ii) stable phenotypic evolution, and (iii) genomic evolution. We show that epigenetic responses vary along a continuum, ranging from short-term acclimatory responses in variable environments within a generation to long-term modifications in populations and species. DNA methylation thus unlocks additional potential for organisms to rapidly acclimate to their environment over short timeframes. If these changes affect fitness, they can circumvent the need for adaptive changes at the genome level. However, methylation has a complex reciprocal relationship with genetic variation as it can be genetically controlled, yet it can also induce point mutations and contribute to genomic evolution. When habitats remain constant over many generations, or populations are separated across habitats, initially plastic phenotypes can become hardwired through epigenetically facilitated mutagenesis. It remains unclear under what circumstances plasticity contributes to evolutionary outcomes, and when plastic changes will become permanently encoded into genotype. We highlight how studies investigating the evolution of epigenetic plasticity need to carefully consider how plasticity in methylation state could evolve among different evolutionary scenarios, the possible phenotypic outcomes, its effects on genomic evolution, and the proximate energetic and ultimate fitness costs of methylation. We argue that accumulating evidence suggests that DNA methylation can contribute toward evolution on various timescales, spanning a continuum from acclimatory plasticity to genomic evolution.

Epigenetic modification of the hypothalamic-pituitary-adrenal (HPA) axis during development in the house sparrow (Passer domesticus)

Epigenetic modifications such as DNA methylation are important mechanisms for mediating developmental plasticity, where ontogenetic processes and their phenotypic outcomes are shaped by early environments. In particular, changes in DNA methylation of genes within the hypothalamic-pituitary-adrenal (HPA) axis can impact offspring growth and development. This relationship has been well documented in mammals but is less understood in other taxa. Here, we use target-enriched enzymatic methyl sequencing (TEEM-seq) to assess how DNA methylation in a suite of 25 genes changes over development, how these modifications relate to the early environment, and how they predict differential growth trajectories in the house sparrow (Passer domesticus). We found that DNA methylation changes dynamically over the postnatal developmental period: genes with initially low DNA methylation tended to decline in methylation over development, whereas genes with initially high DNA methylation tended to increase in methylation. However, sex-specific differentially methylated regions (DMRs) were maintained across the developmental period. We also found significant differences in post-hatching DNA methylation in relation to hatch date, with higher levels of DNA methylation in nestlings hatched earlier in the season. Although these differences were largely absent by the end of development, a number of DMRs in HPA-related genes (CRH, MC2R, NR3C1, NR3C2, POMC)-and to a lesser degree HPG-related genes (GNRHR2)-predicted nestling growth trajectories over development. These findings provide insight into the mechanisms by which the early environment shapes DNA methylation in the HPA axis, and how these changes subsequently influence growth and potentially mediate developmental plasticity.

The Mutagenic Consequences of DNA Methylation within and across Generations

DNA methylation is an epigenetic modification with wide-ranging consequences across the life of an organism. This modification can be stable, persisting through development despite changing environmental conditions. However, in other contexts, DNA methylation can also be flexible, underlying organismal phenotypic plasticity. One underappreciated aspect of DNA methylation is that it is a potent mutagen; methylated cytosines mutate at a much faster rate than other genetic motifs. This mutagenic property of DNA methylation has been largely ignored in eco-evolutionary literature, despite its prevalence. Here, we explore how DNA methylation induced by environmental and other factors could promote mutation and lead to evolutionary change at a more rapid rate and in a more directed manner than through stochastic genetic mutations alone. We argue for future research on the evolutionary implications of DNA methylation driven mutations both within the lifetime of organisms, as well as across timescales.