Continent-wide Multi-Segment Phylogenomic Study of Host and Geographic Structuring of African Hantaviruses

기술 요약: 아프리카 한타바이러스의 숙주 및 지리적 구조화를 규명하기 위한 대륙 규모의 다중 세그먼트 계통유전체 연구

문제 제기
아프리카 한타바이러스는 Hantaviridae 과 내에서 가장 특징이 덜 규명된 구성원 중 하나로 남아 있다. 이전의 분자 감시를 통해 아프리카 대륙 전역의 설치류, 땃쥐류, 박쥐에서 다양한 바이러스 계통이 확인되었으나, 단편적인 게놈 데이터로 인해 이들의 진화적 관계, 숙주 연관성 및 지리적 구조화에 대한 연구는 제한적이었다. 기존의 대부분의 연구는 단일 바이러스 종, 개별 게놈 세그먼트 또는 제한된 지리적 영역에 집중되었다. 결과적으로, 세 가지 코딩 세그먼트(L, M, S) 모두에 걸쳐 진화적 관계, 숙주 연관성 및 게놈 일치성을 통합하는 대륙 규모의 비교 프레임워크가 결여되어 있었다. 또한, 서로 다른 유전 영역에서 숙주 및 지리적 연관된 계통 패턴이 얼마나 일관되게 회복되는지에 대해서도 제대로 규명되지 않았다.

연구 방법론
본 연구는 281개의 전체 길이 및 부분 코딩 서열(L n=166, M n=48, S n=67)을 사용하여 대륙 규모의 비교 계통유전체 분석을 수행하였다. 이 서열들은 2004년부터 2024년 사이에 13개 아프리카 국가(이집트, 세네갈, 베냉, 말리, 기니, 코트디부아르, 시에라리온, 가봉, DRC, 에티오피아, 케냐, 탄자니아, 마다가스카르)에서 수집된 것들과 대표적인 글로벌 참조 균주들을 포함하여 NCBI GenBank에서 추출되었다.

분석 워크플로우는 다음과 같다:

• 서열 정렬(Sequence Alignment): 각 세그먼트의 뉴클레오타이드 서열을 MAFFT v7을 사용하여 독립적으로 정렬한 후, 위치적 상동성을 확보하고 모호하게 정렬된 영역을 제거하기 위해 AliView에서 수동 편집을 수행하였다.
• 계통 추론(Phylogenetic Inference): ModelFinder를 통해 최적의 치환 모델을 선택한 후, IQ-TREE v2.1.3을 사용하여 각 세그먼트에 대한 최대 우도(Maximum-likelihood, ML) 계통수를 구축하였다. 노드 지지도는 1,000회의 울트라 패스트 부트스트랩 반복(ultrafast bootstrap replicates)을 통해 평가되었다.
• 시간적 신호 평가(Temporal Signal Assessment): 시간적 신호를 평가하기 위해 Clockor2를 이용한 root-to-tip 회귀 분석을 수행하였다. 분석 결과 약한 신호(L, M, S에 대해 각각 R² 값 0.069, 0.021, 0.046)가 나타났으며, 이는 신뢰할 수 있는 분자 시계 추론 및 시간 보정 베이지안 분석을 불가능하게 하였다.
• 비교 분석(Comparative Analysis): 유전적 재조합(reassortment)의 신호를 식별하고 위상적 일치성(topological congruence)을 평가하기 위해 Dendroscope를 사용하여 쌍별 탠글그램(pairwise tanglegram) 비교(L–M, L–S, M–S)를 생성하였다.

주요 결과

• 숙주 연관 구조화(Host-Associated Structuring): 세 가지 세그먼트 모두에서 계통 분석 결과 뚜렷한 숙주 연관 구조화가 나타났다. 설치류, 땃쥐류, 박쥐 관련 한타바이러스는 별개의 진화적 클레이드(clade)를 형성하였다. 박쥐 및 땃쥐 매개 계통은 설치류 관련 바이러스에 비해 기저의 분기된 위치를 차지했다. 설치류 관련 그룹 내에서는 숙주 종(예: 기니의 Hylomyscus simus 대 케냐의 Hylomyscus endorobae)과 지리에 의해 차별화되었다.
• 지리적 구조화(Geographic Structuring): 아프리카 토착 한타바이러스는 뚜렷한 지역적 클러스터링을 보이며, 별개의 마다가스카르, 동아프ic, 서아프리카 집단을 형성하였다.
  • 서아프리카: 토착 계통(예: 기니의 Sangassou 바이러스)과 전 세계적으로 분포하는 설치류 관련 서울 한타바이러스(SEOV) 계통을 모두 포함한다.
  • 동아프리카: 케냐, 탄자니아, 에티오피아의 바이러스는 상당한 유전적 다양성을 가진 별도의 지역 집단을 형성하였다.
  • 마다가스카르: 바이러스는 별개의 잘 지지되는 단계통(monophyletic) 계통을 형성하였으며, 이는 아프리카 본토 바이러스와 유의미한 유전적 차이가 있음을 나타낸다.
  • 북아프리카: 이집트 서열은 서부 팔레아틱(Western Palearctic) 및 유라시아 무리네(Murinae) 관련 계통과 클러스터를 형성하였다.
  • SEOV 분포: 세네갈, 베냉, 말리의 서열은 인간과 관련된 유럽 및 아시아의 격리 균주를 포함한 전 세계적으로 분포하는 SEOV 계통 내에 클러스터링되었으며, 이는 인간 거주지 근처의 설치류를 통한 인위적 확산을 시사한다.
• 세그먼트별 해상도(Segment-Specific Resolution): 세 가지 게놈 세그gment는 상호 보완적인 계통 신호를 제공하였다:
  • L 세그먼트(RNA 의존 RNA 중합효소 인코딩)는 깊은 진화적 관계를 해결하고 주요 숙주 관련 그룹을 분리하는 강력한 프레임워크를 제공하였다.
  • M 세그먼트(당단백질 전구체 인코딩)는 밀접하게 관련된 숙주 및 지리적 그룹 간의 명확한 구분을 제공하였다.
  • S 세그먼트(뉴클레오캡시드 단백질 인코딩)는 깊은 노드에서의 해상도는 낮았으나, 주요 계통 경계를 일관되게 회복하였다.
• 위상적 일치성 및 재조합(Topological Congruence and Reassortment): 쌍별 탠글그램 분석 결과 세그먼트 간에 상당한 위상적 일치성이 나타났으며, L-M 사이에서 가장 높은 일치도를 보였다. 불일치는 주로 L-S 비교의 말단 가지(terminal branches)에 국한되어 나타났다. 이러한 국소적 불일치는 유전적 재조한(reassortment)의 신호와 일치하지만 결정적인 증거는 아니다. 주요 계통 구조가 유지되고 있다는 점은 광범적인 재조합이 샘플링된 바이러스의 진화적 구조를 주도하는 주요 동력은 아님을 시사한다.

의의 및 주장
저자들은 본 연구가 아프리카 한타바이러스에 대한 최초의 통합된 대륙 규모 계통유전체 프레임워크를 구축했다고 주장한다. 주요 의의는 아프리카 한타바이러스의 다양성이 숙주 매개체 연관성과 지리에 의해 강력하게 구조화되어 있으며, 주요 진화적 관계가 세 가지 게놈 세그먼트 모두에서 폭넓게 보존되어 있다는 점을 입증한 데 있다.

본 논문은 이러한 발견이 다음의 기초를 제공한다고 상정한다:

1. 미래의 게놈 감시: 바이러스 다양성 및 진화적 변화를 모니터링하기 위한 기준점 확립.
2. 원 헬스(One Health) 전략: 전통적인 설치류 숙주를 넘어 땃쥐와 박쥐까지 감시 범위를 확대해야 할 필요성을 강조하고, 도시화와 무역이 SEOV 확산에 미치는 역할을 인식함으로써 신종 인수공통 감염병에 대한 감시 및 대비를 강화함.
3. 진단 개선: 아프리카 토착 계통의 유전적 차이를 강조함으로써, 유럽이나 미국 바이러스를 대상으로 개발된 진단법이 민감도가 부족할 수 있음을 시사하며, 지역적으로 적응된 접근 방식의 필요성을 제시함.

저자들은 연구의 한계에 대해 신중한 어조를 유지하며, 약한 시간적 신호로 인해 분자 시계 측정이 불가능했다는 점과 중앙 및 남부 아프리카의 불균형한 지리적/숙주 샘플링(특히 이 지역의 경우)으로 인해 현재의 데이터셋이 대륙의 실제 한타바이러스 다양성의 일부만을 나타낼 가능성이 높다는 점을 인정하였다. 그들은 본 연구가 알려진 계통의 진화적 구조를 명확히 하지만, 아프리카 한타바이러스 다양성에 대한 포괄적인 이해를 위해서는 확장된 샘플링과 완전한 게놈 가용성이 필수적이라고 결론지었다.