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

PyPop: a mature open-source software pipeline for population genomics

Python for Population Genomics (PyPop) is a software package that processes genotype and allele data and performs large-scale population genetic analyses on highly polymorphic multi-locus genotype data. In particular, PyPop tests data conformity to Hardy-Weinberg equilibrium expectations, performs Ewens-Watterson tests for selection, estimates haplotype frequencies, measures linkage disequilibrium, and tests significance. Standardized means of performing these tests is key for contemporary studies of evolutionary biology and population genetics, and these tests are central to genetic studies of disease association as well. Here, we present PyPop 1.0.0, a new major release of the package, which implements new features using the more robust infrastructure of GitHub, and is distributed via the industry-standard Python Package Index. New features include implementation of the asymmetric linkage disequilibrium measures and, of particular interest to the immunogenetics research communities, support for modern nomenclature, including colon-delimited allele names, and improvements to meta-analysis features for aggregating outputs for multiple populations. Code available at: https://zenodo.org/records/10080668 and https://github.com/alexlancaster/pypop.

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.

Malaria Resilience in South America: Epidemiology, Vector Biology, and Immunology Insights from the Amazonian International Center of Excellence in Malaria Research Network in Peru and Brazil

The 1990s saw the rapid reemergence of malaria in Amazonia, where it remains an important public health priority in South America. The Amazonian International Center of Excellence in Malaria Research (ICEMR) was designed to take a multidisciplinary approach toward identifying novel malaria control and elimination strategies. Based on geographically and epidemiologically distinct sites in the Northeastern Peruvian and Western Brazilian Amazon regions, synergistic projects integrate malaria epidemiology, vector biology, and immunology. The Amazonian ICEMR's overarching goal is to understand how human behavior and other sociodemographic features of human reservoirs of transmission-predominantly asymptomatically parasitemic people-interact with the major Amazonian malaria vector, Nyssorhynchus (formerly Anopheles) darlingi, and with human immune responses to maintain malaria resilience and continued endemicity in a hypoendemic setting. Here, we will review Amazonian ICEMR's achievements on the synergies among malaria epidemiology, Plasmodium-vector interactions, and immune response, and how those provide a roadmap for further research, and, most importantly, point toward how to achieve malaria control and elimination in the Americas.

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.

Changing transmission dynamics among migrant, indigenous and mining populations in a malaria hotspot in Northern Brazil: 2016 to 2020

Background

Roraima state is the northernmost state in Brazil and the primary border-crossing point between Brazil and Venezuela. The uncontrolled surge of malaria in Venezuela, coupled with mass migration of Venezuelans to neighbouring countries and the upward trend in informal mining in the state, pose a serious threat to the broader region, especially to migrant, indigenous and mining populations, jeopardizing malaria elimination efforts. This study describes changes in the epidemiological profile of malaria in Roraima state related to time, place and populations at risk from 2016 to 2020.

Methods

De-identified malaria surveillance data were obtained from the Malaria Epidemiological Surveillance System from 2016 to 2020. Pearson’s chi-square tested differences between imported and autochthonous cases. Multivariable logistic regression was used to identify risk factors for imported versus autochthonous cases by demographic characteristics.

Results

Odds of being an imported case were higher for Plasmodium falciparum cases (AOR = 2.08). However, as the number of cases from Venezuela decreased in 2020 following closure of the border, the proportion of P. falciparum cases increased markedly, from 6.24% in 2019 to 18.50% in 2020. Over the 5-year period, the odds of being an imported case among miners were about nine times higher than the general population (AOR = 8.99). The proportion of total malaria cases that were among indigenous people increased from 33.09% in 2016 to 54.83% in 2020. Indigenous children had a higher burden of malaria with over 40% of cases in children 0 to 9 years old, compared to 8% in non-indigenous children 0 to 9 years old. In some municipalities, place of infection differed from place of notification, with a large proportion of cases in these municipalities reporting in Boa Vista.

Conclusions

Malaria remains a serious threat in Roraima state, especially among high-risk populations, such as miners, migrants, and indigenous people. As malaria cases have increased among indigenous people and miners, and the proportion of P. falciparum cases has increased, elimination efforts require understanding of these risk factors to tailor interventions appropriately. Furthermore, cross-border surveillance systems need to be urgently strengthened at formal and unofficial border points, especially since the border with Venezuela reopened in July 2021.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-022-04141-6.

Bias-corrected maximum-likelihood estimation of multiplicity of infection and lineage frequencies

Background

The UN’s Sustainable Development Goals are devoted to eradicate a range of infectious diseases to achieve global well-being. These efforts require monitoring disease transmission at a level that differentiates between pathogen variants at the genetic/molecular level. In fact, the advantages of genetic (molecular) measures like multiplicity of infection (MOI) over traditional metrics, e.g., R₀, are being increasingly recognized. MOI refers to the presence of multiple pathogen variants within an infection due to multiple infective contacts. Maximum-likelihood (ML) methods have been proposed to derive MOI and pathogen-lineage frequencies from molecular data. However, these methods are biased.

Methods and findings

Based on a single molecular marker, we derive a bias-corrected ML estimator for MOI and pathogen-lineage frequencies. We further improve these estimators by heuristical adjustments that compensate shortcomings in the derivation of the bias correction, which implicitly assumes that data lies in the interior of the observational space. The finite sample properties of the different variants of the bias-corrected estimators are investigated by a systematic simulation study. In particular, we investigate the performance of the estimator in terms of bias, variance, and robustness against model violations. The corrections successfully remove bias except for extreme parameters that likely yield uninformative data, which cannot sustain accurate parameter estimation. Heuristic adjustments further improve the bias correction, particularly for small sample sizes. The bias corrections also reduce the estimators’ variances, which coincide with the Cramér-Rao lower bound. The estimators are reasonably robust against model violations.

Conclusions

Applying bias corrections can substantially improve the quality of MOI estimates, particularly in areas of low as well as areas of high transmission—in both cases estimates tend to be biased. The bias-corrected estimators are (almost) unbiased and their variance coincides with the Cramér-Rao lower bound, suggesting that no further improvements are possible unless additional information is provided. Additional information can be obtained by combining data from several molecular markers, or by including information that allows stratifying the data into heterogeneous groups.

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

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

Malaria in Southern Venezuela: The hottest hotspot in Latin America

Malaria elimination in Latin America is becoming an elusive goal. Malaria cases reached a historical ~1 million in 2017 and 2018, with Venezuela contributing 53% and 51% of those cases, respectively. Historically, malaria incidence in southern Venezuela has accounted for most of the country's total number of cases. The efficient deployment of disease prevention measures and prediction of disease spread to new regions requires an in-depth understanding of spatial heterogeneity on malaria transmission dynamics. Herein, we characterized the spatial epidemiology of malaria in southern Venezuela from 2007 through 2017 and described the extent to which malaria distribution has changed country-wide over the recent years. We found that disease transmission was focal and more prevalent in the southeast region of southern Venezuela where two persistent hotspots of Plasmodium vivax (76%) and P. falciparum (18%) accounted for ~60% of the total number of cases. Such hotspots are linked to deforestation as a consequence of illegal gold mining activities. Incidence has increased nearly tenfold over the last decade, showing an explosive epidemic growth due to a significant lack of disease control programs. Our findings highlight the importance of spatially oriented interventions to contain the ongoing malaria epidemic in Venezuela. This work also provides baseline epidemiological data to assess cross-border malaria dynamics and advocates for innovative control efforts in the Latin American region.

Measuring the contribution of human mobility to malaria persistence

BACKGROUND

To achieve malaria elimination, it is important to determine the role of human mobility in parasite transmission maintenance. The Alto Juruá basin (Brazil) exhibits one of the largest vivax and falciparum malaria prevalence in the Amazon. The goal of this study was to estimate the contribution of human commutes to malaria persistence in this region, using data from an origin-destination survey.

