Prospective malaria control using entomopathogenic fungi: comparative evaluation of impact on transmission and selection for resistance
Malaria Journalvolume 11, Article number: 383 (2012)
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.
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, F i , 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, I i , is calculated as
i > D
The average probability of survival from start of cycle i to start of cycle (i + 1) is S i , 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.
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 .
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.
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.
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/j.pt.2010.08.004.
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.
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.
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.
World Health Organisation: World malaria report. 2011, World Health Organisation website. Available at http://www.who.int/malaria/world_malaria_report_2011/en/
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.
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.
Beier JC: Malaria parasite develoment in mosquitoes. Annu Rev Entomol. 1998, 43: 519-543. 10.1146/annurev.ento.43.1.519.
Thomas MB, Read AF: Can fungal biopesticides control malaria?. Nat Revs Microbiol. 2007, 5: 377-383. 10.1038/nrmicro1638.
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.
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.
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.
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.
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.
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.
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. 10.1016/j.jip.2009.01.006.
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.
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.
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.
Wang C, St. Leger RJ: A scorpion neurotoxin increases the potency of a fungal insecticide. Nat Biotechnol. 2007, 25: 1455-1456. 10.1038/nbt1357.
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.
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.
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.
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.
Mandal S, Sarkar R, Sinha S: Mathematical models of malaria - a review. Malar J. 2011, 10: 109-10.1186/1475-2875-10-109.
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.
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.
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.
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.
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.
Smith DL, McKenzie FE: Statics and dynamics of malaria infections in Anopheles mosquitoes. Malar J. 2004, 3: 13-10.1186/1475-2875-3-13.
Hancock P: Combining fungal biopesticides and insecticide-treated bednets to enhance malaria control. PLoS Comput Biol. 2009, 5: e1000525-10.1371/journal.pcbi.1000525.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
The authors declare that they have no competing interests.
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.