WBSA: Web Service for Bisulfite Sequencing Data Analysis

Whole-Genome Bisulfite Sequencing (WGBS) and genome-wide Reduced Representation Bisulfite Sequencing (RRBS) are widely used to study DNA methylation. However, data analysis is complicated, lengthy, and hampered by a lack of seamless analytical pipelines. To address these issues, we developed a convenient, stable, and efficient web service called Web Service for Bisulfite Sequencing Data Analysis (WBSA) to analyze bisulfate sequencing data. WBSA focuses on not only CpG methylation, which is the most common biochemical modification in eukaryotic DNA, but also non-CG methylation, which have been observed in plants, iPS cells, oocytes, neurons and stem cells of human. WBSA comprises three main modules as follows: WGBS data analysis, RRBS data analysis, and differentially methylated region (DMR) identification. The WGBS and RRBS modules execute read mapping, methylation site identification, annotation, and advanced analysis, whereas the DMR module identifies actual DMRs and annotates their correlations to genes. WBSA can be accessed and used without charge either online or local version. WBSA also includes the executables of the Portable Batch System (PBS) and standalone versions that can be downloaded from the website together with the installation instructions. WBSA is available at no charge for academic users at http://wbsa.big.ac.cn.     

Detection of differentially methylated regions from whole-genome bisulfite sequencing data without replicates

DNA methylation is an important epigenetic modification involved in many biological processes and diseases. Recent developments in whole genome bisulfite sequencing (WGBS) technology have enabled genome-wide measurements of DNA methylation at single base pair resolution. Many experiments have been conducted to compare DNA methylation profiles under different biological contexts, with the goal of identifying differentially methylated regions (DMRs). Due to the high cost of WGBS experiments, many studies are still conducted without biological replicates. Methods and tools available for analyzing such data are very limited.We develop a statistical method, DSS-single, for detecting DMRs from WGBS data without replicates. We characterize the count data using a rigorous model that accounts for the spatial correlation of methylation levels, sequence depth and biological variation. We demonstrate that using information from neighboring CG sites, biological variation can be estimated accurately even without replicates. DMR detection is then carried out via a Wald test procedure. Simulations demonstrate that DSS-single has greater sensitivity and accuracy than existing methods, and an analysis of H1 versus IMR90 cell lines suggests that it also yields the most biologically meaningful results. DSS-single is implemented in the Bioconductor package DSS.     

Calling differentially methylated regions from whole genome bisulphite sequencing with DMRcate

Whole genome bisulphite sequencing (WGBS) permits the genome-wide study of single molecule methylation patterns. One of the key goals of mammalian cell-type identity studies, in both normal differentiation and disease, is to locate differential methylation patterns across the genome. We discuss the most desirable characteristics for DML (differentially methylated locus) and DMR (differentially methylated region) detection tools in a genome-wide context and choose a set of statistical methods that fully or partially satisfy these considerations to compare for benchmarking. Our data simulation strategy is both biologically informed—employing distribution parameters derived from large-scale consortium datasets—and thorough. We report DML detection ability with respect to coverage, group methylation difference, sample size, variability and covariate size, both marginally and jointly, and exhaustively with respect to parameter combination. We also benchmark these methods on FDR control and computational time. We use this result to backend and introduce an expanded version of DMRcate: an existing DMR detection tool for microarray data that we have extended to now call DMRs from WGBS data. We compare DMRcate to a set of alternative DMR callers using a similarly realistic simulation strategy. We find DMRcate and RADmeth are the best predictors of DMRs, and conclusively find DMRcate the fastest.     

DMRFusion: A differentially methylated region detection tool based on the ranked fusion method

DNA methylation is an important epigenetic modification involved in many biological processes and diseases. Computational analysis of differentially methylated regions (DMRs) could explore the underlying reasons of methylation. DMRFusion is presented as a useful tool for comprehensive DNA methylation analysis of DMRs on methylation sequencing data. This tool is designed base on the integration of several ranking methods; Information gain, Between versus within Class scatter ratio, Fisher ratio, Z-score and Welch's t-test. In this study, DMRFusion on reduced representation bisulfite sequencing (RRBS) data in chronic lymphocytic leukemia cancer displayed 30 nominated regions and CpG sites with a maximum methylation difference detected in the hypermethylation DMRs. We realized that DMRFusion is able to process methylation sequencing data in an efficient and accurate manner and to provide annotation and visualization for DMRs with high fold difference score (p-value and FDR < 0.05 and type I error: 0.04).     

Bisulfighter: accurate detection of methylated cytosines and differentially methylated regions

Analysis of bisulfite sequencing data usually requires two tasks: to call methylated cytosines (mCs) in a sample, and to detect differentially methylated regions (DMRs) between paired samples. Although numerous tools have been proposed for mC calling, methods for DMR detection have been largely limited. Here, we present Bisulfighter, a new software package for detecting mCs and DMRs from bisulfite sequencing data. Bisulfighter combines the LAST alignment tool for mC calling, and a novel framework for DMR detection based on hidden Markov models (HMMs). Unlike previous attempts that depend on empirical parameters, Bisulfighter can use the expectation-maximization algorithm for HMMs to adjust parameters for each data set. We conduct extensive experiments in which accuracy of mC calling and DMR detection is evaluated on simulated data with various mC contexts, read qualities, sequencing depths and DMR lengths, as well as on real data from a wide range of biological processes. We demonstrate that Bisulfighter consistently achieves better accuracy than other published tools, providing greater sensitivity for mCs with fewer false positives, more precise estimates of mC levels, more exact locations of DMRs and better agreement of DMRs with gene expression and DNase I hypersensitivity. The source code is available at http://epigenome.cbrc.jp/bisulfighter.     

Bisulfite sequencing and dinucleotide content analysis of 15 imprinted mouse differentially methylated regions (DMRs): paternally methylated DMRs contain less CpGs than maternally methylated DMRs

Imprinted genes in mammals show monoallelic expression dependent on parental origin and are often associated with differentially methylated regions (DMRs). There are two classes of DMR: primary DMRs acquire gamete-specific methylation in either spermatogenesis or oogenesis and maintain the allelic methylation differences throughout development; secondary DMRs establish differential methylation patterns after fertilization. Targeted disruption of some primary DMRs showed that they dictate the allelic expression of nearby imprinted genes and the establishment of the allelic methylation of secondary DMRs. However, how primary DMRs are recognized by the imprinting machinery is unknown. As a step toward elucidating the sequence features of the primary DMRs, we have determined the extents and boundaries of 15 primary mouse DMRs (including 12 maternally methylated and three paternally methylated DMRs) in 12.5-dpc embryos by bisulfite sequencing. We found that the average size of the DMRs was 3.2 kb and that their average G+C content was 54%. Dinucleotide content analysis of the DMR sequences revealed that, although they are generally CpG rich, the paternally methylated DMRs contain less CpGs than the maternally methylated DMRs. Our findings provide a basis for the further characterization of DMRs.     

MIRA-SNuPE,a quantitative,multiplex method for measuring allele-specific DNA methylation

5-methyl-C (5mC) and 5-hydroxymethyl-C (5hmC) are epigenetic marks with well-known and putative roles in gene regulation, respectively. These two DNA covalent modifications cannot be distinguished by bisulfite sequencing or restriction digestion, the standard methods of 5mC detection. The methylated CpG island recovery assay (MIRA), however, specifically detects 5mC but not 5hmC. We further developed MIRA for the analysis of allele-specific CpG methylation at differentially methylated regions (DMRs) of imprinted genes. MIRA specifically distinguished between the parental alleles by capturing the paternally methylated H19/Igf2 DMR and maternally methylated KvDMR1 in mouse embryo fibroblasts (MEFs) carrying paternal and maternal duplication of mouse distal Chr7, respectively. MIRA in combination with multiplex single nucleotide primer extension (SNuPE) assays specifically captured the methylated parental allele from normal cells at a set of maternally and paternally methylated DMRs. The assay correctly recognized aberrant biallelic methylation in a case of loss of imprinting. The MIRA-SNuPE assays revealed that placenta exhibited less DNA methylation bias at DMRs compared to yolk sac, amnion, brain, heart, kidney, liver and muscle. This method should be useful for the analysis of allele-specific methylation events related to genomic imprinting, X chromosome inactivation and for verifying and screening haplotype-associated methylation differences in the human population.Key words: epigenetics, imprinting, DMR, MIRA, MBD, DNA methylation, SNuPE     

Genome-wide DNA methylation analyses in the brain reveal four differentially methylated regions between humans and non-human primates

ABSTRACT: BACKGROUND: The highly improved cognitive function is the most significant change in human evolutionary history. Recently, several large-scale studies reported the evolutionary roles of DNA methylation; however, the role of DNA methylation on brain evolution is largely unknown. RESULTS: To test if DNA methylation has contributed to the evolution of human brain, with the use of MeDIP-Chip and SEQUENOM MassARRAY, we conducted a genome-wide analysis to identify differentially methylated regions (DMRs) in the brain between humans and rhesus macaques. We first identified a total of 150 candidate DMRs by the MeDIP-Chip method, among which 4DMRs were confirmed by the MassARRAY analysis. All 4 DMRs are within or close to the CpG islands, and a MIR3 repeat element was identified in one DMR, but no repeat sequence was observed in the other 3 DMRs. For the 4 DMR genes, their proteins tend to be conserved and two genes have neural related functions. Bisulfite sequencing and phylogenetic comparison among human, chimpanzee, rhesus macaque and rat suggested several regions of lineage specific DNA methylation, including a human specific hypomethylated region in the promoter of K6IRS2 gene. CONCLUSIONS: Our study provides a new angle of studying human brain evolution and understanding the evolutionary role of DNA methylation in the central nervous system. The results suggest that the patterns of DNA methylation in the brain are in general similar between humans and nonhuman primates, and only a few DMRs were identified.     

Extensive tissue-specific variation of allelic methylation in the Igf2 gene during mouse fetal development: relation to expression and imprinting

The imprinted Igf2 gene is active only on the paternal allele in most tissues. Its imprinting involves a cis-acting imprinting-control region (ICR) located upstream of the neighboring and maternally expressed H19 gene. It is thought that differential methylation of the parental alleles at the ICR is crucial for parental imprinting of both genes. Differentially methylated regions (DMRs) have also been identified within the Igf2 gene and their differential methylation is thought to be established during early development. To gain further insight into the function of these DMRs, we performed a quantitative analysis of their allelic methylation levels in different tissues during fetal development and the postnatal period in the mouse. Surprisingly, we found that the methylation levels of Igf2 DMRs vary extensively during fetal development, mostly on the expressed paternal allele. In particular, in skeletal muscle, differential allelic methylation in both DMR 1 and DMR 2 occurs only after birth, whereas correct paternal monoallelic expression is always observed, including in the embryonic stages. This suggests that differential methylation in the DMR 1 and DMR 2 of the Igf2 gene is dispensable for its imprinting in skeletal muscle. Furthermore, progressive methylation of the Igf2 paternal allele appears to be correlated with concomitant postnatal down-regulation and silencing of the gene. We discuss possible relations between Igf2 allelic methylation and expression during fetal development.     

Fetal-specific DNA methylation ratio permits noninvasive prenatal diagnosis of trisomy 21

The trials performed worldwide toward noninvasive prenatal diagnosis (NIPD) of Down's syndrome (or trisomy 21) have shown the commercial and medical potential of NIPD compared to the currently used invasive prenatal diagnostic procedures. Extensive investigation of methylation differences between the mother and the fetus has led to the identification of differentially methylated regions (DMRs). In this study, we present a strategy using the methylated DNA immunoprecipitation (MeDiP) methodology in combination with real-time quantitative PCR (qPCR) to achieve fetal chromosome dosage assessment, which can be performed noninvasively through the analysis of fetal-specific DMRs. We achieved noninvasive prenatal detection of trisomy 21 by determining the methylation ratio of normal and trisomy 21 cases for each tested fetal-specific DMR present in maternal peripheral blood, followed by further statistical analysis. The application of this fetal-specific methylation ratio approach provided correct diagnosis of 14 trisomy 21 and 26 normal cases.     

Genome-Wide DNA Methylation Patterns and Transcription Analysis in Sheep Muscle

DNA methylation plays a central role in regulating many aspects of growth and development in mammals through regulating gene expression. The development of next generation sequencing technologies have paved the way for genome-wide, high resolution analysis of DNA methylation landscapes using methodology known as reduced representation bisulfite sequencing (RRBS). While RRBS has proven to be effective in understanding DNA methylation landscapes in humans, mice, and rats, to date, few studies have utilised this powerful method for investigating DNA methylation in agricultural animals. Here we describe the utilisation of RRBS to investigate DNA methylation in sheep Longissimus dorsi muscles. RRBS analysis of ∼1% of the genome from Longissimus dorsi muscles provided data of suitably high precision and accuracy for DNA methylation analysis, at all levels of resolution from genome-wide to individual nucleotides. Combining RRBS data with mRNAseq data allowed the sheep Longissimus dorsi muscle methylome to be compared with methylomes from other species. While some species differences were identified, many similarities were observed between DNA methylation patterns in sheep and other more commonly studied species. The RRBS data presented here highlights the complexity of epigenetic regulation of genes. However, the similarities observed across species are promising, in that knowledge gained from epigenetic studies in human and mice may be applied, with caution, to agricultural species. The ability to accurately measure DNA methylation in agricultural animals will contribute an additional layer of information to the genetic analyses currently being used to maximise production gains in these species.     

Allele-specific histone modifications regulate expression of the Dlk1-Gtl2 imprinted domain

Dlk1 and Gtl2 are reciprocally expressed imprinted genes located on mouse chromosome 12. The Dlk1-Gtl2 locus carries three differentially methylated regions (DMRs), which are methylated only on the paternal allele. Of these, the intergenic (IG) DMR, located 12 kb upstream of Gtl2, is required for proper imprinting of linked genes on the maternal chromosome, while the Gtl2 DMR, located across the promoter of the Gtl2 gene, is implicated in imprinting on both parental chromosomes. In addition to DNA methylation, modification of histone proteins is also an important regulator of imprinted gene expression. Chromatin immunoprecipitation was therefore used to examine the pattern of histone modifications across the IG and Gtl2 DMRs. The data show maternal-specific histone acetylation at the Gtl2 DMR, but not at the IG DMR. In contrast, only low levels of histone methylation were observed throughout the region, and there was no difference between the two parental alleles. An existing mouse line carrying a deletion/insertion upstream of Gtl2 is unable to imprint the Dlk1-Gtl2 locus properly and demonstrates loss of allele-specific methylation at the Gtl2 DMR. Further analysis of these animals now shows that the loss of allele-specific methylation is accompanied by increased paternal histone acetylation at the Gtl2 DMR, with the activated paternal allele adopting a maternal acetylation pattern. These data indicate that interactions between DNA methylation and histone acetylation are involved in regulating the imprinting of the Dlk1-Gtl2 locus.     

Differentially methylated regions of imprinted genes in prenatal, perinatal and postnatal human tissues

Epigenetic plasticity in relation to in utero exposures may mechanistically explain observed differences in the likelihood of developing common complex diseases including hypertension, diabetes and cardiovascular disease through the cumulative effects of subtle alterations in gene expression. Imprinted genes are essential mediators of growth and development and are characterized by differentially methylated regulatory regions (DMRs) that carry parental allele-specific methylation profiles. This theoretical 50% level of methylation provides a baseline from which endogenously- or exogenously-induced deviations in methylation can be detected. We quantified DNA methylation at imprinted gene DMRs in a large panel of human conceptal tissues, in matched buccal cell specimens collected at birth and at one year of age, and in the major cell fractions of umbilical cord blood to assess the stability of methylation at these regions. DNA methylation was measured using validated pyrosequencing assays at seven DMRs regulating the IGF2/H19, DLK1/MEG3, MEST, NNAT and SGCE/PEG10 imprinted domains. DMR methylation did not significantly differ for the H19, MEST and SGCE/PEG10 DMRs across all conceptal tissues analyzed (ANOVA p>0.10). Methylation differences at several DMRs were observed in tissues from brain (IGF2 and MEG3-IG DMRs), liver (IGF2 and MEG3 DMRs) and placenta (both DLK1/MEG3 DMRs and NNAT DMR). In most infants, methylation profiles in buccal cells at birth and at one year of age were comparable, as was methylation in the major cell fractions of umbilical cord blood. Several infants showed temporal deviations in methylation at multiple DMRs. Similarity of inter-individual and intra-individual methylation at some, but not all of the DMRs analyzed supports the possibility that methylation of these regions can serve as useful biosensors of exposure.     

Systematic assessment of reduced representation bisulfite sequencing to human blood samples: A promising method for large-sample-scale epigenomic studies

Complementary to the time- and cost-intensive direct bisulfite sequencing, we applied reduced representation bisulfite sequencing (RRBS) to the human peripheral blood mononuclear cells (PBMC) from YH, the Asian individual whose genome and epigenome has been deciphered in the YH project and systematically assessed the genomic coverage, coverage depth and reproducibility of this technology as well as the concordance of DNA methylation levels measured by RRBS and direct bisulfite sequencing for the detected CpG sites. Our result suggests that RRBS can cover more than half of CpG islands and promoter regions with a good coverage depth and the proportion of the CpG sites covered by the biological replicates reaches 80-90%, indicating good reproducibility. Given a smaller data quantity, RRBS enjoys much better coverage depth than direct bisulfite sequencing and the concordance of DNA methylation levels between the two methods is high. It can be concluded that RRBS is a time and cost-effective sequencing method for unbiased DNA methylation profiling of CpG islands and promoter regions in a genome-wide scale and it is the method of choice to assay certain genomic regions for multiple samples in a rapid way.     

Influence of the Prader-Willi syndrome imprinting center on the DNA methylation landscape in the mouse brain

Reduced representation bisulfite sequencing (RRBS) was used to analyze DNA methylation patterns across the mouse brain genome in mice carrying a deletion of the Prader-Willi syndrome imprinting center (PWS-IC) on either the maternally- or paternally-inherited chromosome. Within the ∼3.7 Mb imprinted Angelman/Prader-Willi syndrome (AS/PWS) domain, 254 CpG sites were interrogated for changes in methylation due to PWS-IC deletion. Paternally-inherited deletion of the PWS-IC increased methylation levels ∼2-fold at each CpG site (compared to wild-type controls) at differentially methylated regions (DMRs) associated with 5′ CpG island promoters of paternally-expressed genes; these methylation changes extended, to a variable degree, into the adjacent CpG island shores. Maternal PWS-IC deletion yielded little or no changes in methylation at these DMRs, and methylation of CpG sites outside of promoter DMRs also was unchanged upon maternal or paternal PWS-IC deletion. Using stringent ascertainment criteria, ∼750,000 additional CpG sites were also interrogated across the entire mouse genome. This analysis identified 26 loci outside of the imprinted AS/PWS domain showing altered DNA methylation levels of ≥25% upon PWS-IC deletion. Curiously, altered methylation at 9 of these loci was a consequence of maternal PWS-IC deletion (maternal PWS-IC deletion by itself is not known to be associated with a phenotype in either humans or mice), and 10 of these loci exhibited the same changes in methylation irrespective of the parental origin of the PWS-IC deletion. These results suggest that the PWS-IC may affect DNA methylation at these loci by directly interacting with them, or may affect methylation at these loci through indirect downstream effects due to PWS-IC deletion. They further suggest the PWS-IC may have a previously uncharacterized function outside of the imprinted AS/PWS domain.     

Influence of the Prader-Willi syndrome imprinting center on the DNA methylation landscape in the mouse brain

Reduced representation bisulfite sequencing (RRBS) was used to analyze DNA methylation patterns across the mouse brain genome in mice carrying a deletion of the Prader-Willi syndrome imprinting center (PWS-IC) on either the maternally- or paternally-inherited chromosome. Within the ～3.7 Mb imprinted Angelman/Prader-Willi syndrome (AS/PWS) domain, 254 CpG sites were interrogated for changes in methylation due to PWS-IC deletion. Paternally-inherited deletion of the PWS-IC increased methylation levels ～2-fold at each CpG site (compared to wild-type controls) at differentially methylated regions (DMRs) associated with 5′ CpG island promoters of paternally-expressed genes; these methylation changes extended, to a variable degree, into the adjacent CpG island shores. Maternal PWS-IC deletion yielded little or no changes in methylation at these DMRs, and methylation of CpG sites outside of promoter DMRs also was unchanged upon maternal or paternal PWS-IC deletion. Using stringent ascertainment criteria, ～750,000 additional CpG sites were also interrogated across the entire mouse genome. This analysis identified 26 loci outside of the imprinted AS/PWS domain showing altered DNA methylation levels of ≥25% upon PWS-IC deletion. Curiously, altered methylation at 9 of these loci was a consequence of maternal PWS-IC deletion (maternal PWS-IC deletion by itself is not known to be associated with a phenotype in either humans or mice), and 10 of these loci exhibited the same changes in methylation irrespective of the parental origin of the PWS-IC deletion. These results suggest that the PWS-IC may affect DNA methylation at these loci by directly interacting with them, or may affect methylation at these loci through indirect downstream effects due to PWS-IC deletion. They further suggest the PWS-IC may have a previously uncharacterized function outside of the imprinted AS/PWS domain.     

Reduced representation bisulfite sequencing design for assessing the methylation of human CpG islands in large samples

The reduced representation bisulfite sequencing (RRBS) method has been developed for the high-throughput analysis of DNA methylation based on the sequencing of genomic libraries treated with sodium bisulfite by next-generation approaches. In contrast to whole-genome sequencing, the RRBS approach elaborates specific endonucleases to prepare libraries in order to produce pools of CpG-rich DNA fragments. The original RRBS technology based on the use of the MspI libraries allows one to increase the relative number of CpG islands in the pools of genomic fragments compared to whole-genome bisulfite sequencing. Nevertheless, this technology is rarely used due to the high cost compared with bisulfite methylation analysis with hybridization microarrays and significant residual amount of data represented by the sequences of genomic repeats that complicates the alignment and is not of particular interest for developing DNA methylation markers, which is often the main goal of biomedical research. We have developed an algorithm for estimating the likelihood that recognition sites of restriction endonucleases will be represented in CpG islands and present a method of reducing the effective size of the RRBS library without a significant loss of the CpG islands based on the use of the XmaI endonuclease for library preparation. In silico analysis demonstrates that the optimum range of the XmaI-RRBS fragment lengths is 110–200 base pairs. The sequencing of this library allows one to assess the methylation status of over 125000 CpG dinucleotides, of which over 90000 belong to CpG islands.

     

Bovine DNA methylation imprints are established in an oocyte size-specific manner, which are coordinated with the expression of the DNMT3 family proteins

A subset of genes, known as imprinted genes, is present in the mammalian genome. Genomic imprinting governs the monoallelic expression of these genes, depending on whether the gene was inherited from the sperm or the egg. This parent-of-origin specific gene expression is generally dependent on the epigenetic modification, DNA methylation, and the DNA methylation status of CpG dinucleotides residing in loci known as differentially methylated regions (DMRs). The enzymatic machinery responsible for the addition of methyl (-CH(3)) groups to the cytosine residue in the CpG dinucleotides are known as DNA methyltransferases (DNMTs). Correct establishment and maintenance of methylation patterns at imprinted genes has been associated with placental function and regulation of embryonic/fetal development. Much work has been carried out on imprinted genes in mouse and human; however, little is known about the methylation dynamics in the bovine oocyte. The primary objective of the present study was to characterize the establishment of methylation at maternally imprinted genes in bovine growing oocytes and to determine if the expression of the bovine DNMTs-DNMT3A, DNMT3B, and DNMT3L-was coordinated with DNA methylation during oocyte development. To this end, a panel of maternally imprinted genes was selected (SNRPN, MEST, IGF2R, PEG10, and PLAGL1) and putative DMRs for MEST, IGF2R, PEG10, and PLAGL1 were identified within the 5' regions for each gene; the SNRPN DMR has been reported previously. Conventional bisulfite sequencing revealed that methylation marks were acquired at all five DMRs investigated in an oocyte size-dependent fashion. This was confirmed for a selection of genes using pyrosequencing analysis. Furthermore, mRNA expression and protein analysis revealed that DNMT3A, DNMT3B, and DNMT3L are also present in the bovine oocyte during its growth phase. This study demonstrates for the first time that an increase in bovine imprinted gene DMR methylation occurs during oocyte growth, as is observed in mouse.     

Autonomous silencing of the imprinted Cdkn1c gene in stem cells

Parent-of-origin specific expression of imprinted genes relies on the differential DNA methylation of specific genomic regions. Differentially methylated regions (DMRs) acquire DNA methylation either during gametogenesis (primary DMR) or after fertilisation when allele-specific expression is established (secondary DMR). Little is known about the function of these secondary DMRs. We investigated the DMR spanning Cdkn1c in mouse embryonic stem cells, androgenetic stem cells and embryonic germ stem cells. In all cases, expression of Cdkn1c was appropriately repressed in in vitro differentiated cells. However, stem cells failed to de novo methylate the silenced gene even after sustained differentiation. In the absence of maintained DNA methylation (Dnmt1^-/-), Cdkn1c escapes silencing demonstrating the requirement for DNA methylation in long term silencing in vivo. We propose that postfertilisation differential methylation reflects the importance of retaining single gene dosage of a subset of imprinted loci in the adult.