A Bayesian deconvolution strategy for immunoprecipitation-based DNA methylome analysis

DNA methylation is an indispensible epigenetic modification required for regulating the expression of mammalian genomes. Immunoprecipitation-based methods for DNA methylome analysis are rapidly shifting the bottleneck in this field from data generation to data analysis, necessitating the development of better analytical tools. In particular, an inability to estimate absolute methylation levels remains a major analytical difficulty associated with immunoprecipitation-based DNA methylation profiling. To address this issue, we developed a cross-platform algorithm-Bayesian tool for methylation analysis (Batman)-for analyzing methylated DNA immunoprecipitation (MeDIP) profiles generated using oligonucleotide arrays (MeDIP-chip) or next-generation sequencing (MeDIP-seq). We developed the latter approach to provide a high-resolution whole-genome DNA methylation profile (DNA methylome) of a mammalian genome. Strong correlation of our data, obtained using mature human spermatozoa, with those obtained using bisulfite sequencing suggest that combining MeDIP-seq or MeDIP-chip with Batman provides a robust, quantitative and cost-effective functional genomic strategy for elucidating the function of DNA methylation.     

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.     

MethylAction: detecting differentially methylated regions that distinguish biological subtypes

DNA methylation differences capture substantial information about the molecular and gene-regulatory states among biological subtypes. Enrichment-based next generation sequencing methods such as MBD-isolated genome sequencing (MiGS) and MeDIP-seq are appealing for studying DNA methylation genome-wide in order to distinguish between biological subtypes. However, current analytic tools do not provide optimal features for analyzing three-group or larger study designs. MethylAction addresses this need by detecting all possible patterns of statistically significant hyper- and hypo- methylation in comparisons involving any number of groups. Crucially, significance is established at the level of differentially methylated regions (DMRs), and bootstrapping determines false discovery rates (FDRs) associated with each pattern. We demonstrate this functionality in a four-group comparison among benign prostate and three clinical subtypes of prostate cancer and show that the bootstrap FDRs are highly useful in selecting the most robust patterns of DMRs. Compared to existing tools that are limited to two-group comparisons, MethylAction detects more DMRs with strong differential methylation measurements confirmed by whole genome bisulfite sequencing and offers a better balance between precision and recall in cross-cohort comparisons. MethylAction is available as an R package at http://jeffbhasin.github.io/methylaction.     

Genome wide analysis of acute myeloid leukemia reveal leukemia specific methylome and subtype specific hypomethylation of repeats

Methylated DNA immunoprecipitation followed by high-throughput sequencing (MeDIP-seq) has the potential to identify changes in DNA methylation important in cancer development. In order to understand the role of epigenetic modulation in the development of acute myeloid leukemia (AML) we have applied MeDIP-seq to the DNA of 12 AML patients and 4 normal bone marrows. This analysis revealed leukemia-associated differentially methylated regions that included gene promoters, gene bodies, CpG islands and CpG island shores. Two genes (SPHKAP and DPP6) with significantly methylated promoters were of interest and further analysis of their expression showed them to be repressed in AML. We also demonstrated considerable cytogenetic subtype specificity in the methylomes affecting different genomic features. Significantly distinct patterns of hypomethylation of certain interspersed repeat elements were associated with cytogenetic subtypes. The methylation patterns of members of the SINE family tightly clustered all leukemic patients with an enrichment of Alu repeats with a high CpG density (P<0.0001). We were able to demonstrate significant inverse correlation between intragenic interspersed repeat sequence methylation and gene expression with SINEs showing the strongest inverse correlation (R(2) = 0.7). We conclude that the alterations in DNA methylation that accompany the development of AML affect not only the promoters, but also the non-promoter genomic features, with significant demethylation of certain interspersed repeat DNA elements being associated with AML cytogenetic subtypes. MeDIP-seq data were validated using bisulfite pyrosequencing and the Infinium array.     

Folic acid supplementation alters the DNA methylation profile and improves insulin resistance in high-fat-diet-fed mice

