Malaria parasite diversity and transmission intensity affect development of parasitological immunity in a mathematical model
© Eckhoff; licensee BioMed Central Ltd. 2012
Received: 23 September 2012
Accepted: 13 December 2012
Published: 15 December 2012
The development of parasitological immunity against malaria affects the ability to detect infection, the efficiency of the local human parasite reservoir at infecting mosquitoes, and the response to reintroduction of parasites to previously cleared areas. Observations of similar age-trends in detected prevalence and mean parasitaemia across more than an order-of-magnitude of variation in baseline transmission complicate simple exposure-driven explanations.
Mathematical models often employ age-dependent immune factors to match the observed trends, while the present model uses a new detailed mechanistic model of parasite transmission dynamics to explain age-trends through the mechanism of parasite diversity. Illustrative simulations are performed for multiple field sites in Tanzania and Nigeria, and observed age-trends and seasonality in parasite prevalence are recreated in silico, proffering possible mechanistic explanations of the observational data.
Observed temporal dynamics in measured parasitaemia are recreated for each location and age-prevalence outputs are studied. Increasing population-level diversity in malaria surface antigens delays development of broad parasitological immunity. A local parasite population with high diversity can recreate the observed trends in age-prevalence across more than an order of magnitude of variation in transmission intensities.
Mechanistic models of human immunity and parasite antigen diversity can recreate the observed temporal patterns for the development of parasitological immunity across a wide range of transmission intensities. This has implications for the distribution of disease burden across the population, the human transmission reservoir, design of elimination campaigns, and development and roll-out of potential vaccines.
Notable progress has been made against malaria over the past seven years, with the burden of malaria decreasing from an estimated 500 million cases and 1 million deaths in 2005  to an estimated 200 million cases and 650,000 deaths in 2010 . Despite this progress, the burden of malaria remains intolerably high, and current gains may be fragile. Interventions may wear out and require replacement or become less useful due to resistance. Malaria is thus the focus of a global effort aimed at its eradication.
Design of elimination campaigns should not follow a one-size-fits-all pattern due to geographic heterogeneity in malaria epidemiology. High levels of geographic heterogeneity exist in transmission intensity , local vector species , and population densities and housing conditions . Not only is the epidemiology highly heterogeneous geographically, but it is changing rapidly through urbanization in Africa  and the vast scale-up in intervention coverage . Success in the ongoing eradication campaign will not be due to a single intervention, but to combinations of interventions tailored to local circumstances. In some cases, attainment of local elimination may require development of new tools such as vaccines or novel means of vector control. The challenges facing malaria eradication define an ambitious research agenda as outlined by the MalERA framework [7, 8].
Several key questions concern the development of immunity to severe disease, symptomatic disease, and the parasite itself due to previous exposure and infection. Prior exposure can influence the clinical outcome of a new infection , the extent of symptoms , and the probability of transmission [11, 12], all of which have programmatic implications. Different time intervals are observed to be involved in acquisition of immunity to severe malaria, to symptomatic malaria, and to detectable malaria , but sterilizing immunity is not observed in natural conditions. Immunity to severe malaria can develop in a few years, but rates of detected parasitaemia tend to remain high until the age of ten and then decrease slowly towards adult levels beyond the age of twenty . Malaria exhibits tremendous antigenic variation, both within single infections  and across the parasite population , and variant-specific exposure is a strong predictor of future responses . Adaptive immunity is also observed to Plasmodium merozoites , and the innate inflammatory immune response is important for parasite control and patient symptoms . These human immune responses and parasite antigenic diversity interact to create the patterns of immune acquisition, but modelling these dynamics can be difficult.
Mathematical modelling of malaria has a rich history dating back to the work of Ross and Macdonald and their contemporaries, which established key principles for the first global eradication campaign [19–21]. Mathematical modelling of malaria made major leaps forward through the Garki Project and the model built to study it , as well as modern models for population dynamics which followed [23–27]. In addition to population-level models, recent research has advanced the state of detailed models of the within-host dynamics of Plasmodium falciparum[28–30]. Models have been constructed to study antigenic variation at the single infection level [31, 32]. Other work has advanced the state of modelling of vector population dynamics and intensity of vector-borne disease transmission [33–36]. An important new class of models combines population-level transmission with an ever richer representation of individual-level infection and immune dynamics [14, 25, 26, 37, 38].
Mathematical models have been used to understand the development of immunity over time and exposure at the population level [14, 37]. These approaches tend to use combinations of age, integrated exposure, or total infections to set overall immune levels as a proxy for the specific-immune responses acquired over time. The OpenMalaria model [14, 25], for example, implements pre-erythrocytic immunity, maternal antibodies, changes in biting exposure with age, cumulative number of infections, and lifetime integrated parasite density to achieve good fits to age-prevalence data without distinguishing among the specific antigens for each infection. Multiple other models also use combinations of age and antigen-agnostic exposure to recreate the acquisition of immunity [22, 26, 37]. At another end of the modelling spectrum, models based on strain theory  tend to represent acquisition of immunity as the gradual development of quasi-sterilizing immunity to a large repertoire of separate strains with distinct antigenic profiles for each infection. The lack of observed sterilizing immunity to primary infection is one of multiple issues with strain theory, but the observed roles of immune response to specific antigens  and antigen-specific acquisition of immunity  illustrate possible shortcomings of antigen-agnostic models that rely only on age and cumulative exposure.
This work presents a new model which combines previous work on individual-based models that contain detailed local vector population dynamics  with a more detailed representation of within-host parasite dynamics and immune responses . It aims to develop and exercise a detailed model of transmission with a parasite population exhibiting multiple forms of clonal and intraclonal antigenic variation to which individual human agents in the simulation acquire immunity over the course of repeated exposures. Each new simulated infection may have a mix of previously-experienced and novel antigens. High transient immune responses eventually clear ongoing infections when the parasite fails to introduce an antigenic variant which can propagate. After infections are cleared, antibody responses decay to lower levels of specific immunological memory. Immunological memory reduces the parasite density of future infections presenting that antigen and may aid in the eventual clearance of such infections, although sterilizing immunity is not attained.
The present model is used to investigate and illustrate questions about the acquisition of immunity and its dependence on age and exposure. The acquisition of immunity as a function of age and exposure is compared to field data from the Garki Project [22, 44] and other field sites in Namawala, Tanzania . The diversity in parasite antigens across the entire parasite population in a region affects the time for the development of parasitological immunity and the population dynamics. By implementing a repertoire of specific antigens to which non-sterilizing immunity is acquired over time, the present model provides a platform for studying the roles of antigen-specific immunity and diversity in the parasite population without introducing less-substantiated assumptions of strain theory. As such, this platform can be used to explore more mechanistic interpretations of well-established relationships of immunity with age and exposure.
Model structure overview
Parameters and variables
Value in simulations
Time remaining in current asexual cycle
Number of IRBCs of type i
Density of IRBCs of type i per μL
Non-dimensional capacity to mount antibody response to variant f(i)
Non-dimensional antibody response to variant f(i)
Non-dimensional level of inflammatory cytokines
Fever in degrees C
Adaptive immune response to MSP variant j
Number of IRBCs of type i with progeny expressing type i + j
Success probability of merozoite invasion
Number of gametocytes of stage i
Mature female gametocytes per μL
Probability of clearing IRBC of type i
Probability of successfully infecting mosquito
Varies up to 0.3
Multiple Anopheline species can be tracked simultaneously, each with potentially different behaviors and ecologies, and thus possibly responding distinctly to different combinations of interventions. Within a single species population, adult mosquitoes can be represented as individual agents in the simulation or grouped as identical-state cohorts. The number of required cohorts corresponds to the number of occupied unique states for mosquitoes, and more detail and complexity requires more unique states. The individual-agent version of the vector population model scales with the number of local vectors at a given time.
Mosquitoes begin the life cycle as eggs oviposited by existing adult females following successful feeds. Following hatching, cohorts progress through larval and pupal development at temperature-dependent rates to emergence. Progress through this phase is tracked by a cohort-specific variable which increases from zero by a temperature-dependent amount each time step until it reaches completion, and thus this representation allows varying temperatures and varying lengths of larval development in a simple and memory-sparing structure. All mosquitoes in the cohort are assumed to experience the same temperatures on a given time step. Following emergence, there is a brief delay during which sugar feeding and mating occur. Following mating, the adult females begin a cycle of feeding and oviposition. Feeds upon gametocyte-carrying humans carry a chance of infection, and mosquitoes that are successfully infected by the parasite change state to infected and progress at a temperature-dependent rate towards infectiousness. The structure and implementation of the progress variable is similar to the implementation for larval progress, but with a different temperature-dependent rate function. Successful blood meals lead to oviposition, and this closed-loop egg laying can thus capture the local elimination of certain species following intense vector control.
Within-host parasite dynamics
These vector population and transmission models are coupled to microsolvers for the natural history of the parasite within a human host . The vector population and transmission models are simulated with a one-day time step, to correspond to the interval between nights of vectors feeding. The intrahost microsolver uses a combination of discrete and continuous processes to represent parasite dynamics and immune responses . The parasite begins an infection in the human liver and emerges after a latency of seven days. Each infected red blood cell (IRBC) expresses one Plasmodium falciparum Erythrocyte Membrane Protein 1 (PfEMP-1) surface antigen out of 50 for that infection and one minor epitope out of a set of five. The 50 PfEMP-1 antigens are drawn from the overall population of possible antigens n PfEMP1 (varying between 300 and 1,000 in the present results), as is done for each of the sets of five minor epitopes n minor (varying from 20 to 100 total minor epitopes) and for the merozoite surface antigen n MSP (varying from 4 to 100 in the present results). An indexing function f j (i) relates PfEMP-1 antigen i in infection j to the overall list of antigens which extends from 1 to n PfEMP1 . A related indexing function maps which minor epitope in n minor corresponds to the specific variant i. Superinfection is supported in this model, and each concurrent infection will be simulated individually with its own set of indexing functions. These infections all interact with the same immune system at the same time, and the immune microsolver governs the dynamics of responses to all antigens, not just the antigens of a single infection.
Asexual cycle parasite dynamics and immune response
Gametocyte production and transmission
The human immune response within the microsolvers described by the above equations consists of a broad innate inflammatory response as well as specific antibody responses for each antigenic variant of each type. Increasing concentrations of novel antigens stimulate both an inflammatory response and a delayed antibody response. As the antibody response increases on a dimensionless scale from 0 to 1, the initial inflammatory response is suppressed. After the antibody response appears, continued antigenic stimulation will drive up the adapted response until that antigenic variant is cleared, at which point both the antibody levels and the capacity to respond will decay over time. The capacity to generate specific antibodies upon challenge can be set to decay to a nonzero value Y memory , allowing a shorter latency to response upon reinfection. Higher values of Y memory yield stronger immune response and lower rates of detected parasitaemia upon reinfection, and this results in lower parasite densities in the adult human population. The larger the antigenic population, the more infections and thus time it takes to acquire broad parasitological immunity.
The development of anti-parasite immunity over age and exposure can be modelled with functions of age and exposure , or with exposure and infections with an age parameter in biting exposure . Part of the challenge facing mathematical models is that over a wide range of transmission intensities there is a high level of detected parasitaemia through the age of 10 followed by a decline to a lower level after the age of 20 [14, 22, 45]. The total exposure experienced varies by over an order of magnitude, and yet the development of immunity limiting the rate of detection of parasite infection varies much less. These dynamics can be fit by incorporating a strong age effect in the immune functions , but it is unclear whether such a strong mechanism with these age-linked characteristics exists. In addition, many existing models treat every infection the same, and ignore antigenic diversity across clones.
Figures 3 and 4 show the change in detected parasite prevalence and geometric mean parasitaemia smoothed with a one-year window for Namawala seasonality and either baseline vector populations, 0.5 baseline, or 0.2 baseline. The “0.5 x Namawala” simulation exhibits approximately half of the adult vectors per person as experienced in the baseline Namawala simulation. The actual EIR can vary depending on the parasite reservoir in the human population. The number of available merozoite and minor epitope combinations sets the asymptote for older age groups, and increasing the number of PfEMP-1 variants in the population increases the time for approaching this asymptote. All of the (100, 100, N) trajectories end with the same prevalence, but the one with N = 1,000 PfEMP-1 variants takes much longer to reach this level and begins to approach patterns seen in the Namawala data , which is displayed in panel 3C as point data. Simulations with too little parasite population diversity, for example (4, 20, 300), reach age-prevalence equilibrium after just 3–5 years, which is not realistic. Increasing the number of local parasite antigenic variants drives the age-prevalence curves to converge towards the actual data. This provides a simple, parsimonious, and mechanistic explanation for these curves.
Figures 5 and 6 repeat the analysis for Garki weather and EIRs ranging from a value characteristic of Rafin Marke (panel A) to one characteristic of Sugungum (panel C). Mean detected prevalence by age is plotted for Rafin Marke (panel A) and Sugungum (panel C) for comparison [14, 22], using a detection sensitivity of 10 parasites per microliter . Mean parasitaemia decreases from approximately 1,000 to below 100 parasites per microliter, as seen in the parasitological data [14, 22]. Immune memory levels were set to 0.3, which represents high levels of immune memory within the low range of memory levels permitting observed reinfection dynamics, but such an immune memory with high parasite diversity produces an excellent approximation of observed patterns. With the (4, 20, 300) setting, an immune memory level of 0.225 can recreate the prevalence asymptote seen in field studies, but prevalence by age approaches the adult asymptote with a time-constant under four years. In addition, with a lower memory level of 0.225 the higher parasite diversity simulations asymptote at unrealistically high levels of detected parasitaemia. Thus the higher memory and diversity is preferred.
In addition to the temporally smoothed results for parasite prevalence and parasite density averaged by age, it is important to examine the temporal dynamics over the course of a year by age group. The magnitude of oscillations in detected prevalence over the year can change with transmission intensity, depth of the dry season, and age of the cohort. As seen in Figure 1, transmission is highly seasonal in Namawala, and yet field studies showed remarkably low seasonal variation in detected parasitaemia . In Figure seventy-eight of the Garki Project report , temporal outputs from the Garki model are compared to the temporal pattern exhibited for each age range in the villages of Sugungum (highest transmission) and Rafin Marke (lowest transmission). Not much seasonality is observed for the younger age groups in Sugungum (although the 1–4 age group oscillates from 0.75 to 0.97 in Rafin Marke), but in the older age groups, a strong seasonal oscillation develops in detected parasitaemia. In Sugungum, in the 15–24 year age group, the low in April is 0.35 and the high in November is 0.55. In the 25–44 age group, the April low is 0.25 and the high of 0.45 was observed in August. In Rafin Marke, the 15–24 year age group has an oscillation from 0.38 in March to 0.8 in December, a dramatically high magnitude oscillation. The 25–44 age group ranges from 0.18 in March to 0.5 in October, a degree of oscillation underestimated by models with slower decay constants, including the Garki model itself .
The simulations presented in the foregoing demonstrate that parasite population diversity together with non-sterilizing immunological memory can explain the development of parasitological immunity with minimal ad hoc assumptions. The current model does not have an increase in biting rate exposure as age and body surface area increase, and thus modelled parasite prevalences for infants under one year of age are higher than observed, but incorporating age-dependent biting and maternal immunity as has been done elsewhere [14, 50] would effectively shift the Rafin Marke, Sugungum, and Namawala curves rightward for the first year and produce an even better match to observations than reported above. Introduction of new parasite variants to the local population over time rather than all at the start of the simulation would extend the development of immunity to older age groups as well. These features were left out of present simulations to illustrate what can be achieved merely through the mechanism of immunity to an ever-growing repertoire of specific antigens. The role of such model platforms is not to produce a single best-fit set of parameters but to provide a simple and transparent way to understand the roles of specific biological and physiological mechanisms. The observed shift in age-prevalence curves with increasing parasite population diversity provides a mechanistic bridge between current antigen-agnostic models [14, 22, 26, 34] and detailed antigen-specific field studies [10, 17, 41, 51, 52].
The presented simulations provide a simple and mechanistic explanation for a variety of field observations. First of all, the age-prevalence curve for Rafin Marke achieves the field-observed peak prevalence, although this peak is attained at too early an age. In addition, the convergence in the twenty years and older age group is an excellent match. In contrast, previous models tend to have too low a peak age-prevalence for Rafin Marke and too slow a decline to the old-age asymptote . The present results achieve a more accurate peak and convergence to old-age asymptotic behaviour for Sugungum as well. The results for both villages have better matches on a linear age scale than many pre-existing models. The present simulations fall short on the under-one year old age group, but incorporation of maternal antibodies and decreased biting rates would resolve this issue. These features are left out to provide a simpler illustration of the acquisition of immunity to an antigenic repertoire over time. Finally, the presented simulations perform very well on temporal dynamics. The high oscillation in detected parasitaemia in older age groups in Garki villages is matched more accurately than by the Garki model itself . And the present model provides a clear mechanistic interpretation of the low levels of seasonal variability in detected parasitaemia in Namawala as observed  and predicts an increase in seasonal variability as transmission levels decrease.
It is important to note that the model was not fitted to the parasitological data presented: the intrahost dynamics were parameterized from malariatherapy data and laboratory studies , and the vector population dynamics were constrained by entomological data [22, 45]. Thus, although even better fits could be achieved by allowing all parameters to vary, the goal of the present work is to fix parameters that are better understood from clinical research datasets and then to discern the effects of less understood parameters. The effects of parasite population size are compared to field data which were not used to generate or fit the model. The success of parasite diversity in explaining age-prevalence patterns does not rule out age-dependent effects, and models which use age factors as a proxy for immunity will work very well and will be substantially computationally faster. The present model has a relatively high dimensionality and a correspondingly slow computational rate which make it suboptimal for some approaches, although the level of biomedical detail it admits supports study of questions that are difficult to examine in simpler frameworks.
The present model provides a useful mechanistic framework for recreating and comparing field observations of antibody profiles and responses. In general, the repertoire of antibody responses to PfEMP-1 variants increases with age in malaria endemic regions , and clinical disease tends to correspond to gaps in the antibody repertoire . The role of transmission in the acquisition of this immunity can be seen in the similarity of specific antibody profiles in young children from the same region, having been exposed to contemporaneous parasites . The present model incorporates known mechanisms such as immunological memory, which is observed for malaria antigens , while having antibody responses that decrease rapidly from hyper-immune levels post infection and exhibit longer-term decay on faster time scales than measles . The power of these model structures is that they can examine these effects and implement known antigen-specific concepts and thus provide a mechanistic explanation for age-prevalence patterns of existing antigen-agnostic models.
Mechanistic models of vector dynamics can be exercised to study how sensitivities of local vector population dynamics to weather and other environmental variables can drive sensitivities to weather in parasite dynamics measured at the human population level. Environmental parameters such as temperature, rainfall, ground moisture and humidity can be used to construct static risk maps in the absence of interventions [54, 55], and the impact of seasonality has been examined with careful mathematical modelling . The present simulations demonstrate the impact of dry-season rainfall, vector refugia, and the mix of species on seasonal population dynamics. In the Garki weather data, no rainfall was recorded in the data for over six months. In a rainfall-driven model, drought of this depth and duration induces local extinction of the An. gambiae sl. populations. As a result, a baseline precipitation level of 1 mm per two weeks was artificially added to the data. In reality, this could be provided by low levels of rain that would not register in a daily rain measurement or by semi-permanent standing-water sources. Alternatively, there may be local fadeouts of the vector population requiring influx from refugia at the start of the rainy season when larger areas of habitat will become available. In this case, such refugia become very important in near-elimination conditions. The mix of species has strong effects on the seasonality as well. An. funestus is modelled with a slow integration of rainfall, with populations peaking at the end of the rainy season and decaying slowing into the dry season. This provides a low level of transmission in the first part of the dry season and has a substantial effect on the time-trough in parasitaemia. A combination of these factors in the present model can be used to explain both the lack of a strong seasonal oscillation in Namawala parasitaemia and the much more pronounced oscillation observed for adults in the Garki district.
Other important factors when comparing model simulation outputs to field data include the diagnostic sensitivity for parasitaemia. Each day of the simulation, each in silico individual has a simulated blood smear test which finds any parasites in 0.1 μL of blood, corresponding to a sensitivity of 10 parasites/μL as described in the malariatherapy analyses and the Garki report . These were specialized studies and are not characteristic of typical slide sensitivities. Reducing the sensitivity to more typical values for field microscopy, such as 50 parasites/μL, reduces the measured prevalence and increases the mean detected parasitaemia, because all detected parasite densities are at least at the sensitivity threshold.
Other sensitive and uncertain parameters concern human infectivity to mosquitoes, which depends on a variety of factors [11, 12]. Such infectivity has been studied in previous models [37, 57, 58], and the present model can be used to understand the effect of the development of immunity to the asexual stages upon the degree of transmission. The effects of decreased gametocyte densities are balanced by those of decreases in the level of inflammatory cytokines, which have been shown to limit transmission success of gametocytes [46, 59]. Results from new field studies will advance understanding of the roles of different components of the human population comprising the human transmission reservoir, e.g., when studied by age and symptom-structures.
It is important to understand system sensitivities to parameters and model assumptions governing within-host and mosquito-borne transmission dynamics. Mechanistic models can aid in exploring the effects of assumptions and the weights of each alternative in open questions. The present model recreates baseline transmission intensity and transmission patterns. Modifications to local parasite population antigenic diversity can match changes in temporal dynamics by age and smoothed detected prevalence by age. The mechanistic framework of antigen-specific immune responses and parasite population antigenic diversity provides a simple and parsimonious explanation for observed prevalence trends.
Key areas for further investigation include factors related to human immunity and its development over the course of exposure [14, 32, 37] and the duration of infections [24, 60]. The malariatherapy data comprise an important resource for understanding within-host dynamics and informing population-level models , but they are limited in the questions they can be used to answer. New data will be required to support development of a fuller understanding of the effects of parasite diversity upon population dynamics, and models can be exercised to demonstrate the importance of such studies. Understanding present levels of parasite diversity and how they change with transmission comprises a potentially powerful new synergy between modelling and field studies.
In summary, mathematical models can demonstrate the implications of basic parasitological research and vector ecosystematics for the design and operation of the global campaign for malaria eradication [7, 61–63]. Levels of population immunity affect the ability to measure parasitaemia, to correctly estimate the scales of the human infectious reservoir, and to detect, contain and eliminate parasite reintroductions to previously-cleared areas. Modern computer-based mathematical modelling provides a rich set of opportunities to explore model spaces that were too expensive to work with even a decade ago, and thus models such as the present one will become increasingly useful tools to the broader malaria research community.
The author thanks Bill and Melinda Gates for their active support of this work and their sponsorship through the Global Good Fund. Useful discussions with colleagues at Intellectual Ventures Laboratory, the Malaria Program of the Bill and Melinda Gates Foundation, VECNet, and MalERA are likewise greatly appreciated.
Financial support was provided by Bill and Melinda Gates through the Global Good Fund.
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- World Health Organization: World malaria report: 2011. 2011, Geneva: World Health OrganizationGoogle Scholar
- Hay SI, Guerra CA, Gething PW, Patil AP, Tatem AJ, Noor AM, Kabaria CW, Manh BH, Elyazar IR, Brooker S, Smith DL, Moyeed RA, Snow RW: A world malaria map: Plasmodium falciparum endemicity in 2007. PLoS Med. 2009, 6: e1000048.PubMed CentralView ArticlePubMedGoogle Scholar
- Hay SI, Sinka ME, Okara RM, Kabaria CW, Mbithi PM, Tago CC, Benz D, Gething PW, Howes RE, Patil AP, Temperley WH, Bangs MJ, Chareonviriyaphap T, Elyazar IRF, Harbach RE, Hemingway J, Manguin S, Mbogo CM, Rubio-Palis Y, Godfray HCJ: Developing global maps of the dominant Anopheles vectors of human malaria. PLoS Med. 2010, 7: e1000209-10.1371/journal.pmed.1000209.PubMed CentralView ArticlePubMedGoogle Scholar
- Global Rural–urban Mapping Project (GRUMP): Alpha Version.http://sedac.ciesin.columbia.edu/gpw,
- De Castro MC, Yamagata Y, Mtasiwa D, Tanner M, Utzinger J, Keiser J, Singer BH: Integrated urban malaria control: a case study in Dar Es Salaam Tanzania. Am J Trop Med Hyg. 2004, 71: 103-117.Google Scholar
- Alonso PL, Brown G, Arevalo-Herrera M, Binka F, Chitnis C, Collins F, Doumbo OK, Greenwood B, Hall BF, Levine MM, Mendis K, Newman RD, Plowe CV, Rodríguez MH, Sinden R, Slutsker L, Tanner M: A research agenda to underpin malaria eradication. PLoS Med. 2011, 8: e1000406-10.1371/journal.pmed.1000406.PubMed CentralView ArticlePubMedGoogle Scholar
- The mal ERACGoM: A research agenda for malaria eradication: modeling. PLoS Med. 2011, 8: e1000403.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Bull PC, Lowe BS, Kortok M, Molyneux CS, Newbold CI, Marsh K: Parasite antigens on the infected red cell surface are targets for naturally acquired immunity to malaria. Nat Med. 1998, 4: 358-360. 10.1038/nm0398-358.PubMed CentralView ArticlePubMedGoogle Scholar
- Sinden RE, Butcher GA, Billker O, Fleck SL: Regulation of Infectivity of Plasmodium to the mosquito vector. Adv Parasitol. 1996, 38: 53-117.View ArticlePubMedGoogle Scholar
- Drakeley C, Sutherland C, Bousema JT, Sauerwein RW, Targett GAT: The epidemiology of Plasmodium falciparum gametocytes: weapons of mass dispersion. Trends Parasitol. 2006, 22: 424-430. 10.1016/j.pt.2006.07.001.View ArticlePubMedGoogle Scholar
- Langhorne J, Ndungu FM, Sponaas AM, Marsh K: Immunity to malaria: more questions than answers. Nat Immunol. 2008, 9: 725-732.View ArticlePubMedGoogle Scholar
- Maire N, Smith T, Ross A, Owusu-Agyei S, Dietz K, Molineaux L: A model for natural immunity to asexual blood stages of Plasmodium falciparum malaria in endemic areas. Am J Trop Med Hyg. 2006, 75: 19-31.PubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Bull PC, Lowe BS, Kortok M, Marsh K: Antibody recognition of Plasmodium falciparum erythrocyte surface antigens in Kenya: evidence for rare and prevalent variants. Infect Immun. 1999, 67: 733-739.PubMed CentralPubMedGoogle Scholar
- Fowkes FJI, Richards JS, Simpson JA, Beeson JG: The relationship between anti-merozoite antibodies and incidence of Plasmodium falciparum malaria: a systematic review and meta-analysis. PLoS Med. 2010, 7: e1000218-10.1371/journal.pmed.1000218.PubMed CentralView ArticlePubMedGoogle Scholar
- Stevenson MM, Riley EM: Innate immunity to malaria. Nat Rev Immunol. 2004, 4: 169-180. 10.1038/nri1311.View ArticlePubMedGoogle Scholar
- Ross R: The prevention of malaria. 1910, London: J. MurrayGoogle Scholar
- Macdonald G: Theory of the eradication of malaria. Bull World Health Organ. 1956, 15: 369-387.PubMed CentralPubMedGoogle Scholar
- Macdonald G: The epidemiology and control of malaria. 1957, London: Oxford University PressGoogle Scholar
- Molineaux L, Gramiccia G: The Garki Project: Research on the Epidemiology and Control of Malaria in the Sudan Savanna of West Africa. 1980, Geneva: World Health OrganizationGoogle Scholar
- McKenzie FE, Baird JK, Beier JC, Lal AA, Bossert WH: A biologic basis for integrated malaria control. Am J Trop Med Hyg. 2002, 67: 571-577.PubMed CentralPubMedGoogle Scholar
- Smith DL, McKenzie FE, Snow RW, Hay SI: Revisiting the basic reproductive number for malaria and its implications for malaria control. PLoS Biol. 2007, 5: e42-10.1371/journal.pbio.0050042.PubMed CentralView ArticlePubMedGoogle Scholar
- 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: 1-10.PubMedGoogle Scholar
- Griffin JT, Hollingsworth TD, Okell LC, Churcher TS, White M, Hinsley W, Bousema T, Drakeley CJ, Ferguson NM, Basáñez M-G, Ghani AC: Reducing Plasmodium falciparum malaria transmission in Africa: a model-based evaluation of intervention strategies. PLoS Med. 2010, 7: e1000324-10.1371/journal.pmed.1000324.PubMed CentralView ArticlePubMedGoogle Scholar
- White L, Maude R, Pongtavornpinyo W, Saralamba S, Aguas R, Van Effelterre T, Day N, White N: The role of simple mathematical models in malaria elimination strategy design. Malar J. 2009, 8: 212-10.1186/1475-2875-8-212.PubMed CentralView ArticlePubMedGoogle Scholar
- Molineaux L, Diebner HH, Eichner M, Collins WE, Jeffery GM, Dietz K: Plasmodium falciparum parasitaemia described by a new mathematical model. Parasitology. 2001, 122: 379-391.View ArticlePubMedGoogle Scholar
- McKenzie FE, Bossert WH: The dynamics of Plasmodium falciparum blood-stage infection. J Theor Biol. 1997, 188: 127-140. 10.1006/jtbi.1997.0478.View ArticlePubMedGoogle Scholar
- Metcalf CJE, Graham AL, Huijben S, Barclay VC, Long GH, Grenfell BT, Read AF, Bjørnstad ON: Partitioning regulatory mechanisms of within-host malaria dynamics using the effective propagation number. Science. 2011, 333: 984-988. 10.1126/science.1204588.PubMed CentralView ArticlePubMedGoogle Scholar
- Recker M, Nee S, Bull PC, Kinyanjui S, Marsh K, Newbold C, Gupta S: Transient cross-reactive immune responses can orchestrate antigenic variation in malaria. Nature. 2004, 429: 555-558. 10.1038/nature02486.View ArticlePubMedGoogle Scholar
- Paget-McNicol S, Gatton M, Hastings I, Saul A: The Plasmodium falciparum var gene switching rate, switching mechanism and patterns of parasite recrudescence described by mathematical modelling. Parasitology. 2002, 124: 225-235.View ArticlePubMedGoogle Scholar
- Saul AJ, Graves PM, Kay BH: A cyclical feeding model for pathogen transmission and its application to determine vectorial capacity from vector infection rates. J Applied Ecology. 1990, 27: 123-133. 10.2307/2403572.View ArticleGoogle Scholar
- Smith DL, McKenzie FE: Statics and dynamics of malaria infection in Anopheles mosquitoes. Malar J. 2004, 3: 13-10.1186/1475-2875-3-13.PubMed CentralView ArticlePubMedGoogle Scholar
- Depinay J-M, Mbogo C, Killeen G, Knols B, Beier J, Carlson J, Dushoff J, Billingsley P, Mwambi H, Githure J, Toure A, Ellis McKenzie F: A simulation model of African Anopheles ecology and population dynamics for the analysis of malaria transmission. Malar J. 2004, 3: 29-10.1186/1475-2875-3-29.PubMed CentralView ArticlePubMedGoogle Scholar
- Bomblies A, Duchemin J-B, Eltahir E: A mechanistic approach for accurate simulation of village scale malaria transmission. Malar J. 2009, 8: 223-10.1186/1475-2875-8-223.PubMed CentralView ArticlePubMedGoogle Scholar
- Filipe JA, Riley EM, Drakeley CJ, Sutherland CJ, Ghani AC: Determination of the processes driving the acquisition of immunity to malaria using a mathematical transmission model. PLoS Comput Biol. 2007, 3: e255-10.1371/journal.pcbi.0030255.PubMed CentralView ArticlePubMedGoogle Scholar
- Okell LC, Drakeley CJ, Bousema T, Whitty CJ, Ghani AC: Modelling the impact of artemisinin combination therapy and long-acting treatments on malaria transmission intensity. PLoS Med. 2008, 5: e226-10.1371/journal.pmed.0050226. discussion e226PubMed CentralView ArticlePubMedGoogle Scholar
- Gupta S, Day KP: A strain theory of malaria transmission. Parasitology Today. 1994, 10: 476-481. 10.1016/0169-4758(94)90160-0.View ArticlePubMedGoogle Scholar
- Collins WE, Jeffery GM: A retrospective examinationof secondary sporozoite- and trophozoite-induced infections with Plasmodium falciparum: development of parasitologic and clinical immunity following secondary infection. Am J Trop Med Hyg. 1999, 61: 20-35.View ArticlePubMedGoogle Scholar
- Barry AE, Trieu A, Fowkes FJI, Pablo J, Kalantari-Dehaghi M, Jasinskas A, Tan X, Kayala MA, Tavul L, Siba PM, Day KP, Baldi P, Felgner PL, Doolan DL: The stability and complexity of antibody responses to the major surface antigen of Plasmodium falciparum are associated with age in a malaria endemic area. Mol Cell Proteomics. 2011, 10:Google Scholar
- Eckhoff P: A malaria transmission-directed model of mosquito life cycle and ecology. Malar J. 2011, 10: 303-10.1186/1475-2875-10-303.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Molineaux L, Storey J, Cohen JE, Thomas A: A longitudinal study of human malaria in the West African Savanna in the absence of control measures: relationships between different Plasmodium species, in particular P. falciparum and P. malariae. Am J Trop Med Hyg. 1980, 29: 725-737.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Naotunne TS, Karunaweera ND, Mendis KN, Carter R: Cytokine-mediated inactivation of malarial gametocytes is dependent on the presence of white blood cells and involves reactive nitrogen intermediates. Immunology. 1993, 78: 555-562.PubMed CentralPubMedGoogle Scholar
- Full Data Reanalysis.http://gpcc.dwd.de,
- Holder AA, Patino JAG, 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.PubMedGoogle Scholar
- Bull PC, Marsh K: The role of antibodies to Plasmodium falciparum-infected-erythrocyte surface antigens in naturally acquired immunity to malaria. Trends Microbiol. 2002, 10: 55-58. 10.1016/S0966-842X(01)02278-8.View ArticlePubMedGoogle Scholar
- 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: 11-18.PubMedGoogle Scholar
- Fowkes FJI, McGready R, Cross NJ, Hommel M, Simpson JA, Elliott SR, Richards JS, Lackovic K, Viladpai-nguen J, Narum D, Tsuboi T, Anders RF, Nosten F, Beeson JG: New insights into acquisition, boosting and longevity of immunity to malaria in pregnant women. J Infect Dis. 2012, 206: 1612-1621. 10.1093/infdis/jis566.PubMed CentralView ArticlePubMedGoogle Scholar
- McCallum FJ, Persson KEM, Mugyenyi CK, Fowkes FJI, Simpson JA, Richards JS, Williams TN, Marsh K, Beeson JG: Acquisition of growth-inhibitory antibodies against blood-stage Plasmodium falciparum. PLoS ONE. 2008, 3: e3571-10.1371/journal.pone.0003571.PubMed CentralView ArticlePubMedGoogle Scholar
- Collins WE, Jeffery GM: A retrospective examination of sporozoite- and trophozoite-induced infections with Plasmodium falciparum: development of parasitologic and clinical immunity during primary infection. Am J Trop Med Hyg. 1999, 61: 4-19.View ArticlePubMedGoogle Scholar
- Craig MH, Snow RW, le Sueur D: A climate-based distribution model of malaria transmission in sub-Saharan Africa. Parasitol Today. 1999, 15: 105-111. 10.1016/S0169-4758(99)01396-4.View ArticlePubMedGoogle Scholar
- Kiszewski A, Mellinger A, Spielman A, Malaney P, Sachs SE, Sachs J: A global index representing the stability of malaria transmission. Am J Trop Med Hyg. 2004, 70: 486-498.PubMedGoogle Scholar
- McKenzie FE, Killeen GF, Beier JC, Bossert WH: Seasonality, parasite diversity, and local extinctions in Plasmodium falciparum malaria. Ecology. 2001, 82: 2673-2681.PubMed CentralView ArticlePubMedGoogle Scholar
- Killeen GF, Ross A, Smith T: Infectiousness of malaria-endemic human populations to vectors. Am J Trop Med Hyg. 2006, 75: 38-45.PubMedGoogle Scholar
- Ross A, Killeen G, Smith T: Relationships between host infectivity to mosquitoes and asexual parasite density in Plasmodium falciparum. Am J Trop Med Hyg. 2006, 75: 32-37.PubMedGoogle Scholar
- Naotunne TS, Karunaweera ND, Del Giudice G, Kularatne MU, Grau GE, Carter R, Mendis KN: Cytokines kill malaria parasites during infection crisis: extracellular complementary factors are essential. J Exp Med. 1991, 173: 523-529. 10.1084/jem.173.3.523.View ArticlePubMedGoogle Scholar
- Sama W, Dietz K, Smith T: Distribution of survival times of deliberate Plasmodium falciparum infections in tertiary syphilis patients. Trans R Soc Trop Med Hyg. 2006, 100: 811-816. 10.1016/j.trstmh.2005.11.001.View ArticlePubMedGoogle Scholar
- Cohen JM, Smith DL, Vallely A, Taleo G, Malefoasi G, Sabot O: Holding the Line. Shrinking the Malaria Map: A Prospectus on Malaria Elimination. Edited by: Feachem RGA, Phillips AA, Targett GA. 2009, San Francisco: The Global Health Group: UCSF Global Health Sciences, 40-60.Google Scholar
- Moonen B, Barrett S, Tulloch J, Jamison DT, Moonen B, Barrett S, Tulloch J, Jamison DT: Making the Decision. Shrinking the Malaria Map: A Prospectus on Malaria Elimination. Edited by: Feachem RGA, Phillips AA, Targett GA, Feachem RGA, Phillips AA, Targett GA. 2009, San Francisco: The Global Health Group: UCSF Global Health Sciences, 1-18.Google Scholar
- Smith DL, Smith TA, Hay SI: Measuring Malaria for Elimination. Shrinking the Malaria Map: A Prospectus on Malaria Elimination. Edited by: Feachem RGA, Phillips AA, Targett GA. 2009, San Francisco: The Global Health Group: UCSF Global Health Sciences, 108-126.Google Scholar
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 (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.