Using evolution to generate sustainable malaria control with spatial repellents
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
Article 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.
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

 2,088
 views

 343
 downloads

 15
 citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading

 Epidemiology and Global Health
Artificially sweetened beverages containing noncaloric monosaccharides were suggested as healthier alternatives to sugarsweetened beverages. Nevertheless, the potential detrimental effects of these noncaloric monosaccharides on blood vessel function remain inadequately understood. We have established a zebrafish model that exhibits significant excessive angiogenesis induced by high glucose, resembling the hyperangiogenic characteristics observed in proliferative diabetic retinopathy (PDR). Utilizing this model, we observed that glucose and noncaloric monosaccharides could induce excessive formation of blood vessels, especially intersegmental vessels (ISVs). The excessively branched vessels were observed to be formed by ectopic activation of quiescent endothelial cells (ECs) into tip cells. Singlecell transcriptomic sequencing analysis of the ECs in the embryos exposed to high glucose revealed an augmented ratio of capillary ECs, proliferating ECs, and a series of upregulated proangiogenic genes. Further analysis and experiments validated that reduced foxo1a mediated the excessive angiogenesis induced by monosaccharides via upregulating the expression of marcksl1a. This study has provided new evidence showing the negative effects of noncaloric monosaccharides on the vascular system and the underlying mechanisms.

 Epidemiology and Global Health
 Microbiology and Infectious Disease
Influenza viruses continually evolve new antigenic variants, through mutations in epitopes of their major surface proteins, hemagglutinin (HA) and neuraminidase (NA). Antigenic drift potentiates the reinfection of previously infected individuals, but the contribution of this process to variability in annual epidemics is not well understood. Here, we link influenza A(H3N2) virus evolution to regional epidemic dynamics in the United States during 1997—2019. We integrate phenotypic measures of HA antigenic drift and sequencebased measures of HA and NA fitness to infer antigenic and genetic distances between viruses circulating in successive seasons. We estimate the magnitude, severity, timing, transmission rate, agespecific patterns, and subtype dominance of each regional outbreak and find that genetic distance based on broad sets of epitope sites is the strongest evolutionary predictor of A(H3N2) virus epidemiology. Increased HA and NA epitope distance between seasons correlates with larger, more intense epidemics, higher transmission, greater A(H3N2) subtype dominance, and a greater proportion of cases in adults relative to children, consistent with increased population susceptibility. Based on random forest models, A(H1N1) incidence impacts A(H3N2) epidemics to a greater extent than viral evolution, suggesting that subtype interference is a major driver of influenza A virus infection ynamics, presumably via heterosubtypic crossimmunity.