[논문 리뷰] Metagenome assembly of high-fidelity long reads with hifiasm-meta

Xiaowen Feng, Haoyu Cheng외 2명|arXiv (Cornell University)|2021. 10. 16.

Genomics and Phylogenetic Studies인용 수 9

한 줄 요약

hifiasm-meta는 고정밀도 PacBio HiFi 리드를 위한 새로운 메타게놈 어셈블러로, 리드 선택, 히알로타입 인식 그래프 구축, 커버리지 기반 정제를 통해 어셈블리의 연속성과 정확성을 향상시킨다. 이는 각 데이터셋에서 수십에서 수백 개의 완전한 원형 박테리아 게놈을 재구성했으며, 실존적 및 모의 메타게놈 데이터셋에서 metaFlye와 HiCanu를 지속적으로 능가했다.

ABSTRACT

Current metagenome assemblers developed for short sequence reads or noisy long readswere not optimized for accurate long reads. Here we describe hifiasm-meta, a new metagenome assembler that exploits the high accuracy of recent data. Evaluated on seven empirical datasets, hifiasm-meta reconstructed tens to hundreds of complete circular bacterial genomes per dataset, consistently outperforming other metagenome assemblers.

연구 동기 및 목표

• 짧은 리드나 노이지한 장거리 리드에 최적화된 기존 메타게놈 어셈블러의 한계를 해결하기 위해, 고정밀도 장거리 리드에 특화된 어셈블러 개발
• PacBio HiFi 리드의 높은 정확도와 긴 길이 특성을 활용하여 복잡한 메타게놈 샘플의 *de novo* 어셈블리 향상
• 다양한 종이 혼합된 공동체에서 고형성, 변동 커버리지, 차이 읽기 포함 등으로 인한 어셈블리 오류와 분할을 줄이기
• 다양한 마이크로바이옴에서 완전, 거의 완전, 고품질의 메타게놈-어셈블드 게놈(MAGs)을 회수할 수 있는 견고하고 확장 가능한 어셈블러 개발

제안 방법

• 매우 높은 유전자 빈도를 가진 균주를 줄이면서도 저빈도 종의 리드를 유지하기 위해 선택적 리드 선택 단계 도입
• 차이 리드로 간주되어 폐기될 수 있는 저커버리지 게놈을 어셈블리 그래프 구축 과정에서 보호하기 위해, 해당 리드 유지
• 리드 제거 로직 수정: 포함된 리드가 모두 동일한 히알로타입에서 유래된 것으로 추정될 경우에만 제거
• 커버리지 기반 정제 적용: 한 단위티그의 커버리지 비율이 현저히 낮은 경우(예: rAB > 0.7, rAC < 0.7)에 오버랩을 제거하여 모호성 해결
• 커버리지 정보를 활용해 서로 다른 히알로타입의 단위티그를 연결하여 어셈블리 내 잔여 간극 메우기
• 추가로 포함된 리드로 인한 부가적 중복 서브그래프를 해결하기 위해 수정된 버블 터널링 전략 사용

실험 결과

연구 질문

• RQ1 고정밀도 장거리 리드를 대상으로 설계된 메타게놈 어셈블러가 노이지한 장거리 리드나 짧은 리드에 최적화된 도구들보다 뛰어난 연속성과 완전성을 달성할 수 있는가?
• RQ2 리드 선택이 메타게놈 어셈블리에서 저빈도 종의 회수에 어떤 영향을 미치는가?
• RQ3 히알로타입 인식 그래프 구축과 커버리지 기반 정제가 복잡한 메타게놈 샘플에서 분할과 오염을 어느 정도 줄일 수 있는가?
• RQ4 hifiasm-meta는 인간 및 동물 장내 마이크로바이옴을 포함한 다양한 메타게놈 데이터셋에서 일관되게 완전한 원형 박테리아 게놈을 재구성할 수 있는가?

주요 결과

• ATCC 모의 공동체에서 hifiasm-meta는 14개의 완전한 원형 박테리아 게놈을 재구성했으며, metaFlye와 HiCanu를 능가했고, 18% 빈도의 P. gingivalis 균주를 단일 컨티그로 완전히 복원했다.
• Zymo 모의 공동체에서 hifiasm-meta는 E. coli 균주(B766)를 완전한 원형 컨티그로 어셈블링했으며, metaFlye는 이를 실패했고, HiCanu는 두 균주에서 동일한 결과를 달성했다.
• 저빈도 종(예: M. smithii, 0.04%)의 경우 hifiasm-meta가 다른 도구들보다 더 연속적인 어셈블리 결과를 도출했지만, 여전히 커버리지 부족으로 인해 완전한 재구성에 한계가 있었다.
• 양A 장내 마이크로바이옴 데이터셋에서 hifiasm-meta는 14개의 거의 완전한 MAGs(완전성 ≥90%, 오염 ≤5%)와 23개의 고품질 MAGs(완전성 ≥70%, 오염 ≤10%)를 식별했으며, DAS Tool은 CheckM보다 22% 더 많은 거의 완전한 MAGs를 식별했다.
• 네 명의 인간 장내 샘플의 공통 어셈블리에서 hifiasm-meta는 16개의 ≥1Mb 컨티그를 생성했고, 그 중 7개는 원형이었으며, 계통발생 분석 결과 GTDB-Tk 분류 기준과 일치하는 명확한 고위계 분류군 클러스터링을 보였다.
• 모든 실존적 데이터셋에서 N50 및 품질 점수 측면에서 hifiasm-meta는 뛰어난 성능을 보였으며, 알고리즘 복잡도가 높아도 메모리 사용량과 런타임이 다른 도구들과 유사했다.