Population structure and spatio-temporal transmission dynamics of Plasmodium vivax after radical cure treatment in a rural village of the Peruvian Amazon

Spatiotemporal dynamics of malaria and climate influence on its incidence in Condorcanqui Province, 2005-2022

BACKGROUND

Amazonas is a region in northern Peru with the second-highest incidence of malaria. Approximately 95% of the cases are reported in the Condorcanqui province, where native communities living along the banks of Santiago River lack access to potable water, sewage, and electricity. This study aimed to analyse malaria's spatial, temporal, and climatic characteristics in Condorcanqui to guide future studies and prevention strategies.

METHODS

A database provided by DIRESA-Amazonas was evaluated. Database included cases from 44 health facilities serving 112 native communities. According to the malaria control programs implemented in Peru, the study was divided into three periods: 2005-2010, 2011-2016, and 2017-2022. A Spearman correlation analysis was also conducted to assess the relationship between malaria incidence and climate variables.

RESULTS

During the study periods, 10,632 cases were reported, including Plasmodium vivax (84.87%), Plasmodium falciparum (14.91%) and Plasmodium malariae (0.23%) infections. Annual incidence rates (AIRs) significantly varied across the study periods (p < 0.001). A significant reduction in malaria incidence occurred during the first period, largely attributed to PAMAFRO programme interventions. Subsequent periods, showed a gradual increase in cases, with a peak of P. vivax in 2019 and the reintroduction of P. falciparum. Males and individuals aged 0-11 years presented the greatest number of cases. Significant correlations were found between malaria incidence and the Oceanic Niño Index (ONI) at lag0 (ρ = 0.14, p = 0.037), corrected precipitation at lag1 (ρ = 0.16, p = 0.020), and minimum wind speed at lag1 (ρ = 0.15, p = 0.024).

CONCLUSIONS

Malaria incidence in Condorcanqui has increased over the last 5 years, driven by climatic influences such as the ONI, precipitation, and low wind speeds. Without immediate preventive efforts, cases are expected to continue rising. Effective control strategies must tackle the social, economic, and political issues that heighten vulnerability, such as poverty and limited healthcare access. Maintaining control initiatives and tailoring them to local needs will be essential for achieving long-term reductions of malaria in Peru.

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[identity by descent]>0.99), copy number variation in Pvdbp and multiple Pvrbps, mutations in antimalarial resistance orthologs, 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.

Detection through the use of RT-MqPCR of asymptomatic reservoirs of malaria in samples of patients from the indigenous Comarca of Guna Yala, Panama: Essential method to achieve the elimination of malaria

BACKGROUND

Plasmodium vivax is the main causative agent of malaria in Panama. However, the prevalence of asymptomatic infections in the different endemic regions remains unknown. Understanding the epidemiological behavior of asymptomatic infections is essential for the elimination of malaria. This study aimed to determine the prevalence of asymptomatic malarial infections in one of the main endemic regions of Panama using multiplex real-time reverse transcription RT-MqPCR.

METHODS

A cross-sectional study was conducted in three communities in the Guna Yala Comarca. A total of 551 thick blood smears and their respective samples on filter paper were collected from volunteers of different ages and sexes from June 20 to 25, 2016. Infections by the Plasmodium spp. were diagnosed using microscopy and RT-MqPCR. All statistical analyses were performed using the R software.

RESULTS

The average prevalence of asymptomatic infections by P. vivax in the three communities detected by RT-MqPCR was 9.3%, with Ukupa having the highest prevalence (13.4%), followed by Aidirgandi (11.1%) and Irgandi (3.3%). A total of 74 samples were diagnosed as asymptomatic infections using RT-MqPCR. Light microscopy (LM) detected that 17.6% (13/74) of the asymptomatic samples and 82.4% (61/74) were diagnosed as false negatives. A 100% correlation was observed between samples diagnosed using LM and RT-MqPCR. A total of 52.7% (39/74) of the asymptomatic patients were female and 85.1% (63/74) were registered between the ages of 1 and 21 years. Factors associated with asymptomatic infection were community (aOR = 0.38 (95% CI 0.17-0.83), p < 0.001) and age aOR = 0.98 (95% CI 0.97-1.00), p < 0.05); F = 5.38; p < 0.05).

