Genome-scale DNA methylation mapping of clinical samples at single-nucleotide resolution

Highly efficient derivation of ventricular cardiomyocytes from induced pluripotent stem cells with a distinct epigenetic signature

Cardiomyocytes derived from pluripotent stem cells can be applied in drug testing, disease modeling and cell-based therapy. However, without procardiogenic growth factors, the efficiency of cardiomyogenesis from pluripotent stem cells is usually low and the resulting cardiomyocyte population is heterogeneous. Here, we demonstrate that induced pluripotent stem cells (iPSCs) can be derived from murine ventricular myocytes (VMs), and consistent with other reports of iPSCs derived from various somatic cell types, VM-derived iPSCs (ViPSCs) exhibit a markedly higher propensity to spontaneously differentiate into beating cardiomyocytes as compared to genetically matched embryonic stem cells (ESCs) or iPSCs derived from tail-tip fibroblasts. Strikingly, the majority of ViPSC-derived cardiomyocytes display a ventricular phenotype. The enhanced ventricular myogenesis in ViPSCs is mediated via increased numbers of cardiovascular progenitors at early stages of differentiation. In order to investigate the mechanism of enhanced ventricular myogenesis from ViPSCs, we performed global gene expression and DNA methylation analysis, which revealed a distinct epigenetic signature that may be involved in specifying the VM fate in pluripotent stem cells.

Deep sequencing reveals distinct patterns of DNA methylation in prostate cancer

Beginning with precursor lesions, aberrant DNA methylation marks the entire spectrum of prostate cancer progression. We mapped the global DNA methylation patterns in select prostate tissues and cell lines using MethylPlex-next-generation sequencing (M-NGS). Hidden Markov model-based next-generation sequence analysis identified ∼68,000 methylated regions per sample. While global CpG island (CGI) methylation was not differential between benign adjacent and cancer samples, overall promoter CGI methylation significantly increased from ~12.6% in benign samples to 19.3% and 21.8% in localized and metastatic cancer tissues, respectively (P-value < 2 × 10(-16)). We found distinct patterns of promoter methylation around transcription start sites, where methylation occurred not only on the CGIs, but also on flanking regions and CGI sparse promoters. Among the 6691 methylated promoters in prostate tissues, 2481 differentially methylated regions (DMRs) are cancer-specific, including numerous novel DMRs. A novel cancer-specific DMR in the WFDC2 promoter showed frequent methylation in cancer (17/22 tissues, 6/6 cell lines), but not in the benign tissues (0/10) and normal PrEC cells. Integration of LNCaP DNA methylation and H3K4me3 data suggested an epigenetic mechanism for alternate transcription start site utilization, and these modifications segregated into distinct regions when present on the same promoter. Finally, we observed differences in repeat element methylation, particularly LINE-1, between ERG gene fusion-positive and -negative cancers, and we confirmed this observation using pyrosequencing on a tissue panel. This comprehensive methylome map will further our understanding of epigenetic regulation in prostate cancer progression.

De novo DNA methylation: a germ cell perspective

DNA methylation is a fundamentally important epigenetic modification of the mammalian genome that has widespread influences on gene expression. During germ-cell specification and maturation, epigenetic reprogramming occurs and the DNA methylation landscape is profoundly remodelled. Defects in this process have major consequences for embryonic development and are associated with several genetic disorders. In this review we report our current understanding of the molecular mechanisms associated with de novo DNA methylation in germ cells. We discuss recent discoveries connecting histone modifications, transcription and the DNA methylation machinery, and consider how these new findings could lead to a model for methylation establishment. Elucidating how DNA methylation marks are established in the germline has been a challenge for nearly 20 years, but represents a key step towards a full understanding of several biological processes including genomic imprinting, epigenetic reprogramming and the establishment of the pluripotent state in early embryos.

Targeted bisulfite sequencing by solution hybrid selection and massively parallel sequencing

We applied a solution hybrid selection approach to the enrichment of CpG islands (CGIs) and promoter sequences from the human genome for targeted high-throughput bisulfite sequencing. A single lane of Illumina sequences allowed accurate and quantitative analysis of ~1 million CpGs in more than 21 408 CGIs and more than 15 946 transcriptional regulatory regions. Of the CpGs analyzed, 77–84% fell on or near capture probe sequences; 69–75% fell within CGIs. More than 85% of capture probes successfully yielded quantitative DNA methylation information of targeted regions. Differentially methylated regions (DMRs) were identified in the 5′-end regulatory regions, as well as the intra- and intergenic regions, particularly in the X-chromosome among the three breast cancer cell lines analyzed. We chose 46 candidate loci (762 CpGs) for confirmation with PCR-based bisulfite sequencing and demonstrated excellent correlation between two data sets. Targeted bisulfite sequencing of three DNA methyltransferase (DNMT) knockout cell lines and the wild-type HCT116 colon cancer cell line revealed a significant decrease in CpG methylation for the DNMT1 knockout and DNMT1, 3B double knockout cell lines, but not in DNMT3B knockout cell line. We demonstrated the targeted bisulfite sequencing approach to be a powerful method to uncover novel aberrant methylation in the cancer epigenome. Since all targets were captured and sequenced as a pool through a series of single-tube reactions, this method can be easily scaled up to deal with a large number of samples.

