Abstract:MotivationThe rapid development of next-generation sequencing technology provides an opportunity to study genome-wide DNA methylation at single-base resolution. However, depletion of unmethylated cytosines brings challenges for aligning bisulfite-converted sequencing reads to a large reference. Software tools for aligning methylation reads have not yet been comprehensively evaluated, especially for the widely used reduced representation bisulfite sequencing (RRBS) that involves enrichment for CpG islands (CGIs).ResultsWe specially developed a simulator, RRBSsim, for benchmarking analysis of RRBS data. We performed extensive comparison of seven mapping algorithms for methylation analysis in both real and simulated RRBS data. Eighteen lung tumors and matched adjacent tissues were sequenced by the RRBS protocols. Our empirical evaluation found that methylation results were less consistent between software tools for CpG sites with low sequencing depth, medium methylation level, on CGI shores or gene body. These observations were further confirmed by simulations that indicated software tools generally had lower recall of detecting these vulnerable CpG sites and lower precision of estimating methylation levels in these CpG sites. Among the software tools tested, bwa-meth and BS-Seeker2 (bowtie2) are currently our preferred aligners for RRBS data in terms of recall, precision and speed. Existing aligners cannot efficiently handle moderately methylated CpG sites and those CpG sites on CGI shores or gene body. Interpretation of methylation results from these vulnerable CpG sites should be treated with caution. Our study reveals several important features inherent in methylation data, and RRBSsim provides guidance to advance sequence-based methylation data analysis and methodological development.Availability and implementationRRBSsim is a simulator for benchmarking analysis of RRBS data and its source code is available at https://github.com/xwBio/RRBSsim or https://github.com/xwBio/Docker-RRBSsim.Supplementary informationSupplementary data are available at Bioinformatics online.

Detection of differentially methylated regions from bisulfite-seq data by hidden Markov models incorporating genome-wide methylation level distributions

HMM-DM: identifying differentially methylated regions using a hidden Markov model

HMM-Fisher: identifying differential methylation using a hidden Markov model and Fisher’s exact test

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

Calling differentially methylated regions from whole genome bisulphite sequencing with DMRcate

DMseg: a Python algorithm for de novo detection of differentially or variably methylated regions

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

A family of mixture models for beta valued DNA methylation data

A comprehensive evaluation of computational tools to identify differential methylation regions using RRBS data

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

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

Detection of differentially methylated CpG sites between tumor samples with uneven tumor purities.

Evaluation of epitranscriptome-wide N6-methyladenosine differential analysis methods

A Bayesian framework for genome-wide inference of DNA methylation levels

Spectral Learning of Binomial HMMs for DNA Methylation Data

Identifying Differential Methylation in Cancer Epigenetics via a Bayesian Functional Regression Model

Gbdmr: identifying differentially methylated CpG regions in the human genome via generalized beta regressions

High accuracy methylation identification tools on single molecular level for PacBio HiFi data

Estimation of the methylation pattern distribution from deep sequencing data

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

A data-driven, high-throughput methodology to determine tissue-specific differentially methylated regions able to discriminate body fluids