[논문 리뷰] XORRO: Rapid Paired-End Read Overlapper

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

Genomics and Phylogenetic Studies참고 문헌 9인용 수 3

한 줄 요약

XORRO는 64비트 레지스터에서 비트 수준 연산을 사용하여 이편 쌍 Illumina 리드를 매우 빠르게 겹치는 최적화된 저수준 소프트웨어 도구이다. 염기서열을 2비트 값으로 인코딩하고 XOR, 비트 시프트, popcount 명령어를 활용함으로써 XORRO는 수백만 개의 리드를 분당 이하로 겹치며, 기존 도구인 FLASH에 비해 최대 10배 빠르고 SHE-RA에 비해 수 개체 수준 빠른 성능을 보여줌.

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/

연구 동기 및 목표

• 다음세대 시퀀싱 파이프라인에서 리드 겹침 처리 속도가 데이터 생성 속도를 능가하는 계산적 병목 현상을 해결하기 위해.
• 낮은 인Del 비율을 가지는 Illumina 플랫폼에서 생성된 이편 쌍 리드에 특화된 빠르고 효율적인 알고리즘을 개발하기 위해.
• 기존의 앨리어먼트 알고리즘인 Needleman-Wunsch와 k-mer 방법이 리드 겹침에서 성능 한계를 극복하기 위해.
• 고성능 컴퓨팅(HPC) 인프라 없이도 일반 단일 CPU 시스템에서 실시간으로 이편 쌍 리드 겹침을 수행할 수 있도록 하기 위해.
• 이편 쌍 데이터로부터 더 긴, 더 정확한 공통 서열을 확보할 수 있도록 연구자들이 사용할 수 있는 경량, 빠르고 접근성 있는 솔루션을 제공하기 위해.

제안 방법

• 각 염기서열을 64비트 정수 내에서 2비트 값으로 표현함(A=00, C=01, G=10, T=11)으로써 한 연산당 32개 염기서열을 처리함.
• XOR 연산을 사용해 32염기 세그먼트 단위로 병렬적으로 겹치는 리드 쌍 간의 불일치 위치를 식별함.
• SSE4.2의 popcount 명령어를 활용해 한 CPU 사이클 내에서 불일치 수를 세어 유효한 겹침을 신속하게 필터링함.
• x86_64 어셈블리어를 사용해 루프 카운터 없이도 모든 가능한 겹침 정렬을 빠르게 검사함.
• 가장 긴 겹침에서 시작하여 사용자가 정의한 불일치 수 이하로 떨어지는 첫 번째 유효한 매칭을 찾을 때까지 내림차순으로 검색함.
• 원래의 이편 쌍 품질 점수의 합(또는 불일치 시 차이)을 기반으로 겹친 리드의 품질 점수를 재구성함.

실험 결과

연구 질문

• RQ1 저수준 비트 수준 알고리즘이 기존의 문자열 정렬 방법에 비해 이편 쌍 리드 겹침에서 뚜렷한 성능 향상을 이룰 수 있는가?
• RQ2 XOR 및 popcount과 같은 비트 수준 연산이 짧은 리드 시퀀싱에서 겹침 계산을 얼마나 빠르게 가속화할 수 있는가?
• RQ3 표준 64비트 단일 CPU 머신에서 근접한 실시간 이편 쌍 리드 겹침을 달성하는 것이 가능한가?
• RQ4 속도와 정확도 측면에서 XORRO의 성능은 FLASH 및 SHE-RA와 같은 기존 도구와 비교해 어떻게 되는가?
• RQ5 2비트 염기서열 인코딩과 어셈블리 최적화된 비트 조작을 사용하면 겹침 탐지의 계산 오버헤드를 줄일 수 있는가?

주요 결과

• XORRO는 표준 랩탑에서 3GB의 이편 쌍 리드를 단 4분 만에 겹쳐 높은 처리량을 입증함.
• 기준 비교에서 XORRO는 FLASH에 비해 10배 빠른 성능을 기록했고, SHE-RA에 비해 수 개체 수준의 빠른 성능을 보였음.
• 64비트 비트 수준 연산을 사용함으로써 비교당 필요한 명령어 수가 감소하여 CPU 오버헤드가 최소화됨.
• popcount 명령어 덕분에 한 CPU 사이클 내에서 불일치 수를 효율적으로 세어 유효한 겹침을 필터링하는 데 핵심적인 역할을 함.
• 사용자가 정의한 수준 이하의 불일치를 허용함으로써 높은 정확도를 유지하면서도 빠른 속도를 확보함.
• 옵션으로 제공되는 Perl 래퍼는 64비트 레지스터 한계를 초과하는 대규모 데이터셋에 대한 자동 매개변수 조정 및 배치 처리를 가능하게 함.