Hostname: page-component-7bb8b95d7b-dtkg6 Total loading time: 0 Render date: 2024-10-03T18:02:31.926Z Has data issue: false hasContentIssue false

Transmission pattern of shigellosis in Wuhan City, China: a modelling study

Published online by Cambridge University Press:  02 November 2021

Zeyu Zhao
Affiliation:
State Key Laboratory of Molecular Vaccinology and Molecular Diagnostics, School of Public Health, Xiamen University, Fujian, People's Republic of China
Qi Chen
Affiliation:
Hubei Provincial Center for Disease Control and Prevention, Wuhan City, Hubei Province, People's Republic of China
Bin Zhao
Affiliation:
Laboratory Department, Xiang'an Hospital of Xiamen University, State Key Laboratory of Molecular Vaccinology and Molecular Diagnostics, Xiamen, Fujian, People's Republic of China
Qingqing Hu
Affiliation:
Division of Public Health, School of Medicine, University of Utah, 201 Presidents Circle, Salt Lake City 84112, Utah, USA
Jia Rui
Affiliation:
State Key Laboratory of Molecular Vaccinology and Molecular Diagnostics, School of Public Health, Xiamen University, Fujian, People's Republic of China
Yao Wang
Affiliation:
State Key Laboratory of Molecular Vaccinology and Molecular Diagnostics, School of Public Health, Xiamen University, Fujian, People's Republic of China
Yuanzhao Zhu
Affiliation:
State Key Laboratory of Molecular Vaccinology and Molecular Diagnostics, School of Public Health, Xiamen University, Fujian, People's Republic of China
Xingchun Liu
Affiliation:
State Key Laboratory of Molecular Vaccinology and Molecular Diagnostics, School of Public Health, Xiamen University, Fujian, People's Republic of China
Jingwen Xu
Affiliation:
State Key Laboratory of Molecular Vaccinology and Molecular Diagnostics, School of Public Health, Xiamen University, Fujian, People's Republic of China
Meng Yang
Affiliation:
State Key Laboratory of Molecular Vaccinology and Molecular Diagnostics, School of Public Health, Xiamen University, Fujian, People's Republic of China
Meijie Chu
Affiliation:
State Key Laboratory of Molecular Vaccinology and Molecular Diagnostics, School of Public Health, Xiamen University, Fujian, People's Republic of China
Yanhua Su*
Affiliation:
State Key Laboratory of Molecular Vaccinology and Molecular Diagnostics, School of Public Health, Xiamen University, Fujian, People's Republic of China
Benhua Zhao*
Affiliation:
State Key Laboratory of Molecular Vaccinology and Molecular Diagnostics, School of Public Health, Xiamen University, Fujian, People's Republic of China
Tianmu Chen*
Affiliation:
State Key Laboratory of Molecular Vaccinology and Molecular Diagnostics, School of Public Health, Xiamen University, Fujian, People's Republic of China
*
Author for correspondence: Tianmu Chen, Email: 13698665@qq.com; Yanhua Su, Email: suyanhua813@xmu.edu.cn; Benhua Zhao, Email: benhuazhao@xmu.edu.cn
Author for correspondence: Tianmu Chen, Email: 13698665@qq.com; Yanhua Su, Email: suyanhua813@xmu.edu.cn; Benhua Zhao, Email: benhuazhao@xmu.edu.cn
Author for correspondence: Tianmu Chen, Email: 13698665@qq.com; Yanhua Su, Email: suyanhua813@xmu.edu.cn; Benhua Zhao, Email: benhuazhao@xmu.edu.cn
Rights & Permissions [Opens in a new window]

Abstract

The article aims to estimate and forecast the transmissibility of shigellosis and explore the association of meteorological factors with shigellosis. The mathematical model named Susceptible–Exposed–Symptomatic/Asymptomatic–Recovered–Water/Food (SEIARW) was used to explore the feature of shigellosis transmission based on the data of Wuhan City, China, from 2005 to 2017. The study applied effective reproduction number (Reff) to estimate the transmissibility. Daily meteorological data from 2008 to 2017 were used to determine Spearman's correlation with reported new cases and Reff. The SEIARW model fit the data well (χ2 = 0.00046, p > 0.999). The simulation results showed that the reservoir-to-person transmission of the shigellosis route has been interrupted. The Reff would be reduced to a transmission threshold of 1.00 (95% confidence interval (CI) 0.82–1.19) in 2035. Reducing the infectious period to 11.25 days would also decrease the value of Reff to 0.99. There was a significant correlation between new cases of shigellosis and atmospheric pressure, temperature, wind speed and sun hours per day. The correlation coefficients, although statistically significant, were very low (<0.3). In Wuhan, China, the main transmission pattern of shigellosis is person-to-person. Meteorological factors, especially daily atmospheric pressure and temperature, may influence the epidemic of shigellosis.

Type
Research article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
Copyright © The Author(s), 2021. Published by Cambridge University Press

Introduction

