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.

oxBS-MLE: an efficient method to estimate 5-methylcytosine and 5-hydroxymethylcytosine in paired bisulfite and oxidative bisulfite treated DNA

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

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

A Strategy for Accurate Quantification of 5-Methylcytosine and 5-Hydroxymethylcytosine at CpG Sites Within Gene Promoter.

Ultrafast Bisulfite Sequencing Detection of 5-Methylcytosine in DNA and RNA

Direct enzymatic sequencing of 5-methylcytosine at single-base resolution

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

Base-resolution M5c 5 C Profiling Across the Mammalian Transcriptome by Bisulfite-Free Enzyme-Assisted Chemical Labeling Approach

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

DNA4mC-LIP: a Linear Integration Method to Identify N4-methylcytosine Site in Multiple Species

Base-resolution m5C profiling across the mammalian transcriptome by bisulfite-free enzyme-assisted chemical labeling approach

Identification and Quantification of 5-Methylcytosine and 5-Hydroxymethylcytosine on Random DNA Sequences by a Nanoconfined Electrochemiluminescence Platform

Estimation of the methylation pattern distribution from deep sequencing data

Simultaneous single-cell analysis of 5mC and 5hmC with SIMPLE-seq

Enzymatic methyl sequencing detects DNA methylation at single-base resolution from picograms of DNA

Modular Oxidation of Cytosine Modifications and Their Application in Direct and Quantitative Sequencing of 5-Hydroxymethylcytosine

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

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

Base-Resolution Analysis of 5-Hydroxymethylcytosine in the Mammalian Genome

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

A Bayesian Assignment Method for Ambiguous Bisulfite Short Reads