Fun with maths: exploring implications of mathematical models for malaria eradication
© Eckhoff et al.; licensee BioMed Central. 2014
Received: 15 September 2014
Accepted: 6 December 2014
Published: 11 December 2014
Mathematical analyses and modelling have an important role informing malaria eradication strategies. Simple mathematical approaches can answer many questions, but it is important to investigate their assumptions and to test whether simple assumptions affect the results. In this note, four examples demonstrate both the effects of model structures and assumptions and also the benefits of using a diversity of model approaches. These examples include the time to eradication, the impact of vaccine efficacy and coverage, drug programs and the effects of duration of infections and delays to treatment, and the influence of seasonality and migration coupling on disease fadeout. An excessively simple structure can miss key results, but simple mathematical approaches can still achieve key results for eradication strategy and define areas for investigation by more complex models.
Mathematical models are a helpful tool for testing assumptions and elucidating the quantitative implications of disease features. Most of the time, results of a mathematical modelling study match mental models of expert opinion, which have been shaped by years of research and experience. Occasionally, surprising or counterintuitive results emerge from models. This can happen with simple compartmental models, but sometimes the surprising nature of modelling results is due to mistaken intuition based on mental or simplified mathematical models taken too far. Such findings must be examined closely, but they enhance understanding when validated. This note examines four examples from malaria eradication of how modelling can produce useful results that may initially strike one as unanticipated, especially if one relies too heavily on simple modelling concepts or mental pictures.
Time to eradication
The classic epidemiological model concept of the basic reproductive number R0 nominally represents the secondary infections resulting from a single primary infection in an immunologically naïve population. Macdonald developed the original formulation for malaria [1, 2], and others have extended the theory over time . This concept has become ubiquitous in epidemiological modelling, and while it can be helpful, it can also cause confusion and even misguided planning if used incorrectly. The basic theory is that if R0 is reduced to a controlled reproductive number Rc, the number of secondary infections resulting from a single primary infection under control measures, below 1, then the disease will disappear eventually. Each infection is replaced by less than one infection, and prevalence exponentially decays to zero. While this simplified equilibrium view is true, planning eradication efforts requires investigation of the dynamics and timelines as described by Macdonald and others [4, 5], not just the equilibrium conditions.
When human infections are not actively cleared and elimination is driven by attaining Rc <1 solely through transmission reductions, this slow decay has serious programme implications that limit the chances of elimination. For many diseases the durations of interventions that reduce transmission, such as vaccines, are much longer than the durations of infection. As such, transmission reductions that reduce R0 below 1 can succeed in driving local eliminations of the pathogen while the interventions are still efficacious without replacement. For malaria however, the average infection duration of six months is only a factor of six away from three-year lifetimes of insecticide-treated bed nets (ITNs). Indoor residual spraying (IRS) campaigns need to be repeated every two to six months depending on the active ingredient  in areas with year-round transmission and need to be repeated yearly in highly seasonal settings. These frequencies are considerably faster than the timelines to elimination based on transmission reductions alone. The current lead vaccine candidate, RTS,S , has a duration of efficacy that is not longer than ITNs. As seen in Figure 1, the infectious reservoir may still be quite robust by the time these interventions need to be replaced, or else the reproductive number will rise back above 1 without having achieved elimination.
The result of this matching of decay constants and intervention durability drives elimination programmes down one of two routes. The first is that elimination campaigns need to reduce transmission for many years to maintain Rc <1, with repeated distributions. The second possible route is to actively clear the infectious reservoir with drug-based campaigns in order to reduce the duration of infections and to achieve rapid drops in prevalence. Reducing the duration of infections in fact also drives further decreases in Rc.This makes malaria more like diseases of shorter infectious duration and more susceptible to transient transmission-driven reductions of Rc <1.
Coverage and R0
In contrast to many common childhood diseases with R0 ranging from 1.5 to 12, malaria can exhibit an R0 over 100 in well-suited settings . Traditional epidemiological modelling shows that the vaccine coverage required to interrupt transmission with a perfect 100% effective vaccine without decay can be calculated as 1-1/R0. The higher the value of R0, the higher the covered fraction needs to be. If R0 = 2, then exceeding 50% coverage will result in Rc <1. If R0 = 10, then 90% coverage is required.
This relationship also has implications for the impact of vaccine efficacy. If a vaccine instead of preventing all infections prevents only a fraction Veff of infections, with Veff <1, then the coverage required to achieve Rc <1 becomes (1-1/R0)/Veff. At some combinations of R0 and Veff, this coverage rises above 100% and elimination becomes impossible even at full coverage. So if the vaccine is leaky and exhibits 80% efficacy at preventing infections, the relationship with maximum R0 shifts dramatically at the high end of coverage, as seen in Figure 2. Below 80% coverage, the leaky vaccine can eliminate R0’s close to the R0’s eliminated by the perfect vaccine. At 90% coverage, however, where the perfect vaccine can eliminate an R0 of 10, the leaky vaccine can only eliminate R0 up to 4. Even perfect coverage only approaches a maximum R0 of 5. The situation becomes even more difficult if vaccine efficacy decays over time, requiring repeated boosting vaccination, higher coverage, higher initial efficacy, or all three.
The result of this relationship is to quantify the value of achieving higher intervention coverage. It is commonly accepted that it becomes increasingly costly to increase coverage at higher baseline coverages. To achieve greater than 90% coverage might be prohibitively expensive. An important counterpoint to these dramatically rising costs of achieving coverage is that one actually buys a much higher eliminate-able R0. The high baseline R0 exhibited by malaria in some pockets also drives non-linearly increasing value to higher vaccine efficacies. Increasing vaccine efficacies has a challenging set of scientific hurdles and increasing coverage has logistical hurdles, but understanding what is required for success can quantify the value of surmounting either or both sets of hurdles. A similar effect of non-linearly increasing benefits of coverage has been demonstrated in model studies of drug impacts for elimination . In summary, although coverage increases become increasingly expensive as baseline coverage increases, the benefits correspondingly grow.
Case detection delays and transmission-blocking drugs
There has recently been significant debate about the possible role of low-dose primaquine in elimination scenarios . The basic argument builds from achievable coverage and time constants. If the average duration of a falciparum infection is approximately six months, and the duration of gametocytaemia following clearance with schizonticidal anti-malarial drugs approximately two weeks, then clearing an infection reduces 12 times more onward transmission than adding low-dose primaquine to the schizonticidal regimen. Thus, if MSAT campaigns are reaching 50% of the population with an artemisinin-combination therapy (ACT) to clear infections, adding primaquine to the MSAT regimen for all 50% would only reduce the same amount of transmission as increasing ACT coverage to 54%. Limitations in who could receive primaquine would further reduce this marginal impact. As has previously been pointed out , primaquine only begins to have a discernable impact at high coverages, similar to that discussed in the previous section.
For drug-based clearance of the infectious reservoir, the most important factors remain achieving high coverage with ACT and reducing the delay in case finding, but the impacts of realistic delays in case finding and data-driven distributions of infection duration constitute an important, although minor, effect. This is an example in which excessively simple assumptions may get the basic result correct but miss an important nuance, while somewhat more detailed mathematical approaches can capture realistic effects and improve accuracy. Realistic delays in case finding and data-driven distributions of infection duration together reduce the target ACT coverage at which adding a transmission blocking drug (not necessarily primaquine) to the ACT regimen begins to show equivalent epidemiological impact to a fixed fraction increase in coverage with the schizonticidal regimen alone. Compared to the previous two examples, this impact of infection senescence is a minor effect. However, this effect can still be elucidated with mathematical analyses.
Metapopulations and synchrony
While malaria transmission is very much a local phenomenon driven by mixing of mosquitoes and humans on short spatial scales, malaria eradication is a global issue. Areas are connected spatially by human migration and elimination is not possible without considering the spatial context. Excellent work has been done to show the impact of human mobility on malaria , and how the coupling of regions raises the difficulty of achieving elimination. Human movements superimposed on a landscape of heterogeneous transmission intensities and campaign coverages have the potential to reintroduce infection to temporarily cleared areas before the benefits of elimination ‘stickiness’ can be achieved, contributing to malaria resurgence [14, 15]. As such, initiating an elimination campaign is largely a regional, rather than a local decision .
This benefit of synchrony for achieving elimination has several implications for malaria eradication. The campaign benefits from seasonal forcing being similar in phase over locally connected areas, and the chance of elimination increases if anti-malarial treatments are used to shorten the length of infections sufficiently to have rapid decays of the infectious reservoir during the low transmission season. In addition, locations with more stable year-round transmission, even if they do not have the high spikes in case counts of highly seasonal settings, will be more difficult to eliminate and may require higher campaign coverages for a given annual average R0. Finally, achieving elimination across connected geographic regions will be easier if the external forcing on transmission, namely campaign pressures such as vector control spraying and anti-malarial drug mass screen-and-treating, are synchronized across the region of interest.
In summary, each of these four examples demonstrates the power of data-driven mathematical analyses for informing malaria eradication strategy. Mathematical modelling forces one to make assumptions explicit and reveals the implications of those assumptions. These examples also show how model structure can impact results and conclusions for policy while illustrating the importance of comparing results from simple and complex models. Dynamics and small populations matter in an eradication context, so checking equilibrium analyses with both stochastic and deterministic dynamic models is essential.
Each of the four examples in this note can be examined in the relatively simple frameworks here or in more complex microsimulation models. Simple mathematical analyses can often arrive at the correct big picture, while small modifications can capture nuances that improve accuracy, as in the third example. More complicated models can then confirm results while investigating whether the effects of realistic details affect policy conclusions or not. Since it would be possible to investigate irrelevant details in complex models forever, simple models have an additional important role in focusing attention on areas where the capacities of complex models could influence the answer. Starting with a simple framework is thus strategically useful.
It is important to make sure that seemingly straightforward conclusions from a simple model structure be rigorously examined with model structures that can accommodate improved realism [11, 18–20]. If the same result for a given question is obtained from simple and complex models, then one gains confidence in the results from the simpler structures. More complex microsimulation models can recreate and test results obtained in simple analyses and facilitate extensions of questions that would have been otherwise less tractable. Such questions raised by the examples above could include the interrelationships among decay of vaccine efficacy, heterogeneity in vaccine coverage, delivery strategy and schedule, transmission context, and coverage targets or the effects on disease fadeout of spatial variation in seasonality and local heterogeneity in risk. When used properly in concert, a diversity of mathematical models and approaches will play an essential role in eradication planning.
- 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, New York: Oxford University PressGoogle 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 D, Hay S: Endemicity response timelines for Plasmodium falciparum elimination. Malar J. 2009, 8: 87-10.1186/1475-2875-8-87.PubMed CentralView ArticlePubMedGoogle Scholar
- Macdonald G, Goeckel GW: The Malaria Parasite Rate and the Interruption of Transmission. Bull World Health Organ. 1964, 31: 365-377.PubMed CentralPubMedGoogle 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
- Sama W, Killeen G, Smith T: Estimating the duration of Plasmodium falciparum infection from trials of indoor residual spraying. Am J Trop Med Hyg. 2004, 70: 625-634.PubMedGoogle 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
- Pluess B, Tanser FC, Lengeler C, Sharp BL: Indoor residual spraying for preventing malaria. Cochrane Database of Systematic Reviews. 2010, Issue 4. Art. No.: CD006657. doi:10.1002/14651858.CD006657.pub2Google Scholar
- RTS SCTP: A Phase 3 trial of RTS, S/AS01 malaria vaccine in African infants. N Engl J Med. 2012, 367: 2284-2295.View ArticleGoogle Scholar
- Johnston GL, Gething PW, Hay SI, Smith DL, Fidock DA: Modeling within-host effects of drugs on Plasmodium falciparum transmission and prospects for malaria elimination. PLoS Comput Biol. 2014, 10: e1003434-10.1371/journal.pcbi.1003434.PubMed CentralView ArticlePubMedGoogle Scholar
- Gillespie DT: Exact stochastic simulation of coupled chemical reactions. J Phys Chem. 1977, 81: 2340-2361. 10.1021/j100540a008.View ArticleGoogle Scholar
- Wesolowski A, Eagle N, Tatem AJ, Smith DL, Noor AM, Snow RW, Buckee CO: Quantifying the impact of human mobility on malaria. Science. 2012, 338: 267-270. 10.1126/science.1223467.PubMed CentralView ArticlePubMedGoogle Scholar
- Smith DL, Cohen JM, Chiyaka C, Johnston G, Gething PW, Gosling R, Buckee CO, Laxminarayan R, Hay SI, Tatem AJ: A sticky situation: the unexpected stability of malaria elimination. Philos Trans R Soc London B: Biol Sci. 2013, 368:Google Scholar
- Cohen J, Smith D, Cotter C, Ward A, Yamey G, Sabot O, Moonen B: Malaria resurgence: a systematic review and assessment of its causes. Malar J. 2012, 11: 122-10.1186/1475-2875-11-122.PubMed CentralView ArticlePubMedGoogle Scholar
- Moonen B, Barrett S, Tulloch J, Jamison DT: Making the decision. Shrinking the Malaria Map: A Prospectus on Malaria Elimination. Edited by: Group ME, Feachem RGA, Phillips AA, Targett GA. 2009, San Francisco: The Global Health Group: UCSF Global Health Sciences, 1-18.Google Scholar
- Grenfell BT, Bolker BM, Kleczkowski A: Seasonality and extinction in chaotic metapopulations. Proc Biol Sci. 1995, 259: 97-103. 10.1098/rspb.1995.0015.View ArticleGoogle 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
- Wenger E, Eckhoff P: A mathematical model of the impact of present and future malaria vaccines. Malar J. 2013, 12: 126-10.1186/1475-2875-12-126.PubMed CentralView ArticlePubMedGoogle 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
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/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.