Short-Term Association Between Air Pollution and Infectious Disease Spectrum in Shanghai, China: A Time-Series Study From 2013 To 2019

Research Article

Austin J Public Health Epidemiol. 2024; 11(3): 1169.

Short-Term Association Between Air Pollution and Infectious Disease Spectrum in Shanghai, China: A Time-Series Study From 2013 To 2019

Yihan Lin1#; Hao Meng2#; Yong He2; Wenzhuo Liang1; Yiran Niu1; Zhenliang Liu1; Ziying Wang2; Yuan Lei1; Yangyang Tian1; Shiyang Chang1*

¹Department of Histology and Embryology, College of Basic Medical, Hebei Medical University, China

²Department of Pathogenic Biology, College of Basic Medicine, Hebei Medical University, China

*Corresponding author: Shiyang Chang, Department of Histology and Embryology, College of Basic Medical, Hebei Medical University, Shijiazhuang 050017, China. Tel./Fax: +86-0311-86266084 Email: 18901410@hebmu.edu.cn

#These authors have been equally contributed to this article.

Received: August 13, 2024 Accepted: August 30, 2024 Published: September 09, 2024

Abstract

Background: Epidemiological evidence on the association between air pollution and the risks of infectious diseases remained largely lacking. We aimed to examine associations of exposures to fine Particulate Matter (PM2.5) and ozone (O3) with risks of national notifiable infectious diseases in a mega city, shanghai in China.

Methods: We constructed a double-pollutant model for each air pollutant, applying a time-series analysis incorporating both single and Distributed Lag Model (DLM) separately to model the exposure-lag-response relationship with a total of 43 national Notifiable Infectious Diseases (NNIDs) during 2013 to 2019. The model was adjusted for seasonality and log-term trend, mean temperature, relative humidity, and other air pollutant. Analysis was further conducted for NNIDs categories and specific diseases.

Results: The study included 661,267 NNIDs cases. Exposures to PM2.5 and O3 were associated with increased risks of NNIDs but were not associated with the same categories. Each 10 μg/m3 increase in O3 was associated with an increased risk of total NNIDs (Relative Risk [RR] lag 1 month: 1.29, 95% Confidence Interval [CI]: 1.02 to 1.65), vaccine preventable disease (RR lag 1: 1.75, 95% CI: 1.02 to 3.01) and sexually transmitted and bloodborne diseases (RR at lag 2: 1.12, 95% CI: 1.00 to 1.26), while the association for PM2.5 remained inconclusive.

Conclusion: These findings suggested substantial infectious disease burden was associated with exposures to ambient air pollutants, emphasizing the urgent need a complete picture of association between air pollution and notifiable infectious diseases and comprehensive evaluation of the relevant disparity among spectrum of disease.

Keywords: Infectious diseases; Air pollution; Ozone; Fine particulate matter; Time-series study; Distributed lag model

Introduction

China has undergone a rapid epidemiological transition with remarkable progress in the control of infectious diseases. Yet, possibly due to significant reductions in disease burden, infectious diseases are often overlooked as causes of morbidity and mortality in China, and the assessment of climate and air pollution effects and its effect disparity related to infectious diseases [6,18]. In the past few years, the outbreak of a variety of emerging infectious diseases, such as COVID-19 and monkeypox, has raised a new challenge for human health and emerging attention on the climate impact on infectious diseases [3,11,21,22,28].

Ambient air pollutants, such as fine Particulate Matter (PM2.5) and Ozone (O3), could potentially elevate the presence of bacteria, viruses, or other pathogens within the atmosphere. They may also function as an immunosuppressive agent, thereby compromising the typical immune defenses of the human body [23]. However, current epidemiological evidence of the relationship between air pollutants and infectious diseases remains limited and inconclusive, posing a great challenging to draw a reliable conclusion from existing research. Moreover, comprehensive reports of comparison and disparity in the as sociation between air pollution and spectrum of infectious diseases were not identified in China and other countries [7,23]. In this context, the successive infectious disease surveillance system is an opportunity to provide a complete picture of association between air pollution and notifiable infectious diseases and comprehensive evaluation of the relevant disparity among spectrum of disease in the past decade.

To our knowledge, this is the first study to report comprehensively the short-term effect of air pollution on a wide range of notifiable infectious diseases due to 43 causes and evaluate the disparity in association by specific category. The identification of the potential disparity in infectious disease burdens due to air pollution would provide direction for the precise implementation of prevention and control measures.

Materials and Methods

Study Design and Infectious Diseases Data

The study is a time-series analysis using secondary data conducted in Shanghai, a megacity in China, during 2013 and 2019. Ethnical approval was not applicable to our research since the data collected in this study is secondary data without any personal information.

Monthly National Notifiable Infectious Diseases (NNIDs) data was collected from the surveillance system, detailed description has been published elsewhere [6]. A total of 43 National Notifiable Infectious Diseases (NNIDs) were included in this study, and were divided into seven categories following the previous categorization approach [9]. Specifically, I. Vaccine Preventable Diseases (11 diseases): This category encompasses seasonal influenza, rubella, pertussis, mumps, measles, hepatitis A, B and D, neonatal tetanus, poliomyelitis and diphtheria. II. Bacterial Diseases (4 diseases): including tuberculosis, scarlet fever, meningococcal meningitis, and leprosy. III. Gastrointestinal and Enterovirus Diseases (5 diseases): This group consists of diseases primarily affecting the gastrointestinal system, such as typhoid and paratyphoid, infectious diarrhea, Hand, Foot, and Mouth Disease (HFMD), dysentery, and acute hemorrhagic conjunctivitis. IV. Sexually Transmitted and Bloodborne Diseases (4 diseases): This category includes syphilis, gonorrhea, HIV/AIDS, and hepatitis C. V. Vectorborne Diseases (7 diseases): This group covers typhus, schistosomiasis, malaria, kala-azar, Japanese encephalitis, dengue, and filariasis. VI. Zoonotic Diseases (9 diseases): brucellosis, hepatitis E, hydatid disease, rabies, anthrax, leptospirosis, H5N1, H7N9, and Severe Acute Respiratory Syndrome (SARS). VII. Quarantinable Diseases (3 diseases): This category encompasses hemorrhagic fever, cholera, and plague.

Air Pollutants and Weather Variables

The two air pollutants included in this study are fine particulate matter (PM2.5) and ozone (O3). The unit of PM2.5 and O3 prediction is μg/m3, aligning with the China ambient air quality standards (GB3095-2012). Monthly average PM2.5 concentrations at surface level were downloaded from a nationwide PM2.5 dataset, with a spatial resolution of 10 km [3,26]. The data set is part of the Tracking Air Pollution in China (TAP, http://tapdata.org.cn/) project. The details of the PM2.5 prediction model have been documented elsewhere [13,15]. Briefly, TAP is a publicly available database with a high temporal and spatial resolution, and the enhanced performance through the integration of multisource-fusion data and machine learning algorithms. Initially, a comprehensive dataset comprising ground PM2.5 measurements from monitoring stations, satellite-derived Aerosol Optical Depth (AOD), meteorological parameters, land use characteristics, population figures, and elevation data, along with information from the Weather Research and Forecasting/Community Multiscale Air Quality Modeling System (WRF/CMAQ), was harmonized into a unified 10 km grid. Subsequently, the estimation of PM2.5 concentration in TAP products was predicted through a two-stage machine learning framework employing a synthetic minority oversampling technique in conjunction with a tree-based gap-filling method. The cross-validation of the prediction model yielded a range of 0.80 to 0.88, indicating comparable performance with existing studies [26].

Maximum 8 h average O3 concentration predictions were collected from TAP dataset, predicted using the three-stage random forest model [27. The three-stage O3 prediction model incorporates a comprehensive set of data sources, including ground measurements in the reference state, CMAQ simulations, Ozone Monitoring Instrument (OMI) satellite O3 profiles (PROFOZ), MERRA-2 meteorological parameters, MODIS Normalized Difference Vegetation Index (NDVI), and National Centers for Environmental Information (NCEI) annual night light data. In the initial stage, two sets of maximum 8-hour average O3 concentration predictions were generated, one incorporating satellite data and the other without, aimed at addressing gaps arising from missing satellite retrievals in the subsequent model. The O3 prediction model that excluded satellite data provided full spatial coverage. In the second stage, we employed an elastic-net regression model to merge random forest predictions from both datasets, ensuring comprehensive O3 predictions. To enhance prediction accuracy, a third-stage model was developed to predict the spatiotemporal distribution of the difference between maximum 8-hour average O3 measurements and random forest predictions, utilizing kriging interpolations. These predicted residuals were then integrated into the second-stage predictions, yielding the final predictions. The 5-fold cross-validation predictions of the O3 prediction model demonstrated an R2 value of 0.70 when compared to ground measurements.

Monthly meteorological data during 2013 and 2019, including average temperature (°C), relative humidity (%), was obtained from the National Meteorological Data Sharing Center (http://data.cma.cn/).

Statistical Analysis

Time Series Regression (TSR) with generalized linear model was applied to explore the short-term effect of individual weather variable on NNIDs categories [16]. To allow for the overdispersion of the NNIDs counts, a Quasi-Poisson model is selected. TSR is widely used for mathematical modelling in environmental epidemiology, as it measures short-term effect (which is the association between monthly variation in exposure and outcome in this study), socioeconomic and demographic levels therefore are assumed to be constant over neighboring months. To control for the long-time trend and seasonality of NNIDs, alternative choices of time adjustment have been performed and compared, including linear trend, time interactions, Fourier terms, and different splines with varied degrees of freedom (df). A Natural Cubic B-Spline (NCS) of time with 8 df per year was selected as the best-fit approach in our analysis (Supplementary Figure 1). For modelling the relationship between air pollutants and NNIDs, two modeling approaches including a single lag model and a Distributed Lag Model (DLM) were employed separately, as the health effect of pollutant variability is usually linear and delayed [1,4]. In particular, the DLM accounts for the impacts of other lag periods and provides cumulative exposure-lag-response associations over multiple lag units. Specifically, a cross-basis function for both exposure and lag dimensions assuming linear relations were introduced into the model [5. For non-infectious disease, a maximum lag of 21 day is commonly used (which approximates to lag1 in this study) [12]. However, given the more complicated casual pathway for infectious disease, wide range of lag months of 2 has been considered, also based on recommendation from previous literatures on infectious immune period, i.e., maximum lag up to 6 months [16,17]. We did the analysis for NNIDs by different categories as well as by specific causes. However, only subgroups with a sample size exceeding 5000 were considered eligible for analysis, ensuring adequate statistical power [2].

Double-pollutant models are applied for each pollutant, adjusting for another pollutant as well as other time-varying weather variables (mean temperature, relative humidity) [5]. The adjustment for each confounding variable is implemented via a NCS with 3 df of moving average of the covariate over the lag period. The effect estimate is reported as Relative Risk (RR) with its 95% Confidence Interval (CI), representing the morbidity risk changes per 10 μg/m3 in PM2.5 or O3 at each lag month, after adjusting for other lagged exposures.

Sensitivity Analysis

Sensitivity analyses were performed to check the robustness of the analysis. An extended 3-month lag period to explore a wider range of relationship pattern and lag durations. Additionally, single-pollutant models were employed for comparison with results adjusted for other pollutants. All the analyses were performed in R software (version 4.2.1; https://www.rproject.org/) with “dlnm” package. The significance level was set at 0.05.

Results

Characteristics of NNIDs and Air Pollutants

There were 661,267 incident NNIDs reported in Shanghai from January 2013 to December 2019 (Table 1). The incident cases predominately consisted of vaccine preventable diseases (93,134 cases, 14%), bacteria diseases (73,851, 11%), gastrointestinal and enterovirus diseases (351,464, 53%) and sexually transmitted and bloodborne diseases (137,036 cases, 20%). There were also 447 incident cases for vector borne diseases, 5,300 for zoonotic diseases, and 35 for quarantinable diseases. The time-series plots showed seasonal patterns for each NNIDs category, along with a decreasing trend over time except for vector borne and zoonotic diseases (Supplementary Figure 2).