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

Sintesi Tecnica: Caratterizzazione dell'intero genoma delle infezioni da Plasmodium vivax in viaggiatori locali e non viaggiatori nella foresta amazzonica peruviana

Definizione del problema
Mentre le regioni progrediscono verso l'eliminazione della malaria, distinguere tra infezioni da Plasmodium vivax acquisite localmente e importate è fondamentale per una sorveglianza mirata. Tuttavia, in contesti caratterizzati da alta mobilità umana, bassa diversità genetica dei parassiti e colli di bottiglia storici della popolazione — come la foresta amazzonica peruviana — distinguere l'origine delle infezioni a scale spaziali fini rimane difficile. Gli strumenti molecolari tradizionali (es. microsatelliti, barcode SNP) possono mancare della risoluzione necessaria per differenziare comunità vicine o per identificare segnali recenti di importazione quando le popolazioni di parassiti presentano un'alta parentela. Inoltre, la biologia di P. vivax, specificamente il potenziale di ricaduta da ipnozoiti dormienti, complica l'interpretazione delle infezioni rilevate dopo un viaggio, poiché queste potrebbero rappresentare riattivazioni di precedenti infezioni locali piuttosto che nuove importazioni. Questo studio affronta il vuoto nella comprensione se il sequenziamento dell'intero genoma (WGS) possa distinguere efficacemente le infezioni associate ai viaggiatori dalle trasmissioni locali in un contesto di alta somiglianza genetica e dati di mobilità incompleti.

Metodologia
Lo studio ha utilizzato un disegno trasversale che ha coinvolto 56 infezioni confermate da P. vivax raccolte tra il 2021 e il 2022 in tre comunità del distretto di Mazán, regione di Loreto: Libertad (n=24), Gamitanacocha (n=23) e Urco Miraño (n=9). I partecipanti sono stati classificati come "viaggiatori" se avevano riferito di aver soggiornato fuori dalla propria comunità per almeno una notte nell'ultimo mese (n=16) o come "non viaggiatori" (n=40).

• Processamento dei campioni: Il DNA è stato estratto da sangue intero e la specie del parassita è stata confermata tramite qPCR bersagliando il gene mitocondriale cox1. È stata impiegata l'amplificazione selettiva dell'intero genoma (sWGA) per arricchire il DNA del parassita prima del sequenziamento.
• Sequenziamento: Le librerie sono state preparate utilizzando un protocollo PCR-free e sequenziate sulla piattaforma Illumina NovaSeq, generando read appaiate (paired-end) di 150 bp con una copertura media di ~100×.
• Analisi dei dati:
  • Variant Calling: Le read sono state allineate al riferimento di P. vivax (PvPAM) e le varianti sono state identificate utilizzando GATK. È stato mantenuto un set finale di 16.485 SNP biallelici di alta qualità per l'analisi.
  • Diversità genetica: L'eterozigosi attesa ($H_e$) e la diversità nucleotidica ($\pi$) sono state calcolate per valutare la diversità all'interno e tra le comunità e i gruppi di viaggio.
  • Struttura della popolazione: L'analisi delle componenti principali (PCA), l'analisi discriminante delle componenti principali (DAPC) e ADMIXTURE sono stati utilizzati per identificare il raggruppamento genetico e le proporzioni di ascendenza.
  • Parentela: L'analisi dell'identità per discendenza (IBD) è stata eseguita utilizzando il pacchetto isoRelate per rilevare la discendenza condivisa e i recenti collegamenti di trasmissione.
  • Analisi di rete: Alberi di espansione minima (MST) e clustering gerarchico sono stati generati per visualizzare le distanze genetiche e la connettività.

Risultati chiave

