Skip to main content

Prospective malaria control using entomopathogenic fungi: comparative evaluation of impact on transmission and selection for resistance



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 [16], just as they were last century [7]. 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) [8]. 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 [9]. 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 [10]. 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 [1116].

A wide variety of mortality schedules can be induced in Anopheles by entomopathogenic fungi [17]. 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 [11], dose [18] and malaria-infection status [15, 18], see also [19]. Lethality can also be increased by genetically modifying fungal isolates [2022].

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 [2325]. There is evidence that fungal spores do not repel mosquitoes [26], raising the prospect that, for IRS, fungal biopesticides might more effectively reduce transmission than pyrethroid-based technologies currently in use.


The model

Many malaria transmission models already exist [27], 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 [8]. Other modelling frameworks used to assess the LLA approach are heuristically useful but lack sufficient detail to define target virulence schedules [2830].

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 [3134].

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.

Figure 1
figure 1

FCM model structure for one feeding cycle. Calculation of possible outcomes for feeding cycle i+1 of the model, assuming status when commencing cycle i+1 is age = i, malaria status = m i , fungal infection status = l i , with l0=0 and m0=0. Each arrow represents a probability calculated by the model. If the malaria infection status (number of cycles since infection) of a mosquito when biting is greater than the development time of the malaria parasite in the mosquito, then ‘bites host’ after attacking a human host is recorded as an infectious bite.

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

RAIB = 1 AIB with treatment AIB without treatment

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 Selection Coefficient = 1 LRS for specified mosquito type with intervention LRS for susceptible mosquitoes without intervention . 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.

Table 1 Variables and parameters for the feeding cycle model
Table 2 Values used in FCM for this analysis

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 [36], 76% infection in experimental huts with fungus-impregnated eave curtains [13], and laboratory trials showing >95% infection from treated clay pots [14] or exposure to treated clay tiles [11].

The average number of eggs laid in a given cycle, by mosquitoes surviving to the start of that cycle, F i , is calculated as

F i = m = 0 i 1 l = 0 i 1 f i , m , l v i , m , l V i

This provides the basis for the evaluation of relative fitness using a comparison of values for φ, lifetime egg production, representing LRS, φ = i = 1 λ F i V i .

Comparative levels of transmission control are assessed using u, the average number of infectious bites per mosquito lifetime, u = i = 1 λ I i V i .

The average number of infectious bites during cycle i per mosquito surviving to the beginning of cycle i, I i , is calculated as

I i = m = D i 1 l = 0 i 1 q i , m , l , 2 v i , m , l + q i , m , l , 3 v i , m , l V i i > D

The average probability of survival from start of cycle i to start of cycle (i + 1) is S i , with

S i = m = 0 i 1 l = 0 i 1 s i , m , l v i , m , l V i

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.

Table 3 Variables and parameters for the population model
Table 4 Values used in the population model for this analysis

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 (G 3,n ) and heterozygous (G 2,n ) genotypes, from which R n is calculated, the proportion of the population with a resistant phenotype in period n, with R n = 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 Q n , 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 Q n = M n W n J q .


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).

Figure 2
figure 2

Illustrative survival curves for a range of simple virulence mortality assumptions. Survival curves illustrating mortalities defined by two simple virulence parameters. For each illustrated pair of values, mortality is zero until the specified initiation day, and is thereafter maintained at the indicated fixed daily mortality rate. Initiation day and mortality rate are the two parameters used to define the assumed incremental mortality generated by a given biopesticide infection.

Fungal biopesticides can also impact mosquito feeding propensity and flight capacity in the days before mosquito death [11]. 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.

Figure 3
figure 3

Comparison of virulence characteristics and fitness costs associated with given reductions in EIR. Top panels show different combinations of values for initiation day (x-axis) and daily mortality rate (y-axis) which achieve the denoted reductions in EIR (RAIB). The mortality rate required to achieve a given RAIB increases for later initiation days, up to an initiation day beyond which the target RAIB cannot be achieved, at which point the plots stop. With 50% biopesticide coverage, no virulence parameter combinations achieve 99% RAIB. Bottom panels show the selection coefficients corresponding to the same set of virulence parameter values, eg, the 99% RAIB value plotted for initiation day 2 gives the fitness cost for susceptibility to a biopesticide with initiation day 2 combined with the mortality rate required to achieve a 99% RAIB. Higher selection coefficients indicate stronger selection pressure for resistance.

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.

Figure 4
figure 4

