PLoS Medicine
Home The changing epidemiology of Plasmodium vivax: Insights from conventional and novel surveillance tools
The changing epidemiology of <i>Plasmodium vivax</i>: Insights from conventional and novel surveillance tools
The changing epidemiology of Plasmodium vivax: Insights from conventional and novel surveillance tools

The authors have declared that no competing interests exist.

Article Type: Review Article Article History
Abstract

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

Auburn,Cheng,Marfurt,and Price: The changing epidemiology of Plasmodium vivax: Insights from conventional and novel surveillance tools

Summary points

    Renewed efforts to eliminate malaria have had greater impact on Plasmodium falciparum than Plasmodium vivax, a reflection of the fundamental differences in the biology of the parasite, its transmission dynamics, and ability to form dormant liver stages.

    The main burden of P. vivax malaria is in young children residing in remote communities with poor access to healthcare services.

    The decline in P. vivax malaria has led to an increasing proportion of the parasite reservoir occurring in asymptomatic and low-density P. vivax infections and heterogeneous patterns of parasite transmission.

    Genetic tools to study the spatial and temporal patterns of P. vivax transmission in different endemicities and ultrasensitive PCR (uPCR)-based techniques are expanding our knowledge of the magnitude and biology of low-density P. vivax infections.

    Serology offers alternative ways of detecting recent P. vivax infections and monitoring of the impact of public health interventions at very low endemicity.

Introduction

Plasmodium vivax remains an important public health burden affecting the poorest and most vulnerable communities of more than 49 endemic countries [1]. Over the past 2 decades, investment in malaria control and elimination programmes has helped to shrink the global malaria map considerably; however, in areas where P. vivax and Plasmodium falciparum are co-endemic, P. vivax is becoming the predominant cause of malaria [2]. The rising proportion of P. vivax highlights the greater transmission potential of this species attributable to several biological properties that differ significantly from P. falciparum (Table 1). Key features of P. vivax confounding its control and elimination include the parasite’s ability to form dormant liver stages (hypnozoites) that can cause relapse weeks to months after an initial infection, the ability to circulate at low peripheral parasite densities while remaining transmissible to the mosquito vector, and the early production of sexual stages (gametocytes) that allow the parasite to transmit prior to clinical presentation and gametocytocidal treatment. P. vivax has also evolved to survive in a greater variety of Anopheles vectors. Collectively, these properties sustain parasite survival and transmission in more extreme climates and confer greater resilience than P. falciparum against conventional parasite and vector control activities [3].

Table 1
Biological properties of P. vivax that enhance the capacity for transmission and spread.
Biological propertyChallengePotential solutions
Low-density infectionsClinical illness and transmission occur at low-level parasitaemia, which is hard to diagnoseHigh-sensitivity molecular diagnostics such as PCR and large volume uPCR
Liver-stage reservoirHepatic reservoir is difficult to diagnose and results in recurrent infections that sustain ongoing transmission and confound the assessment of TESSerological markers to identify individuals with recent exposure and high risk of relapse who can be offered radical cure. Genotyping of infections initially and at the time of recurrence to refine classifications (recrudescence, reinfection, or relapse) in clinical and epidemiological studies
Early gametocytaemia prior to clinical presentationEnhances transmission prior to antimalarial treatmentEarly diagnosis and treatment of clinical illness, prevention of relapses, and transmission interventions such as indoor residual spraying with insecticides and LLINs
Capacity to develop in mosquitoes at low ambient temperaturesWide range of ecological receptivity, enhancing transmissionTransmission interventions such as indoor residual spraying with insecticides and LLINs

LLIN, long-lasting insecticide treated net; TES, therapeutic efficacy studies; uPCR, ultrasensitive PCR.

Despite the relative resilience of P. vivax, its incidence has declined in most endemic regions, falling from an estimated global burden of 24.5 million clinical cases in 2000 to 14.3 million in 2017 [2]. This decline has been associated with changing patterns in the epidemiology of P. vivax as endemicity shifts from high and stable to intermediate and pre- and post-elimination levels (Fig 1). Examples of the spectrum of national endemicity across 49 P. vivax–endemic countries are presented in Table 2. However, the mean national incidence measures presented in the table mask marked heterogeneity in endemicity at the micro-epidemiological level. In this review, we discuss the complex dynamics of parasite transmission and spread across this spectrum of endemicity, the epidemiological challenges that they present, and the tools that are being used to monitor them.

Trends in parasite transmission and diversity, with the spectrum of clinical illness and demographics varying with endemicity.
Fig 1

Trends in parasite transmission and diversity, with the spectrum of clinical illness and demographics varying with endemicity.

The schematic was inspired by Fig 1 in Wesolowski et al. [36]. As transmission declines, parasite populations become increasingly fragmented into distinct foci and enhanced inbreeding leads to clonal and highly related infections. An increasing proportion of infections are acquired from importations, with high risk of outbreaks. Cases become harder to detect and treat owing to increasingly asymptomatic outcome and SM level. A shift is also observed towards high-risk occupational groups and hard-to-reach populations. SM, submicroscopic.

Table 2
National incidence of P. vivax malaria in 49 endemic countries.
Incidence of malariaHighMediumLowVery low
API>10 per 1,0001–10 per 1,0000.1–1 per 1,000<0.1 per 1,000
Afghanistan

Guyana

Pakistan

Papua New Guinea

Solomon Islands

South Korea

Sudan

Vanuatu

Venezuela
Belize

Djibouti

Ecuador

Guatemala

Honduras

Nepal

North Korea

Panama

Somalia

Suriname

Thailand

Yemen
Bolivia

Brazil

Cambodia

Colombia

Costa Rica

Eritrea

Ethiopia

French Guiana

India

Indonesia

Laos

Madagascar

Myanmar

Nicaragua

Peru

Philippines
Bangladesh

Bhutan

China*

El Salvador

Iran

Kenya

Malaysia*

Mexico

Oman*

Saudi Arabia

Timor-Leste*

Vietnam

Incidences generated from Battle et al. [2], from the estimated number of clinical cases of P. vivax per annum for the population at risk within each country. These classifications provide only a rough approximation since there is marked heterogeneity in endemicity at subnational level. Quantifying endemicity and its heterogeneity at a local level is critical to guide malaria control activities.

* No indigenous cases since 2018.

API, annual parasite incidence.

Characterising P. vivax incidence: The complexity of recurrence

Whereas the incidence of P. falciparum is driven by recurrent inoculation from infected mosquitoes (reinfections), the incidence of P. vivax is attributable to reactivation of hypnozoites (relapses). In some areas, antimalarial drug resistance, predominantly to chloroquine, results in schizontocidal treatment failure that contributes to early parasite recurrence (recrudescence) [4]. Radical cure of P. vivax refers to killing both the blood and liver stages of the parasites to ensure not only clinical recovery from acute infection, but also prevention of subsequent relapses. The frequency and timing of these relapses vary considerably with geographical region, host immunity, and the number of sporozoites inoculated by the mosquito the absolute risk over 12 months ranging from less than 10% in temperate areas to more than 80% in tropical climates [5]. Recurrent parasitaemia in individuals continuing to reside in malaria endemic settings can also be due to reinfection; however, longitudinal cohorts of patients treated with and without primaquine radical cure suggest that these account for less than 30% of all recurrences [6].

High and stable endemicity

In areas of high P. vivax transmission and relapse periodicity, nonimmune individuals can have recurrent episodes of malaria every few weeks, a similar force of infection as that seen in areas hyperendemic for P. falciparum. Recurrent bouts of malaria result in a cumulative risk of anaemia, childhood malnutrition, and a rise in associated morbidity and mortality [7]. However, recurrent P. vivax infections in early life result in a faster acquisition of immunity than that following P. falciparum, non-sterilising immunity suppressing clinical disease, and a lower risk of illness in older individuals [8]. In these areas, the primary goal of malaria control programmes is to reduce clinical cases and ongoing transmission of the parasite, largely through vector control measures, widespread access to early diagnosis of patients with clinical illness, and provision of effective antimalarial treatment. Effective surveillance tools to monitor treatment efficacy and progress in transmission reduction are critical.

Monitoring the impact of public health interventions

Conventionally, the surveillance of P. vivax malaria has relied upon reporting of the prevalence and incidence of P. vivax infection, identified by light microscopy (LM) or rapid diagnostic tests (RDTs). In highly endemic settings, case detection is generally passive, focused on diagnosing individuals with clinical illness presenting to health facilities. This approach is constrained by a lack of sensitivity in detecting low-level parasitaemia and failure to detect individuals with asymptomatic infections. Although molecular laboratory methods, targeting parasite DNA, offer greater sensitivity in quantifying the extent of the parasite reservoir, these methods are not feasible or cost-effective in highly endemic settings. An alternative approach to gauge transmission intensity is quantifying parasite genetic diversity from population-representative samples.

