De novo identification of differentially methylated regions in the human genome

Abstract

BackgroundThe identification and characterisation of differentially methylated regions (DMRs) between phenotypes in the human genome is of prime interest in epigenetics. We present a novel method, DMRcate, that fits replicated methylation measurements from the Illumina HM450K BeadChip (or 450K array) spatially across the genome using a Gaussian kernel. DMRcate identifies and ranks the most differentially methylated regions across the genome based on tunable kernel smoothing of the differential methylation (DM) signal. The method is agnostic to both genomic annotation and local change in the direction of the DM signal, removes the bias incurred from irregularly spaced methylation sites, and assigns significance to each DMR called via comparison to a null model.ResultsWe show that, for both simulated and real data, the predictive performance of DMRcate is superior to those of Bumphunter and Probe Lasso, and commensurate with that of comb-p. For the real data, we validate all array-derived DMRs from the candidate methods on a suite of DMRs derived from whole-genome bisulfite sequencing called from the same DNA samples, using two separate phenotype comparisons.ConclusionsThe agglomeration of genomically localised individual methylation sites into discrete DMRs is currently best served by a combination of DM-signal smoothing and subsequent threshold specification. The findings also suggest the design of the 450K array shows preference for CpG sites that are more likely to be differentially methylated, but its overall coverage does not adequately reflect the depth and complexity of methylation signatures afforded by sequencing.For the convenience of the research community we have created a user-friendly R software package called DMRcate, downloadable from Bioconductor and compatible with existing preprocessing packages, which allows others to apply the same DMR-finding method on 450K array data.Electronic supplementary materialThe online version of this article (doi:10.1186/1756-8935-8-6) contains supplementary material, which is available to authorized users.

Germline BRCA2 mutations drive prostate cancers with distinct evolutionary trajectories

“…We also filtered probes based on SNPs (65 probes) and non-CpG methylation probes (2,944 probes). Next, we used the DMRcate package (v1.2.0)51 to further filter out 27,296 probes that are known to cross-hybridize to multiple locations in the genome and 17,079 probes that contain a SNP with an annotated minor allele frequency of greater than 5% with a maximum distance of two nucleotides to the nearest CpG site. Average intensity levels were taken for technical replicates.…”

Unique DNA methylation signature in HPV-positive head and neck squamous cell carcinomas

“…While this threshold may lead to an underestimation of HPV-related methylation changes, we estimated this to be a better conservative approach that may avoid the identification of (or reduce the likelihood of identifying) false positive hits due to undetectable and/or not adjustable batch effects by the surrogate variable analysis applied in the study. Identification of DMRs was performed using the DMRcate R package [17]. Statistically significant DMRs (FDR <0.05) with at least three consecutive CpGs included in a bookend of 1000 nucleotides were retained for further analysis.…”

Effect of prenatal DHA supplementation on the infant epigenome: results from a randomized controlled trial

“…Testing for VMRs, independent of treatment group, and for differently methylated regions (DMRs) between the treatment groups was carried out using the Bioconductor package DMRcate [57], which extracts these regions via kernel density modelling. Since there were no limma-significant probes in the treatment group comparison (FDR-adjusted P  < 0.05) to search for DMRs, we relaxed the DMRcate default parameters that threshold significance based on the number of limma-significant probes.…”