Folic acid (FA) supplementation may protect from obesity and insulin resistance, the effects and mechanism of FA on chronic high-fat-diet-induced obesity-related metabolic disorders are not well elucidated. We adopted a genome-wide approach to directly examine whether FA supplementation affects the DNA methylation profile of mouse adipose tissue and identify the functional consequences of these changes. Mice were fed a high-fat diet (HFD), normal diet (ND) or an HFD supplemented with folic acid (20 μg/ml in drinking water) for 10 weeks, epididymal fat was harvested, and genome-wide DNA methylation analyses were performed using methylated DNA immunoprecipitation sequencing (MeDIP-seq). Mice exposed to the HFD expanded their adipose mass, which was accompanied by a significant increase in circulating glucose and insulin levels. FA supplementation reduced the fat mass and serum glucose levels and improved insulin resistance in HFD-fed mice. MeDIP-seq revealed distribution of differentially methylated regions (DMRs) throughout the adipocyte genome, with more hypermethylated regions in HFD mice. Methylome profiling identified DMRs associated with 3787 annotated genes from HFD mice in response to FA supplementation. Pathway analyses showed novel DNA methylation changes in adipose genes associated with insulin secretion, pancreatic secretion and type 2 diabetes. The differential DNA methylation corresponded to changes in the adipose tissue gene expression of Adcy3 and Rapgef4 in mice exposed to a diet containing FA. FA supplementation improved insulin resistance, decreased the fat mass, and induced DNA methylation and gene expression changes in genes associated with obesity and insulin secretion in obese mice fed a HFD.     

Comparative DNA methylome analysis of endometrial carcinoma reveals complex and distinct deregulation of cancer promoters and enhancers

Background

Aberrant DNA methylation is a hallmark of many cancers. Classically there are two types of endometrial cancer, endometrioid adenocarcinoma (EAC), or Type I, and uterine papillary serous carcinoma (UPSC), or Type II. However, the whole genome DNA methylation changes in these two classical types of endometrial cancer is still unknown.

Results

Here we described complete genome-wide DNA methylome maps of EAC, UPSC, and normal endometrium by applying a combined strategy of methylated DNA immunoprecipitation sequencing (MeDIP-seq) and methylation-sensitive restriction enzyme digestion sequencing (MRE-seq). We discovered distinct genome-wide DNA methylation patterns in EAC and UPSC: 27,009 and 15,676 recurrent differentially methylated regions (DMRs) were identified respectively, compared with normal endometrium. Over 80% of DMRs were in intergenic and intronic regions. The majority of these DMRs were not interrogated on the commonly used Infinium 450K array platform. Large-scale demethylation of chromosome X was detected in UPSC, accompanied by decreased XIST expression. Importantly, we discovered that the majority of the DMRs harbored promoter or enhancer functions and are specifically associated with genes related to uterine development and disease. Among these, abnormal methylation of transposable elements (TEs) may provide a novel mechanism to deregulate normal endometrium-specific enhancers derived from specific TEs.

Conclusions

DNA methylation changes are an important signature of endometrial cancer and regulate gene expression by affecting not only proximal promoters but also distal enhancers.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-868) contains supplementary material, which is available to authorized users.     

A streamlined method for analysing genome-wide DNA methylation patterns from low amounts of FFPE DNA

Background

Formalin fixed paraffin embedded (FFPE) tumor samples are a major source of DNA from patients in cancer research. However, FFPE is a challenging material to work with due to macromolecular fragmentation and nucleic acid crosslinking. FFPE tissue particularly possesses challenges for methylation analysis and for preparing sequencing-based libraries relying on bisulfite conversion. Successful bisulfite conversion is a key requirement for sequencing-based methylation analysis.

Methods

Here we describe a complete and streamlined workflow for preparing next generation sequencing libraries for methylation analysis from FFPE tissues. This includes, counting cells from FFPE blocks and extracting DNA from FFPE slides, testing bisulfite conversion efficiency with a polymerase chain reaction (PCR) based test, preparing reduced representation bisulfite sequencing libraries and massively parallel sequencing.

