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.

BASAL: a universal mapping algorithm for nucleotide base-conversion sequencing

A Bayesian Assignment Method for Ambiguous Bisulfite Short Reads

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

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

Leveraging Basecaller's Move Table to Generate a Lightweight k-mer Model

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

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

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

Valibs: A Visual Aligner for Bisulfite Sequences

De novo basecalling of m6A modifications at single molecule and single nucleotide resolution

Objective and Comprehensive Evaluation of Bisulfite Short Read Mapping Tools.

Biological Nanopore Approach for Single‐molecule Analysis of Nucleobase Modifications

An MCMC algorithm for base calling in sequencing-by-synthesis

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

Locating, tracing and sequencing multiple expanded genetic letters in complex DNA context via a bridge-base approach

BaseFormer: Transformer Based Base-Caller for Fast and Accurate Next Generation Sequencing

A Novel Multi-Alignment Pipeline for High-Throughput Sequencing Data.

An efficient Burrows-Wheeler transform-based aligner for short read mapping

Improving adenine and dual base editors through introduction of TadA-8e and Rad51DBD

A/C Simultaneous Conversion Using the Dual Base Editor in Human Cells.

Training Data Diversity Enhances the Basecalling of Novel RNA Modification-Induced Nanopore Sequencing Readouts