Population Genetics of Plasmodium vivax in the Peruvian Amazon

Comparing newly developed SNP barcode panels with microsatellites to explore population genetics of malaria parasites in the Peruvian Amazon

Malaria molecular surveillance (MMS) can provide insights into transmission dynamics, guiding national control/elimination programs. Considering the genetic differences among parasites from different areas in the Peruvian Amazon, we previously designed SNP barcode panels for Plasmodium vivax (Pv) and P. falciparum (Pf), integrated into AmpliSeq assays, to provide population genetics estimates of malaria parasites. These AmpliSeq assays are ideal for MMS: multiplexing different traits of interest, applicable to many use cases, and high throughput for large numbers of samples. The present study compares the genetic resolution of the SNP barcode panels in the AmpliSeq assays with widely used microsatellite (MS) panels to investigate Amazonian malaria parasites. Malaria samples collected in remote areas of the Peruvian Amazon (51 Pv & 80 Pf samples) were characterized using the Ampliseq assays and MS. Population genetics estimates (complexity of infection, genetic diversity and differentiation, and population structure) were compared using the SNP barcodes (Pv: 40 SNPs & Pf: 28 SNPs) and MS panels (Pv: 16 MS & Pf: 7 MS). The genetic diversity of Pv (expected heterozygosity, He ) was similar across the subpopulations for both makers: He MS = 0.68 - 0.78 (p = 0.23) and He SNP = 0.36 - 0.38 (p = 0.80). Pairwise genetic differentiation (fixation index, F ST ) was also comparable: F ST-MS = 0.04 - 0.14 and F ST-SNP = 0.03 - 0.12 (p = 0.34 - 0.85). No geographic clustering was observed with any panel. In addition, Pf genetic diversity trends ( He MS = 0 - 0.48 p = 0.03 - 1; He SNP = 0 - 0.09, p = 0.03 - 1) and pairwise F ST comparisons (F ST-MS = 0.14 - 0.65, F ST-SNP = 0.19 - 0.61, p = 0.24 - 0.83) were concordant between the panels. Similar population structure clustering was observed with both SNP and MS, highlighting one Pf subpopulation in an indigenous community. The SNP barcodes in the Pv AmpliSeq v2 Peru and Pf AmpliSeq v1 Peru assays offer comparable results to MS panels when investigating population genetics in Pv and Pv populations. Therefore, the AmpliSeq assays can efficiently characterize malaria transmission dynamics and population structure and support malaria elimination efforts in Peru.

Lineage-informative microhaplotypes for recurrence classification and spatio-temporal surveillance of Plasmodium vivax malaria parasites

Challenges in classifying recurrent Plasmodium vivax infections constrain surveillance of antimalarial efficacy and transmission. Recurrent infections may arise from activation of dormant liver stages (relapse), blood-stage treatment failure (recrudescence) or reinfection. Molecular inference of familial relatedness (identity-by-descent or IBD) can help resolve the probable origin of recurrences. As whole genome sequencing of P. vivax remains challenging, targeted genotyping methods are needed for scalability. We describe a P. vivax marker discovery framework to identify and select panels of microhaplotypes (multi-allelic markers within small, amplifiable segments of the genome) that can accurately capture IBD. We evaluate panels of 50-250 microhaplotypes discovered in a global set of 615 P. vivax genomes. A candidate global 100-microhaplotype panel exhibits high marker diversity in the Asia-Pacific, Latin America and horn of Africa (median HE = 0.70-0.81) and identifies 89% of the polyclonal infections detected with genome-wide datasets. Data simulations reveal lower error in estimating pairwise IBD using microhaplotypes relative to traditional biallelic SNP barcodes. The candidate global panel also exhibits high accuracy in predicting geographic origin and captures local infection outbreak and bottlenecking events. Our framework is open-source enabling customised microhaplotype discovery and selection, with potential for porting to other species or data resources.

Genomic surveillance of malaria parasites in an indigenous community in the Peruvian Amazon

Hard-to-reach communities represent Peru's main challenge for malaria elimination, but information about transmission in these areas is scarce. Here, we assessed Plasmodium vivax (Pv) and P. falciparum (Pf) transmission dynamics, resistance markers, and Pf hrp2/3 deletions in Nueva Jerusalén (NJ), a remote, indigenous community in the Peruvian Amazon with high population mobility. We collected samples from November 2019 to May 2020 by active (ACD) and passive case detection (PCD) in NJ. Parasites were identified with microscopy and PCR. Then, we analyzed a representative set of positive-PCR samples (Pv = 68, Pf = 58) using highly-multiplexed deep sequencing assays (AmpliSeq) and compared NJ parasites with ones from other remote Peruvian areas using population genetics indexes. The ACD intervention did not reduce malaria cases in the short term, and persistent malaria transmission was observed (at least one Pv infection was detected in 96% of the study days). In Nueva Jerusalen, the Pv population had modest genetic diversity (He = 0.27). Pf population had lower diversity (He = 0.08) and presented temporal clustering, one of these clusters linked to an outbreak in February 2020. Moreover, Pv and Pf parasites from NJ exhibited variable levels of differentiation (Pv Fst = 0.07-0.52 and Pf Fst = 0.11-0.58) with parasites from other remote areas. No artemisin resistance mutations but chloroquine (57%) and sulfadoxine-pyrimethamine (35-67%) were detected in NJ's Pf parasites. Moreover, pfhrp2/3 gene deletions were common (32-50% of parasites with one or both genes deleted). The persistent Pv transmission and the detection of a Pf outbreak with parasites genetically distinct from the local ones highlight the need for tailored interventions focusing on mobility patterns and imported infections in remote areas to eliminate malaria in the Peruvian Amazon.

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.

High-Throughput Genotyping of Plasmodium vivax in the Peruvian Amazon via Molecular Inversion Probes

Plasmodium vivax transmission occurs throughout the tropics and is an emerging threat in areas of Plasmodium falciparum decline, causing relapse infections that complicate treatment and control. Targeted sequencing for P. falciparum has been widely deployed to detect population structure and the geographic spread of antimalarial and diagnostic resistance. However, there are fewer such tools for P. vivax . Leveraging global variation data, we designed four molecular inversion probe (MIP) genotyping panels targeting geographically differentiating SNPs, neutral SNPs, putative antimalarial resistance genes, and vaccine candidate genes. We deployed these MIP panels on 866 infections from the Peruvian Amazon and identified transmission networks with clonality (IBD>0.99), copy number variation in Pvdbp and multiple Pvrbps , fixation of putative antimalarial resistance, and balancing selection in 13 vaccine candidate genes. Our MIP panels are the broadest genotyping panel currently available and are poised for successful deployment in other regions of P. vivax transmission.

Study protocol: improving response to malaria in the Amazon through identification of inter-community networks and human mobility in border regions of Ecuador, Peru and Brazil

INTRODUCTION

Understanding human mobility's role in malaria transmission is critical to successful control and elimination. However, common approaches to measuring mobility are ill-equipped for remote regions such as the Amazon. This study develops a network survey to quantify the effect of community connectivity and mobility on malaria transmission.

METHODS

We measure community connectivity across the study area using a respondent driven sampling design among key informants who are at least 18 years of age. 45 initial communities will be selected: 10 in Brazil, 10 in Ecuador and 25 in Peru. Participants will be recruited in each initial node and administered a survey to obtain data on each community's mobility patterns. Survey responses will be ranked and the 2-3 most connected communities will then be selected and surveyed. This process will be repeated for a third round of data collection. Community network matrices will be linked with each country's malaria surveillance system to test the effects of mobility on disease risk.

ETHICS AND DISSEMINATION

This study protocol has been approved by the institutional review boards of Duke University (USA), Universidad San Francisco de Quito (Ecuador), Universidad Peruana Cayetano Heredia (Peru) and Universidade Federal Minas Gerais (Brazil). Results will be disseminated in communities by the end of the study.

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.

Genomic surveillance of malaria parasites in an indigenous community in the Peruvian Amazon

Hard-to-reach communities represent Peru's main challenge for malaria elimination, but information about transmission in these areas is scarce. Here, we assessed Plasmodium vivax (Pv) and P. falciparum (Pf) transmission dynamics, resistance markers, and Pf hrp2/3 deletions in Nueva Jerusalén (NJ), a remote, indigenous community in the Peruvian Amazon with high population mobility. We collected samples from November 2019 to May 2020 by active (ACD) and passive case detection (PCD) in NJ. Parasites were identified with microscopy and PCR. Then, we analyzed a representative set of positive-PCR samples (Pv = 68, Pf = 58) using highly-multiplexed deep sequencing assays (AmpliSeq) and compared NJ parasites with ones from other remote Peruvian areas using population genetics indexes. The ACD intervention did not reduce malaria cases in the short term, and persistent malaria transmission was observed (at least one Pv infection was detected in 96% of the study days). In Nueva Jerusalen, the Pv population had modest genetic diversity (He = 0.27). Pf population had lower diversity (He = 0.08) and presented temporal clustering, one of these clusters linked to an outbreak in February 2020. Moreover, Pv and Pf parasites from NJ exhibited variable levels of differentiation (Pv Fst = -0.52 & Pf Fst = 0.11-0.58) with parasites from other remote areas. No artemisin resistance mutations but chloroquine (57%) and sulfadoxine-pyrimethamine (35-67%) were detected in NJ's Pf parasites. Moreover, pfhrp2/3 gene deletions were common (32-50% of parasites with one or both genes deleted). The persistent Pv transmission and the detection of a Pf outbreak with parasites genetically distinct from the local ones highlight the need for tailored interventions focusing on mobility patterns and imported infections in remote areas to eliminate malaria in the Peruvian Amazon.

Network Profile: Improving Response to Malaria in the Amazon through Identification of Inter-Community Networks and Human Mobility in Border Regions of Ecuador, Peru, and Brazil

Objectives

Understanding human mobility's role on malaria transmission is critical to successful control and elimination. However, common approaches to measuring mobility are ill-equipped for remote regions such as the Amazon. This study develops a network survey to quantify the effect of community connectivity and mobility on malaria transmission.

Design

A community-level network survey.

Setting

We collect data on community connectivity along three river systems in the Amazon basin: the Pastaza river corridor spanning the Ecuador-Peru border; and the Amazon and Javari river corridors spanning the Brazil-Peru border.

Participants

We interviewed key informants in Brazil, Ecuador, and Peru, including from indigenous communities: Shuar, Achuar, Shiwiar, Kichwa, Ticuna, and Yagua. Key informants are at least 18 years of age and are considered community leaders.

Primary outcome

Weekly, community-level malaria incidence during the study period.

Methods

We measure community connectivity across the study area using a respondent driven sampling design. Forty-five communities were initially selected: 10 in Brazil, 10 in Ecuador, and 25 in Peru. Participants were recruited in each initial node and administered a survey to obtain data on each community's mobility patterns. Survey responses were ranked and the 2-3 most connected communities were then selected and surveyed. This process was repeated for a third round of data collection. Community network matrices will be linked with eadch country's malaria surveillance system to test the effects of mobility on disease risk.

Findings

To date, 586 key informants were surveyed from 126 communities along the Pastaza river corridor. Data collection along the Amazon and Javari river corridors is ongoing. Initial results indicate that network sampling is a superior method to delineate migration flows between communities.

Conclusions

Our study provides measures of mobility and connectivity in rural settings where traditional approaches are insufficient, and will allow us to understand mobility's effect on malaria transmission.

A new Plasmodium vivax reference genome for South American isolates

BACKGROUND

Plasmodium vivax is the second most important cause of human malaria worldwide, and accounts for the majority of malaria cases in South America. A high-quality reference genome exists for Papua Indonesia (PvP01) and Thailand (PvW1), but is lacking for South America. A reference genome specifically for South America would be beneficial though, as P. vivax is a genetically diverse parasite with geographical clustering.

RESULTS

This study presents a new high-quality assembly of a South American P. vivax isolate, referred to as PvPAM (P. vivax Peruvian AMazon). The genome was obtained from a low input patient sample from the Peruvian Amazon and sequenced using PacBio technology, resulting in a highly complete assembly with 6497 functional genes. Telomeric ends were present in 17 out of 28 chromosomal ends, and additional (sub)telomeric regions are present in 12 unassigned contigs. A comparison of multigene families between PvPAM and the PvP01 genome revealed remarkable variation in vir genes, and the presence of merozoite surface proteins (MSP) 3.6 and 3.7. Three dhfr and dhps drug resistance associated mutations are present in PvPAM, similar to those found in other Peruvian isolates. Mapping of publicly available South American whole genome sequencing (WGS) data to PvPAM resulted in significantly fewer variants and truncated reads compared to the use of PvP01 or PvW1 as reference genomes. To minimize the number of core genome variants in non-South American samples, PvW1 is most suited for Southeast Asian isolates, both PvPAM and PvW1 are suited for South Asian isolates, and PvPAM is recommended for African isolates. Interestingly, non-South American samples still contained the least subtelomeric variants when mapped to PvPAM, indicating high quality of the PvPAM subtelomeric regions.

CONCLUSIONS

Our findings show that the PvPAM reference genome more accurately represents South American P. vivax isolates in comparison to PvP01 and PvW1. In addition, PvPAM has a high level of completeness, and contains a similar number of annotated genes as PvP01 or PvW1. The PvPAM genome therefore will be a valuable resource to improve future genomic analyses on P. vivax isolates from the South American continent.

The impact of sustained malaria control in the Loreto region of Peru: a retrospective, observational, spatially-varying interrupted time series analysis of the PAMAFRO program

Background

Although malaria control investments worldwide have resulted in dramatic declines in transmission since 2000, progress has stalled. In the Amazon, malaria resurgence has followed withdrawal of Global Fund support of the Project for Malaria Control in Andean Border Areas (PAMAFRO). We estimate intervention-specific and spatially-explicit effects of the PAMAFRO program on malaria incidence across the Loreto region of Peru, and consider the influence of the environmental risk factors in the presence of interventions.

Methods

We conducted a retrospective, observational, spatial interrupted time series analysis of malaria incidence rates among people reporting to health posts across Loreto, Peru between the first epidemiological week of January 2001 and the last epidemiological week of December 2016. Model inference is at the smallest administrative unit (district), where the weekly number of diagnosed cases of Plasmodium vivax and Plasmodium falciparum were determined by microscopy. Census data provided population at risk. We include as covariates weekly estimates of minimum temperature and cumulative precipitation in each district, as well as spatially- and temporally-lagged malaria incidence rates. Environmental data were derived from a hydrometeorological model designed for the Amazon. We used Bayesian spatiotemporal modeling techniques to estimate the impact of the PAMAFRO program, variability in environmental effects, and the role of climate anomalies on transmission after PAMAFRO withdrawal.

Findings

During the PAMAFRO program, incidence of P. vivax declined from 42.8 to 10.1 cases/1000 people/year. Incidence for P. falciparum declined from 14.3 to 2.5 cases/1000 people/year over this same period. The effects of PAMAFRO-supported interventions varied both by geography and species of malaria. Interventions were only effective in districts where interventions were also deployed in surrounding districts. Further, interventions diminished the effects of other prevailing demographic and environmental risk factors. Withdrawal of the program led to a resurgence in transmission. Increasing minimum temperatures and variability and intensity of rainfall events from 2011 onward and accompanying population displacements contributed to this resurgence.

Interpretation

Malaria control programs must consider the climate and environmental scope of interventions to maximize effectiveness. They must also ensure financial sustainability to maintain local progress and commitment to malaria prevention and elimination efforts, as well as to offset the effects of environmental change that increase transmission risk.

Funding

National Aeronautics and Space Administration, National Institutes of Health, Bill and Melinda Gates Foundation.

Genetic diversity of Plasmodium vivax isolates from pregnant women in the Western Brazilian Amazon: a prospective cohort study

Background

Each year, 92 million pregnant women are at risk of contracting malaria during pregnancy, with the underestimation of the mortality and morbidity burden associated with Plasmodium vivax. During pregnancy, P. vivax infection is associated with low birth weight, maternal anaemia, premature delivery, and stillbirth. In the State of Acre (Brazil), high transmission leaves pregnant women at greater risk of contracting malaria and having a greater number of recurrences. The study of genetic diversity and the association of haplotypes with adverse pregnancy effects is of great importance for the control of the disease. Here we investigate the genetic diversity of P. vivax parasites infecting pregnant women across their pregnancies.

Methods

P. vivax DNA was extracted from 330 samples from 177 women followed during pregnancy, collected in the State of Acre, Brazil. All samples were negative for Plasmodium falciparum DNA. Sequence data for the Pvmsp1 gene was analysed alongside data from six microsatellite (MS) markers. Allelic frequencies, haplotype frequencies, expected heterozygosity (HE) were calculated. Whole genome sequencing (WGS) was conducted on four samples from pregnant women and phylogenetic analysis performed with other samples from South American regions.

Findings

Initially, the pregnant women were stratified into two groups-1 recurrence and 2 or more recurrences-in which no differences were observed in clinical gestational outcomes or in placental histological changes between the two groups. Then we evaluated the parasites genetically. An average of 18.5 distinct alleles were found at each of the MS loci, and the HE calculated for each marker indicates a high genetic diversity occurring within the population. There was a high percentage of polyclonal infections (61.7%, 108/175), and one haplotype (H1) occurred frequently (20%), with only 9 of the haplotypes appearing in more than one patient.

Interpretation

Most pregnant women had polyclonal infections that could be the result of relapses and/or re-infections. The high percentage of H1 parasites, along with the low frequency of many other haplotypes are suggestive of a clonal expansion. Phylogenetic analysis shows that P. vivax population within pregnant women clustered with other Brazilian samples in the region.

Funding

FAPESP and CNPq - Brazil.

Effect of out-of-village working activities on recent malaria exposure in the Peruvian Amazon using parametric g-formula

In the Amazon Region of Peru, occupational activities are important drivers of human mobility and may increase the individual risk of being infected while contributing to increasing malaria community-level transmission. Even though out-of-village working activities and other mobility patterns have been identified as determinants of malaria transmission, no studies have quantified the effect of out-of-village working activities on recent malaria exposure and proposed plausible intervention scenarios. Using two population-based cross-sectional studies in the Loreto Department in Peru, and the parametric g-formula method, we simulated various hypothetical scenarios intervening in out-of-village working activities to reflect their potential health benefits. This study estimated that the standardized mean outcome (malaria seroprevalence) in the unexposed population (no out-of-village workers) was 44.6% (95% CI: 41.7%-47.5%) and 66.7% (95% CI: 61.6%-71.8%) in the exposed population resulting in a risk difference of 22.1% (95% CI: 16.3%-27.9%). However, heterogeneous patterns in the effects of interest were observed between peri-urban and rural areas (Cochran's Q test = 15.5, p < 0.001). Heterogeneous patterns were also observed in scenarios of increased prevalence of out-of-village working activities and restriction scenarios by gender (male vs. female) and age (18 and under vs. 19 and older) that inform possible occupational interventions targetting population subgroups. The findings of this study support the hypothesis that targeting out-of-village workers will considerably benefit current malaria elimination strategies in the Amazon Region. Particularly, males and adult populations that carried out out-of-village working activities in rural areas contribute the most to the malaria seropositivity (recent exposure to the parasite) in the Peruvian Amazon.

Drug resistance and population structure of Plasmodium falciparum and Plasmodium vivax in the Peruvian Amazon

Malaria is a major health problem in Peru despite substantial progress achieved by the ongoing malaria elimination program. This study explored the population genetics of 63 Plasmodium falciparum and 170 P. vivax cases collected in the Peruvian Amazon Basin between 2015 and 2019. Microscopy and PCR were used for malaria detection and positive samples were genotyped at neutral and drug resistance-associated regions. The P. falciparum population exhibited a low nucleotide diversity (π = 0.02) whereas the P. vivax population presented a higher genetic diversity (π = 0.34). All P. falciparum samples (n = 63) carried chloroquine (CQ) resistant mutations on Pfcrt. Most P. falciparum samples (53 out of 54) carried sulfadoxine (SD) resistant mutations on Pfdhfr and Pfdhps. No evidence was found of artemisinin resistance mutations on kelch13. Population structure showed that a single cluster accounted for 93.4% of the P. falciparum samples whereas three clusters were found for P. vivax. Our study shows a low genetic diversity for both species with significant differences in genetic sub-structuring. The high prevalence of CQ-resistance mutations could be a result of indirect selection pressures driven by the P. vivax treatment scheme. These results could be useful for public health authorities to safeguard the progress that Peru has achieved towards malaria elimination.

Population genomics in neglected malaria parasites

Malaria elimination includes neglected human malaria parasites Plasmodium vivax, Plasmodium ovale spp., and Plasmodium malariae. Biological features such as association with low-density infection and the formation of hypnozoites responsible for relapse make their elimination challenging. Studies on these parasites rely primarily on clinical samples due to the lack of long-term culture techniques. With improved methods to enrich parasite DNA from clinical samples, whole-genome sequencing of the neglected malaria parasites has gained increasing popularity. Population genomics of more than 2200 P. vivax global isolates has improved our knowledge of parasite biology and host-parasite interactions, identified vaccine targets and potential drug resistance markers, and provided a new way to track parasite migration and introduction and monitor the evolutionary response of local populations to elimination efforts. Here, we review advances in population genomics for neglected malaria parasites, discuss how the rich genomic information is being used to understand parasite biology and epidemiology, and explore opportunities for the applications of malaria genomic data in malaria elimination practice.

Untangling population structure and genetic diversity of reticulocyte binding protein 2b (PvRBP2b) erythrocytic stage vaccine candidate in worldwide Plasmodium vivax isolates

BACKGROUNDS

Plasmodium vivax is the predominant Plasmodium species distributed extensively in the Americas and Asia-Pacific areas. Encoded protein by Plasmodium vivax Reticulocyte Binding Proteins (PvRBPs) family member are of critical prominence to parasite invasion and have been considered the significant targets in development of malaria vaccine for the blood stage. As high genetic polymorphism of parasites may impede the effectiveness of vaccine development, more research to unraveling genetic polymorphism of pvrbp2b from various geographical regions seems indispensable to map the exact pattern of field isolates.

METHODOLOGY/PRINCIPAL FINDINGS

The aim of this study was to determine the sequences of Iranian pvrbp2b (nt: 502-1896) gene and then, to ascertain polymorphism of pvrbp2b gene, recombination, the level of genetic distances, evaluation of natural selection, and the prediction of B-cell epitopes of Iranian and global P. vivax isolates. Pvrbp2b partial gene was amplified and sequenced from 60 Iranian P. vivax isolates. Iranian pvrbp2b sequences as well as 95 published sequences from five countries were used to evaluate the genetic diversity and neutral evolution signature in worldwide scale. A total of 38 SNPs were identified among 60 Iranian pvrbp2b sequences (32 non-synonymous and 6 synonymous mutations), and 32 amino acid substitutions were observed in 29 positions as compared to Sal-1 sequence. Worldwide sequence analysis showed that 44 amino acid changes had occurred in 37 positions of which seven polymorphic sites had trimorphic mutations while the rest was dimorphic. The overall nucleotide diversity for Iranian isolates was 0.00431 ± 0.00091 while the level of nucleotide diversity was ranged from 0.00337 ± 0.00076 (Peru) to 0.00452 ± 0.00092 (Thailand) in global scale.

CONCLUSIONS/SIGNIFICANCE

Of amino acid substitutions, 12 replacements were located in the B-cell epitopes in which nine polymorphic sites were positioned in N-terminal and three polymorphic sites in predicted B-cell epitopes of C-terminal, signifying both variable and conserved epitopes for vaccine designing. Using the achieved outcome of the current investigation interrogate questions to the selection of conserved regions of pvrbp2b and understanding polymorphism and immune system pressure to pave a way for developing a vaccine based on PvRBP2b candidate antigen.

Population genomic evidence of Plasmodium vivax Southeast Asian origin

Plasmodium vivax is the most common and widespread human malaria parasite. It was recently proposed that P. vivax originates from sub-Saharan Africa based on the circulation of its closest genetic relatives (P. vivax-like) among African great apes. However, the limited number of genetic markers and samples investigated questions the robustness of this hypothesis. Here, we extensively characterized the genomic variations of 447 human P. vivax strains and 19 ape P. vivax-like strains collected worldwide. Phylogenetic relationships between human and ape Plasmodium strains revealed that P. vivax is a sister clade of P. vivax-like, not included within the radiation of P. vivax-like By investigating various aspects of P. vivax genetic variation, we identified several notable geographical patterns in summary statistics in function of the increasing geographic distance from Southeast Asia, suggesting that P. vivax may have derived from a single area in Asia through serial founder effects.

Towards one standard treatment for uncomplicated Plasmodium falciparum and Plasmodium vivax malaria: Perspectives from and for the Peruvian Amazon

Malaria continues to wreak havoc in the Peruvian Amazon. Lengthy research efforts have brought important lessons on its particular epidemiology: the heterogeneous levels of transmission, the large reservoir of both asymptomatic and submicroscopic infections, the co-transmission of Plasmodium vivax and Plasmodium falciparum in the same areas, and the limitations of current diagnostics. Based on these features, the national elimination program could greatly benefit from simplified standard treatment, with the use of artemisinin-based combination therapy and even shorter schemes of primaquine maintaing the total dosing. It is acknowledged that there is some uncertainty regarding the true prevalence of glucose-6-phosphate dehydrogenase deficiency (G6PD) and genetic polymorphisms related to cytochrome P-450 isozyme 2D6 functioning. Once we have a better understanding, tafenoquine, whether or not in combination with a rapid G6PD enzyme test, may become a future pathway to eliminate the otherwise hidden reservoir of the P. vivax hypnozoite through one standard Plasmodium treatment.

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: R_s  = 7.1-6.4, H_E  = 0.77-0.71; Madang: R_s  = 8.2-6.1, H_E  = 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: R_s  = 11.4-9.3, H_E  = 0.83-0.80; Madang: R_s  = 12.2-14.5, H_E  = 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.

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.

Open-Source 3D Printable GPS Tracker to Characterize the Role of Human Population Movement on Malaria Epidemiology in River Networks: A Proof-of-Concept Study in the Peruvian Amazon

Human movement affects malaria epidemiology at multiple geographical levels; however, few studies measure the role of human movement in the Amazon Region due to the challenging conditions and cost of movement tracking technologies. We developed an open-source low-cost 3D printable GPS-tracker and used this technology in a cohort study to characterize the role of human population movement in malaria epidemiology in a rural riverine village in the Peruvian Amazon. In this pilot study of 20 participants (mean age = 40 years old), 45,980 GPS coordinates were recorded over 1 month. Characteristic movement patterns were observed relative to the infection status and occupation of the participants. Applying two analytical animal movement ecology methods, utilization distributions (UDs) and integrated step selection functions (iSSF), we showed contrasting environmental selection and space use patterns according to infection status. These data suggested an important role of human movement in the epidemiology of malaria in the Peruvian Amazon due to high connectivity between villages of the same riverine network, suggesting limitations of current community-based control strategies. We additionally demonstrate the utility of this low-cost technology with movement ecology analysis to characterize human movement in resource-poor environments.

Use of a Plasmodium vivax genetic barcode for genomic surveillance and parasite tracking in Sri Lanka

Background

Sri Lanka was certified as a malaria-free nation in 2016; however, imported malaria cases continue to be reported. Evidence-based information on the genetic structure/diversity of the parasite populations is useful to understand the population history, assess the trends in transmission patterns, as well as to predict threatening phenotypes that may be introduced and spread in parasite populations disrupting elimination programmes. This study used a previously developed Plasmodium vivax single nucleotide polymorphism (SNP) barcode to evaluate the population dynamics of P. vivax parasite isolates from Sri Lanka and to assess the ability of the SNP barcode for tracking the parasites to its origin.

Methods

A total of 51 P. vivax samples collected during 2005–2011, mainly from three provinces of the country, were genotyped for 40 previously identified P. vivax SNPs using a high-resolution melting (HRM), single-nucleotide barcode method. Minor allele frequencies, linkage disequilibrium, pair-wise F_ST values, and complexity of infection (COI) were evaluated to determine the genetic diversity. Structure analysis was carried out using STRUCTURE software (Version 2.3.4) and SNP barcode was used to identify the genetic diversity of the local parasite populations collected from different years. Principal component analysis (PCA) was used to determine the clustering according to global geographic regions.

Results

The proportion of multi-clone infections was significantly higher in isolates collected during an infection outbreak in year 2007. The minor allele frequencies of the SNPs changed dramatically from year to year. Significant linkage was observed in sample sub-sets from years 2005 and 2007. The majority of the isolates from 2007 consisted of at least two genetically distinct parasite strains. The overall percentage of multi-clone infections for the entire parasite sample was 39.21%. Analysis using STRUCTURE software (Version 2.3.4) revealed the high genetic diversity of the sample sub-set from year 2007. In-silico analysis of these data with those available from other global geographical regions using PCA showed distinct clustering of parasite isolates according to geography, demonstrating the usefulness of the barcode in determining an isolate to be indigenous.

Conclusions

Plasmodium vivax parasite isolates collected during a disease outbreak in year 2007 were more genetically diverse compared to those collected from other years. In-silico analysis using the 40 SNP barcode is a useful tool to track the origin of an isolate of uncertain origin, especially to differentiate indigenous from imported cases. However, an extended barcode with more SNPs may be needed to distinguish highly clonal populations within the country.

Temporal and Microspatial Heterogeneity in Transmission Dynamics of Coendemic Plasmodium vivax and Plasmodium falciparum in Two Rural Cohort Populations in the Peruvian Amazon

BACKGROUND

Malaria is highly heterogeneous: its changing malaria microepidemiology needs to be addressed to support malaria elimination efforts at the regional level.

METHODS

A 3-year, population-based cohort study in 2 settings in the Peruvian Amazon (Lupuna, Cahuide) followed participants by passive and active case detection from January 2013 to December 2015. Incidence and prevalence rates were estimated using microscopy and polymerase chain reaction (PCR).

RESULTS

Lupuna registered 1828 infections (1708 Plasmodium vivax, 120 Plasmodium falciparum; incidence was 80.7 infections/100 person-years (95% confidence interval [CI] , 77.1-84.5). Cahuide detected 1046 infections (1024 P vivax, 20 P falciparum, 2 mixed); incidence was 40.2 infections/100 person-years (95% CI, 37.9-42.7). Recurrent P vivax infections predominated onwards from 2013. According to PCR data, submicroscopic predominated over microscopic infections, especially in periods of low transmission. The integration of parasitological, entomological, and environmental observations evidenced an intense and seasonal transmission resilient to standard control measures in Lupuna and a persistent residual transmission after severe outbreaks were intensively handled in Cahuide.

CONCLUSIONS

In 2 exemplars of complex local malaria transmission, standard control strategies failed to eliminate submicroscopic and hypnozoite reservoirs, enabling persistent transmission.

Dynamics of Plasmodium vivax populations in border areas of the Greater Mekong sub-region during malaria elimination

Background

Countries within the Greater Mekong Sub-region (GMS) of Southeast Asia have committed to eliminating malaria by 2030. Although the malaria situation has greatly improved, malaria transmission remains at international border regions. In some areas, Plasmodium vivax has become the predominant parasite. To gain a better understanding of transmission dynamics, knowledge on the changes of P. vivax populations after the scale-up of control interventions will guide more effective targeted control efforts.

Methods

This study investigated genetic diversity and population structures in 206 P. vivax clinical samples collected at two time points in two international border areas: the China-Myanmar border (CMB) (n = 50 in 2004 and n = 52 in 2016) and Thailand-Myanmar border (TMB) (n = 50 in 2012 and n = 54 in 2015). Parasites were genotyped using 10 microsatellite markers.

Results

Despite intensified control efforts, genetic diversity remained high (H_E = 0.66–0.86) and was not significantly different among the four populations (P > 0.05). Specifically, H_E slightly decreased from 0.76 in 2004 to 0.66 in 2016 at the CMB and increased from 0.80 in 2012 to 0.86 in 2015 at the TMB. The proportions of polyclonal infections varied significantly among the four populations (P < 0.05), and showed substantial decreases from 48.0% in 2004 to 23.7 at the CMB and from 40.0% in 2012 to 30.7% in 2015 at the TMB, with corresponding decreases in the multiplicity of infection. Consistent with the continuous decline of malaria incidence in the GMS over time, there were also increases in multilocus linkage disequilibrium, suggesting more fragmented and increasingly inbred parasite populations. There were considerable genetic differentiation and sub-division among the four tested populations. Temporal genetic differentiation was observed at each site (F_ST = 0.081 at the CMB and F_ST = 0.133 at the TMB). Various degrees of clustering were evident between the older parasite samples collected in 2004 at the CMB and the 2016 CMB and 2012 TMB populations, suggesting some of these parasites had shared ancestry. In contrast, the 2015 TMB population was genetically distinctive, which may reflect a process of population replacement. Whereas the effective population size (N_e) at the CMB showed a decrease from 4979 in 2004 to 3052 in 2016 with the infinite allele model, the N_e at the TMB experienced an increase from 6289 to 10,259.

Conclusions

With enhanced control efforts on malaria, P. vivax at the TMB and CMB showed considerable spatial and temporal differentiation, but the presence of large P. vivax reservoirs still sustained genetic diversity and transmission. These findings provide new insights into P. vivax transmission dynamics and population structure in these border areas of the GMS. Coordinated and integrated control efforts on both sides of international borders are essential to reach the goal of regional malaria elimination.

Microsatellite analysis reveals connectivity among geographically distant transmission zones of Plasmodium vivax in the Peruvian Amazon: A critical barrier to regional malaria elimination

Despite efforts made over decades by the Peruvian government to eliminate malaria, Plasmodium vivax remains a challenge for public health decision-makers in the country. The uneven distribution of its incidence, plus its complex pattern of dispersion, has made ineffective control measures based on global information that lack the necessary detail to understand transmission fully. In this sense, population genetic tools can complement current surveillance. This study describes the genetic diversity and population structure from September 2012 to March 2015 in three geographically distant settlements, Cahuide (CAH), Lupuna (LUP) and Santa Emilia (STE), located in the Peruvian Amazon. A total 777 P. vivax mono-infections, out of 3264, were genotyped. Among study areas, LUP showed 19.7% of polyclonal infections, and its genetic diversity (H_exp) was 0.544. Temporal analysis showed a significant increment of polyclonal infections and H_exp, and the introduction and persistence of a new parasite population since March 2013. In STE, 40.1% of infections were polyclonal, with H_exp = 0.596. The presence of four genetic clusters without signals of clonal expansion and infections with lower parasite densities compared against the other two areas were also found. At least four parasite populations were present in CAH in 2012, where, after June 2014, malaria cases decreased from 213 to 61, concomitant with a decrease in polyclonal infections (from 0.286 to 0.18), and expectedly variable H_exp. Strong signals of gene flow were present in the study areas and wide geographic distribution of highly diverse parasite populations were found. This study suggests that movement of malaria parasites by human reservoirs connects geographically distant malaria transmission areas in the Peruvian Amazon. The maintenance of high levels of parasite genetic diversity through human mobility is a critical barrier to malaria elimination in this region.

Plasmodium vivax transmission is heterogeneous and discontinuous in the Peruvian Amazon. Such heterogeneity is the result of factors that include, but are not restricted to, the environment, public policies, and characteristics of the parasite, the vector, and human activities. All these factors make P. vivax transmission resilient to interventions. In order to achieve the goals of control and local elimination, P. vivax surveillance must inform how those factors sustain disease transmission in order to focalize and synchronize control strategies. In this study, we implemented molecular surveillance complemented with population genetic tools in the areas of Cahuide, Lupuna, and Santa Emilia located in the Peruvian Amazon. In particular, we characterize the transmission and the parasite genetic variation in these sites from September 2012 to March 2015. The changes in parasite diversity, the wide geographic dispersion of parasite subpopulation and the introduction of a new parasite clone or subpopulation in Lupuna documented in this study suggest that connectivity among the different endemic areas, likely due to human mobility, sustains disease transmission in the region hindering the success of control measures. This information must be considered in the design of current control strategies.

Limited differentiation among Plasmodium vivax populations from the northwest and to the south Pacific Coast of Colombia: A malaria corridor?

Background

Malaria remains endemic in several countries of South America with low to moderate transmission intensity. Regional human migration through underserved endemic areas may be responsible for significant parasite dispersion making the disease resilient to interventions. Thus, the genetic characterization of malarial parasites is an important tool to assess how endemic areas may connect via the movement of infected individuals. Here, four sites in geographically separated areas reporting 80% of the malaria morbidity in Colombia were studied. The sites are located on an imaginary transect line of 1,500 km from the northwest to the south Pacific Coast of Colombia with a minimal distance of 500 km between populations that display noticeable ethnic, economic, epidemiological, and ecological differences.

Methodology/Principal findings

A total of 624 Plasmodium vivax samples from the four populations were genotyped by using eight microsatellite loci. Although a strong geographic structure was expected between these populations, only moderate evidence of genetic differentiation was observed using a suite of population genetic analyses. High genetic diversity, shared alleles, and low linkage disequilibrium were also found in these P. vivax populations providing no evidence for a bottleneck or clonal expansions as expected from recent reductions in the transmission that could have been the result of scaling up interventions or environmental changes. These patterns are consistent with a disease that is not only endemic in each site but also imply that there is gene flow among these populations across 1,500 km.

Conclusion /Significance

The observed patterns in P. vivax are consistent with a “corridor” where connected endemic areas can sustain a high level of genetic diversity locally and can restore parasite-subdivided populations via migration of infected individuals even after local interventions achieved a substantial reduction of clinical cases. The consequences of these findings in terms of control and elimination are discussed.

The regional movements of infected individuals that connect suitable transmission areas make malaria resilient to control efforts. Those movements are expected to leave genetic signatures in the parasite populations that can be detected using analytical tools. In this study, the genetic makeups of Plasmodium vivax populations were characterized to assess whether the most endemic areas in Colombia were connected. Samples were collected from passive surveillance studies in four locations across an imaginary transect line of 1,500 km from the northwest to the south Pacific Coast of Colombia (South America). Considering the distance, and contrary to expectations, we found weak levels of genetic differentiation between these parasite populations with no evidence indicating that their genetic diversity has been eroded as expected whenever the prevalence of the disease is successfully reduced, e.g., through control programs or environmental changes. Although the sampling lacks the geographic and temporal detail to describe how the dispersion of parasite lineages occurred, the observed patterns are consistent with a series of infected populations that are connected in space by human movements allowing the parasite to diffuse across this 1,500 km transect. This malaria corridor needs to be characterized to achieve elimination.

Use of open mobile mapping tool to assess human mobility traceability in rural offline populations with contrasting malaria dynamics

Infectious disease dynamics are affected by human mobility more powerfully than previously thought, and thus reliable traceability data are essential. In rural riverine settings, lack of infrastructure and dense tree coverage deter the implementation of cutting-edge technology to collect human mobility data. To overcome this challenge, this study proposed the use of a novel open mobile mapping tool, GeoODK. This study consists of a purposive sampling of 33 participants in six villages with contrasting patterns of malaria transmission that demonstrates a feasible approach to map human mobility. The self-reported traceability data allowed the construction of the first human mobility framework in rural riverine villages in the Peruvian Amazon. The mobility spectrum in these areas resulted in travel profiles ranging from 2 hours to 19 days; and distances between 10 to 167 km. Most Importantly, occupational-related mobility profiles with the highest displacements (in terms of time and distance) were observed in commercial, logging, and hunting activities. These data are consistent with malaria transmission studies in the area that show villages in watersheds with higher human movement are concurrently those with greater malaria risk. The approach we describe represents a potential tool to gather critical information that can facilitate malaria control activities.

Intra-host dynamics of co-infecting parasite genotypes in asymptomatic malaria patients

Malaria-infected individuals often harbor mixtures of genetically distinct parasite genotypes. We studied intra-host dynamics of parasite genotypes co-infecting asymptomatic adults in an area of intense malaria transmission in Chikhwawa, Malawi. Serial blood samples (5 ml) were collected over seven consecutive days from 25 adults with asymptomatic Plasmodium falciparum malaria and analyzed to determine whether a single peripheral blood sample accurately captures within-host parasite diversity. Blood samples from three of the participants were also analyzed by limiting dilution cloning and SNP genotyping of the parasite clones isolated to examine both the number and relatedness of co-infecting parasite haplotypes. We observed rapid turnover of co-infecting parasite genotypes in 88% of the individuals sampled (n = 22) such that the genetic composition of parasites infecting these individuals changed dramatically over the course of seven days of follow up. Nineteen of the 25 individuals sampled (76%) carried multiple parasite genotypes at baseline. Analysis of serial blood samples from three of the individuals revealed that they harbored 6, 12 and 17 distinct parasite haplotypes respectively. Approximately 70% of parasite haplotypes recovered from the three extensively sampled individuals were unrelated (proportion of shared alleles <83.3%) and were deemed to have primarily arisen from superinfection (inoculation of unrelated parasite haplotypes through multiple mosquito bites). The rest were related at the half-sib level or greater and were deemed to have been inoculated into individual human hosts via parasite co-transmission from single mosquito bites. These findings add further to the growing weight of evidence indicating that a single blood sample poorly captures within-host parasite diversity and underscore the importance of repeated blood sampling to accurately capture within-host parasite ecology. Our data also demonstrate a more pronounced role for parasite co-transmission in generating within-host parasite diversity in high transmission settings than previously assumed. Taken together, these findings have important implications for understanding the evolution of drug resistance, malaria transmission, parasite virulence, allocation of gametocyte sex ratios and acquisition of malaria immunity.

Exploration of Plasmodium vivax transmission dynamics and recurrent infections in the Peruvian Amazon using whole genome sequencing

BACKGROUND

Plasmodium vivax poses a significant challenge to malaria elimination due to its ability to cause relapsed infections from reactivation of dormant liver parasites called hypnozoites. We analyzed 69 P. vivax whole genome sequences obtained from subjects residing in three different villages along the Peruvian Amazon. This included 23 paired P. vivax samples from subjects who experienced recurrent P. vivax parasitemia following observed treatment with chloroquine and primaquine.

METHODS

Genomic DNA was extracted from whole blood samples collected from subjects. P. vivax DNA was enriched using selective whole genome amplification and whole genome sequencing. We used single nucleotide polymorphisms (SNPs) from the core P. vivax genome to determine characteristics of the parasite population using discriminant analysis of principal components, maximum likelihood estimation of individual ancestries, and phylogenetic analysis. We estimated the relatedness of the paired samples by calculating the number of segregating sites and using a hidden Markov model approach to estimate identity by descent.

RESULTS

We present a comprehensive dataset of population genetics of Plasmodium vivax in the Peruvian Amazonian. We define the parasite population structure in this region and demonstrate a novel method for distinguishing homologous relapses from reinfections or heterologous relapses with improved accuracy. The parasite population in this area was quite diverse with an estimated five subpopulations and evidence of a highly heterogeneous ancestry of some of the isolates, similar to previous analyses of P. vivax in this region. Pairwise comparison of recurrent infections determined that there were 12 homologous relapses and 3 likely heterologous relapses with highly related parasites. To the best of our knowledge, this is the first large-scale study to evaluate recurrent P. vivax infections using whole genome sequencing.

CONCLUSIONS

Whole genome sequencing is a high-resolution tool that can identify P. vivax homologous relapses with increased sensitivity, while also providing data about drug resistance and parasite population genetics. This information is important for evaluating the efficacy of known and novel antirelapse medications in endemic areas and thus advancing the campaign to eliminate malaria.

Genetic Variability of Plasmodium vivax in the North Coast of Peru and the Ecuadorian Amazon Basin

In the Peruvian North Coast (PNC), the number of Plasmodium vivax malaria cases increased steadily from 2007 to 2010 despite a significant decline in the overall number of cases in Peru during the same period. To better understand the transmission dynamics of P. vivax populations in the PNC and the neighboring Ecuadorian Amazon Basin (EAB), we studied the genetic variability and population structure of P. vivax in these areas. One hundred and twenty P. vivax isolates (58 from Piura and 37 from Tumbes in the PNC collected from 2008 to 2010 and 25 from the EAB collected in Pastaza from 2001 to 2004) were assessed by five polymorphic microsatellite markers. Genetic variability was determined by expected heterozygosity (He) and population structure by Bayesian inference cluster analysis. We found very low genetic diversity in the PNC (He = 0-0.32) but high genetic diversity in the EAB (He = 0.43-0.70). Population structure analysis revealed three distinct populations in the three locations. Six of 37 (16%) isolates from Tumbes had an identical haplotype to that found in Piura, suggesting unidirectional flow from Piura to Tumbes. In addition, one haplotype from Tumbes showed similarity to a haplotype found in Pastaza, suggesting that this could be an imported case from EAB. These findings strongly suggest a minimal population flow and different levels of genetic variability between these two areas divided by the Andes Mountains. This work presents molecular markers that could be used to increase our understanding of regional malaria transmission dynamics, which has implications for the development of strategies for P. vivax control.

Molecular approaches to determine the multiplicity of Plasmodium infections

Multiplicity of infection (MOI), also termed complexity of infection (COI), is defined as the number of genetically distinct parasite strains co-infecting a single host, which is an important indicator of malaria epidemiology. PCR-based genotyping often underestimates MOI. Next generation sequencing technologies provide much more accurate and genome-wide characterization of polyclonal infections. However, complete haplotype characterization of multiclonal infections remains a challenge due to PCR artifacts and sequencing errors, and requires efficient computational tools. In this review, the advantages and limitations of current molecular approaches to determine multiplicity of malaria parasite infection are discussed.

Evolution of life in urban environments

Our planet is an increasingly urbanized landscape, with over half of the human population residing in cities. Despite advances in urban ecology, we do not adequately understand how urbanization affects the evolution of organisms, nor how this evolution may affect ecosystems and human health. Here, we review evidence for the effects of urbanization on the evolution of microbes, plants, and animals that inhabit cities. Urbanization affects adaptive and nonadaptive evolutionary processes that shape the genetic diversity within and between populations. Rapid adaptation has facilitated the success of some native species in urban areas, but it has also allowed human pests and disease to spread more rapidly. The nascent field of urban evolution brings together efforts to understand evolution in response to environmental change while developing new hypotheses concerning adaptation to urban infrastructure and human socioeconomic activity. The next generation of research on urban evolution will provide critical insight into the importance of evolution for sustainable interactions between humans and our city environments.

Increasingly inbred and fragmented populations of Plasmodium vivax associated with the eastward decline in malaria transmission across the Southwest Pacific

The human malaria parasite Plasmodium vivax is more resistant to malaria control strategies than Plasmodium falciparum, and maintains high genetic diversity even when transmission is low. To investigate whether declining P. vivax transmission leads to increasing population structure that would facilitate elimination, we genotyped samples from across the Southwest Pacific region, which experiences an eastward decline in malaria transmission, as well as samples from two time points at one site (Tetere, Solomon Islands) during intensified malaria control. Analysis of 887 P. vivax microsatellite haplotypes from hyperendemic Papua New Guinea (PNG, n = 443), meso-hyperendemic Solomon Islands (n = 420), and hypoendemic Vanuatu (n = 24) revealed increasing population structure and multilocus linkage disequilibrium yet a modest decline in diversity as transmission decreases over space and time. In Solomon Islands, which has had sustained control efforts for 20 years, and Vanuatu, which has experienced sustained low transmission for many years, significant population structure was observed at different spatial scales. We conclude that control efforts will eventually impact P. vivax population structure and with sustained pressure, populations may eventually fragment into a limited number of clustered foci that could be targeted for elimination.

Plasmodium vivax is a major human malaria parasite, common in endemic areas outside sub-Saharan Africa, and more difficult to control than other malaria parasite species. The distinct lifecycle biology of P. vivax is thought to contribute to its more stable and efficient transmission allowing the maintenance of high diversity and potentially, gene flow. Independent studies are therefore needed to understand how P. vivax populations respond to changing transmission levels, in order to inform malaria control and elimination efforts. Here we have determined parasite population genetic structure in three countries of the Southwest Pacific, an island chain with a natural west to east decline in transmission intensity (Papua New Guinea > Solomon Islands > Vanuatu). With declining transmission, P. vivax populations experience only a modest decline in diversity but a significant increase in multilocus linkage disequilibrium and population structure, indicating that parasite populations become more inbred and begin to fragment into clustered foci. Analysis of two time points in one study area (Tetere, Solomon Islands) also show similar changes in association with intensifying malaria control. The results indicate that with long term sustained malaria control P. vivax populations will eventually fracture into population clusters that could be targeted for elimination.

High prevalence of very-low Plasmodium falciparum and Plasmodium vivax parasitaemia carriers in the Peruvian Amazon: insights into local and occupational mobility-related transmission

BACKGROUND

The incidence of malaria due both to Plasmodium falciparum and Plasmodium vivax in the Peruvian Amazon has risen in the past 5 years. This study tested the hypothesis that the maintenance and emergence of malaria in hypoendemic regions such as Amazonia is determined by submicroscopic and asymptomatic Plasmodium parasitaemia carriers. The present study aimed to precisely quantify the rate of very-low parasitaemia carriers in two sites of the Peruvian Amazon in relation to transmission patterns of P. vivax and P. falciparum in this area.

METHODS

This study was carried out within the Amazonian-ICEMR longitudinal cohort. Blood samples were collected for light microscopy diagnosis and packed red blood cell (PRBC) samples were analysed by qPCR. Plasma samples were tested for total IgG reactivity against recombinant PvMSP-10 and PfMSP-10 antigens by ELISA. Occupation and age 10 years and greater were considered surrogates of occupation-related mobility. Risk factors for P. falciparum and P. vivax infections detected by PRBC-qPCR were assessed by multilevel logistic regression models.

RESULTS

Among 450 subjects, the prevalence of P. vivax by PRBC-PCR (25.1%) was sixfold higher than that determined by microscopy (3.6%). The prevalence of P. falciparum infection was 4.9% by PRBC-PCR and 0.2% by microscopy. More than 40% of infections had parasitaemia under 5 parasites/μL. Multivariate analysis for infections detected by PRBC-PCR showed that participants with recent settlement in the study area (AOR 2.1; 95% CI 1.03:4.2), age ≥ 30 years (AOR 3.3; 95% CI 1.6:6.9) and seropositivity to P. vivax (AOR 1.8; 95% CI 1.0:3.2) had significantly higher likelihood of P. vivax infection, while the odds of P. falciparum infection was higher for participants between 10 and 29 years (AOR 10.7; 95% CI 1.3:91.1) and with a previous P. falciparum infection (AOR 10.4; 95% CI 1.5:71.1).

CONCLUSIONS

This study confirms the contrasting transmission patterns of P. vivax and P. falciparum in the Peruvian Amazon, with stable local transmission for P. vivax and the source of P. falciparum to the study villages dominated by very low parasitaemia carriers, age 10 years and older, who had travelled away from home for work and brought P. falciparum infection with them.

The Biology of Plasmodium vivax

Plasmodium vivax is the second most prevalent cause of malaria worldwide and the leading cause of malaria outside of Africa. Although infections are seldom fatal clinical disease can be debilitating and imposes significant health and economic impacts on affected populations. Estimates of transmission and prevalence intensity can be problematic because many episodes of vivax originate from hypnozoite stages in the liver that have remained dormant from previous infections by an unknown mechanism. Lack of treatment options to clear hypnozoites and the ability to infect mosquitoes before disease symptoms present represent major challenges for control and eradication of vivax malaria. Compounding these challenges is the unique biology of P. vivax and limited progress in development of experimental research tools, thereby hindering development of new drugs and vaccines. Renewed emphasis on vivax malaria research is beginning to make progress in overcoming some of these challenges.

Micro-epidemiology and spatial heterogeneity of P. vivax parasitaemia in riverine communities of the Peruvian Amazon: A multilevel analysis

Malaria has steadily increased in the Peruvian Amazon over the last five years. This study aimed to determine the parasite prevalence and micro-geographical heterogeneity of Plasmodium vivax parasitaemia in communities of the Peruvian Amazon. Four cross-sectional active case detection surveys were conducted between May and July 2015 in four riverine communities in Mazan district. Analysis of 2785 samples of 820 individuals nested within 154 households for Plasmodium parasitaemia was carried out using light microscopy and qPCR. The spatio-temporal distribution of Plasmodium parasitaemia, dominated by P. vivax, was shown to cluster at both household and community levels. Of enrolled individuals, 47% had at least one P. vivax parasitaemia and 10% P. falciparum, by qPCR, both of which were predominantly sub-microscopic and asymptomatic. Spatial analysis detected significant clustering in three communities. Our findings showed that communities at small-to-moderate spatial scales differed in P. vivax parasite prevalence, and multilevel Poisson regression models showed that such differences were influenced by factors such as age, education, and location of households within high-risk clusters, as well as factors linked to a local micro-geographic context, such as travel and occupation. Complex transmission patterns were found to be related to human mobility among communities in the same micro-basin.

Predominance of asymptomatic and sub-microscopic infections characterizes the Plasmodium gametocyte reservoir in the Peruvian Amazon

Malaria transmission requires that Anopheles mosquitoes ingest Plasmodium gametocyte stages circulating in the human bloodstream. In the context of malaria elimination, understanding the epidemiology of gametocytes relative to all Plasmodium infections and the contribution of asymptomatic and sub-microscopic parasite carriers to the gametocyte reservoir is necessary, especially in low endemic settings with predominance of P.vivax. A 13-month longitudinal study was conducted in two communities (n = 1935 individuals) of Loreto Department, Peru, with five active screenings for Plasmodium infections and gametocyte stages by quantitative real-time PCR (qPCR) and reverse transcription (RT)-qPCR, respectively. Parasite prevalence by qPCR was 7.2% for P.vivax (n = 520/7235; range by survey 6.0%-8.1%) and 3.2% for P.falciparum (n = 235/7235; range by survey 0.4%-7.7%). Sub-microscopic infections accounted for 73.5% of P.vivax (range by survey 60%-89%) and almost the totality of P.falciparum cases. Gametocytes were found in 28.4% P.vivax infections (range by survey 18.7%-34.1%), with a peak of 61.5% in one community at the start of the transmission season. About 59.8% of all P.vivax gametocyte carriers were asymptomatic and 31.9% were sub-microscopic. Age patterns for gametocyte prevalence paralleled asexual stage infections and peaked among >15–25 year old individuals. Asexual parasite density was found to be the strongest predictor for P.vivax gametocyte presence in longitudinal multivariate analysis (odds ratio 2.33 [95% confidence interval 1.96, 2.78]; P<0.001). Despite significant differences in seasonality patterns and P.vivax prevalence found at the local scale, sub-microscopic and asymptomatic infections predominate and contribute significantly to the gametocyte reservoir in different communities of the Peruvian Amazon. Control and elimination campaigns need sensitive tools to detect all infections that escape routine malaria surveillance, which may contribute to maintain transmission in the region.

Malaria elimination, i.e. the complete interruption of parasite transmission in a region, is in the agenda of health authorities in countries that achieved substantial reduction of the disease burden in the past decade. However, our understanding of transmission epidemiology for low transmission areas where Plasmodium vivax is endemic, like the Amazon basin, is still limited. In this study, we describe the prevalence and risk factors for carrying the parasite stages that are transmitted to the mosquito vectors, named gametocytes, in 1935 individuals from two communities of the Peruvian Amazon that were regularly screened during 1 year. We report that malaria infections with no clinical symptoms and those with parasite levels below microscopy detection threshold, account for two thirds of all P.vivax infections with gametocytes, and that the highest infection rate is found among young adults. In addition, almost the totality of P.falciparum infections detected was sub-microscopic. Because all these infections escape current malaria surveillance systems -based on passive case detection and/or microscopy diagnosis-, new approaches are necessary to target all infections in order to eliminate the malaria transmission reservoir in Peru.

Different origin and dispersal of sulfadoxine-resistant Plasmodium falciparum haplotypes between Eastern Africa and Democratic Republic of Congo

Sulfadoxine/pyrimethamine (SP) is still used for malaria control in sub-Saharan Africa; however, widespread resistance is a major concern. This study aimed to determine the dispersal and origin of sulfadoxine resistance lineages in the Democratic Republic of the Congo compared with East African Plasmodium falciparum dihydropteroate synthetase (Pfdhps) haplotypes. The analysis involved 264 isolates collected from patients with uncomplicated malaria from Tanzania, Uganda and DR Congo. Isolates were genotyped for Pfdhps mutations at codons 436, 437, 540, 581 and 613. Three microsatellite loci (0.8, 4.3 and 7.7 kb) flanking the Pfdhps gene were assayed. Evolutionary analysis revealed a shared origin of Pfdhps haplotypes in East Africa, with a distinct population clustering in DR Congo. Furthermore, in Tanzania there was an independent distinct origin of Pfdhps SGEGA resistant haplotype. In Uganda and Tanzania, gene flow patterns contribute to the dispersal and shared origin of parasites carrying double- and triple-mutant Pfdhps haplotypes associated with poor outcomes of intermittent preventive treatment during pregnancy using SP (IPTp-SP). However, the origins of the Pfdhps haplotypes in DR Congo and Eastern Africa sites are different. The genetic structure demonstrated a divergent and distinct population cluster predominated by single-mutant Pfdhps haplotypes at the DR Congo site. This reflects the limited dispersal of double- and triple-mutant Pfdhps haplotypes in DR Congo. This study highlights the current genetic structure and dispersal of high-grade Pfdhps resistant haplotypes, which is important to guide implementation of SP in malaria chemoprevention strategies in the region.

Dissecting malaria biology and epidemiology using population genetics and genomics

Molecular approaches have an increasingly recognized utility in surveillance of malaria parasite populations, not only in defining prevalence and incidence with higher sensitivity than traditional methods, but also in monitoring local and regional parasite transmission patterns. In this review, we provide an overview of population genetic and genomic studies of human-infecting Plasmodium species, highlighting recent advances in the field. In accordance with the renewed impetus for malaria eradication, many studies are now using genetic and genomic epidemiology to support local evidence-based intervention strategies. Microsatellite genotyping remains a popular approach for both Plasmodium falciparum and Plasmodium vivax. However, with the increasing availability of whole genome sequencing data enabling effective single nucleotide polymorphism-based panels tailored to a given study question and setting, this approach is gaining popularity. The availability of new reference genomes for Plasmodium malariae and Plasmodium ovale should see a surge in similar molecular studies on these currently neglected species. Genomic studies are revealing new insights into important adaptive mechanisms of the parasite including antimalarial drug resistance. The advent of new methodologies such as selective whole genome amplification for dealing with extensive human DNA in low density field isolates should see genome-wide approaches becoming routine for parasite surveillance once the economic costs outweigh the current cost benefits of targeted approaches.

Epidemiology of Plasmodium vivax Malaria in Peru

Malaria in Peru, dominated by Plasmodium vivax, remains a public health problem. The 1990s saw newly epidemic malaria emerge, primarily in the Loreto Department in the Amazon region, including areas near to Iquitos, the capital city, but sporadic malaria transmission also occurred in the 1990s–2000s in both north-coastal Peru and the gold mining regions of southeastern Peru. Although a Global Fund-supported intervention (PAMAFRO, 2005–2010) was temporally associated with a decrease of malaria transmission, from 2012 to the present, both P. vivax and Plasmodium falciparum malaria cases have rapidly increased. The Peruvian Ministry of Health continues to provide artemesinin-based combination therapy for microscopy-confirmed cases of P. falciparum and chloroquine–primaquine for P. vivax. Malaria transmission continues in remote areas nonetheless, where the mobility of humans and parasites facilitates continued reintroduction outside of ongoing surveillance activities, which is critical to address for future malaria control and elimination efforts. Ongoing P. vivax research gaps in Peru include the following: identification of asymptomatic parasitemics, quantification of the contribution of patent and subpatent parasitemics to mosquito transmission, diagnosis of nonparasitemic hypnozoite carriers, and implementation of surveillance for potential emergence of chloroquine- and 8-aminoquinoline-resistant P. vivax. Clinical trials of tafenoquine in Peru have been promising, and glucose-6-phosphate dehydrogenase deficiency in the region has not been observed to be a limitation to its use. Larger-scale challenges for P. vivax (and malaria in general) in Peru include logistical difficulties in accessing remote riverine populations, consequences of government policy and poverty trends, and obtaining international funding for malaria control and elimination.