Genetic and epigenetic divergence between disturbed and undisturbed subpopulations of a Mediterranean shrub: a 20-year field experiment

Comparative epigenetic and genetic spatial structure in Mediterranean mountain plants: a multispecies study

Changes in epigenetic states can allow individuals to cope with environmental changes. If such changes are heritable, this may lead to epigenetic adaptation. Thus, it is likely that in sessile organisms such as plants, part of the spatial epigenetic variation found across individuals will reflect the environmental heterogeneity within populations. The departure of the spatial epigenetic structure from the baseline genetic variation can help in understanding the value of epigenetic regulation in species with different breadth of optimal environmental requirements. Here, we hypothesise that in plants with narrow environmental requirements, epigenetic variability should be less structured in space given the lower variability in suitable environmental conditions. We performed a multispecies study that considered seven pairs of congeneric plant species, each encompassing a narrow endemic with habitat specialisation and a widespread species. In three populations per species we used AFLP and methylation-sensitive AFLP markers to characterise the spatial genetic and epigenetic structures. Narrow endemics showed a significantly lower epigenetic than genetic differentiation between populations. Within populations, epigenetic variation was less spatially structured than genetic variation, mainly in narrow endemics. In these species, structural equation models revealed that such pattern was associated to a lack of correlation between epigenetic and genetic information. Altogether, these results show a greater decoupling of the spatial epigenetic variation from the baseline spatial genetic pattern in endemic species. These findings highlight the value of studying genetic and epigenetic spatial variation to better understand habitat specialisation in plants.

The Contribution of Epigenetics to Evolutionary Adaptation in Zingiber kawagoii Hayata (Zingiberaceae) Endemic to Taiwan

We epigenotyped 211 individuals from 17 Zingiber kawagoii populations using methylation-sensitive amplification polymorphism (MSAP) and investigated the associations of methylated (mMSAP) and unmethylated (uMSAP) loci with 16 environmental variables. Data regarding genetic variation based on amplified fragment length polymorphism (AFLP) were obtained from an earlier study. We found a significant positive correlation between genetic and epigenetic variation. Significantly higher mean mMSAP and uMSAP uH_E (unbiased expected heterozygosity: 0.223 and 0.131, respectively, p < 0.001) per locus than that estimated based on AFLP (uH_E = 0.104) were found. Genome scans detected 10 mMSAP and 9 uMSAP F_ST outliers associated with various environmental variables. A significant linear fit for 11 and 12 environmental variables with outlier mMSAP and uMSAP ordination, respectively, generated using full model redundancy analysis (RDA) was found. When conditioned on geography, partial RDA revealed that five and six environmental variables, respectively, were the most important variables influencing outlier mMSAP and uMSAP variation. We found higher genetic (average F_ST = 0.298) than epigenetic (mMSAP and uMSAP average F_ST = 0.044 and 0.106, respectively) differentiation and higher genetic isolation-by-distance (IBD) than epigenetic IBD. Strong epigenetic isolation-by-environment (IBE) was found, particularly based on the outlier data, controlling either for geography (mMSAP and uMSAP β_E = 0.128 and 0.132, respectively, p = 0.001) or for genetic structure (mMSAP and uMSAP β_E = 0.105 and 0.136, respectively, p = 0.001). Our results suggest that epigenetic variants can be substrates for natural selection linked to environmental variables and complement genetic changes in the adaptive evolution of Z. kawagoii populations.

The impact of environment on genetic and epigenetic variation in Trifolium pratense populations from two contrasting semi-natural grasslands

Central European grasslands, such as calcareous grasslands and oat-grass meadows, are characterized by diverse environmental conditions and management regimes. Therefore, we aimed to determine potential differences in genetic and epigenetic variation patterns between the contrasting habitats and to identify the drivers of genetic and epigenetic variation. We investigated the genetic and epigenetic variation of the ecologically variable plant species Trifolium pratense L. applying amplified fragment length polymorphism and methylation-sensitive amplification polymorphism analyses. We observed low levels of genetic and epigenetic differentiation among populations and between habitat types. Genetic and epigenetic variations were not interdependent. Thus, genetic variation was significantly isolated by habitat dissimilarity, whereas epigenetic variation was affected by environment. More specifically, we observed a significant correlation of epigenetic diversity with soil moisture and soil pH (the latter potentially resulting in phosphorus limitation). Genetic variation was, therefore, affected more strongly by habitat-specific environmental conditions induced by land use-related disturbance and gene flow patterns, while epigenetic variation was driven by challenging environmental conditions.

Lifetime genealogical divergence within plants leads to epigenetic mosaicism in the shrub Lavandula latifolia (Lamiaceae)

Epigenetic mosaicism is a possible source of within-plant phenotypic heterogeneity, yet its frequency and developmental origin remain unexplored. This study examines whether extant epigenetic heterogeneity within Lavandula latifolia (Lamiaceae) shrubs reflects recent epigenetic modifications experienced independently by different plant parts or, alternatively, it is the cumulative outcome of a steady lifetime process. Leaf samples from different architectural modules (branch tips) were collected from three L. latifolia plants and characterized epigenetically by global DNA cytosine methylation and methylation state of methylation-sensitive amplified fragment-length polymorphism (MS-AFLP) markers. Epigenetic characteristics of modules were then assembled with information on the branching history of plants. Methods borrowed from phylogenetic research were used to assess genealogical signal of extant epigenetic variation and reconstruct within-plant genealogical trajectory of epigenetic traits. Plants were epigenetically heterogeneous, as shown by differences among modules in global DNA methylation and variation in the methylation states of 6 to 8% of MS-AFLP markers. All epigenetic features exhibited significant genealogical signal within plants. Events of epigenetic divergence occurred throughout the lifespan of individuals and were subsequently propagated by branch divisions. Internal epigenetic diversification of L. latifolia individuals took place steadily during their development, a process which eventually led to persistent epigenetic mosaicism.

Metabolomic Variation Aligns with Two Geographically Distinct Subpopulations of Brachypodium Distachyon before and after Drought Stress

Brachypodium distachyon (Brachypodium) is a non-domesticated model grass that has been used to assess population level genomic variation. We have previously established a collection of 55 Brachypodium accessions that were sampled to reflect five different climatic regions of Turkey; designated 1a, 1c, 2, 3 and 4. Genomic and methylomic variation differentiated the collection into two subpopulations designated as coastal and central (respectively from regions 1a, 1c and the other from 2, 3 and 4) which were linked to environmental variables such as relative precipitation. Here, we assessed how far genomic variation would be reflected in the metabolomes and if this could be linked to an adaptive trait. Metabolites were extracted from eight-week-old seedlings from each accession and assessed using flow infusion high-resolution mass spectrometry (FIE-HRMS). Principal Component Analysis (PCA) of the derived metabolomes differentiated between samples from coastal and central subpopulations. The major sources of variation between seedling from the coastal and central subpopulations were identified. The central subpopulation was typified by significant increases in alanine, aspartate and glutamate metabolism and the tricarboxylic acid (TCA) cycle. Coastal subpopulation exhibited elevated levels of the auxin, indolacetic acid and rhamnose. The metabolomes of the seedling were also determined following the imposition of drought stress for seven days. The central subpopulation exhibited a metabolomic shift in response to drought, but no significant changes were seen in the coastal one. The drought responses in the central subpopulation were typified by changes in amino acids, increasing the glutamine that could be functioning as a stress signal. There were also changes in sugars that were likely to be an osmotic counter to drought, and changes in bioenergetic metabolism. These data indicate that genomic variation in our Turkish Brachypodium collection is largely reflected as distinctive metabolomes (“metabolotypes”) through which drought tolerance might be mediated.

Sensitivity analysis based on the random forest machine learning algorithm identifies candidate genes for regulation of innate and adaptive immune response of chicken

Two categories of immune responses—innate and adaptive immunity—have both polygenic backgrounds and a significant environmental component. The goal of the reported study was to define candidate genes and mutations for the immune traits of interest in chickens using machine learning–based sensitivity analysis for single-nucleotide polymorphisms (SNPs) located in candidate genes defined in quantitative trait loci regions. Here the adaptive immunity is represented by the specific antibody response toward keyhole limpet hemocyanin (KLH), whereas the innate immunity was represented by natural antibodies toward lipopolysaccharide (LPS) and lipoteichoic acid (LTA). The analysis consisted of 3 basic steps: an identification of candidate SNPs via feature selection, an optimisation of the feature set using recursive feature elimination, and finally a gene-level sensitivity analysis for final selection of models. The predictive model based on 5 genes (MAPK8IP3 CRLF3, UNC13D, ILR9, and PRCKB) explains 14.9% of variance for KLH adaptive response. The models obtained for LTA and LPS use more genes and have lower predictive power, explaining respectively 7.8 and 4.5% of total variance. In comparison, the linear models built on genes identified by a standard statistical analysis explain 1.5, 0.5, and 0.3% of variance for KLH, LTA, and LPS response, respectively. The present study shows that machine learning methods applied to systems with a complex interaction network can discover phenotype-genotype associations with much higher sensitivity than traditional statistical models. It adds contribution to evidence suggesting a role of MAPK8IP3 in the adaptive immune response. It also indicates that CRLF3 is involved in this process as well. Both findings need additional verification.

Manipulation of cytosine methylation does not remove latitudinal clines in two invasive goldenrod species in Central Europe

Invasive species frequently differentiate phenotypically in novel environments within a few generations, often even with limited genetic variation. For the invasive plants Solidago canadensis and S. gigantea, we tested whether such differentiation might have occurred through heritable epigenetic changes in cytosine methylation. In a 2-year common-garden experiment, we grew plants from seeds collected along a latitudinal gradient in their non-native Central European range to test for trait differentiation and whether differentiation disappeared when seeds were treated with the demethylation agent zebularine. Microsatellite markers revealed no population structure along the latitudinal gradient in S. canadensis, but three genetic clusters in S. gigantea. Solidago canadensis showed latitudinal clines in flowering phenology and growth. In S. gigantea, the number of clonal offspring decreased with latitude. Although zebularine had a significant effect on early growth, probably through effects on cytosine methylation, latitudinal clines remained (or even got stronger) in plants raised from seeds treated with zebularine. Thus, our experiment provides no evidence that epigenetic mechanisms by selective cytosine methylation contribute to the observed phenotypic differentiation in invasive goldenrods in Central Europe.

Genetic and Methylome Variation in Turkish Brachypodium Distachyon Accessions Differentiate Two Geographically Distinct Subpopulations

Brachypodium distachyon (Brachypodium) is a non-domesticated model grass species that can be used to test if variation in genetic sequence or methylation are linked to environmental differences. To assess this, we collected seeds from 12 sites within five climatically distinct regions of Turkey. Seeds from each region were grown under standardized growth conditions in the UK to preserve methylated sequence variation. At six weeks following germination, leaves were sampled and assessed for genomic and DNA methylation variation. In a follow-up experiment, phenomic approaches were used to describe plant growth and drought responses. Genome sequencing and population structure analysis suggested three ancestral clusters across the Mediterranean, two of which were geographically separated in Turkey into coastal and central subpopulations. Phenotypic analyses showed that the coastal subpopulation tended to exhibit relatively delayed flowering and the central, increased drought tolerance as indicated by reduced yellowing. Genome-wide methylation analyses in GpC, CHG and CHH contexts also showed variation which aligned with the separation into coastal and central subpopulations. The climate niche modelling of both subpopulations showed a significant influence from the “Precipitation in the Driest Quarter” on the central subpopulation and “Temperature of the Coldest Month” on the coastal subpopulation. Our work demonstrates genetic diversity and variation in DNA methylation in Turkish accessions of Brachypodium that may be associated with climate variables and the molecular basis of which will feature in ongoing analyses.

Small spaces, big impacts: contributions of micro-environmental variation to population persistence under climate change

Individuals within natural populations can experience very different abiotic and biotic conditions across small spatial scales owing to microtopography and other micro-environmental gradients. Ecological and evolutionary studies often ignore the effects of micro-environment on plant population and community dynamics. Here, we explore the extent to which fine-grained variation in abiotic and biotic conditions contributes to within-population variation in trait expression and genetic diversity in natural plant populations. Furthermore, we consider whether benign microhabitats could buffer local populations of some plant species from abiotic stresses imposed by rapid anthropogenic climate change. If microrefugia sustain local populations and communities in the short term, other eco-evolutionary processes, such as gene flow and adaptation, could enhance population stability in the longer term. We caution, however, that local populations may still decline in size as they contract into rare microhabitats and microrefugia. We encourage future research that explicitly examines the role of the micro-environment in maintaining genetic variation within local populations, favouring the evolution of phenotypic plasticity at local scales and enhancing population persistence under global change.

No evidence for local adaptation and an epigenetic underpinning in native and non-native ruderal plant species in Germany

Many invasive species have rapidly adapted to different environments in their new ranges. This is surprising, as colonization is usually associated with reduced genetic variation. Heritable phenotypic variation with an epigenetic basis may explain this paradox.Here, we assessed the contribution of DNA methylation to local adaptation in native and naturalized non-native ruderal plant species in Germany. We reciprocally transplanted offspring from natural populations of seven native and five non-native plant species between the Konstanz region in the south and the Potsdam region in the north of Germany. Before the transplant, half of the seeds were treated with the demethylation agent zebularine. We recorded survival, flowering probability, and biomass production as fitness estimates.Contrary to our expectations, we found little evidence for local adaptation, both among the native and among the non-native plant species. Zebularine treatment had mostly negative effects on overall plant performance, regardless of whether plants were local or not, and regardless of whether they were native or non-native. Synthesis. We conclude that local adaptation, at least at the scale of our study, plays no major role in the success of non-native and native ruderal plants. Consequently, we found no evidence yet for an epigenetic basis of local adaptation.

Epigenetic Differentiation of Natural Populations of Lilium bosniacum Associated with Contrasting Habitat Conditions

Epigenetic variation in natural populations with contrasting habitats might be an important element, in addition to the genetic variation, in plant adaptation to environmental stress. Here, we assessed genetic, epigenetic, and cytogenetic structure of the three Lilium bosniacum populations growing on distinct habitats. One population was growing under habitual ecological conditions for this species and the other two were growing under stress associated with high altitude and serpentine soil. Amplified fragment length polymorphism and methylation-sensitive amplification polymorphism analyses revealed that the three populations did not differentiate genetically, but were clearly separated in three distinct clusters according to DNA methylation profiles. Principal coordinate analysis showed that overall epigenetic variation was closely related to habitat conditions. A new methylation-sensitive amplification polymorphism scoring approach allowed identification of mainly unmethylated (φST = 0.190) and fully CpG methylated (φST = 0.118) subepiloci playing a role in overall population differentiation, in comparison with hemimethylated sites (φST = 0.073). In addition, unusual rDNA repatterning and the presence of B chromosomes bearing 5S rDNA loci were recorded in the population growing on serpentine soil, suggesting dynamic chromosome rearrangements probably linked to global genome demethylation, which might have reactivated some mobile elements. We discuss our results considering our earlier data on morphology and leaf anatomy of several L. bosniacum populations, and suggest a possible role of epigenetics as a key element in population differentiation associated with environmental stress in these particular lily populations.

Genetic and epigenetic changes during the invasion of a cosmopolitan species (Phragmites australis)

