Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Spatio-Temporal Pattern and Influencing Factors of Hemorrhagic Fever with Renal Syndrome (HFRS) in Hubei Province (China) between 2005 and 2014

  • Liang Ge ,

    geliang0021@126.com

    Affiliations State Key Laboratory of Information Engineering in Surveying, Mapping and Remote Sensing, Wuhan University, Wuhan city, Hubei Province, PR China, Tianjin Institute of Surveying and Mapping, Tianjin city, PR China

  • Youlin Zhao,

    Affiliation Business School of Hohai University, Nanjing city, Jiangsu Province, PR China

  • Kui Zhou,

    Affiliation Tianjin Institute of Surveying and Mapping, Tianjin city, PR China

  • Xiangming Mu,

    Affiliation School of Information Studies in University of Wisconsin-Milwaukee 2025 E Newpot Ave #NWQB, Milwaukee, WI, United States of America

  • Haibo Yu,

    Affiliation Tianjin Institute of Surveying and Mapping, Tianjin city, PR China

  • Yongfeng Wang,

    Affiliation Tianjin Institute of Surveying and Mapping, Tianjin city, PR China

  • Ning Wang,

    Affiliation First Crust Deformation Monitoring and Application Center, China Earthquake administration, Tianjin city, PR China

  • Hong Fan,

    Affiliation State Key Laboratory of Information Engineering in Surveying, Mapping and Remote Sensing, Wuhan University, Wuhan city, Hubei Province, PR China

  • Liqiang Guo,

    Affiliation Tianjin Institute of Surveying and Mapping, Tianjin city, PR China

  • XiXiang Huo

    Affiliation Hubei Provincial Center for Disease Control and Prevention, Wuhan, China

Abstract

Hemorrhagic Fever with Renal Syndrome (HFRS) is considered as a globally distributed infectious disease, which results in many deaths annually in Hubei Province, China. The outbreak of HFRS is usually characterized with spatio-temporal heterogeneity and is seasonally distributed. Further, it might also be impacted by the influencing factors such as socio-economic and geographical environment. To better understand and predict the outbreak of HFRS in the Hubei Province, the spatio-temporal pattern and influencing factors were investigated in this study. Moran’s I Index value was adopted in spatial global autocorrelation analysis to identify the overall spatio-temporal pattern of HFRS outbreak. Kulldorff scan statistical analysis was performed to further identify the changing trends of the clustering patterns of HFRS outbreak. Spearman's rank correlation analysis was used to explore the possible influencing factors on HFRS epidemics such as climate and geographic. The results demonstrated that HFRS outbreak in Hubei Province decreased from 2005 to 2012 in general while increasing slightly from 2012 to 2014. The spatial and temporal scan statistical analysis indicated that HFRS epidemic was temporally clustered in summer and autumn from 2005 to 2014 except 2008 and 2011. The seasonal epidemic pattern of HFRS in Hubei Province was characterized by a bimodal pattern (March to May and September to November) while peaks often occurring in the spring time. SEOV-type HFRS was presumed to influence more on the total number of HFRS incidence than HTNV-type HFRS do. The average humidity and human population density were the main influencing factors during these years. HFRS outbreaks were more in plains than in other areas of Hubei Province. We did not find that whether the terrain of the wetland (water system) plays a significant role in the outbreak of HFRS incidence. With a better understanding of rodent infection rate, socio-economic status and ecological environment characteristics, this study may help to reduce the outbreak of HFRS disease.

Introduction

In China, HFRS was mainly caused by two types of Hantaviruses, named Hantaan virus (HTNV) and Seoul virus (SEOV), each associated with a unique rodent host [13]. In China, the first case of HFRS was discovered in 1935. After that, data showed that the number of HFRS cases in China accounts for 90% of the globally reported cases over the past 20 years [47].

Generally, human activities and natural factors were related to the occurrence and epidemic of hantaviruses [810]. In Europe, in order to discover the regular pattern and feature of HFRS, hank vole dynamics, soil contamination dynamics and human contamination dynamics were used to explain the spatial variation of HFRS outbreak [11]. In China, some researchers used cluster analysis to study the relationships between the spatial distribution and the influencing factors of the HFRS outbreaks to explore the degree of clustering. Poisson regression analysis was performed by Zhang et al. to identify the HFRS transmission pattern in north-eastern China from year 1997 to 2007 [12]. Climate factors (e.g. monthly rainfall, relative humidity, and land surface temperature) were found to be the determinants to HFRS transmission in this research. Moran’s I spatial autocorrelation statistical and retrospective spatio-temporal clustering methods were used by Wu et al. to analyze the spatio-temporal distribution in Liaoning Province from 1988 to 2001. They demonstrated that the outbreak of HFRS had homogeneous spatio-temporal characteristics [5]. Lin et al performed Spatial smoothing and Martin Kulldorff’s spatial scan test to study the spatial distribution and variation of the HFRS outbreak in Liaoning Province between 2000 and 2005 [4]. They found that the clusters of HFRS cases were consistently influencing by humidity and the amount of forestation. In another study, ARIMA model and historical time series data were utilized by Liu et al. to simulate the temporal distribution tendency of HFRS in China from 1975 to 2008. The results demonstrated that ARIMA model had a good feasibility to forecast the HFRS outbreak [13].

Although a variety of methods have been implemented to reduce the occurrence of HFRS, the HFRS cases were still more than 20,000 annually in China from 1980 to 2009 according to the report [12]. In Hubei Province, the number of HFRS cases reached 23,943 cases in 1983 [12,14,15], and get 104,467 cases in total between 1980 and 2009 [11,16,17]. The relationships between natural environment and number of HFRS outbreaks were explored in Russia and Korea in the 1990s [18]. Since 2000, for the outbreak of HFRS in northeast China, most researches began to focus on a short term (five years or four years) relationships between dynamic spatial distribution and possible influencing factors [12,15]. In this study, we will explore a ten year appeared relationships between dynamic spatial distribution and possible influencing factors in Hubei Province.

Liu et al. found that temporal and geographic factors both played important roles on the outbreaks of HFRS in central-south China between 2000 and 2009 [19]. Zhang et al. studied the spatio-temporal distribution pattern of HFRS in Hubei Province over 1980–2009, they found that the number of HFRS cases was decreasing in general during this period [14]. Through the data from the experiment analysis in Hubei Province, this study will explore the distribution tendencies of clustering areas and the dominated influencing factors on a year-by-year base. The significant of this study is to characterize the overall spreading tendency and moving patterns of HFRS in recent years in Hubei Province.

Research Methodology

Study area and data source

Fig 1 displays the location of Hubei Province in China. Hubei Province is located in the central south part of China. In the middle of the Yangtze River with an area of 186,000 square kilometres, it is located at 108'21"-116'07" east longitude and 29'05"-33'20" north latitude. The population was about 59.37 million in 2005. Up to 2013, it had one autonomous prefecture, 12 prefecture-level cities, 24 county-level cities, 38 counties, two autonomous counties and one forest district. Administrative region of counties in Hubei Province have changed for several times. In order to organize the data effectively during data processing, 76 counties have been integrated according to the latest Chinese Administrative Coding Rules.

thumbnail
Fig 1. The location of Hubei Province in China.

The highlighted block in the above map indicates its location. The underneath map indicates the county boundary of Hubei Province.

https://doi.org/10.1371/journal.pone.0167836.g001

Four data sources were used in this study. First, the number of HFRS cases of 76 counties in Hubei Province was provided by Hubei Provincial Centre for Disease Control and Prevention (CDC) between 2005 and 2014. Secondly, human population data was provided by Hubei Provincial Bureau of Statistics based on the 2009 censuses. Thirdly, the map at a scale of 1:4000, 000 was obtained by National Geomatics Centre of China. Fourthly, climatic data from 2005 to 2014 was provided partly by Hubei Provincial Meteorological Bureau and partly from online open sources. These four sources of data were stored in the geographical information database and attribute database. Geographical information database was linked to administrative region by using ArcGIS 10.2. The extent contour and central points for each county were also stored in the geographical information database. Data fields included name, centre point coordinates, county code and etc. Attribute database consists of HFRS cases data per month and per year for each county, demographics and climate data. Data fields for attributes data included county code, average temperature value, average humidity value, average rainfall value, cases number for each county and death number for each county.

Methods

Spatial autocorrelation analysis.

Spatial autocorrelation analysis is applied to describe the similarity of geographically proximate units. In this study, spatial autocorrelation analysis was used to clarify the overall autocorrelation pattern of HFRS incidence (clustered, dispersed, or random). Global Moran’s I statistic, developed by Patrick Alfred Pierce Moran [20], is one of the most commonly used methods in spatial autocorrelation analysis. It calculates the autocorrelation (I) as follows: (1) Where n is the total number of counties, xi is the number of cases for county i, Xj is the number of cases for county j. zi and zj are the deviation values for county i and county j from their mean and respectively. ωi,j is the spatial weight which used to quantify the spatial relationships between county i and county j. Spatial weight is structured into a matrix, the distance between county i and county j is calculated with K nearest neighbours.

Moran’s I index values for each year are calculated to measure spatial autocorrelation characteristics. Moran’s I index value ranges from −1.0 to +1.0. A zero Moran’s I Index value indicates a random spatial pattern [17]. Negative values indicate negative spatial autocorrelation while positive values indicate positive spatial autocorrelation.

Spatial autocorrelation analysis measures the globe spatial autocorrelation characteristics based on both county locations and number of HFRS cases simultaneously. It evaluates the spatial distribution pattern by a series number of county locations and the associated number values of HFRS cases.

Spatio-temporal scan statistical analysis.

It is helpful to identify the moving tendencies of HFRS incidence in Hubei Province by investigating its spatial and temporal clustering areas. After the first step of spatial autocorrelation analysis, if HFRS occurred in a clustered way (which can be treated as an epidemic), the spatio-temporal scan statistical analysis could be used to identify the spatial and temporal characteristics of this epidemic [5,11,15,21,22]. Kulldorff scan statistic method was widely used to measure the degree of clustering [2326]. This method imposed a cylindrical window to detect spatio-temporal clusters by constantly changing the size of the window, it has a circular (or elliptic) geographic base and a height corresponding to the disease onset time [14,20,27,28].

Using this method, the radius of the cylindrical window varies continuously from zero to a specified maximum size. The maximum-size specified the percentage of the maximum total population at risk within the scanning window [29]. In this research, considering the human population density and the size of each county, the maximum cluster size was set to 50%. Discrete Poisson Model was used to identify the cluster locations of HFRS cases in a specific window. Moreover, SaTScan program was used to implement the scan statistical analysis [2932].

Spearman's rank correlation analysis.

Previous studies have showed that the transmission of HFRS in Hubei/China was influenced by environmental factors such as temperature, humidity and rainfall [10,26,33,34]. Climate variables have some relationships with dynamic rodent population, which can be used as possible influencing factors for the prediction of the HFRS epidemic [12]. In order to better understand influencing factors of the epidemic HFRS incidence, average temperature, average humidity, total rainfall and human population density were studied in this research. Spearman's rank correlation coefficients were used to measure the associated degree between variables [35,36].

Results

Features of HFRS outbreaks in Hubei Province

First, the ratio of HFRS cases in each month over the whole year, we define it as the HFRS outbreak value. Fig 2 demonstrates the HFRS outbreak rates for each month from 2005 to 2014 in Hubei Province. We observed the value of the rate fluctuated month by month, and had the maximum value of 0.6534 (June 2006) and the minimum value of 0.05 (September 2012). In general, the simulated trend line shows that the HFRS incidence rate had a decreasing tendency. It indicates that HFRS incidence decreased in general while increasing slightly from 2012 to 2014.

thumbnail
Fig 2. Incidence and trends in Hubei Province from 2005 to 2014.

Each point represents the HFRS incidence rate in a specific year. All of the points are lined to indicate the trend of the HFRS incidence in Hubei Province. Trend line is simulated according to the incidence rates.

https://doi.org/10.1371/journal.pone.0167836.g002

Spatial distribution of HFRS cases

Fig 3 indicates the Moran’s I index values of the spatial autocorrelation results for HFRS cases from 2005 to 2014. Fig 3 and Table 1 indicate the spatial distribution of HFRS cases happened in 76 counties. Except 2005 and 2011, HFRS outbreak generally presented spatially clustered distribution with z-score>1.65 and p-value<0.05. The clustering degree of HFRS cases was gradually increased from 2006 to 2009 and reached its peak in 2010 (Moran’s I = 0.2374 with z-score = 4.1336 and p-value = 0.0001), while it happening in a high discrete pattern with the lowest point in 2011 (Moran’s I = 0.0550 with z-score = 1.2442 and p-value = 0.2134). Since 2012, it recovered to the clustering pattern till 2014 (Moran’s I = 0.2604 with z-score = 4.4143 and p-value = 0.0001).

thumbnail
Fig 3. Distribution pattern of HFRS according to Moran's Index from 2005–2014.

Each point represents the Moran’s I value in a specific year. All of the points are lined to indicate the trend of the Moran’s I in Hubei Province.

https://doi.org/10.1371/journal.pone.0167836.g003

thumbnail
Table 1. Spatial Autocorrelation results of average annual incidence rate of each year for HFRS in Hubei Province.

https://doi.org/10.1371/journal.pone.0167836.t001

Clustering pattern of HFRS cases

The spatio-temporal scan statistical analysis showed that HFRS was not randomly distributed in space during these ten years, which achieved the same conclusion as our previous studies did [37]. In Table 2, the number of clusters decreases from 5 to 2 from 2005 to 2009 while total cases of HFRS decreasing from 824 to 217. Although the variation from 2010 to 2014 is relatively small, the overall trend is still decreasing in terms of the spatio-temporal clusters of HFRS cases from 2005 to 2014.

thumbnail
Table 2. Spatio-temporal Clusters of HFRS cases from 2005–2014 in Hubei Province.

https://doi.org/10.1371/journal.pone.0167836.t002

Fig 4 illustrates the exact locations of the clusters area and clustering levels. The Most Likely Cluster locates in cluster one. The clustering level decreases from cluster one (in red) to cluster five (in blue). The clusters distributed in various counties from west to east in 2005 and 2007. From 2007 to 2012, the clusters gradually moved towards the central area of Hubei Province. In 2014, only one cluster existed in the central area.

