Skip to main content

Push by a net, pull by a cow: can zooprophylaxis enhance the impact of insecticide treated bed nets on malaria control?



Mass insecticide treated bed net (ITN) deployment, and its associated coverage of populations at risk, had “pushed” a decline in malaria transmission. However, it is unknown whether malaria control is being enhanced by zooprophylaxis, i.e., mosquitoes diverted to feed on hosts different from humans, a phenomenon that could further reduce malaria entomological transmission risk in areas where livestock herding is common.


Between May and July 2009, we collected mosquitoes in 104 houses from three neighboring villages with high ITN coverage (over 80%), along Lake Victoria. We also performed a census of livestock in the area and georeferenced tethering points for all herds, as well as, mosquito larval habitats. Bloodmeal contents from sampled mosquitoes were analyzed, and each mosquito was individually tested for malaria sporozoite infections. We then evaluated the association of human density, ITN use, livestock abundance and larval habitats with mosquito abundance, bloodfeeding on humans and malaria sporozoite rate using generalized linear mixed effects models.


We collected a total of 8123 mosquitoes, of which 1664 were Anopheles spp. malaria vectors over 295 household spray catches. We found that vector household abundance was mainly driven by the number of householders (P < 0.05), goats/sheep tethered around the house (P < 0.05) and ITNs, which halved mosquito abundance (P < 0.05). In general, similar patterns were observed for Anopheles arabiensis, but not An. gambiae s.s. and An. funestus s.s., whose density did not increase with the presence of livestock animals. Feeding on humans significantly increased in all species with the number of householders (P < 0.05), and only significantly decreased for An. arabiensis in the presence of cattle (P < 0.05). Only 26 Anopheles spp. vectors had malaria sporozoites with the sporozoite rate significantly decreasing as the proportion of cattle feeding mosquitoes increased (P < 0.05).


Our data suggest that cattle, in settings with large ITN coverage, have the potential to drive an unexpected “push-pull” malaria control system, where An. arabiensis mosquitoes “pushed” out of human contact by ITNs are likely being further “pulled” by cattle.


The dominant Anopheles spp. (Diptera: Culicidae) malaria vectors in East Africa are An. gambiae s.s. Giles and An. funestus s.s. Giles, which are regarded as highly anthropophagic (and endophagic), and An. arabiensis Patton, which is regarded as zoophagic (and exophagic) [1]. Insecticide treated nets (ITNs) have significantly improved malaria control in endemic areas such as Sub-Saharan Africa (SSA) [2], but also elsewhere, e.g., Vanuatu in Oceania [3, 4]. ITNs alone have significantly reduced morbidity and mortality due to malaria in SSA [5], a fact that is probably connected with the observed decline of vector populations [6], especially An. gambiae s.s.[7, 8]. Moreover, the larger the coverage, i.e., the more individuals sleeping under ITNs, the more effective ITNs seem to be, mainly because of the emergence of community effects [9, 10].

However, the widespread emergence of insecticide resistance in African malaria vectors may hamper malaria control programs using ITNs [11, 12]. Shifting biting hours and locations by An. gambiae s.s. after an increase of ITN coverage has become a strong concern [13, 14]. This species is also known to become opportunistic by feeding on abundant hosts over “innately” preferred ones after an increase of ITNs [15]. Since ITNs are partially inefficient in the control of exophagic vectors like An. arabiensis and early-feeding vectors like An. rivulorum Leeson [16, 17], these new phenomena call for integration of ITNs with a more robust, multi-faceted malaria control strategy [18].

ITNs are the major malaria control tool in Lake Victoria basin (LVB), western Kenya [19]. Nevertheless, despite a high ITN coverage in LVB [20], malaria prevalence remains hyperendemic (over 40% prevalence) [21]. This high prevalence could have been influenced by mosquito-mediated transmission happening outside households protected by ITNs [2224]. The non-domiciliary transmission may ultimately reflect patterns of entomological risk for malaria transmission that could be further driven by the heterogeneity of house Anopheles spp. abundance [2528].

An important factor likely underpinning adult mosquito abundance heterogeneity is their blood-feeding behavior [29]. For example, An. arabiensis, a species whose bloodmeals are not strongly biased towards humans [17], is currently the most abundant vector species within the An. gambiae complex in LVB, a trend that followed ITN use in western Kenya [7]. Within the limited distance that mosquitoes travel for bloodfeeding [30, 31], mosquitoes may be attracted by any available vertebrate host [32], especially if their feeding is not biased [33]. Livestock breeding is a major economic activity in rural East Africa [34], providing abundant hosts that can serve as an alternative bloodmeal source for malaria vectors, especially as humans become unavailable blood sources because of ITN use [17, 35]. More than a century ago, this phenomenon was called “zooprophylaxis” [36, 37] a malaria control strategy where transmission is interrupted by attracting mosquitoes to dead-end hosts [38, 39], and which has been long recommended by WHO as a protective measure against malaria [40].

Here we explored the potential of zooprophylaxis as an additional tool that could improve malaria control in endemic areas of East Africa under large ITN coverage. We hypothesized that zooprophylaxis can enhance the impacts of ITNs on malaria control by acting in a similar way to the “pull” component of a “push-pull” integrated pest management system [41]. Briefly, we consider that ITNs have worked as the “push” in malaria control, as supported by historical trends of decreased household mosquito abundance following the expansion of ITN coverage [68]. Then alternative hosts, livestock, kept at a certain distance from the households could serve as a “pull” to both keep mosquitoes away from humans and to waste their infective bites in dead end-hosts (Figure 1). We tested our hypothesis in three villages along the Kenyan shore of Lake Victoria, with a high ITN coverage.

Figure 1

Insecticide Treated Net (ITN) – Zooprophylaxis as a “Push – Pull” integrated malaria control strategy. Mosquitoes entering the houses will either be killed or “pushed” out by ITNs, then they may be pulled further away by the presence of alternative hosts such as livestock.


Study area

Our study site (Figure 2) comprised three villages, totaling 455 households over a 7.7 km2 surface in Lake Victoria shore, Mbita district, Nyanza province, Kenya (0°28’S and 34°11’E at the approximate center of the study area). Mbita’s rainfall pattern is bimodal, with a long (March - May) and a short (November - December) rainy seasons, all other months being relatively dry. Regardless of the season, the coastal lake environment maintains a high malaria vector abundance across seasons [42]. The high vector abundance in this area seems to drive year-round malaria transmission, with over 40% of the human population harboring Plasmodium falciparum (Welch) malaria parasites [21, 43]. Most houses are built using a stick framework, which is plastered with a mixture of mud and cow dung and commonly covered with a corrugated iron roof or, in few instances, by thatched roofs [44]. Most residents belong to the Luo ethnic group, and depend on small-scale farming, fishing and livestock herding for subsistence. Cattle, goats and sheep are the main livestock species herd in the study area, while domestic animals are mainly chickens and dogs. The Luos sleep inside their houses with all domestic animals and livestock outside the house during the night-time, i.e., the time when Anopheles spp. malaria vectors commonly bloodfeed [29].

Figure 2

Study Site. This map shows the location of our study site in Kenya and the location of all the houses enrolled in our study and other neighboring houses enrolled in the Nagasaki University Health Demographic Survey System.

In 2006, Kenya’s National Malaria Control Programme began to introduce ITNs in our study area through government health facilities [45]. ITNs remain as the main malaria control tool in the area [46]. Currently, there is an estimated 80% ITN household coverage, where a household is considered to be ITN covered if there is at least one ITN for each two residents [44].

Household survey

In May 2009, out of the 455 households in the study area, we chose 104 that were: (i) enrolled in the Health Demographic Surveillance System (HDSS) of Nagasaki University and Kenya Medical Research Institute and (ii) had iron or thatched roofs and eaves, and (iii) household heads of the houses were livestock breeders or took care of livestock belonging to people outside the study area (iv) where household heads provided informed consent for mosquito collection inside the houses. (ii) and (iii) were selection criteria included to minimize differences in socio-economic status among households, provided the relevant role that income differences have had in previous comprehensive studies on zooprophylaxis in Africa [47, 48]. For the subsequent analysis, the floor area (m2) of each house was estimated using a metric tape. We used the floor area as a measurement of house size.

Livestock census

We performed a census of all livestock in the study area, in order to estimate the availability of potentially zooprophilatic hosts. Specifically, we counted the number of cattle, goats and sheep, and georeferenced the night-time tethering location of each herd using a hand-held global positioning system (GPS; GPSmap 60CSx, Garmin International Inc., Olathe, USA). For large herds, those where animals were tethered in groups separated by more than 20 m, we separately recorded the location of each group. For the analysis we considered the combined number of goats and sheep, because they are herd and tether together. We also considered calf, juvenile goat and sheep to be equivalent to half adult, based on domestic animal traits (e.g., body surface area, body mass and CO2) that are relevant for mosquito bloodfeeding [49].