Application of microdroplet PCR for large-scale targeted bisulfite sequencing

Cytosine methylation of DNA CpG dinucleotides in gene promoters is an epigenetic modification that regulates gene transcription. While many methods exist to interrogate methylation states, few current methods offer large-scale, targeted, single CpG resolution. We report an approach combining bisulfite treatment followed by microdroplet PCR with next-generation sequencing to assay the methylation state of 50 genes in the regions 1 kb upstream of and downstream from their transcription start sites. This method yielded 96% coverage of the targeted CpGs and demonstrated high correlation between CpG island (CGI) DNA methylation and transcriptional regulation. The method was scaled to interrogate the methylation status of 77,674 CpGs in the promoter regions of 2100 genes in primary CD4 T cells. The 2100 gene library yielded 97% coverage of all targeted CpGs and 99% of the target amplicons.

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.

Reprogramming factor expression initiates widespread targeted chromatin remodeling

Despite rapid progress in characterizing transcription factor-driven reprogramming of somatic cells to an induced pluripotent stem cell (iPSC) state, many mechanistic questions still remain. To gain insight into the earliest events in the reprogramming process, we systematically analyzed the transcriptional and epigenetic changes that occur during early factor induction after discrete numbers of divisions. We observed rapid, genome-wide changes in the euchromatic histone modification, H3K4me2, at more than a thousand loci including large subsets of pluripotency-related or developmentally regulated gene promoters and enhancers. In contrast, patterns of the repressive H3K27me3 modification remained largely unchanged except for focused depletion specifically at positions where H3K4 methylation is gained. These chromatin regulatory events precede transcriptional changes within the corresponding loci. Our data provide evidence for an early, organized, and population-wide epigenetic response to ectopic reprogramming factors that clarify the temporal order through which somatic identity is reset during reprogramming.

Tissue-specific distribution and dynamic changes of 5-hydroxymethylcytosine in mammalian genomes

Cytosine residues in the vertebrate genome are enzymatically modified to 5-methylcytosine, which participates in transcriptional repression of genes during development and disease progression. 5-Methylcytosine can be further enzymatically modified to 5-hydroxymethylcytosine by the TET family of methylcytosine dioxygenases. Analysis of 5-methylcytosine and 5-hydroxymethylcytosine is confounded, as these modifications are indistinguishable by traditional sequencing methods even when supplemented by bisulfite conversion. Here we demonstrate a simple enzymatic approach that involves cloning, identification, and quantification of 5-hydroxymethylcytosine in various CCGG loci within murine and human genomes. 5-Hydroxymethylcytosine was prevalent in human and murine brain and heart genomic DNAs at several regions. The cultured cell lines NIH3T3 and HeLa both displayed very low or undetectable amounts of 5-hydroxymethylcytosine at the examined loci. Interestingly, 5-hydroxymethylcytosine levels in mouse embryonic stem cell DNA first increased then slowly decreased upon differentiation to embryoid bodies, whereas 5-methylcytosine levels increased gradually over time. Finally, using a quantitative PCR approach, we established that a portion of VANGL1 and EGFR gene body methylation in human tissue DNA samples is indeed hydroxymethylation.

DNA methylation in cancer development, diagnosis and therapy--multiple opportunities for genotoxic agents to act as methylome disruptors or remediators

The role of DNA methylation and recently discovered hydroxymethylation in the function of the human epigenome is currently one of the hottest topics in the life sciences. Progress in this field of research has been further accelerated by the discovery that alterations in the methylome are not only associated with key functions of cells and organisms, such as development, differentiation and gene expression, but may underlie a number of human diseases, including cancer. This review describes both well established and more recent observations concerning alterations in the methylome, i.e. the global and local distribution of 5-methylcytosines, involved in its normal functions. Then, the changes in DNA methylation pattern seen in cancer cells are discussed in the context of their utilisation in cancer diagnostics and treatment. On this basis, comparisons are made between natural covalent DNA modification and that induced by genotoxic agents, chemical carcinogens and antitumour drugs as regards their impact on epigenetic mechanisms. The available data suggest that DNA damage by genotoxins can mimic epigenetic markers and in consequence disrupt the proper function of the epigenome. On the other hand, the same processes in cancer cells, e.g. DNA demethylation as a result of DNA methyltransferase blocking or the induction of DNA repair by DNA adducts, may restore the activity of hypermethylated anticancer genes. The observed multiple mechanisms by which genotoxic agents directly affect methylome function suggest that chemical carcinogens act primarily as epigenome disruptors, whereas mutations are secondary events that occur at later stages of cancer development when genome-protecting mechanisms have already been deregulated.

