|
|
||||||||
| ABSTRACT |
|
|
|---|
| INTRODUCTION |
|
|
|---|
In recent years, Kenya has experienced very rapid human population growth. The population of this country has almost doubled since 1980.12 Population growth in the western Kenya highlands is a particularly severe problem, due in part to a lack of family planning and to migration from other areas. As a consequence of dramatic human population increase, there have been unprecedented land-use changes in the highlands. For example, since 1965 the Malava Forest in the Kaka-mega district has been reduced from 600 to less than 100 hectares.13 Most rain forests have been cleared for crop planting, cattle grazing, commercial logging, firewood collection, and housing construction.14
The drastic land-use alteration can promote vector-borne disease transmission in several ways.15,16 First, more human-made aquatic habitats become available for the Anopheles gambiae complex, which is the primary malaria vector in Africa. For example, increased cattle grazing creates more open temporary habitats that may serve as mosquito breeding habitats.17,18 Second, the physical and chemical properties of mosquito larval habitats may change. Third, the microclimate of mosquito larval habitats may change. Habitats well exposed to sunlight may have a higher temperature than those in shaded areas, and thus mosquito larval development rate may be increased. For example, a higher adult density of the An. gambiae complex was found in houses near cultivated swamps than in those near natural swamps in a Ugandan highland site.9 Natural swamps are primarily populated with tall grasses such as papyrus, whereas cultivated swamps are planted with agricultural crops such as maize.18 Since cultivated swamps generally receive more exposure to sunlight than those in natural swamps, the ambient air temperature in the cultivated swamp area was significantly higher than that in natural swamps.9 Moreover, mosquito larval predators may be more prevalent in natural swamps than in cultivated swamps.19 These differences in aquatic habitat characteristics may have important effects on anopheline larval development and survivorship, and on the distribution of suitable larval habitats.
The aim of the present study is to examine the effects of land use/land cover types on the distribution of anopheline larval habitats in a western Kenya highland site. We determined the spatial distribution pattern of larval habitats and analyzed several physical and chemical components in the aquatic habitats that may be affected by land use and may be relevant to mosquito larval development. We also examined the effects of topography on the spatial distribution of larval habitats. The east African highland region contains numerous small valleys and basin-like depressions in a plateau.7 Stagnant water bodies are more likely to form in valley bottoms and basin-like depressions than in the hills.18 Thus topography may be an important factor for the distribution of aquatic habitats in the highlands. Knowledge of the effects of land use/land cover types and topography on mosquito larval habitat distribution is useful for predicting the distribution and abundance of malaria vectors and the impact of land-use practice on malaria transmission, and thus for designing novel strategies for malaria intervention.
| MATERIALS AND METHODS |
|
|
|---|
Land cover types. Prior to the larval habitat survey, we estimated the surface area (meters2) of each land cover type within the study area. Land cover types were classified into five categories: farmlands, forests, pastures, roads, and natural swamps. Farmlands were characterized by the presence of any agricultural crops and bare ground that had been prepared for planting crops. Pastures were grasslands used for grazing. An area with a mixture of grass and shrubs was also classified as pasture. Natural swamps were characterized by the presence of emergent aquatic plants. Forests and natural swamps were the categories with the least modification by anthropogenic activities. A multispectral one-meter ground resolution Ikonos image (Space Imaging, Atlanta, GA) of the study area was acquired on April 12, 2001. Each land cover type was visually digitized using the geographic information system (GIS) software package ArcView 3.3 (Environmental Systems Research Institute, Redlands, CA). The digitized land covers were validated by ground truthing. The surface area of each land cover type was estimated by first projecting the land cover layer into Universal Transverse Mercator 1983 zone 36N and then using an Avenue Script that calculates the surface area in ArcView. All streams (including the Yala River) and roads in the study area were also mapped and digitized using the satellite image in the GIS and validated by ground truthing.
Habitat characterization. We first counted and categorized all aquatic habitats (excluding running water and water in containers in houses and in tree holes) with respect to land cover types in the long rainy season (May 17 to May 30, 2002) and dry season (August 13 to August 21, 2002). Land cover types, where aquatic habitats were located, were classified into the above five categories. The location and elevation of each habitat was recorded using the global positioning system in differential mode.21 Occurrence of anopheline larvae was examined at each habitat using a standard 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.
We also characterized 130 and 50 randomly selected aquatic habitats (including both anopheline-positive and -negative habitats) with respect to physicochemical variables in the long rainy season and dry season, respectively. The selected habitats were characterized using the following variables: 1) canopy cover, 2) presence of filamentous algae, 3) turbidity, 4) habitat size (volume), 5) pH, and 6) concentrations of ammonium (NH4+), nitrate (NO3), and orthophosphate (PO43). Canopy cover was estimated with a spherical densiometer.22 Surface water was sampled and analyzed for the above chemical variables and turbidity following the standard methods of the U.S. Environmental Protection Agency.23,24 The volume (meters3) of an aquatic habitat of 10 evenly spaced grids was established by measuring its length, mean width, and mean depth, and calculated as the product of these three indices.
In the habitat characterization with the physicochemical variables, anopheline larvae were collected and preserved in 95% ethanol from each positive habitat, and third and fourth instar larvae were identified to the An. gambiae complex, the An. funestus complex, and the other anopheline species based on a morphologic key using a microscope.25 Larvae of the An. gambiae complex and the An. funestus complex were further analyzed with the ribosomal DNA (rDNA) polymerase chain reaction (PCR) method.26,27 First and second instar larvae were subjected to rDNA-PCR analysis for species identification directly.
Data analysis. We determined whether the distribution pattern of anopheline larval habitats was clustered in the study area in each season using point pattern spatial analysis (Ripleys K-function).28 In this analysis, we only used the data from the anopheline-positive habitats. When it was clustered, we determined whether topography influenced the distribution pattern of anopheline larval habitats using two geographic variables, the distance of each habitat to the nearest stream and its elevation. The Kolmogorov-Smirnov one-sample test was used to determine whether anopheline larval habitats were confined toward valley bottoms in each season. The distance of each larval habitat to the nearest stream was estimated in the GIS. Similarly, we examined the relationships between elevation and the distribution of anopheline habitats. We further determined whether the occurrence (presence/absence) of anopheline larvae was significantly correlated with these geographic variables and season using logistic regression analysis.
Since anopheline breeding sites are generally reduced during the dry season,29,30 we suspected that the reduction is reflected in the difference in habitat distribution patterns between the dry season and rainy season. We used bivariate point pattern analysis to determine whether the distribution of anopheline-positive habitats in the dry season was more confined than in the rainy season using the statistical software package Programita.28,31 When it was significantly clustered, we used a t-test to compare the average distances from the nearest streams and the average elevations of anpheline-positive habitats between the seasons.
Since occurrence of the An. funestus complex larvae is generally limited to larger, semipermanent or permanent habitats with aquatic vegetation,30,32,33 its habitat distribution pattern may be different from the An. gambiae complex. We determined whether the distribution of the An. funestus complex-positive habitats was significantly more clustered than that of the An. gambiae complex using bivariate K-function analysis. When it was significantly clustered, we compared the average distances of larval habitats to the nearest streams and the average elevations of larval habitats between the two species complexes. In the analysis for the species complexes, the data from both seasons were included because the number of anopheline habitats during the dry season was not large enough to perform a separate analysis.
Associations between the occurrence of anopheline larval habitats and the land cover types were examined for each season using the chi-square test. The Tukey-type multiple comparison test was used to test differences in proportion of larval occurrence among the land cover types.34 The association was also examined for each species complex. In the analysis for each species complex, the data from both seasons were included together and road was excluded because few positive sites were found within this land cover type. Logistic regression analysis was used to examine the association between occurrence of each species complex and environmental variables. We further examined whether the characteristics of larval habitats (including both positive and negative sites) varied significantly among the land cover types using two-way analysis of variance for continuous variables and chi-square analysis for categorical variables. Tukey multiple comparison tests were used for post-hoc analyses. The statistical significance was 0.05 and adjusted with the Bonferroni correction if necessary. All continuous variables were log-transformed.
| RESULTS |
|
|
|---|
Relationships with topography.
During the rainy season survey, 767 aquatic habitats were identified in the study area, and anopheline larvae were found at 277 habitats (36.1%; Table 1
and Figure 1
). The number of aquatic habitats was reduced to 190, and 48 habitats (25.3%) were anopheline-positive during the dry season. Figure 2
shows measures of the observed L(d) and the 95% confidence interval plotted for values of interpoint distance ranging between 0.1 and 1.4 km. For clustering to be significant at the P < 0.05 level, the observed L(d) must lie above or below the 95% confidence interval lines. We found that anopheline larval habitats were significantly clustered up to 1.4 km during the rainy season, and clustering peaked at 0.70.9 km (Figure 2A
). During the dry season, larval habitat clustering occurred up to 1.0 km and peaked at 0.7 km (Figure 2B
). Bivariate K-function analysis found that the magnitude of the K-function estimate of larval habitats for the dry season was significantly greater than that for the rainy season at distances up to 0.3 km and approximately 0.6 km (Figure 2C
). That is, larval habitat during the dry season was significantly more clustered compared with the rainy season at the distances. Similarly, An. funestus larval habitats were significantly more clustered than the An. gambiae complex habitats in the rainy season (Figure 2D
).
|
|
|
|
|
|
|
|
2 = 36.69, df = 4, P < 0.001), but not significant during the dry season (
2 = 5.13, df = 4, P = 0.274; Table 1
We found that the larvae of the An. gambiae complex were found more frequently in farmlands and pastures than in larval habitats in the forests (
2 = 10.53, df = 3, P = 0.015), while An. funestus larvae were found more frequently in swamps and pastures than in farmlands and forests (
2 = 9.85, df = 3. P = 0.020; Table 4
). Among the eight physicochemical variables examined, canopy cover was the only factor significantly associated with the occurrence of the An. gambiae complex and An. funestus larvae (Table 5
). The average canopy cover was significantly less in the An. gambiae complexpositive habitats and in the An. funestuspositive habitats than in the negative habitats.
|
|
| DISCUSSION |
|
|
|---|
The results from this study suggest that land cover types affect availability and suitability of aquatic habitats for anopheline larvae. Among the physicochemical variables of aquatic habitats that we examined, we found that canopy cover was the only factor significantly associated with the occurrence of the An. gambiae complex and An. funestus larvae. Unshaded aquatic habitats occurred more often in farmlands and pastures than forests. Therefore, the An. gambiae complex larval habitats were more frequently observed in farmlands and pastures. Larval habitats of the An. gambiae complex also commonly occur in temporary habitats in swamp margins in valley bottoms.18 It is interesting to note that the number of anopheline larval habitats in farmlands and pastures was greatly reduced in the dry season, mostly due to low availability of aquatic habitats. Farmers stop drawing water to farmlands from streams through irrigation ditches during the dry season, and there is little rainfall. Anopheles funestus also prefer habitats in unforested areas, but they tend to occur more often in larger, semipermanent or permanent habitats with emergent plants.30,32,33 The confined distribution of An. funestus larval habitats observed in the present study is probably due to limited distribution of longstanding aquatic habitats.
Frequent occurrence of the An. gambiae complex larvae in temporary and sunlit habitats may result from a combination of several factors. At least three plausible explanations for this phenomenon are suggested. First, the An. gambiae complex females preferentially select small, open habitats for oviposition.35 Second, larval predation is less prevalent in temporary habitats than in large, permanent habitats.19,3638 Third, the An. gambiae complex is a typical r-strategist (a species that is opportunistic and reproduces rapidly when density-independent limiting factors are not present), exploiting the increased resources of warmer, open habitats that tend to produce more algae (the main food source for the An. gambiae complex) than do shaded habitats.39 Since maize is a dominant crop in this study area, maize pollen, which is an excellent food source, may also enhance survivorship of the An. gambiae complex in open habitats in farmlands.40,41 Moreover, warmer temperatures encountered in small and open habitats during daytime hours shorten larvae-to-pupae development time, and subsequently larval mortality due to desiccation is reduced.32 The An. gambiae complex may have evolved to exploit these favorable conditions by selecting small and open habitats for oviposition.
Our results suggest that land use practices affect the growing conditions of anopheline larvae and habitat availability in the African highlands. Patches of forest still remain along streams in the Kenyan highlands such as our study area, and clearing these riparian forests will create more larval habitats for the malaria vectors. Assessed from the perspective of larval ecology, our results support the hypothesis that clearing forests leads to an increased incidence of malaria.15 This study also demonstrated that mapping the distribution of anopheline larval habitats may be used as a foundation for larval control. Our finding on the aggregated distribution of larval habitats in the highland sites suggests that larval control can be targeted to aquatic habitats in farmlands and pastures near streams or valley bottoms. Larval control may be more effective if implemented during the dry season when the distribution of larval habitats is more strictly confined toward valley bottoms.
Received July 24, 2004. Accepted for publication January 5, 2005.
Acknowledgments: We thank H. Atieli, B. Ndenga, M. Okonji, S. Ajuma, W. Miheso, and P. Lutiali for field assistance, U. Kitron and J. Clennon for assistance with statistical analyses, M. Ombok for assistance with the global positioning system in differential mode, L. Carson for critical comments, and T. Wiegand for providing the statistical software package Programita.
Financial support: This work was supported by National Institutes of Health grants R01 AI 50243 and D43 TW01505.
* Address correspondence to Noboru Minakawa, Faculty of Medicine, Saga University, 5-1-1 Nabeshima, Saga 849-8501, Japan. E-mail: minakawa{at}post.saga-med.ac.jp ![]()
Authors addresses: Noboru Minakawa, Faculty of Medicine, Saga University, 5-1-1 Nabeshima, Saga 849-8501, Japan, E-mail: minakawa{at}post.saga-med.ac.jp. Stephen Munga, Francis Atieli, Andrew K. Githeko, Kenya Medical Research Institute, Kisumu, Kenya, E-mails: SMunga{at}kisian.mimcom.net, Fatieli{at}kisian.mimcom.net, and agitheko{at}kisian.mimcom.net. Emmanuel Mushinzimana, Guofa Zhou, and Guiyun Yan, Department of Biological Sciences, State University of New York, Buffalo, NY 14260, E-mails: Emushinzimana{at}kisian.mimcom.net, gzhou2{at}acsu.buffalo.edu, and gyan{at}buffalo.edu.
Reprint requests: Dr. Noboru Minakawa, Faculty of Medicine, Saga University, 5-1-1 Nabeshima, Saga 849-8501, Japan.
| REFERENCES |
|
|
|---|
This article has been cited by other articles:
![]() |
D. Zhong, Y. Afrane, A. Githeko, Z. Yang, L. Cui, D. M. Menge, E. A. Temu, and G. Yan Plasmodium falciparum Genetic Diversity in Western Kenya Highlands Am J Trop Med Hyg, December 1, 2007; 77(6): 1043 - 1050. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Zhou, S. Munga, N. Minakawa, A. K. Githeko, and G. Yan Spatial Relationship between Adult Malaria Vector Abundance and Environmental Factors in Western Kenya Highlands Am J Trop Med Hyg, July 1, 2007; 77(1): 29 - 35. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. MINAKAWA, E. OMUKUNDA, G. ZHOU, A. GITHEKO, and G. YAN MALARIA VECTOR PRODUCTIVITY IN RELATION TO THE HIGHLAND ENVIRONMENT IN KENYA. Am J Trop Med Hyg, September 1, 2006; 75(3): 448 - 453. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. A. AFRANE, G. ZHOU, B. W. LAWSON, A. K. GITHEKO, and G. YAN 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, May 1, 2006; 74(5): 772 - 778. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. MUNGA, N. MINAKAWA, G. ZHOU, E. MUSHINZIMANA, O.-O. J. BARRACK, A. K. GITHEKO, and G. YAN ASSOCIATION BETWEEN LAND COVER AND HABITAT PRODUCTIVITY OF MALARIA VECTORS IN WESTERN KENYAN HIGHLANDS Am J Trop Med Hyg, January 1, 2006; 74(1): 69 - 75. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |