- Descriptive analysis of BD and environment factors
During the four-year period from 1 January, 2014 to 31 December, 2017, there were a total of 7,102 BD cases, with an average daily BD cases of 4.84. The cases of BD among male was approximately 1.17 times as high as that among female. Children under 6 years old were at the highest proportion(38.2%) of BD. The average daily mean temperature, PM2.5 and PM10 were 7.73℃, 53.44μg/m3 and 124.47μg/m3, respectively(Table 1).
- Correlation analysis between BD and environment factors
The correlation between the incidence of bacillary dysentery and PM2.5, PM10, CO, NO2 and O3 was statistically significant (p<0.05), respectively. The incidence of bacillary dysentery was negatively correlated with PM2.5, PM10, CO and NO2, while was positively correlated with O3. The correlation coefficients between the mean temperature, PM2.5, PM10 and the cases of BD were 0.58,-0.39 and -0.33, respectively.
- Exploratory analysis of the relationship between temperature and BD
The maximum temperature, minimum temperature, and mean temperature were respectively used to fit the GAM, in which the model fitted by the mean temperature was generally associated with the lowest generalized cross-validation value(GCV, means that the model fitted by mean temperature had the best predictive ability for incidence), therefore we only reported results for associations with mean temperature.
After adjusting the confounding factor such as air pressure, relative humidity, wind speed PM2.5 and PM10, the effect of mean temperature on bacterial dysentery is shown below. As shown in Figure 2a, the three-dimensional plot shows the entire surface between mean temperature and daily cases of bacillary dysentery at all lag days. The estimated effect of mean temperature was non-linear for daily cases of bacillary dysentery, with higher relative risk at hot temperature. For example, hot temperature (19℃) and cold temperature (-16℃) were positively associated with daily cases of bacillary dysentery on the current day. Cold effects (i.e., significant increases in incidence associated with cold temperature) was not apparent after a 6-day lag, with relative risk close to 1 when the temperature below 0℃. Figure 2b shows the estimated effect of mean temperature over 14 days on daily cases of bacillary dysentery. There was J-shaped relationship between mean temperature and daily cases of bacillary dysentery, with a large comfortable temperature range where the relative risk of cases of bacillary dysentery was less than 1 when mean temperature was below 10℃. In contrast, during all lag period, when mean temperature above 10℃, the risk of bacterial dysentery increased with increasing temperature.
The single-lag effects of temperature on bacillary dysentery in each subgroup was shown in Table 2, except for people aged 41 years and over, suggesting that temperature effects on BD appeared to be acute and greatest on the day of exposure(lag of 0 days). Significant temperature effect appeared after a 3-day lag for people aged 41 years and over, whereas that occurred within 0-6 days, 0-3 days, 0-11 days and 0-6 days for male, female, people aged 0-6 years and aged 7-40 years, respectively. Association between temperature and BD for people aged 0-6 years lasted longer than other subgroup, the effect value decreased with increasing lag days, and the effect on the day of exposure was the greatest, a 1℃ increase in mean temperature was associated with a 2.0%(RR=1.02, 95%CI:1.013, 1.028) increase in cases of BD. The effect of temperature on BD for people aged 41 years and over occurred within 9-14 days and the effect values were less than 1.
The cumulative effects of mean temperature on BD stratified by age and sex were illustrated in Figure 3. This suggested that mean temperature was significantly associated with bacillary dysentery in male and female, and people aged 0-6 years, 7-40 years, but not with bacillary dysentery in people aged 41 years and over. The mean temperature-bacillary dysentery patterns for female and people aged 0-6 years were similar, with the relative risk increasing significantly when mean temperature was more than 10℃. In general, hot effect of lag 0-14 days was greater than cold effect of lag 0-14 days.
- Exploratory analysis of the relationship between PM2.5 and BD
After adjusting the confounding factors such as temperature, air pressure, relative humidity, wind speed and PM10, the effect of PM2.5 on bacterial dysentery is shown below. Figure 4 depicts the effect of PM2.5 on BD. It is illustrated by 3D image and overall effect plot of the relative risk of BD. The 3D plot showed a very strong immediate positive effect of the PM2.5 concentration above around 100μg/m3 at lag 0-2 days, indicating that a high number of BD cases would occurred within the first 0-2 days of any polluted weather crossing the PM2.5 of 100μg/m3. Higher concentration of PM2.5 also seemed to have a moderate effect on BD at around 8-14 days lag period, while lower concentration of PM2.5 seemed to have a moderate positive effect on BD at around 0-6 days lag period.. Figure 4b shows the estimated effect of PM2.5 over 14 days on daily cases of bacillary dysentery. There was an inverted S-shaped curve, indicating that PM2.5 was associated with a decrease in BD occurrence with concentration of PM2.5 increasing when it was below 50ug/m3 and above 150ug/m3 , whereas PM2.5 was associated with a increase in BD occurrence with concentration of PM2.5 increase when it ranged 50ug/m3 and 150ug/m3.
Table 3 shows the single-lag effects of PM2.5 on bacillary dysentery in each subgroup. Except for people aged 7-40 years, the effects of temperature on BD appeared to be acute and greatest on the day of exposure(lag of 0 days). However, the lag time of PM2.5 on BD for people aged 0-6 years and 41 years and over, and male was nine days, which differed from that in female (14 days). The effect of PM2.5 for people aged 7-40 years was not statistically significant during the whole lag period. In particular, the effect of PM2.5 on BD for children aged 0-6 years was greater than other subgroup, each 5ug/m3 increase in PM2.5 corresponded to an increase of 2.1% (RR=1.021, 95%CI:1.014, 1.027) in cases of BD.
As shown in Figure 5, the cumulative effects of PM2.5 on different genders and age groups were nonlinear. This suggested that PM2.5 was significantly associated with BD in male and female, and people aged 0-6 years and 41 years and over when it,s concentration below 50ug/m3, but not with BD in people aged 7-40 years. The effects of PM2.5 on BD for male, and people aged 0-6 years were similar, with the relative risk increasing significantly when PM2.5 was above 50ug/m3, while the relative risk decreasing significant when PM2.5 was below 50ug/m3. The effects of PM2.5 on BD for female, and people aged 41 years and over were increasing significantly when PM2.5 was between 50ug/m3 and 150ug/m3.
- Exploratory analysis of the relationship between PM10 and BD
After adjusting the confounding factors such as temperature, air pressure, relative humidity, wind speed and PM2.5, the effect of PM10 on bacterial dysentery is shown below. Figure 6a showed the nonlinear relationship and lag effect between PM10 and BD. When the concentration of PM10 was lower than 110μg/m3, there seemed to have a moderately positive effect on BD at around 10-14 day lag period, which indicated that the lower concentration of PM10 could significantly increase the risk of BD when delayed 10-14 days. On the contrary, when the concentration of PM10 was greater than 110μg/m3, there seemed to have a strong positive effect on BD at lag 10-14 days. There were also some shorter lag periods (0-3 days) but negative effects when the concentration of PM10 was higher than around 110μg/m3. Figure 6b shows the estimated effect of PM10 over 14 days on BD. There was U-shaped relationships between PM10 and BD, means the effect of low or high concentration of PM10 was associated with a increase in BD occurrences.
As shown in Table 4, taking 110μg/m3 as reference value, at lag 14 days, each 10μg/m3 increase in PM10 corresponded to an increase of 0.3% (RR=1.003, 95%CI: 1.001, 1.006) and 0.1% (RR=1.001, 95%CI: 0.999, 1.005) in cases of BD for male and ≥41 years old people, respectively. At the lag 0 day, the increased risk for people aged 7-40 years was 0.1% (RR=1.001, 95%CI: 0.996, 1.003). The risk of BD in women and people aged 0-6 years decreased during the whole lag period, and with the number of lag days increasing, the amount of female decline increases gradually, while the reduction in people aged 0-6 years gradually decreased.
As shown in Figure 6, the cumulative effects of PM10 on different genders and age groups were nonlinear. The effects of PM10 on BD for male and female, and people aged 7-40 years and aged 41 years and over were similar, with the relative risk increasing significantly when PM10 was below 200ug/m3.