INTRODUCTION
Haemorrhagic fever with renal syndrome (HFRS) is a zoonosis caused by different species of Hantavirus, characterized by high fever and varying degrees of renal insufficiency and haemorrhage. It has been reported throughout Europe and Asia with 90% of cases occurring in China, where HFRS is an extremely dangerous viral disease second only to viral hepatitis [Reference Shi1]. The annual incidence of HFRS in China peaked at 11·06 new cases/100 000 persons in 1986, declined in the late 1990s, but has rebounded in recent years. It remains a public health problem with 20 000–50 000 human cases annually in mainland China [Reference Lei2, Reference Song3]. HFRS incidence distinctly shows spatial clustering and temporal variability [Reference Fang4, Reference Lin5]. For example, annualized average incidence at the county-level ranged from 0 to 122·117/100 000 during 1994–1998 [Reference Fang4]. In the past decades, HFRS has occurred in about 10-year endemic cycles in east and northeast China [Reference Song6].
Various rodent species are natural reservoir hosts and serve as sources of hantavirus infection. In China, HFRS is predominantly caused by Hantaan virus and Seoul virus, for which Apodemus agrarius and Rattus norvegicus are the leading rodent hosts, respectively [Reference Lei2]. Transmission of hantaviruses from rodents to humans is usually through inhalation of aerosols generated by contaminated excreta, saliva, or urine and possibly through contaminated food or rodent bites [Reference Shi1, Reference Song3, Reference Clement7]. Factors related to rodent habitats, breeding and survival – such as natural environment including elevation, soil type, vegetation, climate change and human behaviours including farming, cultivating and migration – can influence the transmission of hantaviruses and the epidemic process of HFRS [Reference Viel8]. Large populations of rodents can potentially increase the transmission of hantaviruses to humans [Reference Engelthaler9, Reference Hjelle and Glass10]. Excess proliferation of rodent populations related to climate change is of considerable international public health concern [Reference Kausrud11, Reference Tersago12]. Previous studies have shown that HFRS epidemics are sensitive to climatic conditions [Reference Engelthaler9, Reference Tersago12, Reference Zhang13]. However, the quantitative relationship between climate variation and HFRS incidence remains to be determined, especially in high-incidence areas.
A case-crossover design was proposed as an alternative to spatial-ecological modelling to make causal inferences in environmental epidemiology [Reference Carracedo-Martinez14]. It is a method for studying the effect of a risk factor on a health event using only cases with the idea of comparing a case's exposure immediately prior to or during the case-defining event with that same person's exposure at otherwise similar ‘reference’ times [Reference Maclure15]. Each case's exposure values comprise a matched set so that conditional logistic regression is typically used to estimate an odds ratio as the measure of association [Reference Bateson and Schwartz16]. Case-crossover design has the advantage of controlling time-invariant factors because it only involves cases and each case is compared to him/herself [Reference Nitta17].
This study was aimed at exploring the association of atmosphere temperature and precipitation with HFRS incidence in Junan County of Shandong Province, a region with a high prevalence of HFRS, by the method of case-crossover analysis.
MATERIALS AND METHODS
Study site
Junan County is located in the southeast of Shandong Province and on the border with Jiangsu Province, and has a population of about 970 000 people and a total area of 1752 km2. It has a warm temperate and semi-humid continental climate with four distinct seasons. The average altitude of this county is 200 m (range 19·9–662·2 m) above sea level. The county reported the highest epidemic strength of HFRS in Shandong Province with an annual average incidence reaching 55·20/100 000 over the past decades and therefore it was selected as one of the 39 national monitoring spots of HFRS.
The data of the HFRS cases for this study were from the National Notifiable Disease Surveillance System. Monthly climatic data of this county were obtained from the Meteorological Information Center of Shandong Province.
HFRS surveillance in Junan County
The HFRS surveillance system has been established in Junan County since 1975 to examine the general situation of this disease, to monitor its temporal trend, and to assess the effectiveness of interventions. HFRS case diagnosis and reporting follow the national protocols proposed by the Ministry of Health [Reference Song18]. Briefly, healthcare providers and clinics are required to report clinically diagnosed/suspected cases of HFRS to the county's Center for Disease Control and Prevention. Upon receipt of a case report, trained health professionals check the report and review medical records to confirm the clinical diagnosis. Blood samples from the case patient might be collected and a positive serum test of specific antibody to Hantavirus was used to confirm the disease by identifying the Hantavirus type.
Case-crossover design
The latent period from hantavirus infection to the onset of clinical symptoms ranges from 4 to 45 days, usually 2–4 weeks [Reference Song18]. It usually takes about 3 months for rodents, the reservoirs of hantaviruses, to reach the adult stage of life from birth [Reference Shi and Zhang19, Reference Liu20]. Therefore, we assumed that for a case the hazard period was during the most recent 4 months, just before the onset of disease. We defined an exposure window as a time interval in which the climatic condition is assessed as a potential risk factor for HFRS. In accordance with the available climatic data for calendar month, we defined four case exposure windows: the calendar month when the case was diagnosed (exposure window 0), and each of the three calendar months preceding the month when the case was diagnosed (exposure windows -1, -2, -3, respectively). For example, for a HFRS case diagnosed in November 1985, the case exposure windows -3, -2, -1 and 0 would be August, September, October and November 1985, respectively. Symmetric bidirectional case-crossover design [Reference Bateson and Schwartz16] was used in this study in order to control for confounding due to possible secular trend that was not related to the climatic conditions of interest. Each case exposure window was matched with two referent exposure windows symmetric to the case exposure window, one a year before and the other a year after the case exposure window. For example, if the case exposure window was November 1985, then the referent exposure windows were November 1984 and November 1986. A case exposure window and the corresponding two referent exposure windows resulted in a 1:2 match set analogous to a case-control study.
Data analyses
Because the effect of climatic variables on the incidence of diseases is not linear [Reference Hu21], air temperature and precipitation were transformed to categorical variables. The mean temperature for each month was grouped into three levels: low (<10 °C), moderate (10–25 °C), and high (>25 °C). The monthly amount of precipitation was grouped into <10 mm, 10–120 mm, and >120 mm. Such categorization of temperature and precipitation was based on the reports of the climatic conditions thought to be suitable for rodents' breeding [Reference Zhang22, Reference Li, Pu and Mu23].
The climatic conditions of the exposure windows were compared with the climatic conditions of the referent windows with the application of conditional logistic regression. The odds ratio (OR) of HFRS and its 95% confidence interval (CI) for a climatic condition compared to the referent condition were estimated. All the data analyses were performed using SAS v. 9.13 (SAS Institute Inc., USA).
RESULTS
There were 13 722 clinically confirmed cases of HFRS reported to Junan Center for Disease Control and Prevention from January 1977 to December 2001. The yearly incidence varied from 1·59 to 185·93/100 000 with a mean of 61·72/100 000. The monthly incidence during the period is shown in Figure 1.
(a) Monthly haemorrhagic fever with renal syndrome (HFRS) incidence (1/100 000) in Junan County from January 1977 to December 2001. (b) Annual HFRS incidence rate (1/100 000) in Junan County from 1977 to 2001.
According to the meteorological data, the annual average temperature and accumulative precipitation for the period from 1976 to 2002 were 13·24 °C and 787·99 mm, respectively. January had the lowest temperature with a mean temperature of −0·99 °C and July was the warmest month with a mean monthly temperature of 25·79 °C. Average monthly precipitation from 1976 to 2002 varied from the lowest of 10·35 mm in December to the highest of 226·32 mm in July (Table 1).
Mean temperature and precipitation by calendar month from 1976 to 2002 in Junan County, China
Atmosphere temperature and HFRS risk
The relationship between temperature and HFRS risk is shown in Table 2. During the previous 3 months of a HFRS case being diagnosed, a monthly mean temperature level >25 °C, compared to the monthly mean temperature of 10–25 °C, was significantly associated with a lower risk of HFRS at any exposure window consistently. The estimations of ORs were 0·802 (95% CI 0·672–0·958), 0·289 (95% CI 0·233–0·358), and 0·662 (95% CI 0·563–0·779) for exposure windows 0, -1 and -2, respectively. In contrast, there is no consistent relationship of a temperature <10 °C with HFRS risk over the different exposure windows. The cold weather with mean temperature <10 °C only at exposure window -2 (2 months before the occurrence of HFRS cases) significantly decreased the number of cases (OR 0·339, 95% CI 0·183–0·629). The mean temperature of exposure window -3 (3 months before the occurrence of HFRS cases) was not found to be associated with HFRS incidence (P = 0·4428 and 0·1421 for cold weather and hot weather, respectively).
Monthly mean temperature levels at different exposure windows and HFRS risk: results of conditional logistic regression
HFRS, Haemorrhagic fever with renal syndrome; OR, odds ratio; CI, confidence interval.
* Exposure windows 0, -1, -2, -3 represent the calendar month when the case was diagnosed and the three calendar months preceding the month when the case was diagnosed.
† The same calendar month of the year before the hazard period.
‡ The same calendar month of the year after the hazard period.
Amount of precipitation and HFRS risk
The results of conditional logistic regressions of HFRS regarding the amount of precipitation assessed by the different exposure windows are shown in Table 3. Lack of precipitation (monthly precipitation <10 mm) and excessive precipitation (monthly precipitation >120 mm) in the month that the case occurred were both associated with lower risk of HFRS compared to a moderate amount of precipitation (monthly precipitation between 10 and 120 mm). The estimated ORs were 0·821 (95% CI 0·774–0·871) and 0·860 (95% CI 0·774–0·955), respectively. Other significant associations included excessive precipitation at exposure windows -2 and -3 with increasing risk of HFRS (OR 1·355, 95% CI 1·210–1·472 and OR 1·217, 95% CI 1·108–1·336, respectively), which means that plentiful rain during the preceding 2–3 months would be a risk factor for an HFRS epidemic.
Monthly precipitation levels at different exposure windows and HFRS risk: results of conditional logistic regression
HFRS, Haemorrhagic fever with renal syndrome; OR, odds ratio; CI, confidence interval.
* Exposure windows 0, -1, -2, -3 represent the calendar month when the case was diagnosed and the three calendar months preceding the month when the case was diagnosed.
† The same calendar month of the year before the hazard period.
‡ The same calendar month of the year after the hazard period.
Combinations of temperature and precipitation categories and HFRS risk
In addition to the independent analysis for temperature and precipitation (Tables 2 and 3), we analysed how the weather conditions of cross-classification by temperature and precipitation affected the incidence of HFRS using conditional logistic regressions. There were eight cross-groups, excluding the hot and dry weather condition (mean temperature >25 °C and monthly precipitation <10 mm) classified since no month during the study period had hot and dry weather. The condition with appropriate temperature (10–25 °C) and moderate precipitation (10–120 mm) was treated as the referent group, odds ratios and their confidence intervals are displayed for other conditions of the four respective exposure windows (see Table 4). This shows that compared to the referent climatic condition, any combination with mean temperature and monthly precipitation beyond its appropriate range in the preceding 2 months (exposure windows 0 and -1) tended to be associated with a lower risk of HFRS (all ORs <1), which means that appropriate temperature and moderate precipitation is the most favourable condition for occurrence of HFRS. On the contrary, abundant precipitation combined with appropriate air temperature 2 or 3 months before was a risk factor for HFRS. The estimated ORs of precipitation >120 mm with mean temperature between 10 °C and 25 °C were 1·193 (95% CI 1·063–1·339) and 1·348 (95% CI 1·193–1·524) at exposure windows -2 and -3, respectively. The estimated OR of precipitation >120 mm with mean temperature >25 °C was 1·171 (95% CI 1·005–1·334).
Associations of weather conditions and HFRS incidence: results of conditional logistic regression
HFRS, Haemorrhagic fever with renal syndrome; OR, odds ratio; CI, confidence interval.
Bold values indicate the higher risk weather condition for HFRS incidence compared to the moderate weather condition (temperature of 10–25 °C combined with precipitation of 10–20 mm).
DISCUSSION
HFRS is epidemic in many provinces with high incidence in mainland China, although integrated intervention measures including rodent control, environment management and vaccination have been implemented for more than 10 years [Reference Zhang24]. Geographical distribution of HFRS cases in mainland China show that HFRS mainly occurs in areas where the annual average temperature is <20 °C, the annual average relative humidity is 50–80%, and the annual total rainfall is 400–1600 mm [Reference Fang25]. Although the endemic areas are mainly in zones with mild temperature, HFRS is epidemic from the northernmost provinces to the southernmost provinces in east China. In contrast, there was no case reported in west China with an average altitude >2000 m, where precipitation was so rare that the relative humidity was usually <50% and rainy days were fewer than 60 each year [Reference Bi26]. Therefore, most researchers believe that precipitation, followed by air temperature, are the most important elements among all the influencing meteorological factors [Reference Chen27]. This is why we focused on these two variables in order to investigate the effects of meteorological factors on the occurrence of HFRS. Although humidity was also found to be associated with HFRS incidence in previous studies, it was not included in the present study because of the strong correlation between monthly mean humidity and monthly accumulative precipitation during the study period (r = 0·70, n = 324, P < 0·0001).
The results of our study show that the risk of HFRS cases increasing is positively associated with the amount of precipitation in the preceding 2–3 months, while moderate precipitation during the preceding 1–2 months might be necessary to facilitate the prevalence of HFRS. Plenty of rainfall is beneficial for growth of vegetation, which provides abundant food for rodents. Meanwhile, atmospheric moisture makes the food soft which encourages the rodents to consume more [Reference Chen27]. On the contrary, drought usually results in the number of rodents decreasing, which has a detrimental effect on the prevalence of HFRS [Reference Luo28]. For example, a severe drought lasting for three seasons from October 1994 in Xi'an City, China resulted in an obvious decline in HFRS incidence by 1995, falling to the lowest level in 1996 [Reference Yang, Xia and Wu29]. It is natural that the effect on variation of HFRS incidence lags behind climate change since the growth of vegetation, rodents' breeding and virus transmission are a series of processes. The lag period of rainfall for HFRS was about 2–3 months as shown by the results of this study. Nevertheless, excessive rain might work against an HFRS epidemic probably through curtailing the activities of rodents and humans, thereby reducing their contact opportunities. This is supported by our result that accumulative precipitation >120 mm in the same month, compared to moderate precipitation of 10–120 mm, was a disadvantage for occurrence of HFRS.
Regarding atmosphere temperature, our results show that extremes of weather (too cold or too hot) do not favour HFRS prevalence, while the most appropriate mean temperature was shown to be between 10 °C and 25 °C. It was reported that both Apodemus agrarius and Rattus norvegicus have their highest breeding rate at temperatures between 10 °C and 25 °C, while the breeding rate decreased with temperatures <10 °C or >25 °C [Reference Zhang22, Reference Li, Pu and Mu23, Reference Yang and Guo30].
The results of the combinations of temperature and precipitation suggest that the study factors appear to be associated with higher numbers of HFRS cases, which will be valuable for predicting increasing incidence. For example, abundance of precipitation >120 mm or moderate precipitation and hot weather with average temperature >25 °C might indicate an epidemic of HFRS in 3 months. Therefore, taking preventive measures in advance, e.g. deratization, avoiding contaminated food or water, intensifying individual protection would be necessary if such climatic conditions occurred.
In this study, a case-crossover design was used to analyse the effect of mean air temperature and precipitation on the risk of HFRS onset. This has advantages in controlling time-invariant confounding factors, e.g. geographical and geological factors, socioeconomic factors, and people's lifestyle [Reference Holly, Lianne and Thomas31]. Moreover, the time lag of the effect was effectively detected through setting different matching groups month by month and the nonlinear problem was avoided by defining dichotomous dummy variables in conditional logistic regressions. Compared to complex methods such as a generalized additive Poisson regression model in time-series analyses, the data analysis based on the case-crossover design is much easier to perform and the results are more explicable for epidemiologists.
Nevertheless, there are still limitations of this study. First, the long-term effect of climate on the epidemic trend of HFRS was not investigated since the case-crossover design is used to estimate the risk of a rare event associated with a short-term exposure. Second, the potential interaction effect of other meteorological variables, e.g. humidity, duration of sunshine, and atmospheric pressure was not considered. Therefore, further investigation of modelling the associations between climate and occurrence of HFRS is warranted.
CONCLUSIONS
Atmospheric temperature and amount of precipitation during the preceding 3 months can impact the endemic intensity of HFRS. The climatic conditions with relatively high risk for a HFRS epidemic included: condition with moderate mean air temperature (10–25 °C) and abundant precipitation (>120 mm) 3 months before (OR 1·348, 95% CI 1·193–1·524) or 2 months before (OR 1·193, 95% CI 1·063–1·339); condition with temperature >25 °C and abundant precipitation (>120 mm) 3 months before (OR 1·171, 95% CI 1·005–1·364); condition with appropriate temperature (10–25 °C) and moderate precipitation (10–120 mm) during the same month. A case-crossover study is an alternative method for analysing the time-series data of meteorological epidemiology.
ACKNOWLEDGEMENTS
We thank the Center for Disease Control of Junan County and the Meteorology Information Center of Shandong Province for sharing with us the data needed for this study. This study was supported by the Shandong Outstanding Young and Middle-aged Scientists Research Award Fund (2007BSB02046) and the National Basic Research Program of China (973 Program) (2012CB955500).
DECLARATION OF INTEREST
None.