Ethical clearance

This study was approved by the Ethical Committee of The Graduate School of International Health Development, Nagasaki University, and Kenya Medical Research Institute (KEMRI: SSC No. 1310). Informed consent was obtained from all heads of households after the study was explained in the local language.

Mosquito larval habitat survey

Simultaneous to the household survey and livestock census, we examined occurrence of anopheline larvae in potential permanent and ephemeral habitats throughout the study area [42]. Two co-authors searched for anopheline larvae at each potential habitat for 10 minutes using standard mosquito dippers (350 ml; BioQuip Products, Rancho Dominguez, USA). A location of each positive habitat was recorded using a GPS. The locations of permanent habitats were obtained from a previous study [42]. We considered “permanent” habitats swamps or lagoons occurring along Lake Victoria shore and “ephemeral” as those generated by the impact of rainfall on the study area landscape topography. We also want to clarify that larval habitat separation into permanent and ephemeral could be relevant to understand differences in the potential recruitment of adult mosquito vectors in the area. While An. arabiensis and An. gambiae s.s. can thrive mainly in ephemeral habitat, An funestus s.s. is only successful at colonizing stable water bodies with aquatic vegetation [42].

Adult mosquito sampling and ITN use patterns

Following the mosquito larval habitat survey, adult mosquito abundance and ITN use in the 104 study houses were surveyed once per month from May to July 2009. Indoor-resting mosquitoes were sampled using pyrethrum spray catches (PSC). For each PSC we laid down white sheets on the house floor and sprayed permethrin 0.5%, a synthetic pyrethroid insecticide (Doom©, Mortain Inc, Sharjah, UAE), inside the houses, from outside through the house eaves, a procedure done to minimize the possibility of mosquito escape via the eaves. After 10 minutes knocked down mosquitoes were collected from the white sheets and transported with ice to the Nagasaki University laboratory at ICIPE (International Centre of Insect Physiology and Ecology, Mbita, Kenya). The PSC was performed during the early morning hours (4:30 ~ 6:00 hours), by a team including the first author and 3 ~ 5 local assistants.

We estimated ITN use by counting the number of people who slept in each household and the number of people who slept under an ITN by direct observation during the mosquito survey. For subsequent analyses we considered people older than 15 years as adults, between 5 and 15 years as children and below 5 years as infants. For the statistical analysis, we considered children as 1/2 of an adult, and infants as 1/3 of an adult [50]. On the same day, a co-author, not involved in the direct observations on ITN use, interviewed household heads about ITN use following the adult mosquito sampling (i.e., after 10 am). Since we found no major differences between interviews and direct observations (assuming direct observation as a gold standard, the sensitivity of the interviews was 0.93, the specificity was 0.85 and the kappa coefficient for agreement between the two methods was 0.69, for further details see Additional file 1: Supplement S1). Thus, given the intrusive nature of waking up householders to check their ITN use, during June and July, ITN use was assessed via interviews, and for the statistical analysis, we only employed the data from the interviews.

Mosquito species identification, malaria infection diagnostic and bloodmeal identification

Sampled adult mosquitoes were killed in a freezer (−20°C) in Mbita after their collection. Anopheline mosquitoes were sorted into An. gambiae s.l. and An. funestus s.l. using the morphological criteria of Gillies and Coetzee [51]. Female mosquitoes were classified as fed, gravid, or unfed by examining their abdomen under a dissecting microscope. Chilled mosquitoes were transported to Nairobi, and upon arrival kept at −40°C to avoid the digestion of bloodmeals. In Nairobi (Nagasaki University laboratory at Kenya Medical Research Institute, KEMRI), mosquito samples were dissected transversely between the thorax and the abdomen, and legs were removed.

The head and thorax of each female mosquito were tested for malaria sporozoite infection using a P. falciparum circumsporozoite protein enzyme-linked immunosorbent assay (ELISA) [52]. DNA extracted from legs of each female mosquito [53] was employed for PCR based species identification within An. gambiae s.l.[54] and An. funestus s.l.[55]. Briefly, the PCR technique allows the identification of mosquitoes morphologically identified as An.gambiae s.l. into An.gambiae s.s. or An. arabiensis. Similarly, An.funestus s.l. mosquitoes can be separated into An. funestus s.s., An rivulorum, An. leesoni Evans, or An. parensis Gillies. Mosquitoes whose morphological and molecular identification was impossible were categorized as “other Anopheles mosquitoes”.

The abdomen of each blood-fed female mosquito was used for blood meal identification by ELISA and PCR. We employed these two techniques in order to use results from the most sensitive technique, i.e., the technique able to identify most of the bloodmeals in our mosquito samples, for the subsequent statistical analysis. Our ELISA was based in the protocol developed by Beier et al. [56]. Briefly, each female mosquito abdomen was homogenized in 1 ml of phosphate-buffered saline (PBS; pH7.4). Then, we identified blood meals, humans, cattle, goats/sheep, dogs or chicken, using IgG peroxidase antibodies (Sigma-Aldrich, Steinheim, Germany). For the tests, serum from a target species was employed as a positive control, while sera from the other species were negative controls. A reaction was positive when its absorbance was at least 2 times above the mean absorbance of the highest cross-reacting serum from hosts different from the target species.

For our PCR bloodmeal identification we also employed the PBS abdomen homogenate. From the homogenate we extracted DNA following the method described in Collins et al. [53]. This was also done to have positive controls for the multiplex PCR blood-meal identification protocol developed by Kent and Norris [57]. Briefly, we employed cytochrome-B primers for human (Human741F, ggcttacttctcttcattctctcct), cattle (Cow121F, catcggcacaaatttagtcg), dog (Dog368F, ggaattgtactattattcgcaaccat), a common primer for goat and sheep (Goat894F, cctaatcttagtacttgtacccttcctc) [58], and the universal reverse primer (UNREV1025, ggttgtcctccaattcatgtta) [57]. The PCR employed 10 ng of extracted DNA from the samples using an AccuPower™ Premix (Bioneer, Daejeon, Korea) under the following amplification conditions: 95°C for 5 min; 35 cycles of template denaturing at 95°C for 1 min, primer annealing at 54°C for 1 min, and amplicon extension at 72°C for 1 min; and a final extension at 72°C for 7 min. Each animal’s blood was detected by agarose gel electrophoresis (2% TAE) as 334 bp (Human) and 561 bp (Cow), 132 bp (Goat/Sheep), 680 bp (Dog) bands, respectively.

Data analysis

We employed maximum likelihood statistical models to quantify the impacts of ITN use, presence of larval habitats, human and livestock abundance on: (i) mosquito abundance, (ii) mosquito bloodfeeding on humans and (iii) mosquito sporozoite rate. In each case the analysis began with a full model that considered all relevant covariates that we measured, which then was simplified until the Akaike Information Criterion (AIC) was minimized following the stepwise removal of covariates [59]. For (i) and (ii), we performed separate analyses for each dominant malaria vector species in the area, i.e., An. arabiensis, An. gambiae s.s. and An. funestus s.s., and also for all the Anopheles spp. vectors combined. For (iii) we could only analyze An. funestus s.s. and the combination of all vector species we collected. All statistical analyses were performed with the Statistical Package R version 3.0, employing the library lme4 to fit generalized linear mixed effect models [60].

Mosquito abundance: for this analysis we employed Poisson generalized linear mixed models (Pois-GLMM). The response (i.e., independent variable) was the total number of mosquitoes, i.e., including gravid, unfed and bloodfed, collected by each PSC effort, i.e., per household and collection date. Fixed factors in the full model included: number of ITNs in use, distance to the closest permanent larval habitat from each household, abundance of ephemeral larval habitats in a given buffer area around each household, house area (as a proxy for adult mosquito resting habitat size), the month when each sample was collected (to control for seasonal mosquito abundance variability), the adjusted abundance of human residents in each household, cattle and goats/sheep in a given radius around each household, and number of neighboring households in a given radius around each household to account for the impact that humans living in the neighborhood could have in attracting mosquitoes to the focal house. Since some permanent larval habitats were too large to quantify in a given radius, we considered a distance from each household to the closest habitat as a proxy of their magnitude on potential adult mosquito productivity.

The radius, within which the abundance of ephemeral larval habitats, cattle, combined goats/sheep and neighboring households were counted, was chosen prior to the specification of the full model by a two step process. First we counted their abundance within a set of concentric circular areas with the following radii (with each radius originating from the center of each focal household): 20 m, 50 m, 100 m, 150 m, 200 m, 250 m, 300 m, 350 m, 400 m, 450 m, 500 m; and second we chose the best radius based on AIC minimization for Pois-GLMMs only considering mosquito abundance as a function of each factor in a given radius. We needed to perform this selection given the collinearity that emerges by considering the different radii simultaneously, since abundance in each radius is a linear function of a nested shorter radius, and collinear covariates can cause unidentifiability in parameter estimation by GLMMs [59]. The radii range was chosen based on the distances considered in previous zooprophylaxis studies in Russia [40], Pakistan [61], and western Kenya [34], to better account for the dispersal of anopheline mosquitoes, which on average is around 200 ~ 400 m [62]. Random factors included: the collection date to account for potential lack of temporal independence in our samples and the household identity to account for the potential lack of spatial independence [63].

Human bloodfeeding: for this analysis we employed binomial generalized linear mixed models (Bin-GLMMs). The response was the odds of mosquitoes having human bloodmeals versus bloodmeals from other hosts in each household and collection date. The odds of mosquitoes having human bloodmeals is defined as the ratio of the probability of mosquitoes having human bloodmeals divided by the probability of mosquitoes having a non-human bloodmeal. To compute the probability of human bloodmeals, we considered mosquitoes whose bloodmeal came exclusively from humans or humans and a second host. For the probability of a different host bloodmeal, we considered mosquitoes whose bloodmeal was identified as coming from a livestock species (cattle or goats/sheep) or was unidentified by the primers we employed. The full model included the same fixed and random factors used in the full Pois-GLMM used to study mosquito abundance.

Mosquito sporozoite rate: for this analysis we employed Poisson rate generalized linear mixed models (PoisR-GLMM). The response was the number of sporozoite infected mosquitoes. Briefly, the rate model considers the maximum number of sporozoite positive mosquitoes to be constrained by the number of collected mosquitoes that could have been infected (gravid or bloodfed mosquitoes are more likely infected than unfed mosquitoes) [64]. By contrast, a Pois-GLMM assumes this number as unconstrained [59]. As covariates we considered bednet coverage (i.e., the proportion of people in a household sleeping under an ITN) to control for differences in access to human hosts, the distance to permanent mosquito larval habitats, the abundance of ephemeral mosquito larval habitats in a given radius, house size and collection month. For the same reasons, we included these variables when analyzing mosquito abundance. We also considered cattle and goats/sheep blood indices, which are the proportion of mosquitoes with bloodmeals from each of these two hosts (excluding mixed meals with other hosts) divided by the total number of bloodfed mosquitoes. We employed the blood feeding indices as they allowed us to test if differential increases in the proportion of mosquitoes feeding on hosts other than humans decreases sporozoite infection (sporozoite rate of the vectors), or if other hosts provide a zooprophylactic effect.


Mapping of livestock tethering points, mosquito larval habitats and ITN use

A total of 100 herds were present, comprising 850 cattle heads (adult: 746, calf: 104, adjusted abundance by size: 798.0) (Figure 3A). For goats/sheep, a total of 132 herds were present, and were comprised of 1301 goats (adult: 1094, young: 207, adjusted abundance by size: 1199.5), and 82 sheep (adult: 71, young: 11, adjusted abundance by size: 76.5) (Figure 3B).

Figure 3

Location, in relation to the study houses, of: tethering points for (A) Cattle (B) Goats and Sheep (C) Anophele s spp larval Habitats (D) Insecticide Treated Net (ITN) use. For further details see the inset legends.

All major permanent larval habitats (7 sites) with aquatic vegetation were located in Lake Victoria Coast, and a total of 173 ephemeral larval habitats were found (Figure 3C). Nearly 70% (118 sites) of the ephemeral habitats in the study area were natural pools such as puddles created in gullies. The remaining ones (55 sites) were mainly ditches or drainage channels associated with farms and pit holes.

Due to resident absence, of the 104 houses enrolled in the study, only 95 houses were surveyed in May, 100 in June and 100 in July. The average floor area of the houses was 17.37 ± 6.40 m2 (mean ± S.D.). The average number of people sleeping in a household was 4.19 ± 1.95 (mean ± S.D.), with minimum variation across the months of our study (Figure 3D). The average number of ITNs in use per house was 1.48 ± 0.63 (mean ± S.D.) with non-significant variation through our study, and the average number of residents sleeping under ITNs per house was 3.47 ± 1.90 (mean ± S.D.). Throughout our study, an average of 82% of the people sleeping in our study houses used ITNs. There were some heterogeneities regarding the use by age group. About 90% of the adults were covered, but only 65% of the children and 83% of the infants slept under ITNs.

Mosquitoes, bloodmeal hosts and sporozoite rate

In total we collected 8123 mosquitoes over 295 sessions of PSC, and 1664 of them were female mosquitoes belonging to either An. gambiae s.l. or An. funestus s.l. (811 in May, 416 in June, and 437 in July). According to the morphological species identification, 928 Anopheles mosquitoes were An. gambiae s.l., which, using PCR were further identified as: 726 An. arabiensis, 196 An. gambiae s.s.. Belonging to An. funestus s.l., we morphologically identified 736 individuals, with 711 being An. funestus s.s., and 1 An. rivulorum. We could not molecularly identify 6 An. gambiae s.l. and 24 An. funestus s.l. with the PCR method. The abundance of Anopheles vectors was nearly constant during the study period (Figure 4A). Anopheles arabiensis was mainly present in the early months of our study (Figure 4B), and An. funestus s.s. in July (Figure 4D). Meanwhile, An. gambiae s.s. had its maximum abundance in May, almost disappearing in July (Figure 4C).

Figure 4

Mosquito Abundance per household spray catch and month (A) All Anopheles spp (B) Anopheles arabiensis (C) Anopheles gambiae s.s. (D) Anopheles funestus s.s..

We tried to identify host bloodmeals in all the 1204 fed Anopheles vectors, using the Kent & Norris multiplex PCR (Additional file 2: Table S1) and ELISA (Additional file 3: Table S2). With the PCR method we identified bloodmeal hosts for 1031 (85.63%) samples, while the ELISA only allowed the identification of blood sources in 846 (70.27%) samples. Thus, the rest of our results and analysis will be exclusively based on the PCR data. Main vectors where blood meals were identified were: An. arabiensis (n = 440), An. gambiae s.s. (n = 110), An. funestus s.s. (n = 461). Most bloodmeals for An. gambiae s.s. and An.funestus s.s. came from humans (See Additional file 2: Table S1). The proportions of human bloodmeal for An. gambiae s.s., An. funestus s.s. and An. arabiensis were 49.5% (CI: 42.3-56.7), 63.9% (CI: 60.2-67.4) and 12.0% (CI: 9.7-14.6), respectively, and the proportion for An. arabiensis was significantly less compared with the other species. Meanwhile, the major blood source was cattle for An. arabiensis. The proportions of cattle bloodmeal for An. gambiae s.s., An. funestus s.s. and An. arabiensis were 1.2% (CI: 0.6-2.3), 2.0% (CI: 0.6-5.1) and 40.5% (CI: 36.9-44.2), respectively, and the proportion for An. arabiensis was significantly greater.

Two An. arabiensis, three An.gambiae s.s. and 21 An. funestus s.s. had P. falciparum sporozoites. Sporozoite rate was 0.28% (CI: 0.03-0.99) for An. arabiensis, 1.53% (CI: 0.32-4.40) for An. gambiae s.s. and 2.95% (CI: 1.83-4.48) for An. funestus s.s., and the rate for An. funestus s.s. was significantly higher compared with An. arabiensis.

Statistical models

For all models we indicated the number of observations in the different tables. We ignored factors that showed no statistical significance (P > 0.05) in the description of the best models results, since, in principle, the impacts of those factors are not different from what is expected by random. Here, we also want to mention that for some models spatial or temporal variability was not significant, thus explaining their absence in the final models and during the process of model selection. All the assumptions of the statistical models were not violated.

The impacts of ITN use, presence of larval habitats, human and livestock abundance on mosquito abundance

After a process of model selection (Additional file 4: Table S3), we found that the best model explaining the abundance of all the Anopheles vectors included ITN use, goats/sheep abundance in 20 m around a household, and the adjusted number of residents in a household (Table 1). We found the number of mosquitoes per PSC session, increased 9% by each additional goat or sheep tethered within a 20 m radius from a household, and 36% per each additional person sleeping in a household. By contrast mosquito abundance was halved by each ITN use. We also found a higher spatial variability (Household variance in Table 1) than temporal variability (Collection date variance in Table 1).

Table 1 Parameter estimates for the best Poisson GLMM explaining the abundance of Anopheles vectors

We also performed a process of model selection for each one of the dominant vector species (Additional file 5: Table S4). We found that An. arabiensis abundance was increased by 10% with each additional goats/sheep tethered within 20 m of a household, 6% by each square meter of household basal area, that its abundance increased around 10 times in May when compared with July, while ITNs nearly halved its abundance, and each neighbouring household within a 50 m radius decreased mosquito abundance by 5% (Table 2). The impacts of ITNs on An. gambiae s.s. abundance were stronger. For each ITN, mosquito abundance decreased by 58%, its abundance also decreased with the number of neighbouring houses, and it increased by 36% per each additional household resident (Table 2). Further, the abundance was increased by 92% in May when compared with July. As observed in the raw data, An. funestus s.s. abundance was maximum in July. Each additional resident increased its abundance by 68%, while each ITN in use reduced it by 59% (Table 2). Interestingly, An. funestus s.s. abundance increased by 3% with each additional ephemeral habitat in a 500 m radius around the house. For all the dominant vector species, spatial variability (Household variance in Table 2) was larger than temporal variability (Date variance in Table 2).

Table 2 Parameter estimates for the best Poisson GLMM explaining the abundance of three vector species

The impacts of ITN use, presence of larval habitats, human and livestock abundance on mosquito human bloodfeeding

Model selection (Additional file 6: Table S5) showed that when all Anopheles vector species were combined, bloodfeeding on humans was significantly associated with the number of people sleeping in a household, which increased the odds of an Anopheles mosquito feeding on humans by 1.53 times, and odds of human bloodfeeding were decreased 0.99 times by each goat or sheep tethered within 500 m from the household (Table 3). By contrast odds of human bloodfeeding in May was decreased 0.20 times when compared with July, and additional square meter of household basal area decreased it 0.96 times.

Table 3 Parameter estimates for the best binomial GLMM explaining feeding on humans by Anopheles vectors

Regarding each dominant vector species in the area, we found that factors influencing human bloodfeeding were different (Additional file 7: Table S6). We found that An. gambiae s.s. was significantly sensitive only to the presence of humans, where each additional human increased human bloodfeeding odds two times (Table 4). A similar pattern was also observed for An. funestus s.s.. Its human bloodfeeding significantly increased 1.4 times with each additional person sleeping in a house, and was 4.4 times more likely to happen in June (Table 4). By contrast, human bloodfeeding in An. funestus s.s. decreased 0.93 times with each additional ephemeral habitat in a 200 m radius around a house. While human bloodfeeding in An. arabienesis was 1.54 times most likely by the presence of each human in the household, it also increased 1.01 times for each goats/sheep tethered within 500 m from a household and decreased with the number of cattle heads tethered in the same radius as the goats/sheep (Table 4). Anopheles arabiensis feeding on humans was 1.07 times more likely to occur as the number of houses in 150 m around a focal household increased and human bloodfeeding was 0.69 times less likely to occur during May compared with July. The temporal variances for An. arabiensis and An. gambiae s.s. were around 0.30 (Date variance in Table 4) which, respectively, represent 0.23% and 0.36% of the total deviance (An. arabiensis deviance = 129.7 and An gambiae s.s. deviance =83.1). The spatial variance for An. funestus s.s. was very low (Household variance in Table 4) representing only 0.25% of the total deviance (An funestus s.s. deviance =121.0).

Table 4 Parameter estimates for the best Binomial GLMM explaining feeding on humans by three vector species

The impact of ITN use, presence of larval habitats, human and livestock abundance on mosquito sporozoite rates

The sporozoite rate in all Anopheles vector species combined was significantly associated with the number of ephemeral habitats in 250 m radius around the houses, where each habitat increased the sporozoite rate by 13%. In contrast the rate reduced by 9.2% as the proportion of mosquitoes feeding on cattle increased by 10% (Table 5, for variable selection see Additional file 8: Table S7). Regarding An. funestus s.s., the species with the largest number of sporozoite positive individuals, we found that sporozoite rate increased by 15% with each ephemeral habitat in a 250 m radius around the house (Table 5, for variable selection see Additional file 9: Table S8). The random factors had a very low variability, less than 0.1% of the total deviance (Table 5).

Table 5 Parameter estimates for Poisson rate GLMM explaining the sporozoite rate of malaria vectors


The role of zooprophylaxis, more generally of bloodmeal source diversity, on vector-borne disease has been controversial since early times [32, 37, 47, 65, 66]. An outstanding conclusion from field studies [32, 47, 48] and mathematical modeling [38, 39, 67, 68] is that diversion of feeding to alternative hosts needs to outweigh the potential increase in mosquito population size and attraction to humans by “zooprophylactic” animals [31]. To be more specific, mosquitoes may increase their contact with humans after they are attracted to the proximity of a household, or mosquito population size may increase, because of increases in livestock presence and hoofprint habitats for mosquito larvae [69, 70]. All these factors seem to critically rely on the location and abundance of potential alternative hosts.

Nevertheless, previous field studies with negative results for zooprophylaxis neglected the exact location of overnight tethering points for livestock. They mainly looked at heterogeneities in livestock ownership via interviews [47, 48, 61]. Moreover, the information gathered by the interviews may provide imprecise information about alternative host availability, for example, the cattle belonging to a household is not necessarily the one tethered close to it. To overcome this limitation, we performed a census of all livestock present in our study area, georeferencing overnight tethering points, herd size and the age structure of herds.

Zooprophylaxis in the community of Malaria vector species

Unlike a previous comprehensive study on passive zooprophylaxis in West Africa [47, 48], we found that increases in bloodfeeding in cattle can reduce the likelihood of vector infection, i.e., they might render a zooprophylactic effect. Thus, our data suggest, that for the overall community of malaria vectors, cattle presence might be acting as a “pulling” passive malaria control strategy in areas where high ITN coverage have “pushed” changes in mosquito access to humans.

In our study area, although the abundance of goats/sheep tethered in close proximity (20 m or less) to a household is associated with an increase in the number of Anopheles mosquitoes that can be caught in such house, a phenomenon already described [8], the main driver of domiciliary mosquito abundance seems to be the number of people sleeping in a house. Nevertheless, it is well documented that people sleeping under ITNs have a significant reduction in their malaria infection risk [71], because, as shown by this study and others, the number of mosquitoes is significantly smaller in houses where people use ITNs [7, 72]. In fact, our parameter estimates suggest that ITN use can outweigh the increase in mosquito density by the presence of one human and up to 4 goats/sheep. For example, if we estimate the proportional change in mosquito abundance as a function of these factors ((1.09)4 × 1.36 × 0.48 = 0.92, Table 1), we still can expect an 8% reduction in mosquito abundance because of ITN use. If we consider the subset of bloodfed Anopheles mosquitoes, we also see that odds of bloodfeeding on humans are primarily associated with human abundance and ITN use. Similar phenomena have been observed in previous studies about bloodfeeding in African settings with ITNs [14, 73].

Nevertheless, the likelihood of vector infection can also be increased by the presence of ephemeral habitats. Since most of the malaria infected mosquitoes were An. funestus s.s., whose malaria infection was only increased in the presence of ephemeral habitats, the association of ephemeral habitats can be related to two factors. The first is that, An. funestus s.s., whose immature stages mainly develop in permanent swamps along Lake Victoria shore [42], may be blown into the area where infected mosquitoes were found. We have observed that Anopheles mosquitoes emerge in the evening when the west wind from the lake is strong, and then given the small body size and relative poor flight ability of this vector species, they may get stuck in the area. In fact, several infected mosquitoes were found in houses east of two large permanent habitats (Figure 3C and Figure 5). This explanation may also be applied to the association of An. funestus s.s. abundance with ephemeral habitats that were also recorded east of the permanent habitats (Figure 3C).

Figure 5

House location and abundance of sporozoite positive Anopheles spp mosquitoes. For further details see the inset legend.

The second factor, although uncommon, An. funestus s.s. larvae can colonize what we classified as ephemeral habitats [74], probably reflects some sub-optimal oviposition habitat selection, a phenomenon commonly observed in mosquitoes [75, 76]. Furthermore, the suboptimal oviposition might be also related to the limited dispersal ability of An. funestus s.s. mosquitoes and the fitness costs of not ovipositing when it is preferable to do so in a suboptimal habitat and bloodfeed again. This can increase malaria transmission [67, 68] and might be common, since it has been observed in other mosquito vectors [77]. Although this explanation seems counterintuitive given that human bloodfeeding in An. funestus s.s., slightly decreased with an increase of ephemeral habitats in a 200 m radius, infection still can happen as long as mosquitoes bloodfeed on infected humans.

Zooprophylaxis on the dominant malaria vectors

Regarding zooprophylactic patterns in each dominant vector species, we found that there were species-specific differences. Anopheles gambiae s.s., the less common species after ITN coverage has been scaled in the study area [7], had abundance patterns that were insensitive to livestock abundance, and its feeding on humans decreased (yet not significantly) in the presence of goats and sheep [78]. Nevertheless, its household abundance greatly decreased in the presence of ITNs, and the number of sporozoite positive An. gambiae s.s. was only three, meaning that ITNs might be controlling malaria transmission by this species, as suggested by previous studies [7].

Although not statistically significant in the best model selected by AIC, the presence of cattle and goats/sheep reduced the odds of human bloodfeeding in An. funestus s.s. by magnitudes (8% and 1%, respectively) that prevail over the potential positive impacts of cattle on malaria transmission. This species accounted for most of the sporozoite infected mosquitoes and infections might reflect differential patterns of human exposure to this species because of human activities.

Anopheles arabiensis, showed the most complex relationships with livestock abundance, first it increased about 10% for each goat or sheep tethered in close proximity (20 m or less) to the household. This species also increased its bloodfeeding on humans with the overall numbers of goats/sheep in the landscape (500 m or less). Goats and sheep often spend time outside house walls under extended eaves. They may release semiochemicals, which could be associated with the “goat” fragrance familiar to anybody working with goats and sheep, and goat associated semiochemicals could potentially attract this vector close to the houses. Nevertheless, its feeding on humans decreased with cattle abundance in the landscape (500 m or less), which had no impact on its abundance. The positive impacts of cattle might reflect the low number of mosquitoes with malaria sporozoites. In synthesis, cattle seem to provide a zooprophylactic effect in An. arabiensis and An. funestus s.s.. Goats/sheep abundance is nevertheless associated with an increase in An. arabiensis abundance and likelihood of feeding on humans.

Study limitations

The limited span of our study can only provide inferences about the rainy season. Ideally, a longer study would provide a better picture of the magnitude of zooprophylactic effects by different livestock species on the community of malaria vectors year-round. The increased sampling could also improve the power of our estimates, for example, telling if the non-significance of factors in the models selected as best were artifacts of the sample size, and to better understand any impacts of weather seasonality on feeding behavior and abundance of household resting mosquitoes [79].

Our sampling method only focused on collecting indoor resting mosquitoes by PSC without collecting mosquitoes outside households. We could have underestimated the number of mosquitoes, especially the species that have a preference to feed indoors (endophagy) on humans (anthropophagy), but mainly rest outdoors (exophily) [17]. Moreover, our sampling method might totally ignore exophilic and exophagic species. Many of those exo-exo mosquito species are vectorially competent (i.e., get infected and develop sporozoites) for human malaria parasites and might be responsible for some malaria transmission in the area.

We could also have collected additional data on income, similar to previous studies [48], or measure specific details about house conditions such as eave openness, which affect malaria vector household entrance behavior [80]; that is, we limited our inferences to a relatively homogenous group of households.

Perspectives and recommendations for future studies

Beyond improving the limitations of our study, we believe there is also room to explore the origin of unidentified bloodmeals, and the impact of the bloodmeal source hosts in zooprophylaxis. For example, Kawada et al. [17] found An. rivolorum feeding on hippopotamus, a common vertebrate in Lake Victoria, suggesting that malaria vectors can have unappreciated hosts [78]. These hosts may contribute to zooprophylaxis, and the impacts of unappreciated hosts could be related to their diversity, a hypothesis originally envisioned by Celli when talking about zooprophylaxis [32, 36] and refereed as “dilution effect” by some ecologists working with ticks, who have seen decreases in disease risk/burden when biodiversity increases [32, 65, 66]. Zooprophylaxis impacts can also be potentially enhanced by the use of livestock vermicides [81, 82], which can kill mosquitoes after bloodfeeding on livestock hosts. In addition, it will be worth exploring the use of insecticide collars [8385] in goats tethered around houses, which potentially can further “push” malaria vectors out of human contact.


Our data suggest that malaria control by high ITN coverage might be benefiting from zooprophylactic effects derived from the abundance of cattle in our study area. The synergy between zooprophylaxis and ITNs resulted in an integrated vector control strategy for the dominant vector species in our study area, especially for An. arabiensis. This integrated “push-pull” strategy may also be effective for reducing malaria transmission in the other parts of LVB and similar sites where An. arabiensis is a dominant vector species. Our data also suggest that tethering and keeping goats as far as possible from households could improve zooprophylaxis in An. arabiensis, whose abundance and feeding on humans increases with the density of goats tethered in close proximity to households.


  1. 1.

    Okara RM, Sinka ME, Minakawa N, Mbogo CM, Hay SI, Snow RW: Distribution of the main malaria vectors in Kenya. Malar J. 2010, 9: 69-10.1186/1475-2875-9-69.

    PubMed Central  Article  PubMed  Google Scholar 

  2. 2.

    Steketee RW, Campbell CC: Impact of national malaria control scale-up programmes in Africa: magnitude and attribution of effects. Malar J. 2010, 9: 299-10.1186/1475-2875-9-299.

    PubMed Central  Article  PubMed  Google Scholar 

  3. 3.

    Kaneko A, Taleo G, Kalkoa M, Yamar S, Kobayakawa T, Björkman A: Malaria eradication on islands. Lancet. 2000, 356 (9241): 1560-1564. 10.1016/S0140-6736(00)03127-5.

    CAS  Article  PubMed  Google Scholar 

  4. 4.

    Chaves LF, Kaneko A, Taleo G, Pascual M, Wilson ML: Malaria transmission pattern resilience to climatic variability is mediated by insecticide-treated nets. Malar J. 2008, 7: 100-10.1186/1475-2875-7-100.

    PubMed Central  Article  PubMed  Google Scholar 

  5. 5.

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

    Google Scholar 

  6. 6.

    Magesa SM, Wilkes TJ, Mnzava AE, Njunwa KJ, Myamba J, Kivuyo MD, Hill N, Lines JD, Curtis CF: Trial of pyrethroid impregnated bednets in an area of Tanzania holoendemic for malaria. Part 2. Effects on the malaria vector population. Acta Trop. 1991, 49: 97-108. 10.1016/0001-706X(91)90057-Q.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    Bayoh MN, Mathias DK, Odiere MR, Mutuku FM, Kamau L, Gimnig JE, Vulule JM, Hawley WA, Hamel MJ, Walker ED: Anopheles gambiae: historical population decline associated with regional distribution of insecticide-treated bed nets in western Nyanza Province, Kenya. Malar J. 2010, 9: 62-10.1186/1475-2875-9-62.

    PubMed Central  Article  PubMed  Google Scholar 

  8. 8.

    Mutuku FM, King CH, Mungai P, Mbogo C, Mwangangi J, Muchiri EM, Walker ED, Kitron U: Impact of insecticide-treated bed nets on malaria transmission indices on the south coast of Kenya. Malar J. 2011, 10: 356-10.1186/1475-2875-10-356.

    PubMed Central  Article  PubMed  Google Scholar 

  9. 9.

    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 bed nets on child mortality and malaria morbidity in western Kenya. Am J Trop Med Hyg. 2003, 68 (4 Suppl): 121-127.

    PubMed  Google Scholar 

  10. 10.

    Komazawa O, Kaneko S, K’Opiyo J, Kiche I, Wanyua S, Shimada M, Karama M: Are long-lasting insecticidal nets effective for preventing childhood deaths among non-net users? A community-based cohort study in western Kenya. PLoS One. 2012, 7 (11): e49604-10.1371/journal.pone.0049604.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  11. 11.

    Kawada H, Dida GO, Ohashi K, Komagata O, Kasai S, Tomita T, Sonye G, Maekawa Y, Mwatele C, Njenga SM, Mwandawiro C, Minakawa N, Takagi M: Multimodal pyrethroid resistance in malaria vectors, Anopheles gambiae s.s., Anopheles arabiensis, and Anopheles funestus s.s. in western Kenya. PLoS One. 2011, 6 (8): e22574-10.1371/journal.pone.0022574.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  12. 12.

    Trape JF, Tall A, Diagne N, Ndiath O, Ly AB, Faye J, Dieye-Ba F, Roucher C, Bouganali C, Badiane A, Sarr FD, Mazenot C, Toure-Balde A, Raoult D, Druilhe P, Mercereau-Puijalon O, Rogier C, Sokhna C: Malaria morbidity and pyrethroid resistance after the introduction of insecticide-treated bednets and artemisinin-based combination therapies: a longitudinal study. Lancet Infect Dis. 2011, 11 (12): 925-932. 10.1016/S1473-3099(11)70194-3.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    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.

    PubMed Central  Article  PubMed  Google Scholar 

  14. 14.

    Gatton ML, Chitnis N, Churcher T, Donnelly MJ, Ghani AC, Godfray HC, Gould F, Hastings I, Marshall J, Ranson H, Rowland M, Shaman J, Lindsay SW: The importance of mosquito behavioural adaptations to malaria control in Africa. Evolution. 2013, 67 (4): 1218-1230. 10.1111/evo.12063.

    PubMed Central  Article  PubMed  Google Scholar 

  15. 15.

    Lefevre T, Gouagna LC, Dabire KR, Elguero E, Fontenille D, Renaud F, Costantini C, Thomas F: Beyond nature and nurture: phenotypic plasticity in blood-feeding behavior of Anopheles gambiae s.s. when humans are not readily accessible. Am J Trop Med Hyg. 2009, 81 (6): 1023-1029. 10.4269/ajtmh.2009.09-0124.

    Article  PubMed  Google Scholar 

  16. 16.

    Kawada H, Dida GO, Sonye G, Njenga SM, Mwandawiro C, Minakawa N: Reconsideration of Anopheles rivulorum as a vector of Plasmodium falciparum in western Kenya: some evidence from biting time, blood preference, sporozoite positive rate, and pyrethroid resistance. Parasit Vectors. 2012, 5: 230-10.1186/1756-3305-5-230.

    PubMed Central  Article  PubMed  Google Scholar 

  17. 17.

    Takken W, Verhulst NO: Host preferences of blood-feeding mosquitoes. Annu Rev Entomol. 2013, 58: 433-453. 10.1146/annurev-ento-120811-153618.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Chaves LF, Koenraadt CJM: Climate change and highland malaria: fresh air for a hot debate. Q Rev Biol. 2010, 85 (1): 27-55. 10.1086/650284.

    Article  PubMed  Google Scholar 

  19. 19.

    Githinji S, Herbst S, Kistemann T, Noor AM: Mosquito nets in a rural area of Western Kenya: ownership, use and quality. Malar J. 2010, 9: 250-10.1186/1475-2875-9-250.

    PubMed Central  Article  PubMed  Google Scholar 

  20. 20.

    Atieli HE, Zhou G, Afrane Y, Lee MC, Mwanzo I, Githeko AK, Yan G: Insecticide-treated net (ITN) ownership, usage, and malaria transmission in the highlands of western Kenya. Parasit Vectors. 2011, 4: 113-10.1186/1756-3305-4-113.

    PubMed Central  Article  PubMed  Google Scholar 

  21. 21.

    Noor AM, Gething PW, Alegana VA, Patil AP, Hay SI, Muchiri E, Juma E, Snow RW: The risks of malaria infection in Kenya in 2009. BMC Infect Dis. 2009, 9: 180-10.1186/1471-2334-9-180.

    PubMed Central  Article  PubMed  Google Scholar 

  22. 22.

    Reddy MR, Overgaard HJ, Abaga S, Reddy VP, Caccone A, Kiszewski AE, Slotman MA: Outdoor host seeking behaviour of Anopheles gambiae mosquitoes following initiation of malaria vector control on Bioko Island, Equatorial Guinea. Malar J. 2011, 10: 184-10.1186/1475-2875-10-184.

    PubMed Central  Article  PubMed  Google Scholar 

  23. 23.

    Bousema T, Griffin JT, Sauerwein RW, Smith DL, Churcher TS, Takken W, Ghani A, Drakeley C, Gosling R: Hitting hotspots: spatial targeting of malaria for control and elimination. PLoS Med. 2012, 9 (1): e1001165-10.1371/journal.pmed.1001165.

    PubMed Central  Article  PubMed  Google Scholar 

  24. 24.

    Govella NJ, Ferguson H: Why use of interventions targeting outdoor biting mosquitoes will be necessary to achieve malaria elimination. Front Physiol. 2012, 3: 199-

    PubMed Central  Article  PubMed  Google Scholar 

  25. 25.

    Mbogo CM, Mwangangi JM, Nzovu J, Gu W, Yan G, Gunter JT, Swalm C, Keating J, Regens JL, Shililu JI, Githure JI, Beier JC: Spatial and temporal heterogeneity of Anopheles mosquitoes and Plasmodium falciparum transmission along the Kenyan coast. Am J Trop Med Hyg. 2003, 68 (6): 734-742.

    PubMed  Google Scholar 

  26. 26.

    Zhou G, Munga S, Minakawa N, Githeko AK, Yan G: Spatial relationship between adult malaria vector abundance and environmental factors in western Kenya highlands. Am J Trop Med Hyg. 2007, 77 (1): 29-35.

    PubMed  Google Scholar 

  27. 27.

    Grillet ME, Barrera R, Martinez JE, Berti J, Fortin MJ: Disentangling the effect of local and global spatial variation on a mosquito-borne infection in a neotropical heterogeneous environment. Am J Trop Med Hyg. 2010, 82 (2): 194-201. 10.4269/ajtmh.2010.09-0040.

    PubMed Central  Article  PubMed  Google Scholar 

  28. 28.

    Norris LC, Norris DE: Heterogeneity and changes in inequality of malaria risk after introduction of insecticide-treated bed nets in Macha, Zambia. Am J Trop Med Hyg. 2013, 88 (4): 710-717. 10.4269/ajtmh.11-0595.

    PubMed Central  Article  PubMed  Google Scholar 

  29. 29.

    Silver JB: Mosquito Ecology: Field Sampling Methods. 2008, New York: Springer, 3

    Google Scholar 

  30. 30.

    Edman J, Day J, Walker E: Vector-Host Interplay - Factors Affecting Disease Transmission. Ecology of Mosquitoes: Proceedings of a Workshop. 1985, New Beach, Florida: Florida Medical Entomology Laboratory, 273-285.

    Google Scholar 

  31. 31.

    Edman JD: Disease Control Through Manipulation of Vector-Host Interaction: Some Historical and Evolutionary Perspectives. Proceedings of a Symposium: The Role of Vector-Host Interactions in Disease Tranmission. Edited by: Scott TW, Grumstrup-Scott J. 1988, Washington, D.C: Entomological Society of America, 43-50.

    Google Scholar 

  32. 32.

    Chaves LF, Harrington LC, Keogh CL, Nguyen AM, Kitron UD: Blood feeding patterns of mosquitoes: random or structured?. Front Zool. 2010, 7: 3-10.1186/1742-9994-7-3.

    PubMed Central  Article  PubMed  Google Scholar 

  33. 33.

    Mahande A, Mosha F, Mahande J, Kweka E: Feeding and resting behaviour of malaria vector. Anopheles arabiensis with reference to zooprophylaxis. Malar J. 2007, 6: 100-10.1186/1475-2875-6-100.

    PubMed Central  Article  PubMed  Google Scholar 

  34. 34.

    Minakawa N, Seda P, Yan G: Influence of host and larval habitat distribution on the abundance of African malaria vectors in western Kenya. Am J Trop Med Hyg. 2002, 67 (1): 32-38.

    PubMed  Google Scholar 

  35. 35.

    Port G, Boreham P: The effect of bed nets on feeding by Anopheles gambiae Giles (Diptera, Culicidae). Bull Entomol Res. 1982, 72: 483-488. 10.1017/S0007485300013663.

    Article  Google Scholar 

  36. 36.

    Celli A: La Malaria in Italia Durante il 1907. Ricerche Epidemiologiche e profilatiche. Atti della Societa per gli studi della Malaria. 1908, 9: 675-729.

    Google Scholar 

  37. 37.

    Hackett L, Missiroli A: The natural disappareance of malaria in certain regions of Europe. Am J Hyg. 1931, 13: 57-78.

    Google Scholar 

  38. 38.

    Sota T, Mogi M: Effectiveness of zooprophylaxis in malaria control: a theoretical inquiry, with a model for mosquito populations with two bloodmeal hosts. Med Vet Entomol. 1989, 3 (4): 337-345. 10.1111/j.1365-2915.1989.tb00240.x.

    CAS  Article  PubMed  Google Scholar 

  39. 39.

    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.

    PubMed Central  Article  PubMed  Google Scholar 

  40. 40.

    World Health Organization: Mannal on Environmental Management for Mosquito Control with Special Emphasis on Malaria Cectors. 1982

    Google Scholar 

  41. 41.

    Cook SM, Khan ZR, Pickett JA: The use of push-pull strategies in integrated pest management. Annu Rev Entomol. 2007, 52: 375-400. 10.1146/annurev.ento.52.110405.091407.

    CAS  Article  PubMed  Google Scholar 

  42. 42.

    Minakawa N, Dida GO, Sonye GO, Futami K, Njenga SM: Malaria vectors in Lake Victoria and adjacent habitats in western Kenya. PLoS One. 2012, 7 (3): e32725-10.1371/journal.pone.0032725.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  43. 43.

    ter Kuile FO, Terlouw DJ, Kariuki SK, Phillips-Howard PA, Mirel LB, Hawley WA, Friedman JF, Shi YP, Kolczak MS, Lal AA, Vulule JM, Nahlen BL: Impact of permethrin-treated bed nets on malaria, anemia, and growth in infants in an area of intense perennial malaria transmission in western Kenya. Am J Trop Med Hyg. 2003, 68 (4 Suppl): 68-77.

    PubMed  Google Scholar 

  44. 44.

    Iwashita H, Dida G, Futami K, Sonye G, Kaneko S, Horio M, Kawada H, Maekawa Y, Aoki Y, Minakawa N: Sleeping arrangement and house structure affect bed net use in villages along Lake Victoria. Malar J. 2010, 9: 176-10.1186/1475-2875-9-176.

    PubMed Central  Article  PubMed  Google Scholar 

  45. 45.

    Hightower A, Kiptui R, Manya A, Wolkon A, Vanden Eng JL, Hamel M, Noor A, Sharif SK, Buluma R, Vulule J, Laserson K, Slutsker L, Akhwale W: Bed net ownership in Kenya: the impact of 3.4 million free bed nets. Malar J. 2010, 9: 183-10.1186/1475-2875-9-183.

    PubMed Central  Article  PubMed  Google Scholar 

  46. 46.

    Dye TD, Apondi R, Lugada ES, Kahn JG, Smith J, Othoro C: “Before we used to get sick all the time”: perceptions of malaria and use of long-lasting insecticide-treated bed nets (LLINs) in a rural Kenyan community. Malar J. 2010, 9: 345-10.1186/1475-2875-9-345.

    PubMed Central  Article  PubMed  Google Scholar 

  47. 47.

    Bøgh C, Clarke SE, Pinder M, Sanyang F, Lindsay SW: Effect of passive zooprophylaxis on malaria transmission in the Gambia. J Med Entomol. 2001, 38 (6): 822-828. 10.1603/0022-2585-38.6.822.

    Article  PubMed  Google Scholar 

  48. 48.

    Bogh C, Clarke SE, Walraven GE, Lindsay SW: Zooprophylaxis, artefact or reality? A paired-cohort study of the effect of passive zooprophylaxis on malaria in The Gambia. Trans R Soc Trop Med Hyg. 2002, 96 (6): 593-596. 10.1016/S0035-9203(02)90320-2.

    Article  PubMed  Google Scholar 

  49. 49.

    Dekker T, Takken W: Differential responses of mosquito sibling species Anopheles arabiensis and An. quadriannulatus to carbon dioxide, a man or a calf. Med Vet Entomol. 1998, 12 (2): 136-140. 10.1046/j.1365-2915.1998.00073.x.

    CAS  Article  PubMed  Google Scholar 

  50. 50.

    Port G, Boreham P: The relationship of host size to feeding by mosquitoes of the Anopheles gambiae Giles complex (Diptera: Culicidae). Bull Entomol Res. 1980, 70: 133-144. 10.1017/S0007485300009834.

    Article  Google Scholar 

  51. 51.

    Gilles MT, Coetzee M: A Supplement to the Anophelinae of Africa South of the Sahara (Afrotoropical Region). 1987, Johannesburg, South Africa: Publications of the South African Institute for Medical Research

    Google Scholar 

  52. 52.

    Wirtz RA, Zavala F, Charoenvit Y, Campbell GH, Burkot TR, Schneider I, Esser KM, Beaudoin RL, Andre RG: Comparative testing of monoclonal antibodies against Plasmodium falciparum sporozoites for ELISA development. Bull World Health Organ. 1987, 65 (1): 39-45.

    PubMed Central  CAS  PubMed  Google Scholar 

  53. 53.

    Collins FH, Mendez MA, Rasmussen MO, Mehaffey PC, Besansky NJ, Finnerty V: A ribosomal RNA gene probe differentiates member species of the Anopheles gambiae complex. Am J Trop Med Hyg. 1987, 37 (1): 37-41.

    CAS  PubMed  Google Scholar 

  54. 54.

    Scott JA, Brogdon WG, Collins FH: Identification of single specimens of the Anopheles gambiae complex by the polymerase chain reaction. Am J Trop Med Hyg. 1993, 49 (4): 520-529.

    CAS  PubMed  Google Scholar 

  55. 55.

    Koekemoer LL, Kamau L, Hunt RH, Coetzee M: A cocktail polymerase chain reaction assay to identify members of the Anopheles funestus (Diptera: Culicidae) group. Am J Trop Med Hyg. 2002, 66 (6): 804-811.

    CAS  PubMed  Google Scholar 

  56. 56.

    Beier JC, Perkins PV, Wirtz RA, Koros J, Diggs D, Gargan TP, Koech DK: Bloodmeal identification by direct enzyme-linked immunosorbent assay (ELISA), tested on Anopheles (Diptera: Culicidae) in Kenya. J Med Entomol. 1988, 25 (1): 9-16.

    CAS  Article  PubMed  Google Scholar 

  57. 57.

    Kent RJ, Norris DE: Identification of mammalian blood meals in mosquitoes by a multiplexed polymerase chain reaction targeting cytochrome B. Am J Trop Med Hyg. 2005, 73 (2): 336-342.

    PubMed Central  CAS  PubMed  Google Scholar 

  58. 58.

    Kent RJ: Molecular methods for arthropod bloodmeal identification and applications to ecological and vector-borne disease studies. Mol Ecol Resour. 2009, 9 (1): 4-18. 10.1111/j.1755-0998.2008.02469.x.

    CAS  Article  PubMed  Google Scholar 

  59. 59.

    Faraway JJ: Extending the Linear Model with R: Generalized Linear, Mixed Effects and Nonparametric Regression Models. 2006, Boca Raton: CRC Press

    Google Scholar 

  60. 60.

    Development Core Team: A Language and Environment for Statistical Computing. 2011, Vienna, Austria: R Foundation for Statistical Computing,,

    Google Scholar 

  61. 61.

    Hewitt S, Kamal M, Muhammad N, Rowland M: An entomological investigation of the likely impact of cattle ownership on malaria in an Afghan refugee camp in the North West Frontier Province of Pakistan. Med Vet Entomol. 1994, 8 (2): 160-164. 10.1111/j.1365-2915.1994.tb00156.x.

    CAS  Article  PubMed  Google Scholar 

  62. 62.

    Midega JT, Mbogo CM, Mwnambi H, Wilson MD, Ojwang G, Mwangangi JM, Nzovu JG, Githure JI, Yan G, Beier JC: Estimating dispersal and survival of Anopheles gambiae and Anopheles funestus along the Kenyan coast by using mark-release-recapture methods. J Med Entomol. 2007, 44 (6): 923-929. 10.1603/0022-2585(2007)44[923:EDASOA]2.0.CO;2.

    PubMed Central  Article  PubMed  Google Scholar 

  63. 63.

    Chaves LF: An entomologist guide to demystify pseudoreplication: data analysis of field studies with design constraints. J Med Entomol. 2010, 47 (3): 291-298. 10.1603/ME09250.

    Article  PubMed  Google Scholar 

  64. 64.

    Davis JR, Hall T, Chee EM, Majala A, Minjas J, Shiff CJ: Comparison of sampling anopheline mosquitoes by light-trap and human-bait collections indoors at Bagamoyo, Tanzania. Med Vet Entomol. 1995, 9 (3): 249-255. 10.1111/j.1365-2915.1995.tb00130.x.

    CAS  Article  PubMed  Google Scholar 

  65. 65.

    Randolph SE, Dobson ADM: Pangloss revisited: a critique of the dilution effect and the biodiversity-buffers-disease paradigm. Parasitology. 2012, 139 (7): 847-863. 10.1017/S0031182012000200.

    CAS  Article  PubMed  Google Scholar 

  66. 66.

    Ostfeld RS: A Candide response to Panglossian accusations by Randolph and Dobson: biodiversity buffers disease. Parasitology. 2013, 140: 1196-1198. 10.1017/S0031182013000541.

    Article  PubMed  Google Scholar 

  67. 67.

    Killeen GF, McKenzie FE, Foy BD, Bøgh 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 (5): 469-476. 10.1016/S0035-9203(01)90005-7.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  68. 68.

    Killeen GF, Smith TA: Exploring the contributions of bed nets, cattle, insecticides and excitorepellency to malaria control: a deterministic model of mosquito host-seeking behaviour and mortality. Trans R Soc Trop Med Hyg. 2007, 101 (9): 867-880. 10.1016/j.trstmh.2007.04.022.

    PubMed Central  Article  PubMed  Google Scholar 

  69. 69.

    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 (1): 54-61.

    PubMed  Google Scholar 

  70. 70.

    Mutuku FM, Alaii JA, Bayoh MN, Gimnig JE, Vulule JM, Walker ED, Kabiru E, Hawley WA: Distribution, Description, and local knowledge of larval habitats of Anopheles gambiae s.l. in a village in western Kenya. Am J Trop Med Hyg. 2006, 74 (1): 44-53.

    PubMed  Google Scholar 

  71. 71.

    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 (7): e229-10.1371/journal.pmed.0040229.

    PubMed Central  Article  PubMed  Google Scholar 

  72. 72.

    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 (4 Suppl): 16-22.

    PubMed  Google Scholar 

  73. 73.

    Bøgh C, Pedersen EM, Mukoko DA, Ouma JH: Permethrin-impregnated bednet effects on resting and feeding behaviour of lymphatic filariasis vector mosquitoes in Kenya. Med Vet Entomol. 1998, 12 (1): 52-59. 10.1046/j.1365-2915.1998.00091.x.

    Article  PubMed  Google Scholar 

  74. 74.

    Gimnig JE, Ombok M, Kamau L, Hawley WA: Characteristics of larval anopheline (Diptera: Culicidae) habitats in Western Kenya. J Med Entomol. 2001, 38 (2): 282-288. 10.1603/0022-2585-38.2.282.

    CAS  Article  PubMed  Google Scholar 

  75. 75.

    Chaves LF, Keogh CL, Vazquez-Prokopec GM, Kitron UD: Combined sewage overflow enhances oviposition of Culex quinquefasciatus (Diptera: Culicidae) in Urban Areas. J Med Entomol. 2009, 46 (2): 220-226. 10.1603/033.046.0206.

    Article  PubMed  Google Scholar 

  76. 76.

    Nguyen AT, Williams-Newkirk AJ, Kitron UD, Chaves LF: Seasonal weather, nutrient dynamics and conspecific presence impacts on the southern house mosquito oviposition dynamics in combined sewage overflows. J Med Entomol. 2012, 49 (6): 1328-1338. 10.1603/ME12090.

    Article  PubMed  Google Scholar 

  77. 77.

    Edman JD, Scott TW, Costero A, Morrison AC, Harrington LC, Clark GG: Aedes aegypti (Diptera : Culicidae) movement influenced by availability of oviposition sites. J Med Entomol. 1998, 35 (4): 578-583.

    CAS  Article  PubMed  Google Scholar 

  78. 78.

    Bruce-Chwatt LJ, Garrett-Jones C, Weitz B: Ten years’ study (1955–64) of host selection by anopheline mosquitos. Bull World Health Organ. 1966, 35 (3): 405-439.

    PubMed Central  CAS  PubMed  Google Scholar 

  79. 79.

    Paaijmans KP, Thomas MB: The influence of mosquito resting behaviour and associated microclimate for malaria risk. Malar J. 2011, 10: 183-10.1186/1475-2875-10-183.

    PubMed Central  Article  PubMed  Google Scholar 

  80. 80.

    Lwetoijera DW, Kiware SS, Mageni ZD, Dongus S, Harris C, Devine GJ, Majambere S: A need for better housing to further reduce indoor malaria transmission in areas with high bed net coverage. Parasit Vectors. 2013, 6: 57-10.1186/1756-3305-6-57.

    PubMed Central  Article  PubMed  Google Scholar 

  81. 81.

    Fritz ML, Siegert PY, Walker ED, Bayoh MN, Vulule JR, Miller JR: Toxicity of bloodmeals from ivermectin-treated cattle to Anopheles gambiae s.l. Ann Trop Med Parasitol. 2009, 103 (6): 539-547. 10.1179/000349809X12459740922138.

    CAS  Article  PubMed  Google Scholar 

  82. 82.

    Fritz ML, Walker ED, Miller JR: Lethal and sublethal effects of Avermectin/Milbemycin parasiticides on the African Malaria vector, Anopheles arabiensis. J Med Entomol. 2012, 49 (2): 326-331. 10.1603/ME11098.

    CAS  Article  PubMed  Google Scholar 

  83. 83.

    Bett B, Randolph TF, Irungu P, Nyamwaro SO, Kitala P, Gathuma J, Grace D, Vale G, Hargrove J, McDermott J: Field trial of a synthetic tsetse-repellent technology developed for the control of bovine trypanosomosis in Kenya. Prev Vet Med. 2010, 97 (3–4): 220-227.

    CAS  Article  PubMed  Google Scholar 

  84. 84.

    Molina R, Espinosa-Góngora C, Gálvez R, Montoya A, Descalzo MA, Jiménez MI, Dado D, Miró G: Efficacy of 65% permethrin applied to dogs as a spot-on against Phlebotomus perniciosus. Vet Parasitol. 2012, 187 (3–4): 529-533.

    CAS  Article  PubMed  Google Scholar 

  85. 85.

    Pound JM, Lohmeyer KH, Davey RB, Miller JA, George JE: Efficacy of amitraz-impregnated collars on white-tailed deer (Artiodactyla: Cervidae) in reducing free-living populations of lone star ticks (Acari: Ixodidae). J Econ Entomol. 2012, 105 (6): 2207-2212. 10.1603/EC12219.

    CAS  Article  PubMed  Google Scholar 

