- Open Access
Target product profile choices for intra-domiciliary malaria vector control pesticide products: repel or kill?
Malaria Journalvolume 10, Article number: 207 (2011)
The most common pesticide products for controlling malaria-transmitting mosquitoes combine two distinct modes of action: 1) conventional insecticidal activity which kills mosquitoes exposed to the pesticide and 2) deterrence of mosquitoes away from protected humans. While deterrence enhances personal or household protection of long-lasting insecticidal nets and indoor residual sprays, it may also attenuate or even reverse communal protection if it diverts mosquitoes to non-users rather than killing them outright.
A process-explicit model of malaria transmission is described which captures the sequential interaction between deterrent and toxic actions of vector control pesticides and accounts for the distinctive impacts of toxic activities which kill mosquitoes before or after they have fed upon the occupant of a covered house or sleeping space.
Increasing deterrency increases personal protection but consistently reduces communal protection because deterrent sub-lethal exposure inevitably reduces the proportion subsequently exposed to higher lethal doses. If the high coverage targets of the World Health Organization are achieved, purely toxic products with no deterrence are predicted to generally provide superior protection to non-users and even users, especially where vectors feed exclusively on humans and a substantial amount of transmission occurs outdoors. Remarkably, this is even the case if that product confers no personal protection and only kills mosquitoes after they have fed.
Products with purely mosquito-toxic profiles may, therefore, be preferable for programmes with universal coverage targets, rather than those with equivalent toxicity but which also have higher deterrence. However, if purely mosquito-toxic products confer little personal protection because they do not deter mosquitoes and only kill them after they have fed, then they will require aggressive "catch up" campaigns, with behaviour change communication strategies that emphasize the communal nature of protection, to achieve high coverage rapidly.
The most important front line vector control strategies for malaria prevention rely on killing mosquitoes that enter human houses by delivering insecticidal products to these domestic targets in the form of indoor residual spray (IRS) or long-lasting insecticidal nets (LLINs) [1, 2]. The common rationale underpinning these strategies is based on two well-established biological phenomena: 1) that the most important malaria vectors prefer to feed on humans and rest inside houses and 2) that a mosquito must feed several times on humans and, therefore repeatedly risk exposure to such insecticidal measures, before it is old enough to acquire, incubate and then transmit malaria to any human [3, 4]. As the most common and important species of Plasmodium that cause human malaria infections are strict anthroponoses, malaria vectorial capacity of a given mosquito species is directly and closely related to its human-feeding propensity so these two phenomena obviously co-occur in the most important vector populations . This is particularly true in sub-Saharan Africa where, with some interesting exceptions, the bulk of human exposure to Anopheles gambiae and Anopheles funestus has occurred inside houses and these species feed almost exlusively upon humans [6–8]. As a result, even coverage of only half of the human population with LLINs or IRS can deliver huge reductions of transmission and substantial alleviation of malaria burden in settings where the challenge of eliminating malaria is greatest [9, 10]. Few public health interventions achieve such massive positive externality in the form of protecting those not directly covered [9–11] and the elegant way in which these measures exploit the biology of both the parasite and the vector is both intuitive and appealing [3, 4, 12]. The potential for community-level impact that is far greater than what can be achieved with personal protection alone is obviously hugely attractive [2, 11, 12], but this simple rationale and impressive recent progress with implementation masks a complex set of important product profile choices, which have thus far been made in the absence of decisive evidence or clear evaluation criteria.
However, the two most commonly used pesticides for controlling adult malaria vector mosquitoes, namely the synthetic pyrethroids and dichlorodiphenyltrichloroethane (DDT), combine two very distinct modes of action: 1) conventional toxicity which kills mosquitoes exposed to the pesticide while feeding or attempting to feed upon covered humans, 2) deterrence of mosquitoes away from those humans resulting from either irritation upon direct contact with the treated surface or even through spatial repellence from a distance of several meters [13–15]. Pyrethroids exhibit a strong combination of both contact irritant and spatial repellent properties, so that IRS and LLIN using these compounds often deter as many mosquitoes as they kill [16–20]. DDT is the only commonly used alternative to the pyrethroids for IRS and clearly has strong spatial repellency, as well as strong insecticidal effects upon mosquitoes that are not deterred and actually make contact [13, 14].
While high levels of deterrence enhance the personal protection afforded by a pesticide product and, therefore, uptake by the public, it may also attenuate or even reverse communal protection  because it diverts mosquitoes to non-users  rather than killing them outright. Theoretical analysis suggests that where vectors have a strict preference for human hosts, or their preferred alternative hosts are absent, such deterrent properties may be counterproductive or even dangerous . In principle, diversion of mosquitoes away from protected individuals might cancel out the community-level benefits to non-users arising from decreased mosquito survival and infection rates and could even result in increased exposure because bites are increasingly focused on the unprotected portion of the population .
Numerous large scale field trials of insecticidal nets or IRS have produced overwhelming encouraging results [9, 10] but it is critical to note that these impacts result from products with a combination of deterrent and insecticidal properties. Even larger studies will be required to conclusively distinguish the community-level impacts of alternative profiles with deliberately and distinctly formulated toxic versus deterrent product profiles. It is therefore perhaps unsurprising that no such field trial has been conducted. While current guidelines for evaluating LLIN and IRS products in experimental huts  provide clear instructions on how to quantify personal protection and overall mortality rates of mosquitoes, it is not explicitly required to distinguish between toxic effects that kill mosquitoes before or after they feed and, with one exception , trials following these guidelines report only combined total mortality rates. Furthermore, consensus has yet to be attained regarding which of these evaluation criteria should be considered as primary and secondary or how the relative merits of these properties should be compared when evaluating existing products or designing new ones.
A process-based mathematical model of malaria transmission is outlined here, which captures the sequential interaction between deterrent and toxic actions of vector control pesticides and which accounts for the distinctive impacts of slow and fast-acting toxicity upon mosquitoes (Figure 1). This model is applied to explore how the interaction of deterrent and toxic actions affects both overall transmission intensity and its distribution across user and non-user groups in malarious communities. Furthermore, the consequent influence of alternative and hybrid product profiles upon the choice of optimal delivery system strategy is outlined and further potential applications for this model are discussed.
Initially, a recently published deterministic model  was applied to elucidate how interactions between deterrent and insecticidal properties of hypothetical LLIN or IRS products might affect their impact upon malaria control when applied at high coverage across large populations. This exercise revealed that neither this formulation nor any of its predecessors [12, 15, 25] produced plausible, internally consistent outcomes for the probabilities of a mosquito attacking an encountered LLIN user and of successfully obtaining a blood meal when the proportion of human exposure that occurs at times when LLINs are used (π i ) was set to values less than 1. The uncoupling of the impacts of π i upon repellence and insecticidal activity became particularly obvious when the hypothetical LLIN was defined as being 100% repellent (θ Δ = 1) and 100% insecticidal (θ μ = 1): such simulations indicated that mosquitoes were directly killed by these nets, despite the expectation that coupled and complete repellency should prevent any such fatal contact. Furthermore, this implausible exposure of mosquitoes to direct mortality risk despite complete diversion away from such hazard increased as the proportion of exposure the LLIN can potentially prevent (π i ) decreased. Close examination of equations 6 and 7 of the original formulation  reveals how the previous approach caused the uncoupling of this conditionality to produce increasingly unrealistic outcomes as the fraction of exposure of indoor interventions for which the repellency does not apply (1- π i ) increases, namely increasing estimated exposure of mosquitoes to the insecticidal activity and consequently nonsensically increasing insecticide-related mortality.
These flaws arise from inconsistent definition of protection, which was sometimes, but not always, considered to be synonymous with simply using a net. In simple terms, using a net is something that covered individuals only do for approximately one third of a typical day so protection must be assumed to be partial, even for the most nocturnal, indoor-biting vectors, regardless of net efficacy . Such interactions between mosquito and human behaviours are best summarized for indoor interventions such as LLINs or IRS in terms of the proportion of human exposure that would otherwise occur indoors (π i ) . Published field estimates of this parameter for African malaria vector populations indicate that this proportion may fall far short of its optimal maximum value of 1 and may well be dropping in response to increasing selection pressure as ITN coverage increases [25–27]. Here these components of previously published formulations [12, 15, 24, 25] are harmonized so that this increasingly important de facto gap in coverage is treated with far greater clarity and internal consistency (See Table 1 for parameter definitions). In the interests of brevity and simplicity of language, the model description below refers consistently to an LLIN product but relates equally to an IRS product. Here, the essential changes to the existing model are described in detail and a brief but comprehensive description of the overall model is provided.
Coverage, protection and host availability to mosquitoes
Protection is defined as being conditional upon both using a net and, more specifically, using a net at times when transmission occurs . The de facto protective coverage of humans (C h,p ) is therefore defined as being the product of crude coverage (C h ) and the proportion of human exposure that occurs indoors while asleep at times when LLINs are used (π i ) .
The total availability for attack by mosquitoes  of protected (A h,p ) and unprotected humans (A h,u ) in the community is redefined so that individual users of nets exposed at times when they do not use them are considered to be unprotected. Thus, the effect of π i upon host availability is applied as a conditional probability that affects population-level parameters in a coupled manner, rather than a probability which is independently applied to each of distinct individual-scale processes it influences in an uncoupled manner. The total availability of hosts protected against attack by using a net is therefore adjusted for this fraction of exposure which is directly preventable (π i ): The availability for attack of net users at times when those nets are used and therefore protect them is calculated as follows:
Where ah,p is the availability for attack of an individual protected human, N h is the number of humans and Ch is the crude coverage, estimated as the reported nightly usage rate.
The availability of the remaining fraction of humans which are unprotected (A h,u ) because either they do not use a net (Ah,0,u) or because they are exposed during times when the net is not used (A h,net,u ) can be calculated as follows where a h,u is the attack availability of an unprotected individual.
Which can also be expressed simply as follows in manner consistent with equation 2:
Similarly, to estimate the total availability of blood (Z) from these same categories of human hosts, equivalent formulae based on the availability of blood from individual protected (z h,p ) and unprotected (z h,u ) human hosts are applied:
By redefining protection and thus allowing for attenuated reductions of impact of insecticidal protection by human behaviours  at this population level the consistency and simplicity of parameters describing individual-level processes is improved. Individual mean (a h,p and z h,p ) and population total availability parameters (A h,p and Z h,p ) of the model are specified and calculated separately for protect and unprotected users and derived directly from the simpler respective un-weighted terms γ h,p and ϕ h,p , respectively. For diversion, this is achieved directly, similar to some previous formulations :
Where Δ h,p is the probability that a mosquito will divert away from an encountered, protected human host. However, the probability of feeding is expressed more explicitly than before, to consider only mortality which occurs before the mosquito feeds (μ h,p,pre ) rather than total mortality (μ h,p ) including those which feed but die soon afterwards:
Where μ h,p,pre is the probability that a mosquito will die before feeding if it attacks a protected host.
These terms are calculated as follows based on the probabilities of diversion (Δ h,u ) and death (μ h,u ) for unprotected humans, combined with the additional probability of diversion (θΔ) and death before feeding (θ μ,pre ) caused by the deterrent and insecticidal properties of the net:
This distinction, between toxic activities that act fast enough to prevent blood feeding and those that do not, necessitates that the total excess attack-related mosquito mortality resulting from using an LLIN (θ μ ) is specified as the sum of the excess mortality which occurs before (θ μ,pre ) or after (θ μ,post ) obtaining a blood meal:
While insecticide-related mosquito mortality occurring after the mosquito has fed on the protected host does not contribute to personal protection, it does contribute to community-level suppression of malaria transmission by reducing population mean mosquito survival. The term μ h,p is therefore calculated separately as follows:
This distinction between killing mosquitoes before or after feeding on the protected host allows the proportion of blood meals derived from humans (Qh) to be calculated as previously described  based on this revised feeding probability term. Note, however, that this parameter therefore includes fatal blood meals obtained from insecticide-protected humans which mosquitoes never live long enough to digest. The meaning of parameters depending on the availabilities various categories of attackable hosts (A), rather than blood sources per se (Z) described above, such as the duration of the host-seeking interval (ηov) and the probability of surviving host attack per feeding cycle (Pγ)  are unaffected. Note also that, as described below in equation 14, the latter logically remains based on μh,p rather than the new μh,p,pre term.
Implications of redefining coverage, protection and host availability for mosquito population parameters
Previous versions of this model incorporated the lack of an effect of an LLIN on outdoor malaria transmission π i by either treating it as a weighting term for calculating population mean values for feeding probability and attack-related mortality [12, 25] or by applying directly to the individual level diversion and mortality processes [15, 24]. The changed manner in which protection, coverage and availability are conceptually distributed (equations 1 to 7 and associated text), means that population-level parameters such as the proportion of blood meals obtained from humans (Q h ) and mean host-seeking interval (η v ), can all be simply calculated in terms of total host attack (A) and blood (Z) availability parameters exactly as previously described [15, 24]. Note, however, that this means that the published breakdowns of these population-level parameters into functions of the products of numbers of hosts (N) and mean individual availabilities (a and z, respectively) [15, 24] are no longer valid.
For other population-level parameters, simpler, more direct and intuitively satisfying derivations are implied. For example, this approach allows ready estimation of the probability of surviving host attack per feeding cycle (P γ ) based on the mosquito mortality rates (μ) and corresponding community-wide total attack availabilities (A) of protected humans (h,p), unprotected humans (h,u) and cattle (c).
Otherwise, all the mosquito population parameters are calculated exactly as previously described, and outlined as follows.
The mean seeking interval for vertebrate hosts (η v ) can be calculated as the reciprocal of total host availability (A), using estimates of these feeding probabilities and their corresponding encounter rates [24, 28]:
The feeding cycle length (g) is calculated as the sum of the durations of the gestation period (g), the oviposition site-seeking interval (η v ) and the vertebrate host-seeking interval (η v ):
Survival across all phases of the gonotrophic cycle is calculated as the distinct daily survival probability during each phase to the power of the respective time intervals, namely the host-seeking interval (η v ), gestation period (g) and oviposition site-seeking interval (η o ). The daily survival probability of a resting mosquito is defined as P and the survival probabilities during host-seeking and oviposition site-seeking are assumed to be equal and are both defined using the term P ov . The survival rate per feeding cycle (P f ) was estimated as the combined probability that a vector survives gestation (Pg ), oviposition site-seeking , vertebrate host-seeking and the eventual attack of a host :
Similarly, the human blood index is calculated as the proportion of total blood availability accounted for by humans :
The biodemography component of the model is adapted to a daily cycle and cumulative survival up to each age (x) is estimated as follows :
Similarly, the sporozoite infection prevalence of mosquitoes at each age is considered in days, accounting for superinfection:
where κ denotes the mean infectiousness of the human population to vector mosquitoes  and n is the duration of the sporogonic development period of the parasite from ingestion to infective sporozoite stages . Survival and infectveness probabilities are calculated up to 40 days, after which the contributions of mosquitoes in these age classes to transmission become negligible. Note that P x is multiplied by S x to obtain the corresponding probability of being both alive and infective (I x ) on each day
The following mosquito lifetime biodemographic parameters are calculated by summing these three age-specific outcomes as previously described [15, 30]. The number of human bites the average mosquito takes in a lifetime (b h ) is defined as the sum of the probabilities of surviving and feeding on a human at each age (x):
Note that to enable incorporation of survival-dependent emergence rates, the number of human bites on all hosts, rather than just humans, per mosquito lifetime (b) is similarly calculated:
Accounting for superinfection, the number of infectious bites on humans per mosquito lifetime (β h ) is calculated as the product of the human blood index and sum of the products of the probabilities of biting and being infectious at each age [15, 30]:
Again, the number of sporozoite-infected bites on all hosts per mosquito lifetime (β), regardless of whether that host is susceptible to infection or not, is calculated similarly but ignoring the human blood index term:
The overall sporozoite prevalence in the vector population (S) can then be calculated as β h divided by b h :
Epidemiological outcomes: dealing with partially covered, partially protected humans
Also, the entomologic inoculation rate (EIR) for non-users (EIRh,0) can be directly estimated based on the share of all available blood sources which a single non-user represents (z h,u /Z) multiplied by the total number of infectious bites on all hosts (β ; equation 24) by all emerging mosquitoes (E):
Alternatively, this parameter may be estimated by considering only infectious bites on human hosts (β h ; equation 23) and therefore considering only the share of available human blood which such an individual represents:
Nevertheless, it is essential to retain the protection-weighted mean terms for parameters which reflect the properties of individual net users who are only covered with the protective LLIN for proportion of their normal exposure (π i ) and uncovered and unprotected for the remained (1 - π i ). These terms are therefore retained but annotated more distinctly than previously  so that the attack probability (γ h,net rather than ) and feeding probability (ϕ h,net rather than ) reflect the mean of protected and unprotected periods for net users, but cannot be confused with the corresponding probabilities for net users during the specific periods when they are protected (γ h,p and ϕ h,p , respectively).
Consequently, derived terms such as attack availability (a net rather than ) and blood availability (z h,net rather than ), as well as corresponding terms for relative attack availability (λ h,net rather than λ h,p ) and exposure to bites (ψ h,net rather than ψ h,p ) compared with non-users, can be calculated as previously described.
Consequently, the EIR experienced by net users can be calculated by five different but consistent means:
Additionally, the mean EIR experienced in scenario Ω by the mixture of net users and non-users which comprise the community (ψh,Ω) can be independently calculated in three distinct ways which yield consistent results. Consistent with equation 22 of Okumu et al., this parameter can be estimated by simply weighting the EIR parameters for net users and non-users according to crude coverage and the gap in coverage, respectively:
However, it is also possible to calculate exactly the same values with a simpler formula derived from first principles, based on the assumptions of the very first of this family of models :
Reassuringly, identical values can also be calculated as described above by weighting the availability of blood from protected and unprotected individuals according to de facto protective coverage (C h,p ) rather than crude coverage (C h ).
Similarly, the relative exposure of non-users and users of nets (ψh,0,Ωand ψh,net,Ωrather than ψh,p,Ωand ψh,p,Ω, respectively) and community-wide mean relative exposure (ψh,Ω) in a given intervention scenario (Ω) is calculated exactly as previously described except that the terms EIRh,0,0, EIRh,0,Ωand EIRh, net,Ωreplace EIRh, u,0EIRh, u,Ωand EIRh, p,Ωto denote the EIR experienced by non-users in a scenario with no intervention and that of non-users and users under intervention scenario Ω, respectively:
Survival-dependent mosquito proliferation
Previous formulations of this model have assumed that larval habitats are always at their carrying capacity so the annual emergence rate of mosquitoes (E) is fixed, regardless of vector survival rates. In reality, vector populations experience dramatic seasonal fluctuations in larval habitat availability so while this assumption is largely true during drier times of the year when the quantity of habitat is static or contracting, it is rarely limiting during the onset or peak of the rains when vector populations can grow at their maximum reproduction rate [31, 32]. Furthermore, observations of the differential impact of insecticide-treated nets upon sibling species composition within the An. gambiae complex [33, 34] and impact of indoor-residual spraying upon inter-species competition within the An. funestus group [35, 36], both confirm that oviposition input into larval habitats does limit vector population sizes. These simulations were, therefore, executed both with and without allowing for adult survival-dependent emergence rates which were calculated as follows.
Emergence rate was assumed to vary simply and linearly with mean number of successfully-completed feeding cycles by adult mosquitoes (b; Equation 22). Emergence rate in a given vector control scenario (EΩ) was therefore calculated as the product of the maximum emergence rate expected in the absence of any adult mosquito control (E0) and the relative value of the mean number of feeding cycles per mosquito lifetime in that scenario (bΩ), compared with such baseline conditions (b0):
The calculations for the feeding cycle duration itself (f) as the sum of the gestation (g), oviposition site-seeking (η o ) and vertebrate blood-seeking (η v ) intervals are exactly as previously described :
Consistent with the previously published definition of host availability , it is assumed that protecting hosts does not alter their location, or the rate at which they are encountered by kinesis, but rather extends the spatial distribution of locations to which mosquitoes must disperse to in order to obtain blood. As hosts are increasingly protected, a greater mean number of hosts must be encountered before a blood meal can be successfully obtained. Longer host-seeking intervals, that include a greater mean number of unsuccessful host encounters, will inevitably result in a mean increase in the distance and duration of subsequent return journeys to oviposition sites. Calculation of the oviposition site-seeking interval (η o ) is adapted to account for the expectation that mosquitoes forced to fly further and longer in search of blood will also have to fly proportionally further and longer in search of oviposition sites once the blood meal has been digested and eggs are matured. This term is calculated as the reciprocal of aquatic habitat availability, termed Z a rather than A a , as previously described , to maintain consistency with the separate definitions of rates of initiation and completion of resource utilization processes here and elsewhere :
However, here this term (Z a ) is assumed to vary proportionally with vertebrate blood availability (Z) as it changes from baseline (0) to intervention (Ω) scenarios, reflecting the intrinsically endogenous relationship between host and aquatic habitat availability:
Parameterization of the model
The term π i is set at a values of 0.90, consistent with published reports from areas with high coverage of untreated nets [25, 37] and historical field observations for African vector populations from across Africa (Huho et al., Unpublished) or at 0.50, reflecting more recent observations from vector populations exposed to high coverage of LLINs, IRS or house screening [25–27, 38].
Previous modelling investigations [15, 39] have illustrated that the eventual impact of deterrent pesticide products upon malaria transmission exposure for non-users is very sensitive to the assumed value for mosquito survival while foraging for vertebrate blood or oviposition site resources (P ov ), parameter for which no field estimates exist to the authors knowledge. A range of values of were examined in the absence of any intervention measure (C h = 0) to determine an approximate value that is most compatible with the observed biodemographic profiles of real populations of vectors and sporogonic parasites in the field. Implausibly low values for the proportion of mosquitoes surviving each feeding cycle (P f ) except at high assumed values of P ov , approaching the likely upper limit of 0.90 defined by the estimated survival rate of resting mosquitoes (Figure 2).
Furthermore, surprizingly low sporozoite prevalence (S) rates were predicted for both species, especially at the lower end of the range of assumed P ov , suggesting that values of the latter are high in nature. However, actual field estimates for survival per feeding cycle (P f = 0.62) and sporozoite prevalence (S = 0.016) for the village of Namwawala in the 1990s, where the crucial human population size (N h ) and availability parameters (a h ) were obtained from, were quite low by the standards African vector populations in the absence of LLINs or IRS and compare reasonably well with the An. gambiae scenario simulated here where P ov > 0.85. Note that although transmission in this village was dominated by An. arabiensis at this time, no significant cattle population existed so the An. gambiae scenario assuming no alternative hosts is most representative of this setting. While daily survival rates for actively foraging mosquitoes (P ov ) must be somewhat lower than for resting mosquitoes, normal parity and sporozoite rates of African vector populations can only be plausibly explained if this difference is small, so P ov was set at 0.85 for all subsequent simulations.
All other parameter settings for the two vector population scenarios (An. arabiensis representing a mosquito that can exploit non-human hosts compared with An. gambiae which is almost exclusively dependent on humans for blood) are as previously described for a village with 1,000 people and an equal number of cattle .
Specifically, the mean individual attack availability of unprotected humans (a h,u ) to An. arabiensis in this particular Tanzanian village in the 1990s was calculated as the reciprocal of the estimate of the mean vertebrate host-seeking interval (η v ), based on the distribution of ovariolar stalks dilation status among host-seeking specimens , divided by the number of humans present at the time [24, 28]. This approach to estimating this parameter was first described  before clear distinction between the availability of individual hosts for attack (a h,u ) and the availability of individual blood sources per se (z h,u ) had been explicitly outlined  but is even more appropriate when the former is specified. The same a h,u value of 1.2 × 10-3 attacks per host per night per host-seeking mosquito was assumed for An. gambiae. The mean individual attack availability of cattle (a c ) for each species was calculated by multiplying the equivalent parameter for humans (a h,u ) by field estimates  of the relative availability of cattle blood, compared to that of humans (ψ c ), for both vector species , yielding estimates of 1.9 × 10-3 and 2.5 × 10-5 attacks per host per night per host-seeking mosquito, for An. arabiensis and An. gambiae, respectively. Note that this calculation assumes that for unprotected hosts, the probability of successfully feeding upon an attacked host is equivalent for the two host types (μ h,u = μ c ) so that the relative availability of cattle for attack is equivalent to the relative availability of cattle blood (λ c = ψ c ).
Consistent with previous simulations, the maximum emergence rate of mosquitoes in the absence of adult mosquito control measures (E0) was set at 2 × 107 adult mosquitoes per year. Except where stated otherwise, crude coverage of humans was set at 80% (C h = 0.8) in line with the Roll Back Malaria targets for coverage of all age groups with LLINs which represents an ambitious but realistically achievable target for most malaria afflicted developing nations.
The fundamental trade-off between toxic and deterrent actions (Figure 1) is clearly illustrated by the simulation results presented in Figure 3, all of which are based on the assumption that 80% of humans use LLINs (C h = 0.8). Predictions for toxic-deterrent hybrid product profiles (θ μ,pre = 0.5, θ μ,post = 0, θ Δ > 0) converge with those for purely deterrent product profiles (θ μ,pre = θ μ,post = 0, θ Δ > 0) once deterrence reaches 100% efficacy and prevents any fatal contact with the active ingredient (θ Δ = 1 so that γ h,p = 0). This is to say that given maximum diversion, the probability that a mosquito would attack a protected host becomes zero. A number of further observations suggest this trade-off should be carefully considered when defining a target product profile for developing or selecting a malaria vector control pesticide formulation.
A partially efficacious but purely fast-acting toxic product ((θ μ,pre = 0.5, θ μ,post = 0, θ Δ > 0) consistently delivers better protection of non-users than a completely efficacious but purely deterrent ((θ μ,pre = 0, θ μ,post = 0, θ Δ = 1.0) product (Figure 3). A reasonable degree of community-level protection for non-users is accrued where attractive, non-human hosts exist for diverted mosquitoes to feed upon. However, in the absence of such alternative blood sources, the unprotected minority of the human population could experience greater exposure and this negative externality increases with increasing deterrence (Figure 3b). Furthermore, the consistently strong community-level benefits obtained by non-users when their neighbours use pesticide products with purely toxic activity profiles are undermined in all scenarios by supplementing these lethal effects with increasing levels of deterrence (Figure 3).
Where alternative animal hosts exist, benefits for users of toxic nets are modestly enhanced as high levels of personal protection provided by strong deterrent properties (θ Δ > 0.5) are realized (Figure 3). However, this results in an approximate break-even scenario, in terms of mean relative exposure across the entire community because increased benefit for users is offset by reduced benefit for non-users (Figure 3). Where alternative sources of blood are absent, increasing deterrence actually progressively undermines protection of users because the increased personal protection conferred is more than counterbalanced by dramatically attenuated community-level impact (Figure 3).
Note that for all of these conclusions, the model which includes survival-dependent emergence rates (Figure 3c and 3d versus 3a and 3b) improves the predicted outcomes for purely deterrent products and toxic-deterrent hybrids but in no case does so dramatically enough to alter the overall trend or conclusions reached (Figure 3). These simulations suggest that purely toxic products are preferable to purely deterrent ones and that enhancing the personal protection afforded by a toxic product by increasing its repellent or irritant properties will consistently undermine or even reverse communal protection of non-users. In fact, where vectors lack alternative non-human hosts, increasing deterrence may even undermine benefits for users because the degree of community-level protection obtained with purely toxic products is far greater than personal protection at the high levels of coverage now considered as healthy targets for any malaria control programme [1, 2].
Figure 4 illustrates how such counterintuitive predictions may be rationalized by examining the underlying biodemographic parameters describing the vector and sporogonic-stage parasite populations, which ultimately determine impact on malaria transmission. Vector survival per feeding cycle (P f ) is the most important single determinant of malaria transmission intensity other than temperature and is substantially reduced by toxic, deterrent and toxic-deterrent hybrid products where no alternative blood hosts exist (Figure 4a). Where alternative hosts occur, only toxic products with little or no deterrence are predicted to usefully reduce vector survival (P f ). Regardless of whether alternative hosts are present, increasing deterrence of toxic products consistently weakens impact upon this most important target for adult malaria vector control, modest reductions of which result in quasi-exponential suppression of transmission [4, 15, 39]. As the impact upon vector reproduction (E) has been modelled as a linear function of the number of bites taken per lifetime (b h ), itself a simple function of survival (P f ) , it is unsurprising that the impact of these various product profiles mirrors that upon survival (Figure 4b). Being a squared term in all malaria transmission models [4, 30, 39], the proportion of blood meals that the vector population obtains from humans is the next most important determinant of malaria transmission intensity at global  and local level [15, 39, 42]. Where alternative sources of blood are available, deterrence can dramatically reduce this target parameter in its own right and also enhances the impact of toxic products when added as a supplementary activity (Figure 4c). In the absence of alternative hosts, no toxic, deterrent or hybrid product has any meaningful impact on this target parameter. Consistent with outputs of previous formulations , increasing deterrence can greatly extend the feeding cycle length (f) of the vector where no alternative non-human hosts exits but has a very modest effect where they are present (Figure 4d). Consistent with the recently revised, distinct definitions of host and blood availabilities , toxicity has no influence on this determinant of mosquito survival (P f ), feeding frequency (1/ f ), reproduction (E) and transmission potential (b h , S).
In summary, toxic products consistently reduce vector survival and reproduction rates, especially in the absence of alternative blood sources. In places where such non-human preferred hosts exist, toxic products only reduce the proportion of blood meals that are human but have no impact on vector feeding cycle length. In contrast, purely deterrent products only have useful impacts upon the proportion of blood meals obtained from humans where alternative hosts exist and upon feeding cycles length where they are absent. Deterrent products, therefore, impact one of these two target parameters or the other and it is notable that neither has as strong an influence upon transmission as survival, particularly when further impact upon mosquito reproduction rates is considered.
By definition (Figure 1), increasing deterrence of a product inevitably increases the proportion of available blood that non-users constitute at any given coverage level (Figure 5) and therefore the share of mosquito bites they experience, regardless of whether that product is toxic or not. When the predicted extent of this inequitable redistribution of biting mosquitoes (Figure 5) is combined with the predicted impacts upon the biodemographic properties of the vector population (Figure 4a to 4d), the overall impact is to increase biting rates for non-users (Figure 4e) even where alternative blood sources are absent so vector survival (Figure 4a) and reproduction (Figure 4b) are reduced because the availability of blood becomes limiting. This effect is so dramatic that, even for toxic products, increased exposure of non-users to bites can occur at high levels of deterrence (θ Δ > 0.8). While such negative externality in the form of diverting biting mosquitoes to unprotected non-users has been envisaged and discussed previously, the simulated impact of increasing deterrence of toxic products upon the sporozoite infection prevalence are perhaps more interesting. Consistent with previous simulations , purely deterrent products consistently reduce sporozoite prevalence (Figure 4f) by either lowering human blood indices where alternative hosts are available (Figure 4c) or by reducing survival (Figure 4a) and extending feeding cycle length (Figure 4d) where they are not. More surprising is the prediction that increasing the deterrence of a toxic product can attenuate impact upon sporozoite prevalence. In the case of vector populations lacking an alternative non-human host, this rebound of sporozoite infection prevalence arising from enhancing the personal protection provided by the product, by increasing irritant of repellent properties, is quite substantial. In fact this weakening of impact upon sporozoite prevalence may be as important a contributor to the dramatic attenuation of overall impact upon transmission (Figure 3b and 3d) as redistribution of bites to unprotected non-users (Figure 5).
Figure 6 illustrates just how much more efficacious a purely toxic product can be. In both vector-host scenarios, toxic (Figure 6c and 6d) or toxic-deterrent hybrids (Figure 6e and 6f) are clearly superior to non-toxic deterrent products (Figure 6a and 6b). Obviously, the toxic but not deterrent product confers less personal protection than the toxic-deterrent hybrid but correspondingly provides the best communal protection for non-users as coverage increases. Even in the Anopheles arabiensis scenario where alternative hosts are available, the benefit to users of a purely toxic product arising from combined personal and community-level protection exceeds that of a toxic-deterrent hybrid at 57% coverage where baseline transmission primarily occurs indoors (π i = 0.9) and only 27% coverage where an equal amount of baseline transmission occurs outdoors (π i = 0.5). For An. gambiae-dominated transmission systems without alternative blood hosts, the advantage of purely toxic products conferring less protection than those supplemented with deterrence is even more dramatic and obvious, with almost three orders of reduction of transmission possible within feasible coverage targets and the purely toxic product providing greater protection than the hybrid at 22 and 12% coverage, respectively, where most (π i = 0.9) and half (π i = 0.5) of baseline transmission occurs indoors. Not only do purely toxic products have greater efficacy at reasonable coverage levels, they are also more robust to attenuation by outdoor-feeding behaviours in the target vector population (π i = 0.5) because, under such conditions, deterrent products simply divert mosquitoes to feeding on users at times when they are unprotected, especially when no alternative non-human hosts are available.
With the exception of the two bottom panels of Figure 6, all toxic actions simulated thus far are assumed to kill mosquitoes before they can bite the occupant of the house or net. This kind of scenario is best reflected in reality by LLINs with which the pyrethroid insecticide activity is specifically applied to a physical barrier between the attacking mosquito and the protected host so that most dead mosquitoes collected in experimental hut trials are unfed. However, in the case of IRS with non-deterrent insecticides, such as entomopathogenic fungi , bendiocarb , chlorpyrifos methyl , and even pyrethroid-based LLINs that have been depleted of insecticide after several years of use , most mosquitoes killed succeed in feeding before dying so little, if any, personal protection is conferred. Figure 6g and 6h represent such a scenario and this is reflected in the fact that the predicted degree of protection of users and non-users is identical because this is exclusively mediated by community-level suppression of transmission. Obviously, a purely insecticidal product which kills mosquitoes fast enough to prevent blood feeding and therefore also confers personal protection (Figure 6e and 6f) is preferable to one that kills them afterwards and does not (Figure 6g and 6h). Nevertheless, even a purely toxic product, which confers no personal protection because it only kills mosquitoes after they have fed (Figure 6g and 6h), is a consistently better option in terms of protection of non-users than products with deterrent properties, regardless of whether (Figure 6c and 6d) or not (Figure 6a and 6b) that product also has insecticidal activity that kills mosquitoes before feeding. Comparing the residual transmission levels achieved with products that confer only community-level protection through purely post-feeding toxicity with that attained by more conventional products with purely deterrent or deterrent plus pre-feeding insecticidal activities (Figure 7), shows that the non-user is always better off with the former.
For zoophagic vectors with alternative hosts available that predominantly feed indoors (π i = 0.9), deterrent plus pre-feeding insecticidal activity attains lower residual transmission for users than purely post-feeding insecticidal activity. However, such a scenario with most of the alternative hosts and mosquito feeding activity occurring indoors is probably unusual and occurs in a limited number of settings across the tropics. In all other scenarios, especially where half of transmission occurs outdoors (π i = 0.5), the purely post-feeding insecticide confers superior overall protection to users, despite complete lack of personal protection, once a minimum coverage threshold is surpassed. Compared with pure deterrents, overall protection of users becomes greater for the purely post-feeding insecticidal product at quite modest crude coverage levels (49 and 20% for An. arabiensis with alternative hosts and An. gambiae without them, respectively) where most transmission occurs indoors (π i = 0.9) and even lower thresholds (35 and 14%, respectively) where outdoor feeding and/or resting is common (π i = 0.5). Compared with products combining deterrent with pre-feeding insecticidal activity analogous to LLINs, similar patterns were observed, with the consistent disadvantage of purely post-feeding toxicity where alternative hosts exist and most transmission occurs indoors being reversed when outdoor transmission becomes important and crude coverage exceeds 65%, while it becomes consistently advantageous for vectors lacking alternative non-human hosts at remarkably low coverage thresholds of 39% for predominantly indoor transmission and only 22% where half of transmission occurs outdoors.
The idea that deterrency reduces the impact of toxic activities of pesticides upon mosquito survival is long-established  and was discussed extensively during the previous global campaign to eradicate malaria [47–49] as well as the beginning of the more recent drive to promote scale up of LLINs and IRS for control purposes . Deliberate design of pesticide-based vector control products to match ideal target product profiles has recently been reprioritized as an important issue  now that more ambitious programmes to control, eliminate or even eradicate of malaria are back on the global agenda [51, 52]. The process-explicit model of malaria transmission described here captures the sequential interaction between deterrent and toxic actions of vector control pesticides. In simple terms, it is not realistic to expect that one can discourage mosquitoes from making contact with an active ingredient without compromising the ability of that pesticide to kill them (Figure 1). Sub-lethal exposure that deters mosquitoes inevitably reduces the proportion which is subsequently exposed to higher, lethal doses. In fact, the extreme example outlined on the right hand side of all the panels in Figures 3 and 4, wherein the predicted impacts of products with and without toxic activities converge once 100% deterrency is achieved, clearly demonstrates that this is a choice which must be made: increasing deterrency and personal protection must always be traded off against reduced toxicity-mediated mosquito mortality and potent communal level protection where high coverage is achieved.
The assumptions and definitions of this model (Figure 1 and Methods) are also fully compatible with recent recommendations that toxic activities and both forms of deterrence, namely contact irritance and spatial repellence, are distinct and that each pesticide-affected mosquito collected in an experimental hut trial should be classified as having either responded in a manner characteristic of only one of these possible outcomes . While parameter estimates from published studies have been deliberately avoided to minimize any appearance of recommending for or against specific commercial product choices, this model can be readily and directly parameterized from existing, standardized experimental hut evaluations. The diversion term θ Δ is estimated directly as the proportional change in the number of mosquitoes which either do not enter the hut (deterrence) or which leave unfed (excito-repellency) but do not subsequently die. The mortality terms θ μ,pre and θ μ,post are estimated as the increased proportion of all mosquitoes caught in a hut with a given LLIN or IRS product which were found dead or that subsequently died which were either unfed or fed respectively. However, to enable the application of this model to such experimental hut study outcomes, published summaries will need to explicitly distinguish between pre- and post-feeding mortality  and will ideally include the raw data as supplementary online material. The model described also accounts for the distinctive impacts of toxic activities, which kill mosquitoes before or after they have fed upon the occupant of a covered house or sleeping space. A variety of well-established domestic vector control products and emerging new technologies only kill mosquitoes after they have fed because they are applied as IRS formulations or because they are slow acting. Such alternatives to DDT or pyrethroids include entomopathogenic fungi , bendiocarb , chlorpyrifos methyl , and even pyrethroid-based LLINs that have been depleted of insecticide after several years of use , can take hours or days to kill most of the exposed mosquitoes but clearly can deliver massive levels of malaria control if sufficient coverage can be achieved. Only two previous models distinguish between the effects of pesticidal products that kill mosquitoes before and after they feed upon humans [14, 53]. While one only considers processes that occur in houses and does not capture the community-level effects of different product profiles upon transmission , the other does not account for outdoor biting and, like previous versions of this model , inaccurately treats diversion and mortality as independent, rather than sequentially coupled, events .
The specific results presented suggest that if high coverage levels can be achieved that are consistent with current World Health Organization targets [1, 2], purely toxic products with no deterrence are predicted to generally provide superior protection to non-users and even users, especially where vectors feed exclusively on humans and a substantial amount of transmission occurs outdoors. Remarkably, this is even the case if that product confers no personal protection and only kills mosquitoes after they have fed. Products with purely mosquito-toxic profiles may be preferable to those with equivalent toxicity but which confer superior personal protection because of higher deterrence for programmes with universal coverage targets. Purely mosquito-toxic products which confer modest personal protection because they lack deterrence, or which confer none because they only kill mosquitoes after they have fed, will therefore require aggressive "catch up" campaigns to achieve high coverage rapidly and behaviour change communication strategies that emphasize the communal nature of protection.
As with all mathematical predictions, these predictions should only be considered as evidence of plausible hypotheses based on simplifying assumptions and imprecise parameterization. Lessons from learned from historic mistakes, specifically setting malaria prevention policy based on overconfident interpretation of malaria transmission models [3, 4], are as relevant today as they ever were .
For example, one notable simplification to keep in mind is that complete gonotrophic concordance, meaning that each egg batch requires one and only one blood meal, has been assumed. In reality, the first blood meal typically requires at least one additional pre-gravid blood meal to achieve mature phase II development of the ovaries [55–57] and additional blood meals may even be taken during subsequent gonotrophic cycles . While such increased feeding frequency would undoubtedly increase malaria transmission intensity in the absence of interventions such as LLINs or IRS, it would also be expected to increase the frequency of contact with such measures that mosquitoes would be exposed to early in their lives. Correspondingly, incorporating these subtle aspects of mosquito behaviour would most probably enhance the predicted impact of these measures upon transmission and therefore strengthen, rather than weaken, the contrasts between alternative target product profiles suggested here.
The only potentially major inaccuracy that seems obvious from the outputs of this model lies in the prediction that purely deterrent products will provide weak communal protection for non-users and may even increase their exposure. While this phenomenon appears plausible in theory and has been documented by field trials of some topical repellents , the experimental design of that study define situations in which only single users were protected, equivalent to negligible community level coverage (C h ≈ 0) so community-level effects were neither realized nor evaluated. Furthermore, these predictions seem slightly at odds with observations from field trials of community-wide use of essentially untreated mosquito nets in both Tanzania  and Papua New Guinea . In both cases, high coverage of nets lacking meaningful pesticidal properties but deterring mosquitoes through simple physical barrier effects successfully reduced malaria transmission. Combined with the anecdotal but reasonable attribution of reduced malaria transmission in many settings to housing improvements conferring similarly direct protection through similar physical barriers , these net trials suggest that the disappointing predictions for purely deterrent products described here should be interpreted with a degree of caution. The most obvious possible explanation for such possible discrepancies is the uncertainty associated with survival rate of mosquitoes foraging for blood or aquatic habitat and the extreme sensitivity of predictions to this parameter value and to baseline total availabilities of these resources [15, 28, 39, 62–64]. To go beyond speculation based on sensitivity analysis of these critical but, as yet, unmeasured parameters, will clearly require the development of robust field methods, notably trapping of gravid Anopheles seeking oviposition sites .
With some notable exceptions, these simulations compare well with recent, less generalized, modelling analyses which examine choices between specific product types and combinations thereof [65, 66]. Deliberately, no specific product has been named, nor has any measured parameter value for any specific product been set in any of these simulations. Instead, the product parameters have been tuned them across the full range of possible values so that ideal target product profiles can be objectively outlined for manufacturers and their clients to aim for prospectively rather than restrict discussion to the relative merits of currently available products and product combinations. Nevertheless, the parameter space explored here encompasses all the specific examples of product types evaluated in recent modelling analyses [65, 66], resulting in predictions which are readily comparable in broad terms (Figures 3 and 6). Both these complementary recent studies [65, 66] also conclude that IRS with a highly deterrent product such as DDT will have less impact than a predominantly insecticidal product such as IRS with bendiocarb or pyrethroid-based LLINs. However, their conclusions regarding combining such product types differ somewhat and the existing evidence base is insufficient to inform which of these three formulations appears most accurate. Chitnis et al. predict that supplementing a predominantly insecticidal LLIN products  with a highly deterrent one such as IRS with DDT will have a larger impact upon transmission than either one as a stand-alone measure. In contrast, the simulations of Yakob et al., suggests the opposite: that placing a deterrent product in the same house as a predominantly insecticidal one will undermine the superior impact of the latter for exactly the reasons outlined here and captured in the convergence of outcomes with toxic and non-toxic products in Figure 3.
Perhaps the most important observation about the lack of consensus between these three model formulations is that sufficient field data do not exist to reliably compare them in terms of their predictive value. Recent reviews of the impact of IRS , and specifically IRS combined with insecticidal nets [10, 67], both conclude that rigorous, large-scale, randomized controlled trials are conspicuous by their absence. An abundance of descriptive studies unambiguously demonstrate that IRS has massive overall impact and that combining with ITNs gives generally improved personal protection [10, 49]. To the knowledge of the authors, however, no study yet exists in which the exclusively communal protection afforded to residents of unsprayed houses in IRS programmes has been measured as rigorously as it has for non-users of insecticidal nets in communities with high coverage levels [11, 60, 68–72]. Given these limitations in the evidence base for IRS as a stand-alone prevention strategy, it is perhaps unsurprising that the evidence base to support decisions about combining LLINs and IRS is insufficient and has become a common point of discussion for both theoreticians and practitioners [10, 65–67].
Despite these limitations in both the consistency of outputs from alternative existing models and the empirical evidence base from the field, important lessons can be learned from these simulations which are intuitive and for which no caveats seem obvious. Although no evidence, based on rigorously randomized trials, for the probability of the deterrence-related attenuation of insecticidal impact have been reported, the existing descriptive evidence base presents a strong case for the plausibility of this phenomenon.
The effect of insecticidal attack was enhanced by the use of non-irritant insecticides
In fact, the ideal target product profile outlined here was already suggested during the previous malaria eradication era, when the impact of DDT which has a mixed deterrent-plus-toxic profile, was contrasted with that of Dieldrin which acts by contact toxicity only :
In many instances, Dieldrin proved to be more effective than DDT, but its higher cost, its toxicity to mammalians, and the fast-spreading resistance of A[nopheles] gambiae to this insecticide limited its further use in Africa
This model presented herein simply strengthens, explains and generalizes the plausibility of this argument, highlighting the lack of affordable, safe alternatives to Dieldrin with similarly non-deterrent properties. Three decades later, with insecticide resistance on the rise  and increasing levels of exophagy being reported for residual vector populations in Africa  and Asia , it is likely that several such active ingredients with distinct, complementary mechanisms will be required to prevent and manage insecticide resistance in the long term. These simulation results, therefore, serve as a timely reminder of the need for increased investment in development and evaluation of insecticidal products with purely toxic modes of action to achieve improved and sustained malaria vector control.
Even if the worst-case scenarios predicted here are confirmed through large-scale trials, it is important to remember that this analysis is restricted to typical LLIN or IRS products that are used indoors. One of the most interesting phenomena that this model captures, which is increasingly relevant as the importance of outdoor-biting vectors is recognized [26, 27, 54, 75], is that the advantage of purely toxic products becomes greater where vector mosquitoes tend to feed outdoors (Figures 6 and 7). This suggests that deterrent activities can not only divert mosquitoes to animals or to humans lacking such products but also to the users themselves at times of the day when they are outside of the house and unprotected. This new insight arises directly and intuitively from the reformulation of how coverage and protection have been conceptualized and expressed mathematically. Further extensions of this approach may be useful for examining a wider diversity of possible pesticidal vector control products that target mosquitoes outside of houses [76–78] and even away from humans [24, 79–81]. This conceptual and mathematical formulation represents a useful new tool for rational design of malaria vector control products. Furthermore, the way in which coverage and protection are conceptualized in this formulation represents a substantive change in thinking that may also enable more lucid re-examination of what these terms really mean in practice [82, 83].
RBM: Roll Back Malaria Global Strategic Plan 2005-2015. 2005, Geneva: World Health Organization
WHO: Insecticide treated mosquito nets: A position statement. 2007, Geneva: Global Malaria Programme; World Health Organization, 10-
Garrett-Jones C: Prognosis for interruption of malaria transmission through assessment of the mosquito's vectorial capacity. Nature. 1964, 204: 1173-1175. 10.1038/2041173a0.
MacDonald G: The epidemiology and control of malaria. 1957, London: Oxford University Press
Kiswewski AE, Mellinger A, Spielman A, Malaney P, Sachs SE, Sachs J: A global index representing the stability of malaria transmission. Am J Trop Med Hyg. 2004, 70: 486-498.
Gillies MT, Coetzee M: A supplement to the Anophelinae of Africa South of the Sahara (Afrotropical region). 1987, Johannesburg: South African Medical Research Institute
Gillies MT, DeMeillon B: The Anophelinae of Africa South of the Sahara (Ethiopian zoogeographical region). 1968, Johannesburg: South African Institute for Medical Research
White GB: Anopheles gambiae complex and disease transmission in Africa. Trans R Soc Trop Med Hyg. 1974, 68: 279-301.
Lengeler C: Insecticide-treated bed nets and curtains for preventing malaria. Cochrane Database Syst Rev. 2004, CD000363-
Pluess B, Tanser FC, Lengeler C, Sharp BL: Indoor residual spraying for preventing malaria. Cochrane Database of Systematic Reviews. 2010, 4: CD006657-
Hawley WA, Phillips-Howard PA, ter Kuile FO, Terlouw DJ, Vulule JM, Ombok M, Nahlen BL, Gimnig JE, Kariuki SK, Kolczak MS, Hightower AW: Community-wide effects of permethrin-treated bednets on child mortality and malaria morbidity in western Kenya. Am J Trop Med Hyg. 2003, 68 (Supplement 4): 121-127.
Killeen GF, Smith TA, Ferguson HM, Abdulla S, Mshinda H, Lengeler C, Kachur SP: Preventing childhood malaria in Africa by protecting adults from mosquitoes with insecticide-treated nets. PLoS Med. 2007, 4: e229-10.1371/journal.pmed.0040229.
Grieco JP, Achee NL, Chareonviriyaphap T, Suwonkerd W, Chauhan K, Sardelis MR, Roberts DR: A new classification system for the actions of IRS chemicals traditionally used for malaria control. PLoS One. 2007, 2: e716-10.1371/journal.pone.0000716.
Roberts DR, Alecrim WD, Hshieh P, Grieco JP, Bangs M, Andre RG, Chareonviriphap T: A probability model of vector behavior: effects of DDT repellency, irritancy, and toxicity in malaria control. J Vector Ecol. 2000, 25: 48-61.
Killeen GF, Smith TA: Exploring the contributions of bednets, cattle, insecticides and excito-repellency 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.
Malima RC, Magesa SM, Tungu PK, Mwingira V, Magogo FS, Sudi W, Mosha FW, Curtis CF, Maxwell C, Rowland M: An experimental hut evaluation of Olyset nets against anopheline mosquitoes after seven years use in Tanzanian villages. Malar J. 2008, 7: 38-10.1186/1475-2875-7-38.
N'Guessan R, Asidi A, Boko P, Odjo A, Akogbeto M, Pigeon O, Rowland M: An experimental hut evaluation of PermaNet® 3.0, a deltamethrin-piperonyl butoxide combination net, against pyrethroid-resistant Anopheles gambiae and Culex quinquefasciatus mosquitoes in southern Benin. Trans R Soc Trop Med Hyg. 2010, 104: 758-765. 10.1016/j.trstmh.2010.08.008.
Corbel V, Chabi J, Dabire RD, Etang J, Mwane P, Pigeon O, Akogbeto M, Hougard JM: Field efficacy of a new mosaic long-lasting mosquito net (PermaNet 3.0) against pyrethroid-resistant malaria vectors: a multi centre study in Western and Central Africa. Malar J. 2010, 9: 113-
Akogbeto M, Padonou GG, Gbenou G, Irish S, Yadouleton A: Bendiocarb, a potential alternative against pyrethroid resistant Anopheles gambiae in Benin, West Africa. Malar J. 2010, 9: 204-10.1186/1475-2875-9-204.
Tungu P, Magesa S, Maxwell C, Malima R, Massue D, Sudi W, Myamba J, Pigeon O, Rowland M: Evaluation of PermaNet 3.0 a deltamethrin-PBO combination net against Anopheles gambiae and pyrethroid resistant Culex quinquefasciatus mosquitoes: an experimental hut trial in Tanzania. Malar J. 2010, 9: 21-10.1186/1475-2875-9-21.
Moore SJ, Davies C, Cameron MM: Are mosquitoes diverted from repellent-using individuals to non-users? Results of a field study in Bolivia. Trop Med Int Health. 2007, 12: 532-539. 10.1111/j.1365-3156.2006.01811.x.
WHO: Guidelines for testing mosquito adulticides for indoor residual spraying and treatment of mosquito nets, vol. WHO/CDS/NTD/WHOPES/GCDPP/2006.3. 2006, Geneva, Switzerland: World Health Organization
Van Bortel W, Chinh VD, Berkvens D, Speybroeck N, Trung HD, Coosemans M: Impact of insecticide-treated nets on wild pyrethroid resistant Anopheles epiroticus population from southern Vietnam tested in experimental huts. Malar J. 2009, 8: 248-10.1186/1475-2875-8-248.
Okumu FO, Moore SJ, Govella NJ, Chitnis N, Killeen GF: Potential benefits, limitations and target product-profiles of odor-baited mosquito traps as a means of malaria control. PLoS One. 2010, 5: e11573-10.1371/journal.pone.0011573.
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.
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.
Bugoro H, Cooper RD, Butafa C, Iro'ofa C, Mackenzie DO, Chen C-C, Russell TL: Bionomics of the malaria vector Anopheles farauti in Temotu Province, Solomon Islands: issues for malaria elimination. Malar J. 2011, 10: 133-10.1186/1475-2875-10-133.
Killeen GF, Seyoum A, Knols BGJ: Rationalizing historical successes of malaria control in Africa in terms of mosquito resource availability management. Am J Trop Med Hyg. 2004, 71 (Supplement 2): :87-93.
Killeen GF, Ross A, Smith TA: Infectiousness of malaria-endemic human populations to vector mosquitoes. Am J Trop Med Hyg. 2006, 76 (Suppl 2): 38-45.
Killeen GF, McKenzie FE, Foy BD, Schieffelin C, Billingsley PF, Beier JC: A simplified model for predicting malaria entomologic inoculation rates based on entomologic and parasitologic parameters relevant to control. Am J Trop Med Hyg. 2000, 62: 535-544.
Briet OJ: A simple method for calculating mosquito mortality rates, correcting for seasonal variations in recruitment. Med Vet Entomol. 2002, 16: 22-27. 10.1046/j.0269-283x.2002.00335.x.
Russell TL, Lwetiojera D, Knols BGJ, Takken W, Killeen GF, Ferguson HM: Linking individual phenotype to density-dependent population growth: the influence of body size on the population dynamics of malaria vectors. Proc R Soc B. 2011, E-pub ahead of print
Bayoh MN, Mathias DK, Odiere MR, Mutuku FM, Kamau L, Gimnig JE, Vulule JM, Hawley WA, Hamel MJ, Walker ED: Anopheles gambiae: historical population decline associated with regional distribution of insecticide-treated bed nets in western Nyanza Province, Kenya. Malar J. 2010, 9: 62-10.1186/1475-2875-9-62.
Russell TL, Lwetoijera DW, Maliti D, Chipwaza B, Kihonda J, Charlwood JD, Smith TA, Lengeler C, Mwanyangala MA, Nathan R, Knols BGJ, Takken W, Killeen GF: Impact of promoting longer-lasting insecticide treatment of bednets upon malaria transmission in a rural Tanzanian setting with pre-existing high coverage of untreated nets. Malar J. 2010, 9: 187-10.1186/1475-2875-9-187.
Gillies MT, Smith A: Effect of a residual house-spraying campagn on species balance in the Anopheles funestus group: The replacement of Anopheles gambiae Giles with Anopheles rivulorum Leeson. Bull Entomol Res. 1960, 51: 248-252.
Gillies MT, Furlong M: An investigation into the behaviour of Anopheles parensis Gillies at Malindi on the coast of Kenya. Bull Entomol Res. 1964, 55: 1-16. 10.1017/S0007485300049221.
Killeen GF, Kihonda J, Lyimo E, Okech FR, Kotas ME, Mathenge E, Schellenberg J, Lengeler C, Smith TA, Drakeley C: 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.
Ogoma SB, Kannady K, Sikulu M, Chaki PP, Govella NJ, Mukabana WR, Killeen GF: Window screening, ceilings and closed eaves as sustainable ways to control malaria in Dar es Salaam, Tanzania. Malar J. 2009, 8: 221-10.1186/1475-2875-8-221.
Saul A: Zooprophylaxis or zoopotentiation: the outcome of introducing animals on vector transmission is highly dependent on the mosquito mortality while searching. Malar J. 2003, 2: 32-10.1186/1475-2875-2-32.
Charlwood JD, Smith T, Billingsley PF, Takken W, Lyimo EOL, Meuwissen JHET: Survival and infection probabilities of anthropophagic anophelines from an area of high prevalence of Plasmodium falciparum in humans. Bull Entomol Res. 1997, 87: 445-453. 10.1017/S0007485300041304.
White GB, Magayuka SA, Boreham PFL: Comparative studies on sibling species of the Anopheles gambiae Giles complex (Dipt., Culicidae): bionomics and vectorial activity of species A and species B at Segera, Tanzania. Bull Entomol Res. 1972, 62: 295-317. 10.1017/S0007485300047738.
Killeen GF, McKenzie FE, Foy BD, Bogh C, Beier JC: The availability of potential hosts as a determinant of feeding behaviours and malaria transmission by mosquito populations. Trans R Soc Trop Med Hyg. 2001, 95: 469-476. 10.1016/S0035-9203(01)90005-7.
Kiszewski A, Mellinger A, Spielman A, Malaney P, Sachs SE, Sachs J: A global index representing the stabililty of malaria transmission. Am J Trop Med Hyg. 2004, 70: 486-498.
Scholte EJ, Ng'habi K, Kihonda J, Takken W, Paaijmans K, Abdulla S, Killeen GF, Knols BG: An entomopathogenic fungus for control of adult African malaria mosquitoes. Science. 2005, 308: 1641-1642. 10.1126/science.1108639.
N'Guessan R, Boko P, Odjo A, Chabi J, Akogbeto M, Rowland M: Control of pyrethroid and DDT- resistant Anopheles gambiae by application of indoor residual spraying or mosquito nets treated with a long-lasting organophosphate insecticide, chlorpyrifos-methyl. Malar J. 2010, 9: 44-10.1186/1475-2875-9-44.
Muirhead-Thomson RC: Mosquito behaviour in relation to malaria transmission and control in the tropics. 1951, London: Edward Arnold & Co
Muirhead-Thomson RC: The significance of irritability, behaviouristic avoidance and allied phenomena in malaria eradication. Bull World Health Organ. 1960, 22: 721-734.
Carnevale P, Mouchet J: La lutte antivectorielle au Cameroun. Passé-présent-avenir. Réflexions. Bull Soc Pathol Exot. 2001, 94: 202-209.
Kouznetsov RL: Malaria control by application of indoor spraying of residual insecticides in tropical Africa and its impact on community health. Tropical Doctor. 1977, 7: 81-93.
Pates H, Curtis C: Mosquito behavior and vector control. Annu Rev Entomol. 2005, 50: 53-70. 10.1146/annurev.ento.50.071803.130439.
Feachem R, Sabot O: A new global malaria eradication strategy. Lancet. 2008, 371: 1633-1635. 10.1016/S0140-6736(08)60424-9.
Roberts L, Enserink M: Did they really say...eradication?. Science. 2007, 318: 1544-1545. 10.1126/science.318.5856.1544.
Chitnis N, Smith TA, Steketee R: A mathematical model for the dynamics of malaria in mosquitoes feeding on a heterogeneous host population. J Biol Dynamics. 2008, 2: 259-285. 10.1080/17513750701769857.
Ferguson HM, Dornhaus A, Beeche A, Borgemeister C, Gottlieb M, Mulla MS, Gimnig JE, Fish D, Killeen GF: Ecology: a prerequisite for malaria elimination and eradication. PLoS Med. 2010, 7: e1000303-10.1371/journal.pmed.1000303.
Gillies MT: A modified technique for age grading populations of Anopheles gambiae. Ann Trop Med Parasitol. 1958, 58: 261-273.
Gillies MT: Studies on the dispersion and survival of Anopheles gambiae in East Africa, by means of marking and release experiments. Bull Entomol Res. 1961, 52: 99-127. 10.1017/S0007485300055309.
Gillies MT, Wilkes TJ: A study of the age-composition of populations of Anopheles gambiae Giles and A. funestus Giles in North-Eastern Tanzania. Bull Entomol Res. 1965, 56: 237-262. 10.1017/S0007485300056339.
Beier JC: Frequent blood-feeding and restrictive sugar-feeding behavior enhance the malaria vector potential of Anopheles gambiae s.l. and An. funestus (Diptera:Culicidae) in western Kenya. J Med Entomol. 1996, 33: 613-618.
Killeen GF, Tami A, Kihonda J, Okumu FO, Kotas ME, Grundmann H, Kasigudi N, Ngonyani H, Mayagaya V, Nathan R, Abdulla S, Charlwood JD, Smith TA, Lengeler C: Cost-sharing strategies combining targeted public subsidies with private-sector delivery achieve high bednet coverage and reduced malaria transmission in Kilombero Valley, southern Tanzania. BMC Infect Dis. 2007, 7: 121-10.1186/1471-2334-7-121.
Hii JLK, Smith T, Vounatsou P, Alexander N, Mai A, Ibam E, Alpers MP: Area effects of bednet use in a malaria-endemic area in Papua New Guinea. Trans R Soc Trop Med Hyg. 2001, 95: 7-13. 10.1016/S0035-9203(01)90315-3.
Lindsay SW, Emerson PM, Charlwood JD: Reducing malaria transmission by mosquito-proofing homes. Trends Parasitol. 2002, 18: 510-514. 10.1016/S1471-4922(02)02382-6.
Sota T, Mogi M: Effectiveness of zooprophylaxis in malaria control: a theoretical inquiry with a model for mosquito populations with two bloodmeal hosts. Med Vet Entomol. 1989, 3: 337-345. 10.1111/j.1365-2915.1989.tb00240.x.
Yakob L, Yan G: Modelling the effects of integrating larval habitat source reduction and insecticide treated nets for malaria control. PLoS One. 2009, 4: e6921-10.1371/journal.pone.0006921.
Gu W, Regens JL, Beier JC, Novak RJ: Source reduction of mosquito larval habitats has unexpected consequences on malaria transmission. Proc Natl Acad Sci USA. 2006, 103: 17560-17563. 10.1073/pnas.0608452103.
Chitnis N, Schapira A, Smith T, Steketee R: Comparing the effectiveness of malaria vector-control interventions through a mathematical model. Am J Trop Med Hyg. 2010, 83: 230-240. 10.4269/ajtmh.2010.09-0179.
Yakob L, Dunning R, Yan G: Indoor residual spray and insecticide-treated bednets for malaria control: theoretical synergisms and antagonisms. J R Soc Interface. 2010, 8: 799-806.
Kleinschmidt I, Schwabe C, Shiva M, Segura JL, Sima V, Mabunda SJA, Coleman M: Combining Indoor Residual Spraying and Insecticide-Treated Net Interventions. Am J Trop Med Hyg. 2009, 81: 519-524.
Binka FN, Indome F, Smith T: Impact of spatial distribution of permethrin-impregnated bed nets on child mortality in rural Northern Ghana. Am J Trop Med Hyg. 1998, 59: 80-85.
Gimnig JE, Kolczak MS, Hightower AW, Vulule JM, Schoute E, Kamau L, Phillips-Howard PA, ter Kuile FO, Nahlen BL, Hawley WA: Effect of permethrin-treated bed nets on the spatial distribution of malaria vectors in western Kenya. Am J Trop Med Hyg. 2003, 68 (Suppl 4): 115-120.
Gimnig JE, Vulule JM, Lo TQ, Kamau L, Kolczak MS, Phillips-Howard PA, Mathenge EM, ter Kuile FO, Nahlen BL, Hightower AW, Hawley WA: Impact of permethrin-treated bed nets on entomologic indices in an area of intense year-round malaria transmission. Am J Trop Med Hyg. 2003, 68 (Suppl 4): :16-22.
Howard SC, Omumbo J, Nevill CG, Some ES, Donnelly CA, Snow RW: Evidence for a mass community effect of insecticide treated bednets on the incidence of malaria on the Kenyan coast. Trans R Soc Trop Med Hyg. 2000, 94: 357-360. 10.1016/S0035-9203(00)90103-2.
Maxwell CA, Msuya E, Sudi M, Njunwa KJ, Carneiro IA, Curtis CF: Effect of community-wide use of insecticide-treated nets for 3-4 years on malarial morbidity in Tanzania. Trop Med Int Health. 2002, 7: 1003-1008. 10.1046/j.1365-3156.2002.00966.x.
Habicht JP, Victora CG, Vaughan JP: Evaluation designs for adequacy, plausibility and probability of public health programme performance and impact. Int J Epidemiol. 1999, 28: 10-18. 10.1093/ije/28.1.10.
Ranson H, Abdallah H, Badolo A, Guelbeogo WM, Kerah-Hinzoumbe C, Yangalbe-Kalnone E, Sagnon N, Simard F, Coetzee M: Insecticide resistance in Anopheles gambiae: data from the first year of a multi-country study highlight the extent of the problem. Malar J. 2009, 8: 299-10.1186/1475-2875-8-299.
Griffin JT, Hollingsworth TD, Okell LC, Churcher TS, White M, Hinsley W, Bousema T, Drakeley CJ, Ferguson NM, Basáñez MG, Ghani AC: Strategies towards Plasmodium falciparum malaria elimination in Africa using currently available tools. PLoS Medicine. 2010, 7: e1000324-10.1371/journal.pmed.1000324.
Pates HV, Line JD, Keto AJ, Miller JE: Personal protection against mosquitoes in Dar es Salaam, Tanzania, by using a kerosene oil lamp to vaporize transfluthrin. Med Vet Entomol. 2002, 16: 277-284. 10.1046/j.1365-2915.2002.00375.x.
Kimani EW, Vulule JM, Kuria IW, Mugisha F: Use of insecticide-treated clothes for personal protection against malaria: a community trial. Malar J. 2006, 5: 63-10.1186/1475-2875-5-63.
Macintyre K, Sosler S, Letipila F, Lochigan M, Hassig S, Omar SA, Githure J: A new tool for malaria prevention? Results of a trial of permethrin-impregnated bedsheets (shukas) in an area of unstable transmission. Int J Epidemiol. 2003, 32: 157-160. 10.1093/ije/dyg019.
Rowland M, Durrani N, Kenward M, Mohammed N, Urahman H, Hewitt S: Control of malaria in Pakistan by applying deltamethrin insecticide to cattle: a community-randomised trial. Lancet. 2001, 357: 1837-1841. 10.1016/S0140-6736(00)04955-2.
Lwetoijera DW, Sumaye RD, Madumla EP, Kavishe DR, Mnyone LL, Russell TL, Okumu FO: An extra-domiciliary method of delivering entomopathogenic fungus, Metharizium anisopliae IP 46 for controlling adult populations of the malaria vector, Anopheles arabiensis. Parasit Vectors. 2010, 3: 18-10.1186/1756-3305-3-18.
Okumu FO, Madumla EP, John AN, Lwetoijera DW, Sumaye RD: Attracting, trapping and killing disease-transmitting mosquitoes using odor-baited stations - The Ifakara Odor-Baited Stations. Parasit Vectors. 2010, 3: 12-10.1186/1756-3305-3-12.
Korenromp EL, Miller J, Cibulskis RE, Kabir Cham M, Alnwick D, Dye C: Monitoring mosquito net coverage for malaria control in Africa: possession vs. use by children under 5 years. Trop Med Int Health. 2003, 8: 693-703. 10.1046/j.1365-3156.2003.01084.x.
Miller JM, Korenromp EL, Nahlen BL, R WS: Estimating the number of insecticide-treated nets required by African households to reach continent-wide malaria coverage targets. JAMA. 2007, 297: 2241-2250. 10.1001/jama.297.20.2241.
Beier JC, Killeen GF, Githure J: Short report: Entomologic inoculation rates and Plasmodium falciparum malaria prevalence in Africa. Am J Trop Med Hyg. 1999, 61: 109-113.
Hay SI, Rogers DJ, Toomer JF, Snow RW: Annual Plasmodium falciparum entomological inoculation rates across Africa: literature survey, internet access and review. Trans R Soc Trop Med Hyg. 2000, 94: :113-127. 10.1016/S0035-9203(00)90246-3.
Smith DL, Dushoff J, Snow RW, Hay SI: The entomological inoculation rate and Plasmodium falciparum infection in African children. Nature. 2005, 438: 492-495. 10.1038/nature04024.
Smith TA, Maire N, Dietz K, Killeen GF, Vounatsou P, Molineaux L, Tanner M: Relationship between entomologic inoculation rate and the force of infection for Plasmodium falciparum malaria. Am J Trop Med Hyg. 2006, 75 (Supplement 2): 11-18.
Ross A, Killeen GF, Smith TA: Relationships of host infectivity to mosquitoes and asexual parasite density in Plasmodium falciparum. Am J Trop Med Hyg. 2006, 75 (Suppl 2): 32-37.
Saul AJ, Graves PM, Kay BH: A cyclical feeding model for pathogen transmission and its application to determine vectorial capacity from vector infection rates. J Appl Ecol. 1990, 27: 123-133. 10.2307/2403572.
Guillet P, Alnwick D, Cham MK, Neira M, Zim M, Heymann D, Mukelebai K: Long-lasting treated mosquito nets: A breakthrough in malaria prevention. Bull World Health Organ. 2001, 79: 998-
Gillies MT: Studies in house-leaving and outside resting of Anopheles gambiae Giles and Anopheles funestus Giles in East Africa. Bull Entomol Res. 1954, 45: 375-387. 10.1017/S000748530002719X.
Charlwood JD, Smith T, Kihonda J, Heiz B, Billingsley PF, Takken W: Density independent feeding success of malaria vectors (Diptera: Culicidae) in Tanzania. Bull Entomol Res. 1995, 85: 29-35. 10.1017/S0007485300051981.
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.
We thank Mr S Kiware for his careful review of the model. We also thank Dr N Achee, Dr J Grieco, Dr R N'Guessan, Dr K Aultman and Prof SW Lindsay for valuable discussions which have contributed to the formulation of the model and the interpretation of the results. This work was funded in part by the Bill & Melinda Gates Foundation through award numbers 45114 (Malaria Transmission Consortium), 51431 (Replacing DDT: Rigorous Evaluation of Spatial Repellents for the Control of Vector Borne Diseases) and 39777.01 (A stochastic simulation platform for predicting the effects of different malaria intervention strategies). This work is made possible by the generous support of the American people through the United States Agency for International Development (USAID award number 621-A-00-08-0007-00). The contents are the responsibility of GFK, NC, SJM and FOO and do not necessarily reflect the views of USAID or the United States Government.
The authors declare that they have no competing interests. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
All authors formulated the research question and developed the conceptual basis of the model. GFK drafted the model formulation and manuscript in consultation with FOO, SJM and NC. The contents are the responsibility of GFK, NC, SJM and FOO and do not necessarily reflect the views of USAID, the United States Government or the Bill & Melinda Gates Foundation. All authors have read and approved the final version of the manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.