Population level infectious bite rate and proportion resistant for populations exposed to different biopesticides. Top panels show the population reduction in infectious bites per unit of time for each of five different virulence combinations, and the change in this value over time with the spread of resistance to the treatments, shown in bottom panels. 80% coverage assumed throughout.

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).

Figure 5
figure 5

Comparison of four interventions providing a 90% initial reduction in infectious bites. Plots show the change over time in the proportion of resistant individuals (bottom panel) and the percentage reduction in population level infectious bites (top panel) for a mosquito population consistently exposed to one of four vector control treatments, all chosen to give the same 90% initial reduction in EIR.

Clearly, the probability that a mosquito contacts and is affected by a vector-control treatment has a significant impact on both the reduction in EIR and reproductive success. Reductions in EIR improve as coverage is increased, but the strength of selection for resistance also increases (Figure 6, left panels). This illustrates the predictable trade-off between the best transmission control, obtained at high coverage, and the best resistance management, obtained at low coverage. When compared to the currently available alternative, a conventional instant-kill chemical insecticide, however, the relative values for EIR reduction and resistance management with the biopesticides are maximized at the high coverage values which correspond to the best transmission control and the strongest selection pressures for resistance (Figure 6, right panels). Even a biopesticide with sufficiently high virulence to match the initial EIR reduction of instant-kill insecticides at the same coverage levels offers some benefit in terms of useful life (Figure 7). This is because fungus-infected mosquitoes are still able to achieve some reproduction before being killed, thus somewhat reducing the selection for resistance.

Figure 6
figure 6

Comparison of conventional instant-kill chemical insecticide and four biopesticides across a range of coverage values. Lifetime reproductive success with interventions as a proportion of LRS for untreated mosquitoes (top left panel) and as a proportion of LRS for mosquitoes treated with an instant-kill insecticide (top right panel). Reduction in average infectious bites per mosquito lifetime with interventions, compared to the value for untreated mosquitoes (bottom left panel), 0 = no reduction in infectious bites, 1.00 = no infectious bites. Reduction in infectious bites with interventions vs untreated mosquitoes, compared to the reduction achieved using a conventional instant kill insecticide (bottom right panel), 1.00 means reduction equal to that achieved by instant-kill insecticide.

Figure 7
figure 7

Comparison of resistance spread and consequent increases in infectious bites with instant-kill and fungal biopesticides. Biopesticide virulence selected to give pre-resistance EIR reduction matching instant-kill pesticides at 80% or 30% coverage. Plots show the proportion of the population with resistant phenotypes, and the corresponding values for population-level reduction in infectious bites per unit of time compared to an untreated population.


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 [2325, 38]. There is no indication of any repellency effects for the fungal biopesticides [26]. 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).

Figure 8
figure 8

Comparison between biopesticides and instant-kill insecticide with 50% contact irritancy, across range of coverage values. Lifetime reproductive success with interventions as a proportion of LRS for untreated mosquitoes (top left panel) and as a proportion of LRS for mosquitoes treated with an instant-kill insecticide with 50% contact irritancy (top right panel). Reduction in average infectious bites per mosquito lifetimewith interventions, compared to the value for untreated mosquitoes (bottom left panel), 0 = no reduction in infectious bites, 1.00 = no infectious bites. Reduction in infectious bites with interventions vs untreated mosquitoes, compared to the reduction achieved using a conventional instant kill insecticide with 50% contact irritancy (bottom right panel), 1.00 means reduction in AIB equal to that achieved by instant-kill insecticide with 50% contact irritancy.

Malaria interactions

Some fungal strains have been shown to have higher virulence in malaria-infected mosquitoes than in those without malaria infection [15]. 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).

Figure 9
figure 9

Differential mortality in malaria-infected and malaria-free mosquitoes. Comparison of speed of spread of resistance and consequent loss of transmission control for populations treated with one of five fungal biopesticides with differential mortality rates in malaria-infected mosquitoes. Plots show the proportion of the population with resistant phenotypes, and the corresponding values for population-level reductions in infectious bites per unit of time compared to an untreated population. The biopesticides all have day 3 initiation of a 72% daily mortality rate for malaria infected mosquitoes, giving an initial 99% reduction in infectious bites per time period.


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).

Figure 10
figure 10

Schematic illustration of the trade-off between initial levels of control and product useful life. The green and red areas highlight the differences in infectious bites resulting from a choice between a high-virulence intervention with high initial infectious bite reduction and high selection for resistance versus a lower-virulence intervention with lower initial benefits and lower selection for resistance. The red area highlights the additional infectious bites associated with choosing the lower virulence option, before resistance spreads. The green area shows the additional infectious bites associated with choosing the high-virulence option in the long term once resistance spreads.

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 [8].

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 [2325, 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 [24]. 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 [4144]. 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


Late-life acting


Lifetime reproductive success


Population model


Proportionate reduction in average number of infectious bites per mosquito per lifetime, compared to the number for untreated susceptible mosquitoes.


  1. 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-10.1371/journal.pone.0016066.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  2. 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. 10.1016/

    Article  CAS  PubMed  Google Scholar 

  3. 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-10.1371/journal.pone.0024336.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. 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. 10.1016/S1473-3099(11)70194-3.

    Article  CAS  PubMed  Google Scholar 

  5. 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. 10.1289/ehp.1104340.

    Article  PubMed Central  PubMed  Google Scholar 

  6. World Health Organisation: World malaria report. 2011, World Health Organisation website. Available at

    Google Scholar 

  7. Cohen J, Smith D, Cotter C, Ward A, Yamey G, Sabot O, Moonen B: Malaria resurgence: a systematic review and assessment of its causes. Malar J. 2012, 11: 122-10.1186/1475-2875-11-122.

    Article  PubMed Central  PubMed  Google Scholar 

  8. Read AF, Lynch PA, Thomas MB: How to make evolution-proof insecticides for malaria control. PLoS Biol. 2009, 7: e58-10.1371/journal.pbio.1000058.

    Article  Google Scholar 

  9. Beier JC: Malaria parasite develoment in mosquitoes. Annu Rev Entomol. 1998, 43: 519-543. 10.1146/annurev.ento.43.1.519.

    Article  CAS  PubMed  Google Scholar 

  10. Thomas MB, Read AF: Can fungal biopesticides control malaria?. Nat Revs Microbiol. 2007, 5: 377-383. 10.1038/nrmicro1638.

    Article  CAS  Google Scholar 

  11. 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-10.1371/journal.pone.0023591.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  12. 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. 10.1603/ME10134.

    Article  PubMed  Google Scholar 

  13. 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-10.1186/1475-2875-11-87.

    Article  PubMed Central  PubMed  Google Scholar 

  14. 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 fungus Metarhizium anisopliae to the malaria vectors Anopheles gambiae s.s. and Anopheles funestus. Am J Trop Med Hyg. 2008, 78: 910-916.

    PubMed  Google Scholar 

  15. 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. 10.1126/science.1108423.

    Article  CAS  PubMed  Google Scholar 

  16. 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. 10.1126/science.1108639.

    Article  CAS  PubMed  Google Scholar 

  17. 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

    Google Scholar 

  18. 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. 10.1016/j.jip.2009.01.006.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. 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. 10.1080/09583157.2011.597913.

    Article  Google Scholar 

  20. Lu D, Pava-Ripoll M, Li Z, Wang C: Insecticidal evaluation of Beauveria bassiana engineered to express a scorpion neurotoxin and a cuticle degrading protease. Appl Microbiol Biotechnol. 2008, 81: 515-522. 10.1007/s00253-008-1695-8.

    Article  CAS  PubMed  Google Scholar 

  21. 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. 10.1007/s00253-009-2306-z.

    Article  CAS  PubMed  Google Scholar 

  22. Wang C, St. Leger RJ: A scorpion neurotoxin increases the potency of a fungal insecticide. Nat Biotechnol. 2007, 25: 1455-1456. 10.1038/nbt1357.

    Article  CAS  PubMed  Google Scholar 

  23. Mongkalangoon P, Grieco JP, Achee NL, Suwonkerd W, Chareonviriyaphap T: Irritability and repellency of synthetic pyrethroids on an Aedes aegypti population from Thailand. J Vector Ecol. 2009, 34: 217-224.

    Article  PubMed  Google Scholar 

  24. Grieco JP, Achee NL, Chareonviriyaphap T, Suwonkerd W, Chauhan K, Sardelis MR, Roberts DR: A new classification system for the actions of IRS chemicals traditionally used for malaria control. PLoS One. 2007, 2: e716-10.1371/journal.pone.0000716.

    Article  PubMed Central  PubMed  Google Scholar 

  25. Chandre F, Darriet F, Duchon S, Finot L, Manguin S, Carnevale P, Guillet P: Modifications of pyrethroid effects associated with kdr mutation in Anopheles gambiae. Med Vet Entomol. 2000, 14: 81-88. 10.1046/j.1365-2915.2000.00212.x.

    Article  CAS  PubMed  Google Scholar 

  26. 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 fungi Metarhizium anisopliae and Beauveria bassiana. Parasit Vectors. 2010, 3: 80-10.1186/1756-3305-3-80.

    Article  PubMed Central  PubMed  Google Scholar 

  27. Mandal S, Sarkar R, Sinha S: Mathematical models of malaria - a review. Malar J. 2011, 10: 109-10.1186/1475-2875-10-109.

    Article  Google Scholar 

  28. Koella JC, Lynch PA, Thomas MB, Read AF: Towards evolution-proof malaria control with insecticides. Evol App. 2009, 2: 469-480. 10.1111/j.1752-4571.2009.00072.x.

    Article  Google Scholar 

  29. 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. 10.1098/rspa.2010.0413.

    Article  Google Scholar 

  30. 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. 10.1098/rspb.2008.0689.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. 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 Central  CAS  PubMed  Google Scholar 

  32. Saul AJ, Graves PM, Kay BH: A cyclical feeding model for pathogen transmission and its application to determine vectorial capacity from vector infection-rates. J Appl Ecol. 1990, 27: 123-133. 10.2307/2403572.

    Article  Google Scholar 

  33. Smith DL, McKenzie FE: Statics and dynamics of malaria infections in Anopheles mosquitoes. Malar J. 2004, 3: 13-10.1186/1475-2875-3-13.

    Article  PubMed Central  PubMed  Google Scholar 

  34. Hancock P: Combining fungal biopesticides and insecticide-treated bednets to enhance malaria control. PLoS Comput Biol. 2009, 5: e1000525-10.1371/journal.pcbi.1000525.

    Article  PubMed Central  PubMed  Google Scholar 

  35. 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. 10.1016/j.trstmh.2007.04.022.

    Article  PubMed Central  PubMed  Google Scholar 

  36. Lwetoijera DW, Sumaye RD, Madumla EP, Kavishe DR, Mnyone LL, Russell TL, Okumu FO: An extra-domiciliary method of delivering entomopathogenic fungus, Metharizium anisopliae IP 46 for controlling adult populations of the malaria vector. Anopheles arabiensis. Parasit Vectors. 2010, 3: 18-10.1186/1756-3305-3-18.

    Article  PubMed  Google Scholar 

  37. 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. 10.1126/science.1199115.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  38. Kongmee M, Boonyuan W, Achee NL, Prabaripai A, Lerdthusnee K, Chareonviriyaphap T: Irritant and repellent responses of Anopheles harrisoni and Anopheles minimus upon exposure to bifenthrin or deltamethrin using an excito-repellency system and a live host. J Am Mosq Control Assoc. 2012, 28: 20-29. 10.2987/11-6197.1.

    Article  CAS  PubMed  Google Scholar 

  39. 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. 10.1073/pnas.0908530106.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  40. 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-10.1371/journal.pmed.1001262.

    Article  PubMed Central  PubMed  Google Scholar 

  41. Ferguson HM, Read AF: Why is the effect of malaria parasites on mosquito survival still unresolved?. Trends Parasitol. 2002, 18: 256-261. 10.1016/S1471-4922(02)02281-X.

    Article  PubMed  Google Scholar 

  42. Rivero A, Ferguson HM: The energetic budget of Anopheles stephensi infected with Plasmodium chabaudi: is energy depletion a mechanism for virulence?. Proc R Soc B. 2003, 270: 1365-1371. 10.1098/rspb.2003.2389.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  43. 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. 10.1098/rspb.2006.3483.

    Article  PubMed Central  PubMed  Google Scholar 

  44. Ahmed AM, Hurd H: Immune stimulation and malaria infection impose reproductive costs in Anopheles gambiae via follicular apoptosis. Microbes Infect. 2006, 8: 308-315. 10.1016/j.micinf.2005.06.026.

    Article  CAS  PubMed  Google Scholar 

Download references


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.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Penelope A Lynch.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

PL developed and applied the model and contributed to study inception, interpretation of results and drafting of the manuscript. UG reviewed the model and contributed to drafting of the manuscript. MT and AR contributed to study inception, interpretation of results and drafting of the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material


Additional file 1: Appendix A. Derivation of the feeding cycle model. Detailed description of variables and calculations used in the feeding cycle model. (PDF 173 KB)


Additional file 2: Appendix B. Derivation of the population model. Detailed description of variables and calculations used in the population model. (PDF 163 KB)

Authors’ original submitted files for images

Rights and permissions

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 (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Lynch, P.A., Grimm, U., Thomas, M.B. et al. Prospective malaria control using entomopathogenic fungi: comparative evaluation of impact on transmission and selection for resistance. Malar J 11, 383 (2012).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: