Uncovering the transmission dynamics of Plasmodium vivax using population genetics

Characterizing the genetic diversity and population structure of Plasmodium knowlesi in Aceh Province, Indonesia

As in other parts of Southeast Asia, efforts to achieve or sustain malaria elimination in Indonesia have been threatened by the emergence of human infection with the primate species P. knowlesi. To understand the transmission dynamics of this species, investigation of P. knowlesi genetic diversity and population structure is needed. A molecular surveillance study was conducted in two phases between June 2014 and September 2018 at five primary health facilities in Aceh Province, Indonesia, an area nearing malaria elimination. Dried blood spot samples were collected from patients presenting with suspected malaria and testing positive for malaria by microscopy. PCR was performed for molecular confirmation and species identification. Forty-six samples were confirmed to be P. knowlesi, of which 41 were amplified with genotyping targeting ten known P. knowlesi microsatellite markers. For samples within a site, nearly all (9 of 10 loci) or all loci were polymorphic. Across sites, multiple identical haplotypes were observed, though linkage distribution in the population was low (index of association (IAS) = 0.008). The parasite population was indicative of low diversity (expected heterozygosity [HE] =  0.63) and low complexity demonstrated by 92.7% monoclonal infections, a mean multiplicity of infection of 1.06, and a mean within-host infection fixation index (FST) of 0.05. Principal coordinate and neighbour-joining tree analyses indicated that P. knowlesi strains from Aceh were distinct from those reported in Malaysia. In a near-elimination setting in Indonesia, we demonstrate the first evidence that P. knowlesi strains were minimally diverse and were genetically distinct from Malaysian strains, suggesting highly localized transmission and limited connectivity to Malaysia. Ongoing genetic surveillance of P. knowlesi in Indonesia can inform tracking and planning of malaria control and elimination efforts.

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

Introduction

Malaria molecular surveillance (MMS) can provide insights into transmission dynamics, guiding national control programs. We previously designed AmpliSeq assays for MMS, which include different traits of interest (resistance markers and pfhrp2/3 deletions), and SNP barcodes to provide population genetics estimates of Plasmodium vivax and Plasmodium falciparum parasites in the Peruvian Amazon. The present study compares the genetic resolution of the barcodes in the AmpliSeq assays with widely used microsatellite (MS) panels to investigate population genetics of Amazonian malaria parasites.

Methods

We analyzed 51 P. vivax and 80 P. falciparum samples from three distinct areas in the Loreto region of the Peruvian Amazon: Nueva Jerusalén (NJ), Mazan (MZ), and Santa Emilia (SE). Population genetics estimates and costs were compared using the SNP barcodes (P. vivax: 40 SNPs and P. falciparum: 28 SNPs) and MS panels (P. vivax: 16 MS and P. falciparum: 7 MS).

Results

The P. vivax genetic diversity (expected heterozygosity, He) trends were similar for both markers: He MS = 0.68-0.78 (p > 0.05) and He SNP = 0.36-0.38 (p > 0.05). P. vivax pairwise genetic differentiation (fixation index, FST) was also comparable: FST-MS = 0.04-0.14 and FST-SNP = 0.03-0.12 (pairwise p > 0.05). In addition, P. falciparum genetic diversity trends (He MS = 0-0.48, p < 0.05; He SNP = 0-0.09, p < 0.05) and pairwise FST comparisons (FST-MS = 0.14-0.65, FST-SNP = 0.19-0.61, pairwise p > 0.05) were concordant between both panels. For P. vivax, no geographic clustering was observed with any panel, whereas for P. falciparum, similar population structure clustering was observed with both markers, assigning most parasites from NJ to a distinct subpopulation from MZ and SE. We found significant differences in detecting polyclonal infections: for P. vivax, MS identified a higher proportion of polyclonal infections than SNP (69% vs. 33%, p = 3.3 × 10-5), while for P. falciparum, SNP and MS detected similar rates (46% vs. 31%, p = 0.21). The AmpliSeq assay had a higher estimated per-sample cost compared to MS ($183 vs. $27-49).

Discussion

The SNP barcodes in the AmpliSeq assays offered comparable results to MS for investigating population genetics in P. vivax and P. falciparum populations, despite some discrepancies in determining polyclonality. Given both panels have their respective advantages and limitations, the choice between both should be guided by research objectives, costs, and resource availability.

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.

Extensive genetic diversity in Plasmodium vivax from Sudan and its genetic relationships with other geographical isolates

Plasmodium vivax, traditionally overlooked has experienced a notable increase in cases in East Africa. This study investigated the geographical origin and genetic diversity of P. vivax in Sudan using 14 microsatellite markers. A total of 113 clinical P. vivax samples were collected from two different ecogeographical zones, New Halfa and Khartoum, in Sudan. Additionally, 841 geographical samples from the database were incorporated for a global genetic analysis to discern genetic relationships among P. vivax isolates on regional and worldwide scales. On the regional scale, our findings revealed 91 unique and 8 shared haplotypes among the Sudan samples, showcasing a remarkable genetic diversity compared to other geographical isolates and supporting the hypothesis that P. vivax originated from Africa. On a global scale, distinct genetic clustering of P. vivax isolates from Africa, South America, and Asia (including Papua New Guinea and Solomon Island) was observed, with limited admixture among the three clusters. Principal component analysis emphasized the substantial contribution of African isolates to the observed global genetic variation. The Sudanese populations displayed extensive genetic diversity, marked by significant multi-locus linkage disequilibrium, suggesting an ancestral source of P. vivax variation globally and frequent recombination among the isolates. Notably, the East African P. vivax exhibited similarity with some Asian isolates, indicating potential recent introductions. Overall, our results underscore the effectiveness of utilizing microsatellite markers for implementing robust control measures, given their ability to capture extensive genetic diversity and linkage disequilibrium patterns.

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 populations in the western Greater Mekong Subregion evaluated using a genetic barcode

