Using evolution to generate sustainable malaria control with spatial repellents

Abstract
eLife digest
Introduction
Results
Discussion
Appendix 1
Appendix 2
References
Article and author information
Metrics

Abstract

Evolution persistently undermines vector control programs through insecticide resistance. Here we propose a novel strategy which instead exploits evolution to generate and sustain new control tools. Effective spatial repellents are needed to keep vectors out of houses. Our approach generates such new repellents by combining a high-toxicity insecticide with a candidate repellent initially effective against only part of the vector population. By killing mosquitoes that enter treated properties the insecticide selects for vector phenotypes deflected by the repellent, increasing efficacy of the repellent against the target vector population and in turn protecting the insecticide against the spread of insecticide resistance. Using such evolved spatial repellents offers an evolutionarily sustainable, ‘double-dip’ system of disease control combining mortality and repellence. We formalize this idea using models which explore vector population genetics and disease transmission probabilities and show that using evolved spatial repellents is theoretically achievable, effective and sustainable.

https://doi.org/10.7554/eLife.15416.001

eLife digest

Many of the mosquito species that transmit malaria have evolved to bite humans indoors at night, and therefore health programs target them using insecticides sprayed on surfaces inside people’s homes. This strategy, however, stops working when mosquito populations evolve to resist the insecticide used, either because they are immune to its poisonous effects or because they change their behaviour to avoid it. Consequently, there is now a need to develop alternative strategies to control mosquitoes that are more sustainable in the face of evolution. One possibility is repellents that keep mosquitoes out of homes.

Lynch and Boots have now asked whether evolution could be used to create effective repellents from substances that initially repel only part of the mosquito population by pairing them with lethal insecticides sprayed inside people’s homes. Mathematical models showed that, before insecticide resistance becomes widespread, this “evolved repellence” approach could reduce the spread of malaria by a similar amount to using insecticides alone. This was particularly true if the models considered that, as well as surviving to give fewer infectious bites, repelled infectious mosquitoes may be less likely to transmit malaria with each feed, for example if they feed more on livestock rather than humans.

The models of Lynch and Boots also show that that the success of the evolved repellence concept in a given location depends on a number of factors. The proportion of the starting mosquito population that is repelled or resistant can have a large effect. Similarly, success will also depend on how likely normal, repelled and insecticide-resistant mosquitoes are to reproduce successfully. These values can be influenced by the choice of insecticide and repellent and how the chemicals are applied. Lynch and Boots show that swapping insecticides can allow an evolved repellent to be established where it would otherwise not succeed. Also, the spread of resistance to the paired insecticide is slowed or prevented when the mosquito population evolves to be repelled.

Practical laboratory and field- work is now needed to build on this theoretical groundwork and to determine suitable locations and application strategies to exploit this concept as a way to sustainably reduce the spread of malaria in the future.

https://doi.org/10.7554/eLife.15416.002

Introduction

Vector-borne diseases including malaria and dengue remain a major burden to human health despite decades of funding targeting their eradication. One child per minute still dies of malaria globally (WHO, 2015a), and dengue infects an estimated four hundred million individuals per year (Bhatt et al., 2013; Murray et al., 2013). One of the key reasons for this has been the impact of evolution on our control tools, with both theory and experience pointing to the inevitability of resistance to conventional insecticides (Brogdon and McAllister, 1998; Hemingway et al., 2002; Nauen, 2007; Asidi et al., 2012; Hemingway and Ranson, 2000; Ndiath et al., 2012). It is clear that mosquito populations have considerable potential to evolve and that to achieve sustainable control, it would be prudent to develop approaches that take evolution into account. Very fundamentally, there is a need for evolutionarily rational disease intervention (Boots, 2015). Here, we highlight a potential intervention strategy in which evolution in the vector may be utilized to improve rather than reduce disease control. Specifically, we argue that combining initially limited repellents with toxic insecticides has the potential to select for effective repellence and to lead to evolutionarily sustainable control using evolved spatial repellents (‘ESR’s).

It is important to understand that substances can affect vectors in multiple ways (Kennedy, 1947; Miller et al., 2009) and there is considerable overlap in the terms used to describe their different actions. We adopt the terminology of Greico et al. (Grieco et al., 2007) who characterize three modes of action for the insecticides used in vector control programs; toxicity, contact-irritance and spatial repellence. Toxicity is simply the capacity to generate mortality in mosquitoes that contact the insecticide. Contact-irritance is the repellence of mosquitoes when they contact the insecticide, and spatial repellence acts at a distance, deflecting mosquitoes before they contact a treated surface. Specifically, we are considering repellence that acts to protect people indoors by preventing malaria vectors from entering treated properties, rather than more localized repellent effects that might deflect them away from individuals outdoors or from treated surfaces indoors. The two approaches are distinct, localized repellence from treated bednets, for example, provides protection only for individuals with nets, and only at times of day when they are under the nets.

In the absence of resistance, DDT, the most effective chemical used to date for the control of malaria vectors, has both high toxicity and high spatial repellence (Kennedy, 1947; Grieco et al., 2007). Work inspired by the historic efficacy of DDT as a tool in combating malaria has highlighted the potential role of spatial repellence as a means of reducing transmission (Achee et al., 2012; Loyola et al., 1990; Roberts et al., 2000; Roberts, 2010; Roberts et al., 2000; Curtis, 2002; Roberts and Alecrim, 1991). The malaria vectors targeted by public health campaigns using bednets and indoor residual spraying of insecticides (IRS) typically feed indoors between dusk and dawn. Vector control has therefore focused on exploiting this behaviour to deliver lethal control measures against indoor-feeding mosquitoes (Yakob et al., 2011; WHO, 2015b), with outdoor biting viewed as unwanted behavioural resistance (Bradley et al., 2012; Russell et al., 2011; Reddy et al., 2011; Cooke et al., 2015; Mouchet et al., 1963; Gatton et al., 2013; Killeen et al., 2011). However, it is increasingly being acknowledged that a switch away from indoor biting, by reducing vector fitness and/or by creating a less favourable context for transmission, may offer public health benefits as an end in itself, consistent with the idea that the efficacy of DDT may be partly derived from its repellent rather than solely its toxic effects (Roberts et al., 2000; Roberts, 1993; Grieco et al., 2000). The use of spatial repellents as a tool to maintain vector-free homes is therefore now being actively investigated as a means of reducing transmission of malaria and dengue (Achee et al., 2012; N'Guessan et al., 2006; Ogoma et al., 2014; University of Notre Dame, 2016). This can be viewed as an alternative means to achieve the same benefits as the application of physical barriers to prevent vector access to homes (Roll Back Malaria, 2015; Patrick, 1900; Menger et al., 2016). Syafruddin et al. have conducted a proof-of-concept double-blind placebo-controlled field trial using spatial repellents disbursed from burning coils. This trial recorded an approximately 50% reduction in malaria infections in repellent-protected properties (Syafruddin et al., 2014). They emphasize the need for repellents able to achieve the same effects without combustion in order to develop practical public health interventions. Given the benefits of insect repellence in terms of nuisance reduction as well as disease control, and the consequent long-standing interest in identifying safe and effective repellents, the limited choice and efficacy of currently available compounds suitable for this purpose suggests that finding suitable new repellent substances is a particularly challenging task. Although many promising potential repellents are being evaluated for personal protection, for localized outdoor and indoor protection, and for the property-scale protection relevant for ESRs (Syafruddin et al., 2014; Govella et al., 2015; Abiy et al., 2015; Revay et al., 2013; Müller et al., 2009), sustained, effective action at a distance is still challenging to achieve, and still often requires active dispersal through combustion or powered devices (Revay et al., 2013; Müller et al., 2009). DDT itself carries a reputational burden from its historic overuse, such that environmental concerns make its continued use problematic and its continued availability uncertain.

We propose that for the vectors currently targeted by IRS campaigns, there is an opportunity to exploit the evolutionary processes generated by using lethal insecticides to create effective new spatial repellents from compounds that initially repel only a fraction of the vector population. Specifically, we propose a combination of a spatial repellent, which deters mosquitoes from entering buildings, with a low-contact-repellence high-toxicity insecticide, which kills those which do enter, leading to a ‘double-dip’ system of disease control. If failure to enter buildings is seen as a method of transmission reduction in its own right, then the repellent provides transmission reduction by deflecting mosquitoes. For mosquitoes that are not deflected, the mortality imposed by contact with the insecticide will provide transmission reduction in the same manner as current insecticide-only IRS control methods. Provided that, in a given context, the fitness cost of being deflected is less than the fitness cost of being susceptible to the insecticide, there is potential for selection to favour an increase in the proportion of deflected individuals in the treated population. If deflection is viewed as a form of ‘behavioural resistance’, this system actively exploits ‘resistance’ evolution, as mortality generated by the insecticide serves to select for phenotypes that are deflected by the repellent. Candidates for use as spatial repellents thus only need initially to repel a small proportion of a mosquito population, since in this instance, for once, evolution will work to enhance the efficacy of a disease-control measure. Here, we explore the feasibility of this approach using evolutionary models that explore a wide range of possible vector characteristics and disease parameters.

Modeling

If deflecting mosquitoes from accessing humans indoors is in itself an effective means of reducing transmission, as proposed by Achee et al. (2012), then selection for mosquitoes which are repelled from treated houses could serve to generate new public health tools. Critically, selection would depend on the relative fitness of mosquitoes that are deflected away from buildings compared to those entering buildings. This in turn would depend on the proportion of properties treated with insecticides, the susceptibility or resistance (physiological) of mosquitoes contacting indoor insecticides and the fitness costs associated with being deflected away from the sleeping indoor hosts which vector species have evolved to exploit. However, the fitness costs of being deflected from human dwellings are difficult to determine directly. Following a well-established history of mathematical modeling to explore issues relating to the evolution of resistance in malaria vector populations (for example [Rosenheim and Tabashnik, 1990; Le Menach et al., 2007; Mandal et al., 2011; Georghiou and Taylor, 1977; White et al., 2014]), we have therefore developed an analysis to explore the possible outcome across a range of fitness scenarios, using a feeding-cycle based, two-locus, bi-allelic population genetics model, capturing the mosquito life-history characteristics of one initial mating and pre-adult development period as well as the Plasmodium development period in infected mosquitoes before transmission is possible.

As well as the spread of resistance and deflection alleles over time, our model also tracks adult population size and the expected number of infectious bites given by the population during each modeled time unit. Each modeled time unit corresponds to the length and reproductive outcomes of a single feeding cycle, an approach used in previous models of vector population genetics (Lynch et al., 2012; Read et al., 2009). For details of the model see Appendix 1. To minimize the sensitivity of the model results to specific parameter values, we frame our key disease control results in terms of the proportionate difference between the model’s calculated infectious bite values for a given set of intervention assumptions and those assuming no intervention, minimising the impact of parameter values that are unaffected by the intervention. The impact of deflection on malaria prevalence is determined by the proportion of mosquitoes deflected by a repellent and the probability (compared to non-deflected mosquitoes) that they will then acquire and transmit a Plasmodium infection. The effect on Plasmodium transmission of deflecting vectors to outdoor biting has not been definitively measured in the field, we therefore consider a wide range of per-feed probabilities of Plasmodium transmission to deflected mosquitoes compared to the 4% probability assumed for indoor feeds. Reductions in this parameter are intended to represent the effects of all potential sources of reduced transmission to the vector, including deflection to non-human and therefore non-infectious hosts. In terms of model results, proportionate reductions in this parameter will have the same effect as proportionate changes in the probability that an infectious mosquito which survives to feed will give an infectious bite to a human host. As such, the reductions explored can be interpreted as the product of the proportionate reductions in transmission to and from outdoor-feeding vectors.

Our key model assumptions. (1) Physiological insecticide resistance is controlled by a single locus bi-allelic autosomal gene, with the resistance allele being completely dominant to the susceptibility allele. (2) Deflection by a given spatial repellent is controlled by a single locus bi-allelic autosomal gene, with the deflection allele being completely dominant to the non-deflection allele. (3) Deflected vectors are assumed not to come into contact with the insecticide used in association with the ESR; therefore, mosquitoes that have phenotypes which combine deflection and resistance will not experience any of the fitness benefits associated with resistance if the ESR and insecticide are always present together. (4) The resistance and deflection loci are not linked and re-assort randomly. (5) The genotypes determining adult resistance and deflection phenotypes do not affect the probability of juvenile survival from egg to adult. (6) Mating is random and females mate once, as newly emerged adults, with males in their cohort. (7) Juvenile density dependence means that variation in the absolute number of eggs produced by the adult population does not materially change the rate at which new adults join the population.

Whilst we cannot predict the form that the genetic determinants of resistance, behavioural or otherwise, may take, there are examples of single-locus insecticide-resistance genes, including the knockdown resistance (kdr) alleles that provide resistance to DDT and pyrethroids (Ndiath et al., 2012; Chandre et al., 2000; Jones et al., 2012; Dabiré et al., 2012). Whilst the genetic basis of deflection behaviour is unknown and may often be more complex than that of insecticide resistance, it is parsimonious to model this process initially by assuming simple single genes that determine the likelihood of such responses. The impact on the predictions of continuous traits, or, should suitable data become available, of specific more complex genetic assumptions, can be incorporated into future work. The model considers four possible phenotypes, as shown in Table 1, with associated genotypes (resistance alleles represented by R, and deflection alleles by D).

Table 1

Phenotype definitions and characteristics.

https://doi.org/10.7554/eLife.15416.003

Phenotype	Fitness	Resistant	Deflected	Genotypes
Susceptible	$F_{S}$	No	No	rr\dd
Resistant	$F_{R}$	Yes	No	Rr\dd RR\dd
Deflected and not resistant	$F_{D}$	No	Yes	rr\Dd rr\DD
Deflected and resistant	$F_{R D}$	Yes	Yes	Rr\Dd Rr\DD RR\Dd RR\DD

The average fitness of offspring into which deflection alleles are inherited, ${\bar{F}}_{D}$ , is

{\bar{F}}_{D} = F_{S} + \frac{[D r] [- R] + [D R]}{[D -]} (F_{R D} - F_{S}) + \frac{[D r] [- r]}{[D -]} (F_{D} - F_{S})

With $[d r], [d R], [D r], [D R], [d -]$ and $[D -]$ representing, in the zygote genotypes for the population at a given time point, the proportion of alleles at the deflection locus which are non-deflection alleles paired with susceptible alleles, non-deflection alleles paired with resistant alleles, deflection alleles paired with susceptible alleles, deflection alleles paired with resistant alleles, non-deflection alleles paired with any resistance allele, and deflection alleles paired with any resistance allele, respectively, assuming the same proportions in gametes of mating males and newly emerged females.

The average fitness of offspring into which non-deflection alleles are inherited, ${\bar{F}}_{d}$ , is

{\bar{F}}_{d} = F_{S} + [d R] (1 + \frac{[d r]}{[d -]}) (F_{R} - F_{S}) + \frac{[d r] [D r]}{[d -]} (F_{D} - F_{S}) + ([D R] + \frac{[d R] [D r]}{[d -]}) (F_{R D} - F_{S})

In order for the proportion of deflection alleles in the population to increase, we need the average fitness of the offspring into which deflection alleles are inherited to be greater than the average fitness of the offspring into which non-deflection alleles are inherited. This is true when the following inequality applies:

\begin{array}{l} {\bar{F}}_{D} > {\bar{F}}_{d} \leftrightarrow \frac{[d r]^{2}}{[d -]} (F_{D} - F_{S}) > [d R] (1 + \frac{[d r]}{[d -]}) (F_{R} - F_{D}) \\ + ([D R] (1 - \frac{1}{[D -]}) + [D r] (\frac{[d R]}{[d -]} - \frac{[- R]}{[D -]})) (F_{R D} - F_{D}) \end{array}

The equivalent expression for spread of the resistance allele is:

\begin{aligned} {\bar{F}}_{R} & > {\bar{F}}_{r} \leftrightarrow [r d] (F_{R} - F_{S}) + [r D] (F_{R} - F_{D}) + ([r D] + \frac{[R D] [- d]}{[R -]} - \frac{[R d] [r D]}{[r -]}) (F_{R D} - F_{R}) \\ > \frac{[r d] [r D]}{[r -]} (F_{D} - F_{S}) \end{aligned}

See Appendix 2 for derivation of Equations 1 and 2 and Expressions 1 and 2.

We assume that ESR is only relevant where the fitness cost of being susceptible to an IRS insecticide is greater than the fitness cost of being deflected by an ESR, so $F_{D} > F_{S}$ , requiring that the mortality associated with a mosquito that has a susceptible phenotype entering an insecticide-treated property is greater than that associated with a mosquito being deflected from a property. From Expression 1, it can be seen that the spread of the deflection allele in the vector population will be favoured by maximising the fitness difference between susceptible and deflected phenotypes (increasing the value of the left-hand side of the expression), and by maximising the fitness of deflected relative to resistant phenotypes (reducing the value of the right-hand side of the expression). Avoiding the use of ESR in properties without insecticide, minimising $Y_{3}$ (Table 2), and using ESR in all insecticide-treated properties, minimising $Y_{2}$ , improves the survival of deflected phenotypes (given $B < I$ ) without affecting the survival probability of susceptible or resistant mosquitoes, enhancing the desired fitness relationships and hence favouring the spread of deflection alleles in the vector population and the initial establishment of a new evolved spatial repellent.

Table 2

Feeding related survival probabilities.

https://doi.org/10.7554/eLife.15416.004

		Feeding-related survival probabilities
Proportion of properties		Y₁	Y₂	Y₃	Y₄
Phenotype	Baseline fitness adjustment	Untreated property	Insecticide only	ESR only	Insecticide and ESR	Average survival
Susceptible		$U$	$U - I$	$U$	$U - I$	$U - I (Y_{2} + Y_{4})$
Resistant	COR₁	$U$	$U$	$U$	$U$	$U$
Resistant and deflected	COR₂	$U$	$U$	$U - B$	$U - B$	$U - B (Y_{3} + Y_{4})$
Deflected		$U$	$U - I$	$U - B$	$U - B$	$U - I Y_{2} - B (Y_{3} + Y_{4})$

$U$ =no-treatment survival, $I$ =survival reduction caused by insecticide in susceptible mosquitoes, $B$ =survival reduction caused by deflection from protected building, $Y_{i}$ =proportion of properties in each treatment category, COR₁=fitness cost of resistance experienced by resistant non-deflected phenotypes, COR₂=fitness cost of resistance experienced by resistant deflected phenotypes.

From Expression 1, it can be seen that the spread of deflection alleles is dependent not only upon the relative fitness values of the different phenotypes but also upon there being sufficiently low initial levels of resistance alleles in the population. Since the genotype proportions will change over time, this is a dynamic relationship. In order for the $D$ allele to spread at all, initial allele proportions and fitness relationships must comply with the inequality in Expression 1, but the spread of the resistance allele over time may eventually reverse the relationship, so that the proportion of non-deflection alleles will begin to increase instead. The rate of spread of the resistance allele will also determine whether the deflection allele will spread and be sustained in the population. From Expression 2, which shows the conditions necessary for the resistance allele to spread, it can be seen that the fitness differential between resistant and susceptible phenotypes and that between resistant and deflected-resistant phenotypes, by helping to determine whether the resistance alleles spread, are also determinants of whether the deflection allele, and hence deflected phenotypes, will spread and be maintained in the population. The spread of the deflection allele when resistance alleles are present in the population is also critically determined by the initial proportion of deflection alleles at the deflection locus and by the proportions of deflection and non-deflection alleles that are paired with resistance alleles. These interactions are explored in Figures 1–4.

Figure 1

Download asset Open asset

Spread of deflection in a population over time.

(i) Phenotype and (ii) genotype proportions over time for a population subject to insecticide interventions applied in combination with an ESR. Illustrating (a) long-term establishment of ESR, (b) transient establishment of ESR, and (c) failure to establish ESR. The parameter values used to generate the plots in panels (a) are: 20% per cycle survival of susceptible phenotypes, 60% per cycle survival of resistant phenotypes 45% per cycle survival of deflected phenotypes, 0.5% initial proportion of resistance alleles and 25% initial proportion of deflection alleles. For the panels in (b) and (c) , the per cycle survival of deflected phenotypes is reduced to 40%. For panel (c), other parameters are also amended to 30% per cycle survival for susceptible phenotypes, 5% initial prevalence of resistance alleles and 10% initial prevalence of deflection alleles.

https://doi.org/10.7554/eLife.15416.005

In the absence of resistance, Expression 1 reduces to ${\bar{F}}_{D} > {\bar{F}}_{d} \leftrightarrow [d -] (F_{D} - F_{S}) > 0$ , so deflection would be expected to spread to fixation provided that deflected phenotypes have greater fitness than non-deflected phenotypes. If deflection reaches fixation, then Expression 2 reduces to ${\bar{F}}_{R} > {\bar{F}}_{r} \leftrightarrow [r D] (F_{R D} - F_{D}) > 0$ , and resistance will only spread if resistant-deflected phenotypes are fitter than deflected phenotypes. In this instance, the strategy of using ESR in all insecticide-treated properties would prevent further spread of resistance because deflected mosquitoes would not enter any insecticide-treated properties and would experience no benefits from resistance. hence the fitness of resistant deflected phenotypes will be the same as or lower than that of non-resistant deflected phenotypes (depending upon cost of resistance), giving $F_{R D} \leq F_{D}$ .

Results

We carried out a numerical analysis using the model to explore the establishment over time of a new ESR and the associated population-level changes in infectious bite rate. Assumed baseline parameter values are: (i) time from egg laying to adult emergence equivalent to the length of three gonotrophic cycles, (ii) probability per feed that a non-deflected mosquito will acquire a Plasmodium infection is 4%, (iii) probability per feed that an infectious mosquito gives an infectious bite on a human host is 80%, (iv) time to infectiousness of Plasmodium infection in the vector is approximately equivalent to the length of three gonotrophic cycles. Where not stated otherwise, we use 20% survival of susceptible phenotypes per cycle, 60% survival of resistant phenotypes per cycle, 45% survival of deflected phenotypes per cycle, 0.5% initial proportion of resistance alleles and 25% initial proportion of deflection alleles.

The first question considered is whether, and under what circumstances, selection could generate an effective spatial repellent from a substance that initially repelled only a part of the population. Consistent with Expression 1, we found that the spread of a deflection allele through the population depended on the initial proportions of deflection and resistance alleles in the population and the fitness differentials between susceptible, resistant, deflected and resistant deflected phenotypes.

As illustrated in Figure 1, for some combinations of fitness values and initial allele proportions, the deflection allele spreads rapidly to near-fixation, and remains consistently at that level for at least 300 cycles (panels labelled ‘a’). In other cases, the deflection allele spreads initially, but falls away within 300 cycles as the resistance allele spreads (panels labelled ‘b’), and in some cases the deflection allele shows only minimal spread before being lost as resistance spreads (panels labelled ‘c’). For cases like that in ‘b’, where deflection spreads initially but then falls away, we considered the outcome if the insecticide used is swapped for an alternative, for which resistance alleles are still relatively rare, whilst the deflection allele is close to its peak prevalence. In some cases this allows a ‘ratchet’ effect, whereby the deflection allele is able to spread and reach sustained high levels.

There is presently little or no direct information available about the fitness costs of a switch to outdoor biting. Furthermore, this would be expected to vary with mosquito species, the degree of anthropophilly, the type and accessibility of outdoor hosts, and various other factors. The initial proportion of deflection and resistance alleles in the population will depend upon the choice of insecticide and ESR, but will also be expected to vary between specific populations. We therefore carried out analyses for a range of parameter values, the results of which are summarized in Figure 2.

Figure 2 with 1 supplement see all

Download asset Open asset

Combinations of per cycle survival values for deflected and non-deflected resistant phenotypes, which support the spread and maintenance of deflected phenotypes in the population.

Grid plots indicating which combinations of resistant (x-axis) and deflected (y-axis) phenotype per-cycle survival values (in 1% increments) give rise to a population comprising at least 80% deflected phenotypes after 300 modeled time periods. When this is achieved directly, the applicable square is bright green. Dark green squares indicate combinations for which the required outcome can be achieved via a ‘ratchet’ where the initial paired insecticide is swapped once for a new insecticide, with allele proportions at the time of the swap assumed to be 0.5% resistance alleles and the maximum percentage of deflection alleles achieved whilst using the first insecticide. Results are calculated for 1% increments in each survival value. Gridlines and diagonals are to aid visual location of results on the grid. The baseline parameters (panel i) are: 20% per cyclesurvival of susceptible phenotypes; resistant deflected phenotypes have the same survival probability as non-resistant deflected phenotypes; 0.5% initial prevalence of resistance alleles; and 25% initial prevalence of deflection alleles. Parameter values for panels (ii) to (iv) differ from the baseline values as follows: panel (ii) 10% initial prevalence of deflection alleles; panel (**iii**) 2% initial prevalence of resistance alleles; and panel (iv) 60% per cycle survival of susceptibles .

https://doi.org/10.7554/eLife.15416.006

A comparison of panels (i) to (iv) in Figure 2 shows that the fitness values for all phenotypes, and the initial prevalence of deflection and resistance alleles all affect the potential for deflection alleles to spread and be maintained in the population, consistent with the relationships shown in Expression 1. For example, with the baseline parameter values, in a context in which insecticide-resistant phenotypes have average per-cycle survival of 60% and deflected phenotypes have per-cycle survival of 45%, an ESR introduced at a time when the prevalence of resistance and deflection alleles are 0.5% and 25% will become established, with the deflection allele spreading and deflected phenotypes comprising more than 80% of the population 300 cycles after introduction (panel i), as shown by the light green square for 45% deflected survival with 60% resistant survival. However, if the initial prevalence of deflection alleles is only 10%, then with the same per-cycle survival rates, deflection alleles spread but are not sustained. Deflection can, however, be established by replacement of the initial insecticide whilst deflection alleles are at their maximum prevalence (panel ii), as shown by the dark green square for 45% deflected survival with 60% resistant survival.

From comparison of Figure 2 panel (i) and Figure 3, it can be seen that if a cost of resistance affects resistant deflected phenotypes this serves to increase the range of deflection and resistance fitness combinations for which deflection alleles can spread and be maintained.

Figure 3

Download asset Open asset

Effect of cost of resistance on resistant and deflected per-cycle survival combinations which support the spread and maintenance of deflected phenotypes in the population.

Combinations of phenotype survival values which can result in more than 80% of the population having deflected phenotypes after 300 cycles, assuming various costs of resistance. Colours indicate the lowest cost of resistance (COR) incurred by deflected resistant phenotypes which achieves the threshold 80% deflection phenotypes in the population after 300 cycles (without assuming any ratchet). COR here represents the reduction in per cycle survival of deflected resistant phenotypes arising as result of having a resistant phenotype, so the resistant deflected phenotype has per-cycle survival equal to that for the deflected phenotype (y-axis) less the applicable COR. The survival values shown for the resistant phenotype (x-axis) are those after taking account of any cost of resistance that affects non-deflected resistant phenotypes. The baseline parameters are: 20% per cycle survival of susceptible phenotypes; 0.5% initial prevalence of resistance alleles; and 25% initial prevalence of deflection alleles.

https://doi.org/10.7554/eLife.15416.008

Once the ESR is established, the impact of the ESR–insecticide combination treatment on the transmission of Plasmodium will critically depend on the reduced vectorial capacity arising from the exclusion of mosquitoes from treated properties. This will have two components. Reduced survival of deflected mosquitoes will reduce the population of adult mosquitoes, the probability of infected mosquitoes surviving to give an infectious bite, and the number of bites that an infectious mosquito will survive to give. In addition, anything which reduces the probability of transmission from human host to a feeding mosquito, or from a feeding mosquito to a human host, will enhance the reduction in transmission arising from deflection away from human dwellings. The level of transmission of Plasmodium to or from mosquitoes that have transitioned to outdoor feeding is not yet well-explored. In calculating the levels of infectious bites that correspond to given levels of deflected mosquitoes in the population, we therefore represent all the possible sources of reduced transmission by deflected vectors as a range of possible parameter values for their probability per feed of acquiring a Plasmodium infection, assuming probabilities of 4%, 2%, 1% and 0% (100%, 50%, 25% or 0% of the value assumed in the absence of deflection). From Figure 4, it can be seen that for a context in which transmission is as efficient for deflected outdoor-feeding vectors as for vectors exposed to no intervention, the initial reduction in infectious bites achieved when using an ESR (blue line with squares) is comparable with, but not quite equal to, the reduction in infectious bites achievable using an insecticide alone (pink line with crosses). However, the reduction in bites achieved using the ESR is maintained at a high level in the long term, while resistance rapidly eliminates the effectiveness of the unpaired insecticide. When assuming some reduction in transmission for deflected mosquitoes (lines with slash, circle and triangular markers), this trade-off between immediate and long-term benefits is reduced, with the ESR offering an initial reduction in infectious bites very similar to that achieved initially with insecticide alone, with the benefit again maintained or improved over the long term.

Figure 4

Download asset Open asset

Effect of evolved spatial repellent on infectious bites from a vector population over time.

The plots represent the infectious bites from the vector population per unit of time as a proportion of that with no intervention, assuming use of insecticide with and without ESR and probabilities per feed that ESR-deflected mosquitoes will become infected with *Plasmodium* of 4%, 2%, 1% or 0%. Plots otherwise use baseline parameter values. Probability per feed that non-deflected mosquito acquires *Plasmodium* infection is assumed to be 4%. If the size of the human population is assumed to be the same for all treatments and time periods, this equates to the entomological inoculation rate (EIR) as a percentage of the EIR with no intervention.

https://doi.org/10.7554/eLife.15416.009

Although the results summarized in Figures 1, 2 and 4 assume that resistant deflected phenotypes have the same fitness as susceptible deflected phenotypes, consistent with a context in which deflected phenotypes never enter a property treated with insecticide, and hence never experience any fitness benefit from resistance, note that this is not a requirement for successful establishment of an ESR, as illustrated in Figure 2—figure supplement 1.

Discussion

The use of DDT for indoor residual spray (IRS) programs proved an outstanding success in the history of public health campaigns against malaria. Its withdrawal in the light of environmental concerns reversed successes that approached elimination in some regions (Roberts, 2010; Curtis, 2002; Roberts et al., 1997). By the time of its withdrawal,resistance to the toxic effects of DDT was already widely observed, but there is empirical evidence to suggest that it may nonetheless have maintained efficacy as a transmission-reduction agent through its action as an effective spatial repellent (Roberts et al., 2000; Roberts and Alecrim, 1991). In part inspired by this, the deflection of malaria-vector mosquitoes from indoor feeding at night to outdoor feeding is being actively investigated as a means to reduce malaria transmission (Achee et al., 2012; University of Notre Dame, 2016). Here, we show that combining a spatial repellent which initially repels only a small proportion of a target vector population with indoor residual spraying of a high-toxicity insecticide can serve both to create a highly effective spatial repellent and to protect the companion insecticide from the rapid evolution of direct resistance to its toxic effects.

Our analysis provides the initial theoretical framework to spur empirical testing of this concept. There is, however, already a body of empirical evidence consistent with our proposal that selection can act to increase the efficacy of a repellent paired with a toxic insecticide. Our work predicts that, given heritable behavioural traits in the exposed population, toxic substances encountered with sufficient frequency and with sufficient volatility to be detectable at a distance would be observed to show repellent action against exposed populations. In fact, many substances investigated because of their repellent properties are also toxic on contact with target species. These include DEET, probably the most effective and well-known of the available repellents, as well as many of the naturally repellent plant-based volatiles currently in use or under investigation for personal protection (Xue et al., 2003; Licciardi et al., 2006). In fact, experiments that have combined repellents with insecticides on bed nets have often found that the benefits arise from the additive or synergistic toxicity of the repellents (Pennetier et al., 2007; Faulde et al., 2010; N'guessan et al., 2008) to an equivalent or greater degree than through their localized repellent effects. A comparison of the behavioural and insecticidal effects of three public-health insecticides provides more direct evidence of toxicity giving rise to repellency. The most toxic of the three, carbosulfan, also demonstrated high spatial repellence in a population of susceptible mosquitoes but not in a resistant population (Malima et al., 2009), an observation wholly consistent with our theoretical predictions.

The two elements defining the public-health benefits of an ESR used to deflect vectors away from human dwellings are: (1) the proportion of the mosquito population affected by the ESR; and (2) the reduction in infectious bites resulting from deflection. The former is dependent on the spread of the deflection allele in the population. This is determined by the relative fitness of susceptible, resistant, deflected and resistant-plus-deflected phenotypes, and by the initial prevalence of deflection and resistance alleles when the combined ESR-insecticide intervention is introduced. Choices about how the ESR and its associated insecticide are selected and deployed can influence or determine these values. The potential for establishing an ESR is maximized by lower fitness values for susceptibles and by higher fitness values for deflected phenotypes, indicating that control programs should target high coverage with high-efficacy insecticides to minimize the fitness of non-deflected susceptibles. Critically, all treated properties should ideally be treated with both an insecticide and ESR, rather than with either alone, in order to maximize selection for deflection and to minimize selection for physiological resistance. Minimising the entry of deflected phenotypes to insecticide-treated properties avoids exposing deflected phenotypes to insecticide and hence reducing the average fitness of deflected phenotypes relative to that of resistant phenotypes. Minimising the use of ESR on properties without insecticide coverage, particularly during the establishment stage of an ESR, avoids generating a fitness cost for deflected mosquitoes but not for non-deflected susceptible mosquitoes which can enter and feed. The use of a well-established ESR in a small number of properties for which insecticide use is for some reason impossible may, however, comprise one of ESR’s potential benefits. A high initial proportion of deflection and a low proportion of resistance alleles in the treated population also supports the spread of deflection alleles and the establishment of an ESR. Although these values cannot be directly controlled, desirable values may be targeted through careful choice of insecticide and ESR for each vector population. The candidate ESR can be chosen to target an initial threshold level of deflection genotype in the population. The initial companion insecticide can be chosen to target low initial levels of physiological resistance in the target population, and could exploit high-toxicity insecticides that would normally show a very rapid loss to resistance, serving both to maximize selection for deflected phenotypes and to give maximum immediate transmission reduction benefits.

An ESR needs to be partnered with a suitable IRS insecticide, which should have low contact repellence to maximize the mortality produced in non-deflected mosquitoes. Existing and new chemical insecticides with low contact repellence and high toxicity, such as bendiocarb (Evans, 1993), are therefore suitable as potential partner insecticides for an ESR program, potentially transforming them from short-term solutions to sustainable tools. Furthermore fungal biopesticides being developed for vector control potentially offer an ideal partnership with a spatial repellency treatment, whether DDT or a novel chemical, as they appear to have no contact repellency and have inherent resistance management benefits (Thomas and Read, 2007; Mnyone et al., 2010). These novel biopesticides have previously been proposed as ‘evolution proof’ late-life-acting insecticides (Read et al., 2009; Thomas and Read, 2007), offering respite from the treadmill of insecticide loss to resistance. For late-life-acting insecticides, relatively low-virulence fungal strains are ideal (Lynch et al., 2012), but much work on these organisms has focussed on generating high-virulence strains (Fang et al., 2012). Given their low or absent contact repellence, such high-virulence strains would provide ideal candidates for combination with an evolved spatial repellent, provided action can be taken before resistance becomes established in the target mosquito populations.

Our focus is on repellence that prevents indoor-feeding malaria vectors from entering properties, rather than on more localized repellence away from bed nets once vectors have entered properties. The excito-repellent properties of some of the pyrethroid insecticides currently used on long-lasting insecticide-treated nets (LLINs) reduce the mortality they generate by pushing vectors away before they acquire a lethal dose (Grieco et al., 2007; Tananchai et al., 2012). Such reduced mortality would diminish their likely efficacy as partner insecticides for ESRs. They could nonetheless enhance the establishment of an ESR when deployed in combination with suitable IRS, by reducing the relative fitness both of non-deflected susceptible mosquitoes and of mosquitoes resistant to the IRS insecticide. Such insecticides would also benefit from the resistance protection provided by an established ESR. As one of the key public-health tools currently deployed, and one already showing signs of succumbing to resistance (N'Guessan et al., 2007; Ochomo et al., 2014; Chandre et al., 1999), protecting the efficacy of these compounds could be hugely beneficial. Further, in the search for replacements for pyrethroids on bednets, alternative actives with low contact repellence could make excellent partner insecticides for an ESR and again benefit from protection against a rapid loss to resistance. If the partner insecticide is deployed via bednets then, as for IRS, the ESR can be deployed separately, allowing it to be applied, refreshed, removed or replaced without requiring any changes to the manufacture or maintenance of the LLINs.

The fitness cost of deflection away from indoor human hosts will vary according to vector species and degree of anthropophilly. Clearly for species that are readily zoophagous, such as Anopheles arabiensis, the fitness cost of diverting to a livestock host would be relatively low, provided that the alternative host is present and accessible. For more anthropophilic species, such as Anopheles gambiae, the fitness cost of deflection must be higher, and may be expected to vary between local populations depending on the feeding alternatives adopted. In assessing where ESRs offer high potential to contribute to public health campaigns, therefore, consideration should be given to the detail of the species mix in the local population, and the availability of alternative hosts. Immigration into the population will slow the spread of deflection, and so the size of a treated area is also a likely determinant of success since a large enough treatment area can minimize immigration by effectively including a whole breeding population within the treatment area. Equally, an isolated settlement may provide a closed vector population over a relatively small treatment area.

The separation of repellence and toxicity provides additional benefits. The companion insecticide can be changed whilst maintaining use of an established ESR. In contexts where an ESR can only be established on a transient basis, for instance because resistance alleles are already relatively common in the population when the ESR is introduced or because the population is subject to sustained immigration of non-deflected phenotypes, the ESR may still be established using a ‘ratchet’ approach in which the insecticide used in combination with the ESR can be changed when deflection phenotypes reach their peak, so that alleles for resistance to the new product will only offer a fitness benefit if paired with the relatively small proportion of non-deflection genotypes remaining in the population, allowing the deflection allele to spread. Where an ESR can become well established in populations in which the fitness of resistant phenotypes is higher than that of deflected phenotypes, unless the deflection allele reaches fixation, insecticide resistance will still spread eventually and deflection will eventually disappear. However, replacement of the partner insecticide at a suitable time will preserve the benefits of the established ESR, and the spread of resistance to the new partner insecticide may be wholly suppressed. This might open public-health opportunities such as the short-term use of a more expensive insecticide to generate ESR protection for the long-term use of a cheaper alternative.

When the conditions for establishment of an ESR are met, then differential transmission for outdoor/early biting vs indoor biting is potentially very important for the outcome in terms of reducing infectious bites. Because it is easiest to establish an ESR where deflection has little impact on fitness, if the only mechanism by which deflection to outdoor biting has an effect on infectious bites is through the incremental mortality associated with outdoor biting, then ESRs are most useful in the circumstances where they are hardest to establish. However, many other factors may reduce the probability per feed that a vector will acquire Plasmodium, or that, once infectious, it will transmit the parasite to a human host.

There is empirical evidence of reduced Plasmodium infection in at least one vector population which has transferred from indoor to outdoor biting (Ndiath et al., 2014). Ndiath et al. (Ndiath et al., 2014) found no sporozooites in a population of insecticide-susceptible An. gambiae pushed to outdoor feeding by the deployment of LLINs in Dielmo, whilst transmission continued at a high level in an indoor-feeding resistant population (Ndiath et al., 2014). Whilst this specific example cannot be unambiguously attributed to reduced transmission to/from outdoor feeding vectors rather than, for example, to differential mortality rates, it is reasonable to consider other potential causes of reduced outdoor transmission rates. These might include, for example, increased probabilities of taking feeds from non-human hosts (Lefèvre et al., 2009), differential availability of infectious human hosts indoors and outdoors, reduced effectiveness of transmission outside normal biting times (Gautret and Motard, 1999; Mideo et al., 2013; O'Donnell et al., 2011), and increased active host response to early/outdoor feeding attempts. Deflection to outdoor biting could therefore have dramatic effects on transmission that are independent of direct vector mortality effects, and there is potential for easily established ESRs which are as effective as conventional insecticides in reducing transmission, but on a much more sustainable basis.

There is still some controversy about the desirability of deflecting mosquitoes away from indoor biting and the well-established control methods which exploit this behaviour. The spatial repellent concept (whether evolved or conventional) is predicated on the idea that forcing vectors and Plasmodium into behavioural options that offer them lower fitness outcomes will inherently provide a new and potentially sustainable means of reducing transmission (Achee et al., 2012). However, outdoor biting is commonly viewed as a route to increased vector activity in a context where personal protection and anti-vector measures are hard to action (Bradley et al., 2012). In some contexts, for some compounds, the ESR concept may help to generate and sustain repellents that are able to protect against outdoor biting in the vicinity of treated properties, but in some locations it may nonetheless be the case that producing a move to outdoor biting simply does not generate overall transmission reductions equivalent to those achieved by control measures applied directly to indoor-feeding vectors. Transmission reductions may also change over time; for example, if outdoor biting initially generates low transmission because the accessible outdoor human hosts are all adults who have partial immunity, successful protection of children would eventually generate increasing numbers of susceptible adults who will in turn be exposed to outdoor feeding vectors, eliminating this benefit. It might also be argued that, whether beneficial or detrimental, outdoor biting will inevitably evolve as a behavioural resistance mechanism, and therefore that the use of an ESR is irrelevant. However, we would argue that ESRs offer benefits even in contexts where a transition to outdoor biting is ultimately found to be detrimental. For vectors that preferentially feed indoors, there is an expected fitness cost to a switch to outdoor feeding. Selection should therefore favour the use of a ‘cue’ to determine host choice, allowing the fitness benefits of indoor feeding to be enjoyed wherever safe, whilst avoiding the costs of entering insecticide-treated properties. Where ESRs are used in such contexts, either they will have no effect at all (i.e. they will do no harm nor good) or at least part of the vector population selected to respond to a ‘cue’ will be ‘cued’ by the ESR. This means that, should outdoor feeding prove undesirable at some point, the transition to outdoor feeding will be at least partly reversible since, unlike direct responses to an insecticide or to entering buildings, the response to an ESR will cease to affect vector behaviour if the ESR is withdrawn.

As well as being reversible in its effects, unlike many resistance-management strategies, ESR does not require the withdrawal or reduction of existing control measures in order to be effective. All novel interventions carry potential for harm, however, if they divert limited resources away from more effective control measures. In order to minimize potential harm as well as maximising potential benefits therefore, candidate ESRs should be cheap to obtain, and cheap and easy to deploy in the field. Unlike IRS, there is no imperative to cover interior surfaces with an ESR as they are intended to act at a distance, deployment as a single-spot event should therefore be possible. If this is formatted as a simple physical item, for example a small disc of impregnated paper, it can be deployed without the need for any special equipment. For ease of acceptance, candidate ESRs should have no detectable unpleasant odour for humans. To minimize unintended fitness costs for deflection, candidate substances should not be common in the natural environment of the target vector population. Any carrying material used for ESR distribution should be cheap, readily available, not amenable to other practical use, resistant to degradation and inedible for animals with indoor access.

The evolution of behavioural change in response to pesticides is also a recognized issue for agricultural insecticide use (Castillo-Chavez et al., 1988; Kennedy et al., 1987; Jongsma et al., 2010) and it is interesting to note that in the early 1980s Gould (1984) suggested that some insect populations might be expected to evolve behavioural avoidance of agricultural insecticides which were not originally repellent to them, and that in many situations, physiological resistance would evolve more slowly to insecticide formulations with high repellency than to non-repellent formulations. Whereas use in agriculture is commonly highly detrimental to the sustained utility of public-health tools, the parallel potential for ESRs in human health and agriculture may in some contexts allow agricultural use to enhance their public health role. In areas where vectors rest in insecticide-treated crops, using a vector-control spatial repellent on crops in combination with agricultural insecticides may serve to enhance selection for deflection alleles without generating any selection pressure for physiological resistance to public health insecticides (whilst incidentally providing some resistance-management for agricultural pesticides). This might also be a means to apply selection to nectar-feeding males, or to enhance selection for deflection in contexts where there is some existing resistance to all available public health insecticides but not to the agricultural insecticides being used locally.

Work carried out since the widescale use of DDT was terminated has clarified the role of spatial repellence in its contribution to malaria control (Roberts et al., 2000). Given its highly effective spatial repellence combined with high toxicity and the ubiquity of its use, it is perhaps interesting to consider whether DDT was always an effective repellent, or whether it actually provides the first empirical example of the ESR concept in action?

Our models suggest that there is clear potential to use evolution to create better control through evolved spatial repellence, and we hope to stimulate empirical work to test both our assumptions and the application of this novel approach. More generally, we demonstrate the importance of taking into account and modeling the evolutionary implications of different methods of insect control and medical interventions more broadly. We provide an example in which the inevitable evolution of the target insects can be used to improve rather than to reduce the effectiveness of the intervention. When evaluating a new intervention detailed mathematical modeling of the evolutionary outcomes alongside the more commonly considered epidemiological outcomes has considerable potential to improve infectious disease control. We recommend bringing together theoretical and empirical work to explore fully the potential of the ESR concept.

Appendix 1

Definition of the PM model

Appendix 1—table 1

Variable and parameter definitions for the population genetics model.

https://doi.org/10.7554/eLife.15416.010

Variable and parameter labels
Label	Description	Units
$t$	Index of model time periods length of one time period equivalent to duration of one gonotrophic cycle	Time
$h$	Index of phenotypes 1= susceptible non-deflected 2 = resistant non-deflected 3 = resistant and deflected 4 = susceptible deflected
$g$	Index of genotypes 1=rrdd, 2=rRdd, 3=RRdd, 4=rrDd, 5=rRDd, 6=RRDd, 7=rrDD, 8=rRDD, 9=RRDD
$D_{w, t}$	The proportion of male gametes with allele pair $w$ in time period $t$ . Pairs are: 1=rd, 2=rD, 3=Rd, 4=RD	Proportion of male gametes
$T_{t}$	Infectious bites from population in period $t$ as a proportion of baseline (pre-treatment) infectious bites per period	Proportion of baseline bites/time
$I_{t}$	Females giving infectious bites in period $t$ as a proportion of females in baseline population	Proportion of baseline population
$M_{h, t}$	Females of phenotype $h$ which are infectious with Plasmodium in period $t$ as a proportion of baseline population	Proportion of baseline population
$c_{0}$	Baseline probability that an infectious female in the no-intervention baseline population will give an infectious bite when feeding	Probability
$c_{h}$	Probability that an infectious female of phenotype $h$ will give an infectious bite when feeding	Probability
$m_{h, t}$	Females newly infected with Plasmodium in period $t$ as a proportion of the baseline population	Proportion of baseline population
$β_{0}$	Baseline average probability susceptible individuals survive one model time period (in the absence of any intervention)	Probability
$β_{h}$	Average probability that individuals with phenotype $h$ will survive one time period in an environment including a given intervention	Probability
$γ$	Number of model time periods between start of Plasmodium infection in mosquito, and point from which infection has matured and mosquito can give infectious bites	Model time periods
$b_{0}$	Probability that an uninfected mosquito in a baseline population subject to no interventions will acquire a Plasmodium infection during a single model time-period	Probability
$b_{h}$	Probability that an uninfected mosquito with phenotype $h$ will acquire a Plasmodium infection during a single time-period	Probability
$U_{h, t}$	Females with phenotype $h$ which are not infected with Plasmodium at start of time period $t$ , as a proportion of the baseline population	Proportion of baseline population
$J_{h, t}$	New females of phenotype $h$ in period $t$ , as a proportion of the baseline population	Proportion of baseline population
$B_{g, t}$	Proportion of new adults in time period $t$ with genotype $g$	Proportion of new adults
$τ_{1, t}$	New (female) adults in period $t$ as proportion of baseline population	Proportion of baseline population
$τ_{2, t}$	Adults surviving from preceding period to start of period $t$ as proportion of baseline population	Proportion of baseline ppopulation
$τ_{3, t}$	New (female) adults in period $t$ as a proportion of total (female) population in period $t$	Proportion of current population
$τ_{4, t}$	Adults surviving from preceding period to start of period $t$ as proportion of population in period $t$	Proportion of current population
$τ_{5, t}$	Population in period $t$ as proportion of baseline population	Proportion of baseline population
$S_{g}$	Probability a female with genotype $g$ will survive one model time period	Probability
$E_{g, t}$	Proportion of all eggs laid in time period $t$ which have genotype $g$	Proportion of eggs laid in period
$K_{g, t}$	Proportion of adult females surviving at start of time period $t$ which have genotype $g$	Proportion of adult females
$ϖ$	Number of time periods between egg-laying and emergence of new adults	Model time periods
$W_{g, t}$	The proportion of all females (including new adult females) which have genotype $g$ at start of time period $t$	Proportion of adult females
$ψ_{t}$	Proportion of the population alive at start of time period $t - 1$ which survive to start of time period $t$	Proportion of population
$θ_{g, t}$	Number of eggs with genotype $g$ laid in period $t$ , relative to a normalized number of eggs assumed per mosquito per successful lay	Relative number of eggs
$δ_{t}$	Total number of eggs laid in period $t$ , relative to a normalized number of eggs assumed per mosquito per successful lay	Relative number of eggs
$Z_{g, j, t}$	The proportion of eggs produced by females which are new adults in time period $t$ and have genotype $j$ , which are genotype $g$	Proportion of eggs
$A_{g, t}$	Average number of eggs laid during time period $t$ by females with genotype $g$ alive at start of time period $t$	Normalized number of eggs
$L$	Normalized number of eggs assumed per mosquito per successful lay	Eggs/mosquito/lay
$X_{t}$	Newly mated females in time period $t$ as a proportion of the baseline population	Proportion of baseline population
$G_{1, t}$	Proportion of alleles at the resistance locus which are resistance alleles in model time period $t$	Proportion of alleles
$G_{2, t}$	Proportion of alleles at the deflection locus which are deflection alleles in model time period $t$	Proportion of alleles
$P_{1, t}$	Proportion of population with resistant phenotypes in model time period $t$	Proportion of current population
$P_{2, t}$	Proportion of population with deflected phenotypes in model time period $t$	Proportion of current population
$P_{3, t}$	Proportion of population with susceptible non-deflected phenotypes in model time period $t$	Proportion of current population
$V_{y, t}$	Proportion of population with genotype $y$ in model time period $t$ genotypes; 1=rr, 2=rR or RR, 3=dd, 4=dD or DD	Proportion of current population

The model tracks the population of adult females and assumes a ‘baseline population’ of stable size (new adult recruitment rate is constant and equal to adult mortality rate) in the absence of any insecticide or ESR interventions. It is assumed that larval habitats are over-subscribed, and variations in the adult population do not materially change the number of juveniles reaching maturity and joining the adult population in each time period. Females are assumed to mate once, as newly emerged adults, mating with males having the same genotype proportions as newly emerged adult females.

$t$ represents the index of model time periods, with $t = 0$ representing the periods prior to application of an intervention when the population is the baseline population.

$S_{g}$ , the probability that a female with genotype $g$ will survive one time period, is the primary input value driving the spread of different alleles in the popgen model. Each genotype is associated with one of the four possible phenotypes. Each phenotype has an associated per period survival probability.

Infectious bites from population in period $t$ as a proportion of baseline infectious bites, $T_{t} = \frac{I_{t}}{I_{0}}$ .

Females giving infectious bites in period t as a proportion of females in the baseline population, $I_{t} = \sum_{h = 1}^{4} M_{h, t} c_{h}$ .

The females of genotype $h$ that are newly infected with Plasmodium in period $t$ , as a proportion of the baseline population, $m_{h, t}$ , is calculated as $m_{h, 0} = U_{h, 0} β_{0} b_{0}$ and $m_{h, t} = U_{h, t} β_{h} b_{h}$ with $t > 0$ .

$M_{h, t}$ females of phenotype $h$ that are infectious with Plasmodium in period $t$ as a proportion of baseline population, calculated for the baseline population as the total proportion of newly infected females in each period, for all periods up to $γ$ periods (number of periods between Plasmodium infection and infectiousness) ago, multiplied by the probability of surviving to $γ$ periods ago, multiplied by the probability of surviving $γ$ periods, $M_{h, 0} = \frac{m_{h, 0}}{1 - β_{0}} {β_{0}}^{γ}$ .

After baseline period, $M_{h, t}$ is calculated for each phenotype as surviving infectious females from the previous period, plus survivors from newly infected females with that phenotype from $γ$ periods ago. $M_{h, t} = M_{h, (t - 1)} β_{h} + m_{h, (t - m i n (t, γ))} {β_{h}}^{γ} t > 0$

Females with phenotype $h$ that are not infected with Plasmodium at start of time period $t$ , as a proportion of the baseline population.

For the baseline population, $U_{h, 0}$ , the non-infected females at the start of each time period are calculated for each phenotype as the sum of the surviving, uninfected new females from all past time periods, assuming a constant value for new females of each phenotype in each period, $U_{h, 0} = \frac{J_{h, 0}}{1 - β_{0} (1 - b_{0})}$ . Post-introduction of the intervention, the non-infected females for each genotype at the start of each time period are calculated as the surviving non-infected females from the previous period plus new females of that genotype, giving $U_{h, t} = U_{h, t - 1} β_{h} (1 - b_{h}) + J_{h, t}$ with $t > 0$ .

$J_{h, t}$ , new females of phenotype $h$ in period $t$ , expressed as a proportion of the baseline population, is calculated as the proportion of the new adults in period $t$ that have genotypes corresponding to phenotype $h$ .

New susceptible females in time period t as a proportion of the baseline population,

J_{1, t} = B_{1, t} τ_{1, t} .

New deflected females in time period t as a proportion of the baseline population,

J_{2, t} = (B_{4, t} + B_{7, t}) τ_{1, t} .

New resistant females in time period t as a proportion of the baseline population,

J_{3, t} = (B_{2, t} + B_{3, t}) τ_{1, t} .

New deflected-resistant females in time period t as a proportion of the baseline population,

J_{4, t} = (B_{5, t} + B_{6, t} + B_{8, t} + B_{9, t}) τ_{1, t} .

For this analysis, it is assumed that larval habitats are over-subscribed and that the number of new adults recruited to the population in each time period remains constant with and without interventions. In all time periods therefore, the number of new (female) adults joining the population is the same as for the baseline population. Since the baseline population is assumed to be constant, the proportion of the population comprising new adults in each time period in the baseline population can be calculated as $τ_{1, t} = 1 - β_{0}$ .

$τ_{2, t}$ , adult females surviving from the preceding period to the start of period $t$ as a proportion of the baseline population, is calculated as $τ_{2, 0} = 1 - τ_{1, 0}$ and $τ_{2, t} = \sum_{g = 1}^{9} (τ_{1, t - 1} B_{g, t - 1} + τ_{2, t - 1} K_{g, t - 1}) S_{g}$ with $t > 0$ .

New (female) adults in period $t$ as a proportion of the total of population in period $t$ , $τ_{3, t} = \frac{τ_{1, t}}{τ_{1, t} + τ_{2, t}} .$

Adult females surviving from the preceding period to the start of period $t$ as a proportion of the population in period $t$ , $τ_{4, t} = \frac{τ_{2, t}}{τ_{1, t} + τ_{2, t}} .$

Population of adult females in period $t$ as a proportion of the baseline population, $τ_{5, t} = τ_{1, t} + τ_{2, t}$ .

The proportion of new adults in each cycle that have genotype $g$ in the baseline population, $B_{g, 0}$ is calculated based on specified initial allele proportions, assuming Hardy-Weinberg equilibrium (independently for each locus).

The proportion of new adults in time period $t$ (with $t > 0$ ) that have genotype $g$ is calculated as the proportion of eggs with that genotype from $ϖ$ periods ago, where $ϖ$ represents the number of model time periods required for development from egg to new adults. For early time periods, before $t > ϖ$ , the baseline egg genotype proportions are used ( $t = 0$ ), giving $B_{g, t} = E_{g, t - m i n (t, ϖ)}$ .

The proportion of surviving adult females at start of time period $t$ that have genotype $g$ , $K_{g, t}$ , for the baseline population is the same as the proportion of eggs with that genotype, giving $K_{g, 0} = B_{g, 0}$ , since in the absence of an intervention there is assumed to be no difference in survival or fecundity between the phenotypes associated with each genotype. Note that for simplicity, this assumption is applied to the baseline population even when a cost of resistance is considered once an intervention is introduced. Following introduction of an intervention, the proportion of surviving adult females at start of time period $t$ that have genotype $g$ is calculated as the proportion of females with genotype $g$ that survive from the previous period, divided by the total proportion surviving, with all genotypes, giving $K_{g, t} = \frac{W_{g, t - 1} S_{g}}{ψ_{t}} t > 0$ .

$W_{g, t}$ , the proportion of females with genotype $g$ at start of time period $t$ is calculated as the sum of new adults and surviving adults with that genotype at that time, $W_{g, t} = B_{g, t} τ_{3, t} + K_{g, t} τ_{4, t} .$

$ψ_{t}$ , the proportion of the population alive at start of time period $t - 1$ that survive to start of time period $t$ is calculated as the sum of the proportions with each genotype alive at start of time period $t - 1$ multiplied by the one-period survival probability associated with that genotype $ψ_{t} = \sum_{i = 1}^{9} (W_{i, t - 1} S_{i})$ .

The model assumes that females mate once, as new adults, with newly emerged adult males. The normalized number of eggs with genotype $g$ laid in period $t$ , $θ_{g, t}$ , is therefore calculated as the total for all previous periods of the proportion of eggs of genotype $g$ resulting from mating with new males in that period, for females of each genotype, multiplied by the probability for females of that genotype of surviving to period $t$ : $t θ_{g, t} = \sum_{i = 1}^{t} ⟮ \sum_{j = 1}^{9} (B_{j, i} Z_{g, j, i} A_{j, i} {S_{j}}^{t - i}) τ_{1, i} ⟯$ .

The proportion of the eggs laid in time period $t$ that have genotype $g$ is calculated as $E_{g, t} = \frac{θ_{g, t}}{_{δ_{t}}} .$

The total (normalized) number of eggs laid in period $t$ is calculated as $δ_{t} = \sum_{g = 1}^{9} θ_{g . t} .$

Some combinations of maternal and offspring genotypes are not possible, so the proportion of eggs of genotype $g$ produced by females with genotype $j$ that are new adults in time period $t$ , $Z_{g, j, t}$ ,is zero for all time periods for some genotypes.

\begin{array}{ll} Z_{1, i, t} = 0 & {i = 3, 6, 7, 8, 9} \\ Z_{2, i, t} = 0 & {i = 7, 8, 9} \\ Z_{3, i, t} = 0 & {i = 1, 7, 8, 9} \\ Z_{4, i, t} = 0 & {i = 3, 6, 9} \\ Z_{6, i, t} = 0 & {i = 1, 4, 7} \\ Z_{7, i, t} = 0 & {i = 1, 2, 3, 7, 9} \\ Z_{8, i, t} = 0 & {i = 1, 2, 3} \\ Z_{9, i, t} = 0 & {i = 1, 2, 3, 4, 7} \end{array}

For other combinations values are calculated as follows:

$Z_{1, 1, t} = D_{1, t}$	% eggs which are genotype rrdd in time period $t$ for new females with genotype rrdd
$Z_{1, 2, t} = 0.5 D_{1, t}$	% eggs which are genotype rrdd in time period $t$ for new females with genotype rRdd
$Z_{1, 4, t} = 0.5 D_{1, t}$	% eggs which are genotype rrdd in time period $t$ for new females with genotype rrdD
$Z_{1, 5, t} = 0.25 D_{1, t}$	% eggs which are genotype rrdd in time period $t$ for new females with genotype rRdD
$Z_{2, 1, t} = D_{3, t}$	% eggs which are genotype rRdd in time period $t$ for new females with genotype rrdd
$Z_{2, 2, t} = 0.5 (D_{1, t} + D_{3, t})$	% eggs which are genotype rRdd in time period $t$ for new females with genotype rRdd
$Z_{2, 3, t} = D_{1, t}$	% eggs which are genotype rRdd in time period $t$ for new females with genotype RRdd
$Z_{2, 4, t} = 0.5 D_{3, t}$	% eggs which are genotype rRdd in time period $t$ for new females with genotype rrdD
$Z_{2, 5, t} = 0.25 (D_{1, t} + D_{3, t})$	% eggs which are genotype rRdd in time period $t$ for new females with genotype rRdD
$Z_{2, 6, t} = 0.5 D_{1, t}$	% eggs which are genotype rRdd in time period $t$ for new females with genotype RRdD
$Z_{3, 2, t} = 0.5 D_{3, t}$	% eggs which are genotype RRdd in time period $t$ for new females with genotype rRdd
$Z_{3, 3, t} = D_{3, t}$	% eggs which are genotype RRdd in time period $t$ for new females with genotype RRdd
$Z_{3, 5, t} = 0.25 D_{3, t}$	% eggs which are genotype RRdd in time period $t$ for new females with genotype rRdD
$Z_{3, 6, t} = 0.5 D_{3, t}$	% eggs which are genotype RRdd in time period $t$ for new females with genotype RRdD
$Z_{4, 1, t} = D_{2, t}$	% eggs which are genotype rrdD in time period $t$ for new females with genotype rrdd
$Z_{4, 2, t} = 0.5 D_{2, t}$	% eggs which are genotype rrdD in time period $t$ for new females with genotype rRdd
$Z_{4, 4, t} = 0.5 (D_{1, t} + D_{2, t})$	% eggs which are genotype rrdD in time period $t$ for new females with genotype rrdD
$Z_{4, 5, t} = 0.25 (D_{1, t} + D_{2, t})$	% eggs which are genotype rrdD in time period $t$ for new females with genotype rRdD
$Z_{4, 7, t} = D_{1, t}$	% eggs which are genotype rrdD in time period $t$ for new females with genotype rrDD
$Z_{4, 8, t} = 0.5 D_{1, t}$	% eggs which are genotype rrdD in time period $t$ for new females with genotype rRDD
$Z_{5, 1, t} = D_{4, t}$	% eggs which are genotype rRdD in time period $t$ for new females with genotype rrdd
$Z_{5, 2, t} = 0.5 (D_{2, t} + D_{4, t})$	% eggs which are genotype rRdD in time period $t$ for new females with genotype rRdd
$Z_{5, 3, t} = D_{2, t}$	% eggs which are genotype rRdD in time period $t$ for new females with genotype RRdd
$Z_{5, 4, t} = 0.5 (D_{3, t} + D_{4, t})$	% eggs which are genotype rRdD in time period $t$ for new females with genotype rrdD
$Z_{5, 5, t} = 0.25 (D_{1, t} + D_{2, t} + D_{3, t} + D_{4, t})$	% eggs which are genotype rRdD in time period $t$ for new females with genotype rRdD
$Z_{5, 6, t} = 0.5 (D_{1, t} + D_{2, t})$	% eggs which are genotype rRdD in time period $t$ for new females with genotype RRdD
$Z_{5, 7, t} = D_{3, t}$	% eggs which are genotype rRdD in time period $t$ for new females with genotype rrDD
$Z_{5, 8, t} = 0.5 (D_{1, t} + D_{3, t})$	% eggs which are genotype rRdD in time period $t$ for new females with genotype rRDD
$Z_{5, 9, t} = D_{1, t}$	% eggs which are genotype rRdD in time period $t$ for new females with genotype RRDD
$Z_{6, 2, t} = 0.5 D_{4, t}$	% eggs which are genotype RRdD in time period $t$ for new females with genotype rRdd
$Z_{6, 3, t} = D_{4, t}$	% eggs which are genotype RRdD in time period $t$ for new females with genotype RRdd
$Z_{6, 5, t} = 0.25 (D_{3, t} + D_{4, t})$	% eggs which are genotype RRdD in time period $t$ for new females with genotype rRdD
$Z_{6, 6, t} = 0.5 (D_{3, t} + D_{4, t})$	% eggs which are genotype RRdD in time period $t$ for new females with genotype RRdD
$Z_{6, 8, t} = 0.5 D_{3, t}$	% eggs which are genotype RRdD in time period $t$ for new females with genotype rRDD
$Z_{6, 9, t} = D_{3, t}$	% eggs which are genotype RRdD in time period $t$ for new females with genotype RRDD
$Z_{7, 4, t} = 0.5 D_{2, t}$	% eggs which are genotype rrDD in time period $t$ for new females with genotype rrdD
$Z_{7, 5, t} = 0.25 D_{2, t}$	% eggs which are genotype rrDD in time period $t$ for new females with genotype rRdD
$Z_{7, 7, t} = D_{2, t}$	% eggs which are genotype rrDD in time period $t$ for new females with genotype rrDD
$Z_{7, 8, t} = 0.5 D_{2, t}$	% eggs which are genotype rrDD in time period $t$ for new females with genotype rRDD
$Z_{8, 4, t} = 0.5 D_{4, t}$	% eggs which are genotype rRDD in time period $t$ for new females with genotype rrDD
$Z_{8, 5, t} = 0.25 (D_{2, t} + D_{4, t})$	% eggs which are genotype rRDD in time period $t$ for new females with genotype rRdD
$Z_{8, 6, t} = 0.5 D_{2, t}$	% eggs which are genotype rRDD in time period $t$ for new females with genotype RRdD
$Z_{8, 7, t} = D_{4, t}$	% eggs which are genotype rRDD in time period $t$ for new females with genotype rrDD
$Z_{8, 8, t} = 0.5 (D_{2, t} + D_{4, t})$	% eggs which are genotype rRDD in time period $t$ for new females with genotype rRDD
$Z_{8, 9, t} = D_{2, t}$	% eggs which are genotype rRDD in time period $t$ for new females with genotype RRDD
$Z_{9, 5, t} = 0.25 D_{4, t}$	% eggs which are genotype RRDD in time period $t$ for new females with genotype rRdD
$Z_{9, 6, t} = 0.5 D_{4, t}$	% eggs which are genotype RRDD in time period $t$ for new females with genotype RRdD
$Z_{9, 8, t} = 0.5 D_{4, t}$	% eggs which are genotype RRDD in time period $t$ for new females with genotype rRDD
$Z_{9, 9, t} = D_{4, t}$	% eggs which are genotype RRDD in time period $t$ for new females with genotype RRDD

The average number of eggs laid during time period $t$ by females with genotype $g$ that are alive at the start of time period $t$ , $A_{g, t}$ , is calculated as the normalized number of eggs per mosquito per lay, multiplied by the probability that females of genotype $g$ will survive one time period. For the baseline population, this gives $A_{g, 0} = β_{0} L$ , and for the population post introduction of an intervention, $A_{g, t} = S_{g} L t > 0$ .

The model is interested in comparative and proportionate rather than absolute values for egg production, for convenience therefore, we use a normalized value for the number of eggs assumed per mosquito per successful lay, represented by $L$ .

$X_{t}$ represents the newly mated females in time period $t$ as % baseline population. For the baseline population, the genotype proportions in females and their eggs is assumed to be constant over time. Females surviving from the baseline population will therefore contribute eggs to subsequent periods with the same genotype proportions in the eggs for each maternal genotype, irrespective of age. For this purpose, therefore, as all females will have the same characteristics as newly mated females, we have $X_{0} = 1$ . Thereafter $X_{t} = τ_{1, t}$ , with $t > 0$ .

Mating males in time period $t$ are assumed to have the genotype proportions of newly emerged adults in period $t$ . With allele pairs indexed as 1=rd, 2=rD, 3=Rd, 4=RD, the proportion of male gametes with allele pair $w$ in time period $t$ , $D_{w, t}$ , is therefore calculated as:

\begin{array}{ll} D_{1, t} = & B_{1, t} + 0.5 (B_{2, t} + B_{4, t}) + 0.25 B_{5, t} \\ D_{2, t} = & B_{7, t} + 0.5 (B_{4, t} + B_{8, t}) + 0.25 B_{5, t} \\ D_{3, t} = & B_{3, t} + 0.5 (B_{2, t} + B_{6, t}) + 0.25 B_{5, t} \\ D_{4, t} = & B_{9, t} + 0.5 (B_{8, t} + B_{6, t}) + 0.25 B_{5, t} \end{array}

The proportion of alleles at the resistance locus in the population at time period $t$ that are resistance alleles, $G_{1, t}$ , is calculated as $G_{1, t} = W_{3, t} + W_{6, t} + W_{9, t} + 0.5 (W_{2, t} + W_{5, t} + W_{8, t})$ .

The proportion of alleles at the deflection locus in the population at time period $t$ that are deflection alleles, $G_{2, t}$ , is calculated as $G_{2, t} = 0.5 (W_{4, t} + W_{5, t} + W_{6, t}) + W_{7, t} + W_{8, t} + W_{9, t}$ .

The proportion of the adult female population in time period $t$ that have resistant (non-deflected) phenotypes, $P_{1, t}$ , is the sum of those with genotypes rRdd or RRdd, $P_{1, t} = W_{2, t} + W_{3, t}$ .

The proportion of the adult female population in time period $t$ that have deflected phenotypes, $P_{2, t}$ , is the sum of those with genotypes having Dd or DD in combination with any allele combination at the resistance locus, $P_{2, t} = \sum_{i = 4}^{9} W_{i, t}$ .

The proportion of the adult female population in time period $t$ that have susceptible, non-deflected phenotypes, $P_{3, t}$ , is equal to the proportion with rrdd genotypes: $P_{3, t} = 1 - P_{1, t} - P_{2, t}$ .

The proportion of the adult female population in time period $t$ that have two susceptible alleles (rr) at the resistance locus, $V_{1, t}$ , is calculated as $V_{1, t} = W_{1, t} + W_{4, t} + W_{7, t}$ .

The proportion of the adult female population in time period $t$ that have at least one resistance allele (rR, RR) at the resistance locus, $V_{2, t}$ , is calculated as $V_{2, t} = W_{2, t} + W_{3, t} + W_{5, t} + W_{6, t} + W_{8, t} + W_{9, t}$ .

The proportion of the adult female population in time period $t$ that have two non-deflection alleles (dd) at the deflection locus, $V_{3, t}$ , is calculated as $V_{3, t} = W_{1, t} + W_{2, t} + W_{3, t}$ .

The proportion of the adult female population in time period $t$ that have at least one deflection allele (rD, DD) at the resistance locus, $V_{4, t}$ , is calculated as $V_{4, t} = \sum_{i = 4}^{9} W_{i, t}$ .

Appendix 2

Derivation of Equations 1 and 2 and Expressions 1 and 2

Because the fitness of the individuals carrying deflection and non-deflection alleles depends on whether they are also carrying resistance alleles, and because resistance alleles may originate from the same parent as the relevant deflection allele or from the paired zygote, we consider the average fitness of the phenotypes into which deflection and non-deflection alleles will be inherited, based on the proportions of each possible allele pair in the population of zygotes, assuming the same proportions in gametes of mating males and newly emerged females.

Consider the allele pairs in gametes in the mating adult population:

\begin{array}{ll} [d r] & = p r o p o r t i o n o f d e f l e c t i o n l o c u s a l l e l e s w h i c h a r e n o n - d e f l e c t i o n a l l e l e s p a i r e d w i t h \\ n o n - r e s i s t a n c e a l l e l e s i n g a m e t e s \\ [d R] & = p r o p o r t i o n o f d e f l e c t i o n l o c u s a l l e l e s w h i c h a r e n o n - d e f l e c t i o n a l l e l e s p a i r e d w i t h \\ r e s i s t a n c e a l l e l e s i n g a m e t e s \\ [D r] & = p r o p o r t i o n o f d e f l e c t i o n l o c u s a l l e l e s w h i c h a r e d e f l e c t i o n a l l e l e s p a i r e d w i t h \\ n o n - r e s i s t a n c e a l l e l e s i n g a m e t e s \\ [D R] & = p r o p o r t i o n o f d e f l e c t i o n l o c u s a l l e l e s w h i c h a r e d e f l e c t i o n a l l e l e s p a i r e d w i t h \\ r e s i s t a n c e a l l e l e s i n g a m e t e s \\ [d -] & = p r o p o r t i o n o f d e f l e c t i o n l o c u s a l l e l e s i n g a m e t e s w h i c h a r e n o n - d e f l e c t i o n a l l e l e s \\ [D -] & = p r o p o r t i o n o f d e f l e c t i o n l o c u s a l l e l e s i n g a m e t e s w h i c h a r e d e f l e c t i o n a l l e l e s \\ [- r] & = p r o p o r t i o n o f d e f l e c t i o n l o c u s a l l e l e s w h i c h a r e p a i r e d w i t h s u s c e p t i b l e a l l e l e s i n g a m e t e s \\ [- R] & = p r o p o r t i o n o f d e f l e c t i o n l o c u s a l l e l e s w h i c h a r e p a i r e d w i t h r e s i s t a n c e a l l e l e s i n g a m e t e s \\ [d -] + [D -] = [- r] + [- R] = 1 \end{array}

Let ${[S]}_{d}$ represent the proportion of non-deflection alleles inherited into offspring with genotypes corresponding to susceptible, non-deflected phenotypes:

{[S]}_{d} = \frac{{[d r]}^{2}}{[d -]} .

Let ${[R]}_{d}$ represent the proportion of non-deflection alleles inherited into offspring with genotypes corresponding to resistant, non-deflected phenotypes:

{[R]}_{d} = \frac{[d r] [d R] + [d R] [d r] + [d R] [d R]}{[d -]} = \frac{[d -] [d R] + [d R] [d r]}{[d -]} = [d R] + \frac{[d R] [d r]}{[d -]} .

Let ${[D R]}_{d}$ represent the proportion of non-deflection alleles inherited into offspring with genotypes corresponding to resistant, deflected phenotypes:

{[D R]}_{d} = \frac{[d r] [D R] + [d R] [D r] + [d R] [D R]}{[d -]} = \frac{[d -] [D R] + [d R] [D r]}{[d -]} = [D R] + \frac{[d R] [D r]}{[d -]} .

Let ${[D]}_{d}$ represent the proportion of non-deflection alleles inherited into offspring with genotypes corresponding to non-resistant, deflected phenotypes:

{[D]}_{d} = \frac{[d r] [D r] \}{[d -]} .

Let $F_{S}$ represent the fitness of susceptible, non-deflected phenotypes.

Let $F_{R}$ represent the fitness of resistant, non-deflected phenotypes.

Let $F_{R D}$ represent the fitness of resistant, deflected phenotypes.

Let $F_{D}$ represent the fitness of non-resistant, deflected phenotypes.

Let ${\bar{F}}_{D}$ represent the average fitness of offspring into which deflection alleles are inherited.

Let ${[D R]}_{D}$ represent the proportion of deflection alleles inherited into offspring with genotypes corresponding to resistant, deflected phenotypes:

\begin{array}{ll} {[D R]}_{D} & = \frac{[D r] [d R] + [D r] [D R] + [D R] [D R] + [D R] [D r] + [D R] [d R] + [D R] [d r]}{[D -]} \\ = [- R] + \frac{[D R] [- r]}{[D -]} \end{array}

Let ${[D]}_{D}$ represent the proportion of deflection alleles inherited into offspring with genotypes corresponding to non-resistant, deflected phenotypes:

\begin{array}{ll} {[D]}_{D} & = \frac{[D r] [d r] + [D r] [D r]}{[D -]} \\ = \frac{[D r] [- r]}{[D -]} \\ {\bar{F}}_{D} & = {[D R]}_{D} F_{R D} + {[D]}_{D} F_{D} \\ = ([- R] + \frac{[D R] [- r]}{[D -]}) F_{R D} + \frac{[D r] [- r]}{[D -]} F_{D} \\ = F_{S} + \frac{[D r] [- R] + [D R]}{[D -]} (F_{R D} - F_{S}) + \frac{[D r] [- r]}{[D -]} (F_{D} - F_{S}) \end{array}

This is Equation 1.

Let ${\bar{F}}_{d}$ represent the average fitness of offspring into which non-deflection alleles are inherited:

\begin{array}{ll} {\bar{F}}_{d} & = [S]_{d} F_{S} + [R]_{d} F_{R} + {[R D]}_{d} F_{R D} + {[D]}_{d} F_{D} \\ {\bar{F}}_{d} & = \frac{{[d r]}^{2}}{[d -]} F_{S} + ([d R] + \frac{[d R] [d r]}{[d -]}) F_{R} + ([D R] + \frac{[d R] [D r]}{[d -]}) F_{R D} + \frac{[d r] [D r]}{[d -]} F_{D} \\ {\bar{F}}_{d} & = F_{S} + ([d R] + \frac{[d R] [d r]}{[d -]}) (F_{R} - F_{S}) + \frac{[d r] [D r]}{[d -]} (F_{D} - F_{S}) + ([D R] + \frac{[d R] [D r]}{[d -]}) (F_{R D} - F_{S}) \end{array}

This is Equation 2.

Deflection alleles will spread when the average fitness of offspring into which deflection alleles are inherited is higher than that of offspring into which non-deflection alleles are inherited:

\begin{array}{ll} {\bar{F}}_{D} & > {\bar{F}}_{d} \leftrightarrow F_{S} + \frac{[D r] [- R] + [D R]}{[D -]} (F_{R D} - F_{S}) + \frac{[D r] [- r]}{[D -]} (F_{D} - F_{S}) \\ > F_{S} + ([d R] + \frac{[d R] [d r]}{[d -]}) (F_{R} - F_{S}) + ([D R] + \frac{[d R] [D r]}{[d -]}) (F_{R D} - F_{S}) + \frac{[d r] [D r]}{[d -]} (F_{D} - F_{S}) \end{array}

This can be rearranged as follows to give Expression 1.

\begin{array}{ll} {\bar{F}}_{D} & > {\bar{F}}_{d} \leftrightarrow \frac{[D r] [- R] + [D R]}{[D -]} (F_{R D} - F_{S}) + \frac{[D r] [- r]}{[D -]} (F_{D} - F_{S}) \\ > ([d R] + \frac{[d R] [d r]}{[d -]}) (F_{R} - F_{S}) + ([D R] + \frac{[d R] [D r]}{[d -]}) (F_{R D} - F_{S}) + \frac{[d r] [D r]}{[d -]} (F_{D} - F_{S}) \\ {\bar{F}}_{D} & > {\bar{F}}_{d} \leftrightarrow \frac{[D r] [- R] + [D R]}{[D -]} (F_{R D} - F_{D}) + \frac{[D r] [- r] + [D r] [- R] + [D R]}{[D -]} (F_{D} - F_{S}) \\ > ([d R] + \frac{[d R] [d r]}{[d -]}) (F_{R} - F_{S}) + ([D R] + \frac{[d R] [D r]}{[d -]}) (F_{R D} - F_{D}) \\ + (\frac{[d r] [D r]}{[d -]} + [D R] + \frac{[d R] [D r]}{[d -]}) (F_{D} - F_{S}) \\ {\bar{F}}_{D} & > {\bar{F}}_{d} \leftrightarrow \frac{[D r] [- R] + [D R]}{[D -]} (F_{R D} - F_{D}) + (F_{D} - F_{S}) \\ > ([d R] + \frac{[d R] [d r]}{[d -]}) (F_{R} - F_{S}) + [D -] (F_{D} - F_{S}) + ([D R] + \frac{[D r] [d R]}{[d -]}) (F_{R D} - F_{D}) \\ {\bar{F}}_{D} & > {\bar{F}}_{d} \leftrightarrow (F_{D} - F_{S}) \\ > ([d R] + \frac{[d R] [d r]}{[d -]}) (F_{R} - F_{S}) + [D -] (F_{D} - F_{S}) + ([D R] + \frac{[D r] [d R]}{[d -]} - \frac{[D r] [- R] + [D R]}{[D -]}) (F_{R D} - F_{D}) \end{array}

\begin{aligned} {\bar{F}}_{D} & > {\bar{F}}_{d} \leftrightarrow (F_{D} - F_{S}) > \\ [D -] (F_{D} - F_{S}) + [d R] (1 + \frac{[d r]}{[d -]}) (F_{R} - F_{S}) \\ + ([D R] (1 - \frac{1}{[D -]}) + [D r] (\frac{[d R]}{[d -]} - \frac{[- R]}{[D -]})) (F_{R D} - F_{D}) \\ {\bar{F}}_{D} & > {\bar{F}}_{d} \leftrightarrow (F_{D} - F_{S}) > ([D -] + [d R] + \frac{[d R] [d r]}{[d -]}) (F_{D} - F_{S}) + [d R] (1 + \frac{[d r]}{[d -]}) (F_{R} - F_{D}) \\ + ([D R] (1 - \frac{1}{[D -]}) + [D r] (\frac{[d R]}{[d -]} - \frac{[- R]}{[D -]})) (F_{R D} - F_{D}) \\ [d R] + \frac{[d R] [d r]}{[d -]} = [d R] + \frac{[d R] [d r]}{[d R] + [d r]} = \frac{[d R] ([d R] + [d r]) + [d R] [d r]}{[d R] + [d r]} \\ = \frac{{[d R]}^{2} + [d R] [d r] + [d R] [d r]}{[d R] + [d r]} = \frac{{([d R] + [d r])}^{2} - {[d r]}^{2}}{[d R] + [d r]} = [d -] - \frac{{[d r]}^{2}}{[d -]} \\ ∴ {\bar{F}}_{D} & > {\bar{F}}_{d} \leftrightarrow (F_{D} - F_{S}) > ([D -] + [d -] - \frac{{[d r]}^{2}}{[d -]}) (F_{D} - F_{S}) + [d R] (1 + \frac{[d r]}{[d -]}) (F_{R} - F_{D}) \\ + ([D R] (1 - \frac{1}{[D -]}) + [D r] (\frac{[d R]}{[d -]} - \frac{[- R]}{[D -]})) (F_{R D} - F_{D}) \\ {\bar{F}}_{D} & > {\bar{F}}_{d} \leftrightarrow (F_{D} - F_{S}) > (1 - \frac{{[d r]}^{2}}{[d -]}) (F_{D} - F_{S}) + [d R] (1 + \frac{[d r]}{[d -]}) (F_{R} - F_{D}) \\ + ([D R] (1 - \frac{1}{[D -]}) + [D r] (\frac{[d R]}{[d -]} - \frac{[- R]}{[D -]})) (F_{R D} - F_{D}) \\ {\bar{F}}_{D} & > {\bar{F}}_{d} \leftrightarrow \frac{{[d r]}^{2}}{[d -]} (F_{D} - F_{S}) > [d R] (1 + \frac{[d r]}{[d -]}) (F_{R} - F_{D}) \\ + ([D R] (1 - \frac{1}{[D -]}) + [D r] (\frac{[d R]}{[d -]} - \frac{[- R]}{[D -]})) (F_{R D} - F_{D}) \end{aligned}

Which is Expression 1

The notation and assumptions used to develop Expression 2, the inequality indicating when resistance alleles will be inherited into offspring with an average fitness higher than susceptible alleles, are similar to those used for the equivalent expression with respect to deflection alleles.

Consider the allele pairs in gametes in the mating adult population:

\begin{array}{ll} [R d] & = p r o p o r t i o n o f r e s i s t a n c e l o c u s a l l e l e s w h i c h a r e r e s i s t a n c e a l l e l e s p a i r e d w i t h \\ n o n - d e f l e c t i o n a l l e l e s i n g a m e t e s \\ [R D] & = p r o p o r t i o n o f r e s i s t a n c e l o c u s a l l e l e s w h i c h a r e r e s i s t a n c e a l l e l e s p a i r e d w i t h \\ d e f l e c t i o n a l l e l e s i n g a m e t e s \\ [r d] & = p r o p o r t i o n o f r e s i s t a n c e l o c u s a l l e l e s w h i c h a r e s u s c e p t i b l e a l l e l e s p a i r e d w i t h \\ n o n - d e f l e c t i o n a l l e l e s i n g a m e t e s \\ [r D] & = p r o p o r t i o n o f r e s i s t a n c e l o c u s a l l e l e s w h i c h a r e s u s c e p t i b l e a l l e l e s p a i r e d w i t h \\ d e f l e c t i o n a l l e l e s i n g a m e t e s \\ [R -] & = p r o p o r t i o n o f r e s i s t a n c e l o c u s a l l e l e s i n g a m e t e s w h i c h a r e r e s i s t a n c e a l l e l e s \\ [r -] & = p r o p o r t i o n o f r e s i s t a n c e l o c u s a l l e l e s i n g a m e t e s w h i c h a r e s u s c e p t i b l e a l l e l e s \\ [- d] & = p r o p o r t i o n o f r e s i s t a n c e l o c u s a l l e l e s w h i c h a r e p a i r e d w i t h n o n - d e f l e c t i o n a l l e l e s i n g a m e t e s \\ [- D] & = p r o p o r t i o n o f r e s i s t a n c e l o c u s a l l e l e s w h i c h a r e p a i r e d w i t h d e f l e c t i o n a l l e l e s i n g a m e t e s \\ [r -] + [R -] = [- d] + [- D] = 1 \end{array}

Let ${\bar{F}}_{R}$ represent the average fitness of offspring into which resistance alleles are inherited:

\begin{aligned} {\bar{F}}_{R} & = \frac{1}{[R -]} (\begin{matrix} [R d] [R d] F_{R} + [R d] [R D] F_{R D} + [R d] [r d] F_{R} + [R d] [r D] F_{R D} \\ + [R D] [R d] F_{R D} + [R D] [R D] F_{R D} + [R D] [r d] F_{R D} + [R D] [r D] F_{R D} \end{matrix}) \\ {\bar{F}}_{R} & = \frac{1}{[R -]} (\begin{matrix} ([R d] [R d] + [R d] [r d]) F_{R} \\ + ([R d] [R D] + [R d] [r D] + [R D] [R d] + [R D] [R D] + [R D] [r d] + [R D] [r D]) F_{R D} \end{matrix}) \\ {\bar{F}}_{R} & = \frac{1}{[R -]} ([R d] [- d] F_{R} + ([R d] [- D] + [R D]) F_{R D}) \\ {\bar{F}}_{R} & = \frac{1}{[R -]} ([R d] [- d] F_{R} + ([R -] [- D] + [R D] [- d]) F_{R D}) \\ {\bar{F}}_{R} & = \frac{[R d] [- d]}{[R -]} (F_{R} - F_{S}) + ([- D] + \frac{[R D] [- d]}{[R -]}) (F_{R D} - F_{S}) + F_{S} \end{aligned}

The average fitness of offspring into which susceptible alleles are inherited, ${\bar{F}}_{r}$ , is:

\begin{aligned} {\bar{F}}_{r} & = \frac{1}{[r -]} (\begin{matrix} [r d] [R d] F_{R} + [r d] [R D] F_{R D} + [r d] [r d] F_{S} + [r d] [r D] F_{D} \\ + [r D] [R d] F_{R D} + [r D] [R D] F_{R D} + [r D] [r d] F_{D} + [r D] [r D] F_{D} \end{matrix}) \\ {\bar{F}}_{r} & = \frac{1}{[r -]} (\begin{matrix} [r d] [R d] F_{R} + ([r d] [R D] + [r D] [R d] + [r D] [R D]) F_{R D} \\ + [r d] [r d] F_{S} + ([r d] [r D] + [r D] [r d] + [r D] [r D]) F_{D} \end{matrix}) \\ {\bar{F}}_{r} & = \frac{[r d] [R d]}{[r -]} F_{R} + ([R D] + \frac{[r D] [R d]}{[r -]}) F_{R D} + \frac{[r d] [r d]}{[r -]} F_{S} + ([r D] + \frac{([r d] [r D])}{[r -]}) F_{D} \\ {\bar{F}}_{r} & = \frac{[r d] [R d]}{[r -]} (F_{R} - F_{S}) + ([R D] + \frac{[r D] [R d]}{[r -]}) (F_{R D} - F_{S}) + ([r D] + \frac{([r d] [r D])}{[r -]}) (F_{D} - F_{S}) + F_{S} \end{aligned}

In order for the proportion of resistance alleles in the population to increase, we need the average fitness of the offspring into which resistance alleles are inherited to be greater than the average fitness of the offspring into which susceptible alleles are inherited. This is true when the following inequality applies:

\begin{aligned} {\bar{F}}_{R} & > {\bar{F}}_{r} \leftrightarrow \frac{[R d] [- d]}{[R -]} (F_{R} - F_{S}) + ([- D] + \frac{[R D] [- d]}{[R -]}) (F_{R D} - F_{S}) + F_{S} \\ - \frac{[r d] [R d]}{[r -]} (F_{R} - F_{S}) - ([R D] + \frac{[r D] [R d]}{[r -]}) (F_{R D} - F_{S}) - ([r D] + \frac{([r d] [r D])}{[r -]}) (F_{D} - F_{S}) - F_{S} > 0 \\ \leftrightarrow (\frac{[R d] [- d]}{[R -]} - \frac{[R d] [r d]}{[r -]}) (F_{R} - F_{D}) + (\frac{[R d] [- d]}{[R -]} - \frac{[R d] [r d]}{[r -]} - [r D] - \frac{([r d] [r D])}{[r -]}) (F_{D} - F_{S}) \\ + ([r D] + \frac{[R D] [- d]}{[R -]} - \frac{[r D] [R d]}{[r -]}) (F_{R D} - F_{R}) + ([r D] + \frac{[R D] [- d]}{[R -]} - \frac{[r D] [R d]}{[r -]}) (F_{R} - F_{D}) \\ + ([r D] + \frac{[R D] [- d]}{[R -]} - \frac{[r D] [R d]}{[r -]}) (F_{D} - F_{S}) > 0 \\ \leftrightarrow (\frac{[R d] [- d]}{[R -]} - \frac{[R d] [r d]}{[r -]} + [r D] + \frac{[R D] [- d]}{[R -]} - \frac{[r D] [R d]}{[r -]}) (F_{R} - F_{D}) \\ + (\frac{[R d] [- d]}{[R -]} - \frac{[R d] [r d]}{[r -]} - [r D] - \frac{([r d] [r D])}{[r -]} + [r D] + \frac{[R D] [- d]}{[R -]} - \frac{[r D] [R d]}{[r -]}) (F_{D} - F_{S}) \\ + ([r D] + \frac{[R D] [- d]}{[R -]} - \frac{[r D] [R d]}{[r -]}) (F_{R D} - F_{R}) > 0 \\ \leftrightarrow ([r -]) (F_{R} - F_{D}) + ([r d] - \frac{([r d] [r D])}{[r -]}) (F_{D} - F_{S}) \\ + ([r D] + \frac{[R D] [- d]}{[R -]} - \frac{[r D] [R d]}{[r -]}) (F_{R D} - F_{R}) > 0 \\ \leftrightarrow [r d] (F_{R} - F_{S}) + [r D] (F_{R} - F_{D}) \\ + ([r D] + \frac{[R D] [- d]}{[R -]} - \frac{[r D] [R d]}{[r -]}) (F_{R D} - F_{R}) > \frac{([r d] [r D])}{[r -]} (F_{D} - F_{S}) \end{aligned}

This is Expression 2.

References

(2015) Repellent efficacy of DEET, MyggA, neem (Azedirachta indica) oil and chinaberry (Melia azedarach) oil against Anopheles arabiensis, the principal malaria vector in Ethiopia
Malaria Journal 14:187.

https://doi.org/10.1186/s12936-015-0705-4
- Google Scholar
1. Achee NL
2. Bangs MJ
3. Farlow R
4. Killeen GF
5. Lindsay S
6. Logan JG
7. Moore SJ
8. Rowland M
9. Sweeney K
10. Torr SJ
11. Zwiebel LJ
12. Grieco JP
(2012) Spatial repellents: from discovery and development to evidence-based validation
Malaria Journal 11:164.

https://doi.org/10.1186/1475-2875-11-164
- Google Scholar
1. Asidi A
2. N'Guessan R
3. Akogbeto M
4. Curtis C
5. Rowland M
(2012) Loss of household protection from use of insecticide-treated nets against pyrethroid-resistant mosquitoes, benin
Emerging Infectious Diseases 18:1101–1106.

https://doi.org/10.3201/eid1807.120218
- Google Scholar
1. Bhatt S
2. Gething PW
3. Brady OJ
4. Messina JP
5. Farlow AW
6. Moyes CL
7. Drake JM
8. Brownstein JS
9. Hoen AG
10. Sankoh O
11. Myers MF
12. George DB
13. Jaenisch T
14. Wint GR
15. Simmons CP
16. Scott TW
17. Farrar JJ
18. Hay SI
(2013) The global distribution and burden of dengue
Nature 496:504–507.

https://doi.org/10.1038/nature12060
- Google Scholar
1. Boots M
(2015) The need for evolutionarily rational disease interventions: vaccination can select for higher virulence
PLoS Biology 13:e1002236.

https://doi.org/10.1371/journal.pbio.1002236
- Google Scholar
1. Bradley J
2. Matias A
3. Schwabe C
4. Vargas D
5. Monti F
6. Nseng G
7. Kleinschmidt I
(2012) Increased risks of malaria due to limited residual life of insecticide and outdoor biting versus protection by combined use of nets and indoor residual spraying on Bioko Island, Equatorial Guinea
Malaria Journal 11:242.

https://doi.org/10.1186/1475-2875-11-242
- Google Scholar
1. Brogdon WG
2. McAllister JC
(1998) Insecticide resistance and vector control
Emerging Infectious Diseases 4:605.

https://doi.org/10.3201/eid0404.980410
- Google Scholar
(1988) Physiological and behavioral adaptation to varying environments: a mathematical model
Evolution 42:986–994.

https://doi.org/10.2307/2408914
- Google Scholar
1. Chandre F
2. Darrier F
3. Manga L
4. Akogbeto M
5. Faye O
6. Mouchet J
7. Guillet P
(1999)
Status of pyrethroid resistance in Anopheles gambiae sensu lato

Bulletin of the World Health Organization 77:230–234.
- Google Scholar
1. Chandre F
2. Darriet F
3. Duchon S
4. Finot L
5. Manguin S
6. Carnevale P
7. Guillet P
(2000) Modifications of pyrethroid effects associated with kdr mutation in Anopheles gambiae
Medical and Veterinary Entomology 14:81–88.

https://doi.org/10.1046/j.1365-2915.2000.00212.x
- Google Scholar
1. Cooke MK
2. Kahindi SC
3. Oriango RM
4. Owaga C
5. Ayoma E
6. Mabuka D
7. Nyangau D
8. Abel L
9. Atieno E
10. Awuor S
11. Drakeley C
12. Cox J
13. Stevenson J
(2015) 'A bite before bed': exposure to malaria vectors outside the times of net use in the highlands of western Kenya
Malaria Journal 14:1–15.

https://doi.org/10.1186/s12936-015-0766-4
- Google Scholar
1. Curtis CF
(2002)
Restoration of malaria control in the Madagascar highlands by DDT spraying

The American Journal of Tropical Medicine and Hygiene 66:1.
- Google Scholar
1. Dabiré RK
2. Namountougou M
3. Sawadogo SP
4. Yaro LB
5. Toé HK
6. Ouari A
7. Gouagna LC
8. Simard F
9. Chandre F
10. Baldet T
11. Bass C
12. Diabaté A
(2012) Population dynamics of Anopheles gambiae s.l. in Bobo-Dioulasso city: bionomics, infection rate and susceptibility to insecticides
Parasites & Vectors 5:127.

https://doi.org/10.1186/1756-3305-5-127
- Google Scholar
1. Evans RG
(1993)
Laboratory evaluation of the irritancy of bendiocarb, lambda-cyhalothrin and DDT to Anopheles gambiae

Journal of the American Mosquito Control Association 9:285–293.
- Google Scholar
(2012) Strain improvement of fungal insecticides for controlling insect pests and vector-borne diseases
Current Opinion in Microbiology 15:232–238.

https://doi.org/10.1016/j.mib.2011.12.012
- Google Scholar
(2010) Insecticidal, acaricidal and repellent effects of DEET- and IR3535-impregnated bed nets using a novel long-lasting polymer-coating technique
Parasitology Research 106:957–965.

https://doi.org/10.1007/s00436-010-1749-6
- Google Scholar
1. Gatton ML
2. Chitnis N
3. Churcher T
4. Donnelly MJ
5. Ghani AC
6. Godfray HCJ
7. Gould F
8. Hastings I
9. Marshall J
10. Ranson H
11. Rowland M
12. Shaman J
13. Lindsay SW
(2013) The importance of mosquito behavioural adaptations to malaria control in Africa
Evolution 67:1218–1230.

https://doi.org/10.1111/evo.12063
- Google Scholar
1. Gautret P
2. Motard A
(1999) Periodic infectivity of Plasmodium gametocytes to the vector. A review
Parasite 6:103–111.

https://doi.org/10.1051/parasite/1999062103
- Google Scholar
1. Georghiou GP
2. Taylor CE
(1977) Genetic and biological influences in the evolution of insecticide resistance
Journal of Economic Entomology 70:319–323.

https://doi.org/10.1093/jee/70.3.319
- Google Scholar
1. Gould F
(1984) Role of behavior in the evolution of insect adaptation to insecticides and resistant host plants
Bulletin of the Entomological Society of America 30:34–41.

https://doi.org/10.1093/besa/30.4.34
- Google Scholar
1. Govella NJ
2. Ogoma SB
3. Paliga J
4. Chaki PP
5. Killeen G
(2015) Impregnating hessian strips with the volatile pyrethroid transfluthrin prevents outdoor exposure to vectors of malaria and lymphatic filariasis in urban Dar es Salaam, Tanzania
Parasites & Vectors 8:1–9.

https://doi.org/10.1186/s13071-015-0937-8
- Google Scholar
(2000)
A comparison study of house entering and exiting behavior of Anopheles vestitipennis (Diptera: Culicidae) using experimental huts sprayed with DDT or deltamethrin in the southern district of Toledo, Belize, C.A

Journal of Vector Ecology 25:62–73.
- Google Scholar
(2007) A new classification system for the actions of IRS chemicals traditionally used for malaria control
PLoS One 2:e716.

https://doi.org/10.1371/journal.pone.0000716
- Google Scholar
(2002) An overview of insecticide resistance
Science 298:96–97.

https://doi.org/10.1126/science.1078052
- Google Scholar
1. Hemingway J
2. Ranson H
(2000) Insecticide resistance in insect vectors of human disease
Annual Review of Entomology 45:371–391.

https://doi.org/10.1146/annurev.ento.45.1.371
- Google Scholar
1. Jones CM
2. Machin C
3. Mohammed K
4. Majambere S
5. Ali AS
6. Khatib BO
7. McHa J
8. Ranson H
9. Kelly-Hope LA
(2012) Insecticide resistance in Culex quinquefasciatus from Zanzibar: implications for vector control programmes
Parasites & Vectors 5:78.

https://doi.org/10.1186/1756-3305-5-78
- Google Scholar
1. Jongsma MA
2. Gould F
3. Legros M
4. Yang L
5. van Loon JJA
6. Dicke M
(2010) Insect oviposition behavior affects the evolution of adaptation to Bt crops: consequences for refuge policies
Evolutionary Ecology 24:1017–1030.

https://doi.org/10.1007/s10682-010-9368-3
- Google Scholar
(1987) Ecological, agricultural, genetic, and commercial considerations in the deployment of insect-resistant germplasm
Environmental Entomology 16:327–338.

https://doi.org/10.1093/ee/16.2.327
- Google Scholar
1. Kennedy JS
(1947) The excitant and repellent effects on mosquitos of sub-lethal contacts with DDT
Bulletin of Entomological Research 37:593–607.

https://doi.org/10.1017/S0007485300030091
- Google Scholar
(2011) Target product profile choices for intra-domiciliary malaria vector control pesticide products: repel or kill?
Malaria Journal 10:207.

https://doi.org/10.1186/1475-2875-10-207
- Google Scholar
1. Le Menach A
2. Takala S
3. McKenzie FE
4. Perisse A
5. Harris A
6. Flahault A
7. Smith DL
(2007) An elaborated feeding cycle model for reductions in vectorial capacity of night-biting mosquitoes by insecticide-treated nets
Malaria Journal 6:10.

https://doi.org/10.1186/1475-2875-6-10
- Google Scholar
1. Lefèvre T
2. Gouagna LC
3. Dabiré KR
4. Elguero E
5. Fontenille D
6. Renaud F
7. Costantini C
8. Thomas F
(2009) Beyond nature and nurture: phenotypic plasticity in blood-feeding behavior of Anopheles gambiae s.s. when humans are not readily accessible
American Journal of Tropical Medicine and Hygiene 81:1023–1029.

https://doi.org/10.4269/ajtmh.2009.09-0124
- Google Scholar
(2006) Lethal and behavioural effects of three synthetic repellents (DEET, IR3535 and KBR 3023) on Aedes aegypti mosquitoes in laboratory assays
Medical and Veterinary Entomology 20:288–293.

https://doi.org/10.1111/j.1365-2915.2006.00630.x
- Google Scholar
(1990)
Effect of indoor residual spraying of DDT and bendiocarb on the feeding patterns of Anopheles pseudopunctipennis in Mexico

Journal of the American Mosquito Control Association 6:635–640.
- Google Scholar
1. Lynch PA
2. Grimm U
3. Thomas MB
4. Read AF
(2012) Prospective malaria control using entomopathogenic fungi: comparative evaluation of impact on transmission and selection for resistance
Malaria Journal 11:383.

https://doi.org/10.1186/1475-2875-11-383
- Google Scholar
1. Malima RC
2. Oxborough RM
3. Tungu PK
4. Maxwell C
5. Lyimo I
6. Mwingira V
7. Mosha FW
8. Matowo J
9. Magesa SM
10. Rowland MW
(2009) Behavioural and insecticidal effects of organophosphate-, carbamate- and pyrethroid-treated mosquito nets against African malaria vectors
Medical and Veterinary Entomology 23:317–325.

https://doi.org/10.1111/j.1365-2915.2009.00837.x
- Google Scholar
(2011) Mathematical models of malaria--a review
Malaria Journal 10:202.

https://doi.org/10.1186/1475-2875-10-202
- Google Scholar
1. Menger DJ
2. Omusula P
3. Wouters K
4. Oketch C
5. Carreira AS
6. Durka M
7. Derycke JL
8. Loy DE
9. Hahn BH
10. Mukabana WR
11. Mweresa CK
12. van Loon JJ
13. Takken W
14. Hiscox A
(2016) Eave screening and push-pull tactics to reduce house entry by vectors of malaria
American Journal of Tropical Medicine and Hygiene 94:868–878.

https://doi.org/10.4269/ajtmh.15-0632
- Google Scholar
1. Mideo N
2. Reece SE
3. Smith AL
4. Metcalf CJ
(2013) The Cinderella syndrome: why do malaria-infected cells burst at midnight?
Trends in Parasitology 29:10–16.

https://doi.org/10.1016/j.pt.2012.10.006
- Google Scholar
(2009) Designation of chemicals in terms of the locomotor responses they elicit from insects: an update of Dethier et al. (1960)
Journal of Economic Entomology 102:2056–2060.

https://doi.org/10.1603/029.102.0606
- Google Scholar
(2010) Anopheline and culicine mosquitoes are not repelled by surfaces treated with the entomopathogenic fungi Metarhizium anisopliae and Beauveria bassiana
Parasites & Vectors 3:80.

https://doi.org/10.1186/1756-3305-3-80
- Google Scholar
Book
1. Mouchet J
2. Hamon J
(1963)
Difficulties in Malaria Eradication Campaigns Due to the Behaviour of the Vectors

Geneva: World Health Organization.
- Google Scholar
(2013) Epidemiology of dengue: past, present and future prospects
Clinical Epidemiology 5:299–309.

https://doi.org/10.2147/CLEP.S34440
- Google Scholar
1. Müller GC
2. Junnila A
3. Butler J
4. Kravchenko VD
5. Revay EE
6. Weiss RW
7. Schlein Y
(2009) Efficacy of the botanical repellents geraniol, linalool, and citronella against mosquitoes
Journal of Vector Ecology 34:2–8.

https://doi.org/10.1111/j.1948-7134.2009.00002.x
- Google Scholar
(2007) Reduced efficacy of insecticide-treated nets and indoor residual spraying for malaria control in pyrethroid resistance area, Benin
Emerging Infectious Diseases 13:199–206.

https://doi.org/10.3201/eid1302.060631
- Google Scholar
(2008) DEET microencapsulation: a slow-release formulation enhancing the residual efficacy of bed nets against malaria vectors
Transactions of the Royal Society of Tropical Medicine and Hygiene 102:259–262.

https://doi.org/10.1016/j.trstmh.2007.10.012
- Google Scholar
(2006) Evaluation of synthetic repellents on mosquito nets in experimental huts against insecticide-resistant anopheles gambiae and culex quinquefasciatus mosquitoes
Transactions of the Royal Society of Tropical Medicine and Hygiene 100:1091–1097.

https://doi.org/10.1016/j.trstmh.2006.04.005
- Google Scholar
1. Nauen R
(2007) Insecticide resistance in disease vectors of public health importance
Pest Management Science 63:628–633.

https://doi.org/10.1002/ps.1406
- Google Scholar
1. Ndiath MO
2. Mazenot C
3. Sokhna C
4. Trape JF
(2014) How the malaria vector anopheles gambiae adapts to the use of insecticide-treated nets by african populations
PLoS One 9:e97700.

https://doi.org/10.1371/journal.pone.0097700
- Google Scholar
1. Ndiath MO
2. Sougoufara S
3. Gaye A
4. Mazenot C
5. Konate L
6. Faye O
7. Faye O
8. Sokhna C
9. Trape JF
(2012) Resistance to DDT and pyrethroids and increased kdr mutation frequency in an. gambiae after the implementation of permethrin-treated nets in senegal
PLoS One 7:e31943.

https://doi.org/10.1371/annotation/855f784c-b8db-449c-9de8-83b6404a049d
- Google Scholar
(2011) Fitness costs of disrupting circadian rhythms in malaria parasites
Proceedings of the Royal Society B 278:2429–2436.

https://doi.org/10.1098/rspb.2010.2457
- Google Scholar
1. Ochomo E
2. Bayoh NM
3. Kamau L
4. Atieli F
5. Vulule J
6. Ouma C
7. Ombok M
8. Njagi K
9. Soti D
10. Mathenge E
11. Muthami L
12. Kinyari T
13. Subramaniam K
14. Kleinschmidt I
15. Donnelly MJ
16. Mbogo C
(2014) Pyrethroid susceptibility of malaria vectors in four Districts of western Kenya
Parasites & Vectors 7:310.

https://doi.org/10.1186/1756-3305-7-310
- Google Scholar
1. Ogoma SB
2. Ngonyani H
3. Simfukwe ET
4. Mseka A
5. Moore J
6. Maia MF
7. Moore SJ
8. Lorenz LM
(2014) The mode of action of spatial repellents and their impact on vectorial capacity of Anopheles gambiae sensu stricto
PLoS One 9:e110433.

https://doi.org/10.1371/journal.pone.0110433
- Google Scholar
1. Patrick M
(1900) Experimental proof of the mosquitomalaria theory
BMJ 2:949.

https://doi.org/10.1136/bmj.2.2074.949
- Google Scholar
1. Pennetier C
2. Corbel V
3. Boko P
4. Odjo A
5. N'Guessan R
6. Lapied B
7. Hougard JM
(2007) Synergy between repellents and non-pyrethroid insecticides strongly extends the efficacy of treated nets against Anopheles gambiae
Malaria Journal 6:1–7.

https://doi.org/10.1186/1475-2875-6-38
- Google Scholar
(2009) How to make evolution-proof insecticides for malaria control
PLoS Biology 7:e58.

https://doi.org/10.1371/journal.pbio.1000058
- Google Scholar
1. Reddy MR
2. Overgaard HJ
3. Abaga S
4. Reddy VP
5. Caccone A
6. Kiszewski AE
7. Slotman MA
(2011) Outdoor host seeking behaviour of Anopheles gambiae mosquitoes following initiation of malaria vector control on bioko island, equatorial guinea
Malaria Journal 10:184.

https://doi.org/10.1186/1475-2875-10-184
- Google Scholar
1. Revay EE
2. Kline DL
3. Xue RD
4. Qualls WA
5. Bernier UR
6. Kravchenko VD
7. Ghattas N
8. Pstygo I
9. Müller GC
(2013) Reduction of mosquito biting-pressure: spatial repellents or mosquito traps? A field comparison of seven commercially available products in Israel
Acta Tropica 127:63–68.

https://doi.org/10.1016/j.actatropica.2013.03.011
- Google Scholar
Book
1. Roberts D
(1993)
Insecticide Repellency in Malaria Vector Control: A Position Paper

VBC Project report No 81131.
- Google Scholar
(2000)
A probability model of vector behavior: effects of DDT repellency, irritancy, and toxicity in malaria control

Journal of Vector Ecology 25:48–61.
- Google Scholar
1. Roberts DR
2. Alecrim WD
(1991)
Behavioral response of Anopheles darlingi to DDT-sprayed house walls in Amazonia

Bulletin of the Pan American Health Organization 25:210–217.
- Google Scholar
(1997) DDT, global strategies, and a malaria control crisis in South America
Emerging Infectious Diseases 3:295–302.

https://doi.org/10.3201/eid0303.970305
- Google Scholar
(2000) DDT house spraying and re-emerging malaria
The Lancet 356:330–332.

https://doi.org/10.1016/S0140-6736(00)02516-2
- Google Scholar
1. Roberts DR
(2010) Impact of anti-DDT campaigns on malaria control
Outlooks on Pest Management 21:4–11.

https://doi.org/10.1564/21feb02
- Google Scholar
1. Roll Back Malaria Vector Control Working Group
(2015)
Housing and malaria consensus statement
- Google Scholar
1. Rosenheim J
2. Tabashnik BE
(1990) Evolution of pesticide resistance: interactions between generation time and genetic, ecological, and operational factors
Journal of Economic Entomology 83:1184–1193.

https://doi.org/10.1093/jee/83.4.1184
- Google Scholar
(2011) Increased proportions of outdoor feeding among residual malaria vector populations following increased use of insecticide-treated nets in rural Tanzania
Malaria Journal 10:80.

https://doi.org/10.1186/1475-2875-10-80
- Google Scholar
1. Syafruddin D
2. Bangs MJ
3. Sidik D
4. Elyazar I
5. Asih PB
6. Chan K
7. Nurleila S
8. Nixon C
9. Hendarto J
10. Wahid I
11. Ishak H
12. Bøgh C
13. Grieco JP
14. Achee NL
15. Baird JK
(2014) Impact of a spatial repellent on malaria incidence in two villages in Sumba, Indonesia
American Journal of Tropical Medicine and Hygiene 91:1079–1087.

https://doi.org/10.4269/ajtmh.13-0735
- Google Scholar
(2012) Pyrethroid induced behavioral responses of Anopheles dirus, a vector of malaria in Thailand
Journal of Vector Ecology 37:187–196.

https://doi.org/10.1111/j.1948-7134.2012.00216.x
- Google Scholar
1. Thomas MB
2. Read AF
(2007) Can fungal biopesticides control malaria?
Nature Reviews Microbiology 5:377–383.

https://doi.org/10.1038/nrmicro1638
- Google Scholar
Website
1. University of Notre Dame
(2016) Spatial repellents for control of vector-borne disease

http://ndigd.nd.edu/projects/data/spatial-repellents-for-control-of-vector-borne-disease/
(2014) Negative cross resistance mediated by co-treated bed nets: a potential means of restoring pyrethroid-susceptibility to malaria vectors
PLoS One 9:e95640.

https://doi.org/10.1371/journal.pone.0095640
- Google Scholar
1. WHO
(2015)
Global Technical Strategy for Malaria 2016–2030

Global Technical Strategy for Malaria 2016–2030.
- Google Scholar
1. WHO
(2015)
Vector control and insecticide resistance
- Google Scholar
1. Xue R
2. Ali A
3. Barnard DR
(2003)
Scientific notes-laboratory evaluation of toxicity of 16 insect repellents in aerosol sprays to adult mosquitoes

Journal of the American Mosquito Control Association. 19:271–274.
- Google Scholar
1. Yakob L
2. Dunning R
3. Yan G
(2011) Indoor residual spray and insecticide-treated bednets for malaria control: theoretical synergisms and antagonisms
Journal of the Royal Society Interface 8:799–806.

https://doi.org/10.1098/rsif.2010.0537
- Google Scholar

Article and author information

Author details

Penelope Anne Lynch

Department of Biosciences, University of Exeter, Cornwall Campus, Penryn, United Kingdom

Contribution
PAL, Conception, design and modelling, Drafting or revising the article

For correspondence
P.A.Lynch@Exeter.ac.uk

Competing interests
The authors declare that no competing interests exist.

"This ORCID iD identifies the author of this article:" 0000-0002-0486-8507
Mike Boots
1. Department of Biosciences, University of Exeter, Cornwall Campus, Penryn, United Kingdom
2. Department of Integrative Biology, University of California, Berkeley, United States
Contribution
MB, Design, Drafting or revising the article

Competing interests
The authors declare that no competing interests exist.

Funding

Biotechnology and Biological Sciences Research Council (BB/L010879/1)

Penelope Anne Lynch
Mike Boots

NERC Environmental Bioinformatics Centre (NE/J009784/1)

Penelope Anne Lynch
Mike Boots

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

The authors would like to thank Andrew Read for early discussions of the ESR concept. We would also like to thank the editors and reviewers at eLife who were extremely helpful and supportive during the process of bringing the manuscript to publication. They helped to improve the manuscript substantially and they introduced a number of interesting ideas, including the potential for ESR to support repellents for use outdoors.

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

2,106

views
343

downloads
16

citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Mendeley

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Penelope Anne Lynch
Mike Boots

(2016)

Using evolution to generate sustainable malaria control with spatial repellents

eLife 5:e15416.

https://doi.org/10.7554/eLife.15416

Abstract

eLife digest

Introduction

Modeling

Spread of deflection in a population over time.

Results

Combinations of per cycle survival values for deflected and non-deflected resistant phenotypes, which support the spread and maintenance of deflected phenotypes in the population.

Effect of cost of resistance on resistant and deflected per-cycle survival combinations which support the spread and maintenance of deflected phenotypes in the population.

Effect of evolved spatial repellent on infectious bites from a vector population over time.

Discussion

Appendix 1

Definition of the PM model

Appendix 2

Derivation of Equations 1 and 2 and Expressions 1 and 2

References

Article and author information

Author details

Penelope Anne Lynch

Contribution

For correspondence

Competing interests

Mike Boots

Contribution

Competing interests

Funding

Biotechnology and Biological Sciences Research Council (BB/L010879/1)

NERC Environmental Bioinformatics Centre (NE/J009784/1)

Acknowledgements

Copyright

Metrics

Download links

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Categories and tags

Further reading

Livestock abortion surveillance in Tanzania reveals disease priorities and importance of timely collection of vaginal swab samples for attribution

Systematic evaluation of multifactorial causal associations for Alzheimer’s disease and an interactive platform MRAD developed based on Mendelian randomization analysis

Be the first to read new articles from eLife

Share this article

Cite this article

Spread of deflection in a population over time.

Combinations of per cycle survival values for deflected and non-deflected resistant phenotypes, which support the spread and maintenance of deflected phenotypes in the population.

Effect of cost of resistance on resistant and deflected per-cycle survival combinations which support the spread and maintenance of deflected phenotypes in the population.

Effect of evolved spatial repellent on infectious bites from a vector population over time.

Author details

Penelope Anne Lynch

Contribution

For correspondence

Competing interests

Mike Boots

Contribution

Competing interests

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Categories and tags

Further reading