Preparation of reduced representation bisulfite sequencing libraries for genome-scale DNA methylation profiling

Genome-wide mapping of 5-methylcytosine is of broad interest to many fields of biology and medicine. A variety of methods have been developed, and several have recently been advanced to genome-wide scale using arrays and next-generation sequencing approaches. We have previously reported reduced representation bisulfite sequencing (RRBS), a bisulfite-based protocol that enriches CG-rich parts of the genome, thereby reducing the amount of sequencing required while capturing the majority of promoters and other relevant genomic regions. The approach provides single-nucleotide resolution, is highly sensitive and provides quantitative DNA methylation measurements. This protocol should enable any standard molecular biology laboratory to generate RRBS libraries of high quality. Briefly, purified genomic DNA is digested by the methylation-insensitive restriction enzyme MspI to generate short fragments that contain CpG dinucleotides at the ends. After end-repair, A-tailing and ligation to methylated Illumina adapters, the CpG-rich DNA fragments (40-220 bp) are size selected, subjected to bisulfite conversion, PCR amplified and end sequenced on an Illumina Genome Analyzer. Note that alignment and analysis of RRBS sequencing reads are not covered in this protocol. The extremely low input requirements (10-300 ng), the applicability of the protocol to formalin-fixed and paraffin-embedded samples, and the technique's single-nucleotide resolution extends RRBS to a wide range of biological and clinical samples and research applications. The entire process of RRBS library construction takes ∼9 d.

Reference Maps of human ES and iPS cell variation enable high-throughput characterization of pluripotent cell lines

The developmental potential of human pluripotent stem cells suggests that they can produce disease-relevant cell types for biomedical research. However, substantial variation has been reported among pluripotent cell lines, which could affect their utility and clinical safety. Such cell-line-specific differences must be better understood before one can confidently use embryonic stem (ES) or induced pluripotent stem (iPS) cells in translational research. Toward this goal we have established genome-wide reference maps of DNA methylation and gene expression for 20 previously derived human ES lines and 12 human iPS cell lines, and we have measured the in vitro differentiation propensity of these cell lines. This resource enabled us to assess the epigenetic and transcriptional similarity of ES and iPS cells and to predict the differentiation efficiency of individual cell lines. The combination of assays yields a scorecard for quick and comprehensive characterization of pluripotent cell lines.

Comparative methylome analysis of benign and malignant peripheral nerve sheath tumors

Aberrant DNA methylation (DNAm) was first linked to cancer over 25 yr ago. Since then, many studies have associated hypermethylation of tumor suppressor genes and hypomethylation of oncogenes to the tumorigenic process. However, most of these studies have been limited to the analysis of promoters and CpG islands (CGIs). Recently, new technologies for whole-genome DNAm (methylome) analysis have been developed, enabling unbiased analysis of cancer methylomes. By using MeDIP-seq, we report a sequencing-based comparative methylome analysis of malignant peripheral nerve sheath tumors (MPNSTs), benign neurofibromas, and normal Schwann cells. Analysis of these methylomes revealed a complex landscape of DNAm alterations. In contrast to what has been reported for other tumor types, no significant global hypomethylation was observed in MPNSTs using methylome analysis by MeDIP-seq. However, a highly significant (P < 10(-100)) directional difference in DNAm was found in satellite repeats, suggesting these repeats to be the main target for hypomethylation in MPNSTs. Comparative analysis of the MPNST and Schwann cell methylomes identified 101,466 cancer-associated differentially methylated regions (cDMRs). Analysis showed these cDMRs to be significantly enriched for two satellite repeat types (SATR1 and ARLα) and suggests an association between aberrant DNAm of these sequences and transition from healthy cells to malignant disease. Significant enrichment of hypermethylated cDMRs in CGI shores (P < 10(-60)), non-CGI-associated promoters (P < 10(-4)) and hypomethylated cDMRs in SINE repeats (P < 10(-100)) was also identified. Integration of DNAm and gene expression data showed that the expression pattern of genes associated with CGI shore cDMRs was able to discriminate between disease phenotypes. This study establishes MeDIP-seq as an effective method to analyze cancer methylomes.