While many introduced invasive species can increase genetic diversity through multiple introductions and/or hybridization to colonize successfully in new environments, others with low genetic diversity have to persist by alternative mechanisms such as epigenetic variation. Given that Phragmites australis is a cosmopolitan reed growing in a wide range of habitats and its invasion history, especially in North America, has been relatively well studied, it provides an ideal system for studying the role and relationship of genetic and epigenetic variation in biological invasions. We used amplified fragment length polymorphism (AFLP) and methylation-sensitive (MS) AFLP methods to evaluate genetic and epigenetic diversity and structure in groups of the common reed across its range in the world. Evidence from analysis of molecular variance (AMOVA) based on AFLP and MS-AFLP data supported the previous conclusion that the invasive introduced populations of P. australis in North America were from European and Mediterranean regions. In the Gulf Coast region, the introduced group harbored a high level of genetic variation relative to originating group from its native location, and it showed epigenetic diversity equal to that of the native group, if not higher, while the introduced group held lower genetic diversity than the native. In the Great Lakes region, the native group displayed very low genetic and epigenetic variation, and the introduced one showed slightly lower genetic and epigenetic diversity than the original one. Unexpectedly, AMOVA and principal component analysis did not demonstrate any epigenetic convergence between native and introduced groups before genetic convergence. Our results suggested that intertwined changes in genetic and epigenetic variation were involved in the invasion success in North America. Although our study did not provide strong evidence proving the importance of epigenetic variation prior to genetic, it implied the similar role of stable epigenetic diversity to genetic diversity in the adaptation of P. australis to local environment.

Transgenerational epigenetics: Inheritance of global cytosine methylation and methylation-related epigenetic markers in the shrub Lavandula latifolia

PREMISE OF THE STUDY

The ecological and evolutionary significance of natural epigenetic variation (i.e., not based on DNA sequence variants) variation will depend critically on whether epigenetic states are transmitted from parents to offspring, but little is known on epigenetic inheritance in nonmodel plants.

METHODS

We present a quantitative analysis of transgenerational transmission of global DNA cytosine methylation (= proportion of all genomic cytosines that are methylated) and individual epigenetic markers (= methylation status of anonymous MSAP markers) in the shrub Lavandula latifolia. Methods based on parent-offspring correlations and parental variance component estimation were applied to epigenetic features of field-growing plants ('maternal parents') and greenhouse-grown progenies. Transmission of genetic markers (AFLP) was also assessed for reference.

KEY RESULTS

Maternal parents differed significantly in global DNA cytosine methylation (range = 21.7-36.7%). Greenhouse-grown maternal families differed significantly in global methylation, and their differences were significantly related to maternal origin. Methylation-sensitive amplified polymorphism (MSAP) markers exhibited significant transgenerational transmission, as denoted by significant maternal variance component of marker scores in greenhouse families and significant mother-offspring correlations of marker scores.

CONCLUSIONS

Although transmission-related measurements for global methylation and MSAP markers were quantitatively lower than those for AFLP markers taken as reference, this study has revealed extensive transgenerational transmission of genome-wide global cytosine methylation and anonymous epigenetic markers in L. latifolia. Similarity of results for global cytosine methylation and epigenetic markers lends robustness to this conclusion, and stresses the value of considering both types of information in epigenetic studies of nonmodel plants.

Genetic and epigenetic variations associated with adaptation to heterogeneous habitat conditions in a deciduous shrub

Environmentally induced phenotypic plasticity is thought to play an important role in the adaption of plant populations to heterogeneous habitat conditions, and yet the importance of epigenetic variation as a mechanism of adaptive plasticity in natural plant populations still merits further research. In this study, we investigated populations of Vitex negundo var. heterophylla (Chinese chastetree) from adjacent habitat types at seven sampling sites. Using several functional traits, we detected a significant differentiation between habitat types. With amplified fragment length polymorphisms (AFLP) and methylation-sensitive AFLP (MSAP), we found relatively high levels of genetic and epigenetic diversity but very low genetic and epigenetic differences between habitats within sites. Bayesian clustering showed a remarkable habitat-related differentiation and more genetic loci associated with the habitat type than epigenetic, suggesting that the adaptation to the habitat is genetically based. However, we did not find any significant correlation between genetic or epigenetic variation and habitat using simple and partial Mantel tests. Moreover, we found no correlation between genetic and ecologically relevant phenotypic variation and a significant correlation between epigenetic and phenotypic variation. Although we did not find any direct relationship between epigenetic variation and habitat environment, our findings suggest that epigenetic variation may complement genetic variation as a source of functional phenotypic diversity associated with adaptation to the heterogeneous habitat in natural plant populations.