METHODS

Data from an origin-destination survey were used to describe the intensity and motivation for commutations between rural and urban areas in two Alto Juruá basin (Brazil) municipalities, Mâncio Lima and Rodrigues Alves. The relative time-person spent in each locality per household was estimated. A logistic model was developed to estimate the effect of commuting on the probability of contracting malaria for a certain residence zone inhabitant commuting to another zone.

RESULTS

The main results suggest that the assessed population is not very mobile. A total of [Formula: see text] households reported spending over [Formula: see text] of their annual person-hour in areas within the same residence zone. Study and work were the most prevalent commuting motivations, calculated at [Formula: see text] and [Formula: see text] respectively. Spending person-hours in urban Rodrigues Alves conferred relative protection to urban Mâncio Lima residents. The opposite effect was observed for those spending time in rural areas of both municipalities.

CONCLUSION

Residence area is a stronger determinant for contracting malaria than commuting zones in the Alto Juruá region. As these municipalities are a hotspot for Plasmodium transmission, understanding the main local human fluxes is essential for planning control strategies, since the probability of contracting malaria is dependent on the transmission intensity of both the origin and the displacement area. The natural conditions for the circulation of certain pathogens, such as Plasmodium spp., combined with the Amazon human mobility pattern indicate the need for disease control perspective changes. Therefore, intersectoral public policies should become the basis for health mitigation actions.

Balancing selection and high genetic diversity of Plasmodium vivax circumsporozoite central region in parasites from Brazilian Amazon and Rio de Janeiro Atlantic Forest

Circumsporozoite protein (CSP) is the primary pre-erythrocytic vaccine target in Plasmodium species. Knowledge about their genetic diversity can help predict vaccine efficacy and the spread of novel parasite variants. Thus, we investigated pvcsp gene polymorphisms in 219 isolates (136 from Brazilian Amazon [BA], 71 from Rio de Janeiro Atlantic Forest [AF], and 12 from non-Brazilian countries [NB]). Forty-eight polymorphic sites were detected, 46 in the central repeat region (CR), and two in the C-terminal region. Also, the CR presents InDels and a variable number of repeats. All samples correspond to the VK210 variant, and 24 VK210 subtypes based on CR. Nucleotide diversity (π = 0.0135) generated a significant number of haplotypes (168) with low genetic differentiation between the Brazilian regions (Fst = 0.208). The haplotype network revealed similar distances among the BA and AF regions. The linkage disequilibrium indicates that recombination does not seem to be acting in diversity, reinforcing natural selection's role in accelerating adaptive evolution. The high diversity (low Fst) and polymorphism frequencies could be indicators of balancing selection. Although malaria in BA and AF have distinct vector species and different host immune pressures, consistent genetic signature was found in two regions. The immunodominant B-cell epitope mapped in the CR varies from seven to 19 repeats. The CR T-cell epitope is conserved only in 39 samples. Concerning to C-terminal region, the Th2R epitope presented nonsynonymous SNP only in 6% of Brazilian samples, and the Th3R epitope remained conserved in all studied regions. We conclude that, although the uneven distribution of alleles may jeopardize the deployment of vaccines directed to a specific variable locus, a unique vaccine formulation could protect populations in all Brazilian regions.

Identity-by-descent with uncertainty characterises connectivity of Plasmodium falciparum populations on the Colombian-Pacific coast

Characterising connectivity between geographically separated biological populations is a common goal in many fields. Recent approaches to understanding connectivity between malaria parasite populations, with implications for disease control efforts, have used estimates of relatedness based on identity-by-descent (IBD). However, uncertainty around estimated relatedness has not been accounted for. IBD-based relatedness estimates with uncertainty were computed for pairs of monoclonal Plasmodium falciparum samples collected from five cities on the Colombian-Pacific coast where long-term clonal propagation of P. falciparum is frequent. The cities include two official ports, Buenaventura and Tumaco, that are separated geographically but connected by frequent marine traffic. Fractions of highly-related sample pairs (whose classification using a threshold accounts for uncertainty) were greater within cities versus between. However, based on both highly-related fractions and on a threshold-free approach (Wasserstein distances between parasite populations) connectivity between Buenaventura and Tumaco was disproportionally high. Buenaventura-Tumaco connectivity was consistent with transmission events involving parasites from five clonal components (groups of statistically indistinguishable parasites identified under a graph theoretic framework). To conclude, P. falciparum population connectivity on the Colombian-Pacific coast abides by accessibility not isolation-by-distance, potentially implicating marine traffic in malaria transmission with opportunities for targeted intervention. Further investigations are required to test this hypothesis. For the first time in malaria epidemiology (and to our knowledge in ecological and epidemiological studies more generally), we account for uncertainty around estimated relatedness (an important consideration for studies that plan to use genotype versus whole genome sequence data to estimate IBD-based relatedness); we also use threshold-free methods to compare parasite populations and identify clonal components. Threshold-free methods are especially important in analyses of malaria parasites and other recombining organisms with mixed mating systems where thresholds do not have clear interpretation (e.g. due to clonal propagation) and thus undermine the cross-comparison of studies.

Malaria in Venezuela: changes in the complexity of infection reflects the increment in transmission intensity

Background

Malaria incidence has reached staggering numbers in Venezuela. Commonly, Bolívar State accounted for approximately 70% of the country cases every year. Most cases cluster in the Sifontes municipality, a region characterized by an extractive economy, including gold mining. An increase in migration to Sifontes, driven by gold mining, fueled a malaria spillover to the rest of the country and the region. Here samples collected in 2018 were compared with a previous study of 2003/2004 to describe changes in the parasites population structures and the frequency of point mutations linked to anti-malarial drugs.

Methods

A total of 88 Plasmodium falciparum and 94 Plasmodium vivax isolates were collected in 2018 and compared with samples from 2003/2004 (106 P. falciparum and 104 P. vivax). For P. falciparum, mutations linked to drug resistance (Pfdhfr, Pfdhps, and Pfcrt) and the Pfk13 gene associated with artemisinin delayed parasite clearance, were analysed. To estimate the multiplicity of infection (MOI), and perform P. falciparum and P. vivax population genetic analyses, the parasites were genotyped by using eight standardized microsatellite loci.

Results

The P. falciparum parasites are still harbouring drug-resistant mutations in Pfdhfr, Pfdhps, and Pfcrt. However, there was a decrease in the frequency of highly resistant Pfdhps alleles. Mutations associated with artemisinin delayed parasite clearance in the Pfk13 gene were not found. Consistent with the increase in transmission, polyclonal infections raised from 1.9% in 2003/2004 to 39% in 2018 in P. falciparum and from 16.3 to 68% in P. vivax. There is also a decrease in linkage disequilibrium. Bayesian clustering yields two populations linked to the time of sampling, showing that the parasite populations temporarily changed. However, the samples from 2003/2004 and 2018 have several alleles per locus in common without sharing multi-locus genotypes.

Conclusions

The frequency of mutations linked with drug resistance in P. falciparum shows only changes in Pfdhps. Observations presented here are consistent with an increase in transmission from the previously circulating parasites. Following populations longitudinally, using molecular surveillance, provides valuable information in cases such as Venezuela with a fluid malaria situation that is affecting the regional goals toward elimination.

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.

Malaria Molecular Epidemiology: An Evolutionary Genetics Perspective

Malaria is a vector-borne disease that involves multiple parasite species in a variety of ecological settings. However, the parasite species causing the disease, the prevalence of subclinical infections, the emergence of drug resistance, the scale-up of interventions, and the ecological factors affecting malaria transmission, among others, are aspects that vary across areas where malaria is endemic. Such complexities have propelled the study of parasite genetic diversity patterns in the context of epidemiologic investigations. Importantly, molecular studies indicate that the time and spatial distribution of malaria cases reflect epidemiologic processes that cannot be fully understood without characterizing the evolutionary forces shaping parasite population genetic patterns. Although broad in scope, this review in the Microbiology Spectrum Curated Collection: Advances in Molecular Epidemiology highlights the need for understanding population genetic concepts when interpreting parasite molecular data. First, we discuss malaria complexity in terms of the parasite species involved. Second, we describe how molecular data are changing our understanding of malaria incidence and infectiousness. Third, we compare different approaches to generate parasite genetic information in the context of epidemiologically relevant questions related to malaria control. Finally, we describe a few Plasmodium genomic studies as evidence of how these approaches will provide new insights into the malaria disease dynamics. *This article is part of a curated collection.

Ookinete-Specific Genes and 18S SSU rRNA Evidenced in Plasmodium vivax Selection and Adaptation by Sympatric Vectors

In the southern Pacific coast of Chiapas, Mexico (SM), the two most abundant vector species, Nyssorhynchus albimanus and Anopheles pseudopunctipennis, were susceptible to different Plasmodium vivax Pvs25/28 haplotypes. To broaden our understanding of the existing P. vivax in the area, genes encoding proteins relevant for ookinete development and the 18S rRNA were studied. P. vivax infectivity (percentage of infected mosquitoes and oocyst numbers) was evaluated by simultaneously feeding infected blood samples from patients to Ny. albimanus and An. pseudopunctipennis female mosquitoes. Three infectivity patterns were identified: one group of parasites were more infective to An. pseudopunctipennis than to Ny. albimanus, another group was more infective to Ny. albimanus, while a third group infected both vectors similarly. In 29 parasite isolates, the molecular variations of ookinete-specific genes and the 18S rRNA-type S were analyzed. Using concatenated sequences, phylogenetic trees, and Structure analysis, parasite clustering within SM isolates and between these and those from other geographical origins were investigated. A ML phylogenetic tree resolved two parasite lineages: PvSM-A and PvSM-B. They were associated to a different 18S rRNA variant. PvSM-A parasites had 18S rRNA variant rV2 and correspond to parasites causing high oocyst infection in Ny. albimanus. A new ML tree and Structure analysis, both comprising global sequences, showed PvSM-A clustered with Latin American parasites. Meanwhile, all isolates of PvSM-B had 18S rRNA variant rV1 and remained as unique genetic cluster comprising two subgroups: PvSM-Ba, producing high infection in An. pseudopunctipennis, and PvSM-Bb, causing similar oocyst infection in both vector species. PvSM-A parasites were genetically similar to parasites from South America. Meanwhile, PvSM-B were exclusive to southern Mexico and share ancestry with Asian parasites. The results suggest that these lineages evolved separately, likely by geographic and vector restriction.

High-resolution micro-epidemiology of parasite spatial and temporal dynamics in a high malaria transmission setting in Kenya

Novel interventions that leverage the heterogeneity of parasite transmission are needed to achieve malaria elimination. To better understand spatial and temporal dynamics of transmission, we applied amplicon next-generation sequencing of two polymorphic gene regions (csp and ama1) to a cohort identified via reactive case detection in a high-transmission setting in western Kenya. From April 2013 to July 2014, we enrolled 442 symptomatic children with malaria, 442 matched controls, and all household members of both groups. Here, we evaluate genetic similarity between infected individuals using three indices: sharing of parasite haplotypes on binary and proportional scales and the L1 norm. Symptomatic children more commonly share haplotypes with their own household members. Furthermore, we observe robust temporal structuring of parasite genetic similarity and identify the unique molecular signature of an outbreak. These findings of both micro- and macro-scale organization of parasite populations might be harnessed to inform next-generation malaria control measures.

Here, Nelson et al. use amplicon next-generation sequencing of two P. falciparum polymorphic gene regions to investigate the genetic similarity of parasite populations across time and space in a pediatric cohort in Kenya. They identify both micro- and macro-scale structuring of malaria parasites in this high-transmission setting, which could inform future intervention strategies.

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.