Large scale loss of data in low-diversity illumina sequencing libraries can be recovered by deferred cluster calling

Massively parallel DNA sequencing is capable of sequencing tens of millions of DNA fragments at the same time. However, sequence bias in the initial cycles, which are used to determine the coordinates of individual clusters, causes a loss of fidelity in cluster identification on Illumina Genome Analysers. This can result in a significant reduction in the numbers of clusters that can be analysed. Such low sample diversity is an intrinsic problem of sequencing libraries that are generated by restriction enzyme digestion, such as e4C-seq or reduced-representation libraries. Similarly, this problem can also arise through the combined sequencing of barcoded, multiplexed libraries. We describe a procedure to defer the mapping of cluster coordinates until low-diversity sequences have been passed. This simple procedure can recover substantial amounts of next generation sequencing data that would otherwise be lost.

Epigenetic modifications in pluripotent and differentiated cells

Epigenetic modifications constitute a complex regulatory layer on top of the genome sequence. Pluripotent and differentiated cells provide a powerful system for investigating how the epigenetic code influences cellular fate. High-throughput sequencing of these cell types has yielded DNA methylation maps at single-nucleotide resolution and many genome-wide chromatin maps. In parallel to epigenome mapping efforts, remarkable progress has been made in our ability to manipulate cell states; ectopic expression of transcription factors has been shown to override developmentally established epigenetic marks and to enable routine generation of induced pluripotent stem (iPS) cells. Despite these advances, many fundamental questions remain. The roles of epigenetic marks and, in particular, of epigenetic modifiers in development and in disease states are not well understood. Although iPS cells appear molecularly and functionally similar to embryonic stem cells, more genome-wide studies are needed to define the extent and functions of epigenetic remodeling during reprogramming.

Tackling the epigenome: challenges and opportunities for collaboration

What are the key considerations to take into account when large-scale epigenomics projects are being implemented?

Epigenetic mechanisms in diabetic vascular complications

There has been a rapid increase in the incidence of diabetes as well the associated vascular complications. Both genetic and environmental factors have been implicated in these pathologies. Increasing evidence suggests that epigenetic factors play a key role in the complex interplay between genes and the environment. Actions of major pathological mediators of diabetes and its complications such as hyperglycaemia, oxidant stress, and inflammatory factors can lead to dysregulated epigenetic mechanisms that affect chromatin structure and gene expression. Furthermore, persistence of this altered state of the epigenome may be the underlying mechanism contributing to a 'metabolic memory' that results in chronic inflammation and vascular dysfunction in diabetes even after achieving glycaemic control. Further examination of epigenetic mechanisms by also taking advantage of recently developed next-generation sequencing technologies can provide novel insights into the pathology of diabetes and its complications and lead to the discovery of much needed new drug targets for these diseases. In this review, we highlight the role of epigenetics in diabetes and its vascular complications, and recent technological advances that have significantly accelerated the field.

Two groups of DNA methylation markers to classify colorectal cancer into three epigenotypes

A subgroup of colorectal cancer (CRC) shows non-random accumulation of aberrant DNA methylation, so-called CpG island methylator phenotype (CIMP), which was associated with microsatellite instability and BRAF mutation. As just one group of methylation markers was suitable to extract CIMP+/CIMP-high, and had been commonly used in the "one-panel method", it had been unclear whether another cluster of CRC with DNA methylation accumulation exists in microsatellite-stable CRC. We therefore epigenotyped CRC by a comprehensive approach, that is, the two-way unsupervised hierarchical clustering method using highly quantitative methylation data by a single detection method, MALDI-TOF mass spectrometry, on novel regions selected genome-widely through methylated DNA immunoprecipitation on array analysis. CRC was clearly clustered into three DNA methylation epigenotypes, high-, intermediate- and low-methylation epigenotypes (HME, IME, and LME, respectively). Methylation markers are clustered into two distinct groups: Group-1 methylated specifically in HME and including most reported CIMP-related markers; and Group-2 methylated both in HME and IME. While suitable markers to detect a subgroup of CRC with intermediate methylation and correlation to KRAS mutation have been expected to be developed, our data indicated that a "two-panel method" is necessary to properly classify CRC into three epigenotypes, the first panel to extract HME using Group-1 markers, and the second panel to divide the remaining into IME and LME using Group-2 markers. Here we review and compare our recent study and reported CRC classification methods by DNA methylation information, and propose the use of two panels of methylation markers as CRC classifiers.

Genome-wide distribution of DNA methylation at single-nucleotide resolution

