Next Article in Journal
Gender Division of Labor, Burnout, and Intention to Leave Work Among Young Female Nurses in Japan: A Cross-Sectional Study
Previous Article in Journal
Efficiency of Health Care Risk Waste Management in Rural Healthcare Facilities of South Africa: An Assessment of Selected Facilities in Vhembe District, Limpopo Province
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatiotemporal Heterogeneity in Human Schistosoma japonicum Infection at Village Level in Hubei Province, China

1
Key Laboratory of Monitoring and Estimate for Environment and Disaster of Hubei Province, Institute of Geodesy and Geophysics, Chinese Academy of Sciences, Wuhan 430077, China
2
Hubei Center for Disease Control and Prevention, Hubei Provincial Academy of Preventive Medicine, Wuhan 430079, China
3
College of Earth and Planetary Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
4
State Key Laboratory of Resources and Environmental Information Systems, Institute of Geographical Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China
*
Authors to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2019, 16(12), 2198; https://doi.org/10.3390/ijerph16122198
Submission received: 29 May 2019 / Revised: 19 June 2019 / Accepted: 20 June 2019 / Published: 21 June 2019

Abstract

:
The spatiotemporal dynamics of Schistosoma japonicum, combined with temporal heterogeneity among regions of different epidemic areal-types from a microscale viewpoint might capture the local change dynamics and thus aid in optimizing the combinations of precise schistosomiasis control measures. The prevalence data on schistosomiasis infection from 2007 to 2012 in the 30 most endemic counties of Hubei Province, Central China, were appended to the village-level administrative division polygon layer. Anselin local Moran’s I, a retrospective space–time scan statistic and a multilevel-growth model analysis framework, was used to investigate the spatiotemporal pattern of schistosomiasis resident infection rate (RIR) at the village level and how natural geographical environment influence the schistosomiasis RIR over time. Two spatiotemporal high-risk clusters and continuous high-rate clusters were identified mainly in the embankment region across flooding areas of lakes connected with the Yangze and Hanjiang Rivers. Moreover, 12 other clusters and outlier evolution modes were detected to be scattered across the continuous high-rate clusters. Villages located in embankment region had the highest initial values and most rapidly reduced RIRs over time, followed by villages located in marshland-and-lake regions and finally by villages located in hilly region. Moreover, initial RIR values and rates of change did significantly vary (p < 0.001 and p < 0.001, respectively) irrespective of their epidemic areal-type. These local spatiotemporal heterogeneities could contribute to the formulation of distinct control strategies based on local transmission dynamics and be applied in other endemic areas of schistosomiasis.

1. Introduction

Schistosomiasis japonica is caused by Schistosoma japonicum, one of the five trematode species [1] responsible for human and animal infections. It has considerable public health and economic significance in the People’s Republic of China [2,3]. Knowledge of the spatiotemporal heterogeneity of Schistosomiasis japonica on a small scale can provide important insights into local change dynamics and can thus aid in planning schistosomiasis control and elimination efforts. With the development of national and local surveillance systems, long-term temporal and micro-geographic level parasitological data are available which enable the characterization of changing patterns of disease risk in both temporal and spatial dimensions [4]. Different spatial analytical approaches can be used, including global Moran’s I and local Moran’s I statistics [5,6,7,8,9], spatial or space–time scan statistics [6,7,8,10,11,12], and the standard deviational ellipse [9,11]. Almost all studies on spatiotemporal patterns of schistosomiasis japonica focus on the county level because of the lack of village and township based vector data, and few studies have focused on the micro-scale level, especially the village-level [6]. The precision of spatial analyses is usually determined by the resolution of the underlying data (e.g., pixel size, time interval), and data at a coarse spatial resolution would fail to capture local heterogeneity [13].
The Schistosomiasis japonica endemic areas in China are divided into three types, based on the geographical characteristics of schistosomiasis endemic areas and the habitat environment of Oncomelania hupensis, the sole intermediate host of S. japonicum: (1) plain regions with waterway networks, (2) hilly and mountainous regions, and (3) marshland-and-lake regions and many subtypes such as embankment regions, hilly regions [14]. The comprehensive control strategy and prioritization of interventions should be varied between different endemic areas [15,16,17]. However, little is known about the temporal heterogeneity of Schistosomiasis japonica in different epidemic areal-types.
This study aimed to investigate the spatiotemporal patterns of human S. japonicum infection in most endemic counties in the Hubei Province of China at the village level. Anselin local Moran’s I statistics and a retrospective space–time scan statistic were used to determine the spatial distribution of S. japonicum infection over time. A multilevel-growth model was further used to explore the relationship between S. japonicum and epidemic areal-types and time factors.

2. Materials and Methods

2.1. Study Area

Hubei Province in Central China and the middle reaches of the Yangtze River (Figure 1) have included some of the most serious schistosomiasis outbreak regions of all time [12], and they had the highest estimated incidence of schistosomiasis in China up to 2014 [18]. Major rivers include the Yangtze River in the south and the Han River in the north. Though plains dominate half of Jianghan Province, as well as Southcentral Hubei, a series of mountains and hills dominate Western, Northeastern, and Southeastern Hubei. This study was carried out in the 30 most endemic counties with epidemic areal-types, including embankment, marshland-and-lake, and hilly regions in Hubei Province (Figure 1) [18,19].

