Whole genome bisulfite sequencing of cell-free DNA and its cellular contributors uncovers placenta hypomethylated domains

BackgroundCirculating cell-free fetal DNA has enabled non-invasive prenatal fetal aneuploidy testing without direct discrimination of the maternal and fetal DNA. Testing may be improved by specifically enriching the sample material for fetal DNA. DNA methylation may allow for such a separation of DNA; however, this depends on knowledge of the methylomes of circulating cell-free DNA and its cellular contributors.ResultsWe perform whole genome bisulfite sequencing on a set of unmatched samples including circulating cell-free DNA from non-pregnant and pregnant female donors and genomic DNA from maternal buffy coat and placenta samples. We find CpG cytosines within longer fragments are more likely to be methylated. Comparison of the methylomes of placenta and non-pregnant circulating cell-free DNA reveal many of the 51,259 identified differentially methylated regions are located in domains exhibiting consistent placenta hypomethylation across millions of consecutive bases. We find these placenta hypomethylated domains are consistently located within regions exhibiting low CpG and gene density. Differentially methylated regions identified when comparing placenta to non-pregnant circulating cell-free DNA are recapitulated in pregnant circulating cell-free DNA, confirming the ability to detect differential methylation in circulating cell-free DNA mixtures.ConclusionsWe generate methylome maps for four sample types at single-base resolution, identify a link between DNA methylation and fragment length in circulating cell-free DNA, identify differentially methylated regions between sample groups, and uncover the presence of megabase-size placenta hypomethylated domains.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0645-x) contains supplementary material, which is available to authorized users.     

Dnmt3a Protects Active Chromosome Domains against Cancer-Associated Hypomethylation

Changes in genomic DNA methylation patterns are generally assumed to play an important role in the etiology of human cancers. The Dnmt3a enzyme is required for the establishment of normal methylation patterns, and mutations in Dnmt3a have been described in leukemias. Deletion of Dnmt3a in a K-ras–dependent mouse lung cancer model has been shown to promote tumor progression, which suggested that the enzyme might suppress tumor development by stabilizing DNA methylation patterns. We have used whole-genome bisulfite sequencing to comprehensively characterize the methylomes from Dnmt3a wildtype and Dnmt3a-deficient mouse lung tumors. Our results show that profound global methylation changes can occur in K-ras–induced lung cancer. Dnmt3a wild-type tumors were characterized by large hypomethylated domains that correspond to nuclear lamina-associated domains. In contrast, Dnmt3a-deficient tumors showed a uniformly hypomethylated genome. Further data analysis revealed that Dnmt3a is required for efficient maintenance methylation of active chromosome domains and that Dnmt3a-deficient tumors show moderate levels of gene deregulation in these domains. In summary, our results uncover conserved features of cancer methylomes and define the role of Dnmt3a in maintaining DNA methylation patterns in cancer.     

Reducing the risk of false discovery enabling identification of biologically significant genome-wide methylation status using the HumanMethylation450 array

Background

The Illumina HumanMethylation450 BeadChip (HM450K) measures the DNA methylation of 485,512 CpGs in the human genome. The technology relies on hybridization of genomic fragments to probes on the chip. However, certain genomic factors may compromise the ability to measure methylation using the array such as single nucleotide polymorphisms (SNPs), small insertions and deletions (INDELs), repetitive DNA, and regions with reduced genomic complexity. Currently, there is no clear method or pipeline for determining which of the probes on the HM450K bead array should be retained for subsequent analysis in light of these issues.

Results

We comprehensively assessed the effects of SNPs, INDELs, repeats and bisulfite induced reduced genomic complexity by comparing HM450K bead array results with whole genome bisulfite sequencing. We determined which CpG probes provided accurate or noisy signals. From this, we derived a set of high-quality probes that provide unadulterated measurements of DNA methylation.

Conclusions

Our method significantly reduces the risk of false discoveries when using the HM450K bead array, while maximising the power of the array to detect methylation status genome-wide. Additionally, we demonstrate the utility of our method through extraction of biologically relevant epigenetic changes in prostate cancer.

Electronic supplementary material

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

Bisulfite sequencing reveals that Aspergillus flavus holds a hollow in DNA methylation

Aspergillus flavus first gained scientific attention for its production of aflatoxin. The underlying regulation of aflatoxin biosynthesis has been serving as a theoretical model for biosynthesis of other microbial secondary metabolites. Nevertheless, for several decades, the DNA methylation status, one of the important epigenomic modifications involved in gene regulation, in A. flavus remains to be controversial. Here, we applied bisulfite sequencing in conjunction with a biological replicate strategy to investigate the DNA methylation profiling of A. flavus genome. Both the bisulfite sequencing data and the methylome comparisons with other fungi confirm that the DNA methylation level of this fungus is negligible. Further investigation into the DNA methyltransferase of Aspergillus uncovers its close relationship with RID-like enzymes as well as its divergence with the methyltransferase of species with validated DNA methylation. The lack of repeat contents of the A. flavus' genome and the high RIP-index of the small amount of remanent repeat potentially support our speculation that DNA methylation may be absent in A. flavus or that it may possess de novo DNA methylation which occurs very transiently during the obscure sexual stage of this fungal species. This work contributes to our understanding on the DNA methylation status of A. flavus, as well as reinforces our views on the DNA methylation in fungal species. In addition, our strategy of applying bisulfite sequencing to DNA methylation detection in species with low DNA methylation may serve as a reference for later scientific investigations in other hypomethylated species.     

BSMAP: whole genome bisulfite sequence MAPping program

Background  

Bisulfite sequencing is a powerful technique to study DNA cytosine methylation. Bisulfite treatment followed by PCR amplification specifically converts unmethylated cytosines to thymine. Coupled with next generation sequencing technology, it is able to detect the methylation status of every cytosine in the genome. However, mapping high-throughput bisulfite reads to the reference genome remains a great challenge due to the increased searching space, reduced complexity of bisulfite sequence, asymmetric cytosine to thymine alignments, and multiple CpG heterogeneous methylation.     

Short‐term calorie restriction ameliorates genomewide,age‐related alterations in DNA methylation

DNA methylation plays major roles in many biological processes, including aging, carcinogenesis, and development. Analyses of DNA methylation using next‐generation sequencing offer a new way to profile and compare methylomes across the genome in the context of aging. We explored genomewide DNA methylation and the effects of short‐term calorie restriction (CR) on the methylome of aged rat kidney. Whole‐genome methylation of kidney in young (6 months old), old (25 months old), and OCR (old with 4‐week, short‐term CR) rats was analyzed by methylated DNA immunoprecipitation and next‐generation sequencing (MeDIP‐Seq). CpG islands and repetitive regions were hypomethylated, but 5′‐UTR, exon, and 3′‐UTR hypermethylated in old and OCR rats. The methylation in the promoter and intron regions was decreased in old rats, but increased in OCR rats. Pathway enrichment analysis showed that the hypermethylated promoters in old rats were associated with degenerative phenotypes such as cancer and diabetes. The hypomethylated promoters in old rats related significantly to the chemokine signaling pathway. However, the pathways significantly enriched in old rats were not observed from the differentially methylated promoters in OCR rats. Thus, these findings suggest that short‐term CR could partially ameliorate age‐related methylation changes in promoters in old rats. From the epigenomic data, we propose that the hypermethylation found in the promoter regions of disease‐related genes during aging may indicate increases in susceptibility to age‐related diseases. Therefore, the CR‐induced epigenetic changes that ameliorate age‐dependent aberrant methylation may be important to CR's health‐ and life‐prolonging effects.     

Hierarchical Clustering of Breast Cancer Methylomes Revealed Differentially Methylated and Expressed Breast Cancer Genes

Oncogenic transformation of normal cells often involves epigenetic alterations, including histone modification and DNA methylation. We conducted whole-genome bisulfite sequencing to determine the DNA methylomes of normal breast, fibroadenoma, invasive ductal carcinomas and MCF7. The emergence, disappearance, expansion and contraction of kilobase-sized hypomethylated regions (HMRs) and the hypomethylation of the megabase-sized partially methylated domains (PMDs) are the major forms of methylation changes observed in breast tumor samples. Hierarchical clustering of HMR revealed tumor-specific hypermethylated clusters and differential methylated enhancers specific to normal or breast cancer cell lines. Joint analysis of gene expression and DNA methylation data of normal breast and breast cancer cells identified differentially methylated and expressed genes associated with breast and/or ovarian cancers in cancer-specific HMR clusters. Furthermore, aberrant patterns of X-chromosome inactivation (XCI) was found in breast cancer cell lines as well as breast tumor samples in the TCGA BRCA (breast invasive carcinoma) dataset. They were characterized with differentially hypermethylated XIST promoter, reduced expression of XIST, and over-expression of hypomethylated X-linked genes. High expressions of these genes were significantly associated with lower survival rates in breast cancer patients. Comprehensive analysis of the normal and breast tumor methylomes suggests selective targeting of DNA methylation changes during breast cancer progression. The weak causal relationship between DNA methylation and gene expression observed in this study is evident of more complex role of DNA methylation in the regulation of gene expression in human epigenetics that deserves further investigation.     

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

Background

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

Methods

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

Results

The main features and advantages of this protocol are:

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

Conclusions

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

     

Reconstructing A/B compartments as revealed by Hi-C using long-range correlations in epigenetic data

Analysis of Hi-C data has shown that the genome can be divided into two compartments called A/B compartments. These compartments are cell-type specific and are associated with open and closed chromatin. We show that A/B compartments can reliably be estimated using epigenetic data from several different platforms: the Illumina 450 k DNA methylation microarray, DNase hypersensitivity sequencing, single-cell ATAC sequencing and single-cell whole-genome bisulfite sequencing. We do this by exploiting that the structure of long-range correlations differs between open and closed compartments. This work makes A/B compartment assignment readily available in a wide variety of cell types, including many human cancers.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0741-y) contains supplementary material, which is available to authorized users.     

Bloom syndrome patients and mice display accelerated epigenetic aging

Bloom syndrome (BSyn) is an autosomal recessive disorder caused by variants in the BLM gene, which is involved in genome stability. Patients with BSyn present with poor growth, sun sensitivity, mild immunodeficiency, diabetes, and increased risk of cancer, most commonly leukemias. Interestingly, patients with BSyn do not have other signs of premature aging such as early, progressive hair loss and cataracts. We set out to determine epigenetic age in BSyn, which can be a better predictor of health and disease over chronological age. Our results show for the first time that patients with BSyn have evidence of accelerated epigenetic aging across several measures in blood lymphocytes, as compared to carriers. Additionally, homozygous Blm mice exhibit accelerated methylation age in multiple tissues, including brain, blood, kidney, heart, and skin, according to the brain methylation clock. Overall, we find that Bloom syndrome is associated with accelerated epigenetic aging effects in multiple tissues and more generally a strong effect on CpG methylation levels.     

Genome-Wide Divergence of DNA Methylation Marks in Cerebral and Cerebellar Cortices

Background

Emerging evidence suggests that DNA methylation plays an expansive role in the central nervous system (CNS). Large-scale whole genome DNA methylation profiling of the normal human brain offers tremendous potential in understanding the role of DNA methylation in brain development and function.

Methodology/Significant Findings

Using methylation-sensitive SNP chip analysis (MSNP), we performed whole genome DNA methylation profiling of the prefrontal, occipital, and temporal regions of cerebral cortex, as well as cerebellum. These data provide an unbiased representation of CpG sites comprising 377,509 CpG dinucleotides within both the genic and intergenic euchromatic region of the genome. Our large-scale genome DNA methylation profiling reveals that the prefrontal, occipital, and temporal regions of the cerebral cortex compared to cerebellum have markedly different DNA methylation signatures, with the cerebral cortex being hypermethylated and cerebellum being hypomethylated. Such differences were observed in distinct genomic regions, including genes involved in CNS function. The MSNP data were validated for a subset of these genes, by performing bisulfite cloning and sequencing and confirming that prefrontal, occipital, and temporal cortices are significantly more methylated as compared to the cerebellum.

Conclusions

These findings are consistent with known developmental differences in nucleosome repeat lengths in cerebral and cerebellar cortices, with cerebrum exhibiting shorter repeat lengths than cerebellum. Our observed differences in DNA methylation profiles in these regions underscores the potential role of DNA methylation in chromatin structure and organization in CNS, reflecting functional specialization within cortical regions.     

Genome-wide DNA methylation changes in skeletal muscle between young and middle-aged pigs

Background

Age-related physiological, biochemical and functional changes in mammalian skeletal muscle have been shown to begin at the mid-point of the lifespan. However, the underlying changes in DNA methylation that occur during this turning point of the muscle aging process have not been clarified. To explore age-related genomic methylation changes in skeletal muscle, we employed young (0.5 years old) and middle-aged (7 years old) pigs as models to survey genome-wide DNA methylation in the longissimus dorsi muscle using a methylated DNA immunoprecipitation sequencing approach.

Results

We observed a tendency toward a global loss of DNA methylation in the gene-body region of the skeletal muscle of the middle-aged pigs compared with the young group. We determined the genome-wide gene expression pattern in the longissimus dorsi muscle using microarray analysis and performed a correlation analysis using DMR (differentially methylated region)-mRNA pairs, and we found a significant negative correlation between the changes in methylation levels within gene bodies and gene expression. Furthermore, we identified numerous genes that show age-related methylation changes that are potentially involved in the aging process. The methylation status of these genes was confirmed using bisulfite sequencing PCR. The genes that exhibited a hypomethylated gene body in middle-aged pigs were over-represented in various proteolysis and protein catabolic processes, suggesting an important role for these genes in age-related muscle atrophy. In addition, genes associated with tumorigenesis exhibited aged-related differences in methylation and expression levels, suggesting an increased risk of disease associated with increased age.

Conclusions

This study provides a comprehensive analysis of genome-wide DNA methylation patterns in aging pig skeletal muscle. Our findings will serve as a valuable resource in aging studies, promoting the pig as a model organism for human aging research and accelerating the development of comparative animal models in aging research.

Electronic supplementary material

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

High density DNA methylation array with single CpG site resolution

We have developed a new generation of genome-wide DNA methylation BeadChip which allows high-throughput methylation profiling of the human genome. The new high density BeadChip can assay over 480K CpG sites and analyze twelve samples in parallel. The innovative content includes coverage of 99% of RefSeq genes with multiple probes per gene, 96% of CpG islands from the UCSC database, CpG island shores and additional content selected from whole-genome bisulfite sequencing data and input from DNA methylation experts. The well-characterized Infinium® Assay is used for analysis of CpG methylation using bisulfite-converted genomic DNA. We applied this technology to analyze DNA methylation in normal and tumor DNA samples and compared results with whole-genome bisulfite sequencing (WGBS) data obtained for the same samples. Highly comparable DNA methylation profiles were generated by the array and sequencing methods (average R² of 0.95). The ability to determine genome-wide methylation patterns will rapidly advance methylation research.     

Quantification of epigenetic biomarkers: an evaluation of established and emerging methods for DNA methylation analysis

Background

DNA methylation is an important epigenetic mechanism in several human diseases, most notably cancer. The quantitative analysis of DNA methylation patterns has the potential to serve as diagnostic and prognostic biomarkers, however, there is currently a lack of consensus regarding the optimal methodologies to quantify methylation status. To address this issue we compared five analytical methods: (i) MethyLight qPCR, (ii) MethyLight digital PCR (dPCR), methylation-sensitive and -dependent restriction enzyme (MSRE/MDRE) digestion followed by (iii) qPCR or (iv) dPCR, and (v) bisulfite amplicon next generation sequencing (NGS). The techniques were evaluated for linearity, accuracy and precision.

Results

MethyLight qPCR displayed the best linearity across the range of tested samples. Observed methylation measured by MethyLight- and MSRE/MDRE-qPCR and -dPCR were not significantly different to expected values whilst bisulfite amplicon NGS analysis over-estimated methylation content. Bisulfite amplicon NGS showed good precision, whilst the lower precision of qPCR and dPCR analysis precluded discrimination of differences of < 25% in methylation status. A novel dPCR MethyLight assay is also described as a potential method for absolute quantification that simultaneously measures both sense and antisense DNA strands following bisulfite treatment.

Conclusions

Our findings comprise a comprehensive benchmark for the quantitative accuracy of key methods for methylation analysis and demonstrate their applicability to the quantification of circulating tumour DNA biomarkers by using sample concentrations that are representative of typical clinical isolates.

Electronic supplementary material

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