Mo-Zhu Wang
Hong-Li Li
Min Tang1
Fei-Hai Yu- 1Institute of Wetland Ecology & Clone Ecology/Zhejiang Provincial Key Laboratory of Plant Evolutionary Ecology and Conservation, Taizhou University, Taizhou, China
- 2State Key Laboratory of Systematic and Evolutionary Botany, Institute of Botany, Chinese Academy of Sciences, Beijing, China
- 3School of Ecology and Nature Conservation, Beijing Forestry University, Beijing, China
Epigenetic mechanisms such as DNA methylation are considered as an important pathway responsible for phenotypic responses and rapid acclimation of plants to different environments. To search for empirical evidence that DNA methylation is implicated in stress-responses of non-model species, we exposed genetically uniform, experimental populations of the wetland clonal plant Hydrocotyle vulgaris to two manipulated flood regimes, i.e., semi-submergence vs. submergence, measured phenotypic traits, and quantified different types of DNA methylation using MSAP (methylation-sensitive amplified polymorphism). We found different epi-phenotypes and significant epigenetic differentiation between semi-submerged and submerged populations. Compared to subepiloci (denoting DNA methylation conditions) for the CG-methylated state, unmethylation and CHG-hemimethylation subepiloci types contribute more prominently to the epigenetic structure of experimental populations. Moreover, we detected some epimarker outliers potentially facilitate population divergence between two flood regimes. Some phenotypic variation was associated with flood-induced DNA methylation variation through different types of subepiloci. Our study provides the indication that DNA methylation might be involved in plant responses to environmental variation without altering DNA sequences.
Introduction
Plants exposed to environmental changes often exhibit plastic phenotypes (Putnam et al., 2016; Colicchio et al., 2018). The classical view advocates that an individual phenotype is determined by its environment, genotype, and their interaction, including both plasticity and evolutionary adaptation (Richards et al., 2010; Bian et al., 2013). However, epigenetic regulation (e.g., DNA methylation, histone modifications, chromatin remodeling, and expression of non-coding RNAs) without changing DNA sequence has been widely considered as another candidate mechanism accounting for plant phenotypic variation (Bossdorf et al., 2008; Marfil et al., 2009). Nowadays, the best-studied hallmark of epigenetic modification is DNA methylation, which is mainly the addition of a methyl group to the C5 position of a cytosine residue in three different sequence contexts (CG, CHG, and CHH sites, H = A, C, T), via catalysis of several DNA methyltransferase enzymes (Bossdorf et al., 2008; Zoldoš et al., 2018).
Alteration of epigenetic markers mainly originate from genetic variation, environmental induction, or spontaneous epimutations (Dubin et al., 2015; Trucchi et al., 2016; Richards et al., 2017). According to dependence degree on genetic context, epigenetic variation could be classified into three categories: obligate (fully dependent), facilitated (semi-independent), or pure (completely independent) (Richards, 2006; Robertson and Richards, 2015). Unlike genetic variation, DNA methylation patterns are sensitive to changing environments and commonly possess a much higher variation rate (Schulz et al., 2014; Jueterbock et al., 2020). Such variation can be reversibly transient within one generation or stably heritable to several generations (Angers et al., 2010; Paun et al., 2010; Colicchio et al., 2018). DNA methylation alters gene expression through transcriptional repression or remodeling chromatin, further affecting plant phenotypes (Boyko and Kovalchuk, 2010; Grativol et al., 2012; Griffin et al., 2016; Colicchio and Herman, 2020). Therefore, environment-induced epigenetic regulation could not only offer a rapid pathway for phenotypic plasticity, but also underlie plant adaptive evolution when across-generational plasticity confers fitness benefits in predictable environments (Putnam et al., 2016; Huang et al., 2017; Groot et al., 2018; Colicchio and Herman, 2020).
Depending on sequence context, different DNA methylation types (i.e., CG/CHG/CHH) vary in their responses to environmental factors, associations with genetic variation, or functions in gene expression, etc. (Schulz et al., 2014; Colicchio et al., 2018). For instance, in Arabidopsis thaliana accessions, CHH methylation of transposable elements (TEs) was sensitive to growth temperature and under cis- and trans-acting genetic control, whereas CG methylation on the gene coding regions was independent of genetic effects and instead strongly correlated with the latitude of origin (Dubin et al., 2015). In general, CG methylation in gene bodies (GbM) usually activates gene expression, while that in promoters and TEs is associated with gene silencing (Dubin et al., 2015; Colicchio et al., 2018; Jueterbock et al., 2020). Non-CG methylation (i.e., CHH or CHG methylation) mostly occurring in transposons or repeat regions seems to regulate transcriptional repression through chromatin remodeling (Grativol et al., 2012; Schulz et al., 2013; Dubin et al., 2015; Colicchio et al., 2018).
To explore ecological and evolutionary significance of DNA methylation, the first step is to find evidence that at least part of the epigenome changes correlate with plant stress-responses. In recent years, there are accumulating ecological studies exploring different DNA methylation types responding to environmental stresses and/or their relations to plant phenotypic characteristics. Most of these studies focused on plants with the genome reference at the individual level (e.g., Colicchio et al., 2018; Jueterbock et al., 2020). However, roles of epigenetic variation in plastic responses of natural populations to specific environment changes is still largely unknown, especially for non-traditional model organisms lacking the genome reference (Paun et al., 2010; Abratowska et al., 2012; Rico et al., 2014; Watson et al., 2018).
Given the sensitivity of DNA methylation to environmental variation, a direct test for the epigenetic contribution in natural systems could be confounded by complex and dynamic natural conditions (Schulz et al., 2014). Moreover, some previous studies used genetically diverse plant materials, which might hardly disclose the pure epigenetic effects with the presence of genetic variation (Lira-Medeiros et al., 2010; Schulz et al., 2014; Robertson et al., 2017). Therefore, manipulated experimental populations without any genetic variation, such as those consisting of genetically identical asexual individuals (ramets) vegetatively propagated by a single genet (clone) of clonal plants, are better materials to strictly assess roles of epigenetic variation in plastic responses (Rapp and Wendel, 2005; Verhoeven et al., 2010; Verhoeven and Preite, 2014; Huang et al., 2017).
The wetland clonal plant Hydrocotyle vulgaris L. (Araliaceae) is considered potentially invasive in China due to high phenotypic plasticity, rapid clonal growth, strong adaptability, and exclusion of other native species (Miao et al., 2011; Liu et al., 2014; Dong et al., 2015). Our previous study showed that the natural H. vulgaris populations in southern China possessed low genetic variation but high epigenetic variation, and that their phenotypic variation was largely correlated with epigenetic variation rather than genetic variation (Wang et al., 2020). H. vulgaris often experiences water depth changes, which may represent a strong selective force for its population diversification (Wang et al., 2020). In this study, we explored roles of different DNA methylation types in phenotypic responses of H. vulgaris to flood variation by exposing its genetically uniform experimental populations to two different flood regimes and by evaluating phenotypic and DNA methylation consequences using MSAP (methylation-sensitive amplified polymorphism). Specifically, we addressed the following questions. (1) What are the phenotypic responses of experimental H. vulgaris populations to different flood regimes? (2) What DNA methylation patterns can be generated in different flood regimes? (3) Are environmentally induced alterations in different methylation types related to phenotypic variation?
Materials and Methods
Material Propagation
From June to August 2016, 128 plants of H. vulgaris were collected from 10 natural populations in southern China (Wang et al., 2020). Using AFLP (amplified fragment length polymorphism), we distinguished 20 genotypes from the 128 individuals, among which a single wide spread genotype accounted for 82% of the total samples and dominated in all 10 populations (Wang et al., 2020). Plants of the most dominant genotype were mixed cultivated and vegetatively propagated under the same condition in a greenhouse at Taizhou University. In early July 2017, we selected more than 576 newly generated similar-sized ramets at the same developmental stage. Each ramet consisted of one node, one leaf and some adventitious roots (petiole length: 22.5 ± 0.2 cm, mean ± SE, n = 30). To ensure that all ramets are epigenetically uniform at the start of the experiment, thirty ramets of them were randomly selected for detection of DNA methylation patterns by MSAP, and were identified to be assigned to the same epigenotype.
Experimental Design
The remaining 546 ramets were used for the experiment described below. We constructed experimental populations of H. vulgaris in six big plastic tanks (1.38 m in bottom diameter, 1.60 m in top diameter and 0.89 m in height) filled with a 30-cm-deep mixture of sand and local soil at a 1:1 volume ratio. In each tank, the 91 similar-sized, genetically identical ramets were evenly planted in the range of a 50-cm-edged hexagon from the center point of the soil surface, with two adjacent ramets spacing 10 cm apart (Supplementary Appendix 1A). The soil in the tanks was always kept moist after planting. After 20 days recovery of the six established experimental populations, two flooding treatments were applied. No ramet died before the flooding treatments.
The two flooding treatments were semi-submergence and submergence, each with three replicate tanks (experimental populations). For the semi-submergence treatment, the tank was filled with tap water to a depth of 10 cm above the soil surface, so that the ramets could protrude from the water surface. By contrast, for the submergence treatment, the water level in the tank was maintained 30 cm above the soil surface, so that the ramets were submerged (under the water surface). The experiment was conducted in an open area at Taizhou University, and all the six tanks were placed closely and randomly to avoid potential confounding effects of micro-environmental differences.
Harvest and Measurements
The experiment lasted from 10 September to 20 December 2017. At harvest, we uniformly set 19 sampling points in each tank, similar to the planting approach, and the difference was that the sampling points were at 25 cm intervals (Supplementary Appendix 1B). At each sampling point, we took two connected mature ramets with fully expanded leaves: one ramet was randomly selected for phenotypic measurement and the other for epigenetic analysis.
For epigenetic analysis, the leaf of the ramet was dried in silica gel. Total genomic DNA from 30 mg of the dry leaf was extracted using Dingguo Plant Genomic DNA Kit (Beijing, China), and quantified spectrophotometrically. After verifying integrity and purity by 1% agarose gel electrophoresis, DNA was diluted to 20 ng/μL as the starting material for epigenetic analysis. The MSAP protocol and scoring method were exactly the same as our previous study, with five selective primer combinations, i.e., E-AGT/H-TAT, E-AGT/H-TTC, E-ATC/H-TGA, E-AAC/H-TCG, and E-ATG/H-TGA (Wang et al., 2020). The error rates for HpaII and MspI scores were about 0.65 and 0.46%, respectively (Wang et al., 2020). Only the repeatable markers were involved in the following molecular analyses.
To quantify phenotypic responses, we first measured leaf petiole length, leaf area and stem internode length of the ramet. Then, the petiole, leaf blade and stem internode were dried at 90°C for 48 h and weighed. Specific petiole length was calculated as petiole length per unit petiole dry mass, specific leaf area as leaf area per unit leaf dry mass, and specific internode length as internode length per unit internode dry mass.
Data Analysis
We used nested ANOVA to test the effect of flooding treatments on each of the six phenotypic traits (petiole length, specific petiole length, leaf area, specific leaf area, internode length, and specific internode length) at the population level. Experimental populations were nested within the treatment. Before analyses, specific petiole length, internode length and leaf area were log-transformed to improve homoscedasticity (Supplementary Appendix 2).
For MSAP data, the presence or absence of the bands from specific isoschizomer digestions (EcoRI/HpaII and EcoRI/MspI) results in four conditions of a particular fragment: (I) bands present in both enzyme combinations (1/1), indicating an unmethylated state; (II) bands absent in both enzyme combinations (0/0), indicating an uninformative state; (III) bands present only in EcoRI/MspI profiles (0/1), indicating hemi- or fully methylated CG-sites; (IV) bands present only in EcoRI/HpaII profiles (1/0), indicating hemimethylated CHG-sites. Due to the fact that different methylation states participate in different regulating processes, considering them separately would give the most comprehensive picture of DNA methylation (Schulz et al., 2013). Therefore, we used the “Mixed-Scoring 2” approach implemented in R script “MSAP_calc.r” (Schulz et al., 2013) to transform the three discernible methylation status represented by combination of EcoRI/HpaII and EcoRI/MspI banding patterns at each epilocus into binary matrices of different types of subepiloci. Thus, for each epilocus, up to three subepiloci can be generated: u-subepilocus (denoting the unmethylated loci where type I is scored as 1 and other types were scored as 0), m-subepilocus (denoting the CG-methylated loci where type III is scored as 1), and h-subepilocus (denoting the CHG-hemimethylated loci where type IV is scored as 1) (Schulz et al., 2014).
Based on the binary matrices of all MSAP subepiloci and each subepiloci type, we conducted the following analyses. Epigenetic diversity of each experimental population in terms of the percentage of polymorphic loci (PLP) and Shannon’s information index (H) were assessed by “MSAP_calc.r” (Schulz et al., 2013). To visualize population epigenetic structure, principal coordinate analyses (PCoA) were performed with GenALEx 6.5 based on the matrix of Nei’s distances (Peakall and Smouse, 2012). Epigenetic differentiation at different hierarchical components, that is, between treatments (ϕRT), among experimental populations within treatments (ϕPR) and within experimental populations (ϕPT), was calculated using analysis of molecular variance (AMOVA) with Genalex 6.5. Significance levels were determined after 9,999 permutations.
To identify putatively adaptive epiloci that may facilitate shaping population epigenetic responses to flood variation, we performed outlier detection based on individuals from different experimental populations by using the BayeScan 2.1 (Foll and Gaggiotti, 2008) for estimating the posterior odds (PO) of each epilocus. The analyses were run for 100,000 iterations, with a burn-in of 50,000 iterations, a sample size of 5,000 and a thinning interval of 10. An additional burn-in was carried out by 20 short pilot runs of 5,000 iterations. Only loci exceeding a “strong” detection level [log10 (PO) > 1] were considered as putative outliers.
To establish the relationships between environmental, epigenetic and phenotypic variation, structural equation modeling (SEM) was conducted in AMOS 24.0, by relating flood regime and epigenetic variation on phenotypic variation. For each subepiloci type, we examined the direct effects of methylation variation (first three PCoA axis for corresponding subepiloci) and the environmental factor (two flood regimes; semi-submergence was coded as “0,” while submergence treatment was coded as “1”) on phenotypic variation (six phenotypic traits), and indirect effects of flood regimes on phenotypic traits through methylation variation.
Results
Phenotypic Responses to Different Flood Regimes
Flooding significantly affected phenotypic traits of H. vulgaris (Figure 1). Submerged populations exhibited significantly shorter petiole length, internode length and smaller leaf area, but higher specific petiole length, specific internode length and specific leaf area than semi-submerged populations (Figure 1).
Figure 1. Phenotypic traits of the Hydrocotyle vulgaris at population level under semi-submergence (Semi.) and submergence (Sub.). Phenotypic traits are (A) petiole length, (B) leaf area, (C) internode length, (D) specific petiole length, (E) specific leaf are, and (F) specific internode length. Mean ± SE are shown (n = 3). F, P, and degree of freedom of nested ANOVAs are also given.
Epigenetic Responses to Different Flood Regimes
Epigenetic Diversity
The MSAP analysis resulted in 345 scorable epiloci, of which 75 (21.7%) were polymorphic. Mixed scoring 2 detected 144 polymorphic subepiloci, including 66 u-, 42 m-, and 36 h-subepiloci. Epigenetic diversity of submerged populations was significantly higher than that of semi-submerged populations, as quantified by percentage of polymorphic loci and Shannon’s information index based on all subepiloci, u-subepiloci and m-subepiloci (Table 1 and Supplementary Appendix 3).
Table 1. Epigenetic diversity of the six clonally propagated populations of Hydrocotyle vulgaris under semi-submergence (Semi.) and submergence (Sub.) as quantified by (A) percentage of polymorphic loci and (B) Shannon’s information index based on all subepiloci and u-, m-, and h-subepiloci of MSAP.
Epigenetic Structure
Principal coordinates analysis revealed that the epigenetic structure of experimental populations differed among different types of subepiloci (Figure 2). Based on all MSAP subepiloci, epigenetic distances separated semi-submergence and submergence treatments along the first axis, forming two separated clusters (Figure 2A). Moreover, the experimental populations 1 and 3 in the submergence treatment fell apart, whereas the experimental populations in the semi-submergence treatment were much closer, with higher convergence degree. A similar differentiation pattern was also found in u-subepiloci, with individuals more clumped in each treatment (Figure 2B). For m-subepiloci, most individuals from different flood regimes grouped together without clear population differentiation, and only some individuals of the submerged experimental population 2 scattered from the cluster (Figure 2C). For h-subepiloci, the two treatments were mainly separated along the second coordinate, forming two clusters, with several individuals separated from the group of the submerged populations (Figure 2D).
Figure 2. Results of principal coordinate analyses (PCoA) based on (A) all subepiloci, and (B) u-, (C) m-, and (D) h-subepiloci of MSAP. Semi. and Sub. stand for the semi-submergence and submergence treatment, respectively.
For the combined epigenetic dataset (Table 2), AMOVA showed that 36% of epigenetic variance occurred between treatments, 50% within experimental populations, and only 14% among experimental populations within treatments. Similarly, for u-subepiloci and h-subepiloci, most variance occurred between treatments and within experimental populations (for u-subepiloci, 46 and 41%, respectively; for h-subepiloci, 39 and 55%, respectively). However, for m-subepiloci, variation mainly existed within experimental populations (73%).
Table 2. Results of hierarchical AMOVA based on (A) all subepiloci and (B) u-, (C) m-, and (D) h-subepiloci of MSAP.
Outlier Detection
For the complete set of the 144 MSAP subepiloci, BayeScan identified 10 (6.9%) outliers (Figure 3), among which seven were u-subepiloci, one was m-subepiloci and two were h-subepiloci, accounting for 10.61, 2.38, and 5.56% of the corresponding type of outliers, respectively. Based on PCoA analysis, outliers clearly separated the semi-submergence and submergence populations along the first axis, while there was no clear differentiation between the two treatments for neutral subepiloci (Supplementary Appendix 4).
Figure 3. MSAP outliers identified by Bayescan. The posterior odds (PO) for a locus to be under divergent selection are shown on a log scale on the x-axis. The vertical dash line indicates the threshold for strong evidence for selection. Outliers on the right side of the vertical line are marked with their names.
Relationships Among Environmental, Epigenetic and Phenotypic Variation
The SEM linked the two flood treatments, variation of the different types of subepiloci and the six measured phenotypic traits. The treatments directly affected all traits; however, petiole length and internode length were only significantly correlated with different flood regimes, with no relationship with subepiloci variation (Figure 4). Leaf area and specific leaf area were related to flood-independent u-subepiloci variation (Figure 4), which may arise from spontaneous epimutation. Flood-induced epigenetic variation affected leaf area by m-subepiloci, specific petiole length by all subepiloci types, specific internode length and specific leaf area by h-subepiloci (Figure 4).
Figure 4. Relationships among environmental, epigenetic and phenotypic variation of Hydrocotyle vulgaris. The epigenetic variables are the first three axes from principal coordinate analyses (PCoA) of different types of subepiloci, i.e., (A) u-subepiloci, (B) m-subepiloci, and (C) h-subepiloci. The lines with arrows indicate significant correlations between variables, with coefficients and P-values showing near the lines.
Discussion
Phenotypic Responses to Different Flood Regimes
Submergence inhibited growth of H. vulgaris, possibly due to that decreased irradiance, sediment anoxia, and osmotic stress in this severe environment restrained plant carbohydrate storage, oxygen transport, and nutrient acquisition (Vretare et al., 2001; Santamaría, 2002). However, in response to semi-submergence, H. vulgaris may develop flood-tolerant responses and soil-oxygen deficiency resistance, such as elongating stout petiole to extend above the water surface and enlarging thick leaf to capture light and increase gas exchange (Vretare et al., 2001; Luo and Xie, 2009). Moreover, the oxygen transported to the node may drive the length extension of internode for further dispersal (Vretare et al., 2001). These changes could confer a fitness benefit for H. vulgaris under semi-submergence, with significant higher aboveground biomass and population density [for ramet aboveground biomass (mean ± SE), submergence = 0.038 ± 0.005 g, semi-submergence = 0.188 ± 0.009 g, F1, 4 = 160.791, P < 0.001).
Epigenetic Responses to Different Flood Regimes
Principal coordinate analyses showed a clear epigenetic differentiation between the semi-submergence and the submergence experimental populations of H. vulgaris, indicating that environmental conditions could shape DNA methylation patterns of plant populations (Note that if DNA methylation changes largely arise from random epimutation, the presence/absence of private bands would be observed in many loci and such epiloci could be neutral so that the treatment-induced epigenetic differentiation would not occur) (Boyko and Kovalchuk, 2010; Schulz et al., 2014; Zhang et al., 2016). Consistent to previous studies (e.g., Gao et al., 2010; Lira-Medeiros et al., 2010; Li et al., 2013; Zoldoš et al., 2018), our results also suggest that experimental H. vulgaris populations can not only respond differently in phenotypic traits, but also undergo a genome-wide epigenetic reprogramming under divergent pressures from contrasting treatments. Therefore, such environment-directed DNA methylation mechanism may be involved in plant adaptation to stress (Boyko and Kovalchuk, 2010; González et al., 2016).
Moreover, the submerged experimental populations of H. vulgaris exhibited greater epigenetic diversity and differentiation than the semi-submerged populations. Some of previous studies reported that severe stress could trigger epigenome variability, providing a possible mechanism for fine-tuning short-term adaptive benefits (Boyko et al., 2010; Verhoeven and Preite, 2014; Colicchio et al., 2018). However, long periods of constant stress can fix the allelic variant that confers tolerance to stress via strong directional selection, leading to the constrained epigenetic diversity and differentiation (Lira-Medeiros et al., 2010; Grativol et al., 2012; Rico et al., 2014).
Population epigenetic differentiation based on u-subepiloci was highly similar to that based on all MSAP subepiloci. Also, h-subepiloci revealed the semi-submergence and the submergence population cluster. This could indicate a functional difference of subepiloci types, with the additive contribution of u-subepiloci and h-subepiloci to population divergence between the two flood regimes. Moreover, AMOVA results showed that variation mainly existed between treatments and within experimental populations based on both u- and h-subepiloci, similar to that based on all subepiloci, whereas most variation existed only within experimental populations based on m-subepiloci. Therefore, the hemimethylation or demethylation in the CHG-context may play a more important role in habitat adjustment in plants than changes of CG-context. Several previous studies have revealed that u-subepiloci and m-subepiloci are more significant in shaping epigenetic structure of natural populations from different habitats (e.g., Schulz et al., 2014; Zoldoš et al., 2018). Such inconsistency suggests that the function of CG- and CHG-methylated states in response to environmental factors is species- and/or environment-specific (Rico et al., 2014; Putnam et al., 2016). We identified ten outlier epiloci facilitated separation of H. vulgaris experimental populations between semi-submergence and submergence (Supplementary Appendixes 4, 5), which may contribute to plastic responses of populations to the flood variation.
Relationships Among Environmental, Epigenetic and Phenotypic Variation
Structural equation modeling analyses showed that petiole length and internode length of H. vulgaris were only significantly correlated with flood, but not with epiloci variation. These results may arise from effects of nutritional or physiological activities, or the low-resolution of MSAP technique (Zhang et al., 2016). Leaf area and specific leaf area are partially affected by u-subepiloci without environmental induction, possibly due to the spontaneous epigenetic variation, arising from imperfect action of enzymes that ensure proper maintenance of epigenetic information through cell division (Verhoeven and Preite, 2014). Stochastic DNA methylation variation is a source for phenotypic diversity in plants, which may mediate phenotypes for several generations that could affect subsequent selection and contribute to adaptive processes (Verhoeven and Preite, 2014; van der Graaf et al., 2015; Groot et al., 2018).
Some phenotypic variation was associated with environment-induced DNA methylation variation through different types of subepiloci, possibly due to their functional differences in regulating gene expression. However, it provides no direct causal information about the region or gene influenced by DNA methylation, as MSAP epiloci are anonymous markers. Our results support the emerging three-way link among flood regimes, DNA methylation and phenotypic changes, suggesting that epigenetic variation might be involved in plastic responses to environmental variation (Wu et al., 2013; Putnam et al., 2016).
Conclusions
We conclude that plants can exhibit significant phenotypic differences between flood regimes, with clear DNA methylation differentiation associated with phenotypes. Moreover, by using the mixed scoring approach, we find the different contributions of methylation types to epigenetic processes in habitat-related responses. Our study potentially adds to the knowledge base of DNA methylation-environmental interactions. However, we did not demonstrate heritability of the epigenetic changes in later-generation and their long-term adaptive and evolutionary implications. Moreover, information on the mechanistic link between methylation and phenotype is still limited. Therefore, more profound studies are needed to deeply uncover the epigenetic role in plant ecological and evolutionary processes.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Author Contributions
M-ZW and F-HY designed the research. M-ZW and H-LL performed the research. M-ZW contributed new reagents or analytical tools, analyzed the data, and wrote the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by the National Natural Science Foundation of China (32071527, 31870610, and 31761123001), Ten Thousand Talent Program of Zhejiang Province (2018R52016), and the Joint Fund of Zhejiang Provincial Natural Science Foundation (LTZ20C030001).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We thank Jing-Xia Shen and Qian Zhang for assistance with the experiment.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.831175/full#supplementary-material
References
Abratowska, A., Wasowicz, P., Bednarek, P. T., Telka, J., and Wierzbicka, M. (2012). Morphological and genetic distinctiveness of metallicolous and non-metallicolous populations of Armeria maritime s.l. (Plumbaginaceae) in Poland. Plant Biol. 14, 586–595. doi: 10.1111/j.1438-8677.2011.00536.x
Angers, B., Castonguay, E., and Massicotte, R. (2010). Environmentally induced phenotypes and DNA methylation: how to deal with unpredictable conditions until the next generation and after. Mol. Ecol. 19, 1283–1295. doi: 10.1111/j.1365-294X.2010.04580.x
Bian, R. J., Nie, D. D., Fu, X., Zhou, X. L., Gao, Y., Bai, Z. J., et al. (2013). Adaptational significance of variations in DNA methylation in clonal plant Hierochloe glabra (Poaceae) in heterogeneous habitats. Aust. J. Bot. 61, 274–282.
Bossdorf, O., Richards, C. L., and Pigliucci, M. (2008). Epigenetics for ecologists. Ecol. Lett. 11, 106–115.
Boyko, A., and Kovalchuk, I. (2010). Genome instability and epigenetic modification-heritable responses to environmental stress? Curr. Opin. Plant Biol. 14, 260–266. doi: 10.1016/j.pbi.2011.03.003
Boyko, A., Blevins, T., Yao, Y., Golubov, A., Bilichak, A., Ilnytskyy, Y., et al. (2010). Transgenerational adaptation of Arabidopsis to stress requires DNA methylation and the function of Dicer-like proteins. PLoS One 5:e9514. doi: 10.1371/journal.pone.0009514
Colicchio, J. M., and Herman, J. (2020). Empirical patterns of environmental variation favor adaptive transgenerational plasticity. Ecol. Evol. 10, 1648–1665. doi: 10.1002/ece3.6022
Colicchio, J. M., Kelly, J. K., and Hileman, L. C. (2018). Parental experience modifies the Mimulus methylome. BMC Genomics 19:746. doi: 10.1186/s12864-018-5087-x
Dong, B. C., Wang, J. Z., Liu, R. H., Zhang, M. X., Luo, F. L., and Yu, F. H. (2015). Soil heterogeneity affects ramet placement of Hydrocotyle vulgaris. J. Plant Ecol. 8, 91–100. doi: 10.1093/jpe/rtu003
Dubin, M. J., Zhang, P., Meng, D. Z., Remigereau, M. S., Osborne, E. J., Casale, F. P., et al. (2015). DNA methylation variation in Arabidopsis has a genetic basis and shows evidence of local adaptation. ELife 4:e05255. doi: 10.7554/eLife.05255
Foll, M., and Gaggiotti, O. (2008). A genome-scan method to identify selected loci appropriate for both dominant and codominant markers, a Bayesian perspective. Genetics 180, 977–993. doi: 10.1534/genetics.108.092221
Gao, L. X., Geng, Y. P., Li, B., Chen, J. K., and Yang, J. (2010). Genome-wide DNA methylation alterations of Alternanthera philoxeroides in natural and manipulated habitats: implications for epigenetic regulation of rapid responses to environmental fluctuation and phenotypic variation. Plant Cell Environ. 33, 1820–1827. doi: 10.1111/j.1365-3040.2010.02186.x
González, A. P. R., Chrtek, J., Dobrev, P. I., Dumalasová, V., Fehrer, J., Mráz, P., et al. (2016). Stress-induced memory alters growth of clonal offspring of white clover (Trifolium repens). Am. J. Bot. 103, 1567–1574. doi: 10.3732/ajb.1500526
Grativol, C., Hemerly, A. S., and Ferreira, P. C. G. (2012). Genetic and epigenetic regulation of stress responses in natural plant populations. Biochim. Biophys. Acta 1819, 176–185. doi: 10.1016/j.bbagrm.2011.08.010
Griffin, P. T., Niederhuth, C. E., and Schmitz, R. J. (2016). A comparative analysis of 5-azacytidine and zebularine induced DNA demethylation. G3 6, 2773–2780. doi: 10.1534/g3.116.030262
Groot, M. P., Wagemaker, N., Ouborg, N. J., Verhoeven, K. J. F., and Vergeer, P. (2018). Epigenetic population differentiation in field- and common garden-grown Scabiosa columbaria plants. Ecol. Evol. 8, 3505–3517. doi: 10.1002/ece3.3931
Huang, X. N., Li, S. G., Ni, P., Gao, Y. C., Jiang, B., Zhou, Z. C., et al. (2017). Rapid response to changing environments during biological invasions: DNA methylation perspectives. Mol. Ecol. 26, 6621–6633. doi: 10.1111/mec.14382
Jueterbock, A., Boström, C., Coyer, J. A., Olsen, J. L., Kopp, M., Dhanasiri, A. A., et al. (2020). The seagrass methylome memorizes heat stress and is associated with variation in stress performance among clonal shoots. Front. Plant Sci. 11:571646. doi: 10.1101/787754
Li, W. G., Chen, W. B., Qi, X. F., Wang, Q. F., and Chen, J. M. (2013). Variation of cytosine methylation in response to water availability in two contrasting growth types of an amphibious plant Alternanthera philoxeroides. Biochem. Syst. Ecol. 50, 175–181.
Lira-Medeiros, C. F., Parisod, C., Fernandes, R. A., Mata, C. S., Cardoso, M. A., and Ferreira, P. C. G. (2010). Epigenetic variation in mangrove plants occurring in contrasting natural environment. PLoS One 5:e10326. doi: 10.1371/journal.pone.0010326
Liu, R. H., Chen, Q. W., Dong, B. C., and Yu, F. H. (2014). Effects of vegetative propagule pressure on the establishment of an introduced clonal plant, Hydrocotyle vulgaris. Sci. Rep. 4:5507. doi: 10.1038/srep05507
Luo, W., and Xie, Y. (2009). Growth and morphological responses to water level and nutrient supply in three emergent macrophyte species. Hydrobiologia 624, 151–160.
Marfil, C. F., Camadro, E. L., and Masuelli, R. W. (2009). Phenotypic instability and epigenetic variability in a diploid potato of hybrid origin, Solanumruiz-lealii. BMC Plant Biol. 9:21. doi: 10.1186/1471-2229-9-21
Miao, L. H., Ji, M. C., Wang, Y. Y., Qiao, D. D., and Chen, Y. C. (2011). Study on invasion risk of Hydrocotyle vulgaris as an alien species in wetlands. J. Zhejiang Univ. 37, 425–431.
Paun, O., Bateman, R. M., Fay, M. F., Hedrén, M., Civeyrel, L., and Chase, M. W. (2010). Stable epigenetic effects impact adaptation in Allopolyploid Orchids (Dactylorhiza: Orchidaceae). Mol. Biol. Evol. 27, 2465–2473. doi: 10.1093/molbev/msq150
Peakall, R., and Smouse, P. E. (2012). GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research-an update. Bioinformatics 28, 2537–2539. doi: 10.1093/bioinformatics/bts460
Putnam, H. M., Davidson, J. M., and Gates, R. D. (2016). Ocean acidification influences host DNA methylation and phenotypic plasticity in environmentally susceptible corals. Evol. Appl. 9, 1165–1178. doi: 10.1111/eva.12408
Rapp, R. A., and Wendel, J. F. (2005). Epigenetics and plant evolution. New Phytol. 168, 81–91. doi: 10.1111/j.1469-8137.2005.01491.x
Richards, C. L., Alonso, C., Becker, C., Bossdorf, O., Bucher, E., Colomé-Tatché, M., et al. (2017). Ecological plant epigenetics: evidence from model and non-model species, and the way forward. Ecol. Lett. 20, 1576–1590. doi: 10.1111/ele.12858
Richards, C. L., Bossdorf, O., and Pigliucci, M. (2010). What role does heritable epigenetic variation play in phenotypic evolution? BioScience 60, 232–237. doi: 10.1525/bio.2010.60.3.9
Richards, E. J. (2006). Inherited epigenetic variation-revisiting soft inheritance. Nat. Rev. Genet. 7, 395–401. doi: 10.1038/nrg1834
Rico, L., Ogaya, R., Barbeta, A., and Penuelas, J. (2014). Changes in DNA methylation fingerprint of Quercus ilex trees in response to experimental field drought simulating projected climate change. Plant Biol. 16, 419–427. doi: 10.1111/plb.12049
Robertson, M., Aaron, S., Ashley, S., Christina, J. M., and Christina, L. R. (2017). Genetic and epigenetic variation in Spartina alterniflora following the Deepwater Horizonoil spill. Evol. Appl. 10, 792–801. doi: 10.1111/eva.12482
Robertson, M., and Richards, C. (2015). Opportunities and challenges of next-generation sequencing applications in ecological epigenetics. Mol. Ecol. 24, 3799–3801. doi: 10.1111/mec.13277
Santamaría, L. (2002). Why are most aquatic plants widely distributed? Dispersal, clonal growth and small-scale heterogeneity in a stressful environment. Acta Oecol. 23, 137–154. doi: 10.1016/s1146-609x(02)01146-3
Schulz, B., Eckstein, R. L., and Durka, W. (2013). Scoring and analysis of methylation-sensitive amplification polymorphisms for epigenetic population studies. Mol. Ecol. Resour. 13, 642–653. doi: 10.1111/1755-0998.12100
Schulz, B., Eckstein, R. L., and Durka, W. (2014). Epigenetic variation reflects dynamic habitat conditions in a rare floodplain herb. Mol. Ecol. 23, 3523–3537. doi: 10.1111/mec.12835
Trucchi, E., Mazzarella, A. B., Gilfillan, G. D., Lorenzo, M. T., Schonswetter, P., and Paun, O. (2016). BsRADseq: screening DNA methylation in natural populations of non-model species. Mol. Ecol. 25, 1697–1713. doi: 10.1111/mec.13550
van der Graaf, A., Wardenaar, R., Neumann, D. A., Taudt, A., Shaw, R. G., Jansen, R. C., et al. (2015). Rate, spectrum, and evolutionary dynamics of spontaneous epimutations. Proc. Natl. Acad. Sci. U.S.A. 112, 6676–6681. doi: 10.1073/pnas.1424254112
Verhoeven, K. J. F., and Preite, V. (2014). Epigenetic variation in asexually reproducing organisms. Evolution 68, 644–655. doi: 10.1111/evo.12320
Verhoeven, K. J. F., Jansen, J. J., van Dijk, P. J., and Biere, A. (2010). Stress-induced DNA methylation changes and their heritability in asexual dandelions. New. Phytol. 185, 1108–1118. doi: 10.1111/j.1469-8137.2009.03121.x
Vretare, V., Weisner, S. E. B., Strand, J. A., and Granéli, W. (2001). Phenotypic plasticity in Phragmites australis as a functional response to water depth. Aquat. Bot. 69, 127–145. doi: 10.1016/s0304-3770(01)00134-6
Wang, M. Z., Li, H. L., Li, J. M., and Yu, F. H. (2020). Correlations between genetic, epigenetic and phenotypic variation of an introduced clonal herb. Heredity 124, 146–155. doi: 10.1038/s41437-019-0261-8
Watson, R. G. A., Baldanzi, S., Pérez-Figueroa, A., Gouws, G., and Porri, F. (2018). Morphological and epigenetic variation in mussels from contrasting environments. Mar. Biol. 165:50.
Wu, W., Yi, M. R., Wang, X., Ma, L., Jiang, L., Li, X., et al. (2013). Genetic and epigenetic differentiation between natural Betulaermanii (Betulaceae) populations inhabiting contrasting habitats. Tree Genet. Genomes 9, 1321–1328. doi: 10.1007/s11295-013-0641-9
Zhang, Y. Y., Parepa, M., Fischer, M., and Bossdorf, O. (2016). “Epigenetics of colonizing species? A study of Japanese Knotweed in central Europe,” in Invasion Genetics: The Baker and Stebbins Legacy, eds S. C. H. Barrett, R. I. Colautti, K. M. Dlugosch, and L. H. Rieseberg (Oxford: John Wiley & Sons, Ltd).
Keywords: artificial populations, epigenetic variation, clonal plant, phenotypes, Hydrocotyle vulgaris, flooding
Citation: Wang M-Z, Li H-L, Tang M and Yu F-H (2022) DNA Methylation Correlates With Responses of Experimental Hydrocotyle vulgaris Populations to Different Flood Regimes. Front. Plant Sci. 13:831175. doi: 10.3389/fpls.2022.831175
Received: 08 December 2021; Accepted: 11 February 2022;
Published: 07 March 2022.
Edited by:
Mangal Singh Rathore, Council of Scientific and Industrial Research (CSIR), IndiaReviewed by:
Vinod Singh Gour, Amity University Jaipur, IndiaMichael James Wilkinson, Aberystwyth University, United Kingdom
Copyright © 2022 Wang, Li, Tang and Yu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Fei-Hai Yu, feihaiyu@126.com