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 hightoxicity 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, ‘doubledip’ 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.001eLife 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 insecticideresistant 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.002Introduction
Vectorborne 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, contactirritance and spatial repellence. Toxicity is simply the capacity to generate mortality in mosquitoes that contact the insecticide. Contactirritance 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 indoorfeeding 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 vectorfree 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 proofofconcept doubleblind placebocontrolled field trial using spatial repellents disbursed from burning coils. This trial recorded an approximately 50% reduction in malaria infections in repellentprotected 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 longstanding 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 propertyscale 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 lowcontactrepellence hightoxicity insecticide, which kills those which do enter, leading to a ‘doubledip’ 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 insecticideonly 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 diseasecontrol 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 wellestablished 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 feedingcycle based, twolocus, biallelic population genetics model, capturing the mosquito lifehistory characteristics of one initial mating and preadult 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 nondeflected 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 perfeed 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 nonhuman and therefore noninfectious 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 outdoorfeeding vectors.
Our key model assumptions. (1) Physiological insecticide resistance is controlled by a single locus biallelic 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 biallelic autosomal gene, with the deflection allele being completely dominant to the nondeflection 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 reassort 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 singlelocus insecticideresistance 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).
The average fitness of offspring into which deflection alleles are inherited, $\overline{F}}_{D$, is
With $[dr],\phantom{\rule{thinmathspace}{0ex}}\phantom{\rule{thinmathspace}{0ex}}[dR],\phantom{\rule{thinmathspace}{0ex}}\phantom{\rule{thinmathspace}{0ex}}[Dr],\phantom{\rule{thinmathspace}{0ex}}\phantom{\rule{thinmathspace}{0ex}}[DR],\phantom{\rule{thinmathspace}{0ex}}\phantom{\rule{thinmathspace}{0ex}}[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 nondeflection alleles paired with susceptible alleles, nondeflection alleles paired with resistant alleles, deflection alleles paired with susceptible alleles, deflection alleles paired with resistant alleles, nondeflection 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 nondeflection alleles are inherited, $\overline{F}}_{d$, is
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 nondeflection alleles are inherited. This is true when the following inequality applies:
The equivalent expression for spread of the resistance allele is:
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 insecticidetreated 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 lefthand side of the expression), and by maximising the fitness of deflected relative to resistant phenotypes (reducing the value of the righthand side of the expression). Avoiding the use of ESR in properties without insecticide, minimising $Y}_{3$ (Table 2), and using ESR in all insecticidetreated 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.
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 nondeflection 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 deflectedresistant 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 nondeflection alleles that are paired with resistance alleles. These interactions are explored in Figures 1–4.
In the absence of resistance, Expression 1 reduces to ${\overline{F}}_{D}>{\overline{F}}_{d}\leftrightarrow \left[d\right]\left({F}_{D}{F}_{S}\right)>0$, so deflection would be expected to spread to fixation provided that deflected phenotypes have greater fitness than nondeflected phenotypes. If deflection reaches fixation, then Expression 2 reduces to ${\overline{F}}_{R}>{\overline{F}}_{r}\leftrightarrow [rD]({F}_{RD}{F}_{D})>0$, and resistance will only spread if resistantdeflected phenotypes are fitter than deflected phenotypes. In this instance, the strategy of using ESR in all insecticidetreated properties would prevent further spread of resistance because deflected mosquitoes would not enter any insecticidetreated 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 nonresistant deflected phenotypes (depending upon cost of resistance), giving $F}_{RD}\le {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 populationlevel 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 nondeflected 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 nearfixation, 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.
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 insecticideresistant phenotypes have average percycle survival of 60% and deflected phenotypes have percycle 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 percycle 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.
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 wellexplored. 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 outdoorfeeding 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 tradeoff between immediate and longterm 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.
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 transmissionreduction 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 malariavector 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 hightoxicity 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 wellknown of the available repellents, as well as many of the naturally repellent plantbased 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 publichealth 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 publichealth 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 resistantplusdeflected phenotypes, and by the initial prevalence of deflection and resistance alleles when the combined ESRinsecticide 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 highefficacy insecticides to minimize the fitness of nondeflected 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 insecticidetreated 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 nondeflected susceptible mosquitoes which can enter and feed. The use of a wellestablished 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 hightoxicity 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 nondeflected 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 shortterm 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’ latelifeacting insecticides (Read et al., 2009; Thomas and Read, 2007), offering respite from the treadmill of insecticide loss to resistance. For latelifeacting insecticides, relatively lowvirulence fungal strains are ideal (Lynch et al., 2012), but much work on these organisms has focussed on generating highvirulence strains (Fang et al., 2012). Given their low or absent contact repellence, such highvirulence 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 indoorfeeding malaria vectors from entering properties, rather than on more localized repellence away from bed nets once vectors have entered properties. The excitorepellent properties of some of the pyrethroid insecticides currently used on longlasting insecticidetreated 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 nondeflected 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 publichealth 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 nondeflected 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 nondeflection 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 publichealth opportunities such as the shortterm use of a more expensive insecticide to generate ESR protection for the longterm 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 insecticidesusceptible An. gambiae pushed to outdoor feeding by the deployment of LLINs in Dielmo, whilst transmission continued at a high level in an indoorfeeding 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 nonhuman 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 wellestablished 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 antivector 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 indoorfeeding 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 insecticidetreated 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 resistancemanagement 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 singlespot 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 (CastilloChavez 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 nonrepellent formulations. Whereas use in agriculture is commonly highly detrimental to the sustained utility of publichealth 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 insecticidetreated crops, using a vectorcontrol 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 resistancemanagement for agricultural pesticides). This might also be a means to apply selection to nectarfeeding 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
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 oversubscribed, 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}={\displaystyle \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}{\beta}_{0}{b}_{0}$ and $m}_{h,t}={U}_{h,t}{\beta}_{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 $\gamma $ periods (number of periods between Plasmodium infection and infectiousness) ago, multiplied by the probability of surviving to $\gamma $ periods ago, multiplied by the probability of surviving $\gamma $ periods, $M}_{h,0}={\displaystyle \frac{{m}_{h,0}}{1{\beta}_{0}}{{\beta}_{0}}^{\gamma}$.
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 $\gamma $ periods ago. ${M}_{h,t}={M}_{h,(t1)}{\beta}_{h}+{m}_{h,(tmin(t,\gamma ))}{{\beta}_{h}}^{\gamma}\phantom{\rule{1em}{0ex}}\phantom{\rule{1em}{0ex}}\phantom{\rule{1em}{0ex}}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 noninfected 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}={\displaystyle \frac{{J}_{h,0}}{1{\beta}_{0}(1{b}_{0})}$. Postintroduction of the intervention, the noninfected females for each genotype at the start of each time period are calculated as the surviving noninfected females from the previous period plus new females of that genotype, giving $U}_{h,t}={U}_{h,t1}{\beta}_{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,
New deflected females in time period t as a proportion of the baseline population,
New resistant females in time period t as a proportion of the baseline population,
New deflectedresistant females in time period t as a proportion of the baseline population,
For this analysis, it is assumed that larval habitats are oversubscribed 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 $\tau}_{1,t}=1{\beta}_{0$.
${\tau}_{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 ${\tau}_{2,0}=1{\tau}_{1,0}$ and $\tau}_{2,t}=\sum _{g=1}^{9}({\tau}_{1,t1}{B}_{g,t1}+{\tau}_{2,t1}{K}_{g,t1}){S}_{g$ with $t\text{}\text{}0$.
New (female) adults in period $t$ as a proportion of the total of population in period $t$, $\tau}_{3,t}={\displaystyle \frac{{\tau}_{1,t}}{{\tau}_{1,t}+{\tau}_{2,t}}.$
Adult females surviving from the preceding period to the start of period $t$ as a proportion of the population in period $t$, $\tau}_{4,t}={\displaystyle \frac{{\tau}_{2,t}}{{\tau}_{1,t}+{\tau}_{2,t}}.$
Population of adult females in period $t$ as a proportion of the baseline population, $\tau}_{5,t}={\tau}_{1,t}+{\tau}_{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 HardyWeinberg equilibrium (independently for each locus).
The proportion of new adults in time period $t$ (with $t\text{}\text{}0$) that have genotype $g$ is calculated as the proportion of eggs with that genotype from $\varpi $ periods ago, where $\varpi $ represents the number of model time periods required for development from egg to new adults. For early time periods, before $t\text{}\text{}\varpi$, the baseline egg genotype proportions are used ($t=0$), giving $B}_{g,t}={E}_{g,t\mathrm{m}\mathrm{i}\mathrm{n}(t,\varpi )$.
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}={\displaystyle \frac{{W}_{g,t1}{S}_{g}}{{\psi}_{t}}\phantom{\rule{1em}{0ex}}\phantom{\rule{1em}{0ex}}t>\text{}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}{\tau}_{3,t}+{K}_{g,t}{\tau}_{4,t}.$
${\psi}_{t}$, the proportion of the population alive at start of time period $t1$ 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 $t1$ multiplied by the oneperiod survival probability associated with that genotype ${\psi}_{t}=\sum _{i=1}^{9}({W}_{i,t1}{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$, ${\theta}_{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\phantom{\rule{1em}{0ex}}{\theta}_{g,t}=\sum _{i=1}^{t}\u27ee\sum _{j=1}^{9}({B}_{j,i}{Z}_{g,j,i}{A}_{j,i}{{S}_{j}}^{ti}){\tau}_{1,i}\u27ef$.
The proportion of the eggs laid in time period $t$ that have genotype $g$ is calculated as $E}_{g,t}={\displaystyle \frac{{\theta}_{g,t}}{{}_{{\delta}_{t}}}.$
The total (normalized) number of eggs laid in period $t$ is calculated as ${\delta}_{t}=\sum _{g=1}^{9}{\theta}_{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.
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}={\beta}_{0}L$, and for the population post introduction of an intervention, ${A}_{g,t}={S}_{g}L\phantom{\rule{1em}{0ex}}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}={\tau}_{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:
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 (nondeflected) 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}={\displaystyle \sum _{i=4}^{9}{W}_{i,t}}$.
The proportion of the adult female population in time period $t$ that have susceptible, nondeflected 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 nondeflection 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}={\displaystyle \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 nondeflection 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 nondeflection 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:
Let ${[S]}_{d}$ represent the proportion of nondeflection alleles inherited into offspring with genotypes corresponding to susceptible, nondeflected phenotypes:
Let ${[R]}_{d}$ represent the proportion of nondeflection alleles inherited into offspring with genotypes corresponding to resistant, nondeflected phenotypes:
Let ${[DR]}_{d}$ represent the proportion of nondeflection alleles inherited into offspring with genotypes corresponding to resistant, deflected phenotypes:
Let ${[D]}_{d}$ represent the proportion of nondeflection alleles inherited into offspring with genotypes corresponding to nonresistant, deflected phenotypes:
Let ${F}_{S}$ represent the fitness of susceptible, nondeflected phenotypes.
Let ${F}_{R}$ represent the fitness of resistant, nondeflected phenotypes.
Let ${F}_{RD}$ represent the fitness of resistant, deflected phenotypes.
Let ${F}_{D}$ represent the fitness of nonresistant, deflected phenotypes.
Let $\overline{F}}_{D$ represent the average fitness of offspring into which deflection alleles are inherited.
Let ${[DR]}_{D}$ represent the proportion of deflection alleles inherited into offspring with genotypes corresponding to resistant, deflected phenotypes:
Let ${[D]}_{D}$ represent the proportion of deflection alleles inherited into offspring with genotypes corresponding to nonresistant, deflected phenotypes:
This is Equation 1.
Let $\overline{F}}_{d$ represent the average fitness of offspring into which nondeflection alleles are inherited:
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 nondeflection alleles are inherited:
This can be rearranged as follows to give Expression 1.
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:
Let $\overline{F}}_{R$ represent the average fitness of offspring into which resistance alleles are inherited:
The average fitness of offspring into which susceptible alleles are inherited, $\overline{F}}_{r$, is:
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:
This is Expression 2.
References

Loss of household protection from use of insecticidetreated nets against pyrethroidresistant mosquitoes, beninEmerging Infectious Diseases 18:1101–1106.https://doi.org/10.3201/eid1807.120218

Insecticide resistance and vector controlEmerging Infectious Diseases 4:605.https://doi.org/10.3201/eid0404.980410

Status of pyrethroid resistance in Anopheles gambiae sensu latoBulletin of the World Health Organization 77:230–234.

Modifications of pyrethroid effects associated with kdr mutation in Anopheles gambiaeMedical and Veterinary Entomology 14:81–88.https://doi.org/10.1046/j.13652915.2000.00212.x

Restoration of malaria control in the Madagascar highlands by DDT sprayingThe American Journal of Tropical Medicine and Hygiene 66:1.

Laboratory evaluation of the irritancy of bendiocarb, lambdacyhalothrin and DDT to Anopheles gambiaeJournal of the American Mosquito Control Association 9:285–293.

Strain improvement of fungal insecticides for controlling insect pests and vectorborne diseasesCurrent Opinion in Microbiology 15:232–238.https://doi.org/10.1016/j.mib.2011.12.012

Genetic and biological influences in the evolution of insecticide resistanceJournal of Economic Entomology 70:319–323.https://doi.org/10.1093/jee/70.3.319

Role of behavior in the evolution of insect adaptation to insecticides and resistant host plantsBulletin of the Entomological Society of America 30:34–41.https://doi.org/10.1093/besa/30.4.34

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.AJournal of Vector Ecology 25:62–73.

Insecticide resistance in insect vectors of human diseaseAnnual Review of Entomology 45:371–391.https://doi.org/10.1146/annurev.ento.45.1.371

The excitant and repellent effects on mosquitos of sublethal contacts with DDTBulletin of Entomological Research 37:593–607.https://doi.org/10.1017/S0007485300030091

Beyond nature and nurture: phenotypic plasticity in bloodfeeding behavior of Anopheles gambiae s.s. when humans are not readily accessibleAmerican Journal of Tropical Medicine and Hygiene 81:1023–1029.https://doi.org/10.4269/ajtmh.2009.090124

Effect of indoor residual spraying of DDT and bendiocarb on the feeding patterns of Anopheles pseudopunctipennis in MexicoJournal of the American Mosquito Control Association 6:635–640.

Mathematical models of malariaa reviewMalaria Journal 10:202.https://doi.org/10.1186/1475287510202

Eave screening and pushpull tactics to reduce house entry by vectors of malariaAmerican Journal of Tropical Medicine and Hygiene 94:868–878.https://doi.org/10.4269/ajtmh.150632

The Cinderella syndrome: why do malariainfected cells burst at midnight?Trends in Parasitology 29:10–16.https://doi.org/10.1016/j.pt.2012.10.006

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

BookDifficulties in Malaria Eradication Campaigns Due to the Behaviour of the VectorsGeneva: World Health Organization.

Epidemiology of dengue: past, present and future prospectsClinical Epidemiology 5:299–309.https://doi.org/10.2147/CLEP.S34440

DEET microencapsulation: a slowrelease formulation enhancing the residual efficacy of bed nets against malaria vectorsTransactions of the Royal Society of Tropical Medicine and Hygiene 102:259–262.https://doi.org/10.1016/j.trstmh.2007.10.012

Evaluation of synthetic repellents on mosquito nets in experimental huts against insecticideresistant anopheles gambiae and culex quinquefasciatus mosquitoesTransactions of the Royal Society of Tropical Medicine and Hygiene 100:1091–1097.https://doi.org/10.1016/j.trstmh.2006.04.005

Insecticide resistance in disease vectors of public health importancePest Management Science 63:628–633.https://doi.org/10.1002/ps.1406

Fitness costs of disrupting circadian rhythms in malaria parasitesProceedings of the Royal Society B 278:2429–2436.https://doi.org/10.1098/rspb.2010.2457

BookInsecticide Repellency in Malaria Vector Control: A Position PaperVBC Project report No 81131.

A probability model of vector behavior: effects of DDT repellency, irritancy, and toxicity in malaria controlJournal of Vector Ecology 25:48–61.

Behavioral response of Anopheles darlingi to DDTsprayed house walls in AmazoniaBulletin of the Pan American Health Organization 25:210–217.

DDT, global strategies, and a malaria control crisis in South AmericaEmerging Infectious Diseases 3:295–302.https://doi.org/10.3201/eid0303.970305

DDT house spraying and reemerging malariaThe Lancet 356:330–332.https://doi.org/10.1016/S01406736(00)025162

Impact of antiDDT campaigns on malaria controlOutlooks on Pest Management 21:4–11.https://doi.org/10.1564/21feb02

Housing and malaria consensus statement

Evolution of pesticide resistance: interactions between generation time and genetic, ecological, and operational factorsJournal of Economic Entomology 83:1184–1193.https://doi.org/10.1093/jee/83.4.1184

Impact of a spatial repellent on malaria incidence in two villages in Sumba, IndonesiaAmerican Journal of Tropical Medicine and Hygiene 91:1079–1087.https://doi.org/10.4269/ajtmh.130735

Pyrethroid induced behavioral responses of Anopheles dirus, a vector of malaria in ThailandJournal of Vector Ecology 37:187–196.https://doi.org/10.1111/j.19487134.2012.00216.x

Can fungal biopesticides control malaria?Nature Reviews Microbiology 5:377–383.https://doi.org/10.1038/nrmicro1638

Global Technical Strategy for Malaria 2016–2030Global Technical Strategy for Malaria 2016–2030.

Vector control and insecticide resistance

Scientific noteslaboratory evaluation of toxicity of 16 insect repellents in aerosol sprays to adult mosquitoesJournal of the American Mosquito Control Association. 19:271–274.

Indoor residual spray and insecticidetreated bednets for malaria control: theoretical synergisms and antagonismsJournal of the Royal Society Interface 8:799–806.https://doi.org/10.1098/rsif.2010.0537
Decision letter

Ben CooperReviewing Editor; Mahidol Oxford Tropical Medicine Research Unit, Thailand
In the interests of transparency, eLife includes the editorial decision letter and accompanying author responses. A lightly edited version of the letter sent to the authors after peer review is shown, indicating the most substantive concerns; minor comments are not usually included.
Thank you for submitting your article "Evolution leading to sustainable malaria control through evolved spatial repellents" for consideration by eLife. Your article has been favorably evaluated by Diethard Tautz as the Senior Editor and three reviewers: Fred Gould (Reviewer #2), Gerry Kileen (Reviewer #3), and Ben Cooper (Reviewer #1), who is a member of our Board of Reviewing Editors.
The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.
Summary:
This paper proposes a novel strategy for sustained malaria control: starting with a partially effective repellent (effective only against a fraction of the local mosquito population) and using evolution to generate a highly effective spatial repellent. The idea is to combine use of this repellent (referred to as the ESR – the evolved spatial repellent) to stop some proportion of mosquitoes entering human dwellings where a highly effective insecticide (ideally with low contact repellence) is used inside. Those mosquitoes which are repelled, it is argued, will under plausible conditions have selective advantages (because they won't encounter the insecticide) and genes encoding the repellence behaviour will be likely to spread through the mosquito population, leading to an increasingly effective spatial repellent.
All three reviewers found this work interesting and potentially important and felt that, if appropriate revisions are made, it should be published. Reviewer 3 commented that "this piece brings an invaluable idea to the table at a great time, exactly when several new technologies with the right properties to act on these principles are becoming available" but felt more needs to be done to ensure that this work reaches out to the fieldbased community in language they can relate to more easily.
Essential revisions:
1) The work needs to more fully acknowledge related work and to place in historical context. In particular, it was suggested that the authors look over other papers that have built on the Gould 1984 paper by providing more general analysis and by adapting it to specific cases in agriculture.
Some relevant references are:
Kennedy, G. G., F. Gould, O. M. B. de Ponti, and R. E. Stinner. 1987. Ecological, agricultural, genetic, and commercial consideration in the deployment of insectresistant germ plasm. Environ. Entomol. 16: 327338.
CastilloChavez, C., S. A. Levin, and F. Gould. 1988. Physiological and behavioral adaptation to varying environments: a mathematical model. Evolution 42: 986994.
Jongsma, M. A., F. Gould, M. Legros, Y. M. Yang, J. J. A. van Loon, M. Dicke. 2010. Insect oviposition behavior affects the evolution of adaptation to Bt crops: consequences for refuge policies. Evol. Ecol. 24:10171030.
Gatton, M.L., N. Chitnis, T. Churcher, M. J. Donnelly, A. C. Ghani, H. C. J. Godfray, F. Gould, I. Hastings, J. Marshall, H. Ranson, M. Rowland, J. Shaman, S. W. Lindsay. 2013. The Importance of Mosquito Behavioural Adaptations to Malaria Control in Africa. Evolution. 67:12181230.
It was also noted that this work extends the arguments made by others in relation to combinations of irritants and lethal insecticides (PLoS One 9: e95640), to include vaporphase repellents applied inside houses. This preceding work should be overtly and fairly acknowledged in the background and Discussion sections.
2) The manuscript should be revised to increase accessibility and broad appeal, particularly to the fieldbased community. Changes suggested include:
Revising the title to something that is more likely to reach out to a wider audience including nontheoreticians. In particular, the term "evolved repellents" appears to be a term that many potential readers may not be familiar with.
A punchier Conclusions section and a slightly more direct Abstract in language accessible to nontheoreticians would be of use.
Consistent use of common specific wording used widely by epidemiologists and entomologists alike. For example, instead of "infectious bite values" use "entomological inoculation rate (EIR)" (assuming that is what is meant) or the "rate of human exposure to infectious bites".
3) In the second paragraph of the Modelling section: There have been numerous studies demonstrating shifting of mosquitoes to feeding outdoors by indoor protective measures like bed nets, and the reduced sporozoite rates observed (See Lengeler et al. Cochrane review as just one example) cannot be unambiguously attributed to reduced acquisition of infections, and are more likely to arise from reduced survival or feeding frequency (PLoS Med 4: e229, Malaria Journal 20: 207). There is no biological reason to believe individuals would be any less infectious outdoors than indoors, so this is an implausible set of scenarios to explore that should be removed in the interest of brevity. It would be preferable to see the impact of deflection on other important life history parameters like survival, human blood index and feeding frequency explicitly outlined here, as they have been by other preceding models, such as Trans Roy Soc Trop Med 110: 107.
4) The third paragraph of the Introduction describes a view of malaria vector behaviour and transmission system properties based on the peculiar characteristics of a small number (albeit disproportionately important) of very humanspecialized vectors from Africa in particular. See Durnez & Coosemans book chapter on residual malaria transmission, Malaria Journal 13: 330, and the 2014 guidance note on residual malaria transmission by WHO for a more inclusive view of the diversity of malaria transmission. This partial (but substantive) limitation to their relevance needs to be overtly and explicitly outlined in the Discussion, conclusions and Abstract. Also, the fact that this specific subset of transmission scenarios has been considered needs to be explained and justified in the Modelling section. This should be very easy as restricting consideration to the world's most important vectors, which enter a houses a lot, still encompasses a really a big chunk of the world's malaria and is exactly where this strategy is most likely to work.
5) Given the reliance of this piece on the assumption that sufficiently efficacious spatial repellent products exist to address these needs, the Discussion should consider emerging reports of new prototypes with improved product profiles (Parasites and Vectors 8: 322, which contrary to the third paragraph of the Introduction requires no combustion or electricity). Further encouraging updates are likely to be in the literature in the months ahead. The conclusions could apply equally to a single vaporphase active ingredient that is only lethal if the mosquito ignores the repellent effect at sublethal doses for long enough to expose itself to higher lethal doses. The fluorobenzylated pyrethroids seem like an obvious example of compounds already available to which this hypothesis might well apply (US Patent 2015/0289513 A1). It was also suggested that the Discussion could expand on the broader implications by noting that an equally important implication of this hypothesis is that such indoor lethal insecticides could help preserve the efficacy of repellents to be used outdoors where contact insecticides are not feasible to use (Malaria Journal 11: 17 and 13: 146), by also using them indoors.
6) In endemic settings human infectiousness to mosquitoes spans a much narrower range than evaluated here, typically between 2 and 7% (Trans Roy Soc Hyg Trop Med 94: 472, Nature Communications 1: 108, Am J Trop Med Hyg 75 (Suppl. 2): 32 and 38). Surely this is a range exploration exercise that could be removed, to shorten a very long, sometimes repetitive paper. Use of these implausible values may have caused correspondingly implausible estimates for EIR (is this why these have been anonymized by presenting all impact estimates relative to a zero intervention scenario, rather than in absolute terms?). Please present all graphs with commonly used terms like EIR for all outputs, so that their plausibility can be transparently assessed by field specialists.
7) The authors should compare their output to other two locus models that address similar questions with specified parameters, matching the parameters and see if the results are the same. This will help to confirm that there have not been errors in implementation.
[Editors' note: further revisions were requested prior to acceptance, as described below.]
Thank you for resubmitting your work entitled "Using evolution to generate sustainable malaria control with spatial repellents" for further consideration at eLife. Your revised article has been favorably evaluated by Diethard Tautz as the Senior editor and a Reviewing editor.
The manuscript has been improved but there are some remaining issues that need to be addressed before acceptance, as outlined below:
1) In the second paragraph of the Modelling section it should say "see Appendix 1".
2) Appendix 2 makes it clear that [dr], [dR] etc. refer to proportions of alleles in gametes. Shouldn't the Modelling section also make this clear?
3) Near the end of the Modelling section it contains a couple of "it can be seen" statements that seem either incomplete, imprecise or require further justification. In the first it claims that spread of deflection will by favoured by maximising the "fitness difference between susceptible and deflected phenotypes", but this won't be true if F_{D} < F_{S}. Also not clear why there is no reference here to the F_{R} – F_{D} term. The second "it can be seen" is even less clear as there are terms on the rhs of the inequality that both increase and decrease with increasing frequency of resistance alleles. I'm prepared to believe it's generally true with plausible parameters, but I don't think it's necessarily easily seen.
4) In the fourth paragraph of the Modelling section "data become" rather than "data becomes" (also some slightly erratic punctuation throughout which no doubt will be fixed in the production process).
5) In the fifth paragraph of the Results "0.05%" – the Figure 2 caption says 0.5% – presumably the latter is correct.
https://doi.org/10.7554/eLife.15416.011Author response
Essential revisions:
1) The work needs to more fully acknowledge related work and to place in historical context. In particular, it was suggested that the authors look over other papers that have built on the Gould 1984 paper by providing more general analysis and by adapting it to specific cases in agriculture.
Some relevant references are:
Kennedy, G. G., F. Gould, O. M. B. de Ponti, and R. E. Stinner. 1987. Ecological, agricultural, genetic, and commercial consideration in the deployment of insectresistant germ plasm. Environ. Entomol. 16: 327338.
CastilloChavez, C., S. A. Levin, and F. Gould. 1988. Physiological and behavioral adaptation to varying environments: a mathematical model. Evolution 42: 986994.
Jongsma, M. A., F. Gould, M. Legros, Y. M. Yang, J. J. A. van Loon, M. Dicke. 2010. Insect oviposition behavior affects the evolution of adaptation to Bt crops: consequences for refuge policies. Evol. Ecol. 24:10171030.
Gatton, M.L., N. Chitnis, T. Churcher, M. J. Donnelly, A. C. Ghani, H. C. J. Godfray, F. Gould, I. Hastings, J. Marshall, H. Ranson, M. Rowland, J. Shaman, S. W. Lindsay. 2013. The Importance of Mosquito Behavioural Adaptations to Malaria Control in Africa. Evolution. 67:12181230.
It was also noted that this work extends the arguments made by others in relation to combinations of irritants and lethal insecticides (PLoS One 9: e95640), to include vaporphase repellents applied inside houses. This preceding work should be overtly and fairly acknowledged in the background and Discussion sections.
We thank the reviewers for bringing the above papers to our attention, and have included citations for them in the manuscript. We have also included citations for an illustrative range of previous modelling in relevant areas, and cited a review of mathematical modelling for malaria,[1] to provide additional context for interested readers.
2) The manuscript should be revised to increase accessibility and broad appeal, particularly to the fieldbased community. Changes suggested include:
Revising the title to something that is more likely to reach out to a wider audience including nontheoreticians. In particular, the term "evolved repellents" appears to be a term that many potential readers may not be familiar with.
A punchier Conclusions section and a slightly more direct Abstract in language accessible to nontheoreticians would be of use.
We have revised the title to eliminate the unfamiliar term ‘Evolved repellents’, as suggested. We have also made revisions to the Abstract and changes throughout the text in response to this and other comments which have hopefully helped to clarify and simplify some key sections of the manuscript.
Consistent use of common specific wording used widely by epidemiologists and entomologists alike. For example, instead of "infectious bite values" use "entomological inoculation rate (EIR)" (assuming that is what is meant) or the "rate of human exposure to infectious bites".
We take the reviewer’s point that our message needs to be in the most accessible form possible, and have made changes to the text to try to achieve this. If we make the simple assumption that human population size is not materially affected by the chosen intervention, and is stable over time, we can equate our infectious bite values to EIR and we have now explained this in the manuscript.
3) In the second paragraph of the Modelling section: There have been numerous studies demonstrating shifting of mosquitoes to feeding outdoors by indoor protective measures like bed nets, and the reduced sporozoite rates observed (See Lengeler et al. Cochrane review as just one example) cannot be unambiguously attributed to reduced acquisition of infections, and are more likely to arise from reduced survival or feeding frequency (PLoS Med 4: e229, Malaria Journal 20: 207).
We have amended our wording about the Ndiath result to ensure that we do not misrepresent its significance.
“The impact of deflection on malaria prevalence is determined by the proportion of mosquitoes deflected by a repellent and the probability, compared to nondeflected mosquitoes, that they will then acquire and transmit a Plasmodium infection. […] We have not explicitly assessed the additional effects of reduced transmission from outdoorfeeding infectious mosquitoes resulting from, for example, transfer to feeding on nonhuman hosts, but any such effects would clearly serve to further enhance the public health benefits of deflection.”
There is no biological reason to believe individuals would be any less infectious outdoors than indoors, so this is an implausible set of scenarios to explore that should be removed in the interest of brevity. It would be preferable to see the impact of deflection on other important life history parameters like survival, human blood index and feeding frequency explicitly outlined here, as they have been by other preceding models, such as Trans Roy Soc Trop Med 110: 107.
With apologies, this and other comments from the reviewers demonstrate that our explanation of this part of our analysis was confusing, and we have amended the manuscript to improve this. The impact of deflection on important life history parameters such as human blood index is included in the ‘transmission probability’ parameter in our model, which is the probability of transmission to the vector per feed and encompasses any difference between the probability of acquiring an infection when feeding outside rather than inside a human dwelling. As well as different likelihoods of taking a nonhuman feed, other possible causes of differential per feed transmission probabilities are not yet wellexplored, but there is for example, interesting evidence to suggest that vector feeding times may affect the availability or transmissibility of gametocytes, and we have added some additional references regarding this. It is because of the uncertainty regarding such effects, and the potential for them to vary between species (vector and parasite) and context, that we have chosen to consider a wide range of possible transmission probabilities in our modelling.
Assumed differential vector survival with and without deflection is intrinsic to our modelling, and its effects on survival to feed, survival to reach infectiousness and survival to give each infectious bite are reflected directly through the survival parameters in the model rather than the per feed probabilities of transmission to/from the vector. We have amended the text to ensure this is clear.
4) The third paragraph of the Introduction describes a view of malaria vector behaviour and transmission system properties based on the peculiar characteristics of a small number (albeit disproportionately important) of very humanspecialized vectors from Africa in particular. See Durnez & Coosemans book chapter on residual malaria transmission, Malaria Journal 13: 330, and the 2014 guidance note on residual malaria transmission by WHO for a more inclusive view of the diversity of malaria transmission. This partial (but substantive) limitation to their relevance needs to be overtly and explicitly outlined in the Discussion, conclusions and Abstract. Also, the fact that this specific subset of transmission scenarios has been considered needs to be explained and justified in the Modelling section. This should be very easy as restricting consideration to the world's most important vectors, which enter a houses a lot, still encompasses a really a big chunk of the world's malaria and is exactly where this strategy is most likely to work.
We have amended the text to explicitly identify that ESR is a tool with relevance in the context of public health campaigns based on indoor residual spraying (IRS), which inherently target only indoorfeeding vectors. We thank the reviewer for highlighting that this was not sufficiently clear in the manuscript.
“Malaria vectors typically feed indoors between dusk and dawn and vector control has focused on exploiting this behaviour to deliver lethal control measures against indoorfeeding mosquitoes [Yakob, Dunning and Yan, 2011], 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, Hamon and World Health Organization, 1963].”
The reviewer identifies targeting of indoorfeeding vectors as exclusively affecting highly anthropophilic vectors, however, we would argue that other important vectors which have a high propensity to feed on nonhuman hosts, such as An. arabiensis, also show indoor feeding on human hosts [2] and are targeted by IRS. These are potentially very important targets for ESR, since the fitness costs of deflection would be expected to be lower for a vector more likely to deflect to an alternative host, making it easier to establish and maintain an ESR, whilst increasing the proportion of feeds taken on nonhuman hosts also reduces the probability of transmission to/from vectors, since the vectors can neither acquire nor transmit a Plasmodium infection during a feeding cycle utilising a nonhuman host.
5) Given the reliance of this piece on the assumption that sufficiently efficacious spatial repellent products exist to address these needs, the Discussion should consider emerging reports of new prototypes with improved product profiles (Parasites and Vectors 8: 322, which contrary to the third paragraph of the Introduction requires no combustion or electricity). Further encouraging updates are likely to be in the literature in the months ahead.
Thanks to the reviewer for bringing to our attention to this interesting and positive article by Govella et al. [3], which assesses the performance of close proximity transfluthrinimpregnated hessian strips in providing protection from mosquito bites outdoors. As discussed by Govella et al., there are still a number of issues to consider before translating this into a method for keeping mosquitoes out of properties (and see below), but we agree that this is nonetheless interesting and relevant to a consideration of currently available repellents. Transfluthrin has also been used for short periods in another study as a spatial repellent (in the sense used in our manuscript) [4]. We have therefore added references to this and some other recent articles on research into malaria vector repellents with various modes of repellent action.
“Effective action at a distance is still challenging to achieve, and commonly requires active dispersal through combustion or powered devices [Menger et al., 2016; Syafruddin et al., 2014].”
The conclusions could apply equally to a single vaporphase active ingredient that is only lethal if the mosquito ignores the repellent effect at sublethal doses for long enough to expose itself to higher lethal doses. The fluorobenzylated pyrethroids seem like an obvious example of compounds already available to which this hypothesis might well apply (US Patent 2015/0289513 A1).
Thanks to the reviewer for this interesting comment. Review of specific compounds with potential for use as ESRs is clearly a key part of the next stage of this research, but is beyond the intended scope of the current paper.
It was also suggested that the Discussion could expand on the broader implications by noting that an equally important implication of this hypothesis is that such indoor lethal insecticides could help preserve the efficacy of repellents to be used outdoors where contact insecticides are not feasible to use (Malaria Journal 11: 17 and 13: 146), by also using them indoors.
We think this is a very interesting idea. The context and the specific characteristics of the chosen ESR would be critical to whether this would work. For example, naturally outdoorfeeding vectors will experience little or no selection from indoor insecticides, so ESR will have little influence on their behaviour. If a vector population is feeding both indoors and outdoors, then widespread use of an ESR outdoors would increase the relative fitness costs of deflection so might delay establishment of an ESR and/or accelerate its loss to the spread of insecticide resistance. If the overall fitness cost is low though, say a zoophagic vector with plenty of available nonhuman hosts, or use of outdoor repellent for personal protection by only a small number of atrisk people, then the ESR concept might indeed generate and/or maintain functionality of a compound as an outdoor repellent. This would however also be conditional on finding a compound which has the required characteristics both to protect a contained indoor space and to protect an open outdoor area and/or be used on skin or clothing. We have tried to include a suitably caveated but succinct reference to this idea in the manuscript.
6) In endemic settings human infectiousness to mosquitoes spans a much narrower range than evaluated here, typically between 2 and 7% (Trans Roy Soc Hyg Trop Med 94: 472, Nature Communications 1: 108, Am J Trop Med Hyg 75 (Suppl. 2): 32 and 38). Surely this is a range exploration exercise that could be removed, to shorten a very long, sometimes repetitive paper. Use of these implausible values may have caused correspondingly implausible estimates for EIR (is this why these have been anonymized by presenting all impact estimates relative to a zero intervention scenario, rather than in absolute terms?).
We particularly appreciate this comment as it has brought to our attention that we have at some point lost from the manuscript the values assumed for some base parameter values, including the assumed probability of a vector acquiring a Plasmodium infection when feeding indoors, which is 4%, well placed in the 2%−7% range indicated by the reviewer. With apologies for this omission, we have reinstated this and the other base parameter values in the document.
The equivalent vector infection probabilities assumed for deflected vectors feeding outdoors are 100% x 4% = 4%, 50% x 4% = 2%, 25% x 4% = 1% and 0% x 4% = 0%. Because our results are formulated in comparative rather than absolute terms, the relationship between the assumed indoor and outdoor transmission values is key. However, the way we have given this information has clearly generated confusion, so we have amended both the main text and the Figure 4 labels to make clearer the actual rather than just the relative transmission probabilities assumed.
We agree that our presentation of this issue involved some unnecessary repetition and have amended the manuscript to improve this.
As indicated in the manuscript, we formulate our results in comparative rather than absolute terms to minimize the impact of poorly quantified or location/species specific parameters on the results. The intention is to consider whether or not the proposed innovation can offer a reduction in transmission, and for this purpose comparative metrics are highly suitable, and have previously been used for evaluation of novel antimalaria interventions and assessment of the significance of possible vectorparasite interactions [57]. Our assumption here is that this analysis considers the potential benefits of adding ESR to existing public health strategies, and that detailed, sitespecific analysis would be carried out to establish its suitability for each given setting.
Please present all graphs with commonly used terms like EIR for all outputs, so that their plausibility can be transparently assessed by field specialists.
For Figure 1 the graphs show the proportion of the vector population with given genotypes and phenotypes, plotted against time. We have changed the headings slightly to ensure that this is clear. We have also corrected the xaxis labels for this figure, which had mysteriously acquired unwanted ‘%’ signs during the last stage of figure preparation.
For Figures 2 and 3 the grid plots indicate the combinations of deflected and resistant percycle survival values which generate a vector population with at least 80% deflected phenotypes after 300 time periods. We have amended the figure captions to try to make the meaning of the figures as clear as possible, but we are not aware of any standard terminology that could be used here.
Figure 4 shows the number of infectious bites from the vector population per unit of time as a percentage of that without any intervention. With an additional assumption that the human population size does not change over time this metric equates to EIR expressed as a% of the EIR with no intervention. We have stated this additional assumption and changed the graph labels accordingly.
7) The authors should compare their output to other two locus models that address similar questions with specified parameters, matching the parameters and see if the results are the same. This will help to confirm that there have not been errors in implementation.
With the agreement of the editors we have addressed this point using a separate ‘proof of concept’ (‘POC’) model, which we created prior to building the complete model documented in Appendix 1. Although both enacted in Excel, the two models were built separately from scratch, so would not be expected to include matching typos or other corresponding errors in implementation. The POC model does not reflect the mating of females only at the start of their adult lives, and does not track Plasmodium infection or transmission in the vector population. We compared results from the two models for the three cases illustrated in Figure 1 and for one COR case from Figure 3. The change in genotype proportions was slower using the PM than the POC, consistent with the damping effect of reflecting the male genotypes applicable at the time of initial mating for females surviving from period to period. Otherwise the results were qualitatively consistent in all cases.
References
1) Sandip Mandal RRS, Somdatta Sinha. Mathematical models of malaria  a review. Malaria Journal. 2011;10(202).
2) Sinka ME, Bangs MJ, Manguin S, Coetzee M, Mbogo CM, Hemingway J, et al. The dominant Anopheles vectors of human malaria in Africa, Europe and the Middle East: occurrence data, distribution maps and bionomic précis. Parasites & Vectors. 2010;3(1):134.
3) Govella NJ, Ogoma SB, Paliga J, Chaki PP, Killeen G. 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. 2015;8(1):19.
4) Wagman JM, Grieco JP, Bautista K, Polanco J, Briceño I, King R, et al. The field evaluation of a pushpull system to control malaria vectors in Northern Belize, Central America. Malaria Journal. 2015;14(1):111.
5) Lynch P, Grimm U, Thomas M, Read A. Prospective malaria control using entomopathogenic fungi: comparative evaluation of impact on transmission and selection for resistance. Malaria Journal. 2012;11(1):383.
6) Read AF, Lynch PA, Thomas MB. How to Make EvolutionProof Insecticides for Malaria Control. PLoS Biology. 2009;7(4):e58.
7) Cator L, Lynch P, Thomas M, Read A. Alterations in mosquito behaviour by malaria parasites: potential impact on force of infection. Malaria Journal. 2014;13(1):164.
[Editors' note: further revisions were requested prior to acceptance, as described below.]
1) In the second paragraph of the Modelling section it should say "see Appendix 1".
Corrected.
2) Appendix 2 makes it clear that [dr], [dR] etc. refer to proportions of alleles in gametes. Shouldn't the Modelling section also make this clear?
Main text amended, now reads “With $\left[dr\right]$, $\left[dR\right]$, $\left[Dr\right]$, $\left[DR\right]$, $\left[d\u2013\right]$ and $\left[D\u2013\right]$ representing, in the zygote genotypes for the population at a given time point, the proportion of alleles at the deflection locus which are nondeflection alleles paired with susceptible alleles, nondeflection alleles paired with resistant alleles, deflection alleles paired with susceptible alleles, deflection alleles paired with resistant alleles, nondeflection 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 newlyemerged females.”
3) Near the end of the Modelling section it contains a couple of "it can be seen" statements that seem either incomplete, imprecise or require further justification. In the first it claims that spread of deflection will by favoured by maximising the "fitness difference between susceptible and deflected phenotypes", but this won't be true if F_{D} < F_{S}.
The ESR concept inherently depends on deflected phenotypes having higher fitness than susceptible (not deflected not resistant) phenotypes, i.e. the fitness cost of entering properties and being killed in properties treated with insecticide must be higher than the fitness costs of being deflected away from properties treated with repellent. We had not stated this clearly in the text, and have now amended the text to correct this. (Modelling section, third paragraph from the end).
Also not clear why there is no reference here to the F_{R} – F_{D} term.
We have amended the text to explicitly state that maximising the fitness of deflected relative to resistant phenotypes will support establishment of an ESR.
The second "it can be seen" is even less clear as there are terms on the rhs of the inequality that both increase and decrease with increasing frequency of resistance alleles. I'm prepared to believe it's generally true with plausible parameters, but I don't think it's necessarily easily seen.
We have reviewed the ‘it can be seen’ statements and made amendments to the text to try to make things clearer and more precise. We have relocated an explanation linking the recommended application strategy to the fitness of deflected resistant phenotypes to a later point in the text where this has a clear and simple impact on the spread of resistance if deflection reaches fixation.
4) In the fourth paragraph of the Modelling section "data become" rather than "data becomes" (also some slightly erratic punctuation throughout which no doubt will be fixed in the production process).
Corrected.
5) In the fifth paragraph of the Results "0.05%" – the Figure 2 caption says 0.5% – presumably the latter is correct.
Corrected, with thanks to the reviewer for spotting this.
https://doi.org/10.7554/eLife.15416.012Article and author information
Author details
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.
Reviewing Editor
 Ben Cooper, Mahidol Oxford Tropical Medicine Research Unit, Thailand
Publication history
 Received: February 21, 2016
 Accepted: September 13, 2016
 Version of Record published: October 25, 2016 (version 1)
 Version of Record updated: October 27, 2016 (version 2)
Copyright
© 2016, Lynch et al.
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

 1,911
 Page views

 321
 Downloads

 8
 Citations
Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.