Skip to main content

Embryo CHH hypermethylation is mediated by RdDM and is autonomously directed in Brassica rapa

Abstract

Background

RNA-directed DNA methylation (RdDM) initiates cytosine methylation in all contexts and maintains asymmetric CHH methylation. Mature plant embryos show one of the highest levels of CHH methylation, and it has been suggested that RdDM is responsible for this hypermethylation. Because loss of RdDM in Brassica rapa causes seed abortion, embryo methylation might play a role in seed development. RdDM is required in the maternal sporophyte, suggesting that small RNAs from the maternal sporophyte might translocate to the developing embryo, triggering DNA methylation that prevents seed abortion. This raises the question of whether embryo hypermethylation is autonomously regulated by the embryo itself or influenced by the maternal sporophyte.

Results

Here, we demonstrate that B. rapa embryos are hypermethylated in both euchromatin and heterochromatin and that this process requires RdDM. Contrary to the current models, B. rapa embryo hypermethylation is not correlated with demethylation of the endosperm. We also show that maternal somatic RdDM is not sufficient for global embryo hypermethylation, and we find no compelling evidence for maternal somatic influence over embryo methylation at any locus. Decoupling of maternal and zygotic RdDM leads to successful seed development despite the loss of embryo CHH hypermethylation.

Conclusions

We conclude that embryo CHH hypermethylation is conserved, autonomously controlled, and not required for embryo development. Furthermore, maternal somatic RdDM, while required for seed development, does not directly influence embryo methylation patterns.

Background

DNA methylation is an epigenetic modification that can modulate chromatin structure and gene expression [1]. Plants methylate cytosines in all sequence contexts (CG, CHG, and CHH, where H is any base other than G) and use specific methyltransferases to maintain each context after replication [2]. In addition, the RNA-directed DNA methylation (RdDM) pathway is responsible for de novo methylation, a process that is most clearly observed at CHH positions [3]. RdDM functions primarily at the edges of euchromatin transposons, where constant re-establishment of methylation might be necessary [4, 5].

RdDM can be divided into siRNA production and DNA methylation stages. During siRNA production, RNA polymerase Pol IV produces single-stranded RNA transcripts which are copied into double-stranded RNA by RNA-DEPENDENT RNA POLYMERASE 2 (RDR2) and cut into 24-nucleotide small interfering (si)RNAs by DICER LIKE 3 (DCL3) [6,7,8]. To mediate DNA methylation, these 24-nt siRNAs are loaded onto ARGONAUTE 4 (AGO4), which interacts with a non-coding scaffold transcript produced by RNA polymerase V and recruits DOMAINS REARRANGED METHYLTRANSFERASE 2 (DRM2) to institute methylation marks on cytosine bases [9,10,11]. These two stages of RdDM frequently occur in cis but can also function in trans due to siRNA-AGO4 loading in the cytoplasm [12]. siRNAs can act in trans to trigger DNA methylation at allelic sites [13, 14] or at homologous non-allelic sites [15] or might move between cells to act non-cell autonomously [16].

With the exception of Arabidopsis, which has only a small reduction in seed size, loss of RdDM in most species results in disruption of reproductive development, indicating that RdDM is necessary for successful sexual reproduction [17,18,19,20,21]. Mature embryos accumulate high levels of CHH methylation in Arabidopsis, soybean, and chickpea, suggesting that RdDM might enable reproduction through hypermethylation of the mature embryo [22,23,24,25,26,27]. In Arabidopsis, the developing endosperm is demethylated at sequences that show hypermethylation in the embryo, leading to the hypothesis that siRNAs produced in the endosperm might move to the embryo to direct methylation [22, 28,29,30]. The movement of siRNAs between the maternal integuments and the filial tissues has also been proposed [31]. However, embryos produced through somatic embryogenesis also display hypermethylation, despite a lack of association with either endosperm or maternal integuments [27] and torpedo-stage Arabidopsis embryos accumulate roughly equal maternal and paternal siRNAs [32].

Here, we show that Brassica rapa mature embryos are hypermethylated in the CHH context in both euchromatin and heterochromatin, and we demonstrate that this process requires RdDM. Although maternal RdDM is required for seed development, it is not sufficient for embryo hypermethylation, and methylation in the CHH context is not necessary for proper seed development. Furthermore, we find no evidence that hypermethylation of the embryo is driven by siRNAs produced in adjacent tissues, suggesting that embryo CHH hypermethylation is entirely autonomous.

Results

Brassica rapa embryos are hypermethylated in the CHH context

To analyze global methylation levels in mature embryos, we performed whole-genome bisulfite sequencing on embryos dissected from dry seeds and compared the resulting data with other reproductive tissues (ovule, endosperm, and early embryo) and a non-reproductive control (leaf). Bisulfite conversion was greater than 99% in all samples, with a read depth coverage of > 9 (Additional file 1: Fig. S1). We calculated methylation in 300-bp non-overlapping windows for all three sequence contexts (Fig. 1a). CG methylation was largely unchanged across the tissues, with the exception of endosperm, which was demethylated in CG and CHG contexts, consistent with the observations in Arabidopsis and rice [33,34,35]. As in Arabidopsis, soybean, and chickpea [22,23,24,25,26,27], we also observed elevated global CHH methylation in B. rapa mature embryos, with moderately increased CHH methylation in torpedo-stage embryos (Fig. 1a and Additional file 1: Fig. S1). The increased CHH methylation in torpedo-stage embryos was correlated with CHH hypermethylation in mature embryos (Fig. 1b, correlation coefficient = 0.6), indicating that hypermethylation is a gradual process during embryogenesis.