DNA methylation, a well-known epigenetic modification in mammalian genomes, is important for development and health. Dysregulation of DNA methylation can cause abnormal gene regulation, leading to anomalous development and diseases. Until recently, the ability to understand the functions and dynamics of DNA methylation was limited by the availability of technologies for comprehensively characterizing methylation on a genome-wide scale. Rapid advances in high-throughput approaches (particularly next-generation sequencing), coupled with molecular techniques, have enabled unbiased genome-wide profiling of DNA modifications at single-base resolution and helped to elucidate their impact on gene regulation. Here, we discuss the development of genomic approaches to decipher the global methylome at single-base resolution, the challenges faced, and the emerging new insights. Our ability to decipher this important epigenetic modification and how it impacts gene expression will provide a framework for understanding numerous disease mechanisms, and suggest means to treat or prevent them in the future.

Evaluation of affinity-based genome-wide DNA methylation data: effects of CpG density, amplification bias, and copy number variation

DNA methylation is an essential epigenetic modification that plays a key role associated with the regulation of gene expression during differentiation, but in disease states such as cancer, the DNA methylation landscape is often deregulated. There are now numerous technologies available to interrogate the DNA methylation status of CpG sites in a targeted or genome-wide fashion, but each method, due to intrinsic biases, potentially interrogates different fractions of the genome. In this study, we compare the affinity-purification of methylated DNA between two popular genome-wide techniques, methylated DNA immunoprecipitation (MeDIP) and methyl-CpG binding domain-based capture (MBDCap), and show that each technique operates in a different domain of the CpG density landscape. We explored the effect of whole-genome amplification and illustrate that it can reduce sensitivity for detecting DNA methylation in GC-rich regions of the genome. By using MBDCap, we compare and contrast microarray- and sequencing-based readouts and highlight the impact that copy number variation (CNV) can make in differential comparisons of methylomes. These studies reveal that the analysis of DNA methylation data and genome coverage is highly dependent on the method employed, and consideration must be made in light of the GC content, the extent of DNA amplification, and the copy number.

Aberrant global methylation patterns affect the molecular pathogenesis and prognosis of multiple myeloma

We used genome-wide methylation microarrays to analyze differences in CpG methylation patterns in cells relevant to the pathogenesis of myeloma plasma cells (B cells, normal plasma cells, monoclonal gammopathy of undetermined significance [MGUS], presentation myeloma, and plasma cell leukemia). We show that methylation patterns in these cell types are capable of distinguishing nonmalignant from malignant cells and the main reason for this difference is hypomethylation of the genome at the transition from MGUS to presentation myeloma. In addition, gene-specific hypermethylation was evident at the myeloma stage. Differential methylation was also evident at the transition from myeloma to plasma cell leukemia with remethylation of the genome, particularly of genes involved in cell-cell signaling and cell adhesion, which may contribute to independence from the bone marrow microenvironment. There was a high degree of methylation variability within presentation myeloma samples, which was associated with cytogenetic differences between samples. More specifically, we found methylation subgroups were defined by translocations and hyperdiploidy, with t(4;14) myeloma having the greatest impact on DNA methylation. Two groups of hyperdiploid samples were identified, on the basis of unsupervised clustering, which had an impact on overall survival. Overall, DNA methylation changes significantly during disease progression and between cytogenetic subgroups.

Thirteenth AACR Special Conference on Cancer Epigenetics

The 13th American Association of Cancer Research (AACR) Special Conference on Cancer Epigenetics was held in San Juan, Puerto Rico, on January 20–23, 2010. This event heralded insights arising from the large sequencing efforts within the NIH Epigenetics Roadmap Project and the Cancer Genome Atlas Project, as well as important new discoveries in the basic biology of epigenetics and results of epigenetic drug clinical trials. A summary of the recent findings is presented here, with particular emphasis on the major themes of the conference. Cancer Res; 70(19); 7372–8. ©2010 AACR.

Genome-wide DNA methylation maps in follicular lymphoma cells determined by methylation-enriched bisulfite sequencing

Background

Follicular lymphoma (FL) is a form of non-Hodgkin's lymphoma (NHL) that arises from germinal center (GC) B-cells. Despite the significant advances in immunotherapy, FL is still not curable. Beyond transcriptional profiling and genomics datasets, there currently is no epigenome-scale dataset or integrative biology approach that can adequately model this disease and therefore identify novel mechanisms and targets for successful prevention and treatment of FL.

Methodology/Principal Findings

