Bowtie2:快速空位读段比对

Fast Gapped-read Alignment with Bowtie 2
B. Langmead, S.L. Salzberg · Nat. Methods 9: 357-359, 2012 · ← 导航

1. 背景

Bowtie 1(2009)是首批 BWT 短读段比对器,但它有两个限制:(1) 不支持空位比对——无法处理 indel;(2) 默认 DFS 找到第一个命中即停,可能错过最佳比对。Bowtie 2 为解决这些限制而设计。

2. 核心算法

  1. FM-index 搜索: 通过反向搜索产生精确匹配种子
  2. 多种子策略: 找多个短种子(默认 20bp)而非一个长种子
  3. 动态规划拼接: 种子间用 banded DP 连接,允许空位
  4. 双方向搜索: 同时正向和反向互补搜索,避免 DFS 过早截断

3. 与 BWA-MEM 的关键对比

Bowtie 2 (2012)BWA-MEM (2013)
种子固定长度精确匹配(~20bp)变长 SMEM
索引标准 FM-indexFMD-index(双向)
比对模式默认端到端自动局部/端到端切换
嵌合读段不支持支持
>1kb 读段速度明显下降性能稳定

测试(BWA-MEM 论文):650bp 读段上 BWA-MEM 约快 6 倍。100bp 读段两者相当。BWA-MEM 的快来自 SMEM 一次遍历找所有种子 + banded DP 与 query 长度线性相关。

4. 历史定位

Bowtie 2 是从 Bowtie 1 到长读段比对的关键过渡产品。它在 TopHat 2(RNA-seq 分析)中作为底层比对器广泛使用。BWA-MEM 发表后,在多数场景下速度、准确率和功能(嵌合比对、超长读段)均超越 Bowtie 2,逐步成为行业标准。