Modelling the cost-effectiveness of mass screening and treatment for reducing Plasmodium falciparum malaria burden.

BACKGROUND: Past experience and modelling suggest that, in most cases, mass treatment strategies are not likely to succeed in interrupting Plasmodium falciparum malaria transmission. However, this does not preclude their use to reduce disease burden. Mass screening and treatment (MSAT) is preferred to mass drug administration (MDA), as the latter involves massive over-use of drugs. This paper reports simulations of the incremental cost-effectiveness of well-conducted MSAT campaigns as a strategy for P. falciparum malaria disease-burden reduction in settings with varying receptivity (ability of the combined vector population in a setting to transmit disease) and access to case management. METHODS: MSAT incremental cost-effectiveness ratios (ICERs) were estimated in different sub-Saharan African settings using simulation models of the dynamics of malaria and a literature-based MSAT cost estimate. Imported infections were simulated at a rate of two per 1,000 population per annum. These estimates were compared to the ICERs of scaling up case management or insecticide-treated net (ITN) coverage in each baseline health system, in the absence of MSAT. RESULTS: MSAT averted most episodes, and resulted in the lowest ICERs, in settings with a moderate level of disease burden. At a low pre-intervention entomological inoculation rate (EIR) of two infectious bites per adult per annum (IBPAPA) MSAT was never more cost-effective than scaling up ITNs or case management coverage. However, at pre-intervention entomological inoculation rates (EIRs) of 20 and 50 IBPAPA and ITN coverage levels of 40 or 60%, respectively, the ICER of MSAT was similar to that of scaling up ITN coverage further. CONCLUSIONS: In all the transmission settings considered, achieving a minimal level of ITN coverage is a "best buy". At low transmission, MSAT probably is not worth considering. Instead, MSAT may be suitable at medium to high levels of transmission and at moderate ITN coverage. If undertaken as a burden-reducing intervention, MSAT should be continued indefinitely and should complement, not replace, case management and vector control interventions.

Effects of pyrethroid resistance on the cost effectiveness of a mass distribution of long-lasting insecticidal nets: a modelling study.

BACKGROUND: The effectiveness of insecticide-treated nets in preventing malaria is threatened by developing resistance against pyrethroids. Little is known about how strongly this affects the effectiveness of vector control programmes. METHODS: Data from experimental hut studies on the effects of long-lasting, insecticidal nets (LLINs) on nine anopheline mosquito populations, with varying levels of mortality in World Health Organization susceptibility tests, were used to parameterize malaria models. Both simple static models predicting population-level insecticidal effectiveness and protection against blood feeding, and complex dynamic epidemiological models, where LLINs decayed over time, were used. The epidemiological models, implemented in OpenMalaria, were employed to study the impact of a single mass distribution of LLINs on malaria, both in terms of episodes prevented during the effective lifetime of the batch of LLINs, and in terms of net health benefits (NHB) expressed in disability-adjusted life years (DALYs) averted during that period, depending on net type (standard pyrethroid-only LLIN or pyrethroid-piperonyl butoxide combination LLIN), resistance status, coverage and pre-intervention transmission level. RESULTS: There were strong positive correlations between insecticide susceptibility status and predicted population level insecticidal effectiveness of and protection against blood feeding by LLIN intervention programmes. With the most resistant mosquito population, the LLIN mass distribution averted up to about 40% fewer episodes and DALYs during the effective lifetime of the batch than with fully susceptible populations. However, cost effectiveness of LLINs was more sensitive to the pre-intervention transmission level and coverage than to susceptibility status. For four out of the six Anopheles gambiae sensu lato populations where direct comparisons between standard LLINs and combination LLINs were possible, combination nets were more cost effective, despite being more expensive. With one resistant population, both net types were equally effective, and with one of the two susceptible populations, standard LLINs were more cost effective. CONCLUSION: Despite being less effective when compared to areas with susceptible mosquito populations, standard and combination LLINs are likely to (still) be cost effective against malaria even in areas with strong pyrethroid resistance. Combination nets are likely to be more cost effective than standard nets in areas with resistant mosquito populations.

Ensemble modeling of the likely public health impact of a pre-erythrocytic malaria vaccine.

BACKGROUND: The RTS,S malaria vaccine may soon be licensed. Models of impact of such vaccines have mainly considered deployment via the World Health Organization's Expanded Programme on Immunization (EPI) in areas of stable endemic transmission of Plasmodium falciparum, and have been calibrated for such settings. Their applicability to low transmission settings is unclear. Evaluations of the efficiency of different deployment strategies in diverse settings should consider uncertainties in model structure. METHODS AND FINDINGS: An ensemble of 14 individual-based stochastic simulation models of P. falciparum dynamics, with differing assumptions about immune decay, transmission heterogeneity, and treatment access, was constructed. After fitting to an extensive library of field data, each model was used to predict the likely health benefits of RTS,S deployment, via EPI (with or without catch-up vaccinations), supplementary vaccination of school-age children, or mass vaccination every 5 y. Settings with seasonally varying transmission, with overall pre-intervention entomological inoculation rates (EIRs) of two, 11, and 20 infectious bites per person per annum, were considered. Predicted benefits of EPI vaccination programs over the simulated 14-y time horizon were dependent on duration of protection. Nevertheless, EPI strategies (with an initial catch-up phase) averted the most deaths per dose at the higher EIRs, although model uncertainty increased with EIR. At two infectious bites per person per annum, mass vaccination strategies substantially reduced transmission, leading to much greater health effects per dose, even at modest coverage. CONCLUSIONS: In higher transmission settings, EPI strategies will be most efficient, but vaccination additional to the EPI in targeted low transmission settings, even at modest coverage, might be more efficient than national-level vaccination of infants. The feasibility and economics of mass vaccination, and the circumstances under which vaccination will avert epidemics, remain unclear. The approach of using an ensemble of models provides more secure conclusions than a single-model approach, and suggests greater confidence in predictions of health effects for lower transmission settings than for higher ones.

Importance of factors determining the effective lifetime of a mass, long-lasting, insecticidal net distribution: a sensitivity analysis.

BACKGROUND: Long-lasting insecticidal nets (LLINs) reduce malaria transmission by protecting individuals from infectious bites, and by reducing mosquito survival. In recent years, millions of LLINs have been distributed across sub-Saharan Africa (SSA). Over time, LLINs decay physically and chemically and are destroyed, making repeated interventions necessary to prevent a resurgence of malaria. Because its effects on transmission are important (more so than the effects of individual protection), estimates of the lifetime of mass distribution rounds should be based on the effective length of epidemiological protection. METHODS: Simulation models, parameterised using available field data, were used to analyse how the distribution's effective lifetime depends on the transmission setting and on LLIN characteristics. Factors considered were the pre-intervention transmission level, initial coverage, net attrition, and both physical and chemical decay. An ensemble of 14 stochastic individual-based model variants for malaria in humans was used, combined with a deterministic model for malaria in mosquitoes. RESULTS: The effective lifetime was most sensitive to the pre-intervention transmission level, with a lifetime of almost 10 years at an entomological inoculation rate of two infectious bites per adult per annum (ibpapa), but of little more than 2 years at 256 ibpapa. The LLIN attrition rate and the insecticide decay rate were the next most important parameters. The lifetime was surprisingly insensitive to physical decay parameters, but this could change as physical integrity gains importance with the emergence and spread of pyrethroid resistance. CONCLUSIONS: The strong dependency of the effective lifetime on the pre-intervention transmission level indicated that the required distribution frequency may vary more with the local entomological situation than with LLIN quality or the characteristics of the distribution system. This highlights the need for malaria monitoring both before and during intervention programmes, particularly since there are likely to be strong variations between years and over short distances. The majority of SSA's population falls into exposure categories where the lifetime is relatively long, but because exposure estimates are highly uncertain, it is necessary to consider subsequent interventions before the end of the expected effective lifetime based on an imprecise transmission measure.

Simulation of malaria epidemiology and control in the highlands of Western Kenya.

BACKGROUND: Models of Plasmodium falciparum malaria epidemiology that provide realistic quantitative predictions of likely epidemiological outcomes of existing vector control strategies have the potential to assist in planning for the control and elimination of malaria. This work investigates the applicability of mathematical modelling of malaria transmission dynamics in Rachuonyo South, a district with low, unstable transmission in the highlands of western Kenya. METHODS: Individual-based stochastic simulation models of malaria in humans and a deterministic model of malaria in mosquitoes as part of the OpenMalaria platform were parameterized to create a scenario for the study area based on data from ongoing field studies and available literature. The scenario was simulated for a period of two years with a population of 10,000 individuals and validated against malaria survey data from Rachuonyo South. Simulations were repeated with multiple random seeds and an ensemble of 14 model variants to address stochasticity and model uncertainty. A one-dimensional sensitivity analysis was conducted to address parameter uncertainty. RESULTS: The scenario was able to reproduce the seasonal pattern of the entomological inoculation rate (EIR) and patent infections observed in an all-age cohort of individuals sampled monthly for one year. Using an EIR estimated from serology to parameterize the scenario resulted in a closer fit to parasite prevalence than an EIR estimated using entomological methods. The scenario parameterization was most sensitive to changes in the timing and effectiveness of indoor residual spraying (IRS) and the method used to detect P. falciparum in humans. It was less sensitive than expected to changes in vector biting behaviour and climatic patterns. CONCLUSIONS: The OpenMalaria model of P. falciparum transmission can be used to simulate the impact of different combinations of current and potential control interventions to help plan malaria control in this low transmission setting. In this setting and for these scenarios, results were highly sensitive to transmission, vector exophagy, exophily and susceptibility to IRS, and the detection method used for surveillance. The level of accuracy of the results will thus depend upon the precision of estimates for each. New methods for analysing and evaluating uncertainty in simulation results will enhance the usefulness of simulations for malaria control decision-making. Improved measurement tools and increased primary data collection will enhance model parameterization and epidemiological monitoring. Further research is needed on the relationship between malaria indices to identify the best way to quantify transmission in low transmission settings. Measuring EIR through mosquito collection may not be the optimal way to estimate transmission intensity in areas with low, unstable transmission.

A periodically-forced mathematical model for the seasonal dynamics of malaria in mosquitoes.

We describe and analyze a periodically-forced difference equation model for malaria in mosquitoes that captures the effects of seasonality and allows the mosquitoes to feed on a heterogeneous population of hosts. We numerically show the existence of a unique globally asymptotically stable periodic orbit and calculate periodic orbits of field-measurable quantities that measure malaria transmission. We integrate this model with an individual-based stochastic simulation model for malaria in humans to compare the effects of insecticide-treated nets (ITNs) and indoor residual spraying (IRS) in reducing malaria transmission, prevalence, and incidence. We show that ITNs are more effective than IRS in reducing transmission and prevalence though IRS would achieve its maximal effects within 2 years while ITNs would need two mass distribution campaigns over several years to do so. Furthermore, the combination of both interventions is more effective than either intervention alone. However, although these interventions reduce transmission and prevalence, they can lead to increased clinical malaria; and all three malaria indicators return to preintervention levels within 3 years after the interventions are withdrawn.

Incorporating genetic selection into individual-based models of malaria and other infectious diseases.