thumbnail
Fig 4. Clusters of HFRS cases from 2005 to 2014.

Spatial and temporal clusters were identified each year from 2005 to 2014. Different colors represent different levels of clusters. The cluster one for each year was regarded as the Most Likely Cluster.

https://doi.org/10.1371/journal.pone.0167836.g004

Table 3 represents the Most Likely Clusters of HFRS cases from 2005–2014 in Hubei Province, from Table 3, we can find that HFRS cases mostly occurred in the second half of the year (from summer to autumn during 2005, 2006, 2007, 2009, 2012, 2013 and 2014). Further, the number of the counties included in cluster one varied from 2005 to 2014. For example, cluster one includes only one county (Huang Gang) in 2005 but 5 counties (Macheng, Luotian, Yingshan, Huanggang and Zuoshui) in 2006. Interesting, Yicheng County, as a special epidemic area, has been included in cluster one all the time from 2007 to 2014.

thumbnail
Table 3. The Most Likely Clusters (Cluster one in Fig 4) of HFRS cases from 2005–2014 in Hubei Province.

https://doi.org/10.1371/journal.pone.0167836.t003

Influencing factors of HFRS cases

Monthly HFRS outbreaks cases are presented in Fig 5 to explore the seasonal variation of HFRS outbreaks. Bimodal pattern is obviously observed noted that two-peak points appeared each year, one is in spring and the other is in autumn. In 2005, 2006 and 2007, one more peak points appeared during November.

thumbnail
Fig 5. Monthly HFRS cases from 2005 to 2014.

(A) Average monthly HRFS cases from 2005 to 2009. (B) Average monthly HRFS cases from 2009 to 2014.

https://doi.org/10.1371/journal.pone.0167836.g005

Related studies have demonstrated that climate factors had a great contribution to the outbreak of HFRS [3,7,10,12,38]. From Table 4, it seems that the average humidity influencing factor shows positive associated with HFRS cases (Spearman Correlation = 0.041 at the 0.001 level). Taking Yicheng County as an example, the value of its humidity is much higher than the average humidity level in the second half year (Fig 6). Combining the result of clusters in Fig 4 and the Most Likely Clusters in Table 3, results show that high humidity associate with the HFRS outbreaks consistently in Yicheng County.

thumbnail
Fig 6. Comparison of the humidity of Yicheng County with the average humidity of Hubei Province.

Average monthly humidity values for Hubei Province and Yicheng County are depicted by dotted line and solid line, respectively.

https://doi.org/10.1371/journal.pone.0167836.g006

thumbnail
Table 4. Correlation results for HFRS cases from 2005–2014 in Hubei Province.

https://doi.org/10.1371/journal.pone.0167836.t004

HFRS is an infectious disease. The related studies have indicated that the density of the virus carriers, human immunity and immigration played significant roles in the transmission of HFRS [10,39]. In Table 4, the correlation analysis result, which values are with Spearman Correlation = 0.397 and with p-value≤0.001, means that the human population density factor is positive associate with the outbreak of HFRS epidemic in Hubei Province.

Discussion

In this study, the relationships between the HFRS cases, spatio-temporal distribution and clustering patterns have been studied. We also explored some influencing factors for the HFRS cases in Hubei Province.

Our previous study of the past three decades of HFRS fatalities data revealed a general trend of decreased in Hubei Province [40]. No further outbreaks of HFRS occurred since 2000 after the two peaks happened in the 1980s and the in the 1990s, respectively [37]. In this study, we found the HFRS case increased from 2012 to 2014 while the general trend showed the HFRS cases decreased from 2005 to 2012. This phenomenon could possibly be explained the following thirdly factors.

First, Hantaviruses were carried and transmitted by rodents [39]. Xiao et al. found that rodent population density and the types of the rodent species significantly impact on the occurrence of HFRS [41]. An experiment from State Key Laboratory of Virology, Institute of Medical Virology, School of Medicine in Wuhan University, China has shown that the mice positive rate dropped from 9.16% to 0.28% from 2005 to 2011. But in 2012 the rate suddenly rose to 19%. The more frequently the aerosolized urine, droppings or saliva of infected rodents people touched, the higher infection possibility people will have. In addition, individuals working with living rodents could get a higher probability of infection. The effect of a rapidly growing rodent population on the environmental load of hantaviruses has been studied by Sauvage et al. [42] Their results demonstrated that the increasing of the rodent population and the concomitant in the environmental contamination is linked to an increasing risk of human infection. In this study, we inferred that the increase of HFRS incidence in Hubei Province from 2012 might be caused by the suddenly increasing of infected rodent population.

Secondly, human migration (from rural areas to cities) and urbanization have both reduced the risk of exposing to rodent excreta. As shown in Fig 7, the ratio of rural population over total population in Hubei Province is decreasing year-by-year. It may affect the incidence of HFRS as the decreased trend in the whole period.

thumbnail
Fig 7. The ratio of rural population over the total population of Hubei Province from 2005–2014.

Each point represents the ratio in a specific year. All of the points are lined to indicate the trend of the rural population ratio in Hubei Province from 2005 to 2014.

https://doi.org/10.1371/journal.pone.0167836.g007

Thirdly, Fig 8 shows the continuous rising of GDP in Hubei Province from 2005 to 2014. The growth of GDP improves housing conditions and workplace conditions. These will be beneficial to the reduce of rodent-to-human contact. According to the statistical data acquired from Hubei Provincial Bureau of Statistics, medical and health institutions in Hubei Province increased from 9983 to 36089 during 2005–2014.

thumbnail
Fig 8. GDP of Hubei Province from 2005–2014.

Each point represents the GDP value in a specific year.

https://doi.org/10.1371/journal.pone.0167836.g008

In the study, we also found that the clustering degree was increasing and gathering to the central area of Hubei Province. In addition, from geographic characteristics perspective, HFRS has not been endemic near abundant water resources until 2005. Combining with the topography, geomorphology and vegetation covered environments, this finding might explain by the following.

The previous studies have demonstrated that HFRS in Hubei Province was mainly affected by SEOV virus which was transmitted by Rattus norvegicus during 1980–2012 [37,39]. Hubei Province is surrounded on three sides (east, west and north) by mountains. Its low and flat middle part is a piece of incomplete basin slightly open towards south. Counties located in the central part possess the subtropical humid monsoon and monsoon climate. This can provide a better living environment for the survival and reproduction of the rodents (like Rattus norvegicus), and the Rattus norvegicus could easily carry SEOV virus.

The industrial development often means high-speed economy development and the concentration of population. As listed in Table 3, Yicheng County appeared in the Most Likely Clusters consistently. Industry and agriculture data, collected from Hubei Provincial Bureau of Statistics, indicated that the industrial production from 2007 to 2014 increased from 3.001 billion to 16.717 billion by 457% percentage in Yicheng County. While for the none cluster area like Shennongjia County, for example, the number of industrial production from 2007 to 2014 increased from 42.062 million to 75.322 billion by only 79% percentage. Compared with the high human population density of 192 per square kilometre for Yicheng County, the population density in Shennongjia County only covers 18 per square kilometre. We believe that the industrial development and human population density may affect the number of HFRS cases.

Peng et al. found that HFRS cases were distributed along large water systems (wetland) between 1983 and 1995, such as Yangtze river and Huai river in Anhui Province [26]. The surrounding water system is helpful to the cultivation of crops. It will also create a better environment for rodents to survive and thrive. In Fig 9, the blue area represents water system that mainly distributes in the east part. In Fig 4, the cluster gathering area is in the western region during 2005 and 2006. And then, a cluster gathering area mainly concentrated in the northern central area in the following years (2007–2014). As a result, the cluster aggregation of HFRS has not showed a significant correlation with the water system. For instance, Yicheng County has about 115 square kilometres water area, the average annual HFRS cases rate during this period is 8.11. However, Wuhan City, Honghu City and Danjiangkou City have a water area of 1040, 622 and 400 square kilometres respectively, which rank among the top three largest water area cities. The HFRS cases rate corresponding to the three cities are only 0.13,0.77 and 0. In Fig 10, elevation values from DEM data are presented. The central, southern and eastern region are mainly plain areas or small hilly areas with relatively low elevation values (their average elevation is between -142 to 277 meters). Combined with Fig 4, HFRS case clusters mostly gathered in the central and southern part. One possible explanation may be that HFRS outbreak more easily in plains rather than in mountains.

thumbnail
Fig 9. Water distribution in Hubei Province.

The colored areas indicate the water boundaries in Hubei Province.

https://doi.org/10.1371/journal.pone.0167836.g009

thumbnail
Fig 10. Digital Elevation Model(DEM) Data for Hubei Province.

The color gradient presents the elevation value of Hubei Province in thematic maps. It is grouped in 10 categories based on the elevation value for each location.

https://doi.org/10.1371/journal.pone.0167836.g010

As displayed in Fig 5, the number of HFRS cases shows that there are two peaks annually. One happens in spring time (March to May). It is corresponded to the epidemic peak of SEOV-type HFRS emerged in spring time. Another peak appears between September and November, the peak of HTNV-type HFRS was often emerged during these time period [39,43]. This phenomenon indicates that the causes of HFRS in Hubei Province might be associated with both types of virus (HTNV and SEOV). Compared with the peaks happened in each year, and data revealed from Fig 5, we could infer that if peaks mostly emerged in spring time during a year, it could be presumed that SEOV-type HFRS have more influences vice versa.

Pervious study has concluded that there are many environmental factors that impact the spatio-temporal dynamics of HFRS, such as the changes of precipitation, temperature, land use and vegetation community dynamics[44]. For example, temperatures may affect the dynamics and activities of rodent hosts as well as the infectivity of hantavirus [40]. In this research, it was demonstrated that the factors of human population density and average humidity are significantly related to HFRS epidemic in Hubei Province.

Li et al. found that the increase of human population density could make the transmission of the virus much easier [43]. As shown in Table 4, human population density factor is positively associated with the outbreaks of HFRS. It can be inferred that the congregation of population may provide more opportunities for the transmission of HFRS virus.

Climate factors have been long considered as important impact factors for the spreading of HFRS [10,12,26]. In our research, average humidity was positively associated with the outbreak of HFRS. Other environmental factors (such as average temperature, total rainfall) did not display significant relations with the outbreak of HFRS. The result was different from the other researches. One possible explanation is that we used data from different regions. For example, Liu et al. believed that rainfall, land surface temperature and relative humidity were significantly associated with monthly HFRS cases in Elunchun and Molidawahaner counties (China) for 1997–2007 [12]. However, Hubei Province is categorized as subtropical humid climate. It means the average humidity in Hubei Province is higher than the north part of China. Hubei Province is known as the "land of a thousand lakes". The climate and ecological conditions in Hubei Province provide an ideal habitat for rodents [45]. Xiao et al. demonstrated that moist environments are help for the vitality and infectivity of Hantaan virus and the distribution of rodents [41]. The humidity factor may affect several biological traits of the rodent hosts, from individual life histories to population dynamics. Corresponded with the correlation results from this research, it could be inferred that humid climate may contribute to the survival and reproduction for the virus. And this may contribute to the outbreaks of HFRS.

Conclusions

The findings in our research presented some new features regarding the trend and influencing factors for the HFRS outbreaks in Hubei Province. First, SEOV and HTNV are both associated the HFRS outbreaks. HFRS seasonal distribution in Hubei Province was characterized by a bimodal pattern (March to May, September to November). This bimodal pattern was corresponding to the peaks of SEOV and HTNV virus respectively. Compared with the peaks happened in each year, we could infer that if peaks mostly emerged in spring time during a year, it could be presumed that SEOV-type HFRS have more influences vice versa. Secondly, the average humidity and human population density were associated with the HFRS epidemic. It can be inferred that the congregation of population may provide more opportunities for the transmission of HFRS virus. And humid climate may contribute to the survival and reproduction for the virus. Thirdly, HFRS outbreaks were more in plains than in other areas of Hubei Province. HFRS case clusters mostly gathered in the central and southern part, the central and southern regions are mainly plain areas or small hilly areas. We did not find that whether the terrain of the wetland (water system) plays a significant role in the outbreak of HFRS incidence. We believe these findings can better understand the spatio-temporal pattern and influencing factors, and can further help to improve the reducing cases of HFRS in the future.

Supporting Information

S1 File. Population Data of Hubei Province in 2005

https://doi.org/10.1371/journal.pone.0167836.s001

(RAR)

S2 File. Map at a scale of 1:4000,000 of Hubei Province

https://doi.org/10.1371/journal.pone.0167836.s002

(RAR)

Acknowledgments

This work was supported by the National Natural Science Foundation of China, No. 71503068; Fundamental Research Funds in Key Research Areas for the Central Universities, No. 2015B09614; Fundamental Research Funds for the Central Universities, No. 2014B15114.

Author Contributions

  1. Conceptualization: L. Ge YZ.
  2. Data curation: L. Ge KZ.
  3. Formal analysis: L. Ge YZ.
  4. Funding acquisition: YZ.
  5. Investigation: L. Ge KZ YW.
  6. Methodology: L. Ge KZ.
  7. Project administration: L. Ge YZ.
  8. Resources: YZ XH.
  9. Software: L. Ge NW L. Guo.
  10. Supervision: L. Ge YZ.
  11. Validation: XH HY.
  12. Visualization: L. Ge YZ HF.
  13. Writing – original draft: L. Ge YZ XM.
  14. Writing – review & editing: L. Ge YZ.

References

  1. 1. Xu ZY, Guo CS, Wu YL, Zhang XW, Liu K. Epidemiological studies of hemorrhagic fever with renal syndrome: analysis of risk factors and mode of transmission. J Infect Dis. 1985; 1: 137–144.
  2. 2. Zou Y, Wang JB, Gaowa HS, Yao LS, Hu GW, Li MH, et al. Isolation and genetic characterization of hantaviruses carried by Microtus voles in China. J Med Virol. 2008; 4: 680–688.
  3. 3. Fang L, Wang X, Liang S, Li Y, Song S, Zhang W, et al. Spatio temporal trends and climatic factors of hemorrhagic fever with renal syndrome epidemic in Shandong Province, China. Plos Neglect Trop D. 2010; 8: e789.
  4. 4. Lin H, Liu Q, Guo J, Zhang J, Wang J, Chen H. Analysis of the geographic distribution of HFRS in Liaoning Province between 2000 and 2005. Bmc Public Health. 2007: 207. pmid:17697362
  5. 5. Wu W, Guo J, Guan P, Sun Y, Zhou B. Clusters of spatial, temporal, and space-time distribution of hemorrhagic fever with renal syndrome in Liaoning Province, Northeastern China. Bmc Infect Dis. 2011.
  6. 6. Zuo SQ, Fang LQ, Zhan L, Zhang PH, Jiang JF, Wang LP, et al. Geo-spatial hotspots of hemorrhagic fever with renal syndrome and genetic characterization of Seoul variants in Beijing, China. PLoS Negl Trop Dis. 2011; 1: e945.
  7. 7. Guan P, Huang D, He M, Shen T, Guo J, Zhou B. Investigating the effects of climatic variables and reservoir on the incidence of hemorrhagic fever with renal syndrome in Huludao City, China: a 17-year data analysis based on structure equation model. Bmc Infect Dis. 2009.
  8. 8. Zeier M, Handermann M, Bahr U, Rensch B, Muller S, Kehm R, et al. New ecological aspects of hantavirus infection: a change of a paradigm and a challenge of prevention—a review. Virus Genes. 2005; 2: 157–180.
  9. 9. Li Q, Yang ZQ, Qu H, Xiao H. Variability of G1 gene of hantaviruses occurring in the Hubei Province, P.R. China from 1985 to 2000. Acta Virol. 2005; 1: 51–57.
  10. 10. Bi P, Tong S, Donald K, Parton K, Ni J. Climatic, reservoir and occupational variables and the transmission of haemorrhagic fever with renal syndrome in China. Int J Epidemiol. 2002; 1: 189–193.
  11. 11. WU JQGZ W.. GIS-based spatial, temporal, and space–time analysis of haemorrhagic fever with renal syndrome.: Epidemiol. Infect; 2009.
  12. 12. Zhang WY, Guo WD, Fang LQ, Li CP, Bi P, Glass GE, et al. Climate variability and hemorrhagic fever with renal syndrome transmission in Northeastern China. Environ Health Perspect. 2010; 7: 915–920.
  13. 13. Liu Q, Liu X, Jiang B, Yang W. Forecasting incidence of hemorrhagic fever with renal syndrome in China using ARIMA model. Bmc Infect Dis. 2011.
  14. 14. Chen J RRNA. Geovisual analytics to enhance spatial scan statistic interpretations: an analysis of US cervical cancer mortality. Int J Health Geogr. 2008; 7.
  15. 15. Cui F, Wang T, Wang L, Yang S, Zhang L, Cao H, et al. Spatial analysis of hemorrhagic fever with renal syndrome in Zibo City, China, 2009–2012. Plos One. 2013; 6: e67490.
  16. 16. MORAN PA. Notes on continuous stochastic phenomena. Biometrika. 1950; 1–2: 17–23.
  17. 17. Wikipedia. Wikipedia Vandals. 2011; 3: 35.
  18. 18. Schmaljohn C, Hjelle B. Hantaviruses: a global disease problem. Emerg Infect Dis. 1997; 2: 95–104.
  19. 19. Liu J, Liu DY, Chen W, Li JL, Luo F, Li Q, et al. Genetic analysis of hantaviruses and their rodent hosts in central-south China. Virus Res. 2012; 2: 439–447.
  20. 20. Tiwari N, Adhikari CM, Tewari A, Kandpal V. Investigation of geo-spatial hotspots for the occurrence of tuberculosis in Almora district, India, using GIS and spatial scan statistic. Int J Health Geogr. 2006: 33. pmid:16901341
  21. 21. Fang L, Yan L, Liang S, de Vlas SJ, Feng D, Han X, et al. Spatial analysis of hemorrhagic fever with renal syndrome in China. Bmc Infect Dis. 2006: 77. pmid:16638156
  22. 22. Sugumaran R, Larson SR, Degroote JP. Spatio-temporal cluster analysis of county-based human West Nile virus incidence in the continental United States. Int J Health Geogr. 2009: 43. pmid:19594928
  23. 23. Zheng L YHLB. Epidemic characteristics and spatio-temporal patterns of scrub typhus during 2006–2013 in Tai'an, Northern China. Epidemiol Infect. 2015: 2451–2458. pmid:25543665
  24. 24. Ge E, Zhang X, Wang X, Wei X. Spatial and temporal analysis of tuberculosis in Zhejiang Province, China, 2009–2012. Infect Dis Poverty. 2016; 1: 11.
  25. 25. Zhang WY, Wang LY, Zhang XS, Han ZH, Hu WB, Qian Q, et al. Spatiotemporal Clustering Analysis and Risk Assessments of Human Cutaneous Anthrax in China, 2005–2012. Plos One. 2015; 7: e133736.
  26. 26. Bi P, Wu X, Zhang F, Parton KA, Tong S. Seasonal rainfall variability, the incidence of hemorrhagic fever with renal syndrome, and prediction of the disease in low-lying areas of China. Am J Epidemiol. 1998; 3: 276–281.
  27. 27. Li XZ, Wang JF, Yang WZ, Li ZJ, Lai SJ. A spatial scan statistic for nonisotropic two-level risk cluster. Stat Med. 2012; 2: 177–187.
  28. 28. 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 Health Geogr. 2004; 1: 11.
  29. 29. Azage M, Kumie A, Worku A, Bagtzoglou AC. Childhood Diarrhea Exhibits Spatiotemporal Variation in Northwest Ethiopia: A SaTScan Spatial Statistical Analysis. Plos One. 2015; 12: e144690.
  30. 30. Lefebvre A, Bertrand X, Vanhems P, Lucet JC, Chavanet P, Astruc K, et al. Detection of Temporal Clusters of Healthcare-Associated Infections or Colonizations with Pseudomonas aeruginosa in Two Hospitals: Comparison of SaTScan and WHONET Software Packages. Plos One. 2015; 10: e139920.
  31. 31. Alvarez J, Goede D, Morrison R, Perez A. Spatial and temporal epidemiology of porcine epidemic diarrhea (PED) in the Midwest and Southeast regions of the United States. Prev Vet Med. 2016: 155–160. pmid:26586344
  32. 32. SaTScan software for the spatial and space-time scan statistic [Type of Medium]. City;. Available from: http://www.satscan.org/.
  33. 33. Landscape elements and Hantaan virus–related hemorrhagic fever with renal syndrome, People's Republic of China.
  34. 34. WY Z, LQ F, JF J. Predicting the Risk of Hantavirus Infection in Beijing, People's Republic of China. The American journal of tropical medicine and hygiene. 2009: 678–683. pmid:19346399
  35. 35. Myers JL. Research Design and Statistical Analysis. In: Well AD, editor: Lawrence Erlbaum; 2003. pp. 508.
  36. 36. Puth M, Neuh User M, Ruxton GD. Effective use of Spearman's and Kendall's correlation coefficients for?association between two measured traits. Anim Behav. 2015: 77–84.
  37. 37. Zhang YH, Ge L, Liu L, Huo XX, Xiong HR, Liu YY, et al. The epidemic characteristics and changing trend of hemorrhagic fever with renal syndrome in Hubei Province, China. Plos One. 2014; 3: e92700.
  38. 38. Thomson MC, Garcia-Herrera R, Beniston M, Kelly-Hope L, Thomson MC. Climate and Infectious Diseases. In: Thomson MC, Garcia-Herrera R, Beniston M, ^editors. Seasonal Forecasts, Climatic Change and Human Health: Springer Netherlands; 2008. p. 31–70.
  39. 39. Zhang S, Wang S, Yin W, Liang M, Li J, Zhang Q, et al. Epidemic characteristics of hemorrhagic fever with renal syndrome in China, 2006–2012. Bmc Infect Dis. 2014; 1: 384.
  40. 40. Hardestam J, Simon M, Hedlund KO, Vaheri A, Klingstrom J, Lundkvist A. Ex vivo stability of the rodent-borne Hantaan virus in comparison to that of arthropod-borne members of the Bunyaviridae family. Appl Environ Microbiol. 2007; 8: 2547–2551.
  41. 41. Xiao H, Lin X, Gao L, Huang C, Tian H, Li N, et al. Ecology and geography of hemorrhagic fever with renal syndrome in Changsha, China. Bmc Infect Dis. 2013: 305. pmid:23819824
  42. 42. Piechotowski I, Brockmann SO, Schwarz C, Winter CH, Ranft U, Pfaff G. Emergence of hantavirus in South Germany: rodents, climate and human infections. Parasitol Res. 2008; 1 Supplement: 131–137.
  43. 43. Li S, Ren H, Hu W, Lu L, Xu X, Zhuang D, et al. Spatio temporal Heterogeneity Analysis of Hemorrhagic Fever with Renal Syndrome in China Using Geographically Weighted Regression Models. Int J Env Res Pub He. 2014.
  44. 44. Fang L, Zhao W, De Vlas SJ, Zhang W, Liang S, Looma CW, et al. Spatiotemporal dynamics of hemorrhagic fever with renal syndrome, Beijing, People's Republic of China. Emerg Infect Dis. 2009; 12: 2043–2045.
  45. 45. Du H, Wang PZ, Li J, Bai L, Li H. Clinical characteristics and outcomes in critical patients with hemorrhagic fever with renal syndrome. Bmc Infect Dis. 2014.