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.

BISCUIT: an efficient, standards-compliant tool suite for simultaneous genetic and epigenetic inference in bulk and single-cell studies

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

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

Single- Cell Bisulfite-Free 5mc and 5hmc Sequencing with High Sensitivity and Scalability

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

The SEQC2 epigenomics quality control (EpiQC) study

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

Analyzing single-cell bisulfite sequencing data with

High-coverage allele-resolved single-cell DNA methylation profiling by scDEEP-mC reveals cell lineage, X-inactivation state, and replication dynamics

Shortcut barcoding and early pooling for scalable multiplex single-cell reduced-representation CpG methylation sequencing at single nucleotide resolution

Genome-wide base-resolution mapping of DNA methylation in single cells using single-cell bisulfite sequencing (scBS-seq)

EPIC-TABSAT: analysis tool for targeted bisulfite sequencing experiments and array-based methylation studies

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

Droplet-based combinatorial indexing for massive scale single-cell epigenomics

Droplet-based bisulfite sequencing for high-throughput profiling of single-cell DNA methylomes

Assessment of batch-correction methods for scRNA-seq data with a new test metric

IMPLICON: an ultra-deep sequencing method to uncover DNA methylation at imprinted regions

mLiftOver: harmonizing data across Infinium DNA methylation platforms

LiBis: an ultrasensitive alignment augmentation for low-input bisulfite sequencing

Comprehensive evaluation of the Infinium human MethylationEPIC v2 BeadChip

User-friendly exploration of epigenomic data in single cells using sincei