CONCLUSIONS

This study provides novel evidence of the considerable prevalence of asymptomatic P. vivax infections in the endemic region of Kuna Yala, representing a new challenge that requires immediate attention from the National Malaria Program. The results of this study provide essential information for the health authorities responsible for developing new policies. Furthermore, it will allow program administrators to reorient and design effective malaria control strategies that consider asymptomatic infections as a fundamental part of malaria control and move towards fulfilling their commitment to eliminate it.

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.

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.

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.

Effectiveness of a Malaria Surveillance Strategy Based on Active Case Detection during High Transmission Season in the Peruvian Amazon

Background: Faced with the resurgence of malaria, malaria surveillance in the Peruvian Amazon incorporated consecutive active case detection (ACD) interventions using light microscopy (LM) as reactive measure in communities with an unusual high number of cases during high transmission season (HTS). We assessed the effectiveness in malaria detection of this local ACD-based strategy. Methods: A cohort study was conducted in June–July 2015 in Mazan, Loreto. Four consecutive ACD interventions at intervals of 10 days were conducted in four riverine communities (Gamitanacocha, Primero de Enero, Libertad and Urco Miraño). In each intervention, all inhabitants were visited at home, and finger-prick blood samples collected for immediate diagnosis by LM and on filter paper for later analysis by quantitative real-time polymerase chain reaction (qPCR). Effectiveness was calculated by dividing the number of malaria infections detected using LM by the number of malaria infections detected by delayed qPCR. Results: Most community inhabitants (88.1%, 822/933) were present in at least one of the four ACD interventions. A total of 451 infections were detected by qPCR in 446 participants (54.3% of total participants); five individuals had two infections. Plasmodium vivax was the predominant species (79.8%), followed by P. falciparum (15.3%) and P. vivax-P. falciparum co-infections (4.9%). Most qPCR-positive infections were asymptomatic (255/448, 56.9%). The ACD-strategy using LM had an effectiveness of 22.8% (detection of 103 of the total qPCR-positive infections). Children aged 5–14 years, and farming as main economic activity were associated with P. vivax infections. Conclusions: Although the ACD-strategy using LM increased the opportunity of detecting and treating malaria infections during HTS, the number of detected infections was considerably lower than the real burden of infections (those detected by qPCR).

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.

Nationwide genetic surveillance of Plasmodium vivax in Papua New Guinea reveals heterogeneous transmission dynamics and routes of migration amongst subdivided populations

The Asia Pacific Leaders in Malaria Alliance (APLMA) have committed to eliminate malaria from the region by 2030. Papua New Guinea (PNG) has the highest malaria burden in the Asia-Pacific region but with the intensification of control efforts since 2005, transmission has been dramatically reduced and Plasmodium vivax is now the dominant malaria infection in some parts of the country. To gain a better understanding of the transmission dynamics and migration patterns of P. vivax in PNG, here we investigate population structure in eight geographically and ecologically distinct regions of the country. A total of 219 P. vivax isolates (16-30 per population) were successfully haplotyped using 10 microsatellite markers. A wide range of genetic diversity (H_e=0.37-0.87, R_s=3.60-7.58) and significant multilocus linkage disequilibrium (LD) was observed in six of the eight populations (I_A^S=0.08-0.15 p-value<0.05) reflecting a spectrum of transmission intensities across the country. Genetic differentiation between regions was evident (Jost's D=0.07-0.72), with increasing divergence of populations with geographic distance. Overall, P. vivax isolates clustered into three major genetic populations subdividing the Mainland lowland and coastal regions, the Islands and the Highlands. P. vivax gene flow follows major human migration routes, and there was higher gene flow amongst Mainland parasite populations than among Island populations. The Central Province (samples collected in villages close to the capital city, Port Moresby), acts as a sink for imported infections from the three major endemic areas. These insights into P. vivax transmission dynamics and population networks will inform targeted strategies to contain malaria infections and to prevent the spread of drug resistance in PNG.

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.

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.

Temporal genetic changes in Plasmodium vivax apical membrane antigen 1 over 19 years of transmission in southern Mexico

Background

Mexico advanced to the pre-elimination phase in 2009 due to a significant reduction in malaria cases, and since 2000, Plasmodium vivax is the only species transmitted. During the last two decades, malaria transmission has been mostly local and isolated to a few regions. It is important to gain further insights into the impact of control measures on the parasite population structure. Hence, the aim of the current study was to determine detailed changes in P. vivax genetic diversity and population structure based on analysing the gene that encodes the apical membrane antigen 1 (pvama1). This analysis covered from control to pre-elimination (1993–2011) in a hypo-endemic region in southern Mexico.

Results

The 213 pvama1_I-II sequences presently analysed were grouped into six periods of three years each. They showed low genetic diversity, with 15 haplotypes resolved. Among the DNA sequences, there was a gradual decrease in genetic diversity, the number of mixed genotype infections and the intensity of positive selection, in agreement with the parallel decline in malaria cases. At the same time, linkage disequilibrium (R²) increased. The three-dimensional haplotype network revealed that pvama1_I-II haplotypes were separated by 1–11 mutational steps, and between one another by 0–3 unsampled haplotypes. In the temporal network, seven haplotypes were detected in at least two of the six-time layers, and only four distinct haplotypes were evidenced in the pre-elimination phase. Structure analysis indicated that three subpopulations fluctuated over time. Only 8.5% of the samples had mixed ancestry. In the pre-elimination phase, subpopulation P1 was drastically reduced, and the admixture was absent.

Conclusions

The results suggest that P. vivax in southern Mexico evolved based on local adaptation into three “pseudoclonal” subpopulations that diversified at the regional level and persisted over time, although with varying frequency. Control measures and climate events influenced the number of malaria cases and the genetic structure. The sharp decrease in parasite diversity and other related genetic parameters during the pre-elimination phase suggests that malaria elimination is possible in the near future. These results are useful for epidemiological surveillance.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-017-2156-y) contains supplementary material, which is available to authorized users.

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.

Molecular evidence of high rates of asymptomatic P. vivax infection and very low P. falciparum malaria in Botswana

BACKGROUND

Botswana is one of eight SADC countries targeting malaria elimination by 2018. Through spirited upscaling of control activities and passive surveillance, significant reductions in case incidence of Plasmodium falciparum (0.96 - 0.01) was achieved between 2008 and 2012. As part of the elimination campaign, active detection of asymptomatic Plasmodium species by a highly sensitive method was deemed necessary. This study was carried out to determine asymptomatic Plasmodium species carriage by nested PCR in the country, in 2012.

METHOD

A cross-sectional study involving 3924 apparently healthy participants were screened for Plasmodium species in 14 districts (5 endemic: Okavango, Ngami, Tutume, Boteti and Bobirwa; and 9 epidemic: North East, Francistown, Serowe-Palapye, Ghanzi, Kweneng West, Kweneng East, Kgatleng, South East, and Good Hope). Venous blood was taken from each participant for a nested PCR detection of Plasmodium species.

RESULTS

The parasite rates of asymptomatic Plasmodium species detected were as follows: Plasmodium falciparum, 0.16 %; Plasmodium vivax, 4.66 %; Plasmodium malariae, (Pm) 0.16 %; Plasmodium ovale, 0 %, mixed infections (P. falciparum and P. vivax), 0.055 %; and (P. vivax and P. malariae), 0.027 %, (total: 5.062 %). The high proportion of asymptomatic reservoir of P. vivax was clustered in the East, South Eastern and Central districts of the country. There appeared to be a correlation between the occurrence of P. malariae infection with P. vivax infection, with the former only occurring in districts that had substantial P. vivax circulation. The median age among 2-12 year olds for P. vivax infection was 5 years (Mean 5.13 years, interquartile range 3-7 years). The odds of being infected with P. vivax decreased by 7 % for each year increase in age (OR 0.93, 95 % CI 0.87-1.00, p = 0.056).

CONCLUSION

We have confirmed low parasite rate of asymptomatic Plasmodium species in Botswana, with the exception of P.vivax which was unexpectedly high. This has implication for the elimination campaign so a follow up study is warranted to inform decisions on new strategies that take this evidence into account in the elimination campaign.

Prospective Study of Plasmodium vivax Malaria Recurrence after Radical Treatment with a Chloroquine-Primaquine Standard Regimen in Turbo, Colombia

Plasmodium vivax recurrences help maintain malaria transmission. They are caused by recrudescence, reinfection, or relapse, which are not easily differentiated. A longitudinal observational study took place in Turbo municipality, Colombia. Participants with uncomplicated P. vivax infection received supervised treatment concomitantly with 25 mg/kg chloroquine and 0.25 mg/kg/day primaquine for 14 days. Incidence of recurrence was assessed over 180 days. Samples were genotyped, and origins of recurrences were established. A total of 134 participants were enrolled between February 2012 and July 2013, and 87 were followed for 180 days, during which 29 recurrences were detected. The cumulative incidence of first recurrence was 24.1% (21/87) (95% confidence interval [CI], 14.6 to 33.7%), and 86% (18/21) of these events occurred between days 51 and 110. High genetic diversity of P. vivax strains was found, and 12.5% (16/128) of the infections were polyclonal. Among detected recurrences, 93.1% (27/29) of strains were genotyped as genetically identical to the strain from the previous infection episode, and 65.5% (19/29) of infections were classified as relapses. Our results indicate that there is a high incidence of P. vivax malaria recurrence after treatment in Turbo municipality, Colombia, and that a large majority of these episodes are likely relapses from the previous infection. We attribute this to the primaquine regimen currently used in Colombia, which may be insufficient to eliminate hypnozoites.

Uncomplicated Clinical Malaria Features, the Efficacy of Artesunate-Amodiaquine and Their Relation with Multiplicity of Infection in the Democratic Republic of Congo

Background

In the Democratic Republic of Congo, artesunate-amodiaquine (ASAQ) is the first-line medication recommended for uncomplicated malaria treatment. We conducted a study in Kinshasa to describe the clinical features of the disease and assess the efficacy of ASAQ and its impact on the multiplicity of infection in children with uncomplicated malaria.

Methods

Children aged 12 to 59 months with uncomplicated P. falciparum malaria were treated with ASAQ and followed up passively for 42 days. To distinguish new infections from recrudescent parasites, samples were genotyped using a stepwise strategy with three molecular markers (GLURP, MSP2 and MSP1). We then assessed PCR-corrected and -uncorrected day-42 cure rates and multiplicity of infection (MOI).

Results

In total, 2,796 patients were screened and 865 enrolled in the study. Clinical features were characterized by history of fever (100%), coryza (59.9%) and weakness (59.4%). The crude and PCR-corrected efficacies of ASAQ were 55.3% (95%CI: 51.8–58.8) and 92.8% (95%CI: 91.0–94.6) respectively, as 83.6% (95%CI: 79.1–87.2) of the recurrences were new infections. Compared to monoclonal infections, polyclonal infections were more frequent at enrollment (88.1%) and in recurrences (80.1%; p = 0.005; OR: 1.8, 95%CI: 1.20–2.8). The median MOI at enrollment (MOI = 3.7; IQR: 0.7–6.7) decreased to 3 (IQR: 1–5) in the recurrent samples (p<0.001). Patients infected with a single haplotype on day 0 had no recrudescence; the risk of recrudescence increased by 28% with each additional haplotype (HR: 1.3, 95%CI: 1.24–1.44).

Conclusion

The PCR-corrected efficacy of ASAQ at day 42 was 92.8%, but crude efficacy was relatively poor due to high reinfection rates. Treatment outcomes were positively correlated with MOI. Continued monitoring of the efficacy of ACTs—ASAQ, in this case—is paramount.

Trial Registration

ClinicalTrials.gov NCT01374581

Accelerating to Zero: Strategies to Eliminate Malaria in the Peruvian Amazon

AbstractIn February 2014, the Malaria Elimination Working Group, in partnership with the Peruvian Ministry of Health (MoH), hosted its first international conference on malaria elimination in Iquitos, Peru. The 2-day meeting gathered 85 malaria experts, including 18 international panelists, 23 stakeholders from different malaria-endemic regions of Peru, and 11 MoH authorities. The main outcome was consensus that implementing a malaria elimination project in the Amazon region is achievable, but would require: 1) a comprehensive strategic plan, 2) the altering of current programmatic guidelines from control toward elimination by including symptomatic as well as asymptomatic individuals for antimalarial therapy and transmission-blocking interventions, and 3) the prioritization of community-based active case detection with proper rapid diagnostic tests to interrupt transmission. Elimination efforts must involve key stakeholders and experts at every level of government and include integrated research activities to evaluate, implement, and tailor sustainable interventions appropriate to the region.

Population Genetics of Plasmodium vivax in Four Rural Communities in Central Vietnam

Background

The burden of malaria in Vietnam has drastically reduced, prompting the National Malaria Control Program to officially engage in elimination efforts. Plasmodium vivax is becoming increasingly prevalent, remaining a major problem in the country's central and southern provinces. A better understanding of P. vivax genetic diversity and structure of local parasite populations will provide baseline data for the evaluation and improvement of current efforts for control and elimination. The aim of this study was to examine the population genetics and structure of P. vivax isolates from four communities in Tra Leng commune, Nam Tra My district in Quang Nam, Central Vietnam.

Methodology/Principal Findings

P. vivax mono infections collected from 234 individuals between April 2009 and December 2010 were successfully analyzed using a panel of 14 microsatellite markers. Isolates displayed moderate genetic diversity (He = 0.68), with no significant differences between study communities. Polyclonal infections were frequent (71.4%) with a mean multiplicity of infection of 1.91 isolates/person. Low but significant genetic differentiation (F_ST value from -0.05 to 0.18) was observed between the community across the river and the other communities. Strong linkage disequilibrium (IAS = 0.113, p < 0.001) was detected across all communities, suggesting gene flow within and among them. Using multiple approaches, 101 haplotypes were grouped into two genetic clusters, while 60.4% of haplotypes were admixed.

Conclusions/Significance

In this area of Central Vietnam, where malaria transmission has decreased significantly over the past decade, there was moderate genetic diversity and high occurrence of polyclonal infections. Local human populations have frequent social and economic interactions that facilitate gene flow and inbreeding among parasite populations, while decreasing population structure. Findings provide important information on parasites populations circulating in the study area and are relevant to current malaria elimination efforts.

In Vietnam, Plasmodium vivax (P. vivax) is the second most frequent human malaria parasite and a major obstacle to countrywide malaria elimination. Knowing the local parasite structure is useful for elimination efforts. Therefore, we analyzed, with a panel of 14 microsatellite markers, 234 P. vivax mono infections in blood samples collected from 4 communities in central Vietnam. Genetic diversity in the population was moderate; a high occurrence of polyclonal infections and significant linkage disequilibrium were detected, suggesting inbreeding or recombination between highly related haplotypes. In addition, both genetic differentiation and population structure was low and only detected between communities at each side of the river. Those results suggest gene flow between study communities with the river defining a moderate geographical barrier. Future studies should determine how this genetic variation is maintained in an area of extremely low transmission.

Population Genetics of Plasmodium vivax in the Peruvian Amazon

Background

Characterizing the parasite dynamics and population structure provides useful information to understand the dynamic of transmission and to better target control interventions. Despite considerable efforts for its control, vivax malaria remains a major health problem in Peru. In this study, we have explored the population genetics of Plasmodium vivax isolates from Iquitos, the main city in the Peruvian Amazon, and 25 neighbouring peri-urban as well as rural villages along the Iquitos-Nauta Road.

Methodology/ Results

From April to December 2008, 292 P. vivax isolates were collected and successfully genotyped using 14 neutral microsatellites. Analysis of the molecular data revealed a similar proportion of monoclonal and polyclonal infections in urban areas, while in rural areas monoclonal infections were predominant (p = 0.002). Multiplicity of infection was higher in urban (MOI = 1.5–2) compared to rural areas (MOI = 1) (p = 0.003). The level of genetic diversity was similar in all areas (He = 0.66–0.76, p = 0.32) though genetic differentiation between areas was substantial (PHI_PT = 0.17, p<0.0001). Principal coordinate analysis showed a marked differentiation between parasites from urban and rural areas. Linkage disequilibrium was detected in all the areas (IAs = 0.08–0.49, for all p<0.0001). Gene flow among the areas was stablished through Bayesian analysis of migration models. Recent bottleneck events were detected in 4 areas and a recent parasite expansion in one of the isolated areas. In total, 87 unique haplotypes grouped in 2 or 3 genetic clusters described a sub-structured parasite population.

Conclusion/Significance

Our study shows a sub-structured parasite population with clonal propagation, with most of its components recently affected by bottleneck events. Iquitos city is the main source of parasite spreading for all the peripheral study areas. The routes of transmission and gene flow and the reduction of the parasite population described are important from the public health perspective as well for the formulation of future control policies.

We present the population genetics of malaria vivax parasites in a large area of the Peruvian Amazon. Our results showed that the parasite population had a predominant clonal propagation, reproducing themselves with identically or closely related parasites; therefore, the same genetic characteristics are maintained in the offspring. The clonal propagation may favour the higher levels of genetic differentiation among the parasites from isolated areas compared to areas where human migration is common. The patterns of gene flow have been established, finding Iquitos city as a reservoir of parasite genetic variability. Moreover, a recent reduction of the parasite population was observed in areas where recent control activities were performed. This research provides a picture of the nature and dynamics of the parasite population which have a significant impact in the malaria epidemiology; therefore, this knowledge is crucial for the development of efficient control policies.

Therapeutic efficacy of alternative primaquine regimens to standard treatment in preventing relapses by Plasmodium vivax: A systematic review and meta-analysis

OBJECTIVE

To compare efficacy and safety of primaquine regimens currently used to prevent relapses by P. vivax.

METHODS

A systematic review was carried out to identify clinical trials evaluating efficacy and safety to prevent malaria recurrences by P. vivax of primaquine regimen 0.5 mg/kg/ day for 7 or 14 days compared to standard regimen of 0.25 mg/kg/day for 14 days. Efficacy of primaquine according to cumulative incidence of recurrences after 28 days was determined. The overall relative risk with fixed-effects meta-analysis was estimated.

RESULTS

For the regimen 0.5 mg/kg/day/7 days were identified 7 studies, which showed an incidence of recurrence between 0% and 20% with follow-up 60-210 days; only 4 studies comparing with the standard regimen 0.25 mg/kg/day/14 days and no difference in recurrences between both regimens (RR= 0.977, 95% CI= 0.670 to 1.423) were found. 3 clinical trials using regimen 0.5 mg/kg/day/14 days with an incidence of recurrences between 1.8% and 18.0% during 330-365 days were identified; only one study comparing with the standard regimen (RR= 0.846, 95% CI= 0.484 to 1.477). High risk of bias and differences in handling of included studies were found.

CONCLUSION

Available evidence is insufficient to determine whether currently PQ regimens used as alternative rather than standard treatment have better efficacy and safety in preventing relapse of P. vivax. Clinical trials are required to guide changes in treatment regimen of malaria vivax.

Effectiveness of combined chloroquine and primaquine treatment in 14 days versus intermittent single dose regimen, in an open, non-randomized, clinical trial, to eliminate Plasmodium vivax in southern Mexico

Background

In Mexico, combined chloroquine (CQ) and primaquine (PQ) treatment has been used since the late 1950s to treat Plasmodium vivax infections. Although malaria transmission has declined, current treatment strategies must be evaluated to advance towards malaria elimination.

Methods

The clinical and parasitological outcome of treating symptomatic P. vivax with the 14-day (T14) treatment or intermittent single dose (ISD) regimen was evaluated in southern Mexico between February 2008 and September 2010. Patients over 12 months old with P. vivax mono-infection and asexual parasitaemia ≥500 parasites/µl were treated under supervision. After diagnosis (day 0), treatment began immediately. T14 patients received CQ for 3 days (10, 10 and 5 mg/kg) and PQ daily for 14 days (0.25 mg/kg), while ISD patients received a single dose of CQ (10 mg/kg) and PQ (0.75 mg/kg) on days 0, 30, 60, 180, 210, and 240. Follow-up was done by observing clinical and laboratory (by microscopy, serology and PCR) outcome, considering two endpoints: primary blood infection clearance and clinical response at ~28 days, and the incidence of recurrent blood infection during 12 months. Parasite genotypes of primary/recurrent blood infections were analysed.

Results

During the first 28 days, no differences in parasite clearance or clinical outcome were observed between T14 (86 patients) and ISD (67 patients). On day 3, 95 % of patients in both groups showed no blood parasites, and no recurrences were detected on days 7–28. Contrarily, the therapeutic effectiveness (absence of recurrent parasitaemia) was distinct for T14 versus ISD at 12 months: 83.7 versus 50 %, respectively (p = 0.000). Symptomatic and asymptomatic infections were recorded on days 31–352. Some parasite recurrences were detected by PCR and/or serological testing.

Conclusions

T14 was effective for opportune elimination of the primary blood infection and preventing relapse episodes. The first single dose of CQ-PQ eliminated primary blood infection as efficiently as the initial three-dose scheme of T14, but the ISD regimen should be abandoned. A single combined dose administered to symptomatic patients in remote areas while awaiting parasitological diagnosis may contribute to halting P. vivax transmission. Alternatives for meeting the challenge of T14 supervision are discussed.

Trial registration: NIH-USA, ClinicalTrial.gov Identifier: NCT02394197

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-015-0938-2) contains supplementary material, which is available to authorized users.

Assessment of an automated capillary system for Plasmodium vivax microsatellite genotyping

BACKGROUND

Several platforms have been used to generate the primary data for microsatellite analysis of malaria parasite genotypes. Each has relative advantages but share a limitation of being time- and cost-intensive. A commercially available automated capillary gel cartridge system was assessed in the microsatellite analysis of Plasmodium vivax diversity in the Peruvian Amazon.

METHODS

The reproducibility and accuracy of a commercially-available automated capillary system, QIAxcel, was assessed using a sequenced PCR product of 227 base pairs. This product was measured 42 times, then 27 P. vivax samples from Peruvian Amazon subjects were analyzed with this instrument using five informative microsatellites. Results from the QIAxcel system were compared with a Sanger-type sequencing machine, the ABI PRISM(®) 3100 Genetic Analyzer.

RESULTS

Significant differences were seen between the sequenced amplicons and the results from the QIAxcel instrument. Different runs, plates and cartridges yielded significantly different results. Additionally, allele size decreased with each run by 0.045, or 1 bp, every three plates. QIAxcel and ABI PRISM systems differed in giving different values than those obtained by ABI PRISM, and too many (i.e. inaccurate) alleles per locus were also seen with the automated instrument.

CONCLUSIONS

While P. vivax diversity could generally be estimated using an automated capillary gel cartridge system, the data demonstrate that this system is not sufficiently precise for reliably identifying parasite strains via microsatellite analysis. This conclusion reached after systematic analysis was due both to inadequate precision and poor reproducibility in measuring PCR product size.

Plasmodium vivax Diversity and Population Structure across Four Continents

Plasmodium vivax is the geographically most widespread human malaria parasite. To analyze patterns of microsatellite diversity and population structure across countries of different transmission intensity, genotyping data from 11 microsatellite markers was either generated or compiled from 841 isolates from four continents collected in 1999–2008. Diversity was highest in South-East Asia (mean allelic richness 10.0–12.8), intermediate in the South Pacific (8.1–9.9) Madagascar and Sudan (7.9–8.4), and lowest in South America and Central Asia (5.5–7.2). A reduced panel of only 3 markers was sufficient to identify approx. 90% of all haplotypes in South Pacific, African and SE-Asian populations, but only 60–80% in Latin American populations, suggesting that typing of 2–6 markers, depending on the level of endemicity, is sufficient for epidemiological studies. Clustering analysis showed distinct clusters in Peru and Brazil, but little sub-structuring was observed within Africa, SE-Asia or the South Pacific. Isolates from Uzbekistan were exceptional, as a near-clonal parasite population was observed that was clearly separated from all other populations (F_ST>0.2). Outside Central Asia F_ST values were highest (0.11–0.16) between South American and all other populations, and lowest (0.04–0.07) between populations from South-East Asia and the South Pacific. These comparisons between P. vivax populations from four continents indicated that not only transmission intensity, but also geographical isolation affect diversity and population structure. However, the high effective population size results in slow changes of these parameters. This persistency must be taken into account when assessing the impact of control programs on the genetic structure of parasite populations.

Plasmodium vivax is the predominant malaria parasite in Latin America, Asia and the South Pacific. Different factors are expected to shape diversity and population structure across continents, e.g. transmission intensity which is much lower in South America as compared to Southeast-Asia and the South Pacific, or geographical isolation of P. vivax populations in the South Pacific. We have compiled data from 841 isolates from South and Central America, Africa, Central Asia, Southeast-Asia and the South Pacific typed with a panel of 11 microsatellite markers. Diversity was highest in Southeast-Asia, where transmission is intermediate-high and migration of infected hosts is high, and lowest in South America and Central Asia where malaria transmission is low and focal. Reducing the panel of microsatellites showed that 2–6 markers are sufficient for genotyping for most drug trials and epidemiological studies, as these markers can identify >90% of all haplotypes. Parasites clustered according to continental origin, with high population differentiation between South American and Central Asian populations and the other populations, and lowest differences between Southeast-Asia and the South Pacific. Current attempts to reduce malaria transmission might change this pattern, but only after transmission is reduced for an extended period of time.

Development of a single nucleotide polymorphism barcode to genotype Plasmodium vivax infections

Plasmodium vivax, one of the five species of Plasmodium parasites that cause human malaria, is responsible for 25–40% of malaria cases worldwide. Malaria global elimination efforts will benefit from accurate and effective genotyping tools that will provide insight into the population genetics and diversity of this parasite. The recent sequencing of P. vivax isolates from South America, Africa, and Asia presents a new opportunity by uncovering thousands of novel single nucleotide polymorphisms (SNPs). Genotyping a selection of these SNPs provides a robust, low-cost method of identifying parasite infections through their unique genetic signature or barcode. Based on our experience in generating a SNP barcode for P. falciparum using High Resolution Melting (HRM), we have developed a similar tool for P. vivax. We selected globally polymorphic SNPs from available P. vivax genome sequence data that were located in putatively selectively neutral sites (i.e., intergenic, intronic, or 4-fold degenerate coding). From these candidate SNPs we defined a barcode consisting of 42 SNPs. We analyzed the performance of the 42-SNP barcode on 87 P. vivax clinical samples from parasite populations in South America (Brazil, French Guiana), Africa (Ethiopia) and Asia (Sri Lanka). We found that the P. vivax barcode is robust, as it requires only a small quantity of DNA (limit of detection 0.3 ng/μl) to yield reproducible genotype calls, and detects polymorphic genotypes with high sensitivity. The markers are informative across all clinical samples evaluated (average minor allele frequency > 0.1). Population genetic and statistical analyses show the barcode captures high degrees of population diversity and differentiates geographically distinct populations. Our 42-SNP barcode provides a robust, informative, and standardized genetic marker set that accurately identifies a genomic signature for P. vivax infections.

Plasmodium vivax malaria is a major global public health problem, with nearly 2.5 billion people at risk for infection and approximately 132–391 million clinical infections annually. It has a wide geographical range, with a high disease burden in Asia, Central and South America, the Middle East, Oceania, and East Africa. Advances in sequencing technology and sample processing have made it possible to characterize the genetic diversity of P. vivax populations. This genetic variation provides a means to identify parasites by unique genetic signatures, or “barcodes.” We developed such a genetic barcode for P. vivax, composed of 42 robust and informative variants. Here we report its development and validation based on 87 clinical samples identified by microscopy to contain P. vivax from geographically diverse parasite populations from South America (Brazil, French Guiana), Africa (Ethiopia) and Asia (Sri Lanka). We show that the SNP barcode provides a genotyping tool that can be performed at low cost, providing a means to uniquely identify parasite infections and distinguish geographic origins, and that barcode data may offer new insights into P. vivax population structure and diversity.