[論文レビュー] XORRO: Rapid Paired-End Read Overlapper

Russell J. Dickson, Gregory B. Gloor|arXiv (Cornell University)|Apr 16, 2013

Genomics and Phylogenetic Studies参考文献 9被引用数 3

ひとこと要約

XORRO は、64ビットレジスタ上のビット演算を用いて、ペアエンド Illumina リードの重ね合わせを高速に実行する高度に最適化された低レベルソフトウェアツールである。ヌクレオチドを2ビット値にエンコードし、XOR、ビットシフト、popcount 指令を活用することで、数百万本のリードのサブミニット未満での重ね合わせを実現した。FLASH や SHE-RA といった既存ツールと比較して最大10倍の高速化を達成した。

ABSTRACT

Background: Computational analysis of next-generation sequencing data is outpaced by data generation in many cases. In one such case, paired-end reads can be produced from the Illumina sequencing method faster than they can be overlapped by downstream analysis. The advantages in read length and accuracy provided by overlapping paired-end reads demonstrates the necessity for software to efficiently solve this problem. Results: XORRO is an extremely efficient paired-end read overlapping program. XORRO can overlap millions of short paired-end reads in a few minutes. It uses 64-bit registers with a two bit alphabet to represent sequences and does comparisons using low-level logical operations like XOR, AND, bitshifting and popcount. Conclusions: As of the writing of this manuscript, XORRO provides the fastest solution to the paired-end read overlap problem. XORRO is available for download at: sourceforge.net/projects/xorro-overlap/

研究の動機と目的

• 次世代シーケンシングパイプラインにおける計算ボトルネック、特にリードの重ね合わせがデータ生成速度を上回る問題に対処する。
• Illumina プラットフォーム由来のペアエンドリードに特化した、高速で効率的なアルゴリズムを開発する。このリードはインDEL率が低いことが特徴である。
• 従来のアラインメントアルゴリズム（例：Needleman-Wunsch や k-mer 法）がリードの重ね合わせにおいて性能に限界を示す問題を克服する。
• HPC インfrastrucure を必要としない、コンsumer レベルの単一CPUシステムでも、実時間に近いペアエンドリードの重ね合わせを可能にする。
• 研究者がペアエンドデータからより長い、より正確なコンSENSUS配列を得るために、軽量で高速かつアクセスしやすいソリューションを提供する。

提案手法

• 1つのヌクレオチドを64ビット整数内に2ビット値として表現（A=00、C=01、G=10、T=11）し、1回の演算で32塩基を処理する。
• XOR 演算を用いて、32塩基セグメントごとに並列に重複するリードペア間の不一致位置を特定する。
• SSE4.2 で提供される popcount 指令を活用し、1つのCPUサイクルで不一致数をカウントすることで、有効な重ね合わせの高速フィルタリングを実現する。
• x86_64 アセンブリを用いて、ループカウンタを用いずに、ビットシフトおよびローテートを高速に実行し、すべての可能な重ね合わせアラインメントをテストする。
• 最も長い重ね合わせから開始し、ユーザーが定義した不一致数未満の有効な一致が得られるまで、順次短くしていく。
• 元のペアエンド品質スコアの合計（または不一致の場合の差）に基づいて、重ね合わせ済みリードの品質スコアを再構築する。

実験結果

リサーチクエスチョン

• RQ1 低レベルのビット演算アルゴリズムは、従来の文字列アラインメント手法に比べて、ペアエンドリードの重ね合わせで顕著な性能向上を達成できるか？
• RQ2 XOR や popcount といったビットレベル演算は、短鎖リードシーケンシングにおける重ね合わせ計算をどの程度高速化できるか？
• RQ3 標準的な64ビット単一CPUマシン上で、ほぼリアルタイムのペアエンドリードの重ね合わせが実現可能か？
• RQ4 XORRO の性能は、FLASH や SHE-RA といった既存ツールと比較して、速度および正確性の観点でどの程度優れているか？
• RQ5 2ビットヌクレオチドエンコードとアセンブリ最適化されたビット操作を用いることで、重ね合わせ検出の計算オーバーヘッドを低減できるか？

主な発見

• XORRO は標準のラップトップで3GBのペアエンドリードをわずか4分で重ね合わせ、高いスループットを実証した。
• ベンチマーク比較において、FLASH よりも10倍の高速化を達成し、SHE-RA よりも数個のオーダーも高速であった。
• 64ビットビット演算の活用により、比較処理に必要な命令数が削減され、CPU オーバーヘッドが最小限に抑えられた。
• popcount 指令により、1つのCPUサイクルで不一致数の効率的カウントが可能となり、有効な重ね合わせのフィルタリングに不可欠であった。
• ユーザーが定義した不一致数以内を許容することで、高速性を維持しつつ高い正確性を達成した。
• オプションのPerlラッパーにより、自動パrameterチューニングと、64ビットレジスタの上限を超える大規模データセットのバッチ処理が可能になった。