- Open Access
Assessment of mosquito larval productivity among different land use types for targeted malaria vector control in the western Kenya highlands
Parasites & Vectorsvolume 8, Article number: 356 (2015)
Mosquito larval source management (LSM) is likely to be more effective when adequate information such as dominant species, seasonal abundance, type of productive habitat, and land use type are available for targeted sites. LSM has been an effective strategy for reducing malaria morbidity in both urban and rural areas in Africa where sufficient proportions of larval habitats can be targeted. In this study, we conducted longitudinal larval source surveillance in the western Kenya highlands, generating data which can be used to establish cost-effective targeted intervention tools.
One hundred and twenty-four (124) positive larval habitats were monitored weekly and sampled for mosquito larvae over the 85-week period from 28 July 2009 to 3 March 2011. Two villages in the western Kenya highlands, Mbale and Iguhu, were included in the study.
After preliminary sampling, habitats were classified into four types: hoof prints (n = 21; 17 % of total), swamps (n = 32; 26 %), abandoned goldmines (n = 35; 28 %) and drainage ditches (n = 36; 29 %). Positive habitats occurred in two land use types: farmland (66) and pasture (58). No positive larval habitats occurred in shrub land or forest.
A total of 46,846 larvae were sampled, of which 44.1 % (20,907) were from abandoned goldmines, 30.9 % (14,469) from drainage ditches, 22.4 % (10,499) from swamps and 2.1 % (971) from hoof prints. In terms of land use types, 57.2 % (26,799) of the sampled larvae were from pasture and 42.8 % (20,047) were from farmland. Of the specimens identified morphologically, 24,583 (52.5 %) were Anopheles gambiae s.l., 11,901 (25.4 %) were Culex quinquefasciatus, 5628 (12 %) were An. funestus s.l. and 4734 (10.1 %) were other anopheline species (An. coustani, An. squamosus, An. ziemanni or An. implexus). Malaria vector dynamics varied seasonally, with An.gambiae s.s. dominating during wet season and An.arabiensis during dry season. An increased proportion of An. arabiensis was observed compared to previous studies.
These results suggest that long-term monitoring of larval habitats can establish effective surveillance systems and tools. Additionally, the results suggest that larval control is most effective in the dry season due to habitat restriction, with abandoned goldmines, drainage ditches and swamps being the best habitats to target. Both farmland and pasture should be targeted for effective larval control. An increased proportion of An. arabiensis in the An. gambiae complex was noticed in this study for the very first time in the western Kenya highlands; hence, further control tools should be in place for effective control of An. arabiensis.
Anopheles gambiae sibling species have become dominant across Africa, even in areas where these species previously did not exist [1, 2]. New species recently have been identified in the An. gambiae complex . Currently, An. gambiae s.l. comprises eight species, which have been identified from different geographical locations in sub-Saharan Africa [3, 4]. Previous studies have found An. gambiae s.l. increasing in proportion or occurring in areas where it did not previously exist [1, 2].
Land use changes and topography have provided increased exposure to sunlight, which has contributed to the increased availability of potential breeding habitats in the African highlands [5–10]. Reclamation of swamps for agriculture has resulted in more sunlight reaching the swamps, which have then become potential larval habitats [7–9, 11]. Deforestation for agriculture and timber has led to increased availability of new productive An. gambiae s.l. breeding sites . Farm and pasture habitats are open and exposed to sunlight, thereby attracting gravid female mosquitoes for oviposition [1, 8, 9, 11, 13, 14]. Increased human population across sub-Saharan Africa has led to increased demand for land, resulting in accelerated land use changes [7, 12]. These land use changes result in microclimatic changes in breeding habitats due to increased temperatures [12, 15–17].
The rise in temperatures within cooler regions of Africa has expanded the distribution of malaria vectors and caused epidemics in highland areas where people are non-immune and have not been exposed to parasite infection [18–20]. In these areas, vector control efforts should be emphasized and implemented in a targeted manner. In the recent past, these areas have received wide coverage of insecticide-treated bed nets (ITNs) together with education on bed net usage . These highland sites have also received wide coverage of indoor residual spray (IRS) . While these strategies may be effective in adult vector control, more emphasis should be placed on controlling the aquatic stages, particularly since the protection offered by current tools in adult vector control has been compromised by intensified insecticide resistance within malaria vector populations [23, 24]. Targeting immature stages and developing monitoring information systems is of paramount importance in effective larval control. In different parts of Africa where larval source management has been practiced, reductions in adult mosquito incidence have been observed [25–29].
The current study monitored larval population dynamics in different land use settings and between different seasons (dry and rainy) for a period of 85 weeks to develop a larval surveillance system in the western Kenya highlands. Additionally, this study assessed the best time of year to conduct efficient larval control in the region.
Study area description
This study was conducted in the Mbale and Iguhu constituencies in the western Kenya highlands. Mbale (00. 04′57″N, 340. 43′16″E, 1,620m) and Iguhu (00. 09′41″N, 340. 44′36″E) lie at altitudes of 1,620m and 1,450m above sea level, respectively (Fig. 1). The dry seasons are January to February and temperature ranges between 25 and 30 °C. The main economic activity in these areas is small-scale food and cash crop farming, while a small proportion of inhabitants practice small-scale grazing.
Definitions of land use and habitat types
Land use types were categorized based on natural vegetation and the activities taking place on the land. The land with the highest tree coverage was considered forest; cultivated land was considered farmland; grazing areas were considered pastureland; and areas of bush land were considered shrub land. All positive habitats (having larvae of either anopheline or culicine species) were categorized as either drainage ditches (canals used to drain water from farms during rainy season), abandoned goldmines (pits left uncovered after mining activities have ceased), swamps (shallow bodies of water with vegetation) or hoof prints (cattle prints on wet ground), according to categorization from previous studies [1, 30]. The majority of these habitats were found in the valley bottom, due to the hilly topography of western Kenya. A total of one hundred and twenty-four (124) habitats were categorized. Positive larval habitats were found in farmland and pasture, while no positive habitats occurred in shrub land or forest. The distribution of habitat types was as follows: drainage ditches (36), hoof prints (21), swamps (32) and abandoned goldmines (35). Hoof prints were the only semi-permanent (short term) larval habitats; all other habitats were permanent. Larval surveys were conducted on all positive habitat types in all land use types.
Larval surveys, sampling and identification
Preliminary surveys were conducted in all habitats in all land use types using a standard dipper. All positive habitats were considered for follow-up in the main study. Each habitat was geo-referenced using a portable geographical positioning system (GPS). The 124 positive habitats which were selected were sampled for mosquito larvae once a week from 28 July 2009 through 3 March 2011. In habitats with sufficient water volume, a total of twenty (20) dips were made using a standard dipper (350mL, BioQuip Products, Inc. California, USA); for smaller habitats (mostly hoof prints) fewer dips were made. Larval abundance was calculated as the number of larvae per number of dips made in each habitat. Larval surveys were conducted between 10:00 and 13:00 h. All larvae (stages 1 to 4) sampled from each habitat were identified immediately in the field using morphological keys developed by Gillies and Coetzee . Stage 1 and 2 larval instars were returned to their respective habitats while a small number of specimens of stage 3 and 4 larval instars were taken for molecular identification in the laboratory. The larval specimens were preserved in absolute alcohol (70 %) and kept in a freezer at –20 °C until needed for molecular identification. Some larval specimens belonging to An.gambiae s.l. and the An. funestus group were taken for molecular identification of sibling species by polymerase chain reaction (PCR), following protocols developed by Scott et al. for An.gambiae s.l.  and by Koekemoer et al. for An. funestus .
Species abundance over time for each habitat type was analysed using generalized linear mixed models (GLMM) with multiple samples collected from each site over time. Larval abundance by land use and habitat type was analysed with GLMM using site as a random effect. Larval habitat and land use types were considered as fixed. Models for each species were developed separately. The abundance of identified An.gambiae s.s. and An. arabiensis mosquitoes between seasons was compared using one-way analysis of variance (ANOVA) with a significance level of 5 %. Data was analysed using PASW Statistics version 18.0 for Windows (SPSS Inc., Chicago, IL).
Ethical approval for this study was granted by the National Ethical Review Committee at the Kenya Medical Research Institute under the main project “Ecology of African highland malaria (II), SSC No. 1382”. Prior to implementation of the study, village leaders and elders were called to a meeting where the essence of the study was explained. Written consent to visit the habitats was obtained from all landowners in the selected sites.
Larval abundance and species identification
Among the specimens identified morphologically, 24,583 (52.5 %) were An. gambiae s.l., 5628 (12 %) were An. funestus s.l., 11,901 (25.4 %) were Cx. quinquefasciatus and 4734 (10.1 %) were other anopheline species (An. implexus, An. squamosus, An. ziemanni, or An. coustani). Among the 2350 specimens of An.gambiae s.l. identified using molecular techniques, 1445 (61.5 %) were An.gambiae s.s., 898 (38.2 %) were An. arabiensis and 7 (0.3 %) had no PCR product amplification. Of the 540 specimens in the An.funestus group, 172 (31.9 %) were An. funestus s.s., 161 (29.8 %) were An. leesoni, 68 (12.6 %) were An. rivulorum and 137 (25.7 %) were An. vaneedeni. Overall larval dynamics of the four main mosquito species (An.gambiae s.l., An. funestus, Cx. quinquefasciatus and other anopheline) showed An.gambiae s.l. to dominate over the other species, while An. funestus was the least abundant. The dynamics of the identified An.gambiae s.s. and An. arabiensis displayed similar trends to An.gambiae s.l.
In repeated measure analysis using sampling site as a random effect for 85 weeks, larval abundance was found to differ significantly between habitat types (DF = 3, F = 10.117, P < 0.0001).
Habitat and land use types in relation to larval abundance
In the current study, a total of 46,846 larvae were sampled, of which 30.9 % (14,469) were from drainage ditches, 44.6 % (20,907) from abandoned goldmines, 2.1 % (971) from hoof prints and 22.4 % (10,499) from swamps. Eighty-one percent (37,945) were anopheline species and 19 % (8901) were Cx. quinquefasciatus. Larval abundance varied by species and habitat type. An.gambiae s.l. occurrence showed no significant difference between habitats (DF = 3, F = 0.885, P = 0.448). Other anopheline species differed significantly in occurrence between habitats, with the highest occurrence in swamps (DF = 3, F = 14.460, P < 0.001). Cx.quinquefasciatus occurrence also differed significantly among habitats, again with the highest occurrence in swamps (DF = 3, F = 5.214, P = 0.001). An. funestus occurrence differed significantly between habitat types, with the highest abundance occurring in drainage ditches (DF = 3, F = 9.765, P < 0.001). No An. funestus larvae were recorded from hoof prints.
When analysed by land use type, 42.8 % (20,047) of the sampled specimens were from farmland and 57.2 % (26,799) from pasture. An. gambiae s.l. was statistically more abundant in pasture than in farmland (DF = 1, F = 4.824, P = 0.028), as was An. funestus (DF = 1, F = 5.133, P = 0.024). Other anopheline species had no significant difference in occurrence between pasture and farmland (DF = 1, F = 0.001, P = 0.979); nor did Cx. quinquefasciatus (DF = 1, F = 4.824, P = 0.143).
Anopheline larval abundance by seasonality
Among anopheline species sampled, there were significant differences in larval abundance by season, with An.gambiae s.l larvae more abundant in dry season than in rainy season (F = 17.76, df = 1, P ≥ 0.001) and An.funestus larvae more abundant in rainy season than in dry season (F = 4.16, df =1, P = 0.045). Other anopheline larvae were also significantly more abundant in rainy season than in dry season (F = 5.64, DF = 1, P = 0.020) (Fig. 2). Larval density varied significantly among habitats between weeks (DF = 3, F = 27.18, P < 0.0001).
This study has shown that extending larval habitat follow-up for longer periods of time can provide critical information for larval source management. Such information includes habitat type, larval abundance, species composition and dynamics, seasonality, and land use type. It is interesting to note that, compared to historical findings, this study found a higher proportion of An. arabiensis abundance (38.2 %) in the study area. In a previous study by Wamae et al. at the same sites, the proportion of An.arabiensis was reported to be 2.6 % among all An.gambiae s.l. sampled . It is possible that changes in land cover and topography in the western Kenya highlands have increased the suitability of breeding sites for An.arabiensis, as was found in previous studies [7, 12, 35]. A higher proportion of An.arabiensis (38.2 %) has been found in larval habitat compared to that found in indoor sampling using CDC light traps and pyrethrum spray catches (2.6 %) . Our current results show that we can assess mosquito species composition with greater accuracy by establishing larval habitat monitoring surveillance systems. This is important, as sampling methods such as CDC light traps, pyrethrum spray catches and window/eave traps have been compromised by the wide coverage of ITN and IRS programmes .
In this study, Anopheles and Culex species were found in all four habitat types (hoof prints, abandoned goldmines, swamps and drainage ditches) with the exception of An. funestus, which was not found in hoof prints. Naturally, An. funestus breeds in large water bodies which are permanent with shade and vegetation [38, 39]. Other studies conducted at the same sites have shown similar species composition across habitat types, with An. funestus also not recorded from hoof prints [6–8, 11, 13, 40]. In designing an effective control tool for malaria vectors, the major target habitats should be abandoned goldmines, swamps and drainage ditches which are open to sunlight. These habitats attract more An. gambiae s.l. mosquitoes to oviposit, leading to a higher abundance of larvae than in other habitat types. Also, An.gambiae s.l. was observed in all habitats except swamps during the dry season, a result of restricted water sources. Cx.quinquefasciatus colonized swamps, drainage ditches and abandoned goldmines during dry seasons when these habitats had high rates of organic matter and plant decomposition; by contrast, such characteristics have been found to cause An.gambiae s.l. to avoid breeding in matured habitats . However, Cx. quinquefasciatus has been found to be only a nuisance vector in the study area, not a vector for disease parasites. An. gambiae s.l. has been found to colonize fresh shallow and temporary habitats with relatively low grass cover and high sunlight exposure [1, 11]. This information about species composition and habitat preference is of paramount importance in developing cost-effective and efficient larval control. In Brazil, Killeen and others showed that a well-designed larval control programme can contribute to the elimination of An.gambiae . Long-term monitoring of larval habitats has been found to have a positive impact on our understanding of habitat productivity and larval abundance [1, 11].
Malaria vector and non-vector species were found in both of the land use types that had positive breeding sites (i.e., farmland and pasture). Comparison of larval abundance between the two land use types shows a higher abundance in pasture than in farmland. Previous studies have found higher larval abundance in farmland than in pasture [9, 11]. The differences in larval abundance between the two land use types might have contributed to the increased proportion of An.arabiensis, which can take blood meals from grazing cattle, rest in bushes and, soon after digestion, oviposit in nearby habitats. This maybe a survival strategy and adaptation behaviour for An.arabiensis in this highland site in western Kenya. The coverage of pyrethroids used in IRS and ITNs in western Kenya may also shift species composition in favour of An.arabiensis, since the availability of an alternative blood meal allows An. arabiensis mosquitoes to avoid exposure to the current intervention tools. In other parts of Kenya, similar scenarios have been found. After massive coverage with intervention tools there was a tremendous shift in malaria vector species composition from An.gambiae s.s. to An.arabiensis . The same trend was observed in Tanzania in areas of intensive bed net coverage . Since An.gambiae s.s. mosquitoes feed indoors (mostly on humans) and rest indoors where intervention coverage exists, their blood meal seeking time is increased and survivorship is reduced ; hence, the population of zoophilic and exophilic An.arabiensis overlaps and dominates. This suggests the possibility of outdoor malaria transmission in this area of western Kenya where the population of An.arabiensis is escalating. The only option for reducing outdoor vector populations is to target larval source management, which is possible when efficient and effective surveillance systems are established. It is plausible that outdoor malaria vector populations have increased as a result of induced exophily due to IRS and ITN coverage [46, 47].
In larval habitat monitoring and control, understanding seasonal abundance and species composition is of great importance when designing efficient, cost-effective tools. Seasonality has always been associated with the abundance of immature stages and adult mosquitoes . The abundance of both larvae and adult malaria vectors has been observed to increase soon after the rainy season with the stabilization of habitats . In this study, seasonality has been demonstrated to influence larval abundance for all Anopheles species found in the study site, with the dry season having statistically significantly higher larval abundance than the rainy season. This can be attributed to the flush effect during the rainy season, which may result in unexpected larval loss in habitats [1, 11, 49]. However, the population can proliferate over time, reaching its peak during the beginning of dry season when weather is favourable and the shortened gonotrophic cycle leads to high larval density. This in turn is likely to increase the population of adults, which are then able to lay higher densities of eggs when habitats are completely well established. The topography of western Kenya is undulating, and larval habitats occur mostly in the valley bottoms, which are easily flooded during the rainy season [1, 49]. The dry season has been shown to have more larval abundance for all species. This can be attributed to stable habitats, which favour oviposition by gravid females and provide enough food for larvae survivorship [14, 50]. During dry season, habitats are few and exposed to sunlight, hence facilitating microbial decomposition activity and photosynthetic algal growth, which in turn provides food to mosquito larvae [1, 11, 30, 51]. Additionally, during dry season habitats maintain stable temperature, which has been found to shorten larval immature stages . In this study larval habitats in dry season were few, but each site contained high densities of larvae of all species found in that site. Thus this long-term study has shown that dry season is the best time to target for effective larval habitat control. Habitats are restricted and few, so follow-up management will be effective with regard to time, space, cost and expected outputs. In previous studies, the best time for larval control in this study area was suggested to be between the short rainy season and the dry season [1, 11]. The current study suggests the same but with more detailed habitats, land use patterns and larval abundance dynamics. If control programmes are coordinated in these areas during dry season, larvicidal efficacy will be higher, as there is no larvicidal dilution effect and larval density is at its peak. Monitoring is likely to be easier, as habitats are few and restricted. From different studies conducted across Africa, it is evident that if larval source management is targeted properly, malaria control and elimination can be achieved . Community willingness to be involved in malaria control has also shown a positive impact and should continue to be cultivated and extended throughout sub-Saharan Africa [54–56].
These results suggest that long-term monitoring of larval habitats can establish effective surveillance systems and tools. Additionally, the results suggest that larval control is most effective in the dry season due to habitat restriction, with abandoned goldmines, drainage ditches and swamps being the best habitats to target. The land use type that should be targeted are both pasture and farmlands. An increased proportion of An. arabiensis in the An. gambiae complex was noticed in this study for the very first time in the western Kenya highlands; hence, further control tools should be in place for effective control of An. arabiensis.
Kweka E, Zhou G, Lee M-C, Gilbreath T, Mosha F, Munga S, et al. Evaluation of two methods of estimating larval habitat productivity in western Kenya highlands. Parasit Vectors. 2011;4:110.
Chen H, Githeko A, Zhou G, Githure J, Yan G. New records of Anopheles arabiensis breeding on the Mount Kenya highlands indicate indigenous malaria transmission. Malar J. 2006;5:17.
Coetzee M, Hunt RH, Wilkerson R, Della Torre A, Coulibaly MB, Besansky N. Anopheles coluzzii and Anopheles amharicus, new members of the Anopheles gambiae complex. Zootaxa. 2013;3619:246–74.
Coetzee M, Craig M, le Sueur D. Distribution of African malaria mosquitoes belonging to the anopheles gambiae complex. Parasitol Today. 2000;16:74–7.
Githeko A, Ayisi J, Odada P, Atieli F, Ndenga B, Githure J, et al. Topography and malaria transmission heterogeneity in western Kenya highlands: prospects for focal vector control. Malar J. 2006;5:107.
Atieli H, Zhou G, Lee M-C, Kweka E, Afrane Y, Mwanzo I, et al. Topography as a modifier of breeding habitats and concurrent vulnerability to malaria risk in the western Kenya highlands. Parasit Vectors. 2011;4:241.
Minakawa N, Munga S, Atieli F, Mushinzimana E, Zhou G, Githeko AK, et al. Spatial distribution of anopheline larval habitats in Western Kenyan highlands: effects of land cover types and topography. Am J Trop Med Hyg. 2005;73:157–65.
Munga S, Minakawa N, Zhou G, Mushinzimana E, Barrack OO, Githeko AK, et al. Association between land cover and habitat productivity of malaria vectors in Western Kenyan highlands. Am J Trop Med Hyg. 2006;74:69–75.
Munga S, Yakob L, Mushinzimana E, Zhou G, Ouna T, Minakawa N, et al. Land use and land cover changes and spatiotemporal dynamics of anopheline larval habitats during a four-year period in a highland community of Africa. Am J Trop Med Hyg. 2009;81:1079–84.
Sovi A, Govoetchan R, Tokponnon F, Hounkonnou H, Aikpon R, Agossa F, et al. Impact of land-use on malaria transmission in the Plateau region, Southeastern Benin. Parasit Vectors. 2013;6:352.
Kweka EJ, Zhou G, Munga S, Lee M-C, Atieli HE, Nyindo M, et al. Anopheline larval habitats seasonality and species distribution: a prerequisite for effective targeted larval habitats control programmes. PLoS One. 2012;7:e52084.
Tuno N, Okeka W, Minakawa N, Takagi M, Yan G. Survivorship of anopheles gambiae sensu stricto (Diptera: Culicidae) Larvae in Western Kenya highland forest. J Med Entomol. 2005;42:270–7.
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:29–35.
McCrae AW. Oviposition by African malaria vector mosquitoes. II. Effects of site tone, water type and conspecific immatures on target selection by freshwater Anopheles gambiae Giles, sensu lato. Ann Trop Med Parasitol. 1984;78:307–18.
Zhou G, Minakawa N, Githeko AK, Yan G. Climate variability and malaria epidemics in the highlands of East Africa. Trends Parasitol. 2005;21:54–6.
Himeidan Y, Zhou G, Yakob L, Afrane Y, Munga S, Atieli H, et al. Habitat stability and occurrences of malaria vector larvae in western Kenya highlands. Malar J. 2009;8:234.
Omumbo J, Lyon B, Waweru S, Connor S, Thomson M. Raised temperatures over the Kericho tea estates: revisiting the climate in the East African highlands malaria debate. Malar J. 2011;10:12.
Wanjala C, Waitumbi J, Zhou G, Githeko A. Identification of malaria transmission and epidemic hotspots in the western Kenya highlands: its application to malaria epidemic prediction. Parasit Vectors. 2011;4:81.
Badu K, Siangla J, Larbi J, Lawson B, Afrane Y, Ong'echa J, et al. Variation in exposure to Anopheles gambiae salivary gland peptide (gSG6-P1) across different malaria transmission settings in the western Kenya highlands. Malar J. 2012;11:318.
Badu K, Afrane Y, Larbi J, Stewart V, Waitumbi J, Angov E, et al. Marked variation in MSP-119 antibody responses to malaria in western Kenyan highlands. BMC Infect Dis. 2012;12:50.
Atieli H, Zhou G, Afrane Y, Lee M-C, Mwanzo I, Githeko A, et al. Insecticide-treated net (ITN) ownership, usage, and malaria transmission in the highlands of western Kenya. Parasit Vectors. 2011;4:113.
Zhou G, Afrane Y, Dixit A, Atieli H, Lee M-C, Wanjala C, et al. Modest additive effects of integrated vector control measures on malaria prevalence and transmission in western Kenya. Malar. 2013;12:256.
Asale A, Getachew Y, Hailesilassie W, Speybroeck N, Duchateau L, Yewhalaw D. Evaluation of the efficacy of DDT indoor residual spraying and long-lasting insecticidal nets against insecticide resistant populations of Anopheles arabiensis Patton (Diptera: Culicidae) from Ethiopia using experimental huts. Parasit Vectors. 2014;7:131.
Kitau J, Oxborough R, Kaye A, Chen-Hussey V, Isaacs E, Matowo J, et al. Laboratory and experimental hut evaluation of a long-lasting insecticide treated blanket for protection against mosquitoes. Parasit Vectors. 2014;7:129.
Fillinger U, Lindsay S. Larval source management for malaria control in Africa: myths and reality. Malar J. 2011;10:353.
Worrall E, Fillinger U. Large-scale use of mosquito larval source management for malaria control in Africa: a cost analysis. Malar J. 2011;10:338.
Smith DL, Perkins TA, Tusting LS, Scott TW, Lindsay SW. Mosquito population regulation and larval source management in heterogeneous environments. PLoS One. 2013;8:e71247.
Tusting LS, Thwing J, Sinclair D, Fillinger U, Gimnig J, Bonner KE, et al. Mosquito larval source management for controlling malaria. Cochrane Database Syst Rev. 2013;8:CD008923.
Geissbühler Y, Kannady K, Chaki PP, Emidi B, Govella NJ, Mayagaya V, et al. Microbial larvicide application by a large-scale, community-based program reduces malaria infection prevalence in urban Dar Es Salaam, Tanzania. PLoS ONE. 2009;4:e5107.
Mwangangi J, Kahindi S, Kibe L, Nzovu J, Luethy P, Githure J, et al. Wide-scale application of Bti/Bs biolarvicide in different aquatic habitat types in urban and peri-urban Malindi, Kenya. Parasitol Res. 2011;108:1355–63.
Gillies MT, Coetzee M. A Supplement to the Anophelinae of Africa, South of the Sahara, Publication of the South African Institute of International Medicine and Research No. 55. 1987.
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:520–9.
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:804–11.
Wamae P, Githeko A, Menya D, Takken W. Shading by napier grass reduces malaria vector larvae in natural habitats in Western Kenya highlands. EcoHealth. 2010;7:485–97.
Gouagna L, Rakotondranary M, Boyer S, Lemperiere G, Dehecq J-S, Fontenille D. Abiotic and biotic factors associated with the presence of Anopheles arabiensis immatures and their abundance in naturally occurring and man-made aquatic habitats. Parasit Vectors. 2012;5:96.
Ndenga B, Githeko A, Omukunda E, Munyekenye G, Atieli H, Wamai P, et al. Population dynamics of malaria vectors in Western Kenya highlands. J Med Entomol. 2006;43:200–6.
Sikaala C, Killeen G, Chanda J, Chinula D, Miller J, Russell T, et al. Evaluation of alternative mosquito sampling methods for malaria vectors in Lowland South - East Zambia. Parasit Vectors. 2013;6:91.
Sinka M, Bangs M, Manguin S, Coetzee M, Mbogo C, Hemingway J, et al. The dominant Anopheles vectors of human malaria in Africa, Europe and the Middle East: occurrence data, distribution maps and bionomic precis. Parasit Vectors. 2010;3:117.
Juliano SA. Species interactions among larval mosquitoes: context dependence across habitat gradients. Ann Rev Entomol. 2009;54:37–56.
Nmor J, Sunahara T, Goto K, Futami K, Sonye G, Akweywa P, et al. Topographic models for predicting malaria vector breeding habitats: potential tools for vector control managers. Parasit Vectors. 2013;6:14.
Munga S, Vulule J, Kweka E. Response of anopheles gambiae s.l. (Diptera: Culicidae) to larval habitat age in western Kenya highlands. Parasit Vectors. 2013;6:13.
Killeen GF, Fillinger U, Kiche I, Gouagna LC, Knols BGJ. Eradication of Anopheles gambiae from Brazil: lessons for malaria control in Africa? Lancet Infect Dis. 2002;2:618–27.
Bayoh MN, Mathias D, Odiere M, Mutuku F, Kamau L, Gimnig J, et al. Anopheles gambiae: historical population decline associated with regional distribution of insecticide-treated bed nets in western Nyanza Province, Kenya. Malar J. 2010;9:62.
Kitau J, Oxborough RM, Tungu PK, Matowo J, Malima RC, Magesa SM, et al. Species shifts in the Anopheles gambiae complex: do LLINs successfully control Anopheles arabiensis? PLoS One. 2012;7:e31481.
Mejia P, Teklehaimanot H, Tesfaye Y, Teklehaimanot A. Physical condition of Olyset(R) nets after five years of utilization in rural western Kenya. Malar J. 2013;12:158.
Antonio-Nkondjio C, Demanou M, Etang J, Bouchite B. Impact of cyfluthrin (Solfac EW050) impregnated bed nets on malaria transmission in the city of Mbandjock: lessons for the nationwide distribution of long-lasting insecticidal nets (LLINs) in Cameroon. Parasit Vectors. 2013;6:10.
Futami K, Dida G, Sonye G, Lutiali P, Mwania M, Wagalla S, et al. Impacts of insecticide treated bed nets on Anopheles gambiae s.l. populations in Mbita district and Suba district, Western Kenya. Parasit Vectors. 2014;7:63.
Imbahale S, Paaijmans K, Mukabana W, van Lammeren R, Githeko A, Takken W. A longitudinal study on Anopheles mosquito larval abundance in distinct geographical and environmental settings in western Kenya. Malar J. 2011;10:81.
Paaijmans KP, Wandago MO, Githeko AK, Takken W. Unexpected high losses of Anopheles gambiae larvae due to rainfall. PLoS One. 2007;2:e1146.
McCrae AW. Oviposition by African malaria vector mosquitoes. I. Temporal activity patterns of caged, wild-caught, freshwater Anopheles gambiae Giles sensu lato. Ann Trop Med Parasitol. 1983;77:615–25.
Gilbreath T, Kweka E, Afrane Y, Githeko A, Yan G. Evaluating larval mosquito resource partitioning in western Kenya using stable isotopes of carbon and nitrogen. Parasit Vectors. 2013;6:353.
Lyons C, Coetzee M, Chown S. Stable and fluctuating temperature effects on the development rate and survival of two malaria vectors: anopheles arabiensis and Anopheles funestus. Parasit Vectors. 2013;6:104.
Walker K, Lynch M. Contributions of Anopheles larval control to malaria suppression in tropical Africa: review of achievements and potential. Med Vet Entomol. 2007;21:2–21.
Chaki P, Dongus S, Fillinger U, Kelly A, Killeen G. Community-owned resource persons for malaria vector control: enabling factors and challenges in an operational programme in Dar es Salaam, United Republic of Tanzania. Hum Resour Health. 2011;9:21.
Mukabana W, Kannady K, Kiama GM, Ijumba J, Mathenge E, Kiche I, et al. Ecologists can enable communities to implement malaria vector control in Africa. Malar J. 2006;5:9.
Chaki P, Govella N, Shoo B, Hemed A, Tanner M, Fillinger U, et al. Achieving high coverage of larval-stage mosquito surveillance: challenges for a community-based mosquito control programme in urban Dar es Salaam, Tanzania. Malar J. 2009;8:311.
We thank Wilberforce Muyeso, Amos Wabwile, Geofrey Bullemi and Josephine Shikholwa (RIP) for technical assistance during field work. Dr. Ming-Chieh Lee is acknowledged for providing Fig. 1 in this manuscript. Dr. Guofa Zhou is acknowledged for repeated measure data analysis. This work was supported by grants from the National Institutes of Health (NIH) D43 TW01505 and R01 A150243. We thank the western Kenya community for supporting the study throughout.
The authors declare that they have no competing interests.
EJK, YEH and SM conceived and designed the study and analysed the data. EJK wrote the manuscript. EJK, YEH, SM, AKG and GY revised the manuscript. All authors reviewed the work and approved the final version for submission.