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

MeQA: a pipeline for MeDIP-seq data quality assessment and analysis

MOTIVATION

We present a pipeline for the pre-processing, quality assessment, read distribution and methylation estimation for methylated DNA immunoprecipitation (MeDIP)-sequence datasets. This is the first MeDIP-seq-specific analytic pipeline that starts at the output of the sequencers. This pipeline will reduce the data analysis load on staff and allows the easy and straightforward analysis of sequencing data for DNA methylation. The pipeline integrates customized scripting and several existing tools, which can deal with both paired and single end data.

AVAILABILITY

The package and extensive documentation, and comparison to public data is available at http://life.tongji.edu.cn/meqa/.

Detecting Chemically Modified DNA Bases Using Surface Enhanced Raman Spectroscopy

Post-translational modifications of DNA- changes in the chemical structure of individual bases that occur without changes in the DNA sequence- are known to alter gene expression. They are believed to result in frequently deleterious phenotypic changes, such as cancer. Methylation of adenine, methylation and hydroxymethylation of cytosine, and guanine oxidation are the primary DNA base modifications identified to date. Here we show it is possible to use surface enhanced Raman spectroscopy (SERS) to detect these primary DNA base modifications. SERS detection of modified DNA bases is label-free and requires minimal additional sample preparation, reducing the possibility of additional chemical modifications induced prior to measurement. This approach shows the feasibility of DNA base modification assessment as a potentially routine analysis that may be further developed for clinical diagnostics.

Bio-CAP: a versatile and highly sensitive technique to purify and characterise regions of non-methylated DNA

Across vertebrate genomes methylation of cytosine residues within the context of CpG dinucleotides is a pervasive epigenetic mark that can impact gene expression and has been implicated in various developmental and disease-associated processes. Several biochemical approaches exist to profile DNA methylation, but recently an alternative approach based on profiling non-methylated CpGs was developed. This technique, called CxxC affinity purification (CAP), uses a ZF-CxxC (CxxC) domain to specifically capture DNA containing clusters of non-methylated CpGs. Here we describe a new CAP approach, called biotinylated CAP (Bio-CAP), which eliminates the requirement for specialized equipment while dramatically improving and simplifying the CxxC-based DNA affinity purification. Importantly, this approach isolates non-methylated DNA in a manner that is directly proportional to the density of non-methylated CpGs, and discriminates non-methylated CpGs from both methylated and hydroxymethylated CpGs. Unlike conventional CAP, Bio-CAP can be applied to nanogram quantities of genomic DNA and in a magnetic format is amenable to efficient parallel processing of samples. Furthermore, Bio-CAP can be applied to genome-wide profiling of non-methylated DNA with relatively small amounts of input material. Therefore, Bio-CAP is a simple and streamlined approach for characterizing regions of the non-methylated DNA, whether at specific target regions or genome wide.

High resolution profiling of human exon methylation by liquid hybridization capture-based bisulfite sequencing

Background

DNA methylation plays important roles in gene regulation during both normal developmental and disease states. In the past decade, a number of methods have been developed and applied to characterize the genome-wide distribution of DNA methylation. Most of these methods endeavored to screen whole genome and turned to be enormously costly and time consuming for studies of the complex mammalian genome. Thus, they are not practical for researchers to study multiple clinical samples in biomarker research.

Results

Here, we display a novel strategy that relies on the selective capture of target regions by liquid hybridization followed by bisulfite conversion and deep sequencing, which is referred to as liquid hybridization capture-based bisulfite sequencing (LHC-BS). To estimate this method, we utilized about 2 μg of native genomic DNA from YanHuang (YH) whole blood samples and a mature dendritic cell (mDC) line, respectively, to evaluate their methylation statuses of target regions of exome. The results indicated that the LHC-BS system was able to cover more than 97% of the exome regions and detect their methylation statuses with acceptable allele dropouts. Most of the regions that couldn't provide accurate methylation information were distributed in chromosomes 6 and Y because of multiple mapping to those regions. The accuracy of this strategy was evaluated by pair-wise comparisons using the results from whole genome bisulfite sequencing and validated by bisulfite specific PCR sequencing.

Conclusions

In the present study, we employed a liquid hybridisation capture system to enrich for exon regions and then combined with bisulfite sequencing to examine the methylation statuses for the first time. This technique is highly sensitive and flexible and can be applied to identify differentially methylated regions (DMRs) at specific genomic locations of interest, such as regulatory elements or promoters.

RRBSMAP: a fast, accurate and user-friendly alignment tool for reduced representation bisulfite sequencing

SUMMARY

Reduced representation bisulfite sequencing (RRBS) is a powerful yet cost-efficient method for studying DNA methylation on a genomic scale. RRBS involves restriction-enzyme digestion, bisulfite conversion and size selection, resulting in DNA sequencing data that require special bioinformatic handling. Here, we describe RRBSMAP, a short-read alignment tool that is designed for handling RRBS data in a user-friendly and scalable way. RRBSMAP uses wildcard alignment, and avoids the need for any preprocessing or post-processing steps. We benchmarked RRBSMAP against a well-validated MAQ-based pipeline for RRBS read alignment and observed similar accuracy but much improved runtime performance, easier handling and better scaling to large sample sets. In summary, RRBSMAP removes bioinformatic hurdles and reduces the computational burden of large-scale epigenome association studies performed with RRBS.

