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.

BSReadSim: a versatile and efficient simulator to generate realistic bisulfite sequencing reads

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

Bs-Seeker2: A Versatile Aligning Pipeline for Bisulfite Sequencing Data

BitMapperBS: a fast and accurate read aligner for whole-genome bisulfite sequencing

SCSsim: an Integrated Tool for Simulating Single-Cell Genome Sequencing Data.

BS-RNA: An efficient mapping and annotation tool for RNA bisulfite sequencing data

Objective and Comprehensive Evaluation of Bisulfite Short Read Mapping Tools.

scReadSim: a single-cell RNA-seq and ATAC-seq read simulator

SimuSCoP: Reliably Simulate Illumina Sequencing Data Based on Position and Context Dependent Profiles

Pirs: Profile-Based Illumina Pair-End Reads Simulator

NPBSS: a New PacBio Sequencing Simulator for Generating the Continuous Long Reads with an Empirical Model

TAMeBS: A Sensitive Bisulfite-Sequencing Read Mapping Tool for DNA Methylation Analysis

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

Systematic Review of Next-Generation Sequencing Simulators: Computational Tools, Features and Perspectives

A Bayesian Assignment Method for Ambiguous Bisulfite Short Reads

BSXplorer: analytical framework for exploratory analysis of BS-seq data

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

InSilicoSeq 2.0: Simulating realistic amplicon-based sequence reads

Using Local Alignment to Enhance Single-Cell Bisulfite Sequencing Data Efficiency.

PCBS: an R package for fast and accurate analysis of bisulfite sequencing data

BEERS2: RNA-Seq simulation through high fidelity in silico modeling