A long-term retrospective analysis of the haemorrhagic fever with renal syndrome epidemic from 2005 to 2021 in Jiangxi Province, China

Jiangxi is one of the provinces in China most seriously affected by the haemorrhagic fever with renal syndrome (HFRS) epidemic. The aim of this paper was to systematically explore the HFRS epidemic in Jiangxi from the perspective of Hantavirus (HV) prevalence in rodents and humans and virus molecular characteristics. Individual information on all HFRS cases in Jiangxi from 2005 to 2021 was extracted from the China Information System for Disease Control and Prevention. All S and M fragment sequences of the Seoul virus and Hantan virus strains uploaded by Jiangxi and its neighbouring provinces and some representative sequences from provinces in China or some countries of Southeast Asia with the highest HV prevalence were retrieved and downloaded from NCBI GenBank. Periodogram and spatial autocorrelation were adopted for temporal periodicity and spatial clustering analysis of the HFRS epidemic. Joinpoint regression was utilized to explore the changing morbidity trend patterns of HFRS. Multiple sequence alignment and amino acid variation analysis were used to explore the homology and variation of strain prevalence in Jiangxi. Based on monthly morbidity time series, the periodogram analysis showed that the prevalence of HFRS had periodicities of 6 months and 12 months. Spatial autocorrelation analysis showed that HFRS distributed in Jiangxi was not random, with a “High-High” clustering area around Gaoan County. HFRS morbidity among the 0 ~ 15-year-old and ~ 61-year-old or older populations in Jiangxi increased significantly during the period of 2008–2015. Generally, HFRS morbidity was significantly positively correlated with the index of rat with virus (IRV) (r = 0.742) in the counties surrounding Gaoan from 2005 to 2019. HTNV strains in Jiangxi were in one independent branch, while the SEOV strains in Jiangxi were relatively more diverse. Both the YW89-15 and GAW30/2021 strains shared approximately 85% nucleotide homology and approximately 97% amino acid homology with their corresponding standard strains and vaccine strains. GAW30/2021 and YW89-15 had some amino acid site variations in nucleoprotein, glycoprotein precursor and RNA-dependent polymerase with their corresponding vaccine strains Z10 (HTNV) and Z37 (SEOV). The HFRS epidemic in Jiangxi has obvious temporal periodicity and spatial clustering, and the significant increase in the non-Immunization Expanded Program (EPI) targeted population (children and elderly) suggests that HFRS vaccination in this population needs to be considered. Although applying the EPI played a certain role in curbing the incidence of HFRS in Jiangxi from the perspective of ecological epidemiology, HTNV and SEOV strains prevalent in Jiangxi have some amino acid site variations compared to their corresponding vaccine strains, suggesting that HV variation needs to be continuously monitored in the future to observe vaccine protective efficiency.