AVAILABILITY

http://rrbsmap.computational-epigenetics.org/ http://code.google.com/p/bsmap/

CONTACT

wl1@bcm.tmc.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genomic distribution and inter-sample variation of non-CpG methylation across human cell types

DNA methylation plays an important role in development and disease. The primary sites of DNA methylation in vertebrates are cytosines in the CpG dinucleotide context, which account for roughly three quarters of the total DNA methylation content in human and mouse cells. While the genomic distribution, inter-individual stability, and functional role of CpG methylation are reasonably well understood, little is known about DNA methylation targeting CpA, CpT, and CpC (non-CpG) dinucleotides. Here we report a comprehensive analysis of non-CpG methylation in 76 genome-scale DNA methylation maps across pluripotent and differentiated human cell types. We confirm non-CpG methylation to be predominantly present in pluripotent cell types and observe a decrease upon differentiation and near complete absence in various somatic cell types. Although no function has been assigned to it in pluripotency, our data highlight that non-CpG methylation patterns reappear upon iPS cell reprogramming. Intriguingly, the patterns are highly variable and show little conservation between different pluripotent cell lines. We find a strong correlation of non-CpG methylation and DNMT3 expression levels while showing statistical independence of non-CpG methylation from pluripotency associated gene expression. In line with these findings, we show that knockdown of DNMTA and DNMT3B in hESCs results in a global reduction of non-CpG methylation. Finally, non-CpG methylation appears to be spatially correlated with CpG methylation. In summary these results contribute further to our understanding of cytosine methylation patterns in human cells using a large representative sample set.

A concise review on epigenetic regulation: insight into molecular mechanisms

Epigenetic mechanisms are responsible for the regulation of transcription of imprinted genes and those that induce a totipotent state. Starting just after fertilization, DNA methylation pattern undergoes establishment, reestablishment and maintenance. These modifications are important for normal embryo and placental developments. Throughout life and passing to the next generation, epigenetic events establish, maintain, erase and reestablish. In the context of differentiated cell reprogramming, demethylation and activation of genes whose expressions contribute to the pluripotent state is the crux of the matter. In this review, firstly, regulatory epigenetic mechanisms related to somatic cell nuclear transfer (SCNT) reprogramming are discussed, followed by embryonic development, and placental epigenetic issues.

Physiological regulation of transgene expression by a lentiviral vector containing the A2UCOE linked to a myeloid promoter

Protection against epigenetic silencing is a desirable feature of future gene therapy vectors, in particular for those applications in which transgene expression will not confer growth advantage to gene-transduced cells. The ubiquitous chromatin opening element (UCOE) consisting of the methylation-free CpG island encompassing the dual divergently transcribed promoters of the human HNRPA2B1-CBX3 housekeeping genes (A2UCOE) has been shown to shield constitutive active heterologous promoters from epigenetic modifications and chromosomal position effects. However, it is unclear if this element can be used to improve expression from tissue-specific enhancer/promoters, while maintaining tissue specificity in hematopoietic cells. Here, we evaluated the potential of the A2UCOE in combination with the myeloid-specific myeloid related protein 8 (MRP8) promoter to target transgene expression specifically to myeloid cells in vitro and in vivo from a self-inactivating lentiviral vector. The inclusion of the A2UCOE did not interfere with specific upregulation of MRP8 promoter activity during myeloid differentiation and mediated sustained and vector copy-dependent expression in myeloid cells. Notably, the A2UCOE did not protect the MRP8 promoter from methylation in the P19 embryonal carcinoma cell line, suggesting that this element maintains the inherent epigenetic state and transcriptional activity of cellular promoters in their native configuration. Thus, the A2UCOE could represent a useful protective genetic element in gene therapy vectors, ensuring physiological transcriptional regulation of tissue-specific promoters independent of the chromosomal integration site.

Epigenetics in breast cancer: what's new?

Epigenetic changes are critical for development and progression of cancers, including breast cancer. Significant progress has been made in the basic understanding of how various epigenetic changes such as DNA methylation, histone modification, miRNA expression, and higher order chromatin structure affect gene expression. The present review will focus on methylation and demethylation of histones. While the acetylation of histones has been at the forefront of well-characterized post-translational modifications of histones, including the development of inhibitors targeting de-acetylating enzymes, the past few years have witnessed a dramatic increase in knowledge regarding the role of histone methylation/demethylation. This is an exciting and rapidly evolving area of research, with much promise for potential clinical intervention in several cancers including breast cancer. We also summarize efforts to identity DNA methylation signatures that could be prognostic and/or predictive markers in breast cancer, focusing on recent studies using genome-wide approaches. Finally, we briefly review the efforts made by both the National Institutes of Health Epigenome Project and The Cancer Genome Atlas, especially highlighting the study of breast cancer epigenetics, exciting technological advances, potential roadblocks, and future directions.

Epigenetic-based companion diagnostics

In current medical practice, when a patient is diagnosed with cancer the treating physician generally follows a standard protocol, assigning the treatment that gives a favorable response in the largest proportion of patients. However, in many individual instances this approach may not be the most effective solution and, typically, treatment is only initiated or altered once the cancer has actually started progressing. During this process, patients will lose treatment time waiting to start chemotherapy or will endure severe side effects associated with toxic chemotherapeutic treatments. While some patients are undertreated because current diagnostic methods cannot provide accurate enough information regarding the aggressiveness or drug response of their disease, others with nonaggressive forms of cancer are overtreated and unnecessarily undergo the side effects associated with chemotherapeutic treatment. Epigenetic markers have been widely investigated and are considered key regulators of cellular transcription. Histone modifications and DNA methylation have been demonstrated to play key roles in maintaining stem-cell-like states, cellular differentiation and cancer. In particular, DNA methylation is a frequent, abundant and stable cancer mark, with an inherent role in oncogenesis and tumor progression. In this article, the potential of DNA methylation as a companion diagnostic is assessed, illustrated by exploring some development paths. Epigenetic silencing of MGMT is a key example of how biomarker development, biological pathways and clinical utility come together, serving as a hallmark of epigenetic companion diagnostics.

A new frontier in personalized cancer therapy: mapping molecular changes

Mutations in the genome of a normal cell can affect the function of its many genes and pathways. These alterations could eventually transform the cell from a normal to a malignant state by allowing an uncontrolled proliferation of the cell and formation of a cancer tumor. Each tumor in an individual patient can have hundreds of mutated genes and perturbed pathways. Cancers clinically presenting as the same type or subtype could potentially be very different at the molecular level and thus behave differently in response to therapy. The challenge is to distinguish the key mutations driving the cancer from the background of mutational noise and find ways to effectively target them. The promise is that such a molecular approach to classifying cancer will lead to better diagnostic, prognostic and personalized treatment strategies. This article provides an overview of advances in the molecular characterization of cancers and their applications in therapy.

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.

DNA methylation: a timeline of methods and applications

DNA methylation is a biochemical process where a DNA base, usually cytosine, is enzymatically methylated at the 5-carbon position. An epigenetic modification associated with gene regulation, DNA methylation is of paramount importance to biological health and disease. Recently, the quest to unravel the Human Epigenome commenced, calling for a modernization of previous DNA methylation profiling techniques. Here, we describe the major developments in the methodologies used over the past three decades to examine the elusive epigenome (or methylome). The earliest techniques were based on the separation of methylated and unmethylated cytosines via chromatography. The following years would see molecular techniques being employed to indirectly examine DNA methylation levels at both a genome-wide and locus-specific context, notably immunoprecipitation via anti-5'methylcytosine and selective digestion with methylation-sensitive restriction endonucleases. With the advent of sodium bisulfite treatment of DNA, a deamination reaction that converts cytosine to uracil only when unmethylated, the epigenetic modification can now be identified in the same manner as a DNA base-pair change. More recently, these three techniques have been applied to more technically advanced systems such as DNA microarrays and next-generation sequencing platforms, bringing us closer to unveiling a complete human epigenetic profile.

The state of genome-wide association studies in pulmonary disease: a new perspective

With rapid advances in our knowledge of the human genome and increasing availability of high-throughput investigative technology, genome-wide association (GWA) studies have recently gained marked popularity. As an unbiased approach to identifying genomic regions of importance in complex human disease, the results of such studies have the potential to illuminate novel causal pathways, guide mechanistic research, and aid in prediction of disease risk. The use of a genome-wide approach presents considerable methodological and statistical challenges, and properly conducted studies are essential to avoid false-positive results. A total of 22 GWA studies have been published in pulmonary medicine thus far, implicating several intriguing genomic regions in the determination of pulmonary function measures, onset of asthma, and susceptibility to chronic obstructive pulmonary disease. Many questions remain, however, as most identified genetic variants contribute only nominally to overall disease risk, genetic disease mechanisms remain uncertain, and disease-associated variants are not consistent across studies. Perhaps most fundamentally, the association signals identified have not yet been traced to causal variants. This perspective will review the current state of GWA studies in pulmonary disease. We begin with an introduction to the hypothesis, principles, and limitations of this type of genome-wide approach, highlight key points from available studies, and conclude by addressing future approaches to better understand the genetics of complex pulmonary disease.

Unique DNA methylome profiles in CpG island methylator phenotype colon cancers

A subset of colorectal cancers was postulated to have the CpG island methylator phenotype (CIMP), a higher propensity for CpG island DNA methylation. The validity of CIMP, its molecular basis, and its prognostic value remain highly controversial. Using MBD-isolated genome sequencing, we mapped and compared genome-wide DNA methylation profiles of normal, non-CIMP, and CIMP colon specimens. Multidimensional scaling analysis revealed that each specimen could be clearly classified as normal, non-CIMP, and CIMP, thus signifying that these three groups have distinctly different global methylation patterns. We discovered 3780 sites in various genomic contexts that were hypermethylated in both non-CIMP and CIMP colon cancers when compared with normal colon. An additional 2026 sites were found to be hypermethylated in CIMP tumors only; and importantly, 80% of these sites were located in CpG islands. These data demonstrate on a genome-wide level that the additional hypermethylation seen in CIMP tumors occurs almost exclusively at CpG islands and support definitively that these tumors were appropriately named. When these sites were examined more closely, we found that 25% were adjacent to sites that were also hypermethylated in non-CIMP tumors. Thus, CIMP is also characterized by more extensive methylation of sites that are already prone to be hypermethylated in colon cancer. These observations indicate that CIMP tumors have specific defects in controlling both DNA methylation seeding and spreading and serve as an important first step in delineating molecular mechanisms that control these processes.

Retrotransposon-induced heterochromatin spreading in the mouse revealed by insertional polymorphisms

The "arms race" relationship between transposable elements (TEs) and their host has promoted a series of epigenetic silencing mechanisms directed against TEs. Retrotransposons, a class of TEs, are often located in repressed regions and are thought to induce heterochromatin formation and spreading. However, direct evidence for TE-induced local heterochromatin in mammals is surprisingly scarce. To examine this phenomenon, we chose two mouse embryonic stem (ES) cell lines that possess insertionally polymorphic retrotransposons (IAP, ETn/MusD, and LINE elements) at specific loci in one cell line but not the other. Employing ChIP-seq data for these cell lines, we show that IAP elements robustly induce H3K9me3 and H4K20me3 marks in flanking genomic DNA. In contrast, such heterochromatin is not induced by LINE copies and only by a minority of polymorphic ETn/MusD copies. DNA methylation is independent of the presence of IAP copies, since it is present in flanking regions of both full and empty sites. Finally, such spreading into genes appears to be rare, since the transcriptional start sites of very few genes are less than one Kb from an IAP. However, the B3galtl gene is subject to transcriptional silencing via IAP-induced heterochromatin. Hence, although rare, IAP-induced local heterochromatin spreading into nearby genes may influence expression and, in turn, host fitness.

Genome-wide analysis distinguishes hyperglycemia regulated epigenetic signatures of primary vascular cells

Emerging evidence suggests that poor glycemic control mediates post-translational modifications to the H3 histone tail. We are only beginning to understand the dynamic role of some of the diverse epigenetic changes mediated by hyperglycemia at single loci, yet elevated glucose levels are thought to regulate genome-wide changes, and this still remains poorly understood. In this article we describe genome-wide histone H3K9/K14 hyperacetylation and DNA methylation maps conferred by hyperglycemia in primary human vascular cells. Chromatin immunoprecipitation (ChIP) as well as CpG methylation (CpG) assays, followed by massive parallel sequencing (ChIP-seq and CpG-seq) identified unique hyperacetylation and CpG methylation signatures with proximal and distal patterns of regionalization associative with gene expression. Ingenuity knowledge-based pathway and gene ontology analyses indicate that hyperglycemia significantly affects human vascular chromatin with the transcriptional up-regulation of genes involved in metabolic and cardiovascular disease. We have generated the first installment of a reference collection of hyperglycemia-induced chromatin modifications using robust and reproducible platforms that allow parallel sequencing-by-synthesis of immunopurified content. We uncover that hyperglycemia-mediated induction of genes and pathways associated with endothelial dysfunction occur through modulation of acetylated H3K9/K14 inversely correlated with methyl-CpG content.

Long term transcriptional reactivation of epigenetically silenced genes in colorectal cancer cells requires DNA hypomethylation and histone acetylation

Epigenetic regulation of genes involves the coordination of DNA methylation and histone modifications to maintain transcriptional status. These two features are frequently disrupted in malignancy such that critical genes succumb to inactivation. 5-aza-2'-deoxycytidine (5-aza-dC) is an agent which inhibits DNA methyltransferase, and holds great potential as a treatment for cancer, yet the extent of its effectiveness varies greatly between tumour types. Previous evidence suggests expression status after 5-aza-dC exposure cannot be explained by the DNA methylation status alone.

AIM

We sought to identify chromatin changes involved with short and long term gene reactivation following 5-aza-dC exposure. Two colorectal cancer cell lines, HCT116 and SW480, were treated with 5-aza-dC and then grown in drug-free media to allow DNA re-methylation. DNA methylation and chromatin modifications were assessed with bisulfite sequencing and Chromatin Immuno-Precipitation analysis.

RESULTS

Increased H3 acetylation, H3K4 tri-methylation and loss of H3K27 tri-methylation were associated with reactivation. Hypermethylated genes that did not show increased acetylation were transiently expressed with 5-aza-dC treatment before reverting to an inactive state. Three reactivated genes, CDO1, HSPC105 and MAGEA3, were still expressed 10 days post 5-aza-dC treatment and displayed localised hypomethylation at the transcriptional start site, and also an increased enrichment of histone H3 acetylation.

CONCLUSIONS

These observations suggest that hypomethylation alone is insufficient to reactivate silenced genes and that increased Histone H3 acetylation in unison with localised hypomethylation allows long term reversion of these epigenetically silenced genes. This study suggests that combined DNA methyltransferase and histone deacetylase inhibitors may aid long term reactivation of silenced genes.

Large-scale methylation domains mark a functional subset of neuronally expressed genes

DNA methylation is essential for embryonic and neuronal differentiation, but the function of most genomic DNA methylation marks is poorly understood. Generally the human genome is highly methylated (>70%) except for CpG islands and gene promoters. However, it was recently shown that the IMR90 human fetal lung fibroblast cells have large regions of the genome with partially methylated domains (PMDs, <70% average methylation), in contrast to the rest of the genome which is in highly methylated domains (HMDs, >70% average methylation). Using bisulfite conversion followed by high-throughput sequencing (MethylC-seq), we discovered that human SH-SY5Y neuronal cells also contain PMDs. We developed a novel hidden Markov model (HMM) to computationally map the genomic locations of PMDs in both cell types and found that autosomal PMDs can be >9 Mb in length and cover 41% of the IMR90 genome and 19% of the SH-SY5Y genome. Genomic regions marked by cell line specific PMDs contain genes that are expressed in a tissue-specific manner, with PMDs being a mark of repressed transcription. Genes contained within N-HMDs (neuronal HMDs, defined as a PMD in IMR90 but HMD in SH-SY5Y) were significantly enriched for calcium signaling, synaptic transmission, and neuron differentiation functions. Autism candidate genes were enriched within PMDs and the largest PMD observed in SH-SY5Y cells marked a 10 Mb cluster of cadherin genes with strong genetic association to autism. Our results suggest that these large-scale methylation domain maps could be relevant to interpreting and directing future investigations into the elusive etiology of autism.

Epigenome-wide association studies for common human diseases

Despite the success of genome-wide association studies (GWASs) in identifying loci associated with common diseases, a substantial proportion of the causality remains unexplained. Recent advances in genomic technologies have placed us in a position to initiate large-scale studies of human disease-associated epigenetic variation, specifically variation in DNA methylation. Such epigenome-wide association studies (EWASs) present novel opportunities but also create new challenges that are not encountered in GWASs. We discuss EWAS design, cohort and sample selections, statistical significance and power, confounding factors and follow-up studies. We also discuss how integration of EWASs with GWASs can help to dissect complex GWAS haplotypes for functional analysis.

High resolution detection and analysis of CpG dinucleotides methylation using MBD-Seq technology

Methyl-CpG binding domain protein sequencing (MBD-seq) is widely used to survey DNA methylation patterns. However, the optimal experimental parameters for MBD-seq remain unclear and the data analysis remains challenging. In this study, we generated high depth MBD-seq data in MCF-7 cell and developed a bi-asymmetric-Laplace model (BALM) to perform data analysis. We found that optimal efficiency of MBD-seq experiments was achieved by sequencing ∼100 million unique mapped tags from a combination of 500 mM and 1000 mM salt concentration elution in MCF-7 cells. Clonal bisulfite sequencing results showed that the methylation status of each CpG dinucleotides in the tested regions was accurately detected with high resolution using the proposed model. These results demonstrated the combination of MBD-seq and BALM could serve as a useful tool to investigate DNA methylome due to its low cost, high specificity, efficiency and resolution.

A comparative analysis of DNA methylation across human embryonic stem cell lines

Background

We performed a comparative analysis of the genome-wide DNA methylation profiles from three human embryonic stem cell (HESC) lines. It had previously been shown that HESC lines had significantly higher non-CG methylation than differentiated cells, and we therefore asked whether these sites were conserved across cell lines.

Results

We find that heavily methylated non-CG sites are strongly conserved, especially when found within the motif TACAG. They are enriched in splice sites and are more methylated than other non-CG sites in genes. We next studied the relationship between allele-specific expression and allele-specific methylation. By combining bisulfite sequencing and whole transcriptome shotgun sequencing (RNA-seq) data we identified 1,020 genes that show allele-specific expression, and 14% of CG sites genome-wide have allele-specific methylation. Finally, we asked whether the methylation state of transcription factor binding sites affects the binding of transcription factors. We identified variations in methylation levels at binding sites and found that for several transcription factors the correlation between the methylation at binding sites and gene expression is generally stronger than in the neighboring sequences.

Conclusions

These results suggest a possible but as yet unknown functional role for the highly methylated conserved non-CG sites in the regulation of HESCs. We also identified a novel set of genes that are likely transcriptionally regulated by methylation in an allele-specific manner. The analysis of transcription factor binding sites suggests that the methylation state of cis-regulatory elements impacts the ability of factors to bind and regulate transcription.

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.

Pipeline for large-scale microdroplet bisulfite PCR-based sequencing allows the tracking of hepitype evolution in tumors

Cytosine methylation provides an epigenetic level of cellular plasticity that is important for development, differentiation and cancerogenesis. We adopted microdroplet PCR to bisulfite treated target DNA in combination with second generation sequencing to simultaneously assess DNA sequence and methylation. We show measurement of methylation status in a wide range of target sequences (total 34 kb) with an average coverage of 95% (median 100%) and good correlation to the opposite strand (rho = 0.96) and to pyrosequencing (rho = 0.87). Data from lymphoma and colorectal cancer samples for SNRPN (imprinted gene), FGF6 (demethylated in the cancer samples) and HS3ST2 (methylated in the cancer samples) serve as a proof of principle showing the integration of SNP data and phased DNA-methylation information into “hepitypes” and thus the analysis of DNA methylation phylogeny in the somatic evolution of cancer.

