[Virus] 말라리아 분자 감시를 위한 선택적 전장유전체 증폭(sWGA) 전략 개선 연구 정리

Enhancing whole genome DNA amplification of Plasmodium falciparum for advanced molecular surveillance in malaria control

 

저농도(저기생충혈증) 임상 샘플에서 사람 DNA가 압도적으로 많은 상황에서도, sWGA(SWGA) 조건을 임상 워크플로에 맞게 단축하고(16h→8h), 더 나아가 EquiPhi29 기반 “3시간 프로토콜”로 당일 시퀀싱이 가능한 수준까지 끌어올린 최적화 연구입니다(비용도 크게 절감).

---

왜 이 연구가 필요한가? (배경)

말라리아 분자감시는 약제내성(예: pfcrt, pfmdr1, pfdhps/pfdhfr, pfkelch13) 추적, 유행 집단의 이동/전파 분석에 핵심인데, 실제 현장에서는 무증상/저농도 감염 샘플이 많아 기생충 DNA가 너무 적고 사람 DNA가 너무 많아 WGS가 어렵다는 문제가 큽니다. 그래서 Selective Whole Genome Amplification(SWGA)가 등장합니다. 짧은 올리고 프라이머가 Plasmodium DNA에 더 자주/선호적으로 결합하도록 설계하고, phi29 기반 MDA로 기생충 DNA를 선택적으로 증폭해 WGS를 가능하게 합니다. 다만 기존 프로토콜은 16시간 이상 등 시간이 길어 임상/루틴 환경에 부담이 큽니다.

---

연구 목표 (무엇을 최적화했나)

1. 기존 SWGA(step-down + 30°C 장시간 반응)에서 반응 시간을 줄여도 WGS 품질(커버리지, 변이콜, 내성 프로파일)이 유지되는지 평가
2. 더 빠른 효소(EquiPhi29)를 활용한 “3시간 SWGA”가 실제로 커버리지/깊이/비용에서 이점이 있는지 검증

---

샘플 구성 (임상 현실 반영)

• 총 13개 P. falciparum 양성 DNA (UKHSA Malaria Reference Lab 제공)
• 여행자(2019–2020) 유래: 서아프리카 8, 중/동아프리카 5
• qPCR Ct 범위 13.01–36.43, 이 중 7개는 Ct >35(매우 낮은 기생충 DNA)

---

전체 실험 워크플로 (한눈에 보기)

DNA 추출/정량(qPCR 포함) → SWGA(조건별) → bead 정제 → (분지 DNA 처리) T7 Endonuclease 처리 → ONT 라이브러리(바코딩) → MinION 시퀀싱(R10.4.1) → 맵핑/분류/품질 → 변이콜(FreeBayes) → 약제내성 프로파일(Malaria-Profiler)

---

실험 방법

1) DNA 추출 및 정량

• Whole blood에서 QIAsymphony SP + DSP DNA kit(Qiagen)로 추출
• 사용 전 Qubit HS로 dsDNA 정량 + qPCR로 P. falciparum 양/ Ct 확인

---

2) SWGA 조건 1: “기존(phi29, NEB) + step-down 프로토콜” (2h / 8h / 16h 비교)

• 입력 DNA: gDNA 120 ng, 물로 최대 30 µL로 맞춤
• 50 µL 반응 조성(핵심 구성만):
  • 10× phi29 buffer 5 µL
  • BSA 0.35 µL
  • P. falciparum SWGA primer mix (250 µM) 0.5 µL
  • dNTP(10 mM) 5 µL
  • phi29 DNA polymerase 30U (3 µL)
• step-down 온도: 35°C 5min → 34°C 10min → 33°C 15min → 32°C 20min → 31°C 30min → 30°C에서 2h 또는 8h 또는 16h
• 반응 종료: 65°C 15min 불활성화
• 정제: KAPA beads 1:1 비율 bead clean-up, EB 31 µL 용출
• 시약비(이 SWGA 반응만): £6.59 / sample

포인트: 2h는 “당일 처리”를 노린 조건, 8h는 “overnight”을 노린 조건으로 설정했다는 점이 임상/루틴 감시 목적에 현실적입니다.

---

3) SWGA 조건 2: “EquiPhi29 기반 3시간 프로토콜” (신규/단축)

• 대상: 3개 샘플(C1, IC1, M1; 서로 다른 기생충 농도 대표)
• 핵심 흐름:
  1. gDNA 120 ng + primer mix 포함(총 5 µL)로 만든 뒤
  2. 95°C 3min 변성 → 얼음 5min
  3. EquiPhi29 반응 mix 추가 후
  4. 45°C 3h 반응 → 65°C 10min 불활성화