Studies of P. falciparum populations reveal increasing genetic diversity and prevalence of polyclonal infections with increasing transmission [9]. This trend can be explained by the higher frequency of individuals with multiple infectious mosquito bites in high transmission settings, leading to superinfection (multiple inoculations from different mosquitoes) or coinfection events (a single inoculation with multiple clones) and outcrossing between different clones in the sexual stage of development in the mosquito. However, the correlation between population diversity, polyclonality, and endemicity is less apparent in P. vivax, with several studies revealing high polyclonality and population diversity even at low endemicity [10,11]. As an example, Fig 2 illustrates patterns in polyclonality and diversity across a range of endemic settings in a selection of studies that all used the same set of microsatellite markers, revealing several low endemic sites with high diversity. Despite low endemicity, high diversity and polyclonality have been documented in Iran (65% polyclonal) and Sri Lanka (69% polyclonal), potentially reflecting sustained infection diversity from imported cases [12,13]. Genome-wide studies from South American populations have revealed similar trends of extensive genetic diversity in low transmission settings, potentially reflecting a complex composite of historical events and more recent epidemiological patterns [14,15]. Conversely, longitudinal studies reveal a positive correlation between P. vivax diversity and endemicity. In Sabah, Malaysia, where P. vivax was eliminated in 2020, there has been a marked decline in both population diversity and polyclonality (50% to 8% from 2010 to 2015) [16]. Whereas in the highly endemic setting of Papua, Indonesia, where the incidence of P. vivax fell by almost 30% over a 9-year period, polyclonality fell from 67% to 35%, despite no change in parasite diversity [17,18]. A more consistent finding across studies is increasing genetic structure with declining incidence, but this is apparent only after a substantial reduction in transmission [16,19,20]. Collectively, these studies highlight the potential utility of genetic metrices to identify moderate to large changes in transmission reduction that can complement conventional malariometric surveillance. However, a better understanding of how historical and more recent migratory and epidemiological factors may have shaped local population diversity will be needed. In addition, further research is needed to assess how genetic indices can be made readily available and provide timely, actionable, and cost-effective information to malaria control programmes [21].

Trends in P. vivax polyclonality and population diversity across endemic settings.
Fig 2

Trends in P. vivax polyclonality and population diversity across endemic settings.

The data were derived from published, open-access online vivaxGEN-MS platform [113]. Sites represent first- or second-level administrative divisions and include areas with a minimum of 15 clinical isolates. Measures were restricted to 3 microsatellite markers: Pvmsp1f3, PvMS16, and Pv3.27 [114]. Polyclonality reflects the percentage of P. vivax infections with 1 or more markers exhibiting more than 1 allele. The population diversity was approximated using the mean expected heterozygosity. Approximations of the API per 1,000 population were derived from the original publications [16,18,41,91,115119]. The results demonstrate a general trend of increasing polyclonality with increasing API, with several deviations potentially reflecting factors such as imported cases. In contrast, population-level diversity remains high across endemic settings, potentially reflecting sustained diversity from the liver-stage reservoir. API, annual parasite incidence.

Monitoring treatment efficacy

Chloroquine remains the most widely used first-line treatment of P. vivax, but increasing reports of drug failure highlight the importance of surveillance of its clinical efficacy [4]. Several candidate markers of chloroquine resistance in P. vivax have been derived from orthologues of P. falciparum resistance loci, population genomic signatures of selection, and, most recently, a genetic cross; however, none of these markers have been validated against clinical outcomes [16,2226]. A crucial advance in therapeutic efficacy surveys (TES) of P. falciparum was the ability to differentiate between reinfections (which tend to be genetically different from the initial infection) and recrudescence (which are genetically the same) using genotyping [27]. The challenge is far greater in P. vivax since recurrent infections can also arise from dormant hypnozoites (relapses), which may be genetically similar or different from the initial infection [28,29]. Recent studies have demonstrated the utility of genetics to disentangle reinfections from relapses by identifying paired isolates before and after treatment that, while being genetically different, demonstrate familial relationships [3032]. The rationale is that genetically related infections, such as siblings, are more likely to have derived from the same mosquito inoculation than from different inoculations (reinfections) and are therefore more likely to be relapse events. A mathematical model incorporating genetic and time-to-event data in a TES of primaquine anti-relapse efficacy in Thailand revealed only 2.9% treatment failures compared with an unadjusted clinical estimate of 12% [33]. These genetic tools will greatly enhance the interpretation of TES outcomes for both schizontocidal and hypnozonticidal agents.

Intermediate endemicity

As transmission intensity declines, the parasite population contracts, transmission becomes increasingly heterogeneous even at small spatial scales, and focalised networks of infection are apparent [34]. In these areas, low-level parasitaemia and asymptomatic infections become increasingly important reservoirs of infection that can sustain ongoing transmission, despite falling numbers of clinical cases. Capturing the burden of these low-level and asymptomatic infections is a major priority for control programmes. A better understanding of the spatial heterogeneity in parasite transmission will inform targeted public health interventions, by identifying residual hotspots of transmission and defining the boundaries and connectivity of these foci.

Identifying foci of infection and their connectivity

Currently, foci of infection are identified using geospatial surveillance and case investigations [35]. Information on the genetic relatedness between parasites provides a more direct link between infections, enabling greater accuracy in characterising spatial patterns of transmission [36]. Spatial studies of malaria that integrate genetic data with information on patient mobility, such as travel history and mobile phone data, can help characterise transmission networks. Such approaches are relatively advanced for P. falciparum [37,38], but far more complex for P. vivax due to confounding from multiple and delayed relapses. Recently, studies using genetic data have shed some light on the connectivity between P. vivax populations at a subnational level, revealing a spectrum of patterns [11]. In Indonesia and the Pacific region, the connectivity between P. vivax populations in different islands was surprisingly high despite the geographic barrier to gene flow [3941]. Comparisons between co-endemic P. vivax and P. falciparum populations in these and other studies reveal consistently greater connectivity among the P. vivax populations, suggesting that the hypnozoite reservoir enhances the spread of P. vivax infections through human travel [3943]. In parts of the Peruvian Amazon and along the Pacific Coast of Colombia, connected endemic populations referred to as malaria “corridors” appear to enable P. vivax infections to spread over long distances, some as far as 1,500 km apart [44,45]. Whereas in other areas of Latin America, P. vivax populations have much lower connectivity [42,46,47]; restricted human movement and local adaptative pressures leading to the frequent replacement of parasite strains increase the geographic differentiation between these populations. Given the complex and often unpredictable patterns of parasite connectivity, comprehensive genetic surveillance frameworks will be needed to understand how P. vivax parasites are spreading within and between communities [21]. The methods used to characterise parasite connectivity are also an important consideration for malaria parasites, which have an obligate stage of sexual reproduction in the mosquito midgut. Phylogenetic methods are effective for quantifying genetic relatedness between rapidly mutating organisms that do not recombine but are less suitable for malaria parasites [48]. Furthermore the fixation index (FST), a classical approach used to quantify connectedness between populations, does not account for recombination and is therefore constrained in its ability to infer fine-resolution connectedness [49]. Recent studies have employed measures of the probability of interindividual identity by descent (IBD) to characterise the relatedness between malaria parasites, revealing more fine-scale insights into malaria parasite transmission networks [50,51]. However, all these methods are challenged by polyclonal infections, an obstacle that will require new statistical innovations.

The hidden burden of submicroscopic and asymptomatic infections

The peripheral parasitaemia of P. vivax, at which an individual develops symptoms, is estimated to be between 40 and 2,000 parasites per μl [5254]. While the pyrogenic threshold of P. vivax has been proposed to be significantly lower for P. vivax compared to P. falciparum, these estimates do not account for a significant noncirculating biomass. While sequestration of P. falciparum is well recognised, emerging data suggest a hidden reservoir of P. vivax [55,56].

Expert microscopic examination of a blood film can detect peripheral parasitaemias as low as 25 parasites μl−1, a level at which there are approximately 108 parasites in the peripheral circulation of an adult. However, in routine practice field, microscopy can usually only detect >100 parasites μl−1 reliably. Over the past 2 decades, molecular-based diagnostic assays have enabled the threshold of detection to be lowered to 0.02 μl−1, unveiling a large hidden burden of submicroscopic (SM) and asymptomatic infections [57]. Since individuals with asymptomatic, low-level, or SM P. vivax parasitaemia are unlikely to seek medical attention and their infections are likely to be missed by LM or RDTs, they generally remain undetected and untreated.

The scale of the challenge for malaria elimination has been highlighted by the application of molecular approaches to cross-sectional surveys and clinical cohorts which reveal a substantial “hidden” reservoir of infection across the spectrum of transmission settings. Key questions pertinent to malaria elimination strategies include the magnitude of this hidden burden and the degree to which individuals with asymptomatic and SM P. vivax infections are contributing to ongoing transmission.

In areas of intermediate endemicity, individuals repeatedly exposed to P. vivax malaria, either by reinfection or relapse, develop a degree of premunition that can suppress clinical symptoms despite rising levels of parasitaemia, hence some asymptomatic P. vivax carriage can be detected by either LM or PCR. Cross-sectional surveys in the Asia-Pacific region have shown that up to 97% of individuals with P. vivax parasitaemia detected by LM and 100% of SM infections are asymptomatic [58,59]. In a large cross-sectional study in an area of low and seasonal transmission, ultrasensitive PCR (uPCR, with a level of detection of 0.022 parasites μl−1) detected asymptomatic infections in almost a third of individuals of whom 70% had parasite densities below 100 μl−1 [57]. Importantly, asymptomatic P. vivax infections have been shown to include gametocyte carriage and be equally infectious to mosquitoes as symptomatic individuals, suggesting significant potential to contribute to parasite transmission [60,61].

The relative burden of SM P. vivax infections has been estimated by comparing P. vivax prevalence detected by LM and PCR. Overall, SM parasitaemia accounted for almost 70% of infection, although there was marked heterogeneity in different transmission settings [59,62]. The relative size of this hidden reservoir increases with increasing sensitivity of parasite detection method used and is greatest in areas with lower transmission intensity [59,63].

The contribution of SM and asymptomatic P. vivax infections to transmission

The contribution of SM P. vivax infections to transmission is determined by several factors including the duration of infection, variation in parasite densities, gametocyte carriage, infectivity to mosquitoes, and propensity to relapse. P. vivax gametocytes have been shown to transmit more efficiently to Anopheles than those of P. falciparum and to be more transmissible at lower parasite densities [64]. However, transmission dynamics may vary significantly between different vector species and thus these comparisons between parasite species cannot be taken as absolute across the spectrum of malaria endemic environments.

In naturally acquired P. vivax infections, low-level peripheral parasitaemia fluctuates. Although the duration of infection and dynamics of parasitaemia are likely to determine transmissibility, these are difficult to quantify due to the challenges in conducting longitudinal cohorts with frequent molecular sampling and difficulties in distinguishing persisting infections from relapses and new infections.

In an area of low and seasonal malaria in western Cambodia, 24 asymptomatic individuals with P. vivax monoinfection were followed for 11 months and assessed monthly using uPCR [65]. All of those assessed tested positive for P. vivax at 3 or more monthly visits, and almost 20% tested positive every month. Interestingly, of the individuals who had parasite densities <100 μl−1 at the start of the cohort, two-thirds never had densities >100 μl−1 over the subsequent 11 months. Similar trends were observed in a clinical cohort conducted in Lao PDR [66]. These findings suggest that a significant proportion of individuals with SM P. vivax control their infection, with parasite densities rarely fluctuating above microscopically detectable levels during recurrent or persisting infection; if confirmed, there may be limited transmission of these SM infections to mosquitoes. Importantly, strategies to target low-level P. vivax infections by identifying individuals with transient rises in microscopic parasitaemia will not detect most P. vivax infections. Further studies are required to ascertain the dynamics of parasite densities, clinical outcomes, and transmissibility of these P. vivax infections over time.

P. vivax gametocyte carriage can be detected by LM and by reverse transcription quantitative PCR (RT-qPCR) of the Pvs25 transcript [67]. In acute symptomatic illness, P. falciparum gametocytaemia tends to occur 7 to 14 days after the appearance of asexual stages, whereas individuals infected with P. vivax who develop clinical illness often present with both asexual and sexual parasite stages in the peripheral circulation. In both symptomatic and asymptomatic infection, there is a positive correlation between asexual parasitaemia and gametocytaemia, although this correlation varies significantly [68]. Studies in rural Amazonia and Ethiopia where more than half of P. vivax–infected individuals were asymptomatic revealed that 80% to 90% of asymptomatic and SM-infected individuals had detectable gametocytaemia [60,68]. In contrast, gametocytes were detected in only 20% of asymptomatic and SM P. vivax infections in Peru [69]. In the Solomon Islands and Papua New Guinea, Pvs25 transcripts were detected in under 20% of individuals with SM infections compared to 41% to 80% of those with microscopic infection [70,71].

Blood from asymptomatic P. vivax individuals, mostly detected only by PCR, has been shown to infect 1.2% Anopheles darlingi, compared to 22% of blood from symptomatic individuals [72]. Similarly, membrane feeding experiments demonstrate that 3.3% asymptomatic/PCR-detected P. vivax individuals were infectious compared to 32% of asymptomatic, microscopically positive individuals [68]. Based on these data, the authors estimate that in Ethiopia, SM P. vivax infections contribute about 14% of the infectious reservoir compared to 56% by the asymptomatic, microscopic detected individuals. However, these studies do not demonstrate the degree to which individuals with low-level infection contribute to transmission over a prolonged period of time. The contribution of low-level P. vivax infections to transmission is likely to be amplified by the number and frequency of relapses in infected individuals with hypnozoites. While individuals with SM P. vivax infection may control parasite densities at low levels, a Plasmodium cynomolgi model suggests that during relapses, there may be a higher proportion of gametocytes [73]. The evidence highlights that to achieve timely elimination of the parasite, greater efforts should be made to identify all individual harbouring parasites and offering them radical cure of all parasite stages.

It should be noted that even though both asexual and sexual parasites are detected by PCR from peripheral blood, this may not reflect the total P. vivax biomass. A significant number of trophozoites were seen sequestered in the spleen of one, and both asexual stages and gametocytes were observed in the bone marrow of another P. vivax–infected individual at similar densities as peripheral blood [74,75]. These sequestered parasites may help to sustain parasitemia in peripheral blood for transmission.

Low endemicity: Pre- and post-malaria elimination

As countries approach the end stages of malaria elimination, an increasing proportion of infections are imported, either across porous borders or from long-distance importations [76]. There is a shift in the populations at greatest risk of malaria infection, for instance, a relative rise of infections in adult males in occupations that increase their exposure to the local vectors or travel to endemic areas [7781]. Residual pockets of infection are also observed in hard-to-reach communities such as those in mountainous or densely forested regions, with poor road access [81,82]. New technologies, such as mobile mapping tools, have potential to facilitate tracking human mobility and associated malaria transmission networks in these hard-to-reach populations [83]. At this stage of malaria control, the major goals are to diagnose and treat individuals in the residual pockets of endemicity, identify imported cases, and monitor for evidence of local transmission.

Surveillance of imported cases

Traditionally, imported malaria has been identified and mapped using information from patients’ travel history. However, persistent blood stage infections and long-latency liver stages constrain the accuracy of this approach in P. vivax. Genetic studies of P. vivax have shown that infections from different countries can be differentiated, in some areas with just a handful of single nucleotide polymorphism (SNP) or microsatellite markers [8488]. A recent genomic study across 20 countries identified 28 SNPs that could distinguish imported from local infections in most areas and provides the first online tool with which investigators can map infections to their country of origin [88]. These genetic tools can complement traditional epidemiologic approaches to inform on imported P. vivax cases and facilitate surveillance and final certification of malaria elimination status.

While progress has been made in enhancing the surveillance of long-distance P. vivax importations, the resolution of local from imported cases is still constrained in porous border regions. In regions such as Yunnan Province, China, which shares hundreds of kilometres of border with highly endemic Myanmar, border malaria threatens efforts to prevent community transmission [89]. In large stretches of this border region, parasite gene flow is high, with most cases likely to be imported from Myanmar into China [90]. Similar challenges are faced in several other countries, including southern Bhutan, which shares a border with India, where genetic analyses highlight high homology between local and imported infections [91]. In porous border regions, cross-country collegiate strategies are needed to identify and contain shared reservoirs of P. vivax. Innovative, robust new tools, and multicountry and regional networks bringing countries together to share knowledge and tools will be critical to strengthen these efforts.

Detecting local transmission at very low case incidence: The role of seroepidemiology

Seroepidemiology is an alternative approach that can inform parasite surveillance beyond that offered by LM, RDT, and genetic-based methods [92]. Antibodies are biomarkers for exposure to infection, with the advantage that their response profiles last beyond a resolved acute infection when parasites are no longer detectable by LM, RDTs, or nucleic acid amplification–based tests (NAATs) [93]. To date, most studies have focused on P. falciparum, where seroepidemiological data have been correlated with malaria endemicity and been used to identify transmission hotspots, monitor the impact of control interventions, and evaluate elimination status [9497]. The advent of the new enzyme-linked immunosorbent assays (ELISAs) technology, which allows assessment of antibody responses to single, or multiple antigens at higher throughput, paralleled by the development and application of more sophisticated analytical methods and models, has expanded our seroepidemiological knowledge and offers opportunities to apply this approach to national and regional surveillance. The rate at which conversion to seropositivity occurs and the longevity of the antibody response to a given antigen vary for different targets and is largely dependent on transmission intensity and booster effects by consecutive reinfections [98]. Age-specific prevalence and seroconversion rates of antibodies against stage-specific antigenic targets such as PfAMA-1, PfMSP-119, and PfMSP-2 have potential to define seasonal and long-term transmission trends [99]. Depending on the stage specificity of antigens, estimates of disease exposure and transmission potential can be further refined: Antibody responses to preerythrocytic antigen targets such as circumsporozoite protein (CSP) and liver-stage antigens (LSAs) have been used to predict recent exposure and monitor temporal trends [100]. Furthermore, antibody responses against gametocyte-expressed targets such as Pfmdv1 and Pfs16 have been used as indicators of gametocyte carriage and hence, transmission intensity [101]. Antibody responses against the Anopheles gambiae salivary gland protein gSG6 or the synthetic peptide gSG6-P1 have been correlated with mosquito exposure [94,102]. While these serological markers have been successfully validated at multiple sites in Africa, Asia, and Southeast Asia, other studies suggest that anti-gSG6 targets of local anopheline vectors are more reliable for estimating transmission intensity [103106].

The use of protein microarrays and/or fluorescent microsphere-based technologies has advanced the field of seroepidemiology by allowing multiplexing a larger number of antigens at high throughput and thus, enabling screening and validating a multitude of candidate antigens. In a recent study of patients from Thailand, Brazil and the Solomon Islands were screened for more than 300 P. vivax antigens as potential markers of recent exposure; 8 antigens were identified that detected recent exposure with 80% sensitivity and specificity [107]. Since relapsing infections are key to sustaining P. vivax transmission, it is crucial to identify individuals at high risk of carrying hypnozoites, so that they can be offered radical cure. Currently, this approach focuses predominantly on patients presenting to clinics with vivax malaria but could be expanded to individuals with asymptomatic parasitaemia detected by active screening, detecting carriers of occult hypnozoites. Further work is underway to confirm these findings in different locations with differing transmission intensity and test the utility of simple RDTs in detecting individuals at risk of relapse.

Is the last parasite standing the most resistant?

While the first reports of chloroquine resistance in P. falciparum were in the late 1950s, chloroquine-resistant P. vivax was not reported until 1989 on the island of New Guinea [108,109]. The delayed emergence in P. vivax has been attributed in part to the early, presymptomatic gametocytogenesis in this species [110]. Early gametocytogenesis may also enhance the persistence in the population of drug-resistant P. vivax infections once established locally. In the malaria pre-elimination setting of Sabah, Malaysia, a clinical efficacy study of chloroquine revealed that over 60% of patients failed treatment with chloroquine within 28 days [111]. Interestingly, genomic evaluation of P. vivax infections collected in the same study demonstrated large clonal outbreaks and a high prevalence of drug resistance variants, including the multidrug resistance 1 gene (pvmdr1) Y976F mutation. The latter has been proposed as a marker of chloroquine resistance, although is at best a minor determinant [16,112]. The combination of unstable transmission and high prevalence of drug-resistant parasites in the late stages of malaria elimination presents a significant risk of resurgence, highlighting the need for diligent surveillance of clinical efficacy. Since clinical trials become logistically challenging as case numbers fall, there is a reliance on molecular surveillance, although for this to be applicable to P. vivax, better markers of antimalarial resistance will need to be identified.

Summary

As the map of P. vivax malaria shrinks, there are marked changes in the epidemiology of the parasite that impact upon malaria control activities. Symptomatic malaria is no longer a disease predominantly of children, and there is an increasing proportion of individuals with asymptomatic and SM infection. Transmission patterns become heterogeneous, and imported cases become important obstacles to complete parasite elimination. Several biological properties of P. vivax that enhance the transmissibility of this species relative to P. falciparum present major obstacles to elimination. A key operational challenge is identifying the hypnozoite reservoirs of P. vivax that define the transmission dynamics of the parasite and providing radical cure of all stages of the parasite to infected individuals. In addition to conventional methods, new surveillance tools offer unique insights into the complex and evolving micro-epidemiology of the parasite, although they are currently largely constrained to research-based activities. Standard and uPCR-based methods can help to quantify the hidden P. vivax reservoir, beyond patients with symptomatic illness. Genetic analyses can reveal the diversity and connectivity between individual infections and populations that can inform malaria control activities. In areas of very low incidence, novel serological methods offer promise for measuring populations exposed to recent infection, at high risk of relapse, so that individuals and populations can be targeted for test and treat strategies and even mass drug administration. Many of these tools are in their infancy, and their cost-effectiveness and operational utility have yet to be defined (Table 3), but their development and application to P. vivax will be critical in achieving the ultimate and timely elimination of the parasite.

Table 3
Challenges with conventional and new surveillance methods.
ObjectiveApproachChallenges
Monitoring the impact of public health interventionsPassive case detection with RDTs or LMLack of sensitivity in detecting low-level parasitaemia, failure to detect individuals with asymptomatic infections, and reliance on treatment-seeking behaviour
PCR-based parasite detectionNot feasible or cost-effective in most endemic settings
Parasite genotyping to gauge transmission intensityUtility and cost effectiveness of genetic indices need to be defined in different endemic settings [21]
SerologyRequires validation of responses in different endemic settings to predict individuals at increased risk of recurrent P. vivax parasitaemia
Detecting submicroscopic and asymptomatic infectionsCommunity surveys using ultrasenstive diagnostics such PCR-based parasite detection and next generation RDTsHigh demand on resources and challenges to produce timely results
Identifying foci of infectionGeospatial surveillance and case investigationThe liver-stage reservoir constrains accurate linkage between P. vivax infections
Parasite genotypingGenetic data provide a direct link between infections, but the cost efficacy of this approach at large scale needs to be evaluated [21]
Monitoring treatment efficacyTESLogistically challenging and difficult to classify the origin of recurrences (recrudescence, reinfection, or relapse)
Genotyping parasite drug resistance markersNo validated molecular markers of the most widely use blood- or liver-stage antimalarial drugs for P. vivax [112]
Parasite genotyping to improve recurrence classificationLimited number of clinical isolates are amenable to whole genome sequencing. Parsimonious genotyping and modelling approaches are being developed to characterise the IBD between infections [33]
Identifying imported casesTravel historyPersistent blood stage infections and long-latency liver stages constrain the accuracy of assumptions based on travel history
GenotypingScope for detecting long-distance imported cases, but currently unable to distinguish local from imported cases in porous border regions [88]

IBD, identity by descent; LM, light microscopy; RDT, rapid diagnostic test; TES, therapeutic efficacy studies.

Acknowledgements

We thank our colleagues for helpful discussions.

Abbreviations

CSPcircumsporozoite protein
ELISAenzyme-linked immunosorbent assay
IBDidentity by descent
LMlight microscopy
LSAliver-stage antigen
NAATnucleic acid amplification–based test
RDTrapid diagnostic test
RT-qPCRreverse transcription quantitative PCR
SMsubmicroscopic
SNPsingle nucleotide polymorphism
TEStherapeutic efficacy surveys
uPCRultrasensitive PCR

References

WHO. World Malaria Report 2019. World Health Organization: Geneva; 2019 [cited 2021 Feb 7]. Available from: https://www.who.int/publications/i/item/9789241565721

KEBattle, TCDLucas, MNguyen, REHowes, AKNandi, KATwohig, et al. Mapping the global endemicity and clinical burden of Plasmodium vivax, 2000–17: a spatial and temporal modelling study. Lancet. 2019;394(10195):33243. Epub 2019/06/24. 10.1016/S0140-6736(19)31096-7

