BSMAP: whole genome bisulfite sequence MAPping program.
BACKGROUND
Bisulfite sequencing is a powerful technique to study DNA cytosine methylation. Bisulfite treatment followed by PCR amplification specifically converts unmethylated cytosines to thymine. Coupled with next generation sequencing technology, it is able to detect the methylation status of every cytosine in the genome. However, mapping high-throughput bisulfite reads to the reference genome remains a great challenge due to the increased searching space, reduced complexity of bisulfite sequence, asymmetric cytosine to thymine alignments, and multiple CpG heterogeneous methylation.


RESULTS
We developed an efficient bisulfite reads mapping algorithm BSMAP to address the above issues. BSMAP combines genome hashing and bitwise masking to achieve fast and accurate bisulfite mapping. Compared with existing bisulfite mapping approaches, BSMAP is faster, more sensitive and more flexible.


CONCLUSIONS
BSMAP is the first general-purpose bisulfite mapping software. It is able to map high-throughput bisulfite reads at whole genome level with feasible memory and CPU usage. It is freely available under GPL v3 license at http://code.google.com/p/bsmap/.
Xi
Yuanxin
Y
Division of Biostatistics, Dan L Duncan Cancer Center, Department of Molecular and Cellular Biology, Baylor College of Medicine, One Baylor Plaza, Houston, TX 77030, USA. yxi@bcm.edu
Li
Wei
W
D001483:Base Sequence; D004247:DNA / Q000737:chemistry; D019175:DNA Methylation; D023281:Genomics / Q000379:methods*; D016415:Sequence Alignment; D017422:Sequence Analysis, DNA / Q000379:methods*; D012984:Software; D013447:Sulfites / Q000737:chemistry
D013447:Sulfites; D004247:DNA; C042345:hydrogen sulfite