Program in Bioinformatics and Integrative Biology
Cytosine; *DNA Methylation; Sequence Alignment; Sequence Analysis, DNA; *Software; Sulfites
Bioinformatics | Computational Biology
Analysis of bisulfite sequencing data usually requires two tasks: to call methylated cytosines (mCs) in a sample, and to detect differentially methylated regions (DMRs) between paired samples. Although numerous tools have been proposed for mC calling, methods for DMR detection have been largely limited. Here, we present Bisulfighter, a new software package for detecting mCs and DMRs from bisulfite sequencing data. Bisulfighter combines the LAST alignment tool for mC calling, and a novel framework for DMR detection based on hidden Markov models (HMMs). Unlike previous attempts that depend on empirical parameters, Bisulfighter can use the expectation-maximization algorithm for HMMs to adjust parameters for each data set. We conduct extensive experiments in which accuracy of mC calling and DMR detection is evaluated on simulated data with various mC contexts, read qualities, sequencing depths and DMR lengths, as well as on real data from a wide range of biological processes. We demonstrate that Bisulfighter consistently achieves better accuracy than other published tools, providing greater sensitivity for mCs with fewer false positives, more precise estimates of mC levels, more exact locations of DMRs and better agreement of DMRs with gene expression and DNase I hypersensitivity. The source code is available at http://epigenome.cbrc.jp/bisulfighter.
Rights and Permissions
Citation: Nucleic Acids Res. 2014 Apr;42(6):e45. doi: 10.1093/nar/gkt1373. Epub 2014 Jan 13. Link to article on publisher's site.
DOI of Published Version
Nucleic acids research
Saito, Yutaka; Tsuji, Junko; and Mituyama, Toutai, "Bisulfighter: accurate detection of methylated cytosines and differentially methylated regions" (2014). Open Access Articles. 2583.
Creative Commons License
This work is licensed under a Creative Commons Attribution-Noncommercial 3.0 License