Fig. 1
figure 1

B. rapa embryos are hypermethylated at CHH sites. a Distribution of methylation levels in 300-nt windows across the B. rapa genome. Ridge plots display the density of average methylation in each context, while background box plots enclose the 10th to 90th percentiles of the data. The black bar marks the median for each tissue/context combination. Only windows with a read depth ≥ 5 over all cytosines were included (approximately 1 million windows per tissue/context combination). b Increased CHH methylation (log2 fold change relative to leaves) is correlated in torpedo and mature embryos; 752,405 300-nt windows with a read depth of at least 5 in both tissues are plotted

To assess the types of chromatin responsible for embryo hypermethylation, we analyzed methylation levels in 25-kb windows across each chromosome (Fig. 2a and Additional file 1: Fig. S2). Pericentromeric heterochromatin, which has a denser accumulation of transposons, was strongly methylated in CG and CHG contexts in all tissues. The sole exception was endosperm, which as expected showed a small reduction in CG methylation and stronger loss of CHG methylation. In comparison with these heterochromatic marks, CHH methylation was distributed more equally across the length of the chromosome. Increased CHH methylation in mature embryos relative to other tissues was readily apparent. In Arabidopsis, embryo hypermethylation occurred primarily in pericentromeric heterochromatin [22, 23], while in soybean somatic embryos, CHH hypermethylation was seen across the genome [27]. We observed a similar degree of CHH hypermethylation in both the CG-dense pericentromeric regions and the chromosome arms (Fig. 2b), indicating that both heterochromatin and euchromatin are targets of CHH methylation during B. rapa embryo development.

Fig. 2
figure 2

Embryo CHH hypermethylation is not restricted to the pericentromere. a Heatmaps of transposon density or methylation level in 25-kb windows across chromosome 10. Each methylation context has its own scale bar to visualize changes across tissues. Other chromosomes are presented in Additional file 1: Fig. S2. b CHH hypermethylation in mature embryos (green line) is not correlated with the amount of CG methylation in leaves (purple line). Five-window rolling average of 25-kb windows across chromosome 10 are plotted

Finally, we assessed embryo methylation relative to leaf for each cytosine context within each 300-nt window. Most windows were unchanged with respect to CG and CHG methylation, while CHH methylation showed a pronounced shift toward hypermethylation (Fig. 3a). These changes were highly significant (Fig. 3b), but we selected only windows with the strongest changes in methylation for further analysis. We defined differentially methylated windows (DMWs) as those with at least 5-fold (log2 = 2.32) increase or decrease in methylation in mature embryo compared to leaf, and an FDR adjusted p value less than 0.005 (Fig. 3a, b, Additional file 1: Fig. S3). Hypermethylated windows were more abundant than hypomethylated windows for each sequence context, with CHH hyper-DMWs vastly outnumbering other DMWs in other sequence contexts (Fig. 3c).

Fig. 3
figure 3

Identification of embryo differentially methylated windows. Density distributions (a) and volcano plots (b) of methylation fold change in mature embryos versus leaf for three cytosine contexts; 300-nt windows with a read depth of at least 5 are plotted (650,665 CG, 686,741 CHG, and 869,526 CHH windows). The dashed blue line marks 5-fold hypomethylation, and the dashed green line marks 5-fold hypermethylation. Windows above this threshold with an FDR-adjusted p value <0.005 were collected for subsequent DMW analysis. c Number of differentially methylated windows (DMWs) passing the above thresholds in each methylation context

Together, our observations demonstrate that mature B. rapa embryos are extensively hypermethylated at CHH sites across the genome, and this hypermethylation is the primary difference between the leaf and embryo methylation patterns. This hypermethylation is widespread, not limited to pericentromeric heterochromatin, and progressive throughout embryogenesis.

Embryo CHH hypermethylation is dependent on RdDM

RdDM is the major pathway for de novo methylation in all sequence contexts, and its activity is frequently observed through the accumulation of CHH methylation. However, most of the CHH methylation in the genome is instead placed by CHROMOMETHYLTRANSFERASE 2 (CMT2) [18, 36]. Kawakatsu and colleagues [23] demonstrated that both RdDM and CMT2 contribute to CHH methylation in the embryo, but their analysis did not determine which process was responsible for the hypermethylation relative to non-embryonic tissues. Small RNA accumulation at hypermethylated regions is correlated with embryo hypermethylation [22, 24, 26], but it is not clear whether siRNA accumulation is required for increased methylation, or whether these two processes occur independently.

Our analysis also implicates RdDM in embryo CHH hypermethylation. Firstly, compared to all genomic windows with sufficient WGBS read depth, embryo CHH hyper-DMWs are significantly enriched for class II DNA elements and class III Helitrons (Fig. 4a). CHH hyper-DMWs are also significantly depleted at genes and at long terminal repeat (LTR) retroelements, following the characteristic pattern of RdDM loci in B. rapa [19]. Most importantly, CHH hyper-DMWs are enriched for loci previously shown to produce 24-nt siRNAs (Fig. 4a).

Fig. 4
figure 4

small RNA production is correlated with embryo hypermethylation, but not hypomethylation in leaves. a Enrichment or depletion of genomic features at CHH hyper-DMWs. The percentage of CHH hyper-DMWs overlapping annotated genomic features is plotted compared to the percentage of overlap for all genomic windows with comparable read depth. All differences are significant at p<2.2e−16. b, d CHH hyper-DMWs and all genomic windows were binned based on the number of sRNAs mapping to them in torpedo embryos or leaves, and the fraction of windows in each bin is shown. The embryo sRNA library has 39.1 million mapped reads, while the leaf library has 10.8 million reads. c Absolute CHH methylation in mature embryos is plotted as a function of the number of mapped sRNAs in torpedo embryos at all CHH hyper-DMWs. Box plots circumscribe the 10th–90th percentiles, and the black bar marks the median

To further investigate the association between CHH hyper-DMWs and siRNAs, we analyzed siRNA accumulation at CHH hyper-DMWs in torpedo-stage embryos (Additional file 1: Fig. S2). CHH hyper-DMWs have a significantly higher accumulation of siRNAs than the genome as a whole (Fig. 4b), and windows with greater siRNA accumulation in torpedo embryos have higher CHH methylation in mature embryos (Fig. 4c). However, many CHH hyper-DMWs lack substantial siRNA accumulation, a pattern also detected in chickpea [26]. We also observed a similar enrichment of siRNAs at CHH hyper-DMWs in leaves despite the 5-fold or greater difference in methylation between these tissues (Fig. 4d). This suggests that while production of 24-nt siRNAs is associated with embryo hypermethylation, similar 24-nt siRNA accumulation in leaves does not result in similar hypermethylation.

To directly test whether RdDM is required for embryo hypermethylation in B. rapa, we assayed the differences in methylation levels at CHH hyper-DMWs between wild type and an RdDM-deficient mutant, braA.rdr2-2 (rdr2 hereafter) (Grover et al. [19]). Mature rdr2 embryos have a clear reduction in CHH and CHG methylation compared to wild-type embryos (Fig. 5a, b), with CHH methylation levels similar to wild-type leaves. In contrast, CG methylation at these CHH hyper-DMWs is rdr2-independent. We also compared WT and rdr2 methylation at individual loci (Fig. 5c). Over 80% of CHH hyper-DMWs lose at least half of their methylation in the rdr2 mutant, demonstrating that embryo hypermethylation in the CHH context is primarily due to RdDM, and that other methylation pathways have only a minor contribution to embryo CHH hypermethylation.

Fig. 5
figure 5

Embryo hypermethylation is determined by filial genotype. a Distribution of CHH methylation at CHH hyper-DMWs in leaves and mature embryos. rdr2 embryos were derived either from rdr2 homozygous mothers (maroon) or from rdr2/+ heterozygous mothers (purple). b Cartoon and images of representative seeds measured in (a). Colored tissues have functional RdDM; gray tissues are deficient in RDR2. Scale bar is 5 mm. c Hex plots of mature embryo CHH methylation change by torpedo-stage endosperm CG (left) or CHG (right) methylation. d siRNA accumulation in endosperm at CHH hyper-DMWs. The endosperm siRNA library has 19.6 million mapped reads

No evidence for endosperm-directed hypermethylation of the embryo

It has been suggested that demethylation of the endosperm allows the production of siRNAs that target DNA methylation in the embryo [22, 28,29,30]. To determine whether there is an association between endosperm demethylation and embryo hypermethylation, we compared the changes in methylation levels in these tissues. Compared to leaf samples, endosperm is demethylated for both CG and CHG, while embryos are hypermethylated for CHH (Fig. 1a). However, there is no correlation between CG or CHG demethylation in the endosperm and CHH hypermethylation in the embryo, whether we assessed all genomic windows (Fig. 5d) or only the CHH hyper-DMWs (Additional file 1: Fig. S4). We measured correlations between embryo and endosperm methylation in multiple ways, both globally and at CHH hyper-DMWs (Additional file 1: Fig. S3). Absolute CHH methylation in the embryo is positively correlated with all methylation contexts in the endosperm, indicating that the embryo hypermethylated loci tend to have higher absolute methylation in the endosperm. When fold change in methylation relative to leaf samples are measured, only CHH methylation is correlated between these tissues, indicating that regardless of their absolute methylation level, the filial tissues are coordinately increasing CHH methylation. Only when we use the difference in absolute methylation as a metric do we find a correlation between endosperm CG/CHG demethylation and embryo CHH hypermethylation. However, this correlation is not strong and becomes weaker when only the CHH hyper-DMWs are assessed. Together, these results suggest that while similar loci may be demethylated in endosperm and hypermethylated in embryos, there is no evidence that demethylation of endosperm causes hypermethylation of embryos in B. rapa.

Because the presumed mechanism whereby endosperm demethylation triggers embryo hypermethylation is transport of siRNAs, we also assessed siRNA production at CHH hyper-DMWs in developing endosperm (Fig. 5e). CHH hyper-DMWs produce more siRNAs than the genome as a whole, but at a level that is comparable to developing embryos or leaves (Fig. 4b, d), suggesting that siRNA production occurs at these windows consistently across tissues and is not a response to endosperm demethylation. On the whole, we find no evidence that demethylation of the endosperm triggers siRNA production to cause hypermethylation of the mature embryo.

Maternal sporophytic RdDM is not sufficient for embryo hypermethylation

RdDM mutants in Capsella rubella and Brassica rapa show a high rate of seed abortion that is dependent on the maternal somatic genotype rather than the filial genotype [19, 37]. The few seeds that are produced from rdr2 plants are smaller and irregular in size and shape (Fig. 5b). Because a functional RdDM pathway in the maternal sporophyte is required for seed development, we hypothesized that maternal sporophytic RdDM might drive hypermethylation of the developing embryo. To test this hypothesis, we pollinated heterozygous (rdr2/RDR2) pistils with homozygous (rdr2/rdr2) pollen and identified rdr2 embryos that developed in the presence of functional maternal sporophytic RdDM. Compared to methylation levels of rdr2 mutant embryos from homozygous mutant mothers (rdr2/rdr2), we did not observe the restoration of embryo CHH hypermethylation (Fig. 5a, b).

