Skip to main content

Comparative assessment of diverse strategies for malaria vector population control based on measured rates at which mosquitoes utilize targeted resource subsets

Abstract

Background

Eliminating malaria requires vector control interventions that dramatically reduce adult mosquito population densities and survival rates. Indoor applications of insecticidal nets and sprays are effective against an important minority of mosquito species that rely heavily upon human blood and habitations for survival. However, complementary approaches are needed to tackle a broader diversity of less human-specialized vectors by killing them at other resource targets.

Methods

Impacts of strategies that target insecticides to humans or animals can be rationalized in terms of biological coverage of blood resources, quantified as proportional coverage of all blood resources mosquito vectors utilize. Here, this concept is adapted to enable impact prediction for diverse vector control strategies based on measurements of utilization rate s for any definable, targetable resource subset, even if that overall resource is not quantifiable.

Results

The usefulness of this approach is illustrated by deriving utilization rate estimates for various blood, resting site, and sugar resource subsets from existing entomological survey data. Reported impacts of insecticidal nets upon human-feeding vectors, and insecticide-treated livestock upon animal-feeding vectors, are approximately consistent with model predictions based on measured utilization rates for those human and animal blood resource subsets. Utilization rates for artificial sugar baits compare well with blood resources, and are consistent with observed impact when insecticide is added. While existing data was used to indirectly measure utilization rates for a variety of resting site subsets, by comparison with measured rates of blood resource utilization in the same settings, current techniques for capturing resting mosquitoes underestimate this quantity, and reliance upon complex models with numerous input parameters may limit the applicability of this approach.

Conclusions

While blood and sugar consumption can be readily quantified using existing methods for detecting natural markers or artificial tracers, improved techniques for labelling mosquitoes, or other arthropod pathogen vectors, will be required to assess vector control measures which target them when they utilize non-nutritional resources such as resting, oviposition, and mating sites.

Background

While antiparasitic drugs and vaccines will be essential to the final stages of malaria elimination, their effectiveness as transmission control interventions will rely heavily upon unprecedented levels of vector control in highly endemic settings[1–4]. It will not be possible to eliminate malaria transmission from most of the tropics without developing scalable vector control intervention options which complement long-lasting insecticidal nets (LLINs) and indoor residual spraying (IRS) by targeting adult mosquitoes when they use resources other than human blood indoors, and indoor resting sites[5–9]. Mosquitoes usually also exploit non-human blood and human blood sources outdoors, as well as sugar, outdoor resting sites, oviposition sites, and mating sites, so all of these other biological and environmental resources represent alternative targets for vector control interventions[5].

Faced with such an array of resource target options, the challenge is to define exactly which of these intervention targets are optimal in each of the diverse vectorial systems that exist[10], and to classify these settings into limited numbers of distinct categories where specific intervention combinations maximize impact. Recent analyses indicate that the impact of vector control measures targeting the blood hosts upon which mosquitoes depend can be rationalized in terms of measurements of the biological coverage of all available blood sources that is achieved, rather than merely high demographic coverage of a targeted subset, such as humans while asleep indoors[9, 11]. Blood resources are perhaps the best understood, and most readily quantified, of all the resources used by mosquitoes that could potentially be targeted with vector control measures[5]. However, many established or emerging vector control strategies only target a specific subset of the resource that can be readily identified and treated in the field. Examples include targeting insecticides to indoor resting sites only with IRS[12], and artificially introduced resting sites[13–17] or sugar baits[18–20]. Furthermore, the other available forms of these resources that cannot be targeted with insecticides but compete with these subsets for the attentions of mosquitoes, such as naturally occurring sugar sources[21] or outdoor resting sites[22, 23], are often impossible to identify or quantify[5]. It is therefore not possible to estimate biological coverage as a fraction of all available forms of that resource. Here, the concept of biological coverage is extended beyond blood resources, and adapted to enable impact prediction for more diverse vector control strategies, based on direct measurements of coverage and utilization rates for definable, targetable subsets of less readily quantified resources that are equally important to mosquito survival, and therefore equally valid as potential targets for vector control interventions.

Methods

All symbols used are listed and defined in Table 1.

Table 1 Parameter symbols and definitions

Defining biological coverage based on the example of blood resources

Biological coverage (C A,p ) of all the blood resources upon which mosquitoes rely, with long-lasting insecticidal nets (LLINs) or any other personal protection measure, has been previously defined as the proportion of all mosquito attacks upon all available hosts for which those hosts were covered with a protective (p) intervention at that time and place[11, 24]:

C A , p = A p A
(1)

where the total attack availabilities of the all hosts (A), and covered hosts at times and places when they are actually protected (A p ), are defined kinetically[25, 42] as the rates per night at which an individual host-seeking mosquito respectively encounters and attacks[26] either all hosts or all hosts that are protected at the time of the encounter and attack events.

However, to allow simplified notation for generalization of this approach to a greater diversity of distinct resources, here the symbol A is replaced by R to denote the total kinetic availability rate of a specific given resource, which may be specified as blood (v), resting sites (r), sugar (s) or any other resource (R ∈ v, r, s …). Also, the terms attack and protected which were previously used to define availability in kinetic terms for models of blood resource seeking and acquisition[26, 27] are not entirely appropriate for non-blood resources, so these are replaced with more generally applicable terms attempt to utilize and covered, respectively. Furthermore, now that biological coverage has been defined to explicitly consider only protection that is in place at the times and places when that resource is utilized, de facto coverage (c) and protection (p) are equivalent to each other (c ≈ p), so the former is used to simplify and harmonize the notation. Expressing Equation 1 in terms of this revised notation yields:

C v = v c v
(2)

where C v is the proportion of all mosquito attacks upon real (live vertebrate hosts) or perceived (artificial odor-baited traps, sometimes referred to as pseudo-hosts[27]) blood resources to which effective coverage with a vector control intervention applies at that time and place, v is the total rate at which individual mosquitoes encounter and attack all hosts and pseudo hosts, and v c is the total rate at which individual mosquitoes encounter and attack all hosts and pseudo hosts at times and places when they are effectively covered with a vector control intervention.

In the case of interventions such as LLINs, which only protect humans while they use them indoors, biological coverage can be calculated as the product of the proportion of all bloodmeals (v) that originate from the human (h) host species subset Q v h , the proportion of human exposure to mosquito bites that would otherwise occur indoors (i) without an LLIN π v h , i , and the proportional demographic coverage of humans, measured as the proportion of humans who reported using a net while asleep indoors the previous night C v h , i [11]:

C v = Q v h π v h , i C v h , i
(3)

where all three terms on the right hand side of Equation 3 are defined as sequentially nested fractions and sub-fractions of the total availability of all blood hosts (v) that are represented by humans (h), those humans while indoors (i), and those humans while indoors and protected by coverage with LLIN use at the time (c):

Q v h = v h /v
(4)
Ï€ v h , i = v h , i / v h
(5)
C v h , i = v h , i , c / v h , i
(6)

Adapting the concept of biological coverage to rationalize vector control impact based on utilization rates of diverse resource targets

Expressing Equation 2 in more general terms that may be applied to any given resource (R), rather than blood specifically (v), yields the following formula:

C R = R c R
(7)

LLINs that directly kill mosquitoes when they encounter and attack protected human blood sources are the best established[43] and easiest resource targets to conceptualize and model, but previous formulations predicted their impact upon vector survival by assuming mosquitoes feed once and only once per gonotrophic cycle[26]. However, resting sites, oviposition sites, and even blood resources themselves, may be utilized more than once per gonotrophic cycle[44], while sugar sources may be used less than once[21, 45]. To adapt the concept of biological coverage to more diverse resource targets which are used more than once per gonotrophic cycle, the term resource utilization is defined as the mean rate at which mosquitoes utilize any given resource (R) per gonotrophic cycle (α R ). This definition of resource utilization rate can be expressed mathematically as the product of the duration of the gonotrophic cycle, expressed as nights per gonotrophic cycle (g), and the rate per night at which a mosquito population utilizes that resource (m R ), divided by the size of the mosquito population (M):

α R = g m R M
(8)

where g is expressed as units of nights, m R in units of utilization attempt events per night, and M as the number of individual adult mosquitoes present in the population. For any targetable, intervention-covered (c) proportion of that resource (R c ), the corresponding rate at which mosquitoes encounter and attempt to utilize that covered fraction α R c , by definition, varies in proportion to the fraction of the kinetic availability of that resource that it represents:

α R c = g m R c M = g m R M R c R
(9)

Hence the quotient of the availability or utilization rates for the total resource, divided into those for the intervention-covered fraction, are equivalent to biological coverage of that resource:

α R c α R = m R c m R = R c R = C R
(10)

Most vector control strategies only target a specific subset (x) of the resource that they are delivered to, which is practically definable, identifiable, accessible, and treatable in the field. Similarly to resource coverage (Equation 10), the proportion of all available forms of a specific resource (R) that is accounted for by a given subset (R x ) of that resource Q R x , can be defined and measured in terms of the rate at which mosquitoes encounter and attempt to utilize it[26, 27, 42] by generalizing Equation 4 for subsets of any possible resource, rather than specifying blood:

Q R x = R x R = α R x α R
(11)

Similarly, the proportion of that subset (x) that is effectively protected at times when it is utilized by mosquitoes (y) can be expressed in terms of the proportion of resource utilization attempt events it accounts for in that resource subset:

π R x , y = α R x , y / α R x
(12)

Hence, the biological coverage of all forms of that resource (C R ) can be expressed more explicitly than in Equation 10 as the product of the proportion of that resource represented by that subset Q R x , the proportion of utilization attempt events for that subset to which protection effectively and conditionally applies Ï€ R x , y and measured intervention coverage of that resource subset at times and places when it may be effectively protected C R x , y :

C R = Q R x π R x , y C R x , y = α R x , y α R C R x , y = α R x , y , c α R
(13)

Note that α R x , y , c is the utilization rate of the covered fraction of the targeted resource subset, equivalent to α R c in Equation 10, because all covered forms of the resource occur amongst the intervention-targeted subset of that resource (x) at the times and places at which they mosquitoes actually attempt to utilize them (y):

α R x , y , c = α R c because R c ∈ R x , y
(14)

Predicting intervention impact based on resource subset coverage and utilization rates

Interventions targeting adult mosquitoes may have quite complex modes of action, repelling mosquitoes away from humans[46] or contaminating them with agents that affect their longevity[47, 48], competence[47, 48] or fecundity[49, 50]. Biological agents may be transmitted horizontally or vertically through the population[47, 48, 50], while coverage amplification of chemicals may be achieved by mosquito-mediated transfer between resources[49]. Regardless of their complexity, these diverse strategies can all be enhanced by maximizing biological coverage of the resource targeted to ensure maximum contact with the mosquito population, and this is a critically important determinant of success in its own right. Previous formulations describing mosquito survival and mortality as a function of exposure to LLINs or IRS[26] are therefore adapted and simplified as follows to allow for a range of utilization rates ranging from zero to several times per gonotrophic cycle, rather than the previously assumed utilization rate of once per gonotrophic cycle for all blood resources (α v =1). All predictions of impact upon mosquito survival, and the entomologic inoculation rates they mediate, were implemented and parameterized exactly as previously described[51], except that Equation 14 of the original formulation[26] was adapted to incorporate the mortality risks of utilizing all covered and uncovered resources in a more generally applicable manner:

P γ = P γ , 0 P α Rc
(15)

where P γ is the probability of surviving all utilization attempt events for all resources per gonotrophic cycle, Pγ,0 is the probability of surviving all utilization attempt events for all resources per gonotrophic cycle in the absence of any intervention, and P α Rc is the probability of surviving all attempts to utilize intervention-covered forms of the targeted resources per gonotrophic cycle. The probability of surviving all attempts to utilize intervention-covered forms of the targeted resource per gonotrophic cycle P α Rc can be calculated as an exponential decay function of the product of the mortality probability associated with exposure to a covered form of the resource through a single utilization attempt event μ R c , and the mean utilization rate for all covered forms of that resource α R c :

P α Rc = e - μ R c α R c
(16)

By substituting rearranged forms of Equation 10 and then Equation 13 into Equation 16, a solution with two field measurable parameters for the targetable, quantifiable, surveyable, subset is derived:

P α R c = e - μ R c α R C R = e - μ R c α R x , y C R x , y
(17)

It is therefore not necessary to know the proportion of that total resource which the targeted subset represents, or the coverage (C R ) or utilization rate (α R ) for all available forms of a resource. Impact can be predicted directly so long as the coverage of the targeted subset itself C R x , y , and utilization rates for that subset under conditions that enable the intervention to protect it against safe utilization by the mosquito α R x , y , can be measured.

This approach to predicting the survival probability assumes that utilization attempt events are randomly, and independently, distributed across all resource units and mosquitoes. Specifically, the number of times one mosquito utilizes a resource (or resource subset) in one gonotrophic cycle is assumed to be a non-negative integer valued random variable (0, 1, 2, 3…) since the mosquito may not necessarily use the resource or, alternatively, may access it multiple times. Hence, the utilization rate of these resources should be understood as an expected value depending on random events that may be expressed as a mean. This is clearly not the case in relation to obligate utilization of blood from one of all available blood resources (R = v). Each mosquito must utilize one of these available resources to complete the gonotrophic cycle, so complete coverage of all blood resources (C R  = 1) that are utilized at a mean rate of once per gonotrophic cycle (α R  = 1) with an insecticide which induces comprehensive fatality μ R c = 1 would deterministically result in reduction of survival probability to zero P α R , c = 1 - μ R c α R C R = 0 , rather than merely reduced to the minor proportion of mosquitoes that are inaccurately assumed by Equation 16 to have completed a gonotrophic cycle without taking any bloodmeal. However, for a covered subset of a resource (Equation 17), rather than for all available forms of that resource (Equation 16), it is realistic to assume that the number of utilization attempt events per gonotrophic cycle is a random variable for individual mosquitoes and utilization rates per gonotrophic cycle are expected values (expressed as a mean), even for obligate blood resource utilization behaviours.

Measuring utilization rates for subsets of undefined resources by comparison with those for quantifiable blood resources

Adult mosquitoes use many distinct resources during their lifetimes, including several that they need afresh every time they complete a gonotrophic cycle: blood, resting sites, and oviposition sites. Most of these resources are difficult to quantify directly, so the same is true of the rates at which mosquitoes utilize them, thereby making contact with them. However, measurements of feeding rates upon humans or livestock allow ready quantification of absolute mosquito population size or recruitment rate[31]. This is because the size of these mammalian host populations can be conveniently measured by direct census, and blood acquisition occurs at a measurable rate per host[52] with a measurable probability for a given blood host species[22, 23]. Also, blood acquisition usually occurs at a utilization rate of only once per gonotrophic cycle (α v  ≈ 1) where sugar availability is not limiting[45], except for the first gonotrophic cycle where two bloodmeals may be required[53, 54]. For example, the emergence or recruitment rate of mosquitoes (E) in a given setting can be calculated as function of the measured mean biting rate experienced by individual humans (B h ), the number of humans living there (N h ), the proportion of bloodmeals obtained from humans Q v h , and survival probability per feeding cycle (P f )[31]:

E= B h N h Q v h b
(18)

where the mean lifetime number of bloodmeals per emerging mosquito (b) is calculated as the sum of the probabilities of surviving to all plausible gonotrophic ages, expressed as the number of gonotrophic cycles completed (j)[26, 31]:

b= ∑ j = 1 ∞ P f j
(19)

Similarly, for a very zoophagic (predominantly animal-feeding) vector with a strong preference for a known, accessible, manageable non-human host species such as cattle, goats, sheep, pigs or other livestock (l), it may be easier to accurately measure biting rates upon such livestock (B l ) so the equivalent calculation can be made if the proportion of blood obtained from that host species Q v l , and the population size of that host species (N l ) can be determined:

E= B l N l Q v l b
(20)

The key to applying Equations 18 and 20 to estimate absolute vector population sizes is the assumption that the fraction of all available sources of blood that each entomologically surveyed host represents Q v h / N h or Q v l / N l can be readily estimated by host census and bloodmeal identification from a sample of resting, blood-fed mosquitoes, so it is not necessary to directly detect all biting events on all hosts. Generalizing this principle, the emergence rate of mosquitoes (E) can be estimated based on the rate at which mosquitoes are trapped or observed utilizing m R x , z a surveyed sample subset (z) of any targetable subset (x) of a given resource (R x,z ) if the proportion of all available forms of that resource which that surveyed subset represents (R x,z /R), and the rate at which individual mosquitoes utilize all available forms of that resource per gonotrophic cycle (α R ), are both known. Note also that equation 18 and 20 both also implicitly include a term in the denominator for the utilization rate of all blood sources, that was negated by assumed a value approximating unity (α v ≈ 1) but can be explicitly reintroduced for the purposes of generalization. Substituting m R x , z for B h or B l , R x,z /R for Q v h / N h or Q v l / N l , and α R for α v in Equations 18 and 20, respectively, yields the following general formula:

E= m R x , z R R x , z / α R b
(21)

Blood resources can be readily identified as discrete units and their total numbers can be quantified by head-count census. However, units of sugar, resting site, oviposition site, and mating site resources are difficult to define unambiguously, except where these are introduced artificially (sugar baits, houses, boxes, pots, barrier screens, water containers, or swarming markers), and the total quantity of these resources available in the environment is even more difficult, if not impossible, to ascertain. Therefore, it is not obvious how R x,z /R can be estimated for these resources with existing field survey methods. However, this is not necessary to know, because intervention impact can be conveniently rationalized in terms of coverage and utilization rates for definable, targetable subsets of those resources (Equation 17), and it is possible to calculate their relative rates of utilization compared with those for blood from a known proportion of all available blood resources. Here the quantifiable total blood resource, and a surveyed sample (z) of an identifiable subset (x) of that blood resource, is specified (R = v) and distinguished from equivalent terms for other resources, such as resting sites (R = r) or sugar (R = s), with the specific terms v and v x,z . Also, the per gonotrophic rates at which mosquitoes attempt to utilize (α) all blood resources (v) or a distinct, identifiable subset of blood resources (v x ), as well as the per night rate at which utilization attempt events occur (m) at a surveyed sample (z) of that blood resource subset (v x ), are distinguished from those for other resources with the specific terms a v , a v x , and m v x , z , respectively. Given that resource utilization rate per gonotrophic cycle is proportional to the total rate at which utilization attempt events occur in the overall population (Equations 8, 9, 10 and 11), the relative rate of utilization of such a resource subset compared to all blood resources can be expressed by dividing Equation 21, which is specified for a given non-blood resource subset and a surveyed sample thereof (R x /Rx,z), by an equivalent formulation specified for all blood resources, and a surveyed sample of hosts from a subset for which bloodmeals recovered from the midguts of recently fed specimens can be identified and distinguished from other sources (v/vx,z), and rearranging:

α R x α v = m R x m v = m R x , z R x v x , z m v x , z R x , z v
(22)

Note that the emergence and mean longevity terms cancel each other out so that estimates of these parameters are not required to estimate the relative rate of utilization of a resource compared with blood, as described below.

The most obvious vertebrate blood resource subsets (v x ) which are readily surveyed, including detection of blood feeding events and identification of blood source in specimens of fed mosquitoes, are humans (x = h) and livestock (x = l)[52]. The proportion of all available blood resources sampled by the field survey can be quantified as the product of the proportion of bloodmeals obtained from humans Q v h or livestock Q v l and the number of humans (N h,z ) or cattle (N l,z ) sampled by the host attack survey, divided by the total number of humans (N h ) or livestock (N l ) present:

v h , z v = N h , z Q v h N h or v l , z v = N l , z Q v l N l
(23)

Fortunately, while gonotrophic discordance beyond the first feeding cycle does occur in Anopheles, it is unusual and can be quantified[45]. In most cases, it is therefore reasonable to explicitly assume that mosquitoes predominantly utilize blood approximately once per feeding cycle (α v  ≈ 1), so substituting the host-specified (v x  = v h or v l ) formula of Equation 23 into Equation 22, and replacing α v with unity, yields a solution for α R x for which all the terms are measurable in the field:

α R x = m R x , z R x N h , z Q v h m v h , z R x , z N h
(24a)

or

α R x = m R x , z R x N l , z Q v l m v l , z R x , z N l
(24b)

where R x,z /R x is the proportion of all available forms of the targeted non-blood resource subset that was surveyed entomologically to measure the rate per night at which the entire mosquito population attempts to utilize it m R x , z , where N h,z /N h and N l,z /N l are the proportions of all humans or livestock that were respectively surveyed to measure the rate at which mosquitoes attempted to utilize their blood, and where Q v h and Q v l are the proportions of bloodmeals the vector population obtains from all available humans and livestock, respectively.

Where two resources co-occur and overlap completely with each other, specifically the example of resting sites (R = r i ) and human blood indoors within houses (v h,i ), the proportion of each resource subset that is sampled is no longer required because these cancel each other out. All that is required is an estimate of the number of persons or person nights sampled by the host attack survey (N h,z ), and the total number people staying in those sampled houses (N h,Ω ), or even their ratio, which is commonly referred to as the mean number of occupants per house (Nh,Ω/Nh,z):

α R x = m R x , z N h , z Q v h m v h , z N h , Ω
(25a)

In some experiments, however, resting events following more than one bloodmeal are represented in surveys of resting sites because those events may last two or more days. Recent standardized trials to compare various techniques for catching host-seeking and resting mosquitoes[55–57] placed these in or immediately outside of different houses within a defined sampling frame each night, so that the former would not compete with the latter by trapping mosquitoes before they can feed and rest. Also, spray catches must be spaced by intervals of several days to allow residual pyrethrum to dissipate. In both cases, mosquitoes gestating over two or more preceding nights are allowed to accumulate from multiple nights of blood feeding (τ) in the surveyed houses, and this must be accounted for when estimating utilization rates:

α R x = m R x , z N h , z Q v h m v h , z N h , Ω τ
(25b)

Literature review and utilization rate estimate extraction

Studies, or sets of studies, were identified which presented sufficient parameter estimate data for utilization rates of specific, intervention-targetable resource subsets to be calculated for specific malaria vector species in specific, distinct locations. In addition to the authors’ archives of literature and unpublished data, the Pubmed database was also queried with the search term ‘Anopheles AND ((pyrethrum spray OR aspirator) OR (insecticide AND (cattle OR livestock)) OR odour-baited OR sugar)’. For utilization of livestock blood and sugar, consideration was limited to studies in settings where trials of insecticide-treated livestock or sugar baits, respectively, have been either implemented or specifically suggested. To avoid cluttering of the second figure presented in the Results section, consideration of studies enabling estimation of indoor resting site utilization was restricted to recent unpublished studies of our own and those published in the last decade. Where results for a species complex or group were reported, these are attributed to the most common sibling species identified in that population.

Utilization rates for blood from humans while indoors, when they can be protected with LLINs, was calculated as the product of the proportion of human exposure to mosquito bites occurring indoors π v h , i and the human blood index Q v h , by assuming a single bloodmeal per gonotrophic cycle:

α v h , i = π v h , i Q v h
(26)

Where local estimates for the proportion of bloodmeals obtained from humans Q v h were not available for the vector in question, the median values for that species from a previous review[39] were applied. Where direct estimates of the proportion of bloodmeals obtained from cattle and other treated livestock Q v l were not available, utilization of blood from other non-human sources was assumed to be negligible, so that this quantity could be calculated as the complement of the proportion obtained from humans Q v h :

α v l = Q v l ≈1- Q v h
(27)

Utilization rates for odour-baited traps are calculated by assuming that the probability of a mosquito attacking a trap, rather than a natural host, per gonotrophic cycle is equivalent to the proportion of all available human hosts, animal hosts, and pseudo-hosts, that they represent[27]. This is calculated on the basis of the ratio of traps to people (N t /N h ), the relative availability of those traps (λ t ), and the proportion of bloodmeals obtained from non-human hosts 1 - Q v h [27]:

α v t = λ t N t / N h 1 + 1 - Q v h / Q v h + λ t N t / N h
(28)

Utilization rates for sugar resource subsets α s x were calculated as follows, based on direct field measurements of utilization rates per day for dye-labelled sugar baits α s x * [18–20], and an assumed mean gonotrophic cycle duration of three days (g = 3):

α s x = α s x * g
(29)

