- Open Access
Asymmetric DNA methylation of CpG dyads is a feature of secondary DMRs associated with the Dlk1/Gtl2 imprinting cluster in mouse
Epigenetics & Chromatin volume 10, Article number: 31 (2017)
Differential DNA methylation plays a critical role in the regulation of imprinted genes. The differentially methylated state of the imprinting control region is inherited via the gametes at fertilization, and is stably maintained in somatic cells throughout development, influencing the expression of genes across the imprinting cluster. In contrast, DNA methylation patterns are more labile at secondary differentially methylated regions which are established at imprinted loci during post-implantation development. To investigate the nature of these more variably methylated secondary differentially methylated regions, we adopted a hairpin linker bisulfite mutagenesis approach to examine CpG dyad methylation at differentially methylated regions associated with the murine Dlk1/Gtl2 imprinting cluster on both complementary strands.
We observed homomethylation at greater than 90% of the methylated CpG dyads at the IG-DMR, which serves as the imprinting control element. In contrast, homomethylation was only observed at 67–78% of the methylated CpG dyads at the secondary differentially methylated regions; the remaining 22–33% of methylated CpG dyads exhibited hemimethylation.
We propose that this high degree of hemimethylation could explain the variability in DNA methylation patterns at secondary differentially methylated regions associated with imprinted loci. We further suggest that the presence of 5-hydroxymethylation at secondary differentially methylated regions may result in hemimethylation and methylation variability as a result of passive and/or active demethylation mechanisms.
Genomic imprinting in mammals results in the parent of origin-specific monoallelic expression of a subset of genes. Achieving the appropriate balance of gene expression from the maternally and paternally contributed genomes via the establishment of parental allele-specific imprinting marks is crucial for normal growth and development. Therefore, it is critical to understand the mechanisms responsible for controlling the expression of imprinted genes. To date, approximately 150 mammalian genes have been identified as imprinted [1, 2]. Most imprinted genes are found within clusters that contain a CpG-rich imprinting control region (ICR) that functions both to specify parental origin and to regulate imprinted expression of the genes within the cluster [3, 4]. Monoallelic expression of imprinted genes is achieved via multiple mechanisms, including epigenetic modifications such DNA methylation and histone modifications, as well as the activity of long noncoding RNAs [3, 4].
Differential DNA methylation at imprinted loci has been shown to play an important role in distinguishing the parental alleles and regulating their expression [5,6,7,8,9]. Differentially methylated regions (DMRs) associated with imprinted genes fall into two categories: primary and secondary DMRs. Primary, or gametic, DMRs serve as imprinting control regions (ICRs), functioning both to specify parental origin and as a shared regulatory element that controls the expression of genes throughout the associated imprinting cluster. Primary DMRs acquire DNA methylation on one of the two parental alleles during gametogenesis and remain differentially methylated from fertilization throughout development, thereby marking parental origin . The differentially methylated state of primary DMRs can affect expression in a variety of ways. For example, primary DMRs can regulate gene expression through their differential ability to bind enhancer blocking proteins, thereby influencing the activity of an insulator [10, 11]. In other cases, primary DMRs are located at promoters, where they have been shown to directly influence the expression of long noncoding RNAs that subsequently regulate the expression of other genes in the imprinting cluster [12,13,14,15]. In contrast, secondary DMRs acquire parent of origin-specific DNA methylation after implantation [16,17,18,19]. Secondary DMRs are generally located at promoters or within gene bodies, and the acquisition of parental allele-specific DNA methylation at these sequences is dependent on differential methylation of the associated ICR, while the converse is not true [8, 9, 20]. While secondary DMRs do not function as primary imprinting marks, methylation of these regions frequently corresponds with gene silencing and may play a role in maintaining imprinted expression [17, 21,22,23].
The DNA methylation associated with primary DMRs is very stable, with the methylated allele displaying 90–100% methylation at the cytosines located in CpG dinucleotides throughout development [5, 19, 24,25,26,27]. DNA methylation at secondary DMRs is less consistent. For example, methylation at secondary DMRs located at the H19 and Gtl2 promoters average 70 and 78.9%, respectively, as compared to methylation at their respective primary DMRs, which average ~90 and 95.8% [5, 28]. We recently illustrated that the highly variable DNA methylation pattern at the secondary DMR associated with the imprinted Dlk1 gene is asymmetric, with 35% of the methylated CpG dyads displaying hemimethylation . The trend that DNA methylation is more stable at primary DMRs than at secondary DMRs associated with imprinted genes has also been observed at human imprinted loci .
Our current study investigates the nature of secondary DMRs associated with imprinted loci and potential causes of methylation instability, such as a failure to maintain DNA methylation and/or active demethylation catalyzed by the TET enzymes [30,31,32,33,34]. To test the hypothesis that variably methylated secondary DMRs display higher levels of hemimethylation than stably methylated primary DMRs, we analyzed DNA methylation at two additional DMRs associated with the Dlk1/Gtl2 imprinting cluster: the IG-DMR, a primary DMR, and the Gtl2-DMR, a secondary DMR [14, 28]. We also quantified the level of 5-hydroxymethylation (5-hmC) at the IG-, Gtl2- and Dlk1-DMRs to determine whether there is a correlation between high levels of hemimethylation and high levels of 5-hydroxymethylation. Our results support the hypothesis that high levels of 5-hmC may contribute to methylation instability at secondary DMRs associated with imprinted genes.
CpG dyads within the Gtl2-DMR display high levels of hemimethylation
To determine whether asymmetric methylation is unique to the Dlk1-DMR or is a feature of other secondary DMRs associated with imprinted loci, we examined CpG dyad methylation at the linked Gtl2-DMR. We had previously analyzed DNA methylation on the coding strand of the Gtl2-DMR and observed moderate variability in the methylation status, with the 5′ half of the region analyzed displaying lower average DNA methylation levels than the 3′ half . We therefore assessed the DNA methylation status of cytosines located in 22 pairs of complementary CpG dinucleotides spanning this region to determine whether these CpG dyads were homomethylated versus hemimethylated (Fig. 1). All of our experiments were conducted using F1 hybrid tissues collected from crosses between C57BL/6 (B6) and a specially derived strain containing Mus musculus castaneus (CAST)-derived sequences from chromosome 12 on an otherwise C57BL/6 genetic background (CAST12) [18, 28], allowing us to distinguish paternally inherited alleles from maternally inherited alleles based on sequence polymorphisms (detailed in the “Methods”).
We analyzed CpG dyad methylation in DNA derived from four developmental stages: 7.5 d.p.c. embryo, 14.5 d.p.c. embryo, 5 d.p.p. liver and adult liver. The DNA methylation patterns on each parental allele were consistent throughout development and were also similar in tissues obtained from reciprocal crosses (Figs. 2, 3). Across all four developmental stages, cytosines in CpG dinucleotides were methylated 80-93% of the time on paternal alleles and 6–10% of the time on maternal alleles (Table 1). We assessed the significance of these results at each developmental stage using a Mann–Whitney U test and found that the median level of DNA methylation was significantly higher on the paternal alleles as compared to the maternal alleles in all of the tissues examined, with P values ranging from <0.0001 to 0.0147 (Table 2; Additional File 1). Furthermore, P values derived from Mann–Whitney U tests illustrate that median DNA methylation levels did not vary significantly across development on either the paternal or the maternal allele (Additional File 1). Average DNA methylation levels did not vary substantially between the 5′ half versus the 3′ half of the analyzed region. These results confirm that the Gtl2-DMR is differentially methylated throughout development.
Homomethylation was observed at 68–78% of the CpG dyads containing methylated cytosine, while hemimethylation was detected at 22–32% of these CpG dyads (Table 1). The levels of homo- and hemimethylation at the Gtl2-DMR were similar to the overall average of 65% homomethylation and 35% hemimethylation observed at the Dlk1-DMR . When we restricted our analysis to the same four developmental stages assessed in this study, we observed 69–74% homomethylation at methylated CpG dyads within the Dlk1-DMR and 26–31% hemimethylation (Table 1). We tested the homo- and hemimethylation levels at the Gtl2- and Dlk1-DMRs for statistical independence by performing Chi-squared analysis and determined that hemimethylation levels are not significantly different at these loci (P = 0.1318). Therefore, we conclude that hemimethylation levels are similar at two distinct secondary DMRs located within the Dlk1/Gtl2 imprinting cluster.
CpG dyads within the IG-DMR display low levels of hemimethylation
We next assessed hemimethylation levels at the IG-DMR, which serves as the imprinting control region for the Dlk1/Gtl2 imprinting cluster [19, 35]. We analyzed 22 CpG dyads located within the IG-DMR (Fig. 1). We had previously analyzed DNA methylation on the coding strand of this region and had found it to lack variability, with paternally inherited alleles showing near 100% DNA methylation and maternally inherited alleles displaying less than 10% DNA methylation . Consistent with our previous findings, we observed methylation at 96 and 12% of paternally versus maternally inherited CpG dinucleotides located within the IG-DMR, respectively (Fig. 4; Table 1). The median levels of DNA methylation were significantly higher on paternally derived alleles as compared to maternally derived alleles for all tissues analyzed, with P values ranging from <0.0001 to <0.01 (Table 2; Additional File 1), confirming that this region is differentially methylated throughout development. There were no significant differences in the DNA methylation profile of maternal alleles across development (Additional File 1). In contrast, while median DNA methylation levels on paternal alleles was not significantly different between the 14.5 d.p.c. embryo, 5 d.p.p. liver or adult liver samples, the distribution of DNA methylation on paternal alleles derived from 7.5 d.p.c. embryos was different from the distribution in 14.5 d.p.c. embryos (P = 0.0021), 5 d.p.p. liver (P = 0.0178) and adult liver (P = 0.0006). The significance of these results may be attributed to the fact that paternal alleles derived from 7.5 d.p.c. embryos contain more unmethylated cytosines than paternal alleles derived from other developmental stages (Fig. 4; Additional File 1).
Of the CpG dyads displaying cytosine methylation within the IG-DMR, 92% were homomethylated, while 8% were hemimethylated (Table 1). The frequency of hemimethylation was higher on maternally inherited alleles, which were much more sparsely methylated than paternally inherited alleles (Table 1). On average, the level of hemimethylation we observed at the IG-DMR at each developmental stage and across development was lower than the level we observed at either the Dlk1-DMR or the Gtl2-DMR (Fig. 5). We assessed the significance of this result using a Chi-squared test of independence. The level of hemimethylation at the primary IG-DMR is significantly different than the level at either of the secondary DMRs (IG- vs. Dlk1-DMR, P-5.93 × 10−65; IG- vs. Gtl2-DMR, P = 8.76 × 10−57), supporting our hypothesis that high levels of hemimethylation are characteristic of secondary, but not primary, DMRs associated with imprinted loci.
For all of our DNA methylation analyses, we employed a conservative approach whereby we grouped subclones that were derived from the same PCR amplification and had identical sequence and DNA methylation patterns, as it is not possible to determine whether these amplicons are derived from the same or different template molecules. As it is possible that some of the grouped subclones actually represent independent samples, we performed the same analyses for ungrouped data sets in which each subclone was treated as an independent sample. Ungrouping the identical subclones resulted in greater significant differences between the parental alleles at both the Gtl2-DMR and the IG-DMR, confirming their differentially methylated status (Table 2; Additional File 1). While ungrouping identical subclones reduced hemimethylation values at both the Gtl2- and IG-DMRs, from 31.8 to 29.8 and 8.4 to 6.1%, respectively (Additional file 3: Table S2, Additional file 4: Table S3), the difference in hemimethylation levels between the primary IG-DMR and the secondary Gtl2-DMR remains highly significant as assessed using a Chi-squared test of independence (P = 2.81 × 10−59).
High levels of hemimethylation correlate with high levels of 5-hydroxymethylcytosine at the Dlk1-DMR, but not at the Gtl2-DMR
We hypothesized that hemimethylation at the Dlk1- and Gtl2-DMRs could arise as a result of sequential oxidation of 5-methylcytosine (5-mC) by the TET enzymes followed by either passive depletion of methylation following DNA replication or thymine DNA glycosylase-mediated base excision repair, ultimately resulting in demethylation of that residue [30,31,32,33,34]. If the TET enzymes are responsible for demethylation of cytosines leading to high levels of hemimethylation at secondary DMRs associated with imprinted loci, we would expect to see 5-hydroxymethylcytosine (5-hmC), an oxidation intermediate in this pathway, at these loci [30, 31]. We therefore assessed the relative levels of 5-mC and 5-hmC at the Dlk1-, Gtl2- and IG-DMRs in genomic DNA isolated from 9.5 and 14.5 d.p.c. embryos and from 5 d.p.p. and adult liver. To conduct this analysis, we glucosylated genomic DNA derived from each of the four developmental stages listed above, digested glucosylated and unglucosylated samples with MspI, HpaII or no enzyme, amplified the resulting products using qPCR and calculated percent 5-hmC based on the method previously described by Magalhães et al. . We used 9.5 d.p.c. embryos as our earliest developmental time point, rather than 7.5 d.p.c., in order to have sufficient DNA for these analyses.
At the IG-DMR, where we observed low levels of hemimethylation, we detected correspondingly low levels of 5-hmC in DNA derived from 9.5 to 14.5 d.p.c. embryos and from 5 d.p.p. liver: 3–4% of the methylation detected at the MspI site located within the IG-DMR was 5-hmC (Table 3). Higher 5-hmC levels were detected in adult liver, which also had more variation between biological replicates (Additional file 5: Table S4). In contrast, our analysis of two MspI sites located within the Dlk1-DMR detected 5-hmC levels ranging from 7 to 40% in embryos and neonatal liver, with highest levels observed in 14.5 d.p.c. embryos (Table 3). This higher level of 5-hmC at the Dlk1-DMR correlates with the high levels of hemimethylation detected at this locus. We therefore anticipated that we would observe similarly high levels of 5-hmC at the Gtl2-DMR, which also displays high hemimethylation levels. However, the level of 5-hmC at the Gtl2-DMR inversely correlated with hemimethylation levels at this locus. While we detected high levels of hemimethylation at the Gtl2-DMR, we observed very low levels of 5-hmC at three MspI sites located within this locus (Table 3).
5-hmC is absent at the Dlk1-, Gtl2- and IG-DMRs in triple TET knockout ES cells
To validate the assay used to assess 5-hmC levels in the Dlk1-, Gtl2- and IG-DMRs, we conducted an analysis of 5-hmC in DNA derived from wild-type and triple TET knockout embryonic stem (ES) cells . In wild-type ES cells, we detected 5-hmC at 12, 14 and 10% of the methylated cytosines located within MspI sites at the Dlk1-, Gtl2- and IG-DMRs, respectively (Table 4). In ES cells with a triple knockout of the TET genes, 5-hmC was undetectable at all three loci (Table 4). These data validate that the assay employed to detect 5-hmC and confirm that the presence of 5-hmC at these loci is dependent on functional TET enzyme activity.
Differential DNA methylation plays an important role in the regulation of imprinted genes, directly affecting the activity of ICRs as well as directly or indirectly regulating the expression of genes within imprinting clusters [3 4 and references therein]. Stably maintaining DNA methylation at imprinted loci is critical for normal growth and development, and aberrant DNA methylation patterns are associated not only with abnormal expression of imprinted genes, but also with multiple imprinting disorders [38, 39]. Therefore, understanding the normal methylation patterns and how they are altered are important for elucidating the regulation of imprinted genes.
Investigation of DNA methylation at imprinted loci has identified a difference in methylation stability at primary versus secondary DMRs: methylation patterns at primary DMRs appear to be very stable and consistent, while methylation at secondary DMRs is more variable in both mice and humans [5, 18, 19, 24,25,26,27,28,29]. Our previous analyses identified a correlation between variable methylation patterns and high levels of hemimethylation at CpG dyads at the secondary DMR associated with the imprinted Dlk1 gene . To test the hypothesis that hemimethylation is a normal characteristic of secondary DMRs associated with imprinted loci, we utilized a hairpin bisulfite mutagenesis approach to examine CpG dyad methylation at another secondary DMR (the Gtl2-DMR) as well as at a primary DMR that serves as an ICR (the IG-DMR associated with the Dlk1/Gtl2 imprinting cluster). We found high levels of hemimethylation at the Gtl2-DMR, but not at the IG-DMR, supporting our hypothesis (Figs. 2, 3, 4).
It is possible that some of the hemimethylation we observed resulted from errors associated with bisulfite mutagenesis, including failed conversion of unmethylated cytosines and inappropriate conversion of methylated cytosines . Either one of these errors could lead to a hemimethylated CpG dyad. It is unlikely that these errors had an appreciable effect on the hemimethylation levels we observed, as the high molarity, high temperature, short reaction time methodology we employed have been shown to result in low bisulfite conversion error rates [40, 41]. To test this assumption, we directly calculated the failed conversion rate by examining 20,868 non-CpG cytosine locations in our bisulfite-treated samples. We identified 291 cytosines, some of which may have arisen as a result of PCR-induced error rather than failed conversion (Additional file 6: Table S5). The failed conversion rate we observed of 1.39% is similar to the error rate reported by Genereux et al.  for hairpin-linked molecules treated under high molarity, high temperature, 90-min reaction time mutagenesis conditions. We therefore conclude that the predicted and observed error rates are unlikely to account for an appreciable amount of the hemimethylation we observed, particularly at the secondary DMRs, where hemimethylation levels range from 22 to 35%. Furthermore, even if hemimethylation levels vary slightly due to conversion errors, it would not affect our overall conclusion that the difference in hemimethylation levels between primary and secondary DMRs is highly significant.
Based on the observation that secondary DMRs have high levels of hemimethylation, we further proposed that oxidation of 5-methylcytosine by the TET enzymes could result in hemimethylation either by impeding the activity of DNMT1, resulting in replication-dependent passive demethylation, or by further processing of the oxidized products 5-fC and 5-caC and their removal via base excision repair, a mechanism of active demethylation [30,31,32,33,34]. In support of this hypothesis, we found low levels of 5-hmC at the IG-DMR, which had correspondingly low levels of hemimethylation, and high levels of 5-hmC at the Dlk1-DMR, which had correspondingly high levels of hemimethylation. However, at the Gtl2-DMR, where we observed high levels of hemimethylation, we detected low levels of 5-hmC.
The fact that we observed low levels of 5-hmC at the Gtl2-DMR may be a consequence of our experimental approach. The earliest developmental stage at which we assessed 5-hmC was embryonic day 9.5, when methylation of the Gtl2-DMR is already relatively stable as compared to methylation at earlier embryonic stages such as 6.5 and 7.5 d.p.c. . While we detected very little 5-hmC in DNA derived from 9.5 to 14.5 d.p.c. embryos and from neonatal and adult tissues, it is possible that the Gtl2-DMR contains high levels of 5-hmC earlier in development, when DNA methylation patterns at this locus are more labile. Alternatively, the absence of 5-hmC at the Gtl2-DMR may point to the relative importance of maintaining DNA methylation at this locus in order to appropriately silence Gtl2 expression on the paternal allele and achieve appropriate imprinted expression patterns across the Dlk1/Gtl2 imprinting cluster. Indeed, studies have shown that loss of methylation on the paternally inherited Gtl2-DMR correlates with expression of Gtl2 from this allele [42, 43]. In contrast to the substantiated role of differential DNA methylation at the Gtl2-DMR, we and others have questioned whether differential methylation at the Dlk1-DMR, which is located in the 5th exon of Dlk1, plays an important regulatory role [18, 19].
To further investigate the hypothesis that 5-hmC contributes to DNA methylation instability at secondary DMRs associated with imprinted genes, we are currently examining CpG dyad methylation patterns and 5-hmC levels at additional imprinted loci. The work described herein focused on analyses of primary and secondary DMRs associated with the Dlk1/Gtl2 imprinting cluster on mouse chromosome 12. Our current inquiries include an examination of CpG dyad methylation patterns at imprinted loci distributed across multiple genomic locations, including both paternally and maternally methylated DMRs, and our preliminary results suggest that the relationship between 5-hmC and DNA methylation variability at secondary DMRs may be generalizable (Davis Laboratory, unpublished data).
The presence of 5-hmC at secondary DMRs associated with imprinted genes suggests that the TET proteins play a role in modulating DNA methylation stability at these loci. Several studies have shown that the TET proteins play a role in epigenetic reprogramming throughout development [44 and references therein]. For example, the presence of 5-hmC at repressed promoters and in the gene bodies of expressed genes in ES cells is concomitant with the association of TET1 and TET2, respectively, with these sequences [45,46,47,48]. Recent reports have also illustrated the role of TET enzymes in epigenetic reprogramming during primordial germ cell development, including the erasure of imprinting marks [49,50,51], and that a double knockout of TET1 and TET2 disrupts normal methylation and expression patterns at imprinted loci . Liu and colleagues further reported that some imprinted loci, such as H19, appear to be more sensitive to TET activity based on methylation patterns in wild-type vs. TET knockout (TET-KO) ES cells, suggesting that the TET proteins may have different effects at different imprinted loci . Interestingly, the study conducted by Liu et al.  also showed that there were no substantial changes in methylation at the IG-DMR in the TET-KO ES cells as compared to wild-type ES cells, consistent with our observation that there is very little 5-hmC at the IG-DMR in embryos and neonatal liver (Table 3). Indeed, the IG-DMR has a very stable allele-specific DNA methylation pattern, in accordance with its role as an imprinting control element critical for normal growth and development [9, 54]. In contrast, the consistently higher levels of 5-hmC across we observed across development at the Dlk1-DMR suggests that 5-hmC may be a more stable epigenetic mark as this locus, as it is in the developing brain [data herein; 55].
We have yet to determine whether the high levels of hemimethylation at secondary DMRs plays a functional role. Arand et al.  illustrated a correlation between the reduction in DNA methylation levels and concomitant increase in hemimethylation levels during embryogenesis and primordial germ cell development, and suggested that hemimethylation plays a role in impairing maintenance methylation in order to keep methylation levels low. Similarly, Jin et al.  suggested that TET1 may act to reduce DNA methylation levels at hypomethylated CpG islands. Therefore, it would be interesting to look at allele-specific distribution of 5-hmC at imprinted loci to see whether TETs play a role in keeping unmethylated DMRs unmethylated. The assay we used for the 5-hmC studies described herein was not allele-specific, as that would have required a strain-specific polymorphism in close proximity to an MspI site within the regions of analysis. The development of methods to assess 5-hmC in an allele-specific way would allow us to address this question in order to determine the significance of this epigenetic modification.
Secondary DMRs associated with imprinted loci have low DNA methylation fidelity as compared to primary DMRs that serve as imprinting control regions. Our current analyses illustrate that the variable DNA methylation patterns at secondary DMRs correlate with high levels of hemimethylation at CpG dyads and high levels of 5-hmC. Our work therefore supports the hypothesis that secondary DMRs have a unique epigenetic profile that distinguishes them from primary DMRs. Our studies further provide insight into the molecular mechanisms responsible for methylation instability at secondary DMRs, as oxidation of 5-mC to 5-hmC by the TET enzymes ultimately leads to a loss of DNA methylation passively due to reduced DNMT1 fidelity and/or actively via further oxidation followed by DNA repair. These results are significant as they highlight the complexities associated with the maintenance of the epigenetic profile at secondary DMRs: Differential DNA methylation is maintained at these loci despite activities that function to reduce methylation levels. Further investigation is warranted to understand how parent of origin-specific DNA methylation is established and maintained at secondary DMRs.
C57BL/6J (B6) and Mus musculus castaneus (CAST) mice were purchased from the Jackson Laboratory. To facilitate the isolation of F1 hybrid mice, a strain of mice that served as the source of the M. m. castaneus allele (CAST12) was constructed as previously described . Natural matings between B6 and CAST were used to generate F1 hybrid males for spermatozoa collection; all other F1 hybrid tissues used for bisulfite analyses were generated from natural matings between B6 and CAST12 mice. For all F1 hybrid tissues, the maternal allele is located on the left.
DNA purification and bisulfite analysis
DNA was isolated from 7.5 d.p.c. embryos using the DNeasy Blood & Tissue Kit (Qiagen, Germantown, MD, cat#69504). DNA was isolated from 9.5 to 14.5 d.p.c. embryos and from 5 d.p.p. and adult liver following proteinase K digestion and a series of phenol/chloroform extractions as described previously . Prior to bisulfite mutagenesis, complementary strands of DNA were covalently attached as follows: for IG-DMR analyses, 0.5 µg of genomic DNA was digested with 1 µl SpeI (NEB, Ipswich, MA, cat#R0133) and ligated to 1 µg of phosphorylated hairpin linker IG-DMR-HP (5′-CTAGAGCGATGCGTTCGAGCATCGCT-3′) ; for Gtl2-DMR analyses, 0.5 µg of genomic DNA was digested with 1 µl BanI (NEB, Ipswich, MA, cat#R0118) and ligated to 1 µg of phosphorylated hairpin linker Gtl2-HP-3 (5′-GTACAGCGATGCGTTCGAGCATCGCT-3′). 0.5 µg of hairpin-linked, ligated DNA was denatured by incubating in freshly prepared 3 M NaOH for 20 min at 42 °C, and then subjected to bisulfite mutagenesis using an EZ DNA Methylation-Direct kit (Zymo Research, Irvine, CA, cat#D5020). All mutagenized DNAs were subjected to multiple independent PCR amplifications to ensure analysis of different strands of DNA; subclones derived from independent PCR amplifications are distinguished by different letters of the alphabet. PCR contamination was ruled out via analysis of no template negative control amplifications for both the first and second rounds of PCR. Data from multiple independent tissue samples derived from the same developmental stage were combined, as we did not detect variation between biological replicates. Primer pairs used for nested amplification of the mutagenized DNA were designed to incorporate at least one SNP as well as CpG dinucleotides within the previously analyzed DMRs . All base pair numbers are from GenBank Accession Number NC_000078.6. Primers and PCR cycling conditions for the IG-DMR and for two adjacent regions within the Gtl2-DMR are detailed in Additional file 7: Table S6. Expected second round PCR products for the IG-DMR and the two regions of the Gtl2-DMR are 412 bp, 721 bp, and 695 bp, respectively. Subcloning of amplified products was achieved using a pGEM-T Easy vector (Promega Corporation, Madison, WI, cat#A1360). Sequencing reactions were outsourced to Genewiz (South Plainfield, NJ) or were performed using a Thermo Sequenase Cycle Sequencing Kit (Affymetrix, Cleveland, OH, cat#78500) and analyzed on a 4300 DNA Analyzer (LI-COR Biosciences, Lincoln, NE). Sequence polymorphisms used to distinguish B6 vs. CAST DNA (B/C): IG-DMR, G/A at bp#109,528,369; Gtl2-DMR, G/A at bp#109,541,531, A/G at bp#109,541,671, AA/GC at bp#109,541,820-109,541,821. Percent methylation was calculated based on data obtained from both complementary strands. Percent hemimethylation was calculated by dividing the number of hemimethylated CpG dyads by the number of hemimethylated plus homomethylated CpG dyads.
For 5-hydroxymethylation analyses, DNA was isolated from 9.5 d.p.c. embryos, 14.5 d.p.c. embryos, 5 d.p.p. liver and adult liver as described above. DNA derived from three different genetic backgrounds [C57BL/6 J, B6x(CAST or CAST12) and (CAST or CAST12)xB] was used to ensure that genetic background did not affect the outcome. In addition, DNA was isolated from wild-type and TET triple knockout embryonic stem cells . 5-hydroxymethylation levels were assessed using an EpiMark 5-hmC and 5-mC Analysis Kit (NEB, Ipswich, MA, cat#E3317). Briefly, 2.5 µg genomic DNA or 2 µg ES cell DNA was glucosylated using 30 units of T4 ß-glucosyltransferase at 37 °C overnight. Glucosylated and unglucosylated control DNAs were treated with MspI, HpaII or no restriction endonuclease at 37 °C overnight. Following treatment with proteinase K, products were amplified via PCR and quantitative PCR (StepOnePlus, Applied Biosystems). Primers and PCR cycling conditions used are detailed in Additional file 8: Table S7. qPCR was performed in triplicate for each of three independent biological samples. Amount of 5-mC and 5-hmC in each sample was calculated according to Magalhães et al. . Genomic coordinates for MspI/HpaII sites are: Dlk1-DMR-A, bp#109,459,830; Dlk1-DMR-B, bp#109,460,017; IG-DMR, bp#109,528,624; Gtl2-DMR-A, bp#109,541,314; Gtl2-DMR-B, bp#109,541,776; Gtl2-DMR-C, bp#109,541,811.
imprinting control region
differentially methylated region
- C or CAST:
Mus musculus castaneus
- C12 or CAST12:
Mus musculus castaneus chromosome 12 on a C57BL/6 background
polymerase chain reaction
Morison IM, Ramsay JP, Spencer HG. A census of mammalian imprinting. Trends Genet. 2005;21:457–65.
Williamson CM, Blake A, Thomas S, Beechey CV, Hancock J, Cattanach BM, Peters J. World wide web site, mouse imprinting data and references. Oxfordshire: MRC Hartwell; 2013.
Barlow D, Bartolomei MS. Genomic imprinting in mammals. Cold Spring Harbor Perspect. Biol. 2014. doi:10.1101/cshperspect.a018382.
Adalsteinsson BT, Ferguson-Smith AC. Epigenetic control of the genome: lessons from genomic imprinting. Genes. 2014;5:635–55. doi:10.3390/genes5030635.
Tremblay KD, Duran KL, Bartolomei MS. A 5′ 2-kilobase-pair region of the imprinted mouse H19 gene exhibits exclusive paternal methylation throughout development. Mol Cell Biol. 1997;17:4322–9.
Shemer R, Birger Y, Riggs AD, Razin A. Structure of the imprinted mouse Snrpn gene and establishment of its parental-specific methylation pattern. Proc Natl Acad Sci USA. 1997;94:10267–72.
Yoon B, Herman H, Hu B, Park YJ, Lindroth A, Bell A, West AG, Chang Y, Stablewski A, Piel JC, Loukinov DI, Lobanenkov VV, Soloway PD. Rasgrf1 imprinting is regulated by a CTCF-dependent methylation-sensitive enhancer blocker. Mol Cell Biol. 2005;25:11184–90.
Lopes S, Lewis A, Hajkova P, Dean W, Oswald J, Forné T, Murrell A, Constância M, Bartolomei M, Walter J, Reik W. Epigenetic modifications in an imprinting cluster are controlled by a hierarchy of DMRs suggesting long-range chromatin interactions. Hum Mol Genet. 2003;12:295–305.
Kagami M, O’Sullivan MJ, Green AJ, Watabe Y, Arisaka O, Masawa N, Matsuoka K, Fukami M, Matsubara K, Kato F, Ferguson-Smith AC, Ogata T. The IG-DMR and the MEG3-DMR at human chromosome 14q32.2: hierarchical interaction and distinct functional properties as imprinting control centers. PLoS Genet. 2010;6:1–13 (PMID: 20585555).
Bell AC, Felsenfeld G. Methylation of a CTCF-dependent boundary controls imprinted expression of the Igf2 gene. Nature. 2000;405:482–5.
Hark AT, Schoenherr CJ, Katz DJ, Ingram RS, Levorse JM, Tilghman SM. CTCF mediates methylation-sensitive enhancer-blocking activity at the H19/Igf2 locus. Nature. 2000;405:486–9.
Sleutels F, Zwart R, Barlow DP. The non-coding Air RNA is required for silencing autosomal imprinted genes. Nature. 2002;415:810–3.
Nagano T, Mitchell JA, Sanz LA, Pauler FM, Ferguson-Smith AC, Feil R, Fraser P. The Air noncoding RNA epigenetically silences transcription by targeting G9a to chromatin. Science. 2008;322:1717–20. doi:10.1126/science.1163802.
Latos PA, Pauler FM, Koerner MV, Şenergin HB, Hudson QJ, Stocsits RR, Allhoff W, Stricker SH, Klement RM, Warczok KE, Aumayr K, Pasierbek P, Barlow DP. Airn transcriptional overlap, but not its lncRNA products, induces imprinted Igf2r silencing. Science. 2012;338:1469–72. doi:10.1126/science.1228110.
Mancini-DiNardo D, Steele SJS, Levorse JM, Ingram RS, Tilghman SM. Elongation of the Kcnq1ot1 transcript is required for genomic imprinting of neighboring genes. Genes Dev. 2006;20:1268–82. doi:10.1101/gad.1416906.
Hanel ML, Wevrick R. Establishment and maintenance of DNA methylation patterns in mouse Ndn: implications for maintenance of imprinting in target genes of the imprinting center. Mol Cell Biol. 2001;21:2384–92.
Bhogal B, Arnaudo A, Dymkowski A, Best A, Davis TL. Methylation at mouse Cdkn1c is acquired during postimplantation development and functions to maintain imprinted expression. Genomics. 2004;84:961–70.
Gagne A, Hochman A, Qureshi M, Tong C, Arbon J, McDaniel K, Davis TL. Analysis of DNA methylation acquisition at the imprinted Dlk1 locus reveals asymmetry at CpG dyads. Epigenet Chromatin. 2014;7:9.
Takada S, Paulsen M, Tevendale M, Tsai C-E, Kelsey G, Cattanach BM, Ferguson-Smith AC. Epigenetic analysis of the Dlk1-Gtl2 imprinted domain on mouse chromosome 12: implications for imprinting control from comparison with Igf2-H19. Hum Mol Genet. 2002;11:77–86.
Beygo J, Elbracht M, de Groot K, Begemann M, Kanber D, Platzer K, Gillessen-Kaesbach G, Vierzig A, Green A, Heller R, Buiting K, Eggermann T. Novel deletions affecting the MEG2-DMR provide further evidence for a hierarchical regulation of imprinting in 14q32. Eur J Hum Genet. 2015;23:180–8. doi:10.1038/ejhg.2014.72.
John RM, Lefebvre L. Developmental regulation of somatic imprints. Differentiation. 2011;81:270–80.
Stöger R, Kubicka P, Liu C-G, Kafri T, Razin A, Cedar H, Barlow DP. Maternal-specific methylation of the imprinted mouse Igf2r locus identifies the expressed locus as carrying the imprinting signal. Cell. 1993;73:61–71.
Kobayashi H, Sakurai T, Sato S, Nakabayashi K, Hata K, Kono T. Imprinted DNA methylation reprogramming during early mouse embryogenesis at the Gpr1-Zdbf2 locus is linked to long cis-intergenic transcription. FEBS Lett. 2012;586:827–33.
Yatsuki H, Joh K, Higashimoto K, Soejima H, Arai Y, Wang Y, Hatada I, Obata Y, Morisaki H, Zhang Z, Nakagawachi T, Satoh Y, Mukai T. Domain regulation of imprinting cluster in Kip2/Lit1 subdomain on mouse chromosome 7F4/F5: large-scale DNA methylation analysis reveals that DMR-Lit1 is a putative imprinting control region. Genome Res. 2002;12:1860–70.
Arnaud P, Monk D, Hichins M, Gordon E, Dean W, Beechey CV, Peters J, Craigen W, Preece M, Stanier P, Moore GE, Kelsey G. Conserved methylation imprints in the human and mouse GRB10 genes with divergent allelic expression suggests differential reading of the same mark. Hum Mol Genet. 2003;12:1005–19.
Ono R, Shiura H, Aburatani H, Kohda T, Kaneko-Ishino T, Ishino F. Identification of a large novel imprinted gene cluster on mouse proximal chromosome 6. Genome Res. 2003;13:1696–705.
Coombes C, Arnaud P, Gordon E, Dean W, Coar EA, Williamson CM, Feil R, Peters J, Kelsey G. Epigenetic properties and identification of an imprint mark in the Nesp-Gnasx1 domain of the mouse Gnas imprinted locus. Mol Cell Biol. 2003;23:5475–88.
Nowak K, Stein G, Powell E, He LM, Naik S, Morris J, Marlow S, Davis TL. Establishment of paternal allele-specific DNA methylation at the imprinted mouse Gtl2 locus. Epigenetics. 2011;6:1012–20.
Woodfine K, Huddleston JE, Murrell A. Quantitative analysis of DNA methylation at all human imprinted regions reveals preservation of epigenetic stability in adult somatic tissue. Epigenet Chromatin. 2011;4:1. doi:10.1186/1756-8935-4-1.
Tahiliani M, Koh KP, Shen Y, Pastor WA, Bandukwala H, Brudno Y, Agarwal S, Iyer LM, Liu DR, Aravind L, Rao A. Conversion of 5-methylcytosine to 5-hydroxymethylcytosine in mammalian DNA by MLL partner TET1. Science. 2009;324:930–5.
Ito S, Shen L, Dai Q, Wu SC, Collins LB, Swenberg JA, He C, Zhang Y. Tet proteins can convert 5-methylcytosine to 5-formylcytosine and 5-carboxylcytosine. Science. 2011;333:1300–3.
He YF, Li BZ, Li Z, Liu P, Wang Y, Tang Q, Ding J, Jia Y, Chen Z, Li L, Sun Y, Li X, Dai Q, Song CX, Zhang K, He C, Xu GL. Tet-mediated formation of 5-carboxylcytosine and its excision by TDG in mammalian DNA. Science. 2011;333:1303–7.
Kohli RM, Zhang Y. TET enzymes, TDG and the dynamics of DNA demethylation. Nature. 2013;502:472–9. doi:10.1038/nature12750.
Valinluck V, Sowers LC. Endogenous cytosine damage products alter the site selectivity of human DNA maintenance methyltransferase DNMT1. Cancer Res. 2007;67:946–50.
Lin S-P, Coan P, Teixeira da Rocha S, Seitz H, Cavaille J, Teng P-W, Takada S, Ferguson-Smith AC. Differential regulation of imprinting in the murine embryo and placenta by the Dlk1-Dio3 imprinting control region. Development. 2007;134:417–26.
Magalhães HR, Leite SBP, de Paz CCP, Duarte G, Ramos ES. Placental hydroxymethylation vs methylation at the imprinting control region 2 on chromosome 11p15.5. Braz J Med Biol Res. 2013;46:916–9. doi:10.1590/1414-431X20133035.
Hu X, Zhang L, Mao S-Q, Li Z, Chen J, Zhang R-R, Wu H-P, Gao J, Guo F, Liu W, Xu G-F, Dai H-Q, Shi YG, Li X, Hu B, Tang F, Pei D, Xu G-L. Tet and TDG mediate DNA demethylation essential for mesenchymal-to-epithelial transition in somatic cell reprogramming. Cell Stem Cell. 2014;14:512–22. doi:10.1016/j.stem.2014.01.001.
Mackay DJG, Callaway JLA, Marks SM, White HE, Acerini CL, Boonen SE, Dayanikli P, Firth HV, Goodship JA, Haemers AP, Hahnemann JMD, Kordonouri O, Masoud AF, Oestergaard E, Storr J, Ellard S, Hattersley AT, Robinson DO, Temple IK. Hypomethylation of multiple imprinted loci in individuals with transient neonatal diabetes is associated with mutations in ZFP57. Nat Genet. 2008;40:949–51. doi:10.1038/ng.187.
Eggermann T, Perez de Nanciares G, Maher ER, Temple IK, Tümer Z, Monk D, Mackay DJ, Grøznskov K, Ricci A, Linglart A, Netchine I. Imprinting disorders: a group of congenital disorders with overlapping patterns of molecular changes affecting imprinted loci. Clin Epigenet. 2015;7:123. doi:10.1186/s13148-015-0143-8.
Genereux DP, Johnson WC, Burden AF, Stöger R, Laird CD. Errors in the bisulfite conversion of DNA: modulating inappropriate- and failed-conversion frequencies. Nucleic Acids Res. 2008;36(22):e150. doi:10.1093/nar/gkn691.
Holmes EE, Jung M, Meller S, Leisse A, Sailer V, Zech J, Mengdehl M, Garbe L-A, Uhl B, Kristiansen G, Dietrich D. Performance evaluation of kits for bisulfite-conversion of DNA from tissues, cell lines, FFPE tissues, aspirates, lavages, effusions, plasma, serum and urine. PLoS ONE. 2014;9(4):e93933. doi:10.1371/journal.pone.0093933.
Steshina EY, Carr MS, Glick EA, Yevtodiyenk A, Appelbe OK, Schmidt JV. Loss of imprinting at the Dlk1-Gtl2 locus caused by insertional mutagenesis in the Gtl2 5′ region. BMC Genet. 2006;7:44.
Sekita Y, Wagatsuma H, Irie M, Kobayashi S, Kohda T, Matsuda J, Yokoyama M, Ogura A, Schuster-Gossler K, Gossler A, Ishino F, Kaneko-Ishino T. Aberrant regulation of imprinted gene expression in Gtl2lacZ mice. Cytogenet Genome Res. 2006;113:223–9.
Li D, Guo B, Wu H, Tan L, Lu Q. TET family of dioxygenases: crucial roles and underlying mechanisms. Cytogenet Genome Res. 2015;146:171–80. doi:10.1159/000438853.
Pastor WA, Pape UJ, Huang Y, Henderson HR, Liste R, Ko M, McLoughlin EM, Brudno Y, Mahapatra S, Kapranov P, Tahiliani M, Daley GQ, Liu XS, Ecker JR, Milos PM, Agarwal S, Rao A. Genome-wide mapping of 5-hydroxymethylcytosine in embryonic stem cells. Nature. 2011;473:394–7. doi:10.1038/nature10102.
Wu H, D’Alessilo AC, Ito S, Xia K, Wang Z, Cui K, Zhao K, Sun E, Zhang Y. Dual functions of Tet1 in transcriptional regulation in mouse embryonic stem cells. Nature. 2011;473:389–93. doi:10.1038/nature09934.
Stroud H, Feng S, Morey Kinney S, Pradhan S, Jacobsent SE. 5-hydroxymethylcytosine is associated with enhancers and gene bodies in human embryonic stem cells. Genome Biol. 2011;12(6):R54. doi:10.1186/gb-2011-12-6-r54.
Huang Y, Chavez L, Chang X, Wang X, Pastor WA, Kang J, Zepeda-Martinez JA, Pape UJ, Jacobsen SE, Peters B, Rao A. Distinct roles of the methylcytosine oxidases Tet1 and Tet2 in mouse embryonic stem cells. Proc Natl Acad Sci USA. 2014;111:1361–6. doi:10.1073/pnas.1322921111.
Yamaguchi S, Hong K, Liu R, Inoue A, Shen L, Zhang K, Zhang Y. Dynamics of 5-methylcytosine and 5-hydroxymethylcytosine during germ cell reprogramming. Cell Res. 2013;23:329–39.
Hackett JA, Sengupta R, Zylicz JJ, Murakami K, Lee C, Down TA, Surani MA. Germline DNA demethylation dynamics and imprint erasure through 5-hydroxymethylcytosine. Science. 2013;339:448–52.
Yamaguchi S, Shen L, Liu Y, Sendler D, Zhang Y. Role of Tet1 in erasure of genomic imprinting. Nature. 2013;504:460–4. doi:10.1038/nature12805.
Dawlaty MM, Breiling A, Le T, Raddatz G, Barrasa MI, Cheng AW, Gao Q, Powell BE, Li Z, Xu M, Faull KF, Lyko F, Jaenisch R. Combined deficiency of Tet1 and Tet2 causes epigenetic abnormalities but is compatible with postnatal development. Dev Cell. 2013;24:310–23. doi:10.1016/j.devcell.2012.12.015.
Liu L, Mao S-Q, Ray C, Zhang Y, Bell FT, Ng S-F, Xu G-L, Li X. Differential regulation of genomic imprinting by TET proteins in embryonic stem cells. Stem Cell Res. 2015;15:435–43. doi:10.1016/j.scr.2015.08.010.
Georgiades P, Watkins M, Surani MA, Ferguson-Smith AC. Parental origin-specific developmental defects in mice with uniparental disomy for chromosome 12. Development. 2000;127:4719–28.
Hahn MA, Qiu R, Wu X, Li AX, Zhang H, Wang J, Jui J, Jin SG, Jiang Y, Pfeifer GP, Lu Q. Dynamics of 5-hydroxymethylcytosine and chromatin marks in mammalian neurogenesis. Cell Rep. 2013;3:291–300.
Arand J, Wossidlo M, Lepikhow K, Peat JR, Reik W, Walte J. Selective impairment of methylation maintenance is the major cause of DNA methylation reprogramming in the early embryo. Epigenet Chromatin. 2015;8:1.
Jin C, Lu Y, Jelinek J, Liang S, Estecio MRH, Barton MC. Issa J-PJ. TET1 is a maintenance DNA demethylase that prevents methylation spreading in differentiated cells. Nucleic Acids Res. 2014;42:6956–71. doi:10.1093/nar/gku372.
Davis TL, Trasler JM, Moss SB, Yang GJ, Bartolomei MS. Acquisition of the H19 methylation imprint occurs differentially on the parental alleles during spermatogenesis. Genomics. 1999;58:18–28.
Laird CD, Pleasant ND, Clark AD, Sneeden JL, Hassan KMA, Manley NC, Vary JC, Morgan T, Hansen RS, Stöger R. Hairpin-bisulfite PCR: assessing epigenetic methylation patterns on complementary strands of individual DNA molecules. Proc Natl Acad Sci USA. 2004;101:204–9.
MG and EV carried out molecular genetic studies. TLD conceived of the study and experimental design, carried out molecular genetic studies and wrote the manuscript. All authors read and approved the final manuscript.
We thank Marisa Bartolomei for providing DNA derived from wild-type and TET triple knockout embryonic stem cells. We thank Jessica Arbon and Chris Pathmanabhan for their contributions toward this work, Joshua Shapiro for advice on statistical analyses and Marisa Bartolomei for careful reading of the manuscript and thoughtful discussion.
The authors declare that they have no competing interests.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its supplementary information files.
Ethical approval for procedures involving animals was granted by the Bryn Mawr College Institutional Animal Care and Use Committee, PHS Welfare Assurance Number A3920-01.
This work was supported by awards from the Bryn Mawr College Faculty Research Fund and National Science Foundation Grant 1514600 to TLD. In addition, MG and KV were supported in part by the Bryn Mawr College Summer Science Research program.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1. Data used for statistical analyses of DNA methylation levels at the Gtl2- and IG-DMRs. The numerical data used to perform Mann-Whitney U tests and the resulting P values are contained in this file. Data from the Gtl2-DMR 5’ region, Gtl2-DMR 3’ region and the IG-DMR are presented in separate sheets. Within each sheet, data from each of the developmental stages are presented in chronological order, as they are in the Results, Figures, and Tables. Each sheet presents the information for a specific locus, tissue, cross (maternal allele x paternal allele), and parental allele analyzed, as indicated in columns A-D. % methylation (column F) was calculated by dividing the number of methylated CpG sites observed in a given subclone (column E) by the total number of CpG sites analyzed within the subclone; the raw data used to make these calculations are found in Figs. 2, 3 and 4. For the Gtl2-DMR 3’ region and the IG-DMR, P values were calculated independently for BxC samples vs. CxB samples. In addition, P values were calculated for the combined BxC + CxB samples, as some of the BxC and CxB sample sizes were too small to accurately perform Mann-Whitney U tests. Data presented in sheets labelled “grouped” combine subclones derived from the same PCR with identical DNA methylation patterns and identical sequences as a single sample, while every subclone analyzed is treated as an independent sample in sheets labelled “ungrouped”.
Additional file 2: Table S1. Percent methylation on paternal and maternal alleles at DMRs across four developmental stages.
Additional file 3: Table S2. Comparison of percent methylation at the Gtl2-DMR with like subclones grouped vs. ungrouped.
Additional file 4: Table S3. Comparison of percent methylation at the IG-DMR with like subclones grouped vs. ungrouped.
Additional file 5: Table S4. Percent 5-hydroxymethylcytosine at the Dlk1-, IG- and Gtl2-DMRs at four developmental stages.
Additional file 7: Table S6. Primer and PCR cycling conditions for amplification of bisulfite-mutagenized DNA.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
About this article
Cite this article
Guntrum, M., Vlasova, E. & Davis, T.L. Asymmetric DNA methylation of CpG dyads is a feature of secondary DMRs associated with the Dlk1/Gtl2 imprinting cluster in mouse. Epigenetics & Chromatin 10, 31 (2017). https://doi.org/10.1186/s13072-017-0138-0
- Genomic imprinting
- DNA methylation
- Gametic DMR
- Somatic DMR