Skip to main content

Impact of outdoor air pollution on the incidence of pertussis in China: a time-series study

Abstract

Background

The increasing number of pertussis cases worldwide over the past two decades has challenged healthcare workers, and the role of environmental factors and climate change cannot be ignored. The incidence of pertussis has increased dramatically in mainland China since 2015, developing into a serious public health problem. The association of meteorological factors on pertussis has attracted attention, but few studies have examined the impact of air pollutants on this respiratory disease.

Methods

In this study, we analyzed the relationship between outdoor air pollution and the pertussis incidence. The study period was from January 2013 to December 2018, and monthly air pollutant data and the monthly incidence of patients in 31 provinces of China were collected. Distributed lag nonlinear model (DLNM) analysis was used to estimate the associations between six air pollutants and monthly pertussis incidence in China.

Results

We found a correlation between elevated pertussis incidence and short-term high monthly CO2 and O3 exposure, with a 10 μg/m3 increase in NO2 and O3 being significantly associated with increased pertussis incidence, with RR values of 1.78 (95% CI: 1.29-2.46) and 1.51 (95% CI: 1.16-1.97) at a lag of 0 months, respectively. Moreover, PM2.5 and SO2 also played key roles in the risk of pertussis surged. These associations remain significant after adjusting for long-term trend, seasonality and collinearity.

Conclusions

Overall, these data reinforce the evidence of a link between incidence and climate identified in regional and local studies. These findings also further support the hypothesis that air pollution is responsible for the global resurgence of pertussis. Based on this we suggest that public health workers should be encouraged to consider the risks of the environment when focusing on pertussis prevention and control.

Peer Review reports

Introduction

Pertussis has been resurging worldwide in the past two decades. Pertussis epidemics have been reported increasingly in East Asia (ie, mainland China, India), Europe (ie, United Kingdom), North America (ie, America and Canada), and Oceania (ie, Australia) since 2008 [1,2,3,4]. One study indicated that in 2012, the rate of pertussis in the U.S. was higher than in any previous year since the vaccine was first implemented [2]. A burst of pertussis incidence also appeared in China in recent years and has lasted till today.

Pertussis is a highly contagious respiratory disease caused by the bacterium Bordetella pertussis [5]. After an incubation period of about 1 to 3 weeks (usually 7 to 10 days), patients develop apparent symptoms, from a runny dry cough to a progressive spasmodic cough, the latter lasting for up to 2 months [6]. The pathophysiologic mechanism involves inflammation and spasm of the respiratory tract. The pathogen adheres to the respiratory epithelium and releases toxins causing an inflammatory response that leads to swelling of the respiratory mucosa, increased secretions and decreased mucus membrane secretion. This eventually causes narrowing and obstruction of the respiratory tract, shortness of breath, and coughing. Respiratory muscle spasm will exacerbate this phenomenon and prolong the onset [7,8,9].

The epidemiology of pertussis appears to be the result of a combination of factors such as un-vaccination, reduced vaccine efficacy, gene mutation in Bordetella pertussis, environmental factor, and climate change [10, 11]. Air pollution has been recognized as a pervasive public health threat [12]. Previous studies have reported the relationship between weather parameters and pertussis [6, 13, 14]. Through regression tree models and DLNMs, Zhang et al. revealed that temperature change, high temperature, and relative humidity increased pertussis risk, temperature found to be a key driver of pertussis infection. Interestingly, the association between pertussis infection and meteorological factors was stronger in people over 10 years old [14]. By studying the relationship between meteorological parameters and the risk of pertussis in Chongqing, China from 2004 to 2018, Wang et al., have shown that for each unit increase in monthly average temperature, the number of monthly pertussis cases increased by 19.496% (95%CI 2.368-39.490%), and meteorological factors such as air pressure and precipitation were also linked to pertussis [6]. However, the small size of the target population, or the assessment of the effect of a single meteorological parameter on pertussis in these studies, was not sufficient to compare the changes of disease incidence and show the long-term role of climate change in the spread of pertussis. What is more, very few studies have explored the impact of air pollutants on this respiratory disease. Therefore, it is necessary to conduct a study that also considered common air pollutants such as (PM2.5, PM10, SO2, CO, NO2, and O3), to comprehensively elucidate the long-term quantitative association between atmospheric pollutants and pertussis incidence to address the growing threat posed by pertussis to public health.

In this observational study, we aimed to assess changes in pertussis epidemiology, using DLNM models, to explain the complex causes of pertussis surges and the associations between six air pollutants and pertussis, and discuss how to manage these factors to prevent future epidemics and provide early warning information and reference for pertussis epidemic.

Materials and methods

Data collection

We obtained pertussis cases data from the China Public Health Science Data Center (https://www.phsciencedata.cn/Share/en/index.jsp.). Population data of the corresponding periods were obtained from the National Bureau of Statistics of the People’s Republic of China (http://data.stats.gov.cn). All notified cases were diagnosed according to detailed instructions issued by the National Health Commission in 2007 (http://www.nhc.gov.cn/wjw/s9491/201410/52040bc16d3b4eecae56ec28b3358666.shtml). The monthly average air pollutant concentrations (PM2.5, PM10, SO2, CO, NO2, and O3) and meteorological data (sunshine, mean wind speed, precipitation, relative humidity, and mean temperature) of 31 provinces in China from 2013 to 2018 were obtained from 393 National Air Quality Monitoring Stations (https://www.epmap.org/) and 839 meteorological monitoring stations (https://data.cma.cn/), respectively (Fig. 1).

Fig. 1
figure 1

Geographical location of air quality monitoring stations (during 2013-2018) and meteorological monitoring stations (2004-2018) in China. Air quality monitoring stations (N=393); meteorological monitoring stations (N=839)

Statistical analysis

We analyzed the association between monthly trends in atmospheric pollutants and monthly incidence of pertussis. Since meteorological conditions can affect the concentration of air pollutants, we incorporated five common meteorological conditions into the DLNM model to improve the accuracy of the results. The effects of air pollutants on pertussis were evaluated using DLNM models.

DLNM were based on traditional model ideas such as generalized additive models and distributed lag linear models [15]. The main idea was to add a lag dimension to the exposure-response relationship through a cross basis function. In this way, the change distribution of its effect in the independent variable dimension and the lag dimension was described at the same time.

Both single-pollutant and multi-pollutant models were used to analyze the effect of air pollutant variables on pertussis incidence. The single-pollutant model considered not only air pollutants, but also long-term trend, meteorological conditions, quartiles of average incidence, and incidence in the previous month. The final model looked like this:

$$\mathrm{Model}\;1:\;\mathrm{Log}\;\lbrack\mathrm E(\mathrm{Yt})\rbrack\;=\mathrm\beta\;+\;\mathrm{cb}\;(\mathrm{Kt},\;2,\;\mathrm\beta1)+\mathrm S1\;(\mathrm x)\;+\;\mathrm S2\;(\mathrm z)\;+\;\mathrm S3\;(\mathrm m)\;+\;\mathrm S5\;(\mathrm{Region})\;+\;\mathrm S6\;(\mathrm{Season})$$

Where t is the observation month; [E(Yt)] is the observed pertussis cases in year Yt, and β is the intercept of the whole model; cb (Kt, 2, and β1) is the cross basis of K, and β1 is an estimate of the effect of K in a given lag month t, with the maximum delay month set to 2; K is one of the four air pollutants. S1-3 are "average sunlight", "average temperature", and "average wind speed" representing meteorological factors associated with pertussis incidence; "Region" and "Season" represents the covariates; month is the ordinal variable of January in a year; and S() is the penalized spline function. In this study, we used cubic spline functions S1-3 to adjust for confounders in the model, and S5 (Region) and S6 (Season) to adjust for monthly confounders.

The multi-pollutant model took air pollutants and meteorological conditions (sunshine, wind speed, and mean temperature) associated with pertussis incidence into account, other adjustment variables in the model were long-term trend, seasonal trend, quartile arrays of mean incidence, and the incidence in the previous month. The final model looked like this:

$$\mathrm{Model}:\;\mathrm{Log}\;\lbrack\mathrm E(\mathrm{Yt})\rbrack=\beta+\mathrm{cb}(\mathrm{Kt},\;2,\;\beta1)+\mathrm S1(\mathrm x)+\mathrm S2(\mathrm z)+\mathrm S3(\mathrm m)+\mathrm S5(\mathrm{Region})+\mathrm S6(\mathrm{Season})+\mathrm S7(\mathrm{lag}.\mathrm{value}1)$$

where t is the observation month; [E(Yt)] is the observed pertussis cases in year Yt, and β is the intercept of the whole model; cb (Kt, 2, and β1) is the cross basis of K, and β1 is the estimated impact of K at a specific lag month t, with the maximum delay month set to 2; K is the sum of the four air pollutants. S1-3 are "average sunshine," "average temperature," and "average wind speed," representing meteorological factors associated with pertussis incidence; "region", "season" and "lag.value1" covariates; "lag.value1" represents the average incidence in the previous month; and S() is a penalized spline function. In this study, we used cubic spline functions S1-3 to adjust for confounders in the model, and S5 (Region) and s6 (Season) to adjust for monthly confounders.

The construction of the basis function used the crossbasis () function in the “dlnm” package. The exposure-response dimension fitting was performed using natural cubic splines with node positions set at the 10th, 75th and 90th percentiles. The exposure-lag dimension fitting also used the natural cubic spline function, the lag period was selected to be seven months, and the node positions were set using the logarithmic equal spacing method provided by the package “dlnm” (version 4.1.3, https://cran.r-project.org/web/package/dlnm/index.html) [16]. Second, according to the AIC+BIC criterion we selected the degree of freedom, and chose the degree of freedom with the smallest AIC+BIC value to build the final model.

Results

Descriptive analysis

From 2013 to 2018, 79 757 cases were reported in China. On average, the monthly mean concentration of PM2.5, PM10, SO2, NO2, O3 and CO were 50.65 μg/m3, 87.67 μg/m3, 24.5 μg/m3, 33.19 μg/m3, 84.36 μg/m3, and 0.98 mg/m3, respectively (Table 1).

Table 1 Descriptive statistics for monthly pertussis cases and air pollution concentrations in China, 2013-2018

DLNM model

Single-pollutant models showed that four air pollutants (PM2.5, SO2, NO2, and O3) were associated with pertussis incidence from 2013 to 2018 (Fig. 2), with relative risks ranging from 0.59-1.27 for PM2.5, 0.46-1.79 for SO2, 0.57-3.72 for NO2, and 0.98-2.61 for O3. In the multi-pollutant models (Fig. 3), the relative risk of PM2.5 was 0.43-1.38, SO2 was 0.21-3.82, NO2 was 0.33-2.97, and O3 was 0.85-1.99. The maximum RR for PM2.5 was 1.38 with exposure to 10 μg/m3 PM2.5 at lag 2 months. The maximum RR for SO2 was 3.82 with exposure to 195μg/m3 SO2 at lag 1 months. The maximum RR for NO2 was 2.97 with exposure to 8 μg/m3 NO2 at lag 0 months. The maximum RR for O3 was 1.99 with exposure to 15 μg/m3 O3 at lag 1.6 months.

Fig. 2
figure 2

Contour plots of exposure-response relationships between pertussis incidence and four air pollutants in the single-pollutant model. A-D Four air pollutants associated with (PM2.5, SO2, NO2 and O3) pertussis incidence from 2013 to 2018. The reference level was set to the median value of the corresponding variable. The Y-axis represents the lag period from 0 to 2 months. The e-axis represents the range of observations for each variable. RR stands for relative risk, red stands for RR>1, white stands for RR=1, and blue stands for RR<1

Fig. 3
figure 3

Contour plots of exposure-response relationships between pertussis incidence and four air pollutants in the multi-pollutant model. A-D Four air pollutants associated with (NO2, O3, PM2.5 and SO2) pertussis incidence from 2013 to 2018. The reference level was set to the median value of the corresponding variable. The Y-axis represents the lag period from 0 to 2 months. The X-axis represents the range of observations for each variable. RR stands for relative risk, red stands for RR>1, white stands for RR=1, and blue stands for RR<1

The lag relationship between air pollutants and the incidence of pertussis

The results of DLNMs were shown in DLNMs, we used median PM2.5, median SO2, median NO2, and median O3 as references, and then RR of pertussis incidence with lag 0-2 months was calculated with the extremely low (2.5th), 25th, 75th and extremely high (97.5th) percentile of PM2.5, SO2, NO2, and O3, respectively.

The 3-D plots for multiple pollutants show that the RR values would increase when the NO2 and SO2 concentrations increase at a lag of 1 and 0 months (Fig. 4). Furthermore, single-pollutant plots presented in Fig. 5 showed that under the 75th quantile of NO2, RR for lag 0 months (RR = 1.15, 95% CI: 1.03-1.28) and RR for 2 months (RR = 1.12, 95% CI: 1.01-1.24) were significantly high. However, under extremely low (2.5th) ,25th and high (97.5th) NO2, no statistical significance in RR was observed; Under the 2.5th quantile of O3, the RR for lag 0 months (RR = 1.37, 95% CI: 1.02-1.83) and the RR for lag 2 months (RR = 1.69, 95% CI: 1.32-2.17) were significantly high. In addition, Under the 97.5th quantile of O3, the RR for lag 0 months (RR = 1.21, 95% CI: 1.00-1.44) and the RR for lag 2 months (RR = 1.32, 95% CI: 1.10-1.58) were significantly high. And under the 25th and 75th percentile of O3, no statistical significance in RR was observed. Under the 2.5th, 25th, and 75th percentile of PM2.5, no statistical significance in RR was observed; under the 2.5th, 25th quantile of SO2, RR for lag 0 months (RR = 1.60, 95% CI: 1.10-2.33), RR for lag 0 months (RR = 1.10, 95% CI: 1.00-1.21), were significantly high, respectively. Under the 75th percentile of SO2, no statistical significance in RR was observed.

Fig. 4
figure 4

3-D plot of relative risk of pertussis onset (RR) versus lag time for exposure to four air pollutants. A-D Four air pollutants associated with (NO2, O3, PM2.5 and SO2) pertussis incidence from 2013 to 2018. The X-axis represents the exposure concentration of the four pollutants,the Y-axis represents the lag period from 0 to 2 months, the Z-axis represents the relative risk of each variable (NO2, O3, PM2.5 and SO2)

Fig. 5
figure 5

Summary of relative risk curves for different lags between air pollutants and pertussis incidence in single-pollutant models from 2013 to 2018. The Y-axis represents the relative risk of each variable (NO2, O3, PM2.5 and SO2), the X-axis represents the lag period from 0 to 2 months. Lines represent means estimated using the DLNM single-pollutant model, shaded areas represent 95% confidence intervals

Multi-pollutant plots presented in Fig. 6 showed at the 2.5th, 25th, 75th and 97.5th percentile of NO2, the RR for lag 2 months (RR = 1.45, 95% CI: 1.07-1.98), RR for lag 2 months (RR = 1.26, 95% CI: 1.07-1.48), RR for lag 2 months (RR = 1.15, 95% CI: 1.04-1.27), and RR for 2 months (RR = 1.43, 95% CI: 1.09-1.89) were significantly high, respectively. In addition, under extremely low (2.5th) NO2, the risk was greatest at a lag of 0 months (RR = 1.78, 95% CI: 1.29-2.46). Under the 2.5th, 75th and 97.5th quantile of O3, the RR for lag 2 months (RR = 1.45, 95% CI: 1.17-1.80), RR for lag 2 months (RR = 1.13, 95% CI: 1.05-1.22) and the RR for lag 2 months (RR = 1.32, 95% CI: 1.22-1.57) were significantly high, respectively. Additionally, under extremely low (2.5th) O3, the risk was greatest at a lag of 1 months (RR = 1.51, 95% CI: 1.16-1.97). Under the 25th percentile of O3, no statistical significance was observed for RR. Lastly, Under the 2.5th, 25th, 75th and 97.5th percentile of PM2.5 and the 2.5th, 25th and 75th percentile of SO2, no statistical significance was observed for RR.

Fig. 6
figure 6

Summary of relative risk curves for different lags between air pollutants and pertussis incidence in multi-pollutant models from 2013 to 2018. The Y-axis represents the relative risk of each variable (NO2, O3, PM2.5 and SO2), the X-axis represents the lag period from 0 to 2 months. Lines represent means estimated using the DLNM multi-pollutant model, shaded areas represent 95% confidence intervals

Discussion

Our study mainly analyzed the association between monthly air pollution and monthly pertussis incidence in 31 provinces and regions in China from 2013 to 2018. Estimates of the percent change in pertussis cases with a 10-unit increase in air pollutants were based on a large dataset. Through the indicator, the epidemic trend of pertussis can be predicted, targeted prevention and control measures can be taken in advance, and health education and public health interventions for susceptible populations can be carried out.

Over the past 71 years, the pertussis epidemic in China has gone through four stages: high incidence, decline, low incidence, and resurgence. In addition, the incidence of some respiratory infectious diseases has increased over the past decade. Finally, because of the strong link between climate change and respiratory infectious diseases, we believe that climate has played an important role in the resurgence of pertussis.

The research focused on the association between air pollutants and pertussis incidence. A novel discovery of this research was the exposure-response relationship between the cumulative risk of NO2, O3, PM2.5 and SO2 and the incidence of pertussis. Considering there was no much literature reporting the association between air pollutants and pertussis, we compared the results of our study with other respiratory infectious diseases.

Environmental particulate matter may play an important role in the occurrence and development of respiratory diseases. Our research results revealed that NO2 increased the risk of pertussis and this was consistent with recent discovery. One possible mechanism of NO2 and other air pollutants was that, regardless of the duration of exposure, high levels of NO2 can irritate airways and increase susceptibility to respiratory infections [17, 18]. Compared with adults, children are at a higher risk because they breathe at higher rates and spend more time outdoors and their lungs and immune systems are immature.

Previous studies have shown that exposure to O3 in early life can impair lung function and growth [19,20,21]. In infant rhesus monkeys, sustained exposure to O3 caused structural abnormalities in bronchioles and increased susceptibility to antigen-induced airway remodeling [22]. We found significant evidence of spatial heterogeneity in the associations between air pollutants and morbidity across provinces and regions. Many factors may contribute to this variability, including different geographic distributions, long-term air pollution levels, population susceptibility, environmental policy and planning, and different study period lengths.

Our study found that short-term (0-15 lag days) exposure to high concentrations of PM2.5 increased the risk of pertussis infection, and some studies reported that short-term exposure to outdoor particulate matter (PM10 or PM2.5) have been found to reduce lung function and increase hospitalization rates [23,24,25,26]. A recent study assessed the association between ambient particulate matter exposure and chronic obstructive pulmonary disease, using logistic regression analyses and multiple linear regression models. The results provided evidence that exposure to higher PM10 or PM2.5 concentrations increased the risk of developing COPD and decreased the respiratory function [27]. Mortality rates from respiratory diseases have gradually decreased as air pollutant levels have improved [28, 29]. However, studies have shown that even lower levels of ambient air pollutants can have harmful effects on the respiratory system and lungs, suggesting that even at lower concentrations of pollutants, there is still a risk of disease [30].

As a potent lung irritant that can reduce lung function, SO2 has been identified as a cause of respiratory symptoms including cough and reduced lung function. Exposure to high levels of SO2 can cause severe airway damage. A study examined the relationship between nitrogen dioxide and sulfur dioxide concentrations and children’s lung function in Greece showed that long-term exposure to higher concentrations of air pollutants was associated with subclinical airway narrowing and slower increase of percent forced vital capacity (FVC%) in 8-10-year-old children [31]. Our results also suggest that higher concentrations of SO2 increase the risk of pertussis infection during the lag period of 0-30 days.

Over the past 30 years, the Chinese government has implemented policies to improve air quality and ameliorate the adverse effects of air pollutants on public health. Since 2013, air quality has been continuing to improve in most regions. However, this study showed that NO2, O3, PM2.5 and SO2 concentrations and their lagged effects were associated with the risk of pertussis, which indicated that high baseline levels of pollution and its impact on pertussis cases would persist.

Using air pollutants and case data from 2013 to 2018, We found a non-negligible association between short-term pollutant exposure and the development of pertussis. Our previous study analyzed the risk of pertussis after long-term (0-15 lag months) exposure to these major air pollutants. However, considering that it has an incubation period of 1-2 weeks, and even the duration of the disease lasts only 4-6 weeks, we adjusted the lag time to 0-2 lag months in the current study. Generally speaking, these findings further support the hypothesis that air pollution is responsible for the global resurgence of pertussis. Finally, we recommend encouraging public health workers to consider NO2 and O3 risks in response to the pertussis resurgence.

This study has several limitations. First, we could not obtain the data of air pollutants before 2013 because the national air quality monitoring network was only established at 2013. Second, on the basis of only monthly data we were unable to work out the acute effects of air pollutants on pertussis. Third, we were unable to extrapolate the effect of individual exposure levels on pertussis risk from air pollution across provinces. This study focused on the relationship between six air pollutants and pertussis incidence in China from the perspective of time lag effect. We need to further investigated the potential role of air pollutants in accelerating pertussis transmission in the future.

Availability of data and materials

All data generated or analysed during this study are included in this published article (and its supplementary information files).

Abbreviations

DLNM:

Distributed lag nonlinear model

References

  1. Almaw L, Bizuneh H. Pertussis outbreak investigation in Janamora district, Amhara Regional State, Ethiopia: a case-control study. Pan Afr Med J. 2019;34:65. https://doi.org/10.11604/pamj.2019.34.65.19612.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Cherry JD. Epidemic pertussis in 2012–the resurgence of a vaccine-preventable disease. N Engl J Med. 2012;367(9):785–7. https://doi.org/10.1056/NEJMp1209051.

    Article  CAS  PubMed  Google Scholar 

  3. Christie CD, Marx ML, Marchant CD, Reising SF. The 1993 epidemic of pertussis in Cincinnati. Resurgence of disease in a highly immunized population of children. N Engl J Med. 1994;331(1):16–21.

    Article  CAS  PubMed  Google Scholar 

  4. Güriş D, Strebel PM, Bardenheier B, Brennan M, Tachdjian R, Finch E, et al. Changing epidemiology of pertussis in the United States: increasing reported incidence among adolescents and adults, 1990–1996. Clin Infect Dis. 1999;28(6):1230–7.

    Article  PubMed  Google Scholar 

  5. Xu Z, Hu D, Luu LDW, Octavia S, Keil AD, Sintchenko V, et al. Genomic dissection of the microevolution of Australian epidemic Bordetella pertussis. Emerg Microbes Infect. 2022;11(1):1460–73. https://doi.org/10.1080/22221751.2022.2077129.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Wang Y, Xu C, Ren J, Zhao Y, Li Y, Wang L, et al. The long-term effects of meteorological parameters on pertussis infections in Chongqing, China, 2004–2018. Sci Rep. 2020;10(1):17235. https://doi.org/10.1038/s41598-020-74363-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Pinto MV, Merkel TJ. Pertussis disease and transmission and host responses: insights from the baboon model of pertussis. J Infect. 2017;74(Suppl 1):S114-s119. https://doi.org/10.1016/s0163-4453(17)30201-3.

    Article  PubMed  Google Scholar 

  8. Jiang W, Wei C, Mou D, Zuo W, Liang J, Ma X, et al. Infant rhesus macaques as a non-human primate model of Bordetella pertussis infection. BMC Infect Dis. 2021;21(1):407. https://doi.org/10.1186/s12879-021-06090-y.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Warfel JM, Beren J, Kelly VK, Lee G, Merkel TJ. Nonhuman primate model of pertussis. Infect Immun. 2012;80(4):1530–6. https://doi.org/10.1128/iai.06310-11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Rane MS, Rohani P, Halloran ME. Durability of protection after 5 doses of acellular pertussis vaccine among 5–9 year old children in King County Washington. Vaccine. 2021;39(41):6144–50. https://doi.org/10.1016/j.vaccine.2021.08.070.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Klein NP, Bartlett J, Rowhani-Rahbar A, Fireman B, Baxter R. Waning protection after fifth dose of acellular pertussis vaccine in children. N Engl J Med. 2012;367(11):1012–9. https://doi.org/10.1056/NEJMoa1200850.

    Article  CAS  PubMed  Google Scholar 

  12. Tainio M, Jovanovic Andersen Z, Nieuwenhuijsen MJ, Hu L, de Nazelle A, An R, et al. Air pollution, physical activity and health: A mapping review of the evidence. Environ Int. 2021;147:105954. https://doi.org/10.1016/j.envint.2020.105954.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Zhang Y, Bambrick H, Mengersen K, Tong S, Feng L, Zhang L et al. Using big data to predict pertussis infections in Jinan city, China: a time series analysis. Int J Biometeorol. 2020; 64(1).https://doi.org/10.1007/s00484-019-01796-w.

  14. Zhang Y, Bambrick H, Mengersen K, Tong S, Feng L, Liu G, et al. Association of weather variability with resurging pertussis infections among different age groups: A non-linear approach. Sci Total Environ. 2020;719:137510. https://doi.org/10.1016/j.scitotenv.2020.137510.

    Article  CAS  PubMed  Google Scholar 

  15. Gasparrini A, Armstrong B, Kenward MG. Distributed lag non-linear models. Stat Med. 2010;29(21):2224–34. https://doi.org/10.1002/sim.3940.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Gasparrini A. Distributed Lag Linear and Non-Linear Models in R: The Package dlnm. J Stat Softw. 2011;43(8):1–20.

  17. Liu Y, Ding H, Chang S-T, Lu R, Zhong H, Zhao N, et al. Exposure to air pollution and scarlet fever resurgence in China: a six-year surveillance study. Nat Commun. 2020;11(1):4229. https://doi.org/10.1038/s41467-020-17987-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Mahara G, Wang C, Yang K, Chen S, Guo J, Gao Q, et al. The Association between Environmental Factors and Scarlet Fever Incidence in Beijing Region: Using GIS and Spatial Regression Models. Int J Environ Res Public Health. 2016;13(11):1083.

  19. Holm SM, Balmes JR. Systematic Review of Ozone Effects on Human Lung Function, 2013 Through 2020. Chest. 2022;161(1):190–201. https://doi.org/10.1016/j.chest.2021.07.2170.

    Article  CAS  PubMed  Google Scholar 

  20. Yen Y-C, Yang C-Y, Ho C-K, Yen P-C, Cheng Y-T, Mena KD, et al. Indoor ozone and particulate matter modify the association between airborne endotoxin and schoolchildren’s lung function. Sci Total Environ. 2020;705:135810. https://doi.org/10.1016/j.scitotenv.2019.135810.

    Article  CAS  PubMed  Google Scholar 

  21. Mortimer KM, Tager IB, Dockery DW, Neas LM, Redline S. The effect of ozone on inner-city children with asthma: identification of susceptible subgroups. Am J Respir Crit Care Med. 2000;162(5):1838–45.

    Article  CAS  PubMed  Google Scholar 

  22. Schelegle ES, Miller LA, Gershwin LJ, Fanucchi MV, Van Winkle LS, Gerriets JE, et al. Repeated episodes of ozone inhalation amplifies the effects of allergen sensitization and inhalation on airway immune and structural development in Rhesus monkeys. Toxicol Appl Pharmacol. 2003;191(1):74–85.

    Article  CAS  PubMed  Google Scholar 

  23. Garcia E, Rice MB, Gold DR. Air pollution and lung function in children. J Allergy Clin Immunol. 2021; 148(1).https://doi.org/10.1016/j.jaci.2021.05.006.

  24. Grigg J. Air Pollution and Suppression of Lung Function Growth: A Triumph for Epidemiology. Am J Respir Crit Care Med. 2020;201(4):400–1. https://doi.org/10.1164/rccm.201911-2219ED.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Yang T, Chen R, Gu X, Xu J, Yang L, Zhao J, et al. Association of fine particulate matter air pollution and its constituents with lung function: The China Pulmonary Health study. Environ Int. 2021;156:106707. https://doi.org/10.1016/j.envint.2021.106707.

    Article  CAS  PubMed  Google Scholar 

  26. Havet A, Hulo S, Cuny D, Riant M, Occelli F, Cherot-Kornobis N, et al. Residential exposure to outdoor air pollution and adult lung function, with focus on small airway obstruction. Environ Res. 2020;183:109161. https://doi.org/10.1016/j.envres.2020.109161.

    Article  CAS  PubMed  Google Scholar 

  27. Liu S, Zhou Y, Liu S, Chen X, Zou W, Zhao D, et al. Association between exposure to ambient particulate matter and chronic obstructive pulmonary disease: results from a cross-sectional study in China. Thorax. 2017;72(9):788–95. https://doi.org/10.1136/thoraxjnl-2016-208910.

    Article  PubMed  Google Scholar 

  28. Jorgenson AK, Thombs RP, Clark B, Givens JE, Hill TD, Huang X, et al. Inequality amplifies the negative association between life expectancy and air pollution: A cross-national longitudinal study. Sci Total Environ. 2021;758:143705. https://doi.org/10.1016/j.scitotenv.2020.143705.

    Article  CAS  PubMed  Google Scholar 

  29. Bennett JE, Tamura-Wicks H, Parks RM, Burnett RT, Pope CA, Bechle MJ, et al. Particulate matter air pollution and national and county life expectancy loss in the USA: A spatiotemporal analysis. PLoS Med. 2019;16(7):e1002856. https://doi.org/10.1371/journal.pmed.1002856.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Kariisa M, Foraker R, Pennell M, Buckley T, Diaz P, Criner GJ, et al. Short- and long-term effects of ambient ozone and fine particulate matter on the respiratory health of chronic obstructive pulmonary disease subjects. Arch Environ Occup Health. 2015;70(1):56–62. https://doi.org/10.1080/19338244.2014.932753.

    Article  CAS  PubMed  Google Scholar 

  31. Priftis KN, Anthracopoulos MB, Paliatsos AG, Tzavelas G, Nikolaou-Papanagiotou A, Douridas P, et al. Different effects of urban and rural environments in the respiratory status of Greek schoolchildren. Respir Med. 2007;101(1):98–106.

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by Jiangsu social development project (Grant number and BK20221196), Independent research project (Grant number 2023YQFH08), and National Natural Science Foundation of China (Grant numbers U1602223 and 81773499).

Author information

Authors and Affiliations

Authors

Contributions

Y.X., Y.L., L.A., C.Z. and W.T. Conceptualize and design the study, Y.X., Y.L., H.L., Y.W., D.N., N.Y., D.H and J.Q. analyze and interpret the data, C.L. and G.W. draft the manuscript and W.T. and J.W. collect the data. All authors reviewed the manuscript.

Corresponding authors

Correspondence to Jiahong Wu or Weilong Tan.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

Fig. S1. AIC+BIC parameters for selecting the main model.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Xu, Y., Luo, Y., Yue, N. et al. Impact of outdoor air pollution on the incidence of pertussis in China: a time-series study. BMC Public Health 23, 2231 (2023). https://doi.org/10.1186/s12889-023-16530-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12889-023-16530-w

Keywords