An improved understanding of the Plasmodium vivax populations in the Great Mekong Subregion (GMS) is needed to monitor the progress of malaria elimination. This study aimed to use a P. vivax single nucleotide polymorphism (SNP) barcode to evaluate the population dynamics and explore the gene flow among P. vivax parasite populations in the western GMS (China, Myanmar and Thailand). A total of 315 P. vivax patient samples collected in 2011 and 2018 from four regions of the western GMS were genotyped for 42 SNPs using the high-throughput MassARRAY SNP genotyping technology. Population genetic analysis was conducted to estimate the genetic diversity, effective population size, and population structure among the P. vivax populations. Overall, 291 samples were successfully genotyped at 39 SNPs. A significant difference was observed in the proportion of polyclonal infections among the five P. vivax populations (P = 0.0012, Pearson Chi-square test, χ2 = 18.1), with western Myanmar having the highest proportion (96.2%, 50/52) in 2018. Likewise, the average complexity of infection was also highest in western Myanmar (1.31) and lowest in northeast Myanmar (1.01) in 2018. The older samples from western China in 2011 had the highest pairwise nucleotide diversity (π, 0.388 ± 0.046), expected heterozygosity (He, 0.363 ± 0.02), and the largest effective population size. In comparison, in the neighboring northeast Myanmar, the more recent samples in 2018 showed the lowest values (π, 0.224 ± 0.036; He, 0.220 ± 0.026). Furthermore, the 2018 northeast Myanmar parasites showed high and moderate genetic differentiation from other populations with FST values of 0.162-0.252, whereas genetic differentiation among other populations was relatively low (FST ≤ 0.059). Principal component analysis, phylogeny, and STRUCTURE analysis showed that the P. vivax population in northeast Myanmar in 2018 substantially diverged from other populations. Although the 42 SNP barcode is a valuable tool for tracking parasite origins of worldwide parasite populations, a more extended barcode with additional SNPs is needed to distinguish the more related parasite populations in the western GMS.

Synthesizing environmental, epidemiological and vector and parasite genetic data to assist decision making for disease elimination

We present a framework for identifying when conditions are favourable for transmission of vector-borne diseases between communities by incorporating predicted disease prevalence mapping with landscape analysis of sociological, environmental and host/parasite genetic data. We explored the relationship between environmental features and gene flow of a filarial parasite of humans, Onchocerca volvulus, and its vector, blackflies in the genus Simulium. We generated a baseline microfilarial prevalence map from point estimates from 47 locations in the ecological transition separating the savannah and forest in Ghana, where transmission of O. volvulus persists despite onchocerciasis control efforts. We generated movement suitability maps based on environmental correlates with mitochondrial population structure of 164 parasites from 15 communities and 93 vectors from only four sampling sites, and compared these to the baseline prevalence map. Parasite genetic distance between sampling locations was significantly associated with elevation (r = .793, p = .005) and soil moisture (r = .507, p = .002), while vector genetic distance was associated with soil moisture (r = .788, p = .0417) and precipitation (r = .835, p = .0417). The correlation between baseline prevalence and parasite resistance surface maps was stronger than that between prevalence and vector resistance surface maps. The centre of the study area had high prevalence and suitability for parasite and vector gene flow, potentially contributing to persistent transmission and suggesting the importance of re-evaluating transmission zone boundaries. With suitably dense sampling, this framework can help delineate transmission zones for onchocerciasis and would be translatable to other vector-borne diseases.

Genetic Structure of Introduced Plasmodium vivax Malaria Isolates in Greece, 2015-2019

Greece has been malaria-free since 1974, after an intense malaria control program. However, as Greece hosts migrant populations from P. vivax malaria-endemic countries, there is a risk of introducing the disease to specific vulnerable and receptive areas of the country. Knowledge of the genetic diversity of P. vivax populations is essential for understanding the dynamics of malaria disease transmission in a given region. We used nine highly polymorphic markers to genotype 124 P. vivax-infected archived DNA samples from human blood specimens referred to the NMRL from all over Greece throughout 2015-2019. The genotypic variability of the samples studied was noted, as they comprised several unique haplotypes, indicative of the importation of a large number of different P. vivax strains in the country. However, only a few events of local transmission were recorded. Genotyping revealed and confirmed the same clusters as those identified through epidemiological investigation. In only one introduction event was the index case found. No sustained/ongoing malaria transmissions in/between the studied regions or during consecutive years or additional foci of local transmission were observed. Genotyping is an important component in assisting malaria surveillance, as it provides information concerning the patterns of introduction and the effectiveness of implemented malaria control and elimination measures.

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.

Molecular markers for malaria genetic epidemiology: progress and pitfalls

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

Vivax malaria: a possible stumbling block for malaria elimination in India

Plasmodium vivax is geographically the most widely dispersed human malaria parasite species. It has shown resilience and a great deal of adaptability. Genomic studies suggest that P. vivax originated from Asia or Africa and moved to the rest of the world. Although P. vivax is evolutionarily an older species than Plasmodium falciparum, its biology, transmission, pathology, and control still require better elucidation. P. vivax poses problems for malaria elimination because of the ability of a single primary infection to produce multiple relapses over months and years. P. vivax malaria elimination program needs early diagnosis, and prompt and complete radical treatment, which is challenging, to simultaneously exterminate the circulating parasites and dormant hypnozoites lodged in the hepatocytes of the host liver. As prompt surveillance and effective treatments are rolled out, preventing primaquine toxicity in the patients having glucose-6-phosphate dehydrogenase (G6PD) deficiency should be a priority for the vivax elimination program. This review sheds light on the burden of P. vivax, changing epidemiological patterns, the hurdles in elimination efforts, and the essential tools needed not just in India but globally. These tools encompass innovative treatments for eliminating dormant parasites, coping with evolving drug resistance, and the development of potential vaccines against the parasite.

Heterogeneity in prevalence of subclinical Plasmodium falciparum and Plasmodium vivax infections but no parasite genomic clustering in the Chittagong Hill Tracts, Bangladesh

BACKGROUND

Malaria remains endemic in Bangladesh, with the majority of cases occurring in forested, mountainous region in the Chittagong Hill Tracts (CHT). This area is home to Bengali and diverse groups of indigenous people (Pahari) residing largely in mono-ethnic villages.

METHODS

1002 individuals of the 9 most prominent Pahari and the Bengali population were randomly selected and screened by RDT and qPCR. Parasites were genotyped by msp2 and deep sequencing of 5 amplicons (ama1-D3, cpmp, cpp, csp, and msp7) for Plasmodium falciparum (n = 20), and by microsatellite (MS) typing of ten loci and amplicon sequencing of msp1 for Plasmodium vivax (n = 21). Population structure was analysed using STRUCTURE software. Identity-by-state (IBS) was calculated as a measure of parasite relatedness and used to generate relatedness networks.

RESULTS

The prevalence of P. falciparum and P. vivax infection was 0.7% by RDT (P. falciparum 6/1002; P. vivax 0/1002, mixed: 1/1002) and 4% by qPCR (P. falciparum 21/1002; P. vivax 16/1002, mixed: 5/1002). Infections were highly clustered, with 64% (27/42) of infections occurring in only two Pahari groups, the Khumi and Mro. Diversity was high; expected heterozygosity was 0.93 for P. falciparum and 0.81 for P. vivax. 85.7% (18/21) of P. vivax and 25% (5/20) of P. falciparum infections were polyclonal. No population structure was evident for either species, suggesting high transmission and gene flow among Pahari groups.

CONCLUSIONS

High subclinical infection prevalence and genetic diversity mirror ongoing transmission. Control activities should be specifically directed to Pahari groups at greatest risk.

Spatiotemporal Changes in Plasmodium vivax msp142 Haplotypes in Southern Mexico: From the Control to the Pre-Elimination Phase

For 20 years, Plasmodium vivax has been the only prevalent malaria species in Mexico, and cases have declined significantly and continuously. Spatiotemporal genetic studies can be helpful for understanding parasite dynamics and developing strategies to weaken malaria transmission, thus facilitating the elimination of the parasite. The aim of the current contribution was to analyze P. vivax-infected blood samples from patients in southern Mexico during the control (1993–2007) and pre-elimination phases (2008–2011). Nucleotide and haplotype changes in the pvmsp142 fragment were evaluated over time. The majority of multiple genotype infections occurred in the 1990s, when the 198 single nucleotide sequences exhibited 57 segregating sites, 64 mutations, and 17 haplotypes. Nucleotide and genetic diversity parameters showed subtle fluctuations from across time, in contrast to the reduced haplotype diversity and the increase in the R2 index and Tajima’s D value from 2008 to 2011. The haplotype network consisted of four haplogroups, the geographical distribution of which varied slightly over time. Haplogroup-specific B-cell epitopes were predicted. Since only high-frequency and divergent haplotypes persisted, there was a contraction of the parasite population. Given that 84% of haplotypes were exclusive to Mesoamerica, P. vivax flow is likely circumscribed to this region, representing important information for parasite surveillance.

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.

Genetic characterization of Plasmodium vivax isolates from Pakistan using circumsporozoite protein (pvcsp) and merozoite surface protein-1 (pvmsp-1) genes as genetic markers

Background

Plasmodium vivax contributes to over 70% malaria burden in Pakistan, but limited data exists on various aspects including genetic diversity of the parasite as compared to other parts of the world. Since the information about the genetic diversity of P. vivax assists to understand the population dynamics of the parasite, the current study was designed to understand population divergence of P. vivax in Pakistan using circumsporozoite protein (pvcsp) and merozoite surface protein-1 (pvmsp-1) genes as molecular markers.

Methods

The PCR for pvcsp and pvmsp-1 genes was carried out for 150 P. vivax isolates, followed by DNA sequencing of 35 and 30, respectively. Genetic diversity and polymorphism were analysed using ChromasPro, ClustalW, MEGA7, DnaSP v.5 and WebLogo programs.

Results

The PCR for pvcsp and pvmsp-1 genes was carried out for 150 P. vivax isolates and resulting the PCR products of 1100 bp for pvcsp and ~ 400 bp for pvmsp-1 genes, respectively. In the central-repeat region (CRR) of pvcsp gene, sequences comprised of four variable repeats of PRMs, out of which GDRADGQPA (PRM1), GDRAAGQPA (PRM2) were more extensively dispersed among the P. vivax isolates. Partial sequences (~ 400 bp) of block 2 of pvmsp-1 gene depicted high level of diversity.

Conclusion

The results revealed the polymorphism and genetic diversity especially at the CRR of pvcsp and block 2 of pvmsp-1 genes, respectively. The base-line data presented here warrants future studies to investigate more into the genetic diversity of P. vivax with large sample size from across the country for better understanding of population dynamics of P. vivax that will help to control malaria at individual and community level.

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.

Whole genome sequencing of Plasmodium vivax isolates reveals frequent sequence and structural polymorphisms in erythrocyte binding genes

Plasmodium vivax malaria is much less common in Africa than the rest of the world because the parasite relies primarily on the Duffy antigen/chemokine receptor (DARC) to invade human erythrocytes, and the majority of Africans are Duffy negative. Recently, there has been a dramatic increase in the reporting of P. vivax cases in Africa, with a high number of them being in Duffy negative individuals, potentially indicating P. vivax has evolved an alternative invasion mechanism that can overcome Duffy negativity. Here, we analyzed single nucleotide polymorphism (SNP) and copy number variation (CNV) in Whole Genome Sequence (WGS) data from 44 P. vivax samples isolated from symptomatic malaria patients in southwestern Ethiopia, where both Duffy positive and Duffy negative individuals are found. A total of 123,711 SNPs were detected, of which 22.7% were nonsynonymous and 77.3% were synonymous mutations. The largest number of SNPs were detected on chromosomes 9 (24,007 SNPs; 19.4% of total) and 10 (16,852 SNPs, 13.6% of total). There were particularly high levels of polymorphism in erythrocyte binding gene candidates including merozoite surface protein 1 (MSP1) and merozoite surface protein 3 (MSP3.5, MSP3.85 and MSP3.9). Two genes, MAEBL and MSP3.8 related to immunogenicity and erythrocyte binding function were detected with significant signals of positive selection. Variation in gene copy number was also concentrated in genes involved in host-parasite interactions, including the expansion of the Duffy binding protein gene (PvDBP) on chromosome 6 and MSP3.11 on chromosome 10. Based on the phylogeny constructed from the whole genome sequences, the expansion of these genes was an independent process among the P. vivax lineages in Ethiopia. We further inferred transmission patterns of P. vivax infections among study sites and showed various levels of gene flow at a small geographical scale. The genomic features of P. vivax provided baseline data for future comparison with those in Duffy-negative individuals and allowed us to develop a panel of informative Single Nucleotide Polymorphic markers diagnostic at a micro-geographical scale.

Molecular surveillance over 14 years confirms reduction of Plasmodium vivax and falciparum transmission after implementation of Artemisinin-based combination therapy in Papua, Indonesia

