Plasmodium vivax genomic surveillance in the Peruvian Amazon with Pv AmpliSeq assay.

BACKGROUND: Plasmodium vivax is the most predominant malaria species in Latin America, constituting 71.5% of malaria cases in 2021. With several countries aiming for malaria elimination, it is crucial to prioritize effectiveness of national control programs by optimizing the utilization of available resources and strategically implementing necessary changes. To support this, there is a need for innovative approaches such as genomic surveillance tools that can investigate changes in transmission intensity, imported cases and sources of reintroduction, and can detect molecular markers associated with drug resistance. METHODOLOGY/PRINCIPAL FINDINGS: Here, we apply a modified highly-multiplexed deep sequencing assay: Pv AmpliSeq v2 Peru. The tool targets a newly developed 41-SNP Peru barcode for parasite population analysis within Peru, the 33-SNP vivaxGEN-geo panel for country-level classification, and 11 putative drug resistance genes. It was applied to 230 samples from the Peruvian Amazon (2007-2020), generating baseline surveillance data. We observed a heterogenous P. vivax population with high diversity and gene flow in peri-urban areas of Maynas province (Loreto region) with a temporal drift using all SNPs detected by the assay (nSNP = 2909). In comparison, in an indigenous isolated area, the parasite population was genetically differentiated (FST = 0.07-0.09) with moderate diversity and high relatedness between isolates in the community. In a remote border community, a clonal P. vivax cluster was identified, with distinct haplotypes in drug resistant genes and ama1, more similar to Brazilian isolates, likely representing an introduction of P. vivax from Brazil at that time. To test its applicability for Latin America, we evaluated the SNP Peru barcode in P. vivax genomes from the region and demonstrated the capacity to capture local population clustering at within-country level. CONCLUSIONS/SIGNIFICANCE: Together this data shows that P. vivax transmission is heterogeneous in different settings within the Peruvian Amazon. Genetic analysis is a key component for regional malaria control, offering valuable insights that should be incorporated into routine surveillance.

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 surveillance of Plasmodium falciparum drug-resistance markers in Vietnam using multiplex amplicon sequencing (2000-2016).

Emergence and spread of Plasmodium falciparum resistance to artemisinin-based combination therapies (ACT) is a major challenge for Greater Mekong Subregion countries in their goal to eliminate malaria by 2030. Tools to efficiently monitor drug resistance beyond resource-demanding therapeutic efficacy studies are necessary. A custom multiplex amplicon sequencing assay based on Illumina technology was designed to target the marker of partial resistance to artemisinin (K13), five candidate modulators of artemisinin resistance, the marker of resistance to chloroquine (crt), and four neutral microsatellite loci. The assay was used to genotype 635 P. falciparum-positive blood samples collected across seven provinces of Vietnam and one of Cambodia between 2000 and 2016. Markers of resistance to artemisinin partner-drugs piperaquine (copy number of plasmepsin-2) and mefloquine (copy number of multidrug-resistance 1) were determined by qPCR. Parasite population structure was further assessed using a 101-SNP barcode. Validated mutations of artemisinin partial resistance in K13 were found in 48.1% of samples, first detection was in 2000, and by 2015 prevalence overcame > 50% in Central Highlands and Binh Phuoc province. K13-C580Y variant became predominant country-wide, quickly replacing an outbreak of K13-I543T in Central Highlands. Mutations in candidate artemisinin resistance modulator genes paralleled the trends of K13 mutants, whereas resistance to piperaquine and mefloquine remained low (≈ 10%) by 2015-2016. Genomic tools applied to malaria surveillance generate comprehensive information on dynamics of drug resistance and population structure and reflect drug efficacy profiles from in vivo studies.

Molecular Surveillance of Malaria Using the PF AmpliSeq Custom Assay for Plasmodium falciparum Parasites from Dried Blood Spot DNA Isolates from Peru.

Malaria molecular surveillance has great potential to support national malaria control programs (NMCPs), informing policy for its control and elimination. Here, we present a new three-day workflow for targeted resequencing of markers in 13 resistance-associated genes, histidine rich protein 2 and 3 (hrp2&3), a country (Peru)-specific 28 SNP-barcode for population genetic analysis, and apical membrane antigen 1 (ama1), using Illumina short-read sequencing technology. The assay applies a multiplex PCR approach to amplify all genomic regions of interest in a rapid and easily standardizable procedure and allows simultaneous amplification of a high number of targets at once, therefore having great potential for implementation into routine surveillance practice by NMCPs. The assay can be performed on routinely collected filter paper blood spots and can be easily adapted to different regions to investigate either regional trends or in-country epidemiological changes.

Malaria Molecular Surveillance in the Peruvian Amazon with a Novel Highly Multiplexed Plasmodium falciparum AmpliSeq Assay.

Molecular surveillance for malaria has great potential to support national malaria control programs (NMCPs). To bridge the gap between research and implementation, several applications (use cases) have been identified to align research, technology development, and public health efforts. For implementation at NMCPs, there is an urgent need for feasible and cost-effective tools. We designed a new highly multiplexed deep sequencing assay (Pf AmpliSeq), which is compatible with benchtop sequencers, that allows high-accuracy sequencing with higher coverage and lower cost than whole-genome sequencing (WGS), targeting genomic regions of interest. The novelty of the assay is its high number of targets multiplexed into one easy workflow, combining population genetic markers with 13 nearly full-length resistance genes, which is applicable for many different use cases. We provide the first proof of principle for hrp2 and hrp3 deletion detection using amplicon sequencing. Initial sequence data processing can be performed automatically, and subsequent variant analysis requires minimal bioinformatic skills using any tabulated data analysis program. The assay was validated using a retrospective sample collection (n = 254) from the Peruvian Amazon between 2003 and 2018. By combining phenotypic markers and a within-country 28-single-nucleotide-polymorphism (SNP) barcode, we were able to distinguish different lineages with multiple resistance haplotypes (in dhfr, dhps, crt and mdr1) and hrp2 and hrp3 deletions, which have been increasing in recent years. We found no evidence to suggest the emergence of artemisinin (ART) resistance in Peru. These findings indicate a parasite population that is under drug pressure but is susceptible to current antimalarials and demonstrate the added value of a highly multiplexed molecular tool to inform malaria strategies and surveillance systems. IMPORTANCE While the power of next-generation sequencing technologies to inform and guide malaria control programs has become broadly recognized, the integration of genomic data for operational incorporation into malaria surveillance remains a challenge in most countries where malaria is endemic. The main obstacles include limited infrastructure, limited access to high-throughput sequencing facilities, and the need for local capacity to run an in-country analysis of genomes at a large-enough scale to be informative for surveillance. In addition, there is a lack of standardized laboratory protocols and automated analysis pipelines to generate reproducible and timely results useful for relevant stakeholders. With our standardized laboratory and bioinformatic workflow, malaria genetic surveillance data can be readily generated by surveillance researchers and malaria control programs in countries of endemicity, increasing ownership and ensuring timely results for informed decision- and policy-making.

Band 3-mediated Plasmodium vivax invasion is associated with transcriptional variation in PvTRAg genes.

The Plasmodium vivax reticulocyte invasion process is still poorly understood, with only a few receptor-ligand interactions identified to date. Individuals with the Southeast Asian ovalocytosis (SAO) phenotype have a deletion in the band 3 protein on the surface of erythrocytes, and are reported to have a lower incidence of clinical P. vivax malaria. Based on this observation, band 3 has been put forward as a receptor for P. vivax invasion, although direct proof is still lacking. In this study, we combined functional ex vivo invasion assays and transcriptome sequencing to uncover a band 3-mediated invasion pathway in P. vivax and potential band 3 ligands. Invasion by P. vivax field isolates was 67%-71% lower in SAO reticulocytes compared with non-SAO reticulocytes. Reticulocyte invasion was decreased by 40% and 27%-31% when blocking with an anti-band 3 polyclonal antibody and a PvTRAg38 peptide, respectively. To identify new band 3 receptor candidates, we mRNA-sequenced schizont-stage isolates used in the invasion assays, and observed high transcriptional variability in multigene and invasion-related families. Transcriptomes of isolates with low or high dependency on band 3 for invasion were compared by differential expression analysis, which produced a list of band 3 ligand candidates with high representation of PvTRAg genes. Our ex vivo invasion assays have demonstrated that band 3 is a P. vivax invasion receptor and confirm previous in vitro studies showing binding between PvTRAg38 and band 3, although the lower and variable inhibition levels observed suggest the involvement of other ligands. By coupling transcriptomes and invasion phenotypes from the same isolates, we identified a list of band 3 ligand candidates, of which the overrepresented PvTRAg genes are the most promising for future research.

Novel highly-multiplexed AmpliSeq targeted assay for Plasmodium vivax genetic surveillance use cases at multiple geographical scales.

Although the power of genetic surveillance tools has been acknowledged widely, there is an urgent need in malaria endemic countries for feasible and cost-effective tools to implement in national malaria control programs (NMCPs) that can generate evidence to guide malaria control and elimination strategies, especially in the case of Plasmodium vivax. Several genetic surveillance applications ('use cases') have been identified to align research, technology development, and public health efforts, requiring different types of molecular markers. Here we present a new highly-multiplexed deep sequencing assay (Pv AmpliSeq). The assay targets the 33-SNP vivaxGEN-geo panel for country-level classification, and a newly designed 42-SNP within-country barcode for analysis of parasite dynamics in Vietnam and 11 putative drug resistance genes in a highly multiplexed NGS protocol with easy workflow, applicable for many different genetic surveillance use cases. The Pv AmpliSeq assay was validated using: 1) isolates from travelers and migrants in Belgium, and 2) routine collections of the national malaria control program at sentinel sites in Vietnam. The assay targets 229 amplicons and achieved a high depth of coverage (mean 595.7 ± 481) and high accuracy (mean error-rate of 0.013 ± 0.007). P. vivax parasites could be characterized from dried blood spots with a minimum of 5 parasites/µL and 10% of minority-clones. The assay achieved good spatial specificity for between-country prediction of origin using the 33-SNP vivaxGEN-geo panel that targets rare alleles specific for certain countries and regions. A high resolution for within-country diversity in Vietnam was achieved using the designed 42-SNP within-country barcode that targets common alleles (median MAF 0.34, range 0.01-0.49. Many variants were detected in (putative) drug resistance genes, with different predominant haplotypes in the pvmdr1 and pvcrt genes in different provinces in Vietnam. The capacity of the assay for high resolution identity-by-descent (IBD) analysis was demonstrated and identified a high rate of shared ancestry within Gia Lai Province in the Central Highlands of Vietnam, as well as between the coastal province of Binh Thuan and Lam Dong. Our approach performed well in geographically differentiating isolates at multiple spatial scales, detecting variants in putative resistance genes, and can be easily adjusted to suit the needs in other settings in a country or region. We prioritize making this tool available to researchers and NMCPs in endemic countries to increase ownership and ensure data usage for decision-making and malaria policy.

Two fatal autochthonous cases of airport malaria, Belgium, 2020.

We report an outbreak investigation of two fatal cases of autochthonous Plasmodium falciparum malaria that occurred in Belgium in September 2020. Various hypotheses of the potential source of infection were investigated. The most likely route of transmission was through an infectious exotic Anopheles mosquito that was imported via the international airport of Brussels or the military airport Melsbroek and infected the cases who lived at 5 km from the airports. Based on genomic analysis of the parasites collected from the two cases, the most likely origin of the Plasmodium was Gabon or Cameroon. Further, the parasites collected from the two Belgian patients were identical by descent, which supports the assumption that the two infections originated from the bite of the same mosquito, during interrupted feeding. Although airport malaria remains a rare event, it has significant implications, particularly for the patient, as delayed or missed diagnosis of the cause of illness often results in complications and mortality. Therefore, to prevent such severe or fatal outcomes, we suggest a number of public health actions including increased awareness among health practitioners, especially those working in the vicinity of airports, and increased surveillance of exotic mosquito species at airports.

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

Monitoring the genetic structure of pathogen populations may be an economical and sensitive approach to quantify the impact of control on transmission dynamics, highlighting the need for a better understanding of changes in population genetic parameters as transmission declines. Here we describe the first population genetic analysis of two major human malaria parasites, Plasmodium falciparum (Pf) and Plasmodium vivax (Pv), following nationwide distribution of long-lasting insecticide-treated nets (LLINs) in Papua New Guinea (PNG). Parasite isolates from pre- (2005-2006) and post-LLIN (2010-2014) were genotyped using microsatellite markers. Despite parasite prevalence declining substantially (East Sepik Province: Pf = 54.9%-8.5%, Pv = 35.7%-5.6%, Madang Province: Pf = 38.0%-9.0%, Pv: 31.8%-19.7%), genetically diverse and intermixing parasite populations remained. Pf diversity declined modestly post-LLIN relative to pre-LLIN (East Sepik: Rs  = 7.1-6.4, HE  = 0.77-0.71; Madang: Rs  = 8.2-6.1, HE  = 0.79-0.71). Unexpectedly, population structure present in pre-LLIN populations was lost post-LLIN, suggesting that more frequent human movement between provinces may have contributed to higher gene flow. Pv prevalence initially declined but increased again in one province, yet diversity remained high throughout the study period (East Sepik: Rs  = 11.4-9.3, HE  = 0.83-0.80; Madang: Rs  = 12.2-14.5, HE  = 0.85-0.88). Although genetic differentiation values increased between provinces over time, no significant population structure was observed at any time point. For both species, a decline in multiple infections and increasing clonal transmission and significant multilocus linkage disequilibrium post-LLIN were positive indicators of impact on the parasite population using microsatellite markers. These parameters may be useful adjuncts to traditional epidemiological tools in the early stages of transmission reduction.

Differential impact of malaria control interventions on P. falciparum and P. vivax infections in young Papua New Guinean children.

INTRODUCTION: As malaria transmission declines, understanding the differential impact of intensified control on Plasmodium falciparum relative to Plasmodium vivax and identifying key drivers of ongoing transmission is essential to guide future interventions. METHODS: Three longitudinal child cohorts were conducted in Papua New Guinea before (2006/2007), during (2008) and after scale-up of control interventions (2013). In each cohort, children aged 1-5 years were actively monitored for infection and illness. Incidence of malaria episodes, molecular force of blood-stage infections (molFOB) and population-averaged prevalence of infections were compared across the cohorts to investigate the impact of intensified control in young children and the key risk factors for malaria infection and illness in 2013. RESULTS: Between 2006 and 2008, P. falciparum infection prevalence, molFOB, and clinical malaria episodes reduced by 47%, 59% and 69%, respectively, and a further 49%, 29% and 75% from 2008 to 2013 (prevalence 41.6% to 22.1% to 11.2%; molFOB: 3.4 to 1.4 to 1.0 clones/child/year; clinical episodes incidence rate (IR) 2.6 to 0.8 to IR 0.2 episodes/child/year). P. vivax clinical episodes declined at rates comparable to P. falciparum between 2006, 2008 and 2013 (IR 2.5 to 1.1 to 0.2), while P. vivax molFOB (2006, 9.8; 2008, 12.1) and prevalence (2006, 59.6%; 2008, 65.0%) remained high in 2008. However, in 2013, P. vivax molFOB (1.2) and prevalence (19.7%) had also substantially declined. In 2013, 89% of P. falciparum and 93% of P. vivax infections were asymptomatic, 62% and 47%, respectively, were sub-microscopic. Area of residence was the major determinant of malaria infection and illness. CONCLUSION: Intensified vector control and routine case management had a differential impact on rates of P. falciparum and P. vivax infections but not clinical malaria episodes in young children. This suggests comparable reductions in new mosquito-derived infections but a delayed impact on P. vivax relapsing infections due to a previously acquired reservoir of hypnozoites. This demonstrates the need to strengthen implementation of P. vivax radical cure to maximise impact of control in co-endemic areas. The high heterogeneity of malaria in 2013 highlights the importance of surveillance and targeted interventions to accelerate towards elimination.

Complement Receptor 1 availability on red blood cell surface modulates Plasmodium vivax invasion of human reticulocytes.

Plasmodium vivax parasites preferentially invade reticulocyte cells in a multistep process that is still poorly understood. In this study, we used ex vivo invasion assays and population genetic analyses to investigate the involvement of complement receptor 1 (CR1) in P. vivax invasion. First, we observed that P. vivax invasion of reticulocytes was consistently reduced when CR1 surface expression was reduced through enzymatic cleavage, in the presence of naturally low-CR1-expressing cells compared with high-CR1-expressing cells, and with the addition of soluble CR1, a known inhibitor of P. falciparum invasion. Immuno-precipitation experiments with P. vivax Reticulocyte Binding Proteins showed no evidence of complex formation. In addition, analysis of CR1 genetic data for worldwide human populations with different exposure to malaria parasites show significantly higher frequency of CR1 alleles associated with low receptor expression on the surface of RBCs and higher linkage disequilibrium in human populations exposed to P. vivax malaria compared with unexposed populations. These results are consistent with a positive selection of low-CR1-expressing alleles in vivax-endemic areas. Collectively, our findings demonstrate that CR1 availability on the surface of RBCs modulates P. vivax invasion. The identification of new molecular interactions is crucial to guiding the rational development of new therapeutic interventions against vivax malaria.

Characterization of Plasmodium falciparum and Plasmodium vivax recent exposure in an area of significantly decreased transmission intensity in Central Vietnam.

BACKGROUND: In Vietnam, malaria transmission has been reduced to very low levels over the past 20 years, and as a consequence, the country aims to eliminate malaria by 2030. This study aimed to characterize the dynamics and extent of the parasite reservoir in Central Vietnam, in order to further target elimination strategies and surveillance. METHODS: A 1-year prospective cohort study (n = 429) was performed in three rural communities in Quang Nam province. Six malaria screenings were conducted between November 2014 and November 2015, including systematic clinical examination and blood sampling for malaria parasite identification, as well as molecular and serological analysis of the study population. Malaria infections were detected by light microscopy (LM) and quantitative real time PCR (qPCR), while exposure to Plasmodium falciparum and Plasmodium vivax was measured in the first and last survey by ELISA for PfAMA1, PfGLURP R2, PvAMA1, and PvMSP1-19. Classification and regression trees were used to define seropositivity and recent exposure. RESULTS: Four malaria infections (2 P. falciparum, 2 P. vivax) were detected in the same village by qPCR and/or LM. No fever cases were attributable to malaria. At the same time, the commune health centre (serving a larger area) reported few cases of confirmed malaria cases. Nevertheless, serological data proved that 13.5% of the surveyed population was exposed to P. falciparum and/or P. vivax parasites during the study period, of which 32.6% were seronegative at the start of the study, indicating ongoing transmission in the area. Risk factor analysis for seroprevalence and exposure to P. falciparum and/or P. vivax identified structural or economic risk factors and activity/behaviour-related factors, as well as spatial heterogeneity at the village level. CONCLUSIONS: Previous studies in Central Vietnam demonstrated high occurrence of asymptomatic and sub-microscopic infections. However, in this study very few asymptomatic infections were detected despite serological evidence of continued transmission. Nonetheless, the factors associated with spatial heterogeneity in transmission could be evaluated using serological classification of recent exposure, which supports the usefulness of serological methods to monitor malaria transmission.

Micro-epidemiology of malaria in an elimination setting in Central Vietnam.

BACKGROUND: In Vietnam, malaria persists in remote forested regions where infections are spatially heterogeneous, mostly asymptomatic and with low parasite density. Previous studies in Vietnam have investigated broad behavioural concepts such as 'engaging in forest activities' as risk factors for malaria infection, which may not explain heterogeneity in malaria risk, especially in malaria elimination settings. METHODS: A mixed methods study combining ethnographic research and a cross-sectional survey was embedded in a 1-year malariometric cohort study in three ethnic minority villages in South Tra My district, Quang Nam Province in Central Vietnam. Qualitative data collection included in-depth interviews, informal conversations and participant observations over a 2-month period, and the findings were used to develop the questionnaire used in the cross-sectional survey. The latter collected data on evening activities, mobility patterns and household characteristics. The primary outcome, recent exposure to malaria, was defined using the classification and regression tree method to determine significant changes in antibody titres during the year preceding the survey. Risk factor analyses for recent exposure to malaria were conducted using logistic regression. RESULTS: 22 in-depth interviews and numerous participant observations were recorded during the ethnographic research (April to June 2015), and 160 adults (86% response rate) responded to the cross-sectional survey (November to December 2015). Recent exposure to Plasmodium falciparum malaria was estimated at 22.9 and at 17.1% for Plasmodium vivax. Ongoing malaria transmission appears to be maintained by activities that delay or disrupt sleeping in a permanent structure in which a bed net could be hung, including evening drinking gatherings, fishing, logging in the forest and outdoor TV watching. CONCLUSIONS: Vector control tools for outdoor evening activities in villages as well as at farms, forest and river locations should be incorporated into current malaria elimination efforts in Central Vietnam. Micro-epidemiology studies using mixed-methods designs can provide a comprehensive understanding of the malaria risk at fine spatial scales and better inform the implementation of targeted interventions for malaria elimination.

Sustained Malaria Control Over an 8-Year Period in Papua New Guinea: The Challenge of Low-Density Asymptomatic Plasmodium Infections.

Background: The scale-up of effective malaria control in the last decade has resulted in a substantial decline in the incidence of clinical malaria in many countries. The effects on the proportions of asymptomatic and submicroscopic infections and on transmission potential are yet poorly understood. Methods: In Papua New Guinea, vector control has been intensified since 2008, and improved diagnosis and treatment was introduced in 2012. Cross-sectional surveys were conducted in Madang Province in 2006 (with 1280 survey participants), 2010 (with 2117 participants), and 2014 (with 2516 participants). Infections were quantified by highly sensitive quantitative polymerase chain reaction (PCR) analysis, and gametocytes were quantified by reverse-transcription qPCR analysis. Results: Plasmodium falciparum prevalence determined by qPCR decreased from 42% in 2006 to 9% in 2014. The P. vivax prevalence decreased from 42% in 2006 to 13% in 2010 but then increased to 20% in 2014. Parasite densities decreased 5-fold from 2006 to 2010; 72% of P. falciparum and 87% of P. vivax infections were submicroscopic in 2014. Gametocyte density and positivity correlated closely with parasitemia, and population gametocyte prevalence decreased 3-fold for P. falciparum and 29% for P. vivax from 2010 to 2014. Conclusions: Sustained control has resulted in reduced malaria transmission potential, but an increasing proportion of gametocyte carriers are asymptomatic and submicroscopic and represent a challenge to malaria control.