Genome-wide high throughput analysis of DNA methylation in eukaryotes

Cytosine methylation is the quintessential epigenetic mark. Two well-established methods, bisulfite sequencing and methyl-DNA immunoprecipitation (MeDIP) lend themselves to the genome-wide analysis of DNA methylation by high throughput sequencing. Here we provide an overview and brief review of these methods. We summarize our experience with MeDIP followed by high throughput Illumina/Solexa sequencing, exemplified by the analysis of the methylated fraction of the Neurospora crassa genome ("methylome"). We provide detailed methods for DNA isolation, processing and the generation of in vitro libraries for Illumina/Solexa sequencing. We discuss potential problems in the generation of sequencing libraries. Finally, we provide an overview of software that is appropriate for the analysis of high throughput sequencing data generated by Illumina/Solexa-type sequencing by synthesis, with a special emphasis on approaches and applications that can generate more accurate depictions of sequence reads that fall in repeated regions of a chosen reference genome.     

Standardization and quality controls for the methylated DNA immunoprecipitation technique

MeDIP (Methylated DNA Immunoprecipitation) is a relatively recent technique aimed to enrich the methylated fraction of DNA with an antibody directed against 5-methyl-cytosine. MeDIP processed samples are suitable for investigation of the methylation status of specific genomic loci and for performing genome-wide screening when hybridized to DNA methylation microarrays or analyzed by deep sequencing. Here, we describe a standardization protocol and quality controls to assess the specificity, reproducibility and efficiency of the MeDIP procedure. These may have utility when comparing results between samples and experiments within laboratories and between laboratories.     

Exome localization of complex disease association signals

Background

Recent progress in high-throughput technologies has greatly contributed to the development of DNA methylation profiling. Although there are several reports that describe methylome detection of whole genome bisulfite sequencing, the high cost and heavy demand on bioinformatics analysis prevents its extensive application. Thus, current strategies for the study of mammalian DNA methylomes is still based primarily on genome-wide methylated DNA enrichment combined with DNA microarray detection or sequencing. Methylated DNA enrichment is a key step in a microarray based genome-wide methylation profiling study, and even for future high-throughput sequencing based methylome analysis.

Results

In order to evaluate the sensitivity and accuracy of methylated DNA enrichment, we investigated and optimized a number of important parameters to improve the performance of several enrichment assays, including differential methylation hybridization (DMH), microarray-based methylation assessment of single samples (MMASS), and methylated DNA immunoprecipitation (MeDIP). With advantages and disadvantages unique to each approach, we found that assays based on methylation-sensitive enzyme digestion and those based on immunoprecipitation detected different methylated DNA fragments, indicating that they are complementary in their relative ability to detect methylation differences.

Conclusions

Our study provides the first comprehensive evaluation for widely used methodologies for methylated DNA enrichment, and could be helpful for developing a cost effective approach for DNA methylation profiling.     

Genome-wide mapping of DNA methylation in Nile Tilapia

Cytosine DNA methylation is crucial for gene regulation and maintenance of genome stability. However, the detailed nile tilapia methylome remains uncharacterized. In this study, we present the first high-resolution methylome of tilapia gonad generated using methylated DNA immunoprecipitation (MeDIP) and high-throughput sequencing. In the ovary, 265 and 56 methylation peaks were identified in the genebody and promoter region of 145 genes, respectively. In the testis, 293 and 80 methylation peaks were identified in the genebody and promoter region of 144 genes. Furthermore, 8 and 49 genes showed differentially higher and lower promoter-region methylation rates, respectively, in the ovary relative to those of the testis. Quantitative PCR results revealed that the expression level of fibroblast growth factor 16 (fgf16), sialidase-3-like, fibroblast growth factor 20, aromatase (cyp19a), estrogen receptor, and gonadotropin receptor II precursor were negatively correlated to their methylation levels in the ovary and testis. The methylated levels of cyp19a and fgf16 were validated by bisulfite sequencing PCR technology, and the results were consistent with the MeDIP results. Thus, apart from generating the first methylation map, this study produced a candidate gene repository that provides additional options to explore the relationship between DNA methylation and sex differentiation or maintenance.     

Clinical and public health research using methylated DNA Immunoprecipitation (MeDIP): A comparison of commercially available kits to examine differential DNA methylation across the genome

The methylated DNA immunoprecipitation method (MeDIP) is a genome-wide, high-resolution approach that detects DNA methylation with oligonucleotide tiling arrays or high throughput sequencing platforms. A simplified high-throughput MeDIP assay will enable translational research studies in clinics and populations, which will greatly enhance our understanding of the human methylome. We compared three commercial kits, MagMeDIP Kit TM (Diagenode), Methylated-DNA IP Kit (Zymo Research) and Methylamp? Methylated DNA Capture Kit (Epigentek), in order to identify which one has better reliability and sensitivity for genomic DNA enrichment. Each kit was used to enrich two samples, one from fresh tissue and one from a cell line, with two different DNA amounts. The enrichment efficiency of each kit was evaluated by agarose gel band intensity after Nco I digestion and by reaction yield of methylated DNA. A successful enrichment is expected to have a 1:4 to 10:1 conversion ratio and a yield of 80% or higher. We also evaluated the hybridization efficiency to genome-wide methylation arrays in a separate cohort of tissue samples. We observed that the MagMeDIP kit had the highest yield for the two DNA amounts and for both the tissue and cell line samples, as well as for the positive control. In addition, the DNA was successfully enriched from a 1:4 to 10:1 ratio. Therefore, the MagMeDIP kit is a useful research tool that will enable clinical and public health genome-wide DNA methylation studies.     

Clinical and public health research using methylated DNA immunoprecipitation (MeDIP): a comparison of commercially available kits to examine differential DNA methylation across the genome

The methylated DNA immunoprecipitation method (MeDIP) is a genome-wide, high-resolution approach that detects DNA methylation with oligonucleotide tiling arrays or high throughput sequencing platforms. A simplified high-throughput MeDIP assay will enable translational research studies in clinics and populations, which will greatly enhance our understanding of the human methylome. We compared three commercial kits, MagMeDIP Kit TM (Diagenode), Methylated-DNA IP Kit (Zymo Research) and Methylamp™ Methylated DNA Capture Kit (Epigentek), in order to identify which one has better reliability and sensitivity for genomic DNA enrichment. Each kit was used to enrich two samples, one from fresh tissue and one from a cell line, with two different DNA amounts. The enrichment efficiency of each kit was evaluated by agarose gel band intensity after Nco I digestion and by reaction yield of methylated DNA. A successful enrichment is expected to have a 1:4 to 10:1 conversion ratio and a yield of 80% or higher. We also evaluated the hybridization efficiency to genome-wide methylation arrays in a separate cohort of tissue samples. We observed that the MagMeDIP kit had the highest yield for the two DNA amounts and for both the tissue and cell line samples, as well as for the positive control. In addition, the DNA was successfully enriched from a 1:4 to 10:1 ratio. Therefore, the MagMeDIP kit is a useful research tool that will enable clinical and public health genome-wide DNA methylation studies.     

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.     

Examination of the specificity of DNA methylation profiling techniques towards 5-methylcytosine and 5-hydroxymethylcytosine

DNA cytosine-5 methylation is a well-studied epigenetic pathway implicated in gene expression control and disease pathogenesis. Different technologies have been developed to examine the distribution of 5-methylcytosine (5mC) in specific sequences of the genome. Recently, substantial amounts of 5-hydroxymethylcytosine (5hmC), most likely derived from enzymatic oxidation of 5mC by TET1, have been detected in certain mammalian tissues. Here, we have examined the ability of several commonly used DNA methylation profiling methods to distinguish between 5mC and 5hmC. We show that techniques based on sodium bisulfite treatment of DNA are incapable of distinguishing between the two modified bases. In contrast, techniques based on immunoprecipitation with anti-5mC antibody (methylated DNA immunoprecipitation, MeDIP) or those based on proteins that bind to methylated CpG sequences (e.g. methylated-CpG island recovery assay, MIRA) do not detect 5hmC and are specific for 5mC unless both modified bases occur in the same DNA fragment. We also report that several methyl-CpG binding proteins including MBD1, MBD2 and MBD4 do not bind to sequences containing 5hmC. Selective mapping of 5hmC will require the development of unique tools for the detection of this modified base.     

Methylated DNA Immunoprecipitation

The identification of DNA methylation patterns is a common procedure in the study of epigenetics, as methylation is known to have significant effects on gene expression, and is involved with normal development as well as disease ^1-4. Thus, the ability to discriminate between methylated DNA and non-methylated DNA is essential for generating methylation profiles for such studies. Methylated DNA immunoprecipitation (MeDIP) is an efficient technique for the extraction of methylated DNA from a sample of interest ^5-7. A sample of as little as 200 ng of DNA is sufficient for the antibody, or immunoprecipitation (IP), reaction. DNA is sonicated into fragments ranging in size from 300-1000 bp, and is divided into immunoprecipitated (IP) and input (IN) portions. IP DNA is subsequently heat denatured and then incubated with anti-5''mC, allowing the monoclonal antibody to bind methylated DNA. After this, magnetic beads containing a secondary antibody with affinity for the primary antibody are added, and incubated. These bead-linked antibodies will bind the monoclonal antibody used in the first step. DNA bound to the antibody complex (methylated DNA) is separated from the rest of the DNA by using a magnet to pull the complexes out of solution. Several washes using IP buffer are then performed to remove the unbound, non-methylated DNA. The methylated DNA/antibody complexes are then digested with Proteinase K to digest the antibodies leaving only the methylated DNA intact. The enriched DNA is purified by phenol:chloroform extraction to remove the protein matter and then precipitated and resuspended in water for later use. PCR techniques can be used to validate the efficiency of the MeDIP procedure by analyzing the amplification products of IP and IN DNA for regions known to lack and known to contain methylated sequences. The purified methylated DNA can then be used for locus-specific (PCR) or genome-wide (microarray and sequencing) methylation studies, and is particularly useful when applied in conjunction with other research tools such as gene expression profiling and array comparative genome hybridization (CGH) ⁸. Further investigation into DNA methylation will lead to the discovery of new epigenetic targets, which in turn, may be useful in developing new therapeutic or prognostic research tools for diseases such as cancer that are characterized by aberrantly methylated DNA ^{2, 4, 9-11}.     

Gestational stress induces depressive-like and anxiety-like phenotypes through epigenetic regulation of BDNF expression in offspring hippocampus

Exposure to stressful life events during pregnancy exerts profound effects on neurodevelopment and increases the risk for several neurodevelopmental disorders including major depression. The mechanisms underlying the consequences of gestational stress are complex and remain to be elucidated. This study investigated the effects of gestational stress on depressive-like behavior and epigenetic modifications in young adult offspring. Gestational stress was induced by a combination of restraint and 24-hour light disturbance to pregnant dams throughout gestation. Depressive-like and anxiety-like behaviors of young adult offspring were examined. The expression and promoter methylation of brain derived neurotrophic factor (BDNF) were measured using RT-qPCR, Western blot, methylated DNA immunoprecipitation (MeDIP) and chromatin immunoprecipitation (ChIP). In addition, the expressions of histone deacetylases (HDACs) and acetylated histone H3 lysine 14 (AcH3K14) were also analyzed. Our results show that offspring from gestational stress dams exhibited depressive-like and anxiety-like behaviors. Biochemically, stress-offspring showed decreased expression of BDNF, increased expression of DNMT1, HDAC1, and HDAC2, and decreased expression of AcH3K14 in the hippocampus as compared to non-stress offspring. Data from MeDIP and ChIP assays revealed an increased methylation as well as decreased binding of AcH3K14 on specific BDNF promoters. Pearson analyses indicated that epigenetic changes induced by gestational stress were correlated with depressive-like and anxiety-like behaviors. These data suggest that gestational stress may be a suitable model for understanding the behavioral and molecular epigenetic changes observed in patients with depression.     

Polymerase chain reaction-based methods of DNA methylation analysis

DNA methylation is the main epigenetic modification in humans, and changes in methylation patterns play an important role in tumorigenesis. Hypermethylation of normally unmethylated CpG islands in the promoter regions often occurs in important tumor suppressor genes, DNA repair genes, and metastasis inhibitor genes. The changes of methylation status of various gene promoters seem to be a common feature of malignant cells and these changes can occur early in the progression process. Therefore detection of aberrant promoter hypermethylation of cancer-related genes may be useful for cancer diagnosis or detection of cancer recurrence. The purpose of this review is to provide a summary of the most commonly used techniques for the study of DNA methylation. Current scientific literature involving methylation detection methods was reviewed with an emphasis on polymerase chain reaction (PCR)-based detection methods. The current methodologies may be broadly classed into PCR-based methylation assays and non-PCR-based methylation assays. The problems and advantages of the different methods for detecting aberrant methylation are discussed. As the number of genes known to be hypermethylated in cancer is growing, the detection of aberrant promoter region methylation will be a promising approach for using DNA-based markers for the early detection of human cancers. Many techniques, especially PCR-based methylation assay techniques, make it practical to use these new methylation biomarkers in early cancer diagnosis.     

Analysis of complex methylation data

The development of multiple DNA methylation analysis techniques, including higher-throughput assays, has resulted in data structures of increasing complexity and diversity. Here, we discuss the general principles of DNA methylation analysis and propose a nomenclature for the various types of methylation analysis. We briefly outline several DNA methylation analysis techniques and discuss how these different technologies affect the structure of the resulting methylation data. We then describe the basic statistics and bioinformatic principles relevant to the analysis of simple and complex methylation data.     

Genome-wide methylation profiling identifies hypermethylated biomarkers in high-grade cervical intraepithelial neoplasia

Epigenetic modifications, such as aberrant DNA promoter methylation, are frequently observed in cervical cancer. Identification of hypermethylated regions allowing discrimination between normal cervical epithelium and high-grade cervical intraepithelial neoplasia (CIN2/3), or worse, may improve current cervical cancer population-based screening programs. In this study, the DNA methylome of high-grade CIN lesions was studied using genome-wide DNA methylation screening to identify potential biomarkers for early diagnosis of cervical neoplasia. Methylated DNA Immunoprecipitation (MeDIP) combined with DNA microarray was used to compare DNA methylation profiles of epithelial cells derived from high-grade CIN lesions with normal cervical epithelium. Hypermethylated differentially methylated regions (DMRs) were identified. Validation of nine selected DMRs using BSP and MSP in cervical tissue revealed methylation in 63.2–94.7% high-grade CIN and in 59.3–100% cervical carcinomas. QMSP for the two most significant high-grade CIN-specific methylation markers was conducted exploring test performance in a large series of cervical scrapings. Frequency and relative level of methylation were significantly different between normal and cancer samples. Clinical validation of both markers in cervical scrapings from patients with an abnormal cervical smear confirmed that frequency and relative level of methylation were related with increasing severity of the underlying CIN lesion and that ROC analysis was discriminative. These markers represent the COL25A1 and KATNAL2 and their observed increased methylation upon progression could intimate the regulatory role in carcinogenesis. In conclusion, our newly identified hypermethylated DMRs represent specific DNA methylation patterns in high-grade CIN lesions and are candidate biomarkers for early detection.     

Genome-wide methylation profiling identifies hypermethylated biomarkers in high-grade cervical intraepithelial neoplasia

Epigenetic modifications, such as aberrant DNA promoter methylation, are frequently observed in cervical cancer. Identification of hypermethylated regions allowing discrimination between normal cervical epithelium and high-grade cervical intraepithelial neoplasia (CIN2/3), or worse, may improve current cervical cancer population-based screening programs. In this study, the DNA methylome of high-grade CIN lesions was studied using genome-wide DNA methylation screening to identify potential biomarkers for early diagnosis of cervical neoplasia. Methylated DNA Immunoprecipitation (MeDIP) combined with DNA microarray was used to compare DNA methylation profiles of epithelial cells derived from high-grade CIN lesions with normal cervical epithelium. Hypermethylated differentially methylated regions (DMRs) were identified. Validation of nine selected DMRs using BSP and MSP in cervical tissue revealed methylation in 63.2–94.7% high-grade CIN and in 59.3–100% cervical carcinomas. QMSP for the two most significant high-grade CIN-specific methylation markers was conducted exploring test performance in a large series of cervical scrapings. Frequency and relative level of methylation were significantly different between normal and cancer samples. Clinical validation of both markers in cervical scrapings from patients with an abnormal cervical smear confirmed that frequency and relative level of methylation were related with increasing severity of the underlying CIN lesion and that ROC analysis was discriminative. These markers represent the COL25A1 and KATNAL2 and their observed increased methylation upon progression could intimate the regulatory role in carcinogenesis. In conclusion, our newly identified hypermethylated DMRs represent specific DNA methylation patterns in high-grade CIN lesions and are candidate biomarkers for early detection.     

Methylome analysis using MeDIP-seq with low DNA concentrations

DNA methylation is an epigenetic mark that has a crucial role in many biological processes. To understand the functional consequences of DNA methylation on phenotypic plasticity, a genome-wide analysis should be embraced. This in turn requires a technique that balances accuracy, genome coverage, resolution and cost, yet is low in DNA input in order to minimize the drain on precious samples. Methylated DNA immunoprecipitation-sequencing (MeDIP-seq) fulfils these criteria, combining MeDIP with massively parallel DNA sequencing. Here we report an improved protocol using 100-fold less genomic DNA than that commonly used. We show comparable results for specificity (>97%) and enrichment (>100-fold) over a wide range of DNA concentrations (5,000-50 ng) and demonstrate the utility of the protocol for the generation of methylomes from rare bone marrow cells using 160-300 ng of starting DNA. The protocol described here, i.e., DNA extraction to generation of MeDIP-seq library, can be completed within 3-5 d.