Partial maintenance of organ-specific epigenetic marks during plant asexual reproduction leads to heritable phenotypic variation

Contributed by Detlef Weigel, July 31, 2018 (sent for review March 28, 2018; reviewed by Xiaofeng Cao and Mary Gehring)
September 10, 2018
115 (39) E9145-E9152

Significance

While clonally propagated individuals should share identical genomes, there is often substantial phenotypic variation among them. Both genetic and epigenetic modifications induced during regeneration have been associated with this phenomenon. Here we investigated the fate of the epigenome after asexual propagation by generating clonal individuals from differentiated somatic cells through the manipulation of a zygotic transcription factor. We found that phenotypic novelty in clonal progeny was linked to epigenetic imprints that reflect the organ used for regeneration. Some of these organ-specific imprints can be maintained during the cloning process and subsequent rounds of meiosis. Our findings are fundamental for understanding the significance of epigenetic variability arising from asexual reproduction and have significant implications for future biotechnological applications.

Abstract

Plants differ from animals in their capability to easily regenerate fertile adult individuals from terminally differentiated cells. This unique developmental plasticity is commonly observed in nature, where many species can reproduce asexually through the ectopic initiation of organogenic or embryogenic developmental programs. While organ-specific epigenetic marks are not passed on during sexual reproduction, the fate of epigenetic marks during asexual reproduction and the implications for clonal progeny remain unclear. Here we report that organ-specific epigenetic imprints in Arabidopsis thaliana can be partially maintained during asexual propagation from somatic cells in which a zygotic program is artificially induced. The altered marks are inherited even over multiple rounds of sexual reproduction, becoming fixed in hybrids and resulting in heritable molecular and physiological phenotypes that depend on the identity of the founder tissue. Consequently, clonal plants display distinct interactions with beneficial and pathogenic microorganisms. Our results demonstrate how novel phenotypic variation in plants can be unlocked through altered inheritance of epigenetic marks upon asexual propagation.
Compared with animals, in plants somatic cells can be much more easily coaxed into regenerating entire individuals (1). A potential reason why differentiated plant cells can rapidly acquire “stemness” is that the epigenome of plants is much more flexible than that of animals. Thus, asexual reproduction is much more common in plants than in animals (2, 3), and this has been traditionally exploited by humans for the clonal propagation and genetic manipulation of many economically important plant species, including grapevines, nearly all tuber and root crops, and fruit and forest trees (4).
Although clonal propagation provides ecological and evolutionary benefits, the resulting restricted genetic variation could be detrimental to fitness. Notably, clonally propagated plants are not always phenotypically identical to their parents, a phenomenon termed somaclonal variation. While somaclonal variation is often attributed to the accumulation of random genetic mutations in form of single-base changes or transposon activation (5, 6), increasing evidence suggests that genetic changes are not solely responsible. For example, genome-wide DNA methylation patterns can be perturbed during clonal propagation through hormone-dependent tissue culture, and these epimutations can be associated with altered expression of protein-coding genes (710). In some cases, such tissue-culture induced epimutations can be stably inherited across multiple sexual generations (7, 8, 10) and can be responsible for phenotypes that distinguish clonal descendants from their parents (11). Increased genetic and epigenetic diversity may be deleterious, but also potentially result in advantageous traits. For all these reasons, we would like to better understand the precise origin and mechanistic basis of the molecular and phenotypic changes created during plant regeneration.
Depending on the species, different types of organs, such as roots, stems, and leaves, as well as entire embryos, can be used for clonal propagation. The starting material thus includes tissue with different patterns of gene expression and epigenetic profiles. For example, 1.6% of methylated regions in the Arabidopsis thaliana genome differ in methylation status between shoots and roots, and more than 2,000 genes differ in expression levels (12).
An attractive hypothesis is that complete or partial maintenance of variant epigenetic landscapes present in the starting material will contribute to phenotypic variation among regenerants. We have tested this hypothesis by regenerating A. thaliana plants from different postembryonic organs. Frequent limitations of plant regeneration are the lengthy periods of tissue culture and the hormone mixtures that need to be optimized for each tissue, which confounds the interpretation of phenotypic differences in regenerated plants. To circumvent this problem, we have exploited the fact that somatic embryogenesis can be induced in different tissues of A. thaliana by ectopic expression of certain zygotic factors (2). Here we used transgenic lines carrying an inducible version of the RWP-RK DOMAIN-CONTAINING 4 (RKD4) transcription factor gene (13). We provide evidence that certain organ-specific epigenetic differences are partially maintained during regeneration of A. thaliana. Whole plants derived from roots inherit many aspects of root-specific methylation and gene expression patterns not just in roots, but also in leaves. These epigenetic profiles and the resulting macroscopic and molecular phenotypes are stably transmitted during meiosis for at least four self-crossing generations. Our findings demonstrate that plants with novel methylation and gene expression patterns, as well as physiological phenotypes, can be created using specific regeneration strategies.

Results

Tissue Origin of Regenerants Affects Activity of Defense-Related Genes.

Plants can reproduce asexually from both belowground and aboveground organs, which are known to be epigenetically distinct (12, 14). We took advantage of this situation to determine the extent to which the organ-specific origin of the epigenome could influence phenotypes of clonal progeny. To mimic naturally occurring events associated with asexual propagation (2, 15), we did not resort to hormone-induced regeneration in tissue culture, but instead produced somatic embryos from distinct root (root origin; RO) and leaf (leaf origin; LO) tissues of A. thaliana (Col-0 strain) by controlled expression of the RKD4 zygotic factor (13) (SI Appendix, Fig. S1). We collected seeds from independently regenerated G0 individuals after self-pollination and further propagated each line by selfing for over three consecutive generations (G1–G3) (Fig. 1A). Visual examination revealed no obvious morphological differences between RO and LO plants.
Fig. 1.
Differential expression of defense-related genes in RO and LO plants. (A) Experimental design. RO and LO regenerated plants (n = 5 independent lines each) were propagated through self-fertilization over three generations. (B) GO analysis of DEGs (FDR <0.01; absolute log-twofold change >1.5) between leaves of RO and LO G2 plants revealing enrichment for defense-related functions. (C) Heatmap of scaled, normalized log-transformed read counts for genes underlying the significant enrichment of GO terms in B. FPKM, fragments per kilobase per million. (D) Interaction gene network of DEGs. To facilitate network analysis, we lowered the FDR threshold to <0.05. Of the resulting 4,585 DEGs, 2,752 were represented in the ANAP protein interaction database. Nodes represent genes; triangles highlight defense-related genes from B; edges indicate evidence for gene interaction. Blue filled circles, high expression in LO plants; orange filled circles, high expression in RO plants; red outlined circles, genes with stress-related GO terms.
To determine any potential differences at the molecular level, we performed whole-genome transcriptome analyses in five randomly selected G2 lines. Using stringent thresholds [false discovery rate (FDR) <0.01; absolute log-twofold change >1.5], only 13 differentially expressed genes (DEGs) differentiated roots of RO plants from roots of LO plants, but almost 20-fold more DEGs (n = 239) were identified between leaves of RO and LO plants (SI Appendix, Fig. S2 and Dataset S1). Gene ontology (GO) analysis revealed that these DEGs were enriched for stress- and defense-related genes (FDR <0.05) (Fig. 1B). The 51 genes in these two categories were mostly up-regulated in leaves of RO plants and primarily involved in cellular signaling (27%) and transcriptional regulation (22%). Clustering of all samples based on DEG expression levels suggested that leaves of RO plants had partial root characteristics, sharing apparent similarity with roots from both RO and LO plants (Fig. 1C); in contrast, leaves of LO plants formed a single distinct group (Fig. 1C). To gain further insight into this transcriptional variation, we performed a network analysis of DEGs that distinguished leaves of RO and LO plants, based on a less stringent threshold (FDR <0.05) (16). This network contained the 239 stringent DEGs (SI Appendix, Fig. S2 and Dataset S1), of which 213 had known interactions with other DEGs and 69 were part of a single functional network (Fig. 1D). This network included five genes (WRKY6, SZF1, PYL5, PUB23, and DRIP2) that had been previously implicated in the negative regulation of abiotic and biotic stress responses (1721).

Tissue of Origin of Regenerants Affects Interaction with Microbes.

We next sought to define whether the gene expression differences between RO and LO plants were functionally meaningful and affected whole-plant or tissue-specific functions. Plant microbiota are known to be affected by developmental factors as well as immune system activity (22, 23), and both the altered expression of tissue-specific genes in RO plants as well as the enriched GO categories suggested that it would be relevant to test for a change in the interaction with the biotic environment. We grew plants in natural soils and after 4 weeks assessed the bacterial communities that became associated with their roots. Bacterial communities from roots of RO plants differed from those of both nonregenerated plants (ANOSIM, r = 0.485, P = 0.007) and LO plants (ANOSIM, r = 0.216, P = 0.056) (Fig. 2A and SI Appendix, Fig. S3). In addition, in vitro root colonization by Bacillus amyloliquefaciens FZB42, a soil-borne plant growth-promoting bacterium (24), differed significantly between RO and nonregenerated plants (Wilcoxon rank-sum test, W = 3862.5, P = 5.9e-05) (SI Appendix, Fig. S4). Finally, we inoculated roots of regenerated and nonregenerated plants with synthetic communities (SynComs) consisting of abundant soil- and root-derived bacterial isolates (25). Again, the root bacterial communities differed between regenerated and nonregenerated individuals (7.3% variance explained by genotype; permutation-based ANOVA, P = 0.07) (SI Appendix, Fig. S5); this was especially obvious for the Alcaligenaceae family (Betaproteobacteria) (Fig. 2B and SI Appendix, Fig. S5). Similarly, when analyzing leaves of SynCom-inoculated plants (9.6% variance explained by genotype; P = 0.023), the most notable differences were observed for Xanthomonadaceae (Gammaproteobacteria) (ANOVA, P < 0.05) (SI Appendix, Figs. S5 and S6), which include several phytopathogenic strains.
Fig. 2.
Plants regenerated from roots or leaves interact differently with beneficial and pathogenic microbes. (A) PCA of Bray–Curtis distances of bacterial communities present in roots of nonregenerated and regenerated plants grown in natural soils (n = 10). (B) Canonical analysis of principal coordinates (based on Bray–Curtis distances) showing different root-associated communities of SynComs colonized on nonregenerated and regenerated plants (n = 12). (C) Susceptibility of nonregenerated and regenerated plants to P. syringae pv. tomato strain DC3000 infection. Bacterial growth was determined at 3 d after inoculation (100 cfu mL−1). Data are mean ± SD values from three independent experiments. Statistical significance according to Fisher’s exact test: **P < 0.01, ***P < 0.001; NS, not significant. (D) Susceptibility of nonregenerated and regenerated plants to H. arabidopsidis (Hpa) Noks1 infection, as indicated by the number of conidiospores on leaves at 3 d after inoculation with a suspension of 300,000 spores/mL. Data are mean ± SD values from two independent experiments. P values were determined using Student’s t test.
These differences in microbiota led us to further assess the response to known pathogens. For this, we inoculated leaves with the bacterium Pseudomonas syringae pv. tomato strain DC3000 and the oomycete Hyaloperonospora arabidopsidis isolate Noks1, as the A. thaliana Col-0 strain lacks gene-for-gene resistance to both of these pathogens (26, 27). We found that RO plants were more sensitive than LO plants to infection by either pathogen, and that these differences were stably inherited for at least three generations (Fig. 2 C and D). This result was in line with our previous observation of the up-regulation of negative regulators of biotic stress in RO plants (Fig. 1D).

Heritable Differences in Genome-Wide DNA Methylation in Regenerants.

Given the unexpected widespread and heritable differences between transcriptional profiles and microbe responses in RO and LO plants, which should have very few, if any, group-wise genetic differences, we investigated the potential epigenomic basis for phenotypic differentiation between RO and LO plants. DNA methylation is an important epigenetic mark, for which excellent statistical methods for genome-wide comparisons are available (28). Moreover, dynamic changes in DNA methylation during sexual reproduction have been well documented (29). We used whole-genome bisulfite sequencing of leaves and roots from RO and LO individuals to monitor methylome changes over three consecutive generations (Dataset S2). Principal component analysis (PCA) of 736,413 differentially methylated positions (DMPs) discovered in pairwise contrasts revealed clear differences between root and leaf samples from regenerated and nonregenerated plants (PC1; Fig. 3A). Regenerated samples clustered according to their tissue of origin before regeneration and not, as might be expected, by tissue identity at the time of DNA extraction (PC2; Fig. 3A). When repeating PCA with the methylation level of positions not classified as DMPs, we found little residual variance, indicating a low number of false negatives (SI Appendix, Fig. S7).
Fig. 3.
DNA methylation variation in regenerated plants is stably inherited after selfing and back-crossing. (A) PCA of DNA methylation levels at DMPs identified from pairwise sample comparisons. Numbers in brackets indicate the fraction of overall variance explained by the respective PC. (B) Clustering of LO and RO leaf samples in generation G2, based on 765 DMRs identified in all-against-all pairwise comparisons. (C) PCA of methylation at 255 DMRs identified in G2 RO vs. LO leaf comparison, divided by cytosine sequence context. (Right) PCA on methylation in all contexts within randomly chosen non-DMRs. Numbers in brackets indicate the amount of variance explained by the respective PC. (D) Gains and losses of DNA methylation in DMRs identified between RO and LO leaves in the G2 generation. Color keys indicate methylation rate differences in relation to leaves of nonregenerated plants. (Right) Differences in a random subset of non-DMRs. (E) Methylation frequencies at DMRs in leaves of nonregenerated Col-0, RO, and reciprocal crosses (F1) between nonregenerated Col-0 and RO plants. (F) Methylation analysis of progeny from F1 reciprocal hybrids. The heatmap shows DMR methylation levels in individual F1 hybrid plants (#7 and #16) and each of four independent descendants. The bar plot shows the frequency of hypomethylation in F2 plants of DMRs that were hypomethylated in the F1 hybrid.
Because the functional relevance of individual DMPs in plants is unclear, we also analyzed differentially methylated regions (DMRs). Cluster analysis based on 765 DMRs identified in pairwise comparisons between all G2 leaf samples validated the use of the independent RO and LO lines as replicate groups (Fig. 3B). We found 255 consistent RO vs. LO DMRs in G2 leaves (Dataset S3). Compared with methylated regions in general, RO vs. LO DMRs were overrepresented in exons of gene coding sequences and in 2-kb regions flanking genes (SI Appendix, Fig. S8). Similar distributions of DMRs across genomic features have been reported in the context of biotic and abiotic stresses (3032), suggesting that gene-proximal DNA methylation changes might be a common feature of induced epimutations.
Methylation patterns are known to differ between shoots and roots in A. thaliana and closely related species (12, 14). Methylation in LO leaves was similar to that in nonregenerated leaves (PC 1 in Fig. 3C); likewise, roots from RO samples showed patterns similar to those in nonregenerated roots. However, leaves from RO populations had methylation patterns closer to those of roots, especially in a symmetric cytosine context (Fig. 3C). In contrast, methylation levels at these DMRs in roots of LO plants were similar to those seen in nonregenerated root samples (Fig. 3C). As expected, methylation levels at non-DMRs grouped samples primarily by their tissue identity, regardless of regenerant origin (Fig. 3C).
Our DMR analysis also revealed that leaves from RO plants had less overall CG and CHG (but not CHH) methylation in the identified DMRs compared with leaves of nonregenerated plants or LO plants (Fig. 3D). To account for sampling bias and stochastic effects, we repeated the analyses using randomly selected non-DMRs, which did not produce any evidence of reduced methylation in roots of RO plants (Fig. 3D). The reduced DNA methylation in RO leaves was stably inherited over at least three generations (Fig. 3 C and D), indicating that root-specific DNA methylation patterns are not lost during subsequent sexual reproduction.
To understand the importance of the methylome status at the time of regeneration, we compared the leaf methylome of RO plants with the methylation profiles of roots from nonregenerated seedlings and mature plants (SI Appendix, Fig. S9). The RO leaf methylome was closer to that of roots from nonregenerated seedlings, implying that regeneration had been induced at the seedling stage and that the cell-specific methylation pattern was maintained throughout the regeneration process and in subsequent sexually reproduced progeny.
Because there is evidence of methylation information being transferable between chromosomes (33, 34), we tested whether such information transfer could occur in trans in our system. To this end, we performed reciprocal crosses between RO and nonregenerated plants. DNA methylation in both F1 hybrids was at midparent values, indicating that the DMRs on chromosomes inherited from the RO parent retained their hypomethylated status (Fig. 3E and SI Appendix, Fig. S10). We assessed the heritability of the observed differential methylation at these loci by sequencing individual F2 plant progenies. More than two-thirds (80%) of DMRs with midparental methylation in the F1 hybrid retained their hypomethylated state in at least one F2 descendant, indicating that allele-specific methylation was stably inherited through both mitotic and meiotic divisions (Fig. 3F).
The establishment and maintenance of DNA methylation in plants rely on a series of partially interconnected pathways, depending on the genomic features that are methylated (35). Consequently, when genome-wide demethylation is induced by various mutations, some regions can be remethylated upon restoration of the methylation machinery, while others cannot, and this is a function of the underlying methylation pathways (3638). To gain insight into the mechanisms that supported the partial maintenance of root-specific methylation patterns in RO plants, we investigated whether methylated regions that became hypomethylated in RO plants were under the control of a specific epigenetic pathway, by comparing the methylation changes in different mutant contexts (39). We found that CG methylation in such regions was affected in the chromatin remodeling mutant ddm1 and in the DNA methylation maintenance mutants met1 and vim1 vim2 vim3, while CHG methylation was altered in the de novo methyltransferase cmt3 mutant and in mutants with a compromised H3K9 methylation machinery (SI Appendix, Fig. S11).
Contrary to our expectations, the regions that were hypo-methylated in RO plants were not affected in the triple ros1 dml2 dml3 (rdd) mutant, which lacks three DNA demethylases, suggesting that DNA hypomethylation in RO plants is due to differences in the establishment and/or maintenance of DNA methylation in root initials during embryogenesis, rather than to subsequent active demethylation. Induced methylation changes are often associated with changes in activity of transposable elements (TEs), which in turn is reflected in altered expression patterns of siRNAs, which are a central component of TE silencing pathways that involve RNA-directed DNA methylation (RdDM) (4042). LO vs. RO DMRs tend to not directly overlap with TEs and 24-nt siRNA loci (43), but are often closer to them than would be expected by chance (SI Appendix, Fig. S12).

Differential Methylation Affects Expression of RSM1.

To establish a connection between the molecular and phenotypic variation generated by plant regeneration, we searched for correlations between changes in DNA methylation (765 DMRs) and gene expression in leaves of regenerated plants, using a combined set of 1,537 DEGs from G1 and G2 generations (FDR <0.01). We found 29, mostly hypomethylated, DMRs that were within 2 kb upstream or downstream of DEGs (SI Appendix, Fig. S13 and Dataset S4). To confirm that such DMRs can indeed affect gene expression, we selected a DMR approximately 1 kb downstream of RADIALIS-LIKE SANT/MYB1 (RSM1)/MATERNAL EFFECT EMBRYO ARREST 3 (MEE3), a regulator of growth and flowering time (44, 45). This DMR was hypomethylated in RO leaves compared with LO leaves, in both G1 and G2 generations (SI Appendix, Fig. S14). To manipulate DNA methylation experimentally, we introduced an inverted repeat (IR) hairpin to force DNA hypermethylation at the RSM1-DMR (RSM1-IR) by RdDM. Bisulfite sequencing confirmed that the targeted genome region became specifically hypermethylated in IR transgenic lines (Fig. 4A), and this was accompanied by reduced RSM1 expression (Fig. 4B). RSM1_IR lines suffered from pleiotropic developmental defects, including accelerated senescence and earlier flowering (Fig. 4C). When we transformed RSM1_IR plants with a synthetic RSM1 construct (synRSM1) resistant to RSM1_IR targeting, the RSM1-IR developmental phenotypes were suppressed (Fig. 4D), indicating that RSM1 is the causative gene responsible for the underlying developmental defects observed in IR lines. In aggregate, these observations are consistent with the DMR containing a regulatory region, the activity of which is influenced by DNA methylation.
Fig. 4.
Activity of an RSM1 regulatory element is affected by DNA methylation. (A) Snapshot of DNA methylation in different sequence contexts at the RSM1 locus in wild-type (WT) Col-0 plants and plants carrying an inverted repeat transgene (RSM1-IR). The black horizontal bar indicates the region targeted by RSM1-IR–induced RdDM for DNA hypermethylation. Green ticks, CG methylation; blue, CHG methylation; orange, CHH methylation. (B) RSM1 expression in leaves and roots of WT plants, RSM1-IR plants, and RSM1-IR plants complemented with synRSM1. White bar, leaf; black bar, root. (C) Induced DNA hypermethylation of RSM1-DMR affects morphology. (D) Genetic complementation of RSM1-IR with a synthetic RSM1 transgene (synRSM1) resistant to RSM1-IR–induced hypermethylation.

Discussion

We have set out to determine whether epigenetic phenomena contribute to phenotypic variation arising from asexual propagation in plants. Our data show that the epigenetic profile of certain somatic cells can be at least partially maintained when clonal regenerants are produced by inducing an embryogenic developmental program in these somatic cells. Inappropriate maintenance of somatic epigenetic profiles has also been observed when mammalian clones are produced by nuclear transfer, where it can result in embryonic lethality and postnatal growth defects (46, 47).
Because somatic embryogenesis closely mimics many aspects of normal embryogenesis (48), the partial retention of tissue-specific epigenetic signatures in the primary regenerants is not so surprising, as these marks have not passed through gametogenesis and early stages of embryogenesis. A less expected finding is that progeny of regenerated plants, which have undergone multiple cycles of sexual reproduction, can continue to retain some organ-specific epigenetic marks, as well as transcriptional and phenotypic signatures typical of the founder tissue used for the initial propagation. The stability of these epigenetic marks and the reduced methylation levels of the affected loci resemble epialleles induced by the inactivation of positive regulators of DNA methylation. Such epialleles have been shown to be heritable over multiple generations in epigenetic recombinant inbred lines (36, 37, 49), where they are also associated with stable phenotypic variation (36, 50).
Plant regeneration induced by phytohormones can lead to heritable genome-wide DNA hypomethylation that is correlated with altered gene expression. In both A. thaliana and rice, such epimutations seem to occur in a stochastic manner (7, 9), but they appear to be more consistent with and to be targeted to specific genomic regions in maize (10). While we also observe some random epigenetic changes after somatic embryogenesis induced by zygotic factors, most methylation changes following tissue-specific regeneration are reproducible and reflect the epigenetic state of the tissue of origin.
In A. thaliana, asymmetric CHH methylation is actively reprogrammed during male gametogenesis and proper methylation is rapidly restored in the embryo after fertilization (5153). On the other hand, symmetric CHG and CG methylation is thought to remain stable throughout both male and female gametogenesis and in embryos (5155). Unfortunately, detailed methylation maps of embryonic lineages and stem cells are not available, but the leaf methylation pattern can be envisioned to resemble a ground state more closely than the root methylation pattern, which is distinguished from the leaf pattern primarily by hypomethylation. With roots being further away from the ground state than leaves, one might then expect RO plants to be epigenetically more distinct from nonregenerated plants compared with LO plants.
Targeted hypomethylation likely takes place already during embryogenesis, when root cell initials are first specified, because hypomethylation is a characteristic of many cell types in the postembryonic root (51). The mechanisms underlying epigenetic reprogramming in roots are unknown, but methylation differences between roots and leaves have been found in other plant species as well (12, 14, 56, 57). One possibility is that such marks are involved in regulating distinct tissue-specific transcriptional responses to environmental factors, such as exposure and interactions with microorganisms. In support of this argument, A. thaliana leaves challenged with bacterial pathogens rapidly remodel their DNA methylation profiles (31), and root resistance to fungal pathogens requires active DNA demethylation (58). However, whether DNA methylation is the primary epigenetic mark that changes upon regeneration and whether it is the most important cause for the expression and phenotypic differences observed remains to be clarified.
In summary, we propose that organ-specific epigenetic marks captured during somatic embryogenesis are also likely to contribute to the phenotypic somaclonal variation observed in plants propagated in vitro (5, 11, 59, 60), as well as in natural asexually reproducing plant populations (61). The partial maintenance of organ-specific epigenetic marks during cloning might not be a unique feature of plants, and may also occur in clonal animals. Our findings thus not only suggest new and exciting possibilities for enhancing or, perhaps more importantly, limiting phenotypic variation in clonally propagated elite lines (11), but also raise pertinent questions regarding the adaptive significance of developmental epigenomic changes captured during asexual reproduction in plants.

Materials and Methods

A. thaliana Col-0 was grown in a controlled environment (16-h light/8-h dark, 22 °C). For the direct plant regeneration of plants from differentiated organs, we used transgenic lines carrying a dexamethasone-inducible transgene to overexpress the GRANDE (GRD)/RWPRK motif-containing (RKD4) transcription factor (13). Seeds from a transgenic indRKD4 line were germinated on Murashige and Skoog (MS) plates. After 4–6 d, plants were transferred to MS plates with 20 µM dexamethasone, on which they were incubated for 7 d. Plants were transferred to dexamethasone-free MS plates to allow the formation of somatic embryos, which were dissected manually by micromanipulation using tungsten needles. Isolated somatic embryos were transferred to MS plates to aid whole plant regeneration from leaves (LO) or roots (RO). We generated 10 independent regenerants from each organ (G0 generation), which were grown in soil to produce seeds. We grew 24 plants from each line (G1 generation) and selected 10 individuals at random for each regenerated line to produce seed. We continued the sexual propagation of each regenerant population for two additional generations (G2 and G3) following the same scheme.
Additional information is provided in SI Appendix, Materials and Methods.

Data Availability

Data deposition: Sequencing reads have been deposited at the European Nucleotide Archive (accession nos. PRJEB26932 and PRJEB14117). DNA methylation data have been uploaded to the EPIC-CoGe epigenome browser and can be accessed at https://genomevolution.org/CoGe/NotebookView.pl?lid=835.

Acknowledgments

We thank C. Lanz and J. Hildebrandt for help with Illumina sequencing, J. Engelhorn for assistance with RNA sampling, Rainer Borriss for providing the B. amyloliquefaciens FZB42 strain, and Liliana M. Costa for discussions and comments on the manuscript. This work was supported by a European Research Council (ERC) Marie Sklodowska-Curie Fellowship (751204-H2020-MSCA-IF-2016, to A.W.); the ERC AdG IMMUNEMESIS Project, the Deutsche Forschungsgemeinschaft SPP1529 Program, and the Max Planck Society (D.W.); the ERC AdG ROOTMICROBIOTA Project, the Cluster of Excellence on Plant Sciences, and the Max Planck Society (P.S.-L.); and the Biotechnology and Biological Sciences Research Council (Grants BB/L025892/1, to G.D.B., and BB/L003023/1, BB/N005279/1, BB/N00194X/1, and BB/P02601X/1, to J.G.-M.).

Supporting Information

Appendix (PDF)
Dataset_S01 (XLSX)
Dataset_S02 (XLSX)
Dataset_S03 (XLSX)
Dataset_S04 (XLSX)

References

1
K Sugimoto, SP Gordon, EM Meyerowitz, Regeneration in plants and animals: Dedifferentiation, transdifferentiation, or just differentiation? Trends Cell Biol 21, 212–218 (2011).
2
M Ikeuchi, Y Ogawa, A Iwase, K Sugimoto, Plant regeneration: Cellular origins and molecular mechanisms. Development 143, 1442–1451 (2016).
3
K Sugimoto, Y Jiao, EM Meyerowitz, Arabidopsis regeneration from multiple tissues occurs via a root development pathway. Dev Cell 18, 463–471 (2010).
4
D McKey, M Elias, B Pujol, A Duputié, The evolutionary ecology of clonally propagated domesticated plants. New Phytol 186, 318–332 (2010).
5
C Jiang, et al., Regenerant Arabidopsis lineages display a distinct genome-wide spectrum of mutations conferring variant phenotypes. Curr Biol 21, 1385–1390 (2011).
6
A Miyao, et al., Molecular spectrum of somaclonal variation in regenerated rice revealed by whole-genome sequencing. Plant Cell Physiol 53, 256–264 (2012).
7
H Stroud, et al., Plants regenerated from tissue culture contain stable epigenome changes in rice. eLife 2, e00354 (2013).
8
SC Stelpflug, SR Eichten, PJ Hermanson, NM Springer, SM Kaeppler, Consistent and heritable alterations of DNA methylation are induced by tissue culture in maize. Genetics 198, 209–218 (2014).
9
M Tanurdzic, et al., Epigenomic consequences of immortalized plant cell suspension culture. PLoS Biol 6, 2880–2895 (2008).
10
Z Han, et al., Heritable epigenomic changes to the maize methylome resulting from tissue culture. Genetics 209, 983–995 (2018).
11
M Ong-Abdullah, et al., Loss of Karma transposon methylation underlies the mantled somaclonal variant of oil palm. Nature 525, 533–537 (2015).
12
N Widman, S Feng, SE Jacobsen, M Pellegrini, Epigenetic differences between shoots and roots in Arabidopsis reveals tissue-specific regulation. Epigenetics 9, 236–242 (2014).
13
T Waki, T Hiki, R Watanabe, T Hashimoto, K Nakajima, The Arabidopsis RWP-RK protein RKD4 triggers gene expression and pattern formation in early embryogenesis. Curr Biol 21, 1277–1281 (2011).
14
DK Seymour, D Koenig, J Hagmann, C Becker, D Weigel, Evolution of DNA methylation patterns in the Brassicaceae is driven by differences in genome organization. PLoS Genet 10, e1004785 (2014).
15
HM Garcês, D Koenig, BT Townsley, M Kim, NR Sinha, Truncation of LEAFY COTYLEDON1 protein is required for asexual reproduction in Kalanchoë daigremontiana. Plant Physiol 165, 196–206 (2014).
16
D Warde-Farley, et al., The GeneMANIA prediction server: Biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res 38, W214–W220 (2010).
17
P Lackman, et al., Jasmonate signaling involves the abscisic acid receptor PYL4 to regulate metabolic reprogramming in Arabidopsis and tobacco. Proc Natl Acad Sci USA 108, 5891–5896 (2011).
18
F Qin, et al., Arabidopsis DREB2A-interacting proteins function as RING E3 ligases and negatively regulate plant drought stress-responsive gene expression. Plant Cell 20, 1693–1707 (2008).
19
S Robatzek, IE Somssich, Targets of AtWRKY6 regulation during plant senescence and pathogen defense. Genes Dev 16, 1139–1149 (2002).
20
J Sun, et al., The CCCH-type zinc finger proteins AtSZF1 and AtSZF2 regulate salt stress responses in Arabidopsis. Plant Cell Physiol 48, 1148–1158 (2007).
21
M Trujillo, K Ichimura, C Casais, K Shirasu, Negative regulation of PAMP-triggered immunity by an E3 ubiquitin ligase triplet in Arabidopsis. Curr Biol 18, 1396–1401 (2008).
22
SL Lebeis, et al., Salicylic acid modulates colonization of the root microbiome by specific bacterial taxa. Science 349, 860–864 (2015).
23
JM Chaparro, DV Badri, JM Vivanco, Rhizosphere microbiome assemblage is affected by plant development. ISME J 8, 790–803 (2014).
24
SP Chowdhury, et al., Effects of Bacillus amyloliquefaciens FZB42 on lettuce growth and health under pathogen pressure and its impact on the rhizosphere bacterial community. PLoS One 8, e68818 (2013).
25
Y Bai, et al., Functional overlap of the Arabidopsis leaf and root microbiota. Nature 528, 364–369 (2015).
26
DF Tomé, J Steinbrenner, JL Beynon, A growth quantification assay for Hyaloperonospora arabidopsidis isolates in Arabidopsis thaliana. Methods Mol Biol 1127, 145–158 (2014).
27
MC Whalen, RW Innes, AF Bent, BJ Staskawicz, Identification of Pseudomonas syringae pathogens of Arabidopsis and a bacterial locus determining avirulence on both Arabidopsis and soybean. Plant Cell 3, 49–59 (1991).
28
J Hagmann, et al., Century-scale methylome stability in a recently diverged Arabidopsis thaliana lineage. PLoS Genet 11, e1004920 (2015).
29
L Quadrana, V Colot, Plant transgenerational epigenetics. Annu Rev Genet 50, 467–491 (2016).
30
A Wibowo, et al., Hyperosmotic stress memory in Arabidopsis is mediated by distinct epigenetically labile sites in the genome and is restricted in the male germline by DNA glycosylase activity. eLife 5, e13546 (2016).
31
RH Dowen, et al., Widespread dynamic DNA methylation in response to biotic stress. Proc Natl Acad Sci USA 109, E2183–E2191 (2012).
32
OV Popova, HQ Dinh, W Aufsatz, C Jonak, The RdDM pathway is required for basal heat tolerance in Arabidopsis. Mol Plant 6, 396–410 (2013).
33
IK Greaves, M Groszmann, A Wang, WJ Peacock, ES Dennis, Inheritance of trans chromosomal methylation patterns from Arabidopsis F1 hybrids. Proc Natl Acad Sci USA 111, 2017–2022 (2014).
34
IK Greaves, et al., Trans chromosomal methylation in Arabidopsis hybrids. Proc Natl Acad Sci USA 109, 3570–3575 (2012).
35
MA Matzke, RA Mosher, RNA-directed DNA methylation: An epigenetic pathway of increasing complexity. Nat Rev Genet 15, 394–408 (2014).
36
S Cortijo, et al., Mapping the epigenetic basis of complex traits. Science 343, 1145–1148 (2014).
37
J Reinders, J Paszkowski, Unlocking the Arabidopsis epigenome. Epigenetics 4, 557–563 (2009).
38
FK Teixeira, et al., A role for RNAi in the selective correction of DNA methylation defects. Science 323, 1600–1604 (2009).
39
H Stroud, MV Greenberg, S Feng, YV Bernatavichute, SE Jacobsen, Comprehensive analysis of silencing mutants reveals complex regulation of the Arabidopsis methylome. Cell 152, 352–364 (2013).
40
L Wei, X Cao, The effect of transposable elements on phenotypic variation: Insights from plants to humans. Sci China Life Sci 59, 24–37 (2016).
41
J Zhai, et al., Small RNA-directed epigenetic natural variation in Arabidopsis thaliana. PLoS Genet 4, e1000056 (2008).
42
F Borges, RA Martienssen, The expanding world of small RNAs in plants. Nat Rev Mol Cell Biol 16, 727–741 (2015).
43
R Rajagopalan, H Vaucheret, J Trejo, DP Bartel, A diverse and evolutionarily fluid set of microRNAs in Arabidopsis thaliana. Genes Dev 20, 3407–3425 (2006).
44
A Hamaguchi, et al., A small subfamily of Arabidopsis RADIALIS-LIKE SANT/MYB genes: A link to HOOKLESS1-mediated signal transduction during early morphogenesis. Biosci Biotechnol Biochem 72, 2687–2696 (2008).
45
C Li, Y Zhou, L-M Fan, A novel repressor of floral transition, MEE3, an abiotic stress regulated protein, functions as an activator of FLC by binding to its promoter in Arabidopsis. Environ Exp Bot 113, 1–10 (2015).
46
AHC Wong, II Gottesman, A Petronis, Phenotypic differences in genetically identical organisms: The epigenetic perspective. Hum Mol Genet 14, R11–R18 (2005).
47
3rd WM Rideout, K Eggan, R Jaenisch, Nuclear cloning and epigenetic reprogramming of the genome. Science 293, 1093–1098 (2001).
48
JL Zimmerman, Somatic embryogenesis: A model for early development in higher plants. Plant Cell 5, 1411–1423 (1993).
49
F Johannes, et al., Assessing the impact of transgenerational epigenetic variation on complex traits. PLoS Genet 5, e1000530 (2009).
50
R Kooke, et al., Epigenetic basis of morphological variation and phenotypic plasticity in Arabidopsis thaliana. Plant Cell 27, 337–348 (2015).
51
T Kawakatsu, JR Nery, R Castanon, JR Ecker, Dynamic DNA methylation reconfiguration during seed development and germination. Genome Biol 18, 171 (2017).
52
D Bouyer, et al., DNA methylation dynamics during early plant life. Genome Biol 18, 179 (2017).
53
J Walker, et al., Sexual-lineage-specific DNA methylation regulates meiosis in Arabidopsis. Nat Genet 50, 130–137 (2018).
54
JP Calarco, et al., Reprogramming of DNA methylation in pollen guides epigenetic inheritance via small RNA. Cell 151, 194–205 (2012).
55
P-H Hsieh, et al., Arabidopsis male sexual lineage exhibits more robust maintenance of CG methylation than somatic tissues. Proc Natl Acad Sci USA 113, 15132–15137 (2016).
56
K Chwialkowska, U Nowakowska, A Mroziewicz, I Szarejko, M Kwasniewski, Water-deficiency conditions differently modulate the methylome of roots and leaves in barley (Hordeum vulgare L.). J Exp Bot 67, 1109–1121 (2016).
57
LJ Ferreira, V Azevedo, J Maroco, MM Oliveira, AP Santos, Salt-tolerant and -sensitive rice varieties display differential methylome flexibility under salt stress. PLoS One 10, e0124060 (2015).
58
TN Le, et al., DNA demethylases target promoter transposable elements to positively regulate stress responsive genes in Arabidopsis. Genome Biol 15, 458 (2014).
59
C Miguel, L Marum, An epigenetic view of plant cells cultured in vitro: Somaclonal variation and beyond. J Exp Bot 62, 3713–3725 (2011).
60
D Zhang, et al., Tissue culture-induced heritable genomic variation in rice, and their phenotypic implications. PLoS One 9, e96879 (2014).
61
Y Piquot, et al., Variation in sexual and asexual reproduction among young and old populations of the perennial macrophyte Sparganium erectum. Oikos 82, 139–148 (1998).

Information & Authors

Information

Published in

The cover image for PNAS Vol.115; No.39
Proceedings of the National Academy of Sciences
Vol. 115 | No. 39
September 25, 2018
PubMed: 30201727

Classifications

Data Availability

Data deposition: Sequencing reads have been deposited at the European Nucleotide Archive (accession nos. PRJEB26932 and PRJEB14117). DNA methylation data have been uploaded to the EPIC-CoGe epigenome browser and can be accessed at https://genomevolution.org/CoGe/NotebookView.pl?lid=835.

Submission history

Published online: September 10, 2018
Published in issue: September 25, 2018

Keywords

  1. epigenetics
  2. Arabidopsis thaliana
  3. DNA methylation
  4. transgenerational inheritance
  5. asexual reproduction

Acknowledgments

We thank C. Lanz and J. Hildebrandt for help with Illumina sequencing, J. Engelhorn for assistance with RNA sampling, Rainer Borriss for providing the B. amyloliquefaciens FZB42 strain, and Liliana M. Costa for discussions and comments on the manuscript. This work was supported by a European Research Council (ERC) Marie Sklodowska-Curie Fellowship (751204-H2020-MSCA-IF-2016, to A.W.); the ERC AdG IMMUNEMESIS Project, the Deutsche Forschungsgemeinschaft SPP1529 Program, and the Max Planck Society (D.W.); the ERC AdG ROOTMICROBIOTA Project, the Cluster of Excellence on Plant Sciences, and the Max Planck Society (P.S.-L.); and the Biotechnology and Biological Sciences Research Council (Grants BB/L025892/1, to G.D.B., and BB/L003023/1, BB/N005279/1, BB/N00194X/1, and BB/P02601X/1, to J.G.-M.).

Authors

Affiliations

Anjar Wibowo1
School of Life Sciences, University of Warwick, CV4 7AL Coventry, United Kingdom;
Department of Molecular Biology, Max Planck Institute for Developmental Biology, 72076 Tübingen, Germany;
Claude Becker1
Department of Molecular Biology, Max Planck Institute for Developmental Biology, 72076 Tübingen, Germany;
Gregor Mendel Institute of Molecular Plant Biology, Austrian Academy of Sciences, Vienna Biocenter, 1030 Vienna, Austria;
Present address: Gregor Mendel Institute of Molecular Plant Biology, Austrian Academy of Sciences, Vienna Biocenter, 1030 Vienna, Austria.
Julius Durr
School of Life Sciences, University of Warwick, CV4 7AL Coventry, United Kingdom;
Jonathan Price
School of Life Sciences, University of Warwick, CV4 7AL Coventry, United Kingdom;
Stijn Spaepen
Department of Plant Microbe Interactions, Max Planck Institute for Plant Breeding Research, 50829 Cologne, Germany;
Cluster of Excellence on Plant Sciences, Max Planck Institute for Plant Breeding Research, 50829 Cologne, Germany;
Sally Hilton
School of Life Sciences, University of Warwick, CV4 7AL Coventry, United Kingdom;
Hadi Putra
School of Life Sciences, University of Warwick, CV4 7AL Coventry, United Kingdom;
Ranjith Papareddy
School of Life Sciences, University of Warwick, CV4 7AL Coventry, United Kingdom;
Present address: Gregor Mendel Institute of Molecular Plant Biology, Austrian Academy of Sciences, Vienna Biocenter, 1030 Vienna, Austria.
Quentin Saintain
School of Life Sciences, University of Warwick, CV4 7AL Coventry, United Kingdom;
Sarah Harvey
School of Life Sciences, University of Warwick, CV4 7AL Coventry, United Kingdom;
Department of Biology, University of York, YO10 5DD York, United Kingdom
Gary D. Bending
School of Life Sciences, University of Warwick, CV4 7AL Coventry, United Kingdom;
Paul Schulze-Lefert
Department of Plant Microbe Interactions, Max Planck Institute for Plant Breeding Research, 50829 Cologne, Germany;
Cluster of Excellence on Plant Sciences, Max Planck Institute for Plant Breeding Research, 50829 Cologne, Germany;
Department of Molecular Biology, Max Planck Institute for Developmental Biology, 72076 Tübingen, Germany;
Jose Gutierrez-Marcos3 [email protected]
School of Life Sciences, University of Warwick, CV4 7AL Coventry, United Kingdom;

Notes

3
To whom correspondence may be addressed. Email: [email protected] or [email protected].
Author contributions: A.W., C.B., G.D.B., P.S.-L., D.W., and J.G.-M. designed research; A.W., C.B., J.D., H.P., R.P., Q.S., and S. Harvey performed research; A.W., C.B., J.P., S.S., S. Hilton, D.W., and J.G.-M. analyzed data; and A.W., C.B., D.W., and J.G.-M. wrote the paper.
Reviewers: X.C., Institute of Genetics and Developmental Biology, Chinese Academy of Sciences; and M.G., Massachusetts Institute of Technology.
1
A.W. and C.B. contributed equally to this work.

Competing Interests

A.W., C.B., D.W., and J.G.-M. are inventors on patent application “Stable epigenetic plant variants” (PCT/EP2016/055377 and WO2016146552A1) filed by University of Warwick and Max Planck Society.

Metrics & Citations

Metrics

Note: The article usage is presented with a three- to four-day delay and will update daily once available. Due to ths delay, usage data will not appear immediately following publication. Citation information is sourced from Crossref Cited-by service.


Citation statements




Altmetrics

Citations

Export the article citation data by selecting a format from the list below and clicking Export.

Cited by

    Loading...

    View Options

    View options

    PDF format

    Download this article as a PDF file

    DOWNLOAD PDF

    Login options

    Check if you have access through your login credentials or your institution to get full access on this article.

    Personal login Institutional Login

    Recommend to a librarian

    Recommend PNAS to a Librarian

    Purchase options

    Purchase this article to access the full text.

    Single Article Purchase

    Partial maintenance of organ-specific epigenetic marks during plant asexual reproduction leads to heritable phenotypic variation
    Proceedings of the National Academy of Sciences
    • Vol. 115
    • No. 39
    • pp. 9641-E9260

    Media

    Figures

    Tables

    Other

    Share

    Share

    Share article link

    Share on social media