• 정제/정량: 기존과 동일하게 bead-wash 후 EB 31 µL 용출
• 시약비(이 SWGA 반응만): £2.28 / sample

---

4) 라이브러리 준비 및 ONT WGS

• SWGA 산물은 가지친/분지 구조가 생길 수 있어, 라이브러리 전 T7 Endonuclease I 처리로 linearize/debranch
  • 예: DNA 25.5 µL + NEBuffer2 3 µL + T7 endonuclease I 1.5 µL
  • 37°C 15min → KAPA beads 0.6× 정제 → NFW 27 µL 용출(50°C 1min + RT 5min)
• 라이브러리: Native Barcoding kit 24 V14 (SQK-NBD114.24)
• 시퀀싱:
  • MinION Mk1B + R10.4.1 flow cell
  • 라이브러리 20 fmol 로딩, 22h 런
  • demux: MinKNOW v23.07.12 (run 중)

---

데이터 분석 파이프라인 (실제 운영 관점에서 중요한 부분)

• basecalling: Guppy v6.5.7
• mapping: minimap2 v2.26 (레퍼런스: Pf3D7 v2015-06-18, PlasmoDB)
• taxonomy: kraken2 v2.1.3
• QC: pycoQC v2.5.2
• coverage: samtools depth
• variant calling: FreeBayes v1.3.2
  • alt allele count ≥10, depth ≥10, baseQ ≥10, clumping off
  • diploid model로 혼합감염(다클론) 가능성 고려
• 읽기 수 차이 보정: Rasusa v2.0.0로 fastq에서 random subsampling하여 read 수 정규화
• 일치도 평가: Cohen’s Kappa + SNP concordance
• 약제내성: Malaria-Profiler(온라인)
  • ONT 옵션, soft cutoff depth 5×, MAF 50% 등 사용
  • 불일치 시 IGV로 수동 확인

---

핵심 결과 정리

1) 16시간을 8시간으로 줄여도 “거의 같다”

• 8h vs 16h는 사람/플라스모디움 맵핑 비율 평균 차이가 매우 작았고, 커버리지도 read 수 정규화 후 거의 동일한 수준으로 수렴
• 2h는 플라스모디움 read 비율(~25.8%) 등에서 성능이 떨어져 이후 분석에서 제외

2) 변이콜(SNP) 관점에서도 8h vs 16h는 높은 일치도

• 평균 SNP concordance 87.45–96.41%, Kappa 0.67–1.00
• 깊이가 올라갈수록 concordance가 좋아져 20× 이상에서 대부분 90% 초과

3) “3시간 EquiPhi29 프로토콜”은 커버리지/깊이/속도/비용에서 확실히 유리

• 3시간 반응으로 당일 처리/시퀀싱 가능
• (3개 대표 샘플에서) 커버리지 최대 +28.8%, depth 2.7–5.3배 증가
• 논문 Discussion에서 요약한 성과: 핵 유전체의 ≥10 reads 커버리지 기준 최소 63.4% 확보, 그리고 비용은 기존 대비 “1/3 조금 넘는 수준”

4) 약제내성 프로파일은 조건 바꿔도 대부분 유지

• 32개 조건 비교 중 29개에서 일치, 불일치 3개는 “돌연변이 부재”가 아니라 해당 구간 커버리지 부족 때문

---

실무 적용 관점의 Take-home message

• 루틴 감시(overnight 가능)라면: 기존 프로토콜을 그대로 쓰되 16h→8h로 단축해도 결과 신뢰성이 크게 흔들리지 않을 가능성이 큼(커버리지/변이콜 일치도 관점).
• “당일 turnaround”가 목표라면: EquiPhi29 3h + T7 endonuclease 처리 + ONT 바코딩/MinION 조합이 현실적인 대안(비용도 낮음).
• 단, ONT 특성상(특히 homopolymer) 변이 불일치가 일부 생길 수 있어 충분한 depth 확보(예: ≥20×)가 중요하다는 점을 논문도 강조합니다.

---

마무리

이 연구는 “저농도 말라리아 샘플에서 WGS를 어떻게 루틴으로 돌릴 것인가?”라는 질문에 대해, 시간(16h→8h→3h)과 비용을 줄이면서도 변이/내성 해석의 일관성을 유지하는 방향을 제시한다는 점에서 의미가 큽니다. 특히 3시간 SWGA는 감시 목적뿐 아니라, 향후 임상 의사결정까지 연결될 수 있는 turnaround 혁신의 후보로 보입니다.

https://link.springer.com/article/10.1186/s12936-025-05643-9