We performed methylation-enriched genome-wide bisulfite sequencing of FL cells and normal CD19⁺ B-cells using 454 sequencing technology. The methylated DNA fragments were enriched with methyl-binding proteins, treated with bisulfite, and sequenced using the Roche-454 GS FLX sequencer. The total number of bases covered in the human genome was 18.2 and 49.3 million including 726,003 and 1.3 million CpGs in FL and CD19⁺ B-cells, respectively. 11,971 and 7,882 methylated regions of interest (MRIs) were identified respectively. The genome-wide distribution of these MRIs displayed significant differences between FL and normal B-cells. A reverse trend in the distribution of MRIs between the promoter and the gene body was observed in FL and CD19⁺ B-cells. The MRIs identified in FL cells also correlated well with transcriptomic data and ChIP-on-Chip analyses of genome-wide histone modifications such as tri-methyl-H3K27, and tri-methyl-H3K4, indicating a concerted epigenetic alteration in FL cells.

Conclusions/Significance

This study is the first to provide a large scale and comprehensive analysis of the DNA methylation sequence composition and distribution in the FL epigenome. These integrated approaches have led to the discovery of novel and frequent targets of aberrant epigenetic alterations. The genome-wide bisulfite sequencing approach developed here can be a useful tool for profiling DNA methylation in clinical samples.

Statistical Issues in the Analysis of ChIP-Seq and RNA-Seq Data

The recent arrival of ultra-high throughput, next generation sequencing (NGS) technologies has revolutionized the genetics and genomics fields by allowing rapid and inexpensive sequencing of billions of bases. The rapid deployment of NGS in a variety of sequencing-based experiments has resulted in fast accumulation of massive amounts of sequencing data. To process this new type of data, a torrent of increasingly sophisticated algorithms and software tools are emerging to help the analysis stage of the NGS applications. In this article, we strive to comprehensively identify the critical challenges that arise from all stages of NGS data analysis and provide an objective overview of what has been achieved in existing works. At the same time, we highlight selected areas that need much further research to improve our current capabilities to delineate the most information possible from NGS data. The article focuses on applications dealing with ChIP-Seq and RNA-Seq.

Bisulfite sequencing of DNA

Exact positions of 5-methylcytosine (m(5)C) on a single strand of DNA can be determined by bisulfite genomic sequencing (BGS). Treatment with bisulfite ion preferentially deaminates unmethylated cytosines, which are then converted to uracil upon desulfonation. Amplifying regions of interest from deaminated DNA and sequencing products cloned from amplicons permits determination of methylation at single-nucleotide resolution along single DNA molecules, which is not possible with other methylation analysis techniques. This unit describes a BGS technique suitable for most DNA sources, including formaldehyde-fixed tissue. Considerations for experimental design and common sources of error are discussed.

Comparison of sequencing-based methods to profile DNA methylation and identification of monoallelic epigenetic modifications

Analysis of DNA methylation patterns relies increasingly on sequencing-based profiling methods. The four most frequently used sequencing-based technologies are the bisulfite-based methods MethylC-seq and reduced representation bisulfite sequencing (RRBS), and the enrichment-based techniques methylated DNA immunoprecipitation sequencing (MeDIP-seq) and methylated DNA binding domain sequencing (MBD-seq). We applied all four methods to biological replicates of human embryonic stem cells to assess their genome-wide CpG coverage, resolution, cost, concordance and the influence of CpG density and genomic context. The methylation levels assessed by the two bisulfite methods were concordant (their difference did not exceed a given threshold) for 82% for CpGs and 99% of the non-CpG cytosines. Using binary methylation calls, the two enrichment methods were 99% concordant and regions assessed by all four methods were 97% concordant. We combined MeDIP-seq with methylation-sensitive restriction enzyme (MRE-seq) sequencing for comprehensive methylome coverage at lower cost. This, along with RNA-seq and ChIP-seq of the ES cells enabled us to detect regions with allele-specific epigenetic states, identifying most known imprinted regions and new loci with monoallelic epigenetic marks and monoallelic expression.

Quantitative comparison of genome-wide DNA methylation mapping technologies

DNA methylation is a key component of mammalian gene regulation and the most classical example of an epigenetic mark. DNA methylation patterns are mitotically heritable and stable over time, but they undergo considerable changes in response to cell differentiation, diseases and environmental influences. Several methods have been developed for DNA methylation profiling on a genomic scale. Here, we benchmark four of these methods on two sample pairs, comparing their accuracy and power to detect DNA methylation differences. The results show that all evaluated methods (MeDIP-seq: methylated DNA immunoprecipitation, MethylCap-seq: methylated DNA capture by affinity purification, RRBS: reduced representation bisulfite sequencing, and the Infinium HumanMethylation27 assay) produce accurate DNA methylation data. However, these methods differ in their ability to detect differentially methylated regions between pairs of samples. We highlight strengths and weaknesses of the four methods and give practical recommendations for the design of epigenomic case-control studies.