RNPrice, RJCommons, KEBattle, KThriemer, KMendis. Plasmodium vivax in the Era of the Shrinking P. falciparum Map. Trends Parasitol. 2020;36(6):56070. 10.1016/j.pt.2020.03.009 .

RNPrice, Lvon Seidlein, NValecha, FNosten, JKBaird, NJWhite. Global extent of chloroquine-resistant Plasmodium vivax: a systematic review and meta-analysis. Lancet Infect Dis. 2014;14(10):98291. 10.1016/S1473-3099(14)70855-2

NJWhite. Why Do Some Primate Malarias Relapse? Trends Parasitol. 2016;32(12):91820. 10.1016/j.pt.2016.08.014

RJCommons, JASimpson, JWatson, NJWhite, RNPrice. Estimating the Proportion of Plasmodium vivax Recurrences Caused by Relapse: A Systematic Review and Meta-Analysis. Am J Trop Med Hyg. 2020;103(3):10949. 10.4269/ajtmh.20-0186

SDini, NMDouglas, JRPoespoprodjo, EKenangalem, PSugiarto, IDPlumb, et al. The risk of morbidity and mortality following recurrent malaria in Papua, Indonesia: a retrospective cohort study. BMC Med. 2020;18(1):28. 10.1186/s12916-020-1497-0

PMichon, JLCole-Tobian, EDabod, SSchoepflin, JIgu, MSusapu, et al. The risk of malarial infections and disease in Papua New Guinean children. Am J Trop Med Hyg. 2007;76(6):9971008.

TJAnderson, BHaubold, JTWilliams, JGEstrada-Franco, LRichardson, RMollinedo, et al. Microsatellite markers reveal a spectrum of population structures in the malaria parasite Plasmodium falciparum. Mol Biol Evol. 2000;17(10):146782. 10.1093/oxfordjournals.molbev.a026247 .

10 

SAuburn, AEBarry. Dissecting malaria biology and epidemiology using population genetics and genomics. Int J Parasitol. 2017;47(2–3):7785. 10.1016/j.ijpara.2016.08.006 .

11 

AEBarry, AWaltmann, CKoepfli, CBarnadas, IMueller. Uncovering the transmission dynamics of Plasmodium vivax using population genetics. Pathog Glob Health. 2015;109(3):14252. 10.1179/2047773215Y.0000000012

12 

YHamedi, KSharifi-Sarasiabi, FDehghan, RSafari, STo, IHandayuni, et al. Molecular Epidemiology of P. vivax in Iran: High Diversity and Complex Sub-Structure Using Neutral Markers, but No Evidence of Y976F Mutation at pvmdr1. PLoS ONE. 2016;11(11):e0166124. 10.1371/journal.pone.0166124

13 

SGunawardena, MUFerreira, GMKapilananda, DFWirth, NDKarunaweera. The Sri Lankan paradox: high genetic diversity in Plasmodium vivax populations despite decreasing levels of malaria transmission. Parasitology. 2014;141(7):88090. 10.1017/S0031182013002278 .

14 

DJWinter, MAPacheco, AFVallejo, RSSchwartz, MArevalo-Herrera, SHerrera, et al. Whole Genome Sequencing of Field Isolates Reveals Extensive Genetic Diversity in Plasmodium vivax from Colombia. PLoS Negl Trop Dis. 2015;9(12):e0004252. 10.1371/journal.pntd.0004252

15 

DNHupalo, ZLuo, AMelnikov, PLSutton, PRogov, AEscalante, et al. Population genomics studies identify signatures of global dispersal and drug resistance in Plasmodium vivax. Nat Genet. 2016;48(8):9538. 10.1038/ng.3588

16 

SAuburn, EDBenavente, OMiotto, RDPearson, RAmato, MJGrigg, et al. Genomic analysis of a pre-elimination Malaysian Plasmodium vivax population reveals selective pressures and changing transmission dynamics. Nat Commun. 2018;9(1):2585. 10.1038/s41467-018-04965-4

17 

EKenangalem, JRPoespoprodjo, NMDouglas, FHBurdam, KGdeumana, FChalfein, et al. Malaria morbidity and mortality following introduction of a universal policy of artemisinin-based treatment for malaria in Papua, Indonesia: A longitudinal surveillance study. PLoS Med. 2019;16(5):e1002815. 10.1371/journal.pmed.1002815

18 

ZPava, AMPuspitasari, ARumaseb, IHandayuni, LTrianty, RASUtami, et al. Molecular surveillance over 14 years confirms reduction of Plasmodium vivax and falciparum transmission after implementation of Artemisinin-based combination therapy in Papua, Indonesia. PLoS Negl Trop Dis. 2020;14(5):e0008295. Epub 2020/05/08. 10.1371/journal.pntd.0008295

19 

AWaltmann, CKoepfli, NTessier, SKarl, AFola, AWDarcy, et al. Increasingly inbred and fragmented populations of Plasmodium vivax associated with the eastward decline in malaria transmission across the Southwest Pacific. PLoS Negl Trop Dis. 2018;12(1):e0006146. 10.1371/journal.pntd.0006146 .

20 

JHKattenberg, ZRazook, RKeo, CKoepfli, CJennison, DLautu-Gumal, et al. Monitoring Plasmodium falciparum and Plasmodium vivax using microsatellite markers indicates limited changes in population structure after substantial transmission decline in Papua New Guinea. Mol Ecol. 2020;29(23):452541. 10.1111/mec.15654 .

21 

RNoviyanti, OMiotto, ABarry, JMarfurt, SSiegel, NThuy-Nhien, et al. Implementing parasite genotyping into national surveillance frameworks: feedback from control programmes and researchers in the Asia-Pacific region. Malar J. 2020;19(1):271. 10.1186/s12936-020-03330-5

22 

JMSa, SRKaslow, RRMoraes Barros, NFBrazeau, CMParobek, DTao, et al. Plasmodium vivax chloroquine resistance links to pvcrt transcription in a genetic cross. Nat Commun. 2019;10(1):4300. 10.1038/s41467-019-12256-9

23 

RDPearson, RAmato, SAuburn, OMiotto, JAlmagro-Garcia, CAmaratunga, et al. Genomic analysis of local variation and recent evolution in Plasmodium vivax. Nat Genet. 2016;48(8):95964. 10.1038/ng.3599

24 

RSuwanarusk, BRussell, MChavchich, FChalfein, EKenangalem, VKosaisavee, et al. Chloroquine resistant Plasmodium vivax: in vitro characterisation and association with molecular polymorphisms. PLoS ONE. 2007;2(10):e1089. 10.1371/journal.pone.0001089

25 

JMSa, MMYamamoto, CFernandez-Becerra, MFde Azevedo, JPapakrivos, BNaude, et al. Expression and function of pvcrt-o, a Plasmodium vivax ortholog of pfcrt, in Plasmodium falciparum and Dictyostelium discoideum. Mol Biochem Parasitol. 2006;150(2):21928. 10.1016/j.molbiopara.2006.08.006 .

26 

GCMelo, WMMonteiro, AMSiqueira, SRSilva, BMMagalhaes, ACAlencar, et al. Expression levels of pvcrt-o and pvmdr-1 are associated with chloroquine resistance and severe Plasmodium vivax malaria in patients of the Brazilian Amazon. PLoS ONE. 2014;9(8):e105922. 10.1371/journal.pone.0105922

27 

GSnounou, HPBeck. The use of PCR genotyping in the assessment of recrudescence or reinfection after antimalarial drug treatment. Parasitol Today. 1998;14(11):4627. 10.1016/s0169-4758(98)01340-4 .

28 

MImwong, GSnounou, SPukrittayakamee, NTanomsing, JRKim, ANandy, et al. Relapses of Plasmodium vivax infection usually result from activation of heterologous hypnozoites. J Infect Dis. 2007;195(7):92733. 10.1086/512241 .

29 

NChen, AAuliff, KRieckmann, MGatton, QCheng. Relapses of Plasmodium vivax infection result from clonal hypnozoites activated at predetermined intervals. J Infect Dis. 2007;195(7):93441. 10.1086/512242 .

30 

ANCowell, HOValdivia, DKBishop, EAWinzeler. Exploration of Plasmodium vivax transmission dynamics and recurrent infections in the Peruvian Amazon using whole genome sequencing. Genome Med. 2018;10(1):52. 10.1186/s13073-018-0563-0

31 

JPopovici, LRFriedrich, SKim, SBin, VRun, DLek, et al. Genomic Analyses Reveal the Common Occurrence and Complexity of Plasmodium vivax Relapses in Cambodia. MBio. 2018;9(1):e0188817. 10.1128/mBio.01888-17(1)

32 

ATBright, MJManary, RTewhey, EMArango, TWang, NJSchork, et al. A high resolution case study of a patient with recurrent Plasmodium vivax infections shows that relapses were caused by meiotic siblings. PLoS Negl Trop Dis. 2014;8(6):e2882. 10.1371/journal.pntd.0002882

33 

ARTaylor, JAWatson, CSChu, KPuaprasert, JDuanguppama, NPJDay, et al. Resolving the cause of recurrent Plasmodium vivax malaria probabilistically. Nat Commun. 2019;10(1):5595. 10.1038/s41467-019-13412-x

34 

TBousema, JTGriffin, RWSauerwein, DLSmith, TSChurcher, WTakken, et al. Hitting hotspots: spatial targeting of malaria for control and elimination. PLoS Med. 2012;9(1):e1001165. 10.1371/journal.pmed.1001165

35 

WHO. Overview of malaria surveillance. [cited 2020 Oct 22]. Available from: https://www.who.int/malaria/areas/surveillance/overview/en/.

36 

AWesolowski, ARTaylor, HHChang, RVerity, STessema, JABailey, et al. Mapping malaria by combining parasite genomic and epidemiologic data. BMC Med. 2018;16(1):190. 10.1186/s12916-018-1181-9

37 

HHChang, AWesolowski, ISinha, CGJacob, AMahmud, DUddin, et al. Mapping imported malaria in Bangladesh using parasite genetic and human mobility data. Elife. 2019;8:e43481. 10.7554/eLife.43481

38 

STessema, AWesolowski, AChen, MMurphy, JWilheim, ARMupiri, et al. Using parasite genetic and human mobility data to infer local and cross-border malaria connectivity in Southern Africa. Elife. 2019;8:e43510. 10.7554/eLife.43510

39 

KAGray, SDowd, LBain, ABobogare, LWini, GDShanks, et al. Population genetics of Plasmodium falciparum and Plasmodium vivax and asymptomatic malaria in Temotu Province, Solomon Islands. Malar J. 2013;12:429. 10.1186/1475-2875-12-429

40 

CKoepfli, LTiminao, TAntao, AEBarry, PSiba, IMueller, et al. A Large Reservoir and Little Population Structure in the South Pacific. PLoS ONE. 2013;8(6):e66041. 10.1371/journal.pone.0066041

41 

RNoviyanti, FCoutrier, RAUtami, HTrimarsanto, YKTirta, LTrianty, et al. Contrasting Transmission Dynamics of Co-endemic Plasmodium vivax and P. falciparum: Implications for Malaria Control and Elimination. PLoS Negl Trop Dis. 2015;9(5):e0003739. 10.1371/journal.pntd.0003739

42 

MUFerreira, NDKarunaweera, Mda Silva-Nunes, NSda Silva, DFWirth, DLHartl. Population structure and transmission dynamics of Plasmodium vivax in rural Amazonia. J Infect Dis. 2007;195(8):121826. 10.1086/512685 .

43 

POrjuela-Sanchez, JMSa, MCBrandi, PTRodrigues, MSBastos, CAmaratunga, et al. Higher microsatellite diversity in Plasmodium vivax than in sympatric Plasmodium falciparum populations in Pursat, Western Cambodia. Exp Parasitol. 2013;134(3):31826. 10.1016/j.exppara.2013.03.029

44 

MAPacheco, KASchneider, NCespedes, SHerrera, MArevalo-Herrera, AAEscalante. Limited differentiation among Plasmodium vivax populations from the northwest and to the south Pacific Coast of Colombia: A malaria corridor? PLoS Negl Trop Dis. 2019;13(3):e0007310. 10.1371/journal.pntd.0007310

45 

PManrique, JMiranda-Alban, JAlarcon-Baldeon, RRamirez, GCarrasco-Escobar, HHerrera, et al. Microsatellite analysis reveals connectivity among geographically distant transmission zones of Plasmodium vivax in the Peruvian Amazon: A critical barrier to regional malaria elimination. PLoS Negl Trop Dis. 2019;13(11):e0007876. 10.1371/journal.pntd.0007876

46 

PVan den Eede, GVan der Auwera, CDelgado, THuyse, VESoto-Calle, DGamboa, et al. Multilocus genotyping reveals high heterogeneity and strong local population structure of the Plasmodium vivax population in the Peruvian Amazon. Malar J. 2010;9:151. 10.1186/1475-2875-9-151

47 

SMChenet, KASchneider, LVillegas, AAEscalante. Local population structure of Plasmodium: impact on malaria control and elimination. Malar J. 2012;11:412. 10.1186/1475-2875-11-412

48 

RBiek, OGPybus, JOLloyd-Smith, XDidelot. Measurably evolving pathogens in the genomic era. Trends Ecol Evol. 2015;30(6):30613. 10.1016/j.tree.2015.03.009

49 

S.Wright Isolation by Distance. Genetics. 1943;28(2):11438.

50 

ARTaylor, SFSchaffner, GCCerqueira, SCNkhoma, TJCAnderson, KSriprawat, et al. Quantifying connectivity between local Plasmodium falciparum malaria parasite populations using identity by descent. PLoS Genet. 2017;13(10):e1007065. 10.1371/journal.pgen.1007065

51 

ARTaylor, DFEcheverry, TJCAnderson, DENeafsey, COBuckee. Identity-by-descent with uncertainty characterises connectivity of Plasmodium falciparum populations on the Colombian-Pacific coast. PLoS Genet. 2020;16(11):e1009101. 10.1371/journal.pgen.1009101

52 

NMAnstey, NMDouglas, JRPoespoprodjo, RNPrice. Plasmodium vivax: clinical spectrum, risk factors and pathogenesis. Adv Parasitol. 2012;80:151201. 10.1016/B978-0-12-397900-1.00003-7 .

53 

MLGatton, QCheng. Evaluation of the pyrogenic threshold for Plasmodium falciparum malaria in naive individuals. Am J Trop Med Hyg. 2002;66(5):46773. 10.4269/ajtmh.2002.66.467 .

54 

S.Kitchen Symptomatology: general considerations. In: MFBoyd, editor. Philadelphia: W.B. Saunders; 1949.

55 

BEBarber, TWilliam, MJGrigg, UParameswaran, KAPiera, RNPrice, et al. Parasite biomass-related inflammation, endothelial activation, microvascular dysfunction and disease severity in vivax malaria. PLoS Pathog. 2015;11(1):e1004558. Epub 2015/01/09. 10.1371/journal.ppat.1004558

56 

SKho, BAndries, JRPoespoprodjo, RJCommons, PAIShanti, EKenangalem, et al. High Risk of Plasmodium vivax Malaria Following Splenectomy in Papua, Indonesia. Clin Infect Dis. 2019;68(1):5160. Epub 2018/05/18. 10.1093/cid/ciy403

57 

MImwong, KStepniewska, RTripura, TJPeto, KMLwin, BVihokhern, et al. Numerical Distributions of Parasite Densities During Asymptomatic Malaria. J Infect Dis. 2016;213(8):13229. 10.1093/infdis/jiv596

58 

IHarris, WWSharrock, LMBain, KAGray, ABobogare, LBoaz, et al. A large proportion of asymptomatic Plasmodium infections with low and sub-microscopic parasite densities in the low transmission setting of Temotu Province, Solomon Islands: challenges for malaria diagnostics in an elimination setting. Malar J. 2010;9:254. Epub 2010/09/09. 1475-2875-9-254 [pii]. 10.1186/1475-2875-9-254

59 

QCheng, JCunningham, MLGatton. Systematic review of sub-microscopic P. vivax infections: prevalence and determining factors. PLoS Negl Trop Dis. 2015;9(1):e3413. 10.1371/journal.pntd.0003413 eCollection 2015 Jan. .

60 

SBarbosa, ABGozze, NFLima, CLBatista, SBastos Mda, VCNicolete, et al. Epidemiology of disappearing Plasmodium vivax malaria: a case study in rural Amazonia. PLoS Negl Trop Dis. 2014;8(8):e3109. 10.1371/journal.pntd.0003109

61 

AFVallejo, JGarcia, ABAmado-Garavito, MArevalo-Herrera, SHerrera. Plasmodium vivax gametocyte infectivity in sub-microscopic infections. Malar J. 2016;15:48. 10.1186/s12936-016-1104-1

62 

CMMoreira, MAbo-Shehada, RNPrice, CJDrakeley. A systematic review of sub-microscopic Plasmodium vivax infection. Malar J. 2015;14:360. 10.1186/s12936-015-0884-z

63 

MImwong, TNNguyen, RTripura, TJPeto, SJLee, KMLwin, et al. The epidemiology of subclinical malaria infections in South-East Asia: findings from cross-sectional surveys in Thailand-Myanmar border areas, Cambodia, and Vietnam. Malar J. 2015;14:381. 10.1186/s12936-015-0906-x

64 

TBousema, CDrakeley. Epidemiology and infectivity of Plasmodium falciparum and Plasmodium vivax gametocytes in relation to malaria control and elimination. Clin Microbiol Rev. 2011;24(2):377410. 10.1128/CMR.00051-10

65 

RTripura, TJPeto, JChalk, SJLee, PSirithiranont, CNguon, et al. Persistent Plasmodium falciparum and Plasmodium vivax infections in a western Cambodian population: implications for prevention, treatment and elimination strategies. Malar J. 2016;15:181. Epub 2016/03/26. 10.1186/s12936-016-1224-7 [pii].

66 

KPhommasone, Fvan Leth, MImwong, GHenriques, TPongvongsa, BAdhikari, et al. The use of ultrasensitive quantitative-PCR to assess the impact of primaquine on asymptomatic relapse of Plasmodium vivax infections: a randomized, controlled trial in Lao PDR. Malar J. 2020;19(1):4. Epub 2020/01/05. 10.1186/s12936-019-3091-5 [pii].

67 

RWampfler, FMwingira, SJavati, LRobinson, IBetuela, PSiba, et al. Strategies for detection of Plasmodium species gametocytes. PLoS ONE. 2013;8(9):e76316. 10.1371/journal.pone.0076316 eCollection 2013. .

68 

FGTadesse, HCSlater, WChali, KTeelen, KLanke, MBelachew, et al. The Relative Contribution of Symptomatic and Asymptomatic Plasmodium vivax and Plasmodium falciparum Infections to the Infectious Reservoir in a Low-Endemic Setting in Ethiopia. Clin Infect Dis. 2018;66(12):188391. 10.1093/cid/cix1123 .

69 

ERovira-Vallbona, JJContreras-Mancilla, RRamirez, MGuzman-Guzman, GCarrasco-Escobar, ALlanos-Cuentas, et al. Predominance of asymptomatic and sub-microscopic infections characterizes the Plasmodium gametocyte reservoir in the Peruvian Amazon. PLoS Negl Trop Dis. 2017;11(7):e0005674. 10.1371/journal.pntd.0005674

70 

AWaltmann, AWDarcy, IHarris, CKoepfli, JLodo, VVahi, et al. High Rates of Asymptomatic, Sub-microscopic Plasmodium vivax Infection and Disappearing Plasmodium falciparum Malaria in an Area of Low Transmission in Solomon Islands. PLoS Negl Trop Dis. 2015;9(5):e0003758. 10.1371/journal.pntd.0003758

71 

CKoepfli, LJRobinson, PRarau, MSalib, NSambale, RWampfler, et al. Blood-Stage Parasitaemia and Age Determine Plasmodium falciparum and P. vivax Gametocytaemia in Papua New Guinea. PLoS ONE. 2015;10(5):e0126747. 10.1371/journal.pone.0126747

72 

FPAlves, LHGil, MTMarrelli, PERibolla, EPCamargo, LHDa Silva. Asymptomatic carriers of Plasmodium spp. as infection source for malaria vector mosquitoes in the Brazilian Amazon. J Med Entomol. 2005;42(5):7779. 10.1093/jmedent/42.5.777 .

73 

CJJoyner, CFABrito, CLSaney, RJoice Cordy, MLSmith, SALapp, et al. Humoral immunity prevents clinical malaria during Plasmodium relapses without eliminating gametocytes. PLoS Pathog. 2019;15(9):e1007974. Epub 2019/09/20. 10.1371/journal.ppat.1007974 PPATHOGENS-D-18-02403 [pii].

74 

AMachado Siqueira, BMLopes Magalhaes, GCardoso Melo, MFerrer, PCastillo, LMartin-Jaular, et al. Spleen rupture in a case of untreated Plasmodium vivax infection. PLoS Negl Trop Dis. 2012;6(12):e1934. 10.1371/journal.pntd.0001934

75 

BBaro, KDeroost, TRaiol, MBrito, ACAlmeida, Ade Menezes-Neto, et al. Plasmodium vivax gametocytes in the bone marrow of an acute malaria patient and changes in the erythroid miRNA profile. PLoS Negl Trop Dis. 2017;11(4):e0005365. 10.1371/journal.pntd.0005365

76 

CCotter, HJSturrock, MSHsiang, JLiu, AAPhillips, JHwang, et al. The changing epidemiology of malaria elimination: new strategies for new challenges. Lancet. 2013;382(9895):90011. 10.1016/S0140-6736(13)60310-4 Epub 2013 Apr 15. .

77 

TYangzom, CSGueye, RNamgay, GNGalappaththy, KThimasarn, RGosling, et al. Malaria control in Bhutan: case study of a country embarking on elimination. Malar J. 2012;11:9. 10.1186/1475-2875-11-9

78 

RRAbeyasinghe, GNGalappaththy, CSmith Gueye, JGKahn, RGFeachem. Malaria control and elimination in Sri Lanka: documenting progress and success factors in a conflict setting. PLoS ONE. 2012;7(8):e43162. 10.1371/journal.pone.0043162

79 

AErhart, DTNgo, VKPhan, TTTa, CVan Overmeir, NSpeybroeck, et al. Epidemiology of forest malaria in central Vietnam: a large scale cross-sectional survey. Malar J. 2005;4:58. 10.1186/1475-2875-4-58

80 

ABhumiratana, PSorosjinda-Nunthawarasilp, WKaewwaen, PManeekan, SPimnon. Malaria-associated rubber plantations in Thailand. Travel Med Infect Dis. 2013;11(1):3750. 10.1016/j.tmaid.2012.11.002 .

81 

RChuquiyauri, MParedes, PPenataro, STorres, SMarin, ATenorio, et al. Socio-demographics and the development of malaria elimination strategies in the low transmission setting. Acta Trop. 2012;121(3):292302. 10.1016/j.actatropica.2011.11.003

82 

HHiwat, LSHardjopawiro, WTakken, LVillegas. Novel strategies lead to pre-elimination of malaria in previously high-risk areas in Suriname, South America. Malar J. 2012;11:10. 10.1186/1475-2875-11-10

83 

GCarrasco-Escobar, MCCastro, JLBarboza, JRuiz-Cabrejos, ALlanos-Cuentas, JMVinetz, et al. Use of open mobile mapping tool to assess human mobility traceability in rural offline populations with contrasting malaria dynamics. PeerJ. 2019;7:e6298. 10.7717/peerj.6298

84 

VRougeron, EElguero, CArnathau, BAcuna Hidalgo, PDurand, SHouze, et al. Human Plasmodium vivax diversity, population structure and evolutionary origin. PLoS Negl Trop Dis. 2020;14(3):e0008072. 10.1371/journal.pntd.0008072

85 

MLBaniecki, ALFaust, SFSchaffner, DJPark, KGalinsky, RFDaniels, et al. Development of a single nucleotide polymorphism barcode to genotype Plasmodium vivax infections. PLoS Negl Trop Dis. 2015;9(3):e0003539. 10.1371/journal.pntd.0003539

86 

CKoepfli, PTRodrigues, TAntao, POrjuela-Sanchez, PVan den Eede, DGamboa, et al. Plasmodium vivax Diversity and Population Structure across Four Continents. PLoS Negl Trop Dis. 2015;9(6):e0003872. 10.1371/journal.pntd.0003872

87 

JHKattenberg, ZRazook, RKeo, CKoepfli, CJennison, DLautu-Ninda, et al. Monitoring of Plasmodium falciparum and Plasmodium vivax using microsatellite markers indicates limited changes in population structure after substantial transmission decline in Papua New Guinea. bioRxiv [Preprint]. 2020. 10.1111/mec.15654

88 

HTrimarsanto, RAmato, RDPearson, ESutanto, RNoviyanti, LTrianty, et al. A molecular barcode and online tool to identify and map imported infection with Plasmodium vivax. bioRxiv [Preprint]. 2019. 10.1101/776781

89 

DWang, SLi, ZCheng, NXiao, CCotter, JHwang, et al. Transmission Risk from Imported Plasmodium vivax Malaria in the China-Myanmar Border Region. Emerg Infect Dis. 2015;21(10):18614. 10.3201/eid2110.150679

90 

ELo, NLam, EHemming-Schroeder, JNguyen, GZhou, MCLee, et al. Frequent Spread of Plasmodium vivax Malaria Maintains High Genetic Diversity at the Myanmar-China Border, Without Distance and Landscape Barriers. J Infect Dis. 2017;216(10):125463. 10.1093/infdis/jix106

91 

SWangchuk, TDrukpa, KPenjor, TPeldon, YDorjey, KDorji, et al. Where chloroquine still works: the genetic make-up and susceptibility of Plasmodium vivax to chloroquine plus primaquine in Bhutan. Malar J. 2016;15(1):277. 10.1186/s12936-016-1320-8

92 

PCorran, PColeman, ERiley, CDrakeley. Serology: a robust indicator of malaria transmission intensity? Trends Parasitol. 2007;23(12):57582. 10.1016/j.pt.2007.08.023 .

93 

CDrakeley, JCook. Chapter 5. Potential contribution of sero-epidemiological analysis for monitoring malaria control and elimination: historical and current perspectives. Adv Parasitol. 2009;69:299352. 10.1016/S0065-308X(09)69005-9 .

94 

ZMIdris, CWChan, MMohammed, MKalkoa, GTaleo, KJunker, et al. Serological measures to assess the efficacy of malaria control programme on Ambae Island, Vanuatu. Parasit Vectors. 2017;10(1):204. 10.1186/s13071-017-2139-z

95 

TBousema, CDrakeley, SGesase, RHashim, SMagesa, FMosha, et al. Identification of hot spots of malaria transmission for targeted malaria control. J Infect Dis. 2010;201(11):176474. 10.1086/652456 .

96 

KKerkhof, VSluydts, LWillen, SKim, LCanier, SHeng, et al. Serological markers to measure recent changes in malaria at population level in Cambodia. Malar J. 2016;15(1):529. 10.1186/s12936-016-1576-z

97 

RLDewasurendra, JNDias, NSepulveda, GSGunawardena, NChandrasekharan, CDrakeley, et al. Effectiveness of a serological tool to predict malaria transmission intensity in an elimination setting. BMC Infect Dis. 2017;17(1):49. 10.1186/s12879-016-2164-0

98 

DLDoolan, CDobano, JKBaird. Acquired immunity to malaria. Clin Microbiol Rev. 2009;22(1):1336. 10.1128/CMR.00025-08

99 

CJDrakeley, PHCorran, PGColeman, JETongren, SLMcDonald, ICarneiro, et al. Estimating medium- and long-term trends in malaria transmission by using serological markers of malaria exposure. Proc Natl Acad Sci U S A. 2005;102(14):510813. 10.1073/pnas.0408725102 .

100 

JBSarr, EOrlandi-Pradines, SFortin, CSow, SCornelie, FRogerie, et al. Assessment of exposure to Plasmodium falciparum transmission in a low endemicity area by using multiplex fluorescent microsphere-based serological assays. Parasit Vectors. 2011;4:212. 10.1186/1756-3305-4-212

101 

JSkinner, CYHuang, MWaisberg, PLFelgner, OKDoumbo, AOngoiba, et al. Plasmodium falciparum Gametocyte-Specific Antibody Profiling Reveals Boosting through Natural Infection and Identifies Potential Markers of Gametocyte Exposure. Infect Immun. 2015;83(11):422936. 10.1128/IAI.00644-15

102 

APoinsignon, SCornelie, MMestres-Simon, ALanfrancotti, MRossignol, DBoulanger, et al. Novel peptide marker corresponding to salivary protein gSG6 potentially identifies exposure to Anopheles bites. PLoS ONE. 2008;3(6):e2472. 10.1371/journal.pone.0002472

103 

ABSagna, JBSarr, LGaayeb, PMDrame, MONdiath, SSenghor, et al. gSG6-P1 salivary biomarker discriminates micro-geographical heterogeneity of human exposure to Anopheles bites in low and seasonal malaria areas. Parasit Vectors. 2013;6:68. 10.1186/1756-3305-6-68

104 

BLondono-Renteria, PMDrame, TWeitzel, RRosas, CGripping, JCCardenas, et al. An. gambiae gSG6-P1 evaluation as a proxy for human-vector contact in the Americas: a pilot study. Parasit Vectors. 2015;8:533. 10.1186/s13071-015-1160-3

105 

PYa-Umphan, DCerqueira, DMParker, GCottrell, APoinsignon, FRemoue, et al. Use of an Anopheles Salivary Biomarker to Assess Malaria Transmission Risk Along the Thailand-Myanmar Border. J Infect Dis 2017;215(3):396404. 10.1093/infdis/jiw543

106 

EJMPollard, CPatterson, TLRussell, AApairamo, JOscar, BArca, et al. Human exposure to Anopheles farauti bites in the Solomon Islands is not associated with IgG antibody response to the gSG6 salivary protein of Anopheles gambiae. Malar J. 2019;18(1):334. 10.1186/s12936-019-2975-8

107 

RJLongley, MTWhite, ETakashima, JBrewster, MMorita, MHarbers, et al. Development and validation of serological markers for detecting recent Plasmodium vivax infection. Nat Med. 2020;26(5):7419. 10.1038/s41591-020-0841-4 .

108 

D.Payne Spread of chloroquine resistance in Plasmodium falciparum. Parasitol Today. 1987;3(8):2416. 10.1016/0169-4758(87)90147-5 .

109 

KHRieckmann, DRDavis, DCHutton. Plasmodium vivax resistance to chloroquine? Lancet. 1989;2(8673):11834. 10.1016/s0140-6736(89)91792-3 .

110 

KASchneider, AAEscalante. Fitness components and natural selection: why are there different patterns on the emergence of drug resistance in Plasmodium falciparum and Plasmodium vivax? Malar J. 2013;12:15. 10.1186/1475-2875-12-15

111 

MJGrigg, TWilliam, JMenon, BEBarber, CSWilkes, GSRajahram, et al. Efficacy of Artesunate-mefloquine for Chloroquine-resistant Plasmodium vivax Malaria in Malaysia: An Open-label, Randomized, Controlled Trial. Clin Infect Dis. 2016;62(11):140311. 10.1093/cid/ciw121

112 

RNPrice, SAuburn, JMarfurt, QCheng. Phenotypic and genotypic characterisation of drug-resistant Plasmodium vivax. Trends Parasitol. 2012;28(11):5229. 10.1016/j.pt.2012.08.005 .

113 

HTrimarsanto, EDBenavente, RNoviyanti, RAUtami, LTrianty, ZPava, et al. VivaxGEN: An open access platform for comparative analysis of short tandem repeat genotyping data in Plasmodium vivax Populations. PLoS Negl Trop Dis. 2017;11(3):e0005465. 10.1371/journal.pntd.0005465 .

114 

CKoepfli, IMueller, JMarfurt, MGoroti, ASie, OOa, et al. Evaluation of Plasmodium vivax genotyping markers for molecular monitoring in clinical trials. J Infect Dis. 2009;199(7):107480. 10.1086/597303 .

115 

NRAbdullah, BEBarber, TWilliam, NANorahmad, URSatsu, PKMuniandy, et al. Plasmodium vivax population structure and transmission dynamics in Sabah Malaysia. PLoS ONE. 2013;8(12):e82553. 10.1371/journal.pone.0082553

116 

SGetachew, STo, HTrimarsanto, KThriemer, TGClark, BPetros, et al. Variation in Complexity of Infection and Transmission Stability between Neighbouring Populations of Plasmodium vivax in Southern Ethiopia. PLoS ONE. 2015;10(10):e0140780. 10.1371/journal.pone.0140780

117 

YHamedi, MNateghpour, PTan-ariya, MTiensuwan, USilachamroon, SLooareesuwan. Plasmodium vivax malaria in Southeast Iran in 1999–2001: establishing the response to chloroquine in vitro and in vivo. Southeast Asian J Trop Med Public Health. 2002;33(3):5128. .

118 

JYKim, YKGoo, YGZo, SYJi, HTrimarsanto, STo, et al. Further Evidence of Increasing Diversity of Plasmodium vivax in the Republic of Korea in Recent Years. PLoS ONE. 2016;11(3):e0151514. 10.1371/journal.pone.0151514

119 

MWHtun, NCNMon, KMAye, CMHlaing, MPKyaw, IHandayuni, et al. Chloroquine efficacy for Plasmodium vivax in Myanmar in populations with high genetic diversity and moderate parasite gene flow. Malar J. 2017;16(1):281. 10.1186/s12936-017-1912-y