Results

The main features and advantages of this protocol are:

• An optimized method for extracting good quality DNA from FFPE tissues.
• An efficient bisulfite conversion and next generation sequencing library preparation protocol that uses 50 ng DNA from FFPE tissue.
• Incorporation of a PCR-based test to assess bisulfite conversion efficiency prior to sequencing.

Conclusions

We provide a complete workflow and an integrated protocol for performing DNA methylation analysis at the genome-scale and we believe this will facilitate clinical epigenetic research that involves the use of FFPE tissue.

     

DNA methylation analysis of murine hematopoietic side population cells during aging

Stem cells have been found in most tissues/organs. These somatic stem cells produce replacements for lost and damaged cells, and it is not completely understood how this regenerative capacity becomes diminished during aging. To study the possible involvement of epigenetic changes in somatic stem cell aging, we used murine hematopoiesis as a model system. Hematopoietic stem cells (HSCs) were enriched for via Hoechst exclusion activity (SP-HSC) from young, medium-aged and old mice and subjected to comprehensive, global methylome (MeDIP-seq) analysis. With age, we observed a global loss of DNA methylation of approximately 5%, but an increase in methylation at some CpG islands. Just over 100 significant (FDR &lt; 0.2) aging-specific differentially methylated regions (aDMRs) were identified, which are surprisingly few considering the profound age-based changes that occur in HSC biology. Interestingly, the polycomb repressive complex -2 (PCRC2) target genes Kiss1r, Nav2 and Hsf4 were hypermethylated with age. The promoter for the Sdpr gene was determined to be progressively hypomethylated with age. This occurred concurrently with an increase in gene expression with age. To explore this relationship further, we cultured isolated SP-HSC in the presence of 5-aza-deoxycytdine and demonstrated a negative correlation between Sdpr promoter methylation and gene expression. We report that DNA methylation patterns are well preserved during hematopoietic stem cell aging, confirm that PCRC2 targets are increasingly methylated with age, and suggest that SDPR expression changes with age in HSCs may be regulated via age-based alterations in DNA methylation.     

The Mechanism of Environmental Endocrine Disruptors (DEHP) Induces Epigenetic Transgenerational Inheritance of Cryptorchidism

Discussion on the role of DEHP in the critical period of gonadal development in pregnant rats (F0), studied the evolution of F1-F4 generation of inter-generational inheritance of cryptorchidism and the alteration of DNA methylation levels in testis. Pregnant SD rats were randomly divided into two groups: normal control group and DEHP experimental group. From pregnancy 7d to 19d, experimental group was sustained to gavage DEHP 750mg/kg bw/day, observed the incidence of cryptorchidism in offspring and examined the pregnancy rate of female rats through mating experiments. Continuous recording the rat’s weight and AGD value, after maturation (PND80) recording testis and epididymis’ size and weight, detected the sperm number and quality. Subsequently, we examined the evolution morphological changes of testicular tissue for 4 generation rats by HE staining and Western Blot. Completed the MeDIP-sequencing analysis of 6 samples (F1 generation, F4 generation and Control). DEHP successfully induced cryptorchidism occurrence in offspring during pregnancy. The incidence of cryptorchidism in F1 was 30%, in F2 was 12.5%, and there was no cryptorchidism coming up in F3 and F4. Mating experiment shows conception rate 50% in F1, F2 generation was 75%, the F3 and F4 generation were 100%. HE staining showed that the seminiferous epithelium of F1 generation was atrophy and with a few spermatogenic cell, F2 generation had improved, F3 and F4 generation were tend to be normal. The DNA methyltransferase expression was up-regulated with the increase of generations by Real Time-PCR, immunohistochemistry and Western Blot. MeDIP-seq Data Analysis Results show many differentially methylated DNA sequences between F1 and F4. DEHP damage male reproductive function in rats, affect expression of DNA methyltransferase enzyme, which in turn leads to genomic imprinting methylation pattern changes and passed on to the next generation, so that the offspring of male reproductive system critical role in the development of imprinted genes imbalances, and eventually lead to producing offspring cryptorchidism. This may be an important mechanism of reproductive system damage.     

Assessment of RainDrop BS-seq as a method for large-scale,targeted bisulfite sequencing

We present a systematic assessment of RainDrop BS-seq, a novel method for large-scale, targeted bisulfite sequencing using microdroplet-based PCR amplification coupled with next-generation sequencing. We compared DNA methylation levels at 498 target loci (1001 PCR amplicons) in human whole blood, osteosarcoma cells and an archived tumor tissue sample. We assessed the ability of RainDrop BS-seq to accurately measure DNA methylation over a range of DNA quantities (from 10 to 1500 ng), both with and without whole-genome amplification (WGA) following bisulfite conversion. DNA methylation profiles generated using at least 100 ng correlated well (median R = 0.92) with those generated on Illumina Infinium HumanMethylation450 BeadChips, currently the platform of choice for epigenome-wide association studies (EWAS). WGA allowed for testing of samples with a starting DNA amount of 10 and 50 ng, although a reduced correlation was observed (median R = 0.79). We conclude that RainDrop BS-seq is suitable for measuring DNA methylation levels using nanogram quantities of DNA, and can be used to study candidate epigenetic biomarker loci in an accurate and high-throughput manner, paving the way for its application to routine clinical diagnostics.     

Assessment of RainDrop BS-seq as a method for large-scale,targeted bisulfite sequencing

We present a systematic assessment of RainDrop BS-seq, a novel method for large-scale, targeted bisulfite sequencing using microdroplet-based PCR amplification coupled with next-generation sequencing. We compared DNA methylation levels at 498 target loci (1001 PCR amplicons) in human whole blood, osteosarcoma cells and an archived tumor tissue sample. We assessed the ability of RainDrop BS-seq to accurately measure DNA methylation over a range of DNA quantities (from 10 to 1500 ng), both with and without whole-genome amplification (WGA) following bisulfite conversion. DNA methylation profiles generated using at least 100 ng correlated well (median R = 0.92) with those generated on Illumina Infinium HumanMethylation450 BeadChips, currently the platform of choice for epigenome-wide association studies (EWAS). WGA allowed for testing of samples with a starting DNA amount of 10 and 50 ng, although a reduced correlation was observed (median R = 0.79). We conclude that RainDrop BS-seq is suitable for measuring DNA methylation levels using nanogram quantities of DNA, and can be used to study candidate epigenetic biomarker loci in an accurate and high-throughput manner, paving the way for its application to routine clinical diagnostics.     

A MBD-seq protocol for large-scale methylome-wide studies with (very) low amounts of DNA

We recently showed that, after optimization, our methyl-CpG binding domain sequencing (MBD-seq) application approximates the methylome-wide coverage obtained with whole-genome bisulfite sequencing (WGB-seq), but at a cost that enables adequately powered large-scale association studies. A prior drawback of MBD-seq is the relatively large amount of genomic DNA (ideally >1 µg) required to obtain high-quality data. Biomaterials are typically expensive to collect, provide a finite amount of DNA, and may simply not yield sufficient starting material. The ability to use low amounts of DNA will increase the breadth and number of studies that can be conducted. Therefore, we further optimized the enrichment step. With this low starting material protocol, MBD-seq performed equally well, or better, than the protocol requiring ample starting material (>1 µg). Using only 15 ng of DNA as input, there is minimal loss in data quality, achieving 93% of the coverage of WGB-seq (with standard amounts of input DNA) at similar false/positive rates. Furthermore, across a large number of genomic features, the MBD-seq methylation profiles closely tracked those observed for WGB-seq with even slightly larger effect sizes. This suggests that MBD-seq provides similar information about the methylome and classifies methylation status somewhat more accurately. Performance decreases with <15 ng DNA as starting material but, even with as little as 5 ng, MBD-seq still achieves 90% of the coverage of WGB-seq with comparable genome-wide methylation profiles. Thus, the proposed protocol is an attractive option for adequately powered and cost-effective methylome-wide investigations using (very) low amounts of DNA.     

Comparison of the Genome-Wide DNA Methylation Profiles between Fast-Growing and Slow-Growing Broilers

Introduction

Growth traits are important in poultry production, however, little is known for its regulatory mechanism at epigenetic level. Therefore, in this study, we aim to compare DNA methylation profiles between fast- and slow-growing broilers in order to identify candidate genes for chicken growth. Methylated DNA immunoprecipitation-sequencing (MeDIP-seq) was used to investigate the genome-wide DNA methylation pattern in high and low tails of Recessive White Rock (WRR_h; WRR_l) and that of Xinhua Chickens (XH_h; XH_l) at 7 weeks of age. The results showed that the average methylation density was the lowest in CGIs followed by promoters. Within the gene body, the methylation density of introns was higher than that of UTRs and exons. Moreover, different methylation levels were observed in different repeat types with the highest in LINE/CR1. Methylated CGIs were prominently distributed in the intergenic regions and were enriched in the size ranging 200–300 bp. In total 13,294 methylated genes were found in four samples, including 4,085 differentially methylated genes of WRR_h Vs. WRR_l, 5,599 of XH_h Vs. XH_l, 4,204 of WRR_h Vs. XH_h, as well as 7,301 of WRR_l Vs. XH_l. Moreover, 132 differentially methylated genes related to growth and metabolism were observed in both inner contrasts (WRR_h Vs. WRR_l and XH_h Vs. XH_l), whereas 129 differentially methylated genes related to growth and metabolism were found in both across-breed contrasts (WRR_h Vs. XH_h and WRR_l Vs. XH_l). Further analysis showed that overall 75 genes exhibited altered DNA methylation in all four contrasts, which included some well-known growth factors of IGF1R, FGF12, FGF14, FGF18, FGFR2, and FGFR3. In addition, we validate the MeDIP-seq results by bisulfite sequencing in some regions.

Conclusions

This study revealed the global DNA methylation pattern of chicken muscle, and identified candidate genes that potentially regulate muscle development at 7 weeks of age at methylation level.     

Profound flanking sequence preference of Dnmt3a and Dnmt3b mammalian DNA methyltransferases shape the human epigenome

Mammalian DNA methyltransferases methylate cytosine residues within CG dinucleotides. By statistical analysis of published data of the Human Epigenome Project we have determined flanking sequences of up to +/-four base-pairs surrounding the central CG site that are characteristic of high (5'-CTTGCGCAAG-3') and low (5'-TGTTCGGTGG-3') levels of methylation in human genomic DNA. We have investigated the influence of flanking sequence on the catalytic activity of the Dnmt3a and Dnmt3b de novo DNA methyltransferases using a set of synthetic oligonucleotide substrates that covers all possible +/-1 flanks in quantitative terms. Methylation kinetics experiments revealed a >13-fold difference between the preferred (RCGY) and disfavored +/-1 flanking base-pairs (YCGR). In addition, AT-rich flanks are preferred over GC-rich ones. These experimental preferences coincide with the genomic methylation patterns. Therefore, we have expanded our experimental analysis and found a >500-fold difference in the methylation rates of the consensus sequences for high and low levels of methylation in the genome. This result demonstrates a very pronounced flanking sequence preference of Dnmt3a and Dnmt3b. It suggests that the methylation pattern of human DNA is due, in part, to the flanking sequence preferences of the de novo DNA MTases and that flanking sequence preferences could be involved in the origin of CG islands. Furthermore, similar flanking sequence preferences have been found for the stimulation of the immune system by unmethylated CGs, suggesting a co-evolution of DNA MTases and the immune system.