Chemical insecticides against adult mosquitoes are a key element in most malaria management programmes, but their efficacy is threatened by the evolution of insecticide-resistant mosquitoes. By killing only older mosquitoes, entomopathogenic fungi can in principle significantly impact parasite transmission while imposing much less selection for resistance. Here an assessment is made as to which of the wide range of possible virulence characteristics for fungal biopesticides best realise this potential.
With mathematical models that capture relevant timings and survival probabilities within successive feeding cycles, transmission and resistance-management metrics are used to compare susceptible and resistant mosquitoes exposed to no intervention, to conventional instant-kill interventions, and to delayed-action biopesticides with a wide range of virulence characteristics.
Fungal biopesticides that generate high rates of mortality at around the time mosquitoes first become able to transmit the malaria parasite offer potential for large reductions in transmission while imposing low fitness costs. The best combinations of control and resistance management are generally accessed at high levels of coverage. Strains which have high virulence in malaria-infected mosquitoes but lower virulence in malaria-free mosquitoes offer the ultimate benefit in terms of minimizing selection pressure whilst maximizing impact on transmission. Exploiting this phenotype should be a target for product development. For indoor residual spray programmes, biopesticides may offer substantial advantages over the widely used pyrethroid-based insecticides. Not only do fungal biopesticides provide substantial resistance management gains in the long term, they may also provide greater reductions in transmission before resistance has evolved. This is because fungal spores do not have contact irritancy, reducing the chances that a blood-fed mosquito can survive an encounter and thus live long enough to transmit malaria.
Delayed-action products, such as fungal biopesticides, have the potential to achieve reductions in transmission comparable with those achieved with existing instant-kill insecticides, and to sustain this control for substantially longer once resistant alleles arise. Given the current insecticide resistance crisis, efforts should continue to fully explore the operational feasibility of this alternative approach.
The impressive reductions in global malaria burden achieved this century by chemical insecticides against adult mosquitoes could be eroded by insecticide-resistant mosquitoes [1–6], just as they were last century . In principle, the evolution of insecticide resistance could be considerably slowed and perhaps prevented altogether by vector control aimed at killing only older mosquitoes, so-called late-life action (LLA) . Malaria parasites in a mosquito host take at least nine days to develop to a stage which can be transmitted to a human via an infectious bite . Since mortality in wild mosquito populations is high, the majority of eggs are produced by young mosquitoes. Thus, a vector-control treatment which kills only older mosquitoes could remove infected mosquitoes before they can transmit malaria whilst only impacting the reproductive success of only the relatively few mosquitoes that survive to old age. This would dramatically reduce transmission while exerting only weak selection for resistance.
One option for an LLA vector-control measure is entomopathogenic fungi . Naturally occurring strains of two fungi, Beauveria bassiana and Metarhizium anisopliae, are already in commercial use for agricultural applications and have been shown to infect and kill mosquitoes in laboratory and field settings. Fungal spores can be picked up by mosquitoes following contact with treated surfaces, and so could be used against mosquitoes in indoor residual spray (IRS) programmes, or delivered via traps, curtains or netting [11–16].
A wide variety of mortality schedules can be induced in Anopheles by entomopathogenic fungi . In some cases, all mosquitoes can be killed within a few days; in others, background mortality rates can be barely altered. This virulence variation depends on isolate , dose  and malaria-infection status [15, 18], see also . Lethality can also be increased by genetically modifying fungal isolates [20–22].
If fungal entomopathogens are to realize the potential of the LLA approach to sustainable malaria control, candidate biopesticides need to be chosen which balance reductions in parasite transmission (maximized by high fungal virulence) with resistance management (maximized by low fungal virulence). Here a mathematical model is used to ask which virulence phenotypes best achieve this balance. The intention is to guide the development of target product profiles. The possible efficacy of fungal biopesticides in IRS campaigns is compared with that of pyrethroid-based insecticides now in widespread use. Pyrethroids are highly lethal if contacted by a mosquito, but they also have a strong excito-repellency effect, which can drive away mosquitoes before they receive a lethal dose [23–25]. There is evidence that fungal spores do not repel mosquitoes , raising the prospect that, for IRS, fungal biopesticides might more effectively reduce transmission than pyrethroid-based technologies currently in use.
Many malaria transmission models already exist , but most do not capture the detailed timings and probabilities of infection, infectiousness, reproduction and mortality over the mosquito lifespan which are key to assessing whether LLAs can provide a useful balance of transmission control and low selection for resistance. In order to encompass these elements, a model has been developed with two separate components, a markovian, deterministic, feeding cycle model (FCM) which calculates survival, egg-laying and infectious bite values during the lifetime of an adult mosquito, and a population model (PM) which tracks the population-level spread of resistance alleles and corresponding loss of transmission control. The model is a development of a simpler version previously used to evaluate putative chemical LLAs . Other modelling frameworks used to assess the LLA approach are heuristically useful but lack sufficient detail to define target virulence schedules [28–30].
The feeding cycle model
The FCM calculates survival, egg-laying and infectious bite values across a series of discrete adult age classes for a specified type of mosquito (e.g., susceptible) subjected to a given intervention (e.g., a particular fungal biopesticide at a particular coverage). Each sequential age class is defined as lasting for the average length of one gonotrophic cycle. Use of the mosquito feeding cycle as the basis for age-structured analyses of mosquito populations is well established [31–34].
The FCM tracks possible states and transitions through each age class (i), applying survival, exposure and infection probabilities (Figure 1). Infection status for a biopesticide (l) or malaria (m), is zero for no infection, otherwise equal to the age of the infection. State changes depend on the preceding state, the passage of time, mortality rates and the probabilities of certain events, such as contacting a biopesticide when resting after a human blood meal. For example, for a case analysing the effects of a fungal biopesticide, a mosquito commencing its fourth cycle with an infectious, three-cycle-old malaria infection, and no fungal infection, will spend a defined period of time searching for a host, with an associated probability of dying from background mortality while it does so. It will then attack a host, with a given probability that the selected host will be a non-infectious human, a malaria-infectious human, or non-human, and a given probability of being killed whilst attacking the host before biting. If it survives to bite, and if the host is human, this is recorded as an infectious bite. There is then a given probability that it is killed by the host after biting. If it is not killed, it begins a period of resting, during which, if the chosen host was human, it has a fixed probability of encountering and being infected by the fungus, as well as a given probability of dying from background mortality before leaving to search for an egg-laying site. During the search for an egg-laying site there is a given probability that the mosquito may die from background mortality, or from the effects of its newly acquired fungal infection. If it survives searching for an egg-laying site it may die before or after laying with given probabilities. If still alive at the end of the cycle, it begins its fifth cycle with an infectious, four-cycle-old malaria infection, and a one-cycle-old fungal infection. For a case analysing the effects of a conventional instant-kill chemical pesticide, the analysis would include a probability of contacting the pesticide after biting the host, and a probability of death, assumed to be instant, resulting from that contact, with zero probability of contacting a biopesticide.
Both conventional instant-kill and delayed-action biopesticides offer public health benefits by reducing the numbers of mosquitoes that survive to give infectious bites in a treated population. Clearly the extent to which a reduction in infectious bites maps to reduced transmission and reduced numbers or severity of malaria cases in a human host population involves many complex, context-specific factors. For comparative purposes, however, it is assumed that in a given context, a given reduction in infectious bites will generate the same reduction in malaria transmission and malaria morbidity and mortality irrespective of the type of intervention from which it results. For generality, therefore, the comparative public health benefits of the insecticides considered in this analysis are all evaluated based on the reduction in infectious bites which they provide. This is quantified in the FCM for a given phenotype by calculating RAIB, the proportionate reduction in the average number of infectious bites per mosquito per lifetime (AIB), defined as
Assuming that the rate at which newly maturing adults join a population is constant through time, and that the size of the human host population is unaffected by the intervention being assessed, RAIB is equal to the proportionate reduction in the entomological inoculation rate (EIR), the number of infectious bites experienced per person per unit of time.
To evaluate mosquito fitness, the average number of eggs produced per mosquito per lifetime is used as a proxy for lifetime reproductive success (LRS). The selection coefficient, the proportionate fitness benefit of resistance to a given intervention, is calculated as . A selection coefficient of zero means no selection pressure in favour of resistance, with higher selection coefficients indicating increasingly strong selection for resistance.
Formulating these key variables in relative terms minimizes the sensitivity of the conclusions to parameter values that are independent of the vector-control treatment or mosquito phenotype being evaluated.
The primary definitions of the FCM are given in Table 1, and its main features are detailed below. A detailed derivation of the FCM is given in Additional file 1: Appendix A. Baseline parameter values used in the analysis are summarized in Table 2.
Variables and parameters for the feeding cycle model
Variable or Parameter
Comments and Constraints
Time, measured in whole units equal to average length of sporogonic cycle, from infection of mosquito by malaria to cycle from which mosquito gives infectious bites
input 0< D
Number of age classes included in analysis
Cycle number (identifies specific cycle in the λ cycles over which probabilities are tracked in the FCM)
0≤ i ≤λ
Malaria status, the number of whole or partial cycles since infection with malaria
0≤ m ≤λ, m = 0 means not infected
Biopesticide infection status, the number of whole or partial cycles since infection with biopesticide
0≤ l≤ λ, l = 0 means not infected
Average number of eggs laid in cycle i by mosquitoes surviving to the start of cycle i
Average lifetime number of eggs laid per mosquito
Average number of eggs laid in cycle i, by mosquitoes starting cycle i with malaria status m and biopesticide status l
Average probability of survival from start of cycle i to start of cycle i+1
Average probability that a mosquito starting cycle i with malaria status m and biopesticide status l, will survive to start of cycle i+1
Average probability of a mosquito being alive at start of cycle i.
Average probability of a mosquito being alive, with malaria status m and biopesticide status l, at start of period i.
Probability that a mosquito alive at start of cycle i with malaria status m and biopesticide status l, survives and bites host type h in cycle i
Type of host attacked
h=2, non-infectious human
h=3, infectious human
Average number of infectious bites in cycle i per mosquito alive at the start of cycle i
Average lifetime number of infectious bites per mosquito
Values used in FCM for this analysis
Variable or Parameter
Background instantaneous mortality rate for mosquito age i
Length of gonotrophic cycle
Time spent host searching and feeding during a cycle
Time spent finding oviposition site and laying during a cycle
Length of resting period (days)
Proportion human population infectious for malaria4
Probability attacks non-human host
Probability killed when attacking host before biting
Probability killed when attacking host after biting (excluding mortality from insecticide treatments)
Probability becomes infected with malaria when biting infectious human host4
Number of eggs laid per successfully laying mosquito per cycle
Time, measured in whole units equal to length of gonotrophic cycle, from infection of mosquito to cycle from which mosquito gives infectious bites
33 Based on 10.78 1 days
Baseline probability that mosquito contacts and is killed by conventional instant-kill chemical insecticide (CC) whilst resting after biting human host
0for cases not assessing use of CC
0.8for cases assessing use of CC
Baseline probability that mosquito contacts and is affected by delayed action pesticide whilst resting after biting human host
0for cases not assessing use of delayed action pesticide
0.8for cases assessing use of delayed action pesticide
Number of age classes included in analysis
1.Averages taken from four geographic locations . Results using individual geographic data sets are expected to give qualitatively equivalent results. Limited sensitivity analysis was consistent with this assumption, so use of the average figures was considered adequate for the present analysis.
2.Since we are only interested in comparative values, the absolute value for the number of eggs per lay is immaterial, 100 has been used as a convenient normalised value.
3.The number of cycles assumed for sporogonic development is calculated from the average number of days for sporogonic development and the average number of days per gonotrophic cycle, rounded down to give a whole number of cycles. This is a conservative assumption with respect to the amount of EIR reduction calculated for given fungal virulence parameters.
4.The data set used provides a total probability of acquiring a malaria infection when biting a human host. This has been used as the value for parameter p, with M=1.00, to give the appropriate combined probability, Mp.
5. Assumes c.11.1% of every cycle is spent resting (8 hours in a 72 hour cycle), with the rest of the gonotrophic cycle divided equally between laying and feeding.
6. Estimated 10% mortality per feeding attempt, divided equally between pre- and post-bite.
The probability that a mosquito contacts and is affected (killed or infected) by a conventional or biological insecticide after biting a human host is input as a single ‘coverage’ value, incorporating the probabilities of being in a treated property, of contacting the pesticide, and of being affected by the pesticide during contact. It is assumed that physical constraints on the proportion of surfaces and internal areas treated will apply equally to conventional and fungal insecticides, and that for mosquitoes contacting treated surfaces, biopesticides can potentially offer rates of infection equivalent to the rates of mortality generated by conventional insecticides. The latter assumption is supported by field trials showing >86% infection of mosquitoes entering outdoor bait boxes , 76% infection in experimental huts with fungus-impregnated eave curtains , and laboratory trials showing >95% infection from treated clay pots  or exposure to treated clay tiles .
The average number of eggs laid in a given cycle, by mosquitoes surviving to the start of that cycle, Fi, is calculated as
This provides the basis for the evaluation of relative fitness using a comparison of values for φ, lifetime egg production, representing LRS, .
Comparative levels of transmission control are assessed using u, the average number of infectious bites per mosquito lifetime, .
The average number of infectious bites during cycle i per mosquito surviving to the beginning of cycle i, Ii, is calculated as
i > D
The average probability of survival from start of cycle i to start of cycle (i + 1) is Si, with
The population model
The PM tracks susceptible and resistant phenotypes over a sequence of time periods for a population subject to a given vector-control treatment. The key outputs, calculated for each time period, are the proportion of the population with resistant and susceptible phenotypes and the overall reduction in infectious bites across the population compared to a susceptible population with no vector-control treatment.
The variables and parameters for the PM are described in Table 3, and baseline values used in the analysis are summarized in Table 4.
Variables and parameters for the population model
Variable or Parameter
Comments & Constraints
Period number (periods over which the population is tracked)*
Dominance of resistance allele
dominant d = 1
recessive d = 0
Genotype (normal allele s, resistant allele r)
(s,s) g = 1
(s,r) g = 2
(r,r) g = 3
Proportion of total population having genotype g at start of period n
Proportion of the population resistant at start of period n
Average number of infectious bites per mosquito in population in period n
Size of initial population (susceptibles in the presence of treatment) as proportion of base population (susceptibles without treatment)
value from FCM
Population size in period n as proportion of initial population size
Average infectious bites during one time period from an untreated population
value from FCM
Number of infectious bites from treated population during time period n, expressed as a % of the number of infectious bites during one time period from a susceptible population without treatment,
Chosen measure of control
Number of periods between egg-laying and adult emergence
*As for the FCM, the duration of one gonotrophic cycle is used as a unit of time. For convenience we use ‘cycles’ to refer to mosquito age and ‘periods’ to refer to the sequential time periods for which values are calculated in the PM.
Values used in the population model for this analysis
Variable or Parameter
Proportion of total population having genotype g at start of period 1
G3,1 = 0
Dominance of resistance allele (0=recessive, 1=dominant)
d = 1
Number of periods between egg-laying and adult emergence
Fitness factor for males with genotype g
All other input values use results calculated by the FCM.
A detailed derivation of the model is given in Additional file 2: Appendix B. In brief, the PM works in discrete time periods, each equivalent to the length of one gonotrophic cycle, with recruitment of newly emerged adult mosquitoes treated as occurring at the start of each time period. For each sequential time period, the proportion of the population comprised by each genotype in each age class is calculated, reflecting the genotypes of new adult recruits and the survival of adults in each age class from the preceding period. This is then used to calculate the proportion of the total population in time period n with homozygous recessive (G3,n) and heterozygous (G2,n) genotypes, from which Rn is calculated, the proportion of the population with a resistant phenotype in period n, with Rn = G3,n + G2,nd. Dominance is actioned by the value of d, which is 0 when resistance is assumed recessive, and 1 when it is assumed to be dominant.
Results from the FCM are used by the PM to calculate the average number of infectious bites per mosquito in the population during each time period. From this Qn, the number of infectious bites given by the population as a whole relative to those given by an untreated population, can be calculated for each time period as .
The model does not attempt to capture the effects of mutational processes or stochastic demographic effects on the origin and initial spread of very low numbers of resistance alleles, and so it is assumed that resistant phenotypes are already established at a low frequency in the population at the start of the analysis. Resistance involves a single gene and a simple dominant/recessive process. Moreover, it is assumed that the size and age structure of the population at the start of the PM analysis is that achieved after sustained use in a susceptible population of the insecticide being evaluated, that there is no immigration or emigration, and the proportion of each genotype in the new adults joining the population matches that in the eggs from which they originate. Density dependence is assumed to occur at the mosquito larval stage and the number of adult mosquitoes recruited to the population per unit of time remains constant.
All model parameters are age-independent, apart from background mortality rates and the action of age-linked pesticides, with incremental mortality from fungal biopesticide infection varying according to the number of days since infection. Conventional insecticides affecting a susceptible individual are assumed to be instantly fatal. Mosquitoes choose human hosts at random, and the model does not capture feedback between numbers of infectious bites and the proportion of human hosts with infectious malaria. Malaria-infected mosquitoes never become uninfected. All feeding cycles are of equal duration and mosquitoes bite once in each cycle. All eggs laid are of equal quality and viability. The analysis assumes that malaria infection produces no effects on behaviour, background mortality or fecundity in infected mosquitoes, and fungus-infected mosquitoes that survive and lay eggs are assumed to lay as many eggs at each laying event as uninfected individuals.
Mosquitoes are assumed to contact the chemical or biopesticide when resting after biting a human host, reflecting an application method essentially consistent with IRS. Avoidance behaviour such as outdoor feeding and outdoor resting is not reflected in the coverage values for susceptible mosquitoes since it comprises a method of resistance.
A number of fungal strains have now been tested in laboratory mosquito populations, and a wide range of mortality characteristics have been observed around the basic pattern of initial fungal growth and development followed by an increase in observable mosquito mortality [11, 15, 16, 18, 37]. This suggests that most virulence profiles are potentially available, and in a search for generalizable results this analysis therefore uses highly simplified virulence mortality characteristics, defined by two parameters, ‘initiation day’, the time from infection to the onset of fungus-induced mortality, and the daily mortality rate from that point (Figure 2).
Fungal biopesticides can also impact mosquito feeding propensity and flight capacity in the days before mosquito death . A mosquito which no longer attempts to feed or to lay eggs is effectively dead from the perspectives of fitness and disease transmission. For the purpose of the model therefore, ‘mortality’ encompasses cessation of feeding and reproduction, as well as actual death.
Coverage and virulence
The proportionate reduction in EIR generated by use of a biopesticide is affected by fungal virulence and coverage (Figure 3). For a given level of coverage, similar levels of EIR reduction are achieved by various combinations of the two parameters used to summarize virulence (initiation day and mortality rate (Figure 2)). Unsurprisingly, the longer a fungus takes to initiate mortality, the greater the subsequent mortality rate has to be to maintain a given level of reduction in EIR. There are limits to the EIR reductions that can be achieved at low virulence and/or low coverage.
For equivalent reductions in EIR, selection for resistance is best minimized by high coverage with late initiation day, high mortality rate biopesticides. For example, the lowest selection coefficient associated with a 90% RAIB at 80% coverage is 21%, with day 9 initiation and a 91% mortality rate. At 50% coverage the lowest selection coefficient available in combination with 90% RAIB is 40%.
The temporal dynamics of EIR reduction and resistance evolution are shown in Figure 4. Predictably, more virulent biopesticides give better population-level reductions in EIR to begin with, but they then drive the evolution of resistance more rapidly. The speed of resistance evolution is more sensitive to the timing of mortality onset than to the incremental mortality rate.
The evolutionary dynamics and resulting pattern of control failure differ markedly for different insecticides even when they give identical reductions in EIR in the pre-evolutionary phase (Figure 5). Conventional instant-kill chemical insecticide (with coverage adjusted to achieve the same initial control) fails first. The longest time to product failure is offered by a fungal biopesticide with relatively late mortality initiation, which then kills at a very high rate (Figure 5).
One of the most commonly used classes of conventional insecticides, pyrethroids, have high contact irritancy (also called excito-repellency), causing approximately 50% of mosquitoes contacting treated surfaces to be repelled without acquiring a harmful dose [23–25, 38]. There is no indication of any repellency effects for the fungal biopesticides . For IRS, if 50% of mosquitoes contacting the instant-kill insecticide are unaffected by it, then, for equivalent spray coverage, fungal biopesticides offer better reductions in EIR at all coverage levels, whilst maintaining selection benefits for all but the most virulent strain at the lowest coverage (Figure 8).
Some fungal strains have been shown to have higher virulence in malaria-infected mosquitoes than in those without malaria infection . The trade-off between reducing EIR and resistance management is greatly reduced where fungal virulence is lower in malaria-free mosquitoes, with selection for resistance virtually eliminated if the fungus induces mortality exclusively in malaria-infected mosquitoes (Figure 9).
Variation in the virulence characteristics of potential biopesticides offers scope for selecting strains targeted to provide desirable combinations of reduced transmission and resistance management. A number of virulence phenotypes can provide equivalent levels of EIR reduction (Figure 3), and in general high biopesticide-induced mortality rates commencing as late as possible offer better resistance management for a given level of EIR reduction (Figures 3 and 5). There is nonetheless a trade-off between extending the time taken for resistance evolution to undermine efficacy of a pesticide, and the initial reductions in transmission (Figure 4). In general terms, more virulent fungal strains better reduce transmission initially, but at the cost of stronger selection for resistance, and consequently a shorter useful life (Figure 10).
Although high coverage offers scope to use less virulent fungal strains to reduce EIR, for given virulence parameters, higher levels of coverage also generate stronger selection for resistance, for both conventional and biopesticide interventions. Remembering that the biopesticides must be considered in relation to the best currently used approaches, it is interesting to note that in relative terms, the benefits of biopesticides versus conventional instant-kill insecticides are maximized at high coverage for both transmission control and resistance management (Figure 6).
The relative importance of initial control versus product lifespan depends on a large number of factors, including the availability of alternative replacement treatments, the meaning in terms of human morbidity and mortality of a smaller reduction in EIR at the outset, and the realities of public budgets and other resources. The relative costs and benefits also change if the biopesticide is being considered for use as part of a combination treatment with other interventions [34, 39, 40]. There is, therefore, no simple mathematical optimum for the many possible virulence schedules; the many possibilities need to be considered in context. In so far as it can be done without compromising transmission control, however, it is clearly beneficial to choose the biopesticide that generates the lowest selection for resistance in a particular context. For resistance management, the aim should be to achieve high levels of coverage, allowing less virulent fungal strains to achieve a given level of control, and maximizing their resistance management benefits over instant-kill insecticides.
Even strains sufficiently virulent to match the transmission-reducing characteristics of conventional instant-kill chemical insecticides at matching coverage levels still offer a small benefit in terms of the rate of spread of resistance (Figure 7). Such a resistance management gain would be enhanced by any fitness costs associated with resistance .
The conclusions presented here are independent of the method of resistance (e.g., metabolic or behavioural), provided resistance is genetically determined. It is assumed however that resistance is a binary quality, with mosquitoes either experiencing the full effects of a control measure, or remaining completely unaffected by it. The analysis of the speed of spread of resistance here thus assumes that susceptible mosquitoes experience infections with the specified virulence characteristics, and that resistant mosquitoes have no fungal mortality. In reality, it is more probable that a resistance/tolerance process would operate, with resistant mosquitoes still becoming infected, but experiencing a lower mortality rate than fully susceptible individuals. The spread of resistance would therefore effectively comprise a reduction in fungal virulence, rather than a complete loss of control. Considering the results presented in Figure 4, for example, this would mean that the spread of resistance to the highest virulence biopesticides, rather than comprising a steep function to complete resistance and total loss of transmission control, would move to the curves calculated for sequentially less virulent strains, as resistance converts high virulence strains to low virulence strains, offering even more beneficial resistance management possibilities. Future analyses could explore the impact of hypothetical resistance mechanisms that might operate with respect to conventional and fungal pesticides. The analyses presented here could also be extended to evaluate the impact of malaria infection on mosquito survival, fecundity and behaviour and variation in fecundity with mosquito age.
Certain widely used pyrethroid insecticides have high contact repellency, with studies suggesting that around 50% of mosquitoes landing on treated surfaces may leave before acquiring a fatal dose [23–25, 38]. Whilst this potentially enhances the impact of pyrethroid-treated bed nets on transmission by deflecting mosquitoes away from protected humans before they bite, for IRS it results in mosquitoes surviving to potentially transmit malaria in later feeding cycles . Thus, for this group of conventional insecticides, the composite ‘coverage’ value at a given level of spray cover, would be half that for biopesticides, and could never be greater than 50%. Comparing biopesticide performance with that of a conventional insecticide, and assuming 50% contact repellency (Figure 8) across a full range of coverage values, fungi better reduce transmission than pyrethroid IRS, while still maintaining some resistance management benefits. This suggests that, for all spray coverage values, suitably virulent fungal strains might provide a better option for IRS-based vector interventions than contact-repellent pyrethroids. If only low levels of spray coverage are achievable, replacing repellent pyrethroids with high-virulence fungal treatments could significantly improve the achievable EIR reduction, without significantly increasing selection for resistance, which is in any case relatively weak at low coverage (Figure 8). Where high spray coverage is achievable, replacing pyrethroids with relatively low-virulence fungal treatments could give improvements in both transmission control and resistance management, since the relative fitness of susceptible mosquitoes would be potentially doubled.
The analysis shows that in all cases, having higher fungal-induced mortality in malaria-infected mosquitoes than in uninfected mosquitoes minimizes the fitness costs associated with a given reduction in transmission (Figure 9). The ideal biopesticide from the resistance management perspective would be one that had little or no impact on mosquitoes not infected with malaria, but was strongly virulent in malaria-infected individuals. This might be possible since malaria infection can impose significant metabolic and immunological challenges to mosquitoes [41–44]. There is only a minimal trade-off between transmission control and resistance management in malaria-linked incremental biopesticide mortality. By changing the fitness cost to the mosquito of malaria infection, pesticides working in this way might also exert selection in favour of vector resistance to malaria, further enhancing the transmission-control benefits from the intervention. Strain selection or genetic modification should ideally target this trait. A further development of this principle would be fungal strains which specifically block development of the malaria parasite in the mosquito, or simply act as a delivery mechanism for anti-malaria interventions in the mosquito host (‘paratransgenesis’ [10, 37]), with minimum survival or fecundity costs to the mosquito. It must be noted, however, that this potentially moves selection for resistance from the mosquito to the malaria parasite, which has so far proved extraordinarily adept at evolving its way out of trouble.
This analysis shows that fungal biopesticides have the potential to significantly reduce EIR while imposing only weak selection for resistance. There is always a trade-off between the magnitude of the initial reductions in transmission and maintaining those reductions in the longer term. Given the severe human and economic consequences of malaria transmission, choosing an intervention which does not maximally reduce transmission at the outset requires very careful justification. However, the analyses presented here show that fungal biopesticides can offer equivalent or better reductions in transmission than existing interventions in both the short and long term. This is especially true where existing conventional chemical pesticides have high contact irritancy or resistance to them has already begun to spread. The theoretical analyses presented here should help define the vector mortality profiles required to maximize the sustained malaria control potential of fungal biopesticides, or indeed other novel biological or chemical insecticides.
Average number of infectious bites per mosquito per lifetime
Entomological inoculation rate, the number of infectious bites per person per unit of time
Feeding cycle model
Indoor residual spraying
Lifetime reproductive success
Proportionate reduction in average number of infectious bites per mosquito per lifetime, compared to the number for untreated susceptible mosquitoes.
We thank S Blanford for beautiful survival curves, members of the Research and Policy in Infectious Disease Dynamics Program of the Science and Technology Directorate, Department of Homeland Security and the Fogarty International Center, National Institutes of Health for discussion, and T Ayerst for patience, encouragement and support. This work was supported by a grant (R21 AI088094) from the National Institute for Allergy and Infectious Disease.
Mathematics and Statistics Department, The Faculty of Mathematics, Computing and Technology, The Open University
Center for Infectious Disease Dynamics, Department of Entomology, Pennsylvania State University
Center for Infectious Disease Dynamics, Departments of Biology and Entomology, Pennsylvania State University
Fogarty International Center, National Institutes of Health
Yewhalaw D, Wassie F, Steurbaut W, Spanoghe P, Van Bortel W, Denis L, Tessema DA, Getachew Y, Coosemans M, Duchateau L, Speybroeck N: Multiple insecticide resistance: an impediment to insecticide-based malaria vector control program.PLoS One 2011, 6:e16066.PubMedView Article
Ranson H, N’Guessan R, Lines J, Moiroux N, Nkuni Z, Corbel V: Pyrethroid resistance in African Anopheline mosquitoes: what are the implications for malaria control?Trends Parasitol 2011, 27:91–98.PubMedView Article
Chanda E, Hemingway J, Kleinschmidt I, Rehman AM, Ramdeen V, Phiri FN, Coetzer S, Mthembu D, Shinondo CJ, Chizema-Kawesha E, Kamuliwo M, Mukonka V, Baboo KS, Coleman M: Insecticide resistance and the future of malaria control in Zambia.PLoS One 2011, 6:e24336.PubMedView Article
Trape J-F, Tall A, Diagne N, Ndiath O, Ly AB, Faye J, Dieye-Ba F, Roucher C, Bouganali C, Badiane A, Sarr FD, Mazenot C, Touré-Baldé A, Raoult D, Druilhe P, Mercereau-Puijalon O, Rogier C, Sokhna C: Malaria morbidity and pyrethroid resistance after the introduction of insecticide-treated bednets and artemisinin-based combination therapies: a longitudinal study.Lancet Infect Dis 2011, 11:925–932.PubMedView Article
Van den Berg H, Zaim M, Yadav RS, Soares A, Ameneshewa B, Mnzava A, Hii J, Dash AP, Ejov M: Global trends in the use of insecticides to control vector-borne diseases.Environ Health Perspect 2012, 120:577–580.PubMedView Article
Thomas MB, Read AF: Can fungal biopesticides control malaria?Nat Revs Microbiol 2007, 5:377–383.View Article
Blanford S, Shi W, Christian R, Marden JH, Koekemoer LL, Brooke BD, Coetzee M, Read AF, Thomas MB: Lethal and pre-lethal effects of a fungal biopesticide contribute to substantial and rapid control of malaria vectors.PLoS One 2011, 6:e23591.PubMedView Article
Farenhorst M, Hilhorst A, Thomas MB, Knols BGJ: Development of fungal applications on netting substrates for malaria vector control.J Med Entomol 2011, 48:305–313.PubMedView Article
Mnyone L, Lyimo I, Lwetoijera D, Mpingwa M, Nchimbi N, Hancock P, Russell T, Kirby M, Takken W, Koenraadt C: Exploiting the behaviour of wild malaria vectors to achieve high infection with fungal biocontrol agents.Malar J 2012, 11:87.PubMedView Article
Farenhorst M, Farina D, Scholte E-J, Takken W, Hunt RH, Coetzee M, Knols BGJ: African water storage pots for the delivery of the entomopathogenic fungusMetarhizium anisopliaeto the malaria vectorsAnopheles gambiaes.s. andAnopheles funestus.Am J Trop Med Hyg 2008, 78:910–916.PubMed
Blanford S, Chan BHK, Jenkins N, Sim D, Turner RJ, Read AF, Thomas MB: Fungal pathogen reduces potential for malaria transmission.Science 2005, 308:1638–1641.PubMedView Article
Scholte E-J, Ng'habi K, Kihonda J, Takken W, Paaijmans K, Abdulla S, Killeen GF, Knols BGJ: An entomopathogenic fungus for control of adult African malaria mosquitoes.Science 2005, 308:1641–1642.PubMedView Article
Blanford S, Jenkins NJ, Read AF, Thomas MB: Evaluating the lethal and pre-lethal effect of a range of fungi against adult mosquitoes.Malar J 2012. in press
Bell AS, Blanford S, Jenkins N, Thomas MB, Read AF: Real-time quantitative PCR for analysis of candidate fungal biopesticides against malaria: technique validation and first applications.J Invertebr Pathol 2009, 100:160–168.PubMedView Article
Darbro JM, Graham RI, Kay BH, Ryan PA, Thomas MB: Evaluation of entomopathogenic fungi as potential biological control agents of the dengue mosquito,Aedes aegypti(Diptera: Culicidae).Biocontrol Sci 2011, 21:1027–1047.View Article
Lu D, Pava-Ripoll M, Li Z, Wang C: Insecticidal evaluation ofBeauveria bassianaengineered to express a scorpion neurotoxin and a cuticle degrading protease.Appl Microbiol Biotechnol 2008, 81:515–522.PubMedView Article
St. Leger R, Wang C: Genetic engineering of fungal biocontrol agents to achieve greater efficacy against insect pests.Appl Microbiol Biotechnol 2010, 85:901–907.PubMedView Article
Wang C, St. Leger RJ: A scorpion neurotoxin increases the potency of a fungal insecticide.Nat Biotechnol 2007, 25:1455–1456.PubMedView Article
Mongkalangoon P, Grieco JP, Achee NL, Suwonkerd W, Chareonviriyaphap T: Irritability and repellency of synthetic pyrethroids on anAedes aegyptipopulation from Thailand.J Vector Ecol 2009, 34:217–224.PubMed
Grieco JP, Achee NL, Chareonviriyaphap T, Suwonkerd W, Chauhan K, Sardelis MR, Roberts DR: A new classification system for the actions of IRS chemicals traditionally used for malaria control.PLoS One 2007, 2:e716.PubMedView Article
Chandre F, Darriet F, Duchon S, Finot L, Manguin S, Carnevale P, Guillet P: Modifications of pyrethroid effects associated with kdr mutation inAnopheles gambiae.Med Vet Entomol 2000, 14:81–88.PubMedView Article
Mnyone L, Koenraadt C, Lyimo I, Mpingwa M, Takken W, Russell T: Anopheline and Culicine mosquitoes are not repelled by surfaces treated with the entomopathogenic fungiMetarhizium anisopliaeandBeauveria bassiana.Parasit Vectors 2010, 3:80.PubMedView Article
Mandal S, Sarkar R, Sinha S: Mathematical models of malaria - a review.Malar J 2011, 10:109.View Article
Koella JC, Lynch PA, Thomas MB, Read AF: Towards evolution-proof malaria control with insecticides.Evol App 2009, 2:469–480.View Article
Gourley SA, Liu R, Wu J: Slowing the evolution of insecticide resistance in mosquitoes: a mathematical model.Proc R Soc A 2011, 467:2127–2148.View Article
Hancock P, Thomas MB, Godfray C: An age-structured model to evaluate the potential of novel malaria-control interventions: a case study of fungal biopesticide sprays.Proc R Soc B 2009, 276:71–80.PubMedView Article
Killeen GF, McKenzie FE, Foy BD, Schieffelin C, Billingsley PF, Beier JC: A simplified model for predicting malaria entomologic inoculation rates based on entomologic and parasitologic parameters relevant to control.Am J Trop Med Hyg 2000, 62:535–544.PubMed
Saul AJ, Graves PM, Kay BH: A cyclical feeding model for pathogen transmission and its application to determine vectorial capacity from vector infection-rates.J Appl Ecol 1990, 27:123–133.View Article
Smith DL, McKenzie FE: Statics and dynamics of malaria infections in Anopheles mosquitoes.Malar J 2004, 3:13.PubMedView Article
Hancock P: Combining fungal biopesticides and insecticide-treated bednets to enhance malaria control.PLoS Comput Biol 2009, 5:e1000525.PubMedView Article
Killeen GF, Smith TA: Exploring the contributions of bed nets, cattle, insecticides and excitorepellency to malaria control: a deterministic model of mosquito host-seeking behaviour and mortality.Trans R Soc Trop Med Hyg 2007, 101:867–880.PubMedView Article
Lwetoijera DW, Sumaye RD, Madumla EP, Kavishe DR, Mnyone LL, Russell TL, Okumu FO: An extra-domiciliary method of delivering entomopathogenic fungus,Metharizium anisopliaeIP 46 for controlling adult populations of the malaria vector.Anopheles arabiensis. Parasit Vectors 2010, 3:18.View Article
Fang W, Vega-Rodríguez J, Ghosh AK, Jacobs-Lorena M, Kang A, St. Leger RJ: Development of transgenic fungi that kill human malaria parasites in mosquitoes.Science 2011, 331:1074–1077.PubMedView Article
Kongmee M, Boonyuan W, Achee NL, Prabaripai A, Lerdthusnee K, Chareonviriyaphap T: Irritant and repellent responses ofAnopheles harrisoniandAnopheles minimusupon exposure to bifenthrin or deltamethrin using an excito-repellency system and a live host.J Am Mosq Control Assoc 2012, 28:20–29.PubMedView Article
Farenhorst M, Mouatcho JC, Kikankie CK, Brooke BD, Hunt RH, Thomas MB, Koekemoer LL, Knols BGJ, Coetzee M: Fungal infection counters insecticide resistance in African malaria mosquitoes.Proc Natl Acad Sci U S A 2009, 106:17443–17447.PubMedView Article
Thomas MB, Godfray HCJ, Read AF, van den Berg H, Tabashnik BE, van Lenteren JC, Waage JK, Takken W: Lessons from agriculture for the sustainable management of malaria vectors.PLoS Med 2012, 9:e1001262.PubMedView Article
Ferguson HM, Read AF: Why is the effect of malaria parasites on mosquito survival still unresolved?Trends Parasitol 2002, 18:256–261.PubMedView Article
Rivero A, Ferguson HM: The energetic budget of Anopheles stephensi infected withPlasmodium chabaudi: is energy depletion a mechanism for virulence?Proc R Soc B 2003, 270:1365–1371.PubMedView Article
Lambrechts L, Chavatte J-M, Snounou G, Koella JC: Environmental influence on the genetic basis of mosquito resistance to malaria parasites.Proc R Soc B 2006, 273:1501–1506.PubMedView Article
Ahmed AM, Hurd H: Immune stimulation and malaria infection impose reproductive costs inAnopheles gambiaevia follicular apoptosis.Microbes Infect 2006, 8:308–315.PubMedView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.