Importance of factors determining the effective lifetime of a mass, long-lasting, insecticidal net distribution: a sensitivity analysis
© Briët et al; licensee BioMed Central Ltd. 2012
Received: 15 November 2011
Accepted: 13 January 2012
Published: 13 January 2012
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.
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.
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.
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.
Over the period 2008-2010, an estimated 290 million long-lasting insecticidal nets (LLINs) were distributed in sub-Saharan Africa . LLINs reduce malaria transmission by protecting individuals from infectious bites, and by reducing the probability that a mosquito survives the extrinsic incubation period. Whereas continuous distribution through antenatal clinics is common, most LLINs are being distributed through mass campaigns, reaching a large proportion of the population at risk of malaria. Over time, after a mass distribution, the proportion of the population sleeping under an LLIN decreases. This is partly due to attrition (the loss of nets available for their intended use, e.g. by alternative use), but also due to new births and user fatigue adding to the unprotected population.
The effective protection of LLINs against mosquito bites also wanes as they decay physically (hole formation) and chemically (insecticide loss).
Many LLIN programmes work with the assumptions that there is little variability in the decay among nets and that they last about 3 years, at which time they need replacement. However, variability in net decay appears to be substantial and the average 'lifespan' could be considerably less than 3 years . The number of LLINs remaining in households does not take the physical and chemical state of the nets into account, and a proportion of those nets may have lost considerable functionality . The World Health Organization (WHO) recommends tracking the physical integrity of nets (number, size and location of holes) and the insecticidal activity, measured by knock down and killing in standard WHO cone and tunnel tests . Unfortunately, little is known about how these quantities, alone or in interaction, affect personal protection, and how they could be used to define when a net is worn out and at the end of its 'useful life'.
Even if individual nets do not adequately prevent mosquitoes from inoculating the user, they may still reduce mosquito survival and thus affect transmission at the population level. This community effect is likely to be more important than personal protection in preventing inoculations [4, 5]. The timing of repeat LLIN distributions may also depend on the characteristics of the human population, in particular the transmission level and immune status. Whereas knowing the 'useful life' of individual LLINs (for which a cut-off minimum functionality would need to be defined, below which an LLIN would be declared 'dead') might facilitate planning in continuous distribution programmes. For a round of mass distributed LLINs, the 'effective lifetime', based on the duration of the malaria preventive effect at population level capturing all the effects described above, might be more useful for planning the timing of subsequent rounds.
This paper describes a simulation experiment to predict the duration of epidemiological protection offered by a mass LLIN distribution targeting the general population and identifies the factors that are important in determining it.
The OpenMalaria modelling platform  is an open source C++ programme and takes scenario specification inputs in eXtensible Markup Language (XML). In this platform, stochastic individual-based models for malaria in humans are combined with a deterministic model for malaria in mosquitoes, which have been fitted to multiple field data sets . For each five-day time step, data on a human population is updated via components representing new infections, parasite densities, acquired immunity, uncomplicated and severe malaria episodes, direct and indirect mortality, infectiousness to mosquitoes, and case management. Each simulated malaria infection has a distinct parasite density that varies by time step, while the malaria transmission level varies seasonally. The models can accommodate multiple mosquito species with varying periodical emergence rates, and non-human hosts [8, 9]. An ensemble of 14 model variants  is currently available, capturing a range of possibilities for the dynamics of malaria in humans.
For this experiment, the existing insecticide treated net (ITN) intervention model component [8, 9] was developed to include capability to model physical and chemical decay of LLINs. The effect of LLINs, depending on their physical and chemical state, on deterrence, pre-prandial and post-prandial killing of malaria vector mosquitoes (see Appendix) was parameterised, using published experimental hut data [11–14].
Also, selected parameter and parameter group combinations that seemed likely to be interdependent or to have multiplicative effects (Figure 1), were varied. For each combination, parameter values were chosen that acted on the outcome in the same direction. This was done both with those parameter values in the combination associated with lower outcomes, and with those values associated with higher outcomes. The parameters outside the combinations were kept at central values.
The relationship between the effective lifetime and pre-intervention EIR was studied in more detail by varying it over a wider range, together with LLIN attrition half-life. The effect of initial coverage was also looked at in more detail. Because the cut-off at half the impact on malaria episodes in the definition of the effective lifetime is arbitrary, sensitivity to this was studied. Finally, whether or not the protective epidemiological impact could be sustained in hypothetical situations where coverage is sustained throughout and where nets do not decay, depending on the pre-intervention EIR and coverage, was explored.
The results of the sensitivity analysis are plotted in a 'skeleton' diagram (Figure 1), which is an adaptation and expansion of a tornado diagram. A skeleton diagram not only displays how strongly the outcome varies with each parameter over its tested range, with all other parameters at their central values; it also indicates, for selected parameter combinations only, how the sensitivity to a parameter is altered if other parameters on which the parameter in question might be contingent are at their extremes.
Effective lifetime appears to be particularly sensitive to three parameters: pre-intervention EIR (annualEIR), attrition of nets (attritionOfNets half-life) and half-life of insecticide decay (insecticideDecay L).
The skeleton diagram shows that, for example, the bar lengths for coverage, deterrency, preprandialKillingEffect, and postprandrialKillingEffect are shorter for the centre-to-low (red) bars connected to the lower extreme for attritionOfNets half-life (red cross), than for those connected to the central level of attritionOfNets half-life (blue cross).
Because the attrition rate determines in part how many nets are left (together with the initial receipt proportion, defined by the parameter coverage), it strongly interacts with parameters that describe the LLIN effects on the vector population. At the lower extremes of both attritionOfNets half-life and preprandialKillingEffect, the centre-to-low (red) bar for the insecticideScalingFactor is only slightly shorter than the centre-to-low (red) bar with both attritionOfNets half-life and preprandialKillingEffect (and all other parameters) at central level. Further, one level deeper, at the lower extreme of the insecticideScalingFactor, the length of the centre-to-low (red) bar for insecticideDecay L is shorter than that with all other parameters at central level, yet still important. From Figure 1, it is not possible to establish how much the variation in each of the parameters for attritionOfNets half-life, preprandialKillingEffect, and insecticideScalingFactor contributes to the shortening of the insecticideDecay L bar. However, it is likely that attritionOfNets half-life has an important role, as insecticide decay only acts on remaining nets and is thus strongly contingent on the attrition rate. At a long attrition half-life, the sensitivity to the mean insecticide decay rate and its 'downstream' parameter (the parameter for variation in insecticide decay is contingent on the mean insecticide decay rate) will be much stronger than at a shorter attrition half-life because most nets will disappear in the latter scenario before the decay rate matters much. This is similar for the mean hole formation rate and its 'downstream' parameters (the parameters for variation in the hole formation rate and the rip factor are contingent on the mean hole formation rate). The holeRate mean is more important at the higher extreme of attritionOfNets half-life, preprandialKillingEffect, and holeScalingFactor than with all parameters at central level. Nevertheless, the outcome appears to be insensitive to parameters determining physical decay (hole formation).
To plan the timing of malaria control interventions after a mass LLIN distribution, it would be useful to understand how long the protective effect of the LLIN distribution will last. Following intervention, the impact on clinical malaria incidence reaches a maximum within a short period and is sustained for some time, after which the number of episodes rises steeply. Clearly, re-intervention (LLIN distributions or other) should take place before case incidence exceeds pre-intervention levels (due to the loss of acquired immunity). The 'effective lifetime', defined as the period from intervention until the reduction in incidence falls to a set cut-off proportion of the maximum impact, can be used to decide when to re-intervene. Because of the sharp rise in incidence it matters little whether the cut-off is set at 40, 50 or 60% of the maximum, and there is little interaction between the cut-off value and the effect of the pre-intervention EIR, indicating that the sensitivity analysis would probably give similar results if a cut-off other than 50% had been used. Even though optimal criteria for determining when to re-intervene would include economic factors, the use of the 'effective lifetime' thus provides a reasonably robust alternative.
Pre-intervention EIR level
LLINs have approximatelya the same proportionate effect on vectorial capacity at different pre-intervention EIRs. However, the apparent rate at which malaria transmission resurges from a level lower than its steady state is strongly positively related to vectorial capacity. High resurgence rates appear to be associated with shallower minima in the annual number of episodes and in the annual number of (asymptomatic) infections, but this observation results mainly from temporal smoothing.
Depending on the pre-intervention EIR, three model variants showed very different effective lifetimes from the base model R0000. These are models R0063, R0065, and R0068, in which the age-specific susceptibility is independent of exposure (contrasting with the sigmoidal relationships with exposure in the other variants), and that include extra-Poisson variation in the probability of being bitten. Model R0063 assigns most of this to inter-host variation; model R0065 is intermediate, and model R0068 assigns the variation predominantly to 'within host variation', assuming that some individuals are always more likely to get bitten than others . An explanation for why these model variants have different results is outside the scope of this paper. The models R0063 and R0065, both with inter-host variation in exposure, showed a stronger relationship between effective lifetime and pre-intervention EIR than did the base model. Heterogeneity in exposure to infectious bites between individuals is highly likely, as there are differences in mosquito access to houses  (depending on house quality and geographical situation) and in individual attractiveness to mosquitoes [17, 18]. The effect of pre-intervention EIR is thus possibly even stronger than shown in the sensitivity analysis, based on averages of 14 model variants, 12 of which ignore inter-host variation in exposure.
For the sake of simplicity, results in the sensitivity analysis were averaged over 14 model variants. This can be problematic if there is strong interaction between the effect and the model over the range studied, as was the case with the pre-intervention EIR. Sensitivity to coverage, however, was similar for all model variants. Although the 14 model variants reflect a plausible range of models, they are not necessarily evenly spread out over this range. In order to calculate a mean effect, models should ideally be weighted based on their overall fit, and on correlations between both structure and parameter values with other variants in the ensemble. Model variants with a poor fit, and/or similar to other variants included in the analysis, should receive a low weight. To perform such a weighting is a challenging task and there is a need for more methodological development in this area to address the problem appropriately.
In all the model variants, pre-intervention EIR was varied by scaling vector emergence rates. The pre-intervention EIR could possibly also be varied by changing other variables influencing vectorial capacity, such as pre-intervention survival rate, and such a pre-intervention EIR may have a somewhat different relationship with the LLIN effective lifetime.
Vector emergence was modelled as a fixed repeating seasonal pattern, independent of the adult vector population size. If the models were to include feedback from the reduced vector population due to interventions, leading to fewer emerging mosquitoes, then longer effective lifetimes would be expected, especially at locations with lower emergence rates (where local extinction of vectors might occur). Therefore, these models give a conservative estimate of effective lifetime and its dependency on emergence rates.
Effective lifetime is highly sensitive to the attrition rate. With a short attrition half-life, nets disappear before net decay can have much impact. Hole formation and insecticidal content decay are therefore of more importance with slower attrition rates. In turn, variability in hole rate and insecticide decay rate, hole size (rip factor) and the effects of decayed nets on mosquito biting and survival, are only important when hole formation and insecticidal content decay are themselves important. Thus, the shorter the attrition half-life, the more important is its accuracy relative to that of the other net decay parameters.
Insecticide decay rate is conditional on slow attrition and one of the most important factors determining the effective lifetime of an LLIN distribution. Lifetime varies strongly between the lower half-life estimate (0.5 years, characteristic for first generation LLINs) and the central value (1.5 years, characteristic for second generation LLINs). Insecticide half-life is strongly product-dependent and the values for specific products are relatively well-defined, so, despite its importance, insecticide decay is not a major contributor to uncertainty in effective lifetime, provided that mosquitoes remain sensitive to the active compound. It remains to be studied how sensitive the effective lifetime of LLINs will be to this parameter in the presence of pyrethroid resistance .
Effective lifetime of an LLIN distribution was surprisingly insensitive to parameters specifying the hole formation process in the nets. It is possible that the guessed values for the holeScalingFactor, not based on any data, were unreasonable, and that the effect of LLINs on mosquitoes wanes much faster than presumed as the hole index increases. Even though insensitive to hole parameters, it would still be useful to have evidence-based estimates as, with the emergence and spread of pyrethroid resistance [20–22], tear resistance is expected to gain in relative importance.
Coverage targets can be varied relatively easily in mass distribution campaigns, thus the effects of coverage were examined in more detail. But, effective lifetime was found to be insensitive to the initial coverage, expressed as the percentage of people that had access to an LLIN at the time of distribution. Since simulated people were not grouped into households, where they could share commonly owned nets, each covered person was simulated independently with a distinct net, used every night. This precluded explicit modelling of the distinction between household ownership of nets and personal net use, net transfer patterns within families, or local protective effects shared between net users and non-users . Thus the simulated coverage is equivalent to the percentage of people that used an LLIN. However, the insensitivity of effective lifetime to coverage implies that it would also be insensitive to measures of usage or familial correlations in ownership or usage.
The simulation experiment with sustained net coverage sheds some light onto what might happen if mass LLIN distributions are repeated at regular intervals, or supplemented by LLIN distribution through continuous delivery channels, such as the Expanded Program on Immunization (EPI). Although the assumption of no net attrition or decay is unrealistic, this experiment allows to distinguish the effects of attrition and decay from the transient dynamics induced by reducing exposure.
In settings with a medium to high pre-intervention EIR, after initial reduction the episode incidence will increase over time even in the absence of decay in number and state of LLINs, and reach a new equilibrium that is higher than its minimum. This is due to a reduction in the equilibrium level of acquired immunity caused by the reduction in exposure. The rate of decline in prevented episodes with sustained coverage is relatively slower than with a mass distribution with a four-year attrition half-life and central values for decay of both net coverage and physical and chemical states of the nets. For instance, at pre-intervention EIR of 128 ibpapa, with 70% coverage, for sustained coverage and for a mass distribution, the effective lifetimes were 15.8 and 3.0 years, respectively. This suggests that a decline in natural immunity has little impact on the effective lifetime of a single mass LLIN distribution.
Over a long period, lower acquired immunity levels will reduce the number of clinical episodes prevented by LLINs and will be accompanied by a shift in the age distribution of the clinical episodes. Some of the simulations suggest that prolonged low coverage levels in areas with a high pre-intervention EIR could lead to an increase in incidence that exceeds pre-intervention levels. However, many of the episodes in these simulations occur in older children or adults, and may well be milder than had they occurred at younger ages, resulting in overall lower case fatality. Nevertheless, these results suggest that a sustained high level of coverage of 80% could continue to suppress episodes at medium pre-intervention EIRs of 128 ibpapa and below, despite (complete) loss of acquired immunity. This supports high coverage targets for long-term sustained intervention planning. However, a separate simulation study involving multiple distribution mechanisms, such as antenatal care and EPI continuous distribution and repeated mass distributions, would be required before answering the question "what target coverage is most cost effective?".
The strong dependency of the effective lifetime of an LLIN mass distribution on pre-intervention transmission indicates 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 monitoring malaria before and during intervention programmes, particularly since there are likely to be strong variations between years and over short distances. The majority of sub-Saharan Africa's population probably falls into exposure categories where the effective lifetime is relatively long , but because exposure estimates are highly uncertain, it is necessary to consider subsequent control measures sooner than at the end of the expected effective lifetime based on an imprecise measure of transmission.
aOnly different death rates might influence LLIN ownership, and thus LLIN effect on the vectorial capacity.
Experiment parameterisation of ITN effects in OpenMalaria and experiment parameter values
In this Appendix, for selected parameters (italicised) and parameter groups (italicised) important for this study, detailed information is given on the choice of the parameter values.
The discussion of the parameters and parameter groups is organised according to the hierarchical organisation of a scenario script. The experiment's 'central' scenario specification in the machine readable language XML is given as Additional file 1, and also on the OpenMalaria site . This XML scenario is richly annotated to explain what function the parameters have. Additional information on the function of the parameters is documented in the wiki section of the project webpage 
The 'Ifakara' demography was used. The population is stationary and approximately stable: individuals move up in age group with time, and because this structure is monotonically decreasing with age, surplus individuals are out-migrated (also above maximumAgeYrs).
A population size (popSize) of 10,000 was used. This is a balance between computational effort (which increases with larger population size) and the level of stochasticity (which decreases with larger population size). At a size of 10,000, the effect of rounding integers in the population demography (which is noticeable below a size of 5,000) is minimal.
For this study, the following output variables are relevant: nets owned: the total number of nets (irrespective of physical and chemical state) present in the population; nUncomp: the number of uncomplicated malaria episodes; nSevere: the number of complicated, severe malaria episodes; nNewInfections: the number of new infections.
interventions > ITN > usage value
Usage represents the proportion of time during the night that a net is used by a simulated individual. It is not the proportion of people that use a net conditional on ownership. Because, in the current model, mosquito species bite homogenously throughout the night, the usage value can also be interpreted as the probability that host searching occurs during the time that people who own a net are using the net. In literature, this is called the 'πi value' . Govella and colleagues  define it as "the proportion of normal exposure of unprotected humans lacking nets that occurs at times and places when net users would be protected by sleeping under them". The current model version (schema 29) allows only one usage value to be set, thus it is not possible to vary the πi value for different species within the same scenario through the usage value, nor can the usage value be varied over net users.
In this experiment, all mosquitoes were assumed to display the same host searching behaviour, with a fixed probability to search in places during times when people are protected by a net (if they own one). This, as opposed to situations where the mosquito population might be, to a degree, divided into sub-populations which either always or never search only during times and in places when people are protected by a net. Such separate behaviour could be caused by genetics, and by learning; repeating the behaviour of whatever happened in the first feeding cycle.
Comparisons of indoor versus outdoor human landing catches throughout the night, combined with studies of human behaviour and the source of blood meals, give insights into the proportion of host searching mosquitoes that would encounter an ITN protected host (given ownership). However, the degree to which sub-populations exist that display different behaviours is unknown.
This parameter was varied for the sensitivity analysis. A usage value of 0.75 was used as the central value, assuming high endophagy (the propensity to bite indoors) and biting peaks after average bed-time of the population. The extreme low parameter value was taken as 0.5 and the extreme high value was taken as 1.0. These values are based on the πi values reported by Govella and colleagues  and Russell and colleagues .
interventions > ITN > holeRate mean
interventions > ITN > holeRate sigma
The value of the hole formation rate is varied among nets by multiplying with a distribution factor which is log normally distributed with mean one and the standard deviation of the log transformed variable sigma). The distribution factor is generated by taking one sample per net from a Gaussian distribution with mean zero and standard deviation one. For each parameter (holeRate, ripRate, insecticideDecay rate), the same sample is multiplied by the respective sigma and a constant (mu) added such that, once exponentiated, the mean of the variable over nets is one. For insecticideDecay rate, this constant can be chosen freely. The transformed sample is then exponentiated to obtain the respective distribution factor. This procedure implies that the distribution of holeRate, ripRate and insecticide decay rate are supposed to be covariant: nets that are heavily used decay fast both chemically and physically, whereas nets that are gently used decay slowly both chemically and physically. There is some evidence that these are indeed associated .
interventions > ITN > ripRate mean
The ripRate mean value was set equal to the value of the holeRate mean. (The ripping process was assumed to be similar to the hole formation process). This parameter was thus varied for the sensitivity analysis, but not independently.
interventions > ITN > ripRate sigma
The riprate sigma value was set equal to the value of the holeRate sigma. (The ripping process was assumed to be similar to the hole formation process). This parameter was thus varied for the sensitivity analysis, but not independently.
interventions > ITN > ripFactor value
interventions > ITN > initialInsecticide mu
The mean insecticide content of new nets (initialInsecticide mu) is set to 68.4 corresponding to the declared deltamethrin content of 68.4 mg.m-2 for long-lasting (incorporated into filaments) insecticidal nets according to WHO interim specification 333/LN/3 .
interventions > ITN > initialInsecticide sigma
The insecticide concentration of new nets is Gaussian distributed. The standard deviation (sigma) was set to 14, based on the interquartile range observed by Kilian and colleagues , for Permanet 2nd generation.
interventions > ITN > insecticideDecay L and function
interventions > ITN > insecticideDecay sigma (and mu)
interventions > ITN > attritionOfNets L and function and k
interventions > ITN > anophelesParams
After encountering a host of type i, they determine whether to attack or not (in which case they continue host tracking). These terms are not separately modelled, and included in PA i and PA, respectively. In this model, once a mosquito is determined to attack, it will either successfully feed, or die in the process. Unfed alive (UA) mosquitoes found in experimental huts are thus regarded as those that encountered a host and entered the hut but decided not to attack. This is a simplification of reality, where a mosquito may survive after unsuccessfully trying an attack. Deterrency acts on the determining phase. Deterrency is defined as one minus the relative number of affected mosquitoes (RA1 vs 2) of a host of type 1, as compared to another host of type 2. The number of affected mosquitoes is calculated as the sum of fed alive (FA), fed dead (FD) and unfed dead (UD) mosquitoes. A host type that is protected by an ITN will likely have a ratio below one, relative to a similar host type without ITN protection.
A mosquito determined to attack host i will either succeed in inserting its proboscis with probability PB i , or die during the process without inserting its proboscis, with probability PB μ i . For transmission from mosquito to human, it is important that the proboscis is inserted, and not if blood feeding was successful. However, in OpenMalaria, for simplicity, only blood-fed mosquitoes are assumed to have potentially inoculated hosts with sporozoites. ITNs will have an effect on the pre-prandial killing probability PB μ i , which can be approximated by the proportion of UD mosquitoes out of the total of determined (unfed dead, and fed dead or alive) mosquitoes.
After proboscis insertion, feeding takes place and the mosquito tries to escape from the host's vicinity, which is successful with probability PC i , or unsuccessful and the mosquito dies in the process. ITNs will have an effect on the probability of successfully escaping the host after a blood meal, called the post-prandial killing effect. This can be approximated by the proportion of FD mosquitoes out of the total of fed (dead or alive) mosquitoes.
Having escaped the host's vicinity, the mosquito will search for an appropriate resting spot. An ITN might interfere with this (through an excito-repellent effect), but this was not modelled specifically.
After a mosquito finds a resting spot it rests, and survives with probability PD i . Whereas indoor residual spraying acts during this phase, ITNs were assumed to have no influence here. Experimental hut study procedures will typically collect mosquitoes in the morning (before the resting is complete), and observe alive mosquitoes for an extended period to account for deaths during the resting phase. This is because contact with insecticide picked up during earlier phases may have a delayed effect on mortality.
Rested mosquitoes will search for an oviposition site, oviposit and start another gonotrophic cycle by host tracking.
ITNs give both personal protection by reducing the number of (infectious) bites, and reduce transmission by reducing the mosquito survival per gonotrophic cycle through increased mortality during attack and escape phases, and more time spent host tracking (with associated mortality) due to deterrency.
Host types used
Person without a net
Person with an intact treated net
Person with an intact untreated net
Person with a holed treated net
Person with a holed untreated net
Person with a treated net, holed to saturation
Person with an untreated net, holed to saturation
Person with a treated net, with undefined hole index
Person with an untreated net, with undefined hole index
interventions > ITN > anophelesParams > deterrency
Relative number of affected mosquitoes by host type
Gokool and colleagues  and Curtis and colleagues  (in their second experiment) found more fed and/or dead Anopheles gambiae s.l. in huts with untreated, holed mosquito nets than without nets. Clearly this is an artefact caused by imperfect hut traps, with more mosquitoes escaping from huts without nets, than from huts with nets, where mosquitoes got trapped under the net. It is extremely unlikely that a holed untreated net makes a person more attractive to mosquitoes than when without a net.
The geometric mean data suggest that 0.5 might be a reasonable central value for the RA of IU and IT host types versus NN, and 0.67 a good central value for an HT host type. As explained above, if the HU host type is considered not to increase the number of affected mosquitoes, a central value of 1 seems reasonable. These values were used to compute the 'medium' deterrency parameter group setting values.
The highest RA values were 0.78, 0.79 and 0.77 for IU, HT, and IT host types, respectively. These values were used to compute the 'low' deterrency parameter group setting values.
The lowest RA values were 0.41, 0.55 and 0.35 for IU, HT, and IT host types, respectively. These values were used to compute the 'high' deterrency parameter group setting values.
Although the holeScalingFactor and the insecticideScalingFactor can be specified for each ITN effect (deterrency, preprandialKillingEffect and postprandialKillingEffect) separately and for each mosquito (sub) population or species separately, the same holeScalingFactor and insecticideScalingFactor were used for ITN effects and species in these simulations.
For holes, very little information was available. Carnevale and colleagues  published data over a range of physical damage to nets (0.5, 1 and 2% of the net surface) and number of bites, but the relationships varied over the experiments. For this work, the same range of values for the holeScalingFactor as used for the insecticideScalingFactor was chosen: 0.05, 0.1 and 0.2.
interventions > ITN > anophelesParams > preprandialKillingEffect
In order to calculate a standardised effect, corrected for different mortality in the NN host type arm, probably due to varying environmental and experimental conditions, the following process was adopted:
Standardised pre-prandial killing
The mean proportion killed for NN host types was lowest (0.09), and the mean proportion for HU host types was somewhat higher (0.23). IU host types had a much larger proportion killed (0.66), followed by HT host types (0.73) and IT host types (0.84).
The fact that the average proportion killed for HU host types (0.23) was higher than that of NN host types (0.09) indicates that holed untreated nets do have a small effect on pre-prandial killing, despite the (severe) damage done to the nets. If the nets were damaged even more, the effect would presumably be smaller and eventually, at a saturation point, such a net (host type HsatU), would no longer impede flight of mosquitoes and the pre-prandial mortality would be the same as for NN. It is of interest to estimate the effect of the insecticide at this saturation point, thus for treated nets holed to saturation (HsatT). Presumably, the effect of an HsatT net would be similar to sleeping next to, not under an HT net. If the effect decays with the same shape function (linear or exponential) for both untreated nets and treated nets, this can be calculated.
Let h denote the hole index in a net in arbitrary units. The host type i can be described by the following letter combinations: NN, HU, IU, IT, HT, HsatT, HsatU, θT or θU (Table 1), depending on the net type the host uses. The letter θ indicates here an undefined hole state. For example, if a linear decay in effect is presumed, if hHsatU = 1 and hIU = 0, then PB μθU = PB μIU - (PB μIU - PB μ NN )hθU = 0.66 - (0.66 - 0.09)hθU. Thus, if PB μHsatU = 0.23, then . Then, if hHsaT = 1, hIT = 0, and hHT = hHU = 0.7544, . Thus the values are 0.66, 0.69, 0.84, and 0.09 for IU, HsatT, IT, and HsatU host types, respectively. These values were used to compute the 'medium' pre-prandial killing parameter group values.
Similarly, the highest values are 0.83, 0.89, 0.93, and 0.09 for IU, HsatT, IT, and HsatU host types, respectively. These values were used to compute the 'high' pre-prandial killing parameter group values.
The lowest values are 0.44, 0.46, 0.72, and 0.09 for IU, HsatT, IT, and HsatU host types, respectively. These values were used to compute the 'low' pre-prandial killing parameter group values.
For the medium preprandialKillingEffect level, this is thus: .
Note that for NN, PB μ NN = 0.09, and that the preprandial killing due to a net is the difference between PB μ 1 (where the subscript 1 indicates a host type with a net) and PB μ NN .
interventions > ITN > anophelesParams > postprandialKillingEffect
Standardised post-prandial killing
The mean post-prandial killing for the HU host type (0.08) was similar to that of the IU host type (0.09), Similarly, the mean post-prandial killing for the HT host type (0.64) was similar to that of the IT host type (0.65). For the mean post-prandial killing effect, the physical net state (intact or badly holed) thus has a minimal effect on post-prandial killing.
The HU and IU host type values were slightly lower than the value for NN, which is counter-intuitive, because nets are not expected to protect against mosquito mortality.
The values 0.10, 0.65, 0.0.65, and 0.10 are proposed for IU, HsatT, IT, and HsatU host types, respectively. These values were used to compute the 'medium' postprandialKillingEffect parameter group setting values.
Similarly, the values 0.36, 0.70, 0.70, and 0.10 for IU, HsatT, IT, and HsatU host types, respectively, were used to compute the 'high' postprandialKillingEffect parameter group setting values.
The values 0.10, 0.50, 0.50, and 0.10 for IU, HsatT, IT, and HsatU host types, respectively, were used to compute the 'low' postprandialKillingEffect parameter group setting values.
with PC μ i the probability of dying (for a mosquito) after biting as a result of the biting process on host type i. For example for the medium postprandialKillingEffect level this is:
interventions > ITN > timed coverage
interventions > importedInfections
From time step zero onwards, 10 infections per 1,000 population per year were imported by stochastically infecting individuals in the population. This was done to ensure that malaria would not be eliminated from the simulated population, which might overestimate the protective effect of an intervention.
Thus, even if an intervention provides full protection to the entire population, 1% of the population will be infected once per year. These infections do not necessarily develop into disease episodes. These could be seen as infections obtained while travelling to a malarious area.
The "Tanzania ACT" health system was used, described elsewhere .
entomology > annualEIR
In the OpenMalaria schema version used (Schema 29), during the warm-up phase of the simulation run, the mosquito emergence rate is scaled such that the average annual malaria transmission to an adult (expressed as the number infectious bites per adult per annum) is approximately equal to this parameter value. Post warm-up, in the absence of interventions, this annual EIR is approximately constant. However, interventions such as ITN distributions will affect the transmission. Hence, the annualEIR is referred to is the pre-intervention EIR. For the sensitivity analysis, the annualEIR was varied between 8, 16, and 32 ibpapa. For the in-depth analysis, also values of 2, 4, 64, 128 and 256 ibpapa were used, covering a low to high range of stable malaria transmission.
entomology > mode and name
The "Namawala" seasonality and relative species abundance was used , with a dynamic transmission mode. Briefly, there are three species (Anopheles gambiae s.s., Anopheles arabiensis and Anopheles funestus. An. gambiae s.s. and An. arabiensis had the same seasonality. The seasonality of the natural logarithm of the observed density was fitted with a third order Fourier series, thus with a mean term, two terms for the annual frequency and two terms for the semi-annual frequency. The other parameter values in this section are described elsewhere .
We acknowledge our colleagues in the Swiss TPH Biostatistics & Computational Sciences Group, particularly Michael Tarantino, for help with the models and simulations. We also thank Christian Lengeler, Don de Savigny, Jo Lines, Sarah Hoibak and Albert Kilian for helpful discussions, and Amena Briët for language editing. Simulations were run with the help of volunteers providing computational power. This work was supported through the Malaria Modelling Project #OPP1032350 funded by the Bill and Melinda Gates Foundation (BMGF) and through the NetWorks project funded by USAID in partnership with John Hopkins Bloomberg School of Public Health.
- WHO Global Malaria Programme: World Malaria Report 2010. 2010, Geneva: World Health OrganizationGoogle Scholar
- Control of Neglected Tropical Diseases WHO Pesticide Evaluation Scheme, Global Malaria Programme Vector Control Unit: Guidelines for monitoring the durability of long-lasting insecticidal mosquito nets under operational conditions. 2011, Geneva: World Health Organization, WHO/HTM/NTD/WHOPES/2011.5Google Scholar
- Kilian A: How long does a long-lasting insecticidal net last in the field?. Public Health Journal. 2010, 21: 43-47.Google Scholar
- Hawley WA, Phillips-Howard P, ter Kuile F, Terlouw DJ, Vulule JM, Ombok M, Nahlen B, Gimnig JE, Kariuki SK, Kolczak MS, Hightower AW: Community-wide effects of permethrin-treated bed nets on child mortality and malaria morbidity in western Kenya. Am J Trop Med Hyg. 2003, 68 (Suppl. 4): 121-127.PubMedGoogle Scholar
- Killeen GF, Smith TA: Exploring the contributions of bed nets, cattle, insecticides and excitorepellency to malaria control: a deterministic model of mosquito host-seeking behaviour and mortality. Trans R Soc Trop Med Hyg. 2007, 101: 867-880. 10.1016/j.trstmh.2007.04.022.PubMed CentralView ArticlePubMedGoogle Scholar
- OpenMalaria: a simulator of malaria epidemiology and control. [http://code.google.com/p/openmalaria/]
- 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 (2 Suppl): 1-10.PubMedGoogle Scholar
- Chitnis N, Smith T, Steketee R: A mathematical model for the dynamics of malaria in mosquitoes feeding on a heterogeneous host population. J Biol Dyn. 2008, 2: 259-285. 10.1080/17513750701769857.View ArticlePubMedGoogle Scholar
- Chitnis N, Hardy D, Smith T: A periodically-forced mathematical model for the seasonal dynamics of malaria in mosquitoes. Bull Math Biol. 2012, [http://www.ncbi.nlm.nih.gov/pubmed/22218880]Google Scholar
- Smith T, Ross A, Maire N, Chitnis N, Studer A, Hardy D, Brooks A, Penny M, Tanner M: Ensemble modeling of the likely public health impact of a pre-erythrocytic malaria vaccine. PLoS Med. 2012, 9: e1001157-10.1371/journal.pmed.1001157.PubMed CentralView ArticlePubMedGoogle Scholar
- Curtis CF, Myamba J, Wilkes TJ: Various pyrethroids on bednets and curtains. Mem Inst Oswaldo Cruz. 1992, 87 (Suppl 3): 363-370.View ArticlePubMedGoogle Scholar
- Lines JD, Curtis CF, Myamba J, Njau R: Tests of repellent or insecticide impregnated curtains, bednets and anklets against malaria vectors in Tanzania. 1985, World Health Organization, WHO/VBC/85.920Google Scholar
- Lines JD, Myamba J, Curtis CF: Experimental hut trials of permethrin-impregnated mosquito nets and eave curtains against malaria vectors in Tanzania. Med Vet Entomol. 1987, 1: 37-51. 10.1111/j.1365-2915.1987.tb00321.x.View ArticlePubMedGoogle Scholar
- Gokool S, Smith DF, Curtis CF: The use of PCR to help quantify the protection provided by impregnated bednets. Parasitol Today. 1992, 8: 347-350. 10.1016/0169-4758(92)90072-A.View ArticlePubMedGoogle Scholar
- Central scenario. [http://openmalaria.googlecode.com/svn/application_deployment/examples/centre_scenario.xml]
- Smith T, Charlwood JD, Takken W, Tanner M, Spiegelhalter DJ: Mapping the densities of malaria vectors within a single village. Acta Trop. 1995, 59: 1-18. 10.1016/0001-706X(94)00082-C.View ArticlePubMedGoogle Scholar
- Qiu YT, Smallegange RC, Van Loon JJ, Ter Braak CJ, Takken W: Interindividual variation in the attractiveness of human odours to the malaria mosquito Anopheles gambiae s.s. Med Vet Entomol. 2006, 20: 280-287. 10.1111/j.1365-2915.2006.00627.x.View ArticlePubMedGoogle Scholar
- Lindsay SW, Adiamah JH, Miller JE, Pleass RJ, Armstrong JR: Variation in attractiveness of human subjects to malaria mosquitoes (Diptera: Culicidae) in The Gambia. J Med Entomol. 1993, 30: 368-373.View ArticlePubMedGoogle Scholar
- Briët O, N'Guessan R, Hardy D, Rowland M: Epidemiological effectiveness of insecticide treated nets in the presence of pyrethroid resistance [abstract]. Trop Med Int Health. 2011, 16 (Suppl. 1): 140-Google Scholar
- Kawada H, Dida GO, Ohashi K, Komagata O, Kasai S, Tomita T, Sonye G, Maekawa Y, Mwatele C, Njenga SM, Mwandawiro C, Minakawa N, Takagi M: Multimodal pyrethroid resistance in malaria vectors, Anopheles gambiae s.s., Anopheles arabiensis, and Anopheles funestus s.s. in western Kenya. PLoS ONE. 2011, 6: e22574-10.1371/journal.pone.0022574.PubMed CentralView ArticlePubMedGoogle Scholar
- Chandre F, Darrier F, Manga L, Akogbeto M, Faye O, Mouchet J, Guillet P: Status of pyrethroid resistance in Anopheles gambiae sensu lato. Bull World Health Organ. 1999, 77: 230-234.PubMed CentralPubMedGoogle Scholar
- Doannio JM, sou-Yovo J, Diarrassouba S, Chauvancy G, Darriet F, Chandre F, Henry MC, Nzeyimana I, Guillet P, Carnevale P: Efficacité des moustiquaires pré-imprégnées de perméthrine Olyset Net® en zone de résistance des vecteurs aux pyréthrinoïdes. I--Evaluation entomologique. Med Trop (Mars). 1999, 59: 349-354.Google Scholar
- Lines JD, Myamba J, Curtis CF: Experimental hut trials of permethrin-impregnated mosquito nets and eave curtains against malaria vectors in Tanzania. Med Vet Entomol. 1987, 1: 37-51. 10.1111/j.1365-2915.1987.tb00321.x.View ArticlePubMedGoogle Scholar
- Maire N, Shillcutt S, Walker DG, Tediosi F, Smith T: Cost effectiveness of the introduction of a pre-erythrocytic malaria vaccine into the Expanded Program on Immunization in sub-Saharan Africa: analysis of uncertainties using a stochastic individual-based simulation model of Plasmodium falciparum malaria. Value Health. 2012, 14 (8): 1028-1038.View ArticleGoogle Scholar
- Maire N, Smith T, Ross A, Owusu-Agyei S, Dietz K, Molineaux L: A model for natural immunity to asexual blood stages of Plasmodium falciparum malaria in endemic areas. Am J Trop Med Hyg. 2006, 75 (2 Suppl): 19-31.PubMedGoogle Scholar
- Killeen GF, Kihonda J, Lyimo E, Oketch FR, Kotas ME, Mathenge E, Schellenberg JA, Lengeler C, Smith TA, Drakeley CJ: Quantifying behavioural interactions between humans and mosquitoes: evaluating the protective efficacy of insecticidal nets against malaria transmission in rural Tanzania. BMC Infect Dis. 2006, 6: 161-10.1186/1471-2334-6-161.PubMed CentralView ArticlePubMedGoogle Scholar
- Govella NJ, Okumu FO, Killeen GF: Insecticide-treated nets can reduce malaria transmission by mosquitoes which feed outdoors. Am J Trop Med Hyg. 2010, 82: 415-419. 10.4269/ajtmh.2010.09-0579.PubMed CentralView ArticlePubMedGoogle Scholar
- Russell TL, Govella NJ, Azizi S, Drakeley CJ, Kachur SP, Killeen GF: Increased proportions of outdoor feeding among residual malaria vector populations following increased use of insecticide-treated nets in rural Tanzania. Malar J. 2011, 10: 80-10.1186/1475-2875-10-80.PubMed CentralView ArticlePubMedGoogle Scholar
- Tami A, Mubyazi G, Talbert A, Mshinda H, Duchon S, Lengeler C: Evaluation of Olyset insecticide-treated nets distributed seven years previously in Tanzania. Malar J. 2004, 3: 19-10.1186/1475-2875-3-19.PubMed CentralView ArticlePubMedGoogle Scholar
- Kilian A, Byamukama W, Pigeon O, Gimnig J, Atieli F, Koekemoer LL, Protopopoff N: Evidence for a useful life of more than three years for a polyester-based long-lasting insecticidal mosquito net in Western Uganda. Malar J. 2011, 10: 299-10.1186/1475-2875-10-299.PubMed CentralView ArticlePubMedGoogle Scholar
- Kilian A, Byamukama W, Pigeon O, Atieli F, Duchon S, Phan C: Long-term field performance of a polyester-based long-lasting insecticidal mosquito net in rural Uganda. Malar J. 2008, 7: 49-10.1186/1475-2875-7-49.PubMed CentralView ArticlePubMedGoogle Scholar
- Deltamethrin long-lasting (incorporated into filaments) insecticidal net. WHO specifications and evaluations for public health pesticides. 2010, World Health OrganizationGoogle Scholar
- Smith SC, Joshi UB, Grabowsky M, Selanikio J, Nobiya T, Aapore T: Evaluation of bednets after 38 months of household use in northwest Ghana. Am J Trop Med Hyg. 2007, 77 (6 Suppl): 243-248.PubMedGoogle Scholar
- Report of the 12th WHOPES working group meeting 8-11 December 2008. 2009, World Health Organization, WHO/HTM/NTD/WHOPES/2009.1Google Scholar
- Lindsay SW, Adiamah JH, Miller JE, Armstrong-Schellenberg JR: Pyrethroid-treated bednet effects on mosquitoes of the Anopheles gambiae complex in The Gambia. Med Vet Entomol. 1991, 5: 477-483. 10.1111/j.1365-2915.1991.tb00576.x.View ArticlePubMedGoogle Scholar
- Carnevale P, Bitsindou P, Diomande L, Robert V: Insecticide impregnation can restore the efficiency of torn bed nets and reduce man-vector contact in malaria endemic areas. Trans R Soc Trop Med Hyg. 1992, 86: 362-364. 10.1016/0035-9203(92)90219-3.View ArticlePubMedGoogle Scholar
- Otten M, Lines JD: Where did the LLINs go?. An analysis of data from 7 countries with the most recent surveys (2008-2009). 2011, Slide 22.Google Scholar
- Tediosi F, Maire N, Smith T, Hutton G, Utzinger J, Ross A, Tanner M: An approach to model the costs and effects of case management of Plasmodium falciparum malaria in sub-Saharan Africa. Am J Trop Med Hyg. 2006, 75 (2 Suppl): 90-103.PubMedGoogle Scholar
- Smith T, Charlwood JD, Kihonda J, Mwankusye S, Billingsley P, Meuwissen J, Lyimo E, Takken W, Teuscher T, Tanner M: Absence of seasonal variation in malaria parasitaemia in an area of intense seasonal transmission. Acta Trop. 1993, 54: 55-72. 10.1016/0001-706X(93)90068-M.View ArticlePubMedGoogle Scholar
- Ross A, Maire N, Molineaux L, Smith T: An epidemiologic model of severe morbidity and mortality caused by Plasmodium falciparum. Am J Trop Med Hyg. 2006, 75 (2 Suppl): 63-73.PubMedGoogle Scholar
- Maire N, Tediosi F, Ross A, Smith T: Predictions of the epidemiologic impact of introducing a pre-erythrocytic vaccine into the expanded program on immunization in sub-Saharan Africa. Am J Trop Med Hyg. 2006, 75 (2 Suppl): 111-118.PubMedGoogle Scholar
- Penny MA, Maire N, Studer A, Schapira A, Smith TA: What should vaccine developers ask? Simulation of the effectiveness of malaria vaccines. PLoS ONE. 2008, 3: e3193-10.1371/journal.pone.0003193.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/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.