PLoS Genetics
Public Library of Science
RNA-directed DNA methylation prevents rapid and heritable reversal of transposon silencing under heat stress in Zea mays

DOI: 10.1371/journal.pgen.1009326, Volume: 17, Issue: 6
Article Type: research-article, Article History

Most plant genomes are mostly transposable elements (TEs), most of which are held in check by modifications of both DNA and histones. The bulk of silenced TEs are associated with methylated DNA and histone H3 lysine 9 dimethylation (H3K9me2). In contrast, epigenetically silenced genes are often associated with histone lysine 27 trimethylation (H3K27me3). Although stress can affect each of these modifications, plants are generally competent to rapidly reset them following that stress. Here we demonstrate that although DNA methylation is not required to maintain silencing of the MuDR element, it is essential for preventing heat-induced, stable and heritable changes in both H3K9me2 and H3K27me3 at this element, and for concomitant changes in transcriptional activity. These finding suggest that RdDM acts to buffer the effects of heat on silenced transposable elements, and that a loss of DNA methylation under conditions of stress can have profound and long-lasting effects on epigenetic silencing in maize.


Transposable elements (TEs) are a ubiquitous feature of all genomes. They survive in large measure because they can out-replicate the rest of the genome [1]. As a consequence of that replication TEs can threaten the integrity of the host genome. In response to this threat, all forms of life have evolved mechanisms by which TEs can be silenced when they are recognized as such and, importantly, maintained in a silenced state over long periods of time, even when the initial trigger for silencing is no longer present [24]. Because plant genomes are largely composed of TEs, the majority of plant DNA is maintained in an epigenetically silent state [5]. Because they are the primary target of epigenetic silencing in plants, TEs are an excellent model for understanding the means by which particular DNA sequences are targeted for silencing, and for understanding the means by which silencing can be maintained from one generation to the next [6]. Finally, because TEs have proved to be exquisitely sensitive to a variety of stresses [79], they can also teach us a great deal about the relationship between stress and epigenetically encoded memory of stress response.

In plants, heritable epigenetic silencing of TEs is almost invariably associated with DNA methylation [1012]. The vast bulk of TEs in plant genomes are methylated and, with some notable exceptions [13], epigenetically silenced [14,15]. DNA methylation has a number of features that makes it an appealing mechanism by which silencing can be heritably propagated, either following cell divisions during somatic development, or transgenerationally, from one generation to the next. Because methylation in both the CG and CHG sequence contexts (where H = A, T or G) are symmetrical, information concerning prior DNA methylation can be easily propagated by methylating newly synthesized DNA strands using the parent strand as a template. For CG methylation, this is achieved by reading the methylated cytosine using VARIANT IN METHYLATION 1–3 (VIM1-3) [16,17] and writing new DNA methylation using the methyl transferase MET1 [1820]. For CHG, methylation is read indirectly by recognition of H3K9 dimethylation (H3K9me2) by CMT3, which catalyzes methylation of newly synthesized DNA, which in turn triggers methylation of H3K9 [2123].

Maintenance methylation of most CHH involves RNA-directed DNA methylation (RdDM). The primary signal for de novo methylation of newly synthesized DNA from previously methylated DNA sequences is thought to be transcription by RNA POLYMERASE IV (POLIV) of short transcripts from previously methylated templates [2426]. This results in the production of small RNAs that are tethered to the target DNA by RNA POLYMERASE V (POLV), which is targeted by SU(VAR)3-9 homologs SUVH2 and SUVH9, which bind to methylated DNA [27]. This in turn triggers de novo methylation of newly synthesized DNA strands using the methyl transferases DRMT1/2 [28,29]. In addition to the RdDM pathway, CHH methylation can also be maintained due to the activity of CHROMOMETHYLASE2 (CMT2), which, similar to CMT3, works in conjunction with H3K9me2 to methylate non-CG cytosines, particularly in deeply heterochromatic regions of the genome [30]. Finally, because both histones and DNA must be accessible in order to be modified, chromatin remodelers such as DDM1 are also often required for successful maintenance of TE silencing [23,31]. In plants, effective silencing of TEs requires coordination between DNA methylation and histone modifications [32]. Together, these pathways can in large part explain heritable propagation of both DNA methylation and histone modification of TEs.

In large genomes such as that of maize, much of RdDM activity is focused not on deeply silenced heterochromatin, which is often concentrated in pericentromeric regions, but on regions immediately adjacent to genes, referred to as “CHH islands” because genes in maize are often immediately adjacent to silenced TEs [15,33]. In maize, mutations in components of the RdDM pathway affect both paramutation and transposon silencing [34]. Mutations in Mediator of Paramutation 1 (Mop1), a homolog of RNA DEPENDENT RNA POLYMERASE2 (RDR2) , result in the loss of nearly all 24 nucleotide small RNAs, as well as the CHH methylation that is associated with them [3537]. Despite this, mop1 has only minimal effects on gene expression in any tissue except the meristem [33,38], and the plants are largely phenotypically normal [39]. This, along with similar observations in Arabidopsis, has led to the suggestion that the primary role of RdDM is to reinforce boundaries between genes and adjacent TEs, rather than to regulate gene expression [33]. However, it should be noted that the mop1 mutation can in some cases have effects on plant phenotype [40]. Further, mop1 mutants can enhance the effects of exogenously applied ABA [41] and mutants of Required to maintain repression6 (Rmr6 ), a homolog of the PolIV subunit DNA-directed RNA polymerase IV subunit 1 (NRPD1) [42], are altered in their response to drought, suggesting that the RdDM pathway may play a role in buffering stress responses in maize [43,44]. Further, even in wild-type backgrounds, there is evidence that the process of heritable paramutation of an allele of R1 , which is known to be dependent on RdDM, is sensitive to changes in temperature and light during specific stages of development [45].

Unlike animals, plants do not experience a global wave of DNA demethylation either in the germinal cells of the gametophyte or in the early embryo [46]. Thus, DNA methylation and associated histone modifications are an attractive mechanism for transgenerationally propagated silencing. Indeed, there is strong evidence that mutants that trigger a global loss of methylation can cause heritable reactivation of previously silenced TEs, although it is worth noting that even in mutants in which the vast majority of DNA methylation has been lost, only a subset of TEs are transcriptionally reactivated [47,48], and DNA methylation of many TEs can be rapidly reestablished at many loci via RdDM in wild-type progenies of mutant plants, suggesting that memory propagated via DNA methylation can be restored due to the presence of small RNAs that can trigger de novo methylation of previously methylated sequences [49,50].

In contrast to TEs, most genes that are silenced during somatic development in plants are associated with H3K27 trimethylation (H3K27me3), which requires the activity of the polycomb complexes PRC2 and PRC1, which together catalyze H3K27 methylation and facilitate its heritable propagation [5153]. In plants, H3K27me3 enrichment is generally associated with genes rather than TEs [54,55], and numerous developmental pathways require the proper deposition and maintenance of this modification [56,57]. The most well explored example of this involves epigenetic setting of FLOWERING LOCUS C (FLC ), a negative regulator of flowering in Arabidopsis [58,59]. In a process known as vernalization, prolonged exposure to cold results in somatically heritable silencing of this gene, which in turn results in flowering under favorable conditions in the spring. Somatically heritable silencing of FLC is initially triggered by non-coding RNAs, which are involved in recruitment of components of PRC2, which catalyze H3K27me3, which in turn mediates a somatically heritable silent state [58]. Importantly, H3K27me3 at genes like FLC is erased each generation, both in pollen and in the early embryo [6062]. The fact that H3K27me3 must be actively reset suggests that in the absence of this resetting, H3K27me3 in plants is competent to mediate transgenerational silencing but is normally prevented from doing so.

Dramatic differences in TE content between even closely related plant species suggest that despite the relative stability of TE silencing under laboratory conditions, TEs frequently escape silencing and proliferate in natural settings [63]. Stress, both biotic and abiotic can often trigger TE transcription and, at least in some cases, transposition [7,6467]. Further, there is evidence that the association of TEs and genes can result in de novo stress induction of adjacent genes [64,68,69].

Because of its dramatic and global effects on both gene expression and protein stability, heat stress has attracted considerable attention, particularly with respect to heritable transmission of TE activity. Although heat stress can trigger somatically heritable changes in gene expression, there appear to be a variety of mechanisms to prevent or gradually ameliorate transgenerational transmission of those changes [70,71]. Thus, for instance, although the ONSEN retrotransposon is sensitive to heat, it is only in mutants in the RdDM pathway that transposed elements are transmitted to the next generation [9,72]. Given that various components of regulatory pathways that have evolved to regulate TEs are up-regulated in germinal lineages, it is not surprising that a defect in one of these pathways would lead to an enhancement in the number of germinally transmitted new insertions [73,74]. The observation that it is the combination of both heat and components of the RdDM pathway results in reactivation of TEs, rather than each by itself has led to the suggestion that a key role of RdDM is to prevent TE activation specifically under conditions of stress [9,75].

Similar experiments using silenced transgenes have demonstrated that double mutants of mom1 and ddm1 cause silenced transgenes as well as several TEs to be highly responsive to heat stress, and the observed reversal of silencing can be passed on to a subsequent generation, but only in mutant progeny [76]. It is also worth noting that in many cases of TE reactivation, silencing is rapidly re-established in wild-type progeny [77,78]. The degree to which this is the case likely depends on a variety of factors, from the copy number of a given element, its position within the genome, its mode of transposition and the presence or absence of trans-acting small RNAs targeting that TE [79].

Our model for epigenetic silencing is the Mutator system of transposons in maize. The Mutator system is a family of related elements that share similar, 200 bp terminal inverted repeats (TIRs) but that contain distinct internal sequences. Nonautonomous Mu elements can only transpose in the presence of the autonomous element, MuDR. MuDR is a member of the MULE superfamily of Class II cut and paste transposons [80,81]. In addition to being required for transposition, the 200 bp TIRs within MuDR elements serve as promoters for the two genes encoded by MuDR, mudrA, which encodes a transposase, and mudrB, which encodes a novel protein that is required for Mu element integration. Both genes are expressed at high levels in rapidly dividing cells, and expression of both of them is required for full activity of the Mutator system [82,83]. MURA, the protein produced by mudrA, is sufficient for somatic excision of Mu elements, which results in characteristically small revertant sectors in somatic tissue. MuDR elements can be heritably silenced when they are in the presence of Mu killer (Muk), a rearranged variant of MuDR whose transcript forms a hairpin that is processed into 21–22 nt small RNAs that directly trigger transcriptional gene silencing (TGS) of mudrA and indirectly trigger silencing of mudrB when it is in trans to mudrA [4,84]. Because Muk can be used to heritably silence MuDR through a simple cross, and because silencing of MuDR can be stably maintained after Muk is segregated away, the MuDR/Muk system is an excellent model for understanding both initiation and maintenance of silencing. Prior to exposure to Muk, MuDR is fully active and is not prone to spontaneous silencing [85]. After exposure, MuDR silencing is exceptionally stable over multiple generations [84].

When mudrA is silenced, DNA methylation in all three sequence contexts accumulates within the 5’ end of the TIR immediately adjacent to mudrA (TIRA) [86]. Methylation at the 5’ and 3’ portions of this TIR have distinctive causes and consequences. The 5’ end of the TIR is readily methylated in the absence of the transposase, but this methylation does not induce transcriptional silencing of mudrA [87]. Methylation in this end of TIRA is readily eliminated in the presence of functional transposase. However, the loss of methylation in a silenced element in this part of the TIRA does not result in heritable reactivation of a silenced element. In contrast, CG and CHG methylation in the 3’ portion of TIRA, which corresponds to the mudrA transcript as well as to Muk-derived 22 nt small RNAs that trigger silencing, is not eliminated in the presence of active transposase and is specifically associated with heritable transcriptional silencing of mudrA.

The second gene encoded by MuDR elements, mudrB, is also silenced by Muk, but the trajectory of silencing of this gene is entirely distinct, despite the fact that the Muk hairpin has near sequence identity to the TIR adjacent to mudrB (TIRB) [4,84]. By the immature ear stage of growth in F1 plants that carry both MuDR and Muk, mudrA is transcriptionally silenced and densely methylated. In contrast, mudrB in intact elements remains transcriptionally active in this tissue, but its transcript is not polyadenylated. It is only in the next generation that steady state levels of transcript become undetectable. Further, experiments using deletion derivatives of MuDR that carry only mudrB are not silenced by Muk when they are on their own, or when they are in trans to an intact MuDR element that is being silenced by Muk. This suggests that heritable silencing of mudrB is triggered by the small RNAs that target mudrA, but the means by which this occurs is indirect and involves spreading of silencing information from mudrA to mudrB.

Silencing of mudrA can be destabilized by the mop1 mutant. MOP1 is homolog of RDR2 that is required for the production of the vast bulk of 24 nt small RNAs in maize, including those targeting Mu TIRs [3537,88]. However, silencing of MuDR by Muk is unimpeded in a mop1 mutant background, likely because Muk-derived small RNAs are not dependent on mop1 [89]. Further, although reversal of silencing of MuDR in a mop1 mutant background does occur, it only occurs gradually, over multiple generations, and only affects mudrA. In contrast, mudrB is not reactivated in this mutant background and, because mudrB is required for insertional activity, although these reactivated elements can excise during somatic development, they cannot insert into new positions.


DNA methylation is not required to maintain silencing of MuDR elements in mop1 mutants

Given that MuDR elements are only activated after multiple generations in a mop1 mutant background, we wanted to understand how silencing of MuDR is maintained in mop1 mutants prior to reactivation. To do this, we examined expression and DNA methylation at TIRA by performing bisulfite sequencing of TIRA of individuals in families that were segregating for a single silenced MuDR element, designated MuDR*, and that were homozygous or heterozygous for mop1 (S1 Fig).

In control plants carrying an active MuDR element, all cytosines in TIRA were unmethylated, which was consistent with our previous results (Fig 1B). Also consistent with previous results, F2 MuDR*/-; mop1/+ plants, whose F1 parent carried both MuDR and Muk, exhibited dense methylation at TIRA. In contrast, DNA methylation in the CG, CHH and CHG contexts at TIRA was absent in mop1 mutant siblings. Interestingly, mop1 had effects on TIRB that are more consistent with the known effects of this mutant specifically on CHH methylation. While F2 MuDR*/-; mop1/+ plants exhibited dense methylation at TIRB in all sequence contexts, mop1 homozygous siblings exhibited a loss of methylation only in the CHH context. Despite the effects of mop1 on MuDR methylation at both TIRA and TIRB, qRT-PCR results demonstrated that these mop1 mutant plants did not exhibit reactivation of mudrA or mudrB (Fig 1A).

DNA methylation patterns at TIRA and TIRB of stably silenced F2 plants.
Fig 1

DNA methylation patterns at TIRA and TIRB of stably silenced F2 plants.

(A) qPCR analysis of mudrA and mudrB expression from MuDR*/-; mop1/+ and MuDR*/-;mop1/mop1 plants. MuDR: active element. MuDR*: inactive element. Tub2 is used as an internal control gene. Six biological replicates are used for each experiment; two of six biological replicates are pooled together for each amplification. Error bars indicate mean ± standard deviation (SD) of three individuals. (B) DNA methylation patterns at TIRA and TIRB. Ten individual clones were sequenced from amplification of bisulfite-treated samples of the indicated genotypes. The cytosines in different sequence contexts are represented by different colors (red, CG; blue, CHG; green, CHH, where H = A, C, or T). For each genotype, DNA from six biological replicates were pooled.

A loss of MOP1 enhances enrichment of H3K9 and H3K27 dimethylation at TIRA

Transposon silencing is often associated with H3K9 and H3K27 dimethylation, two hallmarks of transcriptional silencing in plants [21,55]. DNA methylation, particularly in the CHG context, is linked with H3K9me2 through a self-reinforcing loop, and these two epigenetic marks often colocalize at TEs and associated nearby genes [90]. We had previously demonstrated that these two repressive histone modifications corresponded well with DNA methylation of silenced MuDR elements at TIRA [86]. However, our observation that silencing of mudrA can be maintained in the absence of DNA methylation in mop1 mutants suggests that additional repressive histone modifications may be responsible for maintaining the silenced state of mudrA. To test this hypothesis, we examined the enrichment of H3K9me2 at TIRA in individuals in a family that segregated for silenced MuDR and for mop1 homozygotes and heterozygotes by performing a chromatin immunoprecipitation quantitative PCR (ChIP-qPCR) assay. As controls, we also examined these two histone modifications in leaf tissue from plants carrying active and deeply silenced MuDR elements in a wild-type background. Compared with active MuDR/-; +/+ plants, H3K9me2 and levels were significantly enriched at TIRA in the MuDR* /-; +/+ plants (Fig 2A). The same was true of H3K27me2 (S2 Fig). Surprisingly, a significant increase in H3K9me2 and H3K27me2 at TIRA was observed in mop1 mutants compared with their mop1 heterozygous siblings and with the silenced MuDR*/-; +/+ control plants, suggesting that the loss of DNA methylation that resulted from the loss of MOP1 in these mutants actually resulted in an increase in both of these repressive chromatin marks.

ChIP-qPCR analysis of enrichment of histone marks H3K9me2 and H3K27me3 at TIRA and TIRB in mop1 mutants.
Fig 2

ChIP-qPCR analysis of enrichment of histone marks H3K9me2 and H3K27me3 at TIRA and TIRB in mop1 mutants.

ChIP-qPCR analysis of enrichment of histone marks, H3K9me2 and H3K27me3 at TIRA and TIRB. (A) Relative enrichment of H3K9me2 and H3K27me3 in leaf 3 of plants of the indicated genotypes. MuDR: active element. MuDR*: inactive element. (B) Relative enrichment of H3K9me2 and H3K27me3 in leaf 3 of plants of the indicated genotypes. qPCR signal was normalized to Copia and then to the value of input sample. An unpaired t-test was performed. Error bars indicate mean ± standard deviation (SD) of the three biological replicates. *P<0.05; **P < 0.01

Silencing of TIRB is associated with an increase in H3K27me3

Like mudrA, mudrB is silenced by Muk, but maintenance of mudrB silencing has distinct requirements. Unlike mudrA, which is eventually reactivated in a mop1 mutant background under normal conditions, mudrB remains silenced, suggesting that maintenance of silencing of this gene is independent of MOP1 [36]. ChIP-qPCR revealed that silencing of mudrB is not associated with H3K9me2 methylation. Instead, heritably silenced TIRB is enriched for H3K27me3, a modification normally associated with somatically silenced genes rather than transposable elements (Fig 2B). The mop1 mutant appears to enhance H3K27me3 at TIRB relative to the mop1 heterozygous siblings, although the enrichment is no greater that observed in the MuDR*/-; +/+ controls.

Application of heat stress specifically in the early stage of growth can promote the reactivation of silenced MuDR elements in mop1 mutants

There is ample evidence that a variety of stresses can reactivate epigenetically silenced TEs. One particularly effective treatment is heat stress. Given that a loss of methylation by itself is not sufficient to reactivate silenced MuDR elements, we subjected mop1 mutant and mop1 heterozygous sibling seedlings carrying silenced MuDR elements (MuDR*) to heat stress. Fourteen-day-old MuDR*/-; mop1/mop1 and MuDR*/-; mop1/+ sibling seedlings were heated at 42°C for four hours and leaf samples were collected immediately after that treatment (Fig 3A). qRT-PCR for the heat response factor Hsp90 (Zm00001d024903) confirmed that the seedlings were responding to the heat treatment (S3 Fig). We then examined MuDR transcription by performing qRT-PCR on RNA from leaf three immediately after the plants had been removed from heat and from control plants that had not been subjected to heat stress. In the mop1 mutants, both mudrA and mudrB became transcriptionally reactivated upon heat treatment (Fig 3B). MuDR elements in plants that were mop1 mutant that were not heat stressed and were those that were wild-type and that were heat stressed were not reactivated, demonstrating that both a mutant background and heat stress are required for efficient reactivation. To determine if the application of heat stress at a later stage of plant development can also promote reactivation, we heat-stressed 28-day-old plants and examined MuDR transcription in leaf seven at a similar stage of development (~10 cm) as had been examined in heat stressed leaf three in the previous experiment. In these plants, we saw no evidence of MuDR reactivation although qRT-PCR Hsp90 indicated that these seedlings were responding to the heat treatment (Figs 3B and S3). Taken together, these data suggest that the application of heat stress specifically at an early stage of plant development can promote the reactivation of a silenced TE in a mutant that is deficient in the RdDM pathway.

Expression of mudrA and mudrB in plants under heat stress.
Fig 3

Expression of mudrA and mudrB in plants under heat stress.

(A) Schematic diagram of the heat-reactivation experiment. (B) qRT-PCR of mudrA and mudrB in leaf 3 and leaf 7 in plants of the indicated genotypes. Twelve biological replicates are used for each experiment. Tub2 is an internal control gene. Additional controls for each experiment include MuDR/-, pooled samples from twelve heated MuDR*/-; mop1/+ plants and twelve unheated plants. Red text is used to indicate samples that were subjected to heat stress.

TIRA in a mop1 mutant background already lacks any DNA methylation prior to heat treatment and thus heat would not be expected to reduce TIRA methylation. However, in mop1 mutants TIRB retained CG and CHG methylation and also remained inactive (Fig 1B). To determine if reactivation after heat treatment is associated with a loss of this methylation, we examined DNA methylation at TIRB in mop1 mutants in the presence or absence of heat treatment. This assay was performed on the same tissues that we collected for MuDR expression reactivation analysis. We found that the DNA methylation pattern was the same for both the heat-treated and the control mop1 mutant plants, indicating that heat stress does not alter TIRB methylation and that a further loss of DNA methylation is not the cause of mudrB reactivation in this tissue (S4 Fig).

Heat stress reverses TE silencing by affecting histone modifications at TIRA and TIRB

Under normal conditions, we found that H3K9me2 at TIRA is associated with silencing, and H3K9me2 is actually enriched when TIRA methylation is lost in mop1 mutants (Fig 2A). In contrast, we find that H3K27me3, rather than H3K9me2, is enriched at TIRB and is maintained at similar or slightly elevated levels in mop1 mutant relative to mop1 heterozygous siblings (Fig 2B). Given these observations, we hypothesized that heat stress may reverse H3K9me2 enrichment at TIRA and H3K27me3 enrichment at TIRB. To test this hypothesis, we determined the level of H3K9me2 and H3K27me3 at TIRA and TIRB under normal and stressed conditions using ChIP-qPCR.

Upon heat stress, the level of H3K9me2 at TIRA was significantly decreased in mop1 mutants compared to that of non-treated mop1/mop1 mutant siblings (Fig 4A). Interestingly, however, H3K9me2 enrichment only decreased to the level observed at TIRA in silenced MuDR*/-; +/+ plants, and it remained significantly higher than that of TIRA in the naturally active MuDR/-; +/+ plants. In contrast, we observed no changes in H3K27me3 at TIRA.

ChIP-qPCR analysis of histone marks TIRA and TIRB under heat stress.
Fig 4

ChIP-qPCR analysis of histone marks TIRA and TIRB under heat stress.

Relative enrichment of H3K9me2 and H3K27me3 at TIRA (A) and TIRB (B) in leaf 3 of plants of the indicated genotypes. (Relative enrichment of H3K4me3 at TIRA (C) and TIRB (D) in leaf 3 of plants of the indicated genotypes. qPCR signals were normalized to Copia and then to the value of input samples. MuDR* refers to a silenced MuDR element. MuDR~ refers to a reactivated element. Red text indicates a sample that has been heat-treated. Error bars indicate mean ± standard deviation (SD) of the three biological replicates. **P < 0.01; ***P < 0.001

At TIRB, we observed no changes in H3K9me2 enrichment in any of our samples. Instead, we found that heat treatment reversed previously established H3K27me3 at TIRB, supporting the hypothesis that this modification, rather than H3K9me2, mediates heritable silencing of mudrB (Fig 4B). Consistent with evidence for transcriptional activation of both mudrA and mudrB , we observed enrichment of the active mark histone H3 lysine 4 trimethylation (H3K4me3) in reactivated TIRA and TIRB (Fig 4C and 4D). Taken together, these data demonstrate that heat stress can simultaneously reduce two often mutually exclusive repressive histone modifications, H3K9me2 and H3K27me3 at the two ends of a single TE.

The reactivation state is somatically transmitted to the new emerging tissues

We next sought to determine whether or not the reactivated state can be propagated to cells in somatic tissues after the heat had been removed. We performed quantitative RT-PCR to detect mudrA and mudrB transcripts in mature leaf ten of plants 35 days after the heat stress and in immature tassels ten days after that. At V2, when the heat stress was applied and leaf three was assayed, cells within leaf 10 primordia are present and may have experienced the heat stress. In contrast, because the tassel primordia are not formed until V5, the cells of the tassel could not have experienced the heat stress directly [91,92]. We found that both genes stayed active in both tissues, indicating heat-induced reactivation is stably transmitted to new emerging cells and tissues (Fig 5).

Expression of mudrA and mudrB in new emerging tissues following heat stress.
Fig 5

Expression of mudrA and mudrB in new emerging tissues following heat stress.

(A) Diagram of the experiment. (B) qPCR was performed to measure transcript levels of mudrA and mudrB using expression of Tub2 as an internal control. Expression levels were normalized to that of an active MuDR element. Error bars indicate mean ± standard deviation (SD) of the ten biological replicates.

MuDR activity is stably transmitted to subsequent generations

Our previous work had demonstrated that silenced mudrA (but not mudrB) can be progressively and heritably reactivated only after multiple generations of exposure to the mop1 mutation under normal conditions. Only after eight generations could this activity be stably transmitted to subsequent generations in the absence of the mop1 mutation [36]. To determine if the somatic activity we observed after heat stress can be transmitted to the next generation, we crossed the heat-treated mop1 homozygous plants that carried transcriptionally reactivated MuDR (designated MuDR~) and the sibling mop1 homozygous MuDR* control plants, to a tester that was homozygous wild-type for mop1 and that lacked MuDR (Fig 6A). MURA, the protein encoded by mudrA causes excision of a reporter element at the a1-mum2 allele of the A1 gene, resulting pale kernels with spots of colored revertant tissue. All plants used in these experiments were homozygous for a1-mum2. If mudrA were fully heritably reactivated, a cross between a MuDR~/-; mop1/mop1 plant and a tester would be expected to give rise to 50% spotted kernels, and this phenotype would be expected to cosegregate with the reactivated MuDR element (Fig 6A). The progeny of ten independent heat-reactivated individuals gave a total of 45% spotted kernels. In contrast, ten mop1 homozygous siblings that carried MuDR* and that had not been heat-treated gave rise to an average of only 0.7% spotted kernels after test crossing (Fig 6B and S2 Table). These results show that MuDR activity induced by heat treatment was transmitted to the next generation. qRT-PCR in both endosperms and embryos of the spotted and pale progeny kernels and genotyping for the presence or absence of MuDR at position 1 on chromosome 9L [85] demonstrated that activity was transmitted to both the embryo and the endosperm, and that this activity cosegregated with the single MuDR present in these families (S5 Fig). We employed a similar strategy to test stability of heritability (Fig 6C). We crossed three subsequent generations to testers and counted the spotted kernels. We observed that the progeny of heat-reactivated individuals gave a total of 51%, 48% and 47% spotted kernels in the three subsequent generations. In contrast, subsequent generations of the lineage carrying MuDR* that had not been heat-treated gave rise to only a small number of weakly spotted kernels (Fig 6D, S2 Table). These results demonstrate that heat reactivation is stable over multiple generations in a non-mutant genetic background, as is silencing in the absence of heat stress.

Testing transgenerational inheritance.
Fig 6

Testing transgenerational inheritance.

(A) A schematic diagram showing the crosses used to determine transgenerational inheritance. (B) Ears derived from heat-treated and control individuals. (C) Crosses done in the generations following heat stress. (D) Ratios of spotted kernels in subsequent generations following the heat stress (H1) generation. Red text indicates a sample that has been heat-reactivated.

DNA hypomethylation is not associated with transgenerational inheritance of activity

We have shown that DNA methylation is not reduced under heat stress at TIRB, and that even a complete absence of methylation of TIRA under normal conditions does not result in transcriptional activation. These results suggest that, at least under normal conditions, DNA methylation of MuDR is neither necessary nor sufficient to mediate silencing. However, only plants that were mop1 mutant and whose TIRs were missing either methylation of cytosines in all sequence contexts in the case of TIRA or those in the CHH sequence context in the case of TIRB were reactivated under heat stress. This suggests that a loss of methylation may be a precondition for initiation, and perhaps propagation, of continued activity after that stress. To test the latter possibility, we examined DNA methylation at TIRA and TIRB in the mop1 heterozygous H2 progenies of heat-reactivated mop1 mutant plants and those of their unheated mop1 mutant sibling controls (S1 Fig). Surprisingly, we found that both TIRA and TIRB were extensively methylated in all three sequence contexts in all progenies examined regardless of their activity status (Fig 7). Indeed, their methylation was indistinguishable from that observed at silenced MuDR elements. This suggests that although the restoration of MOP1 function does result in the restoration of methylation at both TIRA and TIRB in these heritably reactivated MuDR elements, this methylation is not sufficient for reestablishment of silencing at either of these TIRs. In order to determine whether DNA methylation we observed in these wild-type H2 plants was stable, we examined TIRA and TIRB methylation in plants three and four generations removed from the initial heat stress. Surprisingly, we found that the observed patterns of methylation in this generation at both TIRs closely resembled that of fully active MuDR elements (Fig 7). This suggests that patterns of methylation consistent with activity are in fact restored in the heat stressed lineage after MOP1 function is restored, but only after changes in H3K9me2 at TIRA and in H3K27me3 at TIRB have already taken place.

DNA methylation patterns at TIRA and TIRB in H2, H4 and H5 progeny of heat-treated plants.
Fig 7

DNA methylation patterns at TIRA and TIRB in H2, H4 and H5 progeny of heat-treated plants.

(A) DNA methylation patterns at TIRA. (B) DNA methylation patterns at TIRB. Ten individual clones were sequenced from each amplification of bisulfite-treated sample. The cytosines in different sequence contexts are represented by different colors (red, CG; blue, CHG; green, CHH, where H = A, C, or T). Red text indicates plants derived from heat-treated plants. MuDR* refers to a silenced MuDR element. MuDR~ refers to a reactivated element. For each assay, six independent samples were pooled together.

Transgenerational heritability of activity is associated with heritability of histone modifications

DNA hypomethylation is not associated with transgenerational inheritance of MuDR activity, and DNA hypermethylation does not result in a restoration of silencing in wild-type progeny of heat reactivated mutants. A plausible alternative is that the observed changes in histone marks mediate heritable propagation of activity of both mudrA and mudrB independent of methylation status. To test this hypothesis, we determined the levels of H3K9me2, H3K27me3 and H3K4me3 at TIRA and TIRB in the mop1 heterozygous H2 progenies of heat-reactivated MuDR~/-; mop1/mop1 plants and those of their sibling untreated MuDR*/-; mop1/mop1 sibling controls. Consistent with the continued activity of mudrB in the progeny of the heat stressed plants, relative levels of H3K27me3 levels remained low and H3K4me3 remained high at TIRB in these plants, suggesting that heritable propagation of H3K27me3 is responsible for that continued activity (Fig 8). Similarly, at TIRA, H3K9me2 remained low and H3K4me3 remained high in these progenies. Interestingly, the increase in DNA methylation in these MuDR active mop1 heterozygous plants was associated with a further decrease in levels of H3K9me2 at TIRA relative to that of their heat stressed mop1 homozygous parents, down to the levels of the active MuDR control. This suggests that an increase in methylation of these active elements in the wild-type background resulted in a concomitant decrease in H3K9me2 at TIRA.

ChIP-qPCR analysis of enrichment of histone marks, H3K9me2, H3K27me3 and H3K4me3 at TIRA and TIRB.
Fig 8

ChIP-qPCR analysis of enrichment of histone marks, H3K9me2, H3K27me3 and H3K4me3 at TIRA and TIRB.

Relative enrichment of H3K9me2, H3K27me3 and H3K4me3 at TIRA and TIRB in leaf 3 of plants of the indicated genotypes. qPCR signals were normalized to Copia and then to the value of input samples. Red text indicates plants derived from heat-treated plants. MuDR* refers to a silenced MuDR element. MuDR~ refers to a reactivated element. An unpaired t-test was performed. Error bars indicate mean ± standard deviation (SD) of the three biological replicates. *P<0.05; **P < 0.01; ***P<0.001.


DNA methylation is neither necessary nor sufficient for the maintenance of silencing at TIRA or TIRB

Our results demonstrating that methylation is not necessary for maintenance of epigenetic silencing in mop1 mutant plants (Fig 1) and is not sufficient to trigger silencing in H2 reactivated plants (Fig 7) suggest that at this particular locus, DNA methylation is not the key determinative factor with respect to either silencing or its reversal. In contrast, changes in H3K9me2 are closely correlated with changes in TIRA activity, suggesting that it is this modification, rather than DNA methylation, that mediates both activity and heritable transmission of silencing of mudrA . Given that H3K9me2 is normally tightly associated with cytosine methylation, particularly in the CHG context [21,93], this result is unexpected. However, our results clearly demonstrate that this modification can be heritably propagated in the absence of DNA methylation and in the absence of the original trigger for silencing, Muk. Even more unexpected is our observation that, once mudrA becomes silenced, in mop1 mutants there appears to be reciprocal relationship between DNA methylation of TIRA and H3K9me2 enrichment. Methylation in all three-sequence-context is eliminated throughout TIRA in mop1 mutants, but this does not result in reactivation of mudrA. Instead, H3K9me2 actually significantly increases in the mop1 mutant. This suggests that silencing at this locus is maintained via a balance between DNA and histone methylation, such that a loss of DNA methylation actually triggers an increase in histone modification. This in turn suggests that the state of activity of mudrA in some way determines the balance between histone and DNA modification, since neither modification by itself appears to be determinative. Our heat experiment supports this hypothesis. Heat rapidly reduces histone modification, but only back down to the level of the silent mop1 heterozygous siblings rather that to the level of TIRA in an active element. In this case, the combination of an absence of DNA methylation with this reduced level of H3K9me2 appears to be sufficient to permit transcription of mudrA, as well as somatic propagation of the reactivated state to daughter cells after the heat is removed. Also supporting a balance hypothesis is the observation that in reactivated mop1 heterozygous progeny of mop1 homozygous heat-treated plants, methylation is restored to that observed in silenced elements and levels of H3K9me2 are then reduced to the level observed in active elements. This again suggests that levels of DNA and histone modification balance each other, such that in increase in methylation in the wild-type progeny of reactivated mop1 mutant plants results in a concomitant decrease in histone modification. Interestingly, however, at some point after leaf 3 of the H2 generation methylation levels are reduced to those of active MuDR elements, suggesting that this reduced methylation level is a consequence, rather than a cause, of maintenance of activity. Collectively, these data suggest that DNA methylation can be a lagging indicator that is responding to a given epigenetic state, rather than determining it.

There are other instances in which silencing can be reversed without a loss of methylation. For instance, mutations in the putative chromatin remodeler MOTHER OF MORPHEOUS1(MOM1 ) can result in activation of silenced transgenes and some endogenous loci in the absence of a loss of DNA methylation [9496]. Similarly, Microrchidia (MORC) ATPase genes, as well the H3K27 monomethyltransferases, TRITHORAX-RELATED PROTEIN 5 (ATXR5) and ATXR6 in Arabidopsis, are required for heterochromatin condensation and TE silencing but not for DNA methylation or histone modification associated with that silencing [9799]. However, unlike reactivated MuDR elements in our experiments, reintroduction of the wild-type MOM1 or MORC alleles result in immediate re-silencing. Finally, mutations in two closely related Arabidopsis genes, MAINTENANCE OF MERISTEMS-LIKE 1 (MAIL1) and MAINTENANCE OF MERISTEMS (MAIN ), can also result in activation of a subset of Arabidopsis TEs in the absence of a loss of methylation [100].

The RdDM pathway buffers the effects of heat stress on silenced MuDR elements

Heat stress rapidly reverses silencing and is associated with a reduction of H3K9me2, but only in a mop1 mutant background. This suggests that although DNA methylation is not required for the maintenance of silencing of mudrA and is not sufficient to trigger de novo silencing of this gene, it is required to prevent a response to heat stress. Thus, we suggest that the primary role of DNA methylation in this instance is to buffer the effects of heat. We note that this observation is similar but distinct from what has been observed for the ONSEN retrotransposon in Arabidopsis. In that case, although heat stress by itself can induce transcription of ONSEN [9,75], it is only when the RdDM pathway is deficient that new insertions are transmitted to the next generation. However, in wild-type progenies of heat stressed mutants, ONSEN elements are rapidly re-silenced [101]. In contrast, reactivated MuDR elements remain active for at least five generations, despite the fact that the RdDM pathway rapidly restores DNA methylation at both TIRA and TIRB. This is likely due to differences between these two elements with respect to the means by which the two elements are maintained in a silenced state. In the absence of Muk, MuDR elements are stably active over multiple generations [85,102]. This suggests that silencing of MuDR requires aberrant transcripts that are distinct from those produced by MuDR that are not present in the minimal Mutator line. Experiments involving some low copy number elements in Arabidopsis that are activated in the DNA methylation deficient ddm1 mutant background suggest that the same is true for these elements as well; once activated, these elements remain active even in wild-type progeny plants [103]. In contrast, evidence from other TEs suggests that transcripts from these elements or their derivatives contribute to their own silencing [47,104,105].

The effects of heat on mop1 mutants are dependent on the stage of development

Our heat stress experiments demonstrated that although heat exposure has a rapid and dramatic effect on MuDR activity in juvenile leaves, heat stress later during adult growth has no effect on this element. Expression analysis of Hsp90, a key marker of heat stress in maize, suggests that the older maize mop1 mutant plants are in fact responding to the heat, but the response does not include reactivation of MuDR . The reason for this difference is not clear. Presumably there are factors expressed later during development that can compensate for the lack of MOP1 in these later leaves. Expression analysis shows dramatic differences between juvenile and adult leaves, including differences in a large number of genes related to stress response [106]. Further, the transition from juvenile to adult growth in maize is associated with a transient loss of mudrA silencing in F1 plants carrying both MuDR and Muk , suggesting that this transition represents an important stage of development with respect to silencing pathways [86]. Future experiments will focus on mutations that affect the time of the juvenile to adult transition that are known to affect the transient loss of MuDR silencing.

Heritably transmitted silencing of TIRB is associated with H3K27me3

Our observation that transgenerationally heritable silencing of mudrB is associated with H3K27me3 was surprising, given that this mark is generally associated with somatic silencing of genes that is reset each generation [107]. However, in the absence of that resetting, silencing can be heritably transmitted to the next generation [60,62]. Our data clearly shows that this is the case for mudrB, whose H3K27me3 enrichment can be heritably transmitted following the loss of Mu killer through at least two rounds of meiosis, and we have evidence that mudrB remains stably silenced for at least eight generations [36]. Given that there is no selective pressure to reset TE silencing mediated by H3K27me3, this is not surprising. Interestingly, although most mutants that affect paramutation are, like Mop1 , components of the RdDM pathway [34], Maintain repression 12 (Rmr12 ), is a gene encoding a protein orthologous to PICKLE, a putative CHD3-type chromatin remodeling factor in Arabidopsis [108]. In that species, PICKLE can either (directly) reduce or (indirectly) enhance H3K27me3 at target genes. Given that Rmr12 is required for stable silencing of the paramutable Pl’ epiallele of the Pl1 gene, perhaps via modification of H3K27me3, it will be interesting to see whether or not Rmr12 is required for stable maintenance of silencing of mudrB.

There is evidence that heat stress can heritably reverse H3K27me3 at specific loci. H3K27 trimethylation can be reversed by the H3K27me3 demethylase RELATIVE OF EARLY FLOWERING 6 (REF6), which acts in conjunction the chromatin remodeler BRAHMA (BRM) to relax silencing at loci containing CTCTGYTY motifs [109]. In Arabidopsis, under heat stress, HEAT SHOCK TRANSCRIPTION FACTOR A2 (HSFA2) activates REF6, which can in turn de-repress HSFA2 by reducing H3K27me3 at this gene. This feedback loop can extend to the progeny of heat stressed plants, resulting in a heritable reduction in levels of H3K27me3 at REF6-targeting genes [110,111]. However, as in the case for all transgenerational shifts in gene expression, the effect is temporary, and both H3K27me3 and gene expression levels are restored to their original state after two generations.

It should also be noted that at both TIRA and TIRB, H3K4me3 is the most reliable indicator of activity. Thus, it is possible that stable maintenance of activity of both mudrA and mudrB is mediated via heritable maintenance of H3K4me3. Indeed, an active epiallele derived from callus culture in rice was associated with decrease in DNA methylation and H3K27me3 and an increase in H3K4me3. Overexpression of histone H3K4 demethylase JUMONJI C 703 (JMJ703 ) resulted in a loss of activity of the reactivated epiallele, a restoration of DNA methylation, an increase in H3K27me3 and a decrease in H3K4me3 [112]. This same H3K4 demethylase has also been shown to be required for stable silencing and DNA methylation of a number of Rice TEs [113]. These data suggest that silencing requires a stable balance between activating and inactivating chromatin marks as well as DNA methylation, changes in any one of which can result in resetting of epigenetic states.


Overall, our data suggests that even when examining a single TE in a single organism, a wide variety of epigenetic processes can be seen to play a role in both silencing and its reversal. At TIRA, a loss of DNA methylation in mop1 mutants is associated with what appears to be a compensatory increase in H3K9me2, which is heritably reversed by a brief exposure to heat. Heritable transmission of a reactivated state of mudrA is refractive to a restoration of DNA methylation, which instead appears to adjust over time to reflect that activity rather than to block it. In contrast to mudrA (and most other TE genes) heritable mudrB silencing is associated with H3K37me3 enrichment, which, like H3K9me2 enrichment at TIRA, is readily and heritably reversed by heat treatment. At both TIRA and TIRB, methylation is neither necessary nor sufficient for silencing, but a lack of MOP1 and an associated loss of DNA methylation at both TIRs does appear to be required to precondition both mudrA and mudrB for responsiveness to heat, consistent with a role for RdDM in buffering the effects of high temperature in maize. Clearly, these results are primarily phenomenological, as the precise mechanism for the reversal of silencing we observe remains a mystery. However, they do suggest that there is a great deal that we do not yet understand about how silenced states can be maintained and how they can be reversed.

Materials and methods

Plant materials

Maize seedlings and adult plants were grown in MetroMix under standard long-day greenhouse conditions at 26°C unless otherwise noted. The minimal Mutator line consists of one full-length functional MuDR element and one nonautonomous Mutator element, Mu1. Mu killer (Muk), a derivative version of the MuDR transposon, can heritably trigger epigenetic silencing of that transposon. Mutator activity is monitored in seeds via excisions of a Mu1 element inserted into the a1-mum2 allele of the A1 gene, resulting in small sectors of revertant tissue, or spots, in the kernels when activity is present. When MuDR activity is absent, the kernels are pale. All plants described in these experiments are homozygous for a1-mum2. Although MuDR can be present in multiple copies, all of the experiments described here have a single copy of MuDR at position 1 on chromosome 2L [102].

All of the crosses used to generate the materials examined in this paper are depicted in S1 Fig. Active MuDR/-; mop1/mop1 plants were crossed to Muk/-; mop1/+ plants. The resulting progeny plants were genotyped using PCR Mix (Syd Labs) to screen for plants that carried MuDR, Muk and that were homozygous for mop1, which were designated F1 plants. F1 plants were then crossed to mop1 heterozygotes. Progeny plants lacking Muk but carrying silenced MuDR elements, designated MuDR*, were designated F2 MuDR* progeny. F2 MuDR* progeny that were homozygous for mop1 were crossed to mop1 heterozygotes. The resulting F3 plants were genotyped for the presence of MuDR. These plants were either homozygous or heterozygous for mop1. These F3 plants were those that were used for the heat stress experiments. H1 refers to the first generation of these F3 plants that were subjected to heat stress, with successive generations designated H2, H3, etc. MuDR was genotyped using primers Ex1 and RLTIR2. Because Ex1 is complementary to sequences flanking MuDR in these families, this primer combination is specific to the single MuDR element segregating in these families. Muk was genotyped using primers TIRAout and 12-4R3. The mop1 mutation was genotyped using primers ZmRDR2F, ZmRDR2R and TIR6. All primer sequences are provided in S1 Table.

Tissue sampling

Plants used in all experiments were genotyped individually. The visible portion of each developing leaf blade, when it was ≈10 cm, was harvested when it emerged from the leaf whorl. Only leaf blades of mature leaves were harvested. For the heat reactivation experiment, seedlings were grown at 26°C for 14 days with a 12–12 light dark cycle. Seedlings were incubated at 42°C for 4 hours and leaf 3 was harvested immediately after stress treatment. As a control, leaf 3 was also collected from sibling seedlings grown at 26°C. For each genotype and treatment, 12 biological replicates were used, all of which were siblings. Samples were stored in -80°C. After sample collection, all seedlings were transferred to a greenhouse at 26°C. In order to determine if reactivation could be propagated to new emerging tissues, leaf 10 at a similar stage of development (~10 cm, as it emerged from the leaf whorl) and the immature tassel (~20 cm) were collected from each individual (Fig 5A). To determine if the application of heat stress at a later stage of plant development can promote reactivation, an independent set of these seedlings from the same family were used. A similar strategy was employed. However, in this case, seedlings were heat stressed for 4 hours after the plants had grown 28 days at 26°C. Leaf 7 was collected instead (Fig 3A). For the bisulfite sequencing experiment, leaf 3 was collected from each individual, when it was ≈10 cm, as it emerged from the leaf whorl. In order to minimize potential variation among different individuals, leaves from 6 individuals with the same genotype and treatment were pooled together. For the ChIP assays, a total of ~ 2 g of leaves from leaf 3 of 6 sibling plants with the indicated genotypes was harvested. Three independent sets of these sample collections were collected and analyzed for each genotype and treatment. Leaf samples were fixed with 1% methanol-free formaldehyde and then stored in -80°C.

RNA isolation and qRT-PCR analysis

Total RNA was extracted using TRIzol reagent (Invitrogen) and purified by Zymo Direct-zol RNA Miniprep Plus kit. 2 μl of total RNA was first loaded on a 1% agarose gel to check for good quality. Then, RNA was quantified by a NanoDrop spectrophotometer (Thermo Fisher Scientific) and reverse transcribed using an oligo-dT primer and GoScript Reverse Transcriptase (Promega). Quantitative RT-PCR was performed by using SYBR Premix Ex Taq (TaKaRa Bio) on an ABI StepOnePlus Real-Time PCR thermocycler (Thermo Fisher Scientific) according to the manufacturer’s instructions. Relative expression values for all experiments were calculated using Tub2 (Zm00001d010275) as an internal control and determined by using the comparative CT method. Sequences for all primers used for qRT-PCR are available in S1 Table and numerical data used for the figures are available in S4 Table.

Genomic bisulfite sequencing

These experiments were performed as previously described [87]. In brief, genomic DNA was isolated and digested with RNase A (Thermo Fisher Scientific). 2 μl of this DNA was loaded on a 1% agarose gel to check for good quality and then quantified using a Qubit fluorometer (Thermo Fisher Scientific). 0.5–1 μg of genomic DNA from each genotype and treatment were used for bisulfite conversion. The EZ DNA Methylation-Gold kit (Zymo Research) was used to perform this conversion. Fragments from TIRA and TIRB were PCR-amplified using EpiMark Hot Start Taq DNA Polymerase (New England BioLabs). For TIRA, the first amplification was for 20 cycles using p1bis2f and TIRAbis2R with an annealing temperature of 48°C, followed by re-amplification for 17 cycles using TIRAbis2R and TIRAmF6 with an annealing temperature of 50°C. Amplicons from TIRB were amplified for 30 cycles using methy_TIRBF and methy_TIRBR with an annealing temperature of 55°C. The resulting fragments were purified and cloned into pGEM-T Easy Vector (Promega). Ligations and transformations were performed as directed by the manufacturer’s instructions. The resulting colonies were screened for the presence of insertions by performing a colony-based PCR using primers of pGEMF and pGEMTR with an annealing temperature of 52°C. The sequences of all primers are provided in S1 Table. Plasmid was extracted from positive colonies using the Zyppy Plasmid Kit (Zymo Research). Plasmid from at least 10 independent clones were sequenced at Purdue Genomics Core Facility. The sequences were analyzed using Kismeth ( [114]. Sequences for all primers are available in S1 Table and numerical data used for figures are available in S3 Table.

Chromatin immunoprecipitation (ChIP)

The ChIP assay was performed as described previously with some modifications [115117]. Briefly, leaf samples were treated with 1% methanol-free formaldehyde for 15 minutes under vacuum. Glycine was added to a final concentration of 125 mM, and incubation was continued for 5 additional minutes. Plant tissues were then washed with distilled water and homogenized in liquid nitrogen. Nuclei were isolated and resuspended in 1 mL nuclei lysis buffer (50 Mm tris-HCl pH8, 10 mM EDTA, 0.25% SDS, protease inhibitor). 50 μl of nuclei lysis was harvested for a quality check. DNA was sheared by sonication (Bioruptor UCD-200 sonicator) sufficiently to produce 300 to 500 bp fragments. After centrifugation, the supernatants were diluted to a volume of 3 mL in dilution buffer (1.1% Triton X-100, 1.2 mM EDTA, 16.7 mM Tris-HCl pH8, 167 mM NaCl). Each sample of supernatant was sufficient to make 6 immunoprecipitation (IP) reactions. Every 500 μl sample was precleared with 25 μl protein A/G magnetic beads (Thermo Fisher Scientific) for 1 hour at 4°C. After the beads were removed using a magnet, the supernatant was removed to a new pre-chilled tube. 50 μl from each sample was used to check for sonication efficiency and set aside to serve as the 10% input control. Antibodies used were anti-H3K9me2 (Millipore), H3K27me2 (Millipore), H3K27me3 (Active Motif), H3K4me3 (Millipore) and H3KAc (Millipore). After incubation overnight with rotation at 4°C, 30 μl of protein A/G magnetic beads was added and incubation continued for 1.5 hours. The beads were then sequentially washed with 0.5 mL of the following: low salt wash buffer (20 mM Tris (pH 8), 150 mM NaCl, 0.1% (wt/vol) SDS, 1% (vol/vol) Triton X-100, 2 mM EDTA), high salt wash buffer (20 mM Tris (pH 8), 500 mM NaCl, 0.1% (wt/vol) SDS, 1% (vol/vol) Triton X-100, 2 mM EDTA), LiCl wash buffer (10 mM Tris (pH 8), 250 mM LiCl, 1% (wt/vol) sodium deoxycholate, 1% (vol/vol) NP-40 substitute, 1 mM EDTA), TE wash buffer (10 mM Tris (pH 8), 1 mM EDTA). After the final wash, the beads were collected using a magnet and resuspended with 200 μl X-ChIP elution buffer (100 mM NaHCO3, 1% (wt/vol) SDS). A total of 20 μl 5M NaCl was then added to each tube including those samples used for quality checks. Cross-links were reversed by incubation at 65°C for 6 hours. Residual protein was digested by incubating with 20 μg protease K (Thermo Fisher Scientific) at 55°C for 1 hour, followed by phenol/chloroform/isoamyl alcohol extraction and DNA precipitation. Final precipitated DNA was dissolved in 50 μl TE. Quantitative RT-PCR was performed by using SYBR Premix Ex Taq (TaKaRa Bio) on an ABI StepOnePlus Real-Time PCR thermocycler (Thermo Fisher Scientific) according to the manufacturer’s instructions. The primers used in this study are listed in S2 Table. The primers used to detect H3K9 and H3K27 dimethylation of Copia retrotransposons and H3K4 trimethylation of actin that were used as internal controls in this study have been validated previously [117]. Primers used for TIRA (TIRAR and TIRAUTRR) and TIRB (Ex1 and RLTIR2) were those used previously to detect changes in chromatin at these TIRs [86]. Expression values were normalized to the input sample that had been collected earlier using the comparative CT method. Sequences for all primers are available in S1 Table and numerical data used for figures are available in S5 Table.


We thank R. Keith Slotkin for critical reading of the manuscript and Anthony Cannon for performing genetic analysis to test the stability of transgenerational heritability.


IR Arkhipova. Neutral Theory, Transposable elements, and eukaryotic genome evolution. Mol Biol Evol. 2018;35(6):13327. Epub 2018/04/25. doi: doi: 10.1093/molbev/msy083 ; PubMed Central PMCID: PMC6455905.

A Hosaka, T Kakutani. Transposable elements, genome evolution and transgenerational epigenetic variation. Curr Opin Genet Dev. 2018;49:438. Epub 2018/03/12. doi: doi: 10.1016/j.gde.2018.02.012 .

CJ Underwood, IR Henderson, RA Martienssen. Genetic and epigenetic variation of transposable elements in Arabidopsis. Curr Opin Plant Biol. 2017;36:13541. Epub 2017/03/28. doi: doi: 10.1016/j.pbi.2017.03.002 ; PubMed Central PMCID: PMC5746046.

RK Slotkin, M Freeling, D Lisch. Heritable transposon silencing initiated by a naturally occurring transposon inverted duplication. Nat Genet. 2005;37(6):6414. Epub 2005/05/24. doi: doi: 10.1038/ng1576 .

JL Bennetzen, H Wang. The contributions of transposable elements to the structure, function, and evolution of plant genomes. Annu Rev Plant Biol. 2014;65:50530. Epub 2014/03/04. doi: doi: 10.1146/annurev-arplant-050213-035811 .

D Lisch. Epigenetic regulation of transposable elements in plants. Annu Rev Plant Biol. 2009;60:4366. Epub 2008/11/15. doi: doi: 10.1146/annurev.arplant.59.032607.092744 .

MA Grandbastien. LTR retrotransposons, handy hitchhikers of plant regulation and stress response. Biochim Biophys Acta. 2015;1849(4):40316. Epub 2014/08/03. doi: doi: 10.1016/j.bbagrm.2014.07.017 .

MA Grandbastien, C Audeon, E Bonnivard, JM Casacuberta, B Chalhoub, AP Costa, et al . Stress activation and genomic impact of Tnt1 retrotransposons in Solanaceae. Cytogenet Genome Res. 2005;110(1–4):22941. Epub 2005/08/12. doi: doi: 10.1159/000084957 .

H Ito, H Gaubert, E Bucher, M Mirouze, I Vaillant, J Paszkowski. An siRNA pathway prevents transgenerational retrotransposition in plants subjected to stress. Nature. 2011;472(7341):1159. Epub 2011/03/15. doi: doi: 10.1038/nature09861 .


MY Kim, D Zilberman. DNA methylation as a system of plant genomic immunity. Trends Plant Sci. 2014;19(5):3206. Epub 2014/03/13. doi: doi: 10.1016/j.tplants.2014.01.014 .


H Zhang, Z Lang, JK Zhu. Dynamics and function of DNA methylation in plants. Nat Rev Mol Cell Biol. 2018;19(8):489506. Epub 2018/05/23. doi: doi: 10.1038/s41580-018-0016-z .


JA Law, SE Jacobsen. Establishing, maintaining and modifying DNA methylation patterns in plants and animals. Nat Rev Genet. 2010;11(3):20420. Epub 2010/02/10. doi: doi: 10.1038/nrg2719 ; PubMed Central PMCID: PMC3034103.


G Martinez, RK Slotkin. Developmental relaxation of transposable element silencing in plants: functional or byproduct? Curr Opin Plant Biol. 2012;15(5):496502. Epub 2012/10/02. doi: doi: 10.1016/j.pbi.2012.09.001 .


JM Noshay, SN Anderson, P Zhou, L Ji, W Ricci, Z Lu, et al . Monitoring the interplay between transposable element families and DNA methylation in maize. PLoS Genet. 2019;15(9):e1008291. Epub 2019/09/10. doi: doi: 10.1371/journal.pgen.1008291 ; PubMed Central PMCID: PMC6752859.


NM Springer, D Lisch, Q Li. Creating order from chaos: epigenome dynamics in plants with complex genomes. Plant Cell. 2016;28(2):31425. Epub 2016/02/13. doi: doi: 10.1105/tpc.15.00911 ; PubMed Central PMCID: PMC4790878.


MS Shook, EJ Richards. VIM proteins regulate transcription exclusively through the MET1 cytosine methylation pathway. Epigenetics. 2014;9(7):9806. Epub 2014/04/26. doi: doi: 10.4161/epi.28906 ; PubMed Central PMCID: PMC4143413.


HR Woo, TA Dittmer, EJ Richards. Three SRA-domain methylcytosine-binding proteins cooperate to maintain global CpG methylation and epigenetic silencing in Arabidopsis. PLoS Genet. 2008;4(8):e1000156. Epub 2008/08/16. doi: doi: 10.1371/journal.pgen.1000156 ; PubMed Central PMCID: PMC2491724.


SW Chan, IR Henderson, SE Jacobsen. Gardening the genome: DNA methylation in Arabidopsis thaliana. Nat Rev Genet. 2005;6(5):35160. Epub 2005/04/30. doi: doi: 10.1038/nrg1601 .


MJ Ronemus, M Galbiati, C Ticknor, J Chen, SL Dellaporta. Demethylation-induced developmental pleiotropy in Arabidopsis. Science. 1996;273(5275):6547. Epub 1996/08/02. doi: doi: 10.1126/science.273.5275.654 .


MW Kankel, DE Ramsey, TL Stokes, SK Flowers, JR Haag, JA Jeddeloh, et al . Arabidopsis MET1 cytosine methyltransferase mutants. Genetics. 2003;163(3):110922. Epub 2003/03/29. ; PubMed Central PMCID: PMC1462485.


J Du, LM Johnson, SE Jacobsen, DJ Patel. DNA methylation pathways and their crosstalk with histone methylation. Nat Rev Mol Cell Biol. 2015;16(9):51932. Epub 2015/08/22. doi: doi: 10.1038/nrm4043 ; PubMed Central PMCID: PMC4672940.


LM Johnson, M Bostick, X Zhang, E Kraft, I Henderson, J Callis, et al . The SRA methyl-cytosine-binding domain links DNA and histone methylation. Curr Biol. 2007;17(4):37984. Epub 2007/01/24. doi: doi: 10.1016/j.cub.2007.01.009 ; PubMed Central PMCID: PMC1850948.


WJ Soppe, Z Jasencakova, A Houben, T Kakutani, A Meister, MS Huang, et al . DNA methylation controls histone H3 lysine 9 methylation and heterochromatin assembly in Arabidopsis. EMBO J. 2002;21(23):654959. Epub 2002/11/29. doi: doi: 10.1093/emboj/cdf657 ; PubMed Central PMCID: PMC136960.


JM Wendte, CS Pikaard. The RNAs of RNA-directed DNA methylation. Biochim Biophys Acta Gene Regul Mech. 2017;1860(1):1408. Epub 2016/08/16. doi: doi: 10.1016/j.bbagrm.2016.08.004 ; PubMed Central PMCID: PMC5203809.


M Matzke, T Kanno, L Daxinger, B Huettel, AJ Matzke. RNA-mediated chromatin-based silencing in plants. Curr Opin Cell Biol. 2009;21(3):36776. Epub 2009/02/27. doi: doi: 10.1016/ .


MA Matzke, RA Mosher. RNA-directed DNA methylation: an epigenetic pathway of increasing complexity. Nat Rev Genet. 2014;15(6):394408. Epub 2014/05/09. doi: doi: 10.1038/nrg3683 .


LM Johnson, J Du, CJ Hale, S Bischof, S Feng, RK Chodavarapu, et al . SRA- and SET-domain-containing proteins link RNA polymerase V occupancy to DNA methylation. Nature. 2014;507(7490):1248. Epub 2014/01/28. doi: doi: 10.1038/nature12931 ; PubMed Central PMCID: PMC3963826.


W Liu, SH Duttke, J Hetzel, M Groth, S Feng, J Gallego-Bartolome, et al . RNA-directed DNA methylation involves co-transcriptional small-RNA-guided slicing of polymerase V transcripts in Arabidopsis. Nat Plants. 2018;4(3):1818. Epub 2018/01/31. doi: doi: 10.1038/s41477-017-0100-y ; PubMed Central PMCID: PMC5832601.


X Cao, W Aufsatz, D Zilberman, MF Mette, MS Huang, M Matzke, et al . Role of the DRM and CMT3 methyltransferases in RNA-directed DNA methylation. Curr Biol. 2003;13(24):22127. Epub 2003/12/19. doi: doi: 10.1016/j.cub.2003.11.052 .


H Stroud, T Do, J Du, X Zhong, S Feng, L Johnson, et al . Non-CG methylation patterns shape the epigenetic landscape in Arabidopsis. Nat Struct Mol Biol. 2014;21(1):6472. Epub 2013/12/18. doi: doi: 10.1038/nsmb.2735 ; PubMed Central PMCID: PMC4103798.


A Zemach, MY Kim, PH Hsieh, D Coleman-Derr, L Eshed-Williams, K Thao, et al . The Arabidopsis nucleosome remodeler DDM1 allows DNA methyltransferases to access H1-containing heterochromatin. Cell. 2013;153(1):193205. Epub 2013/04/02. doi: doi: 10.1016/j.cell.2013.02.033 ; PubMed Central PMCID: PMC4035305.


JM Wendte, RJ Schmitz. Specifications of targeting heterochromatin modifications in plants. Mol Plant. 2018;11(3):3817. Epub 2017/10/17. doi: doi: 10.1016/j.molp.2017.10.002 .


Q Li, JI Gent, G Zynda, J Song, I Makarevitch, CD Hirsch, et al . RNA-directed DNA methylation enforces boundaries between heterochromatin and euchromatin in the maize genome. Proc Natl Acad Sci U S A. 2015;112(47):1472833. Epub 2015/11/11. doi: doi: 10.1073/pnas.1514680112 ; PubMed Central PMCID: PMC4664327.


JB Hollick. Paramutation and related phenomena in diverse species. Nat Rev Genet. 2017;18(1):523. Epub 2016/11/01. doi: doi: 10.1038/nrg.2016.115 .


M Alleman, L Sidorenko, K McGinnis, V Seshadri, JE Dorweiler, J White, et al . An RNA-dependent RNA polymerase is required for paramutation in maize. Nature. 2006;442(7100):2958. Epub 2006/07/21. doi: doi: 10.1038/nature04884 .


MR Woodhouse, M Freeling, D Lisch. The mop1 (mediator of paramutation1) mutant progressively reactivates one of the two genes encoded by the MuDR transposon in maize. Genetics. 2006;172(1):57992. Epub 2005/10/13. doi: doi: 10.1534/genetics.105.051383 ; PubMed Central PMCID: PMC1456185.


K Nobuta, C Lu, R Shrivastava, M Pillay, E De Paoli, M Accerbi, et al . Distinct size distribution of endogeneous siRNAs in maize: Evidence from deep sequencing in the mop1-1 mutant. Proc Natl Acad Sci U S A. 2008;105(39):1495863. Epub 2008/09/26. doi: doi: 10.1073/pnas.0808066105 ; PubMed Central PMCID: PMC2567475.


Y Jia, DR Lisch, K Ohtsu, MJ Scanlon, D Nettleton, PS Schnable. Loss of RNA-dependent RNA polymerase 2 (RDR2) function causes widespread and unexpected changes in the expression of transposons, genes, and 24-nt small RNAs. PLoS Genet. 2009;5(11):e1000737. Epub 2009/11/26. doi: doi: 10.1371/journal.pgen.1000737 ; PubMed Central PMCID: PMC2774947.


JE Dorweiler, CC Carey, KM Kubo, JB Hollick, JL Kermicle, VL Chandler. Mediator of paramutation1 is required for establishment and maintenance of paramutation at multiple maize loci. Plant Cell. 2000;12(11):210118. Epub 2000/11/23. doi: doi: 10.1105/tpc.12.11.2101 ; PubMed Central PMCID: PMC150161.


JF Hultquist, JE Dorweiler. Feminized tassels of maize mop1 and ts1 mutants exhibit altered levels of miR156 and specific SBP-box genes. Planta. 2008;229(1):99113. Epub 2008/09/19. doi: doi: 10.1007/s00425-008-0813-2 .


S Vendramin, J Huang, PA Crisp, TF Madzima, KM McGinnis. Epigenetic regulation of ABA-induced transcriptional responses in maize. G3 (Bethesda). 2020;10(5):172743. Epub 2020/03/18. doi: doi: 10.1534/g3.119.400993 ; PubMed Central PMCID: PMC7202028.


KF Erhard Jr, JL Stonaker, SE Parkinson, JP Lim, CJ Hale, JB Hollick. RNA polymerase IV functions in paramutation in Zea mays. Science. 2009;323(5918):12015. Epub 2009/03/03. doi: doi: 10.1126/science.1164508 .


C Forestan, R Aiese Cigliano, S Farinati, A Lunardon, W Sanseverino, S Varotto. Stress-induced and epigenetic-mediated maize transcriptome regulation study by means of transcriptome reannotation and differential expression analysis. Sci Rep. 2016;6:30446. Epub 2016/07/28. doi: doi: 10.1038/srep30446 ; PubMed Central PMCID: PMC4962059.


A Lunardon, C Forestan, S Farinati, MJ Axtell, S Varotto. Genome-wide characterization of maize small RNA loci and their regulation in the required to maintain repression6-1 (rmr6-1) mutant and long-term abiotic stresses. Plant Physiol. 2016;170(3):153548. Epub 2016/01/10. doi: doi: 10.1104/pp.15.01205 ; PubMed Central PMCID: PMC4775107.


BC Mikula. Environmental programming of heritable epigenetic changes in paramutant r-gene expression using temperature and light at a specific stage of early development in maize seedlings. Genetics. 1995;140(4):137987. Epub 1995/08/01. ; PubMed Central PMCID: PMC1206701.


E Heard, RA Martienssen. Transgenerational epigenetic inheritance: myths and mechanisms. Cell. 2014;157(1):95109. Epub 2014/04/01. doi: doi: 10.1016/j.cell.2014.02.045 ; PubMed Central PMCID: PMC4020004.


S Tsukahara, A Kobayashi, A Kawabe, O Mathieu, A Miura, T Kakutani. Bursts of retrotransposition reproduced in Arabidopsis. Nature. 2009;461(7262):4236. Epub 2009/09/08. doi: doi: 10.1038/nature08351 .


Z Lippman, B May, Yordan, Singer T, Martienssen R. Distinct Mechanisms Determine Transposon Inheritance and Methylation via Small Interfering RNA and Histone Modification. PLoS Biol. 2003;1(3):420. doi: doi: 10.1371/journal.pbio.0000067 ; PMCID: PMC300680.


FK Teixeira, F Heredia, A Sarazin, F Roudier, M Boccara, C Ciaudo, et al . A role for RNAi in the selective correction of DNA methylation defects. Science. 2009;323(5921):16004. Epub 2009/01/31. doi: doi: 10.1126/science.1165313 .


TK To, Y Nishizawa, S Inagaki, Y Tarutani, S Tominaga, A Toyoda, et al . RNA interference-independent reprogramming of DNA methylation in Arabidopsis. Nat Plants. 2020. Epub 2020/12/02. doi: doi: 10.1038/s41477-020-00810-z .


B Schuettengruber, HM Bourbon, L Di Croce, G Cavalli. Genome Regulation by Polycomb and Trithorax: 70 Years and Counting. Cell. 2017;171(1):3457. Epub 2017/09/25. doi: doi: 10.1016/j.cell.2017.08.002 .


I Mozgova, L Hennig. The polycomb group protein regulatory network. Annu Rev Plant Biol. 2015;66:26996. Epub 2015/01/27. doi: doi: 10.1146/annurev-arplant-043014-115627 .


PA Steffen, L Ringrose. What are memories made of? How Polycomb and Trithorax proteins mediate epigenetic memory. Nat Rev Mol Cell Biol. 2014;15(5):34056. Epub 2014/04/24. doi: doi: 10.1038/nrm3789 .


J Xiao, D Wagner. Polycomb repression in the regulation of growth and development in Arabidopsis. Curr Opin Plant Biol. 2015;23:1524. Epub 2014/12/03. doi: doi: 10.1016/j.pbi.2014.10.003 .


F Roudier, I Ahmed, C Berard, A Sarazin, T Mary-Huard, S Cortijo, et al . Integrative epigenomic mapping defines four main chromatin states in Arabidopsis. EMBO J. 2011;30(10):192838. Epub 2011/04/14. doi: doi: 10.1038/emboj.2011.103 ; PubMed Central PMCID: PMC3098477.


Y You, A Sawikowska, M Neumann, D Pose, G Capovilla, T Langenecker, et al . Temporal dynamics of gene expression and histone marks at the Arabidopsis shoot meristem during flowering. Nat Commun. 2017;8:15120. Epub 2017/05/18. doi: doi: 10.1038/ncomms15120 ; PubMed Central PMCID: PMC5442315.


ES Gan, Y Xu, T Ito. Dynamics of H3K27me3 methylation and demethylation in plant development. Plant Signal Behav. 2015;10(9):e1027851. Epub 2015/08/28. doi: doi: 10.1080/15592324.2015.1027851 ; PubMed Central PMCID: PMC4883920.


S Costa, C Dean. Storing memories: the distinct phases of Polycomb-mediated silencing of Arabidopsis FLC. Biochem Soc Trans. 2019;47(4):118796. Epub 2019/07/07. doi: doi: 10.1042/BST20190255 .


I Searle, Y He, F Turck, C Vincent, F Fornara, S Krober, et al . The transcription factor FLC confers a flowering response to vernalization by repressing meristem competence and systemic signaling in Arabidopsis. Genes Dev. 2006;20(7):898912. Epub 2006/04/08. doi: doi: 10.1101/gad.373506 ; PubMed Central PMCID: PMC1472290.


P Crevillen, H Yang, X Cui, C Greeff, M Trick, Q Qiu, et al . Epigenetic reprogramming that prevents transgenerational inheritance of the vernalized state. Nature. 2014;515(7528):58790. Epub 2014/09/16. doi: doi: 10.1038/nature13722 ; PubMed Central PMCID: PMC4247276.


Z Tao, L Shen, X Gu, Y Wang, H Yu, Y He. Embryonic epigenetic reprogramming by a pioneer transcription factor in plants. Nature. 2017;551(7678):1248. Epub 2017/10/27. doi: doi: 10.1038/nature24300 .


M Borg, Y Jacob, D Susaki, C LeBlanc, D Buendia, E Axelsson, et al . Targeted reprogramming of H3K27me3 resets epigenetic memory in plant paternal chromatin. Nat Cell Biol. 2020;22(6):6219. Epub 2020/05/13. doi: doi: 10.1038/s41556-020-0515-y .


TP Michael. Plant genome size variation: bloating and purging DNA. Brief Funct Genomics. 2014;13(4):30817. Epub 2014/03/22. doi: doi: 10.1093/bfgp/elu005 .


K Naito, F Zhang, T Tsukiyama, H Saito, CN Hancock, AO Richardson, et al . Unexpected consequences of a sudden and massive transposon amplification on rice gene expression. Nature. 2009;461(7267):11304. Epub 2009/10/23. doi: doi: 10.1038/nature08479 .


SR Wessler. Turned on by stress. Plant retrotransposons. Curr Biol. 1996;6(8):95961. Epub 1996/08/01. doi: doi: 10.1016/s0960-9822(02)00638-3 .


P Woodrow, G Pontecorvo, LF Ciarmiello, A Fuggi, P Carillo. Ttd1a promoter is involved in DNA-protein binding by salt and light stresses. Mol Biol Rep. 2011;38(6):378794. Epub 2010/11/26. doi: doi: 10.1007/s11033-010-0494-3 .


Y Kimura, Y Tosa, S Shimada, R Sogo, M Kusaba, T Sunaga, et al . OARE-1, a Ty1-copia retrotransposon in oat activated by abiotic and biotic stresses. Plant Cell Physiol. 2001;42(12):134554. Epub 2002/01/05. doi: doi: 10.1093/pcp/pce171 .


I Makarevitch, AJ Waters, PT West, M Stitzer, CN Hirsch, J Ross-Ibarra, et al . Transposable elements contribute to activation of maize genes in response to abiotic stress. PLoS Genet. 2015;11(1):e1004915. Epub 2015/01/09. doi: doi: 10.1371/journal.pgen.1004915 ; PubMed Central PMCID: PMC4287451.


K Kashkush, M Feldman, AA Levy. Transcriptional activation of retrotransposons alters the expression of adjacent genes in wheat. Nat Genet. 2003;33(1):1026. Epub 2002/12/17. doi: doi: 10.1038/ng1063 .


C Lang-Mladek, O Popova, K Kiok, M Berlinger, B Rakic, W Aufsatz, et al . Transgenerational inheritance and resetting of stress-induced loss of epigenetic gene silencing in Arabidopsis. Mol Plant. 2010;3(3):594602. Epub 2010/04/23. doi: doi: 10.1093/mp/ssq014 ; PubMed Central PMCID: PMC2877484.


J Molinier, G Ries, C Zipfel, B Hohn. Transgeneration memory of stress in plants. Nature. 2006;442(7106):10469. Epub 2006/08/08. doi: doi: 10.1038/nature05022 .


VV Cavrak, N Lettner, S Jamge, A Kosarewicz, LM Bayer, O Mittelsten Scheid. How a retrotransposon exploits the plant’s heat stress response for its activation. PLoS Genet. 2014;10(1):e1004115. Epub 2014/02/06. doi: doi: 10.1371/journal.pgen.1004115 ; PubMed Central PMCID: PMC3907296.


K Ohtsu, MB Smith, SJ Emrich, LA Borsuk, R Zhou, T Chen, et al . Global gene expression analysis of the shoot apical meristem of maize (Zea mays L.). Plant J. 2007;52(3):391404. Epub 2007/09/04. doi: doi: 10.1111/j.1365-313X.2007.03244.x ; PubMed Central PMCID: PMC2156186.


T Baubec, A Finke, O Mittelsten Scheid, A Pecinka. Meristem-specific expression of epigenetic regulators safeguards transposon silencing in Arabidopsis. EMBO Rep. 2014;15(4):44652. Epub 2014/02/25. doi: doi: 10.1002/embr.201337915 ; PubMed Central PMCID: PMC3989676.


Y Hayashi, K Takehira, K Nozawa, T Suzuki, Y Masuta, A Kato, et al . ONSEN shows different transposition activities in RdDM pathway mutants. Genes Genet Syst. 2020. Epub 2020/09/08. doi: doi: 10.1266/ggs.20-00019 .


M Iwasaki, J Paszkowski. Identification of genes preventing transgenerational transmission of stress-induced epigenetic states. Proc Natl Acad Sci U S A. 2014;111(23):854752. Epub 2014/06/10. doi: doi: 10.1073/pnas.1402275111 ; PubMed Central PMCID: PMC4060648.


O Mathieu, J Reinders, M Caikovski, C Smathajitt, J Paszkowski. Transgenerational stability of the Arabidopsis epigenome is coordinated by CG methylation. Cell. 2007;130(5):85162. Epub 2007/09/07. doi: doi: 10.1016/j.cell.2007.07.007 .


Z Lippman, B May, C Yordan, T Singer, R Martienssen. Distinct mechanisms determine transposon inheritance and methylation via small interfering RNA and histone modification. PLoS Biol. 2003;1(3):E67. Epub 2003/12/24. doi: doi: 10.1371/journal.pbio.0000067 ; PubMed Central PMCID: PMC300680.


MJ Sigman, RK Slotkin. The first rule of plant transposable element silencing: location, location, location. Plant Cell. 2016;28(2):30413. Epub 2016/02/13. doi: doi: 10.1105/tpc.15.00869 ; PubMed Central PMCID: PMC4790875.


D Lisch. Mutator and MULE transposons. Microbiol Spectr. 2015;3(2):MDNA3-0032-2014. Epub 2015/06/25. doi: doi: 10.1128/microbiolspec.MDNA3-0032-2014 .


D Lisch. Mutator transposons. Trends Plant Sci. 2002;7(11):498504. Epub 2002/11/06. doi: doi: 10.1016/s1360-1385(02)02347-6 .


D Lisch, L Girard, M Donlin, M Freeling. Functional analysis of deletion derivatives of the maize transposon MuDR delineates roles for the MURA and MURB proteins. Genetics. 1999;151(1):33141. Epub 1999/01/05. ; PubMed Central PMCID: PMC1460458.


GN Rudenko, V Walbot. Expression and post-transcriptional regulation of maize transposable element MuDR and its derivatives. Plant Cell. 2001;13(3):55370. Epub 2001/03/17. doi: doi: 10.1105/tpc.13.3.553 ; PubMed Central PMCID: PMC135511.


RK Slotkin, M Freeling, D Lisch. Mu killer causes the heritable inactivation of the Mutator family of transposable elements in Zea mays. Genetics. 2003;165(2):78197. Epub 2003/10/24. ; PubMed Central PMCID: PMC1462800.


D Lisch, P Chomet, M Freeling. Genetic characterization of the Mutator system in maize: behavior and regulation of Mu transposons in a minimal line. Genetics. 1995;139(4):177796. Epub 1995/04/01. ; PubMed Central PMCID: PMC1206502.


H Li, M Freeling, D Lisch. Epigenetic reprogramming during vegetative phase change in maize. Proc Natl Acad Sci U S A. 2010;107(51):221849. Epub 2010/12/08. doi: doi: 10.1073/pnas.1016884108 ; PubMed Central PMCID: PMC3009802.


D Burgess, H Li, M Zhao, SY Kim, D Lisch. Silencing of Mutator elements in maize involves distinct populations of small RNAs and distinct patterns of DNA methylation. Genetics. 2020;215(2):37991. Epub 2020/04/02. doi: doi: 10.1534/genetics.120.303033 ; PubMed Central PMCID: PMC7268996.


D Lisch, CC Carey, JE Dorweiler, VL Chandler. A mutation that prevents paramutation in maize also reverses Mutator transposon methylation and silencing. Proc Natl Acad Sci U S A. 2002;99(9):61305. Epub 2002/04/18. doi: doi: 10.1073/pnas.052152199 ; PubMed Central PMCID: PMC122914.


MR Woodhouse, M Freeling, D Lisch. Initiation, establishment, and maintenance of heritable MuDR transposon silencing in maize are mediated by distinct factors. PLoS Biol. 2006;4(10):e339. Epub 2006/09/14. doi: doi: 10.1371/journal.pbio.0040339 ; PubMed Central PMCID: PMC1563492.


SK Kenchanmane Raju, EJ Ritter, CE Niederhuth. Establishment, maintenance, and biological roles of non-CG methylation in plants. Essays Biochem. 2019;63(6):74355. Epub 2019/10/28. doi: doi: 10.1042/EBC20190032 ; PubMed Central PMCID: PMC6923318.


RS Poethig. The Maize Shoot. In: M Freeling. WV, editor. The Maize Handbook New York, NY.: Springer ; 1994. p. 117.


J Hanway. How a corn plant develops. Ames, IA: Iowa State University Cooperative Extension Service; 1966.


J Du, LM Johnson, M Groth, S Feng, CJ Hale, S Li, et al . Mechanism of DNA methylation-directed histone methylation by KRYPTONITE. Mol Cell. 2014;55(3):495504. Epub 2014/07/16. doi: doi: 10.1016/j.molcel.2014.06.009 ; PubMed Central PMCID: PMC4127122.


P Amedeo, Y Habu, K Afsar, O Mittelsten Scheid, J Paszkowski. Disruption of the plant gene MOM releases transcriptional silencing of methylated genes. Nature. 2000;405(6783):2036. Epub 2000/05/23. doi: doi: 10.1038/35012108 .


M Caikovski, C Yokthongwattana, Y Habu, T Nishimura, O Mathieu, J Paszkowski. Divergent evolution of CHD3 proteins resulted in MOM1 refining epigenetic control in vascular plants. PLoS Genet. 2008;4(8):e1000165. Epub 2008/08/30. doi: doi: 10.1371/journal.pgen.1000165 ; PubMed Central PMCID: PMC2507757.


AV Probst, PF Fransz, J Paszkowski, O Mittelsten Scheid. Two means of transcriptional reactivation within heterochromatin. Plant J. 2003;33(4):7439. Epub 2003/03/01. doi: doi: 10.1046/j.1365-313x.2003.01667.x .


G Moissiard, SJ Cokus, J Cary, S Feng, AC Billi, H Stroud, et al . MORC family ATPases required for heterochromatin condensation and gene silencing. Science. 2012;336(6087):144851. Epub 2012/05/05. doi: doi: 10.1126/science.1221472 ; PubMed Central PMCID: PMC3376212.


Y Jacob, S Feng, CA LeBlanc, YV Bernatavichute, H Stroud, S Cokus, et al . ATXR5 and ATXR6 are H3K27 monomethyltransferases required for chromatin structure and gene silencing. Nat Struct Mol Biol. 2009;16(7):7638. Epub 2009/06/09. doi: doi: 10.1038/nsmb.1611 ; PubMed Central PMCID: PMC2754316.


CJ Harris, D Husmann, W Liu, FE Kasmi, H Wang, A Papikian, et al . Arabidopsis AtMORC4 and AtMORC7 Form Nuclear Bodies and Repress a Large Number of Protein-Coding Genes. PLoS Genet. 2016;12(5):e1005998. Epub 2016/05/14. doi: doi: 10.1371/journal.pgen.1005998 ; PubMed Central PMCID: PMC4865129.


Y Ikeda, T Pelissier, P Bourguet, C Becker, MN Pouch-Pelissier, R Pogorelcnik, et al . Arabidopsis proteins with a transposon-related domain act in gene silencing. Nat Commun. 2017;8:15122. Epub 2017/05/04. doi: doi: 10.1038/ncomms15122 ; PubMed Central PMCID: PMC5418596.


W Matsunaga, A Kobayashi, A Kato, H Ito. The effects of heat induction and the siRNA biogenesis pathway on the transgenerational transposition of ONSEN, a copia-like retrotransposon in Arabidopsis thaliana. Plant Cell Physiol. 2012;53(5):82433. Epub 2011/12/17. doi: doi: 10.1093/pcp/pcr179 .


P Chomet, D Lisch, KJ Hardeman, VL Chandler, M Freeling. Identification of a regulatory transposon that controls the Mutator transposable element system in maize. Genetics. 1991;129(1):26170. Epub 1991/09/01. ; PubMed Central PMCID: PMC1204575.


M Kato, K Takashima, T Kakutani. Epigenetic control of CACTA transposon mobility in Arabidopsis thaliana. Genetics. 2004;168(2):9619. Epub 2004/10/30. doi: doi: 10.1534/genetics.104.029637 ; PubMed Central PMCID: PMC1448851.


A Mari-Ordonez, A Marchais, M Etcheverry, A Martin, V Colot, O Voinnet. Reconstructing de novo silencing of an active plant retrotransposon. Nat Genet. 2013;45(9):102939. Epub 2013/07/16. doi: doi: 10.1038/ng.2703 .


H Hirochika, K Sugimoto, Y Otsuki, H Tsugawa, M Kanda. Retrotransposons of rice involved in mutations induced by tissue culture. Proc Natl Acad Sci U S A. 1996;93(15):77838. Epub 1996/07/23. doi: doi: 10.1073/pnas.93.15.7783 ; PubMed Central PMCID: PMC38825.


B Beydler, K Osadchuk, CL Cheng, JR Manak, EE Irish. The juvenile phase of maize sees upregulation of stress-response genes and is extended by exogenous jasmonic acid. Plant Physiol. 2016;171(4):264858. Epub 2016/06/17. doi: doi: 10.1104/pp.15.01707 ; PubMed Central PMCID: PMC4972259.


CC Sheldon, MJ Hills, C Lister, C Dean, ES Dennis, WJ Peacock. Resetting of FLOWERING LOCUS C expression after epigenetic repression by vernalization. Proc Natl Acad Sci U S A. 2008;105(6):22149. Epub 2008/02/06. doi: doi: 10.1073/pnas.0711453105 ; PubMed Central PMCID: PMC2542874.


NC Deans, BJ Giacopelli, JB Hollick. Locus-specific paramutation in Zea mays is maintained by a PICKLE-like chromodomain helicase DNA-binding 3 protein controlling development and male gametophyte function. PLoS Genet. 2020;16(12):e1009243. Epub 2020/12/16. doi: doi: 10.1371/journal.pgen.1009243 ; PubMed Central PMCID: PMC7837471 The Regents of the University of California.


C Li, L Gu, L Gao, C Chen, CQ Wei, Q Qiu, et al . Concerted genomic targeting of H3K27 demethylase REF6 and chromatin-remodeling ATPase BRM in Arabidopsis. Nat Genet. 2016;48(6):68793. Epub 2016/04/26. doi: doi: 10.1038/ng.3555 ; PubMed Central PMCID: PMC5134324.


J Liu, L Feng, X Gu, X Deng, Q Qiu, Q Li, et al . An H3K27me3 demethylase-HSFA2 regulatory loop orchestrates transgenerational thermomemory in Arabidopsis. Cell Res. 2019;29(5):37990. Epub 2019/02/20. doi: doi: 10.1038/s41422-019-0145-8 ; PubMed Central PMCID: PMC6796840.


X Cui, F Lu, Q Qiu, B Zhou, L Gu, S Zhang, et al . REF6 recognizes a specific DNA sequence to demethylate H3K27me3 and regulate organ boundary formation in Arabidopsis. Nat Genet. 2016;48(6):6949. Epub 2016/04/26. doi: doi: 10.1038/ng.3556 .


X Chen, X Liu, Y Zhao, DX Zhou. Histone H3K4me3 and H3K27me3 regulatory genes control stable transmission of an epimutation in rice. Sci Rep. 2015;5:13251. Epub 2015/08/20. doi: doi: 10.1038/srep13251 ; PubMed Central PMCID: PMC4541256.


X Cui, P Jin, X Cui, L Gu, Z Lu, Y Xue, et al . Control of transposon activity by a histone H3K4 demethylase in rice. Proc Natl Acad Sci U S A. 2013;110(5):19538. Epub 2013/01/16. doi: doi: 10.1073/pnas.1217020110 ; PubMed Central PMCID: PMC3562835.


E Gruntman, Y Qi, RK Slotkin, T Roeder, RA Martienssen, R Sachidanandam. Kismeth: analyzer of plant methylation states through bisulfite sequencing. BMC Bioinformatics. 2008;9:371. Epub 2008/09/13. doi: doi: 10.1186/1471-2105-9-371 ; PubMed Central PMCID: PMC2553349.


L He, W Wu, G Zinta, L Yang, D Wang, R Liu, et al . A naturally occurring epiallele associates with leaf senescence and local climate adaptation in Arabidopsis accessions. Nat Commun. 2018;9(1):460. Epub 2018/02/02. doi: doi: 10.1038/s41467-018-02839-3 ; PubMed Central PMCID: PMC5792623.


B Carter, B Bishop, KK Ho, R Huang, W Jia, H Zhang, et al . The Chromatin Remodelers PKL and PIE1 Act in an Epigenetic Pathway That Determines H3K27me3 Homeostasis in Arabidopsis. Plant Cell. 2018;30(6):133752. Epub 2018/05/29. doi: doi: 10.1105/tpc.17.00867 ; PubMed Central PMCID: PMC6048792.


M Haring, S Offermann, T Danker, I Horst, C Peterhansel, M Stam. Chromatin immunoprecipitation: optimization, quantitative analysis and data normalization. Plant Methods. 2007;3:11. Epub 2007/09/26. doi: doi: 10.1186/1746-4811-3-11 ; PubMed Central PMCID: PMC2077865.
Supplementary materials
  • pgen.1009326.s001.tif info     save_alt  
  • pgen.1009326.s002.tif info     save_alt  
  • pgen.1009326.s003.tif info     save_alt  
  • pgen.1009326.s004.tif info     save_alt  
  • pgen.1009326.s005.tif info     save_alt  
  • pgen.1009326.s006.xlsx info     save_alt  
  • pgen.1009326.s007.xlsx info     save_alt  
  • pgen.1009326.s008.xlsx info     save_alt  
  • pgen.1009326.s009.xlsx info     save_alt  
  • pgen.1009326.s010.xlsx info     save_alt  
  • pgen.1009326.s011.docx info     save_alt  
  • pgen.1009326.s012.docx info     save_alt is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. DNA methylation prevents rapid and heritable reversal of transposon silencing under heat stress in <i>Zea mays</i>&author=Wei Guo,Dafang Wang,Damon Lisch,Ortrun Mittelsten Scheid,&keyword=&subject=Research Article,Biology and life sciences,Cell biology,Chromosome biology,Chromatin,Chromatin modification,DNA methylation,Biology and life sciences,Genetics,Epigenetics,Chromatin,Chromatin modification,DNA methylation,Biology and life sciences,Genetics,Gene expression,Chromatin,Chromatin modification,DNA methylation,Biology and life sciences,Genetics,DNA,DNA modification,DNA methylation,Biology and life sciences,Biochemistry,Nucleic acids,DNA,DNA modification,DNA methylation,Biology and life sciences,Genetics,Epigenetics,DNA modification,DNA methylation,Biology and life sciences,Genetics,Gene expression,DNA modification,DNA methylation,Biology and Life Sciences,Genetics,Gene Expression,Gene Regulation,Silencer Elements,Physical Sciences,Physics,Classical Mechanics,Mechanical Stress,Thermal Stresses,Biology and Life Sciences,Plant Science,Plant Anatomy,Leaves,Biology and Life Sciences,Genetics,Genetic Elements,Mobile Genetic Elements,Transposable Elements,Biology and Life Sciences,Genetics,Genomics,Mobile Genetic Elements,Transposable Elements,Biology and life sciences,Biochemistry,Proteins,DNA-binding proteins,Histones,Research and Analysis Methods,Animal Studies,Experimental Organism Systems,Model Organisms,Maize,Research and Analysis Methods,Model Organisms,Maize,Biology and Life Sciences,Organisms,Eukaryota,Plants,Grasses,Maize,Research and Analysis Methods,Animal Studies,Experimental Organism Systems,Plant and Algal Models,Maize,Biology and Life Sciences,Cell Biology,Chromosome Biology,Chromatin,Chromatin Modification,Histone Modification,Biology and Life Sciences,Genetics,Epigenetics,Chromatin,Chromatin Modification,Histone Modification,Biology and Life Sciences,Genetics,Gene Expression,Chromatin,Chromatin Modification,Histone Modification,Biology and Life Sciences,Genetics,Gene Expression,Histone Modification,