Genetic epidemiology can provide important insights into parasite transmission that can inform public health interventions. The current study compared long-term changes in the genetic diversity and structure of co-endemic Plasmodium falciparum and P. vivax populations. The study was conducted in Papua Indonesia, where high-grade chloroquine resistance in P. falciparum and P. vivax led to a universal policy of Artemisinin-based Combination Therapy (ACT) in 2006. Microsatellite typing and population genetic analyses were undertaken on available isolates collected between 2004 and 2017 from patients with uncomplicated malaria (n = 666 P. falciparum and n = 615 P. vivax). The proportion of polyclonal P. falciparum infections fell from 28% (38/135) before policy change (2004-2006) to 18% (22/125) at the end of the study (2015-2017); p<0.001. Over the same period, polyclonal P. vivax infections fell from 67% (80/119) to 35% (33/93); p<0.001. P. falciparum strains persisted for up to 9 years compared to 3 months for P. vivax, reflecting higher rates of outbreeding in the latter. Sub-structure was observed in the P. falciparum population, but not in P. vivax, confirming different patterns of outbreeding. The P. falciparum population exhibited 4 subpopulations that changed in frequency over time. Notably, a sharp rise was observed in the frequency of a minor subpopulation (K2) in the late post-ACT period, accounting for 100% of infections in late 2016-2017. The results confirm epidemiological evidence of reduced P. falciparum and P. vivax transmission over time. The smaller change in P. vivax population structure is consistent with greater outbreeding associated with relapsing infections and highlights the need for radical cure to reduce recurrent infections. The study emphasizes the challenge in disrupting P. vivax transmission and demonstrates the potential of molecular data to inform on the impact of public health interventions.

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.

Resolving the cause of recurrent Plasmodium vivax malaria probabilistically

Relapses arising from dormant liver-stage Plasmodium vivax parasites (hypnozoites) are a major cause of vivax malaria. However, in endemic areas, a recurrent blood-stage infection following treatment can be hypnozoite-derived (relapse), a blood-stage treatment failure (recrudescence), or a newly acquired infection (reinfection). Each of these requires a different prevention strategy, but it was not previously possible to distinguish between them reliably. We show that individual vivax malaria recurrences can be characterised probabilistically by combined modelling of time-to-event and genetic data within a framework incorporating identity-by-descent. Analysis of pooled patient data on 1441 recurrent P. vivax infections in 1299 patients on the Thailand-Myanmar border observed over 1000 patient follow-up years shows that, without primaquine radical curative treatment, 3 in 4 patients relapse. In contrast, after supervised high-dose primaquine only 1 in 40 relapse. In this region of frequent relapsing P. vivax, failure rates after supervised high-dose primaquine are significantly lower (∼3%) than estimated previously.

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.

Tafenoquine for the prophylaxis, treatment and elimination of malaria: eagerness must meet prudence

Malaria puts more than 3 billion people at risk of infection and causes high morbidity and mortality. Plasmodium vivax forms hypnozoites, which may initiate recurrences, even in the absence of reinfection or superinfection. Until recently, the only drug available for eliminating hypnozoites was primaquine (PQ), which, given its short half-life, requires a relatively long course of treatment. Tafenoquine (TQ) is a PQ analog with a longer half-life. This enables radical cure of malaria with a single dose and overcomes adherence issues associated with PQ, thereby increasing effectiveness in real-life settings. Clinical studies have provided sound evidence for TQ's safety and efficacy against malaria, which recently led to its approval by the US FDA. Here, we review aspects of TQ, including how to avoid hemolytic anemia in G6PD deficient patients. We believe that TQ promises to be a major advance toward malaria elimination.

Molecular detection of drug resistant malaria in Southern Thailand

BACKGROUND

Drug resistance within the major malaria parasites Plasmodium vivax and Plasmodium falciparum threatens malaria control and elimination in Southeast Asia. Plasmodium vivax first-line treatment drug is chloroquine together with primaquine, and the first-line treatment for P. falciparum malaria is artemisinin in combination with a partner drug. Plasmodium vivax and P. falciparum parasites resistant to their respective first-line therapies are now found within Southeast Asia. The resistance perimeters may include high transmission regions of Southern Thailand which are underrepresented in surveillance efforts.

METHODS

This study investigated blood samples from malaria centres in Southern Thailand. Genetic loci associated with drug resistance were amplified and sequenced. Drug resistance associated genes Pvmdr1, Pvcrt-o, Pvdhfr, and Pvdhps were characterized for 145 cases of P. vivax malaria, as well as the artemisinin resistance-associated Pfkelch13 gene from 91 cases of P. falciparum malaria.

RESULTS

Plasmodium vivax samples from Southern Thai provinces showed numerous chloroquine and antifolate resistance-associated mutations, including SNP and Pvcrt-o K10-insertion combinations suggestive of chloroquine resistant P. vivax phenotypes. A high proportion of the C580Y coding mutation (conferring artemisinin resistance) was detected in P. falciparum samples originating from Ranong and Yala (where the mutation was previously unreported).

CONCLUSIONS

The results demonstrate a risk of chloroquine and antifolate resistant P. vivax phenotypes in Southern Thailand, and artemisinin resistant P. falciparum observed as far south as the Thai-Malaysian border region. Ongoing surveillance of antimalarial drug resistance markers is called for in Southern Thailand to inform case management.

Genetic Diversity of Plasmodium vivax Causing Epidemic Malaria in the Republic of Korea

Plasmodium vivax is more challenging to control and eliminate than P. falciparum due to its more asymptomatic infections with low parasite densities making diagnosis more difficult, in addition to its unique biological characteristics. The potential re-introduction of incidence cases, either through borders or via human migrations, is another major hurdle to sustained control and elimination. The Republic of Korea has experienced re-emergence of vivax malaria in 1993 but is one of the 32 malaria-eliminating countries to-date. Despite achieving successful nationwide control and elimination of vivax malaria, the evolutionary characteristics of vivax malaria isolates in the Republic of Korea have not been fully understood. In this review, we present an overview of the genetic variability of such isolates to increase understanding of the epidemiology, diversity, and dynamics of vivax populations in the Republic of Korea.

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.

Paleopathological Considerations on Malaria Infection in Korea before the 20th Century

Malaria, one of the deadliest diseases in human history, still infects many people worldwide. Among the species of the genus Plasmodium, P. vivax is commonly found in temperate-zone countries including South Korea. In this article, we first review the history of malarial infection in Korea by means of studies on Joseon documents and the related scientific data on the evolutionary history of P. vivax in Asia. According to the historical records, malarial infection was not unusual in pre-20th-century Korean society. We also found that certain behaviors of the Joseon people might have affected the host-vector-pathogen relationship, which could explain why malarial infection prevalence was so high in Korea at that time. In our review of genetic studies on P. vivax, we identified substantial geographic differentiation among continents and even between neighboring countries. Based on these, we were able to formulate a strategy for future analysis of ancient Plasmodium strains in Korea.