Transgenerational stress-adaption: an opportunity for ecological epigenetics

In the recent years, there has been considerable interest to investigate the adaptive transgenerational plasticity of plants and how a "stress memory" can be transmitted to the following generation. Although, increasing evidence suggests that transgenerational adaptive responses have widespread ecological relevance, the underlying epigenetic processes have rarely been elucidated. On the other hand, model plant species have been deeply investigated in their genome-wide methylation landscape without connecting this to the ecological reality of the plant. What we need is the combination of an ecological understanding which plant species would benefit from transgenerational epigenetic stress-adaption in their natural habitat, combined with a deeper molecular analysis of non-model organisms. Only such interdisciplinary linkage in an ecological epigenetic study could unravel the full potential that epigenetics could play for the transgenerational stress-adaption of plants.

Uncovering the contribution of epigenetics to plant phenotypic variation in Mediterranean ecosystems

Epigenetic signals can affect plant phenotype and fitness and be stably inherited across multiple generations. Epigenetic regulation plays a key role in the mechanisms of plant response to the environment, without altering DNA sequence. As plants cannot adapt behaviourally or migrate instantly, such dynamic epigenetic responses may be particularly crucial for survival of plants within changing and challenging environments, such as the Mediterranean-Type Ecosystems (MTEs). These ecosystems suffer recurrent stressful events (warm and dry summers with associated fire regimes) that have selected for plants with similar phenotypic complex traits, resulting in similar vegetation growth forms. However, the potential role of epigenetics in plant adaptation to recurrent stressful environments such as the MTEs has generally been ignored. To understand the full spectrum of adaptive processes in such contexts, it is imperative to prompt study of the causes and consequences of epigenetic variation in natural populations. With this purpose, we review here current knowledge on epigenetic variation in natural populations and the genetic and epigenetic basis of some key traits for plants in the MTEs, namely those traits involved in adaptation to drought, fire and oligotrophic soils. We conclude there is still much to be learned about 'plant epigenetics in the wild' and, thus, we propose future research steps in the study of natural epigenetic variation of key traits in the MTEs at different scales.

Comparative epigenetic and genetic spatial structure of the perennial herb Helleborus foetidus: Isolation by environment, isolation by distance, and functional trait divergence

PREMISE OF THE STUDY

Epigenetic variation can play a role in local adaptation; thus, there should be associations among epigenetic variation, environmental variation, and functional trait variation across populations. This study examines these relationships in the perennial herb Helleborus foetidus (Ranunculaceae).

METHODS

Plants from 10 subpopulations were characterized genetically (AFLP, SSR markers), epigenetically (MSAP markers), and phenotypically (20 functional traits). Habitats were characterized using six environmental variables. Isolation-by-distance (IBD) and isolation-by-environment (IBE) patterns of genetic and epigenetic divergence were assessed, as was the comparative explanatory value of geographical and environmental distance as predictors of epigenetic, genetic, and functional differentiation.

KEY RESULTS

Subpopulations were differentiated genetically, epigenetically, and phenotypically. Genetic differentiation was best explained by geographical distance, while epigenetic differentiation was best explained by environmental distance. Divergence in functional traits was correlated with environmental and epigenetic distances, but not with geographical and genetic distances.

CONCLUSIONS

Results are compatible with the hypothesis that epigenetic IBE and functional divergence reflected responses to environmental variation. Spatial analyses simultaneously considering epigenetic, genetic, phenotypic and environmental information provide a useful tool to evaluate the role of environmental features as drivers of natural epigenetic variation between populations.