|
|
||||||||
| ABSTRACT |
|
|
|---|
| INTRODUCTION |
|
|
|---|
Substantial spatial and temporal variation in transmission of schistosomiasis occurs within human populations.4,5 Human infection patterns vary as a function of water contact patterns, immunity, the presence of competent intermediate snail hosts,6 and the availability of suitable aquatic habitats for water use. Along the southern coast of Kenya, human exposure to S. haematobium occurs at snail habitats (ponds and streams) infested with the intermediate host snail Bulinus nasutus. Complicating the effective targeting of control measures in coastal Kenya is the use by residents of a variety of different ponds for bathing, swimming, and washing laundry.7 As Chandiwana and Woolhouse8 and others911 have shown, the frequency and intensity with which people use a diversity of contaminated and uncontaminated water sources affects transmission patterns. Because water use of different ponds is not independent, the association of the human population with various ponds in the area needs to be considered when analyzing spatial patterns of infection.
Clennon and others12 described the spatial clustering of S. haematobium infection around one infested pond in a single rural village in Msambweni. To better understand the influence of multiple transmission sources on infection patterns, retrospective and current spatial patterns were examined at the household level by integrating human infection and snail habitat locational data, and applying spatial analyses in 10 villages throughout Msambweni Division, Kwale District, Kenya.
| METHODS |
|
|
|---|
|
Nine water bodies were monitored for human water contact, snail distribution, and environmental conditions. Monitored water sources included six ponds (Kiziamkala Pond, Maridzani Pond, Bovo Pond, Nimbodze Pond, Mwachiangwa Pond, and Mwamagongo Pond), a rice field, the Mukurumudzi River, and the Lukungwi Stream (Figure 1
). When the ponds contain water, a variety of snail species can be collected, including B. nasutus, the intermediate host of S. haematobium in the study area.16 During the 1980s, B. nasutus snails were found at every water source at least once, with the river being the least favorable habitat. Most snails found shedding S. haematobium cercariae were from Kiziamkala and Maridzani Dams, followed by Bovo, Nimbodze and Marigiza ponds. The percentages of B. nasutus found shedding cercariae were generally less than 2%, even at Kiziamkala and Maridzani Dams. In 2001, Nimbodze Pond had the greatest numbers of B. nasutus (approximately 3%) shedding S. haematobium cercariae, followed by Maridzani, Kiziamkala and Mwamagongo ponds. Bovo and Mwamagongo ponds were also found to have B. nasutus shedding, but in lower numbers. Along the stream, B. nasutus snails were found, but none shed.
Ethical oversight. Prior to collecting urine samples for determining S. haematobium infection status, informed consent was obtained from area residents (or for children their parents). These studies were performed under human investigations protocols approved by the Ethical Review Board of Kenya Medical Research Institute (Nairobi, Kenya) and by the Human Investigations Review Board of University Hospitals (Cleveland, OH).
Infection prevalence. Residents submitted two 10-mL urine specimens at noon a week apart for examination for S. haematobium infection during 20002002. Infections were detected using Nuclepore filtration, and measured up to 1,000 eggs per urine specimen.17 Aggregated household infection levels are reported as geometric mean density (number of parasites per person in a household regardless of infection status), which is a function of both prevalence (the proportion of people infected in a household) and geometric mean intensity (the geometric mean number of parasites for only infected people infected in a household).12,18 Schistosoma haematobium infection patterns in the population of Milalani Village were originally described by Clennon and others.12 These data were compared with those collected in the 1980s by King and others.19
Statistical analysis. Our study population comprised children 617 years of age, the age group that exhibits the highest infection levels, and who unlike older residents were not previously enrolled in a school-based schistosome chemotherapy program.19 Differences in prevalence between age groups were tested using the homogeneity chi-square test. Correlation was used to assess the correspondence of eggs counts (log10 [n + 1]) between urine tests. Logistic regression (forward conditional) was used to appraise the relationship between whether children were infected with sex, age (continuous from 6 to 17 years), age-by-sex interaction, the nearest household distance to an alternative water source (open well, borehole, piped-water), known school attendance, total household population, and number of children per household. A follow-up univariate logistic regression using age as a factor (69, 1013 and 1417 years) was used to calculate odds ratios for each age group. Stepwise linear regression was then used to evaluate the relationship of individual infection levels (log10 [average eggs + 1]) with those variables. These non-spatial statistical analyses were performed using SPSS version 11.5.0 software (SPSS Inc., Chicago, IL).
Geospatial processing. Households and water sources were mapped as described by Clennon and others using a high-resolution satellite image comprised of 1-m2 panchromatic and 4-m2 multi-spectral images of Msambweni.12 During the 2001 house-to-house demographic survey, each household was given a village affiliated household identification (HID) number that was marked on the exterior of the home, and each person was assigned a individual identification (ID) number incorporating the HID. Demographic information collected on each individual included name (common and tribal), sex, year of birth, mothers name and ID, mothers mother name, fathers name and ID, and fathers father name. When urine samples where submitted, the participants name, ID (or at least HID), sex, and year of birth were again recorded so that the infection data could be correctly linked with demographic data. When household locations were mapped, homes were identified based on roof type (thatch, tin, rusty-tin), household relationship to other homes and landmarks (e.g., roads, walking paths, trees by type water sources) in the pan-sharpened Ikonos image. The HID were then marked directly upon the Ikonos image in the field, and digitized in the GIS later the same day. When there was any uncertainty, a GPS reading was taken and a small sketch map of the location was drawn.
To identify 1984 household locations, village workers familiar with the residents in the area went to the participants (identifiable by name, date of birth, sex, and prior school attendance and past participation) current homes and inquired if the home where they currently resided was their home in 1984 or if they have since moved. In cases when the participant had moved, they were asked about the previous homes location (if they knew who was currently residing there or if they could show them the location). We were able to locate 1,014 homes from the 1980s.
Household and water source locations were joined with demographic, parasitologic, malacologic, and environmental data. Locational data were georectified to the Universal Transverse Mercator (UTM) zone 37S projection, 1984 datum, in the GIS software package ArcGIS 9.0 (Environmental Systems Research Institute, Redlands, CA).
Spatial statistics. For spatial analyses, children were grouped by age (69, 1013, and 1417 years). Global, local, focal, and directional spatial analyses were used to examine the spatial structure of S. haematobium infection patterns and identify significant clustering of elevated infection levels in the study area. Global Ripleys K-function,20 global weighted K-function,21 local Gi*(d), and focal Gi(d)2225 were applied to human infection patterns using Point Pattern Analysis (Chen and Getis24) and ClusterSeer (TerraSeer, Ann Arbor, MI) software.
Global second-order (K-function and weighted K-function) spatial analyses that compare the observed pattern of points with a homogenous Poisson process were used to classify spatial patterns over the entire study area as being random, clustered, or uniformly dispersed (tending towards regular), with significance determined using Monte Carlo simulations. The observed point pattern in a K-function analysis is composed of a 0/1 distance matrix that includes points according to distance, over a range of distances. The weighted K-function considers values in each location (e.g., household density of infection), and permutations that randomize these values among household locations are used to determine significance, thus rendering it independent of site locations. We used global K-function to assess the general distribution of household locations, and weighted K-function to evaluate the overall submission distribution and spatial pattern of infection.
The local statistic Gi*(d) was applied to describe local variation in the spatial structure of infection within the study area. Local Gi(d) and Gi*(d) spatial statistics detect significant autocorrelation or clustering around each point by comparing point values (e.g., household infection density) to the overall mean. Mathematically, the values at each point are weighed according to distance (1 = within the considered distance; 0 = outside the considered distance) using a circular extent (Figure 2A
). The Gi(d) and Gi*(d) statistics yield normal standard variables. To account for multiple comparisons, significance levels for local Gi*(d) analyses were determined according to Ord and Getis.23 Clennon and others12 and others26,27 have demonstrated that the local Gi(d) spatial statistic can be applied as a focal statistic by only evaluating the spatial association of point values (infection density, infestation levels) at different distances surrounding unique locations of interest.
|
The 2000 data were compared with data from a 1984 school survey from 1,647 households that could be located19 and for which spatial statistics could be applied to school age children of three age groups (69, 1013, and 1417 years).
| RESULTS |
|
|
|---|
The logistic regression of infection of school age (617 years) children adequately fit the data (Hosmer and Leme-show
2 = 9.7, P > 2.29), but explained only 12% of the variance (Nagelkerke pseudo-R2). Age (as a covariate) was the only significant variable (Wald = 90.6, P < 0.001), and had positive association with infection (ß = 0.21). When age groups (as a factor) were further examined using univariate logistic regression, both adolescents (1417 years of age) (odds ratio [OR] = 2.55, 95% confidence interval [CI] = 2.153.02) and 1013 years-old persons (OR = 1.69, 95% CI = 1.692.32) had higher odds of being infected than 69 year-old children. Although linear regression explained less than 10% of the variance (P < 0.0001) in infection levels, it detected a significant positive relationship between age and infection level (ß = 0.36, t = 9.13, P < 0.0001). Other factors such as sex (t = 0.74, P > 0.46), distance to nearest alternative water source (t = 0.23, P > 0.82), and school attendance (t = 0.85, P > 0.40) did not have significant relationships with degree of infection.
Spatial infection patterns. Household locations were highly aggregated (K-function, P < 0.01) at distances up to 2,500 meters, and submission of urine samples by children was evenly distributed. Household density of infection in each age group was not aggregated (weighted K-function) during both study periods.
During 2000, significant local spatial autocorrelation of high infection levels at 500 meters occurred near Nimbodze Pond for all age groups (Figure 3A
). Additional clustering at 500 meters also occurred near Maridzani Dam among older children (1017 years of age) and Mwachiangwa Dam only for adolescents (1417 years of age). At distances greater than 500 meters, a significant positive autocorrelation was detected throughout the area between Nimbodze, Kiziamkala, and Mwachiangwa ponds, with additional clustering among young children (69 years of age) occurring south of Lukungwi Stream near Nimbodze Pond.
|
Focal clustering of infection.
Focal clustering of infection density around the different water sources varied by age and between the two study periods. In 2000, high levels of infection in children 69 years of age were clustered starting and peaking in close proximity to Nimbodze Pond with clustering persisting past 1,500 meters (Figure 4A
). High infection levels were also clustered among older children (1013 and 1417 years of age) around Nimbodze Pond, but the degrees of clustering were much lower initially and peaked farther from Nimbodze Pond (Figure 4A
). During 2000, the greatest degree of clustering around Maridzani Dam was detected among 1013 year-old persons, with clustering highest close to the dam and significant levels up to a distance of 650 meters (Figure 4B
), as was the clustering of infections in adolescents (1417 years of age). In contrast, clustering among 69 year-old children was significant only between 350 and 450 meters (Figure 4B
). Near the river in 2000, significant clustering of low levels of infection was detected for 1013 and 1417 year-old persons (Figure 4C
).
|
When patterns of recent infection density are compared with data from 1984, differences in the patterns of clustering are apparent. In 1984, clustering of high infection density among all age groups was focused to a greater degree and more extensively around Maridzani Dam, rather than around Nimbodze Pond (Figures 4D and 4E
). In 2000, significant clustering of high infection levels occurred more extensively around Nimbodze Pond (up to a distance of 1,000 meters) than around Maridzani Dam (Figures 4A and 4B
).
Directional clustering near transmission foci.
When directionality was considered in the focal cluster analyses, most of the clustering of high household infection density occurred from east of Nimbodze Pond to west of Maridzani Dam. During 2000, all age groups showed significant clustering extending east from Nimbodze Pond, peaking at 750 meters (Figure 5A
), decreasing until 1,2501,500 meters (the halfway mark between Nimbodze Pond and Maridzani Dam), and then increasing until 2,2502,500 m (the distance to Maridzani Dam). Although the general trend of Gi(d) values in 1984 was similar (significant clustering early with a subsequent decrease in levels that bottom out at approximately 1,250 meters before they increase again) (Figure 5B
). Only at distances of 500750 meters east of Nimbodze Pond was infection significantly clustered for children 69 and 1417 years of age.
|
Some clustering of infection density was detected in other directions around Maridzani Dam. Clustering of high infection levels extended north 500 meters among 10 to 13 year-old persons in both 2000 and 1984. Additional clustering to the east of Maridzani Dam was detected in 1984 through to 2,500 meters in older children (1013 and 1417 years of age) but not in children 69 years of age.
| DISCUSSION |
|
|
|---|
Our findings suggest that infection levels of human urinary schistosomiasis are clustered as a function of the spatial distribution of water sources that are contaminated, as well as those that are not contaminated (e.g., the river) with S. haematobium-infected snails. The effects of main transmission sources on household infection density were significant through 1,500 meters, which is well within the known distance range that people will travel to ponds in the area. Significant directional patterns of clustering can be attributed to most clusters being located between two transmission foci (Nimbodze Pond and Maridzani Dam). The sharp gradient in levels of focal clustering (Figure 4B
) was greatly reduced in the directional focal analysis (Figure 5C
), which demonstrates the strong anisotropy of the data. The clusters of low infection levels are likely associated with the use of uncontaminated water sources that allow people to avoid infection (e.g., rivers and streams).
By considering two study periods 16 years apart and childrens ages, a temporal shift in the primary source of S. haematobium transmission was identified. The spatial patterns of infection in 1984 can be attributed to more intense transmission that was occurring around Maridzani and Bovo Ponds. In 2000, clustering was found primarily around Nimbodze Pond. The change in clustering pattern between 1984 and 2000 suggests a shift in the primary source of transmission in the area, which may be associated with changes in environmental conditions affecting intermediate host snails in the ponds, or with shifts in patterns in human use of ponds. Some spatial heterogeneity in infection patterns during 1984 may not have been recognized because households (linked with infection data) could not be mapped as extensively as the current data.
Although Kiziamkala Dam was found to have many B. nasutus shedding S. haematobium cercariae during the 1980s, people living east of the dam would frequent the river where no B. nasutus have been found shedding cercariae. The clustering of intense infections between Bovo Pond and Maridzani Dam, and not by Kiziamkala Dam, suggests that the use of a non-cercarial infested water source (the river) mitigated infections levels. Only limited 2000 era water contact data was collected in 2001 because of drought conditions and subsequent drying of the primary water bodies used by residents.
Differences in clustering among children of different age groups reflect the effects of exposure to waters contaminated with S. haematobium cercariae, as well as of acquired immunity. The most dramatic shift in local and focal clustering patterns was seen among children 69 years of age, those with limited exposure to Maridzani and Bovo ponds before the contribution of these two ponds to transmission was apparently reduced. Variations in exposure to S. haematobium-contaminated water sources among the different age groups are likely due to differences in water-use activities and in walking distances to water sources.
Maridzani Dam (with interaction with Bovo Pond) appears to have been the primary S. haematobium transmission epicenter in the 1980s, and Nimbodze Pond acted as a main focus in 2000. Our results also suggest that older children in 2000 are carrying infections acquired when Maridzani was still the prominent transmission focus. It is possible that with an increasing number of households built in closer proximity to Nimbodze Pond, the epicenter for S. haematobium transmission shifted from Maridzani Dam to Nimbodze Pond. However, changes in snail distribution patterns suggest that the shift in parasite transmission and the resulting infection clustering pattern are more likely associated with climatic and environmental changes.4,16
During 2000, more snails and more snails shedding cercariae were found at the Nimbodze Pond water contact sites compared with Maridzani Dam, and the reverse was observed in 1984. This may be related to an El Niño event during 19971998. The flooding that occurred in 19971998 was devastating in many parts of Kenya including Kwale District, where the Msambweni study area is located, and was associated with significant environmental degradation (e.g., soil erosion), as well as more small temporary ponds. It is possible that the snail habitats at Maridzani Dam were more negatively affected by soil erosion. During El Niño years, much of the area surrounding Nimbodze Pond, which is more level than the area around Maridzani Dam, was flooded; this likely increased the amount of snail habitat in that part of the study area. It is also possible that the intensity of S. haematobium transmission was decreased at Maridzani Dam because people living in the vicinity of the pond began using small temporary ponds created by the 19971998 El Niño flooding. In addition, the conversion of sugar cane fields to rice fields (which requires an aquatic environment) west of Nimbodze during the 1990s increased the connectivity of the B. nasutus source and sink habitats of Nimbodze Pond, leading to a more stable intermediate host population, which subsequently led to increased S. haematobium transmission.
Applying spatial statistics, Clennon and others12 determined that S. haematobium infection was clustered around a single pond (Nimbodze) known to be contaminated with S. haematobium cercariae. Age-related water contact behaviors and acquired immunity were the most probable determinants underlying the observed pattern of clustering of high density of infection. By moving upscale and expanding our study area to include neighboring villages and additional ponds in the study presented here, we were able to account for the area-wide nature of household use of water sources. By applying the same spatial analysis and a directional focus statistic to data from an earlier study period, we could consider changes in transmission patterns over time. Because of sharp discontinuity in transmission habitats (and infection patterns) and that the process driving infection in Msambweni is not intrinsically stable, a risk map was not constructed.
Based on the complex interactions between study period, age, and spatial associations with water sources, we conclude that a range of spatial and temporal scales needs to be considered for an understanding of S. haematobium transmission patterns in an area. As environmental and social conditions change, spatial analyses allow control programs to better to focus interventions targeting schistosomiasis transmission in space and time.
Received February 8, 2006. Accepted for publication August 21, 2006.
Acknowledgments: We thank the residents of Nganja, Milalani, Vidungeni, Maragiza, Mabatani, Sawa Sawa, Bomani, Mwaembe, Kisimachande, Vingujini and Vidungeni villages in Msambweni Division, Kwale District, Kenya for their involvement in this study. We also thank Idi Masemo, Anthony Chome, Joyce Bongo, Robinson Ireri, and Jackson Muinde for helping locate households, and Anna Schot-thoefer and Gonzalo Vazquez-Prokopec for helpful discussions. The directional Gi(d) was developed in a collaborative effort with Gonzalo Vazquez-Prokopec. Publication of this paper is done with the kind permission of the Director of Medical Services, Ministry of Health, Kenya. This research was conducted in partial fulfillment of Ph.D. degree requirements by Julie A. Clennon.
Financial support: This research was supported by the National Institutes of Health/National Science Foundation Ecology of Infectious Disease program grant (AI45473) and by the Fogarty International Center grant (TW/ES01543) to Charles H. King. Additional support for this study was provided by TerraSeer Inc., (Ann Arbor, MI) through a graduate student research award to Julie A. Clennon.
* Address correspondence to Julie A. Clennon, Department of Molecular Microbiology and Immunology, Johns Hopkins Bloomberg School of Public Health, 615 North Wolfe Street, Baltimore, MD 21205. E-mail: jclennon{at}jhsph.edu ![]()
Authors addresses: Julie A. Clennon, Department of Molecular Microbiology and Immunology, Johns Hopkins Bloomberg School of Public Health, 615 North Wolfe Street, Baltimore, MD 21205, E-mail: jclennon{at}jhsph.edu. Peter L. Mungai, c/o Case Western Reserve University/Division of Vector-Borne Diseases/Kenya Medical Research Institute Filariasis-Schistosomiasis Research Unit, PO Box 8, Msambweni, Kenya, Telephone: 254-40-52267, E-mail: dvbdcwru{at}wananchi.com. Eric M. Muchiri, Division of Vector Borne Diseases, Ministry of Health, PO Box 20750, Nairobi, Kenya, Telephone: 254-20-725833, Fax: 254-20-720030, E-mail: schisto{at}wananchi.com. Charles H. King, Center for Global Health and Diseases, W137, Case Western Reserve University School of Medicine, 10900 Euclid Avenue, Cleveland, OH 44106-4983, Telephone: 216-368-4818, Fax: 216-368-4825, E-mail: chk{at}po.cwru.edu. Uriel Kitron, Division of Epidemiology and Preventive Medicine, Department of Pathobiology, University of Illinois, 2001 South Lincoln Avenue, Urbana, IL 61802, Telephone: 217-265-0714, Fax: 217-244-7421, E-mail: ukitrom{at}uiuc.edu.
Reprint requests: Uriel Kitron, Division of Epidemiology and Preventive Medicine, Department of Pathobiology, University of Illinois, 2001 South Lincoln Avenue, Urbana, IL 61802.
| REFERENCES |
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |