Background

5' methylation of cytosines in DNA molecules is an important epigenetic mark in eukaryotes. Bisulfite sequencing is the gold standard of DNA methylation detection, and whole-genome bisulfite sequencing (WGBS) has been widely used to detect methylation at single-nucleotide resolution on a genome-wide scale. However, sodium bisulfite is known to severely degrade DNA, which, in combination with biases introduced during PCR amplification, leads to unbalanced base representation in the final sequencing libraries. Enzymatic conversion of unmethylated cytosines to uracils can achieve the same end product for sequencing as does bisulfite treatment and does not affect the integrity of the DNA; enzymatic methylation sequencing may, thus, provide advantages over bisulfite sequencing.

Results

Using an enzymatic methyl-seq (EM-seq) technique to selectively deaminate unmethylated cytosines to uracils, we generated and sequenced libraries based on different amounts of Arabidopsis input DNA and different numbers of PCR cycles, and compared these data to results from traditional whole-genome bisulfite sequencing. We found that EM-seq libraries were more consistent between replicates and had higher mapping and lower duplication rates, lower background noise, higher average coverage, and higher coverage of total cytosines. Differential methylation region (DMR) analysis showed that WGBS tended to over-estimate methylation levels especially in CHG and CHH contexts, whereas EM-seq detected higher CG methylation levels in certain highly methylated areas. These phenomena can be mostly explained by a correlation of WGBS methylation estimation with GC content and methylated cytosine density. We used EM-seq to compare methylation between leaves and flowers, and found that CHG methylation level is greatly elevated in flowers, especially in pericentromeric regions.

Conclusion

We suggest that EM-seq is a more accurate and reliable approach than WGBS to detect methylation. Compared to WGBS, the results of EM-seq are less affected by differences in library preparation conditions or by the skewed base composition in the converted DNA. It may therefore be more desirable to use EM-seq in methylation studies.

DNA methylation and nucleosome densities play a critical role in the regulation of gene expression. While much is known about the mechanisms of transcriptional control that are mediated by these, less is known about the degree to which they are tissue-specific. By comparing DNA methylation, nucleosome densities and transcriptional levels in different tissue types we can gain a clearer understanding of the extent to which these mechanisms influence gene expression in a tissue specific manner. We compared DNA methylation in Arabidopsis shoots and roots and found extensive differences across the genome. We computed DNA methylation differences between roots and shoots at single cytosines and found that one in every 173 cytosines was differentially methylated. In addition we compared DNA methylation with tissue specific gene expression and nucleosome density measurements to identify associations between these. We also identified a group of genes that are strongly correlated with these epigenetic marks and are significantly differentially methylated between roots and shoots. These root-specific genes are part of the extensin family, and are preferentially methylated and have at least 10-fold higher expression and lower nucleosome density in roots relative to shoots.

Abstract Background We performed a comparative analysis of the genome-wide DNA methylation profiles from three human embryonic stem cell (HESC) lines. It had previously been shown that HESC lines had significantly higher non-CG methylation than differentiated cells, and we therefore asked whether these sites were conserved across cell lines. Results We find that heavily methylated non-CG sites are strongly conserved, especially when found within the motif TACAG. They are enriched in splice sites and are more methylated than other non-CG sites in genes. We next studied the relationship between allele-specific expression and allele-specific methylation. By combining bisulfite sequencing and whole transcriptome shotgun sequencing (RNA-seq) data we identified 1,020 genes that show allele-specific expression, and 14% of CG sites genome-wide have allele-specific methylation. Finally, we asked whether the methylation state of transcription factor binding sites affects the binding of transcription factors. We identified variations in methylation levels at binding sites and found that for several transcription factors the correlation between the methylation at binding sites and gene expression is generally stronger than in the neighboring sequences. Conclusions These results suggest a possible but as yet unknown functional role for the highly methylated conserved non-CG sites in the regulation of HESCs. We also identified a novel set of genes that are likely transcriptionally regulated by methylation in an allele-specific manner. The analysis of transcription factor binding sites suggests that the methylation state of cis-regulatory elements impacts the ability of factors to bind and regulate transcription.

CRISPR-based targeted modification of epigenetic marks such as DNA cytosine methylation is an important strategy to regulate the expression of genes and their associated phenotypes. Although plants have DNA methylation in all sequence contexts (CG, CHG, CHH, where H = A, T, C), methylation in the symmetric CG context is particularly important for gene silencing and is very efficiently maintained through mitotic and meiotic cell divisions. Tools that can directly add CG methylation to specific loci are therefore highly desirable but are currently lacking in plants. Here we have developed two CRISPR-based CG-specific targeted DNA methylation systems for plants using a variant of the bacterial CG-specific DNA methyltransferase MQ1 with reduced activity but high specificity. We demonstrate that the methylation added by MQ1 is highly target specific and can be heritably maintained in the absence of the effector. These tools should be valuable both in crop engineering and in plant genetic research.

Abstract Background 5-Hydroxymethylcytosine (5hmC) was recently found to be abundantly present in certain cell types, including embryonic stem cells. There is growing evidence that TET proteins, which convert 5-methylcytosine (5mC) to 5hmC, play important biological roles. To further understand the function of 5hmC, an analysis of the genome-wide localization of this mark is required. Results Here, we have generated a genome-wide map of 5hmC in human embryonic stem cells by hmeDIP-seq, in which hydroxymethyl-DNA immunoprecipitation is followed by massively parallel sequencing. We found that 5hmC is enriched in enhancers as well as in gene bodies, suggesting a potential role for 5hmC in gene regulation. Consistent with localization of 5hmC at enhancers, 5hmC was significantly enriched in histone modifications associated with enhancers, such as H3K4me1 and H3K27ac. 5hmC was also enriched in other protein-DNA interaction sites, such as OCT4 and NANOG binding sites. Furthermore, we found that 5hmC regions tend to have an excess of G over C on one strand of DNA. Conclusions Our findings suggest that 5hmC may be targeted to certain genomic regions based both on gene expression and sequence composition.

DNA methylation is an epigenetic mark that is associated with transcriptional repression of transposable elements and protein-coding genes. Conversely, transcriptionally active regulatory regions are strongly correlated with histone 3 lysine 4 di- and trimethylation (H3K4m2/m3). We previously showed that Arabidopsis thaliana plants with mutations in the H3K4m2/m3 demethylase JUMONJI 14 (JMJ14) exhibit a mild reduction in RNA-directed DNA methylation (RdDM) that is associated with an increase in H3K4m2/m3 levels. To determine whether this incomplete RdDM reduction was the result of redundancy with other demethylases, we examined the genetic interaction of JMJ14 with another class of H3K4 demethylases: lysine-specific demethylase 1-like 1 and lysine-specific demethylase 1-like 2 (LDL1 and LDL2). Genome-wide DNA methylation analyses reveal that both families cooperate to maintain RdDM patterns. ChIP-seq experiments show that regions that exhibit an observable DNA methylation decrease are co-incidental with increases in H3K4m2/m3. Interestingly, the impact on DNA methylation was stronger at DNA-methylated regions adjacent to H3K4m2/m3-marked protein-coding genes, suggesting that the activity of H3K4 demethylases may be particularly crucial to prevent spreading of active epigenetic marks. Finally, RNA sequencing analyses indicate that at RdDM targets, the increase of H3K4m2/m3 is not generally associated with transcriptional de-repression. This suggests that the histone mark itself--not transcription--impacts the extent of RdDM.

Chromosomes form 3D structures that are critical to the regulation of cellular and genetic processes. Here, we present a study of global chromatin interaction patterns in Arabidopsis thaliana. Our genome-wide approach confirmed interactions that were previously observed by other methods as well as uncovered long-range interactions such as those among small heterochromatic regions embedded in euchromatic arms. We also found that interactions are correlated with various epigenetic marks that are localized in active or silenced chromatin. Arabidopsis chromosomes do not contain large local interactive domains that resemble the topological domains described in animals but, instead, contain relatively small interactive regions scattered around the genome that contain H3K27me3 or H3K9me2. We generated interaction maps in mutants that are defective in specific epigenetic pathways and found altered interaction patterns that correlate with changes in the epigenome. These analyses provide further insights into molecular mechanisms of epigenetic regulation of the genome.

DNA methylation is a major epigenetic modification found across species and has a profound impact on many biological processes. However, its influence on chromatin accessibility and higher-order genome organization remains unclear, particularly in plants. Here, we present genome-wide chromatin accessibility profiles of 18 Arabidopsis mutants that are deficient in CG, CHG, or CHH DNA methylation. We find that DNA methylation in all three sequence contexts impacts chromatin accessibility in heterochromatin. Many chromatin regions maintain inaccessibility when DNA methylation is lost in only one or two sequence contexts, and signatures of accessibility are particularly affected when DNA methylation is reduced in all contexts, suggesting an interplay between different types of DNA methylation. In addition, we found that increased chromatin accessibility was not always accompanied by increased transcription, suggesting that DNA methylation can directly impact chromatin structure by other mechanisms. We also observed that an increase in chromatin accessibility was accompanied by enhanced long-range chromatin interactions. Together, these results provide a valuable resource for chromatin architecture and DNA methylation analyses and uncover a pivotal role for methylation in the maintenance of heterochromatin inaccessibility.

The plant-specific DNA-dependent RNA polymerase V (Pol V) evolved from Pol II to function in an RNA-directed DNA methylation pathway. Here, we have identified targets of Pol V in Arabidopsis thaliana on a genome-wide scale using ChIP-seq of NRPE1, the largest catalytic subunit of Pol V. We found that Pol V is enriched at promoters and evolutionarily recent transposons. This localization pattern is highly correlated with Pol V-dependent DNA methylation and small RNA accumulation. We also show that genome-wide chromatin association of Pol V is dependent on all members of a putative chromatin-remodeling complex termed DDR. Our study presents a genome-wide view of Pol V occupancy and sheds light on the mechanistic basis of Pol V localization. Furthermore, these findings suggest a role for Pol V and RNA-directed DNA methylation in genome surveillance and in responding to genome evolution.

DNA methylation occurs in CG and non-CG sequence contexts. Non-CG methylation is abundant in plants and is mediated by CHROMOMETHYLASE (CMT) and DOMAINS REARRANGED METHYLTRANSFERASE (DRM) proteins; however, its roles remain poorly understood. Here we characterize the roles of non-CG methylation in Arabidopsis thaliana. We show that a poorly characterized methyltransferase, CMT2, is a functional methyltransferase in vitro and in vivo. CMT2 preferentially binds histone H3 Lys9 (H3K9) dimethylation and methylates non-CG cytosines that are regulated by H3K9 methylation. We revealed the contributions and redundancies between each non-CG methyltransferase in DNA methylation patterning and in regulating transcription. We also demonstrate extensive dependencies of small-RNA accumulation and H3K9 methylation patterning on non-CG methylation, suggesting self-reinforcing mechanisms between these epigenetic factors. The results suggest that non-CG methylation patterns are critical in shaping the landscapes of histone modification and small noncoding RNA.