Habitat stability and occurrences of malaria vector larvae in western Kenya highlands
Malaria Journalvolume 8, Article number: 234 (2009)
Although the occurrence of malaria vector larvae in the valleys of western Kenya highlands is well documented, knowledge of larval habitats in the uphill sites is lacking. Given that most inhabitants of the highlands actually dwell in the uphill regions, it is important to develop understanding of mosquito breeding habitat stability in these sites in order to determine their potential for larval control.
A total of 128 potential larval habitats were identified in hilltops and along the seasonal streams in the Sigalagala area of Kakamega district, western Kenya. Water availability in the habitats was followed up daily from August 3, 2006 to February 23, 2007. A habitat is defined as stable when it remains aquatic continuously for at least 12 d. Mosquito larvae were observed weekly. Frequencies of aquatic, stable and larvae positive habitats were compared between the hilltop and seasonal stream area using χ2-test. Factors affecting the presence/absence of Anopheles gambiae larvae in the highlands were determined using multiple logistic regression analysis.
Topography significantly affected habitat availability and stability. The occurrence of aquatic habitats in the hilltop was more sporadic than in the stream area. The percentage of habitat occurrences that were classified as stable during the rainy season is 48.76% and 80.79% respectively for the hilltop and stream area. Corresponding frequencies of larvae positive habitats were 0% in the hilltop and 5.91% in the stream area. After the rainy season, only 23.42% of habitat occurrences were stable and 0.01% larvae positive habitats were found in the hilltops, whereas 89.75% of occurrences remained stable in the stream area resulting in a frequency of 12.21% larvae positive habitats. The logistic regression analysis confirmed the association between habitat stability and larval occurrence and indicated that habitat surface area was negatively affecting the occurrence of An. gambiae larvae. While An. gambiae and An. funestus larvae occurred throughout the study period along the streams, a total of only 15 An. gambiae larvae were counted in the hilltops, and no An. funestus were found. Moreover, no larvae managed to develop into adults in the hilltops, and the density of adult An. gambiae was consistently low, averaging at 0.06 females per house per survey.
The occurrence of malaria vector larvae in the hilltop area was uncommon as a result of the low availability and high instability of habitats. To optimize the cost-effectiveness of malaria interventions in the western Kenya highlands, larval control should be focused primarily along the streams, as these are likely the only productive habitats at high altitude.
Anopheline larvae occupy habitats that vary widely in their physicochemical properties, surface areas and vegetation. Habitats can be highly species specific and extreme variation in larval production has been shown to occur on a small scale both spatially and temporally [1–5]. Informed larval interventions that target the more prolific habitats have been suggested to have great potential in combating malaria [6–9], especially at a local, rather than nationwide, scale. Intervention therefore needs to be adaptable to, and specific for local communities.
A recent series of malaria epidemics have occurred in the western Kenya highlands [10–12]. Although year-round the valleys constitute the majority of vector breeding sites, potentially important larval habitats can also be found in the uphill areas during the rainy season [13, 14]. Anopheles gambiae is considered a typical pioneer species and it has been shown to be capable of colonizing these uphill habitats within just a few days of the rain falls . This colonization by blood-fed female mosquitoes brings with it the risk of malaria infection to the immunologically naïve, uphill human population . However, extending breeding site interventions to account for this spatial expansion in the vector population is not necessarily justified. Larvae of An. gambiae usually occur in small, temporary and sunlit pools [17–19], which are prone to drying up quickly . In order to justify larval control in the uphill regions during the rainy season, it must be determined whether the transient pools of water in these locations are found are sufficiently stable for the locally recorded threshold of 12 days necessary for larvae to develop into adults . It was the objective of this study to develop the knowledge of the association between habitat stability and vector productivity in the uphill regions of western Kenyan highlands, and to elucidate the role of these habitats in local malaria transmission.
Materials and methods
The study was conducted in Sigalagala (0°12'N, 34°46'E) area in Kakamega district, western Kenya highlands with elevation of 1500-1600 m above sea level. Distinct topographical features of the study villages typically consist of hills and valleys. The hillside is mostly dotted with maize plantations. A few patches of indigenous forests are located along the valleys. Six different seasonal streams flow within the valleys in the study villages and empty into the main Yala river (Figure 1). The inhabitants live in traditional stick and mud houses with thatch or iron sheet roof. Most residences in the area are located uphill, with a few on the hillsides. The areas within 50 m of the nearest stream were considered to be the 'stream area', otherwise they were considered as 'hilltop' (Figure 1). The area experienced a short rainy season from August to November 2006 and the total rainfall was 1050.2 mm during the period of study from August 3, 2006 - February 23, 2007 (Figure 2).
Habitat monitoring and larval sampling
A total number of 128 aquatic habitats were identified within the study area, with 47 in the hilltops and 81 in the stream area (Figure 1). These habitats were used as fixed observation sites throughout the study period and no additional habitats were added into the study after the initial stage. The habitats were mapped using a Global Positioning System (GPS) unit. Individual habitats were numbered and their surface area recorded along with the land use type. The water availability in each habitat was observed daily. Mosquito larval surveys were done weekly in all study sites. The immature stages were sampled using the standard dipping method (dipper size of 350 ml) . Water was dipped up to 20 times. When a habitat was too small to make 20 dips, water was dipped as many times as possible. Mosquito larval species were identified , counted and returned to the habitats. Samples of young larvae of An. gambiae s.l. were identified to species using the polymerase chain reaction (PCR) .
Adult malaria vectors were monitored in the uphill area by collection from a randomly selected sample of twenty-eight houses (Figure 1). These 28 houses were fixed for adult samplings throughout the study period. All selected houses were numbered and geo-referenced using the GPS. Adult mosquitoes were collected biweekly from August 2006 to February 2007 using the indoor pyrethrum spray catches . The mosquitoes collected were identified to species morphologically .
The number of days that a habitat remained aquatic and the number of times that immature stages of An. gambiae were present for each individual habitat over the 205 d study period were recorded. To determine the impact of rainfall on habitat stability and larval occurrence, the study period was divided into two: the rain season lasting from August to the end of November, 2006, and the post-rain season lasting from the 1st of December, 2006 to the 23rd February, 2007. A habitat was defined as stable if it remained aqueous for 12 d or longer, following on from a previous study in the same area that showed the egg-to-adult development cycle of An. gambiae requires at least this length of time . The mean densities of immature stages (larvae/dip) were calculated per habitat. The χ2-test was used to compare the frequencies of aquatic, stable and larvae positive habitats between hilltop and stream areas. Tukey-Kramer HSD test of repeated measure MANOVA was used to compare the mean occurrence of anopheline larvae between different seasons and different topographic areas. Multiple logistic regression with backward stepwise analysis was then performed to identify the factors which significantly influenced the occurrences (presence/absence) of An. gambiae larvae. Any independent variable with a significance level of < 95% was eliminated from the final model. The analysis was conducted using JMP statistical software (JMP SAS Institute Inc. 2003).
Species composition and habitat types of malaria vector larvae
A total of 12,410 mosquito larvae were collected, of which 1,675 (13.5%) were Anopheles spp. Among the anopheline larvae collected, 61.3% (n = 1027) were An. gambiae s.l., 37.2% (n = 623) were Anopheles funestus, 0.8% (n = 14) were Anopheles coustani, 0.4% (n = 7) were Anopheles squamosus and 0.2% (n = 4) were Anopheles implexus. Due to the small number of Anopheles arabiensis (n = 17) identified by PCR, the An. arabiensis and An. gambiae s.s. were pooled together.
Among the larvae collected, about 98.5% (n = 1012) of An. gambiae and all An. funestus were collected from the stream areas. A large proportion of An. gambiae (52.2%) and An. funestus (73.2%) were collected from the drainage ditch habitats. Other habitat types in the stream area where An. gambiae was found include man made pools (22.2%), natural swamps (17.1%) and hoof footprints (8.5%). Only fifteen larvae of An. gambiae were found in the hilltop area, specifically, in four tire track habitats. No An. funestus larvae were collected from hilltops. A total of twenty-seven An. gambiae pupae were found in the stream areas.
Effect of rainfall and topography on habitat stability
In general, the trend in aquatic habitats closely followed the pattern of rainfall; peak habitat frequencies were concurrent with peak rainfall (Figure 2). However, the extent to which habitat frequencies fluctuated between seasons varied greatly between the two topographical regions. While the frequencies of aquatic habitats in the stream area remained relatively static over the study period (fluctuating by no more than 20%), habitat frequency fluctuated by up to 100% in the hilltop areas.
Concordantly, the mean length of days that a habitat remained aquatic was significantly shorter in the hilltop area compared to the stream areas (16.7 d vs 50.2 d during rain and 15.7 d vs 64.0 d after rain) (Table 1). The mean number of times that a habitat dried up in the hilltop area was double that of the stream areas during the rainy season (3.8 vs 2.0), and six times more than in the stream areas after the rainy season (2.6 vs 0.4) (Table 1). Overall, the percentage of habitats that were aquatic on any given day in the hilltop area dropped significantly from 63.77% during the rainy season to 44.62% after the rainy season, whereas these percentages remained unchanged in the stream area regardless of season (91.8% during rain season vs 89.9% after rain season) (Table 1).
In total, habitats remained aquatic for 12 d or more only about 30% of the time in hilltops. The significantly greater longevity of habitats in the stream areas is evident from Figure 3 (χ2 = 856.12, d.f. = 119, P < 0.001 during the rainy season and χ2 = 1300.98, d.f. = 83, P < 0.001 after the rainy season).
Effect of topography and rainfall on the occurrences of malaria vectors
Overall, malaria vector larvae and pupae occurred mostly in the stream areas. Larvae of An. gambiae were only observed four times in the hilltops (in January 2007). Each time we only found one larvae positive habitat and none of these larvae developed into pupae. Pre-adult stages of An funestus were not observed in the hilltop areas throughout the entire study period.
Larvae of An. gambiae were consistently found in the stream areas throughout the study period (Figure 4). However, significantly fewer habitats were found to be larval positive during the rainy season than after (5.91% vs. 12.21%, repeated measure MANOVA Tukey-Kramer HSD test, P < 0.05). Anopheles gambiae pupae occurrence frequencies were much lower and so no significant difference was found between topographies or seasons (0.85% vs. 0.57% during and after rainy season, Tukey-Kramer HSD test, P > 0.05). The mean density of An. gambiae larvae was significantly higher in the stream areas (mean of 0.27 per dip, 95% CI [0.21, 0.32]) compared to the hilltops (mean of 0.02, 95% CI [0.00, 0.10]). The 1.95% prevalence of larval An. funestus in the stream areas during the rainy season was shown to significantly increase to 7.66% after the rainy season (Figure 4, Tukey-Kramer HSD test, P < 0.05). The occurrence of An. funestus pupae was negligible, with only three sightings within the stream areas in total.
The results of logistic regression analyses indicated that the occurrence (presence/absence) of An. gambiae larvae was significantly and positively associated with both the topography (hilltops = 0 and stream area = 1) and the number of days that a habitat remained aquatic, but negatively associated with habitat surface area (χ2 = 84.447. d.f. = 3, P < 0.0001). The effect of habitat types and land use types on the occurrence of mosquito larval was not significant.
Adult mosquito dynamics
During the study period, adult An. gambiae appeared sporadically in the area (Figure 5), and were only found inside one or two houses during each survey. A total of 21 female adult An. gambiae were collected during the sampling period, and the mean density was estimated at 0.06 females per house per survey. Only one adult An. funestus was ever collected (mid-December, 2006).
The present study aimed to investigate the occurrence of malaria vector larvae in relation to topography and habitat stability in order to assess the applicability of a more targeted larval vector control in the highlands of western Kenya. Specifically, the intention was to measure the malaria vector productive potential of breeding sites located in high altitude regions (1,500-1,600 m above sea level). Breeding site productivity studies normally take place at lower altitudes or major river valley area where both pre-adult and adult vectors are found at higher, and more consistent frequencies [2, 4, 14, 16, 24]. However, these studies fail to consider critical aspects of vector ecology at the higher altitudes, where human populations are more crowded and more immunologically naïve to infection.
As expected, we found more vector breeding sites nearby streams than along hilltops in our study. The relatively few breeding sites that were found along hilltops greatly fluctuated in frequency between seasons and were highly dependent on rainfall. The increased stability of breeding sites found within the stream areas is simply a result of the more permissive conditions these areas provide for water accumulation [15, 25]. Not only was the duration of breeding site viability significantly reduced in the hilltops, but the frequency of anopheline larvae positive habitats was also significantly lower. Very few An. gambiae larvae were found in hilltop breeding sites, with a zero pupation rate. Both hilltops and stream areas received equal rainfall, were of identical temperature and had very similar land use. The significant variation between the vector productivity of the two topographically distinct areas is attributed to the significant differences in habitat stability.
Similar observations have been reported in lowlands area of the Lake Victoria Basin in western Kenya, where habitat stability was positively correlated with larval density and pupal productivity . As the Basin is very flat, Anopheles larval habitats are widely dispersed and sufficiently stable to be productive even hundreds of meters from the lake [5, 24]. The hilly landscape at higher altitudes of the western Kenyan highlands only allows for smaller, discrete areas of sustained water accumulation that are mostly located nearby streams [26–28]. These constitute hotspots for productive anopheline habitats which, given their confined spatiality, might justify localized larval control.
The superior capacity of An. gambiae as a colonizer of new territory as illustrated by our study corroborates what is already known for this species [29–32]. Adult malaria vectors were routinely collected from human dwellings of the uphill sites, albeit at low densities. While the uphill habitats appear to be inappropriate targets for vector control, the stream areas of the western Kenyan highlands appear to act as stepping stones between the valleys and the regions of high altitude. There is a tendency for vectors to move to higher altitudes following the rainy season as a result of the increased availability of breeding sites in these areas. Therefore, in eliminating the stepping stone provided by the habitats of the stream areas at high altitude, the extended adult vector distribution might be prevented and the threat of epidemic malaria in the western Kenyan highlands removed.
Shililu J, Ghebremeskel T, Seulu F, Mengistu S, Fekadu H, Zerom M, Ghebregziabiher A, Sintasath D, Bretas G, Mbogo C, Githure J, Brantly E, Novak R, Beier JC: Larval habitat diversity and ecology of anopheline larvae in Eritrea. J Med Entomol. 2003, 40: 921-929. 10.1603/0022-2585-40.6.921.
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.
Sattler MA, Mtasiwa D, Kiama M, Premji Z, Tanner M, Killeen GF, Lengeler C: Habitat characterization and spatial distribution of Anopheles sp. mosquito larvae in Dar es Salaam (Tanzania) during an extended dry period. Malar J. 2005, 4: 4-10.1186/1475-2875-4-4.
Munga S, Minakawa N, Zhou G, Mushinzimana E, Barrrack OO, Githeko AK, Yan G: Association between land cover and habitat productivity of malaria vectors in Western Kenya highlands. Am J Trop Med Hyg. 2006, 74: 69-75.
Mutuku FM, Bayoh MN, Gimnig JE, Vulule JM, Kamau L, Walker ED, Kabiru E, Hawley WA: Pupal habitat productivity of Anopheles gambiae complex mosquitoes in a rural village in western Kenya. Am J Trop Med Hyg. 2006, 74: 54-61.
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-52.
Depinay J-MO, Mbogo CMN, Killeen GF, Knols BGJ, Beier JC, Carlson J, Dushoff J, Billingsely P, Mwambi H, Githure JI, Toure AM, McKenzie FE: A simulation model of African Anopheles ecology and population dynamics for the analysis of malaria transmission. Malar J. 2004, 3: 29-50. 10.1186/1475-2875-3-29.
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.
Killeen GF, Tanner M, Mukabana RW, Kalongolela MS, Kannady K, Lindsay SW, Fillinger U, de Castro MC: Habitat targeting for controlling aquatic stages of malaria vectors in Africa. Am J Trop Med Hyg. 2006, 74: 517-518.
Snow RW, Ikoku A, Omumbo J, Ouma J: The Epidemiology, Politicsand Control of Malaria Epidemics in Kenya: 1900-1998. 1999, Roll Back Malaria report. World Health Organization, Geneva
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.
Zhou G, Minakawa N, Githeko AK, Yan G: Association between climate variability and malaria epidemics in the East African highlands. PNAS. 2004, 101: 2375-2380. 10.1073/pnas.0308714100.
Minakawa N, Sonye G, Mogi M, Yan G: Habitat characteristics of Anopheles gambiae s.s. larvae in a Kenyan highland. Med Vet. 2004, 18: 301-305.
Minakawa N, Munga S, Atieli F, Mushinzimana E, Zhou G, Githeko AK, Yan G: 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-165.
Minakawa N, Sonye G, Yan G: Relationships between occurrences of Anopheles gambiae s.l. (Diptera: Culicidae) and size and stability of larval habitats. J Med Entomol. 2005, 42: 295-300. 10.1603/0022-2585(2005)042[0295:RBOOAG]2.0.CO;2.
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.
Gillies MT, De Meillon B: The Anophelinae of Africa south of the Sahara. 1968, Johannesburg: South African Institute for Medical Research Press
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.
Gimnig JE, Ombok M, Kamau L, Hawley WA: Characteristics of larval anopheline (Diptera:Culicide) habitats in Western Kenya. J Med Entomol. 2001, 38: 282-288. 10.1603/0022-2585-38.2.282.
Minakawa N, Githure JI, Beier JC, Yan G: Anopheline mosquito survival strategies during the dry period in western Kenya. J Med Entomol. 2001, 38: 388-392. 10.1603/0022-2585-38.3.388.
WHO: Manual on Practical Entomology in Malaria. Part II. Methods and Techniques. No. 13. Geneva. 1975
Gillies MT, Coetzee M: A supplement to the Anophelinae of Africa South of the Sahara (Afrotropical Region) (publication no. 55). 1987, Johannesburg: South Africa Institute for Medical Research
Scott JA, Brogdon WG, Collins FH: Identification of single species of the Anopheles gambiae complex by the polymerase chain reaction. Am J Trop Med Hyg. 1993, 49: 520-529.
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.
Carter R, Mendis KN, Roberts D: Spatial targeting of interventions against malaria. Bull World Health Organ. 2000, 78: 1401-1411.
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. 10.1186/1476-072X-7-50.
Li L, Bian L, Yan G: A study of the distribution and abundance of the adult malaria vector in western Kenya highlands. Int J Health Geogr. 2008, 7: 50-10.1186/1475-2875-7-40.
Cohen JM, Ernst KC, Lindblade KA, Vulule JM, John CC, Wilson ML: Topography-derived wetness indices are associated with household-level malaria risk in two communities in the western Kenyan highlands. Malaria J. 2008, 7: 40-10.1016/0035-9203(79)90036-1.
Coluzzi M, Sabatini A, Petrarca V, Di Deco MA: Chromosomal differentiation and adaptation to human environment in the Anopheles gambiae complex. Trans R Soc Trop Med Hyg. 1979, 37: 483-497. 10.1016/0035-9203(79)90036-1.
Toure YT, Petrarca V, Traore SF, Coulibaly A, Maiga HM, Sankare O, Sow M, Di Deco MA, Coluzzi M: The distribution and inversion polymorphism of chromosomally recognized taxa of Anopheles gambiae complex in Mali, West Africa. Parasitologia. 1998, 40: 477-511. 10.1603/0022-2585(2005)042[0993:OSPAEH]2.0.CO;2.
Munga S, Minakawa N, Zhou G, Okey-Owuor J, Barrrack OO, Githeko AK, Yan G: Oviposition site preference and egg hatchability of Anopheles gambiae: effects of land cover types. J Med Entomol. 2005, 42: 993-997. 10.1603/0022-2585(2005)042[0270:SOAGSS]2.0.CO;2.
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.
The authors thank Peter Wamai, Amos Wabwile for technical assistance in the field. This work was supported by NIH grant (D43 TW01505 and R01 AI-50243 to GY).
The authors declare that they have no competing interests.
YH carried out the field surveys, assembled data and drafted the manuscript. GZ performed statistical analysis and revised the manuscript. LY revised the manuscript. YA, SM and HA participated in the field surveys. AR revised the manuscript. AG and GY participated in the study design and coordination. All authors read and approved the final manuscript.