Association and Prediction of Influenza-like Illness with Meteorological Factors in Mississippi, USA

SHI Xue Zhong SONG Zhi Juan YANG Yong Li BAO Jun Zhe ZHU Hui Li LUO Zi Xiao WANG Yu Ping WU Mei Na JIA Xiao Can

SHI Xue Zhong, SONG Zhi Juan, YANG Yong Li, BAO Jun Zhe, ZHU Hui Li, LUO Zi Xiao, WANG Yu Ping, WU Mei Na, JIA Xiao Can. Association and Prediction of Influenza-like Illness with Meteorological Factors in Mississippi, USA[J]. Biomedical and Environmental Sciences, 2022, 35(10): 962-967. doi: 10.3967/bes2022.123
Citation: SHI Xue Zhong, SONG Zhi Juan, YANG Yong Li, BAO Jun Zhe, ZHU Hui Li, LUO Zi Xiao, WANG Yu Ping, WU Mei Na, JIA Xiao Can. Association and Prediction of Influenza-like Illness with Meteorological Factors in Mississippi, USA[J]. Biomedical and Environmental Sciences, 2022, 35(10): 962-967. doi: 10.3967/bes2022.123

doi: 10.3967/bes2022.123

Association and Prediction of Influenza-like Illness with Meteorological Factors in Mississippi, USA

More Information
    Author Bio:

    SHI Xue Zhong, male, born in 1964, Professor, majoring in the method and application of health statistics

    Corresponding author: JIA Xiao Can, PhD, Tel: 15903660307, E-mail: jxc@zzu.edu.cn
图(5) / 表ll (6)
计量
  • 文章访问数:  656
  • HTML全文浏览量:  331
  • PDF下载量:  57
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-05-18
  • 录用日期:  2022-08-29
  • 刊出日期:  2022-10-20

Association and Prediction of Influenza-like Illness with Meteorological Factors in Mississippi, USA

doi: 10.3967/bes2022.123
    作者简介:

    SHI Xue Zhong, male, born in 1964, Professor, majoring in the method and application of health statistics

    通讯作者: JIA Xiao Can, PhD, Tel: 15903660307, E-mail: jxc@zzu.edu.cn

English Abstract

SHI Xue Zhong, SONG Zhi Juan, YANG Yong Li, BAO Jun Zhe, ZHU Hui Li, LUO Zi Xiao, WANG Yu Ping, WU Mei Na, JIA Xiao Can. Association and Prediction of Influenza-like Illness with Meteorological Factors in Mississippi, USA[J]. Biomedical and Environmental Sciences, 2022, 35(10): 962-967. doi: 10.3967/bes2022.123
Citation: SHI Xue Zhong, SONG Zhi Juan, YANG Yong Li, BAO Jun Zhe, ZHU Hui Li, LUO Zi Xiao, WANG Yu Ping, WU Mei Na, JIA Xiao Can. Association and Prediction of Influenza-like Illness with Meteorological Factors in Mississippi, USA[J]. Biomedical and Environmental Sciences, 2022, 35(10): 962-967. doi: 10.3967/bes2022.123
  • Seasonal influenza viruses have caused significant morbidity and mortality worldwide. Especially in the United States, influenza viruses spread during the winter, leading to a sharp increase in acute respiratory infections and hospitalizations. Influenza affects approximately 5% to 20% of the population each season[1]. Influenza is a climate-related acute respiratory infectious disease with a seasonal trend. In temperate regions, influenza infection usually has an obvious annual peak, which mainly occurs in winter months. However, the seasonal patterns in tropical and subtropical regions seem to be more complex. In the tropics, influenza activity can occur throughout the year or as multiple epidemics, sometimes varying from year-to-year[2].

    Several meteorological factors are responsible for the seasonal outbreak of influenza. These include temperature, relative humidity, pressure, wind speed, precipitation, and other factors. The incidence and mortality of influenza rise with the reduction of temperature[3]. The exposure-lag-response relationships between influenza virus activity and environmental factors, such as temperature, are not completely defined. Available evidence is often based on traditional statistical models, which do not reveal the cumulative or lag effects of environmental factors on influenza. Moreover, despite the increasing number of population-based estimates of influenza-associated excess deaths and hospitalizations, few studies of influenza-like illness (ILI) have been published. To address this shortcoming, we used distributed lagged nonlinear models (DLNMs) to explore the influence of meteorological factors on ILI.

    The study site is the state of Mississippi in the southern USA. The state is located in a humid subtropical climate region, with an annual average temperature ranging from 17 °C in the north to 20 °C along the coast. Normal precipitation ranges from approximately 1,270 to 1,650 mm, and is fairly evenly spread throughout the year. ILI epidemics of various magnitudes have plagued Mississippi annually for decades. The study utilized weekly data on ILI incidence from the United States Outpatient Influenza-like Illness Surveillance Network (ILINet). The data were from the first week of 2011 to week 29 of 2020. ILI is defined as “fever (temperature of 100 °F (37.8 °C) or greater) and a cough and/or a sore throat without a known cause other than influenza” (https://www.cdc.gov/flu/weekly/overview.htm accessed on November 24, 2021). Data of weekly meteorological variables, including mean temperature, mean wind speed, accumulated precipitation, mean atmospheric pressure, and mean relative humidity, were obtained from the National Oceanic and Atmospheric Administration. We also downloaded data of particulate matter ≤ 2.5 µm in diameter (PM2.5) from the US Environmental Protection Agency. The data were selected based on previous studies about influenza and other respiratory diseases.

    DLNMs were employed to determine the correlation between meteorological factors and ILI incidence. The Poisson regression was constructed with a quasi-Poisson function that allowed the excess dispersion in the weekly number of ILI cases to combine DLNMs. The weekly ILI incidence was calculated as:

    $$\begin{aligned} E\left({Y}_{t}\right)=\; &\alpha +cb.x+\sum ns\left(meteorological\;factors,df\right)+\\ &ns\left({{PM}}_{2.5},df\right)+ns\left(time,df\right)+factor\left(holiday\right) \end{aligned}$$

    where E(Yt) is the weekly ILI incidence; t is the week during the study period from 2011 through 2020; α is the intercept; cb represents the cross-basis matrix in the distributed lag models; x is the interest study factor; ns() is the natural cubic spline; df is the degree of freedom; time represents seasonality and long term trends with a degree of freedom of 7/year; and holiday refers to a binary variable of public holidays used to control the impact of public holidays. If one of the factors was studied, the remaining five variables were adjusted as covariates. Based on the potential lag effect and influenza incubation period reported in previous studies, the maximum time lag was selected as 4 weeks. To select the best model and evaluate the robustness of the model, df of each factor was determined based on previous studies and combined with Akaike Information Criterion for quasi-Poisson (Q-AIC). The optimal df of each factor was 3, except the temperature (df = 4). The “dlnm” package in R3.6.0 was used to create the DLNMs

    The seasonal autoregressive integrated moving average with exogenous variables (SARIMAX) model is an extension of SARIMA model. The input series includes covariates, such as meteorological factors. In the resulting SARIMAX (p, d, q) (P, D, Q)s (X), p, d, and q denote the order of autoregression, order of difference, and order of moving average, respectively; P, D, and Q denote the order of seasonal autoregression, order of difference, and order of moving average, respectively; and s is the specific value of cycle; s is the cycle of American influenza (52 weeks); and X is the vector of external variables. The data in the training set, from the first week of 2011 to week 52 of 2018, were used to construct the SARIMA model. The data in the validation set from the first week of 2019 to week 29 of 2020 were used to test the effect of model. The autocorrelation function (ACF) and partial autocorrelation function (PACF) used to estimate the model residuals were tested by the Ljung-Box Q test. The minimum Bayesian information criterion (BIC) was taken as the optimal SARIMAX model. The “TSA” package in R3.6.0 was used to create the SARIMAX model.

    From the first week of 2011 to week 29 of 2020, a total of 280, 871 ILI cases were reported in Mississippi. The weekly ILI incidence fluctuated from 3.35 to 81.23 per 100,000 population (Supplementary Table S1, available in www.besjournal.com). ILI occurred throughout the year. The peaks of incidence were in winter months, but with temporal variations in the peaks among different years. The temperature and PM2.5 were negatively correlated with ILI incidence in Mississippi. Wind speed and pressure were positively correlated with ILI incidence. Temperature very strongly correlated with ILI incidence, and pressure moderately strongly correlated. Correlation of the other factors was weakly or absent (Supplementary Table S2, Supplementary Figure S1, available in www.besjournal.com).

    Table S1.  Descriptive summary for ILI incidence, meteorological factors and PM2.5 in Mississippi, 2011–2020

    VariablesMeanStdP25P50P75Min.Max.
    ILI incidence (per 100,000) 18.8113.469.6114.2022.673.3581.23
    Temperature (°C)18.547.2612.1118.9725.730.1330.26
    Wind speed (m/s)2.730.812.052.713.300.825.03
    Precipitation (mm)31.2537.673.3119.5642.930.00283.21
    Pressure (kpa)100.550.35100.32100.50100.7799.75101.97
    Relative humidity (%)71.467.1067.1471.7176.4347.2988.57
    PM2.5 (μg/m3)9.072.827.178.4410.493.4724.20
      Note. ILI: influenza-like illness.

    Table S2.  Spearman correlation of meteorological factors and PM2.5 with ILI incidence in Mississippi, 2011–2020

    TemperatureWind speedPrecipitationPressureRelative humidityPM2.5
    ILI incidence−0.854**0.430**0.0490.539**−0.087−0.360**
      Note. **P < 0.001. ILI: influenza-like illness.

    Figure S1.  Scatter plot of meteorological factors and PM2.5 with ILI incidence. ILI: influenza-like illness.

    The DLNM analysis revealed that the duration of the effect of temperature and pressure on the ILI incidence was 4 weeks. The median temperature and pressure were 18.97 °C and 100.5 kpa, respectively, adjusted for wind speed, precipitation, humidity, and PM2.5. The relationship between temperature and pressure and ILI incidence at different percentiles was analyzed using percentage points set as the 1st, 10th, 90th, and 99th percentiles.

    The relative risk (RR) of ILI incidence was elevated as temperature decreased. RR was highest (1.60, 95% CI: 1.37−1.86) when the temperature decreased to 0.5 °C at 0-week lag (Figure 1A). The cumulative exposure of temperature from 0 to 4 weeks indicated an increased RR with decreased temperature. RR was highest when the temperature dropped to 0.5 °C (4.63, 95% CI: 3.63−5.90; Figure 1C). As shown in Figure 2, lower temperatures (P1 = 2.85 °C, P10 = 8.62 °C) could increase the risk of ILI incidence, while higher temperatures (P90 = 27.52 °C, P99 = 29.21 °C) could increase the risk of ILI incidence for lag 0 to 3 weeks (Supplementary Table S3, available in www.besjournal.com). Consistent with our findings, low temperature is reportedly a risk factor for influenza at different lag periods[4]. Low temperature can promote the spread of influenza by prolonging the survival of the virus, increasing crowd contact, and affecting host immunity[5]. In addition, experiments performed with guinea pigs at different temperatures demonstrated that low temperature enhanced the transmission efficiency of influenza virus in the air. Breathing cold air could slow mucociliary clearance and thereby encourage the transmission of influenza virus in the respiratory tract. Moreover, low temperature may make the host more susceptible to influenza virus infection or prolong the course of disease after infection[6].

    Figure 1.  Associations of influenza-like illness (ILI) incidence with temperature and pressure. Lag-response associations of ILI incidence with temperature (A) and pressure (B). Cumulative associations of ILI incidence with temperature (C) and pressure (D).

    Figure 2.  Effects of temperature (left) and pressure (right) on influenza-like illness (ILI) incidence estimated at the extreme values of P1, P10, P90, and P99.

    Table S3.  The RR values of the effect of temperatures at each percentile node on ILI at lag of 0–4 weeks

    Lag weeksValuesP1 = 2.85 °CP10 = 8.62 °CP90 = 27.52 °CP99 = 29.21 °C
    0RR1.529*1.369*0.667*0.617*
    (95% CI)1.350–1.7321.276–1.4690.601–0.7410.542–0.701
    1RR1.411*1.275*0.763*0.725*
    (95% CI)1.312–1.5181.225–1.3280.721–0.8080.676–0.778
    2RR1.303*1.188*0.873*0.853*
    (95% CI)1.253–1.3541.164–1.2120.852–0.8950.827–0.879
    3RR1.202*1.106*0.9991.003
    (95% CI)1.129–1.2811.069–1.1460.944–1.0560.937–1.074
    4RR1.1101.0311.143*1.180*
    (95% CI)0.991–1.2420.967-1.0991.03–1.2681.039–1.34
    0–4RR3.750*2.366*0.508*0.451*
    (95% CI)3.087–4.5562.139–2.6160.448–0.5750.387–0.526
      Note. *P < 0.05. ILI: influenza-like illness.

    At 1−4 weeks lag, pressure had a positive effect on ILI incidence. The RR was highest in the second week (1.66, 95% CI: 1.36−2.02; Figure 1B). Data of the cumulative exposure to pressure from 0 to 4 weeks showed that the RR increased with the increment of pressure. Pressure > 100.5 kpa increased the risk of ILI incidence. RR was highest when the pressure was 102.0 kpa (6.64, 95% CI: 3.74−11.80; Figure 1D). Lower pressures (P1 = 99.92 kpa, P10 = 100.11 kpa) did not influence ILI incidence, while general high pressure (P90 = 101.06 kpa) could increase the risk of ILI incidence for lag 1 to 2 weeks (Figure 2). Extremely high pressure (P99 = 101.42 kpa) could increase the risk of ILI incidence for lag 1 to 4 weeks (Supplementary Table S4, available in www.besjournal.com). Very few studies have examined the effects of atmospheric pressure on influenza. Some studies found that high pressure is associated with influenza. Pressure is an important determinant of weather. High pressure during winter can result in contraction of local blood vessels in the nasal cavity and small cracks in the mucous membrane. These changes can allow respiratory viruses to easily invade the body. In addition, because the inverse temperature layer of the atmosphere is relatively stable, air pollutants do not easily diffuse and self-purify. This favors the replication and spread of viruses, resulting in numerous viruses in the atmosphere[7]. High pressure often means a sunny day. The resulting increase in outdoor activities increases the risk of virus infection[8].

    Table S4.  The RR values of the effect of pressures at each percentile node on ILI at lag of 0–4 weeks

    Lag weeksValuesP1 = 99.92 kpaP10 = 100.11 kpaP90 = 101.06 kpaP99 = 101.42 kpa
    0RR0.9881.0251.0131.052
    (95% CI)0.816–1.1960.923–1.1370.913–1.1230.923–1.199
    1RR0.9570.9711.126*1.247*
    (95% CI)0.813–1.1270.885–1.0661.035–1.2251.125–1.383
    2RR0.9620.9411.107*1.262*
    (95% CI)0.841–1.10.873–1.0151.04–1.1781.167–1.366
    3RR0.9800.9231.0481.191*
    (95% CI)0.828–1.1590.842–1.0130.962–1.141.075–1.32
    4RR0.9880.9071.0411.146*
    (95% CI)0.813–1.2020.815–1.0080.941–1.1511.018–1.29
    0–4RR0.8800.7841.375*2.261*
    (95% CI)0.564–1.3720.613–1.0021.116–1.6951.783–2.868
      Note. *P < 0.05.

    The optimal SARIMA (1, 0, 0) (1, 1, 0)52 structured on our previous research results were used to build the SARIMAX models with meteorological factors as covariates[9]. The cross-correlation function was used to investigate relationships of meteorological factors and PM2.5 with ILI incidence. All the variables, except precipitation and relative humidity, were significantly correlated with ILI incidence (Supplementary Figure S2, available in www.besjournal.com). These findings were consistent with the preceding Spearman’s correlation analysis results. A total of four SARIMAX model parameters with covariates were significant, with decreased AIC values. The temperature lagged 1-week and 4-week (β = −0.276, P = 0.008; β = 0.216, P = 0.038), wind speed lagged 1-week (β = −0.672, P = 0.013), and pressure lagged 4-week (β = −1.802, P = 0.002) influenced ILI incidence (Table 1). SARIMAX (1, 0, 0) (1, 1, 0)52 with temperature at lag 1-week, wind speed at lag 0-week, and pressure at lag 4-week was the optimal model with the lowest AIC and BIC values. The SARIMAX model displayed a good fitting effect on the prediction of ILI incidence (Supplementary Table S5, Supplementary Figure S3, available in www.besjournal.com).

    Figure S2.  Correlations with pre-whitened meteorological factors and PM2.5. CCF: Cross-correlation Function

    Table 1.  SARIMAX models of ILI incidence with different meteorological factors and PM2.5

    ModelMeteorological factorsAICBIC
    VariablesLagβSE (β)P
    SARIMAX (1, 0, 0) (1, 1, 0)52Temperature00.0420.1050.6882237.4002252.990
    SARIMAX (1, 0, 0) (1, 1, 0)521−0.2760.1040.008*2224.4202240.000
    SARIMAX (1, 0, 0) (1, 1, 0)5220.0230.1040.8262222.7902238.360
    SARIMAX (1, 0, 0) (1, 1, 0)523−0.1880.1040.0702214.5202230.070
    SARIMAX (1, 0, 0) (1, 1, 0)5240.2160.1040.038*2208.3402223.880
    SARIMAX (1, 0, 0) (1, 1, 0)52Wind speed0−0.6720.2700.013*2231.4102247.000
    SARIMAX (1, 0, 0) (1, 1, 0)5210.0570.2720.8352231.3302246.910
    SARIMAX (1, 0, 0) (1, 1, 0)522−0.0670.2710.8052222.7802238.350
    SARIMAX (1, 0, 0) (1, 1, 0)5230.4420.2710.1032215.1502230.700
    SARIMAX (1, 0, 0) (1, 1, 0)524−0.3610.2750.1892210.9102226.450
    SARIMAX (1, 0, 0) (1, 1, 0)52Pressure0−0.0590.5940.9212237.5502253.140
    SARIMAX (1, 0, 0) (1, 1, 0)5210.4240.6020.4812230.8802230.880
    SARIMAX (1, 0, 0) (1, 1, 0)5221.1010.5890.0622219.3702234.930
    SARIMAX (1, 0, 0) (1, 1, 0)523−0.1020.5940.8642217.7602233.310
    SARIMAX (1, 0, 0) (1, 1, 0)524−1.8020.5860.002*2203.3002218.840
    SARIMAX (1, 0, 0) (1, 1, 0)52PM2.500.1220.0830.1432235.4202251.010
    SARIMAX (1, 0, 0) (1, 1, 0)521−0.0160.0830.8492231.3402246.910
    SARIMAX (1, 0, 0) (1, 1, 0)5220.0220.0830.7952222.7702238.340
    SARIMAX (1, 0, 0) (1, 1, 0)5230.0520.0830.5342217.4002232.960
    SARIMAX (1, 0, 0) (1, 1, 0)524−0.1030.0830.2142211.0902226.640
    SARIMAX (1, 0, 0) (1, 1, 0)52Temperature1−0.1950.1090.0742198.7402225.950
    40.1290.1050.222
    Wind speed0−0.6000.2670.024*
    Pressure4−1.5040.5960.012*
    SARIMAX (1, 0, 0) (1, 1, 0)52Temperature1−0.2110.1090.049*2198.2302221.550
    Wind speed0−0.5960.2670.026*
    Pressure4−1.6720.5810.004*
      Note. *: P < 0.05. AIC: Akaike Information Criterion; BIC: Bayesian information criterion.

    Table S5.  Comparison between SARIMA model and SARIMAX model of ILI incidence

    ModelβPLjung-Box Q TestMAERMSEMAPE
    StatisticsP
    SARIMA (1, 0, 0) (1, 1, 0)5221.8220.1495.2307.90533.698
    AR10.886< 0.001
    SAR1−0.607< 0.001
    SARIMAX (1, 0, 0) (1, 1, 0)521.9000.1685.2227.88533.386
    AR10.865< 0.001
    SAR10.0970.065
    Temperature lag1−0.2110.049
    Wind speed lag0−0.5960.026
    Pressure lag4−1.6720.004

    Figure S3.  The observed and predicted value of influenza-like illness (ILI) incidence in Mississippi.

    Wind speed was also included in the model as a covariate, in addition to the temperature and pressure variables, which were most strongly correlated with ILI. Wind speed is a common meteorological variable. However, its correlation with disease has not been studied extensively. We observed that the RR of ILI incidence increased significantly with increased wind speed at the 0-week lag. However, the incidence gradually decreased at wind speeds above 4 m/s. Wind can resuspend bacteria from soil or plant surfaces[10]. This increases the chances that the influenza virus will be redrawn into the respiratory tract. Based on common resuspension, higher wind speeds may dilute local bacterial concentrations[7,10], thus reducing the infection of influenza.

    DLNMs and a 10-year time-series with covariates were respectively adopted in the correlation analysis and prediction. This enabled the simultaneous evaluation of the real-time and delayed effect of meteorological variation by capturing the true incidence of ILI incidence as much as possible. Nonetheless, several limitations need to be acknowledged. First, only the PM2.5 air pollution factor was considered when controlling for confounding factors. In subsequent studies, more air pollution factors should be included in the model and more mixed effects of meteorological factors and air pollution factors should be assessed. Second, the study was performed in a single subtropical city. The findings may not be generalizable to the regions with different climatic conditions.

    In conclusion, lower temperatures and higher pressures were risk factors of ILI incidence in Mississippi. Greater attention should be paid to ILI incidence during periods of cold and in the winter when the air pressure is high. These findings will be helpful for establishing an early warning system, strengthening active prevention and intervention of disease risk factors, and reducing the burden of disease in society.

    Declarations of Interest None

参考文献 (10)
补充材料:
22178Supplementary Materials.pdf

目录

    /

    返回文章
    返回