This article has Open Peer Review reports available.
Analysis of the geographic distribution of HFRS in Liaoning Province between 2000 and 2005
© Lin et al; licensee BioMed Central Ltd. 2007
Received: 17 October 2006
Accepted: 15 August 2007
Published: 15 August 2007
Hemorrhagic fever with renal syndrome (HFRS) is endemic in Liaoning Province, China, and this province was the most serious area affected by HFRS during 2004 to 2005. In this study, we conducted a spatial analysis of HFRS cases with the objective to determine the distribution of HFRS cases and to identify key areas for future public health planning and resource allocation in Liaoning Province.
The annual average incidence at the county level was calculated using HFRS cases reported between 2000 and 2005 in Liaoning Province. GIS-based spatial analyses were conducted to detect spatial distribution and clustering of HFRS incidence at the county level, and the difference of relative humidity and forestation between the cluster areas and non-cluster areas was analyzed.
Spatial distribution of HFRS cases in Liaoning Province from 2000 to 2005 was mapped at the county level to show crude incidence, excess hazard, and spatial smoothed incidence. Spatial cluster analysis suggested 16 and 41 counties were at increased risk for HFRS (p < 0.01) with the maximum spatial cluster sizes at ≤ 50% and ≤ 30% of the total population, respectively, and the analysis showed relative humidity and forestation in the cluster areas were significantly higher than in other areas.
Some clustering of HFRS cases in Liaoning Province may be etiologically linked. There was strong evidence some HFRS cases in Liaoning Province formed clusters, but the mechanism underlying it remains unknown. In this study we found the clustering was consistent with the relative humidity and amount of forestation, and showed data indicating there may be some significant relationships.
HFRS is a zoonosis caused by different species of hantavirus (HV) causing fever, hemorrhage, kidney damage and hypotension. HFRS in China has the most severe epidemics having about 90% of the total reported HFRS cases in the world [1–4]. Although some prevention and control measures such as scientific rodent control, vaccination and environment management have been performed, HFRS remains a serious public health problem with about 20,000~50,000 human cases annually in mainland China . And Liaoning Province is one of the most serious affected areas with the most cases in China and the highest incidence during the years 2004 and 2005. The incidence of HFRS shows high variability at the county level. The environmental factors, climatic factors, low social and economic development, and lack of health consciousness are thought to be related to the high incidence and spatial distribution of HFRS [3, 4]. A better understanding of the spatial distribution patterns of HFRS and the relationship between the high incidence and relevant factors would help to identify areas and population at high risk and may better prevent and control HFRS.
Spatial analysis with spatial smoothing and cluster analysis are commonly used to characterize spatial patterns of disease [6–10]. Spatial smoothing is used to reduce random variation associated with small populations and enables observation of gradients or holes in disease incidence that may not be apparent from direct observation of the raw data. Spatial cluster analysis is conducted to identify whether cases of disease are geographically clustered [11–13]. In this study, we conducted GIS-based spatial analysis involving spatial smoothing and spatial clustering analysis to characterize geographic distribution patterns of HFRS cases. Spatial analysis was used to identify the distribution pattern of HFRS and population at high risk at the county level. The technique corrects for multiple comparisons, adjusts for the heterogeneous population densities among the different areas, detects foci without prior specification of suspected location or size, and thereby overcomes pre-selection bias and allows for adjustment of confounders [13, 14].
Data collection and management
The study site is located in Liaoning Province, one of the three northeastern provinces of China (118°53'~125°46', 38°43'~43°26'N). The maximum distance from east to west is about 550km and from north to south is also about 550km. The area has a population of 41,940,051 (fifth national census in 2000) and encompasses 145,900 square kilometers.
There are 19 counties, 8 autonomic counties, 17 county-level cities, and 14 cities in Liaoning Province. For our study, we considered the city-governed region and others as county level, therefore there are 58 county-level study areas. The information includes the number of HFRS cases per month and per year in every county from 2000 to 2005. For the 6-year period (2000~2005), the average annual incidence was 8.65 cases per 100,000 persons, 21,760 cases were involved in this study.
HFRS cases were first diagnosed using clinical symptoms, then blood samples were collected in the hospitals, serologic identification was performed at the laboratory of Liaoning Provincial Center for Disease Control and Prevention (CDC) to confirm the clinical diagnosis, and then the data was collected by case number according to the sampling results. There might be admission rate bias in the disease report, but this has been reduced as much as possible.
GIS mapping and smoothing for the incidence of HFRS
Records on HFRS cases between 2000 and 2005 were obtained from Liaoning Provincial CDC. To conduct a GIS-based analysis on the spatial distribution of HFRS, a county-level polygon map at a scale of 1:1,000,000 was obtained and on which the county-level point layers containing information regarding latitudes and longitudes of central points of each county were created. Demographic information based on the fifth national census (2000) was integrated in terms of an administration code . All HFRS cases were geocoded and matched to the county-level layers on the polygons and points by administration code using the software ArcGIS9.1.
Based on annual average incidence, all counties were grouped into four categories: non-endemic areas, low endemic areas with an annual average incidence between 0 and 5 per 100,000, medium endemic areas with an incidence between 5 and 30 per 100,000, and high endemic areas with an incidence over 30 per 100,000. The four types of counties were color-coded on the maps.
To reduce variations of incidence in areas with a small population, annual average incidence of HFRS per 100,000 for each county over the 6-year period were calculated and spatial smoothing was used.
Spatial cluster analysis
A "cluster" used here is a real or perceived aggregation of more than the expected number of HFRS cases in a population over a specified time period. In the present study, spatial cluster analysis was also performed on the confirmed HFRS cases to test whether the cases were distributed randomly over space and time, and if not, to evaluate any identified spatial and temporal disease clusters for statistical significance. A spatial scan statistic (SaTScan)  was then applied to identify clusters of HFRS. The spatial scan statistic places a circular window of varying size on the map surface and allows its center to move so that at any given position and size, the window includes different sets of adjacent neighborhoods. If the window contains the center of a neighborhood, then the whole neighborhood is included in the window. As the window is placed at each neighborhood center, its radius is varied continuously from zero up to a maximum radius which never includes more than 50% of the total population. The method creates a large number of distinct circular windows, each containing a distinct set of adjacent neighborhoods, and each is a possible candidate for containing a cluster of HFRS cases. For each window, the method uses a Monte Carlo simulation to test the null hypothesis that there is not an elevated risk of HFRS. Details of how the likelihood function is maximized over all windows under the Poisson assumption have been described elsewhere [12, 13, 15]. In this study, retrospective spatial cluster analysis for higher incidence was used, in which the maximum window radius was set to be smaller than 50% and 30% of the total population, respectively, to find possible sub-clusters. For each window of movable position and size change, the software tested the risk of HFRS within and outside the window using the null hypothesis of the same risk.
The statistical analysis of relative humidity and forestation
Data analysis was conducted using the Statistical Package for the Social Sciences (SPSS). T-test analysis was conducted on the average relative humidity over the six-year period and forestation in every county between the cluster areas and the other areas.
Spatial distribution of HFRS in Liaoning Province
Spatial clustering of HFRS in Liaoning Province
SaTScan statistics for the most likely cluster, Liaoning Province, China, 2000–2005. (Total population = 41,940,051 Total cases =21,760)
Most likely cluster
The 1st iteration1
The 2nd iteration2
The incidence of HFRS compared to relative humidity and forestation
The comparison of relative humidity and forestation between cluster areas and non-cluster areas
1st cluster and others
2nd cluster and non-cluser
Spatial variation in disease distribution occurs where multiple statistical methods have been developed to determine whether patterns of variation occur by chance alone, or whether variation is unlikely to have happened at random [16–18]. Using GIS-based spatial statistics, we investigated the spatial distribution of HFRS cases and identified counties with high endemic HFRS and clustering patterns. The Martin Kulldorff's spatial scan test has several advantages over other local cluster detection tests, namely, a formal likelihood test statistic, avoidance of multiple testing, and not restricting scanning for clusters of a pre-specified size, thereby avoiding pre-selection bias . The spatial statistics analysis yields a nonrandom distribution of HFRS within the 58-county area. We identified the eastern counties as having a cluster in Liaoning Province which was not previously reported. The relative risk in this cluster was at the highest level of 5.121 (50% analysis) and 5.898 (30% analysis). Although much is still unknown about HFRS epidemic initiation and development, there is sufficient data presented here to argue that HFRS is most likely influenced by a complex combination of factors rather than a single foci pathogenic factor. It is likely that many factors contribute to its epidemic and development from natural environmental agents to economic and social factors.
Liaoning Province can be divided into three parts in respect to terrain: the Liaodong mountainous region, the Liaohe plain region, and the Liaoxi mountainous region . We found as shown in Figure 4 the relatively higher risk areas are in the eastern and western mountainous regions where these regions have a humid or semi-humid climate with an annual average precipitation of 650–1,200 mm and 500 – 600 mm, respectively . There is more forest in these regions than on the Liaohe plain, which contribute to a habitat for rodents, the reservoir and vector of HV, and increased opportunities of human contact with rodents and their excreta. Therefore, the terrain, the amount of forestation and the relative high humidity may be the important factors responsible for the epidemic development of HFRS in Liaoning Province.
As mapping and spatial analysis software become easier to use, public health professionals can increasingly employ these tools to aid in disease surveillance and cluster investigations . The spatial scan statistic method approach illustrated here is a useful instrument for HFRS analysis of cluster patterns. Further, space-time disease-surveillance methods have been proposed as a dynamic supplement to purely spatial statistical methods for outbreak detection to detect and predict localized outbreaks before they spread to larger regions. Therefore, this study should be used as a preliminary study showing only spatial patterns. In the areas in which HFRS are highly epidemic, aiming prevention strategies at areas of highest risk can potentially increase the health maintenance program's effectiveness. Persons at highest risk should be informed of the risks and methods for risk reduction. Funds spent on these programs may be better spent in the areas identified here where cost-effectiveness can be maximized. The data presented here will assist state and local health departments for HFRS surveillance systems in the investigation of HFRS clusters by providing definitions and recommendations to aid in the development of a better prevention strategy.
This study analyzed the spatial distribution of HFRS in Liaoning Province, China, during 2000 to 2005 using the spatial smoothing and spatial scan statistic method. SaTScan identified a geographic area in eastern Liaoning Province as the most likely endemic cluster. This has not been previously reported. The data provide evidence there are environment factors including the amount of forestation and climatic factors that may contribute to cause these clusters. The data show practical HFRS control measures, as well as methods for future study of HFRS and other vector-borne diseases. Further, GIS and GIS-based spatial statistical techniques may provide an opportunity to clarify and quantify the epidemic situation of HFRS within highly epidemic areas, and lay a foundation to pursue future investigations into the environmental factors responsible for the increased disease risk. To implement specific and geographically appropriate risk-reduction programs, the use of such spatial analysis tools should become an integral component in epidemiology research and risk assessment of HFRS.
The authors extend appreciation to Zhi Li and Yingwei Sun for providing the data, and Haixia Wu, Liang Lu, Fengxia Meng, Jian Sun, Lingling Lang, Lei Wang, Jimin Sun, Haiqi Wang, and Yilan Liao for help in the process of the study.
We would also like to acknowledge Dr. Weizhong Yang who contributed to the revision of the manuscript and the final manuscript.
This study was supported by the social welfare research special program of the Ministry of Science and Technology, China (grant number: 2003DIA6N009).
- Chen HX, Luo CW: The surveillance of haemorrhagic fever with renal syndrome in China, 1931–95. Chin J Vector Biol Control. 1996, 7: 1-26.Google Scholar
- Bai X, Huang C: Study farther on hemorrhagic fever with renal syndrome. Chin J Infect Dis. 2002, 20: 197-198.Google Scholar
- Chen HX, Qiu FX: Studies on the environment structure of natural nidi and epidemic areas of hemorrhagic fever with renal syndrome in China. Chin Med J. 1994, 107: 107-112.PubMedGoogle Scholar
- Bi P, Tong SL, Donald K: Climatic, reservoir and occupational variables and the transmission of haemorrhagic fever with renal syndrome in China. International Journal of Epidemiology. 2002, 31: 189-193. 10.1093/ije/31.1.189.View ArticlePubMedGoogle Scholar
- Fang L, Yan L, Liang S, de Vlas SJ, Feng D, Han X, Zhao W, Xu B, Bian L, Yang H, Gong P, Richardus JH, Cao W: Spatial analysis of hemorrhagic fever with renal syndrome in China. BMC Infectious Diseases. 2006, 6: 77-10.1186/1471-2334-6-77.View ArticlePubMedPubMed CentralGoogle Scholar
- Frank C, Fix AD, Peña CA, Strickland GT: Mapping Lyme disease for diagnostic and preventive decisions, Maryland. Emerg Infect Dis. 2002, 8: 427-429.View ArticlePubMedPubMed CentralGoogle Scholar
- Odoi A, Martin SW, Michel P, Middleton D, Holt J, Wilson J: Investigation of clusters of giardiasis using GIS and a spatial scan statistic. Int J of Health Geog. 2004, 3: 11-21. 10.1186/1476-072X-3-11.View ArticleGoogle Scholar
- Nkhoma ET, Ed Hsu C, Hunt VI, Harris AM: Detecting spatiotemporal clusters of accidental poisoning mortality among Texas counties, U.S., 1980–2001. Int J Health Geog. 2004, 3: 25-37. 10.1186/1476-072X-3-25.View ArticleGoogle Scholar
- Wu J, Wang J, Meng B, Chen G, Pang L, Song X, Zhang K, Zhang T, Zheng X: Exploratory spatial data analysis for the identification of risk factors to birth defects. BMC Public Health. 2004, 4: 23-10.1186/1471-2458-4-23.View ArticlePubMedPubMed CentralGoogle Scholar
- Zeman P: Objective assessment of risk maps of tick-borne encephalitis and Lyme borreliosis based on spatial patterns of located cases. Int J Epidemiol. 1997, 26: 1121-1129. 10.1093/ije/26.5.1121.View ArticlePubMedGoogle Scholar
- Kulldorff M, Feuer EJ, Miller BA, Freedman LS: Breast cancer clusters in the northeast United States: a geographic analysis. Am J Epidemiol. 1997, 146: 161-170.View ArticlePubMedGoogle Scholar
- Kulldorff M, Nagarwalla N: Spatial disease clusters: detection and inference. Stat Med. 1995, 14: 799-810. 10.1002/sim.4780140809.View ArticlePubMedGoogle Scholar
- Kulldorff M: A spatial scan statistic. Communications in Statistics. Theory and Methods. 1997, 26: 1481-1496.View ArticleGoogle Scholar
- Chaput EK, Meek JI, Heimer R: Spatial analysis of human granulocytic ehrlichiosis near Lyme, Connecticut. Emerg Inf Dis. 2002, 8 (9): 943-948.View ArticleGoogle Scholar
- Kulldorff M: Information Management Services, Inc: SaTScan™ v6.0: Software for the spatial and space-time scan statistics. 2005, [http://www.satscan.org/]Google Scholar
- Fukuda Y, Umezaki M, Nakamura K: Variations in societal characteristics of spatial disease clusters:examples of colon, lung and breast cancer in Japan. International Journal of Health Geographics. 2005, 4: 16-10.1186/1476-072X-4-16.View ArticlePubMedPubMed CentralGoogle Scholar
- Greena C, Robert D, Hoppab T: Geographic analysis of diabetes prevalence in an urban area. Social Science & Medicine. 2003, 57: 551-560. 10.1016/S0277-9536(02)00380-5.View ArticleGoogle Scholar
- Chen Y, Zhang KL, Wang JF, Zheng XY, Bian L: Research of the relation between birth defects and environment in Heshun County of Shanxi Province. Guang Dong Wei Liang Yuan Su Ke Xue. 2004, 11 (12): 21-28.Google Scholar
- Mather FJ, Chen VW, Morgan LH, Correa CN, Shaffer JG, Srivastav SK, Rice JC, Blount G, Swalm CM, Wu X, Scribner RA: Hierarchical modeling and other spatial analyses in prostate cancer incidence data. Am J Prev Med. 2006, 30 (2S): 88-100. 10.1016/j.amepre.2005.09.012.View ArticleGoogle Scholar
- Wang YJ, Liu XF: Research of heamorrhagic fever with renal syndrome foci and the landscape characteristics. Chin J Vector Bio & Control. 2003, 14: 255-268.Google Scholar
- Klassen AC, Kulldorff M, Curriero F: Geographical clustering of prostate cancer grade and stage at diagnosis, before and after adjustment for risk factors. International Journal of Health Geographics. 2005, 4: 1-10.1186/1476-072X-4-1.View ArticlePubMedPubMed CentralGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://0-www.biomedcentral.com.brum.beds.ac.uk/1471-2458/7/207/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.