Abstract:Background: The combination of sodium bisulfite treatment with highly-parallel sequencing is a common method for quantifying DNA methylation across the genome. The power to detect between-group differences in DNA methylation using bisulfite-sequencing approaches is influenced by both experimental (e.g. read depth, missing data and sample size) and biological (e.g. mean level of DNA methylation and difference between groups) parameters. There is, however, no consensus about the optimal thresholds for filtering bisulfite sequencing data with implications for the reproducibility of findings in epigenetic epidemiology. Results: We used a large reduced representation bisulfite sequencing (RRBS) dataset to assess the distribution of read depth across DNA methylation sites and the extent of missing data. To investigate how various study variables influence power to identify DNA methylation differences between groups, we developed a framework for simulating bisulfite sequencing data. As expected, sequencing read depth, group size, and the magnitude of DNA methylation difference between groups all impacted upon statistical power. The influence on power was not dependent on one specific parameter, but reflected the combination of study-specific variables. As a resource to the community, we have developed a tool, POWEREDBiSeq, which utilizes our simulation framework to predict study-specific power for the identification of DNAm differences between groups, taking into account user-defined read depth filtering parameters and the minimum sample size per group. Conclusions: Our data-driven approach highlights the importance of filtering bisulfite-sequencing data by minimum read depth and illustrates how the choice of threshold is influenced by the specific study design and the expected differences between groups being compared. The POWEREDBiSeq tool, which can be applied to different types of bisulfite sequencing data (e.g. RRBS, whole genome bisulfite sequencing (WGBS), targeted bisulfite sequencing and amplicon-based bisulfite sequencing), can help users identify the level of data filtering needed to optimize power and aims to improve the reproducibility of bisulfite sequencing studies.

A Classification Approach for DNA Methylation Profiling with Bisulfite Next-Generation Sequencing Data

Genome-scale sequence data processing and epigenetic analysis of DNA methylation]

A Comprehensive Evaluation of Alignment Software for Reduced Representation Bisulfite Sequencing Data

A Bayesian Framework to Identify Methylcytosines from High-Throughput Bisulfite Sequencing Data

Detecting differentially methylated regions in bisulfite sequencing data using quasi-binomial mixed models with smooth covariate effect estimates

Estimation of the methylation pattern distribution from deep sequencing data

A Hybrid Method Based on Beta-Binomial Bayesian Hierarchical Model and Combination of Ranking Method to Detect Differential Methylation Regions in Bisulfite Sequencing Data

A Framework for Analyzing DNA Methylation Data from Illumina Infinium HumanMethylation450 BeadChip

BayMeth: improved DNA methylation quantification for affinity capture sequencing data using a flexible Bayesian approach

A Bayesian Assignment Method for Ambiguous Bisulfite Short Reads

Data-Driven-Based Approach to Identifying Differentially Methylated Regions Using Modified 1D Ising Model.

Approaches for the Analysis and Interpretation of Whole-Genome Bisulfite Sequencing Data

methylClass : an R package to construct DNA methylation-based classification models

Addressing dispersion in mis‐measured multivariate binomial outcomes: A novel statistical approach for detecting differentially methylated regions in bisulfite sequencing data

Calling differentially methylated regions from whole genome bisulphite sequencing with DMRcate

Direction-aware functional class scoring enrichment analysis of Infinium DNA methylation data

BatMeth: Improved Mapper for Bisulfite Sequencing Reads on DNA Methylation.

Functional DNA methylation differences between tissues, cell types, and across individuals discovered using the M&M algorithm

Characterizing the properties of bisulfite sequencing data: maximizing power and sensitivity to identify between-group differences in DNA methylation

Rapid Quantification of DNA Methylation by Measuring Relative Peak Heights in Direct Bisulfite-Pcr Sequencing Traces

Amplicon-Based Bisulfite Conversion-NGS DNA Methylation Analysis Protocol