INTRODUCTION: Control strategies for human infections are often investigated using individual-based models (IBMs) to quantify their impact in terms of mortality, morbidity and impact on transmission. Genetic selection can be incorporated into the IBMs to track the spread of mutations whose origin and spread are driven by the intervention and which subsequently undermine the control strategy; typical examples are mutations which encode drug resistance or diagnosis- or vaccine-escape phenotypes. METHODS AND RESULTS: We simulated the spread of malaria drug resistance using the IBM OpenMalaria to investigate how the finite sizes of IBMs require strategies to optimally incorporate genetic selection. We make four recommendations. Firstly, calculate and report the selection coefficients, s, of the advantageous allele as the key genetic parameter. Secondly, use these values of "s" to calculate the wait time until a mutation successfully establishes itself in the pathogen population. Thirdly, identify the inherent limits of the IBM to robustly estimate small selection coefficients. Fourthly, optimize computational efficacy: when "s" is small, fewer replicates of larger IBMs may be more efficient than a larger number of replicates of smaller size. DISCUSSION: The OpenMalaria IBM of malaria was an exemplar and the same principles apply to IBMs of other diseases.

Can we depend on case management to prevent re-establishment of P. falciparum malaria, after local interruption of transmission?

Recent declines in malaria burden in many parts of the world have prompted consideration of how interruption of Plasmodium falciparum transmission could be maintained, if achieved, and notably whether large-scale vector control could be replaced with surveillance. This information is essential for elimination feasibility assessments and planning. The risk of re-establishment of transmission depends mainly on vectorial capacity (receptivity), likely to rebound once vector control is removed, the rate of importation of infections (vulnerability), the capacity to detect and treat infections and the level of immunity in infected individuals. Timely detection and removal of new infections is likely to be critical to prevent re-establishment of transmission. We assess, through mathematical modeling and simulation, which levels of case detection and treatment (case management) are required to prevent re-establishment of transmission of P. falciparum after local interruption of transmission has been achieved, in settings with varying receptivity and vulnerability. We find that, even at rather low levels of receptivity, case management alone cannot reliably prevent re-establishment of P. falciparum malaria transmission in the face of medium to high importation rates. Thus, if vector control is to be discontinued, preventing the importations by controlling transmission in source areas will generally be necessary for preventing reintroduction in such settings, and cannot be substituted by very high levels of case management coverage.

Simulated impact of RTS,S/AS01 vaccination programs in the context of changing malaria transmission.

INTRODUCTION: The RTS,S/AS01 pre-erythrocytic malaria vaccine is in phase III clinical trials. It is critical to anticipate where and how it should be implemented if trials are successful. Such planning may be complicated by changing levels of malaria transmission. METHODS/RESULTS: Computer simulations were used to examine RTS,S/AS01 impact, using a vaccine profile based on phase II trial results, and assuming that protection decays only slowly. Settings were simulated in which baseline transmission (in the absence of vaccine) was fixed or varied between 2 and 20 infectious mosquito bites per person per annum (ibpa) over ten years. Four delivery strategies were studied: routine infant immunization (EPI), EPI plus infant catch-up, EPI plus school-based campaigns, and EPI plus mass campaigns. Impacts in changing transmission settings were similar to those in fixed settings. Assuming a persistent effect of vaccination, at 2 ibpa, the vaccine averted approximately 5-7 deaths per 1000 doses of vaccine when delivered via mass campaigns, but the benefit was less at higher transmission levels. EPI, catch-up and school-based strategies averted 2-3 deaths per 1000 doses in settings with 2 ibpa. In settings where transmission was decreasing or increasing, EPI, catch-up and school-based strategies averted approximately 3-4 deaths per 1000 doses. DISCUSSION: Where transmission is changing, it appears to be sufficient to consider simulations of pre-erythrocytic vaccine impact at a range of initial transmission levels. At 2 ibpa, mass campaigns averted the most deaths and reduced transmission, but this requires further study. If delivered via EPI, RTS,S/AS01 could avert approximately 6-11 deaths per 1000 vaccinees in all examined settings, similar to estimates for pneumococcal conjugate vaccine in African infants. These results support RTS,S/AS01 implementation via EPI, for example alongside vector control interventions, providing that the phase III trials provide support for our assumptions about efficacy.