Plant-RRBS, a bisulfite and next-generation sequencing-based methylome profiling method enriching for coverage of cytosine positions

Background: Cytosine methylation in plant genomes is important for the regulation of gene transcription and transposon activity. Genome-wide methylomes are studied upon mutation of the DNA methyltransferases, adaptation to environmental stresses or during development. However, from basic biology to breeding programs, there is a need to monitor multiple samples to determine transgenerational methylation inheritance or differential cytosine methylation. Methylome data obtained by sodium hydrogen sulfite (bisulfite)-conversion and next-generation sequencing (NGS) provide genome- wide information on cytosine methylation. However, a profiling method that detects cytosine methylation state dispersed over the genome would allow high-throughput analysis of multiple plant samples with distinct epigenetic signatures. We use specific restriction endonucleases to enrich for cytosine coverage in a bisulfite and NGS-based profiling method, which was compared to whole-genome bisulfite sequencing of the same plant material. Methods: We established an effective methylome profiling method in plants, termed plant-reduced representation bisulfite sequencing (plant-RRBS), using optimized double restriction endonuclease digestion, fragment end repair, adapter ligation, followed by bisulfite conversion, PCR amplification and NGS. We report a performant laboratory protocol and a straightforward bioinformatics data analysis pipeline for plant-RRBS, applicable for any reference-sequenced plant species. Results: As a proof of concept, methylome profiling was performed using an Oryza sativa ssp. indica pure breeding line and a derived epigenetically altered line (epiline). Plant-RRBS detects methylation levels at tens of millions of cytosine positions deduced from bisulfite conversion in multiple samples. To evaluate the method, the coverage of cytosine positions, the intra-line similarity and the differential cytosine methylation levels between the pure breeding line and the epiline were determined. Plant-RRBS reproducibly covers commonly up to one fourth of the cytosine positions in the rice genome when using MspI-DpnII within a group of five biological replicates of a line. The method predominantly detects cytosine methylation in putative promoter regions and not-annotated regions in rice. Conclusions: Plant-RRBS offers high-throughput and broad, genome- dispersed methylation detection by effective read number generation obtained from reproducibly covered genome fractions using optimized endonuclease combinations, facilitating comparative analyses of multi-sample studies for cytosine methylation and transgenerational stability in experimental material and plant breeding populations

Detection of Oxidation Products of 5-Methyl-2′-Deoxycytidine in Arabidopsis DNA

Epigenetic regulations play important roles in plant development and adaptation to environmental stress. Recent studies from mammalian systems have demonstrated the involvement of ten-eleven translocation (Tet) family of dioxygenases in the generation of a series of oxidized derivatives of 5-methylcytosine (5-mC) in mammalian DNA. In addition, these oxidized 5-mC nucleobases have important roles in epigenetic remodeling and aberrant levels of 5-hydroxymethyl-29-deoxycytidine (5-HmdC) were found to be associated with different types of human cancers. However, there is a lack of evidence supporting the presence of these modified bases in plant DNA. Here we reported the use of a reversed-phase HPLC coupled with tandem mass spectrometry method and stable isotope-labeled standards for assessing the levels of the oxidized 5-mC nucleosides along with two other oxidatively induced DNA modifications in genomic DNA of Arabidopsis. These included 5- HmdC, 5-formyl-29-deoxycytidine (5-FodC), 5-carboxyl-29-deoxycytidine (5-CadC), 5-hydroxymethyl-29-deoxyuridine (5- HmdU), and the (59S) diastereomer of 8,59-cyclo-29-deoxyguanosine (S-cdG). We found that, in Arabidopsis DNA, the levels of 5-HmdC, 5-FodC, and 5-CadC are approximately 0.8 modifications per 106 nucleosides, with the frequency of 5-HmdC (per 5-mdC) being comparable to that of 5-HmdU (per thymidine). The relatively low levels of the 5-mdC oxidation products suggest that they arise likely from reactive oxygen species present in cells, which is in line with the lack of homologous Tetfamily dioxygenase enzymes in Arabidopsis

AGO6 Functions in RNA-Mediated Transcriptional Gene Silencing in Shoot and Root Meristems in Arabidopsis thaliana

RNA-directed DNA methylation (RdDM) is a small interfering RNA (siRNA)-mediated epigenetic modification that contributes to transposon silencing in plants. RdDM requires a complex transcriptional machinery that includes specialized RNA polymerases, named Pol IV and Pol V, as well as chromatin remodelling proteins, transcription factors, RNA binding proteins, and other plant-specific proteins whose functions are not yet clarified. In Arabidopsis thaliana, DICER-LIKE3 and members of the ARGONAUTE4 group of ARGONAUTE (AGO) proteins are involved, respectively, in generating and using 24-nt siRNAs that trigger methylation and transcriptional gene silencing of homologous promoter sequences. AGO4 is the main AGO protein implicated in the RdDM pathway. Here we report the identification of the related AGO6 in a forward genetic screen for mutants defective in RdDM and transcriptional gene silencing in shoot and root apical meristems in Arabidopsis thaliana. The identification of AGO6, and not AGO4, in our screen is consistent with the primary expression of AGO6 in shoot and root growing points

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