2.2. Schistosomiasis Information and Geospatial Processing

The prevalence data on schistosomiasis infections in the 30 most endemic counties were collected from repeated cross-sectional surveys conducted annually between 2007 and 2012 by health professionals from the Hubei Institute of Schistosomiasis Control. This study targeted residents of each village who were 6–65 years old, and over 90% of these residents were screened every September to November using the indirect hemagglutination assay. Stool samples were subsequently collected from over 90% of the individuals with positive serological results, and used to perform the miracidium-hatching test [12]. Resident infection rates (RIRs) were calculated by multiplying the positive rate of serological tests with the positive rate of fecal tests. A total of 4682, 4689, 4710, 4722, 4714, and 4945 epidemic villages were surveyed from 2007 to 2012 in the 30 counties, respectively, and the coverage rate of epidemic villages is about 37%; all of them were appended to the village-level administrative division polygon layer based on the names of county, town, and village. An actual geospatial position in the administrative village layer was used to represent the “village” location through field investigation or by consulting local anti-schistosomiasis employers if a name matching the village-level administrative division was unavailable, for example because of a changed name, or because of not being an administrative village unit.

2.3. Spatial Statistical Approach

We performed cluster and outlier analyses (Anselin local Moran’s I) and a retrospective space–time scan statistic to investigate the spatiotemporal pattern of schistosomiasis RIR at village level.
Anselin local Moran’s I was used to identify spatial clusters and spatial outliers of regions with high or low RIR values by calculating a local Moran’s I value, a z-score, a p-value, and a code representing the cluster/outlier type for each region [20]. The z-scores and p-values represent the statistical significance of the computed index values using a Monte Carlo permutation approach [21,22]. Five cluster/outlier types for local autocorrelation exist: (1) HH (high cluster) indicates that a village and its neighbors had high RIR; (2) LL (low cluster) indicates that a village and its neighbors had low RIR; (3) HL indicates that a village was a high outlier among villages with low RIR; (4) LH indicates that a village was a low outlier among villages with high RIR; and (5) ‘not significant’ indicates no spatial autocorrelation. To determine how spatial autocorrelation may evolve with time and location, cluster and outlier analysis results from 2007 to 2012 were joined by a ‘-’, and then all combined types were classified as “HH,” “HL,”“LH,” “LL,” “H to L,” “L to H,” and “Others” (detailed classification in Supplemental Table S1). For instance, “HH-HH-HH-HH-HH-HL” means HH type from 2007 to 2011 and HL type in 2012 and is classified as “H to L.”
The space–time scan statistic using a moving cylindrical window with circular (or elliptic) geographic base and with height corresponding to time is widely used in disease studies to identify statistically significant high-risk spatiotemporal clusters [8,11,23,24]. SaTScan™ (version 9.4, M. Kulldorff, Harvard Medical School, Boston and Information Management Services Inc., Silver Spring, MD, USA.) was used to identify high-risk clusters defined by circular shapes, constrained to clusters with both RIR values and study periods less than 50% [10]. The RIRs showed skewed distributions and were thus normalized by log (base 10) transformation after adding the mean RIR limit to allow for zeros. The spatiotemporal distribution of log transformed RIR was assessed using the normal model designed for continuous data [25]. Monte Carlo simulations (999 times) were used for significance testing at the 0.05 level. The attributed values for these output circles included X and Y coordinates for the centers and cluster radius.
All spatial processing and mapping were carried out in ArcGIS 10.1 (ESRI Inc., Redlands, CA, USA).

2.4. Multilevel-Growth Model

To evaluate the impact of the natural geographical environment (epidemic areal-type) on differences in schistosomiasis RIR over time, a multilevel-growth model analysis framework was developed with “occasions of measures” assigned to level 1 and “village data” to level 2 [26,27,28]. We first coded the predictors to aid substantive interpretations as follows: time (rectilinear: Year = 0, 1…, 6 representing 2007, 2008, …, 2012, respectively) and epidemic areal-type (0—Embankment region, 1—Marshland-and-lake region, 2—Hilly region). A taxonomy of multilevel-growth models was specified following methods described by Singer and Willett (2003) [28].
The unconditional means model, also referred to as an intercept-only model or a null model [29], is expressed as:
RIR i j = γ 00 + ( 𝒰 0 j + e i j )
where RIR i j represent the village j’s schistosomiasis RIR on occasion i, 𝒰 0 j and e i j represent the village- and occasion-specific random effects, respectively. These random effects were assumed to have a normal distribution.
The unconditional growth model, including predictor Year i j into the unconditional means model, is expressed as:
RIR i j = γ 00 + γ 10 Year i j + ( 𝒰 0 j + 𝒰 1 j Year i j + e i j )
Adding a village-level covariate, Areal _ type j , indicating the schistosomiasis epidemic areal-type of village j, leads to:
RIR i j = γ 00 + γ 10 Year i j + γ 01 Areal _ type j + ( 𝒰 0 j + 𝒰 1 j year i j + e i j )
Cross-level interactions can be assessed as:
RIR i j = γ 00 + γ 10 year i j + γ 01 Areal _ type j + γ 11 year i j * Areal _ type i j + ( 𝒰 0 j + 𝒰 1 j year i j + e i j )
where 𝒰 0 j and 𝒰 1 j are bivariate normally distributed with a mean of 0 and constant variance.
All analyses were conducted using SAS (PROC MIXED) version 9.3 (SAS Institute, Inc., Cary, NC, USA).

3. Results

3.1. Spatiotemporal Pattern Analysis

Figure 2 depicts the heterogeneous geospatial distribution of RIR clusters and outliers identified by local Moran’s I over time, corresponding mainly to villages with high-rate clusters (HH; red) shown in the descriptive maps. Specifically, major high-risk clusters occurring between 2007 and 2012 (HH) comprised 3068 villages (Table 1), accounting for 77% of the total number of significant (0.05 level) villages and included a large geographic area in Southcentral Hubei Province (i.e., Jingzhou, Qianjiang, and Xiantao cities; Figure 2) and two small areas in Hanchuan and Yangxin. A few sporadic villages among the HH villages were identified as LH. Two others unchanged types (HL and LL) throughout the seven years ware mainly located in Caidian and Qujialing. In a small number of cases, comprising 93 villages (Table 1) scattered across Jiangzhou, Qianjiang, Xiantao, and Hanchuang (Figure 2), cluster/outlier types changed over time (L to H, and H to L).
Space–time scan analysis identified two spatiotemporal high-risk clusters (Figure 3 and Table 2). Primary clusters were detected within the period 2007–2009 and represented 2184 villages distributed in 10 counties. These were mainly embankment region epidemic areal-types in the middle-south of Hubei Province. Secondary clusters of embankment region epidemic areal-types, representing 125 villages and also throughout 2007–2009, were located in Northeastern Hanchuan. Both clusters remained statistically significant (p < 0.001). The two spatiotemporal high-risk clusters were basically consistent with the high cluster (HH) of cluster and outlier analysis.

3.2. Multilevel-Growth Model

Overall, 4379 surveyed villages from 2007 to 2012 were eligible for model inclusion. A total of 3501 villages were located in embankment regions, 446 in hilly regions, and 432 in marshland-and-lake regions, respectively (Figure 1). Table 3 summarizes the results of multilevel-growth models.

3.2.1. Unconditional Means Model

The variance components of the unconditional means model showed statistically significant variance associated with villages ( σ ^ 𝓊 0 2   = 1.045; p < 0.001) and statistically significant residual variance ( σ ^ 2   = 1.275; p < 0.001). The intraclass correlation coefficient calculated as 1.045 1.045 + 1.2757 = 0.45 indicated that approximately 45% of RIR variance arose from differences between villages [30,31].

3.2.2. Unconditional Growth Model

The unconditional growth model resulted in significant fit improvement with reduced −2LL, AIC, and BIC compared with the empty model (Table 3). The linear fixed effect of time differed significantly and negatively from zero ( γ ^ 10   = −0.351; p < 0.001), indicating that villages’ RIR values declined by about 0.35% per year on average, irrespective of their epidemic areal-type. Adding the fixed and random effects for time accounted for 73.6% of the explainable variance, as calculated using the Raudenbush and Bryk method [32].

3.2.3. Final Model

Table 3 summarizes the goodness-of-fit of the RIR models. The −2LL, AIC, and BIC for models, including the cross-level interactions model, was smaller than the corresponding value obtained for the other models, indicating that the assessing cross-level interaction model fitted the data better than the other models. After exploring the possibility that time (year), areal-type, and their interactions predicted differences in RIR values, the variance of mean initial status ( σ ^ 𝓊 0 2   = 3.809; p < 0.001) revealed that villages significantly varied in their initial RIR values and that the variance of rate of change ( σ ^ 𝓊 1 2   = 0.136; p < 0.001) indicated that village rates of change varied significantly. The estimated covariance ( σ ^ 𝓊 01 2   = −0.716; p < 0.001) was negative, indicating that villages with low initial RIR values tended to have higher rates of change, with statistical significance. The mean initial RIR value ( γ ^ 00 ) was 3.418%. A significant negative effect of areal-type on initial status ( γ ^ 01   = −0.896; p < 0.001) showed that these initial RIR values differed among the three epidemic areal-types, with higher initial values in embankment regions, medium initial values in marshland-and-lake regions, and lower initial values in hilly regions. The effect of areal-type on rate of change ( γ ^ 11   = 0.143; p < 0.001) was positive, indicating that RIR growth rate among epidemic areal-types significantly differed, on average, by 0.143% per year. Specifically, villages located in embankment regions reduced their RIR most rapidly over time, followed by villages located in marshland-and-lake regions and then by villages located in hilly regions.

4. Discussion

The heterogeneous geographic pattern of schistosomiasis RIR was found based on this spatial statistical approach, including HH clusters detected by local Moran’s I statistics and primary clusters detected by the space–time scan statistics. High-risk clusters for the schistosomiasis RIR were mainly identified in flooding areas of lakes connected to the Yangze and Hanjiang Rivers, particularly in the cities of Jingzhou, Qianjiang, and Xiantao, and in the embankment regions, which were consistent with the cluster regions at county level [12] and provided more details. These areas have natural geographical conditions and socioeconomic factors favorable to schistosomiasis; i.e., an eco-hydrological environment suitable for the survival of the intermediate snail host and its multiplication [33,34]. The plain region with waterway networks evolved from a typical lake and marshland to a schistosomiasis epidemic region connected by rivers and lakes because of farmland reclamation from lakes and large-scale excavation channels after the mid-1970s. Accordingly, the snail distribution evolved linearly from sheet distributions in the bottomland along the river and lakes, and along the ditches, irrigation canals, and rivers with conducive humidity, vegetation, elevation and flow velocity. The broad dispersion of snails, and the frequent flooding, as well as people’s frequent contact with contaminated water due to traditional livelihoods and lifestyles, is a challenge to effective control [35].
The results of the cluster and outlier evolution analyses for the 2007–2012 schistosomiasis RIR had 8 combination modes (including “not significant”) and various spatial distributions of combination modes (Figure 2 and Table 1), indicating that the different risk patterns require distinct control strategies and decision-making processes. This finding may be due to the different infectious sources and populations of schistosomiasis infection in different endemic areas and periods [36]. Additionally worth taking into account, aside from HH clusters and primary clusters, were secondary clusters located in Northeast Hanchuan (Figure 3), and villages with initially low rates that became high rates within a year (L to H, Figure 2), together with adjacent villages.
After exploring the effect of time (year), areal-type, and their interactions on RIR values, village RIR values significantly declined by about 0.537% per year on average (Table 3), which provided some explanation for the periods of two spatiotemporal high-risk clusters (Figure 3 and Table 2) detected by space–time scan analysis was the first three years. We also found temporal heterogeneity of schistosomiasis RIRs between the three epidemic areal-types. The embankment region with the widest area had the highest RIR values in the first year, which was possibly due to more frequent human activity and repeated infections compared to the marshland-and-lake, and hilly regions. Meanwhile, villages located in embankment regions most rapidly reduced their RIR over time, which may be because more attention was paid to integrated control strategies [37]. Significant differences existed in the initial RIR values and rates of change between the three epidemic areal-types, as well as within villages located in the same epidemic areal-type, which reminds us again that suitable comprehensive control measures must take account of the local epidemic situation rather than depend upon a uniform strategy.
In addition, villages in all epidemic regions with higher initial RIR values and faster rates of decline may, with time, have a sustained low prevalence of schistosomiasis. A low transmission of schistosomiasis has been detected in the hilly and mountainous endemic regions of Sichuan [38]. Consequently, new challenges toward the elimination of schistosomiasis japonica transmission in China may ensue because of the low sensitivity/specificity of current diagnostic tools for infections, praziquantel resistance, and climate change [38,39,40,41].
These analyses require more updated data for confirmation, which are not available, which is one limitation of our study. Moreover, two spatial analysis techniques were used to investigate spatial patterns, and although both results are broadly in agreement, no comparative study was conducted among different statistical methods because of their large diversity and different focuses on spatial patterns.
China continues to invest substantially in control activities aimed at eliminating schistosomiasis by 2020 [42]. Currently, the elimination of schistosomiasis depends on the sustained implementation of an integrated control strategy with great emphasis on the control of infectious sources [39,43,44]. For example, aspects to avoid infection in humans include: Improving sanitation by supplying tap water and building lavatories and latrines, providing boats with fecal-matter containers, and implementing an intensive health-education program and vaccination for long-term prevention [45]. Given all these different potential strategies, the importance of suiting the precise targeted control and priority of interventions to the local conditions is of concern [15,16,17].

5. Conclusions

This study systematically and exhaustively investigated the spatiotemporal dynamics of Schistosoma japonicum and temporal heterogeneity among regions of different epidemic areal-types by using spatial statistical analysis and multilevel-growth models at the microscale. Its high spatial resolution captured local heterogeneity, which could aid in regional targeting of schistosomiasis control measures and the spatiotemporal analysis method could be applied in other endemic areas of schistosomiasis.

Supplementary Materials

The following are available online at https://www.mdpi.com/1660-4601/16/12/2198/s1, Table S1: Combined clusters and outliers evolution modes and their categories.

Author Contributions

Conceptualization, J.Q., Y.W. and R.L.; methodology, J.Q.; formal analysis, Y.N., D.H., Q.S., Y.C.; data curation, J.X., Y.X. and H.Z.; writing—original draft preparation, J.Q.; writing—review and editing, Y.W. and R.L.; visualization, J.Q.; funding acquisition, R.L. and J.Q.

Funding

This research was funded by the CRSRI Open Research Program (Program SN: CKWV2017526/KY), and supported by National Nature Science Foundation of China (41571487), (81803297).

Acknowledgments

We acknowledge the support of anti-schistosomiasis stations from Yichang City, Jingzhou City, Qianjiang City, Jingmen City, Tianmen City, Xiaogan City, Xiantao City, Wuhan City, Xianning City, Ezhou City, Huanggang City and Huangshi City.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. King, C.H. Global Health Toward the Elimination of Schistosomiasis. N. Engl. J. Med. 2009, 360, 106–109. [Google Scholar] [CrossRef] [PubMed]
  2. Utzinger, J.; Zhou, X.N.; Chen, M.G.; Bergquist, R. Conquering schistosomiasis in China: The long march. Acta Trop. 2005, 96, 69–96. [Google Scholar] [CrossRef] [PubMed]
  3. Zhou, X.N.; Bergquist, R.; Leonardo, L.; Yang, G.J.; Yang, K.; Sudomo, M.; Olveda, R. Schistosomiasis Japonica: Control and Research Needs. In Important Helminth Infections in Southeast Asia: Diversity and Potential for Control and Elimination; Zhou, X.N., Bergquist, R., Olveda, R., Utzinger, J., Eds.; Academic Press Inc.: London, UK, 2010; Volume 72, pp. 145–178. [Google Scholar]
  4. Hu, Y.; Li, R.; Chen, Y.; Gao, F.H.; Wang, Q.Z.; Zhang, S.Q.; Zhang, Z.J.; Jiang, Q.W. Shifts in the Spatiotemporal Dynamics of Schistosomiasis: A Case Study in Anhui Province, China. PLoS Negl. Trop. Dis. 2015, 9, e0003715. [Google Scholar] [CrossRef] [PubMed]
  5. Chen, Y.Y.; Cai, S.X.; Liu, J.B.; Xiao, Y.; Li, G.; Shan, X.W.; Zhang, J. Epidemic and spatial distribution of Schistosomiasis in Hubei province from 2008 to 2012. Chin. J. Epidemiol. 2014, 35, 1366–1370. (In Chinese) [Google Scholar]
  6. Hu, Y.; Xiong, C.L.; Zhang, Z.J.; Luo, C.; Ward, M.; Gao, J.; Zhang, L.J.; Jiang, Q.W. Dynamics of spatial clustering of schistosomiasis in the Yangtze River Valley at the end of and following the World Bank Loan Project. Parasitol. Int. 2014, 63, 500–505. [Google Scholar] [CrossRef] [PubMed]
  7. Manyangadze, T.; Chimbari, M.J.; Gebreslasie, M.; Mukaratirwa, S. Risk factors and micro-geographical heterogeneity of Schistosoma haematobium in Ndumo area, uMkhanyakude district, KwaZulu-Natal, South Africa. Acta Trop. 2016, 159, 176–184. [Google Scholar] [CrossRef] [PubMed]
  8. Martins-Melo, F.R.; Pinheiro, M.C.C.; Ramos, A.N.; Alencar, C.H.; Bezerra, F.S.D.; Heukelbach, J. Spatiotemporal Patterns of Schistosomiasis-Related Deaths, Brazil, 2000–2011. Emerg. Infect. Dis. 2015, 21, 1820–1823. [Google Scholar] [CrossRef] [PubMed]
  9. Yang, K.; Li, W.; Sun, L.P.; Huang, Y.X.; Zhang, J.F.; Wu, F.; Hang, D.R.; Steinmann, P.; Liang, Y.S. Spatio-temporal analysis to identify determinants of Oncomelania hupensis infection with Schistosoma japonicum in Jiangsu province, China. Parasit. Vectors 2013, 6. [Google Scholar] [CrossRef] [PubMed]
  10. Meurs, L.; Mbow, M.; Boon, N.; van den Broeck, F.; Vereecken, K.; Dieye, T.N.; Abatih, E.; Huyse, T.; Mboup, S.; Polman, K. Micro-Geographical Heterogeneity in Schistosoma mansoni and S. haematobium Infection and Morbidity in a Co-Endemic Community in Northern Senegal. PLoS Negl. Trop. Dis. 2013, 7. [Google Scholar] [CrossRef] [PubMed]
  11. Yang, K.; Sun, L.P.; Liang, Y.S.; Wu, F.; Li, W.; Zhang, J.F.; Huang, Y.X.; Hang, D.R.; Liang, S.; Bergquist, R.; et al. Schistosoma japonicum risk in Jiangsu province, People’s Republic of China: Identification of a spatio-temporal risk pattern along the Yangtze River. Geospat. Health 2013, 8, 133–142. [Google Scholar] [CrossRef]
  12. Zhu, H.; Cai, S.X.; Liu, J.B.; Tu, Z.W.; Xia, J.; Shan, X.W.; Qiu, J.; Jiang, Y.; Xiao, Y.; Tang, L.; et al. A spatial analysis of human Schistosoma japonicum infections in Hubei, China, during 2009–2014. Parasit. Vectors 2016, 9. [Google Scholar] [CrossRef] [PubMed]
  13. Yang, K.; Zhou, X.-N.; Wu, X.-H.; Steinmann, P.; Wang, X.-H.; Yang, G.-J.; Utzinger, J.; Li, H.-J. Landscape Pattern Analysis and Bayesian Modeling for Predicting Oncomelania hupensis Distribution in Eryuan County, People’s Republic of China. Am. J. Trop. Med. Hyg. 2009, 81, 416–423. [Google Scholar] [CrossRef] [PubMed]
  14. Zhou, X.N. Science on Oncomelania Snail; Sciences Press: Beijing, China, 2005. [Google Scholar]
  15. Jia, T.W.; Guo, J.G.; Wu, X.H.; Zhang, L.J.; Wang, L.Y.; Lu, M.; Chen, C.; Zhou, X.N. Surveillance of schistosomiasis in 36 pilot villages implementing comprehensive control strategy with focus on infectious source in China, 2008. Zhongguo Xue Xi Chong Bing Fang Zhi Za Zhi 2010, 22, 333–338. (In Chinese) [Google Scholar]
  16. Qian, Y.L.; Wang, W.; Hong, Q.B.; Liang, Y.S. Bibliometric analysis of literature regarding integrated schistosomiasis control strategy with emphasis on infectious source control. Zhongguo Xue Xi Chong Bing Fang Zhi Za Zhi 2014, 26, 626–631. (In Chinese) [Google Scholar] [PubMed]
  17. Wang, J.S.; He, L.C.; Rong, X.B.; Yuan, M.Z.; He, Z.W.; Liao, C.J.; Hou, Y.Y.; Yuan, J.Y. Comparative evaluation of models for infectious source control of schistosomiasis in lake and marshland regions. Zhongguo Xue Xi Chong Bing Fang Zhi Za Zhi 2011, 23, 262–267. (In Chinese) [Google Scholar] [PubMed]
  18. Lei, Z.L.; Zhang, L.J.; Xu, Z.M.; Dang, H.; Xu, J.; Lv, S.; Cao, C.L.; Li, S.Z.; Zhou, X.N. Endemic status of schistosomiasis in People’s Republic of China in 2014. Zhongguo Xue Xi Chong Bing Fang Zhi Za Zhi 2015, 27, 563–569. (In Chinese) [Google Scholar] [PubMed]
  19. Qiu, J.; Li, R.D.; Xu, X.J. Characteristics analysis of Onocomelania hupensis’s geographical distribution in Hubei Province. Resour. Environ. Yangtze Basin 2012, 21, 100–104. (In Chinese) [Google Scholar]
  20. Anselin, L. Local Indicators of Spatial Association-LISA. Geogr. Anal. 1995, 27, 93–115. [Google Scholar] [CrossRef]
  21. Greene, S.K.; Ionides, E.L.; Wilson, M.L. Patterns of influenza-associated mortality among US elderly by geographic region and virus subtype, 1968–1998. Am. J. Epidemiol. 2006, 163, 316–326. [Google Scholar] [CrossRef] [PubMed]
  22. Pouliou, T.; Elliott, S.J. An exploratory spatial analysis of overweight and obesity in Canada. Prev. Med. 2009, 48, 362–367. [Google Scholar] [CrossRef] [PubMed]
  23. Hjalmars, U.; Kulldorff, M.; Wahlqvist, Y.; Lannering, B. Increased incidence rates but no space-time clustering of childhood astrocytoma in Sweden, 1973–1992: A population-based study of pediatric brain tumors. Cancer 1999, 85, 2077–2090. [Google Scholar] [PubMed]
  24. Kulldorff, M. A spatial scan statistic. Commun. Stat. Theory Methods 1997, 26, 1481–1496. [Google Scholar] [CrossRef]
  25. Kulldorff, M.; Huang, L.; Konty, K. A scan statistic for continuous data based on the normal probability model. Int. J. Health Geogr. 2009, 8. [Google Scholar] [CrossRef] [PubMed]
  26. Geert, V.; Geert, M. Linear Mixed Models for Longitudinal Data; Springer: New York, NY, USA, 2009. [Google Scholar]
  27. Singer, J.D. Using SAS PROC MIXED to fit multilevel models, hierarchical models, and individual growth models. J. Educ. Behav. Stat. 1998, 23, 323–355. [Google Scholar] [CrossRef]
  28. Singer, J.D.; Willett, J.B. Applied Longitudinal Data Analysis: Modeling Change and Event Occurrence; Oxford University Press: Oxford, UK, 2003. [Google Scholar]
  29. Garson, G.D. Hierarchical Linear Modeling: Guide and Applications; Garson, G.D., Ed.; SAGE Publications, Inc.: Thousand Oaks, CA, USA, 2013. [Google Scholar]
  30. Killip, S.; Mahfoud, Z.; Pearce, K. What is an intracluster correlation coefficient? Crucial concepts for primary care researchers. Ann. Fam. Med. 2004, 2, 204–208. [Google Scholar] [CrossRef]
  31. Shrout, P.E.; Fleiss, J.L. Intraclass correlations: Uses in assessing rater reliability. Psychol. Bull. 1979, 86, 420–428. [Google Scholar] [CrossRef] [PubMed]
  32. Raudenbush, S.; Bryk, A. Hierarchical Linear Models: Applications and Data Analysis Methods, 2nd ed.; Sage Publications, Inc.: Thousand Oaks, CA, USA, 2002. [Google Scholar]
  33. Qiu, J.; Li, R.D.; Xu, X.J.; Yu, C.H.; Xia, X.; Hong, X.C.; Chang, B.R.; Yi, F.J.; Shi, Y.Y. Identifying Determinants of Oncomelania hupensis Habitats and Assessing the Effects of Environmental Control Strategies in the Plain Regions with the Waterway Network of China at the Microscale. Int. J. Environ. Res. Public Health 2014, 11, 6571–6585. [Google Scholar] [CrossRef]
  34. Wang, W.L. Schistosomiasis epidemic characteristics after lake beach type region evolving into irrigation network endemic area. Chin. J. Parasitol. Parasit. Dis. 1990, 8, 311. (In Chinese) [Google Scholar]
  35. Peng, X.W.; Wang, J.S.; Rong, X.B.; Wang, S.H.; He, L.C. Schistosomiasis prevalence characteristic and epidemic factors in Four-lake basin of Jianghan Plain, Hubei province. J. Pub. Health Prev. Med. 2007, 18, 66–68. (In Chinese) [Google Scholar]
  36. Ding, Z.J.; Peng, G.Y. Discussion of schistosomiasis infection source control. Chin. J. Schisto Control 2012, 24, 117–118. (In Chinese) [Google Scholar]
  37. Zhu, H.; Huang, X.B.; Cai, S.X.; Tu, Z.W.; Chen, Y.Y.; Li, G.; Xia, J.; Zhang, J.; Wang, W.; Xiang, D.J. Analysis of the economic benefits of integrated control strategies emphasizing the sources of schistosomiasis in different areas where schistosomiasis is prevalent. J. Pathogen Biolog. 2012, 7, 679–681. [Google Scholar]
  38. Spear, R.C.; Seto, E.Y.W.; Carlton, E.J.; Liang, S.; Remais, J.V.; Zhong, B.; Qiu, D. The challenge of effective surveillance in moving from low transmission to elimination of schistosomiasis in China. Int. J. Parasitol. 2011, 41, 1243–1247. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Lei, Z.L.; Zhou, X.N. Eradication of schistosomiasis: A new target and a new task for the National Schistosomiasis Control Porgramme in the People’s Republic of China. Zhongguo Xue Xi Chong Bing Fang Zhi Za Zhi 2015, 27, 1–4. (In Chinese) [Google Scholar] [PubMed]
  40. Wang, W.; Dai, J.R.; Liang, Y.S. Apropos: Factors impacting on progress towards elimination of transmission of Schistosomiasis japonica in China. Parasit. Vectors 2014, 7. [Google Scholar] [CrossRef] [PubMed]
  41. Zhou, Y.B.; Liang, S.; Jiang, Q.W. Factors impacting on progress towards elimination of transmission of Schistosomiasis japonica in China. Parasit. Vectors 2012, 5. [Google Scholar] [CrossRef] [PubMed]
  42. Yang, G.J.; Liu, L.; Zhu, H.R.; Griffiths, S.M.; Tanner, M.; Bergquist, R.; Utzinger, J.; Zhou, X.N. China’s sustained drive to eliminate neglected tropical diseases. Lancet Infect. Dis. 2014, 14, 881–892. [Google Scholar] [CrossRef]
  43. Wang, L.D.; Chen, H.G.; Guo, J.G.; Zeng, X.J.; Hong, X.L.; Xiong, J.J.; Wu, X.H.; Wang, X.H.; Wang, L.Y.; Xia, G.; et al. A Strategy to Control Transmission of Schistosoma japonicum in China. N. Engl. J. Med. 2009, 360, 121–144. [Google Scholar] [CrossRef]
  44. Xiao, Y.; Zhang, J.; Li, G.; Wei, Z.Y.; Jiang, N.; Jun, H.; Cai, S.X. Retrospective investigation on endemic situation of schistosomiasis in Hubei Province. Zhongguo Xue Xi Chong Bing Fang Zhi Za Zhi 2014, 26, 486–490. (In Chinese) [Google Scholar]
  45. Gray, D.J.; Li, Y.S.; Williams, G.M.; Zhao, Z.Y.; Harn, D.A.; Li, S.M.; Ren, M.Y.; Feng, Z.; Guo, F.Y.; Guo, J.G.; et al. A multi-component integrated approach for the elimination of schistosomiasis in the People’s Republic of China: Design and baseline results of a 4-year cluster-randomised intervention trial. Int. J. Parasitol. 2014, 44, 659–668. [Google Scholar] [CrossRef]
Figure 1. Study area and epidemic areal-types of schistosomiasis. The inset shows the location of the 30 most endemic counties in Hubei and in China. The geographical layers of Yangtze and Han Rivers are overlaid.
Figure 1. Study area and epidemic areal-types of schistosomiasis. The inset shows the location of the 30 most endemic counties in Hubei and in China. The geographical layers of Yangtze and Han Rivers are overlaid.
Ijerph 16 02198 g001
Figure 2. Cluster and outlier evolution for schistosomiasis RIR (resident infection rate) from 2007 to 2012 in the 30 most endemic counties, Hubei Province.
Figure 2. Cluster and outlier evolution for schistosomiasis RIR (resident infection rate) from 2007 to 2012 in the 30 most endemic counties, Hubei Province.
Ijerph 16 02198 g002
Figure 3. Mapping of significant spatiotemporal clusters of schistosomiasis RIR, by administrative village.
Figure 3. Mapping of significant spatiotemporal clusters of schistosomiasis RIR, by administrative village.
Ijerph 16 02198 g003
Table 1. Combined cluster/outlier found by the local Moran’s I analysis for schistosomiasis RIR (resident infection rate).
Table 1. Combined cluster/outlier found by the local Moran’s I analysis for schistosomiasis RIR (resident infection rate).
TypeNo. of VillagesDescription of Types
HH3068a high cluster that is statistically significant (0.05 level) in all years or in some years from 2007 to 2012
LH631a low outlier that is statistically significant (0.05 level) in all years or in some years from 2007 to 2012
HL90a high outlier that is statistically significant (0.05 level) in some years from 2007 to 2012
LL97a low cluster that is statistically significant (0.05 level) in some years from 2007 to 2012
H to L52a high-value village or its neighbors converted into a low-value village or its neighbors, i.e., a high cluster converted into a high outlier (HH to HL), or a high cluster converted into a low outlier (HH to LH), or a low outlier converted into a low cluster (LH to LL)
L to H41a low-value village or its neighbors converted into a high value village or its neighbors, i.e., a high outlier converted into a high cluster (HL to HH), or a low outlier converted into a high cluster (LH to HH)
Others8Others (Table S1)
Table 2. Significant spatiotemporal clusters of schistosomiasis RIR as defined by space–time scan statistic in the 30 most endemic counties, Hubei Province, 2007–2012.
Table 2. Significant spatiotemporal clusters of schistosomiasis RIR as defined by space–time scan statistic in the 30 most endemic counties, Hubei Province, 2007–2012.
Cluster 1PeriodCluster LocationCluster Radius (km)No. VillagesCountyLLR 2p Value
12007–200929.485 N, 112.956 E113.942184Qianjian, Xiantao, Honghu, Chibi, Jianli, Shishou, Gong’an, Jiangling, Shashi, Jiayu3764.32<0.001
22007–200930.720 N, 113.766 E13.72125Hanchuan176.65<0.001
1 The most likely or primary clusters (1) and secondary clusters (2) were detected by the LLR. The most likely cluster was defined as the one with the maximum LLR. 2 LLR: loglikelihood ratio test.
Table 3. Results of fitting multilevel-growth models to schistosomiasis RIR in 4379 villages in the 30 most endemic counties, Hubei Province, China.
Table 3. Results of fitting multilevel-growth models to schistosomiasis RIR in 4379 villages in the 30 most endemic counties, Hubei Province, China.
Unconditional Means (Empty Model)Unconditional GrowthAdding a Village-Level CovariateAssessing Cross-Level Interactions
Fixed Effects
Initial Status ( γ ^ 00 )1.3742.2522.5053.418
Year (Rate of Change) ( γ ^ 10 ) −0.351−0.351−0.537
Areal-type (Effect of areal-type on initial status) ( γ ^ 01 ) −0.194−0.896
Areal-type (Effect of areal-type on rate of Change) ( γ ^ 11 ) 0.143
Random Effects
Level 1 (within village)
Residual ( σ ^ 2 )1.2750.3370.3370.337
Level 2 (between village)
Village Mean Initial Status ( σ ^ 𝓊 0 2 )1.0454.1424.0133.809
Village Mean Rate of Change ( σ ^ 𝓊 1 2 ) 0.1450.1450.136
Rate of change covariance ( σ ^ 𝓊 01 2 ) −0.769−0.758−0.716
Fit Statistics
−2LL88,737.361,858.361,586.161,354.3
AIC88,743.361,870.361,600.161,370.3
BIC88,762.461,908.661,644.861,421.4

Share and Cite

MDPI and ACS Style

Qiu, J.; Li, R.; Xiao, Y.; Xia, J.; Zhu, H.; Niu, Y.; Huang, D.; Shao, Q.; Cui, Y.; Wang, Y. Spatiotemporal Heterogeneity in Human Schistosoma japonicum Infection at Village Level in Hubei Province, China. Int. J. Environ. Res. Public Health 2019, 16, 2198. https://doi.org/10.3390/ijerph16122198

AMA Style

Qiu J, Li R, Xiao Y, Xia J, Zhu H, Niu Y, Huang D, Shao Q, Cui Y, Wang Y. Spatiotemporal Heterogeneity in Human Schistosoma japonicum Infection at Village Level in Hubei Province, China. International Journal of Environmental Research and Public Health. 2019; 16(12):2198. https://doi.org/10.3390/ijerph16122198

Chicago/Turabian Style

Qiu, Juan, Rendong Li, Ying Xiao, Jing Xia, Hong Zhu, Yingnan Niu, Duan Huang, Qihui Shao, Ying Cui, and Yong Wang. 2019. "Spatiotemporal Heterogeneity in Human Schistosoma japonicum Infection at Village Level in Hubei Province, China" International Journal of Environmental Research and Public Health 16, no. 12: 2198. https://doi.org/10.3390/ijerph16122198

APA Style

Qiu, J., Li, R., Xiao, Y., Xia, J., Zhu, H., Niu, Y., Huang, D., Shao, Q., Cui, Y., & Wang, Y. (2019). Spatiotemporal Heterogeneity in Human Schistosoma japonicum Infection at Village Level in Hubei Province, China. International Journal of Environmental Research and Public Health, 16(12), 2198. https://doi.org/10.3390/ijerph16122198

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop