Monitoring Plasmodium falciparum and Plasmodium vivax using microsatellite markers indicates limited changes in population structure after substantial transmission decline in Papua New Guinea

Tracing the origins of Plasmodium vivax resurgence after malaria elimination on Aneityum Island in Vanuatu

BACKGROUND

Five years after successful malaria elimination, Aneityum Island in Vanuatu experienced an outbreak of Plasmodium vivax of unknown origin in 2002. Epidemiological investigations revealed several potential sources of P. vivax. We aimed to identify the genetic origin of P. vivax responsible for the resurgence.

METHODS

Five P. vivax microsatellite markers were genotyped using DNA extracted from archived blood samples. A total of 69 samples from four P. vivax populations was included: 29 from the outbreak in 2002, seven from Aneityum in 1999 and 2000, 18 from visitors to Aneityum in 2000, and 15 from nearby Tanna Island in 2002. A neighbour-joining phylogenetic tree was constructed to elucidate the relationships among P. vivax isolates. STRUCTURE and principal component analysis were used to assess patterns of genetic structure.

RESULTS

Here we show distinct genetic origins of P. vivax during the outbreak on Aneityum. While the origin of most P. vivax lineages found during the outbreak remains unidentified, limited genetic diversity among these lineages is consistent with a rapid expansion from a recent common ancestor. Contemporaneous P. vivax from neighboring Tanna and potential relapse of P. vivax acquired from other islands in 1999 and 2000 are also identified as minor contributors to the outbreak.

CONCLUSIONS

Multiple reintroductions of P. vivax after elimination highlight the high receptivity and vulnerability to malaria resurgence in island settings of Vanuatu, despite robust surveillance and high community compliance to control measures.

Imported malaria into Australia: surveillance insights and opportunities

BACKGROUND

Malaria continues to pose a significant burden in endemic countries, many of which lack access to molecular surveillance. Insights from malaria cases in travellers returning to non-endemic areas can provide valuable data to inform endemic country programmes. To evaluate the potential for novel global insights into malaria, we examined epidemiological and molecular data from imported malaria cases to Australia.

METHODS

We analysed malaria cases reported in Australia from 2012 to 2022 using National Notifiable Disease Surveillance System data. Molecular data on imported malaria cases were obtained from literature searches.

RESULTS

Between 2012 and 2022, 3204 malaria cases were reported in Australia. Most cases (69%) were male and 44% occurred in young adults aged 20-39 years. Incidence rates initially declined between 2012 and 2015, then increased until 2019. During 2012-2019, the incidence in travellers ranged from 1.34 to 7.71 per 100 000 trips. Cases were primarily acquired in Sub-Saharan Africa (n = 1433; 45%), Oceania (n = 569; 18%) and Southern and Central Asia (n = 367; 12%). The most common countries of acquisition were Papua New Guinea (n = 474) and India (n = 277). Plasmodium falciparum accounted for 58% (1871/3204) of cases and was predominantly acquired in Sub-Saharan Africa, and Plasmodium vivax accounted for 32% (1016/3204), predominantly from Oceania and Asia. Molecular studies of imported malaria cases to Australia identified genetic mutations and deletions associated with drug resistance and false-negative rapid diagnostic test results, and led to the establishment of reference genomes for P. vivax and Plasmodium malariae.

CONCLUSIONS

Our analysis highlights the continuing burden of imported malaria into Australia. Molecular studies have offered valuable insights into drug resistance and diagnostic limitations, and established reference genomes. Integrating molecular data into national surveillance systems could provide important infectious disease intelligence to optimize treatment guidelines for returning travellers and support endemic country surveillance programmes.

Microsatellites reveal high polymorphism and high potential for use in anti-malarial efficacy studies in areas with different transmission intensities in mainland Tanzania

BACKGROUND

Tanzania is currently implementing therapeutic efficacy studies (TES) in areas of varying malaria transmission intensities as per the World Health Organization (WHO) recommendations. In TES, distinguishing reinfection from recrudescence is critical for the determination of anti-malarial efficacy. Recently, the WHO recommended genotyping polymorphic coding genes, merozoite surface proteins 1 and 2 (msp1 and msp2), and replacing the glutamate-rich protein (glurp) gene with one of the highly polymorphic microsatellites in Plasmodium falciparum to adjust the efficacy of antimalarials in TES. This study assessed the polymorphisms of six neutral microsatellite markers and their potential use in TES, which is routinely performed in Tanzania.

METHODS

Plasmodium falciparum samples were obtained from four TES sentinel sites, Kibaha (Pwani), Mkuzi (Tanga), Mlimba (Morogoro) and Ujiji (Kigoma), between April and September 2016. Parasite genomic DNA was extracted from dried blood spots on filter papers using commercial kits. Genotyping was done using six microsatellites (Poly-α, PfPK2, TA1, C3M69, C2M34 and M2490) by capillary method, and the data were analysed to determine the extent of their polymorphisms and genetic diversity at the four sites.

RESULTS

Overall, 83 (88.3%) of the 94 samples were successfully genotyped (with positive results for ≥ 50.0% of the markers), and > 50.0% of the samples (range = 47.6-59.1%) were polyclonal, with a mean multiplicity of infection (MOI) ranging from 1.68 to 1.88 among the four sites. There was high genetic diversity but limited variability among the four sites based on mean allelic richness (RS = 7.48, range = 7.27-8.03, for an adjusted minimum sample size of 18 per site) and mean expected heterozygosity (He = 0.83, range = 0.80-0.85). Cluster analysis of haplotypes using STRUCTURE, principal component analysis, and pairwise genetic differentiation (FST) did not reveal population structure or clustering of parasites according to geographic origin. Of the six markers, Poly-α was the most polymorphic, followed by C2M34, TA1 and C3M69, while M2490 was the least polymorphic.

CONCLUSION

Microsatellite genotyping revealed high polyclonality and genetic diversity but no significant population structure. Poly-α, C2M34, TA1 and C3M69 were the most polymorphic markers, and Poly-α alone or with any of the other three markers could be adopted for use in TES in Tanzania.

Population genomic evidence of structured and connected Plasmodium vivax populations under host selection in Latin America

Pathogen genomic epidemiology has the potential to provide a deep understanding of population dynamics, facilitating strategic planning of interventions, monitoring their impact, and enabling timely responses, and thereby supporting control and elimination efforts of parasitic tropical diseases. Plasmodium vivax, responsible for most malaria cases outside Africa, shows high genetic diversity at the population level, driven by factors like sub-patent infections, a hidden reservoir of hypnozoites, and early transmission to mosquitoes. While Latin America has made significant progress in controlling Plasmodium falciparum, it faces challenges with residual P. vivax. To characterize genetic diversity and population structure and dynamics, we have analyzed the largest collection of P. vivax genomes to date, including 1474 high-quality genomes from 31 countries across Asia, Africa, Oceania, and America. While P. vivax shows high genetic diversity globally, Latin American isolates form a distinctive population, which is further divided into sub-populations and occasional clonal pockets. Genetic diversity within the continent was associated with the intensity of transmission. Population differentiation exists between Central America and the North Coast of South America, vs. the Amazon Basin, with significant gene flow within the Amazon Basin, but limited connectivity between the Northwest Coast and the Amazon Basin. Shared genomic regions in these parasite populations indicate adaptive evolution, particularly in genes related to DNA replication, RNA processing, invasion, and motility - crucial for the parasite's survival in diverse environments. Understanding these population-level adaptations is crucial for effective control efforts, offering insights into potential mechanisms behind drug resistance, immune evasion, and transmission dynamics.

Molecular markers for malaria genetic epidemiology: progress and pitfalls

Over recent years, progress in molecular markers for genotyping malaria parasites has enabled informative studies of epidemiology and transmission dynamics. Results have highlighted the value of these tools for surveillance to support malaria control and elimination strategies. There are many different types and panels of markers available for malaria parasite genotyping, and for end users, the nuances of these markers with respect to 'use case', resolution, and accuracy, are not well defined. This review clarifies issues surrounding different molecular markers and their application to malaria control and elimination. We describe available marker panels, use cases, implications for different transmission settings, limitations, access, cost, and data accuracy. The information provided can be used as a guide for molecular epidemiology and surveillance of malaria.

Temporal and spatial analysis of Plasmodium falciparum genomics reveals patterns of parasite connectivity in a low-transmission district in Southern Province, Zambia

BACKGROUND

Understanding temporal and spatial dynamics of malaria transmission will help to inform effective interventions and strategies in regions approaching elimination. Parasite genomics are increasingly used to monitor epidemiologic trends, including assessing residual transmission across seasons and importation of malaria into these regions.

METHODS

In a low and seasonal transmission setting of southern Zambia, a total of 441 Plasmodium falciparum samples collected from 8 neighbouring health centres between 2012 and 2018 were genotyped using molecular inversion probes (MIPs n = 1793) targeting a total of 1832 neutral and geographically informative SNPs distributed across the parasite genome. After filtering for quality and missingness, 302 samples and 1410 SNPs were retained and used for downstream population genomic analyses.

RESULTS

The analyses revealed most (67%, n = 202) infections harboured one clone (monogenomic) with some variation at local level suggesting low, but heterogenous malaria transmission. Relatedness identity-by-descent (IBD) analysis revealed variable distribution of IBD segments across the genome and 6% of pairs were highly-related (IBD ≥ 0.25). Some of the highly-related parasite populations persisted across multiple seasons, suggesting that persistence of malaria in this low-transmission region is fueled by parasites "seeding" across the dry season. For recent years, clusters of clonal parasites were identified that were dissimilar to the general parasite population, suggesting parasite populations were increasingly fragmented at small spatial scales due to intensified control efforts. Clustering analysis using PCA and t-SNE showed a lack of substantial parasite population structure.

CONCLUSION

Leveraging both genomic and epidemiological data provided comprehensive picture of fluctuations in parasite populations in this pre-elimination setting of southern Zambia over 7 years.

Asia-Pacific ICEMR: Understanding Malaria Transmission to Accelerate Malaria Elimination in the Asia Pacific Region

Gaining an in-depth understanding of malaria transmission requires integrated, multifaceted research approaches. The Asia-Pacific International Center of Excellence in Malaria Research (ICEMR) is applying specifically developed molecular and immunological assays, in-depth entomological assessments, and advanced statistical and mathematical modeling approaches to a rich series of longitudinal cohort and cross-sectional studies in Papua New Guinea and Cambodia. This is revealing both the essential contribution of forest-based transmission and the particular challenges posed by Plasmodium vivax to malaria elimination in Cambodia. In Papua New Guinea, these studies document the complex host-vector-parasite interactions that are underlying both the stunning reductions in malaria burden from 2006 to 2014 and the significant resurgence in transmission in 2016 to 2018. Here we describe the novel analytical, surveillance, molecular, and immunological tools that are being applied in our ongoing Asia-Pacific ICEMR research program.

The impact of indoor residual spraying on Plasmodium falciparum microsatellite variation in an area of high seasonal malaria transmission in Ghana, West Africa

Here, we report the first population genetic study to examine the impact of indoor residual spraying (IRS) on Plasmodium falciparum in humans. This study was conducted in an area of high seasonal malaria transmission in Bongo District, Ghana. IRS was implemented during the dry season (November-May) in three consecutive years between 2013 and 2015 to reduce transmission and attempt to bottleneck the parasite population in humans towards lower diversity with greater linkage disequilibrium. The study was done against a background of widespread use of long-lasting insecticidal nets, typical for contemporary malaria control in West Africa. Microsatellite genotyping with 10 loci was used to construct 392 P. falciparum multilocus infection haplotypes collected from two age-stratified cross-sectional surveys at the end of the wet seasons pre- and post-IRS. Three-rounds of IRS, under operational conditions, led to a >90% reduction in transmission intensity and a 35.7% reduction in the P. falciparum prevalence (p < .001). Despite these declines, population genetic analysis of the infection haplotypes revealed no dramatic changes with only a slight, but significant increase in genetic diversity (He : pre-IRS = 0.79 vs. post-IRS = 0.81, p = .048). Reduced relatedness of the parasite population (p < .001) was observed post-IRS, probably due to decreased opportunities for outcrossing. Spatiotemporal genetic differentiation between the pre- and post-IRS surveys (D = 0.0329 [95% CI: 0.0209 - 0.0473], p = .034) was identified. These data provide a genetic explanation for the resilience of P. falciparum to short-term IRS programmes in high-transmission settings in sub-Saharan Africa.

The changing epidemiology of Plasmodium vivax: Insights from conventional and novel surveillance tools

Sarah Auburn and co-authors discuss the unique biology and epidemiology of P. vivax and current evidence on conventional and new approaches to surveillance.

SNP barcodes provide higher resolution than microsatellite markers to measure Plasmodium vivax population genetics

BACKGROUND

Genomic surveillance of malaria parasite populations has the potential to inform control strategies and to monitor the impact of interventions. Barcodes comprising large numbers of single nucleotide polymorphism (SNP) markers are accurate and efficient genotyping tools, however may need to be tailored to specific malaria transmission settings, since 'universal' barcodes can lack resolution at the local scale. A SNP barcode was developed that captures the diversity and structure of Plasmodium vivax populations of Papua New Guinea (PNG) for research and surveillance.

METHODS

Using 20 high-quality P. vivax genome sequences from PNG, a total of 178 evenly spaced neutral SNPs were selected for development of an amplicon sequencing assay combining a series of multiplex PCRs and sequencing on the Illumina MiSeq platform. For initial testing, 20 SNPs were amplified in a small number of mono- and polyclonal P. vivax infections. The full barcode was then validated by genotyping and population genetic analyses of 94 P. vivax isolates collected between 2012 and 2014 from four distinct catchment areas on the highly endemic north coast of PNG. Diversity and population structure determined from the SNP barcode data was then benchmarked against that of ten microsatellite markers used in previous population genetics studies.

RESULTS

From a total of 28,934,460 reads generated from the MiSeq Illumina run, 87% mapped to the PvSalI reference genome with deep coverage (median = 563, range 56-7586) per locus across genotyped samples. Of 178 SNPs assayed, 146 produced high-quality genotypes (minimum coverage = 56X) in more than 85% of P. vivax isolates. No amplification bias was introduced due to either polyclonal infection or whole genome amplification (WGA) of samples before genotyping. Compared to the microsatellite panels, the SNP barcode revealed greater variability in genetic diversity between populations and geographical population structure. The SNP barcode also enabled assignment of genotypes according to their geographic origins with a significant association between genetic distance and geographic distance at the sub-provincial level.

CONCLUSIONS

High-throughput SNP barcoding can be used to map variation of malaria transmission dynamics at sub-national resolution. The low cost per sample and genotyping strategy makes the transfer of this technology to field settings highly feasible.