MetMap enables genome-scale Methyltyping for determining methylation states in populations

The ability to assay genome-scale methylation patterns using high-throughput sequencing makes it possible to carry out association studies to determine the relationship between epigenetic variation and phenotype. While bisulfite sequencing can determine a methylome at high resolution, cost inhibits its use in comparative and population studies. MethylSeq, based on sequencing of fragment ends produced by a methylation-sensitive restriction enzyme, is a method for methyltyping (survey of methylation states) and is a site-specific and cost-effective alternative to whole-genome bisulfite sequencing. Despite its advantages, the use of MethylSeq has been restricted by biases in MethylSeq data that complicate the determination of methyltypes. Here we introduce a statistical method, MetMap, that produces corrected site-specific methylation states from MethylSeq experiments and annotates unmethylated islands across the genome. MetMap integrates genome sequence information with experimental data, in a statistically sound and cohesive Bayesian Network. It infers the extent of methylation at individual CGs and across regions, and serves as a framework for comparative methylation analysis within and among species. We validated MetMap's inferences with direct bisulfite sequencing, showing that the methylation status of sites and islands is accurately inferred. We used MetMap to analyze MethylSeq data from four human neutrophil samples, identifying novel, highly unmethylated islands that are invisible to sequence-based annotation strategies. The combination of MethylSeq and MetMap is a powerful and cost-effective tool for determining genome-scale methyltypes suitable for comparative and association studies.

Next generation sequencing: advances in characterizing the methylome

Epigenetic modifications play an important role in lymphoid malignancies. This has been evidenced by the large body of work published using microarray technologies to generate methylation profiles for numerous types and subtypes of lymphoma and leukemia. These studies have shown the importance of defining the epigenome so that we can better understand the biology of lymphoma. Recent advances in DNA sequencing technology have transformed the landscape of epigenomic analysis as we now have the ability to characterize the genome-wide distribution of chromatin modifications and DNA methylation using next-generation sequencing. To take full advantage of the throughput of next-generation sequencing, there are many methodologies that have been developed and many more that are currently being developed. Choosing the appropriate methodology is fundamental to the outcome of next-generation sequencing studies. In this review, published technologies and methodologies applicable to studying the methylome are presented. In addition, progress towards defining the methylome in lymphoma is discussed and prospective directions that have been made possible as a result of next-generation sequencing technology. Finally, methodologies are introduced that have not yet been published but that are being explored in the pursuit of defining the lymphoma methylome.

Whole-genome DNA methylation profiling using MethylCap-seq

MethylCap-seq is a robust procedure for genome-wide profiling of DNA methylation. The approach consists of the capture of methylated DNA using the MBD domain of MeCP2, and subsequent next-generation sequencing of eluted DNA. Elution of the captured methylated DNA is done in steps using a salt gradient, which stratifies the genome into fractions with different CpG density. The enrichment reached within the individual eluates allows for cost-effective deep sequence coverage. The profiles together yield a detailed genome-wide map of methylated regions and readily allows detection of DNA methylation in known and novel regions. Here, we describe principles and details of the MethylCap-seq procedure using different sources of starting material.

The application of next generation sequencing in DNA methylation analysis

DNA methylation is a major form of epigenetic modification and plays essential roles in physiology and disease processes. In the human genome, about 80% of cytosines in the 56 million CpG sites are methylated to 5-methylcytosines. The methylation pattern of DNA is highly variable among cells types and developmental stages and influenced by disease processes and genetic factors, which brings considerable theoretical and technological challenges for its comprehensive mapping. Recently various high-throughput approaches based on bisulfite conversion combined with next generation sequencing have been developed and applied for the genome wide analysis of DNA methylation. These methods provide single base pair resolution, quantitative DNA methylation data with genome wide coverage. We review these methods here and discuss some technical points of special interest like the sequence depth necessary to reach conclusions, the identification of clonal DNA amplification after bisulfite conversion and the detection of non-CpG methylation. Future application of these methods will greatly facilitate the profiling of the DNA methylation in the genomes of different species, individuals and cell types under healthy and disease states.

Introduction into the analysis of high-throughput-sequencing based epigenome data

Sequencing-based approaches now allow high-resolution, genome-scale investigation of cellular epigenetic landscapes. For example, mapping of open chromatin regions, post-translational histone modifications and DNA methylation across a whole genome is now feasible, and new non-coding regulatory RNAs can be sensitively identified via RNA sequencing. The resulting large-scale data sets promise to contribute towards a more precise and complete understanding of gene regulation and to yield insights into the interplay between genomes and the environment. In this article, I review some of the conceptual issues and currently available software tools for the analysis of sequencing-based whole-genome epigenetics data.

Identification and characterization of putative methylation targets in the MAOA locus using bioinformatic approaches

Monoamine oxidase A (MAO A) is an enzyme that catalyzes the oxidation of neurotransmitter amines. A functional polymorphism in the human MAOA gene (high- and low-MAOA) has been associated with distinct behavioral phenotypes. To investigate directly the biological mechanism whereby this polymorphism influences brain function, we recently measured the activity of the MAO A enzyme in healthy volunteers. When found no relationship between the individual's brain MAO A level and the MAOA genotype, we postulated that there are additional regulatory mechanisms that control the MAOA expression. Given that DNA methylation is linked to the regulation of gene expression, we hypothesized that epigenetic mechanisms factor into the MAOA expression. Our underplaying assumption was that the differences in an individual's genotype play a key role in the epigenetic potential of the MAOA locus and, consequently, determine the individual's level of MAO A activity in the brain. As a first step towards experimental validation of the hypothesis, we performed a comprehensive bioinformatic analysis aiming to interrogate genomic features and attributes of the MAOA locus that might modulate its epigenetic sensitivity. Major findings of our analysis are the following: (1) the extended MAOA regulatory region contains two CpG islands (CGIs), one of which overlaps with the canonical MAOA promoter and the other is located further upstream; both CGIs exhibit sensitivity to differential methylation. (2) The uVNTR's effect on the MAOA's transcriptional activity might have epigenetic nature: this polymorphic region resides within the MAOA's CGI and itself contains CpGs, thus, the number of repeating increments effectively changes the number of methylatable cytosines in the MAOA promoter. An array of in silico analyses (the nucleosome positioning, the physical properties of the local DNA, the clustering of transcription-factor binding sites) together with experimental data on histone modifications and Pol 2 sites and data from the RefSeq mRNA library suggest that the MAOA gene might have an alternative promoter. Based on our findings, we propose a regulatory mechanism for the human MAOA according to which the MAOA expression in vivo is executed by the generation of tissue-specific transcripts initiated from the alternative promoters (both CGI-associated) where transcriptional activation of a particular promoter is under epigenetic control.

Identification of driver and passenger DNA methylation in cancer by epigenomic analysis

Human cancer genomes are characterized by widespread aberrations in DNA methylation patterns including DNA hypomethylation of mostly repetitive sequences and hypermethylation of numerous CpG islands. The analysis of DNA methylation patterns in cancer has progressed from single gene studies examining potentially important candidate genes to a more global analysis where all or almost all promoter and CpG island sequences can be analyzed. We provide a brief overview of these genome-scale methylation-profiling techniques, summarize some of the information that has been obtained with these approaches, and discuss what we have learned about the specificity of methylation aberrations in cancer at a genome-wide level. The challenge is now to identify those methylation changes that are thought to be crucial for the processes of tumor initiation, tumor progression, or metastasis and distinguish these from methylation changes that are merely passenger events that accompany the transformation process but have no effect per se on the process of carcinogenesis.

Epigenetic regulation of pluripotency

Epigenetic regulation refers to the mechanisms that alter gene expression patterns in the absence of changes in the nucleotide sequence of the DNA molecule. The best understood epigenetic marks include posttranslational modifications of the histone tails and DNA methylation. Both play central roles in normal development and in diseases. Pluripotent stem cells have great promise for regenerative medicine and recent efforts have focused on identifying molecular networks that govern pluripotency. This chapter provides an overview of epigenetic regulation in embryonic stem cells. We present a brief introduction into epigenetic mechanisms and focus on their role in pluripotent cells.

An Empirical Bayes Approach for Methylation Differentiation at the Single Nucleotide Resolution

DNA methylation is an important epigenetic phenomenon that is associated with a variety of diseases, particularly cancers. Recent development of high throughput sequencing technology has enabled researchers to investigate the methylation rate at a single nucleotide resolution for any given sample. Testing for methylation rate equality or difference between two samples, however, is challenged by the small sample size observed at many sites across the genome. Fisher's exact test is typically used in this situation; however, it is conservative and it cannot be used to test for specific difference in methylation rate between two samples. In this paper, we propose an empirical Bayes approach that utilizes the genome-wide data as prior information for methylation differentiation between two samples. We show that this new approach is more powerful than Fisher's exact test. In addition, it can be used to test for any specific methylation difference while controlling the false discovery rate (FDR). The new method is applied to a real data set from a colon tumor study.