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.     

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.     

Comparison of novel and existing methods for detecting differentially methylated regions

Background

Single-probe analyses in epigenome-wide association studies (EWAS) have identified associations between DNA methylation and many phenotypes, but do not take into account information from neighboring probes. Methods to detect differentially methylated regions (DMRs) (clusters of neighboring probes associated with a phenotype) may provide more power to detect associations between DNA methylation and diseases or phenotypes of interest.

Results

We proposed a novel approach, GlobalP, and perform comparisons with 3 methods—DMRcate, Bumphunter, and comb-p—to identify DMRs associated with log triglycerides (TGs) in real GAW20 data before and after fenofibrate treatment. We applied these methods to the summary statistics from an EWAS performed on the methylation data. Comb-p, DMRcate, and GlobalP detected very similar DMRs near the gene CPT1A on chromosome 11 in both the pre- and posttreatment data. In addition, GlobalP detected 2 DMRs before fenofibrate treatment in the genes ETV6 and ABCG1. Bumphunter identified several DMRs on chromosomes 1 and 20, which did not overlap with DMRs detected by other methods.

Conclusions

Our novel method detected the same DMR identified by two existing methods and detected two additional DMRs not identified by any of the existing methods we compared.

     

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.     

The regulation of LncRNA GTL2 expression by DNA methylation during sheep skeletal muscle development

DNA methylation has crucial roles in regulating the expression of genes involved in skeletal muscle development. However, the DNA methylation pattern of lncRNA during sheep skeletal muscle development remains unclear. This study investigated previous WGBS and LncRNA data in skeletal muscle of sheep (fetus and adult). We then focused on LncRNA GTL2, which is differentially expressed in skeletal muscle and has multiple DMRs. We found that the expression level of GTL2 decreased with age. GTL2 DMRs methylation levels were significantly higher in adult muscle than in fetal muscle. After 5AZA treatment, GTL2 expression was significantly increased in a dose-dependent manner.The dCas9-DNMT3A-sgRNA significantly reduced the expression level of GTL2 in cells, but increased GTL2 DMR methylation levels. The above studies indicate that dCas9-DNMT3A can effectively increase the methylation level in the DMR region of GTL2, the expression level of GTL2 is regulated by DNA methylation during muscle 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.     

Dynamic stage-specific changes in imprinted differentially methylated regions during early mammalian development and prevalence of non-CpG methylation in oocytes

Mammalian imprinted genes are associated with differentially methylated regions (DMRs) that are CpG methylated on one of the two parental chromosomes. In mice, at least 21 DMRs acquire differential methylation in the germline and many of them act as imprint centres. We previously reported the physical extents of differential methylation at 15 DMRs in mouse embryos at 12.5 days postcoitum. To reveal the ontogeny of differential methylation, we determined and compared methylation patterns of the corresponding regions in sperm and oocytes. We found that the extent of the gametic DMRs differs significantly from that of the embryonic DMRs, especially in the case of paternal gametic DMRs. These results suggest that the gametic DMR sequences should be used to extract the features specifying methylation imprint establishment in the germline: from this analysis, we noted that the maternal gametic DMRs appear as unmethylated islands in male germ cells, which suggests a novel component in the mechanism of gamete-specific marking. Analysis of selected DMRs in blastocysts revealed dynamic changes in allelic methylation in early development, indicating that DMRs are not fully protected from the major epigenetic reprogramming events occurring during preimplantation development. Furthermore, we observed non-CpG methylation in oocytes, but not in sperm, which disappeared by the blastocyst stage. Non-CpG methylation was frequently found at maternally methylated DMRs as well as non-DMR regions, suggesting its prevalence in the oocyte genome. These results provide evidence for a unique methylation profile in oocytes and reveal the surprisingly dynamic nature of DMRs in the early embryo.     

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.     

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     

A Comparison of the Whole Genome Approach of MeDIP-Seq to the Targeted Approach of the Infinium HumanMethylation450 BeadChip? for Methylome Profiling

DNA methylation is one of the most studied epigenetic marks in the human genome, with the result that the desire to map the human methylome has driven the development of several methods to map DNA methylation on a genomic scale. Our study presents the first comparison of two of these techniques - the targeted approach of the Infinium HumanMethylation450 BeadChip® with the immunoprecipitation and sequencing-based method, MeDIP-seq. Both methods were initially validated with respect to bisulfite sequencing as the gold standard and then assessed in terms of coverage, resolution and accuracy. The regions of the methylome that can be assayed by both methods and those that can only be assayed by one method were determined and the discovery of differentially methylated regions (DMRs) by both techniques was examined. Our results show that the Infinium HumanMethylation450 BeadChip® and MeDIP-seq show a good positive correlation (Spearman correlation of 0.68) on a genome-wide scale and can both be used successfully to determine differentially methylated loci in RefSeq genes, CpG islands, shores and shelves. MeDIP-seq however, allows a wider interrogation of methylated regions of the human genome, including thousands of non-RefSeq genes and repetitive elements, all of which may be of importance in disease. In our study MeDIP-seq allowed the detection of 15,709 differentially methylated regions, nearly twice as many as the array-based method (8070), which may result in a more comprehensive study of the methylome.     

Novel region discovery method for Infinium 450K DNA methylation data reveals changes associated with aging in muscle and neuronal pathways

We describe a methodology for detecting differentially methylated regions (DMRs) and variably methylated regions (VMRs), in data from Infinium 450K arrays that are very widely used in epigenetic studies. Region detection is more specific than single CpG analysis as it increases the extent of common findings between studies, and is more powerful as it reduces the multiple testing problem inherent in epigenetic whole‐genome association studies (EWAS). In addition, results driven by single erroneous probes are removed. We have used multiple publicly available Infinium 450K data sets to generate a consensus list of DMRs for age, supporting the hypothesis that aging is associated with specific epigenetic modifications. The consensus aging DMRs are significantly enriched for muscle biogenesis pathways. We find a massive increase in VMRs with age and in regions of the genome associated with open chromatin and neurotransmission. Old age VMRs are significantly enriched for neurotransmission pathways. EWAS studies should investigate the role of this interindividual variation in DNA methylation, in the age‐associated diseases of sarcopenia and dementia.     

A tripartite paternally methylated region within the Gpr1-Zdbf2 imprinted domain on mouse chromosome 1 identified by meDIP-on-chip

The parent-of-origin specific expression of imprinted genes relies on DNA methylation of CpG-dinucleotides at differentially methylated regions (DMRs) during gametogenesis. To date, four paternally methylated DMRs have been identified in screens based on conventional approaches. These DMRs are linked to the imprinted genes H19, Gtl2 (IG-DMR), Rasgrf1 and, most recently, Zdbf2 which encodes zinc finger, DBF-type containing 2. In this study, we applied a novel methylated-DNA immunoprecipitation-on-chip (meDIP-on-chip) method to genomic DNA from mouse parthenogenetic- and androgenetic-derived stem cells and sperm and identified 458 putative DMRs. This included the majority of known DMRs. We further characterized the paternally methylated Zdbf2/ZDBF2 DMR. In mice, this extensive germ line DMR spanned 16 kb and possessed an unusual tripartite structure. Methylation was dependent on DNA methyltransferase 3a (Dnmt3a), similar to H19 DMR and IG-DMR. In both humans and mice, the adjacent gene, Gpr1/GPR1, which encodes a G-protein-coupled receptor 1 protein with transmembrane domain, was also imprinted and paternally expressed. The Gpr1-Zdbf2 domain was most similar to the Rasgrf1 domain as both DNA methylation and the actively expressed allele were in cis on the paternal chromosome. This work demonstrates the effectiveness of meDIP-on-chip as a technique for identifying DMRs.     

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.     

Comparative genome-wide DNA methylation analysis of colorectal tumor and matched normal tissues

Aberrant DNA methylation often occurs in colorectal cancer (CRC). In our study we applied a genome-wide DNA methylation analysis approach, MethylCap-seq, to map the differentially methylated regions (DMRs) in 24 tumors and matched normal colon samples. In total, 2687 frequently hypermethylated and 468 frequently hypomethylated regions were identified, which include potential biomarkers for CRC diagnosis. Hypermethylation in the tumor samples was enriched at CpG islands and gene promoters, while hypomethylation was distributed throughout the genome. Using epigenetic data from human embryonic stem cells, we show that frequently hypermethylated regions coincide with bivalent loci in human embryonic stem cells. DNA methylation is commonly thought to lead to gene silencing; however, integration of publically available gene expression data indicates that 75% of the frequently hypermethylated genes were most likely already lowly or not expressed in normal tissue. Collectively, our study provides genome-wide DNA methylation maps of CRC, comprehensive lists of DMRs, and gives insights into the role of aberrant DNA methylation in CRC formation.     

Methylation status of putative differentially methylated regions of porcine IGF2 and H19

This study was designed to identify the putative differentially methylated regions (DMRs) of the porcine imprinted genes insulin-like growth factor 2 and H19 (IGF2-H19), and to assess the genomic imprinting status of IGF2-H19 by identifying the methylation patterns of these regions in germ cells, and in tissues from porcine fetuses, an adult pig, as well as cloned offspring produced by somatic cell nuclear transfer (SCNT). Porcine IGF2-H19 DMRs exhibit a normal monoallelic methylation pattern (i.e., either the paternally- or the maternally derived allele is methylated) similar to the pattern observed for the same genes in the human and mice genomes. Examination of the methylation patterns of the IGF2-H19 DMRs revealed that the zinc finger protein binding sites CTCF1 and 2 did not exhibit differential methylation in both control and cloned offspring. In contrast, the CTCF3 and DMR2 loci of the IGF2 gene showed abnormal methylation in cloned offspring, but a normal differential or moderate methylation pattern in tissues from control offspring and an adult pig. Our data thus suggest that regulation of genomic imprinting at the porcine IGF2-H19 loci is conserved among species, and that the abnormal methylation pattern in the regulatory elements of imprinted genes may lead to an alteration in the coordinated expression of genes required for successful reprogramming, which, in consequence, may contribute to the low efficiency of porcine genome reprogramming induced by nuclear transfer.     

Comparative genome-wide DNA methylation analysis of colorectal tumor and matched normal tissues

Aberrant DNA methylation often occurs in colorectal cancer (CRC). In our study we applied a genome-wide DNA methylation analysis approach, MethylCap-seq, to map the differentially methylated regions (DMRs) in 24 tumors and matched normal colon samples. In total, 2687 frequently hypermethylated and 468 frequently hypomethylated regions were identified, which include potential biomarkers for CRC diagnosis. Hypermethylation in the tumor samples was enriched at CpG islands and gene promoters, while hypomethylation was distributed throughout the genome. Using epigenetic data from human embryonic stem cells, we show that frequently hypermethylated regions coincide with bivalent loci in human embryonic stem cells. DNA methylation is commonly thought to lead to gene silencing; however, integration of publically available gene expression data indicates that 75% of the frequently hypermethylated genes were most likely already lowly or not expressed in normal tissue. Collectively, our study provides genome-wide DNA methylation maps of CRC, comprehensive lists of DMRs, and gives insights into the role of aberrant DNA methylation in CRC formation.