• Diversità genetica: Libertad ha esibito la più alta diversità genetica complessiva ($H_e$ = 0,363), seguita da Urco Miraño (0,198) e Gamitanacocha (0,1). All'interno di Libertad, le infezioni associate ai viaggiatori hanno mostrato una diversità significativamente maggiore rispetto ai non viaggiatori ($H_e$ = 0,432 vs. 0,304; $\pi$ = 0,290 vs. 0,168).
• Differenziazione della popolazione: L'analisi $F_{st}$ pairwise ha rivelato una bassa differenziazione tra Libertad e Gamitanacocha ($F_{st}$ = 0,056), suggerendo una rete di trasmissione condivisa. Al contrario, Urco Miraño era altamente differenziata sia da Libertad ($F_{st}$ = 0,548) che da Gamitanacocha ($F_{st}$ = 0,773).
• Clustering e struttura: PCA e DAPC hanno confermato che i parassiti di Urco Miraño formano un cluster distinto, mentre Libertad e Gamitanacocha mostrano una parziale sovrapposizione. Le infezioni associate ai viaggiatori non hanno formato un cluster globale distinto; invece, sono rimaste ampiamente connesse e intersporse nei cluster della loro comunità di residenza.
• Ascendenza e parentela: ADMIXTURE ha indicato una maggiore eterogeneità tra i viaggiatori in Libertad, con alcuni che mostravano un'ascendenza mista. L'analisi IBD a soglie elevate (50–80%) ha identificato due infezioni associate ai viaggiatori (L1 da Libertad e U1 da Urco Miraño) come singletons geneticamente distinti, potenzialmente rappresentativi di linee evolutive provenienti da località non campionate. Tuttavia, la maggior parte delle infezioni dei viaggiatori era inserita all'interno delle reti di trasmissione locale.
• Connettività: La rete era altamente connessa a soglie IBD inferiori, con una forte parentela osservata all'interno di Gamitanacocha e tra Libertad e Gamitanacocha. Nessuna coppia correlata è stata rilevata tra Urco Miraño e le altre due comunità.

Significato e affermazioni
Gli autori affermano che questo studio dimostra il potenziale e i limiti del WGS per distinguere le infezioni da P. vivax associate ai viaggiatori da quelle acquisite localmente a scale spaziali fini. I risultati supportano l'integrazione di dati genomici ed epidemiologici per interpretare la trasmissione associata alla mobilità. Nello specifico, lo studio conclude che:

1. Il WGS può risolvere la struttura a livello di comunità: Ha identificato con successo cluster genetici distinti in comunità isolate (Urco Miraño) rispetto a quelle connesse (Libertad e Gamitanacocha).
2. I rapporti di viaggio non equivalgono sempre a un'importazione: La maggior parte delle infezioni associate ai viaggiatori si è raggruppata con i parassiti locali, suggerendo che il viaggio riportato non ha necessariamente indicato un'importazione; molte infezioni erano probabilmente acquisite all'interno della rete di trasmissione locale o tramite ricaduta.
3. Risoluzione dipendente dal contesto: La capacità di distinguere le origini delle infezioni dipende fortemente dal livello di differenziazione genetica e dalla disponibilità di dati di riferimento rappresentativi. In contesti con bassa diversità e alta parentela, il solo WGS potrebbe non essere sufficiente per assegnare l'origine senza un campionamento denso e un robusto contesto epidemiologico.
4. Identificazione di outlier: Sebbene la maggior parte dei viaggiatori si sia raggruppata localmente, il WGS ha comunque identificato specifiche infezioni associate ai viaggiatori con profili genetici distinti, suggerendo l'utilità del metodo nel rilevare linee evolutive potenzialmente associate alla mobilità quando esiste una sufficiente differenziazione.

Il documento sottolinea che questi risultati devono essere interpretati con cautela a causa delle piccole dimensioni del campione, in particolare per i gruppi di viaggiatori, e della complessità intrinseca della biologia di P. vivax (ipnozoiti). Gli autori raccomandano studi futuri che incorporino un campionamento regionale più denso, dati temporali e metriche di mobilità più accurate per affinare le strategie di sorveglianza genomica per l'eliminazione della malaria.