To further probe the possibility that maternal sporophytically derived siRNAs might trigger hypermethylation of the embryo, we assessed DNA methylation levels at previously defined “siren” loci [31]. These loci produced over 90% of the 24-nt siRNAs in maternal integuments and are also the most highly expressed siRNA loci in the endosperm. We found no evidence that siRNAs produced from siren loci in the integument were able to direct DNA methylation in rdr2 embryos (Additional file 1: Fig. S5). These results indicate that although siRNA production in the maternal sporophyte is necessary for seed development, it is not sufficient for embryo hypermethylation.

Together, these observations provide no evidence that embryo hypermethylation is directed by siRNAs from another tissue. Combined with the observation that embryos formed through somatic embryogenesis also have elevated DNA methylation despite their lack of interaction with endosperm or integuments, we conclude that embryo hypermethylation is autonomously directed.

Discussion

Seeds form the majority of the world’s food supply, making the development of the seed and interactions between its multiple tissues critically important areas for research. Double fertilization gives rise to the diploid embryo and the triploid endosperm, which are surrounded by the seed coat, a maternal somatic tissue. Communication between maternal and filial tissues, as well as between the embryo and endosperm, is essential to coordinate the development of a seed [38, 39]. Small RNAs have been proposed to move between seed tissues and to establish robust methylation of transposons at this transition between generations [22, 28,29,30,31].

Here, we provide direct evidence that 24-nt siRNAs are responsible for the hypermethylation of mature embryos by demonstrating that rdr2 embryos lose hypermethylation. However, our evidence suggests that these siRNAs are derived autonomously in the embryo and are not transported from other tissues. Maternal sporophytic RDR2 (and hence, siRNA production) is not sufficient for embryo hypermethylation (Fig. 5a, b), clearly indicating that the siRNAs responsible for hypermethylation are produced in the filial tissues. Because the embryo and the endosperm have the same genotype (differing only in maternal ploidy), we cannot separate them genetically. However, we find that endosperm does not produce more siRNAs than embryos from CHH hyper-DMWs (Fig. 5e), nor is there a correlation between endosperm CG/CHG demethylation and embryo CHH hypermethylation (Fig. 5d). Furthermore, somatic soybean embryos produced in tissue culture also display embryo hypermethylation [27]. The most parsimonious explanation for these observations is that embryo CHH hypermethylation is autonomously directed by siRNAs synthesized in the embryo.

Although maternal sporophytic siRNA production is not sufficient for embryo hypermethylation, it remains possible that siRNAs from the maternal integument might trigger the expression of 24-nt siRNAs in the embryo and initiate autonomous methylation in the embryo. This process would be analogous to the production of piRNAs in Drosophila melanogaster, whereby maternally derived small RNAs initiate subsequent filial siRNA production and transposon silencing [40]. It also remains possible that triggering siRNAs are brought to the zygote during fertilization by the sperm nucleus, although this model remains to be tested [41, 42].

Despite a lack of embryo hypermethylation, rdr2 homozygous seeds from heterozygous mothers phenocopy wild-type seeds (Fig. 5b), indicating that embryo hypermethylation is not necessary for seed development in B. rapa. Similarly, Arabidopsis does not require DRM2 methyltransferase for embryo development [25]; however, Arabidopsis does not require RdDM generally, while other species in Brassicaceae have reproductive defects in the absence of RdDM [19]. Decoupling of embryo methylation and seed development in B. rapa supports the hypothesis that embryo hypermethylation is important for seed dormancy or longevity but not for seed development [22,23,24]. We assessed segregating seed populations and observed no difference in germination timing or frequency for unmethylated rdr2 embryos relative to their methylated siblings (data not shown), suggesting that other hypotheses should also be considered.

In Arabidopsis, embryo hypermethylation is preferentially targeted to transposons in the pericentromeric heterochromatin [22, 23], while hypermethylation also occurs at euchromatic transposons in soybean [25, 27]. In B. rapa, we detect hypermethylation in both heterochromatin and euchromatin (Fig. 2b), suggesting that euchromatic embryo hypermethylation might be common among plants. Recent work demonstrates that Arabidopsis heterochromatin is decondensed and produces abundant 24-nt siRNAs during embryogenesis [43], providing an opportunity for the RdDM machinery to access this chromatin for hypermethylation.

Our demonstration that RdDM is responsible for the hypermethylation of embryos leads to the hypothesis that siRNAs would be abundant during embryogenesis. However, we were surprised by the low level of siRNAs at CHH hyper-DMWs in torpedo-stage embryos. Correlation between CHH levels in torpedo-stage and mature embryos (Fig. 1b) indicates that hypermethylation occurs throughout embryogenesis rather than during embryo maturation, and therefore, robust siRNA accumulation would be predicted. The 81,556 CHH hyper-DMWs account for ~10% of windows with sufficient read depth, and they accumulate 16.9% of the mapped siRNAs (20.1% of the mapped 24-nt siRNAs). While this is a substantial enrichment compared to the genome as a whole, these windows account for 13.8% of the mapped siRNAs (15.9% of the mapped 24-nt siRNAs) in leaves. This discrepancy suggests that while siRNA production is required for embryo hypermethylation, developmental-specific factors are required for robust methylation.

Conclusions

Brassica rapa embryos are hypermethylated at both euchromatic and heterochromatic CHH positions. This hypermethylation requires RdDM, and there is no evidence that siRNAs from the endosperm or maternal somatic tissue direct embryo methylation. Successful development of seeds lacking embryo hypermethylation indicates that this methylation is not necessary for embryogenesis, even in species that require RdDM for seed development.

Methods

Plant materials and growth conditions

Brassica rapa ssp trilocularis variety R-o-18 were grown in a greenhouse at 70°/60°F (day/night) under at least 16 h of illumination. Plants were fully dried before seed collection. Dry seeds were soaked in water for no more than 60 min before manual dissection to remove mature embryos. Three wild-type or five rdr2 mutant embryos were pooled before DNA extraction with the GeneJET Plant Genomic DNA Purification Kit (Thermo Fisher Scientific, K0791). Embryos from rdr2/RDR2 heterozygous mothers were individually collected, prepped, and genotyped prior to DNA pooling. rdr2 mutants were used in this study due to their complete loss of 24-nt siRNAs and strong developmental phenotype. Torpedo-stage endosperm and embryo samples were dissected from pistils that were manually pollinated with B. rapa genotype R500. Whole-genome bisulfite sequencing libraries were prepared as previously described [44]. Lambda Phage DNA (Promega D1521) was included as a bisulfite conversion control. Libraries were pooled and sequenced in a single lane of paired-end 76 nt on an Illumina NextSeq500 at the University of Arizona Genetics Core.

Methylation analysis

Whole-genome bisulfite sequencing data from the ovule and leaf were obtained from NCBI (BioProject PRJNA588293, Additional file 1: Table S1). For other tissues, sequencing reads were quality controlled with FastQC [45, 46] and trimmed using Trim Galore (options --trim-n and --quality 20) [47]. Trimmed reads were aligned to Brassica rapa R-o-18 genome (v2.2, a kind gift from G.J. King and the B. rapa sequencing consortium) with bwameth [48]. To mark PCR duplicates and determine properly paired the alignment rate, Picard Tools [49] and Samtools [50] were respectively used, with options -q 10, -c, -F 3840, -f 66 for Samtools. We used Mosdepth [51] with option -x and -Q 10 and a custom Python script developed previously in the lab (bed_coverage_to_x_coverage.py, https://github.com/The-Mosher-Lab/grover_et_al_sirens_2020) to determine genomic coverage. Statistics for all libraries are found in Additional file 1: Table S1.

Percentage methylation per cytosine was extracted with MethylDackel [52] in two successive steps. The first step was to identify inclusion bounds based on methylation bias per read position using MethylDackel mbias, followed by MethylDackel extract. Since the default for MethylDackel is the CG context, we also used --CHG and --CHH options. We determined bisulfite conversion rates by alignment to the bacteriophage lambda (NCBI Genbank accession J02459.1) and Brassica rapa var. pekinensis chloroplast (NCBI Genbank accession NC 015139.1) genomes with a custom Python script developed previously in the lab (bedgraph_bisulfite_conv_calc.py, https://github.com/The-Mosher-Lab/grover_et_al_sirens_2020). Conversion frequencies were all above 99.4% (Additional file 1: Table S1). Replicates were checked for consistency by principal component analysis before pooling to increase read depth (Additional file 1: Fig. S6).

Methylation was calculated for each sample on 300-bp non-overlapping windows, which was made with the help of BEDTools makewindows [53] feature on the Brassica rapa R-o-18 genome. Pairwise methylation differences between tissues were measured using the methylKit package on merged alignment files for each tissue. We considered only those windows which had a q value of < 0.005 when calculating differentially methylated windows.

We calculated the enrichment of genomic features like small RNA loci, TEs, and genes within the CHH hyper-DMWs using BEDTools intersect [53, 54]. Transposable elements were annotated as in [31]. We considered features to be overlapping if there was at least 1 nucleotide shared. Genomic features were annotated onto the 300-bp non-overlapping windows using BEDTools makewindows [53], and the number of overlaps and non-overlaps between the hyper-DMWs and the genomic features were recorded. Fisher’s exact test was performed in R to determine if the number of overlaps indicated significant enrichment or depletion.

Methylation over pre-defined siren loci [31] was determined with BEDTools intersect [53] and a custom Python script (bedgraph_methylation_by_bed.py) developed previously [31] in the lab.

Small RNA analysis

Small RNA sequencing datasets were obtained from NCBI (BioProject PRJNA588293, Additional file 1: Table S2). Small RNA processing (quality checking, non-coding RNA filtering, removal of reads mapping to chloroplast and mitochondrial genomes) was carried out with a publicly available small RNA data processing pipeline [55]. Only 19 to 26-nt reads were retained for further analysis. Replicates were pooled for better read alignment and depth. The genome was divided into 300-bp non-overlapping windows using BEDTools makewindows [53], and ShortStack [56, 57] was used to get read counts on each window (options --mismatches 0, --mmap u, --mincov 0.5 rpm, --pad 75 and --foldsize 1000). The sum of all 19-26nt small RNA reads from genomic windows or CHH hyper-DMWs were low and susceptible to count-based bias when normalized against total library size. Therefore, windows were binned into 5 sRNA expression levels and compared only within the same tissue.

Availability of data and materials

The datasets supporting the conclusions of this article are available in the National Center for Biotechnology Information Sequence Read Archive under Bio Projects PRJNA588293 [58] and PRJNA657007 [59]. Small RNA datasets including leaf, embryo, and endosperm samples were obtained from Bio Project PRJNA588293 [58]. Whole-genome bisulfite datasets (leaf, ovule endosperm, torpedo embryo) were obtained from PRJNA588293 [58]. Whole-genome bisulfite datasets from mature embryo samples have been deposited to Bio Project PRJNA657007 [59].

References

  1. Jones, PA. Functions of DNA methylation: islands, start sites, gene bodies and beyond. Nat Rev Genet 2012;13:484–92. https://doi.org/10.1038/nrg3230.

  2. Zhang, H., Lang, Z. & Zhu, JK. Dynamics and function of DNA methylation in plants. Nat Rev Mol Cell Biol 2018;19:489–506. https://doi.org/10.1038/s41580-018-0016-z.

  3. Matzke MA, Mosher RA. RNA-directed DNA methylation: an epigenetic pathway of increasing complexity. Nature Rev Genet. 2014;15(6):394–408. Nature Publishing Group; Available from:. https://doi.org/10.1038/nrg3683.

    Article  CAS  PubMed  Google Scholar 

  4. Zhong X, Hale CJ, Law JA, Johnson LM, Feng S, Tu A, et al. DDR complex facilitates global association of RNA polymerase v to promoters and evolutionarily young transposons. Nat Struct Mol Biol. Nat Struct Mol Biol. 2012;19(9):870–5. https://doi.org/10.1038/nsmb.2354.

    Article  CAS  PubMed  Google Scholar 

  5. Li Q, Gent JI, Zynda G, Song J, Makarevitch I, Hirsch CD, 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):14728–33. https://doi.org/10.1073/pnas.1514680112 National Academy of Sciences.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Blevins T, Podicheti R, Mishra V, Marasco M, Wang J, Rusch D, et al. Identification of pol IV and RDR2-dependent precursors of 24 nt siRNAs guiding de novo DNA methylation in arabidopsis. Elife. 2015;4:4. eLife Sciences Publications Ltd; Available from:. https://doi.org/10.7554/eLife.09591.

    Article  Google Scholar 

  7. Zhai J, Bischof S, Wang H, Feng S, Lee TF, Teng C, et al. A one precursor one siRNA model for pol IV-dependent siRNA biogenesis. Cell. 2015;163:445–55 Cell Press.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Singh J, Mishra V, Wang F, Huang HY, Pikaard CS. Reaction mechanisms of Pol IV, RDR2, and DCL3 drive RNA channeling in the siRNA-directed DNA methylation pathway. Mol Cell. 2019;75:576–589.e5 Cell Press.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Cao X, Jacobsen SE. Role of the Arabidopsis DRM methyltransferases in de novo DNA methylation and gene silencing. Curr Biol. 2002;12:1138–44 Cell Press.

    Article  CAS  PubMed  Google Scholar 

  10. Wierzbicki AT, Ream TS, Haag JR, Pikaard CS. RNA polymerase v transcription guides ARGONAUTE4 to chromatin. Nat Genet. 2009;41(5):630–4. https://doi.org/10.1038/ng.365 Nature Publishing Group.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Havecker ER, Wallbridge LM, Hardcastle TJ, Bush MS, Kelly KA, Dunn RM, et al. The Arabidopsis RNA-directed DNA methylation argonautes functionally diverge based on their expression and interaction with target loci. Plant Cell. 2010;22:321–34 American Society of Plant Biologists.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Ye R, Wang W, Iki T, Liu C, Wu Y, Ishikawa M, et al. Cytoplasmic assembly and selective nuclear import of Arabidopsis ARGONAUTE4/siRNA complexes. Mol Cell. 2012;46(6):859–70. https://doi.org/10.1016/j.molcel.2012.04.013 Cell Press.

    Article  CAS  PubMed  Google Scholar 

  13. Greaves I, Groszmann M, Dennis ES, Peacock WJ. Trans-chromosomal methylation. Epigenetics. 2012;7:800–5 Taylor and Francis Inc.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Hollick JB. Paramutation: a trans-homolog interaction affecting heritable gene regulation. Curr Opin Plant Biol. 2012;15(5):536–43. Elsevier Current Trends; Available from:. https://doi.org/10.1016/j.pbi.2012.09.003.

    Article  CAS  PubMed  Google Scholar 

  15. Aufsatz W, Mette MF, Van Der Winden J, Matzke AJM, Matzke M. RNA-directed DNA methylation in Arabidopsis. Proc Natl Acad Sci U S A. 2002;99(Supplement 4):16499–506. https://doi.org/10.1073/pnas.162371499 National Academy of Sciences.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Molnar A, Melnyk CW, Bassett A, Hardcastle TJ, Dunn R, Baulcombe DC. Small silencing RNAs in plants are mobile and direct epigenetic modification in recipient cells. Science. 2010;328(5980):872–5. https://doi.org/10.1126/science.1187959 American Association for the Advancement of Science.

    Article  CAS  PubMed  Google Scholar 

  17. Hollick JB. Paramutation and development. Annu Rev Cell Dev Biol. 2010;26(1):557–79. https://doi.org/10.1146/annurev.cellbio.042308.113400 Annual Reviews.

    Article  CAS  PubMed  Google Scholar 

  18. Gouil Q, Baulcombe DC. DNA methylation signatures of the plant chromomethyltransferases. Mittelsten Scheid O, editor. PLoS Genet. 2016;12:e1006526 Public Library of Science.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Grover JW, Kendall T, Baten A, Burgess D, Freeling M, King GJ, et al. Maternal components of RNA-directed DNA methylation are required for seed development in Brassica rapa. Plant J. 2018;94:575–82 Blackwell Publishing Ltd.

    Article  CAS  PubMed  Google Scholar 

  20. Chow HT, Chakraborty T, Mosher RA. RNA-directed DNA Methylation and sexual reproduction: expanding beyond the seed. Curr Opin Plant Biol. 2020; 54. 11–17. Elsevier Ltd; Available from: https://doi.org/10.1016/j.pbi.2019.11.006

  21. Xu L, Yuan K, Yuan M, Meng X, Chen M, Wu J, et al. Regulation of rice tillering by RNA-directed DNA methylation at miniature inverted-repeat transposable elements. Mol Plant. 2020;13(6):851–63. https://doi.org/10.1016/j.molp.2020.02.009 Cell Press.

    Article  CAS  PubMed  Google Scholar 

  22. Bouyer D, Kramdi A, Kassam M, Heese M, Schnittger A, Roudier F, et al. DNA methylation dynamics during early plant life. Genome Biol. 2017;18:179 BioMed Central Ltd.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Kawakatsu T, Nery JR, Castanon R, Ecker JR. Dynamic DNA methylation reconfiguration during seed development and germination. Genome Biol. 2017;18:171 BioMed Central Ltd.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Narsai R, Gouil Q, Secco D, Srivastava A, Karpievitch YV, Liew LC, et al. Extensive transcriptomic and epigenomic remodelling occurs during Arabidopsis thaliana germination. Genome Biol. 2017;18:172 BioMed Central Ltd.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Lin JY, Le BH, Chen M, Henry KF, Hur J, Hsieh TF, et al. Similarity between soybean and Arabidopsis seed methylomes and loss of non-CG methylation does not affect seed development. Proc Natl Acad Sci U S A. 2017;114(45):E9730–9. https://doi.org/10.1073/pnas.1716758114 National Academy of Sciences.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Rajkumar MS, Gupta K, Khemka NK, Garg R, Jain M. DNA methylation reprogramming during seed development and its functional relevance in seed size/weight determination in chickpea. Commun Biol. 2020;3:1–13 Nature Research.

    Article  Google Scholar 

  27. Ji L, Mathioni SM, Johnson S, Tucker D, Bewick AJ, Kim KD, et al. Genome-wide reinforcement of DNA methylation occurs during somatic embryogenesis in soybean. Plant Cell. 2019;31(10):2315–31. https://doi.org/10.1105/tpc.19.00255 American Society of Plant Biologists.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Bauer MJ, Fischer RL. Genome demethylation and imprinting in the endosperm. Curr Opin Plant Biol. 2011;14(2):162–7. https://doi.org/10.1016/j.pbi.2011.02.006.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Lafon-Placette C, Köhler C. Embryo and endosperm, partners in seed development. Curr Opin Plant Biol. 2014;17:64–9. https://doi.org/10.1016/j.pbi.2013.11.008.

    Article  PubMed  Google Scholar 

  30. Mosher RA, Melnyk CW. siRNAs and DNA methylation: seedy epigenetics. Trends Plant Sci. 2010;15(4):204–10. https://doi.org/10.1016/j.tplants.2010.01.002.

    Article  CAS  PubMed  Google Scholar 

  31. Grover JW, Burgess D, Kendall T, Baten A, Pokhrel S, King GJ, et al. Abundant expression of maternal siRNAs is a conserved feature of seed development. Proc Natl Acad Sci U S A. 2020;117:15305–15 NLM (Medline).

    Article  PubMed  PubMed Central  Google Scholar 

  32. Erdmann RM, Satyaki PRV, Klosinska M, Gehring M. A Small RNA pathway mediates allelic dosage in endosperm. Cell Rep. 2017;21(12):3364–72. https://doi.org/10.1016/j.celrep.2017.11.078.

    Article  CAS  PubMed  Google Scholar 

  33. Gehring M, Bubb KL, Henikoff S. Extensive demethylation of repetitive elements during seed development underlies gene imprinting. Science. 2009;324(5933):1447–51. https://doi.org/10.1126/science.1171609 American Association for the Advancement of Science.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Hsieh TF, Ibarra CA, Silva P, Zemach A, Eshed-Williams L, Fischer RL, et al. Genome-wide demethylation of Arabidopsis endosperm. Science. 2009;324(5933):1451–4. https://doi.org/10.1126/science.1172417.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Zemach A, Kim MY, Silva P, Rodrigues JA, Dotson B, Brooks MD, et al. Local DNA hypomethylation activates genes in rice endosperm. Proc Natl Acad Sci U S A. 2010;107(43):18729–34. https://doi.org/10.1073/pnas.1009695107 National Academy of Sciences.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Zemach A, Kim MY, Hsieh PH, Coleman-Derr D, Eshed-Williams L, Thao K, et al. The arabidopsis nucleosome remodeler DDM1 allows DNA methyltransferases to access H1-containing heterochromatin. Cell. 2013;153(1):193–205. https://doi.org/10.1016/j.cell.2013.02.033.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Wang Z, Baulcombe DC. Transposon age and non-CG methylation. Nat Commun. 2020;11:1–9 Nature Research.

    Google Scholar 

  38. Nowack MK, Ungru A, Bjerkan KN, Grini PE, Schnittger A. Reproductive cross-talk: seed development in flowering plants. Biochem Soc Trans. 2010;38(2):604–12.  https://doi.org/10.1042/BST0380604.

  39. Figueiredo DD, Köhler C. Signalling events regulating seed coat development. Biochem Soc Trans. 2014;42(2):358–63. https://doi.org/10.1042/BST20130221.

  40. Iwasaki YW, Siomi MC, Siomi H. PIWI-interacting RNA: its biogenesis and functions. Annu Rev Biochem. 2015;84:405–33 Annual Reviews Inc.

    Article  CAS  PubMed  Google Scholar 

  41. Martínez G, Panda K, Köhler C, Slotkin RK. Silencing in sperm cells is directed by RNA movement from the surrounding nurse cell. Nat Plants. 2016;2(4):16030. https://doi.org/10.1038/nplants.2016.30.

    Article  CAS  PubMed  Google Scholar 

  42. Martinez G, Wolff P, Wang Z, Moreno-Romero J, Santos-González J, Conze LL, et al. Paternal easiRNAs regulate parental genome dosage in Arabidopsis. Nat Genet. 2018;50(2):193–8. https://doi.org/10.1038/s41588-017-0033-4.

    Article  CAS  PubMed  Google Scholar 

  43. Papareddy RK, Páldi K, Paulraj S, Kao P, Nodine MD. Chromatin regulates bipartite-classified small RNA expression to maintain epigenome homeostasis in Arabidopsis; 2020. p. 2020.05.04.076885. Available from: https://www.biorxiv.org/content/10.1101/2020.05.04.076885v1.abstract. [cited 2020 Sep 2]

    Google Scholar 

  44. Urich MA, Nery JR, Lister R, Schmitz RJ, Ecker JR. MethylC-seq library preparation for base-resolution whole-genome bisulfite sequencing. Nat Protoc. 2015;10(3):475–83. https://doi.org/10.1038/nprot.2014.114 Nature Publishing Group.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Babraham Bioinformatics - FastQC a quality control tool for high throughput sequence data. Available from: https://www.bioinformatics.babraham.ac.uk/projects/fastqc/. [cited 2020 Sep 2].

  46. Andrews S. FastQC. Github; Available from: https://github.com/s-andrews/FastQC. [cited 2020 Sep 2].

  47. Babraham Bioinformatics - Trim Galore!. Available from: https://www.bioinformatics.babraham.ac.uk/projects/trim_galore/. [cited 2020 Sep 2].

  48. Pedersen BS, Eyring K, De S, Yang IV, Schwartz DA. Fast and accurate alignment of long bisulfite-seq reads. 2014; Available from: http://arxiv.org/abs/1401.1129

    Google Scholar 

  49. Picard Tools - By Broad Institute. Available from: http://broadinstitute.github.io/picard/. Last accessed 30 Apr 2021.

  50. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25(16):2078–9. https://doi.org/10.1093/bioinformatics/btp352 Oxford Academic.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Pedersen BS, Quinlan AR. Mosdepth: quick coverage calculation for genomes and exomes. Bioinformatics. 2018;34:867–8 Oxford University Press.

    Article  CAS  PubMed  Google Scholar 

  52. dpryan79/MethylDackel: A (mostly) universal methylation extractor for BS-seq experiments. Available from: https://github.com/dpryan79/methyldackel. Last accessed 30 Apr 2021.

  53. Quinlan AR. BEDTools: The Swiss-Army tool for genome feature analysis. Curr Protoc Bioinformatics. 2014;2014:11.12.1–11.12.34 Wiley.

    Google Scholar 

  54. Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26(6):841–2. https://doi.org/10.1093/bioinformatics/btq033 Oxford Academic.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. boseHere/sRNA_snakemake_workflow: a snakemake workflow for an sRNA-seq data analysis pipeline developed by the Mosher Lab. Available from: https://github.com/boseHere/sRNA_snakemake_workflow. Last accessed 30 Apr 2021.

  56. Axtell MJ. ShortStack: comprehensive annotation and quantification of small RNA genes. RNA. 2013;19:740–51 Cold Spring Harbor Laboratory Press.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Johnson NR, Yeoh JM, Coruh C. Axtell MJ. Improved placement of multi-mapping small RNAs. G3: Genes, Genomes, Genetics. 2016;6:2103–11 Genetics Society of America.

    Article  CAS  Google Scholar 

  58. Grover JW, Kendall T, and Mosher RA. Brassica rapa seed epigenome profiling, Accession PRJNA588293, (2019) https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA588293.

    Google Scholar 

  59. Chakraborty T, Kendall T, and Mosher RA. Embryo hypermethylation in Brassica rapa, Accession PRJNA657007, (2020) https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA657007.

    Google Scholar 

Download references

Acknowledgements

The authors are grateful to the Brassica rapa R-o-18 Genome Sequencing Consortium for the access to the B. rapa v2.3 genome and annotations and to Dr. Diane Burgess for transposable element annotation. Illumina sequencing was performed by the University of Arizona Genetics Core.

Peer review information

Kevin Pang was the primary editor of this article and managed its editorial process and peer review in collaboration with the rest of the editorial team.

Review history

The review history is available as Additional file 2.

Funding

The authors gratefully acknowledge support from the National Science Foundation (IOS-1546825 to RAM).

Author information

Authors and Affiliations

Authors

Contributions

RAM conceptualized the project. TC, TK, and JWG performed the experiments. TC and RAM analyzed and visualized the data and wrote the paper. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Rebecca A. Mosher.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

: Supplementary Figure S1-6, Supplementary Table S1-3.

Additional file 2.

Review history.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chakraborty, T., Kendall, T., Grover, J.W. et al. Embryo CHH hypermethylation is mediated by RdDM and is autonomously directed in Brassica rapa. Genome Biol 22, 140 (2021). https://doi.org/10.1186/s13059-021-02358-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13059-021-02358-3

Keywords