Evaluation of two methods of estimating larval habitat productivity in western Kenya highlands
Parasites & Vectorsvolume 4, Article number: 110 (2011)
Malaria vector intervention and control programs require reliable and accurate information about vector abundance and their seasonal distribution. The availability of reliable information on the spatial and temporal productivity of larval vector habitats can improve targeting of larval control interventions and our understanding of local malaria transmission and epidemics. The main objective of this study was to evaluate two methods of estimating larval habitat productivity in the western Kenyan highlands, the aerial sampler and the emergence trap.
The study was conducted during the dry and rainy seasons in 2008, 2009 and 2010. Aerial samplers and emergence traps were set up for sixty days in each season in three habitat types: drainage ditches, natural swamps, and abandoned goldmines. Aerial samplers and emergence traps were set up in eleven places in each habitat type. The success of each in estimating habitat productivity was assessed according to method, habitat type, and season. The effect of other factors including algae cover, grass cover, habitat depth and width, and habitat water volume on species productivity was analysed using stepwise logistic regression
Habitat productivity estimates obtained by the two sampling methods differed significantly for all species except for An. implexus. For for An. gambiae s.l. and An. funestus, aerial samplers performed better, 21.5 and 14.6 folds, than emergence trap respectively, while the emergence trap was shown to be more efficient for culicine species. Seasonality had a significant influence on the productivity of all species monitored. Dry season was most productive season. Overall, drainage ditches had significantly higher productivity in all seasons compared to other habitat types. Algae cover, debris, chlorophyll-a, and habitat depth and size had significant influence with respect to species.
These findings suggest that the aerial sampler is the better of the two methods for estimating the productivity of An. gambiae s.l. and An. funestus in the western Kenya highlands and possibly other malaria endemic parts of Africa. This method has proven to be a useful tool for monitoring malaria vector populations and for control program design, and provides useful means for determining the most suitable sites for targeted interventions.
Although recent studies [1–3], have demonstrated that malaria disease burden is on the decline in several sub-Saharan African countries, the disease still remains a major public health problem. This reduction in the disease burden has parsimoniously been attributed to recent scale-up of control tools in some parts of Africa . To design effective vector control programmes, national malaria control programmes in Africa require accurate information on vector densities and species composition. Further, a clear understanding of the epidemiology of the disease is important if an effective intervention programme is to be developed. Monitoring of habitat productivity and vector population dynamics would provide critical information for vector surveillance in areas where control interventions are implemented. In sub-Saharan African countries, Anopheles gambiae sensu lato (s.l.) and An. funestus groups contain the most efficient malaria transmitting vectors [4–6].
In the recent past, western Kenya highlands have experienced increased frequencies of malaria epidemics and are now reported to have high rates of transmission [7–11]. Significant increases in the human population and subsequent land use changes, such as deforestation and swamp cultivation, have been hypothesized to be among the several mechanisms responsible for this increase in malaria epidemics in the highland regions of western Kenya, [12, 13]. Swamp reclamation for agricultural development has resulted in the creation of potential breeding habitats for An. gambiae s.s and An. arabiensis[4, 14]. High rates of deforestation have led to a rise in local temperatures which has been shown to lead to increased mosquito larvae survivorship and development of parasites in adult mosquitoes [13, 15]. These factors have increased the productivity of malaria vectors , thus increasing the risk of malaria transmission at these highland sites .
The principal malaria vector An. gambiae s.l has been demonstrated to breed in sunlit, temporary water bodies such as hoof prints, goldmines, and drainage ditches while An. funestus is typically associated with permanent water bodies such as swamps and drainage ditches [4, 16]. The habitats of these vectors in these highlands are highly concentrated in valley bottoms in both the rainy and dry seasons [17, 18]. Although adult emergence from habitats is the primary determinant of vector density, the productivity of larval habitats for aquatic stages of the malaria vector Anopheles gambiae (the stage preceding adult metamorphosis) is not clearly understood . This is partly contributed to by methodological problems associated with larval ecology studies of Anopheles gambiae. For example, at the moment, there is no standardized method for accurate estimation of productivity of adult mosquitoes from larval habitats. However, this is essential for effective vector surveillance and planning larval control programmes . Thus, there is need to evaluate an efficient and sensitive method for estimating malaria vector productivity in different habitats so as to translate the information into meaningful epidemiological data for monitoring An. gambiae s.l. and An. funestus. Accurate information on habitat productivity can lead to the determination of the seasons during which habitat source reduction targeting is appropriate, thus leading to vector reduction [22, 23]. Source reduction have shown to be effective in vector reduction in different part of Africa [24–27]
In estimating habitat productivity, three main methods have been used for sampling different mosquito larvae: dipping, netting, and pipetting . These techniques have been used to estimate the larval density in different habitat types. Among them, pipetting is predominately used in estimating the productivity of tree holes in breeding mosquitoes . Use of netting and dipping methods however are limited in habitats with higher grass cover . Due to the inconsistency of these methods, there is need to evaluate emergence traps and aerial samplers in estimating habitat productivity in different seasons and habitats to determine if they would be more effective than methods used previously. Therefore, the objective of this study was to compare the efficiency of the emergence trap and the aerial sampler in estimating larval habitat productivity under different habitat types in highland areas of western Kenya.
Materials and methods
The present study was conducted in six sites at Mbale village, Vihiga district in the western Kenyan highlands (Figure 1). The area has different land use types and topography that have enhanced the availability of potential mosquito breeding habitats. The three main land use types at Mbale are farmland, pasture, and forests. No productive mosquito larval habitats were found in forested areas during the study period. Habitat types used were drainage ditches (these are canals constructed to channel excess water to nearby streams and rivers, from farmland), natural swamps, and abandoned goldmines. Drainage ditches were found in the farmland, while swamps and abandoned goldmines were found in pastures. The main cultivated crops are tea, maize, beans, finger millet, and sorghum.
Larval habitats selection and survey
The different types of larval habitats, i.e., swamps, drainage ditches, and abandoned goldmines, were identified at the beginning of the study in each season. For each survey, the different types of breeding sites were sampled using standard dippers (350 ml). Ten dips were made at each habitat and the presence or absence of mosquito larvae was recorded. In the first rainy season survey, (June - 2008 to August 2008), sixty habitats were selected and each habitat type had ten emergence traps and ten aerial samplers. In the first, (February-2009 to April 2009) and the second dry season (February-2010 to April, 2010), and the second rainy season (June-2009 to August, 2009), 66 habitats were sampled using eleven paired emergence traps and aerial samplers for each habitat type per season. Each trapping method was done during dry and wet season in two different years to avoid yearly weather change effects.
The aerial samplers used in this study had a diameter of 0.28 m and depth of 0.30 m, and were used eight times per habitat (surface area of 0.5 m2) in each sampling day. (Figure 2A). The emergence trap was made of an iron-framed cage, either 1 m × 0.5 m × 1 m (surface area of 0.5 m2) or 1 m × 1 m × 1 m (surface area of 1 m2), covered by fine mosquito polystyrene netting material (Figure 2B). The 0.5 m2 traps were used for small goldmines and drainage ditches while 1 m2 traps were used for large goldmines and swamps. Emergence traps were shifted daily in the same habitat to allow gravid mosquitoes to lay eggs. In all seasons, the samplings were done for 60 days consecutively for both emergence traps and aerial samplers.
Mosquito sampling and processing
Adult mosquitoes were sampled from emergence traps using a mechanical aspirator as described in the World Health Organization (WHO) entomological field manual . Pupae sampled in the aerial sampler were retained in a field insectary until the adults emerged. All adults were anaesthetized using chloroform, and identified using a morphological key in the laboratory as described in Gillies and Coetzee . Mosquito specimens were preserved in absolute ethanol (purity 96%) for molecular species identification using the polymerized chain reaction (PCR) protocols developed by Scott et al.,  for An. gambiae s.l. and Koekemoer et al., for An. funestus group. We utilized a morphological key developed by Darsie et al. for culicine species identification.
Data were recorded daily by date, sampling method, habitat type, and mosquito species. Data analysis was performed using PWAS Statistics 18 (SPSS Inc., Chicago, IL). The most efficient sampling method for estimating mosquito habitat productivity was chosen by comparing sampling methods in all seasons using paired sample t-tests. Assessment of the effects of seasonality, habitat type, and sampling method in estimating habitat productivity were analysed by multi-factorial analysis of variance (MANOVA) using species as a dependent variable. Pairwise comparison of mosquito productivity with method, habitat type, and season was done using Tukey's HSD test of MANOVA with repeated measure. The contribution of algal cover, grass cover, habitat depth and width, and habitat water surface area on species productivity was analysed using stepwise regression analysis (only culicines species, An. gambiae s.l. and An. funestus were included; other species were excluded due to their low densities). The efficiency of the sampling methods in estimating habitat productivity was determined using the three most abundant species (culicine, An. gambiae s.l. and An. funestus) for each season, using cumulative habitat productivity.
A total of 11,514 female mosquitoes were sampled during the study period. Out of these, 5,014 (43.5%) were An. gambiae s.l., 1,247 (10.8%) An. funestus and 5,253 (45.6%) were culicine species mosquito. Drainage ditches produced 4,565 (39.6%), abandoned goldmines 2,887 (25.1%), and swamps 4,062 (35.3%) mosquitoes. Among ten percent (501) of all An. gambiae s.l. identified; 303 (60.5%) were identified as An. gambiae s.s, 169 (33.7%) were An. arabiensis and 29 (5.8%) specimens of mosquitoes had no PCR products. Among identified An. gambiae s.s, 152 (50.2%) were sampled from drainage ditches while 151 (49.8%) were from abandoned goldmines. For An. arabiensis, 109 (64.7%) came from drainage ditches and 60 (35.3%) were from abandoned goldmines.
The 10% (525) of all An. funestus group identified, 514 (97.9%) were found to be An. funestus while the 11 (2.1%) mosquitoes specimens had no PCR products. All An. funestus were sampled, came from swamps. By using a morphological key, it was revealed that all culicines were Culex quinquefasciatus.
Overall estimation of habitat productivity by comparison of sampling method has shown the aerial sampler to be the more efficient method in estimating the productivity of anopheline species while emergence traps did better for culicine species (Table 1). Upon comparison, it was found that for An. gambiae s.l. and An. funestus, the aerial sampler performed better by 21.5 and 14.6 folds than the emergence trap respectively. However, the two sampling methods showed no difference in estimating the productivity of An. implexus (Table 1).
The pairwise analysis using Tukey's HSD test of MANOVA with repeated measures showed that, the productivity of An. gambiae s.l. was significantly higher in drainage ditches than in any other habitat types for all seasons. An. funestus productivity in swamps was significantly higher than in other habitat types throughout all seasons (Tukey HSD test P < 0.05) (Table 2). The culicine species productivity was significantly higher in goldmines when compared with habitat types and sampling methods in all the seasons (Tukey HSD test P < 0.05).
Results of stepwise regression analysis show that the productivity of An. gambiae s.l was significantly associated with high algae cover (P < 0.0001), low grass cover (P < 0.0001) and debris (P = 0.05), but chlorophyll a, habitat width and depth had no significant effect. Only grass coverage (P < 0.01) showed a significant effect on the production of An. funestus. In assessing the impact of environmental factors on productivity of culicine species, only chlorophyll a (P = 0.05) was found to have a significant effect.
The efficiency of the sampling methods was found to vary in all four seasons for all three mosquito species. The aerial sampler was shown to be a better tool than emergence traps for estimating habitat productivity in both dry and rainy seasons (Table 2, Figure 3).
The findings of this study indicate that aerial sampler is the most efficient in estimating habitat productivity of malaria vectors in different habitat types and seasons in the highlands of western Kenya. Mosquito species identified from An. gambiae s.l and An. funestus group have revealed the malaria vectors found in this study area are An. gambiae s.s, An. arabiensis and An. funestus. An. gambiae s.s and An. funestus have been found to be anthropophilic vectors while An. arabiensis is both anthropophilic and zoophilic depend on geographical location and host availability [28, 32]. Similar species composition information was reported from other sites in the western Kenyan highlands [6, 33], where mosquitoes were sampled by other methods. In this study mosquito species productivity was found to be habitat-type dependent as An. gambiae s.s and An. arabiensis were abundant in drainage ditches and abandoned goldmines, An. funestus in swamps, and Culicine species in abandoned goldmines. This pattern was similar to other studies conducted in western Kenya [13, 14, 33, 34].
Aerial samplers performed better than emergence traps in estimating the productivity in larval habitats, presumably due to the netting material covering the emergence trap which inhibited gravid mosquitoes from laying eggs in those habitats. Aerial samplers showed higher efficiency for all species' productivity estimations and did not disrupt mosquito oviposition in habitats, thereby producing more reliable estimates . Aerial samplers yielded more reliable information on mosquito larval habitat productivity in the study areas for control purposes than other methods in previous studies [35, 36]. Unlike other studies conducted in western Kenya [13, 19, 35, 36], in all 60 days of every season, aquatic habitat stability had no effect on the estimation of the productivity because habitats remained stable throughout although the size of the habitats varied with the seasons. The sample size used indicates that the efficiency of these sampling methods in estimating productivity can be useful where any of these species are found. These methods can be utilised in areas of Africa where malaria vectors breeding habitats are prevalent and have increased as a result of environmental degradation and land use changes [17, 36, 37]. The efficiency of these methods provides additional information for planning improved larval habitat interventions intending to reduce the adult density effectively . Larval habitat control by use of larvicides has been a major practice in urban malaria control programmes in Africa [25–27, 39], but the use of these methods evaluated here may significantly reduce the costs of larvicides and increase the efficacy of targeted larval habitat control programmes. Previously used methods for estimating habitats productivity such as dipping, netting, and mobile and static quadrant devices, were not effective in terms time spent in sampling and operational costs as evaluated traps [21, 40]. Belleville traps were structurally better but were not useful as the methods evaluated in this study due to their lower ability to estimate the productivity in habitats with deep water, higher grass cover, water turbidity, algal abundance and debris cover which could not affect current methods . The reported methods in this study are more reliable in estimating vectors population in an area of interest for vector intervention than other methods such as human landing catch, pyrethrum spray catch, pit trap, and odour-baited traps [21, 41–43]. The failure of the above mentioned methods in estimating local vector densities seasonally is based on the fact that they estimate only adult densities that might have flown from other locations for host-seeking. A study conducted in West Africa showed aerial samplers and emergence traps to be efficient in estimating the abundance of anopheline and culicine mosquito species during the rainy season in river beds and rice fields .
Seasonality differences were observed to have contributed small changes in estimating of the vector density. These traps' efficiency in estimating habitat productivity was mostly affected by other natural factors such as dry and rainy seasons' variations. In Figure 3, both trapping methods used had lower density estimation efficiency for all species in second dry season (February 2010 - April 2010) than in the first dry season (February 2009 - March 2009). This might be explained as the effect of over washing which was attributed to El Nino (heavy) rains which took place between November, 2009 and January, 2010. The habitats where larval productivity was most affected were the drainage ditches. These ditches are small in size and therefore, are more susceptibility to instability during long periods of drought or flooding during heavy rains. Since An. gambiae ss and An. arabiensis were found predominantly in these ditches, their productivity was affected by the habitat flooding. An. funestus and culicine species productivity estimation were higher as swamps and abandoned goldmines had larger volume of water that could support their gravid female oviposition sites longer than the drainage ditches. It is probable that the heavy rains experienced during the rainy season from June to August 2009 and subsequently El Niño rains affected the productivity of An. gambiae s.s, An. arabiensis and culicine species as drainage ditches and abandoned goldmines were flooded and washed out respectively. Meanwhile An. funestus productivity was not affected due to the large size of the swamps and their ability to retain large amounts of water (Figure 3). In the rainy season from June to August 2008, rainfalls were short and less intense and therefore, facilitated better productivity of An. gambiae s.l in drainage ditches as there was no flooding of ditches resulting in more oviposition sites for the gravid Anopheles females. Productivity of culicine species in this season was also favoured as no flooding in the goldmines leading to increased larvae survivorship and subsequently a higher density of pupae and adults (Figure 3). For An. funestus, the productivity was found to decrease in swamps during this rainy season with the shorter rains due to the fact that the shorter rains could not reduce the amount of brown and green algal cover in habitats and hence the larval abundance. The algal cover favoured culicine species abundance at high levels of decompositions and An. gambiae s.l abundance at lower levels of decomposition.
Despite the efficiency of the aerial sampler, other habitat ecological factors such as chlorophyll a, debris abundance, grass cover and algae abundance were found to influence productivity of mosquitoes. An. gambiae ss and An. arabiensis were found to have significant relationship with presence of debris, low grass cover and algal abundance. In other studies it was found that, An. gambiae s.l larvae use plant debris and algal biomass as a food source thus increasing productivity of the habitats [13, 45]. Lower grass coverage implies that these habitats are exposed to more sunlight. Previously it has been shown that larval habitats in open areas are more likely to have higher water temperatures, resulting into acceleration of the development rate, increased survival of An. gambiae s.l larvae and decreased larval developmental time into adults [13, 45]. An. funestus was found to occur in large water bodies significantly associated with high grass coverage as reported in previously ecological studies [28, 46]. The increase of chlorophyll-a has been associated with increased habitat pollution, a phenomenon that favoured culicine species larvae resulting in higher densities in polluted habitats [13, 37, 46, 47].
Studies done on mosquito dynamics in the highlands of western Kenya in the past, showed that An. gambiae s.s and An. funestus were the main malaria vectors [18, 48, 49]. Previous studies reported An. arabiensis abundance at 1.3% from larval habitats and 3% from indoor spray catch respectively [13, 24]. Another study by Wamae and colleagues  reported An. arabiensis abundance of 10% in the highlands. This study is the first one of its kind to determine An. arabiensis abundance of 33.7% among all the mosquitoes specimens identified by PCR sampled for two years. This can be justified in three ways; first, this high occurrence of An. arabiensis in habitats has been in contrary to the number of adults collected indoors by pyrethrum spray catches (PSC). The low number of An. arabiensis caught by PSC can be attributed to the duality of the An. arabiensis behaviour in that it is zoophagic (i.e. feeding on animals), exophilic (i.e. prefers to rest outdoors) and exophagic (i.e. prefers feeding outdoors) [28, 32, 51]. Secondly, the high rate of bednet ownership and usage reported among community members in this area might have contributed to the low indoor mosquito densities [52, 53]. Thirdly, environmental degradation (including deforestation) and land use type changes might have contributed to the creation of convenient habitats and raised temperatures for An. arabiensis[54, 55]. This study has also revealed that An. arabiensis has reinvaded and as begun establishing itself within the Kenyan highlands due to high abundance observed ever reported.
The aerial sampler information reported here has the potential to reduce the cost of planning for targeted intervention of malaria vectors species as the preferred habitats of each species is already known. This method allow for improved temporal and spatial planning for the monitoring and control of malaria transmission and epidemics in the highlands of western Kenya and in other highlands and urban areas of Africa. Currently, there are major malaria control efforts underway in Africa, and these have demonstrated remarkable success in malaria vectors reduction in selected countries  where these methods can be of additional value in empowering malaria vector larval habitats control programmes.
These results have given additional information in understanding of habitat productivity for targeted intervention and control of anopheline and culicine species in highland and urban areas of Africa. These sampling methods have been proven to have potential for improved efficiency in the planning for integrated vector management.
Conflict of interest
All authors declare they have no competing interest in this study.
Aregawi MW, Ali AS, Al-Mafazy AW, Molteni F, Katikiti S, Warsame M, Njau RJ, Komatsu R, Korenromp E, Hosseini M, Low-Beer D, Bjorkman A, D'Alessandro U, Marc Coosemans M, Otten M: Reductions in malaria and anaemia case and death burden at hospitals following scale-up of malaria control in Zanzibar, 1999-2008. Malar J. 2011, 10: 46. 10.1186/1475-2875-10-46.
Okiro EA, Hay SI, Gikandi PW, Sharif SK, Noor AM, Peshu N, Marsh K, Snow RW: The decline in paediatric malaria admissions on the coast of Kenya. Malar J. 2007, 6: 151. 10.1186/1475-2875-6-151.
O'Meara WP, Mangeni JN, Steketee R, Greenwood B: Changes in the burden of malaria in sub-Saharan Africa. Lancet Infect Dis. 2010, 10: 545-555. 10.1016/S1473-3099(10)70096-7.
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: 32-38.
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.
Zhou G, Minakawa N, Githeko A, Yan G: Spatial distribution patterns of malaria vectors and sample size determination in spatially heterogeneous environments: a case study in the west Kenyan highland. J Med Entomol. 2004, 41: 1001-1009. 10.1603/0022-2585-41.6.1001.
Balls MJ, Bodker R, Thomas CJ, Kisinza W, Msangeni HA, Lindsay SW: Effect of topography on the risk of malaria infection in the Usambara Mountains, Tanzania. Trans R Soc Trop Med Hyg. 2004, 98: 400-408. 10.1016/j.trstmh.2003.11.005.
Kreuels B, Kobbe R, Adjei S, Kreuzberg C, von Reden C, Bater K, Klug S, Busch W, Adjei O, May J: Spatial variation of malaria incidence in young children from a geographically homogeneous area with high endemicity. J Infect Dis. 2008, 197: 85-93. 10.1086/524066.
Garnham PC: The incidence of malaria at high altitudes. J Natl Malar Soc. 1948, 7: 275-284.
Malakooti MA, Biomndo K, Shanks GD: Reemergence of epidemic malaria in the highlands of western Kenya. Emerg Infect Dis. 1998, 4: 671-676. 10.3201/eid0404.980422.
Githeko AK, Ndegwa W: Predicting malaria epidemics in the Kenyan highlands using climate data: a tool for decision makers. Global Change & Human Health. 2001, 2: 54-63. 10.1023/A:1011943131643.
Munga S, Minakawa N, Zhou G, Githeko AK, Yan G: Survivorship of immature stages of Anopheles gambiae s.l. (Diptera: Culicidae) in natural habitats in western Kenya highlands. J Med Entomol. 2007, 44: 758-764. 10.1603/0022-2585(2007)44[758:SOISOA]2.0.CO;2.
Munga S, Minakawa N, Zhou G, Mushinzimana E, Barrack OO, Githeko AK, Yan G: Association between land cover and habitat productivity of malaria vectors in western Kenyan highlands. Am J Trop Med Hyg. 2006, 74: 69-75.
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-277. 10.1603/0022-2585(2005)042[0270:SOAGSS]2.0.CO;2.
Afrane YA, Zhou G, Lawson BW, Githeko AK, Yan G: Effects of microclimatic changes caused by deforestation on the survivorship and reproductive fitness of Anopheles gambiae in western Kenya highlands. Am J Trop Med Hyg. 2006, 74: 772-778.
Minakawa N, Mutero CM, Githure JI, Beier JC, Yan G: Spatial distribution and habitat characterization of anopheline mosquito larvae in Western Kenya. Am J Trop Med Hyg. 1999, 61: 1010-1016.
Githeko AK, Ayisi JM, Odada PK, Atieli FK, Ndenga BA, Githure JI, Yan G: Topography and malaria transmission heterogeneity in western Kenya highlands: prospects for focal vector control. Malar J. 2006, 5: 107. 10.1186/1475-2875-5-107.
Minakawa N, Sonye G, Mogi M, Githeko A, Yan G: The effects of climatic factors on the distribution and abundance of malaria vectors in Kenya. J Med Entomol. 2002, 39: 833-841. 10.1603/0022-2585-39.6.833.
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.
Gimnig JE, Ombok M, Kamau L, Hawley WA: Characteristics of larval anopheline (Diptera: Culicidae) habitats in Western Kenya. J Med Entomol. 2001, 38: 282-288. 10.1603/0022-2585-38.2.282.
WHO: Manual on practical entomology in malaria. Part II.Methods and Techniques. Division of Malaria and Other Parasitic Diseases. 1975, Geneva
Gu W, Novak RJ: Habitat-based modeling of impacts of mosquito larval interventions on entomological inoculation rates, incidence, and prevalence of malaria. Am J Trop Med Hyg. 2005, 73: 546-552.
Gu W, Utzinger J, Novak RJ: Habitat-based larval interventions: a new perspective for malaria control. Am J Trop Med Hyg. 2008, 78: 2-6.
Fillinger U, Ndenga B, Githeko A, Lindsay SW: Integrated malaria vector control with microbial larvicides and insecticide-treated nets in western Kenya: a controlled trial. Bull World Health Organ. 2009, 87: 655-665. 10.2471/BLT.08.055632.
Castro MC, Tsuruta A, Kanamori S, Kannady K, Mkude S: Community-based environmental management for malaria control: evidence from a small-scale intervention in Dar es Salaam, Tanzania. Malar J. 2009, 8: 57. 10.1186/1475-2875-8-57.
Yohannes M, Haile M, Ghebreyesus TA, Witten KH, Getachew A, Byass P, Lindsay SW: Can source reduction of mosquito larval habitat reduce malaria transmission in Tigray, Ethiopia?. Trop Med Int Health. 2005, 10: 1274-1285. 10.1111/j.1365-3156.2005.01512.x.
Majambere S, Lindsay SW, Green C, Kandeh B, Fillinger U: Microbial larvicides for malaria control in The Gambia. Malar J. 2007, 6: 76. 10.1186/1475-2875-6-76.
Gillies T, Coetzee M: Supplement of the Anopheles of Africa South of Sahara (Afro tropical Region). 1987, Johannesburg, Republic of South Africa: Publication of The South Africa Institute of Medical Research
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-529.
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-811.
Darsie RFJ, Ward RA, CC C: Identification and geographical distribution of the mosquitoes of North America, north of Mexico. 1981, Fresno, California American Mosquito Control Association, 1:
Coetzee M, Craig M, le Sueur D: Distribution of African malaria mosquitoes belonging to the Anopheles gambiae complex. Parasitol Today. 2000, 16: 74-77. 10.1016/S0169-4758(99)01563-X.
Ndenga BA, Simbauni JA, Mbugi JP, Githeko AK, Fillinger U: Productivity of malaria vectors from different habitat types in the Western kenya highlands. PLoS One. 2011, 6: e19473. 10.1371/journal.pone.0019473.
Minakawa N, Sonye G, Mogi M, Yan G: Habitat characteristics of Anopheles gambiae s.s. larvae in a Kenyan highland. Med Vet Entomol. 2004, 18: 301-305. 10.1111/j.0269-283X.2004.00503.x.
Fillinger U, Sonye G, Killeen GF, Knols BG, Becker N: The practical importance of permanent and semipermanent habitats for controlling aquatic stages of Anopheles gambiae sensu lato mosquitoes: operational observations from a rural town in western Kenya. Trop Med Int Health. 2004, 9: 1274-1289. 10.1111/j.1365-3156.2004.01335.x.
Himeidan YE, Zhou G, Yakob L, Afrane Y, Munga S, Atieli H, El-Rayah el A, Githeko AK, Yan G: Habitat stability and occurrences of malaria vector larvae in western Kenya highlands. Malar J. 2009, 8: 234. 10.1186/1475-2875-8-234.
Munga S, Yakob L, Mushinzimana E, Zhou G, Ouna T, Minakawa N, Githeko A, Yan G: 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-1084. 10.4269/ajtmh.2009.09-0156.
John CC, Riedesel MA, Magak NG, Lindblade KA, Menge DM, Hodges JS, Vulule JM, Akhwale W: Possible interruption of malaria transmission, highland Kenya, 2007-2008. Emerg Infect Dis. 2009, 15: 1917-1924.
Fillinger U, Lindsay SW: Suppression of exposure to malaria vectors by an order of magnitude using microbial larvicides in rural Kenya. Trop Med Int Health. 2006, 11: 1629-1642. 10.1111/j.1365-3156.2006.01733.x.
Knight KL: Quantitative Methods for Mosquito Larval Surveys. J Med Entomol. 1964, 1: 109-115.
Kweka EJ, Mwang'onde BJ, Kimaro E, Msangi S, Massenga CP, Mahande AM: A resting box for outdoor sampling of adult Anopheles arabiensis in rice irrigation schemes of lower Moshi, northern Tanzania. Malar J. 2009, 8: 82. 10.1186/1475-2875-8-82.
Govella NJ, Chaki PP, Mpangile JM, Killeen GF: Monitoring mosquitoes in urban Dar es Salaam: Evaluation of resting boxes, window exit traps, CDC light traps, Ifakara tent traps and human landing catches. Parasit Vectors. 2011, 4: 40. 10.1186/1756-3305-4-40.
Kweka EJ, Mwang'onde BJ, Mahande AM: Optimization of odour-baited resting boxes for sampling malaria vector, Anopheles arabiensis Patton, in arid and highland areas of Africa. Parasit Vectors. 2010, 3: 75. 10.1186/1756-3305-3-75.
Fillinger U, Sombroek H, Majambere S, van Loon E, Takken W, Lindsay SW: Identifying the most productive breeding sites for malaria mosquitoes in The Gambia. Malar J. 2009, 8: 62. 10.1186/1475-2875-8-62.
Merritt R, Dadd RH, Walker ED: Feeding behavior, natural food, and nutritional relationships of larval mosquitoes. Annu Rev Entomol. 1992, 37: 349-376. 10.1146/annurev.en.37.010192.002025.
Mwingira V, Mayala BK, Senkoro KP, Rumisha SF, Shayo EH, Mlozi MR, Mboera LE: Mosquito larval productivity in rice-fields infested with Azolla in Mvomero District, Tanzania. Tanzan J Health Res. 2009, 11: 17-22.
Mwangangi JM, Muturi EJ, Shililu J, Muriu SM, Jacob B, Kabiru EW, Mbogo CM, Githure J, Novak R: Contribution of different aquatic habitats to adult Anopheles arabiensis and Culex quinquefasciatus (Diptera: Culicidae) production in a rice agroecosystem in Mwea, Kenya. J Vector Ecol. 2008, 33: 129-138. 10.3376/1081-1710(2008)33[129:CODAHT]2.0.CO;2.
Ndenga B, Githeko A, Omukunda E, Munyekenye G, Atieli H, Wamai P, Mbogo C, Minakawa N, Zhou G, Yan G: Population dynamics of malaria vectors in western Kenya highlands. J Med Entomol. 2006, 43: 200-206. 10.1603/0022-2585(2006)043[0200:PDOMVI]2.0.CO;2.
Shililu JI, Maier WA, Seitz HM, Orago AS: Seasonal density, sporozoite rates and entomological inoculation rates of Anopheles gambiae and Anopheles funestus in a high-altitude sugarcane growing zone in Western Kenya. Trop Med Int Health. 1998, 3: 706-710. 10.1046/j.1365-3156.1998.00282.x.
Wamae PM, Githeko AK, Menya DM, Takken W: Shading by Napier Grass Reduces Malaria Vector Larvae in Natural Habitats in Western Kenya Highlands. Ecohealth. 2010
Sinka ME, Bangs MJ, Manguin S, Coetzee M, Mbogo CM, Hemingway J, Patil AP, Temperley WH, Gething PW, Kabaria CW, Okara RM, Van Boeckel T, Godfray HC, Harbach RE, Hay SI: 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. 10.1186/1756-3305-3-117.
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.
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.
Chen H, Githeko AK, Zhou G, Githure JI, Yan G: New records of Anopheles arabiensis breeding on the Mount Kenya highlands indicate indigenous malaria transmission. Malar J. 2006, 5: 17. 10.1186/1475-2875-5-17.
Githinji S, Herbst S, Kistemann T: The human ecology of malaria in a highland region of South-West Kenya. Methods Inf Med. 2009, 48: 451-453. 10.3414/ME9242.
Mayor S: WHO report shows progress in efforts to reduce malaria incidence. Bmj. 2008, 337: a1678. 10.1136/bmj.a1678.
The help of the field assistants from Mbale village and the highland malaria programme field team at KEMRI is highly appreciated for data collection throughout the study. Dr. Afrane Yaw is acknowledged for field logistics and organisation during this study. This work had financial support from American National Institute of Health (D43 TW01505 and R01 AI-50243 to GY). This paper is published with the permission of the Director General, Kenya Medical Research Institute, (KEMRI).
EJK, AKG, SM and GY conceived, designed and implemented the study; EJK, TMG and SM conducted field experiments. AKG and GY did the supervision of the project activities. EJK, ML and GZ did data analysis and interpretation. EJK wrote this manuscript. AKG, GZ, FM, EJK and GY revised the manuscript. All authors edited and approved the final version.