Skip to main content

ORIGINAL RESEARCH article

Front. Public Health, 11 January 2024
Sec. Infectious Diseases: Epidemiology and Prevention

Epidemiological characteristics and prediction model construction of hemorrhagic fever with renal syndrome in Quzhou City, China, 2005–2022

Qing Gao&#x;Qing Gao1Shuangqing Wang&#x;Shuangqing Wang2Qi WangQi Wang1Guoping CaoGuoping Cao2Chunfu FangChunfu Fang2Bingdong Zhan,
Bingdong Zhan1,2*
  • 1School of Public Health, Zhejiang Chinese Medical University, Hangzhou, Zhejiang, China
  • 2Quzhou Center for Disease Control and Prevention, Quzhou, Zhejiang, China

Background: Hemorrhagic fever with renal syndrome (HFRS) is one of the 10 major infectious diseases that jeopardize human health and is distributed in more than 30 countries around the world. China is the country with the highest number of reported HFRS cases worldwide, accounting for 90% of global cases. The incidence level of HFRS in Quzhou is at the forefront of Zhejiang Province, and there is no specific treatment for it yet. Therefore, it is crucial to grasp the epidemiological characteristics of HFRS in Quzhou and establish a prediction model for HFRS to lay the foundation for early warning of HFRS.

Methods: Descriptive epidemiological methods were used to analyze the epidemic characteristics of HFRS, the incidence map was drawn by ArcGIS software, the Seasonal AutoRegressive Integrated Moving Average (SARIMA) and Prophet model were established by R software. Then, root mean square error (RMSE) and mean absolute error (MAE) were used to evaluate the fitting and prediction performances of the model.

Results: A total of 843 HFRS cases were reported in Quzhou City from 2005 to 2022, with the highest annual incidence rate in 2007 (3.93/100,000) and the lowest in 2022 (1.05/100,000) (P trend<0.001). The incidence is distributed in a seasonal double-peak distribution, with the first peak from October to January and the second peak from May to July. The incidence rate in males (2.87/100,000) was significantly higher than in females (1.32/100,000). Farmers had the highest number of cases, accounting for 79.95% of the total number of cases. The incidence is high in the northwest of Quzhou City, with cases concentrated on cultivated land and artificial land. The RMSE and MAE values of the Prophet model are smaller than those of the SARIMA (1,0,1) (2,1,0)12 model.

Conclusion: From 2005 to 2022, the incidence of HFRS in Quzhou City showed an overall downward trend, but the epidemic in high-incidence areas was still serious. In the future, the dynamics of HFRS outbreaks and host animal surveillance should be continuously strengthened in combination with the Prophet model. During the peak season, HFRS vaccination and health education are promoted with farmers as the key groups.

1 Introduction

Hantavirus (HV) is a single-stranded enveloped RNA virus in the Bunyaviridae family. It primarily causes two infectious diseases in humans, namely Hantavirus Cardiopulmonary Syndrome (HCPS) in the Americas and hemorrhagic fever with renal syndrome (HFRS) in Europe and Asia (1). The morbidity and mortality rates vary from 0 to 43% among various pathogenic virus genotypes (2). Global HV infections show a fluctuating upward trend (3). HFRS is a naturally zoonotic disease, with rodents and other small mammals as the main source of transmission and natural hosts, the Apodemus agrarius and Rattus norvegicus, which carry HVs known as hantaan virus (HTVN) and Seoul virus (SEOV), respectively. The virus is mainly transmitted by aerosol or direct contact (4). The main clinical manifestations include kidney damage, bleeding, and fever (5). At present, HFRS occurs in more than 30 countries in the world, such as China, the United States, Russia, etc. (68). It can report 150,000–200,000 cases annually and is a global public health problem (9). According to statistics, China has the largest number of HFRS cases in the world, with an annual report of 20,000–50,000 cases, accounting for 90% of the global cases (10, 11). It poses a serious threat to the health of residents and socio-economic development and has become one of the major public health problems (12).

HFRS cases were first detected in Zhejiang Province in 1963, and outbreaks occurred in many areas of southern China in the 1980s. Zhejiang Province has a high incidence of HFRS in southern provinces, China. In the past decade, outbreaks in Zhejiang Province have continued to be sporadic. Quzhou, located in the western part of Zhejiang Province, has always been a high-incidence area, and Kaihua County ranks among the top five in the province in terms of cumulative cases from 2005 to 2020 (13). Therefore, it is of great significance to establish an epidemic prediction model for HFRS, explore its epidemic characteristics and development patterns, and seek optimal prevention and control strategies.

The prevalence of HFRS has been found to be influenced by various aspects, such as the accelerated progress of urbanization, and changes in land use types may also increase the risk of HFRS (14). Different populations have different probabilities of contact with rodents, leading to different risks of infection (15). Meanwhile research evidence suggests that a significant decline in the incidence of vector-borne infections occurred as a result of the 2020 policy of non-pharmacological interventions for novel coronavirus pneumonia. And there is no specific treatment for HFRS (16). In recent years, many scholars have used time series analyses to study the temporal changes of HFRS in order to predict future trends, such as the product seasonal autoregressive moving average model (SARIMA) (17). The model has a simple structure, requires only morbidity information for prediction, is highly feasible, and has been widely used in healthcare (18). The Prophet model is a procedure for predicting time series data based on an additive model, which can simulate multiple seasons simultaneously by means of a generalized additive model, and the fitting effect has been validated in time series prediction of infectious diseases such as Hand, Foot and Mouth Disease (HFMD) and AIDS (19, 20), but the effect of the fit in HFRS is not yet conclusive.

In order to better understand the epidemiological characteristics of HFRS in Quzhou, improve the prevention and control measures, and optimize the allocation of resources. We conducted a descriptive analysis on the incident data of HFRS from 2005 to 2022. Using GIS technology, combining with land use data, we analyzed the spatial distribution characteristics of HFRS in Quzhou. The SARIMA and Prophet models were used to establish a prediction model for the reported incidence of HFRS in Quzhou City from 2005 to 2022, to compare the prediction effects of each model, to explore the optimal model, and to provide scientific suggestions for prediction and early warning in the prevention and control of HFRS in Quzhou City.

2 Methods

2.1 Setting

Quzhou is located in the western region of Zhejiang Province, China, the upper reaches of the Qiantang River, between 118°01′15″–119°20′20″ east longitude and 28°15′26″–29°30′00″ north latitude. The city’s terrain is dominated by mountainous hills, which is a subtropical monsoon climate, with four distinct seasons throughout the year, abundant sunshine and precipitation. By 2023, the city has jurisdiction over 2 districts, 3 counties, and 1 county-level city.

2.2 Data sources and collection

Data of confirmed HFRS cases in Quzhou from January 2005 to December 2022 were obtained from the National Infectious Disease Reporting System (NIDRS) and the Chinese Center for Disease Control and Prevention (CDC). Cases were reported by medical and health institutions at all levels in Quzhou and reviewed by CDC at all levels. All cases were confirmed according to the HFRS diagnostic criteria of the Ministry of Health of the People’s Republic of China (WS 278–2008). Confirmed cases must have a history of traveling to an infected area or a history of direct or indirect contact or suspected contact with rats or their excrement (stool, urine) components within 2 months prior to the onset of disease, and have at least two of the following clinical characteristics: fever, chills, bleeding, headache, back pain, abdominal pain, acute renal insufficiency, and hypotension. In addition, a laboratory diagnostic criterion must be met: a positive result for hantavirus-specific immunoglobulin M, a 4-fold rise in titers of hantavirus-specific immunoglobulin G, a positive result for hantavirus-specific ribonucleic acid by reverse transcription polymerase chain reaction in clinical specimens, or having hantavirus isolated from clinical specimens. Land use data from GlobeLand30 2020 edition,1 which was developed by the spatial resolution of 30 meters surface coverage data around the world, it divides land into 10 first-level types: arable land, Woodland, grassland, shrub land, wetland, water body, tundra, artificial surface, unused land, glacier and permanent snow cover. Population data were obtained from Quzhou Statistical Yearbook. The urban (county) boundary maps and township boundary maps of Quzhou are provided by the CDC.

2.3 Statistical analysis

2.3.1 Descriptive analysis

Microsoft Excel (Microsoft Corporation, Redmond, WA, United States) was used to generate a database, organize and analyze the cases, and calculate the incidence rate. R software (version 4.2.1, R Foundation for Statistical Computing, Vienna, Austria) was used to analyze the prevalent characteristics of HFRS (time distribution, age distribution, sex distribution and occupation distribution) in Quzhou city from 2005 to 2022. Qualitative data were described as frequency and percent, and χ2-test was used for inter-group comparison, with p < 0.05 indicating a statistically significant difference.

2.3.2 Spatiotemporal analysis

The incidence data of Quzhou City were divided by township/street, and the incidence rate was calculated. The annual average incidence map and cumulative incidence map from 2005 to 2022 were established on the township/street scale using GIS technology in ArcGIS software (version 10.8; ESRI Inc., Redlands, CA, United States), and the spatial distribution of HFRS in Quzhou City was dynamically displayed. The land use of Quzhou City was collected by GlobeLand30, and the incidence data were linked to draw the HFRS incidence maps of different land use types.

2.3.3 Model construction

Autoregrestic moving average model (ARIMA) is a classical time series prediction method proposed by Geogre Box and Gwilym Jenkins, and has been widely used in the field of public health (2123). ARIMA (p, d, q) model is a hybrid model consisting of auto regressive (AR) model and moving average (MA) model, p is the autoregressive order, q is the moving average order, and d is the difference order. On the basis of ARIMA model, SARIMA adds seasonal parameters, P is the order of seasonal autoregression, D is the order of seasonal difference, Q is the order of seasonal moving average, and S is the time length of seasonal cycle (24). The general form of the model is (p, d, q) × (P, D, Q) S (25), which is suitable for both trend and seasonality. The model construction is divided into four steps: 1. Stationarity test: The application of SARIMA model requires time series to meet the requirements of stationarity, and the stationarity of data is preliminarily observed according to the incidence time series diagram; Augmented Dickey-Fuller (ADF) test is carried out. If p > 0.05 and the sequence is not stable, the original sequence needs to be stabilized by difference and/or seasonal difference and variable transformation. 2. Model identification order: According to the characteristics of the auto correlation function (ACF) graph and partial auto correlation function (PACF) graph, the order is preliminatively determined and multiple models are compared. The optimal model is selected according to Akaike’s information criterion (AIC). 3. Model parameter estimation and diagnosis: parameter estimation and hypothesis test are performed on the identified model, and white noise (Ljung-box, LB) test is performed on the residual sequence, and p > 0.05 passes the test, indicating that it could be used for prediction. 4. Model fitting prediction.

Prophet model is a time series model developed by Facebook in 2017 based on the C++ programming language, which can be used in Python and R (26). The model can handle flexibly the missing values of data, as well as the effects of holidays and special events in time series. The fitting speed of the model is fast, and it is suitable for predicting time series with strong seasonal influence and historical data of multiple seasons (19). The general form of Prophet model is y(t) = g(t) + s(t) + h(t) + ε, g(t) is the trend function, s(t) is the periodic function, h(t) is the holiday function, ε is the error, with the addition model for accumulation. Model construction is mainly divided into four steps: model building, model evaluation, problem presentation, and visual feedback (20).

In this study, the optimal model selected above was used to fit the HFRS case data from 2005 to 2021, and then was used to predict the data of 2022. Then, root mean square error (RMSE) and mean absolute error (MAE) were used to evaluate the fitting and prediction performances of the model, and the smaller the values of the indicators, the better the model performance. The SARIMA model was constructed using the “tseries” and “forecast” packages in R software (version 4.2.1, R Foundation for Statistical Computing, Vienna, Austria), and the “prophet” package was used to construct the prophet model.

3 Results

3.1 General characteristics

A total of 843 cases of HFRS were reported in Quzhou from 2005 to 2022, with the number of cases reported annually fluctuating in the range of 24–84. The annual average incidence was 2.15/100,000, with the lowest in 2022 (1.05/100,000) and the highest in 2007 (3.93/100,000). During the study period, the incidence decreased in a fluctuating manner, with a statistically significant difference across different years (χ2 = 112.791, p < 0.001), as shown in Figure 1. Seasonal bimodal distribution of cases, with the first peak in October–January (48.64%) and the second peak in May–July (24.79%). Among the HFRS cases, 584 were males (69.28%) and 259 were females (30.72%), with a sex ratio of 2.25:1. The incidence of HFRS in males (2.87/100,000) was higher than that in females (1.32/100,000) (χ2 = 114.188, p < 0.001). The age of onset showed a single peak distribution, ranging from 4 to 91 years old, with a median age of 51 years old. The number of cases in the age group of 50 years was the highest, accounting for 29.42% (248/843). Significant difference was observed in the incidence across different age groups, with the highest incidence in the age group of 50 years, followed by the age group of 60 years, and the overall trend showed a picture of rising first and then decreasing. The main occupation included farmers 79.95% (674/843), workers 4.86% (41/843), and students 2.61% (22/843) (Table 1).

FIGURE 1
www.frontiersin.org

Figure 1. Cases and incidence rate of HFRS in Quzhou City from 2005 to 2022.

TABLE 1
www.frontiersin.org

Table 1. Population distribution of HFRS in Quzhou City, 2005–2022.

3.2 Spatio-temporal characteristics

The spatial distribution of the average annual incidence of HFRS in Quzhou from 2005 to 2022 was shown in Figure 2. The disease occurred in 6 districts and counties of the city, and the highest annual incidence was reported in Kaihua County. Among 43 streets, 39 townships, and 18 subdistricts of Quzhou, the incidence rate was highest in Yinkeng Township (32.41/100,000), followed by Cuntou Township (20.65/100,000), Zhongcun Township (17.36/100,000), and Tongcun Township (16.92/100,000). All the above townships belonged to Kaihua County. A map of incidence by year is shown in Supplementary Figure S1.

FIGURE 2
www.frontiersin.org

Figure 2. Distribution of average annual incidence rate of HFRS in each township (street) in Quzhou from 2005 to 2022.

According to the national land resources classification system and the land use type standard, there are seven types of land use in Quzhou, namely artificial surface, woodland, water body, wetland, arable land, grassland and bare land. Combined with the cumulative incidence data from 2005 to 2022, it was found that the distribution of cases was consistent with that of cultivated land and artificial surface, as shown in Figure 3.

FIGURE 3
www.frontiersin.org

Figure 3. Incidence of HFRS in different land use types in Quzhou City, 2005–2022.

3.3 Model prediction

Using the incidence data of HFRS in Quzhou from 2005 to 2021, a time series was established in monthly basis, the time series diagram (Supplementary Figure S2) was drawn, and the component decomposition of the time series (Supplementary Figure S3) showed that the incidence of HFRS in Quzhou tended to be stable. The time series component decomposition diagram showed that the series fluctuated with seasonality, so a seasonal difference was performed on the original series. ADF test (ADF = −4.571, p < 0.05) showed that the sequence was stationary, that is, d = 0, D = 1. A seasonal product model with a period of 12 months was developed: SARIMA (p, 0, q) (P, 1, Q)12. The autocorrelation and partial autocorrelation analysis (Supplementary Figure S4) were carried out, q was determined to be 1, p was 1, and the values of P and Q were generally between 0 and 2. The model with the smallest AIC was selected as the optimal model by exhaustive method. It was concluded that the optimal model was SARIMA (1, 0, 1) (2, 1, 0)12, and the results showed that the model parameters were statistically significant (Table 2). The diagnostic diagram of the model (Supplementary Figure S5) showed the residual error tested by Ljung-Box test (χ2 = 0.241, p = 0.624). It was considered that the model was a white noise model with sufficient information extraction and could be used for prediction.

TABLE 2
www.frontiersin.org

Table 2. Parameter estimation and model verification of SARIMA model.

SARIMA (1,0,1) (2,1,0)12 model was used to fit the trend of HFRS incidence in Quzhou from 2005 to 2021. The overall trend of the actual value was consistent with the fitting value, but the peak value was not fitted well. The predicted number of HFRS cases from January to December in 2022 was compared with the actual number of monthly cases, and all actual values were within the 95% CI of the predicted values (Figure 4).

FIGURE 4
www.frontiersin.org

Figure 4. Fitting and forecast result of SARIMA.

The prophet model was constructed using the monthly incidence data of HFRS in Quzhou from 2005 to 2021 using the prophet package of R software. The interval width was set to 0.95, and other parameters were the default values. The model includes a “trend, year” time-series component (Supplementary Figure S6). Figure 5 shows that the model captures the incidence trend of HFRS well and achieves good results in the prediction of HFRS data.

FIGURE 5
www.frontiersin.org

Figure 5. Fitting and forecast result of Prophet model.

The predicted values and 95% CI of SARIMA and Prophet models are shown in Table 3. The results showed that the RMSE and MAE values of the Prophet model were smaller than those of the SARIMA (1,0,1) (2,1,0)12 model in both the training set and the test set, as shown in Table 4. The Prophet model had the best performance for this time series prediction.

TABLE 3
www.frontiersin.org

Table 3. The monthly number of HFRS cases in the Quzhou City in 2022 predicted by the SARIMA, and Prophet models.

TABLE 4
www.frontiersin.org

Table 4. Comparison of the fitting and prediction performance of the two models.

4 Discussion

Results from the present study showed that the incidence of HFRS in Quzhou fluctuated from 2005 to 2022, the incidence was higher in males than in females, and the occupational distribution was dominated by farmers. The main endemic area is located in the northwest of Quzhou, and cases are concentrated on cultivated land and artificial surfaces. The SARIMA and Prophet models were used for fitting and prediction, respectively, and it was found that the Prophet model performed better.

HFRS is classified as a Class B infectious disease in China according to the Law of the People’s Republic of China on the Prevention and Treatment of Infectious Diseases. In recent decades, it has been considered as one of the top 10 infectious diseases (27). Quzhou has reported HFRS cases every year. From 2005 to 2022, the average annual incidence was 2.15/100,000, which was higher than the average level of Zhejiang Province (0.91/100,000) (13). Differences in incidence due to different topography and climate (28). The city’s terrain is dominated by mountainous hills, with abundant rainfall and plenty of sunshine. It is suitable for the survival of rodents carrying HV. It has been shown that increased rodent populations can accelerate HFRS epidemics (29), and humid environments favor the breeding of mites on rodents, increasing the chances of the host carrying HV (30). From 2004 to 2015, the epidemic area of HFRS in China changed from northeast to southeast China (31). The incidence of HFRS in Quzhou was high from 2005 to 2007 and kept increasing. HV vaccination was carried out in areas such as Liaoning and Zhejiang in 2008 (32), and the incidence of the disease decreased significantly after immunization, so it is assumed that the use of the vaccine can reduce the incidence of HFRS (33). However, the epidemic began to rebound in 2010, which may be caused by changes in natural and social factors (34). In 2020–2022, the incidence of HFRS in the city was significantly lower and below average, and in 2022 it fell to the lowest value (1.05/100,000) in that study period. There are two possible reasons for the analysis: On the one hand, due to the outbreak of novel coronavirus infection at the end of 2019, the state has taken a series of non-pharmacological interventions such as quarantine, work suspension and hygiene improvement to improve people’s awareness of hygiene, develop the hygiene habit of frequent hand washing and wearing masks, and reduce garbage piling, which effectively reduced the phenomenon of mosquitoes and insects, and the rampant presence of insects and rodents (35, 36). On the other hand, the spread of COVID-19 limited the detection of other infectious diseases and reduced the behavior of seeking medical help. As a result, the incidence of HFRS was underestimated (37), which was the same as the change level of HFRS in China (38).

HFRS incidence in Quzhou City has obvious seasonality, showing a bimodal distribution with a large peak in October–January and a small peak in May–July, which is consistent with the seasonal peaks in Shandong Province and Jiangxi Province (11, 17), and a single-peak distribution in February–May in Hebei Province (39). It may be due to the different dominant species of rats in different regions, with the Apodemus agrarius dominated by the HTVN, which peaks in winter, while the Rattus norvegicus carries the SEOV with a high prevalence in the spring, which is the epidemiological season (32). Zhejiang Province is a mixed area of Apodemus agrarius and Rattus norvegicus, and the former accounts for a relatively high proportion (40), creating a dual-peak seasonal pattern that is slightly higher in winter than in spring.

The cases in all age groups were more likely to be male, and the sex difference may be caused by the difference in daily activities (41). Men have more opportunities for agricultural and outdoor labor, which is consistent with the occupational distribution characteristics of farmers (42). The over 50 age group accounted for nearly half of the total prevalence, which is related to its long history of agricultural production and lifestyle, its greater exposure to HFRS, and its greater susceptibility to the disease (43). Therefore, targeted health education on HFRS among the middle-aged and older adult population to enhance their self-protection awareness is an important measure for prevention and control.

According to the distribution map of HFRS epidemic and cumulative incidence rate in Quzhou City from 2005 to 2022, the northwestern part of Quzhou City is a high incidence area, in which Yinkeng Township, Cuntou Township, Zhongcun Township and other townships (streets) bordering Kaihua and Changshan Counties are high incidence areas. All of the above areas are mountainous and most of the population are farmers (44). Furthermore, the incidence of HFRS gradually spread from northwest to the east, which might be attributed to the accelerated urbanization process, urban expansion in the eastern part of Quzhou, and the susceptible population flocking in this area. Due to the change of land use type (45), the opportunity to contact with rodents increased, and the imperfect housing and surrounding infrastructure in the early stage increased the risk of infection (46). It is suggested that the hygiene of living environment should be improved to prevent and control rodents.

The land use distribution map showed that the risk of HFRS was concentrated in arable land and artificial surface near arable land, which was the same as the results in Hunan and Shaanxi Provinces (47, 48). Artificial land surface is the surface formed by artificial construction activities, including all kinds of residential land, industrial and mining, etc. Arable land mainly includes paddy fields, dry land, and vegetable fields. It was found that the type of land use affects the spread of HFRS by altering the environment in which rodents live, that crops are a major food source for rodents, and that farmers are more likely to come into contact with rodents during their farming practices (49). Living in areas near cropland with high rodent density increases the chance of HFRS infection (50).

In this study, based on the incidence data of HFRS in Quzhou from 2005 to 2022, the SARIMA model and Prophet model were constructed to fit and predict the incidence trend of HFRS. We find that Prophet’s optimal model is more effective as the evaluation indexes of RMSE and MAE are lower than those of SARIMA model, and the error of prediction results is smaller. SARIMA model is one of the most commonly used models in time series analysis. It has the advantages of simple structure and easy interpretation of parameters, and is widely used in the field of public health (51). However, it requires the data to be stationary series, and cannot well fit the nonlinear trend (52). HFRS morbidity information has both linear and non-linear characteristics, and SARIMA does not capture the information adequately. Taylor and Letham (26) in their analysis of hand, foot and mouth disease and Mohan et al. (53) in their study of COVID-19 found that the Prophet model had a higher predictive accuracy when using the ARIMA model versus the Prophet model to predict disease incidence. Since Prophet model has more advantages in dealing with outliers and trend migration, it can better predict the development and change of things over time (26). In practice, it is necessary to continuously adjust the model according to the new measured data to further improve the prediction ability of the model.

Our study also has limitations. First of all, the data of HFRS cases were collected by passive surveillance, which may have biases such as underreporting. Secondly, only the inherent characteristics of the time series were included in the construction of the model, and other exogenous variables, such as meteorological factors, economic factors, and rodent density, were not considered. Finally, because the study data came from Quzhou City, China, our conclusions apply only to the population in that area, and extrapolation of the results is limited. In the future, we need to explore the interactions of Prophet’s model in depth, taking into account various influencing factors, in order to obtain forecasts that are closer to reality.

5 Conclusion

The overall incidence of HFRS in Quzhou showed a downward trend from 2005 to 2022, but the epidemic was still serious in Cuntou town and Zhongcun Township of Kaihua County. Male middle-aged and older adult farmers should be the focus of HFRS prevention and control. In the future, we should continue to strengthen the monitoring of the epidemic dynamics and host animals of HFRS, carry out scientific rodent control and control before the onset of the high incidence season of HFRS, and promote HFRS vaccination and health education, especially for farmers. The Prophet model could well fit and predict the incidence trend of HFRS in Quzhou, and provide early warning and prediction for the occurrence of human HFRS and reliable data for the prevention and control of HFRS in Quzhou.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to BZ, bd_zhan@126.com.

Author contributions

QG: Methodology, Visualization, Writing – original draft. SW: Data curation, Software, Writing – review & editing. QW: Writing – review & editing. GC: Writing – review & editing. CF: Writing – review & editing. BZ: Funding acquisition, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research was supported by the Zhejiang Provincial Natural Science Foundation of China (Grant No. LTGG23H260001).

Acknowledgments

The authors would like to sincerely thank Quzhou Center for Disease Control and Prevention for the support of this study, and most importantly, the participants of the study and the members of the survey team.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2023.1333178/full#supplementary-material

Footnotes

References

1. Zou, LX, Chen, MJ, and Sun, L. Haemorrhagic fever with renal syndrome: literature review and distribution analysis in China. Int J Infect Dis. (2016) 43:95–100. doi: 10.1016/j.ijid.2016.01.003

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Gavrilovskaya, IN, Gorbunova, EE, and Mackow, ER. Pathogenic hantaviruses direct the adherence of quiescent platelets to infected endothelial cells. J Virol. (2010) 84:4832–9. doi: 10.1128/JVI.02405-09

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Avšič-Županc, T, Saksida, A, and Korva, M. Hantavirus infections. Clin Microbiol Infect. (2019) 21s:e6–e16. doi: 10.1111/1469-0691.12291

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Sehgal, A, Mehta, S, Sahay, K, Martynova, E, Rizvanov, A, Baranwal, M, et al. Hemorrhagic fever with renal syndrome in Asia: history, pathogenesis, diagnosis, treatment, and prevention. Viruses. (2023) 15:561. doi: 10.3390/v15020561

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Vaheri, A, Henttonen, H, Voutilainen, L, Mustonen, J, Sironen, T, and Vapalahti, O. Hantavirus infections in Europe and their impact on public health. Rev Med Virol. (2013) 23:35–49. doi: 10.1002/rmv.1722

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Tkachenko, EA, Ishmukhametov, AA, Dzagurova, TK, Bernshtein, AD, Morozov, VG, Siniugina, AA, et al. Hemorrhagic fever with renal syndrome, Russia. Emerg Infect Dis. (2019) 25:2325–8. doi: 10.3201/eid2512.181649

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Wang, Y, Wei, X, Jia, R, Peng, X, Zhang, X, Yang, M, et al. The spatiotemporal pattern and its determinants of hemorrhagic fever with renal syndrome in northeastern China: spatiotemporal analysis. JMIR Public Health Surveill. (2023) 9:e42673. doi: 10.2196/42673

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Armién, B, Muñoz, C, Cedeño, H, Salazar, JR, Salinas, TP, González, P, et al. Hantavirus in Panama: twenty years of epidemiological surveillance experience. Viruses. (2023) 15:1395. doi: 10.3390/v15061395

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Jonsson, CB, Figueiredo, LT, and Vapalahti, O. A global perspective on hantavirus ecology, epidemiology, and disease. Clin Microbiol Rev. (2010) 23:412–41. doi: 10.1128/CMR.00062-09

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Zhang, WY, Wang, LY, Liu, YX, Yin, WW, Hu, WB, Magalhaes, RJ, et al. Spatiotemporal transmission dynamics of hemorrhagic fever with renal syndrome in China, 2005-2012. PLoS Negl Trop Dis. (2014) 8:e3344. doi: 10.1371/journal.pntd.0003344

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Liu, SW, Li, JX, Zou, L, Liu, XQ, Xu, G, Xiong, Y, et al. Orthohantavirus infections in humans and rodents in the Yichun region, China, from 2016 to 2021. PLoS Negl Trop Dis. (2023) 17:e0011540. doi: 10.1371/journal.pntd.0011540

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Khan, A, Khan, M, Ullah, S, and Wei, DQ. Hantavirus: the next pandemic we are waiting for? Interdiscip Sci. (2021) 13:147–52. doi: 10.1007/s12539-020-00413-4

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Zhang, R, Zhang, N, Ling, F, Liu, Y, Guo, S, Shi, XG, et al. Study on epidemic trend of hemorrhagic fever with renal syndrome in Zhejiang province, 2005-2020. Zhonghua Liu Xing Bing Xue Za Zhi. (2021) 42:2030–6. doi: 10.3760/cma.j.cn112338-20210528-00435

PubMed Abstract | CrossRef Full Text | Google Scholar

14. She, K, Li, C, Qi, C, Liu, T, Jia, Y, Zhu, Y, et al. Epidemiological characteristics and regional risk prediction of hemorrhagic fever with renal syndrome in Shandong Province, China. Int J Environ Res Public Health. (2021) 18:8495. doi: 10.3390/ijerph18168495

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Riccò, M, Peruzzi, S, Ranzieri, S, and Magnavita, N. Occupational hantavirus infections in agricultural and forestry workers: a systematic review and Metanalysis. Viruses. (2021) 13:2150. doi: 10.3390/v13112150

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Brocato, RL, and Hooper, JW. Progress on the prevention and treatment of hantavirus disease. Viruses. (2019) 11:610. doi: 10.3390/v11070610

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Wang, T, Liu, J, Zhou, Y, Cui, F, Huang, Z, Wang, L, et al. Prevalence of hemorrhagic fever with renal syndrome in Yiyuan County, China, 2005-2014. BMC Infect Dis. (2016) 16:69. doi: 10.1186/s12879-016-1404-7

CrossRef Full Text | Google Scholar

18. Ceylan, Z. Estimation of COVID-19 prevalence in Italy, Spain, and France. Sci Total Environ. (2020) 729:138817. doi: 10.1016/j.scitotenv.2020.138817

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Luo, Z, Jia, X, Bao, J, Song, Z, Zhu, H, Liu, M, et al. A combined model of SARIMA and prophet models in forecasting AIDS incidence in Henan Province, China. Int J Environ Res Public Health. (2022) 19:5910. doi: 10.3390/ijerph19105910

CrossRef Full Text | Google Scholar

20. Xie, C, Wen, H, Yang, W, Cai, J, Zhang, P, Wu, R, et al. Trend analysis and forecast of daily reported incidence of hand, foot and mouth disease in Hubei, China by prophet model. Sci Rep. (2021) 11:1445. doi: 10.1038/s41598-021-81100-2

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Lou, HR, Wang, X, Gao, Y, and Zeng, Q. Comparison of ARIMA model, DNN model and LSTM model in predicting disease burden of occupational pneumoconiosis in Tianjin, China. BMC Public Health. (2022) 22:2167. doi: 10.1186/s12889-022-14642-3

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Alzahrani, SI, Aljamaan, IA, and Al-Fakih, EA. Forecasting the spread of the COVID-19 pandemic in Saudi Arabia using ARIMA prediction model under current public health interventions. J Infect Public Health. (2020) 13:914–9. doi: 10.1016/j.jiph.2020.06.001

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Ge, L, Zhao, Y, Zhou, K, Mu, X, Yu, H, Wang, Y, et al. Spatio-temporal pattern and influencing factors of hemorrhagic fever with renal syndrome (HFRS) in Hubei Province (China) between 2005 and 2014. PLoS One. (2016) 11:e0167836. doi: 10.1371/journal.pone.0167836

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Hernandez-Matamoros, A, Fujita, H, Hayashi, T, and Perez-Meana, H. Forecasting of COVID19 per regions using ARIMA models and polynomial functions. Appl Soft Comput. (2020) 96:106610. doi: 10.1016/j.asoc.2020.106610

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Azeez, A, Obaromi, D, Odeyemi, A, Ndege, J, and Muntabayi, R. Seasonality and trend forecasting of tuberculosis prevalence data in eastern cape, South Africa, using a hybrid model. Int J Environ Res Public Health. (2016) 13:757. doi: 10.3390/ijerph13080757

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Taylor, SJ, and Letham, B. Forecasting at scale. Am Stat. (2018) 72:37–45. doi: 10.1080/00031305.2017.1380080

CrossRef Full Text | Google Scholar

27. Bi, P, Wu, X, Zhang, F, Parton, KA, and 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) 148:276–81. doi: 10.1093/oxfordjournals.aje.a009636

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Tian, H, Yu, P, Cazelles, B, Xu, L, Tan, H, Yang, J, et al. Interannual cycles of Hantaan virus outbreaks at the human-animal interface in Central China are controlled by temperature and rainfall. Proc Natl Acad Sci U S A. (2017) 114:8041–6. doi: 10.1073/pnas.1701777114

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Castel, G, Alburkat, H, Tatard, C, Dutra, L, Criado, M, Bouilloud, M, et al. Puumala orthohantavirus circulation in its wild reservoir, the bank vole, during the 2021 outbreak of hemorrhagic fever with renal syndrome in Jura, France. Infect Dis Now. (2023) 53:104767. doi: 10.1016/j.idnow.2023.104767

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Luo, Y, Lv, H, Yan, H, Zhu, C, Ai, L, Li, W, et al. Meteorological change and hemorrhagic fever with renal syndrome epidemic in China, 2004–2018. Sci Rep. (2022) 12:20037. doi: 10.1038/s41598-022-23945-9

CrossRef Full Text | Google Scholar

31. Sun, L, and Zou, L-X. Spatiotemporal analysis and forecasting model of hemorrhagic fever with renal syndrome in mainland China. Epidemiol Infect. (2018) 146:1680–8. doi: 10.1017/S0950268818002030

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Zhang, S, Fau-Wang, S, Wang, W, Fau-Yin, S, Yin, M, Fau-Liang, W, et al. Epidemic characteristics of hemorrhagic fever with renal syndrome in China, 2006–2012. BMC Infect Dis. (2014) 14:384. doi: 10.1186/1471-2334-14-384

CrossRef Full Text | Google Scholar

33. Jung, J, Ko, SJ, Oh, HS, Moon, SM, Song, JW, and Huh, K. Protective effectiveness of inactivated hantavirus vaccine against hemorrhagic fever with renal syndrome. J Infect Dis. (2018) 217:1417–20. doi: 10.1093/infdis/jiy037

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Li, S, Zhu, L, Zhang, L, Zhang, G, Ren, H, and Lu, L. Urbanization-related environmental factors and hemorrhagic fever with renal syndrome: a review based on studies taken in China. Int J Environ Res Public Health. (2023) 20:3328. doi: 10.3390/ijerph20043328

CrossRef Full Text | Google Scholar

35. Maciel, EA, de Carvalho Al Fau-Nascimento, SF, Matos, RBNS F-d, de Matos Rb Fau-Gouveia, EL, El Fau-Reis, MGG, Fau-Ko, AIRM, et al. Household transmission of leptospira infection in urban slum communities. PLoS Negl Trop Dis. (2008) 2:e154. doi: 10.1371/journal.pntd.0000154

CrossRef Full Text | Google Scholar

36. Estee, AAL, Cramer, Y, and Lover, AA. Forgotten diseases: the neglected tropical diseases and their impact on Global Health and development. Emerg Infect Dis. (2022) 28:2140. doi: 10.3201/eid2810.220740

CrossRef Full Text | Google Scholar

37. Hibiya, K, Iwata, H, Kinjo, T, Shinzato, A, Tateyama, M, Ueda, S, et al. Incidence of common infectious diseases in Japan during the COVID-19 pandemic. PLoS One. (2022) 17:e0261332. doi: 10.1371/journal.pone.0261332

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Bai, BK, Jiang, QY, and Hou, J. The COVID-19 epidemic and other notifiable infectious diseases in China. Microbes Infect. (2022) 24:104881. doi: 10.1016/j.micinf.2021.104881

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Wang, Q, Yue, M, Yao, P, Zhu, C, Ai, L, Hu, D, et al. Epidemic trend and molecular evolution of HV family in the Main hantavirus epidemic areas from 2004 to 2016, in P.R. China. Front Cell Infect Microbiol. (2020) 10:584814. doi: 10.3389/fcimb.2020.584814

CrossRef Full Text | Google Scholar

40. Wu, Y, Wang, J, Liu, Q, Li, T, Luo, M, and Gong, Z. Practice of integrated vector surveillance of arthropod vectors, pathogens and reservoir hosts to monitor the occurrence of tropical vector-borne diseases in 2020 in Zhejiang Province, China. Front Vet Sci. (2022) 9:1003550. doi: 10.3389/fvets.2022.1003550

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Puca, E, Pipero, P, Harxhi, A, Abazaj, E, Gega, A, Puca, E, et al. The role of gender in the prevalence of human leptospirosis in Albania. J Infect Dev Ctries. (2018) 12:150–5. doi: 10.3855/jidc.9805

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Joshi, YP, Kim, EH, and Cheong, HK. The influence of climatic factors on the development of hemorrhagic fever with renal syndrome and leptospirosis during the peak season in Korea: an ecologic study. BMC Infect Dis. (2017) 17:406. doi: 10.1186/s12879-017-2506-6

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Wu, H, Wang, X, Xue, M, Wu, C, Lu, Q, Ding, Z, et al. Spatial-temporal characteristics and the epidemiology of haemorrhagic fever with renal syndrome from 2007 to 2016 in Zhejiang Province, China. Sci Rep. (2018) 8:10244. doi: 10.1038/s41598-018-28610-8

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Zhang, R, Mao, Z, Yang, J, Liu, S, Liu, Y, Qin, S, et al. The changing epidemiology of hemorrhagic fever with renal syndrome in southeastern China during 1963-2020: a retrospective analysis of surveillance data. PLoS Negl Trop Dis. (2021) 15:e0009673. doi: 10.1371/journal.pntd.0009673

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Tian, H, Yu, P, Bjørnstad, ON, Cazelles, B, Yang, J, Tan, H, et al. Anthropogenically driven environmental changes shift the ecological dynamics of hemorrhagic fever with renal syndrome. PLoS Pathog. (2017) 13:e1006198. doi: 10.1371/journal.ppat.1006198

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Tian, H, Hu, S, Cazelles, B, Chowell, G, Gao, L, Laine, M, et al. Urbanization prolongs hantavirus epidemics in cities. Proc Natl Acad Sci U S A. (2018) 115:4707–12. doi: 10.1073/pnas.1712767115

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Chen, Z, Liu, F, Li, B, Peng, X, Fan, L, and Luo, A. Prediction of hot spot areas of hemorrhagic fever with renal syndrome in Hunan Province based on an information quantity model and logistical regression model. PLoS Negl Trop Dis. (2020) 14:e0008939. doi: 10.1371/journal.pntd.0008939

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Li, X, Liu, K, Gu, X, Yuan, XJ, and Shao, ZJ. Epidemiological characteristics and environmental risk factors of hemorrhagic fever with renal syndrome in Wei River basin, China, 2005-2015. Zhonghua Liu Xing Bing Xue Za Zhi. (2018) 39:1159–64. doi: 10.3760/cma.j.issn.0254-6450.2018.09.004

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Xiao, H, Tong, X, Huang, R, Gao, L, Hu, S, Li, Y, et al. Landscape and rodent community composition are associated with risk of hemorrhagic fever with renal syndrome in two cities in China, 2006-2013. BMC Infect Dis. (2018) 18:37. doi: 10.1186/s12879-017-2827-5

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Hansen, A, Cameron, S, Liu, Q, Sun, Y, Weinstein, P, Williams, C, et al. Transmission of haemorrhagic fever with renal syndrome in China and the role of climate factors: a review. Int J Infect Dis. (2015) 33:212–8. doi: 10.1016/j.ijid.2015.02.010

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Wang, Y, Xu, C, Yao, S, and Zhao, Y. Forecasting the epidemiological trends of COVID-19 prevalence and mortality using the advanced α-Sutte Indicator. Epidemiol Infect. (2020) 148:e236. doi: 10.1017/S095026882000237X

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Wang, YW, Shen, ZZ, and Jiang, Y. Comparison of autoregressive integrated moving average model and generalised regression neural network model for prediction of haemorrhagic fever with renal syndrome in China: a time-series study. BMJ Open. (2019) 9:e025773. doi: 10.1136/bmjopen-2018-025773

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Mohan, S, Solanki, AK, Taluja, HK, Anuradha,, and Singh, A. Predicting the impact of the third wave of COVID-19 in India using hybrid statistical machine learning models: a time series forecasting and sentiment analysis approach. Comput Biol Med. (2022) 144:105354. doi: 10.1016/j.compbiomed.2022.105354

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: hemorrhagic fever with renal syndrome, spatiotemporal analysis, land use, SARIMA, Prophet

Citation: Gao Q, Wang S, Wang Q, Cao G, Fang C and Zhan B (2024) Epidemiological characteristics and prediction model construction of hemorrhagic fever with renal syndrome in Quzhou City, China, 2005–2022. Front. Public Health. 11:1333178. doi: 10.3389/fpubh.2023.1333178

Received: 04 November 2023; Accepted: 29 December 2023;
Published: 11 January 2024.

Edited by:

Jian Wu, Suzhou Municipal Hospital, China

Reviewed by:

Edmond Puca, Service of Infection Diseases University Hospital Center, Albania
Samira Sanami, Shahid Beheshti University of Medical Sciences, Iran
Hassan Alhassan, Al Jouf University, Saudi Arabia

Copyright © 2024 Gao, Wang, Wang, Cao, Fang and Zhan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Bingdong Zhan, bd_zhan@126.com

These authors have contributed equally to this work and share first authorship

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.