Utilization rates for resting site subsets α r x , such as the insides of houses or artificial shelters and netting barriers placed in or around them, were estimated using Equation 25a or 25b based on the quotient of the mean rate at which mosquitoes were caught resting inside a sample of them (m r,x,z ) adjusted, where necessary, for an assumed indoor resting period of 2 days (τ = 2), divided by the rate at which they were caught attacking human hosts indoors m v h , i , z , the mean number of exposed occupants per house or room (N h,Ω /N h,z ) adjusted for reported usage rates of LLINs C v h , i , n which were assumed to confer complete protection, the proportion of bloodmeals obtained from humans Q v h , and the proportion of human bloodmeals obtained indoors Ï€ v h , i , with the latter two parameters usually assumed from mean literature values reported for that species[32, 39]:

α r x = m r x , z π v h , i N h , z Q v h m v h , i , z 1 - C v h , i , n N h , Ω τ
(30)

Results

Dependence of impact upon utilization rates of intervention-targeted resource subsets

Figure 1 illustrates how the impact of a vector control intervention with a mosquito-toxic active ingredient depends on the utilization rate of the resource subset to which it is targeted, and on the pre-existing level of transmission mediated by the vector. Utilization rates exceeding, or at least approaching, one event per gonotrophic cycle are required to achieve useful reductions of intense, saturating transmission[4, 58] mediated by human-feeding mosquitoes such as Anopheles arabiensis and Anopheles gambiae. Less efficient vectors that predominantly feed upon animals, such as Anopheles culicifacies, cause far lower levels of transmission and morbidity burden[39] that can respond more sensitively to effective transmission control because baseline exposure levels are not sufficient to saturate the human population with parasite infections[24, 58–61]. Very valuable impact upon transmission by such zoophagic mosquitoes may even be achieved by targeting resource subsets that are used as rarely as once every three feeding cycles because most malaria is transmitted by mosquitoes that are at least four gonotrophic cycles old[31, 54]. However, unless a resource subset is used at least once in every five gonotrophic cycles, it is unlikely to be a useful target for suppressing mosquito survival and population density through vector control, even for weak zoophagic vectors that rarely feed upon humans and mediate modest, manageable levels of transmission.

Figure 1
figure 1

The predicted relationship between the impact of vector control upon malaria transmission and the rate at which mosquitoes utilize the resource subset targeted by the intervention α R x . Anopheles epiroticus, An. arabiensis, and An. gambiae were chosen as examples of vector species that feed primarily upon animals, on both humans and animals, or primarily upon humans. Zero coverage of the human blood resource with LLINs C v h = 0 and high biological coverage of the resource subset which is targeted with an intervention C R x = 0.8 that induces high mortality among mosquitoes that utilize covered forms of that resource μ R c = 0.8 were assumed. Apart from the adaptation of the original formulation to capture intervention impact upon vector survival as a generally applicable function of target resource subset coverage and utilization rate (Equation 15 to 17), all predictions were made by parameterizing and executing deterministic models of malaria transmission and control, exactly as previously described[26, 51, 62].

Field estimates of utilization rates for defined, targetable resource subsets

Figure 2 illustrates the range of utilization rate estimates that could be extracted from the literature for a diversity of resource subsets used by mosquitoes. The predictions of Figure 1, and the generally high utilization rates of human blood while people are indoors (Figure 2), are approximately consistent with the documented impacts of LLINs upon malaria transmission[43]. It is also notable that the lowest estimates for utilization of human blood indoors (Figure 2) are for An. arabiensis, Anopheles farauti, Anopheles darlingi, and Anopheles nuneztovari, all of which have been known to persist and dominate residual transmission systems following scale-up of LLINs or IRS[6–8, 33]. Insecticide-treated livestock also appear to be as promising a target for zoophagic vectors as LLINs are for anthropophagic (predominantly human-feeding) vectors, and the utilization rates estimated for Anopheles culifacies and Anopheles stephensi are consistent (Figure 2) with the proven impact of this approach upon malaria transmission by these species[63]. Odour-baited traps act as pseudo-hosts by mimicking, and even surpassing[64], the taxis stimuli of normal blood sources for mosquitoes, luring them to fatal trap devices or insecticide-treated surfaces[27]. While these can achieve useful utilization rates if their attractiveness and positioning can be optimized, these estimates are somewhat lower than for insecticidal nets and livestock treatments directed at the best-matched vector species (Figure 2) because these devices must compete with natural hosts that, therefore, constitute an inevitable gap in biological coverage[27].

Figure 2
figure 2

Estimated utilization rates for a variety of blood, resting site and sugar resource subsets α R x based on existing field data. See Additional file1 for details of literature, parameter values and equations used to calculate the estimates presented.

Apart from blood, the other important nutrition source that facilitates mosquito survival and malaria transmission is plant-derived sugar[21, 65, 66]. Estimated minimum utilization rates for dye-labelled sugar baits by Anopheles claviger[19], Anopheles sergenti[18], and An. gambiae[20], in three distinct settings appear to be at least comparable with those for utilization of human blood indoors by very anthropophagic vectors, and for utilization of animal blood by zoophagic vectors (Figure 2). The impressive impacts upon all three of these vector populations that have been achieved by adding insecticides to such sugar baits[18–20] are therefore consistent with the predictions outlined in Figure 1, as well as other recent modelling analyses[67]. Given the widespread dependence of mosquitoes upon sugar[21, 65, 66], many important vector populations probably use it at similarly high rates, especially when infected with malaria parasites[68]. Mosquitoes should therefore be at least as amenable to control with this approach as anthropophagic vectors are to LLINs[43], and as zoophagic vectors are to insecticide-treated livestock[63].

Utilization rate estimates for indoor resting sites (Figure 2) are generally lower than that required to explain (Figure 1) the often massive impact of IRS[13] and insecticidal wall linings[69] on many target vector species. However, this is not entirely surprising because even the best techniques for capturing mosquitoes resting indoors, such as pyrethrum spray catch[52] and backpack aspirators[70], are known to consistently under-sample them. Nevertheless, the estimates of resting site utilization presented in Figure 2 are clearly useful for comparison of different potential resting site targets, confirming that comprehensive spraying or lining of entire rooms and houses is probably superior to targeted treatment of pots, boxes (some of which were baited with host odours) or screening barriers, all of which were placed indoors for endophilic vectors or outdoors for exophilic ones[13–17]. However, perhaps the most important observation in relation to these estimates of resting site subset utilization rates is that they rely on upon quite complicated calculations requiring at least five distinct input parameters (Equations 25a, 25b and 30), many of which have to be assumed based on best guesses or literature values from a different setting (Additional file1). In fact, none of the estimates presented in Figure 2 are based entirely upon local estimates for all the input parameters (Additional file1), and are therefore not entirely independent of each other or representative of the full range of values for any of the vector species described.

Discussion

Defining and surveying targetable resource subsets

Estimating coverage and utilization of a resource subset primarily depends upon defining it in a quantifiable manner that can be readily surveyed and targeted, or artificially created in the field. The most obvious and familiar examples are the human populations targeted for universal coverage with LLINs to protect the blood resource they represent to mosquitoes[28, 43, 71]. While wild animals are difficult to survey or deliver interventions to, livestock represent blood resources that can be readily quantified and targeted with interventions[63]. It is even easier to track numbers and functionality of artificial odour-baited traps, which mimic and compete with natural blood sources for the attentions of host-seeking mosquitoes, so their potential impact can also be predicted as a function of biological coverage of all available hosts and pseudo-hosts[27].

The subset of all resting sites represented by the inner surfaces of human dwellings (walls, ceilings and even furniture) are the defined target for IRS[12, 71], as well as insecticide-treated tents[72], shelters[73, 74], and wall linings[69], so coverage can be quantified as the proportion of residential structures treated. On the other hand, artificially introduced pots, boxes, curtains, linings or screening barriers compete with natural resting sites[13–17]. Mosquitoes can be captured relatively efficiently on these well-defined, convenient, standardized surfaces, so it has been suggested that these could also be treated with toxic insecticides to improve control efficiency[13–17]. While it remains difficult to consistently identify and define sugar, oviposition site or mating site resources[5], recent progress with observational[75, 76], trapping[77], tracing and labelling[18–20, 22, 23, 78–81] methods for mosquitoes is encouraging.

Adapting entomological survey techniques to measure resource utilization rates

Comparing the range of utilization rates described in Figure 2 with the predictions of potential impact illustrated in Figure 1 confirms that, despite their known limitations[52, 82], existing entomological field methods may be very useful for designing and evaluating a wide diversity of vector control products[83]. Both blood and sugar meals can be readily identified using a variety of naturally-occurring markers and artificially added tracers[18–20, 22, 23, 78–81], thus enabling very direct, robust measurement of label uptake as a function of time or age. Utilization rates can therefore be estimated directly for subsets of these naturally occurring resources (Equations 26, 27, and 29) or indirectly for artificially introduced subsets such as odour-baited traps (Equation 28).

Utilization rates for resting site subsets (Equation 30), or indeed any other non-blood resource (Equation 25a and 25b), can also be estimated indirectly by calibrating against measurable utilization rates for quantifiable, preferred blood sources. However, the complexity of these models, and their reliance upon local measurements of several entomological input parameters, all of which have limited precision and accuracy, may well limit broader application of this approach beyond the crude application to existing data presented in Figure 2. Recent attempts to rejuvenate and improve existing entomological survey methodology for detecting resource utilization attempt events with electrified grids[84–86], sticky traps[77, 87], mechanized aspirators[70], and high resolution cameras[75, 76], should enable improved sensitivity of utilization attempt event detection at surveyed samples of resource subset targets. However, while such technical advances may well address the inaccuracies of attempts to estimate utilization rates for subsets of resting sites or other non-nutritional resources by improving event detection sensitivity, they are unlikely to improve their precision because considerable uncertainty arises from the need for relatively complex models (Equations 25a, 25b and 30) that require correspondingly numerous measurements of input parameters.

Fortunately, a wide range of more sensitive chemical, biochemical, genetic and biological markers, that could be applied to labelling mosquitoes when they use these other resources, are now available[82] but these remain to be fully exploited. In fact, field studies using artificial tracers to label of both mosquitoes feeding upon sugar[18–20] and sand flies feeding upon rodent blood[88, 89], in which addition of insecticide removed almost all labelled insects from these vector populations, clearly demonstrate the validity of this strategy as a means to estimate biological coverage or utilization rates. The major advantage of labelling mosquitoes when they utilize a resource subset, rather than trapping or observing them, is that the measured proportions of marked insects can be readily analyzed with robust off-the-shelf statistical methods for binary outcomes, and are relatively precise because they have a nominator and denominator which both vary in proportion to population size or event detection sensitivity. The largest caveat to this approach is that essentially all targeted forms of that resource must be labelled on geographic scales large enough to negate the effect that mosquito dispersal has upon measurements of label uptake: immigration of unlabelled mosquitoes into the study area will increase the denominator while emigration will reduce the nominator, so that true local utilization rates will be systematically underestimated[52, 90, 91]. However, this phenomenon could also be exploited to great advantage if multiple distinct labels for various treatment arms were used to measure, and correct for, the effects of mosquito upon impact distribution in large-scale trials of vector control interventions[92].

The conceptual framework and entomological measurement priorities outlined here should be readily and directly applicable to almost any population of mosquitoes, vectors or other pest. It should therefore be possible to simultaneously tackle multiple vectors with integrated vector management programmes[93] that prioritize interventions based on simultaneous, comparative field assessment of respective utilization rates for each potential target species. Recent demonstrations of the usefullness of dummy bait products containing appropriate labels but no insecticide[18–20, 88, 89] illustrate how cost-effective, robust measurements of utilization rates could be used to select and optimize available technologies for immediate use or new prototypes for development.

Conclusions

The concept of biological coverage can be extended to enable prediction of intervention impact for diverse vector control strategies based on estimated utilization rates for any definable, targetable resource subset. Indeed the applicability of this approach has been demonstrated here using existing entomological measurement methods to rationalize the observed impacts of LLINs, insecticide-treated livestock, and attractive toxic sugar baits upon malaria vectors. The development of improved and diversified technologies for controlling transmission of malaria, as well as a diversity of other vector-borne pathogens, could therefore be accelerated, rationalized and streamlined based on field measurements of the rates at which mosquitoes utilize targetable biological resource subsets.

While blood and sugar consumption can be readily quantified using existing methods for detecting natural markers or artificial tracers, improved techniques for labelling mosquitoes will be required to assess and optimize vector control measures which target them when they utilize resting, oviposition and mating sites. All mosquito species need sugar, resting sites, oviposition sites, and mating sites, as indeed do most arthropods of medical and veterinary importance. These resources are therefore important potential targets for the new or improved vector control methods that are clearly needed to eliminate malaria, and also a variety of other vector-borne pathogens. To enable comparative assessment of all potential resource subset targets, including sites which mosquitoes rest, oviposit or mate at, existing tracer technologies need be adapted to enable reliable, non-toxic, non-disruptive labelling of mosquitoes when they utilize these non-nutritional resource subset targets.

References

  1. Wenger EA, Eckhoff PA: A mathematical model of the impact of present and future malaria vaccines. Malar J. 2013, 12: 126-10.1186/1475-2875-12-126.

    Article  PubMed Central  PubMed  Google Scholar 

  2. Eckhoff PA: Mathematical models of within-host and transmission dynamics to determine effects of malaria Interventions in a variety of transmission settings. Am J Trop Med Hyg. 2013, 88: 817-827. 10.4269/ajtmh.12-0007.

    Article  PubMed Central  PubMed  Google Scholar 

  3. Okell LC, Griffin JT, Kleinschmidt I, Hollingsworth TD, Churcher TS, White MT, Bousema T, Drakeley CJ, Ghani AC: The potential contribution of mass treatment to the control of Plasmodium falciparum malaria. PLoS One. 2011, 6: e20179-10.1371/journal.pone.0020179.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. Killeen GF: A second chance to tackle African malaria vector mosquitoes that avoid houses and don’t take drugs. Am J Trop Med Hyg. 2013, 88: 809-816. 10.4269/ajtmh.13-0065.

    Article  PubMed Central  PubMed  Google Scholar 

  5. Ferguson HM, Dornhaus A, Beeche A, Borgemeister C, Gottlieb M, Mulla MS, Gimnig JE, Fish D, Killeen GF: Ecology: a prerequisite for malaria elimination and eradication. PLoS Med. 2010, 7: e1000303-10.1371/journal.pmed.1000303.

    Article  PubMed Central  PubMed  Google Scholar 

  6. Russell TL, Beebe NW, Cooper RD, Lobo NF, Burkot TR: Successful malaria elimination strategies require interventions that target changing vector behaviours. Malar J. 2013, 12: 56-10.1186/1475-2875-12-56.

    Article  PubMed Central  PubMed  Google Scholar 

  7. Durnez L, Coosemans M: Residual Transmission of Malaria: An Old Issue for New Approaches. Anopheles Mosquitoes – New Insights into Malaria Vectors. Edited by: Manguin S. 2013, Rijeka: Intech, 671-704.

    Google Scholar 

  8. Govella NJ, Chaki PP, Killeen GF: Entomological surveillance of behavioural resilience and resistance in residual malaria vector populations. Malar J. 2013, 12: 124-10.1186/1475-2875-12-124.

    Article  PubMed Central  PubMed  Google Scholar 

  9. Killeen GF, Seyoum A, Sikaala CH, Zomboko AS, Gimnig JE, Govella NJ, White MT: Eliminating malaria vectors. Parasit Vectors. 2013, 6: 172-10.1186/1756-3305-6-172.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Sinka ME, Bangs MJ, Manguin S, Rubio-Palis Y, Chareonviriyaphap T, Coetzee M, Mbogo CM, Hemingway J, Patil AP, Temperley WH, Gething PW, Kabaria CW, Burkot TR, Harbach RE, Hay SI: A global map of dominant malaria vectors. Parasit Vectors. 2012, 5: 69-10.1186/1756-3305-5-69.

    Article  PubMed Central  PubMed  Google Scholar 

  11. Kiware SS, Chitnis N, Devine GJ, Moore SJ, Majambere S, Killeen GF: Biologically meaningfull coverage indicators for eliminating malaria transmission. Biol Lett. 2012, 8: 874-877. 10.1098/rsbl.2012.0352.

    Article  PubMed Central  PubMed  Google Scholar 

  12. Pluess B, Tanser FC, Lengeler C, Sharp BL: Indoor residual spraying for preventing malaria. Cochrane Database Syst Rev. 2010, 4: CD006657

    Google Scholar 

  13. Scholte EJ, Ng’habi K, Kihonda J, Takken W, Paaijmans K, Abdulla S, Killeen GF, Knols BG: An entomopathogenic fungus for control of adult African malaria mosquitoes. Science. 2005, 308: 1641-1642. 10.1126/science.1108639.

    Article  CAS  PubMed  Google Scholar 

  14. Odiere M, Bayoh MN, Gimnig J, Vulule J, Irungu L, Walker E: Sampling outdoor, resting Anopheles gambiae and other mosquitoes (Diptera: Culicidae) in Western Kenya with clay pots. J Med Entomol. 2007, 44: 14-22. 10.1603/0022-2585(2007)44[14:SORAGA]2.0.CO;2.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Farenhorst M, Farina D, Scholte EJ, Takken W, Hunt RH, Coetzee M, Knols BGJ: African water storage pots for the delivery of the entomopathogenic fungus Metarhizium anisopliae to the malaria vectors Anopheles gambiae s.s. and Anopheles funestus. Am J Trop Med Hyg. 2008, 78: 910-916.

    PubMed  Google Scholar 

  16. van den Bijllaardt W, Ter Braak R, Shekalaghe S, Otieno S, Mahande A, Sauerwein R, Takken W, Bousema T: The suitability of clay pots for indoor sampling of mosquitoes in an arid area in northern Tanzania. Acta Trop. 2009, 111: 197-199. 10.1016/j.actatropica.2009.04.003.

    Article  PubMed  Google Scholar 

  17. Müller G, Schlein Y: Sugar-questing mosquitoes in arid areas gather on scarce blossoms that can be used for control. Int J Parasitol. 2006, 36: 1077-1080. 10.1016/j.ijpara.2006.06.008.

    Article  PubMed  Google Scholar 

  18. Burkot TR, Russell TL, Reimer LJ, Bugoro H, Beebe NW, Cooper RD, Sukawati S, Collins FH, Lobo NF:Barrier screens: a method to sample blood-fed and host-seeking exophilic mosquitoes. Malar J. 2013, 12: 49-10.1186/1475-2875-12-49.

    Article  PubMed Central  PubMed  Google Scholar 

  19. Müller G, Schlein Y: Efficacy of toxic sugar baits against cistern-dwelling Anopheles claviger. Trans R Soc Trop Med Hyg. 2008, 102: 480-484. 10.1016/j.trstmh.2008.01.008.

    Article  PubMed  Google Scholar 

  20. Müller GC, Beier JC, Traore SF, Toure MB, Traore MM, Bah S, Doumbia S, Schlein Y:Successful field trial of attractive toxic sugar bait (ATSB) plant-spraying methods against malaria vectors in theAnopheles gambiaecomplex in Mali, West Africa.Malar J. 2010, 9: 210-10.1186/1475-2875-9-210.

    Article  PubMed Central  PubMed  Google Scholar 

  21. Foster WA: Mosquito sugar feeding and reproductive energetics. Annu Rev Entomol. 1995, 40: 443-474. 10.1146/annurev.en.40.010195.002303.

    Article  CAS  PubMed  Google Scholar 

  22. Garrett-Jones C: The human blood index of malarial vectors in relationship to epidemiological assessment. Bull World Health Organ. 1964, 30: 241-261.

    PubMed Central  CAS  PubMed  Google Scholar 

  23. Garrett-Jones C, Boreham P, Pant CP: Feeding habits of anophelines (Diptera: Culicidae) in 1971–1978, with reference to the human blood index: a review. Bull Entomol Res. 1980, 70: 165-185. 10.1017/S0007485300007422.

    Article  Google Scholar 

  24. Kiware SS, Chitnis N, Moore SJ, Devine GJ, Majambere S, Killeen GF: Simplified models of vector control impact upon malaria transmission by zoophagic mosquitoes. PLoS One. 2012, 7: e37661-10.1371/journal.pone.0037661.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

  26. Killeen GF, Chitnis N, Moore SJ, Okumu FO: Target product profile choices for intra-domiciliary malaria vector control pesticide products: repel or kill?. Malar J. 2011, 10: 207-10.1186/1475-2875-10-207.

    Article  PubMed Central  PubMed  Google Scholar 

  27. Okumu FO, Moore SJ, Govella NJ, Chitnis N, Killeen GF: Potential benefits, limitations and target product-profiles of odor-baited mosquito traps as a means of malaria control. PLoS One. 2010, 5: e11573-10.1371/journal.pone.0011573.

    Article  PubMed Central  PubMed  Google Scholar 

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

  29. Killeen GF, Smith TA: Exploring the contributions of bednets, cattle, insecticides and excito-repellency to malaria control: A deterministic model of mosquito host-seeking behaviour and mortality. Trans R Soc Trop Med Hyg. 2007, 101: 867-880. 10.1016/j.trstmh.2007.04.022.

    Article  PubMed Central  PubMed  Google Scholar 

  30. Killeen GF, Moore SJ: Target product profiles for protecting against outdoor malaria transmission. Malar J. 2012, 11: 17-10.1186/1475-2875-11-17.

    Article  PubMed Central  PubMed  Google Scholar 

  31. Killeen GF, McKenzie FE, Foy BD, Schieffelin C, Billingsley PF, Beier JC: A simplified model for predicting malaria entomologic inoculation rates based on entomologic and parasitologic parameters relevant to control. Am J Trop Med Hyg. 2000, 62: 535-544.

    PubMed Central  CAS  PubMed  Google Scholar 

  32. Huho BJ, Briët O, Seyoum A, Sikaala CH, Bayoh N, Gimnig JE, Okumu FO, Diallo D, Abdulla S, Smith TA, Killeen GF: Consistently high estimates for the proportion of human exposure to malaria vector populations occurring indoors in rural Africa. Int J Epidemiol. 2013, 42: 235-247. 10.1093/ije/dys214.

    Article  PubMed Central  PubMed  Google Scholar 

  33. Elliott R: The influence of vector behaviour upon malaria transmission. Am J Trop Med Hyg. 1972, 21: 755-763.

    CAS  PubMed  Google Scholar 

  34. Seyoum A, Sikaala CH, Chanda J, Chinula D, Ntamatungiro AJ, Hawela M, Miller JM, Russell TL, Briët OJT, Killeen GF: Most exposure to Anopheles funestus and Anopheles quadriannulatus in Luangwa valley, south-east Zambia occurs indoors, even for users of insecticidal nets. Parasit Vectors. 2012, 5: 101-10.1186/1756-3305-5-101.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Elliott R: Studies on Man-Vector Contact in Some Malarious Areas in Colombia. 1967, Geneva: World Health Organization, 32-

    Google Scholar 

  36. Garrett-Jones C: A Method for Estimating the Man-Biting Rate. 1964, Geneva: World Health Organization, 22-

    Google Scholar 

  37. Bugoro H, Cooper RD, Butafa C, Iro’ofa C, Mackenzie DO, Chen CC, Russell TL: Bionomics of the malaria vector Anopheles farauti in Temotu Province, Solomon Islands: issues for malaria elimination. Malar J. 2011, 10: 133-10.1186/1475-2875-10-133.

    Article  PubMed Central  PubMed  Google Scholar 

  38. Russell TL, Govella NJ, Azizi S, Drakeley CJ, Kachur SP, Killeen GF: Increased proportions of outdoor feeding among residual malaria vector populations following increased use of insecticide-treated nets in rural Tanzania. Malar J. 2011, 10: 80-10.1186/1475-2875-10-80.

    Article  PubMed Central  PubMed  Google Scholar 

  39. Kiswewski AE, 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.

    Google Scholar 

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

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

    Article  Google Scholar 

  42. Killeen GF, Seyoum A, Knols BGJ: Rationalizing historical successes of malaria control in Africa in terms of mosquito resource availability management. Am J Trop Med Hyg. 2004, 71 (Supplement 2): 87-93.

    PubMed  Google Scholar 

  43. Lengeler C: Insecticide-treated bed nets and curtains for preventing malaria. Cochrane Database Syst Rev. 2004, 2: CD000363

    Google Scholar 

  44. Clements AN: Development, Nutrition and Reproduction. The Biology of Mosquitoes. Volume 1. 1992, London: Chapman & Hall

    Google Scholar 

  45. Beier JC: Frequent blood-feeding and restrictive sugar-feeding behavior enhance the malaria vector potential of Anopheles gambiae s.l. and An. funestus (Diptera:Culicidae) in western Kenya. J Med Entomol. 1996, 33: 613-618.

    Article  CAS  PubMed  Google Scholar 

  46. Achee NL, Bangs MJ, Farlow R, Killeen GF, Lindsay S, Logan JG, Moore SJ, Rowland M, Sweeney K, Torr SJ, Zwiebel LJ, Grieco JP: Spatial repellents: from discovery and development to evidence-based validation. Malar J. 2012, 11: 164-10.1186/1475-2875-11-164.

    Article  PubMed Central  PubMed  Google Scholar 

  47. Kambris Z, Blagborough AM, Pinto SB, Blagrove MSC, Godfray HCJ, Sinden RE, Sinkins SP: Wolbachia stimulates immune gene expression and inhibits Plasmodium development in Anopheles gambiae. PLoS Pathog. 2010, 6: e1001143-10.1371/journal.ppat.1001143.

    Article  PubMed Central  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  49. Devine GJ, Perea EZ, Killeen GF, Stancil JD, Clark SJ, Morrison AC: Autodissemination of an insecticide by adult mosquitoes drammatically amplifies lethal coverage of their aquatic habitats. Proc Natl Acad Sci U S A. 2009, 106: 11530-11534. 10.1073/pnas.0901369106.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  50. Harris AF, McKemey AR, Nimmo D, Curtis Z, Black I, Morgan SA, Oviedo MN, Lacroix R, Naish N, Morrison NI, Collado A, Stevenson J, Scaife S, Dafa’alla T, Fu G, Phillips C, Miles A, Raduan N, Kelly N, Beech C, Donnelly CA, Petrie WD, Alphey L: Successful suppression of a field mosquito population by sustained release of engineered male mosquitoes. Nat Biotechnol. 2012, 30: 828-830. 10.1038/nbt.2350.

    Article  CAS  PubMed  Google Scholar 

  51. Killeen GF, Okumu FO, N’Guessan R, Coosemans M, Adeogun A, Awolola S, Etang J, Dabiré RK, Corbel V: The importance of considering community-level effects when selecting insecticidal malaria vector products. Parasit Vectors. 2011, 4: 160-10.1186/1756-3305-4-160.

    Article  PubMed Central  PubMed  Google Scholar 

  52. Silver JB, Service MW: Mosquito Ecology: Field Sampling Methods. 2008, Dordrecht, the Netherlands: Springer, 3

    Book  Google Scholar 

  53. Gillies MT: A modified technique for age grading populations of Anopheles gambiae. Ann Trop Med Parasitol. 1958, 58: 261-273.

    Google Scholar 

  54. Gillies MT, Wilkes TJ: A study of the age-composition of populations of Anopheles gambiae Giles and A. funestus Giles in North-Eastern Tanzania. Bull Entomol Res. 1965, 56: 237-262. 10.1017/S0007485300056339.

    Article  CAS  PubMed  Google Scholar 

  55. Govella NJ, Chaki P, Mpangile JM, Killeen GF: Monitoring mosquitoes in urban Dar es Salaam: Evaluation of resting boxes, window exit traps, CDC light traps,Ifakara tent traps and human landing catches. Parasit Vectors. 2011, 4: 40-10.1186/1756-3305-4-40.

    Article  PubMed Central  PubMed  Google Scholar 

  56. Sikaala CH, Killeen GF, Chanda J, Chinula D, Miller J, Russell TL, Seyoum A: Evaluation of alternative mosquito sampling methods for malaria vectors in lowland south-east Zambia. Parasit Vectors. 2013, 6: 91-10.1186/1756-3305-6-91.

    Article  PubMed Central  PubMed  Google Scholar 

  57. Wong J, Bayoh MN, Olang G, Killeen GF, Hamel MJ, Vulule JM, Gimnig JE: Standardizing operational vector sampling techniques for measuring malaria transmission intensity: Evaluation of six mosquito collection methods in western Kenya. Malar J. 2013, 12: 143-10.1186/1475-2875-12-143.

    Article  PubMed Central  PubMed  Google Scholar 

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

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

  60. Smith DL, Hay SI: Endemicity response timelines for Plasmodium falciparum elimination. Malar J. 2009, 8: 87-10.1186/1475-2875-8-87.

    Article  PubMed Central  PubMed  Google Scholar 

  61. Smith TA, Maire N, Dietz K, Killeen GF, Vounatsou P, Molineaux L, Tanner M: Relationship between entomologic inoculation rate and the force of infection for Plasmodium falciparum malaria. Am J Trop Med Hyg. 2006, 75 (Supplement 2): 11-18.

    PubMed  Google Scholar 

  62. Killeen GF, Seyoum A, Gimnig JE, Stevenson JC, Drakeley CJ, Chitnis N: Made-to-measure malaria vector control strategies: rational design based on insecticide properties and coverage of blood resources for mosquitoes. Malar J. 2014, 13: 146-10.1186/1475-2875-13-146.

    Article  PubMed Central  PubMed  Google Scholar 

  63. Rowland M, Durrani N, Kenward M, Mohammed N, Urahman H, Hewitt S: Control of malaria in Pakistan by applying deltamethrin insecticide to cattle: a community-randomised trial. Lancet. 2001, 357: 1837-1841. 10.1016/S0140-6736(00)04955-2.

    Article  CAS  PubMed  Google Scholar 

  64. Okumu FO, Killeen GF, Ogoma SB, Biswaro L, Smallegange RC, Mbeyela E, Titus E, Munk C, Ngonyani H, Takken W, Mshinda H, Mukabana WR, Moore SJ: Development and field evaluation of a mosquito lure that is more attractive than humans. PLoS One. 2010, 5: e8591-10.1371/journal.pone.0008591.

    Article  Google Scholar 

  65. Gu W, Muller G, Schlein Y, Novak RJ, Beier JC: Natural plant sugar sources of Anopheles mosquitoes strongly impact malaria transmission potential. PLoS One. 2011, 6: e15996-10.1371/journal.pone.0015996.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  66. Huestis DL, Yaro AS, Traore AI, Adamou A, Kassogue Y, Diallo M, Timbine S, Dao A, Lehmann T: Variation in metabolic rate of Anopheles gambiae and A. arabiensis in a Sahelian village. J Exp Biol. 2011, 214: 2345-2353. 10.1242/jeb.054668.

    Article  PubMed Central  PubMed  Google Scholar 

  67. Marshall JM, White MT, Ghani AC, Schlein Y, Muller GC, Beier JC: Quantifying the mosquito’s sweet tooth: modelling the effectiveness of attractive toxic sugar baits (ATSB) for malaria vector control. Malar J. 2013, 12: 291-10.1186/1475-2875-12-291.

    Article  PubMed Central  PubMed  Google Scholar 

  68. Nyasembe VO, Teal PE, Sawa P, Tumlinson JH, Borgemeister C, Torto B: Plasmodium falciparum infection increases Anopheles gambiae attraction to nectar sources and sugar uptake. Curr Biol. 2014, 24: 217-221. 10.1016/j.cub.2013.12.022.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  69. Messenger LA, Matias A, Manana AN, Stiles-Ocran JB, Knowles S, Boakye DA, Coulibaly MB, Larsen M-L, Traoré AS, Diallo B, Konaté M, Guindo A, Traoré SF, Mulder CEG, Le H, Kleinschmidt I, Rowland M: Multicentre studies of insecticide-treated durable wall lining in Africa and South-East Asia: entomological efficacy and household acceptability during one year of field use. Malar J. 2012, 11: 358-10.1186/1475-2875-11-358.

    Article  PubMed Central  PubMed  Google Scholar 

  70. Maia MF, Robinson A, John AN, Mgando J, Simfukwe E, Moore SJ: Comparison of the CDC Backpack aspirator and the Prokopack aspirator for sampling indoor- and outdoor-resting mosquitoes in southern Tanzania. Parasit Vectors. 2011, 4: 124-10.1186/1756-3305-4-124.

    Article  PubMed Central  PubMed  Google Scholar 

  71. WHO: World Malaria Report 2013. 2013, Geneva: World Health Organization

    Google Scholar 

  72. Graham K, Rehman H, Ahmad M, Kamal M, Khan I, Rowland M: Tents pre-treated with insecticide for malaria control in refugee camps: an entomological evaluation. Malar J. 2004, 3: 25-10.1186/1475-2875-3-25.

    Article  PubMed Central  PubMed  Google Scholar 

  73. Burns M, Rowland M, N’guessan R, Carneiro I, Beeche A, Ruiz SS, Kamara S, Takken W, Carnevale P, Allan R: Insecticide-treated plastic sheeting for emergency malaria prevention and shelter among displaced populations: an observational cohort study in a refugee setting in Sierra Leone. Am J Trop Med Hyg. 2012, 87: 242-250. 10.4269/ajtmh.2012.11-0744.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  74. Mittal PK, Sreehari U, Razdan RK, Dash AP: Evaluation of the impact of ZeroFly®, an insecticide incorporated plastic sheeting on malaria incidence in two temporary labour shelters in India. J Vector Borne Dis. 2012, 48: 138-143.

    Google Scholar 

  75. Diabate A, Yaro AS, Dao A, Diallo M, Huestis DL, Lehmann T: Spatial distribution and male mating success of Anopheles gambiae swarms. BMC Evol Biol. 2011, 11: 184-10.1186/1471-2148-11-184.

    Article  PubMed Central  PubMed  Google Scholar 

  76. Butail S, Manoukis N, Diallo M, Ribeiro JM, Lehmann T, Paley DA: Reconstructing the flight kinematics of swarming and mating in wild mosquitoes. J R Soc Interface. 2012, 9: 2624-2638. 10.1098/rsif.2012.0150.

    Article  PubMed Central  PubMed  Google Scholar 

  77. Harris C, Kihonda J, Lwetoijera D, Dongus S, Devine G, Majambere S: A simple and efficient tool for trapping gravid Anopheles at breeding sites. Parasit Vectors. 2011, 4: 125-10.1186/1756-3305-4-125.

    Article  PubMed Central  PubMed  Google Scholar 

  78. Alcaide M, Rico C, Ruiz S, Soriguer R, Munoz J, Figuerola J: Disentangling vector-borne transmission networks: a universal DNA barcoding method to identify vertebrate hosts from arthropod bloodmeals. PLoS One. 2009, 4: e7092-10.1371/journal.pone.0007092.

    Article  PubMed Central  PubMed  Google Scholar 

  79. Fornadel CM, Norris DE: Increased endophily by the malaria vector Anopheles arabiensis in southern Zambia and identification of digested blood meals. Am J Trop Med Hyg. 2008, 79: 876-880.

    PubMed Central  CAS  PubMed  Google Scholar 

  80. Junnila A, Muller GC, Schlein Y: Species identification of plant tissues from the gut of Anopheles sergentii by DNA analysis. Acta Trop. 2010, 115: 227-233. 10.1016/j.actatropica.2010.04.002.

    Article  CAS  PubMed  Google Scholar 

  81. Manda H, Gouagna LC, Nyandat E, Kabiru EW, Jackson RR, Foster WA, Githure JI, Beier JC, Hassanali A: Discriminative feeding behaviour of Anopheles gambiae s.s. on endemic plants in western Kenya. Med Vet Entomol. 2007, 21: 103-111. 10.1111/j.1365-2915.2007.00672.x.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  82. Hagler JR, Jackson CG: Methods for marking insects: current techniques and future prospects. Annu Rev Entomol. 2001, 46: 511-543. 10.1146/annurev.ento.46.1.511.

    Article  CAS  PubMed  Google Scholar 

  83. Vontas J, Moore SJ, Kleinschmidt I, Ranson H, Lindsay SW, Lengeler C, Hamon N, McLean T, Hemingway J: Framework for rapid assessment and adoption of new vector control tools. Trend Parasitol. 2014, 30: 191-204. 10.1016/j.pt.2014.02.005.

    Article  Google Scholar 

  84. Dugassa S, Lindh JM, Torr SJ, Oyieke F, Lindsay SW, Fillinger U: Electric nets and sticky materials for analysing oviposition behaviour of gravid malaria vectors. Malar J. 2012, 11: 374-10.1186/1475-2875-11-374.

    Article  PubMed Central  PubMed  Google Scholar 

  85. Majambere S, Massue DJ, Mlacha Y, Govella NJ, Magesa SM, Killeen GF: Advantages and limitations of commercially available electrocuting grids for studying mosquito behaviour. Parasit Vectors. 2013, 6: 53-10.1186/1756-3305-6-53.

    Article  PubMed Central  PubMed  Google Scholar 

  86. Torr SJ, Della Torre A, Calzetta M, Costantini C, Vale GA:Towards a fuller understanding of mosquito behaviour: use of electrocuting grids to compare the odour-orientated responses ofAnopheles arabiensisandAn. quadriannulatusin the field.Med Vet Entomol. 2008, 22: 93-108. 10.1111/j.1365-2915.2008.00723.x.

    Article  CAS  PubMed  Google Scholar 

  87. Müller GC, Beier JC, Traore SF, Toure MB, Traore MM, Bah S, Doumbia S, Schlein Y: Field experiments of Anopheles gambiae attraction to local fruits/seedpods and flowering plants in Mali to optimize strategies for malaria vector control in Africa using attractive toxic sugar bait methods. Malar J. 2010, 9: 262-10.1186/1475-2875-9-262.

    Article  PubMed Central  PubMed  Google Scholar 

  88. Mascari TM, Stout RW, Clark JW, Gordon SW, Bast JD, Foil LD: Insecticide-treated rodent baits for sand fly control. Pestic Biochem Physiol. 2013, 106: 113-117. 10.1016/j.pestbp.2013.04.006.

    Article  CAS  Google Scholar 

  89. Mascari TM, Hanafi HA, Jackson RE, Ouahabi S, Ameur B, Faraj C, Obenauer PJ, Diclaro JW, Foil LD: Ecological and control techniques for sand flies (Diptera: Psychodidae) associated with rodent reservoirs of leishmaniasis. PLoS Negl Trop Dis. 2013, 7: e2434-10.1371/journal.pntd.0002434.

    Article  PubMed Central  PubMed  Google Scholar 

  90. Killeen GF, Knols BG, Gu W: Taking malaria transmission out of the bottle: implications of mosquito dispersal for vector-control interventions. Lancet Infect Dis. 2003, 3: 297-303. 10.1016/S1473-3099(03)00611-X.

    Article  PubMed  Google Scholar 

  91. Service MW: Mosquito (Diptera: Culicidae) dispersal-the long and short of it. J Med Entomol. 1997, 34: 579-588.

    Article  CAS  PubMed  Google Scholar 

  92. Hawley WA, Phillips-Howard PA, Ter Kuile FO, Terlouw DJ, Vulule JM, Ombok M, Nahlen BL, Gimnig JE, Kariuki SK, Kolczak MS, Hightower AW: Community-wide effects of permethrin-treated bednets on child mortality and malaria morbidity in western Kenya. Am J Trop Med Hyg. 2003, 68 (Supplement 4): 121-127.

    PubMed  Google Scholar 

  93. WHO: Global Strategic Framework for Integrated Vector Management. 2004, Geneva: World Health Organization

    Google Scholar 

Download references

Acknowledgements

We thank Dr. K Aultman and Dr. D Malone for discussions that stimulated and influenced the content of this manuscript. We thank Prof T A Smith for guidance on the probablistic basis of exponential decay models, and Dr. T R Burkot for critical comments on the manuscript, as well as providing population size data for Haleta. We are also grateful to three anonymous reviewers, whose comments had a substantive influence on the final interpretation and conclusions. This work was funded by the Bill & Melinda Gates Foundation (Award numbers 45114, 52644 and OPP1032350).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Gerry F Killeen.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

GFK, AS, JEG, JS, and CJD developed the study hypotheses, contributed data for the analysis, and assisted in the revising the manuscript. GFK conceived the study, developed and applied the model and drafted the manuscript in consultation with the other authors. NC, SSK, and GC refined the mathematical formulation of the model. NC oversaw finalization of the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12936_2014_3388_MOESM1_ESM.xls

Additional file 1: Data sources, publications and calculations supporting Figure 2.(XLS 54 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 credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Killeen, G.F., Kiware, S.S., Seyoum, A. et al. Comparative assessment of diverse strategies for malaria vector population control based on measured rates at which mosquitoes utilize targeted resource subsets. Malar J 13, 338 (2014). https://doi.org/10.1186/1475-2875-13-338

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1475-2875-13-338

Keywords