Mouse Oocyte Methylomes at Base Resolution Reveal Genome-Wide Accumulation of Non-CpG Methylation and Role of DNA Methyltransferases

DNA methylation is an epigenetic modification that plays a crucial role in normal mammalian development, retrotransposon silencing, and cellular reprogramming. Although methylation mainly occurs on the cytosine in a CG site, non-CG methylation is prevalent in pluripotent stem cells, brain, and oocytes. We previously identified non-CG methylation in several CG-rich regions in mouse germinal vesicle oocytes (GVOs), but the overall distribution of non-CG methylation and the enzymes responsible for this modification are unknown. Using amplification-free whole-genome bisulfite sequencing, which can be used with minute amounts of DNA, we constructed the base-resolution methylome maps of GVOs, non-growing oocytes (NGOs), and mutant GVOs lacking the DNA methyltransferase Dnmt1, Dnmt3a, Dnmt3b, or Dnmt3L. We found that nearly two-thirds of all methylcytosines occur in a non-CG context in GVOs. The distribution of non-CG methylation closely resembled that of CG methylation throughout the genome and showed clear enrichment in gene bodies. Compared to NGOs, GVOs were over four times more methylated at non-CG sites, indicating that non-CG methylation accumulates during oocyte growth. Lack of Dnmt3a or Dnmt3L resulted in a global reduction in both CG and non-CG methylation, showing that non-CG methylation depends on the Dnmt3a-Dnmt3L complex. Dnmt3b was dispensable. Of note, lack of Dnmt1 resulted in a slight decrease in CG methylation, suggesting that this maintenance enzyme plays a role in non-dividing oocytes. Dnmt1 may act on CG sites that remain hemimethylated in the de novo methylation process. Our results provide a basis for understanding the mechanisms and significance of non-CG methylation in mammalian oocytes.     

Cooperative activity of DNA methyltransferases for maintenance of symmetrical and non-symmetrical cytosine methylation in Arabidopsis thaliana

Maintenance of cytosine methylation in plants is controlled by three DNA methyltransferases. MET1 maintains CG methylation, and DRM1/2 and CMT3 act redundantly to enforce non-CG methylation. RPS, a repetitive hypermethylated DNA fragment from Petunia hybrida, attracts DNA methylation when transferred into Petunia or other species. In Arabidopsis thaliana, which does not contain any RPS homologues, RPS transgenes are efficiently methylated in all sequence contexts. To test which DNA methylation pathways regulate RPS methylation, we examined maintenance of RPS methylation in various mutant backgrounds. Surprisingly, CG methylation was lost in a drm1/2/cmt3 mutant, and non-CG methylation was almost completely eliminated in a met1 mutant. An unusual cooperative activity of all three DNA methyltransferases is therefore required for maintenance of both CG and non-CG methylation in RPS. Other unusual features of RPS methylation are the independence of its non-CG methylation from the RNA-directed DNA methylation (RdDM) pathway and the exceptional maintenance of methylation at a CC(m)TGG site in some epigenetic mutants. This is indicative of activity of a methylation system in plants that may have evolved from the DCM methylation system that controls CC(m)WGG methylation in bacteria. Our data suggest that strict separation of CG and non-CG methylation pathways does not apply to all target regions, and that caution is required in generalizing methylation data obtained for individual genomic regions.     

DNA methylation in mouse A-repeats in DNA methyltransferase-knockout ES cells and in normal cells determined by bisulfite genomic sequencing

Mouse ES cells with a null mutation of the known DNA methyltransferase retain some residual DNA methylation and can methylate foreign sequences de novo. We have used bisulfite genomic sequencing to examine the sequence specificity and distributions of methylation of a hypermethylated CG island sequence, mouse A-repeats. There were 13 CG dinucleotides in the region examined, 12 of which were methylated to variable extents in all DNAs. We found that: (1) there is considerable residual DNA methylation in ES cells lacking the known DNA methyltransferase (29% of normal methylation in the complete knockout ES DNA); (2) this other activity methylates at exactly the same CG sites as the major methyltransferase; and (3) differences in the distribution of methylated sites between A-repeats in these DNAs are consistent with this other activity methylating in a random de novo fashion. Also, the lack of any methylation in non-CG sites argues that, in other studies where non-CG methylation sites have been found by bisulfite sequencing, detection of such sites of non-CG methylation is not an inherent artifact in this methodology.     

RNAi, DRD1, and histone methylation actively target developmentally important non-CG DNA methylation in arabidopsis

Cytosine DNA methylation protects eukaryotic genomes by silencing transposons and harmful DNAs, but also regulates gene expression during normal development. Loss of CG methylation in the Arabidopsis thaliana met1 and ddm1 mutants causes varied and stochastic developmental defects that are often inherited independently of the original met1 or ddm1 mutation. Loss of non-CG methylation in plants with combined mutations in the DRM and CMT3 genes also causes a suite of developmental defects. We show here that the pleiotropic developmental defects of drm1 drm2 cmt3 triple mutant plants are fully recessive, and unlike phenotypes caused by met1 and ddm1, are not inherited independently of the drm and cmt3 mutations. Developmental phenotypes are also reversed when drm1 drm2 cmt3 plants are transformed with DRM2 or CMT3, implying that non-CG DNA methylation is efficiently re-established by sequence-specific signals. We provide evidence that these signals include RNA silencing though the 24-nucleotide short interfering RNA (siRNA) pathway as well as histone H3K9 methylation, both of which converge on the putative chromatin-remodeling protein DRD1. These signals act in at least three partially intersecting pathways that control the locus-specific patterning of non-CG methylation by the DRM2 and CMT3 methyltransferases. Our results suggest that non-CG DNA methylation that is inherited via a network of persistent targeting signals has been co-opted to regulate developmentally important genes.     

Human non-CG methylation: Are human stem cells plant-like?

Non-CG methylation is well characterized in plants where it appears to play a role in gene silencing and genomic imprinting. Although strong evidence for the presence of non-CG methylation in mammals has been available for some time, both its origin and function remain elusive. In this review we discuss available evidence on non-CG methylation in mammals in light of evidence suggesting that the human stem cell methylome contains significant levels of methylation outside the CG site.Key words: non-CG methylation, stem cells, Dnmt1, Dnmt3a, human methylomeIn plant cells non-CG sites are methylated de novo by Chromomethylase 3, DRM1 and DRM2. Chromomethylase 3, along with DRM1 and DRM2 combine in the maintenance of methylation at symmetric CpHpG as well as asymmetric DNA sites where they appear to prevent reactivation of transposons.¹ DRM1 and DRM2 modify DNA de novo primarily at asymmetric CpH and CpHpH sequences targeted by siRNA.²Much less information is available on non-CG methylation in mammals. In fact, studies on mammalian non-CG methylation form a tiny fraction of those on CG methylation, even though data for cytosine methylation in other dinucleotides, CA, CT and CC, have been available since the late 1980s.³ Strong evidence for non-CG methylation was found by examining either exogenous DNA sequences, such as plasmid and viral integrants in mouse and human cell lines,^4,5 or transposons and repetitive sequences such as the human L1 retrotransposon⁶ in a human embryonic fibroblast cell line. In the latter study, non-CG methylation observed in L1 was found to be consistent with the capacity of Dnmt1 to methylate slippage intermediates de novo.⁶Non-CG methylation has also been reported at origins of replication^7,8 and a region of the human myogenic gene Myf3.⁹ The Myf3 gene is silenced in non-muscle cell lines but it is not methylated at CGs. Instead, it carries several methylated cytosines within the sequence CCTGG. Gene-specific non-CG methylation was also reported in a study of lymphoma and myeloma cell lines not expressing many B lineage-specific genes.¹⁰ The study focused on one specific gene, B29 and found heavy CG promoter methylation of that gene in most cell lines not expressing it. However, in two other cell lines where the gene was silenced, cytosine methylation was found almost exclusively at CCWGG sites. The authors provided evidence suggesting that CCWGG methylation was sufficient for silencing the B29 promoter and that methylated probes based on B29 sequences had unique gel shift patterns compared to non-methylated but otherwise identical sequences.¹⁰ The latter finding suggests that the presence of the non-CG methylation causes changes in the proteins able to bind the promoter, which could be mechanistically related to the silencing seen with this alternate methylation.Non-CG methylation is rarely seen in DNA isolated from cancer patients. However, the p16 promoter region was reported to contain both CG and non-CG methylation in breast tumor specimens but lacked methylation at these sites in normal breast tissue obtained at mammoplasty.¹¹ Moreover, CWG methylation at the CCWGG sites in the calcitonin gene is not found in normal or leukemic lymphocyte DNA obtained from patients.¹² Further, in DNA obtained from breast cancer patients, MspI sites that are refractory to digestion by MspI and thus candidates for CHG methylation were found to carry CpG methylation.¹³ Their resistance to MspI restriction was found to be caused by an unusual secondary structure in the DNA spanning the MspI site that prevents restriction.¹³ This latter observation suggests caution in interpreting EcoRII/BstNI or EcoRII/BstOI restriction differences as due to CWG methylation, since in contrast to the 37°C incubation temperature required for full EcoRII activity, BstNI and BstOI require incubation at 60°C for full activity where many secondary structures are unstable.The recent report by Lister et al.¹⁴ confirmed a much earlier report by Ramsahoye et al.¹⁵ suggesting that non-CG methylation is prevalent in mammalian stem cell lines. Nearest neighbor analysis was used to detect non-CG methylation in the earlier study on the mouse embryonic stem (ES) cell line,¹⁵ thus global methylation patterning was assessed. Lister et al.¹⁴ extend these findings to human stem cell lines at single-base resolution with whole-genome bisulfite sequencing. They report¹⁴ that the methylome of the human H1 stem cell line and the methylome of the induced pluripotent IMR90 (iPS) cell line are stippled with non-CG methylation while that of the human IMR90 fetal fibroblast cell line is not. While the results of the two studies are complementary, the human methylome study addresses locus specific non-CG methylation. Based on that data,¹⁴ one must conclude that non-CG methylation is not carefully maintained at a given site in the human H1 cell line. The average non-CG site is picked up as methylated in about 25% of the reads whereas the average CG methylation site is picked up in 92% of the reads. Moreover, non-CG methylation is not generally present on both strands and is concentrated in the body of actively transcribed genes.¹⁴Even so, the consistent finding that non-CG methylation appears to be confined to stem cell lines,^14,15 raises the possibility that cancer stem cells¹⁶ carry non-CG methylation while their nonstem progeny in the tumor carry only CG methylation. Given the expected paucity of cancer stem cells in a tumor cell population, it is unlikely that bisulfite sequencing would detect non-CG methylation in DNA isolated from tumor cells since the stem cell population is expected to be only a very minor component of tumor DNA. Published sequences obtained by bisulfite sequencing generally report only CG methylation, and to the best of our knowledge bisulfite sequenced tumor DNA specimens have not reported non-CG methylation. On the other hand, when sequences from cell lines have been reported, bisulfite-mediated genomic sequencing⁸ or ligation mediated PCR¹⁷ methylcytosine signals outside the CG site have been observed. In a more recent study plasmid DNAs carrying the Bcl2-major breakpoint cluster¹⁸ or human breast cancer DNA¹³ treated with bisulfite under non-denaturing conditions, cytosines outside the CG side were only partially converted on only one strand¹⁸ or at a symmetrical CWG site.¹³ In the breast cancer DNA study the apparent CWG methylation was not detected when the DNA was fully denatured before bisulfite treatment.¹³In both stem cell studies, non-CG methylation was attributed to the Dnmt3a,^14,15 a DNA methyltransferase with similarities to the plant DRM methyltransferase family¹⁹ and having the capacity to methylate non-CG sites when expressed in Drosophila melanogaster.¹⁵ DRM proteins however, possess a unique permuted domain structure found exclusively in plants¹⁹ and the associated RNA-directed non-CG DNA methylation has not been reproducibly observed in mammals despite considerable published^20–23 and unpublished efforts in that area. Moreover, reports where methylation was studied often infer methylation changes from 5AzaC reactivation studies²⁴ or find that CG methylation seen in plants but not non-CG methylation is detected.^21,22,25,26 In this regard, it is of interest that the level of non-CG methylation reported in stem cells corresponds to background non-CG methylation observed in vitro with human DNA methyltransferase I,²⁷ and is consistent with the recent report that cultured stem cells are epigenetically unstable.²⁸The function of non-CG methylation remains elusive. A role in gene expression has not been ruled out, as the studies above on Myf3 and B29 suggest.^9,10 However, transgene expression of the bacterial methyltransferase M.EcoRII in a human cell line (HK293), did not affect the CG methylation state at the APC and SerpinB5 genes²⁹ even though the promoters were symmetrically de novo methylated at ^mCWGs within each CCWGG sequence in each promoter. This demonstrated that CG and non-CG methylation are not mutually exclusive as had been suggested by earlier reports.^9,10 That observation is now extended to the human stem cell line methylome where CG and non-CG methylation co-exist.¹⁴ Gene expression at the APC locus was likewise unaffected by transgene expression of M.EcoRII. In those experiments genome wide methylation of the CCWGG site was detected by restriction analysis and bisulfite sequencing,²⁹ however stem cell characteristics were not studied.Many alternative functions can be envisioned for non-CG methylation, but the existing data now constrains them to functions that involve low levels of methylation that are primarily asymmetric. Moreover, inheritance of such methylation patterns requires low fidelity methylation. If methylation were maintained with high fidelity at particular CHG sites one would expect that the spontaneous deamination of 5-methylcytosine would diminish the number of such sites, so as to confine the remaining sites to those positions performing an essential function, as is seen in CG methylation.^30–33 However, depletion of CWG sites is not observed in the human genome.³⁴ Since CWG sites account for only about 50% of the non-CG methylation observed in the stem cell methylome¹⁴ where methylated non-CG sites carry only about 25% methylation, the probability of deamination would be about 13% of that for CWG sites that are subject to maintenance methylation in the germ line. Since mutational depletion of methylated cytosines has to have its primary effect on the germ line, if the maintenance of non-CG methylation were more accurate and more widespread, one would have had to argue that stem cells in the human germ lines lack CWG methylation. As it is the data suggests that whatever function non-CG methylation may have in stem cells, it does not involve accurate somatic inheritance in the germ line.The extensive detail on non-CG methylation in the H1 methylome¹⁴ raises interesting questions about the nature of this form of methylation in human cell lines. A key finding in this report is the contrast between the presence of non-CG methylation in the H1 stem cell line and its absence in the IMR90 human fetal lung fibroblast cell line.¹⁴ This suggests that it may have a role in the origin and maintenance of the pluripotent lineage.¹⁴By analogy with the well known methylated DNA binding proteins specific for CG methylation,³⁵ methylated DNA binding proteins that selectively bind sites of non-CG methylation are expected to exist in stem cells. Currently the only protein reported to have this binding specificity is human Dnmt1.^36–38 While Dnmt1 has been proposed to function stoichiometrically³⁹ and could serve a non-CG binding role in stem cells, this possibility and the possibility that other stem-cell specific non-CG binding proteins might exist remain to be been explored.Finally, the nature of the non-CG methylation patterns in human stem cell lines present potentially difficult technical problems in methylation analysis. First, based on the data in the H1 stem cell methylome,⁴⁰ a standard MS-qPCR for non-CG methylation would be impractical because non-CG sites are infrequent, rarely clustered and are generally characterized by partial asymmetric methylation. This means that a PCR primer that senses the 3 adjacent methylation sites usually recommended for MS-qPCR primer design^41,42 cannot be reliably found. For example in the region near Oct4 (Chr6:31,246,431), a potential MS-qPCR site exists with a suboptimal set of two adjacent CHG sites both methylated on the + strand at Chr6:31,252,225 and 31,252,237.^14,40 However these sites were methylated only in 13/45 and 30/52 reads. Thus the probability that they would both be methylated on the same strand is about 17%. Moreover, reverse primer locations containing non-CG methylation sites are generally too far away for practical bisulfite mediated PCR. Considering the losses associated with bisulfite mediated PCR⁴³ the likelihood that such an MS-qPCR system would detect non-CG methylation in the H1 cell line or stem cells present in a cancer stem cell niche^44,45 is very low.The second difficulty is that methods based on the specificity of MeCP2 and similar methylated DNA binding proteins for enriching methylated DNA (e.g., MIRA,⁴⁶ COMPARE-MS⁴⁷) will discard sequences containing non-CG methylation since they require cooperative binding afforded by runs of adjacent methylated CG sites for DNA capture. This latter property of the methylated cytosine capture techniques makes it also unlikely that methods based on 5-methylcytosine antibodies (e.g., meDIP⁴⁸) will capture non-CG methylation patterns accurately since the stem cell methylome shows that adjacent methylated non-CG sites are rare in comparison to methylated CG sites.¹⁴In summary, whether or not mammalian stem cells in general or human stem cells in particular possess functional plant-like methylation patterns is likely to continue to be an interesting and challenging question. At this point we can conclude that the non-CG patterns reported in human cells appear to differ significantly from the non-CG patterns seen in plants, suggesting that they do not have a common origin or function.     

Locus-specific control of DNA methylation by the Arabidopsis SUVH5 histone methyltransferase

In Arabidopsis thaliana, heterochromatin formation is guided by double-stranded RNA (dsRNA), which triggers methylation of histone H3 at Lys-9 (H3 mK9) and CG plus non-CG methylation on identical DNA sequences. At heterochromatin targets including transposons and centromere repeats, H3 mK9 mediated by the Su(var)3-9 homologue 4 (SUVH4)/KYP histone methyltransferase (MTase) is required for the maintenance of non-CG methylation by the CMT3 DNA MTase. Here, we show that although SUVH4 is the major H3 K9 MTase, the SUVH5 protein also has histone MTase activity in vitro and contributes to the maintenance of H3 mK9 and CMT3-mediated non-CG methylation in vivo. Strikingly, the relative contributions of SUVH4, SUVH5, and a third related histone MTase, SUVH6, to non-CG methylation are locus-specific. For example, SUVH4 and SUVH5 together control transposon sequences with only a minor contribution from SUVH6, whereas SUVH4 and SUVH6 together control a transcribed inverted repeat source of dsRNA with only a minor contribution from SUVH5. This locus-specific variation suggests different mechanisms for recruiting or activating SUVH enzymes at different heterochromatic sequences. The suvh4 suvh5 suvh6 triple mutant loses both monomethyl and dimethyl H3 K9 at target loci. The suvh4 suvh5 suvh6 mutant also displays a loss of non-CG methylation similar to a cmt3 mutant, indicating that SUVH4, SUVH5, and SUVH6 together control CMT3 activity.     

Multiplex CRISPR-Cas9 editing of DNA methyltransferases in rice uncovers a class of non-CG methylation specific for GC-rich regions

DNA methylation in the non-CG context is widespread in the plant kingdom and abundant in mammalian tissues such as the brain and pluripotent cells. Non-CG methylation in Arabidopsis thaliana is coordinately regulated by DOMAINS REARRANGED METHYLTRANSFERASE (DRM) and CHROMOMETHYLASE (CMT) proteins but has yet to be systematically studied in major crops due to difficulties in obtaining genetic materials. Here, utilizing the highly efficient multiplex CRISPR-Cas9 genome-editing system, we created single- and multiple-knockout mutants for all the nine DNA methyltransferases in rice (Oryza sativa) and profiled their whole-genome methylation status at single-nucleotide resolution. Surprisingly, the simultaneous loss of DRM2, CHROMOMETHYLASE3 (CMT2), and CMT3 functions, which completely erases all non-CG methylation in Arabidopsis, only partially reduced it in rice. The regions that remained heavily methylated in non-CG contexts in the rice Os-dcc (Osdrm2/cmt2/cmt3a) triple mutant had high GC contents. Furthermore, the residual non-CG methylation in the Os-dcc mutant was eliminated in the Os-ddccc (Osdrm2/drm3/cmt2/cmt3a/cmt3b) quintuple mutant but retained in the Os-ddcc (Osdrm2/drm3/cmt2/cmt3a) quadruple mutant, demonstrating that OsCMT3b maintains non-CG methylation in the absence of other major methyltransferases. Our results showed that OsCMT3b is subfunctionalized to accommodate a distinct cluster of non-CG-methylated sites at highly GC-rich regions in the rice genome.

Examination of knockout mutants reveals that rice methyltransferases have subfunctionalized to accommodate a distinct cluster of non-CG-methylated sites at highly GC-rich regions in the rice genome.     

Base-resolution analyses of sequence and parent-of-origin dependent DNA methylation in the mouse genome

Differential methylation of the two parental genomes in placental mammals is essential for genomic imprinting and embryogenesis. To systematically study this epigenetic process, we have generated a base-resolution, allele-specific DNA methylation (ASM) map in the mouse genome. We find parent-of-origin dependent (imprinted) ASM at 1,952 CG dinucleotides. These imprinted CGs form 55 discrete clusters including virtually all known germline differentially methylated regions (DMRs) and 23 previously unknown DMRs, with some occurring at microRNA genes. We also identify sequence-dependent ASM at 131,765 CGs. Interestingly, methylation at these sites exhibits a strong dependence on the immediate adjacent bases, allowing us to define a conserved sequence preference for the mammalian DNA methylation machinery. Finally, we report a surprising presence of non-CG methylation in the adult mouse brain, with some showing evidence of imprinting. Our results provide a resource for understanding the mechanisms of imprinting and allele-specific gene expression in mammalian cells.     

Regulated AtHKT1 gene expression by a distal enhancer element and DNA methylation in the promoter plays an important role in salt tolerance

Through sos3 (salt overly sensitive 3) suppressor screening, two allelic suppressor mutants that are weak alleles of the strong sos3 suppressor sos3hkt1-1 were recovered. Molecular characterization identified T-DNA insertions in the distal promoter region of the Arabidopsis thaliana HKT1 (AtHKT1, At4g10310) in these two weak sos3 suppressors, which results in physical separation of a tandem repeat from the proximal region of the AtHKT1 promoter. The tandem repeat is approximately 3.9 kb upstream of the ATG start codon and functions as an enhancer element to promote reporter gene expression. A putative small RNA target region about 2.6 kb upstream of the ATG start codon is heavily methylated. CHG and CHH methylation but not CG methylation is significantly reduced in the small RNA biogenesis mutant rdr2, indicating that non-CG methylation in this region is mediated by small RNAs. Analysis of AtHKT1 expression in rdr2 suggests that non-CG methylation in the putative small RNA target region represses AtHKT1 expression in shoots. The DNA methylation-deficient mutant met1-3 has nearly complete loss of total cytosine methylation in the putative small RNA target region and is hypersensitive to salt stress. The putative small RNA target region and the tandem repeat are essential for maintaining AtHKT1 expression patterns crucial for salt tolerance.     

A comparative analysis of DNA methylation across human embryonic stem cell lines

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.     

Involvement of a Jumonji-C domain-containing histone demethylase in DRM2-mediated maintenance of DNA methylation

Histone demethylases-both lysine-specific demethylase 1 (LSD1) and Jumonji-C (JmjC) domain-containing proteins-are broadly implicated in the regulation of chromatin-dependent processes. In Arabidopsis thaliana, histone marks directly affect DNA methylation, and mutations in LSD1 homologues show reduced DNA methylation at some loci. We screened transfer DNA mutations in genes encoding JmjC domains for defects in DNA methylation. Mutations in jmj14 result in reduced DNA methylation in non-CG contexts at targets of DRM2 (domains rearranged methyltransferase 2)-mediated RNA-directed DNA methylation (RdDM), which is associated with an increase in H3K4m3. Unlike other components of RdDM, JMJ14 is not required for de novo methylation of a transgene, suggesting that JMJ14 is specifically involved in the maintenance phase of DRM2-mediated RdDM.     

An Arabidopsis SET domain protein required for maintenance but not establishment of DNA methylation

Cytosine methylation is critical for correct development and genome stability in mammals and plants. In order to elucidate the factors that control genomic DNA methylation patterning, a genetic screen for mutations that disrupt methylation-correlated silencing of the endogenous gene PAI2 was conducted in Arabidopsis: This screen yielded seven loss-of-function alleles in a SET domain protein with histone H3 Lys9 methyltransferase activity, SUVH4. The mutations conferred reduced cytosine methylation on PAI2, especially in non-CG sequence contexts, but did not affect methylation on another PAI locus carrying two genes arranged as an inverted repeat. Moreover, an unmethylated PAI2 gene could be methylated de novo in the suvh4 mutant background. These results suggest that SUVH4 is involved in maintenance but not establishment of methylation at particular genomic regions. In contrast, a heterochromatin protein 1 homolog, LHP1, had no effect on PAI methylation.