Analysis of polymorphisms in the circumsporozoite protein gene of Plasmodium vivax isolates from Henan Province, China

BACKGROUND

Plasmodium vivax malaria has historically been a major source of disease in Henan, China. In the 1970s, the morbidity of malaria was highest in the country. With support from the government and the efforts of healthcare personnel, the reported malaria cases have declined dramatically and a national elimination programme was launched in 2010. To achieve the goal, it is essential to study the diversity of autochthonous malaria and transmission of Plasmodium parasites, which will provide baseline data for disease control and management.

METHODS

Thirty-two P. vivax isolates from Henan province were collected from 2008 to 2011, and circumsporozoite protein (csp) genes were analysed to estimate the genetic diversity of this parasite.

RESULTS

The assessment of csp sequences indicated that all the isolates were the VK210 type, however, none of them was identical to the VK210 strain. The sequences displayed variations in the central region, and eight sub-types were observed. Among the sub-types, HN7 was the most prevalent (37.5%), followed by HN3 (34.4%). A total of 653 repeat units were discovered in 32 Henan isolates. Nucleotide sequences were grouped in 13 unique repeat nucleotide sequence allotypes that coded for 7 different repeated amino acid allotypes. B (GNGAGGQAA) and D (GDRAAGQPA) were more frequent based on the results; they represented 53.9% (352/653) of the total. In comparison to the basic repeat units of VK210, more than 75% of the central repeat units had at least one non-synonymous nucleotide change.

CONCLUSIONS

Recent P. vivax populations in Henan province showed some degree of genetic diversity in csp, with 8 sub-types among 32 samples. Meantime, the results also suggested its relative conserved parasite populations. This could provide interesting baseline data that allow identifying whether potential new cases differ from the parasites already circulating in the area.

Plasmodium vivax Biology: Insights Provided by Genomics, Transcriptomics and Proteomics

During the last decade, the vast omics field has revolutionized biological research, especially the genomics, transcriptomics and proteomics branches, as technological tools become available to the field researcher and allow difficult question-driven studies to be addressed. Parasitology has greatly benefited from next generation sequencing (NGS) projects, which have resulted in a broadened comprehension of basic parasite molecular biology, ecology and epidemiology. Malariology is one example where application of this technology has greatly contributed to a better understanding of Plasmodium spp. biology and host-parasite interactions. Among the several parasite species that cause human malaria, the neglected Plasmodium vivax presents great research challenges, as in vitro culturing is not yet feasible and functional assays are heavily limited. Therefore, there are gaps in our P. vivax biology knowledge that affect decisions for control policies aiming to eradicate vivax malaria in the near future. In this review, we provide a snapshot of key discoveries already achieved in P. vivax sequencing projects, focusing on developments, hurdles, and limitations currently faced by the research community, as well as perspectives on future vivax malaria research.

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.

Passively versus Actively Detected Malaria: Similar Genetic Diversity but Different Complexity of Infection

The surveillance of malaria is generally undertaken on the assumption that samples passively collected at health facilities are comparable to or representative of the broader Plasmodium reservoir circulating in the community. Further characterization and comparability of the hidden asymptomatic parasite reservoir are needed to inform on the potential impact of sampling bias. This study explores the impact of sampling strategy on molecular surveillance by comparing the genetic make-up of Plasmodium falciparum and Plasmodium vivax isolates collected by passive versus active case detection. Sympatric isolates of P. falciparum and P. vivax were collected from a large community survey and ongoing clinical surveillance studies undertaken in the hypomesoendemic setting of Mimika District (Papua, Indonesia). Plasmodium falciparum isolates were genotyped at nine microsatellite loci and P. vivax at eight loci. Measures of diversity and differentiation were used to compare different patient and parasitological sample groups. The results demonstrated that passively detected cases (symptomatic) had comparable population diversity to those circulating in the community (asymptomatic) in both species. In addition, asymptomatic patent infections were as diverse as subpatent infections. However, a significant difference in multiplicity of infection (MOI) and percentage of polyclonal infections was observed between actively and passively detected P. vivax cases (mean MOI: 1.7 ± 0.7 versus 1.4 ± 1.4, respectively; P = 0.001). The study findings infer that, in hypomesoendemic settings, passive sampling is appropriate for molecular parasite surveillance strategies using the predominant clone in any given infection; however, the findings suggest caution when analyzing complexity of infection. Further evaluation is required in other endemic settings.

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.

Plasmodium vivax genetic diversity and heterozygosity in blood samples and resulting oocysts at the Thai-Myanmar border

Background

Polyclonal blood-stage infections of Plasmodium vivax are frequent even in low transmission settings, allowing meiotic recombination between heterologous parasites. Empirical data on meiotic products are however lacking. This study examined microsatellites in oocysts derived by membrane feeding of mosquitoes from blood-stage P. vivax infections at the Thai–Myanmar border.

Methods

Blood samples from patients presenting with vivax malaria were fed to Anopheles cracens by membrane feeding and individual oocysts from midguts were obtained by dissection after 7 days. DNA was extracted from oocysts and parental blood samples and tested by microsatellite analysis.

Results

A focused study of eight microsatellite markers was undertaken for nine blood stage infections from 2013, for which derived oocysts were studied in six cases. One or more alleles were successfully amplified for 131 oocysts, revealing high levels of allelic diversity in both blood and oocyst stages. Based on standard criteria for defining minor alleles, there was evidence of clear deviation from random mating (inbreeding) with relatively few heterozygous oocysts compared to variance across the entire oocyst population (F_IT = 0.89). The main explanation appeared to be natural compartmentalisation at mosquito (F_SC = 0.27) and human stages (F_CT = 0.68). One single human case produced a total of 431 successfully amplified loci (across 70 oocysts) that were homozygous and identical to parental alleles at all markers, indicating clonal infection and transmission. Heterozygous oocyst alleles were found at 15/176 (8.5%) successfully amplified loci in the other five cases. There was apparently reduced oocyst heterozygosity in individual oocysts compared to diversity within individual mosquitoes (F_IS = 0.55), but this may simply reflect the difficulty of detecting minor alleles in oocysts, given the high rate of amplification failure. Inclusion of minor allele peaks (irrespective of height) when matching peaks were found in related blood or oocyst samples, added 11 minor alleles for 9 oocysts, increasing the number of heterozygous loci to 26/176 (14.8%; p = 0.096).

Conclusion

There was an apparently low level of heterozygous oocysts but this can be explained by a combination of factors: relatively low complexity of parental infection, natural compartmentalisation in humans and mosquitoes, and the methodological challenge of detecting minor alleles.

Electronic supplementary material

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

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.

Transmission dynamics of co-endemic Plasmodium vivax and P. falciparum in Ethiopia and prevalence of antimalarial resistant genotypes

Ethiopia is one of the few African countries where Plasmodium vivax is co-endemic with P. falciparum. Malaria transmission is seasonal and transmission intensity varies mainly by landscape and climate. Although the recent emergence of drug resistant parasites presents a major issue to malaria control in Ethiopia, little is known about the transmission pathways of parasite species and prevalence of resistant markers. This study used microsatellites to determine population diversity and gene flow patterns of P. falciparum (N = 226) and P. vivax (N = 205), as well as prevalence of drug resistant markers to infer the impact of gene flow and existing malaria treatment regimes. Plasmodium falciparum indicated a higher rate of polyclonal infections than P. vivax. Both species revealed moderate genetic diversity and similar population structure. Populations in the northern highlands were closely related to the eastern Rift Valley, but slightly distinct from the southern basin area. Gene flow via human migrations between the northern and eastern populations were frequent and mostly bidirectional. Landscape genetic analyses indicated that environmental heterogeneity and geographical distance did not constrain parasite gene flow. This may partly explain similar patterns of resistant marker prevalence. In P. falciparum, a high prevalence of mutant alleles was detected in codons related to chloroquine (pfcrt and pfmdr1) and sulfadoxine-pyrimethamine (pfdhps and pfdhfr) resistance. Over 60% of the samples showed pfmdr1 duplications. Nevertheless, no mutation was detected in pfK13 that relates to artemisinin resistance. In P. vivax, while sequences of pvcrt-o were highly conserved and less than 5% of the samples showed pvmdr duplications, over 50% of the samples had pvmdr1 976F mutation. It remains to be tested if this mutation relates to chloroquine resistance. Monitoring the extent of malaria spread and markers of drug resistance is imperative to inform policy for evidence-based antimalarial choice and interventions. To effectively reduce malaria burden in Ethiopia, control efforts should focus on seasonal migrant populations.

Chloroquine efficacy for Plasmodium vivax in Myanmar in populations with high genetic diversity and moderate parasite gene flow

BACKGROUND

Plasmodium vivax malaria remains a major public health burden in Myanmar. Resistance to chloroquine (CQ), the first-line treatment for P. vivax, has been reported in the country and has potential to undermine local control efforts.

METHODS

Patients over 6 years of age with uncomplicated P. vivax mono-infection were enrolled into clinical efficacy studies in Myawaddy in 2014 and Kawthoung in 2012. Study participants received a standard dose of CQ (25 mg/kg over 3 days) followed by weekly review until day 28. Pvmdr1 copy number (CN) and microsatellite diversity were assessed on samples from the patients enrolled in the clinical study and additional cross-sectional surveys undertaken in Myawaddy and Shwegyin in 2012.

RESULTS

A total of 85 patients were enrolled in the CQ clinical studies, 25 in Myawaddy and 60 in Kawthoung. One patient in Myawaddy (1.2%) had an early treatment failure and two patients (2.3%) in Kawthoung presented with late treatment failures on day 28. The day 28 efficacy was 92.0% (95% CI 71.6-97.9) in Myawaddy and 98.3% (95% CI 88.7-99.8) in Kawthoung. By day 2, 92.2% (23/25) in Myawaddy and 85.0% (51/60) in Kawthoung were aparasitaemic. Genotyping and pvmdr1 CN assessment was undertaken on 43, 52 and 46 clinical isolates from Myawaddy, Kawthoung and Shwegyin respectively. Pvmdr1 amplification was observed in 3.2% (1/31) of isolates in Myawaddy, 0% (0/49) in Kawthoung and 2.5% (1/40) in Shwegyin. Diversity was high in all sites (H _E 0.855-0.876), with low inter-population differentiation (F _ST 0.016-0.026, P < 0.05).

CONCLUSIONS

Treatment failures after chloroquine were observed following chloroquine monotherapy, with pvmdr1 amplification present in both Myawaddy and Shwegyin. The results emphasize the importance of ongoing P. vivax drug resistance surveillance in Myanmar, particularly given the potential connectivity between parasite population at different sites.

Variable number of tandem repeats of 9 Plasmodium vivax genes among Southeast Asian isolates

The variable number of tandem repeats (VNTRs) provides valuable information about both the functional and evolutionary aspects of genetic diversity. Comparative analysis of 3 Plasmodium falciparum genomes has shown that more than 9% of its open reading frames (ORFs) harbor VNTRs. Although microsatellites and VNTR genes of P. vivax were reported, the VNTR polymorphism of genes has not been examined widely. In this study, 230 P. vivax genes were analyzed for VNTRs by SERV, and 33 kinds of TR deletions or insertions from 29 P. vivax genes (12.6%) were found. Of these, 9 VNTR fragments from 8 P. vivax genes were used for PCR amplification and sequence analysis to examine the genetic diversity among 134 isolates from four Southeast Asian countries (China, Republic of Korea, Thailand, and Myanmar) with different malaria endemicity. We confirmed the existence of extensive polymorphism of VNTR fragments in field isolates. This detection provides several suitable markers for analysis of the molecular epidemiology of P. vivax field isolates.

Genetic micro-epidemiology of malaria in Papua Indonesia: Extensive P. vivax diversity and a distinct subpopulation of asymptomatic P. falciparum infections

BACKGROUND

Genetic analyses of Plasmodium have potential to inform on transmission dynamics, but few studies have evaluated this on a local spatial scale. We used microsatellite genotyping to characterise the micro-epidemiology of P. vivax and P. falciparum diversity to inform malaria control strategies in Timika, Papua Indonesia.

METHODS

Genotyping was undertaken on 713 sympatric P. falciparum and P. vivax isolates from a cross-sectional household survey and clinical studies conducted in Timika. Standard population genetic measures were applied, and the data was compared to published data from Kalimantan, Bangka, Sumba and West Timor.

RESULTS

Higher diversity (HE = 0.847 vs 0.625; p = 0.017) and polyclonality (46.2% vs 16.5%, p<0.001) were observed in P. vivax versus P. falciparum. Distinct P. falciparum substructure was observed, with two subpopulations, K1 and K2. K1 was comprised solely of asymptomatic infections and displayed greater relatedness to isolates from Sumba than to K2, possibly reflecting imported infections.

CONCLUSIONS

The results demonstrate the greater refractoriness of P. vivax versus P. falciparum to control measures, and risk of distinct parasite subpopulations persisting in the community undetected by passive surveillance. These findings highlight the need for complimentary new surveillance strategies to identify transmission patterns that cannot be detected with traditional malariometric methods.

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.

Higher Complexity of Infection and Genetic Diversity of Plasmodium vivax Than Plasmodium falciparum Across All Malaria Transmission Zones of Papua New Guinea

Plasmodium falciparum and Plasmodium vivax have varying transmission dynamics that are informed by molecular epidemiology. This study aimed to determine the complexity of infection and genetic diversity of P. vivax and P. falciparum throughout Papua New Guinea (PNG) to evaluate transmission dynamics across the country. In 2008-2009, a nationwide malaria indicator survey collected 8,936 samples from all 16 endemic provinces of PNG. Of these, 892 positive P. vivax samples were genotyped at PvMS16 and PvmspF3, and 758 positive P. falciparum samples were genotyped at Pfmsp2. The data were analyzed for multiplicity of infection (MOI) and genetic diversity. Overall, P. vivax had higher polyclonality (71%) and mean MOI (2.32) than P. falciparum (20%, 1.39). These measures were significantly associated with prevalence for P. falciparum but not for P. vivax. The genetic diversity of P. vivax (PvMS16: expected heterozygosity = 0.95, 0.85-0.98; PvMsp1F3: 0.78, 0.66-0.89) was higher and less variable than that of P. falciparum (Pfmsp2: 0.89, 0.65-0.97). Significant associations of MOI with allelic richness (rho = 0.69, P = 0.009) and expected heterozygosity (rho = 0.87, P < 0.001) were observed for P. falciparum. Conversely, genetic diversity was not correlated with polyclonality nor mean MOI for P. vivax. The results demonstrate higher complexity of infection and genetic diversity of P. vivax across the country. Although P. falciparum shows a strong association of these parameters with prevalence, a lack of association was observed for P. vivax and is consistent with higher potential for outcrossing of this species.

VivaxGEN: An open access platform for comparative analysis of short tandem repeat genotyping data in Plasmodium vivax populations

BACKGROUND

The control and elimination of Plasmodium vivax will require a better understanding of its transmission dynamics, through the application of genotyping and population genetics analyses. This paper describes VivaxGEN (http://vivaxgen.menzies.edu.au), a web-based platform that has been developed to support P. vivax short tandem repeat data sharing and comparative analyses.

RESULTS

The VivaxGEN platform provides a repository for raw data generated by capillary electrophoresis (FSA files), with fragment analysis and standardized allele calling tools. The query system of the platform enables users to filter, select and differentiate samples and alleles based on their specified criteria. Key population genetic analyses are supported including measures of population differentiation (FST), expected heterozygosity (HE), linkage disequilibrium (IAS), neighbor-joining analysis and Principal Coordinate Analysis. Datasets can also be formatted and exported for application in commonly used population genetic software including GENEPOP, Arlequin and STRUCTURE. To date, data from 10 countries, including 5 publicly available data sets have been shared with VivaxGEN.

CONCLUSIONS

VivaxGEN is well placed to facilitate regional overviews of P. vivax transmission dynamics in different endemic settings and capable to be adapted for similar genetic studies of P. falciparum and other organisms.

Overlap Extension Barcoding for the Next Generation Sequencing and Genotyping of Plasmodium falciparum in Individual Patients in Western Kenya

Large-scale molecular epidemiologic studies of Plasmodium falciparum parasites have provided insights into parasite biology and transmission, can identify the spread of drug resistance, and are useful in assessing vaccine targets. The polyclonal nature infections in high transmission settings is problematic for traditional genotyping approaches. Next-generation sequencing (NGS) approaches to parasite genotyping allow sensitive detection of minority variants, disaggregation of complex parasite mixtures, and scalable processing of large samples sets. Therefore, we designed, validated, and applied to field parasites an approach that leverages sequencing of individually barcoded samples in a multiplex manner. We utilize variant barcodes, invariant linker sequences and modular template-specific primers to allow for the simultaneous generation of high-dimensional sequencing data of multiple gene targets. This modularity permits a cost-effective and reproducible way to query many genes at once. In mixtures of reference parasite genomes, we quantitatively detected unique haplotypes comprising as little as 2% of a polyclonal infection. We applied this genotyping approach to field-collected parasites collected in Western Kenya in order to simultaneously obtain parasites genotypes at three unlinked loci. In summary, we present a rapid, scalable, and flexible method for genotyping individual parasites that enables molecular epidemiologic studies of parasite evolution, population structure and transmission.

Molecular Epidemiology of P. vivax in Iran: High Diversity and Complex Sub-Structure Using Neutral Markers, but No Evidence of Y976F Mutation at pvmdr1

BACKGROUND

Malaria remains endemic at low levels in the south-eastern provinces of Iran bordering Afghanistan and Pakistan, with the majority of cases attributable to P. vivax. The national guidelines recommend chloroquine (CQ) as blood-stage treatment for uncomplicated P. vivax, but the large influx of imported cases enhances the risk of introducing CQ resistance (CQR).

METHODOLOGY AND PRINCIPAL FINDINGS

The genetic diversity at pvmdr1, a putative modulator of CQR, and across nine putatively neutral short tandem repeat (STR) markers were assessed in P. vivax clinical isolates collected between April 2007 and January 2013 in Hormozgan Province, south-eastern Iran. One hundred blood samples were collected from patients with microscopy-confirmed P. vivax enrolled at one of five district clinics. In total 73 (73%) were autochthonous cases, 23 (23%) imported cases from Afghanistan or Pakistan, and 4 (4%) with unknown origin. 97% (97/100) isolates carried the F1076L mutation, but none carried the Y976F mutation. STR genotyping was successful in 71 (71%) isolates, including 57(57%) autochthonous and 11 (11%) imported cases. Analysis of population structure revealed 2 major sub-populations, K1 and K2, with further sub-structure within K2. The K1 sub-population had markedly lower diversity than K2 (HE = 0.06 vs HE = 0.82) suggesting that the sub-populations were sustained by distinct reservoirs with differing transmission dynamics, possibly reflecting local versus imported/introduced populations. No notable separation was observed between the local and imported cases although the sample size was limited.

CONCLUSIONS

The contrasting low versus high diversity in the two sub-populations (K1 and K2) infers that a combination of local transmission and cross-border malaria from higher transmission regions shape the genetic make-up of the P. vivax population in south-eastern Iran. There was no molecular evidence of CQR amongst the local or imported cases, but ongoing clinical surveillance is warranted.

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.

Implications of Plasmodium vivax Biology for Control, Elimination, and Research

This paper summarizes our current understanding of the biology of Plasmodium vivax, how it differs from Plasmodium falciparum, and how these differences explain the need for P. vivax-tailored interventions. The article further pinpoints knowledge gaps where investments in research are needed to help identify and develop such specific interventions. The principal obstacles to reduce and eventually eliminate P. vivax reside in 1) its higher vectorial capacity compared with P. falciparum due to its ability to develop at lower temperature and over a shorter sporogonic cycle in the vector, allowing transmission in temperate zones and making it less sensitive to vector control measures that are otherwise effective on P. falciparum; 2) the presence of dormant liver forms (hypnozoites), sustaining multiple relapsing episodes from a single infectious bite that cannot be diagnosed and are not susceptible to any available antimalarial except primaquine, with routine deployment restricted by toxicity; 3) low parasite densities, which are difficult to detect with current diagnostics leading to missed diagnoses and delayed treatments (and protracted transmission), coupled with 4) transmission stages (gametocytes) occurring early in acute infections, before infection is diagnosed.

Genomic analysis of local variation and recent evolution in Plasmodium vivax

The widespread distribution and relapsing nature of Plasmodium vivax infection present major challenges for the elimination of malaria. To characterize the genetic diversity of this parasite in individual infections and across the population, we performed deep genome sequencing of >200 clinical samples collected across the Asia-Pacific region and analyzed data on >300,000 SNPs and nine regions of the genome with large copy number variations. Individual infections showed complex patterns of genetic structure, with variation not only in the number of dominant clones but also in their level of relatedness and inbreeding. At the population level, we observed strong signals of recent evolutionary selection both in known drug resistance genes and at new loci, and these varied markedly between geographical locations. These findings demonstrate a dynamic landscape of local evolutionary adaptation in the parasite population and provide a foundation for genomic surveillance to guide effective strategies for control and elimination of P. vivax.

Where chloroquine still works: the genetic make-up and susceptibility of Plasmodium vivax to chloroquine plus primaquine in Bhutan

BACKGROUND

Bhutan has made substantial progress in reducing malaria incidence. The national guidelines recommend chloroquine (CQ) and primaquine (PQ) for radical cure of uncomplicated Plasmodium vivax, but the local efficacy has not been assessed. The impact of cases imported from India on the genetic make-up of the local vivax populations is currently unknown.

METHODS

Patients over 4 years of age with uncomplicated P. vivax mono-infection were enrolled into a clinical efficacy study and molecular survey. Study participants received a standard dose of CQ (25 mg/kg over 3 days) followed by weekly review until day 28. On day 28 a 14-day regimen of PQ (0.25 mg/kg/day) was commenced under direct observation. After day 42, patients were followed up monthly for a year. The primary and secondary endpoints were risk of treatment failure at day 28 and at 1 year. Parasite genotyping was undertaken at nine tandem repeat markers, and standard population genetic metrics were applied to examine population diversity and structure in infections thought to be acquired inside or outside of Bhutan.

RESULTS

A total of 24 patients were enrolled in the clinical study between April 2013 and October 2015. Eight patients (33.3 %) were lost to follow-up in the first 6 months and another eight patients lost between 6 and 12 months. No (0/24) treatment failures occurred by day 28 and no (0/8) parasitaemia was detected following PQ treatment. Some 95.8 % (23/24) of patients were aparasitaemic by day 2. There were no haemolytic or serious events. Genotyping was undertaken on parasites from 12 autochthonous cases and 16 suspected imported cases. Diversity was high (H E 0.87 and 0.90) in both populations. There was no notable differentiation between the autochthonous and imported populations.

CONCLUSIONS

CQ and PQ remains effective for radical cure of P. vivax in Bhutan. The genetic analyses indicate that imported infections are sustaining the local vivax population, with concomitant risk of introducing drug-resistant strains.

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.

Diversity and evolutionary genetics of the three major Plasmodium vivax merozoite genes participating in reticulocyte invasion in southern Mexico

Background

Reported malaria cases in the Americas had been reduced to about one-half million by 2012. To advance towards elimination of this disease, it is necessary to gain insights into how the malaria parasite is evolving, including the emergence, spread and persistence of new haplotypes in affected regions. In here, the genetic diversity of the three major P. vivax merozoite genes was analyzed.

Methods

From P. vivax-infected blood samples obtained in southern Mexico (SMX) during 2006–2007, nucleotide sequences were achieved for: the 42 kDa carboxyl fragment of the merozoite surface protein-1 (msp1₄₂), domains I-II of the apical membrane antigen-1 (ama1_I-II), and domain II of the Duffy binding protein (dbp_II). Gene polymorphism was examined and haplotype networks were developed to depict parasite relationships in SMX. Then genetic diversity, recombination and natural selection were analyzed and the degree of differentiation was determined as F_ST values.

Results

The diversity of P. vivax merozoite genes in SMX was less than that of parasites from other geographic origins, with dbp_II < ama1_I-II < msp1₄₂. Ama1_I-II and msp1₄₂ exposed the more numerous haplotypes exclusive to SMX. While, all dbp_II haplotypes from SMX were separated from one to three mutational steps, the networks of ama1_I-II and msp1₄₂ were more complex; loops and numerous mutational steps were evidenced, likely due to recombination. Sings of local diversification were more evident for msp1₄₂. Sixteen combined haplotypes were determined; one of these haplotypes not detected in 2006 was highly frequent in 2007. The Rm value was higher for msp1₄₂than for ama1_I-II, being insignificant for dbp_II. The dN-dS value was highly significant for ama1_I-II and lesser so for dbp_II. The F_ST values were higher for dbp_II than msp1₄₂, and very low for ama1_I-II.

Conclusions

In SMX, P. vivax ama1_I-II, dbp_II and msp1₄₂ demonstrated limited diversity, and exhibited a differentiated parasite population. The results suggest that differential intensities of selective forces are operating on these gene fragments, and probably related to their timing, length of exposure and function during reticulocyte adhesion and invasion. Therefore, these finding are essential for mono and multivalent vaccine development and for epidemiological surveillance.

Electronic supplementary material

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