In 1892, dysentery was described as ‘one of the four great epidemic diseases of the world’ by the eminent physician William Osler [Reference Kotloff1]. Shigellosis (or bacterial dysentery) caused by shigella is an acute intestinal infectious disease that often occurs in summer and autumn. High-risk groups in high-income countries include travellers and MSM (men who have sex with men), and in low- and middle-income countries include children aged 1–4 years [Reference Kotloff1]. Shigella spp. are among the causes of diarrhoea in children [Reference Ma3]. Despite reductions in morbidity and mortality of shigellosis in the past 30 years, shigellosis led to approximately 164 000 annual deaths worldwide [Reference Zhang4, Reference Carlton5].

With China's rapid economic growth, the disease once associated with poverty is becoming a chronic illness [Reference Wang2]. Economic level directly influences the incidence of shigellosis [Reference Ma3]. In China, shigellosis incidence dropped by 46.29% from 2005 to 2010, but the disease burden was not evenly distributed [Reference Zhang4]. Although the water and sanitation infrastructure has improved in China, the usage of safe water and proper sanitation varies widely across the country [Reference Carlton5]. From 2006 to 2011, Wuhan City reported 36 487 bacillary dysentery cases, with a mean annual incidence of 77.4 per 100 000 population (ranging from 68.9 to 81.5 per 100 000) [Reference Li6]. Furthermore, the reported incidence varies in Anhui Province, Sichuan Province and Baise City [Reference Ma3, Reference Hao7, Reference Liu8]. A relatively high disease burden was found in Wuhan City compared to Chinese average level [Reference Wang9, Reference Kosek, Bern and Guerrant10].

Kermack and McKendrick [Reference Chen14] built a Susceptible–Infectious–Recovered (SIR) and a Susceptible–Infectious–Susceptible (SIS) model, on which the Susceptible–Infectious–Recovered–Water (SIRW) model is based. The SIRW model is used to examine disease outbreaks caused by waterborne pathogens [Reference Tien and Earn11]. The mathematical model named Susceptible–Exposed–Symptomatic/Asymptomatic–Recovered–Water/Food (SEIARW) is suitable to explore the mechanism of water/food-born disease, such as the transmission process of a small-scale outbreak in a school in Changsha City [Reference Chen12]. Thereafter, a study reported the application of the SEIARW model to the Hubei Province and confirmed that it is mainly concentrated on person-to-person transmission [Reference Chen13]. However, its prediction of transmissibility was at the entire province level. The study also showed that the incidence in the province was heterogeneous. Among its 13 cities or prefectures, Wuhan City had the highest incidence [Reference Heffernan, Smith and Wahl17]. Therefore, the transmission characteristics in Wuhan City might be different from the rest of the province. It is essential to employ the SEIARW model to calculate the transmissibility in a large city, such as Wuhan City. Previous studies always applied the basic reproduction number (R 0) to evaluate infectious diseases [Reference Chen14Reference Heffernan, Smith and Wahl17]. Because R 0 is difficult to quantify, here, the effective reproduction number (R eff) was adopted to evaluate the shigellosis transmission instead [Reference Wang18].

Recently, studies have suggested a correlation between shigellosis incidence and climatic factors, including atmospheric pressure, sun hours, temperature, relative humidity, wind speed and precipitation [Reference Li6, Reference Gao19Reference Xu22]; however, they have not clarified how meteorological factors affect the transmission of the disease. Climate may directly influence the survival of shigellosis in the environment, or indirectly affect the transmission of shigellosis through human behaviours [Reference Zhang, Bi and Hiller23]. However, there is still an epidemiological ‘black box’ [Reference Robert24] (a system of which the input and the output are known but of which nothing is known about intermediate steps) between the climatic factors and the incidence of the disease. It is essential to study the association between shigellosis transmissibility and meteorological variables. That might supply us with a new pathway to explore the transmission mechanism of the disease.

The SEIARW dynamics model was used to simulate shigellosis transmission to further quantify the transmissibility and estimate its correlation with meteorological factors.

Methods

Source of data

Data regarding shigellosis from January 2005 to December 2017 were obtained from Wuhan City through the China Information System for Disease Control and Prevention (CISDCP). Daily meteorological data (average daily value of atmospheric pressure, temperature, relative humidity, precipitation, wind speed and sun hours) from 2008 to 2017 were obtained from the National Population Health Data Centre (https://www.ncmi.cn) and the demographic data including total population, birth rate and death rate were obtained from the Wuhan Statistical Yearbook.

Model development

Based on our previous study [Reference Chen12], the population in the person-to-person route differed into five compartments, and another reservoir route including water and food was also simulated in the model (Fig. 1). As data for the infectious disease model covered more than a decade, we added demographic data such as birth and death rates. The definition and units of measurement of the storage places are shown in Table 1. Model assumptions are listed as follows:

  1. a) Shigellosis is not vertically transmitted, and newly born individuals are all susceptible. The natural change of population was considered in the model, with the natural birth rate denoted by br, and the death rate by dr.

  2. b) Susceptible people may be infected after contact with contaminated water/food or symptomatic/asymptomatic infected people. The infection rate coefficients are βW and β, respectively.

  3. c) The average incubation period of illness is 1/ω. We assumed the transition rate from E to I and E to A is equal. The asymptomatic proportion was denoted by p (0 ⩽ p ⩽ 1). Therefore, the transferred rate from exposed to asymptomatic and symptomatic is denoted by pωE and (1 – p)ωE, respectively.

  4. d) People in I or A will be recovered after a symptomatic infectious period (1/γ) or an asymptomatic infectious period (1/γ'), respectively.

  5. e) Symptomatic and asymptomatic individuals can shed the shigella, with a μI and μ'A shedding rate, respectively.

  6. f) Shigella spp. die in water/food after a period, and the daily rate of decrease of the pathogen is ɛW.