Epigenomics of human embryonic stem cells and induced pluripotent stem cells: insights into pluripotency and implications for disease

Human pluripotent cells such as human embryonic stem cells (hESCs) and induced pluripotent stem cells (iPSCs) and their in vitro differentiation models hold great promise for regenerative medicine as they provide both a model for investigating mechanisms underlying human development and disease and a potential source of replacement cells in cellular transplantation approaches. The remarkable developmental plasticity of pluripotent cells is reflected in their unique chromatin marking and organization patterns, or epigenomes. Pluripotent cell epigenomes must organize genetic information in a way that is compatible with both the maintenance of self-renewal programs and the retention of multilineage differentiation potential. In this review, we give a brief overview of the recent technological advances in genomics that are allowing scientists to characterize and compare epigenomes of different cell types at an unprecedented scale and resolution. We then discuss how utilizing these technologies for studies of hESCs has demonstrated that certain chromatin features, including bivalent promoters, poised enhancers, and unique DNA modification patterns, are particularly pervasive in hESCs compared with differentiated cell types. We outline these unique characteristics and discuss the extent to which they are recapitulated in iPSCs. Finally, we envision broad applications of epigenomics in characterizing the quality and differentiation potential of individual pluripotent lines, and we discuss how epigenomic profiling of regulatory elements in hESCs, iPSCs and their derivatives can improve our understanding of complex human diseases and their underlying genetic variants.

Statistical approaches for the analysis of DNA methylation microarray data

Following the rapid development and adoption in DNA methylation microarray assays, we are now experiencing a growth in the number of statistical tools to analyze the resulting large-scale data sets. As is the case for other microarray applications, biases caused by technical issues are of concern. Some of these issues are old (e.g., two-color dye bias and probe- and array-specific effects), while others are new (e.g., fragment length bias and bisulfite conversion efficiency). Here, I highlight characteristics of DNA methylation that suggest standard statistical tools developed for other data types may not be directly suitable. I then describe the microarray technologies most commonly in use, along with the methods used for preprocessing and obtaining a summary measure. I finish with a section describing downstream analyses of the data, focusing on methods that model percentage DNA methylation as the outcome, and methods for integrating DNA methylation with gene expression or genotype data.

Next generation sequencing based approaches to epigenomics

Next generation sequencing has brought epigenomic studies to the forefront of current research. The power of massively parallel sequencing coupled to innovative molecular and computational techniques has allowed researchers to profile the epigenome at resolutions that were unimaginable only a few years ago. With early proof of concept studies published, the field is now moving into the next phase where the importance of method standardization and rigorous quality control are becoming paramount. In this review we will describe methodologies that have been developed to profile the epigenome using next generation sequencing platforms. We will discuss these in terms of library preparation, sequence platforms and analysis techniques.

Emerging patterns of epigenomic variation

Fuelled by new sequencing technologies, epigenome mapping projects are revealing epigenomic variation at all levels of biological complexity, from species to cells. Comparisons of methylation profiles among species reveal evolutionary conservation of gene body methylation patterns, pointing to the fundamental role of epigenomes in gene regulation. At the human population level, epigenomic changes provide footprints of the effects of genomic variants within the vast nonprotein-coding fraction of the genome, and comparisons of the epigenomes of parents and their offspring point to quantitative epigenomic parent-of-origin effects confounding classical Mendelian genetics. At the organismal level, comparisons of epigenomes from diverse cell types provide insights into cellular differentiation. Finally, comparisons of epigenomes from monozygotic twins help dissect genetic and environmental influences on human phenotypes and longitudinal comparisons reveal aging-associated epigenomic drift. The development of new bioinformatic frameworks for comparative epigenome analysis is putting epigenome maps within the reach of researchers across a wide spectrum of biological disciplines.

Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applications

Summary: A combination of bisulfite treatment of DNA and high-throughput sequencing (BS-Seq) can capture a snapshot of a cell's epigenomic state by revealing its genome-wide cytosine methylation at single base resolution. Bismark is a flexible tool for the time-efficient analysis of BS-Seq data which performs both read mapping and methylation calling in a single convenient step. Its output discriminates between cytosines in CpG, CHG and CHH context and enables bench scientists to visualize and interpret their methylation data soon after the sequencing run is completed.

Availability and implementation: Bismark is released under the GNU GPLv3+ licence. The source code is freely available from www.bioinformatics.bbsrc.ac.uk/projects/bismark/.

Contact:felix.krueger@bbsrc.ac.uk

Supplementary information:Supplementary data are available at Bioinformatics online.

Lsh, chromatin remodeling family member, modulates genome-wide cytosine methylation patterns at nonrepeat sequences

DNA methylation is critical for normal development and plays important roles in genome organization and transcriptional regulation. Although DNA methyltransferases have been identified, the factors that establish and contribute to genome-wide methylation patterns remain elusive. Here, we report a high-resolution cytosine methylation map of the murine genome modulated by Lsh, a chromatin remodeling family member that has previously been shown to regulate CpG methylation at repetitive sequences. We provide evidence that Lsh also controls genome-wide cytosine methylation at nonrepeat sequences and relate those changes to alterations in H4K4me3 modification and gene expression. Deletion of Lsh alters the allocation of cytosine methylation in chromosomal regions of 50 kb to 2 Mb and, in addition, leads to changes in the methylation profile at the 5' end of genes. Furthermore, we demonstrate that loss of Lsh promotes--as well as prevents--cytosine methylation. Our data indicate that Lsh is an epigenetic modulator that is critical for normal distribution of cytosine methylation throughout the murine genome.

Epigenetics and the brain: transcriptome sequencing reveals new depths to genomic imprinting

Transcriptome sequencing has identified more than a thousand potentially imprinted genes in the mouse brain. This comes as a revelation to someone who cut his teeth on the identification of imprinted genes when only a handful was known. Genomic imprinting, an epigenetic mechanism that determines expression of alleles according to sex of transmitting parent, was discovered over 25  years ago in mice but remains an enigmatic phenomenon. Why do these genes disobey the normal Mendelian logic of inheritance, do they function in specific processes, and how is their imprinting conferred? Next generation sequencing technologies are providing an unprecedented opportunity to survey the whole genome for imprinted genes and are beginning to reveal that imprinting may be more pervasive than we had come to believe. Such advances should lay the foundation for a definitive account of imprinting, but may also challenge accepted views on what it means to be imprinted. Editor's suggested further reading in BioEssays RNA as the substrate for epigenome-environment interactions Abstract.

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.

Genomic mapping of 5-hydroxymethylcytosine in the human brain

Methylation at the 5-position of cytosine is a well-studied epigenetic pathway. In addition to 5-methylcytosine (5mC), substantial amounts of 5-hydroxymethylcytosine (5hmC) also referred to as the sixth DNA base have been detected in certain tissues, most notably the brain. However, the genomic distribution of this cytosine modification is unknown. Here, we have used an immunoprecipitation technique (5hmC-IP) to examine the occurrence of 5hmC in DNA from human brain frontal lobe tissue. The distribution of 5hmC was compared to that of 5mC. We show that 5hmC is more selectively targeted to genes than is 5mC. 5hmC is particularly enriched at promoters and in intragenic regions (gene bodies) but is largely absent from non-gene regions. 5hmC peaks at transcription start sites did not correlate with gene expression levels for promoters with intermediate or high CpG content. However, the presence of 5hmC in gene bodies was more positively correlated with gene expression levels than was the presence of 5mC. Promoters of testis-specific genes showed strong 5mC peaks in brain DNA but were almost completely devoid of 5hmC. Our data provide an overview of the genomic distribution of 5hmC in human brain and will set the stage for further functional characterization of this novel DNA modification.

Genome studies at the PAG 2011 conference

The contents of the plenary lectures presented at the Plant and Animal Genome (PAG) meeting in January 2011 are summarized in order to provide some insights into the advances in plant, animal and microbe genome studies as they impact on our understanding of complex biological systems. The areas of biology covered include the dynamics of genome change, biological recognition processes and the new processes that underpin investment in science. This overview does not attempt to summarize the diversity of activities that are covered during the PAG through workshops, posters and the suppliers of cutting-edge technologies, but reviews major advances in specific research areas.

Putting epigenome comparison into practice

Epigenomic analysis efforts have so far focused on the multiple layers of epigenomic information within individual cell types. With the rapidly increasing diversity of epigenomically mapped cell types, unprecedented opportunities for comparative analysis of epigenomes are opening up. One such opportunity is to map the bifurcating tree of cellular differentiation. Another is to understand the epigenomically mediated effects of mutations, environmental influences, and disease processes. Comparative analysis of epigenomes therefore has the potential to provide wide-ranging fresh insights into basic biology and human disease. The realization of this potential will critically depend on availability of a cyberinfrastructure that will scale with the volume of data and diversity of applications and a number of other computational challenges.

A twin approach to unraveling epigenetics

The regulation of gene expression plays a pivotal role in complex phenotypes, and epigenetic mechanisms such as DNA methylation are essential to this process. The availability of next-generation sequencing technologies allows us to study epigenetic variation at an unprecedented level of resolution. Even so, our understanding of the underlying sources of epigenetic variability remains limited. Twin studies have played an essential role in estimating phenotypic heritability, and these now offer an opportunity to study epigenetic variation as a dynamic quantitative trait. High monozygotic twin discordance rates for common diseases suggest that unexplained environmental or epigenetic factors could be involved. Recent genome-wide epigenetic studies in disease-discordant monozygotic twins emphasize the power of this design to successfully identify epigenetic changes associated with complex traits. We describe how large-scale epigenetic studies of twins can improve our understanding of how genetic, environmental and stochastic factors impact upon epigenetics, and how such studies can provide a comprehensive understanding of how epigenetic variation affects complex traits.

Methylated DNA immunoprecipitation genome-wide analysis

DNA methylation is an important and potentially heritable component of the epigenetic machinery that has a major role in the control of gene expression and can be deregulated in many diseases. This modification of genomic DNA can be assessed using the methylated DNA immunoprecipitation (MeDIP) method, based on the quantification of methylated DNA fragments enriched using an antibody specific for methyl-cytosines.The relationship between the enrichment level and the real DNA methylation status is complex, and only a few methods have been developed to evaluate MeDIP enrichment measures to estimate the absolute or relative number of methyl-cytosines in a given sample. Two such methods are MEDME and BATMAN. This chapter focuses on the description and use of the former with a brief discussion of the latter.

From RNA-seq reads to differential expression results

Many methods and tools are available for preprocessing high-throughput RNA sequencing data and detecting differential expression.

A unique chromatin signature uncovers early developmental enhancers in humans

Cell-fate transitions involve the integration of genomic information encoded by regulatory elements, such as enhancers, with the cellular environment. However, identification of genomic sequences that control human embryonic development represents a formidable challenge. Here we show that in human embryonic stem cells (hESCs), unique chromatin signatures identify two distinct classes of genomic elements, both of which are marked by the presence of chromatin regulators p300 and BRG1, monomethylation of histone H3 at lysine 4 (H3K4me1), and low nucleosomal density. In addition, elements of the first class are distinguished by the acetylation of histone H3 at lysine 27 (H3K27ac), overlap with previously characterized hESC enhancers, and are located proximally to genes expressed in hESCs and the epiblast. In contrast, elements of the second class, which we term 'poised enhancers', are distinguished by the absence of H3K27ac, enrichment of histone H3 lysine 27 trimethylation (H3K27me3), and are linked to genes inactive in hESCs and instead are involved in orchestrating early steps in embryogenesis, such as gastrulation, mesoderm formation and neurulation. Consistent with the poised identity, during differentiation of hESCs to neuroepithelium, a neuroectoderm-specific subset of poised enhancers acquires a chromatin signature associated with active enhancers. When assayed in zebrafish embryos, poised enhancers are able to direct cell-type and stage-specific expression characteristic of their proximal developmental gene, even in the absence of sequence conservation in the fish genome. Our data demonstrate that early developmental enhancers are epigenetically pre-marked in hESCs and indicate an unappreciated role of H3K27me3 at distal regulatory elements. Moreover, the wealth of new regulatory sequences identified here provides an invaluable resource for studies and isolation of transient, rare cell populations representing early stages of human embryogenesis.

Systems approaches to molecular cancer diagnostics

The search for improved molecular cancer diagnostics is a challenge for which systems approaches show great promise. As is becoming increasingly clear, cancer is a perpetually-evolving, highly multi-factorial disease. With next generation sequencing providing an ever-increasing amount of high-throughput data, the need for analytical tools that can provide meaningful context is critical. Systems approaches have demonstrated an ability to separate meaningful signal from noise that arises from population heterogeneity, heterogeneity within and across tumors, and multiple sources of technical variation when sufficient sample sizes are obtained and standardized measurement technologies are used. The ability to develop clinically useful molecular cancer diagnostics will be predicated on advancements on two major fronts: 1) more comprehensive and accurate measurements of multiple endpoints, and 2) more sophisticated analytical tools that synthesize high-throughput data into meaningful reflections of cellular states. To this end, systems approaches that have integrated transcriptomic data onto biomolecular networks have shown promise in their ability to classify tumor subtypes, predict clinical progression, and inform treatment options. Ultimately, the success of systems approaches will be measured by their ability to develop molecular cancer diagnostics through distilling complex, systems-wide information into actionable information in the clinic.

Pash 3.0: A versatile software package for read mapping and integrative analysis of genomic and epigenomic variation using massively parallel DNA sequencing

Background

Massively parallel sequencing readouts of epigenomic assays are enabling integrative genome-wide analyses of genomic and epigenomic variation. Pash 3.0 performs sequence comparison and read mapping and can be employed as a module within diverse configurable analysis pipelines, including ChIP-Seq and methylome mapping by whole-genome bisulfite sequencing.

Results

Pash 3.0 generally matches the accuracy and speed of niche programs for fast mapping of short reads, and exceeds their performance on longer reads generated by a new generation of massively parallel sequencing technologies. By exploiting longer read lengths, Pash 3.0 maps reads onto the large fraction of genomic DNA that contains repetitive elements and polymorphic sites, including indel polymorphisms.

Conclusions

We demonstrate the versatility of Pash 3.0 by analyzing the interaction between CpG methylation, CpG SNPs, and imprinting based on publicly available whole-genome shotgun bisulfite sequencing data. Pash 3.0 makes use of gapped k-mer alignment, a non-seed based comparison method, which is implemented using multi-positional hash tables. This allows Pash 3.0 to run on diverse hardware platforms, including individual computers with standard RAM capacity, multi-core hardware architectures and large clusters.

The DNA methylome

Methylation of cytosines is a pervasive feature of eukaryotic genomes and an important epigenetic layer that is fundamental for cellular differentiation processes and control of transcriptional potential. DNA methylation patterns can be inherited and influenced by the environment, diet and aging, and disrupted in diseases. Complete DNA methylomes for several organisms are now available, helping clarify the evolutionary story of this epigenetic mark and its distribution in key genomic elements. Nonetheless, a complete understanding of its role, the mechanisms responsible for its establishment and maintenance, and its cross talk with other components of cellular machinery remains elusive.