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.

wgbstools: A computational suite for DNA methylation sequencing data representation, visualization, and analysis

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

CpGtools: a python package for DNA methylation analysis

Approaches for the Analysis and Interpretation of Whole-Genome Bisulfite Sequencing Data

CGmapTools Improves the Precision of Heterozygous SNV Calls and Supports Allele-Specific Methylation Detection and Visualization in Bisulfite-Sequencing Data.

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

WBSA: Web Service for Bisulfite Sequencing Data Analysis

Generation of Whole-Genome Bisulfite Sequencing Libraries for Comprehensive DNA Methylome Analysis

DMAP2: A Pipeline for Analysis of Whole-Genome-Scale DNA Methylation Sequencing Data

msgbsR: An R package for analysing methylation-sensitive restriction enzyme sequencing data

Whole-Genome Bisulfite Sequencing Protocol for the Analysis of Genome-Wide DNA Methylation and Hydroxymethylation Patterns at Single-Nucleotide Resolution

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

An Integrated Package for Bisulfite DNA Methylation Data Analysis with Indel-sensitive Mapping

Amplicon-Based Bisulfite Conversion-NGS DNA Methylation Analysis Protocol

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

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

methylC Track: visual integration of single-base resolution DNA methylation data on the WashU EpiGenome Browser.

DNA methylation-calling tools for Oxford Nanopore sequencing: a survey and human epigenome-wide evaluation

wgatools: an ultrafast toolkit for manipulating whole genome alignments

RnBeads 2.0: comprehensive analysis of DNA methylation data

GobyWeb: simplified management and analysis of gene expression and DNA methylation sequencing data