Transposable elements (TEs) comprise a substantial portion of many eukaryotic genomes and are typically transcriptionally silenced. RNA–dependent RNA polymerase 2 (RDR2) is a component of the RNA–directed DNA methylation (RdDM) silencing pathway. In maize, loss of mediator of paramutation1 (mop1) encoded RDR2 function results in reactivation of transcriptionally silenced Mu transposons and a substantial reduction in the accumulation of 24 nt short-interfering RNAs (siRNAs) that recruit RNA silencing components. An RNA–seq experiment conducted on shoot apical meristems (SAMs) revealed that, as expected based on a model in which RDR2 generates 24 nt siRNAs that suppress expression, most differentially expressed DNA TEs (78%) were up-regulated in the mop1 mutant. In contrast, most differentially expressed retrotransposons (68%) were down-regulated. This striking difference suggests that distinct silencing mechanisms are applied to different silencing templates. In addition, >6,000 genes (24% of analyzed genes), including nearly 80% (286/361) of genes in chromatin modification pathways, were differentially expressed. Overall, two-thirds of differentially regulated genes were down-regulated in the mop1 mutant. This finding suggests that RDR2 plays a significant role in regulating the expression of not only transposons, but also of genes. A re-analysis of existing small RNA data identified both RDR2–sensitive and RDR2–resistant species of 24 nt siRNAs that we hypothesize may at least partially explain the complex changes in the expression of genes and transposons observed in the mop1 mutant

DNA Methylation Causes Predominant Maternal Controls of Plant Embryo Growth

The parental conflict hypothesis predicts that the mother inhibits embryo growth counteracting growth enhancement by the father. In plants the DNA methyltransferase MET1 is a central regulator of parentally imprinted genes that affect seed growth. However the relation between the role of MET1 in imprinting and its control of seed size has remained unclear. Here we combine cytological, genetic and statistical analyses to study the effect of MET1 on seed growth. We show that the loss of MET1 during male gametogenesis causes a reduction of seed size, presumably linked to silencing of the paternal allele of growth enhancers in the endosperm, which nurtures the embryo. However, we find no evidence for a similar role of MET1 during female gametogenesis. Rather, the reduction of MET1 dosage in the maternal somatic tissues causes seed size increase. MET1 inhibits seed growth by restricting cell division and elongation in the maternal integuments that surround the seed. Our data demonstrate new controls of seed growth linked to the mode of reproduction typical of flowering plants. We conclude that the regulation of embryo growth by MET1 results from a combination of predominant maternal controls, and that DNA methylation maintained by MET1 does not orchestrate a parental conflict

The Tnt1 Retrotransposon Escapes Silencing in Tobacco, Its Natural Host

Retrotransposons' high capacity for mutagenesis is a threat that genomes need to control tightly. Transcriptional gene silencing is a general and highly effective control of retrotransposon expression. Yet, some retrotransposons manage to transpose and proliferate in plant genomes, suggesting that, as shown for plant viruses, retrotransposons can escape silencing. However no evidence of retrotransposon silencing escape has been reported. Here we analyze the silencing control of the tobacco Tnt1 retrotransposon and report that even though constructs driven by the Tnt1 promoter become silenced when stably integrated in tobacco, the endogenous Tnt1 elements remain active. Silencing of Tnt1-containing transgenes correlates with high DNA methylation and the inability to incorporate H2A.Z into their promoters, whereas the endogenous Tnt1 elements remain partially methylated at asymmetrical positions and incorporate H2A.Z upon induction. Our results show that the promoter of Tnt1 is a target of silencing in tobacco, but also that endogenous Tnt1 elements can escape this control and be expressed in their natural host

DNA Methylation and Normal Chromosome Behavior in Neurospora Depend on Five Components of a Histone Methyltransferase Complex, DCDC

Methylation of DNA and of Lysine 9 on histone H3 (H3K9) is associated with gene silencing in many animals, plants, and fungi. In Neurospora crassa, methylation of H3K9 by DIM-5 directs cytosine methylation by recruiting a complex containing Heterochromatin Protein-1 (HP1) and the DIM-2 DNA methyltransferase. We report genetic, proteomic, and biochemical investigations into how DIM-5 is controlled. These studies revealed DCDC, a previously unknown protein complex including DIM-5, DIM-7, DIM-9, CUL4, and DDB1. Components of DCDC are required for H3K9me3, proper chromosome segregation, and DNA methylation. DCDC-defective strains, but not HP1-defective strains, are hypersensitive to MMS, revealing an HP1-independent function of H3K9 methylation. In addition to DDB1, DIM-7, and the WD40 domain protein DIM-9, other presumptive DCAFs (DDB1/CUL4 associated factors) co-purified with CUL4, suggesting that CUL4/DDB1 forms multiple complexes with distinct functions. This conclusion was supported by results of drug sensitivity tests. CUL4, DDB1, and DIM-9 are not required for localization of DIM-5 to incipient heterochromatin domains, indicating that recruitment of DIM-5 to chromatin is not sufficient to direct H3K9me3. DIM-7 is required for DIM-5 localization and mediates interaction of DIM-5 with DDB1/CUL4 through DIM-9. These data support a two-step mechanism for H3K9 methylation in Neurospora

