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

Resumen Técnico: Caracterización de todo el genoma de las infecciones por Plasmodium vivax en viajeros locales y no viajeros de la Amazonía peruana

Planteamiento del problema
A medida que las regiones avanzan hacia la eliminación de la malaria, distinguir entre las infecciones por Plasmodium vivax adquiridas localmente e importadas es fundamental para una vigilancia dirigida. Sin embargo, en entornos caracterizados por una alta movilidad humana, baja diversidad genética del parásito y cuellos de botella poblacionales históricos —como la Amazonía peruana—, diferenciar los orígenes de la infección a escalas espaciales finas sigue siendo un desafío. Las herramientas moleculares tradicionales (por ejemplo, microsatélites, códigos de barras de SNP) pueden carecer de la resolución para diferenciar entre comunidades vecinas o para identificar señales recientes de importación cuando las poblaciones de parásitos exhiben una alta parentela. Además, la biología de P. vivax, específicamente el potencial de recaída a partir de hipnozoitos latentes, complica la interpretación de las infecciones detectadas tras viajar, ya que estas pueden representar reactivaciones de infecciones locales previas en lugar de nuevas importaciones. Este estudio aborda la brecha en la comprensión de si la Secuenciación de Todo el Genoma (WGS, por sus siglas en inglés) puede distinguir eficazmente las infecciones asociadas a viajeros de la transmisión local en un contexto de alta similitud genética e información de movilidad incompleta.

Metodología
El estudio utilizó un diseño transversal que involucró 56 infecciones confirmadas por P. vivax recolectadas entre 2021 y 2022 de tres comunidades en el distrito de Mazán, región de Loreto: Libertad (n=24), Gamitanacocha (n=23) y Urco Miraño (n=9). Los participantes fueron clasificados como "viajeros" si informaron haber permanecido fuera de su comunidad durante al menos una noche en el último mes (n=16) o como "no viajeros" (n=40).

• Procesamiento de muestras: El ADN se extrajo de sangre total, y la especie del parásito se confirmó mediante qPCR dirigido al gen mitocondrial cox1. Se empleó la Amplificación de Genoma Completo selectiva (sWGA) para enriquecer el ADN del parásito antes de la secuenciación.
• Secuenciación: Las librerías se prepararon utilizando un protocolo libre de PCR y se secuenciaron en la plataforma Illumina NovaSeq, generando lecturas de pares de extremos (paired-end) de 150 pb con una cobertura promedio de ~100×.
• Análisis de datos:
  • Identificación de variantes (Variant Calling): Las lecturas se alinearon con la referencia de P. vivax (PvPAM), y las variantes se identificaron utilizando GATK. Se retuvo un conjunto final de 16,485 SNPs bialélicos de alta calidad para el análisis.
  • Diversidad genética: Se calcularon la heterocigosidad esperada ($H_e$) y la diversidad de nucleótidos ($\pi$) para evaluar la diversidad dentro y entre comunidades y grupos de viaje.
  • Estructura poblacional: Se utilizaron el Análisis de Componentes Principales (PCA), el Análisis Discriminante de Componentes Principales (DAPC) y ADMIXTURE para identificar el agrupamiento genético y las proporciones de ancestría.
  • Parentesco (Relatedness): Se realizó un análisis de Identidad por Descenso (IBD) utilizando el paquete isoRelate para detectar ancestría compartida y vínculos de transmisión reciente.
  • Análisis de redes: Se generaron Árboles de Expansión Mínima (MST) y agrupamientos jerárquicos para visualizar las distancias genéticas y la conectividad.

Resultos clave

• Diversidad genética: Libertad exhibió la mayor diversidad genética general ($H_e$ = 0.363), seguida de Urco Miraño (0.198) y Gamitanacocha (0.1). Dentro de Libertad, las infecciones asociadas a viajeros mostraron una diversidad significativamente mayor que las de no viajeros ($H_e$ = 0.432 vs. 0.304; $\pi$ = 0.290 vs. 0.168).
• Diferenciación poblacional: El análisis de $F_{st}$ por pares reveló una baja diferenciación entre Libertad y Gamitanacocha ($F_{st}$ = 0.056), lo que sugiere una red de transmisión compartida. En contraste, Urco Miraño estuvo altamente diferenciado tanto de Libertad ($F_{st}$ = 0.548) como de Gamitanacocha ($F_{st}$ = 0.773).
• Agrupamiento y estructura: El PCA y el DAPC confirmaron que los parásitos de Urco Miraño formaron un grupo distinto, mientras que los de Libertad y Gamitanacocha mostraron un solapamiento parcial. Las infecciones asociadas a viajeros no formaron un grupo global distinto; en cambio, permanecieron mayoritariamente conectadas e interpuestas dentro de los grupos de su comunidad de residencia.
• Ancestría y parentesco: ADMIXTURE indicó una mayor heterogeneidad entre los viajeros en Libertad, donde algunos mostraron una ancestría mixta. El análisis de IBD en umbrales altos (50–80%) identificó dos infecciones asociadas a viajeros (L1 de Libertad y U1 de Urco Miraño) como linajes únicos (singletons) genéticamente distintos, que potencialmente representan linajes de ubicaciones no muestreadas. Sin embargo, la mayoría de las infecciones de viajeros estaban integradas en las redes de transmisión local.
• Conectividad: La red estaba altamente conectada en umbrales de IBD más bajos, con una fuerte parentela observada dentro de Gamitanacocha y entre Libertad y Gamitanacocha. No se detectaron pares relacionados entre Urco Miraño y las otras dos comunidades.

Significancia y afirmaciones
Los autores afirman que este estudio demuestra el potencial y las limitaciones de la WGS para distinguir las infecciones asociadas a viajeros de las adquiridas localmente por P. vivax a escalas espaciales finas. Los hallazgos respaldan la integración de datos genómicos y epidemiológicos para interpretar la transmisión asociada a la movilidad. Específicamente, el estudio concluye que:

1. La WGS puede resolver la estructura a nivel de comunidad: Identificó con éxito grupos genéticos distintos en comunidades aisladas (Urco Miraño) frente a las conectadas (Libertad y Gamitanacocha).
2. Los reportes de viajes no siempre equivalen a importación: La mayoría de las infecciones asociadas a viajeros se agruparon con parásitos locales, lo que sugiere que el viaje reportado no necesariamente indicó una importación; muchas de estas infecciones fueron probablemente adquiridas dentro de la red de transmisión local o vía recaída.
3. Resolución dependiente del contexto: La capacidad de distinguir los orígenes de la infección depende de la diferenciación genética y de la disponibilidad de datos de referencia representativos. En entornos con baja diversidad y alta parentela, la WGS por sí sola puede no ser suficiente para asignar el origen sin un muestreo denso y un contexto epidemiológico robusto.
4. Identificación de valores atípicos (outliers): Aunque la mayoría de los viajeros se agruparon localmente, la WGS también identificó infecciones asociadas a viajeros con perfiles genéticos distintos, lo que sugiere la utilidad del método para detectar linajes potencialmente asociados a la movilidad cuando existe suficiente diferenciación.

El artículo enfatiza que estos hallazgos deben interpretarse con cautela debido al pequeño tamaño de la muestra, particularmente para los grupos de viajeros, y la complejidad inherente de la biología de P. vivax (hipnozoitos). Los autores abogan por futuros estudios que incorporen un muestreo regional más denso, datos temporales y métricas de movilidad más precisas para refinar las estrategias de vigilancia genómica para la eliminación de la malaria.