Developing a BS-Seq Analysis Workflow for Genomic Variation and Methylation Level Calling
Next-Generation Sequencing technologies combined with a bisulfite treatment of DNA (BS-Seq) present an efficient method in the field of epigenomics, allowing for a precise analysis of methylation patterns at single-nucleotide resolution. The expanding field leads to the production of enormous amounts of data, that need to be analysed using precise and efficient software tools. The analysis needs to address the specific challenges arising from the bisulfite treatment, which converts unmethylated Cs to Ts. We present a powerful analysis workflow, implemented using the C++ software library SeqAn, for accurate bisulfite read mapping, the detection of single-nucleotide polymorphisms and methylation level calling at single-nucleotide resolution. Furthermore, we provide a multiple sequence realignment method to improve the alignment accuracy for indel reads. We show that in comparison to existing methods, we achieve comparable and even better results.
 Xi, Yuanxin, and Wei Li. 2009. “BSMAP: whole genome bisulfite sequence MAPping program..” BMC Bioinformatics 10: 232. doi:10.1186/1471-2105-10-232.
 Chen, Pao-Yang, Shawn J Cokus, and Matteo Pellegrini. 2010. “BS Seeker: precise mapping for bisulfite sequencing..” BMC Bioinformatics 11: 203. doi:10.1186/1471-2105-11-203.
 Lyko, Frank, Sylvain Foret, Robert Kucharski, Stephan Wolf, Cassandra Falckenhayn, and Ryszard Maleszka. 2010. “The Honey Bee Epigenomes: Differential Methylation of Brain DNA in Queens and Workers.” PLoS biology 8 (11): –. doi:10.1371/journal.pbio.1000506.