Small RNAs Prevent Transcription-Coupled Loss of Histone H3 Lysine 9 Methylation in Arabidopsis thaliana

In eukaryotes, histone H3 lysine 9 methylation (H3K9me) mediates silencing of invasive sequences to prevent deleterious consequences including the expression of aberrant gene products and mobilization of transposons. In Arabidopsis thaliana, H3K9me maintained by SUVH histone methyltransferases (MTases) is associated with cytosine methylation (5meC) maintained by the CMT3 cytosine MTase. The SUVHs contain a 5meC binding domain and CMT3 contains an H3K9me binding domain, suggesting that the SUVH/CMT3 pathway involves an amplification loop between H3K9me and 5meC. However, at loci subject to read-through transcription, the stability of the H3K9me/5meC loop requires a mechanism to counteract transcription-coupled loss of H3K9me. Here we use the duplicated PAI genes, which stably maintain SUVH-dependent H3K9me and CMT3-dependent 5meC despite read-through transcription, to show that when PAI sRNAs are depleted by dicer ribonuclease mutations, PAI H3K9me and 5meC levels are reduced and remaining PAI 5meC is destabilized upon inbreeding. The dicer mutations confer weaker reductions in PAI 5meC levels but similar or stronger reductions in PAI H3K9me levels compared to a cmt3 mutation. This comparison indicates a connection between sRNAs and maintenance of H3K9me independent of CMT3 function. The dicer mutations reduce PAI H3K9me and 5meC levels through a distinct mechanism from the known role of dicer-dependent sRNAs in guiding the DRM2 cytosine MTase because the PAI genes maintain H3K9me and 5meC at levels similar to wild type in a drm2 mutant. Our results support a new role for sRNAs in plants to prevent transcription-coupled loss of H3K9me

Three SRA-Domain Methylcytosine-Binding Proteins Cooperate to Maintain Global CpG Methylation and Epigenetic Silencing in Arabidopsis

Methylcytosine-binding proteins decipher the epigenetic information encoded by DNA methylation and provide a link between DNA methylation, modification of chromatin structure, and gene silencing. VARIANT IN METHYLATION 1 (VIM1) encodes an SRA (SET- and RING-associated) domain methylcytosine-binding protein in Arabidopsis thaliana, and loss of VIM1 function causes centromere DNA hypomethylation and centromeric heterochromatin decondensation in interphase. In the Arabidopsis genome, there are five VIM genes that share very high sequence similarity and encode proteins containing a PHD domain, two RING domains, and an SRA domain. To gain further insight into the function and potential redundancy among the VIM proteins, we investigated strains combining different vim mutations and transgenic vim knock-down lines that down-regulate multiple VIM family genes. The vim1 vim3 double mutant and the transgenic vim knock-down lines showed decreased DNA methylation primarily at CpG sites in genic regions, as well as repeated sequences in heterochromatic regions. In addition, transcriptional silencing was released in these plants at most heterochromatin regions examined. Interestingly, the vim1 vim3 mutant and vim knock-down lines gained ectopic CpHpH methylation in the 5S rRNA genes against a background of CpG hypomethylation. The vim1 vim2 vim3 triple mutant displayed abnormal morphological phenotypes including late flowering, which is associated with DNA hypomethylation of the 5′ region of FWA and release of FWA gene silencing. Our findings demonstrate that VIM1, VIM2, and VIM3 have overlapping functions in maintenance of global CpG methylation and epigenetic transcriptional silencing

H3K27me3 Profiling of the Endosperm Implies Exclusion of Polycomb Group Protein Targeting by DNA Methylation

Polycomb group (PcG) proteins act as evolutionary conserved epigenetic mediators of cell identity because they repress transcriptional programs that are not required at particular developmental stages. Each tissue is likely to have a specific epigenetic profile, which acts as a blueprint for its developmental fate. A hallmark for Polycomb Repressive Complex 2 (PRC2) activity is trimethylated lysine 27 on histone H3 (H3K27me3). In plants, there are distinct PRC2 complexes for vegetative and reproductive development, and it was unknown so far whether these complexes have target gene specificity. The FERTILIZATION INDEPENDENT SEED (FIS) PRC2 complex is specifically expressed in the endosperm and is required for its development; loss of FIS function causes endosperm hyperproliferation and seed abortion. The endosperm nourishes the embryo, similar to the physiological function of the placenta in mammals. We established the endosperm H3K27me3 profile and identified specific target genes of the FIS complex with functional roles in endosperm cellularization and chromatin architecture, implicating that distinct PRC2 complexes have a subset of specific target genes. Importantly, our study revealed that selected transposable elements and protein coding genes are specifically targeted by the FIS PcG complex in the endosperm, whereas these elements and genes are densely marked by DNA methylation in vegetative tissues, suggesting that DNA methylation prevents targeting by PcG proteins in vegetative tissues