Download references


We are deeply grateful to Dr. Yoshiki Aoki, Dr. Solomon Mpoke, Dr. Sam M Kariuki, Dr. Masaaki Shimada, Dr. Yoshio Ichinose, Dr. Satoshi Kaneko, Dr. Masahiro Horio, Dr. Hitoshi Kawada, Dr. Yoshinori Mitsui, Dr. Shingo Inoue, Dr. Shusuke Nakazawa, Dr. Yukiko Higa, Dr. Kentaro Kato, Dr. Takuya Kubo, Mr. Haruki Kazama, Dr. Yoshihide Maekawa and Mr. Tomonori Hoshi for critical comments and technical support, and Ms. Junko Sakemoto and Ms. Yukie Saito for providing administrative support. Special thanks go to Mr. George Odhiambo Orume, Ms. Beatrice Awuor Agola, Mr. James opodi Okwatch, Ms. Lorine Awuor Oyoo, Ms. Jecinta Odeo Lumumba, Mr. James Omondi Kongere, Ms. Mercy Syombua Mwania, Mr. Peter Lutiali Akweyia and Ms. Scholastica Achieng Wagalla devoted themselves to the fieldwork and experiments. Finally, we wish to express our gratitude to the residents in Koguna, Nyachabe, and Tabla who participated in this study. This study was supported in part by the Global COE Program, Nagasaki University, Japan. GOD and LFC were funded by Nagasaki University (Program for Nurturing Global Leaders in Tropical and Emerging Communicable Diseases).

Author information



Corresponding author

Correspondence to Noboru Minakawa.

Additional information

Competing interests

All authors declare that they have no competing interests.

Authors’ contributions

HI and NM conceived and designed this study. SN helped designing and planning the study in Kenya. GD, GS, and HI collected the field data, and HI and KF organized and conducted the laboratory work. HI, TS, LFC and NM performed the data analyses. HI drafted the first manuscript, and LFC and NM finalized the manuscript. All authors have read and approved the final manuscript.

Electronic supplementary material

Additional file 1: Supplement S1: Comparison of direct observation vs interview to estimate insecticide treated net (ITN) use. (DOCX 19 KB)

Additional file 2: Table S1: Female Anopheles spp. mosquito abundance by month, parity (gravid), feeding status and PCR based bloodmeal source. (DOCX 21 KB)

Additional file 3: Table S2: ELISA based bloodmeal sources in fed Anopheles spp. mosquitoes. (DOCX 21 KB)

Additional file 4: Table S3: Poisson GLMM selection for the abundance of all Anopheles species combined. (DOCX 15 KB)

Additional file 5: Table S4: Poisson GLMM selection for the abundance of each Anopheles species. (DOCX 17 KB)

Additional file 6: Table S5: Binomial GLMM selection for Anopheles spp. mosquito feeding on humans over other bloodmeal sources. (DOCX 15 KB)

Additional file 7: TableS 6: Binomial GLMM selection for each Anopheles species feeding on humans over other bloodmeal sources. (DOCX 17 KB)

Additional file 8: Table S7: Poisson rate GLMM selection for the sporozoite rate of all Anopheles spp. combined. (DOCX 17 KB)

Additional file 9: Table S8: Poisson rate GLMM selection for the sporozoite rate of Anopheles funestus s.s..(DOCX 15 KB)

Authors’ original submitted files for images

Rights and permissions

Reprints and Permissions

About this article

Cite this article

Iwashita, H., Dida, G.O., Sonye, G.O. et al. Push by a net, pull by a cow: can zooprophylaxis enhance the impact of insecticide treated bed nets on malaria control?. Parasites Vectors 7, 52 (2014).

Download citation


  • Anopheles
  • Bed net
  • Bloodfeeding
  • Plasmodium falciparum
  • Zooprophylaxis