 Research
 Open Access
 Published:
Resurgence of malaria infection after mass treatment: a simulation study
Malaria Journal volume 18, Article number: 409 (2019)
Abstract
Background
Field studies are evaluating if mass drug administration (MDA) might shorten the time to elimination of Plasmodium falciparum malaria, when vector control measures and reactive surveillance strategies are scaledup. A concern with this strategy is that there may be resurgence of transmission following MDA.
Methods
A conceptual model was developed to classify possible outcomes of an initial period of MDA, followed by continuously implementing other interventions. The classification considered whether elimination or a new endemic stable state is achieved, and whether changes are rapid, transient, or gradual. These categories were informed by stability analyses of simple models of vector control, case management, and testandtreat interventions. Individualbased stochastic models of malaria transmission (OpenMalaria) were then used to estimate the probability and likely rates of resurgence in realistic settings. Effects of concurrent interventions, including routine case management and testandtreat strategies were investigated.
Results
Analysis of the conceptual models suggest resurgence will occur after MDA unless transmission potential is very low, or the postMDA prevalence falls below a threshold, which depends on both transmission potential and on the induction of bistability. Importation rates are important only when this threshold is very low. In most OpenMalaria simulations the approximately stable state achieved at the end of the simulations was independent of inclusion of MDA and the final state was unaffected by importation of infections at plausible rates. Elimination occurred only with high effective coverage of case management, low initial prevalence, and high intensity testandtreat. High coverage of case management but not by testandtreat induced bistability. Where resurgence occurred, its rate depended mainly on transmission potential (not treatment rates).
Conclusions
A short burst of high impact MDA is likely to be followed by resurgence. To avert resurgence, concomitant interventions need either to substantially reduce average transmission potential or to be differentially effective in averting or clearing infections at low prevalence. Case management at high effective coverage has this differential effect, and should suffice to avert resurgence caused by imported cases at plausible rates of importation. Once resurgence occurs, its rate depends mainly on transmission potential, not on treatment strategies.
Background
Substantial reductions in Plasmodium falciparum malaria transmission in Africa have been achieved by rolling out insecticidetreated nets and scaling up access to artemisininbased combination therapy [1]. In some previously highly endemic areas where transmission levels have been substantially reduced, additional intervention strategies are being trialled with a view to achieving elimination. In particular, in Southern Province of Zambia, the application of 3 dryseason rounds of mass testandtreat with artemether–lumefantrine was trialled, but failed to interrupt transmission [2]. A recent trial aims to test if the achievement of elimination by such a strategy can be accelerated by supplementing it with mass drug administration (MDA) using a longer lasting antimalarial combination (dihydroartemisinin–piperaquine, DHAP) [3, 4]. In an otherwise stable environment, the maximal impact on transmission of a malaria intervention programme generally occurs shortly after maximal scaleup [5] (Fig. 1).
This supports frontloading of drugbased elimination programmes and is the reason why malaria elimination programmes in the mid20th century envisaged an intense initial phase that would break transmission [6, 7]. Consequently, the MDA under investigation is conceived of as a timelimited pulsed intervention, with the objective of shortening timelines, primarily to accelerate the achievement of elimination. Though there have been trials of longterm preventive chemotherapy in endemic settings [8], MDA is not feasible as a longterm recurrent intervention, since it is difficult to organize and obtain consent for recurrent treatment of uninfected people at high coverage; also, it is expensive and potentially selects drugresistant parasites [9, 10]. Even if transmission is interrupted by the MDA, in the absence of effective surveillance reintroductions will occur and both R_{c} and prevalence will return to the initial stable states (blue lines in Fig. 1a–c) [5, 9]. This phenomenon is referred to as resurgence [11, 12].
Many different strategies might be considered as ways of averting resurgence. Most of these can be classified as either (i) vector control approaches, that reduce transmission potential without directly impacting the reservoir of infection in the human population; or (ii) chemoprevention strategies that provide both prophylaxis and clear new infections that have arisen. The latter might include recurrent focal MDA, focal testandtreat, or reactive case detection (RCD) strategies, in which neighbours of passively detected cases are either presumptively treated or screened for parasites and treated if positive.
This paper reports a conceptual model developed for classifying the possible outcomes when an enhanced intervention strategy (consisting of either intensified vector control, or chemoprevention) is initiated alongside a pulse of one or more rounds of highcoverage MDA. The conceptual model draws on the results of stability analyses that have been carried out for a range of malaria models. A largescale simulation experiment using the OpenMalaria (https://github.com/SwissTPH/OpenMalaria/wiki) microsimulation platform was then used to determine which of these outcomes are likely in realworld settings. The parameterizations were based on data from Southern Zambia [4, 13, 14]. In particular, simulated timeseries’ of prevalence were analysed to identify which characteristics of interventions determine whether their effectiveness in averting or slowing down resurgence.
Methods
Conceptual model
The conceptual model was illustrated by graphical representations of possible temporal trends in malariological indices induced by combining initial MDA with additional continual or recurrent interventions (Fig. 1). For simplicity, this figure neglects seasonality. For conciseness, the package of additional interventions, (which in principle may comprise vector control and/or chemoprevention programmes), is referred to as testandtreat. The indices considered were (i) transmission potential, R_{c}, defined as the average number of secondary infections due to each primary infection in the presence of control measures [15], and equal to the basic reproduction number, R_{0}, reduced proportionately by the effect of the interventions on each of the parameters; and (ii) the corresponding prevalence of infection.
The trends in Fig. 1 were then used to develop a classification of the potential outcomes of programmes. This analysis of trends is complemented by bifurcation diagrams [16] showing how the equilibrium prevalence achieved by the programme depends on R_{0} (Fig. 2). These diagrams are intended only to illustrate possible qualitative patterns. Both Figs. 1 and 2 are supported by qualitative arguments advanced in a previous publication [5, 17] and by the results of mathematical models published elsewhere [17,18,19,20].
Microsimulation model
Microsimulations were used to obtain approximate quantitative results applicable for more or less realistic transmission settings, taking into account factors of secondary importance such as the duration of chemoprevention and levels of acquired immunity. Using OpenMalaria (https://github.com/SwissTPH/OpenMalaria/wiki), a single site individualbased simulation model of malaria in humans [21] linked to deterministic model of malaria in mosquitoes [22]. The simulations track parasite densities, corresponding to different infection events, with different submodels for infection of humans [23], and for bloodstage parasite densities with the main immune effects controlling parasite densities and infection duration [24]. Clinical malaria, and hence routine case management and passive case detection, is triggered by hyperparasitaemia, with the models for both parasite densities and clinical incidence calibrating using the data of various field studies [24, 25]. OpenMalaria include various sources of variation between hosts in exposure and susceptibility, while treating the infections in mosquitoes as perfectly mixing.
Simulation of interventions
The simulations considered effects of MDA, of routine case management (treatment of clinical malaria attacks) and of the package of additional interventions (testandtreat). MDA was simulated via a single round of DHAP treatment, parameterized as described previously [26]. This implementation captures both the effects of truncating infections, and also the subsequent prophylactic effect based on pharmacokinetic/dynamic studies. The efficacy of the programme is defined as the relative percentage drop in average parasite prevalence in 2 to 10 year olds from the 3 years preintervention to the average value 7–10 years postMDA. Routine case management was simulated in OpenMalaria as described previously [13, 27]. The set of OpenMalaria simulations of testandtreat were parameterized by simulating a testandtreat strategy, where a positive from a rapid diagnostic test (with a diagnostic cutoff of 50 parasites per μl of blood and a specificity of 0.942, based on fitting to the data of Murray et al. [28] results in cure of bloodstage parasites at the next timestep.
The simulated testandtreat programme was additional to the routine case management system, and comprised treatment of a fixed number (\(\nu\)) of the neighbours of passivelydetected index cases. The number of index cases investigated each week (\(\upiota\)) was set either to a predefined maximum (\(\upiota_{\text{max} }\)), or to the number of cases presenting for treatment during 1 week, whichever was greatest. To simulate targeting, the number of simulated individuals to be tested in the same interval was inflated by the targeting index (\(\tau\)), equivalent to the ratio of the prevalence in all individuals tested (including asymptomatics) to that in the general population [17], which was set to a constant value within a simulation. The application of this approach to account for spatial clustering and heterogeneity in OpenMalaria is described in detail elsewhere [29]. The effective per capita rate of testing was then: \(\Phi = \frac{\iota \nu \tau }{N}\) per week, where \(N\) is the population size. Treatment rates implemented in OpenMalaria were rescaled to allow for the 5day simulation timestep.
Simulation experiment
For analysis of resurgence rates simulations were carried out using the OpenMalaria R0000 base model parameterization [30] with a simulated population size of 10,000 (other exploratory simulations that used the R0670 parameterization, which allows for heterogeneity in susceptibility to comorbidity, gave the same general results). The same pattern of seasonality of transmission based on the pattern for southern Zambia used in previous analyses [13] was used throughout. This was scaled to achieve a predetermined annual entomological inoculation rate during the initial phase of the simulation.
The preexisting intervention programme was assumed to include healthfacility based routine case management with a defined effective coverage, E_{14}, (per 14 day period) [31] to have been in place for long enough for an approximate stable endemic state to have been reached by the time the testandtreat programme is introduced. Values of E_{14} of 16% (corresponding to a weak health care system) and 84% (corresponding to a very strong health care system) were simulated.
It is assumed that in any real situation it would be impracticable to introduce both MDA and testandtreat programmes at the same time, and it would be illogical to commence a testandtreat programme after the MDA was completed, since resurgence might then occur before the programme that might avert it was in place. The MDA was, therefore, carried out 4 years after the testandtreat programme was initiated.
The importation rate is relevant because the intervention package needs to be able to prevent onward transmission from imported cases. To test sensitivity to importation, simulations were carried out with a value similar to a previously published estimate for Zanzibar [32].
To examine how the transmission level, the effective coverage of case management, the importation rate, and the parameters of the testandtreat programme affect the outcome, a factorially designed simulation experiment was carried out with each combination of the levels listed in Table 1 for each factor. Each such scenario was simulated with each of 3 different seeds, both with and without one round of MDA at 90% coverage, leading to a total of 16,200 population simulations. These simulations were performed at the sciCORE (http://scicore.unibas.ch/) scientific computing core facility at University of Basel.
Analysis of elimination and resurgence rate
Point prevalence for the whole population was extracted at an interval of 30 days (simulation time) throughout each simulation. Each simulation with MDA was then matched with the corresponding nonMDA simulation, and each pair of simulations was classified according to whether:
 i.
There was a significant prevalence reduction induced by the testandtreat programme alone. This was evaluated using a paired ttest to compare prevalence in the year before introduction, with that 11 years later at the same time of year. A 5% significance level was used to define statistical significance.
 ii.
Elimination occurred. This was defined as an average prevalence less than a threshold (set to 0.005 except where indicated otherwise). Sensitivity analyses were used to evaluate the effect on the classification of alternative thresholds. A criterion of zero infections would be biased because some simulations included importation and would therefore have nonzero prevalence even if local transmission was minimal.
The resurgence rate was measured by analysis of the rate of approach of the timepoint specific prevalence in the simulation with MDA, \(p_{I} \left( t \right)\), to that of controls without MDA (matched on all other aspects of the scenario), \(p_{C} \left( t \right)\). Justified by the patterns observed in simulations with resurgence, the resurgence process was approximated with an exponential function of the form:
where \(p_{\min} = p_{I} \left( 0 \right)\) denotes the minimum prevalence achieved in the year after MDA, \(\lambda\) is the resurgence rate, and \(\left( {ln2} \right)/\lambda\) is the halflife of resurgence, and where timedependent values \(p_{C} \left( t \right)\) are used as comparators with \(p_{I} \left( t \right)\) in order to allow for the seasonal forcing.
Values for the 4 years after this minimum were analysed by fitting a linear regression model to the logarithmic transform of Eq. (1), in order to estimate \(\lambda\) separately for each matched pair of simulations. Resurgence was said to occur in simulations where there was no elimination and for which the halflife was less than 10 years. Where elimination did not occur, the longterm percentage reduction in prevalence was computed from the ratio of average prevalence during years 9–11 after the round of MDA, to that in the year before the testandtreat began.
Results
Conceptual model
The analysis in the present study concerns settings where routine control measures, such as vector control and casemanagement, have been scaled up maximally but where P. falciparum malaria remains endemic. The situation initially approximates an endemic stable state, with minimal fluctuations in yearly prevalence and incidence. The intervention programme is then enhanced by adding components of testandtreat, MDA or both.
Figure 1a illustrates schematically the time course of average transmission potential in such a population for each strategy as quantified by the controlled reproduction number, R_{c}. In the absence of either testandtreat or MDA, R_{c} remains constant (the yellow line in Fig. 1a). Prevalence also remains constant (yellow lines in Fig. 1b, c). MDA and testandtreat both reduce R_{c}. MDA leads to a transient dip in R_{c} as a result of the prophylactic effect of the intervention (blue line in Fig. 1a). In addition, by clearing parasitaemia, a single round of MDA has an immediate, but transient effect on the infectious reservoir, leading to more extended effects on prevalence. Even if transmission is interrupted by the MDA, in the absence of effective surveillance reintroductions will occur and both R_{c} and prevalence will return to the initial stable states (blue lines in Fig. 1a–c) [5, 9].
The testandtreat programme is envisaged as targeted investigations around passivelydetected cases of existing infections, with resources for these investigations remaining constant over time, once the programme is initiated. The testandtreat programme is assumed to be added to the routine case management system, which serves as a passive surveillance component of a surveillanceresponse system. It functions by testing and treating neighbours of passivelydetected index cases. In Fig. 1a, the introduction of testandtreat is represented by a very rapid decrease in R_{c} (the red line). In response to the reduced value of R_{c}, prevalence decreases gradually, either until a new endemic stable state is reached (red line in Fig. 1b), or until transmission is interrupted (Fig. 1c).
While there is substantial empirical and theoretical support for the broad patterns shown for both MDA and testandtreat in Fig. 1, it is not obvious what may be achieved by combining a pulse of MDA with scaleup of testandtreat. The black lines in Fig. 1b, c illustrate schematically hypothetical trajectories for the prevalence when the two interventions are combined, with the testandtreat introduced shortly before a pulse of MDA. Table 2 lists possible results of adding MDA to a testandtreat programme (as assessed in terms of the timecourse of prevalence), with patterns C and E corresponding to the dashed lines in Fig. 1b, c, respectively, and patterns D and F, corresponding to the solid black lines in Fig. 1b, c, representing the desired outcomes. Patterns A and B correspond to cases where the final prevalence remains similar when MDA is added, to that achieved with testandtreat alone.
Equivalent to the question of whether MDA will be followed by resurgence is that of what will be the subsequent trajectory of prevalence if a system at equilibrium is disturbed by removing a substantial proportion of the infections. Figure 2 (which is derived from a very simple (susceptibleinfectedsusceptible; SIS) model of malaria dynamics in humans [17]) illustrates how the possible outcomes listed in Table 2 arise as functions of the minimum prevalence achieved after the round of MDA, and of the dynamics induced by the intervention programme. Figure 3 indicates how such systems adjust to perturbations at each prevalence value with each panel corresponds to different verticals in Fig. 2, and the vertical axis giving the rate of change (arbitrary units) over time in the prevalence, following a perturbation to the equilibrium; \(dp/dt = 0\), corresponds to an equilibrium.
In systems conforming to the pattern in Fig. 2a the prevalence will resurge to the stable equilibrium (the blue line) (see also Fig. 3a, b). The only exceptions are (i) if there is independently a change in the transmission potential, \(R_{\text{c}}\), so that the prevalence now falls into the lefthand part of the diagram, where the diseasefree equilibrium is stable, or (ii) if infections have been cleared entirely, and there are no reintroductions (in which case the system remains at the unstable, diseasefree equilibrium, corresponding to the red dashed line). If interventions are introduced that reduce the reproduction number to some new value \(R_{\text{c}}\), (corresponding to a lower endemic level) they may change the diagram into one like Fig. 2b. Qualitatively the system has not changed, but there is an extended range of settings over which elimination is occurs (corresponding to outcome E or F in Table 2).
Alternatively, the system may exhibit bistability (Fig. 2c). This occurs when there is a range of \(R_{\text{c}}\) values for which there are two stable (blue line) prevalences, separated by a nonzero unstable equilibrium at some intermediate prevalence, \(p_{uns}^{*}\). This range always has a lower bound but does not necessarily have an upper bound [17]. Bistability arises if infectiousness or persistence of infections is lower at lower prevalence and is an inherent characteristic of some mathematical models of immunity, in particular, those that include hosts that are temporarily resistant [16, 18] or with exposuredependent clearance rates [20]. In such systems, once elimination is achieved it can be maintained even in the face of low levels of introduction, providing that \(p_{uns}^{*}\) is not exceeded [19, 20]. In some such systems the bifurcation occurs in the absence of interventions, and at a value of \(R_{0} < 1 \left[ {16,18} \right]\). However simple malaria models without interventions (such as RossMacdonald, or the even simpler SIS model obtained by reducing the mosquito model to a constant factor [17]), generally give diagrams like Fig. 2a, where there is no bistability, implying that clearing parasites (e.g. with MDA) in an endemic setting, will at best lead to an unstable parasitefree state, where any reintroductions will lead to resurgence.
The endemic equilibrium in OpenMalaria (and in the microsimulation model of Griffin et al. [33]), increases monotonically with \(R_{0}\) and this has been interpreted as implying that there is no bistability [34]. However, while such nonmonotonicity does lead to bistability, it is not a necessary condition (as illustrated by Fig. 2c). Moreover, there may be a qualitative change in the bifurcation diagram from one of the topology of Fig. 2b to that of Fig. 2c depending on the interventions that are in place [17].
Results of microsimulations
Typical time profiles of prevalence in the microsimulations are shown in Fig. 4.
These were analysed to determine the minimum prevalence achieved, the average prevalence during the resurgent period, and the resurgence rates, enabling each pair of simulations to be classified according to which of the patterns listed in Table 2 they match. Each simulation could be classified into one of the six patterns, and examples of simulated timeseries of prevalence corresponding to each pattern are given in Fig. 4. Most fell into patterns A–C where there was no elimination. In 3.2% (255) of the 8100 matched pairs of simulations there was elimination in the control (patterns E and F), but only 3.0% (246) exhibited pattern D, where the MDA made a difference to whether elimination was achieved (Table 2). Even when there was elimination defined using the parasite density cutoff, the simulated EIR (during the 4 years after the MDA) was nonzero (Table 3), though at values of less than 0.1 infectious bites per personyear. These EIR levels are clearly too low to sustain transmission and would be unmeasurable in the field, but do lead to small numbers of secondary cases, which can explain the seasonal oscillations in prevalence in simulations classified as having achieved elimination.
In total 88.3% (7150 of 8100) of simulations with MDA were classified as having resurgence with a median halflife of resurgence of 2.3 years (interquartile range 1.1–2.8). When elimination was not achieved, the MDA had minimal effect on the percentage reduction in final prevalence (median reduction 0.7%, interquartile range − 2%, 5%).
Because of the factorial design, the marginal averages are unbiased estimates of the effects of the individual factors that varied by design. The parameters of the testandtreat programme rarely affected whether elimination occurred (Fig. 5a–c), but had a large effect on the proportion of simulation pairs exhibiting patterns A or B, indicating that intensive testandtreat changes the steadystate prevalence, when sustained over a long period (in this case of approximately a decade). Even in the simulations without MDA there was a substantial reduction in the EIR (from the median value of 5 infectious bites per personyear), resulting from the effect of the testandtreat (Table 3).
The greater the intensity of testandtreat, the greater the longterm effect on the prevalence (as assessed by the approximate steady state in years 9–11 after MDA). At the same time, the parameters of the testandtreat programme had little effect on the resurgence rate in those simulations where there was resurgence (Fig. 6a–c).
Case management coverage remained the same throughout the simulation (it was not scaledup, or introduced partway through the programme, like the testandtreat intervention), so its consequences result from interactions with the other interventions.
These included a very dramatic effect on the proportion of simulations leading to elimination (Fig. 5d). Elimination occurred in none of the simulations at the low case management coverage. At the same time, only 23% (832) of noneliminating MDA simulations with high case management coverage had pattern A (where the final state was not significantly different from the preMDA one) compared with, 32% (1292) of the corresponding simulations with low case management coverage, indicating that the testandtreat strategy is more likely to lead to an accumulated effect on prevalence if routine case management is at high coverage. Related to this, the resurgence rate was lower in those simulations with high case management coverage (Fig. 6d).
The transmission potential, as measured by the initial EIR, also had a dramatic effect on the percentage of simulations leading to elimination (Fig. 5e). 15% (485) of simulations of initial EIR ≤ 4 exhibited elimination, but in only 16/4860 (0.3%) of simulations with higher initial EIR was there elimination. Twelve of these exhibited pattern D (elimination in the control), and only 4 had pattern E or pattern F, indicating that the MDA had minimal effect on the longterm outcome if the EIR exceeded this rather modest value. In the simulations where elimination occurred with MDA, the control simulations were ones with very low control EIR during the postMDA period (Table 3). When elimination did not occur, the longterm effect of the testandtreat programme on prevalence was greater in low (initial) EIR simulations (Fig. 7) and in these simulations the resurgence rate was slower than with a higher EIR (Fig. 6e).
Imported infections decreased the proportion of simulations leading to elimination (Fig. 5f), but did not completely avert elimination when case management coverage was high. Importation, however made no difference to the final steady state in the nonelimination scenarios, and hence did not affect the proportion of simulations with pattern A, or the rate of resurgence when there was no elimination (Fig. 6f).
It is impracticable to draw bifurcation diagrams for the OpenMalaria model, because it is both computationally expensive and stochastic, but the conceptual model (Fig. 2) can account for the occurrence of the different outcome patterns (Fig. 5) in terms of the impact of the interventions on transmission potential, and whether they induce bistability. Where transmission is interrupted in simulations with MDA, but not without MDA this implies that the system is bistable. Since this occurs only in simulations with high casemanagement coverage, it implies that high case management coverage can induce bistability.
Where the treatment rate (considering both treatment of clinical cases and of infections diagnosed via testandtreat strategies), increases less than proportionately with prevalence, the system is bistable. A simple theoretical analysis (Additional file 1) indicates why this arises when the effective coverage of case management is high, even when there is no bistability in the model without case management. The same analysis indicates that testandtreat (as simulated here) reduces transmission by greater proportion at higher levels of prevalence, and so does not induce bistability.
In the simulations, this differential effect of high coverage of case management, compared with the testandtreat programme is reinforced by the relationship between simulated treatment rates and the EIR during the period after MDA (Table 3).
Following MDA there were substantial declines in both the EIR and the rates of case management below the levels in the control simulations. This was found in all outcome groups. However, the decrease in in the numbers of treatments of clinical cases was less than that in the EIR so the ratio increased in the scenarios where transmission was reduced to the lowest levels (Fig. 8a), generally reaching the highest levels in the elimination scenarios with patterns D–F (Table 3) (though at the very lowest levels of residual EIR there was considerable stochasticity in the number of treatments). At low case management coverage (Fig. 8b), the residual EIR never fell into the range where this ratio was substantial.
The number of simulated treatments resulting from the testandtreat programme did not adjust in the same way to the residual EIR. At very low EIR, (and hence low caseincidence) the effectiveness of the testandtreat is limited by the number of index cases, i.e. \(\iota < \iota_{\text{max}}\). At higher residual EIR the effective screening rate \(\Phi = \Phi_{\text{max} } = \frac{{\iota_{\text{max} } \nu \tau }}{N}\) is constant within any given scenario, and is independent of transmission, so the consequent rate of treatment is proportional to the prevalence (Table 3, Fig. 8) (and hence the ratio to the EIR decreases somewhat with increasing EIR, owing to saturation).
Discussion
Field trials of malaria elimination strategies including multiple interventions are challenging undertakings, and it is impossible to trial all the different strategies that might be considered. In silico trials using microsimulation models such as OpenMalaria make it possible to test a wide range of different strategies, calibrated to approximate real settings, so that the analyses allow in a quantitative way for complicating factors such as seasonality, acquired immunity, and case importation. In the simulations reported here, MDA led to elimination in only a small proportion of the simulations, and in most scenarios there was resurgence, despite the very high coverage of the intervention. Testandtreat strategies also only contributed to elimination in only a small number of the simulated scenarios.
Concepts derived from simpler models, such as the RossMacdonald [35] model, or even SIS compartment models that reduce the complexity of the system to a minimum [17] are invaluable for understanding why these rather disappointing results were obtained, and for deducing under which circumstances success might be achieved.
The impact of MDA was simulated as equivalent to that of a single round of MDA with high (90%) coverage, a rather optimistic coverage level. In principle, MDA programmes can consist of any number of rounds, each potentially covering different overlapping subsets of the population. However, previous OpenMalaria simulations suggested that the impact can mostly be explained by the escape probability, defined as the proportion of the population that never receive any MDA [26]. Other factors compromising coverage (such as seasonal migration [36]) could also be captured by effect on the escape probability. A low escape probability (10%) thus errs on the side of the optimistic, while approximating many different patterns of coverage of multiple MDA rounds.
The small proportion of the scenarios in which elimination occurred is readily explained on the basis of the inferred form of the bifurcation diagram. MDA is a pulsed intervention with only a transient effect on \(R_{\text{c}}\) (which arises because of the prophylactic effect and the truncation of existing infections). Only in a subset of settings is there bistability and only in a subset of these does MDA move the system abruptly below \(p_{uns}^{ *}\) and into the dark green zone in Fig. 2c and thus achieve elimination. This is consistent with simple stochastic models for the probability of extinction by MDA, which indicate that this is extremely unlikely unless the population size is very small and the coverage of MDA extremely high [37].
The existence of the bistability that enables elimination depends crucially on the coverage of case management. In various models of clinical malaria, the ratio of the incidence of clinical attacks to the number of inoculations is higher in lower prevalence settings. If there is significant treatment in reaction to these attacks, this leads to the pattern of bifurcation diagram of Fig. 2c, with a saddlenode bifurcation at the lowest value of \(R_{\text{c}}\) at which transmission persists. In such systems, elimination following MDA may be a stable state, even if it is not achieved without MDA. This can arise if the health system functions better when malaria is nearabsent or when prevalence is very low [19] or if new infections are more likely to be treated in low prevalence settings where there is less clinical immunity [19]. With OpenMalaria this pattern is induced by highly effective case management as a result of the selection effect described in Additional file 1, which requires neither of these preconditions. At values of \(R_{\text{c}}\) higher than that at the saddlenode there are two possible nonzero equilibrium prevalence values (Fig. 2c). The lower branch of unstable endemic equilibrium points divides the basins of attraction of the diseasefree equilibrium point and the larger endemic equilibrium points. If the prevalence is greater than the unstable endemic equilibrium point, then the prevalence will resurge until it reaches the larger endemic equilibrium point. If the prevalence after MDA is less than this unstable equilibrium point, \(p_{uns}^{*}\), then the infection dies out and case management will lead to elimination.
At the simulated importation rate and coverages of case management, imported infections did not affect whether resurgence occurred, indicating that the case management alone was sufficient to control importations in those cases where MDA had brought prevalence below \(p_{uns}^{*}\). If prevalence is reduced to zero, then importation is necessary for the system to escape from an unstable diseasefree equilibrium, but this is a highly unlikely scenario since the diseasefree equilibrium is only achieved when there is a high coverage of case management. Since importation has minimal effects on resurgence rates, it is unlikely to be an important determinant of outcomes in such settings.
In contrast, the simulations reported here agree that the transmission potential is both the primary determinant of whether elimination can be achieved by a malaria intervention programme, and also of the rate of resurgence after MDA when it cannot. Because of its modest impact on \(R_{c}\), the coverage of case management has only a small effect on the resurgence rate, despite its critical role in determining whether resurgence occurs at all. If coverage of case management is low, reducing transmission leads to diagrams like that of Fig. 2b, where the threshold value of \(R_{0}\) (where the endemic equilibrium meets the axis) is increased, but where there is no bistability. With this form of bifurcation diagram, MDA will not increase the range of settings where elimination occurs.
Testandtreat strategies do not lead to bistability in the scenarios shown in this paper, and only contribute to elimination via their modest reduction in average transmission. Testandtreat accelerates the approach to elimination only in a subset of those simulations where the testandtreat programme already rapidly reduces prevalence. About the same number of simulations were classified into patterns F (elimination with no resurgence, arguably equivalent to acceleration) and E (elimination following temporary resurgence), but in all these cases, the prevalence in both control and MDA simulations was already very low by the time of the MDA (Fig. 4).
Untargeted testandtreat is an inefficient way of reducing average transmission because its coverage is linearly related to the reproduction number (in contrast to killing adult female mosquitoes which is related via an exponential function of coverage [38]), so it is unsurprising that testandtreat leads to elimination in only a small subset of simulations. In principle, effective targeting can substantially improve the impact of testandtreat. This may also lead to bistability if the targeting ratio, \(\uptau\), decreases with prevalence [17]. The testandtreat simulated here in OpenMalaria treats \(\tau\) as fixed (and also assumes that the resources allocated to the testandtreat system, measured by the product of \(\iota\) and \(\nu\), are constant). These are likely to be reasonable approximations to recurrent focal screenandtreat or recurrent focal MDA, especially where targeting in space and time is rather imprecise, but give a poor representation of highly targeted reactive case detection. There is a need to evaluate microsimulations of reactive case detection approaches where \(\uptau\) may adapt over time [17].
Conclusions
The main determinants of whether resurgence will occur following MDA are the residual transmission potential and the routine treatment of passively detected cases. A high coverage of case management is essential for elimination, but in most simulated scenarios, transmission was not interrupted and the transmission potential proved the main determinant of the resurgence rate. Other reactive interventions, such as high intensity reactive case detection and reactive vector control might also be effective in averting resurgence at low transmission potential and require further analysis.
Availability of data and materials
The underlying code and simulation results are available on request from the authors.
Abbreviations
 MDA:

Mass Drug Administration
References
 1.
Bhatt S, Weiss DJ, Cameron E, Bisanzio D, Mappin B, Dalrymple U, et al. The effect of malaria control on Plasmodium falciparum in Africa between 2000 and 2015. Nature. 2015;526:207–11.
 2.
Larsen DA, Bennett A, Silumbe K, Hamainza B, Yukich JO, Keating J, et al. Populationwide malaria testing and treatment with rapid diagnostic tests and artemetherlumefantrine in southern Zambia: a community randomized stepwedge control trial design. Am J Trop Med Hyg. 2015;92:913–21.
 3.
Eisele TP, Silumbe K, Finn T, Chalwe V, Kamuliwo M, Hamainza B, et al. Assessing the effectiveness of householdlevel focal mass drug administration and communitywide mass drug administration for reducing malaria parasite infection prevalence and incidence in Southern Province, Zambia: study protocol for a community randomized controlled trial. Trials. 2015;16:347.
 4.
Eisele TP, Bennett A, Silumbe K, Finn TP, Chalwe V, Kamuliwo M, et al. Shortterm impact of mass drug administration with dihydroartemisinin plus piperaquine on malaria in Southern Province Zambia: a clusterrandomized controlled trial. J Infect Dis. 2016;214:1831–9.
 5.
Smith T, Schapira A. Reproduction numbers in malaria and their implications. Trends Parasitol. 2012;28:3–8.
 6.
Pampana E. A textbook of malaria eradication. London: Oxford; 1963.
 7.
Macdonald G, Göckel GW. The malaria parasite rate and interruption of transmission. Bull World Health Organ. 1964;31:365–77.
 8.
Greenwood BM, David PH, OtooForbes LN, Allen SJ, Alonso PL, ArmstrongSchellenberg JR, et al. Mortality and morbidity from malaria after stopping malaria chemoprophylaxis. Trans R Soc Trop Med Hyg. 1995;89:629–33.
 9.
von Seidlein L, Greenwood BM. Mass administrations of antimalarial drugs. Trends Parasitol. 2003;19:452–60.
 10.
Greenwood B. The use of antimalarial drugs to prevent malaria in the population of malariaendemic areas. Am J Trop Med Hyg. 2004;70:1–7.
 11.
BruceChwatt LJ. Resurgence of malaria and its control. J Trop Med Hyg. 1974;77:s6.
 12.
Cohen JM, Smith DL, Cotter C, Ward A, Yamey G, Sabot OJ, et al. Malaria resurgence: a systematic review and assessment of its causes. Malar J. 2012;11:122.
 13.
Stuckey EM, Miller JM, Littrell M, Chitnis N, Steketee R. Operational strategies of antimalarial drug campaigns for malaria elimination in Zambia’s southern province: a simulation study. Malar J. 2016;15:148.
 14.
Yukich J, Bennett A, Yukich R, Stuck L, Hamainza B, Silumbe K, et al. Estimation of malaria parasite reservoir coverage using reactive case detection and active community fever screening from census data in Southern Zambia: a resampling approach. Malar J. 2017;16:317.
 15.
Smith DL, Smith T, Hay SI. Measuring malaria for elimination. In: Feachem R, Phillips A, Targett GA, editors. Shrinking the malaria map: a prospectus on malaria elimination. San Francisco: The Global Health Group; 2009. p. 108–26.
 16.
Chitnis N, Cushing JM, Hyman JM. Bifurcation analysis of a mathematical model for malaria transmission. SIAM J Appl Math. 2006;67:24–45.
 17.
Chitnis N, PembertonRoss P, Yukich J, Hamainza B, Miller J, Reiker T, et al. Theory of reactive interventions in the elimination and control of malaria. Malar J. 2019;18:266.
 18.
Aguas R, White LJ, Snow RW, Gomes MG. Prospects for malaria eradication in subSaharan Africa. PLoS ONE. 2008;3:e1767.
 19.
Smith DL, Cohen JM, Chiyaka C, Johnston G, Gething PW, Gosling R, et al. A sticky situation: the unexpected stability of malaria elimination. Philos Trans R Soc Lond B Biol Sci. 2013;368:20120145.
 20.
Keegan LT, Dushoff J. Populationlevel effects of clinical immunity to malaria. BMC Infect Dis. 2013;13:428.
 21.
Smith T, Killeen GF, Maire N, Ross A, Molineaux L, Tediosi F, et al. 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.
 22.
Chitnis N, Hardy D, Smith T. A periodicallyforced mathematical model for the seasonal dynamics of malaria in mosquitoes. Bull Math Biol. 2012;74:1098–124.
 23.
Smith T, Maire N, Dietz K, Killeen GF, Vounatsou P, Molineaux L, et al. Relationships between the entomological inoculation rate and the force of infection for Plasmodium falciparum malaria. Am J Trop Med Hyg. 2006;75(Suppl 2):11–8.
 24.
Maire N, Smith T, Ross A, OwusuAgyei 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.
 25.
Smith T, Ross A, Maire N, Rogier C, Trape JF, Molineaux L. An epidemiologic model of the incidence of acute illness in Plasmodium falciparum malaria. Am J Trop Med Hyg. 2006;75:56–62.
 26.
Brady OJ, Slater HC, PembertonRoss P, Wenger E, Maude RJ, Ghani AC, et al. Role of mass drug administration in elimination of Plasmodium falciparum malaria: a consensus modelling study. Lancet Global Health. 2017;5:E680–7.
 27.
Crowell V, Yukich JO, Briet O, Ross A, Smith T. A novel approach for measuring the burden of uncomplicated Plasmodium falciparum malaria: application to data from Zambia. PLoS ONE. 2013;8:e57297.
 28.
Murray CK, Gasser RA Jr, Magill AJ, Miller RS. Update on rapid diagnostic testing for malaria. Clin Microbiol Rev. 2008;21:97–110.
 29.
Reiker T, Chitnis N, Smith T. Modelling reactive case detection strategies for interrupting transmission of Plasmodium falciparum malaria. Malar J. 2019;18:259.
 30.
Smith T, Ross A, Maire N, Chitnis N, Studer A, Hardy D, et al. Ensemble modeling of the likely public health impact of a preerythrocytic malaria vaccine. PLoS Med. 2012;9:e1001157.
 31.
Galactionova K, Tediosi F, De Savigny D, Smith T, Tanner M. Effective coverage and systems effectiveness for malaria case management in subSaharan African countries. PLoS ONE. 2015;10:e0127818.
 32.
Tatem AJ, Qiu Y, Smith DL, Sabot O, Ali AS, Moonen B. The use of mobile phone data for the estimation of the travel patterns and imported Plasmodium falciparum rates among Zanzibar residents. Malar J. 2009;8:287.
 33.
Griffin JT, Hollingsworth TD, Okell LC, Churcher TS, White M, Hinsley W, et al. Reducing Plasmodium falciparum malaria transmission in Africa: a modelbased evaluation of intervention strategies. PLoS Med. 2010;7:e1000324.
 34.
Griffin JT. Is a reproduction number of one a threshold for Plasmodium falciparum malaria elimination? Malar J. 2016;15:389.
 35.
Smith DL, Battle KE, Hay SI, Barker CM, Scott TW, McKenzie FE. Ross, Macdonald, and a theory for the dynamics and control of mosquitotransmitted pathogens. PLoS Pathog. 2012;8:e1002588.
 36.
Gerardin J, BertozziVilla A, Eckhoff PA, Wenger EA. Impact of mass drug administration campaigns depends on interaction with seasonal human movement. Int Health. 2018;10:252–7.
 37.
PembertonRoss P, Chitnis N, Pothin E, Smith TA. A stochastic model for the probability of malaria extinction by mass drug administration. Malar J. 2017;16:376.
 38.
Macdonald G. Theory of the eradication of malaria. Bull World Health Organ. 1956;15:369–87.
Acknowledgements
Not applicable.
Funding
TS, PPR, and NC acknowledge funding by the Bill and Melinda Gates Foundation (#OPP1032350). MAP is funded by the Swiss National Science Foundation through SNSF Professorship PP00P3_170702. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Author information
Affiliations
Contributions
TAS, PPR, MAP, and NC jointly conceived of the study, and designed the experiments. PPR and TAS analysed the results. TAS drafted the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
About this article
Cite this article
Smith, T.A., PembertonRoss, P., Penny, M.A. et al. Resurgence of malaria infection after mass treatment: a simulation study. Malar J 18, 409 (2019). https://doi.org/10.1186/s1293601930190
Received:
Accepted:
Published:
Keywords
 Plasmodium falciparum
 Simulation
 Elimination