BSMAP

BSMAP maps short reads from bisulfite sequencing to reference genomes to enable genome-wide analysis of DNA cytosine methylation, where unmethylated cytosines are converted to thymine following bisulfite treatment and PCR amplification.

Key Features:

• Bisulfite-aware mapping: Accounts for asymmetric alignments between cytosine and thymine resulting from bisulfite conversion.
• Genome hashing: Uses genome hashing to accelerate the search space for mapping bisulfite-treated reads.
• Bitwise masking: Employs bitwise masking techniques to improve mapping accuracy.
• Reduced-complexity handling: Addresses reduced sequence complexity and increased search space after bisulfite treatment.
• CpG heterogeneity support: Handles multiple CpG sites with heterogeneous methylation patterns.
• Whole-genome scalability: Maps high-throughput bisulfite reads across entire genomes while maintaining feasible memory and CPU usage.
• Performance optimization: Optimized for speed and sensitivity in mapping bisulfite sequencing data.

Scientific Applications:

• DNA methylation profiling: Determine methylation status of cytosines across the genome from next-generation bisulfite sequencing data.
• Epigenetics and disease studies: Support whole-genome analysis of DNA methylation patterns for studies of epigenetic modifications and disease-related methylation changes.
• High-throughput bisulfite sequencing analysis: Enable mapping of high-throughput bisulfite sequencing reads for downstream methylation analysis.

Methodology:

BSMAP implements a bisulfite-aware read mapping algorithm that leverages genome hashing and bitwise masking to handle C→T asymmetry and reduced sequence complexity, improving mapping speed and accuracy.