# Fun with maths: exploring implications of mathematical models for malaria eradication

- Philip A Eckhoff
^{1}Email author, - Caitlin A Bever
^{1}, - Jaline Gerardin
^{1}and - Edward A Wenger
^{1}

**13**:486

**DOI: **10.1186/1475-2875-13-486

© Eckhoff et al.; licensee BioMed Central. 2014

**Received: **15 September 2014

**Accepted: **6 December 2014

**Published: **11 December 2014

## Abstract

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.

## Background

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 R_{0} 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 [3]. 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 R_{0} is reduced to a controlled reproductive number R_{c}, 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.

_{0}tends to be implemented in simple compartmental models or ordinary differential equations, characterized by an exponential distribution for infection durations. A reduction of R

_{0}to R

_{c}<1 corresponds to exponential decay of the infected population. The fastest possible decay of the infected population without active clearance of infections is attained when R

_{c}= 0 and has the inverse of the average duration of infection as its decay constant. For short duration infections, such as flu, the infectious reservoir decays rapidly. Malaria, on the other hand, has average infection durations of up to six months [6–8]. As such, exponential decay of the infectious reservoir can take years under the best case, and over ten years if R

_{c}= 0.9. Figure 1 shows the decay of the malaria infectious reservoir for four values of R

_{c}, with infection duration averaging 180 days. A controlled reproductive number of 0.9, 0.7, 0.5, or even 0 exhibits slow declines in the total number of infected individuals. In contrast, a disease with R

_{c}= 0.9 but an infectious duration of just 3.5 days rapidly depletes its infectious reservoir.

When human infections are not actively cleared and elimination is driven by attaining R_{c} <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 R_{0} 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 [9] 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 [10], 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 R_{c} <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 R_{c}[3].This makes malaria more like diseases of shorter infectious duration and more susceptible to transient transmission-driven reductions of R_{c} <1.

## Coverage and R_{0}

In contrast to many common childhood diseases with R_{0} ranging from 1.5 to 12, malaria can exhibit an R_{0} over 100 in well-suited settings [3]. 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/R_{0}. The higher the value of R_{0}, the higher the covered fraction needs to be. If R_{0} = 2, then exceeding 50% coverage will result in R_{c} <1. If R_{0} = 10, then 90% coverage is required.

_{0}that can be overcome for each increase in coverage. With a perfect infection-preventing vaccine, one that provides complete sterilizing immunity in vaccinated individuals, high levels of coverage provide dramatic increases in the maximum R

_{0}that can be eliminated. Coverage of 50% can eliminate an R

_{0}of 2, coverage of 80% can eliminate an R

_{0}of 5, and coverage of 90% can eliminate an R

_{0}of 10. Above 90%, the gains become dramatic as seen in Figure 2: coverage of 95% will eliminate an R

_{0}of 20 and coverage of 99% will eliminate an R

_{0}of 100.

This relationship also has implications for the impact of vaccine efficacy. If a vaccine instead of preventing all infections prevents only a fraction V_{eff} of infections, with V_{eff} <1, then the coverage required to achieve R_{c} <1 becomes (1-1/R_{0})/V_{eff}. At some combinations of R_{0} and V_{eff}, 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 R_{0} shifts dramatically at the high end of coverage, as seen in Figure 2. Below 80% coverage, the leaky vaccine can eliminate R_{0}’s close to the R_{0}’s eliminated by the perfect vaccine. At 90% coverage, however, where the perfect vaccine can eliminate an R_{0} of 10, the leaky vaccine can only eliminate R_{0} up to 4. Even perfect coverage only approaches a maximum R_{0} 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 R_{0}. The high baseline R_{0} 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 [11]. 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 [11]. 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 [11], 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 [13], 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 [16].

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 R_{0}. 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.

## Conclusions

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.

## Declarations

## Authors’ Affiliations

## References

- 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

## Copyright

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.

## Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.