High genetic diversity of Plasmodium vivax on the north coast of Papua New Guinea

Genetic diversity and molecular evolution of Plasmodium vivax Duffy Binding Protein and Merozoite Surface Protein-1 in northwestern Thailand

The local diversity and population structure of malaria parasites vary across different regions of the world, reflecting variations in transmission intensity, host immunity, and vector species. This study aimed to use amplicon sequencing to investigate the genotypic patterns and population structure of P. vivax isolates from a highly endemic province of Thailand in recent years. Amplicon deep sequencing was performed on 70 samples for the 42-kDa region of pvmsp1 and domain II of pvdbp. Unique haplotypes were identified and a network constructed to illustrate genetic relatedness in northwestern Thailand. Based on this dataset of 70 samples collected between 2015 and 2021, 16 and 40 unique haplotypes were identified in pvdbpII and pvmsp142kDa, respectively. Nucleotide diversity was higher in pvmsp142kDa than in pvdbpII (π = 0.027 and 0.012), as was haplotype diversity (Hd = 0.962 and 0.849). pvmsp142kDa also showed a higher recombination rate and higher levels of genetic differentiation (Fst) in northwestern Thailand versus other regions (0.2761-0.4881). These data together suggested that the genetic diversity of P. vivax in northwestern Thailand at these two studied loci evolved under a balancing selection, most likely host immunity. The lower genetic diversity of pvdbpII may reflect its stronger functional constrain. In addition, despite the balancing selection, a decrease in genetic diversity was observed. Hd of pvdbpII decreased from 0.874 in 2015-2016 to 0.778 in 2018-2021; π of pvmsp142kDa decreased from 0.030 to 0.022 over the same period. Thus, the control activities must have had a strong impact on the parasite population size. The findings from this study provide an understanding of P. vivax population structure and the evolutionary force on vaccine candidates. They also established a new baseline for tracking future changes in P. vivax diversity in the most malarious area of Thailand.

Measurably recombining malaria parasites

Genomic epidemiology has guided research and policy for various viral pathogens and there has been a parallel effort towards using genomic epidemiology to combat diseases that are caused by eukaryotic pathogens, such as the malaria parasite. However, the central concept of viral genomic epidemiology, namely that of measurably mutating pathogens, does not apply easily to sexually recombining parasites. Here we introduce the related but different concept of measurably recombining malaria parasites to promote convergence around a unifying theoretical framework for malaria genomic epidemiology. Akin to viral phylodynamics, we anticipate that an inferential framework developed around recombination will help guide practical research and thus realize the full public health potential of genomic epidemiology for malaria parasites and other sexually recombining pathogens.

Plasmodium vivax Genetic Diversity in Panama: Challenges for Malaria Elimination in Mesoamerica

Panama and all nations within the Mesoamerican region have committed to eliminate malaria within this decade. With more than 90% of the malaria cases in this region caused by Plasmodium vivax, an efficient national/regional elimination plan must include a comprehensive study of this parasite's genetic diversity. Here, we retrospectively analyzed P. vivax genetic diversity in autochthonous and imported field isolates collected in different endemic regions in Panama from 2007 to 2020, using highly polymorphic markers (csp, msp-1, and msp-3α). We did the analysis using molecular techniques that are cost-effective for malaria molecular surveillance within Mesoamerica. Thus, we used molecular analyses that are feasible for malaria molecular surveillance within the region, and that can provide useful information for policy and decision making about malaria elimination. We also evaluated if haplotypes established by combining the genotypes found in these genes were associated with relevant epidemiological variables and showed structure across the transmission foci that have been observed in Panama. Ten different haplotypes were identified, some of them strongly associated with geographical origin, age, and collection year. Phylogenetic analysis of csp (central repeat domain) revealed that both major variant types (vk210 and vk247) were circulating in Panama. Variant vk247 was restricted to the eastern endemic regions, while vk210 was predominant (77.3%) and widespread, displaying higher diversity (14 alleles) and geographically biased alleles. The regional implications of these molecular findings for the control of P. vivax malaria to achieve elimination across Mesoamerica are discussed.

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.

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.

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

Background

The scale-up of effective malaria control in the last decade has resulted in a substantial decline in the incidence of clinical malaria in many countries. The effects on the proportions of asymptomatic and submicroscopic infections and on transmission potential are yet poorly understood.

Methods

In Papua New Guinea, vector control has been intensified since 2008, and improved diagnosis and treatment was introduced in 2012. Cross-sectional surveys were conducted in Madang Province in 2006 (with 1280 survey participants), 2010 (with 2117 participants), and 2014 (with 2516 participants). Infections were quantified by highly sensitive quantitative polymerase chain reaction (PCR) analysis, and gametocytes were quantified by reverse-transcription qPCR analysis.

Results

Plasmodium falciparum prevalence determined by qPCR decreased from 42% in 2006 to 9% in 2014. The P. vivax prevalence decreased from 42% in 2006 to 13% in 2010 but then increased to 20% in 2014. Parasite densities decreased 5-fold from 2006 to 2010; 72% of P. falciparum and 87% of P. vivax infections were submicroscopic in 2014. Gametocyte density and positivity correlated closely with parasitemia, and population gametocyte prevalence decreased 3-fold for P. falciparum and 29% for P. vivax from 2010 to 2014.

Conclusions

Sustained control has resulted in reduced malaria transmission potential, but an increasing proportion of gametocyte carriers are asymptomatic and submicroscopic and represent a challenge to malaria control.

Cost-effectiveness of artemisinin-naphthoquine versus artemether-lumefantrine for the treatment of uncomplicated malaria in Papua New Guinean children

BACKGROUND

A recent randomized trial showed that artemisinin-naphthoquine (AN) was non-inferior to artemether-lumefantrine (AL) for falciparum malaria and superior for vivax malaria in young Papua New Guinean children. The aim of this study was to compare the cost-effectiveness of these two regimens.

METHODS

An incremental cost-effectiveness analysis was performed using data from 231 children with Plasmodium falciparum and/or Plasmodium vivax infections in an open-label, randomized, parallel-group trial. Recruited children were randomized 1:1 to receive once daily AN for 3 days with water or twice daily AL for 3 days given with fat. World Health Organisation (WHO) definitions were used to determine clinical/parasitological outcomes. The cost of transport between the home and clinic, plus direct health-care costs, served as a basis for determining each regimen's incremental cost per incremental treatment success relative to AL by Day 42 and its cost per life year saved.

RESULTS

In the usual care setting, AN was more effective for the treatment of uncomplicated malaria in children aged 0.5-5.9 years. AL and AN were equally efficacious for the treatment of falciparum malaria, however AN had increased anti-malarial treatment costs per patient of $10.46, compared with AL. AN was the most effective regimen for treatment of vivax malaria, but had increased treatment costs of $14.83 per treatment success compared with AL.

CONCLUSIONS

Whilst AN has superior overall efficacy for the treatment of uncomplicated malaria in PNG children, AL was the less costly regimen. An indicative extrapolation estimated the cost per life year saved by using AN instead of AL to treat uncomplicated malaria to be $12,165 for girls and $12,469 for boys (discounted), which means AN may not be cost-effective and affordable for PNG at current cost. However, AN may become acceptable should it become WHO prequalified and/or should donated/subsidized drug supply become available.

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.

Behaviour and molecular identification of Anopheles malaria vectors in Jayapura district, Papua province, Indonesia

BACKGROUND

Members of the Anopheles punctulatus group dominate Papua, Indonesia and Papua New Guinea (PNG), with a geographic range that extends south through Vanuatu. An. farauti and An. punctulatus are the presumed major vectors in this region. Although this group of species has been extensively studied in PNG and the southern archipelagoes within their range, their distribution, ecology and vector behaviours have not been well characterized in eastern Indonesia.

METHODS

Mosquitoes were collected in five villages in Jayapura province, Papua, Indonesia using human-landing collections, animal-baited tents and backpack aspirators. Mosquitoes were morphologically typed and then molecularly distinguished based on ribosomal ITS2 sequences and tested for Plasmodium falciparum and P. vivax infection using circumsporozoite ELISA and PCR.

RESULTS

The presence and vector status of An. farauti 4 in Papua, Indonesia is confirmed here for the first time. The data indicate that this species is entering houses at a rate that increases its potential to come into contact with humans and act as a major malaria vector. An. farauti 4 was also abundant outdoors and biting humans during early evening hours. Other species collected in this area include An. farauti 1, An. hinesorum, An. koliensis, An. punctulatus, and An. tessellatus. Proboscis morphology was highly variable within each species, lending support to the notion that this characteristic is not a reliable indicator to distinguish species within the An. punctulatus group.

CONCLUSIONS

The vector composition in Papua, Indonesia is consistent with certain northern areas of PNG, but the behaviours of anophelines sampled in this region, such as early and indoor human biting of An. farauti 4, may enable them to act as major vectors of malaria. Presumed major vectors An. farauti and An. punctulatus were not abundant among these samples. Morphological identification of anophelines in this sample was often inaccurate, highlighting the importance of using molecular analysis in conjunction with morphological investigations to update keys and training tools.

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.

Significant geographical differences in prevalence of mutations associated with Plasmodium falciparum and Plasmodium vivax drug resistance in two regions from Papua New Guinea

Background

Drug resistance remains a major obstacle to malaria treatment and control. It can arise and spread rapidly, and vary substantially even at sub-national level. National malaria programmes require cost-effective and timely ways of characterizing drug-resistance at multiple sites within their countries.

Methods

An improved multiplexed post-PCR ligase detection reaction—fluorescent microsphere assay (LDR-FMA) was used to simultaneously determine the presence of mutations in chloroquine resistance transporter (crt), multidrug resistance 1 (mdr1), dihydrofolate reductase (dhfr) and dihydropteroate synthase (dhps) genes in Plasmodium falciparum (n = 727) and Plasmodium vivax (n = 574) isolates collected in 2006 from cross-sectional community population surveys in two geographically distinct regions (Madang and East Sepik) of Papua New Guinea (PNG) where strong regional differences in in vivo aminoquinoline and antifolate therapeutic efficacy had previously been observed. Data were compared to those of a follow-up survey conducted in 2010.

Results

Despite some very low parasite densities, the assay successfully amplified all P. falciparum and P. vivax loci in 77 and 69 % of samples, respectively. In 2006, prevalences of pfdhfr (59R-108 N) double mutation/wild type pfdhps haplotype, pfcrt SVMNT haplotype (72S-76T double mutation), and 86Y pfmdr1 mutation all exceeded 90 %. For P. vivax, 65 % carried at least two pvdhfr mutations, 97 % the 647P pvdhps mutation and 54 % the 976F pvmdr1 mutation. Prevalence of mutant haplotypes was higher in Madang than East Sepik for pfcrt SVMNT (97.4 vs 83.3 %, p = 0.001), pfdhfr (59R-108 N) (100 vs 90.6 %, p = 0.001), pvdhfr haplotypes (75.8 vs 47.6 %, p = 0.001) and pvmdr1 976F (71.2 vs 26.2 %, p < 0.001). Data from a subsequent Madang survey in 2010 showed that the prevalence of pfdhps mutations increased significantly from <5 % to >30 % (p < 0.001) as did the prevalence of pvdhfr mutant haplotypes (from 75.8 to 97.4 %, p = 0.012).

Conclusions

This LDR-FMA multiplex platform shows feasibility for low-cost, high-throughput, rapid characterization of a broad range of drug-resistance markers in low parasitaemia infections. Significant geographical differences in mutation prevalence correlate with previous genotyping surveys and in vivo trials and may reflect variable drug pressure and differences in health-care access in these two PNG populations.

Electronic supplementary material

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

Strategies for understanding and reducing the Plasmodium vivax and Plasmodium ovale hypnozoite reservoir in Papua New Guinean children: a randomised placebo-controlled trial and mathematical model

BACKGROUND

The undetectable hypnozoite reservoir for relapsing Plasmodium vivax and P. ovale malarias presents a major challenge for malaria control and elimination in endemic countries. This study aims to directly determine the contribution of relapses to the burden of P. vivax and P. ovale infection, illness, and transmission in Papua New Guinean children.

METHODS AND FINDINGS

From 17 August 2009 to 20 May 2010, 524 children aged 5-10 y from East Sepik Province in Papua New Guinea (PNG) participated in a randomised double-blind placebo-controlled trial of blood- plus liver-stage drugs (chloroquine [CQ], 3 d; artemether-lumefantrine [AL], 3 d; and primaquine [PQ], 20 d, 10 mg/kg total dose) (261 children) or blood-stage drugs only (CQ, 3 d; AL, 3 d; and placebo [PL], 20 d) (263 children). Participants, study staff, and investigators were blinded to the treatment allocation. Twenty children were excluded during the treatment phase (PQ arm: 14, PL arm: 6), and 504 were followed actively for 9 mo. During the follow-up time, 18 children (PQ arm: 7, PL arm: 11) were lost to follow-up. Main primary and secondary outcome measures were time to first P. vivax infection (by qPCR), time to first clinical episode, force of infection, gametocyte positivity, and time to first P. ovale infection (by PCR). A basic stochastic transmission model was developed to estimate the potential effect of mass drug administration (MDA) for the prevention of recurrent P. vivax infections. Targeting hypnozoites through PQ treatment reduced the risk of having at least one qPCR-detectable P. vivax or P. ovale infection during 8 mo of follow-up (P. vivax: PQ arm 0.63/y versus PL arm 2.62/y, HR = 0.18 [95% CI 0.14, 0.25], p < 0.001; P. ovale: 0.06 versus 0.14, HR = 0.31 [95% CI 0.13, 0.77], p = 0.011) and the risk of having at least one clinical P. vivax episode (HR = 0.25 [95% CI 0.11, 0.61], p = 0.002). PQ also reduced the molecular force of P. vivax blood-stage infection in the first 3 mo of follow-up (PQ arm 1.90/y versus PL arm 7.75/y, incidence rate ratio [IRR] = 0.21 [95% CI 0.15, 0.28], p < 0.001). Children who received PQ were less likely to carry P. vivax gametocytes (IRR = 0.27 [95% CI 0.19, 0.38], p < 0.001). PQ had a comparable effect irrespective of the presence of P. vivax blood-stage infection at the time of treatment (p = 0.14). Modelling revealed that mass screening and treatment with highly sensitive quantitative real-time PCR, or MDA with blood-stage treatment alone, would have only a transient effect on P. vivax transmission levels, while MDA that includes liver-stage treatment is predicted to be a highly effective strategy for P. vivax elimination. The inclusion of a directly observed 20-d treatment regime maximises the efficiency of hypnozoite clearance but limits the generalisability of results to real-world MDA programmes.

CONCLUSIONS

These results suggest that relapses cause approximately four of every five P. vivax infections and at least three of every five P. ovale infections in PNG children and are important in sustaining transmission. MDA campaigns combining blood- and liver-stage treatment are predicted to be a highly efficacious intervention for reducing P. vivax and P. ovale transmission.

TRIAL REGISTRATION

ClinicalTrials.gov NCT02143934.

Blood-Stage Parasitaemia and Age Determine Plasmodium falciparum and P. vivax Gametocytaemia in Papua New Guinea

A better understanding of human-to-mosquito transmission is crucial to control malaria. In order to assess factors associated with gametocyte carriage, 2083 samples were collected in a cross-sectional survey in Papua New Guinea. Plasmodium species were detected by light microscopy and qPCR and gametocytes by detection of pfs25 and pvs25 mRNA transcripts by reverse-transcriptase PCR (qRT-PCR). The parasite prevalence by PCR was 18.5% for Plasmodium falciparum and 13.0% for P. vivax. 52.5% of all infections were submicroscopic. Gametocytes were detected in 60% of P. falciparum-positive and 51% of P. vivax-positive samples. Each 10-fold increase in parasite density led to a 1.8-fold and 3.3-fold increase in the odds of carrying P. falciparum and P. vivax gametocytes. Thus the proportion of gametocyte positive and gametocyte densities was highest in young children carrying high asexual parasite densities and in symptomatic individuals. Dilution series of gametocytes allowed absolute quantification of gametocyte densities by qRT-PCR and showed that pvs25 expression is 10-20 fold lower than pfs25 expression. Between 2006 and 2010 parasite prevalence in the study site has decreased by half. 90% of the remaining infections were asymptomatic and likely constitute an important reservoir of transmission. However, mean gametocyte densities were low (approx. 1-2 gametocyte/μL) and it remains to be determined to what extent low-density gametocyte positive individuals are infective to mosquitos.

Uncovering the transmission dynamics of Plasmodium vivax using population genetics

Population genetic analysis of malaria parasites has the power to reveal key insights into malaria epidemiology and transmission dynamics with the potential to deliver tools to support control and elimination efforts. Analyses of parasite genetic diversity have suggested that Plasmodium vivax populations are more genetically diverse and less structured than those of Plasmodium falciparum indicating that P. vivax may be a more ancient parasite of humans and/or less susceptible to population bottlenecks, as well as more efficient at disseminating its genes. These population genetic insights into P. vivax transmission dynamics provide an explanation for its relative resilience to control efforts. Here, we describe current knowledge on P. vivax population genetic structure, its relevance to understanding transmission patterns and relapse and how this information can inform malaria control and elimination programmes.

Plasmodium vivax populations are more genetically diverse and less structured than sympatric Plasmodium falciparum populations

Introduction

The human malaria parasite, Plasmodium vivax, is proving more difficult to control and eliminate than Plasmodium falciparum in areas of co-transmission. Comparisons of the genetic structure of sympatric parasite populations may provide insight into the mechanisms underlying the resilience of P. vivax and can help guide malaria control programs.

Methodology/Principle findings

P. vivax isolates representing the parasite populations of four areas on the north coast of Papua New Guinea (PNG) were genotyped using microsatellite markers and compared with previously published microsatellite data from sympatric P. falciparum isolates. The genetic diversity of P. vivax (H_e = 0.83–0.85) was higher than that of P. falciparum (H_e = 0.64–0.77) in all four populations. Moderate levels of genetic differentiation were found between P. falciparum populations, even over relatively short distances (less than 50 km), with 21–28% private alleles and clear geospatial genetic clustering. Conversely, very low population differentiation was found between P. vivax catchments, with less than 5% private alleles and no genetic clustering observed. In addition, the effective population size of P. vivax (30353; 13043–69142) was larger than that of P. falciparum (18871; 8109–42986).

Conclusions/Significance

Despite comparably high prevalence, P. vivax had higher diversity and a panmictic population structure compared to sympatric P. falciparum populations, which were fragmented into subpopulations. The results suggest that in comparison to P. falciparum, P. vivax has had a long-term large effective population size, consistent with more intense and stable transmission, and limited impact of past control and elimination efforts. This underlines suggestions that more intensive and sustained interventions will be needed to control and eventually eliminate P. vivax. This research clearly demonstrates how population genetic analyses can reveal deeper insight into transmission patterns than traditional surveillance methods.

The neglected human malaria parasite Plasmodium vivax is responsible for a large proportion of the global malaria burden. Efforts to control malaria have revealed that P. vivax is more resilient than the other major human malaria parasite, Plasmodium falciparum. This study utilised population genetics to compare patterns of P. vivax and P. falciparum transmission in Papua New Guinea, a region where infection rates of the two species are similar. The results demonstrated that P. vivax populations are more genetically diverse than those of P. falciparum suggestive of a parasite population that is more resilient to environmental challenges, undergoing higher levels of interbreeding locally and between distant parasite populations. Unique characteristics of P. vivax such as relapse, which allows different strains from past infections to produce subsequent infections, may provide more opportunities for the exchange and dissemination of genetic material. The contrasting patterns observed for the two species may be the result of a differential impact of past elimination attempts and indicate that more rigorous interventions will be needed in efforts to control and eventually eliminate P. vivax.

Predicting antidisease immunity using proteome arrays and sera from children naturally exposed to malaria

Malaria remains one of the most prevalent and lethal human infectious diseases worldwide. A comprehensive characterization of antibody responses to blood stage malaria is essential to support the development of future vaccines, sero-diagnostic tests, and sero-surveillance methods. We constructed a proteome array containing 4441 recombinant proteins expressed by the blood stages of the two most common human malaria parasites, P. falciparum (Pf) and P. vivax (Pv), and used this array to screen sera of Papua New Guinea children infected with Pf, Pv, or both (Pf/Pv) that were either symptomatic (febrile), or asymptomatic but had parasitemia detectable via microscopy or PCR. We hypothesized that asymptomatic children would develop antigen-specific antibody profiles associated with antidisease immunity, as compared with symptomatic children. The sera from these children recognized hundreds of the arrayed recombinant Pf and Pv proteins. In general, responses in asymptomatic children were highest in those with high parasitemia, suggesting that antibody levels are associated with parasite burden. In contrast, symptomatic children carried fewer antibodies than asymptomatic children with infections detectable by microscopy, particularly in Pv and Pf/Pv groups, suggesting that antibody production may be impaired during symptomatic infections. We used machine-learning algorithms to investigate the relationship between antibody responses and symptoms, and we identified antibody responses to sets of Plasmodium proteins that could predict clinical status of the donors. Several of these antibody responses were identified by multiple comparisons, including those against members of the serine enriched repeat antigen family and merozoite protein 4. Interestingly, both P. falciparum serine enriched repeat antigen-5 and merozoite protein 4 have been previously investigated for use in vaccines. This machine learning approach, never previously applied to proteome arrays, can be used to generate a list of potential seroprotective and/or diagnostic antigens candidates that can be further evaluated in longitudinal studies.

Distinct patterns of diversity, population structure and evolution in the AMA1 genes of sympatric Plasmodium falciparum and Plasmodium vivax populations of Papua New Guinea from an area of similarly high transmission

Background

As Plasmodium falciparum and Plasmodium vivax co-exist in most malaria-endemic regions outside sub-Saharan Africa, malaria control strategies in these areas must target both species in order to succeed. Population genetic analyses can predict the effectiveness of interventions including vaccines, by providing insight into patterns of diversity and evolution. The aim of this study was to investigate the population genetics of leading malaria vaccine candidate AMA1 in sympatric P. falciparum and P. vivax populations of Papua New Guinea (PNG), an area of similarly high prevalence (Pf = 22.3 to 38.8%, Pv = 15.3 to 31.8%).

Methods

A total of 72 Pfama1 and 102 Pvama1 sequences were collected from two distinct areas, Madang and Wosera, on the highly endemic PNG north coast.

Results

Despite a greater number of polymorphic sites in the AMA1 genes of P. falciparum (Madang = 52; Wosera = 56) compared to P. vivax (Madang = 36, Wosera = 34), the number of AMA1 haplotypes, haplotype diversity (Hd) and recombination (R) was far lower for P. falciparum (Madang = 12, Wosera = 20; Hd ≤0.92, R ≤45.8) than for P. vivax (Madang = 50, Wosera = 38; Hd = 0.99, R = ≤70.9). Balancing selection was detected only within domain I of AMA1 for P. vivax, and in both domains I and III for P. falciparum.

Conclusions

Higher diversity in the genes encoding P. vivax AMA1 than in P. falciparum AMA1 in this highly endemic area has important implications for development of AMA1-based vaccines in PNG and beyond. These results also suggest a smaller effective population size of P. falciparum compared to P. vivax, a finding that warrants further investigation. Differing patterns of selection on the AMA1 genes indicate that critical antigenic sites may differ between the species, highlighting the need for independent investigations of these two leading vaccine candidates.

Global Population Structure of the Genes Encoding the Malaria Vaccine Candidate, Plasmodium vivax Apical Membrane Antigen 1 (PvAMA1)

Background

The Plasmodium vivax Apical Membrane Antigen 1 (PvAMA1) is a promising malaria vaccine candidate, however it remains unclear which regions are naturally targeted by host immunity and whether its high genetic diversity will preclude coverage by a monovalent vaccine. To assess its feasibility as a vaccine candidate, we investigated the global population structure of PvAMA1.

Methodology and Principal Findings

New sequences from Papua New Guinea (PNG, n = 102) were analysed together with published sequences from Thailand (n = 158), India (n = 8), Sri Lanka (n = 23), Venezuela (n = 74) and a collection of isolates from disparate geographic locations (n = 8). A total of 92 single nucleotide polymorphisms (SNPs) were identified including 22 synonymous SNPs and 70 non-synonymous (NS) SNPs. Polymorphisms and signatures of balancing (positive Tajima's D and low F_ST values) selection were predominantly clustered in domain I, suggesting it is a dominant target of protective immune responses. To estimate global antigenic diversity, haplotypes comprised of (i) non-singleton (n = 40) and (ii) common (≥10% minor allele frequency, n = 23) polymorphic amino acid sites were then analysed revealing a total of 219 and 210 distinct haplotypes, respectively. Although highly diverse, the 210 haplotypes comprised of only common polymorphisms were grouped into eleven clusters, however substantial geographic differentiation was observed, and this may have implications for the efficacy of PvAMA1 vaccines in different malaria-endemic areas. The PNG haplotypes form a distinct group of clusters not found in any other geographic region. Vaccine haplotypes were rare and geographically restricted, suggesting potentially poor efficacy of candidate PvAMA1 vaccines.

Conclusions

It may be possible to cover the existing global PvAMA1 diversity by selection of diverse alleles based on these analyses however it will be important to first define the relationships between the genetic and antigenic diversity of this molecule.

Traditionally misclassified as benign and neglected as a research priority, it is now understood that P. vivax is an increasingly important cause of human malaria. This important human pathogen poses an enormous obstacle to malaria control and elimination efforts due its broad geographic distribution, ability to cause recurring episodes of malaria after long periods of inactivity and extreme biodiversity. Vaccines are an essential component of global malaria control and elimination campaigns but the diversity of malaria antigens is thought to be a major cause of vaccine failure. Furthermore, at present the majority of current vaccine research is directed toward P. falciparum. The aims of this study were to investigate the global diversity of the P. vivax vaccine candidate, Apical Membrane Antigen 1 (PvAMA1), to determine the feasibility of designing a globally effective PvAMA1 vaccine and to determine which region of PvAMA1 is targeted by host immune responses, in order to identify the most promising vaccine candidates. We report that PvAMA1 diversity is extremely high, and that PvAMA1 domain I is a dominant target of host immune responses. These analyses of PvAMA1 diversity from several geographic regions provide a framework to guide development of a broadly efficacious P. vivax vaccine.