A simulation model of the within-host dynamics of Plasmodium vivax infection
© Kerlin and Gatton; licensee BioMed Central. 2015
Received: 21 October 2014
Accepted: 25 January 2015
Published: 5 February 2015
The benign reputation of Plasmodium vivax is at odds with the burden and severity of the disease. This reputation, combined with restricted in vitro techniques, has slowed efforts to gain an understanding of the parasite biology and interaction with its human host.
A simulation model of the within-host dynamics of P. vivax infection is described, incorporating distinctive characteristics of the parasite such as the preferential invasion of reticulocytes and hypnozoite production. The developed model is fitted using digitized time-series’ from historic neurosyphilis studies, and subsequently validated against summary statistics from a larger study of the same population. The Chesson relapse pattern was used to demonstrate the impact of released hypnozoites.
The typical pattern for dynamics of the parasite population is a rapid exponential increase in the first 10 days, followed by a gradual decline. Gametocyte counts follow a similar trend, but are approximately two orders of magnitude lower. The model predicts that, on average, an infected naïve host in the absence of treatment becomes infectious 7.9 days post patency and is infectious for a mean of 34.4 days. In the absence of treatment, the effect of hypnozoite release was not apparent as newly released parasites were obscured by the existing infection.
The results from the model provides useful insights into the dynamics of P. vivax infection in human hosts, in particular the timing of host infectiousness and the role of the hypnozoite in perpetuating infection.
KeywordsPlasmodium vivax Simulation model Preferential invasion Hypnozoite Within-host dynamics
Historically, the high profile and lethal reputation of Plasmodium falciparum has meant that the bulk of malaria research and management has focused on this parasite; comparatively, the importance of Plasmodium vivax (and other species of human malaria) as a human pathogen has been neglected [1,2]. Yet despite a preconception that P. vivax is a benign infection, at least in comparison with P. falciparum, vivax malaria is a significant burden to nearly 40% of the world’s population, with estimates of between 132 and 391 million cases per year ; the greatest burden (52%) observed is in South and East Asia. Multiple reports have linked infection with P. vivax to patient death, often as a result of fatal complications and particularly in pregnant women and young children [3,4].
As the preconception of P. vivax as a benign disease agent is challenged, key questions remain regarding the basic ecology and dynamics of P. vivax infection .
The capacity of researchers to address these questions has been hampered by ongoing difficulties in producing an ex vitro culture system. Such a system was developed for P. falciparum in the 1970s  and drove subsequent research which provided a wealth of information on aspects of the parasite and infection . With limited opportunities to investigate the dynamics of the parasite in the laboratory, modelling provides a useful alternative to explore aspects of P. vivax infection.
This manuscript describes the development of a stochastic simulation model of the within-host dynamics of a P. vivax infection in a naïve host, enabling elucidation of a number of significant questions regarding the interaction between host and parasite. It is intended that the development of this within-host model will subsequently have utility as the basis for a between-host model which incorporates multiple hosts and vectors. As P. falciparum is far better studied than P. vivax, the dynamics of this parasite is used as the basis for the model. However, it is important to recognize that there are several significant features that differentiate P. vivax from P. falciparum. These are detailed below.
Preferential red blood cell invasion
Plasmodium falciparum, while showing a weak preference for younger red blood cells (RBCs), is generally assumed to be capable of invading all RBC age classes, whilst P. vivax shows a strong preference for the youngest RBCs, reticulocytes [7,8]. It has been hypothesized that preferential invasion may play a key role in regulating the abundance of parasites in vivax malaria; increased competition between merozoites for a limited population of susceptible RBCs can slow the growth of an infection, ensuring that the host remains alive sufficiently long to maximize the chances of subsequent generations reinfecting a feeding mosquito . Certainly the less discriminant P. falciparum is the human malaria species predominantly responsible for cases of severe disease and mortalities, while the other human Plasmodium species, which are more discerning in the RBCs they invade, are considered to be more benign disease agents .
Fever is the most common clinical manifestation of malaria infection, and is a primary mechanism by which the host immune system can control infection, with increased body temperatures associated with parasite death. Higher fevers generally lead to greater parasite mortality. Fever is triggered by the rupture of mature schizont stage parasites and associated RBC lysis, with the pyrogenic threshold representing the density of parasites required to stimulate a febrile response to infection. Vivax malaria is generally accepted to have a lower pyrogenic threshold than falciparum malaria [2,11], thus fewer P. vivax parasites are required to trigger a fever in the host.
Production of gametocytes
During the erythrocytic stage of parasite development in humans, a fraction of parasites commit to a sexual development pathway. Different cells will produce male and female gametocytes, which can then be transmitted to a mosquito during a blood meal. Plasmodium vivax gametocytes vary from their P. falciparum equivalents by their presence in the blood at detectable levels during the primary attack of the disease whereas P. falciparum gametocytes cannot usually be seen on microscopy until some 9–12 days post infection . It is unclear whether this difference is due to P. falciparum gametocytes sequestering in the bone marrow, or a longer gametocyte development time. The survival time of P. vivax gametocytes also appears to be significantly shorter than that of P. falciparum gametocytes .
Approximately 12–14 hours after invasion, RBCs parasitized by P. falciparum express a strain-specific variant surface antigen (P. falciparum erythrocyte membrane protein 1 PfEMP1), which mediates agglutination and adhesion of infected erythrocytes to the vascular endothelium and results in parasitized RBCs sequestering in the deep tissues, protected from splenic destruction [13,14]. PfEMP1 is encoded by a highly diverse family of genes, known as var genes, enabling rapid switching of presented antigens to evade the host immune response. While PfEMP1 and var gene switching have been the focus of intense study in P. falciparum, no clinically relevant equivalent has been identified in P. vivax .
Hypnozoites, an exo-erythrocytic parasite life cycle stage derived from the initial inoculum of sporozoites, are characteristic of P. vivax, but not P. falciparum [16,17]. The hypnozoite stage is believed to provide a mechanism for the parasite to survive in environments where the presence of the mosquito vector is intermittent (e.g. in Finland or the Korean Peninsula, where mosquitoes overwinter during the winter months [18,19]). It is generally assumed that the duration of this dormant stage is genetically determined, but it has also been hypothesized that mosquito bites, or a host immune response triggered by other infections may ‘alert’ the hypnozoite to the presence of the vector .
The basis of this work is a previously reported simulation model of preferential invasion of Plasmodia merozoites . In brief, the probability of merozoites successfully invading a new host red blood cell was assumed to be a function of the relative abundance of RBCs (anaemic status of the host or how likely is it that a free merozoite will encounter an RBC) and the relative abundance of susceptible RBC age classes (having encountered an RBC, how likely is it that the RBC is of a susceptible age class). As an infection progresses, the probability of a free merozoite encountering a susceptible RBC is expected to fall. Such a fall reduces the probability of a successful merozoite invasion. As it is widely recognized that P. vivax preferentially invades reticulocytes, merozoite invasion was restricted to 1 and 2 day old RBCs.
If the number of successful merozoites exceeds the number of available susceptible RBCs, the excess are removed from the system.
To develop the model of within-host dynamics this base model was augmented to incorporate the effects of fever, the production of gametocytes and the release of dormant hypnozoite stages of the parasite. These additions were created as separate modules. In the absence of an identifiable mechanism of clinically relevant antigenic variation, it is assumed that the clonal immunity described above is the only specific host immune response that develops during infection. The model uses a 24 hour (1 day) time step and the duration of the erythrocytic cycle in P. vivax was assumed to be 48 hours; variability in the duration of the erythrocytic cycle and differences in the duration of each stage of the erythrocytic lifecycle  were not considered.
Representation of the fever response is largely derived from the work of Gatton and Cheng on P. falciparum , with some changes to account for known differences with P. vivax. Pyrogenic threshold figures for P. vivax have been reported in a number of prior studies; 310 parasites/μL in southern Papua, Indonesia , 500 parasites/μL in analyses of children in Papua New Guinea , and 181 parasites/μL in Thailand . For the current model, a pyrogenic threshold of 300 ± 200 parasites/μL is adopted. This value is considerably lower than the comparative value of 1,644 - 2,500 parasites/μL for P. falciparum [11,27]. Since information on the dynamics of the threshold during an infection is missing for P. vivax, it is assumed this threshold is not reached for a period of 7 days, and introduce a log-linear increase over these first seven days, following work by Gatton and Cheng on the pyrogenic threshold of P. falciparum .
Following a review of the relevant literature [28,29], the model allows fever to impact on parasite survival such that late stage parasites are more susceptible to fever than early stage parasites (remembering that fever is most contemporaneous with early ring stage parasites, occurring shortly after schizont rupture), while gametocytes are considered less susceptible than ring stage parasites. Within the gametocyte sub-population, late stage gametocytes are considered more susceptible to fever than early stage gametocytes. Thus α = 1 for non-rupturing (late stage) asexual parasites, α = 2 for rupturing asexual parasites, α = 0.364 for early-stage gametocytes, and α = 0.667 for late-stage gametocytes.
It is assumed that gametocyte production begins at the onset of infection and that the gametocyte conversion rate (GCR) – the rate at which merozoites commit to the gametocyte development pathway – is log-normally distributed. A mean of the natural logarithm of GCR equal to −1.59, with a standard deviation of 0.13, is used based on published GCRs of 1.4% to 4.7% (mean of 3%)  and 4% . There is also considerable evidence that the gametocyte conversion rate is elevated under stressful conditions that are unfavourable for parasite survival . To incorporate this into the model, the GCR was allowed to rise by 2% on average when fever is present, following a previous model of gametocyte conversion in P. falciparum . All gametocytes are assumed to require two days to mature before they became infectious to mosquitoes [32,33], and all survived for a further three days thereafter (with no stochastic variation) [34,35]. This longevity is considerably shorter than reported values for P. falciparum, however the reduced longevity of P. vivax gametocytes is offset by a higher rate of production .
Given the desire to develop a model that can be subsequently augmented with a between-host model of transmission, the predicted the number of gametocytes required to ensure at least one male and one female is taken up with any mosquito blood meal was also considered. The sex ratio of P. vivax gametocytes is difficult to ascertain in the absence of a reliable culture system. To predict the density of gametocytes required to ensure the carriage of at least one male and one female in a mosquito blood meal, a generic sex ratio for malaria is used - one male gametocyte is assumed to be created for every 10 female gametocytes . Using a negative binomial distribution where the probability of finding a male gametocyte (p) is thus equal to 0.09, the predicted the number of gametocytes required to ensure at least one male and one female with 95% confidence is 31. Two assumptions are made in order to calculate this value: 1) a feeding mosquito will take a blood meal equivalent to 1 μL of blood , and 2) the distribution of gametocytes ingested by a feeding mosquito follows a negative binomial distribution with constant overdispersion, where k = 3.105 .
At the onset of infection, it is assumed that between one and four hypnozoites develop within the liver, each with a pre-programmed relapse date. This assumption was made on the basis of observations of infected hosts demonstrating recurrent relapses in the absence of reinfection; it is uncertain whether this repeated recurrence is due to multiple hypnozoites with different relapse dates developing during the initial infection, or new hypnozoites being created during each subsequent relapse event. Detailed data was available describing relapse patterns in three different regions, allowing the generation of three different relapse scenarios. For the purposes of this study the focus is on a Chesson style relapse pattern. Chloroquine studies from New Guinea and the Pacific describe a rapid pattern of relapse, with a duration of sequestration modelled using a combination of normal (87% - μ = 22.7, σ = 2.14) and binomial (13% - N = 30, p = 0.95) distributions (given a degree of bimodality in the source data), and providing a mean time between primary infection and first relapse, and between subsequent relapses, of 23.8 days . This short relapse time falls within the mean time to relapse predicted for parasites in the wider South-East Asia region . While the Chesson strain is largely confined to South-East Asia, this rapid pattern of relapse has resulted in the Chesson strain of P. vivax gaining a central role in experimental studies of the species .
It is assumed that in 40% of new infections only one hypnozoite is sequestered. A pre-programmed relapse time is assigned by a random draw from the probability distribution described above and converted to a relapse date. Two, three or four hypnozoites are sequestered in 30%, 20% and 10% of infections, respectively; where multiple hypnozoites are sequestered, each is also assigned a relapse date by a random draw, multiplied by the hypnozoite number (e.g. a second hypnozoite, randomly assigned a relapse time of 25 days, will be assigned a relapse date of 25 × 2 = 50 days after the new infection is created).
Once the relapse date is reached, 10,000 merozoites are released from the hypnozoite. These merozoites are treated in the same way as a new infection in the model.
Check to determine if fever will be triggered by the rupturing RBCs,
- 2.If schizonts have matured, or hypnozoites have reached their relapse date then merozoites are released from within rupturing schizonts or from relapsing hypnozoites, or both:
Merozoites compete to find a suitable host cell according to the preferential invasion module,
A fraction of the merozoites commit to the gametocyte pathway,
Any existing gametocytes age one day and gametocyte death occurs if necessary,
If fever has been triggered, parasite death is simulated with a proportion of the parasite burden removed from the each parasite age class.
Fitting the model
The model was fit using digitized copies of six published time series of parasitaemia, taken from neurosyphilis patients treated with P. vivax infection . Time series run for a maximum of 73 days before treatment, so the model simulation period was 80 days. Two variables in the model were not informed by the literature, and were instead fitted using available data. These variables were the strength of the clonal immunity response (β), and strength of the ‘killing’ power of fever (g). Gamma distributions were used to model each variable. After aligning the initial growth of the infection to account for the delay between exo-erythrocytic schizont rupture and patent parasitaemia, simulated annealing was used to fit the model to the mean of the six available time-series.
The fitted model was subsequently used to simulate a set of 1,000 infections, each within a malaria-naïve host, which were used to determine 95% confidence intervals for model output. These numbers were then compared to a more comprehensive set of summary statistics reported for 48 sporozoite induced P. vivax infections in the same neurosyphilis study (henceforth referred to as the ‘McKenzie summary statistics’, or MSS) to validate the model: peak parasitaemia, day of peak parasitaemia, day of peak fever, peak gametocytes and day of peak gametocytes. Statistics recorded for trophozoite-induced infections, while available, were not considered, as the current model was based on a single release of merozoites from a single hepatic schizont, and thus more closely resembled a sporozoite-induced infection.
A version of the fitted model was also run where chemotherapy was used to terminate infections. This model was included to clearly demonstrate the utility of the hypnozoite release component of the model. The treatment was assumed to kill all asexual parasites and gametocytes, and have no residual activity.
All analyses were conducted in the R statistical computing environment .
Comparisons between the model output and summary statistics for 48 sporozoite-induced infections (drawn from McKenzie et. al. )
Day of 1st parasite peak
Peak parasitaemia (/μL)
Day of 1st fever
Day of peak gametocytes
Peak gametocytemia (/μL)
The period of host infectivity can be determined by analysis of the density of gametocytes in the host. It is estimated that a density of 31 gametocytes/μL of host blood will ensure at least one male with 95% confidence. It is assumed that due to the skewed sex ratio in favour of females that the likelihood of a 1 μL sample containing 31 males (no females) is trivial in comparison to the likelihood of 31 females (no males). Given aggregation in the distribution of gametocytes in host blood, a density of 116.1 gametocytes/μL is required to ensure with 95% confidence that any particular blood meal of 1 μL will contain at least 31 gametocytes, and thus at least one male and one female gametocyte. The model predicts that, on average, an infected naïve host in the absence of treatment will exhibit gametocytes densities greater than 116.1/μL for a mean of 34.4 days (95% CI 22, 46 days), first reaching this threshold on average 7.9 days post patency (95% CI 5, 10).
A model of the infection dynamics of P. vivax in a naïve host has been developed. This model is founded on a model of preferential merozoite invasion incorporating non-variant clonal immunity that captures the high affinity displayed by P. vivax merozoites for reticulocytes, and the impact this may have on limiting parasite population growth. The model, augmented by the inclusion of fever, the production of gametocytes, and hypnozoite-mediated relapse, can be used to answer questions regarding the infectivity of hosts and offer suggestions as to the role of hypnozoite-mediated relapse in the evolution of drug resistance.
The model presented is a simple abstraction of a true infection. The model was fit to a small subset of data from neurosyphillis patients as no other time-course data was readily available. Comparisons between model output and available data thus require some caveats. For example, it is unrealistic to expect a perfect match between the model and the summary data: it is not clear how representative the six available time-series are of the neurosyphilis dataset as a whole, there is no way to quantify observation error in the summary data, and timings are dependent on the fitting process. The similarity between model output and the summary MSS data does however suggest that the six available time-series do provide a reasonable representation of the larger data set. Thus it is not surprising that the model is unable to predict the full variation present in patient data. It may be that fitting to the range of values for the six time series, rather than simply the mean, or adding additional parameters would improve model fit, as would having more than six time series to use for fitting. Despite this shortcoming, the model does appear to provide a reasonable description of the infection dynamics, and can subsequently still provide some useful insights into the biology of P. vivax.
The timing of the various peaks of infection are generally earlier in the MSS than in the model output, though in some instances (e.g. maximum day of first parasite peak, median/mean day of peak gametocytes) the estimates are similar. Some part of this is a consequence of the fitting process, but generally the cause appears to be the slow decay in the MSS following the peak. The model had difficulty replicating this decline, and has likely compensated by overshooting the peak, allowing a sharper decline to more closely match the MSS dataset. Thus, the model has overestimated the peak parasitaemia, in order to better fit the slower decline following the peak.
It is unsurprising that estimates of density were generally much larger in the model output than reported in MSS; MSS statistics are a product of microscopy of patient blood, which is subject to the transient nature of parasite detection in peripheral blood. Factors such as parasite sequestration and non-synchronized rupture of parasites impact on the parasite density observed. These factors are not accounted for in the model, therefore it is expected that the model would over-estimate the measured density.
One of the objections for preferential invasion as a mechanism for regulating parasite abundance is that if parasitaemia is suppressed, it follows that gametocyte abundance is also suppressed, reducing the probability of transmission to a new host. The model presented offers evidence that refutes this objection, demonstrating the possibility that preferential invasion can regulate parasite abundance while producing sufficient gametocytes to ensure an infectious host for a period of more than a month. It is further worth noting that the figures quoted here assume a worst case scenario of 1 male gametocyte produced for every 10 female gametocytes, when the generic sex ratio for malaria parasites suggest the proportion of male gametocytes could be as high as 1 in 5, greatly reducing the parasite density required to make a host infectious.
Hypnozoite-mediated relapse provides P. vivax infections a sequestering mechanism, a potential sanctuary where the parasite can effectively lie dormant and await potentially more advantageous circumstances. At this stage it is not definitively known whether the hypnozoite is responsive to signals indicating a beneficial environment to re-activate an infection, or whether the parasite relies on chance, programming an ‘alarm clock’ to promote reactivation at a future time. The most notable conclusion to be drawn from the presented model is that hypnozoite-mediated relapse can be greatly suppressed by an existing infection and associated immune response. It is only when the host immune response, or treatment, sufficiently reduce the pre-existing infection that a hypnozoite-mediated relapse can be detected. Parasites from hypnozoites released prior to this time are absorbed into the existing parasite pool.
Given the short time span of available patient parasitaemia simulated by the model (maximum of 80 days), the relapse pattern typical of the Chesson strain of vivax malaria, a relapse pattern that is extremely short compared to others observed, was demonstrated. If a longer relapse pattern, such as demonstrated in strains from East Timor  or Sri Lanka , were used in the model, the probability that a relapsing infection would be suppressed by an existing infection greatly diminishes due to the reduction in frequency of fever throughout an infection. This raises an interesting question: if hypnozoites evolved in response to environmental conditions where the presence of the mosquito vector is intermittent, what is the purpose of a short relapse pattern such as exhibited by the Chesson strain? The original infection is likely to still be active, and mosquito abundance is unlikely to significantly change in such a short time period. This interesting evolutionary question warrants further consideration.
In the absence of an ex vitro culture system for P. vivax, modelling provides a useful tool for exploring aspects of the biology and ecology of the parasite. Here a model of the dynamics of an infection within a host is presented and used to address questions about host infectivity, the role of preferential invasion in suppressing gametocyte production, and hypnozoite-mediated relapse. While such a model is ultimately an abstraction of a real infection, useful insights can be drawn from the outputs. Future work will allow the development of community-scale, multiple-host models, where parasite infections can be transmitted between hosts by mosquito vectors; such models will facilitate investigations into the suitability and effectiveness of interventions to control disease outbreaks and elimination/eradication strategies for P. vivax.
This work was supported by NHMRC grant #613648.
- Mendis K, Sina BJ, Marchesini P, Carter R. The neglected burden of Plasmodium vivax malaria. Am J Trop Med Hyg. 2001;64:97–106.PubMedGoogle Scholar
- Price RN, Tjitra E, Guerra CA, Yeung S, White NJ, Anstey NM. Vivax malaria: Neglected and not benign. Am J Trop Med Hyg. 2007;77:79–87.PubMed CentralPubMedGoogle Scholar
- Gogia A, Kakar A, Byotra SP. Is benign tertian malaria actually benign? Trop Doct. 2012;42:92–3.View ArticlePubMedGoogle Scholar
- Singh H, Parakh A, Basu S, Rath B. Plasmodium vivax malaria: Is it actually benign? J Infect Public Health. 2011;4:91–5.View ArticlePubMedGoogle Scholar
- Trager W, Jensen JB. Human malaria parasites in continuous culture. Science. 1976;193:673–5.View ArticlePubMedGoogle Scholar
- Trager W, Jensen JB. Continuous culture of Plasmodium falciparum: its impact on malaria research. Int J Parasitol. 1997;27:989–1006.View ArticlePubMedGoogle Scholar
- Garnham PCC. Malaria parasites and other haemosporidia. Oxford: Blackwell; 1966.Google Scholar
- White NJ. Malaria. In: Cook GC, Zumla AI, editors. Manson’s Tropical Diseases. 22nd ed. Edinburgh: Saunders Ltd; 2009.Google Scholar
- Kerlin D, Gatton M. Preferential invasion by Plasmodium merozoites and the self-regulation of parasite burden. PLoS One. 2013;6:e57434.View ArticleGoogle Scholar
- Global Malaria Programme WHO. World Malaria Report 2012. Geneva: World Health Organization; 2012.Google Scholar
- Anstey NM, Russell B, Yeo TW, Price RN. The pathophysiology of vivax malaria. Trends Parasitol. 2009;25:220–7.View ArticlePubMedGoogle Scholar
- Bousema T, Drakeley C. Epidemiology and infectivity of Plasmodium falciparum and Plasmodium vivax gametocytes in relation to malaria control and elimination. Clin Microbiol Rev. 2011;24:377–410.View ArticlePubMed CentralPubMedGoogle Scholar
- Bull PC, Berriman M, Kyes S, Quail MA, Hall N, Kortok MM, et al. Plasmodium falciparum variant surface antigen expression patterns during malaria. PLoS Pathog. 2005;1:202–13.View ArticleGoogle Scholar
- Gatton ML, Cheng Q. Investigating antigenic variation and other parasite–host interactions in Plasmodium falciparum infections in naïve hosts. Parasitology. 2004;128:367–76.View ArticlePubMedGoogle Scholar
- Mueller I, Galinski MR, Baird JK, Carlton JM, Kochar DK, Alonso PL, et al. Key gaps in the knowledge of Plasmodium vivax, a neglected human malaria parasite. Lancet Infect Dis. 2009;9:555–66.View ArticlePubMedGoogle Scholar
- Cogswell FB. The hypnozoite and relapse in primate malaria. Clin Microbiol Rev. 1992;5:26–35.PubMed CentralPubMedGoogle Scholar
- Krotoski WA. Discovery of the hypnozoite and a new theory of malarial relapse. Trans R Soc Trop Med Hyg. 1985;79:1–11.View ArticlePubMedGoogle Scholar
- Huldén L, Huldén L, Heliovaara K. Natural relapses in vivax malaria induced by Anopheles mosquitoes. Malar J. 2008;7:64.View ArticlePubMed CentralPubMedGoogle Scholar
- Nishiura H, Lee HW, Cho SH, Lee WG, In TS, Moon SU, et al. Estimates of short- and long-term incubation periods of Plasmodium vivax malaria in the Republic of Korea. Trans R Soc Trop Med Hyg. 2007;101:338–43.View ArticlePubMedGoogle Scholar
- Huldén L, Huldén L. Activation of the hypnozoite: a part of Plasmodium vivax life cycle and survival. Malar J. 2011;10:90.View ArticlePubMed CentralPubMedGoogle Scholar
- McKenzie FE, Jeffery GM, Collins WE. Plasmodium vivax blood-stage dynamics. J Parasitol. 2002;88:521–35.View ArticlePubMed CentralPubMedGoogle Scholar
- Kerlin DH, Boyce K, Marfurt J, Simpson JA, Kenangalem E, Cheng Q, et al. An analytical method for assessing stage-specific drug activity in Plasmodium vivax malaria: Implications for ex vivo drug susceptibility testing. PLoS Negl Trop Dis. 2012;6:e1772.View ArticlePubMed CentralPubMedGoogle Scholar
- Gatton ML, Cheng Q. Interrupting malaria transmission: quantifying the impact of interventions in regions of low to moderate transmission. PLoS One. 2010;5:e15149.View ArticlePubMed CentralPubMedGoogle Scholar
- Karyana M, Burdarm L, Yeung S, Kenangalem E, Wariker N, Maristela R, et al. Malaria morbidity in Papua Indonesia, an area with multidrug resistant Plasmodium vivax and Plasmodium falciparum. Malar J. 2008;7:148.View ArticlePubMed CentralPubMedGoogle Scholar
- Lin E, Kiniboro B, Gray L, Dobbie S, Robinson L, Laumaea A, et al. Differential patterns of infection and disease with P falciparum and P vivax in young Papua New Guinean children. PLoS One. 2010;5:e9047.View ArticlePubMed CentralPubMedGoogle Scholar
- Luxemburger C, Ricci F, Nosten F, Raimond D, Bathet S, White NJ. The epidemiology of severe malaria in an area of low transmission in Thailand. Trans R Soc Trop Med Hyg. 1997;91:256–62.View ArticlePubMedGoogle Scholar
- Gatton ML, Cheng Q. Evaluation of the pyrogenic threshold for Plasmodium falciparum malaria in naive individuals. Am J Trop Med Hyg. 2002;66:467–73.PubMedGoogle Scholar
- Carter R, Graves PM. Gametocytes. In: Wernssdorfer WH, McGregor I, editors. Malaria Principles and Practice of Malariology, vol. 1. Edinburgh: Churchill Livingstone; 1988. p. 253–305.Google Scholar
- Kwiatkowski D. Febrile temperatures can synchronize the growth of Plasmodium falciparum in vitro. J Exp Med. 1989;169:357–61.View ArticlePubMedGoogle Scholar
- McClean CM, Alvarado HG, Neyra V, Llanos-Cuentas A, Vinetz JM. Optimized in vitro production of Plasmodium vivax ookinetes. Am J Trop Med Hyg. 2010;83:1183–6.View ArticlePubMed CentralPubMedGoogle Scholar
- Mason DP, McKenzie FE. Blood-stage dynamics and clinical implications of mixed Plasmodium vivax-Plasmodium falciparum infections. Am J Trop Med Hyg. 1999;61:367–74.PubMed CentralPubMedGoogle Scholar
- Boyd MF, Kitchen SF. On the infectiousness of patients infected with Plasmodium vivax and Plasmodium falciparum. Am J Trop Med Hyg. 1937;s1-17:253–62.Google Scholar
- Coatney G, Collins W, Warren M, Contacos P. The Primate Malarias. Bethesda MD: Department of Health, Education and Welfare; 1971.Google Scholar
- Mackerras MJ, Ercole QN. Some observations on the action of quinine, atebrin, and plasmoquine on Plasmodium vivax. Trans R Soc Trop Med Hyg. 1949;42:443–54.View ArticlePubMedGoogle Scholar
- Jeffery GM. Infectivity to mosquitoes of Plasmodium vivax following treatment with chloroquine and other antimalarials. Am J Trop Med Hyg. 1958;7:207–11.PubMedGoogle Scholar
- Gatton ML, Cheng Q. Plasmodium falciparum infection dynamics and transmission potential following treatment with sulfadoxine-pyrimethamine. J Antimicrob Chemother. 2006;58:47–51.View ArticlePubMed CentralPubMedGoogle Scholar
- Pichon G, Awono-Amebe HP, Robert V. High heterogeneity in the number of Plasmodium falciparum gametocytes in the bloodmeal of mosquitoes fed on the same host. Parasitology. 2000;121:115–20.View ArticlePubMedGoogle Scholar
- Baird JK, Leksana B, Masbar S, Fryauff DJ, Sutanihardja MA, Suradi, et al. Diagnosis of resistance to chloroquine by Plasmodium vivax: timing of recurrence and whole blood chloroquine levels. Am J Trop Med Hyg. 1997;56:621–6.PubMedGoogle Scholar
- Battle KE, Karhunen MS, Bhatt S, Gething PW, Howes RE, Golding N, et al. Geographical variation in Plasmodium vivax relapse. Malar J. 2014;13:144.View ArticlePubMed CentralPubMedGoogle Scholar
- Ehrman FC, Ellis JM, Young MD. Plasmodium vivax Chesson strain. Science. 1945;101:377.View ArticlePubMedGoogle Scholar
- R Development Core Team. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2011.Google Scholar
- Chen NH, Auliff A, Rieckmann K, Gatton M, Cheng Q. Relapses of Plasmodium vivax infection result from clonal hypnozoites activated at predetermined intervals. J Infect Dis. 2007;195:934–41.View ArticlePubMedGoogle Scholar
- Ishikawa H, Ishii A, Nagai N, Ohmae H, Harada M, Suguri S, et al. A mathematical model for the transmission of Plasmodium vivax malaria. Parasitol Int. 2003;52:81–93.View ArticlePubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.