Genome-wide characterization of Plasmodium vivax infections in local travelers and non-travelers from the Peruvian Amazon

기술 요약: 현지 여행자 및 비여행자의 페루 아마존 지역 내 Plasmodium vivax 감염에 대한 전전체 게놈 특성 규명

문제 제기
지역 사회가 말라리아 퇴치 단계로 나아감에 따라, 국지적으로 획득된 감염과 외부에서 유입된 감염을 구분하는 것은 표적 감시를 위해 매우 중요하다. 그러나 높은 인구 이동성, 낮은 기생충 유전 다양성, 그리고 역사적 인구 병목 현상을 특징으로 하는 지역(예: 페루 아마존)에서는 미세한 공간 척도에서 감염의 기원을 구별하는 것이 여전히 도전적인 과제이다. 기존의 분자 도구(예: 마이크로새틀라이트, SNP 바코드)는 기생충 집단이 높은 친연성을 보이는 경우, 인접한 공동체 간의 차이를 구별하거나 최근의 유입 신호를 식별하는 데 해상도가 부족할 수 있다. 또한, 휴면 상태의 히포조이트(hypnozoite)로부터의 재발 가능성이라는 P. vivax의 생물학적 특성은 여행 후 발견된 감염이 새로운 유입인지, 아니면 이전의 국지적 감염의 재활성화인지를 해석하는 과정을 복잡하게 만든다. 본 연구는 전전체 게놈 시퀀싱(WGS)이 높은 유전적 유사성과 불완전한 이동성 데이터를 가진 맥락에서 여행자 관련 감염과 국지적 전파를 효과적으로 구별할 수 있는지에 대한 이해의 격차를 다룬다.

연구 방법론
본 연구는 2021년에서 2022년 사이에 페루 로레토(Loreto) 지역 마잔(Mazan) 지구의 세 공동체(Libertad, n=24; Gamitanacocha, n=23; Urco Miraño, n=9)에서 수집된 56건의 확진된 P. vivax 감염 사례를 대상으로 단면 연구 설계를 활용하였다. 참가자들은 지난 한 달 동안 거주지 외 지역에서 최소 하룻밤 이상 체류한 경우 "여행자"(n=16)로, 그렇지 않은 경우 "비여행자"(n=40)로 분류되었다.

• 시료 처리: 전혈에서 DNA를 추출하였으며, 미토콘드리아 cox1 유전자를 표적으로 하는 qPCR을 통해 기생충 종을 확인하였다. 시퀀싱 전 기생충 DNA를 농축하기 위해 선택적 전전체 게놈 증폭(sWGA)법을 사용하였다.
• 시퀀싱: 라이브러리는 PCR-free 프로토콜을 사용하여 준비되었으며, Illumina NovaSeq 플랫폼에서 평균 약 100× 커버리지의 150 bp paired-end 리드로 시퀀싱되었다.
• 데이터 분석:
  • 변이 호출(Variant Calling): 리드를 P. vivax 참조 서열(PvPAM)에 정렬하였고, GATK를 사용하여 변이를 호출하였다. 최종적으로 16,485개의 고품질 이형 접합(biallelic) SNP 세트를 분석에 유지하였다.
  • 유전 다양성: 공동체 및 여행 그룹 내/간의 다양성을 평가하기 위해 기대 이형 접합성($H_e$)과 뉴클레오타이드 다양성($\pi$)을 계산하였다.
  • 집단 구조: 주성분 분석(PCA), 판별 분석(DAPC), ADMIXTURE를 사용하여 유전적 클러스터링과 조상 비율을 식별하였다.
  • 친연성: 공유된 조상 및 최근 전파 연결을 탐지하기 위해 isoRelate 패키지를 사용하여 Identity-by-Descent(IBD) 분석을 수행하였다.
  • 네트워크 분석: 유전적 거리와 연결성을 시각화하기 위해 최소 스패닝 트리(MST)와 계층적 클러스터링을 생성하였다.

주요 결과

• 유전 다양성: Libertad가 가장 높은 전체 유전 다양성($H_e$ = 0.363)을 보였으며, Urco Miraño(0.198), Gamitanacocha(0.1) 순이었다. Libertad 내에서는 여행자 관련 감염이 비여행자보다 유의미하게 높은 다양성을 나타냈다($H_e$ = 0.432 vs. 0.304; $\pi$ = 0.290 vs. 0.168).
• 집단 분화: 쌍별 $F_{st}$ 분석 결과, Libertad와 Gamitanacocha 사이의 분화가 낮게 나타났으며($F_{st}$ = 0.056), 이는 공유된 전파 네트워크를 시사한다. 반면, Urco Miraño는 Libertad($F_{st}$ = 0.548) 및 Gamitanacocha($F_{st}$ = 0.773) 모두와 높은 분화를 보였다.
• 클러스터링 및 구조: PCA와 DAPC는 Urco Miraño의 기생충이 별도의 클러스터를 형성함을 확인하였으나, Libertad와 Gamitanacocha는 부분적인 중첩을 보였다. 여행자 관련 감염은 별도의 전역적 클러스터를 형성하지 않았으며, 대신 거주 공동체의 클러스터 내에 연결되어 분포하거나 그 사이에 끼어 있는 양상을 보였다.
• 조상 및 친연성: ADMIXTURE 분석 결과, Libertad의 여행자들 사이에서 더 큰 이질성이 나타났으며 일부는 혼합된 조상 성분을 보였다. 높은 임계값(50–80%)에서의 IBD 분석을 통해 두 건의 여행자 관련 감염(Libertad의 L1 및 Urco Miraño의 U1)이 유전적으로 구별되는 단일 개체(singleton)임을 확인하였으며, 이는 샘플링되지 않은 지역의 계통일 가능성이 있다. 그러나 대부분의 여행자 감염은 국지적 전파 네트워크 내에 포함되어 있었다.
• 연결성: 네트워크는 낮은 IBD 임계값에서 높은 연결성을 보였으며, Gamitanacocha 내부 및 Libertad와 Gamitanacocha 사이에서 강한 친연성이 관찰되었다. Urco Miraño와 나머지 두 공동체 사이에서는 관련된 쌍이 발견되지 않았다.

의의 및 주장
저자들은 본 연구가 미세한 공간 척도에서 여행자 관련 감염과 국지적으로 획득된 P. vivax 감염을 구별하는 데 있어 WGS의 잠재력과 한계를 입증한다고 주장한다. 연구 결과는 이동성과 관련된 전파를 해석하기 위해 유전체 데이터와 역학 데이터를 통합하는 것을 지지한다. 구체적으로, 본 연구는 다음과 같이 결론짓는다:

1. WGS는 공동체 수준의 구조를 해결할 수 있음: 고립된 공동체(Urco Miraño)와 연결된 공동체(Libertad 및 Gamitanacocha) 간의 뚜렷한 유전적 클러스터를 성공적으로 식별하였다.
2. 여행 보고가 항상 유입을 의미하는 것은 아님: 대부분의 여행자 관련 감염이 국지적 기생충과 클러스터를 형성했다는 점은, 보고된 여행이 반드시 유입을 의미하지 않음을 시사한다. 즉, 많은 감염이 국지적 전파 네트워크 내에서 획득되었거나 재발에 의한 것일 수 있다.
3. 맥락 의존적 해상도: 감염 기원을 구별하는 능력은 유전적 분화 정도와 참조 데이터의 가용성에 달려 있다. 다양성이 낮고 친연성이 높은 환경에서는 밀도 높은 샘플링과 강력한 역학적 맥락 없이는 WGS 단독으로는 충분하지 않을 수 있다.
4. 이상치(Outliers) 식별: 대부분의 여행자가 국지적으로 클러스터링되었으나, WGS는 뚜렷한 유전적 프로파일을 가진 특정 여행자 관련 감염을 식별해 냈으며, 이는 충분한 분화가 존재하는 경우 이동성 관련 계통을 탐지하는 데 있어 해당 방법의 유용성을 시사한다.

저자들은 작은 표본 크기, 특히 여행자 그룹의 경우와 P. vivax 생물학의 본질적인 복잡성(히포조이트)으로 인해 이러한 결과 해석에 주의가 필요함을 강조한다. 저자들은 말라리아 퇴치를 위한 유전체 감시 전략을 정교화하기 위해 향후 연구에서 더 밀도 높은 지역 샘격, 시간적 데이터, 그리고 더 정확한 이동성 지표를 포함할 것을 권고한다.