Fig. 1. Flow diagram of SEIARW model of shigellosis.

Table 1. Description and unit of variables in SEIARW model

The differential equations of the model are as follows:

$$\displaystyle{{dS} \over {dt}} = brN-\beta S( I + A) -\beta _WSW-drS$$
$$\displaystyle{{dE} \over {dt}} = \beta S( I + A) + \beta _WSW-\omega E-drE$$
$$\displaystyle{{dI} \over {dt}} = ( 1-p) \omega E-\gamma I-drI$$
$$\displaystyle{{dA} \over {dt}} = p\omega E-{\gamma }^{\prime}A-drA$$
$$\displaystyle{{dR} \over {dt}} = \gamma I + {\gamma }^{\prime}A-drR$$
$$\displaystyle{{dW} \over {dt}} = \mu I + {\mu }^{\prime}A-\varepsilon W$$

Twelve parameters were used to develop the SEIARW model, including β, βW, κ, p, ω, γ', γ, c, ɛ, μ, μ', br and dr, which were defined and valued in Table 2. Let N denotes total number in population, we perform substitution b = βN, bW = μβWN/ɛ, s = S/N, e = E/N, a = A/N, i = I/N, w = ɛW/μN, r = R/N and μ = cμ', and use the dimensionless differential equations of the model (Model 1):

$$\displaystyle{{ds} \over {dt}} = br-bs( i + a) -b_Wsw-drs$$
$$\displaystyle{{de} \over {dt}} = bs( i + a) + b_Wsw-\omega e-dre$$
$$\displaystyle{{di} \over {dt}} = ( 1-p) \omega e-\gamma i-dri$$
$$\displaystyle{{da} \over {dt}} = p\omega e-{\gamma }^{\prime}a-dra$$
$$\displaystyle{{dr} \over {dt}} = \gamma i + {\gamma }^{\prime}a-drr$$
$$\displaystyle{{dw} \over {dt}} = \varepsilon ( i + ca-w) $$

The person-to-person route of shigellosis transmission is described by the following equations (Model 2):

$$\displaystyle{{ds} \over {dt}} = br-bs( i + a) -drs$$
$$\displaystyle{{de} \over {dt}} = bs( i + a) -\omega e-dre$$
$$\displaystyle{{di} \over {dt}} = ( 1-p) \omega e-\gamma i-dri$$
$$\displaystyle{{da} \over {dt}} = p\omega e-{\gamma }^{\prime}a-dra$$
$$\displaystyle{{dr} \over {dt}} = \gamma i + {\gamma }^{\prime}a-drr$$

We calculated the R eff to estimate transmissibility by the value of b simulated for each time segment such as seven parts in 2014 (Fig. 2). The equation is as follows:

$$R_{{\rm eff}} = bs\left({\displaystyle{{1-p} \over \gamma } + \displaystyle{\,p \over {{\gamma }^{\prime}}}} \right)$$

Fig. 2. The division segments according to the prevalence reported per day of Wuhan City in 2014.

Table 2. Descriptive information of parameters in SEIARW model

- means not applicable.

Estimation of parameters

We set κ, ω, p, γ, γ', c and ɛ as 0.3125 (0–1), 1.0000 (0.3333–1), 0.1000 (0.0037–0.27), 0.0741 (0.0477–0.1428), 0.0286 (0–0.0357), 0.3125 (0–1) and 0.6931 (0–1), respectively, based on previous research (Table 2) [Reference Chen12, Reference Chen13, Reference Zhao25]. According to the Wuhan City Yearbook, the annual birth rate and death rate are 0.0063 and 0.0045, respectively. We set br and dr per day as 0.000017 and 0.000012, respectively. The total population from 2005 to 2017 in Wuhan City is shown in Supplementary Table S1. Parameters b and bW were generated by calibration of Models 1 (SEIARW) and 2 (SEIAR) to previously reported shigellosis data.

Simulation method and statistical analysis

To evaluate the contribution of the two transmission routes, we employed a method called ‘knock-out’ simulation whose theory comes from the technique of gene ‘knock-out’, which could make one of an organism's genes inoperative. Furthermore, we performed four scenarios with the ‘knock-out’ simulation (cutting off diverse routes): (a) b and bW were set to 0; (b) b was set to 0; (c) bW was set to 0; and (d) b and bW were taken value from the model fitting (None). The bW = 0 means either μ = 0 or βW = 0 or both are 0 so this is not just setting water/food-to-person contact rate to 0. The total number of infected cases was monitored in each scenario.

We adopted the decreasing the infectious period (DIP) method to estimate the intervention effectiveness in 2017 of Wuhan City. The baseline in 2017 was R eff = 1.16, γ = 0.0741 and infectious period (IP) 13.5 days. The DIP was assumed as follows: (1) an infected individual would go to a hospital immediately after becoming symptomatic, and (2) the hospital can diagnose and treat the patient. We calculated the value of R eff, IP and DIP with a change proportion (Δγ) of 10% and 20%, respectively. The DIP was calculated as follows:

$${\rm DIP} = \displaystyle{1 \over \gamma }-\displaystyle{1 \over {\gamma \times ( {1-\Delta \gamma } ) }} = {-}\displaystyle{1 \over \gamma } \times \displaystyle{{\Delta \gamma } \over {1-\Delta \gamma }}$$

Berkeley Madonna 8.3.18 (Robert Macey and George Oster of the University of California at Berkeley. Copyright ©1993–2001 Robert I. Macey & George F. Oster) was adopted for simulations. Calibration was performed via minimizing the root-mean-square deviation and the Runge–Kutta method of order four with tolerance set to 0.001 was adopted to solve the differential equations [Reference Chen12Reference Chen14, Reference Zhang, Bi and Hiller23, Reference Robert24]. Microsoft Office Excel 2019 (Microsoft, Redmond, WA, USA) was employed to run the data analysis, develop figures and calculate the value of χ 2 to evaluate the goodness of fit of the SEIARW model.

Data in 2014 were divided into seven parts that were called segments and the time step size of the simulation was set to be one day (Fig. 2). We summarised the 3 months with the highest number of new cases each year to estimate the seasonality of illness onset and used the peak R eff to estimate the seasonality of transmission.

The Spearman's correlation coefficient (r s), determined using SPSS 21.0 (IBM Corp., Armonk, NY, USA), was used to evaluate how meteorological variables were related to shigellosis. To consider the lag effect between risk factors and disease, in one study of meteorological variables and shigellosis, it was found to be 3 and 4 days [Reference Li20]. We lagged the meteorological variables by 1, 2, 3 and 4 days and calculated Spearman's correlation coefficient for each lag. Furthermore, we matched the daily meteorological variables with the R eff of the corresponding period, then calculated the mean of the meteorological variables and estimated the correlation.

Sensitivity analysis

We performed a sensitivity analysis for seven parameters which were divided into 1000 values based on their ranges from reference, then the standard deviation (s.d.) and mean of the simulation were obtained. The sensitivity analysis was only performed for data in 2014 to avoid the yearly repeated simulation. Usually, the transmissibility of infectious diseases is different between rising and reduction trends. The year 2014 was divided into seven segments and each was calibrated (Fig. 2).

Results

Epidemiology description of shigellosis in Wuhan City

As shown in Figure 3, Wuhan City reported a total of 51 948 shigellosis cases from 2005 to 2017 with a range of yearly incidence rate from 18.44 to 77.53 per 100 000 persons (median: 66.79 per 100 000 persons). The incidence rate and the number of reported cases had a significantly decreasing trend (test for linear trend χ 2 = 9.735, p = 0.002). From 2005 to 2011, the incidence rate and the number of reported cases were stable. However, from 2011 to 2015, the incidence rate and the number of reported cases decreased quickly, followed by a steady trend after 2015.

Fig. 3. Number of reported cases and incidence rate of shigellosis from 2005 to 2017 in Wuhan City, China.

Curve fitting and estimation of transmissibility

In Figure 4, Model 1 fitted the data well (χ 2 = 0.00046, p > 0.999). The bW was 3.4331 × 10–10 (95% CI 1.4179 × 10–11–6.7245 × 10–10) and b was 0.0892 (95% CI 0.0741–0.1044). The values of b and bW with no obvious seasonality in different segments per year are shown in Supplementary Table S1.

Fig. 4. Fitting effectiveness between SEIARW model and reported data in Wuhan City, China from 2005 to 2017.

We found the person-to-person transmission route is more significant (Fig. 5) according to the ‘knock-out’ simulation. We obtained the same results when cutting off the route from person-to-person (b = 0) as cutting off double routes (b = 0 and bW = 0). Meanwhile, a consistent number of cases was observed between interrupting water/food-to-person route (bW = 0) and control (b and bW ≠ 0).

Fig. 5. Simulated number of cases by four scenarios (b and bW = 0; b = 0; bW = 0; and b and bW ≠ 0 defined as ‘None’).

Although R eff had a relative fluctuation from 2005 to 2017, the trend was gradually decreasing. The mean values (Table 3) were calculated from 1.09 (95% CI 0.96–1.22) to 1.29 (95% CI 0.94–1.65). Figure 6 shows that an exponential model fits the trend of R eff well (χ 2 = 0.03473, p > 0.999). It was forecasted to reach the epidemic threshold of 1.00 (95% CI 0.82–1.19) in the year 2035, according to the exponential model.

Fig. 6. Simulated trends of R eff in Wuhan City, China.

Table 3. The yearly R eff was calculated by the SEIAR model from 2005 to 2017 in Wuhan City, China

We observed that the transmissibility of shigellosis fluctuated within the range of 1–2 (Fig. 7A). The parameter b had the same trend as R eff. Meanwhile, we found that the contribution of water/food-to-person had some effects in the previous 2 years (Fig. 7B), but the transmission from water/food had been interrupted after 2006. In addition, the results showed that there was obvious seasonality of illness onset, mainly concentrated from June to September (Fig. 8). The seasonality of peak transmissibility (R eff) was mainly concentrated in October.

Fig. 7. Simulated trends of R eff, b and bW per segment from 2005 to 2017 in Wuhan City, China. (A) R eff; (B) b and bW. The numbers above the Years represent segments per year. The detailed dates of several numbers above each year are shown in Supplementary Table S1.

Fig. 8. The peak of shigellosis onset (3 months of the top number of cases) and transmission (R eff).

Correlation between shigellosis and meteorological factors

Our results (Fig. 9) showed significant correlation between new cases and atmospheric pressure (r s = −0.282), temperature (r s = 0.299), relative humidity (r s = −0.178), wind speed (r s = −0.076) and sun hours (r s = 0.096); there was no significant correlation between new cases and precipitation (r s = −0.010). However, we just found a positive significant correlation between R eff and sun hours (r s = 0.211). Furthermore, there was no obvious change in r s when the meteorological variables were lagged by 1, 2, 3 and 4 days.

Fig. 9. The correlation between meteorological variables with new cases and R eff. *Significant difference from 0.

Modelling DIP

We further simulated different scenarios to decrease the infectious period (increase the percentage of parameter γ) and to determine the conditions to achieve R eff = 1.00. The results showed that the transmission of shigellosis would probably be interrupted (R eff = 0.99) if the infectious period were shortened by 2.25 days (Table 4).

Table 4. Reducing the R eff down to the transmission threshold (R eff = 1.00) by shortening the symptomatic infectious period

Sensitivity analysis

We obtained the same results in the simulation when setting the parameters c and ɛ (Fig. 10G) to the maximum, minimum and mean values, respectively. The SEIARW model was not sensitive to three parameters, including ω, ɛ and c, but it was sensitive to four parameters, including κ, p, γ and γ' (Fig. 10).

Fig. 10. The sensitivity analysis of parameters in the SEIARW model in 2014. Red curve, black curve and grey shaded area represent the parameters used in our model, mean value and range of standard deviation, respectively. (The upper limit of the shaded area represents the most extreme possibility when taking the maximum value of the parameter.) A narrow grey shaded interval means less sensitivity to this parameter.

Discussion

This study innovatively explored the association between meteorological variables and R eff of shigellosis. We investigated the transmissibility and forecasted the transmission of shigellosis by using the SEIARW model for the whole population in a large city and further explored the relationship between transmissibility and meteorological factors.

Simulation validity

This study adopted the SEIARW dynamics model to explore the feature of shigellosis transmission. Although this transmission model has a disadvantage in that it could fit the data poorly when the disease is sporadic, shigellosis was epidemic in Wuhan City during the study period. Therefore, our model is suitable to perform the simulation.

We added the birth and death rates in the assumptions to the model. Compared with a previous shigellosis model of Hubei Province, we further considered the variation of demographic features [Reference Chen13]. According to the χ 2 test, the SEIARW model fits the reported data well. Our findings suggest the validity of the simulation, and this is consistent with our previous findings [Reference Chen12, Reference Chen13, Reference Zhao25]. Our results indicate that the simulation was more sensitive to parameters κ, p, γ and γ'. We recommend that further study could collect the parameters κ, p, γ and γ' from first-hand data, rather than derived from references.

Description of shigellosis

Several studies indicated a reduction of the incidence of shigellosis in China [Reference Wang2, Reference Zhang4, Reference Wang9, Reference Yan26, Reference Zhao27], but the disease has remained a high burden in Wuhan City. A study indicated the incidence heterogeneity in Hubei Province [Reference Chen13]. Compared with Hubei Province overall, the incidence in Wuhan City was the highest (median: 66.79 per 100 000 persons). Kotloff et al. indicated that water/food transmission may still play a significant role in transmitting the disease, which is often further propagated by person-to-person transmission [Reference Kotloff1]. This is consistent with our result of ‘knock-out’ simulations when setting bW = 0, which means either μ = 0 or βW = 0 or both are 0. This might relate to the following reasons: cleaner water and toilets reducing the frequency of people coming into contact with water polluted by Shigella spp. (it can be considered that βW = 0); or managing shigellosis cases so that they rarely excrete pathogens to drinking water and food (it can be considered that μ = 0). Although transmission has been interrupted when b and bW are set to zero, there is still a section of exposed people in compartment E who will develop cases. The SEIAR model across sex was employed in a study to discuss interpersonal transmissibility; different transmissibility was found in different genders [Reference Zhao25]. Meanwhile, a study of Hubei Province confirmed the interruption of the water/food-to-person route [Reference Chen13], which is closely related to the improvement of water, lavatories and food safety in China. Therefore, the control of shigellosis should be focused on person-to-person transmission, including interventions of case isolation, treating the patients and hand washing.

Transmissibility and seasonality of shigellosis

Compared with influenza, Ebola virus disease and norovirus infection, shigellosis is not highly contagious [Reference Chen14Reference Chen16, Reference Yang28]. We found the range of R eff to be 1.09–1.29 in Wuhan City (higher than the epidemic threshold of 1.00). This finding suggests that one symptomatic/asymptomatic individual can infect at least 1.03 susceptible persons.

According to an exponential model of the transmissibility of shigellosis, R eff has a decreasing trend. We forecast that shigellosis in Wuhan City would not lose epidemic transmission capability until 2035 (R eff = 1.00). On one hand, meteorological factors can influence the transmissibility of shigellosis [Reference Zhang4]. On the other hand, our disease control strategy may be strengthened, which might also affect R eff. However, we found that the transmission of shigellosis would be interrupted in Hubei Province after 2029 [Reference Chen13]. This may be related to the fact that Wuhan City has the highest incidence of shigellosis in Hubei Province. Also, the transmission features in different areas may vary, and a relatively high disease burden may increase the duration of the transmission. The transmissibility in different areas may be influenced by environmental factors, economic factors, education and population density.

Simulated scenarios to shorten the infectious period showed that reducing the infectious period by 2.25 days achieved R eff = 0.99. Therefore, reducing the infectious period is especially important for secondary prevention. We should take measures such as isolation or antibiotic treatment as early as possible.

The obvious seasonality of shigellosis onset was observed in this study (mainly concentrated from June to September). This may be related to the summer vacation of primary and middle school students from July to August. Furthermore, we found that there is a seasonality of peak transmissibility (R eff), which was mainly concentrated in October. We consider that it may be related to the outbreak of shigellosis in October, and the seasonality of transmissibility needs to be further studied.

Correlation between shigellosis and meteorological factors

Several studies have demonstrated the correlation between reported cases of shigellosis and climatic factors over many years, such as temperature (positive), rainfall (positive), relative humidity (positive), wind speed (negative), sunshine duration (negative) and atmospheric pressure (negative) [Reference Li20, Reference Zhang, Bi and Hiller23, Reference Kelly-Hope29, Reference Meng30]. In this study, we found a significant correlation between new cases and atmospheric pressure, temperature, wind speed and sun hours. A previous study indicated a positive correlation between the monthly incidence of shigellosis and both average temperature (r s = 0.878) and average rainfall (r s = 0.931) in Beijing City [Reference Du31]. However, we did not find a significant correlation between new cases and precipitation (r s = 0.010); and we found a weak correlation between new cases and temperature (r s = 0.299). These findings suggest that analysing the correlation with daily data might yield weaker results than with monthly data [Reference Xu22]. A previous study showed a weak negative correlation of sun hours with shigellosis incidence in Chaoyang City of the Liaoning Province (load factor: −0.15) [Reference Zhao32]. In our study, a weak positive correlation was found between new cases and sun hours. This might be because the sun hours affected human behaviours such as increasing travelling.

R eff is an indicator reflecting person-to-person transmissibility. Two previous studies found that person-to-person is an important mode of transmission [Reference Chen13, Reference Zhao25], as did our study. We consider that meteorological factors mainly influence the survival of bacteria and human behaviour. Rainfall and high temperatures can enhance the reproductive capacity of bacteria. Social behaviours such as shaking hands, embracing and kissing can increase the risk of infection. Although the attribution analysis could not be directly employed between disease and meteorological variables, we calculated Spearman's correlation coefficient between meteorological variables and R eff. We found a weak positive correlation between R eff and sun hours. This is different from the results of the correlation between meteorological variables and new cases. Some studies have found no correlation or a small negative correlation between shigellosis cases and sun hours [Reference Du31, Reference Song33]. Sun hours may influence human behaviour. For example, it may reduce social behaviours such as shaking hands if sweating. This variation in results also may be due to the nonlinear relation between R eff and new cases. The incidence of shigellosis could be related to meteorological and other variables.

Limitations

Some studies have indicated an effect of meteorological variables on shigellosis incidence with a lag of 3 or 4 days or even several weeks [Reference Li20, 34]. As the incubation period is typically 1–3 days [Reference Xiao35], a lag effect of more than 4 days was not analysed in this study. In addition, seasonality should be considered in future studies and incorporated in model development. Tien and Earn collected the key data for the SIWR model, including the parameters α (person–reservoir contact rate (‘shedding’) and 1/ξ (pathogen lifetime in water reservoir)), which is very important for modelling results [Reference Chen15]. In this study, the key parameters ɛ and c of our model were from a reference, not first-hand data, which might underestimate the contribution of water/food in shigellosis transmission. We could not consider the population shifting in the compartmental model because we did not collect the internal migrants. Furthermore, further analysis such as Poisson regression should be performed to explore the association between shigellosis onset and several additional factors (economic, population shifting, age and so on). For further work that the pathogen concentration in the water is measured on a regular basis, be considered and the incidence of shigellosis could be related to this and other variables.

Conclusions

In Wuhan City, the incidence of shigellosis was relatively high from 2005 to 2017. The main transmission pattern is person-to-person. The transmission of shigellosis might be stopped in the year 2035 according to our forecast. Meteorological factors, especially daily atmospheric pressure and temperature, may influence the epidemic of shigellosis.

Supplementary material

The supplementary material for this article can be found at https://doi.org/10.1017/S0950268821002363

Acknowledgements

The authors thank all study participants for providing the data and field investigators for collecting the data. The authors also thank Undergraduate Innovation Practice Platform of School of Public Health, Xiamen University, for support.

Author contributions

TC, BZ and YS designed the study. QC collected data. TC, ZZ, QC, BinZ, JR, YW, YZ, XL, JX, MY and MC performed the analysis. TC, ZZ and QH wrote the first draft of this paper. All authors contributed to the writing of the manuscript.

Financial support

This work was partly supported by the Open Research Fund of State Key Laboratory of Molecular Vaccinology and Molecular Diagnostics (SKLVD2018KF001, SKLVD2018KF002 and SKLVD2019KF005), the Bill & Melinda Gates Foundation (INV-005834), the Science and Technology Program of Fujian Province (No: 2020Y0002), and the Xiamen New Coronavirus Prevention and Control Emergency Tackling Special Topic Program (No: 3502Z2020YJ03).

Conflict of interest

None.

Ethical standards

This effort of outbreak investigation and control was part of routine responsibility of Hubei Provincial Center for Disease Control and Prevention; therefore, institutional review and informed consent were not required for this study. All data analysed were anonymised.

Data availability statement

The datasets used and analysed during the current study are available from Qi Chen () on reasonable request.

Footnotes

*

These authors contributed equally to this study.

References

Kotloff, KL et al. (2018) Shigellosis. Lancet (London, England) 391, 801812.CrossRefGoogle ScholarPubMed
Wang, L et al. (2005) Preventing chronic diseases in China. Lancet (London, England) 366, 18211824.CrossRefGoogle ScholarPubMed
Ma, Y et al. (2015) Spatio-temporal pattern and socio-economic factors of bacillary dysentery at county level in Sichuan Province, China. Scientific Reports 5, 15264.CrossRefGoogle ScholarPubMed
Zhang, H et al. (2016) Patterns of bacillary dysentery in China, 2005–2010. International Journal of Environmental Research and Public Health 13, 164.CrossRefGoogle Scholar
Carlton, EJ et al. (2012) Regional disparities in the burden of disease attributable to unsafe water and poor sanitation in China. Bulletin of the World Health Organization 90, 578587.CrossRefGoogle ScholarPubMed
Li, Z et al. (2013) Identifying high-risk areas of bacillary dysentery and associated meteorological factors in Wuhan, China. Scientific Reports 3, 3239.CrossRefGoogle ScholarPubMed
Hao, Y et al. (2019) Effects of ambient temperature on bacillary dysentery: a multi-city analysis in Anhui Province, China. The Science of the Total Environment 671, 12061213.CrossRefGoogle ScholarPubMed
Liu, X et al. (2017) The effects of floods on the incidence of bacillary dysentery in baise (Guangxi Province, China) from 2004 to 2012. International Journal of Environmental Research and Public Health 14, 179.CrossRefGoogle ScholarPubMed
Wang, XY et al. (2006) Trend and disease burden of bacillary dysentery in China (1991–2000). Bulletin of the World Health Organization 84, 561568.CrossRefGoogle Scholar
Kosek, M, Bern, C and Guerrant, RL (2003) The global burden of diarrhoeal disease, as estimated from studies published between 1992 and 2000. Bulletin of the World Health Organization 81, 197204.Google ScholarPubMed
Tien, JH and Earn, DJ (2010) Multiple transmission pathways and disease dynamics in a waterborne pathogen model. Bulletin of Mathematical Biology 72, 15061533.CrossRefGoogle Scholar
Chen, T et al. (2014) Investigation of key interventions for shigellosis outbreak control in China. PLoS ONE 9, e95006.CrossRefGoogle ScholarPubMed
Chen, Q et al. (2020) Epidemiological characteristics and transmissibility of shigellosis in Hubei Province, China, 2005–2017. BMC Infectious Diseases 20, 272.CrossRefGoogle Scholar
Chen, T et al. (2016) Evidence-based interventions of Norovirus outbreaks in China. BMC Public Health 16, 1072.CrossRefGoogle ScholarPubMed
Chen, T et al. (2014) Risk of imported Ebola virus disease in China. Travel Medicine and Infectious Disease 12, 650658.CrossRefGoogle ScholarPubMed
Chen, TM et al. (2017) The transmissibility estimation of influenza with early stage data of small-scale outbreaks in Changsha, China, 2005–2013. Epidemiology & Infection 145, 424433.CrossRefGoogle Scholar
Heffernan, JM, Smith, RJ and Wahl, LM (2005) Perspectives on the basic reproductive ratio. Journal of the Royal Society Interface 2, 281293.CrossRefGoogle ScholarPubMed
Wang, Y et al. (2011) Hand, foot, and mouth disease in China: patterns of spread and transmissibility. Epidemiology 22, 781792.CrossRefGoogle ScholarPubMed
Gao, L et al. (2014) Meteorological variables and bacillary dysentery cases in Changsha City, China. American Journal of Tropical Medicine and Hygiene 90, 697704.CrossRefGoogle ScholarPubMed
Li, ZJ et al. (2015) Nonlinear and threshold of the association between meteorological factors and bacillary dysentery in Beijing. China. Epidemiology & Infection 143, 35103519.CrossRefGoogle ScholarPubMed
Nie, C et al. (2014) Socio-economic factors of bacillary dysentery based on spatial correlation analysis in Guangxi Province, China. PLoS ONE 9, e102020.CrossRefGoogle ScholarPubMed
Xu, C et al. (2017) Spatiotemporal risk of bacillary dysentery and sensitivity to meteorological factors in Hunan Province, China. International Journal of Environmental Research and Public Health 15, 47.CrossRefGoogle ScholarPubMed
Zhang, Y, Bi, P and Hiller, JE (2008) Weather and the transmission of bacillary dysentery in Jinan, northern China: a time-series analysis. Public Health Reports 123, 6166.CrossRefGoogle ScholarPubMed
Robert, M. black box epidemiology. Collins Dictionary of Medicine. 2004; 2005. Available at https://medical-dictionary.thefreedictionary.com/black+box+epidemiology (Accessed 21 October 2020).Google Scholar
Zhao, ZY et al. (2020) Relative transmissibility of shigellosis among male and female individuals: a modeling study in Hubei Province, China. Infectious Diseases of Poverty 9, 39.CrossRefGoogle ScholarPubMed
Yan, C et al. (2018) Spatiotemporal characteristics of bacillary dysentery from 2005 to 2017 in Zhejiang Province, China. International Journal of Environmental Research and Public Health 15, 1826.CrossRefGoogle ScholarPubMed
Zhao, L et al. (2017) An 11-year study of shigellosis and Shigella species in Taiyuan, China: active surveillance, epidemic characteristics, and molecular serotyping. Journal of Infection and Public Health 10, 794798.CrossRefGoogle ScholarPubMed
Yang, Y et al. (2009) The transmissibility and control of pandemic influenza A (H1N1) virus. Science (New York, N.Y.) 326, 729733.CrossRefGoogle ScholarPubMed
Kelly-Hope, LA et al. (2008) Temporal trends and climatic factors associated with bacterial enteric diseases in Vietnam, 1991–2001. Environmental Health Perspectives 116, 712.CrossRefGoogle ScholarPubMed
Meng, Q et al. (2019) Epidemiological characteristics of bacillary dysentery from 2009 to 2016 and its incidence prediction model based on meteorological factors. Environmental Health and Preventive Medicine 24, 82.CrossRefGoogle ScholarPubMed
Du, Z et al. (2018) Association between distribution of bacillary dysentery and meteorological factors in Beijing, 2004–2015. Zhonghua Liu Xing Bing Xue Za Zhi 39, 656660.Google Scholar
Zhao, Y et al. (2016) Association between meteorological factors and bacillary dysentery incidence in Chaoyang city, China: an ecological study. BMJ Open 6, e013376.CrossRefGoogle Scholar
Song, YJ et al. (2018) The epidemiological influence of climatic factors on shigellosis incidence rates in Korea. International Journal of Environmental Research and Public Health 15, 2209.CrossRefGoogle ScholarPubMed
World Health Organization. Foodborne disease outbreaks: guidelines for investigation and control: World Health Organisation. Available at https://apps.who.int/iris/handle/10665/43771 (Accessed 21 October 2020).Google Scholar
Xiao, GG et al. (2012) A school outbreak of Shigella sonnei infection in China: clinical features, antibiotic susceptibility and molecular epidemiology. Indian Pediatrics 49, 287290.CrossRefGoogle ScholarPubMed
Figure 0

Fig. 1. Flow diagram of SEIARW model of shigellosis.

Figure 1

Table 1. Description and unit of variables in SEIARW model

Figure 2

Fig. 2. The division segments according to the prevalence reported per day of Wuhan City in 2014.

Figure 3

Table 2. Descriptive information of parameters in SEIARW model

Figure 4

Fig. 3. Number of reported cases and incidence rate of shigellosis from 2005 to 2017 in Wuhan City, China.

Figure 5

Fig. 4. Fitting effectiveness between SEIARW model and reported data in Wuhan City, China from 2005 to 2017.

Figure 6

Fig. 5. Simulated number of cases by four scenarios (b and bW = 0; b = 0; bW = 0; and b and bW ≠ 0 defined as ‘None’).

Figure 7

Fig. 6. Simulated trends of Reff in Wuhan City, China.

Figure 8

Table 3. The yearly Reffwas calculated by the SEIAR model from 2005 to 2017 in Wuhan City, China

Figure 9

Fig. 7. Simulated trends of Reff, b and bW per segment from 2005 to 2017 in Wuhan City, China. (A) Reff; (B) b and bW. The numbers above the Years represent segments per year. The detailed dates of several numbers above each year are shown in Supplementary Table S1.

Figure 10

Fig. 8. The peak of shigellosis onset (3 months of the top number of cases) and transmission (Reff).

Figure 11

Fig. 9. The correlation between meteorological variables with new cases and Reff. *Significant difference from 0.

Figure 12

Table 4. Reducing the Reff down to the transmission threshold (Reff = 1.00) by shortening the symptomatic infectious period

Figure 13

Fig. 10. The sensitivity analysis of parameters in the SEIARW model in 2014. Red curve, black curve and grey shaded area represent the parameters used in our model, mean value and range of standard deviation, respectively. (The upper limit of the shaded area represents the most extreme possibility when taking the maximum value of the parameter.) A narrow grey shaded interval means less sensitivity to this parameter.

Supplementary material: File

Zhao et al. supplementary material

Zhao et al. supplementary material 1

Download Zhao et al. supplementary material(File)
File 20.4 KB
Supplementary material: Image

Zhao et al. supplementary material

Zhao et al. supplementary material 2

Download Zhao et al. supplementary material(Image)
Image 2.1 MB