Skip to main content

Calibration of an intrahost malaria model and parameter ensemble evaluation of a pre-erythrocytic vaccine



A pre-erythrocytic vaccine could provide a useful tool for burden reduction and eventual eradication of malaria. Mathematical malaria models provide a mechanism for evaluating the effective burden reduction across a range of transmission conditions where such a vaccine might be deployed.


The EMOD model is an individual-based model of malaria transmission dynamics, including vector lifecycles and species-specific behaviour, coupled to a mechanistic intrahost model of malaria parasite and host immune system dynamics. The present work describes the extension of the EMOD model to include diagnoses of severe malaria and iterative calibration of the immune system parameters and parasite antigenic variation to age-stratified prevalence, incidence and severe disease incidence data obtained from multiple regions with broadly varying transmission conditions in Africa. An ensemble of calibrated model parameter sets is then employed to evaluate the potential impact of routine immunization with a pre-erythrocytic vaccine.


The reduction in severe malaria burden exhibits a broad peak at moderate transmission conditions. Under sufficiently intense transmission, a vaccine that reduces but does not eliminate the probability of acquisition from a single challenge bite may delay infections but produces minimal or no net reduction. Conversely, under sufficiently weak transmission conditions, a vaccine can provide a high fractional reduction but avert a relatively low absolute number of cases due to low baseline burden.


Roll-out of routine immunization with pre-erythrocytic malaria vaccines can provide substantial burden reduction across a range of transmission conditions typical to many regions in Africa.


Recent years have seen encouraging progress in the reduction of global malaria burden, with the WHO’s global mortality estimates declining roughly 30% from 2000 to 2013 [1]. However, malaria continues to kill hundreds of thousands a year, with the burden of mortality falling most heavily on children in sub-Saharan Africa. Improvements in vector control, diagnostics and treatment have all contributed to the declining mortality [24], and a preventative vaccine could provide an important additional tool to the malaria control and eradication effort. Mathematical models of malaria transmission and within-host disease progression can aid in understanding the effects of potential vaccine distributions under a broad range of transmission conditions. A variety of potential vaccines are under investigation; depending on which stage of the parasite lifecycle is targeted, a vaccine could act to prevent human acquisition, reduce morbidity post-acquisition, or prevent transmission from human hosts back to the mosquito vectors.

The EMOD model is an individual-based, stochastic simulation framework that couples a mechanistic model of intra-host parasite lifecycle and immune system response dynamics to a detailed model of the mosquito population lifecycle, including responses to changes in climate and descriptions of species-dependent behaviours [59]. Campaign interventions targeting the vector population (insecticide-treated nets, habitat clearance) or the human population (vaccines, anti-malarial drug regimens, etc.) can be flexibly distributed to specified subpopulations, e.g., age-based immunization schedules, calendar-based distribution campaigns or treatment-seeking upon onset of symptoms.

Malaria control efforts rely on multiple tools to reduce disease burden: vector control techniques, chemoprophylaxis, and post-exposure drug treatment [10]. The proliferation of drug-resistant parasites [11] or insecticide-resistant mosquitos [12] remains a concern for these techniques, and malaria vaccines could provide a crucial tool in burden reduction and potential elimination of malaria. Pre-erythrocytic vaccines target the Plasmodium falciparum parasite prior to blood-stage disease, aiming to block liver-stage infection from occurring or from progressing to blood-stage infection, thus preventing the onset of symptomatic malaria as well as onward transmission of gametocytes. Among current pre-erythrocytic vaccine candidates, the RTS,S vaccine is in the most advanced stage of clinical trials [1315]. Human trials of the vaccine (with various adjuvants) have been conducted since the early 1990s, demonstrating safety and immunogenicity in clinical and field settings [16]. A recent Phase 3 trial across eleven African sites [13] reports promising results; a 3-dose infant immunization schedule provides efficacies against clinical malaria of 46% over an 18 month follow-up period, with similarly strong efficacies against severe malaria and malaria hospitalization (34% and 41%, respectively).

An effective pre-erythrocytic vaccine could provide a useful tool for burden reduction and eventual eradication of malaria. However, a vaccine that provides limited protection against infection could delay the development of adaptive immunity to blood-stage disease, increasing the malaria burden in older children whose vaccine-derived immunity has waned. This “rebound effect” has been observed to occur in field studies of intermittent preventive treatment and chemoprophylaxis [17, 18], and should be investigated in the case of pre-erythrocytic vaccines. Evaluation of the effective burden reduction across the broad range of transmission conditions where such a vaccine could be deployed is a natural target for mathematical malaria models [19, 20].


The EMOD malaria model has been described in detail in previous papers [59]. The intrahost model features a mechanistic description of within-host infection dynamics that lead to the acquisition of parasitological and clinical immunities with repeated exposure. This model includes descriptions of both the innate and adaptive human immune responses to antigens presented by merozoite surface proteins (MSP), unique epitopes of P. falciparum erythrocyte membrane proteins (PfEMP), and cross-reactive PfEMP epitopes. Each of these antigenic compartments is also allowed a configurable number of antigenic variants [8, 2123]. The adaptive immune system response to each variant within each antigenic compartment is described by a growth rate, capacity, decay rate, memory level, and parasite clearance efficacy. The unique PfEMP epitopes are also characterized by a switching rate, allowing the presentation of distinct variants over the course of a single infection. Individual body temperatures and red blood cell (RBC) concentrations are tracked throughout the simulation to identify clinical and severe disease incidence (defined in Calibration subsection). Body temperature is governed by innate cytokine production, which is stimulated by parasite density (through a pyrogenic threshold) and schizont rupture events, and mitigated by the production of adaptive antibodies [6]. RBC concentrations are governed by processes of destruction through schizont rupture and production by the host.


Prior to evaluating the efficacy of a pre-erythrocytic vaccine in field conditions, the model was calibrated in two steps. First, age-stratified prevalence and fever incidence curves were used to calibrate parameters governing antigenic dynamics and immune response, and then age-stratified severe disease incidence curves were used to calibrate model parameters governing the probabilistic presentation of a severe malarial incident.

The parameters of the immune system model that are most relevant to acquisition of long-term clinical and parasitological immunity were calibrated to age-stratified prevalence and incidence data from six sites experiencing measured entomological inoculation rates (EIRs) ranging from 18 to over 300. Prevalence data comes from Matsari, Sugungum and Rafin Marke in Nigeria [24], and Namawala, Tanzania [25]; incidence data was obtained from Ndiop and Dielmo, Senegal [2628].

Clinical incidence data present some challenges for model calibration. In a model of malaria transmission with no other pathogens, as is the case here, all fevers are caused by malaria. This is not the case in field studies, and the attribution of clinical symptoms can be problematic, particularly at high prevalence. Studies often use cutoffs on the parasitaemia level concurrent with fever incidence to determine a ‘malaria attributable fraction’ , but this is subject to heterogeneities based on local transmission intensity [29]. Other works have also noted that higher frequencies of active case detection visits tend to produce higher estimates of incidence [30, 31]. It is not obvious whether this correlation arises because of over-attribution of non-malarial fevers caught in frequent visits, or to true malaria-attributable fevers being missed in between visits. At these Senegal sites, patients were visited three times per week, and the present work makes no attempt to ‘correct’ the data for case detection frequency or to censor model-predicted fevers from being attributed to malaria if concurrent parasitaemia is low. Such corrections may become necessary when expanding the set of calibration sites or comparing predicted incidence rates against out-of-sample data.

At the Senegal study sites, patients were also treated with quinine when fever presented with concurrent parasitaemia. The high rate of case detection combined with rapid treatment could exert strong effects on the development of immunity. To account for this, the simulated population undergoes a burn-in of forty years, during which the population is exposed to the expected transmission conditions but does not undergo rapid case detection and treatment; this setup prevents high treatment rates from limiting age-dependent immunity development in the population. After the burn-in period, fever incidence reporting begins, and individuals who present fever begin parasite-clearing treatment with a 30% probability per day from the fever onset, meant to roughly capture the ‘time to next visit’ for fever detection. The rapid parasite clearance after fever onset also likely prevents the dataset from informing the model calibration with respect to the effects of coinfection on morbidity.

The simulated transmission conditions at each site are characterized by a total EIR and a normalized seasonal profile. Simulations with local climate data were necessary to determine the seasonal profile of certain sites, while others were taken from literature. Throughout this study, the annual EIR characterizes the biting rate experienced by grown adults; children experience a reduction in this biting rate described by a piecewise-linear function that approximates their expected surface area as a function of age (i.e., biting rates are assumed to be proportional to skin surface area). Age-dependent biting is well documented [32, 33] and has been used in previous modelling work [34].

Because of the interactions between the human and vector systems in the EMOD model, broad variations in human immune dynamics induce variations in the transmission conditions. This presents a difficulty in the calibration process, as the target distributions are age-specific prevalence and incidence curves determined at measured EIRs. It is computationally expensive to adjust vector habitats at each tested set of intrahost parameters to compare results at equal transmission levels. To remove this complication, the measured transmission conditions were reproduced by removing the vectors from the simulation entirely, and instead subjecting the simulated individuals to periodic sporozoite challenge at a known intensity. The frequency of challenge is varied monthly to reproduce seasonal variation in EIR.

The model parameters included in the calibration are the number of antigenic variants in each of the three immune compartments, the killing strengths of the MSP antigens and shared minor epitopes, and the PfEMP antigenic switching rate; other parameters of the model had been calibrated to reproduce individual time courses of infection and distributions of infection durations and were held fixed during this calibration.

For each set of model input parameters, a set of simulations are run to characterize outputs under the measured transmission conditions at each of the field sites from the data, and the overall likelihood of an input parameter set is the product of likelihoods over the six sites. Simulated prevalence output is modelled as random binomial slide positivity measurements, characterized by a number of slides viewed and an individual’s current blood parasite density. The definition of a clinical malarial incident is configurable in the EMOD model; for this study, a clinical incident begins when an individual’s body temperature is raised by 1.5°C. Parasite density alone does not trigger a clinical incident, though temperature and parasite density are implicitly linked through the innate immune response. The clinical incident continues until the fever remains below 0.5°C for two weeks; this refractory period prevents the multiple recrudescent fever events typical to malaria from being recorded as multiple independent clinical incidents.

Incremental mixture importance sampling [35, 36] was employed for parameter space exploration and calibration. Once the target posterior likelihood has been sufficiently well-sampled, this stage of calibration is finished by producing an acceptable region of parameter space, defined by placing a threshold on the overall log-likelihood (based on Wilks’ theorem) and forming the convex hull of all parameter sets above this threshold. This calibrated parameter volume is used as an acceptance region for resampling model parameter sets. More details of the simulation setup, likelihood functions and algorithm can be found in Additional file 1.

The results of the model calibration to age-stratified measurements of parasite prevalence and clinical incidence are presented in Figure 1. The two left columns present the measured fraction of children testing positive for parasites by slide microscopy at four sites with annual EIRs ranging from 18–329. The data are presented as black points with statistical error bars, and the blue-shaded regions that enclose the results of all simulations passing the total log-likelihood threshold (100/95/68% quantiles are presented in progressively darker shading). The right column similarly presents the rate of fever incidence at two sites, with annual EIRs of 20 and 200.

Figure 1
figure 1

(Two left columns) Age-stratified prevalence data and simulation results from four regions in Nigeria and Tanzania. (Right column) Age-stratified incidence data and simulation results from two regions in Senegal. In both panels, the blue-shaded regions represent the 68/95/100% quantiles of simulations passing the likelihood threshold.

Sites with increasing annual EIR exhibit a more rapid rise in detectable parasitaemia levels at the youngest ages and a higher overall peak prevalence, but all sites exhibit a decrease in prevalence over the second decade of life. Above 20 years of age, the observed prevalence does not vary dramatically over this order-of-magnitude range of transmission intensities. In the two Senegal sites, an order of magnitude increase in EIR produces a substantially higher peak rate of clinical incidence, an earlier peak age (2–3 in Dielmo vs 5–10 in Ndiop), and a lower incidence rate in adults (approximately 0.6/year in Ndiop vs 0.2/year in Dielmo). The more intense transmission in Dielmo confers sufficient immune protection that the cumulative lifetime incidence is lower than in Ndiop [26].

The effects of each of the calibration parameters on the features described above are not simple to describe, as there are strong interactions among the terms. Broadly, increasing the number of antigenic variants in the three immune compartments delays the onset of parasitological and symptomatic immunity, while the antibody killing strengths to MSP and cross-reactive PfEMP epitopes most strongly affect the magnitude of the prevalence and incidence in early ages. The antigenic switching rate also strongly affects the onset of symptomatic immunity; when the switching rate is too slow, immunity is acquired slowly as the host is presented with few variants in any individual infection. The antigenic switching rate has also been constrained by previous work on the duration of infection and infectiousness in naïve individuals [8].

The calibration procedure continues by sampling parameter values from the acceptable regions defined in the first stage and varying a new set of parameters governing severe disease, targeting data on age-specific severe disease rates at five sites in The Gambia and Kenya [37], and proportional rates of cerebral or anemic malaria from seventeen sites across Africa [38]. In the EMOD model, an individual’s current state is mapped onto three probabilities of diagnosis of severe malaria due to three underlying causes. An individual’s current RBC count sets a probability of diagnosis due to anaemia and associated presentations, their current body temperature acts as a proxy for presentation of severe cerebral malaria, and their current parasitaemia level acts as a catch-all proxy for other complications (e.g., respiratory involvement). Similar protections to those described above regarding clinical incidents are in place to prevent a single severe presentation from being recorded as multiple incidents.

The functions translating RBC counts, asexual parasite densities and fever levels to probabilities of diagnosis are sigmoid (logistic) functions characterized by a width and a midpoint. The functional form is:

Where k represents the inverse width, and x 0 is the threshold, at which the probability is equal to 0.5. The calibration includes the widths and midpoints associated with severe cerebral malaria or severe parasitaemia, a multiplier governing how many red blood cells are destroyed when a schizont ruptures, and the level of maternal antibody protection (which varies from site to site based on measured transmission conditions – see “Details of calibration setup” in Additional file 1 for more details). The width and threshold for severe anaemia is not included in the calibration, because this diagnosis in the source data is based on a well-defined threshold at haemoglobin density of 5 g/dl. The frequency of anaemia is governed by the competing processes of erythrocyte destruction by rupturing schizonts and increased erythropoiesis in response to decreased haemoglobin density. Because these two effects produce opposing results on the rate of severe anaemia, the parameter governing erythropoiesis was fixed to a value hand-calibrated to measured individual time courses in naïve patients, and only the erythrocyte destruction factor is varied in this calibration.

The results of calibrating the model to age-stratified severe disease diagnoses are presented in Figure 2. Age-stratified total severe disease incidence rates, measured from hospital records and census data, at five sites with parasite prevalence rates in one to nine years olds ranging from 2 to 83% (mapped through simulation to annual EIRs ranging from 0.1 to 55) are presented on the left. The proportions of severe disease incidence in the under 24 months and five to nine years age bins, and the proportions of cerebral and anemic disease at 17 sites, along with simulation results, are presented on the right. The blue-shaded regions are defined as in Figure 1.

Figure 2
figure 2

(Left) Age-stratified severe disease incidence data and results from simulations from five regions in Africa; note the changing y-axis when comparing. (Right) Proportions of severe disease incidence in various categories. In both panels, the blue-shaded regions represent the 68/95/100% quantiles of simulations passing the likelihood threshold. PfPR2–10 is the proportion of 2–10 years olds testing positive for falciparum parasites.

In the parameter region favored post-calibration, the probability of experiencing a severe cerebral malaria episode is essentially zero for low-grade fevers and increases rapidly over a range of body temperatures from 40-41°C (base body temperature is assumed to be 37°C). Similarly, a severe disease episode caused by hyperparasitaemia is highly unlikely at parasite densities below approximately 200,000 per μl, rapidly increasing to a probability approaching 1 by approximately 400,000 per μl (figures available in Additional file 1).

Sites with increasing transmission intensity exhibit an earlier peak age in severe disease rate, and a corresponding increase in the proportion of severe disease due to anemia. Maternal antibodies play a large role in the initial rise in severe incidence. The degree of protection conferred by maternal antibodies was modelled as an increasing function of transmission intensity, but a mother’s antibody levels against particular antigenic variants were not considered in these simulations. The cumulative severe disease incidence is highest at moderate EIRs. The model predicts that under intense transmission conditions, many children are exposed to the parasite while maternal antibodies still provide sufficient protection to prevent severe symptoms, and these early exposures confer a degree of protection against severe malaria in future infections. This idea has also been hypothesized to explain observations that cumulative paediatric severe malaria incidence is highest under moderate conditions and declines at extremely high levels of transmission [37, 39]. The other five calibration parameters interact strongly in setting the overall rate of severe disease incidence as well as the proportions of severe disease ascribed to the three underlying mechanisms.

The limitations of this calibration procedure are worth addressing here. The prevalence and incidence calibration was performed only to sites experiencing annual EIRs > 18, which represents substantially higher transmission intensity than currently experienced by most of the African population [40]. To provide confidence that the calibrated model can reasonably be applied to study lower transmission intensities, the results of this calibration procedure were checked against data outside of the calibration sample down to annual EIRs of 1. Details and figures are contained in Additional file 1. Individuals in the model are also homogenous in terms of immune system dynamics, response to treatment, and biting risk. Individual- and site-level heterogeneities in these quantities and others have been observed in studies [14, 4145] and can exhibit significant effects on population-level data. These effects are under consideration in extensions to the present work.

Parameter ensemble model vaccine evaluation

After calibration, the effects of a potential rollout of a pre-erythrocytic vaccine were evaluated across a wide range of transmission conditions and for an ensemble of acceptable input parameters. Each parameter set in the ensemble represents an independent set of input parameters sampled from the 12-dimensional calibrated parameter space. One-hundred and seventy parameter sets were included in the ensemble; this number resulted from available computation time, rather than any considerations specific to the analysis. Each of the 170 models was run under a variety of transmission conditions; these were not explicitly based on particular locations, but on a set of four seasonal profiles and 18 total annual EIR magnitudes from 1 to 300. The four seasonal profiles included two realistic profiles based on Sugungum, Nigeria and Mocuba, Mozambique; a non-seasonal profile (constant rate of biting year-round); and a highly seasonal profile in which all biting takes place at a flat rate for six months and no biting occurs for the remaining six months. Eight stochastic repetitions of each model/transmission scenario were performed, beginning at different times of the year to ameliorate potential systematic transient effects. Finally, each setup was run once without vaccine distribution to provide baseline incidence rates, and run twice more with different vaccine efficacy profiles. Simulation outputs are binned into six-month time periods for analysis.

Each individual simulation tracks a cohort of 1,000 individuals from birth until ten years of age. Individuals receive vaccination between six and nine months of age. This cohort-tracking simulation setup necessarily considers only individual-level protection, ignoring any population-level effects of the vaccination. Crosschecks of results from cohort-based and full-population models indicate that vaccinating only newborns with a waning vaccine produces minimal effects on population-level transmission intensity. The pre-erythrocytic vaccine model is characterized by two parameters: an initial degree of protection against infection and a half-life of protection. The degree of protection is defined as the reduction in probability of successful infection in a single challenge bite. Two vaccine profiles were evaluated against baseline – one with an initial degree of protection of 80% and a half-life of nine months (referred to as fast-waning), and another with an initial protection of 60% but a half-life of three years (slow-waning). Figure 3 illustrates the two vaccine profiles evaluated.

Figure 3
figure 3

The two pre-erythrocytic vaccine models evaluated: 80% efficacy with nine-month halflife and 60% efficacy with three-year halflife.


Over essentially the entire range of simulated transmission conditions, the model predicts that a pre-erythrocytic vaccine delivered early in life results in a reduction in the rate of severe malarial incidence for some period of time, followed by a rebound in case counts, when the delayed onset of acquired immunity results in vaccinated children exhibiting a higher rate of severe malaria incidence than unvaccinated children. At all levels of transmission, the median net number of severe cases averted is positive, though it approaches zero at the highest simulated EIR. Figure 4 presents a typical illustrative example of this reduction and rebound behaviour, comparing the median difference in disease incidence in the baseline scenario vs the fast-waning vaccine scenario at an annual EIR of 10. The figure shows a strong reduction in disease incidence over the first two years’ post vaccination (at six months of age), with a rebound effect causing a slightly higher rate of severe disease in the vaccinated group as compared to baseline after two-and-half years of age. The relative magnitudes and timescales of case reduction and rebound are most strongly dependent on the magnitude of the EIR and the vaccine profile, and less so on the model input parameters or the seasonality (though it should be noted that the idea of optimizing the timing of vaccination with respect to local seasonal transmission was not tested). This section will thus focus on case reduction as a function of annual EIR and vaccine profile.

Figure 5 presents the period of case reduction prior to rebound against the annual EIR for both vaccine models. Referring to Figure 4, the period of case reduction is defined as the interpolated crossing of the incidence curves in baseline and vaccine scenarios (marked in black). The eight stochastic repetitions are averaged, and the median/quantiles are computed over the four seasonal profiles and 170 model parameter sets at each annual EIR value. As should be expected, the protective period decreases as the annual EIR increases for both vaccine profiles, and the slowly waning vaccine provides a substantially longer period of case reduction under low transmission intensity. At higher transmission intensities, there is little difference in the period of case reduction between the two vaccine models. Modelling the vaccine as inducing imperfect, probabilistic protection against infection in a single bite implies that vaccination provides minimal effective protection when exposed to highly frequent infectious biting.

The fractional (top) and absolute (bottom) reduction in the cumulative severe (left) and clinical (right) incidence over the first ten years of life are presented against annual EIR for both vaccine profiles in Figure 6. The total severe incidence reduction from distribution of a pre-erythrocytic vaccine exhibits a broad maximum over a range of moderate EIR values from roughly 10 to 40. Under mild transmission conditions, the rate of severe malaria is sufficiently low that, although the vaccine is highly efficacious in terms of relative case reduction, the absolute number of cases averted is quite small. And under intense transmission, transmission is sufficiently strong that a vaccine that acts to randomly block a maximum of 60 or 80% of potential acquisitions will provide only a slight delay in disease acquisition, with insignificant net burden reduction. It should be noted that the uncertainties in the effects of vaccination arise from varying the input parameters within the acceptable region defined in calibration; uncertainties due to model structure (e.g., individual- or site-level heterogeneity) are outside of the scope of the present work.

Figure 4
figure 4

Severe disease incidence per 1,000 children per year at annual EIR = 15. The median baseline incidence is in dark red, and incidence with the fast-waning vaccine is shown in dark blue. The region of incidence reduction at young ages is highlighted in light blue. The disease incidence rebounds at ~2.5 years of age, and the rebound is highlighted in light red.

Figure 5
figure 5

Period of severe (left) and clinical (right) case reduction prior to rebound for the fast-waning and slow-waning vaccine models. The lines follow the median output value from the ensemble of model parameters, and shaded areas outline the 68% quantiles.

Figure 6
figure 6

EIR-dependent reduction in relative (top) and total (bottom) incidence of severe (left) and clinical (right) malaria; incidence, per 1,000 children over the first ten years of life, for both vaccine models. The lines follow the median output value from the ensemble of model parameters, and shaded areas outline the 68% quantiles.

The maximal reduction in cumulative clinical incidence occurs at lower annual EIRs, with a peak around an EIR of 5. At EIRs from 30–100, the rebound effect can actually outweigh the initial case reduction, such that recipients of the vaccine experience slightly increased cumulative rates of clinical incidence. Reductions in transmission intensity have previously been observed to produce increased cumulative clinical incidence, so it is not unexpected that delayed acquisition of protective immunity due to vaccination and increased biting risk at older ages would combine to produce a net increase in cumulative incidence. At the highest EIRs tested, vaccination produces no significant change in total clinical incidence, as the intensity of transmission overwhelms the partially protective vaccination.

The Malaria Atlas Project [40] has produced estimates of transmission conditions on a finely partitioned spatial grid throughout sub-Saharan Africa. These estimates are combined with the results of this parameter ensemble evaluation to produce the visualizations in Figure 7. The left panel maps the estimated EIR values, and the center and right panels represent, for the fast and slow-waning vaccines, respectively, the expected number of cases averted per 1,000 children over a ten-year horizon. The resulting maps project where the rollout of a pre-erythrocytic vaccine would produce optimal benefits. The magnitude of the total cases averted assumes 100% vaccination coverage of newborns, and does not consider coordinated rollout of the vaccine with other transmission-reducing interventions.

Figure 7
figure 7

(left) Malaria Atlas Project estimates of transmission intensity (EIR values interpolated from PfPR 2–10 ). Expected cases averted per 1,000 children vaccinated over a ten-year horizon, 100% infant immunization coverage with (center) fast-waning vaccine and (right) slow-waning vaccine.


Routine immunization of infants has proven to be an incredibly effective tool for reducing the burden of numerous infectious diseases, particularly so for diseases in which protection is long-lasting and naïve individuals dominate the transmission chain. The present work finds that routine immunization with a vaccine providing imperfect protection against acquisition of malaria would provide optimal burden reduction within a range of EIRs from approximately 10–40, where there is sufficient burden for a vaccine to produce large effects, but the frequency of parasite inoculation is not so high as to overwhelm the vaccine’s protective effects. Simulations indicate that, due to malaria’s ability to re-infect human hosts numerous times, infant immunization alone does not strongly perturb population-level transmission dynamics, and efforts towards local elimination must include broader intervention strategies. This evaluation was performed using an ensemble of model parametrizations, and it is found that the results are robust within a broad volume of model parameter space consistent with existing population-level measurements.

Extensions to this work may consider the effects of routine boosters as the vaccine’s protective effects wane, mass vaccinations at sufficient scales to affect malaria’s transmission chain, or correlated distribution of other anti-malarial interventions in regions where the vaccine is deployed. The presented models also do not consider potential differences in overall mortality resulting from shifting a fraction of the severe disease burden onto older children. The calibration framework developed also allows for straightforward addition of new data targets, including similar measurements taken in other regions or different types of measurement (e.g., measured parasite densities), and further work along this direction is also merited.



Entomological inoculation rate

PfPR2–10 :

Proportion of two- to ten-years-olds testing positive for Plasmodium falciparum


P. falciparum erythrocyte membrane protein


Merozoite surface proteins


Red blood cells.


  1. WHO: World Malaria Report: 2013. Geneva: World Health Organization; 2013.

    Google Scholar 

  2. Akachi Y, Atun R: Effect of investment in malaria control on child mortality in sub-Saharan Africa in 2002–2008. PLoS ONE 2011, 6:e21309. 10.1371/journal.pone.0021309

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Lim SS, Fullman N, Stokes A, Ravishankar N, Masiye F, Murray CJL, Gakidou E: Net benefits: a multicountry analysis of observational data examining associations between insecticide-treated mosquito nets and health outcomes. PLoS Med 2011, 8:e1001091. 10.1371/journal.pmed.1001091

    Article  PubMed Central  PubMed  Google Scholar 

  4. Snow RW, Guerra CA, Noor AM, Myint HY, Hay SI: The global distribution of clinical episodes of Plasmodium falciparum malaria. Nature 2005, 434:214–217. 10.1038/nature03342

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Eckhoff PA: A malaria transmission-directed model of mosquito life cycle and ecology. Malar J 2011, 10:303. 10.1186/1475-2875-10-303

    Article  PubMed Central  PubMed  Google Scholar 

  6. Eckhoff PA: Malaria parasite diversity and transmission intensity affect development of parasitological immunity in a mathematical model. Malar J 2012, 11:419. 10.1186/1475-2875-11-419

    Article  PubMed Central  PubMed  Google Scholar 

  7. Eckhoff P: Mathematical models of within-host and transmission dynamics to determine effects of malaria interventions in a variety of transmission settings. Am J Trop Med Hyg 2013, 88:817–827. 10.4269/ajtmh.12-0007

    Article  PubMed Central  PubMed  Google Scholar 

  8. Eckhoff P: P. falciparum infection durations and infectiousness are shaped by antigenic variation and innate and adaptive host immunity in a mathematical model. PLoS ONE 2012, 7:e44950. 10.1371/journal.pone.0044950

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Wenger EA, Eckhoff PA: A mathematical model of the impact of present and future malaria vaccines. Malar J 2013, 12:126. 10.1186/1475-2875-12-126

    Article  PubMed Central  PubMed  Google Scholar 

  10. Girard MP, Reed ZH, Friede M, Kieny MP: A review of human vaccine research and development: malaria. Vaccine 2007, 25:1567–1580. 10.1016/j.vaccine.2006.09.074

    Article  CAS  PubMed  Google Scholar 

  11. White NJ: Antimalarial drug resistance. J Clin Invest 2004, 113:1084–1092. 10.1172/JCI21682

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  12. Hemingway J, Ranson H: Insecticide resistance in insect vectors of human disease. Annu Rev Entomol 2000, 45:371–391. 10.1146/annurev.ento.45.1.371

    Article  CAS  PubMed  Google Scholar 

  13. RTS,S Clinical Trials Partnership: Efficacy and safety of the RTS,S/AS01 malaria vaccine during 18 months after vaccination: a phase 3 randomized, controlled trial in children and young infants at 11 African sites. PLoS Med 2014, 11:e1001685. 10.1371/journal.pmed.1001685

    Article  Google Scholar 

  14. Olotu AI Ph.D. Thesis. In Long Term Efficacy of a Pre-Erythrocytic Malaria Vaccine and Correlates of Protection in Children Residing in a Malaria Endemic Country. University of Oxford; 2013.

    Google Scholar 

  15. White MT, Bejon P, Olotu A, Griffin JT, Bojang K, Lusingu J, Salim N, Abdulla S, Otsyula N, Agnandji ST, Lell B, Asante KP, Owusu-Agyei S, Mahama E, Agbenyega T, Ansong D, Sacarlal J, Aponte JJ, Ghani AC: A combined analysis of immunogenicity, antibody kinetics and vaccine efficacy from phase 2 trials of the RTS,S malaria vaccine. BMC Med 2014, 12:117.

    Article  PubMed Central  PubMed  Google Scholar 

  16. Regules JA, Cummings JF, Ockenhouse CF: The RTS, S vaccine candidate for malaria. Expert Rev Vaccines 2011, 10:589–599. 10.1586/erv.11.57

    Article  CAS  PubMed  Google Scholar 

  17. Konaté AT, Yaro JB, Ouédraogo AZ, Diarra A, Gansané A, Soulama I, Kangoyé DT, Kaboré Y, Ouédraogo E, Ouédraogo A, Tiono AB, Ouédraogo IN, Chandramohan D, Cousens S, Milligan PJ, Sirima SB, Greenwood BM, Diallo DA: Morbidity from malaria in children in the year after they had received intermittent preventive treatment of malaria: a randomised trial. PLoS ONE 2011, 6:e23391. 10.1371/journal.pone.0023391

    Article  PubMed Central  PubMed  Google Scholar 

  18. Aponte JJ, Menendez C, Schellenberg D, Kahigwa E, Mshinda H, Vountasou P, Tanner M, Alonso PL: Age interactions in the development of naturally acquired immunity to Plasmodium falciparum and its clinical presentation. PLoS Med 2007, 4:e242. 10.1371/journal.pmed.0040242

    Article  PubMed Central  PubMed  Google Scholar 

  19. Maire N, Tediosi F, Ross A, Smith T: Predictions of the epidemiologic impact of introducing a pre-erythrocytic vaccine into the expanded program on immunization in sub-Saharan Africa. Am J Trop Med Hyg 2006,75(2 suppl):111–118.

    PubMed  Google Scholar 

  20. Smith T, Killeen GF, Maire N, Ross A, Molineaux L, Tediosi F, Hutton G, Utzinger J, Dietz K, Tanner M: Mathematical modeling of the impact of malaria vaccines on the clinical epidemiology and natural history of Plasmodium falciparum malaria: overview. Am J Trop Med Hyg 2006,75(2 suppl):1–10.

    PubMed  Google Scholar 

  21. Horrocks P, Kyes SA, Bull PC, Deitsch KW: Molecular Aspects of Antigenic Variation in Plasmodium falciparum . In Mol Approaches Malar Edited by: Sherman IW. 2005, 399–415.

    Google Scholar 

  22. Recker M, Buckee CO, Serazin A, Kyes S, Pinches R, Christodoulou Z, Springer AL, Gupta S, Newbold CI: Antigenic variation in Plasmodium falciparum malaria involves a highly structured switching pattern. PLoS Pathog 2011, 7:e1001306. 10.1371/journal.ppat.1001306

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  23. Holder AA, Guevara Patiño JA, Uthaipibull C, Syed SE, Ling IT, Scott-Finnigan T, Blackman MJ: Merozoite surface protein 1, immune evasion, and vaccines against asexual blood stage malaria. Parassitologia 1999, 41:409–414.

    CAS  PubMed  Google Scholar 

  24. Molineaux L, Gramiccia G: The Garki Project : Research on the Epidemiology and Control of Malaria in the Sudan Savanna of West Africa. Geneva: World Health Organization; 1980:311.

    Google Scholar 

  25. Smith T, Charlwood JD, Kihonda J, Mwankusye S, Billingsley P, Meuwissen J, Lyimo E, Takken W, Teuscher T, Tanner M: Absence of seasonal variation in malaria parasitaemia in an area of intense seasonal transmission. Acta Trop 1993, 54:55–72. 10.1016/0001-706X(93)90068-M

    Article  CAS  PubMed  Google Scholar 

  26. Rogier C, Tall A, Diagne N, Fontenille D, Spiegel A, Trape JF: Plasmodium falciparum clinical malaria: lessons from longitudinal studies in Senegal. Parassitologia 1999, 41:255–259.

    CAS  PubMed  Google Scholar 

  27. Trape JF, Rogier C, Konate L, Diagne N, Bouganali H, Canque B, Legros F, Badji A, Ndiaye G, Ndiaye P: The Dielmo project: a longitudinal study of natural malaria infection and the mechanisms of protective immunity in a community living in a holoendemic area of Senegal. Am J Trop Med Hyg 1994, 51:123–137.

    CAS  PubMed  Google Scholar 

  28. Rogier C, Trape J-F: Etude de l’acquisition de la prémunition en zones d’holo- et de méso-endémie palustre et Dielmo et a Ndiop (Sénégal): Résultats préliminaires, 1990–1994. Méd Trop (Marseille) 1995, 55:71–76.

    CAS  Google Scholar 

  29. Smith T, Schellenberg JA, Hayes R: Attributable fraction estimates and case definitions for malaria in endemic areas. Stat Med 1994, 13:2345–2358. 10.1002/sim.4780132206

    Article  CAS  PubMed  Google Scholar 

  30. Snow RW, Craig M, Deichmann U, Marsh K: Estimating mortality, morbidity and disability due to malaria among Africa’s non-pregnant population. Bull World Health Organ 1999, 77:624–640.

    PubMed Central  CAS  PubMed  Google Scholar 

  31. Patil AP, Okiro EA, Gething PW, Guerra CA, Sharma SK, Snow RW, Hay SI: Defining the relationship between Plasmodium falciparum parasite rate and clinical disease: statistical models for disease burden estimation. Malar J 2009, 8:186. 10.1186/1475-2875-8-186

    Article  PubMed Central  PubMed  Google Scholar 

  32. Carnevale P, Frezil J, Bosseno M, Lepont F, Lancien J: Study of agressivity of Anopheles gambiae in relation to age and sex. Bull World Health Organ 1978, 56:147–154.

    PubMed Central  CAS  PubMed  Google Scholar 

  33. Port GR, Boreham PFL, Bryan JH: The relationship of host size to feeding by mosquitoes of the Anopheles gambiae Giles complex (Diptera: Culicidae). Bull Entomol Res 1980, 70:133–144. 10.1017/S0007485300009834

    Article  Google Scholar 

  34. Smith T, Maire N, Dietz K, Killeen GF, Vounatsou P, Molineaux L, Tanner M: Relationship between the entomologic inoculation rate and the force of infection for Plasmodium falciparum malaria. Am J Trop Med Hyg 2006,75(2 Suppl):11–18.

    PubMed  Google Scholar 

  35. Raftery AE, Bao L: Estimating and projecting trends in HIV/AIDS generalized epidemics using incremental mixture importance sampling. Biometrics 2010, 66:1162–1173. 10.1111/j.1541-0420.2010.01399.x

    Article  PubMed Central  PubMed  Google Scholar 

  36. Steele RJ, Raftery AE, Emond MJ: Computing normalizing constants for finite mixture models via incremental mixture importance sampling (IMIS. J Comput Graph Stat 2006, 15:712–734. 10.1198/106186006X132358

    Article  Google Scholar 

  37. Snow RW, Omumbo JA, Lowe B, Molyneux CS, Obiero JO, Palmer A, Weber MW, Pinder M, Nahlen B, Obonyo C, Newbold C, Gupta S, Marsh K: Relation between severe malaria morbidity in children and level of Plasmodium falciparum transmission in Africa. Lancet 1997, 349:1650–1654. 10.1016/S0140-6736(97)02038-2

    Article  CAS  PubMed  Google Scholar 

  38. Okiro EA, Al-Taiar A, Reyburn H, Idro R, Berkley JA, Snow RW: Age patterns of severe paediatric malaria and their relationship to Plasmodium falciparum transmission intensity. Malar J 2009, 8:4. 10.1186/1475-2875-8-4

    Article  PubMed Central  PubMed  Google Scholar 

  39. Gupta S, Snow RW, Donnelly CA, Marsh K, Newbold C: Immunity to non-cerebral severe malaria is acquired after one or two infections. Nat Med 1999, 5:340–343. 10.1038/6560

    Article  CAS  PubMed  Google Scholar 

  40. Gething PW, Patil AP, Smith DL, Guerra CA, Elyazar IR, Johnston GL, Tatem AJ, Hay SI: A new world malaria map: Plasmodium falciparum endemicity in 2010. Malar J 2011, 10:378. 10.1186/1475-2875-10-378

    Article  PubMed Central  PubMed  Google Scholar 

  41. Smith DL, Dushoff J, Snow RW, Hay SI: The entomological inoculation rate and Plasmodium falciparum infection in African children. Nature 2005, 438:492–495. 10.1038/nature04024

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  42. Griffin JT, Ferguson NM, Ghani AC: Estimates of the changing age-burden of Plasmodium falciparum malaria disease in sub-Saharan Africa. Nat Commun 2014, 5:3136.

    Article  PubMed Central  PubMed  Google Scholar 

  43. Woolhouse MEJ, Dye C, Etard J-F, Smith T, Charlwood JD, Garnett GP, Hagan P, Hii JLK, Ndhlovu PD, Quinnell RJ, Watts CH, Chandiwana SK, Anderson RM: Heterogeneities in the transmission of infectious agents: implications for the design of control programs. Proc Natl Acad Sci U S A 1997, 94:338–342. 10.1073/pnas.94.1.338

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  44. Walther M, Woodruff J, Edele F, Jeffries D, Tongren JE, King E, Andrews L, Bejon P, Gilbert SC, Souza JBD, Sinden R, Hill AVS, Riley EM: Innate immune responses to human malaria: heterogeneous cytokine responses to blood-stage Plasmodium falciparum correlate with parasitological and clinical outcomes. J Immunol 2006, 177:5736–5745. 10.4049/jimmunol.177.8.5736

    Article  CAS  PubMed  Google Scholar 

  45. White MT, Griffin JT, Drakeley CJ, Ghani AC: Heterogeneity in malaria exposure and vaccine response: implications for the interpretation of vaccine efficacy trials. Malar J 2010, 9:82. 10.1186/1475-2875-9-82

    Article  PubMed Central  PubMed  Google Scholar 

Download references


The authors would like to thank Bill and Melinda Gates for their active support of this work and their sponsorship through the Global Good Fund. The authors also thank Dan Klein at IDM for his support in developing model calibration tools, and Ewan Cameron and Katie Battle for useful discussions.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Kevin A McCarthy.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

KM commissioned the simulations, analysed the results and drafted the paper. EW, GH and PE designed and constructed various features of the malaria model. All authors have read and approved of the final manuscript.

Electronic supplementary material


Additional file 1: Supplemental material: provides further details of calibration setup, results, and validation. (PDF 2 MB)

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

McCarthy, K.A., Wenger, E.A., Huynh, G.H. et al. Calibration of an intrahost malaria model and parameter ensemble evaluation of a pre-erythrocytic vaccine. Malar J 14, 6 (2015).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: