Research | Open | Published:
Spatial pattern of schistosomiasis in Xingzi, Jiangxi Province, China: the effects of environmental factors
Parasites & Vectorsvolume 6, Article number: 214 (2013)
The recent rebounds of schistosomiasis in the middle and lower reaches of the Yangtze River pose a challenge to the current control strategies. In this study, identification of potential high risk snail habitats was proposed, as an alternative sustainable control strategy, in Xingzi County, China. Parasitological data from standardized surveys were available for 36,208 locals (aged between 6–65 years) from 42 sample villages across the county and used in combination with environmental data to investigate the spatial pattern of schistosomiasis risks.
Environmental factors measured at village level were examined as possible risk factors by fitting a logistic regression model to schsitosomiasis risk. The approach of ordinary kriging was then used to predict the prevalence of schistosomiasis over the whole county.
Risk analysis indicated that distance to snail habitat and wetland, rainfall, land surface temperature, hours of daylight, and vegetation are significantly associated with infection and the residual spatial pattern of infection showed no spatial correlation. The predictive map illustrated that high risk regions were located close to Beng Lake, Liaohuachi Lake, and Shixia Lake.
Those significant environmental factors can perfectly explain spatial variation in infection and the high risk snail habitats delineated by the predicted map of schistosomiasis risks will help local decision-makers to develop a more sustainable control strategy.
Schistosomiasis japonica, caused by Schistosoma japonicum, was one of the most serious parasitic diseases in China with a documented history of over 2,100 years . Over the past six decades, China has made great strides toward reducing the prevalence of schistosomiasis japonica, largely through a strategy based on chemotherapy and snail control . Despite this achievement, schistosomiasis has re-emerged in the past decade due to changes of biological, ecological and socio-economic factors and the termination of the World Bank Loan Project (WBLP) on schistosomiasis control in 2001 . A recent review reported that approximately one million people remained infected and more than 50 million people were at risk, especially those living in the marsh and lake (Dongting Lake and Poyang Lake) regions along the Yangtze River basin . The re-emergence of schistosomiasis poses a challenge to the sustainability of current implemented control strategies [4, 5].
It is imperative to adopt an alternative sustainable control strategy in the lake and marshland regions in China. Chemotherapy was the backbone of previous control programs; however, repeated drug therapy resulted in a reduction in chemotherapy compliance . In addition, large scale chemotherapy is no longer possible because the financial support for schistosomiasis control is generally not enough although it has increased since 2009 . The transmission of S. japonicum is closely related to the distribution of its sole intermediate host - Oncomelania hupensis, which controls the distribution of schistosomiasis. However, large-scale control strategies for snail habitats over a whole endemic area are impractical . The identification of areas at high risk of transmission is of great importance in the efficient control schistosomiasis [8, 9].
In this study conducted in Xingzi County, we firstly mapped the spatial distribution of schistosomiasis infection, then analyzed its relevant environmental determinants, and finally identified potential high risk snail habitats.
Xingzi County is located in the northern part of Jiangxi Province, adjacent to the Poyang Lake, the largest freshwater lake in China. It has a total population of 25 million with 16 million living in the areas affected by schistosomiasis. Ideal snail habitats are common in Xingzi County; these are characterized by habitats described as “land in winter, water in summer”. The water level rises from April to October (wet season) to submerge the whole substratum upon which the snails live (making the area unsuitable for snail activity). During the dry season (mainly from November to the following March) water levels in the marshlands gradually subside and the grassy wet-lands favored by the snails re-emerge. This provides an ideal environment for snail growth and reproduction, which leads Xingzi County to be one of the most affected regions for schistosomisis historically.
S. japonicum infection data were obtained from a cross-sectional survey conducted by health professional of the local anti-schistosomiasis station at 42 randomly selected locations in November 2008. The locations of surveyed villages are shown in Figure 1. A total of 36,208 residents aged 6–65 years participated in the survey. Indirect Hemagglutination Assay (IHA) test was first used to screen the individuals according to the manufacturer’s instructions and then the parasitological test was applied to the seropositive individuals by reading three Kato-Katz thick smears from one stool specimen. All those individuals found to be egg-positive were diagnosed with S. japonicum infection. Approval for oral consent and other aspects of this survey was granted by the Ethics Committee of Fudan University (ID: IRB#2011-03-0295).
Rainfall and hours of daylight
Monthly data of rainfall and hours of daylight in 2008 were obtained from China Meteorological Data Sharing Service System (http://cdc.cma.gov.cn/home.do). With available data at 756 meteorological stations in China, Inverse Distance Weighting (IDW) interpolation was used to derive data at locations of sampled villages. Four indexes (the minimum, the maximum, mean, and standard deviation) were calculated for each village for the consideration that together they capture, albeit crudely, the effects of overall climate condition and seasonal variation in local climate.
Land Surface Temperature (LST) and vegetation index
All 8-day global 1 km products for LST in 2008 and monthly global 1 km products for vegetation index (Normalized Different Vegetation Index, NDVI) that covered Xingzi County were downloaded from the Level 1 and Atmosphere Archive and Distribution System (http://ladsweb.nascom.nasa.gov/data/search.html). Similarly, the four indexes (the minimum, the maximum, mean, and standard deviation) were obtained for each village as explained above.
Altitude of sampled villages were obtained from an interpolated Digital Elevation Model (DEM) from the Global Land Information System (GLIS) of the United States Geological Survey (http://egsc.usgs.gov/isb/pubs/factsheets/fs06994.html).
Distance to snail habitats and water body
Locations of snail habitats were recorded using handhold Global Positioning System (GPS). Data on water body, including lake and wetland, were downloaded from Conservation Science Data Sets of World Wild Foundation at http://worldwildlife.org. For each sample village, the shortest Euclidian distance to snail habitats and water body were calculated respectively.
Information on population density was derived from Center for International Earth Science Information Network (CIESIN) at Columbia University (http://sedac.ciesin.columbia.edu/data/collection/gpw-v3). All the data manipulation and calculation were carried out in the ArcGIS 10 (ESRI, Redlands, CA).
Logistic regression models were fitted to schistosomasis data to identify significant environmental covariates. Initially, univariate analyses were conducted and variables with P>0.2 were excluded. Colinearity was investigated among all possible pairs of potential predictor variables and if any pair had a correlation coefficient >|0.2|, the member of the pair that was thought less likely to be biologically important was excluded. With the remaining variables, backwards-stepwise regression was conducted using P >0.1 as the exit criterion and P ≤0.05 as the entry criterion. Finally, the residuals from this model were tested for spatial autocorrelation using variogram , which will be introduced later in model-based geostatistical model. If the residuals show no spatial autocorrelation over the study area, this traditional (non-spatial) model was chosen as the fitting model; otherwise, the following model-based geostatistics model was used instead.
Model-based geostatistical model
Let Y i denotes the number of positive samples out of n i individuals tested at village i. The model assumes that the Y i are conditionally independent binomial variates given an unobserved spatial stochastic process S(•), and that the mean response at i depends on explanatory variables observed at village i and on S(i). The spatial random-effect (S(i)) was modeled as a stationary Gaussian process (signal) with the mean 0, variance σ2 and correlation function ρ(u), where u is the distance between two villages. The exponential correlation function used here is defined as follows:
in which ϕ >0 is a scale parameter with the dimensions of distance, measuring the rate at which the spatial correlation decays over distance. Under this function, the practical range (defined as the minimum distance at which spatial correlation between locations is below 5%) is approximately 3 ϕ. Hence if p(i) is the probability that a randomly selected person at a village i will test positive for schistosomiasis, then
where μ(i) refers to a trend surface model or covariate information is available and ϵ(i) is the error term which is assumed to be mutually independent N(0, τ2). In geostatistics, ϵ(i) refers to a discontinuity at the origin in the variogram in geostatistical practice . The role of S(i) in the model is to capture residual spatial variation after adjusting for covariates, without which equation (2) becomes traditional (non-spatial) logit model.
Finally, the predicted map of prevalence of infection with Schistosoma japonicum was produced and overlaid with locations of snail habitats. The ordinary kriging  with the Langevin-Hastings conditional simulation algorithm  was used to generate samples from the predicted distribution of the complete surface S(i), given the observed values of the response variable Y i at each sampled village location. The Markov chain was run for 11,000 iterations to perform conditional simulation. After an initial burn-in of 1,000 iterations, the other 10,000 iterations were run and sampled every 10th iteration to generate the predicted distribution of the risk surface and the corresponding variance surface. The convergence was judged by looking at the Markov chain traces of the model parameters. The traces differed in the extent to which they showed good mixing, but they all exhibited a reasonable degree of convergence to a stationary distribution.
All the statistical analysis was implemented using R2.14.1 software (R Development Core Team, Vienna, Austria).
Observed prevalence of S. japonicum infection for the 42 villages ranged from 0% to 9.12%, with the median of 1.08%. Figure 2 displays the spatial variation in schistosomiasis risk and shows that most of the villages were located along the shores of the major lakes, that villages closer to the lakes tended to have a higher prevalence, and that there seems a trend of reduction in risk moving from northeast to southwest.
Table 1 reports the results of fitting the data to a logistic regression model. Statistics for the covariate regression coefficients and spatial component of residuals from this model are summarized. Note that both σ2 and ϕ are less than 0.01, indicating that almost no spatial trend exists of the infection data after adjusting for the environmental factors. Model coefficients show that distance to snail habitats and wetland, maximum and minimum rainfall, and mean hours of sunshine were significantly negatively associated with schistosomiasis risk while the maximum LST at daytime, the maximum and minimum LST at night, and mean NDVI were significantly positively associated with schistosomiasis.
The map of predicted prevalence for S. japonicum infection (Figure 2) shows relatively high schistosomiasis prevalence (>1.5%) in the southeastern part of the county surrounded by Beng Lake, Liaohuachi Lake, and Shixia Lake. Figure 3 represents estimates of the variance of the predictions. As shown in Figure 3, a lower level of uncertainty was apparent in locations close to sampled villages and a higher level of uncertainty was apparent in locations distant from sampled villages. An exception to this was the areas around Beng Lake where the prediction uncertainty is high though some sampled villages are nearby Beng Lake. The prediction uncertainty, however, is generally low over the study area. A comparison of Figure 3 with Figure 2 shows that prediction variances were large at locations where the predictions themselves are high. The map of snail habitat overlaid by the predicted S. japonicum infection risk is illustrated in Figure 4. Within areas with high predicted prevalence, a portion of snail habitats (marked with triangle) were densely located between Beng Lake, Liaohuachi Lake, and Shixia Lake.
This study represents an application of GIS and RS techniques and geostatistics to predict S. japonicum infection risks in Xingzi County, one of the most serious endemic regions of schistosomiasis in China, where re-emergences and transmission of japonicum were still ongoing . This study depicted the spatial pattern of S. japonicum infection and identified the important role of environmental factors in explaining this geographical variation. The predicted map of schistosomiasis risk provides an empirical basis for identifying prior areas when implementing schistosomiasis controls in local regions.
The observed risk factors found to be related to S. japonicum infection are already well known and are consistent and interpretable with the epidemiology of infection and known biology of freshwater snails. Malacological studies confirm that the habitats of O. hupensis snail species in lower reaches of the Yangtze River include lake/marshland regions and hill regions, both of which have extensive physical connections with the Yangtze River through channels or in low floodplains beside the Yangtze River [13, 14]. With frequent floodings of the Poyang Lake due to its connection to the Yangtze River, snails in these habitats can be dispersed and subsequently deposited widely in various localities, such as rivers, lakes, and wetland. Hence, individuals living on or near the shore are more likely to undertake risky water contact behavior as a result of agricultural activities and fishing, increasing their exposure to S. japonicum. The distances to snail habitats and wetland could be regarded as the indicators of exposure because people tend to contact these environments more frequently if they live nearby. This finding is consistent with previous reports [8, 15]. Climatic condition, such as LST, rainfall, and hours of daylight, has been shown to influence the distribution and density of snails and the rate of schistosomal development in the snail host [16–18], which is supported by our results. The NDVI is the most important vegetation index and has been widely and successfully used for prediction of intermediate host snails of schistosomiasis . Higher NDVI indicates higher abundance of green vegetation which could be used as a proxy for suitable O. hupensis habitats, thus facilitating transmission.
Re-emergence of schistosomiasis indicates that maintaining and consolidating the achievements obtained needs a further adjustment. It is well known that the intermediate host --- snail, is an important part of the life cycle of S. japonicum, therefore, controlling snails is an efficient way to keep long-term sustainable control effects [20, 21]. To achieve this, it is important to locate potential high risk snail habitats and take corresponding effective control strategies, such as mollusciciding and environmental modification. The predicted map of schistosomiasis risk (Figure 2) shows that regions with a relatively high prevalence of >1.5% are mainly around lakes, indicating snail habitats within these regions (marked with triangle in Figure 4) are the targets for action. There is a relatively high level of uncertainty in regions with high prevalence however. When disease maps are used in planning control programs, it is important to acknowledge such uncertainties when interpreting maps for disease control .
Our study highlights two limitations which deserve further discussions. One is that our prediction model relies on an assumption of stationarity. This assumption is appropriate when mapping disease over small areas, but questionable over wide areas, such as a country . And geostatistical methods for modeling non-stationarity data have not well developed in spite of some studies [24, 25]. However, to investigate the effects of non-stationarity is interesting and challenging and needs to be considered in further studies. The other potential limitation is that there might be diagnostic errors of infection. Generally low infection with S. japonica in recent years results in uncertainty in both sensitivity and specificity of infection . Wang el at  used a Bayesian method deal with this uncertainty in a case study. The assumption of this method is that the diagnostic sensitivity and specificity for every village are all the same, which are impractical. In addition, this method might introduce new errors because of subjective prior knowledge for model parameters which heavily influences the corresponding posterior distributions . Modeling with diagnostic errors is of consideration in future studies.
In summary, this study investigated a region where schistosomiasis endemic remained an important public health problem. The results revealed the spatial patterns of S. japonicum infection in Xingzi County, which help understand its epidemiological characteristics. More importantly, we identified locations of high risk snail habitats where control efforts should be taken as a strategic priority. These results will help in the development of long-term and sustainable strategies for schistosomiasis control.
Zhou XN, Wang LY, Chen MG, Wu XH, Jiang QW, Chen XY, Zheng J, Utzinger J: The public health significance and control of schistosomiasis in China–then and now. Acta Trop. 2005, 96: 97-105. 10.1016/j.actatropica.2005.07.005.
Zhou XN, Guo JG, Wu XH, Jiang QW, Zheng J, Dang H, Wang XH, Xu J, Zhu HQ, Wu GL: Epidemiology of schistosomiasis in the People’s Republic of China, 2004. Emerg Infect Dis. 2007, 13: 1470-1476. 10.3201/eid1310.061423.
McManus DP, Gray DJ, Li Y, Feng Z, Williams GM, Stewart D, Rey-Ladino J, Ross AG: Schistosomiasis in the People’s Republic of China: the Era of the Three Gorges Dam. Clin Microbiol Rev. 2010, 23: 442-466. 10.1128/CMR.00044-09.
Liang S, Yang C, Zhong B, Qiu D: Re-emerging schistosomiasis in hilly and mountainous areas of Sichuan, China. Bull World Health Organ. 2006, 84: 139-144. 10.2471/BLT.05.025031.
Utzinger J, Bergquist R, Shu-Hua X, Singer BH, Tanner M: Sustainable schistosomiasis control–the way forward. Lancet. 2003, 362: 1932-1934. 10.1016/S0140-6736(03)14968-9.
Zhou YB, Yang MX, Wang QZ, Zhao GM, Wei JG, Peng WX, Jiang QW: Field comparison of immunodiagnostic and parasitological techniques for the detection of Schistosomiasis japonica in the People’s Republic of China. AmJTrop Med Hyg. 2007, 76: 1138-1143.
Chen XY: PhD thesis. Schistosomiasis control strategy, administrative pattern, financial support and technological processes under different background in China. 2005, Shanghai: Fudan University, Department of Epidemiology
Zhang ZJ, Carpenter TE, Lynn HS, Chen Y, Bivand R, Clark AB, Hui FM, Peng WX, Zhou YB, Zhao GM: Location of active transmission sites of Schistosoma japonicum in lake and marshland regions in China. Parasitology. 2009, 136: 737-746. 10.1017/S0031182009005885.
Utzinger J, Mayombana C, Mez K, Tanner M: Evaluation of chemical and physical-morphological factors as potential determinants of Biomphalaria pfeifferi (Krauss, 1848) distribution. Mem Inst Oswaldo Cruz. 1997, 92: 323-328.
Cressie N: Statistics for spatial data. 1993, New York: John Wiley & Sons, Inc
Isaaks EH, Srivastava RM: An Introduction to Applied Geostatistics. 1989, New York: Oxford University Press
Papaspiliopoulus O, Roberts GO, Skold M: Non-centered parameterizations for hierarchical models and data augmentation. Bayesian statistics 7. Edited by: Bernardo M, Bayarri S, Berger JO, Dawid AP, Heckerman D, Smith AFM, West M. 2003, USA: Oxford University Press, 307-326.
Zhao QP, Jiang MS, Littlewood DT, Nie P: Distinct genetic diversity of Oncomelania hupensis, intermediate host of Schistosoma japonicum in mainland China as revealed by ITS sequences. PLoS Negl Trop Dis. 2010, 4: e611-10.1371/journal.pntd.0000611.
Garcia RG: Tolerance of Oncomelania hupensis quadrasi to varying concentrations of dissolved oxygen and organic pollution. Bull World Health Organ. 1972, 47: 59-70.
Clements AC, Moyeed R, Brooker S: Bayesian geostatistical prediction of the intensity of infection with Schistosoma mansoni in East Africa. Parasitology. 2006, 133: 711-719. 10.1017/S0031182006001181.
Appleton C: Review of literature on abiotic factors influencing the distribution and life-cycles of Bilharziasis intermediate host snails. Malacol Rev. 1978, 11: 1-25.
Sturrock RF: The intermediate hosts and host- parasite relationships. Human schistosomiasis. Edited by: Jordan P, Webbe G, Sturrock R. 1993, Wallingford: CAB International, 33-85.
Zheng Y, Qang Q, Zhao G, Zhong J, Zhang S: The function of the overlaying climate data in analysis of Oncomelania snail distribution. Chin. Public Health. 1998, 14: 724-725.
Guo JG, Vounatsou P, Cao CL, Utzinger J, Zhu HQ, Anderegg D, Zhu R, He ZY, Li D, Hu F: A geographic information and remote sensing based model for prediction of Oncomelania hupensis habitats in the Poyang Lake area, China. Acta Trop. 2005, 96: 213-222. 10.1016/j.actatropica.2005.07.029.
Singer BH, de Castro MC: Bridges to sustainable tropical health. Proc Natl Acad Sci. 2007, 104: 16038-16043. 10.1073/pnas.0700900104.
Zhang Z, Carpenter TE, Chen Y, Clark AB, Lynn HS, Peng W, Zhou Y, Zhao G, Jiang Q: Identifying high-risk regions for schistosomiasis in Guichi, China: a spatial analysis. Acta Trop. 2008, 107: 217-223. 10.1016/j.actatropica.2008.04.027.
Clements AC, Firth S, Dembele R, Garba A: Use of Bayesian geostatistical prediction to estimate local variations in Schistosoma haematobium infection in West Africa. B World Health Organ. 2009, 87: 921-929. 10.2471/BLT.08.058933.
Gosoniu L, Vounatsou P, Sogoba N, Smith T: Bayesian modeling of geostatistical malaria risk data. Geospat Health. 2006, 1: 127-139.
Sampson PD, Gottorp P: Nonparametric estimation of nonstationary spatial covariance structure. J Am Stat Assoc. 1992, 87: 108-119. 10.1080/01621459.1992.10475181.
Damian D, Sampson PD, Guttorp P: Bayesian estimation of semi-parametric nonstationary spatial covariance structures. Environmetrics. 2001, 12: 161-178. 10.1002/1099-095X(200103)12:2<161::AID-ENV452>3.0.CO;2-G.
Zhou YB, Zheng HM, Jiang QW: A diagnostic challenge for Schistosomiasis japonica in China: consequences on praziquantel-based morbidity control. Parasit Vectors. 2011, 4: 194-10.1186/1756-3305-4-194.
Wang XH, Zhou XN, Vounatsou P, Chen Z, Utzinger J, Yang K, Steinmann P, Wu XH: Bayesian spatio-temporal modeling of Schistosoma japonicum prevalence data in the absence of a diagnostic ‘gold’ standard. PLoS Negl Trop Dis. 2008, 2: e250-10.1371/journal.pntd.0000250.
Gelman A: Objections to Bayesian statistics. Bayesian Anal. 2008, 3: 445-450.
This work was financially supported by the National Natural Science Foundation of China (grant numbers are 81102167, 81172609, and J1210041), Specialized Research Fund for the Doctoral Program of Higher Education, SRFDP (grant number is 20110071120040), a Foundation for the Author of National Excellent Doctoral Dissertation of PR China (FANEDD) (201186), the National S&T Major Program (2012ZX10004-220, 2008ZX10004-011), and Ecological Environment and Humanities/Social Sciences Interdisciplinary Research Project of Tyndall Center of Fudan University (FTC98503A09). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
All authors declare that they have no competing interests.
ZJ and QJ conceived the study; ZL, JG, BT, and QJ performed the field collections; YH, ZJ, and YC wrote the manuscript; YH performed statistical analyses. All the authors read and approved the final version of the manuscript.