Skip to main content

A malaria transmission-directed model of mosquito life cycle and ecology

Abstract

Background

Malaria is a major public health issue in much of the world, and the mosquito vectors which drive transmission are key targets for interventions. Mathematical models for planning malaria eradication benefit from detailed representations of local mosquito populations, their natural dynamics and their response to campaign pressures.

Methods

A new model is presented for mosquito population dynamics, effects of weather, and impacts of multiple simultaneous interventions. This model is then embedded in a large-scale individual-based simulation and results for local elimination of malaria are discussed. Mosquito population behaviours, such as anthropophily and indoor feeding, are included to study their effect upon the efficacy of vector control-based elimination campaigns.

Results

Results for vector control tools, such as bed nets, indoor spraying, larval control and space spraying, both alone and in combination, are displayed for a single-location simulation with vector species and seasonality characteristic of central Tanzania, varying baseline transmission intensity and vector bionomics. The sensitivities to habitat type, anthropophily, indoor feeding, and baseline transmission intensity are explored.

Conclusions

The ability to model a spectrum of local vector species with different ecologies and behaviours allows local customization of packages of interventions and exploration of the effect of proposed new tools.

Background

Malaria is transmitted by the blood feeding of infectious female Anopheles mosquitoes, and understanding mosquito ecology and population dynamics can inform how best to defeat malaria. Malaria is an important global health issue, causing over half a billion cases and on the order of one million deaths a year [1], and is the focus of a global eradication campaign announced in 2007. Basic vector ecology is a fundamental driver of transmission patterns, and changes in land usage [2] or land modification can dramatically change transmission for better or worse. The growing urbanization in Africa is a powerful current example of such phenomena [3]. Climate and weather affect larval development and parasite maturation within the infected mosquito, and spatial models are able to predict malaria prevalence based primarily on climate details in the absence of interventions [4]. This climate-driven predictability has broken down more recently, possibly due to more widespread interventions such as insecticide-treated bed nets [5], but predictive modelling for global eradication incorporates these geographic effects on baseline transmission. Geographic variation and spatial effects become increasingly important as heterogeneity in transmission allows malaria to persist in some areas while other areas achieve elimination but remain at risk of reintroduction [6].

A successful global eradication campaign will include substantial vector control components, and mathematical models for planning eradication will benefit from accurate and robust representation of the basic vector transmission ecology in each area of interest as well as the ability to incorporate interventions both singly and in combination. Vector population dynamics exhibit latencies such as the time required for sporogony. Spatial processes include vector oviposition, larval habitat, host-seeking, and migration. For aptly modelling eradication, representation of the steady state is not sufficient; elimination-predictive models may need to be accurate at very low prevalence. Finally, models must address sensitivity of results to model parameters and assumptions whenever presenting possible routes to eradication.

Mathematical modelling of the vector-borne transmission of malaria dates back to the early dynamical models of Ross and Macdonald [7, 8], the classical assumptions of which have been clearly exposited [9]. Next steps in vector modelling included cyclical feeding models, which were easier to parameterize from field data and more accurately tracked the mosquito life cycle [10, 11]. Other models emphasize the effect of rainfall and temperature correlations to transmission [12], or compute the Entomological Inoculation Rate (EIR) driven by human infectiousness [13]. Recent work has focused on the effect of hydrology on larval habitat and vector prevalence [14], and on vector population dynamics [15]. Other groups have built comprehensive simulations for both the vector transmission dynamics and within-human parasite dynamics [16–19]. Vector population models have been constructed for other vector-borne diseases, such as dengue [20].

The present work introduces a vector model which has detailed vector population resolution for near elimination phenomena, tracks explicit latencies of larval development and sporogony, implements closed-loop population dynamics, and can implement a wide variety of vector control interventions in combination. Careful attention is given to vector behaviours, such as host preference and feeding locations, and the effects of these parameters on intervention effectiveness are explored. Alternative implementations of this model are discussed, along with parameter sensitivities. The present model is then exercised on several issues of local elimination for simulations based on transmission patterns for a single-location in Tanzania, varying baseline transmission and vector bionomics, and key results are discussed.

Methods: model design

Aquatic habitat

Available larval habitat is a primary driver of local mosquito populations, and different mosquito species can have different habitat preferences, with utilization of an ecological niche driving speciation in some cases [21]. Classifications of larval habitat include temporary, permanent or semi-permanent [22], and some species, such as Anopheles gambiae ss and Anopheles arabiensis will share an ecological niche for larval habitat [23]. Humans can affect available habitat through terrain changes which affect hydrology, through agricultural practices, such as rice cultivation [2], or through creating or eliminating standing water. Remote sensing through satellite imagery is becoming a powerful tool for mapping vector ecology [2, 23], and this trend will most likely continue to increase as eradication planning drives increasing data requirements. Several detailed models already exist of habitat and the impact of rainfall, temperature, humidity, and soil quality [14, 24, 25].

Rainfall and humidity can strongly affect available larval habitat [14, 22, 23], although this depends on the mosquito species and its habitat preference. In fact, preference can be more specific than the species level, as Anopheles funestus exhibits differences in population responses to rainfall which are correlated with chromosomal diversity [26]. Rainfall, rather than habitats with water, is best correlated with numbers of An. gambiae s.l. This effect is not as strong as it is for culicines, nor is it universal, since An. gambiae s.l. have been found in stable aquatic habitats [22]. Even An. funestus, which prefers more semi-permanent larval habitat, has a rainfall dependence in its larval habitat [27], partly due to vegetation on edges of water [26] and the interaction of rainfall with agricultural schedules for crops such as rice.

In the present model, different models for larval habitat are developed for temporary, semi-permanent, permanent, and human-driven habitats. Temporary habitat Htemp in a grid of diameter Dcell increases with rainfall Prain and decays with a rate Ï„temp proportional to the evaporation rate driven by temperature T (K) and humidity RH:

H t e m p + = P r a i n K t e m p D c e l l 2 - H t e m p ( Δ t τ t e m p )
1 τ t e m p = ( 5 . 1 x 1 0 1 1 P a ) e - 5 6 2 8 . 1 K T K k t e m p d e c a y . 0 1 8 k g ∕ m o l 2 π R T K 1 - R H

in which the exponential results from the Clausius-Clayperon relation, the root is from the expression for vapour evaporation rates due to molecular mass given a partial pressure, and the constant is the Clausius-Clayperon integration constant multiplied by a factor ktempdecay to relate mass evaporation per unit area to habitat loss. The value of ktempdecay is initially chosen to set the habitat half-lives near 1 day for hot and dry conditions and 2-3 weeks for more typical tropical conditions. The variation in Ï„temp with temperature T and humidity RH can be seen in Figure 1a. Semi-permanent habitat increases with a constant Ksemi Dcell2Prain and decays with a longer time constant Ï„semi. Permanent habitat is fixed at Kperm Dcell2, and human population-driven habitat is calculated as population N* Kpop. The values of ktempdecay, Ï„semi, Ksemi, and Ktemp can be fit to local data on vector abundance by species over time or to local data on EIR to tailor a simulation to a specific setting.

Figure 1
figure 1

Effects of climate and weather on vector populations. a) Effect of temperature and humidity on time constant Ï„ temp for temporary rainfall-driven larval habitat. The habitat decay is faster for warmer and drier weather. b) Temperature effects on duration larval development, with the functional form from [4] and the present Arrhenius formulation. c) Temperature effects on duration of sporogony. The traditional degree-day formula and the present Arrhenius function are plotted, along with Beier's data from [49].

Larval development and mortality rates are affected by a variety of factors including weather and densities of other larvae. Climate and weather affect not only larval habitat availability but also larval development rates and larval mortality [15]. The duration of larval development is a decreasing function of temperature [4], and the present model replaces earlier mathematical formulations [4] with an Arrhenius temperature-dependent rate a1exp(a2/TK) as seen in Figure 1b. In some cases, this temperature-dependent rate must be modified by local larval density [28], although the presented results do not include such a modification. Rainfall and temperature then combine through habitat creation and larval development to create varying local patterns of distribution by larval instar [23], and larval mortality and development duration determine pupal rates [29].

Heavy rainfall can directly kill larvae by dislodging them from habitat and causing them to dry out [28]. Other factors increasing mortality are cannibalism of 1st instar larvae by 4th instar larvae and overpopulation of larval habitat acting to reduce food availability.

The present model includes this preferential survival of older larvae during overpopulation conditions by only treating as viable those new larvae, which do not cause the larval population to exceed capacity. If capacity shrinks so that the existing population exceeds available capacity, mortality is increased by the degree of overpopulation. With all these factors taken together, about 2-8 percent of larvae typically survive egg to adult [28]. The present model includes a daily larval mortality, which translates into a probability of survival of larval development as a function of temperature and mortality rate presented in Figure 2a. The larval survival plotted is from egg hatch to emergence, not from oviposition to adult maturity, which is significantly less due to egg survival and death during the immature phase.

Figure 2
figure 2

Intermediate outputs which affect vector population or disease transmission dynamics. a) Larval survival of development as a function of temperature and larval mortality. Survival is from successful egg hatch to emergence; survival from oviposition would be significantly lower. b) Adult survival of sporogony as a function of temperature and adult mortality. At lower temperatures, mosquitoes spend longer in each progress queue and the overall effect of a daily mortality is greater. Adult mortality can be artificially increased through interventions such as bed nets, insecticide spraying, and baited traps.

Improved cohort model

There are different possible implementations of the basic model, each with different computational efficiencies, resolutions, and flexibilities. Possible implementations include a modified cohort simulation, a cohort simulation with explicit mosquito ages, a simulation of every individual mosquito in the population, as well as a simulation of a sampled subset of mosquitoes to represent the population as the whole. The basic model is presented in the context of the modified cohort simulation with explanations of the modifications for individual mosquitoes. In the modified cohort simulation, rather than representing the entire population by three compartments for susceptible, latently infected, and infectious mosquitoes, the simulation dynamically allocates a cohort for every distinct state, and the cohort maintains the count of all mosquitoes in that state. This allows temperature-dependent progression through sporogony as described below, even with a different mean temperature each day, with no mosquitoes passing from susceptible to infectious before the full discrete latency. For the cohort simulation with explicit ages to allow modelling of senescence, mosquito age is part of the state definition, and many more cohorts are required to represent the population. The overall progression of cohorts or individual mosquitoes through different states is outlined in Figure 3.

Figure 3
figure 3

Vector development state space. All eggs of a similar state (species, gender, habitat, Wolbachia type) hatching in a time step begin larval development as a cohort. The only changes to this cohort are the population and the progress, and each time step, mortality reduces the population and progress increments by the Arrhenius temperature-driven rate multiplied by the time step. The progress added can vary depending on the daily temperature and is not constrained to be constant or an integer number of total days, so n1 would be the total development period at the mean temperature of the first time step. When progress through development is complete for a cohort, emergence occurs, and the cohort begins the latency to blood feeding as immature emerged adults. This latency can last for several hours up to several days, at which point the cohort begins the cycle of blood feeding. Adults infected in a time step are removed from their cohort and a new cohort is created for newly infected adults. This new cohort then proceeds through the infected development queue, with mortality reducing the population and temperature-dependent incrementing of progress. Once sporogony is complete, the cohort becomes infectious and remains so until the population is reduced to zero, at which point the cohort is de-allocated.

Immature mosquito populations

Upon emergence, there is a period of hours to days before bloodmeal-seeking begins [30]. This period is represented in the model as a fixed latency, during which predators and interventions such as outdoor spraying can still cause mortality. At the end of this immature interval, before the start of host-seeking, female mosquitoes mate. Male mosquitoes are included in the simulation to allow simulation of the mosquito population genetic structure, as well as interventions and phenomena such as release of modified males or mosquitoes with Wolbachia infection. Each female is mated once, with fertility only if the male is not sterile and there is no cytoplasmic incompatibility due to Wolbachia type [31]. Mating outcomes are based on the current distribution of male mosquitoes. Sterile-mated females will blood-feed, but do not produce viable eggs.

Adult mosquito populations

Host-seeking and blood-feeding are key aspects of the reproductive life of an adult female Anopheles, and these are also the key aspects for malaria transmission. After completion of post-emergence maturation, female adults enter a cycle of feeding and egg-laying which will consume the rest of their lives. Female Anopheles mosquitoes bloodfeed every 1/f = 2-4 days [11], and in the model, a fixed proportion (= fΔt) of all mosquitoes in a state cohort attempt to feed during a time step Δt. Subsequent versions of this model include a state with a timer for blood feed processing, which replaces the draw for fraction feeding each time step. The total number of mosquitoes in a cohort that attempt to feed during a time step are then stochastically sorted into a variety of outcomes depending on vector behaviour, host availability, and interventions, such as insecticide-treated nets (ITN) and indoor residual spraying (IRS) as described in Figure 4.

Figure 4
figure 4

Calculation of outcomes for each mosquito every time step in the presence of combined interventions. Each choice has a defined probability, and the conditional probabilities can be summed for each overall possible outcome as described in the Appendix. Bed nets can kill or not, and vector feeding time can be adjusted to change the proportion of bites during the period protected by nets. Indoor feeding and resting can be split by adding in an additional decision fork after indoor and outdoor feeds. After a successful indoor feed, a mosquito must make it to an oviposition site alive to lay eggs and survive. Closed loop egg-laying allows interference by interventions to eventually limit population sizes. Individual resolution of the human population ensures that only those infectious mosquitoes that successfully pass through a gauntlet to get to a human successfully can transmit infection.

Possible outcomes of an attempted feeding include death, survival without feeding, successful feed on a human, or successful feed on an animal. A binary decision tree is created for progress through a feeding cycle as seen in Figure 4. If a feed is attempted, the first branching point is the choice of host type depending on the vector host preference, and if human is selected, the location of feeding is chosen based on the vector indoor feeding probability. Each possible choice is thus conditional on arriving at that stage of the decision tree, allowing simpler definitions of efficacy. Blocking efficacy of nets is specified as the probability that a net blocks a feed, given an attempted indoor feed on a protected human, rather than the reduction in overall successful biting. This binary structure allows logical combination of the effects of ITNs and IRS and makes it simple to add new interventions to the model. An indoor feed only occurs if the net does not block the feed and the treated net does not kill the mosquito. Mosquitoes that complete a feed are eligible to rest on an IRS-treated wall with a specified killing efficacy. Thus the effects of ITNs and IRS in the same house are not independent, and blocked feeds reduce the number of mosquitoes that arrive at the IRS section of the decision tree. Successful feeds on humans have an additional draw for whether the mosquito is infected with Plasmodium or not which depends on human infectiousness, and the conditional probability of surviving a feed on an infectious individual. En route to assembling the distribution of feeding outcomes, human biting rate and entomological inoculation rate (EIR) are calculated, including all human feeds whether or not the mosquito survives, as mosquito death can occur before, during, or after feeding, but transmission can occur even for feeds which the mosquito does not survive. The presented version of the model, with a draw for number feeding each day and without the timer for blood feed processing, calculates the outcomes for the full feeding cycle including blood meal processing and oviposition survival, but assembles these into outcomes in a single calculation with eggs laid that time step for those feeds. The detailed equations for feeding outcomes are contained in Additional File 1.

The effect of a local mosquito species population on disease transmission depends on several species-specific characteristics. Among the most important is anthropophily (phumanfeed), the fraction of bites, which are taken on a human host, or human blood index [32, 33]. Indoor versus outdoor feeding and resting, represented as an probability of feeding indoors (pindoorfeeding), is another important behaviour of the local mosquito population, especially when indoor interventions such as IRS and ITN's are introduced. Vectors which feed predominantly indoors can be decimated by these interventions, while those species which feed outdoors will not experience the same applied mortality. Indoor feeding and resting are not necessarily equal for a mosquito species, and this would be simple to implement in the present outcome calculation structures. In the presented sample simulations, indoor feeds are associated with indoor resting and outdoor feeds with outdoor resting, although a given species may have a mixed proportion of indoor and outdoor feeds.

The most important factor for baseline transmission is the adult mortality [9, 10], which can be calculated per day or per feeding cycle. Mosquitoes also experience additional mortality at high temperatures with low humidity [4]. The formula of Martens for daily survival to a temperature-dependent mortality rate, with T in Celsius, is approximated as (.001e(T-32)), which does not have the mathematical pathologies at the roots of Martens' polynomial. Mosquitoes do exhibit age effects and senescence in laboratory settings, and senescence has been observed in field studies as well [34]. In fact, mosquitoes have not been found in the field having taken more than 14 feeds [11]. These possible effects are studied by adding age to the state space, which results in a much larger number of cohorts, and adding an age-dependent mortality rate to the standard daily or feeding cycle mortality.

In the cohort model, the number of eggs laid per time step is calculated from the number of successful feeds on humans and animals occurring in that time step, with corrections for number of eggs per feeding type. There is no delay currently in the present model for egg-production, and the population growth dynamics are constrained by the days between feeds and full larval-development latencies. However, oviposition timers can be incorporated both in the individual-mosquito based simulations and in the cohort model. Determining the number of eggs from successful blood feeds allows second order effects of interventions on the mosquito population to be captured, which is not possible in models which utilize a pre-determined temporal pattern for emergence rate of mosquitoes.

Infection

A bite on an infected human can result in mosquito infection with Plasmodium, with a probability of infection dependent on a variety of human and mosquito factors [35–44]. In general, human infectivity tends to increase with gametocyte densities in a typical blood meal of 1-3 μl [35], but this can be reduced by high parasitaemia provoking an inflammatory cytokine response [45, 46], by age and immunity of the human host [13], and by gametocyte-killing drug treatments [47]. Once within the mosquito gut, Plasmodium progresses through several stages of development in the mosquito finally resulting in sporozoites within the salivary glands which can infect human hosts [48]. The mosquito attempts to avoid infection through various defenses against Plasmodium gametocytes [49] and melanotic encapsulation of its oocysts [50].

The effect of weather and climate on malaria transmission is seen again in temperature-dependent latencies in sporozoite development [4, 49, 51, 52]. The development times are seen in Figure 1c, and the corresponding survival probabilities are plotted in Figure 2b. These are included in traditional continuous compartmental models as a factor for mosquito survival of this latency, which is multiplied by the rate of change in infectious mosquitoes. This effect can be implemented in cyclical model as a changing probability of surviving incubation with fixed probability of surviving a feeding cycle [53]. This cohort implementation avoids instantaneous transport from susceptible to infected status, even scaled appropriately for steady state. Steady states are rare to non-existent in malaria transmission: seasonality in temperature and rainfall changes vector population sizes and infection rates, monthly rainfall for the same month varies from year to year, and human population infectiousness may not be at the same level at the same time each year, all of which may affect the impact of interventions as a function of their timing. Therefore, it is important that full latencies are enforced, with infectious mosquitoes not appearing until completion of the intervening stages.

Progress towards infectiousness is included as a state variable, and the infection state variable is not changed from infected to infectious until the progress state variable reaches completion. Enforcement of larval and immature latencies similarly captures the dynamics for population growth. Mosquitoes of the same state infected the same evening become a new state cohort in the simulation, and each time step, progress towards infectiousness is incremented by the temperature-dependent rate. At infection, the number infected is subtracted from the population of the uninfected cohort, and a new cohort is allocated with the newly infected population and zero progress towards sporogony. Upon completion, either the cohort is either merged with an identical-state infectious mosquito state cohort, or maintained separately if no identical state likely exists, as in the case of age-tracking.

Outcome probabilities can change in response to infection status. Sporozoite infection of the salivary glands can result in increased feeding mortality when infectious [54, 55]. Fecundity can be affected as smaller egg-batch size is observed due to maturing oocysts [56] but not salivary-gland sporozoites [57]. Once infectious, a mosquito's bites have a probability of infecting a human host, whether the mosquito survives the feed or dies during or after the feed. Probability of human infection from a sporozoite-positive bite can be calculated from field data or laboratory experiments, with a value of approximately 0.5 per bite probably being reasonable [58].

Interventions

Simulations of vector populations in the absence of interventions, such as bed nets, are important, but the purpose of the present model is to evaluate the effects of interventions singly and in combination, especially in the global eradication context. Key issues include incorporation of the effects of each of the full spectrum of interventions to examine possible effects and determining how interventions combine their effects. Some interventions target adult female Anopheles feeding, and these include insecticide treated bed nets [59–61], indoor residual spraying and screening [62]. Figure 4 shows how each intervention affects the feeding cycle as discussed above and how total outcomes can be calculated when interventions are combined, with full equations in Additional File 1.

Other interventions affect the population through the larval stage, either directly with larvicides [63, 64] and larval predators or indirectly through habitat destruction. Land usage modification, either intentional or unintentional, due to urbanization, agriculture, or draining swamps can have powerful effects [62, 63]. Larval control options available in the model include temporary increases in larval mortality through larvicides in a subset of local habitat, repeated treatments with sufficient mortality to render a fraction of local habitat unavailable for longer intervals, or land usage reducing larval habitat long term. Depending on the option, the model either implements a temporary increase in larval mortality in a subset of the habitat carrying capacity, or proportionately reduces the habitat carrying capacity for the duration of effect.

Individual mosquito model

In addition to state cohorts, this basic model can be implemented through simulation of every individual mosquito or simulation of a subset of individual mosquitoes to represent the full population. Each mosquito's state contains the same features as the state cohort model, with status, timers for transition to adult from immature and infected to infectious, mating status and Wolbachia infection, and age. An oviposition timer to enforce a fixed feeding cycle may be included as well. If mosquitoes are sampled and a subset used to represent the local population, each sample mosquito will have an associated sampling weight as well.

Setting up simulations

Vector dynamics are simulated for single human populations well-mixed with multiple vector populations. All simulations are based upon a single-location with temperature [19] and rainfall [65] based upon lat-long (-8.5, 36.5) in Tanzania. Three local vector populations are simulated: An. gambiae s.s., An. arabiensis, and An. funestus. Anopheles gambiae and An. arabiensis are modelled to track the rainfall with the short temperature and humidity-dependent time constant, while An. funestus larval habitat integrates rainfall with a smaller forcing term and decays with a much slower time constant, here set to 100 days to correspond to the length of an agricultural season. The habitat scaling parameters and habitat-specific time constants were obtained by simulation of one species at a time, comparing to measurements of local EIR by species. Parameters which exhibit high uncertainty or geographic variability, such as the host preference of An. arabiensis, are studied over broad numerical ranges for their impact on results. A simplified human disease model is used in all simulations, with a constant latent period of 22 days from bite to infectiousness to mosquitoes, and exponentially-distributed period of infectiousness with mean 180 days. Infectiousness is a constant 0.2, without development of immunity to allow resolution of vector-specific effects. Superinfection is allowed, with a maximum of five simultaneous infections. General model and simulation-specific parameters and their values are summarized in Table 1.

Table 1 Model and simulation parameters

Results and Discussion

Baseline vector population and transmission dynamics are simulated for the single location simulation described in Methods. The habitat scaling parameters are varied to show different baseline EIRs with the same weather-driven seasonality. The total vector population when summed across all three species is seen in Figure 5a, with the rainfall patterns and temperature in 5 b. Figure 5a shows the effect of scaling the time series of available larval habitat, Figure 5c presents the resulting sporozoite rates, and Figure 5d the entomological inoculation rate, which represents the infectious bites received per person per night. If mortality increases as a function of age, the total population numbers do not change greatly, but the sporozoite rate drops due to the suppression of the older part of the mosquito age-distribution. Note that increasing larval habitat has a second order effect on EIR by allowing the human prevalence to rise earlier in the season, which increases the infectivity of the human population to vectors.

Figure 5
figure 5

Baseline population dynamics summed over local populations of An. gambiae ss, An. funestus , and An. arabiensis for different larval habitat multipliers. a) Local weather will drive both the temporary and semi-permanent larval habitats. b) The adult vector population changes as a function of the scaling of the larval habitat carrying capacity, which is driven by local weather. c) The sporozoite rate of mosquito population changes in response to the changing age structure of the vector population over the course of two years. d) Daily EIR combines the adult vector population and the sporozoite rate.

Further simulations demonstrate effects of combined vector control interventions such as insecticide-treated nets (ITN), indoor residual spraying (IRS), larval control, and space spraying. Figures 6 a, c, e show the changes in vector dynamics as coverage with perfect IRS is increased, killing all indoor feeding mosquitoes and maintained at full efficacy for the specified coverage without decay. In addition, all mosquitoes are set to feed indoors, making this an unrealistic scenario, but a useful boundary case showing the maximum possible effect. Larval habitat is set to 3.0 for all three species, and the simple human disease model is used. The detailed model outputs can be used to determine the change in entomologic inoculation rate as described in [9]. These results can be compared to field results for bed net campaigns in the presence of multiple vector species [66]. The size of the local vector population is reduced, but the effects on sporozoite rate and EIR are much more dramatic for several reasons, especially the restructuring of the mosquito population age distribution. The higher mortality results in fewer old mosquitoes in the population, which is the segment of the population with sporozoites. The older cohorts are thus responsible for the major portion of EIR but only minor portions of adult populations, human biting, and fecundity. In many cases, larval habitat remains the limiting factor in determining the number of emerging mosquitoes and interventions primarily act through adult mortality, but at high IRS or ITN coverage levels, it is possible in simulations to reduce emergence rate by limiting successful feeds. This phenomenon has been seen in high-coverage field studies [67], and the disappearance of An. funestus from parts of its earlier habitat as intervention coverage increases is the extreme limit of this phenomenon. The simulations are repeated for IRS with 0.6 killing of post-feeding mosquitoes and results are shown in Figures 6 b, d, f. The effects on adult vectors and sporozoite rate are reduced, and these reductions are compounded in the effect on EIR. Feeds in houses with IRS now have a 40 percent survival probability in contrast to the 0 percent survival in the previous simulations. Thus for 60 percent coverage, the survival for indoor feeds on the population is now 64 percent instead of 40 percent, and the probability of surviving three feeds rises to 26 percent from 6 percent.

Figure 6
figure 6

Effects of increasing coverage with perfect IRS. Effects of increasing coverage with perfect IRS (p kill, IRSpostfeed = 1) on a) Adult population c) Sporozoite rate e) EIR. Effects on sporozoite rate and EIR are much more dramatic than on the adult population because of the restructuring of the age distribution of the mosquito population. For most coverage levels, larval habitat remains the limiting factor in the rate of emerging mosquitoes, and number of young, unfed mosquitoes remains similar as IRS coverage increases up to a point. However, the increased feeding mortality results in a decreased life expectancy for mosquitoes, a moderate reduction in total population, but a strong reduction in mosquitoes older than 10 days. b, d, f) Repeated for IRS with (p kill, IRSpostfeed = 0.6). The effects on sporozoite rate and EIR are not as dramatic due to improved mosquito survival. The larval habitat multiplier is set to 3.0 for these simulations.

The effects of co-varying IRS and ITN coverage are studied and many simulations are run over sections of campaign and parameter space, with each simulation a trial for local elimination success or failure. These trials can be used to estimate the probability of local elimination as a function over campaign or parameter space. In Figure 7, coverage with IRS and ITN is varied, and the trials are assembled into plots, which map the regions of high probability of success and low probability of success. For purposes of this demonstration, both IRS and ITN kill every single relevant mosquito (pkill, ITN = 1 and pkill, IRSpostfeed = 1) and do not decay. All three mosquito species are simulated to feed indoors and take 95 percent of their blood meals on human hosts. The larval habitat scaling is set to 1.0, and the simple human disease model is used. Given these assumptions, it is not surprising that the region of success is large. Decreasing pkill, ITN and pkill, IRSpostfeed to 0.6 and maintaining bed nets at 100% blocking of indoor feeding produces the 2-D plot in Figure 7c. High coverage of bed nets still successfully locally eliminates the disease, since all feeds in this simulation occur indoors at night, but this is a quadratic effect of blocking part of the population from acquisition and transmission, while the exponential effect of increasing local mosquito mortality is drastically reduced with dramatic effects. Even with full IRS coverage and all mosquitoes feeding indoors, a reduced killing efficacy may not permit local elimination without supplemental interventions. The reduced efficacy is intended to represent an effect of insecticide resistance, and the effect on elimination has important implications for campaigns.

Figure 7
figure 7

Effects of combining IRS and ITN. a, b) Probability of Eradication and Estimate Variance for perfect bed nets and IRS, which do not decay, for fully indoor feeding and resting mosquitoes c, d) Bed nets still prevent all nighttime feeds, but only kill 60% of mosquitoes attempting to feed. IRS kills 60% of post-indoor feeding mosquitoes in treated houses. Eradication is no longer possible in many previously possible parameter regions. e, f) Effect of varying indoor feeding and anthropophily of the An. arabiensis population for 90 percent IRS coverage with no decay of insecticide and p kill, IRSpostfeed = 1. The multiplier for larval habitat set to 1.0 for all three sets of simulations, and increasing larval habitat increases the level of coverage required, but not as dramatically as changing the adult mortality. Dark blue regions in a, c, e correspond to parameter regimes in which the estimated probability of local elimination is over 0.9 and dark red less than 0.1. Level sets for mean estimated probability of elimination and for probability estimate variance are labeled.

Changes in vector behaviours such as indoor feeding and anthropophily change the results drastically with several key implications for eradication. The previous system is rerun with IRS at 90 percent coverage with full mosquito post-feeding mortality for indoor feeds in treated-houses, and the human feeding and indoor feeding preferences of the local arabiensis population are varied. Changes in indoor feeding have a dramatic effect on campaign success as seen in Figure 7 e, f, as would be expected for interventions which only affect indoor feeding mosquitoes. Decreasing anthropophily is often viewed as decreasing transmission, but in the presence of high intervention coverage, decreasing anthropophily reduces the mortality during sporogony, allowing a higher proportion of infected mosquitoes to complete sporogony and thereby reducing the probability of campaign success. Feeds on animals during sporogony are safe compared to the protected human feeds and increase the probability of surviving the multiple feeds during sporogony and becoming infectious. At lower anthropophily, fewer mosquitoes become infected and infectious mosquitoes bite humans less frequently, and this effect eventually wins and probability of success increases with decreasing anthropophily below a certain value. Detailed model representations of vector behaviour help explain the failure of elimination campaigns which only targeted indoor feeding. For eradication to succeed, the full transmission cycle must be sufficiently broken, which may involve targeting outdoor-feeding mosquitoes in some areas.

Other available but currently less-used options for vector control include larvicides and space spraying for targeting larval and adult populations, respectively. Figure 8 shows results for simulations of larval control on the left and adult-targeting space spraying on the right. In the presented simulations, larval control is simulated as a temporary decrease in the larval habitat carrying capacity for a specified duration, such as a 30 percent reduction for 180 days. This produces decreases in the adult population, but not in the adult age-structure-driven sporozoite rate, except for transients at the start and conclusion of larval control. Thus the effects on EIR tend to be linear, but this should not be ignored as larval control may be one of the only ways to target outdoor feeding mosquitoes.

Figure 8
figure 8

Larval control and space spraying. Larval control and space spraying introduced to the baseline simulations from Figure 5, with larval habitat multipliers set to 3.0. a, c, e) Larval control on the left reduces adult vector populations and EIR, but does not affect sporozoite rate at these high EIRs as the age structure of the vector population is unaffected. b, d, f) Adult mortality affects all three measures, as it reduces the number of adult vectors, and dramatically changes the age structure so that fewer mosquitoes are old enough to have sporozoites. This creates a compounded effect on EIR.

In order to achieve local elimination of malaria in areas with high rates of outdoor feeding by the local vector populations, some form of control of outdoor mosquitoes may be necessary. Available options can be costly and logistically difficult, and studies of their effects can place constraints on required target efficacy, duration, and frequency of such efforts. In the model, space spraying increases the mortality for all adult mosquitoes, regardless of the stage in the feeding cycle. The artificial daily mortality probability is calculated as 1-exp(-killrate*Δt), so that a killrate of 1.0 will tend to kill approximately 63 percent of adult vectors. As seen in Figure 8, even one day of spraying with a high killrate can produce a large several week drop in EIR as it takes time for the newly emerging mosquitoes to become infected and infectious. Longer duration efforts, such as repeated spraying with a given knockdown of adult vectors each cycle, increase the daily mortality and reduce the number of adult vectors, reshape the age structure and reduce the sporozoite rate and achieve a strong multiplicative decrease on the EIR. Such maintained repetitive space spraying with non-residual insecticides may be logistically difficult, but it could provide leverage on the outdoor feeding population if other options fall through. Simulations such as these can estimate the effects on vector populations for a given input efficacy, duration, and frequency of application, providing inputs to cost-effectiveness analyses.

Conclusions

The present model creates a flexible framework for exploring the effects of combined vector control interventions on vector population dynamics and disease transmission. Campaigns with IRS and ITN are studied and success of elimination is seen to depend on coverage and efficacy as expected, but also on mosquito behaviour. As An. arabiensis outdoor feeding increases, interventions which target indoor feeding become less effective. Decreasing anthropophily from unity at high coverage decreases the rate of killing mosquitoes during sporogony, initially reducing elimination success, but this returns to the expected relationship of less anthropophily increasing elimination success as anthropophily continues to decline to low levels. Climate and weather data with high spatial resolution can help predict spatial and temporal patterns of vector dynamics and assist the rational planning of regional campaigns, especially when included combined with a population transmission model [16, 17, 19].

Further work will exhibit the effect of vector migration and seasonality on interventions such as locally-applied larvicides. Understanding the role and scale of migration is important for estimating effect of adult vector interventions [61, 67] and larval control interventions [9, 59, 63, 64]. The model supports spatially-distributed simulations, and future work will explore the effects of human and vector migration on spatial transmission. Each individual in the simulation has a relative biting rate, which allows study of heterogeneous biting as has been done [37], and link to studies of attractiveness of humans to mosquitoes versus level of malaria infection [68]. Improvements will be made to the density-dependence effects of larval dynamics and models for habitat calculations will be improved. Other future work will model the response of vector populations to applied pressure, such as the lower anthropophily of An. gambiae ss in The Gambia after extended bed net pressure [32]. Changes in behavior or development of insecticide resistance require careful consideration to ensure success of a designed Eradication campaign. The effects of such changes can be seen in the above-presented results, and the present model has the flexibility to incorporate dynamic changes. The modular structure of the model and the implementation of the life and feeding cycles also make it simple to add new potential interventions to the model. Other species are easy to add provided that an applicable habitat model has been developed. Future work can couple this detailed vector transmission model to a more detailed model of human disease and immunity.

Campaigns must address the ecology and behaviour of local mosquito populations in order to ensure that sufficient resources with broad enough effects for all relevant components of the local mosquito populations are introduced. A one-size-fits-all campaign is not optimal, being wasteful in some circumstances and insufficient in others; local tailoring and design are important. Modelling can be used to estimate the risk of disease transmission given reintroduction to areas that had achieved local elimination before their neighbours [6]. Modelling at this level of detail also serves to identify basic data gaps such as local vector ecology and behaviour that must be answered to reduce uncertainty of campaign success. Numerical studies can reveal to which parameters the results of interest are most sensitive, and such parameters which are also poorly constrained by data or are highly geographically-variable can then be highlighted as important data gaps. Modelling studies can also explore the extent to which more through and extensive campaigns can overwhelm data uncertainties and achieve more robust success.

References

  1. Snow RW, Guerra CA, Noor AM, Myint HY, Hay SI: The global distribution of clinical episodes of Plasmodium falciparum malaria. Nature. 2005, 434: 214-217. 10.1038/nature03342.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  2. Jacob BG, Muturi E, Halbig P, Mwangangi J, Wanjogu RK, Mpanga E, Funes J, Shililu J, Githure J, Regens JL, Novak RJ: Environmental abundance of Anopheles (Diptera: Culicidae) larval habitats on land cover change sites in Karima Village, Mwea Rice Scheme, Kenya. Am J Trop Med Hyg. 2007, 76: 73-80.

    PubMed  Google Scholar 

  3. De Castro MC, Yamagata Y, Mtasiwa D, Tanner M, Utzinger J, Keiser J, Singer BH: Integrated urban malaria control: a case study in Dar Es Salaam, Tanzania. Am J Trop Med Hyg. 2004, 71: 103-117.

    Google Scholar 

  4. Craig MH, Snow RW, le Sueur D: A climate-based distribution model of malaria transmission in sub-Saharan Africa. Parasitol Today. 1999, 15: 105-111. 10.1016/S0169-4758(99)01396-4.

    Article  CAS  PubMed  Google Scholar 

  5. Hay SI, Guerra CA, Gething PW, Patil AP, Tatem AJ, Noor AM, Kabaria CW, Manh BH, Elyazar IR, Brooker S, Smith DL, Moyeed RA, Snow RW: A world malaria map: Plasmodium falciparum endemicity in 2007. PLoS Med. 2009, 6: e1000048-

    Article  PubMed Central  PubMed  Google Scholar 

  6. Cohen JM, Smith DL, Vallely A, Taleo G, Malefoasi G, Sabot O: Holding the Line. Shrinking the Malaria Map: A Prospectus on Malaria Elimination. Edited by: Feachem RGA, Phillips AA, Targett GA, Group ME. 2009, San Francisco: The Global Health Group: UCSF Global Health Sciences, 40-60.

    Google Scholar 

  7. Ross R: The prevention of malaria. 1910, New York,: E.P. Dutton & company

    Google Scholar 

  8. Macdonald G: The epidemiology and control of malaria. 1957, London, New York,: Oxford University Press

    Google Scholar 

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

    Article  PubMed Central  PubMed  Google Scholar 

  10. 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 Applied Ecol. 1990, 27: 123-133. 10.2307/2403572.

    Article  Google Scholar 

  11. Killeen G, McKenzie F, Foy B, Schieffelin C, Billingsley P, Beier J: 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 

  12. Yazoume Y, Hoshen M, Kyobutungi C, Louis VR, Sauerborn R: Local scale prediction of Plasmodium falciparum malaria transmission in an endemic region using temperature and rainfall. Global Health Action. 2009, 2:

    Google Scholar 

  13. Killeen GF, Ross A, Smith T: Infectiousness of malaria-endemic human populations to vectors. Am J Trop Med Hyg. 2006, 75: 38-45.

    PubMed  Google Scholar 

  14. Bomblies A, Duchemin J-B, Eltahir EAB: Hydrology of malaria: Model development and application to a Sahelian village. Water Resour Res. 2008, 44: W12445-

    Google Scholar 

  15. Depinay J-M, Mbogo C, Killeen G, Knols B, Beier J, Carlson J, Dushoff J, Billingsley P, Mwambi H, Githure J, Toure AM, McKenzie FE: A simulation model of African Anopheles ecology and population dynamics for the analysis of malaria transmission. Malar J. 2004, 3: 29-10.1186/1475-2875-3-29.

    Article  PubMed Central  PubMed  Google Scholar 

  16. Filipe JA, Riley EM, Drakeley CJ, Sutherland CJ, Ghani AC: Determination of the processes driving the acquisition of immunity to malaria using a mathematical transmission model. PLoS Comput Biol. 2007, 3: e255-10.1371/journal.pcbi.0030255.

    Article  PubMed Central  PubMed  Google Scholar 

  17. Griffin JT, Hollingsworth TD, Okell LC, Churcher TS, White M, Hinsley W, Bousema T, Drakeley CJ, Ferguson NM, Basáñez M-G, Ghani AC: Reducing Plasmodium falciparum malaria transmission in Africa: a model-based evaluation of intervention strategies. PLoS Med. 2010, 7: e1000324-10.1371/journal.pmed.1000324.

    Article  PubMed Central  PubMed  Google Scholar 

  18. Smith T, Maire N, Ross A, Penny M, Chitnis N, Schapira A, Studer A, Genton B, Lengeler C, Tediosi F, de Savigny D, Tanner M: Towards a comprehensive simulation model of malaria epidemiology and control. Parasitology. 2008, 135: 1507-1516. 10.1017/S0031182008000371.

    Article  CAS  PubMed  Google Scholar 

  19. Smith T, Killeen GF, Maire N, Ross A, Molineaux L, Tediosi F, Hutton G, Utzinger J, Dietz K, Tanner M: Mathematical modeling of the impact of malaria vaccines on the clinical epidemiology and natural history of Plasmodium falciparum malaria: Overview. Am J Trop Med Hyg. 2006, 75: 1-10.

    PubMed  Google Scholar 

  20. Magori K, Legros M, Puente ME, Focks DA, Scott TW, Lloyd AL, Gould F: Skeeter Buster: a stochastic, spatially explicit modeling tool for studying Aedes aegypti population replacement and population suppression strategies. PLoS Negl Trop Dis. 2009, 3: e508-10.1371/journal.pntd.0000508.

    Article  PubMed Central  PubMed  Google Scholar 

  21. Simard F, Ayala D, Kamdem GC, Pombi M, Etouna J, Ose K, Fotsing JM, Fontenille D, Besansky NJ, Costantini C: Ecological niche partitioning between Anopheles gambiae molecular forms in Cameroon: the ecological side of speciation. BMC Ecol. 2009, 9: 17-10.1186/1472-6785-9-17.

    Article  PubMed Central  PubMed  Google Scholar 

  22. Fillinger U, Sonye G, Killeen GF, Knols BGJ, Becker N: The practical importance of permanent and semipermanent habitats for controlling aquatic stages of Anopheles gambiae sensu lato mosquitoes: operational observations from a rural town in western Kenya. Trop Med Int Health. 2004, 9: 1274-1289. 10.1111/j.1365-3156.2004.01335.x.

    Article  PubMed  Google Scholar 

  23. Koenraadt CJ, Githeko AK, Takken W: The effects of rainfall and evapotranspiration on the temporal dynamics of Anopheles gambiae s.s. and Anopheles arabiensis in a Kenyan village. Acta Trop. 2004, 90: 141-153. 10.1016/j.actatropica.2003.11.007.

    Article  CAS  PubMed  Google Scholar 

  24. Bomblies A, Duchemin J-B, Eltahir E: A mechanistic approach for accurate simulation of village scale malaria transmission. Malar J. 2009, 8: 223-10.1186/1475-2875-8-223.

    Article  PubMed Central  PubMed  Google Scholar 

  25. Shaman J, Spiegelman M, Cane M, Stieglitz M: A hydrologically driven model of swamp water mosquito population dynamics. Ecological Modelling. 2006, 194: 395-404. 10.1016/j.ecolmodel.2005.10.037.

    Article  Google Scholar 

  26. Guelbeogo W, Sagnon NF, Grushko O, Yameogo M, Boccolini D, Besansky N, Costantini C: Seasonal distribution of Anopheles funestus chromosomal forms from Burkina Faso. Malar J. 2009, 8: 239-10.1186/1475-2875-8-239.

    Article  PubMed Central  PubMed  Google Scholar 

  27. Gillies MT, De Meillon B: The Anophelinae of Africa, South of the Sahara. 1968, Johannesburg, South Africa: South African Institute for Medical Research

    Google Scholar 

  28. Paaijmans KP, Wandago MO, Githeko AK, Takken W: Unexpected high losses of Anopheles gambiae larvae due to rainfall. PLoS One. 2007, 2: e1146-10.1371/journal.pone.0001146.

    Article  PubMed Central  PubMed  Google Scholar 

  29. Mutuku FM, Bayoh MN, Gimnig JE, Vulule JM, Kamau L, Walker ED, Kabiru E, Hawley WA: Pupal habitat productivity of Anopheles gambiae complex mosquitoes in a rural village in western Kenya. Am J Trop Med Hyg. 2006, 74: 54-61.

    PubMed  Google Scholar 

  30. Clements AN: The Biology of Mosquitoes, Volume 1: Development, Nutrition, and Reproduction. 2000, New York: CABI Publishing

    Google Scholar 

  31. Dobson SL, Fox CW, Jiggins FM: The effect of Wolbachia-induced cytoplasmic incompatibility on host population size in natural and manipulated systems. Proc R Soc London Series B: Biol Sciences. 2002, 269: 437-445. 10.1098/rspb.2001.1876.

    Article  Google Scholar 

  32. Killeen GF, McKenzie FE, Foy BD, Bogh C, Beier JC: The availability of potential hosts as a determinant of feeding behaviours and malaria transmission by African mosquito populations. Trans R Soc Trop Med Hyg. 2001, 95: 469-476. 10.1016/S0035-9203(01)90005-7.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. Saul A: Zooprophylaxis or zoopotentiation: the outcome of introducing animals on vector transmission is highly dependent on the mosquito mortality while searching. Malar J. 2003, 2: 32-10.1186/1475-2875-2-32.

    Article  PubMed Central  PubMed  Google Scholar 

  34. Clements AN, Paterson GD: The analysis of mortality and survival rates in wild populations of mosquitoes. J Applied Ecol. 1981, 18: 373-399. 10.2307/2402401.

    Article  Google Scholar 

  35. Jeffery GM, Eyles DE: Infectivity to mosquitoes of Plasmodium falciparum as related to gametocyte density and duration of infection. Am J Trop Med Hyg. 1955, 4: 781-789.

    CAS  PubMed  Google Scholar 

  36. Sinden RE, Butcher GA, Billker O, Fleck SL: Regulation of infectivity of Plasmodium to the mosquito vector. Adv Parasitol. 1996, 38: 53-117.

    Article  CAS  PubMed  Google Scholar 

  37. Smith DL, McKenzie FE, Snow RW, Hay SI: Revisiting the basic reproductive number for malaria and its implications for malaria control. PLoS Biol. 2007, 5: e42-10.1371/journal.pbio.0050042.

    Article  PubMed Central  PubMed  Google Scholar 

  38. Githeko AK, Brandling-Bennett AD, Beier M, Atieli F, Owaga M, Collins FH: The reservoir of Plasmodium falciparum malaria in a holoendemic area of western Kenya. Trans R Soc Trop Med Hyg. 1992, 86: 355-358. 10.1016/0035-9203(92)90216-Y.

    Article  CAS  PubMed  Google Scholar 

  39. Bonnet S, Gouagna C, Safeukui I, Meunier J-Y, Boudin C: Comparison of artificial membrane feeding with direct skin feeding to estimate infectiousness of Plasmodium falciparum gametocyte carriers to mosquitoes. Trans R Soc Trop Med Hyg. 2000, 94: 103-106. 10.1016/S0035-9203(00)90456-5.

    Article  CAS  PubMed  Google Scholar 

  40. Smalley ME, Sinden RE: Plasmodium falciparum gametocytes: their longevity and infectivity. Parasitology. 1977, 74: 1-8. 10.1017/S0031182000047478.

    Article  CAS  PubMed  Google Scholar 

  41. Schneider P, Bousema JT, Gouagna LC, Otieno S, van de Vegte-Bolmer M, Omar SA, Sauerwein RW: Submicroscopic Plasmodium falciparum gametocyte densities frequently result in mosquito infection. Am J Trop Med Hyg. 2007, 76: 470-474.

    PubMed  Google Scholar 

  42. Muirhead-Thomson RC: The malarial infectivity of an African village population to mosquitoes (Anopheles gambiae); a random xenodiagnostic survey. Am J Trop Med Hyg. 1957, 6: 971-979.

    CAS  PubMed  Google Scholar 

  43. Drakeley C, Sutherland C, Bousema JT, Sauerwein RW, Targett GAT: The epidemiology of Plasmodium falciparum gametocytes: weapons of mass dispersion. Trends Parasitol. 2006, 22: 424-430. 10.1016/j.pt.2006.07.001.

    Article  PubMed  Google Scholar 

  44. Bonnet S, Gouagna LC, Paul RE, Safeukui I, Meunier JY, Boudin C: Estimation of malaria transmission from humans to mosquitoes in two neighbouring villages in south Cameroon: evaluation and comparison of several indices. Trans R Soc Trop Med Hyg. 2003, 97: 53-59. 10.1016/S0035-9203(03)90022-8.

    Article  CAS  PubMed  Google Scholar 

  45. Naotunne TS, Karunaweera ND, Del Giudice G, Kularatne MU, Grau GE, Carter R, Mendis KN: Cytokines kill malaria parasites during infection crisis: extracellular complementary factors are essential. J Exp Med. 1991, 173: 523-529. 10.1084/jem.173.3.523.

    Article  CAS  PubMed  Google Scholar 

  46. Naotunne TS, Karunaweera ND, Mendis KN, Carter R: Cytokine-mediated inactivation of malarial gametocytes is dependent on the presence of white blood cells and involves reactive nitrogen intermediates. Immunology. 1993, 78: 555-562.

    PubMed Central  CAS  PubMed  Google Scholar 

  47. Okell L, Drakeley C, Ghani A, Bousema T, Sutherland C: Reduction of transmission from malaria patients by artemisinin combination therapies: a pooled analysis of six randomized trials. Malar J. 2008, 7: 125-10.1186/1475-2875-7-125.

    Article  PubMed Central  PubMed  Google Scholar 

  48. Al-Olayan EM, Beetsma AL, Butcher GA, Sinden RE, Hurd H: Complete Development of mosquito phases of the malaria parasite in vitro. Science. 2002, 295: 677-679. 10.1126/science.1067159.

    Article  CAS  PubMed  Google Scholar 

  49. Beier JC: Malaria parasite development in mosquitoes. Ann Rev Entomol. 1998, 43: 519-543. 10.1146/annurev.ento.43.1.519.

    Article  CAS  Google Scholar 

  50. Gorman MJ, Cornel AJ, Collins FH, Paskewitz SM: A shared genetic mechanism for melanotic encapsulation of CM-Sephadex beads and a malaria parasite, Plasmodium cynomolgi B, in the mosquito, Anopheles gambiae. Exp Parasitol. 1996, 84: 380-386. 10.1006/expr.1996.0126.

    Article  CAS  PubMed  Google Scholar 

  51. Macdonald G: Theory of the eradication of malaria. Bull World Health Organ. 1956, 15: 369-387.

    PubMed Central  CAS  PubMed  Google Scholar 

  52. Afrane YA, Little TJ, Lawson BW, Githeko AK, Yan G: Deforestation and vectorial capacity of Anopheles gambiae Giles mosquitoes in malaria transmission, Kenya. Emerg Infect Dis. 2008, 14: 1533-1538. 10.3201/eid1410.070781.

    Article  PubMed Central  PubMed  Google Scholar 

  53. Burkot TR, Graves PM, Paru R, Battistutta D, Barnes A, Saul A: Variations in malaria transmission rates are not related to anopheline survivorship per feeding cycle. Am J Trop Med Hyg. 1990, 43: 321-327.

    CAS  PubMed  Google Scholar 

  54. Wekesa JW, Copeland RS, Mwangi RW: Effect of Plasmodium falciparum on blood feeding behavior of naturally infected anopheles mosquitoes in western Kenya. Am J Trop Med Hyg. 1992, 47: 484-488.

    CAS  PubMed  Google Scholar 

  55. Anderson RA, Knols BGJ, Kolella JC: Plasmodium falciparum sporozoites increase feeding-associated mortality of their mosquito hosts Anopheles gambiae s.l. Parasitology. 2000, 120: 329-333. 10.1017/S0031182099005570.

    Article  PubMed  Google Scholar 

  56. Carwardine SL, Hurd H: Effects of Plasmodium yoelii nigeriensis infection on Anopheles stephensi egg development and resorption. Med Vet Entomol. 1997, 11: 265-269. 10.1111/j.1365-2915.1997.tb00405.x.

    Article  CAS  PubMed  Google Scholar 

  57. Hogg JC, Hurd H: The effects of natural Plasmodium falciparum infection on the fecundity and mortality of Anopheles gambiae s. l. in north east Tanzania. Parasitology. 1997, 114: 325-331. 10.1017/S0031182096008542.

    Article  PubMed  Google Scholar 

  58. Smith DL, Dushoff J, Snow RW, Hay SI: The entomological inoculation rate and Plasmodium falciparum infection in African children. Nature. 2005, 438: 492-495. 10.1038/nature04024.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  59. Fillinger U, Ndenga B, Githeko A, Lindsay SW: Integrated malaria vector control with microbial larvicides and insecticide-treated nets in western Kenya: a controlled trial. Bull World Health Organ. 2009, 87: 655-665. 10.2471/BLT.08.055632. vvv

    Article  PubMed Central  PubMed  Google Scholar 

  60. Killeen GF, Smith TA, Ferguson HM, Mshinda H, Abdulla S, Lengeler C, Kachur SP: Preventing childhood malaria in Africa by protecting adults from mosquitoes with insecticide-treated nets. PLoS Med. 2007, 4: e229-10.1371/journal.pmed.0040229.

    Article  PubMed Central  PubMed  Google Scholar 

  61. Gimnig JE, Vulule JM, Lo TQ, Kamau L, Kolczak MS, Phillips-Howard PA, Mathenge EM, ter Kuile FO, Nahlen BL, Hightower AW, Hawley WA: Impact of permethrin-treated bed nets on entomologic indices in an area of intense year-round malaria transmission. Am J Trop Med Hyg. 2003, 68: 16-22.

    PubMed  Google Scholar 

  62. Snowden F: The conquest of malaria: Italy, 1900-1962. 2006, New Haven (CT): Yale University Press

    Book  Google Scholar 

  63. Fillinger U, Kannady K, William G, Vanek M, Dongus S, Nyika D, Geissbuhler Y, Chaki P, Govella N, Mathenge E: A tool box for operational mosquito larval control: preliminary results and early lessons from the Urban Malaria Control Programme in Dar es Salaam, Tanzania. Malar J. 2008, 7: 20-10.1186/1475-2875-7-20.

    Article  PubMed Central  PubMed  Google Scholar 

  64. Killeen G, Fillinger U, Knols B: Advantages of larval control for African malaria vectors: Low mobility and behavioural responsiveness of immature mosquito stages allow high effective coverage. Malar J. 2002, 1: 8-10.1186/1475-2875-1-8.

    Article  PubMed Central  PubMed  Google Scholar 

  65. Full Data Reanalysis. [http://gpcc.dwd.de]

  66. Killeen G, McKenzie F, Foy B, Schieffelin C, Billingsley P, Beier J: The potential impact of integrated malaria transmission control on entomologic inoculation rate in highly endemic areas. Am J Trop Med Hyg. 2000, 62: 545-551.

    PubMed Central  CAS  PubMed  Google Scholar 

  67. Gimnig JE, Kolczak MS, Hightower AW, Vulule JM, Schoute E, Kamau L, Phillips-Howard PA, Ter Kuile FO, Nahlen BL, Hawley WA: Effect of permethrin-treated bed nets on the spatial distribution of malaria vectors in western Kenya. Am J Trop Med Hyg. 2003, 68: 115-120.

    PubMed  Google Scholar 

  68. Lacroix R, Mukabana WR, Gouagna LC, Koella JC: Malaria infection increases attractiveness of humans to mosquitoes. PLoS Biol. 2005, 3: e298-10.1371/journal.pbio.0030298.

    Article  PubMed Central  PubMed  Google Scholar 

  69. Kiszewski A, Mellinger A, Spielman A, Malaney P, Sachs SE, Sachs J: A global index representing the stability of malaria transmission. Am J Trop Med Hyg. 2004, 70: 486-498.

    PubMed  Google Scholar 

Download references

Acknowledgements

The author thanks Bill and Melinda Gates for their active support of this work and their sponsorship through the Global Good Fund. This work was performed at Intellectual Ventures Laboratories, and useful discussions with colleagues in this Program, the Malaria Program of the Bill and Melinda Gates Foundation, MalERA, and the Vector Control Development Network are likewise greatly appreciated.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Philip A Eckhoff.

Additional information

Competing interests

The authors declare that they have no competing interests.

Electronic supplementary material

12936_2011_1903_MOESM1_ESM.DOC

Additional file 1: Detailed equations for A malaria transmission-directed model of mosquito life cycle and ecology. The detailed calculations of vector feeding outcomes for feeds on a single individual and on the full local human population are provided. (DOC 62 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 (http://creativecommons.org/licenses/by/2.0), 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

Eckhoff, P.A. A malaria transmission-directed model of mosquito life cycle and ecology. Malar J 10, 303 (2011). https://doi.org/10.1186/1475-2875-10-303

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1475-2875-10-303

Keywords