Skip to main content
Download PDF
Download ePub

Association between short-term exposure to air pollution and ischemic stroke onset: a time-stratified case-crossover analysis using a distributed lag nonlinear model in Shenzhen, China

Environmental Health volume 19, Article number: 1 (2020) Cite this article

Abstract

Background

Stroke, especially ischemic stroke (IS), has been a severe public health problem around the world. However, the association between air pollution and ischemic stroke remains ambiguous.

Methods

A total of 63, 997 IS cases aged 18 years or above in Shenzhen were collected from 2008 to 2014. We used the time-stratified case-crossover design combining with distributed lag nonlinear model (DLNM) to estimate the association between air pollution and IS onset. Furthermore, this study explored the variability across gender and age groups.

Results

The cumulative exposure-response curves were J-shaped for SO2, NO2 and PM10, and V-shaped for O3, and crossed over the relative risk (RR) of one. The 99th, 50th (median) and 1st percentiles of concentration (μg/m3) respectively were 37.86, 10.06, 3.71 for SO2, 116.26, 41.29, 18.51 for NO2, 145.94, 48.29, 16.14 for PM10, and 111.57, 49.82, 16.00 for O3. Extreme high-SO2, high-NO2, high-PM10, high-O3, and low-O3 concentration increased the risk of IS, with the maximum RR values and 95% CIs: 1.50(1.22, 1.84) (99th vs median) at 0–12 lag days, 1.37(1.13, 1.67) (99th vs median) at 0–10 lag days, 1.26(1.04, 1.53) (99th vs median) at 0–12 lag days, 1.25(1.04, 1.49) (99th vs median) at 0–14 lag days, and 1.29(1.03, 1.61) (1st vs median) at 0–14 lag days, respectively. The statistically significant minimal RR value and 95% CI was 0.79(0.66,0.94) at 0–10 lag days for extreme low-PM10. The elderly aged over 65 years were susceptible to extreme pollution conditions. Difference from the vulnerability of males to extreme high-SO2, high-NO2 and low-O3, females were vulnerable to extreme high-PM10 and high-O3. Comparing with the elderly, adults aged 18–64 year were immune to extreme low-NO2 and low-PM10. However, no association between CO and IS onset was found.

Conclusions

SO2, NO2, PM10 and O3 exerted non-linear and delayed influence on IS, and such influence varied with gender and age. These findings may have significant public health implications for the prevention of IS.

Peer Review reports

Introduction

Stroke has been the second major cause of death and DALYs globally [1, 2]. There were about 10.3 million new strokes and 6.5 million deaths from stroke worldwide in 2013 [3]. In recent years, stroke has been the leading cause of death and DALYs in China [1]. A nationwide population-based survey of 480,687 Chinese adults reported that the age-standardized prevalence, incidence and mortality rates of stroke were 1114.8/100000 people, 246.8/100000 and 114.7/100000 person-years, respectively [4]. IS, which comprises 80% of stroke cases, has had a significant impact on healthcare expenditures and the Chinese economy [5]. Therefore, identification of modifiable risk factors for ischemic stroke would better informed public health prevention work. Exposure to air pollution is the largest environmental health risk and a growing global health problem estimated to contribute to as many as 3.1 million all-cause deaths per year [6,7,8]. Nevertheless, the impact of air pollution on morbidity from IS might be important and is less certain. Some studies have investigated the associations between short-time exposure to air pollution and stroke onset but no consensus have been reached [9, 10]. The heterogeneity in results may be partially explained by the difference across studies in regions and populations studied, levels, and constituents of air pollution, and sample size or others [11,12,13]. Besides, most of the studies on the association between air pollution and stroke in China only focused on metropolitan area or heavily polluted regions [14,15,16,17]. whereas few studies have been conducted in cities with relatively low levels of air pollution.

Shenzhen, a southern city of the Pearl River Delta region with a sub-tropical maritime climate, is less polluted than most cities in China. While, it is one of the three largest financial centers in China and one of the first-tier cities all over the world. Furthermore, Shenzhen has ranked first in air quality among China’s top 20 cities in GDP in recent years, achieving the development of both air quality and economy. Therefore, we conducted a study on the association between short-time exposure to air pollution and stroke in Shenzhen, an economically advanced and less polluted city of China.

Methods

Study population

The data on IS was obtained from 62 major municipal and district general hospitals in Shenzhen, including the National Ministry of Health’s Stroke Prevention and Treatment Base Hospital. The information included basic demographics, date of diagnosis, type of report, and the corresponding International Classification of Diseases, 10th Revision (ICD-10) codes and so on. The daily IS stroke count between 1 January 2008 and 31 December 2014 was identified using ICD-10 codes I63. In order to minimize the influence of coding inaccuracy, we used the corresponding diagnosis to check the identified cases. We then reserved the new cases, and individuals aged below 18 years were excluded from this study because of the low incidence.

Air pollution and meteorological data

Data on air pollution, including levels of sulfur dioxide (SO2), nitrogen dioxide (NO2), particulate matter less than 10 μm in aerodynamic diameter (PM10), carbon monoxide (CO), and ozone (O3), were obtained from the Shenzhen Environmental Monitoring Station between 1 January 2008 and 31 December 2014. There were 9 fixed-site air monitoring stations in Shenzhen, located in different districts. Monitoring of air pollution was done in accordance with mandatory quality assurance/quality control (QA/QC) procedures set by the State Environmental Protection Administration of China, ensuring the quality of automatic environmental air monitoring data. We averaged the measurements from all valid monitoring sites. The daily (24 h) mean concentrations of air pollutants were averaged from the available monitoring data across various stations. Data on meteorological factors, including temperature (°C), and relative humidity (%), were obtained from Shenzhen Meteorological Service Center during the period of 1 January 2008 to 31 December 2014.

Statistical analysis

In this study, we applied the time-stratified case-crossover (ts-CCO) design, regarded as a self-matched case-control study, which compares the exposure in the case period when events occurred with exposures in nearby referent periods, to examine the differences in exposure which may contribute to the differences in the daily count of cases [18]. Therefore, a time-stratified case-crossover design was adopted to regulate potential confounders (e.g., age, gender, etc) using self-control and exclude long-term impact of air pollutants (e.g., secular trend, seasonality, etc.) by stratification of time. We used the calendar month as the time stratum, to control the effects of long-term trend, seasonality, and day of the week [17]. We conducted a quasi-Poisson regression, controlling over-dispersion problem, combined with distributed lag non-linear model (DLNM) to estimate the non-linear and delayed influence of air pollution on IS onset. DLNM is based on the definition of “cross-basis”, a bi-dimensional space of functions to reflect the non-linear exposure-responses and lag structure of the association [19, 20]. In this study, consequently, the combination of DLNM with the ts-CCO design was employed, which allows estimating the short-term, non-linear and delayed effect of air pollutant using cross-basis functions for depicting the relationship between air pollutant and IS onset along the dimensions of exposure and lag simultaneously based on removing control confounders and long-term trend by ts-CCO design.

A quasi-Poisson regression model combined with time-stratified case-crossover design and DLNM was built as follows:

$$ {Y}_t\sim Poisson\left({\mu}_t\right) $$
$$ Log\left({\mu}_t/{Population}_{ye}\right)=\alpha +\sum \limits_{i=1}^m cb\left({Pollutant}_{i,t},{df}_{2i-1},{maxlag}_i,{df}_{2i}\right)+ cb\left({Temp}_t,{df}_{2m+1},{maxlag}_{m+1},{df}_{2m+2}\right)+ cb\left({RH}_t,{df}_{2m+3},{maxlag}_{m+2},{df}_{2m+4}\right)+\gamma {Holiday}_t+\lambda\ Stratum $$

where t is the day of observation; Yt is the count of IS cases on t; μt is the expectation of Yt; Populationye is the year-end population size; α is an intercept; Pollutanti, t, Tempt, and RHt are the ith pollutant concentration, temperature and relative humidity on t, respectively; cb() represents the cross-basis function with three pre-specified parameters of maximal lag maxlagi, degree of freedom for lag-response natural spline df2i − 1, and degree of freedom for exposure-response natural spline df2i for pollutant, temperature or relative humidity; Holiday is used to control the effect of public holidays; Stratum is the time stratum in the time-stratified case-crossover design. We defined natural cubic spline function with 3 df for air pollution and meteorological factors to mimic the exposure-response pattern of air pollution-IS onset associations, as well as lag spaces with 3 df to estimate the lag effects. To capture the complete lag-response curve, the maximal lag of air pollutants was set to 14 days; for the sake of simplification and without loss of generalization, meanwhile, this maximal lag was assigned to the length of the case and control periods. In addition, a 3-day duration was specified to be the maximal lag of meteorological factors. The df and maximum lag days for air pollution determination referred to the Akaike information criterion for quasi-Poisson (Q-AIC), which could produce the relatively superior model.

We initially conducted single-pollutant model to evaluate the association between air pollution and IS onset, and then the significant air pollutants were included in multi-pollutant model. Spearman’s correlation tests were used to estimate the associations between air pollution and meteorological factors, and pollutants with correlation coefficient r > 0.60 were not included in multi-pollutant model simultaneously to address the collinearity between air pollutants. In order to identify the high-risk or low-risk air pollution condition, the influence of extreme air pollution was evaluated and presented as relative risk (RR) by comparing the 99th above or 1st below percentiles of air pollution to the median values. We calculated the single day lag influence and the cumulative lag influence (lag0–1, lag0–6, lag0–8, lag0–10, lag0–12, lag0–13, and lag0–14) to effectively depict the characteristics of the association between air pollution and IS onset. In addition, we conducted stratified analysis to investigate the impact of air pollution on subgroups according to gender (male and female) and age groups (adult: 18–64 years; the elderly: ≥ 65 years).

All analyses were conducted using R version 3.5.1 with the dlnm package for fitting the DLNM, the gnm package for conditional quasi-Poisson regression.

Sensitivity analyses were performed to test the robustness of the selected model, which were as following: df [2,3,4,5,6] for air pollution, and df [2,3,4,5,6] for lag space were changed; the maximum lag days (12–21 days) for air pollution were extended.

Results

The characteristics of study population are presented in Table 1. There were 63,997 IS cases averagely aged 63.41 years met the inclusion criteria for the study between 2008 and 2014, of which 59.49% were males and 48.44% were the elderly.

Table 1 Characteristics of the study population in Shenzhen from 2008 to 2014
Full size table

The summary statistics for daily IS cases, air pollution and meteorological factors in Shenzhen are shown in Table 2. On average, 25 (range: 4–95) IS cases were identified each day during the study period. Of these, there were 15 (range: 1–64) male cases and 10 (range: 1–42) cases, 12 (range: 1–54) elderly cases and 13 (range: 0–42) adult cases respectively. The daily average air pollution levels were 11.97 μg/m3 (range: 2.95–70.63 μg/m3) for SO2, 45.89 μg/m3 (range: 13.13–166.14 μg/m3) for NO2, 55.91 μg/m3 (range: 10.86–182.23 μg/m3) for PM10, 1.22 μg/m3 (range: 0.37–3.25 μg/m3) for CO, and 53.49 μg/m3 (range: 6.02–31.04 μg/m3) for O3. The means of meteorological factors were 23.03 °C (range: 5.40–32.00 °C) for temperature and 72.83% (range: 19.00–100.00%) for relative humidity. The time-series analysis of IS cases and air pollution and meteorological factors are shown in Fig. 1.

Table 2 Descriptions of the day ischemic stroke, air pollution and meteorological factors in Shenzhen from 2008 to 2014
Full size table
Fig. 1
figure 1

Time-series results regarding the association of ischemic stroke onset with air pollution indicators and meteorological factors in Shenzhen from 2008 to 2014

Full size image

The correlations between air pollution and meteorological factors are presented in Table 3. The daily concentrations of SO2, NO2, and PM10 were highly and positively correlated with each other (correlation coefficient r > 0.6, P < 0.0001).

Table 3 Spearman correlation coefficients among air pollution and meteorological factors in Shenzhen from 2008 to 2014
Full size table

The results of the initial single-pollutant models indicated that SO2, NO2, PM10 and O3 were associated with IS onset. However, there was no statistical association between CO and IS onset. The single-day and cumulative relative risks of each pollutant in the single-pollutant models are shown in Additional file 1: Figure S1 and Figure S2.

Considering their greater Q-AIC values, the three-pollutant models for the full set combinations of SO2, NO2, PM10 and O3 were cleaned out from this study. In view of Q-AIC values of two-pollutant models, and the correlations among SO2, NO2, and PM10, they were included in two-pollutant model with O3, separately. Different air pollution variables exerted varied extremely influence on IS cases, and this difference were also observed among subgroups of the population.

The cumulative exposure-response curves for total IS cases at lag0–14 are shown in Fig. 2. The cumulative exposure-response curves were J-shaped for SO2, NO2 and PM10. It can be seen intuitively that only when concentration of SO2, NO2 and PM10 reached a certain limit the influence occurred. In addition, the cumulative lag influence of SO2, NO2 and PM10 were enhanced with the increase of concentration. The cumulative exposure-response curve of O3 was V-shaped, which meant that both high and low concentration of O3 may increase the risk of IS onset.

Fig. 2
figure 2

Summary of cumulative exposure-response curves on ischemic stroke for air pollution factors (SO2, NO2, PM10 and O3) for total cases at lag0–14 using two-pollutant model in Shenzhen, 2008–2014

Full size image

The single day lag-response curves for total IS cases are shown in Fig. 3. The curve of extreme high-SO2 was inverted V-shaped under different lag days, and RR value peaked on the eighth lag day and then decreased. No significant association between extreme low-SO2 and IS onset was observed at different lag days. The curve of the extreme high-NO2 was almost straight line, and the largest RR value was observed on the current day and lasted for 6 days. The influence of extreme high-PM10 also presented an inverted V-shape curve through the lag days, and the RR value reached maximum at lag7 and then decreased. The extreme low-PM10 showed a certain protective influence at lag0 to lag6. The extreme high-O3 presented statistically significant hazardous influence from lag4 to lag13.

Fig. 3
figure 3

Summary of single day lag-response curves on ischemic stroke for air pollution factors (SO2, NO2, PM10 and O3) for total cases at different lags using two-pollutant model in Shenzhen, 2008–2014. The extreme-high influence was estimated by the RR of ischemic stroke by comparing the 99th percentile of daily air pollution value to the median value, whereas the extreme-low influence was estimated by comparing the 1st percentile of daily air pollution value to the median value

Full size image

Table 4 summarizes the cumulative influence of extreme air pollution on IS at different lag days for total IS cases. The result showed that extremely high-SO2 could significantly increase the risk of IS, and the cumulative influence almost increased with lag days extending, and the maximum RR value was 1.50 (1.22, 1.84) (99th vs median), appearing at lag0–12. No prominently cumulative influence of extreme low-SO2 on the IS onset was observed at different lag days. Extreme high-NO2 was significant associated with IS, showing a hazardous influence, and the maximum RR value was 1.37 (1.13, 1.67) (99th vs median) at lag0–10. The influence of extreme low -NO2 was not observed in total group. Extreme high-PM10 can increase the risk of IS onset, and the maximum RR value was 1.26(1.04, 1.53) (99th vs median) at lag0–12. Extreme low-PM10 showed a protective influence on IS onset, and the minimal RR value was 0.79(0.66,0.94) at lag0–10. Both extreme high-O3 and low-O3 could increase the risk of IS, regardless of whether it was statistically significant or not, and the maximum RR values were 1.25(1.04, 1.49) (99th vs median) at lag0–14 and 1.29(1.03, 1.61) (1st vs median) at lag0–14, respectively.

Table 4 Extreme influence analysis of different air pollution factors from 2008 to 2014 using two-pollutant models. Relative risk (RR) and 95% confidence interval (CI) were used to estimate the cumulative influence of air pollution factors in total IS cases
Full size table

The cumulative exposure-response curves for subgroups shown in Fig. 4 were similar with total IS cases. It indicated that SO2 showed stronger influence on males and adults. The influence of NO2 and PM10 was a little stronger in males and the elderly than in females and adults. Yet, the impact of O3 shows almost no difference among males and adults, females and the elderly respectively.

Fig. 4
figure 4

Summary of cumulative exposure-response curves on ischemic stroke for air pollution factors (SO2, NO2, PM10 and O3) for subgroups at lag0–14 using two-pollutant model in Shenzhen, 2008–2014. Male and female were subgroups according to gender. The elderly and adult were subgroups according to age (adult: 18–64 years; the elderly: ≥ 65 years)

Full size image

The single day lag-response curves of subgroups presented in Additional file 1: Figure S3 resembled these of total IS cases.

The cumulative influence of extreme air pollution on IS at different lag days for subgroups are shown in Table 5. Males were more susceptible to extreme high-SO2 than females, and the maximum RR value is 40% higher. Extreme high-SO2 had a stronger influence on the elderly and adults. However, no prominently cumulative influence of extreme low-SO2 on subgroups was observed at different lag days. The influence of extreme high-NO2 on males was slightly stronger than females, and the RR value was 0.19 higher. Extreme high-NO2 showed a strong influence on the elderly, but insignificantly on adults. Extreme low-NO2 presented protective influence on males and adults. Females and the elderly were more sensitive for extreme high-PM10. Extreme low-PM10 showed protective influence, especially on females and adults. Extreme high-O3 showed significant risk influence in females and the elderly. Besides, extreme low-O3 showed significant risk influence on males and adults.

Table 5 Extreme influence analysis of different air pollution factors from 2008 to 2014 using two-pollutant models. Relative risk (RR) and 95% confidence interval (CI) were used to estimate the cumulative influence of air pollution factors in subgroups
Full size table

To facilitate comparison with previous studies, the RRs and 95% CIs of air pollution factors on IS onset at lag0, lag1 and lag0–1 are listed in Additional file 1: Table S1. In addition, the cumulative exposure-response curves of air pollutants on IS at lag0–1, lag0–6 and lag0–14 are presented in Additional file 1: Figure S4.

The results of sensitivity analysis, which were conducted by changing the df for air pollution from 2 to 6, the df for the lag space from 2 to 6, and by changing the maximum lag days from 12 to 21 days, were similar to the results obtained above. Results of sensitivity analysis were shown in Additional file 1: Figure S5-S10.

Discussion

In this study, we found significant association between short-term exposure to SO2, NO2, PM10, and O3 and IS onset. The influence of exposing to high-level and low-level air pollution was not concerted for subgroups depending on gender and age.

The present study found that extreme high-SO2, high-NO2, and high-PM10 could increase the risk of IS. A multicity case-crossover study found significantly positive associations between SO2 [RR and 95% confidence interval (CI): 1.016(01.0100, 1.0230) per interquartile range (IQR) increase], NO2 [RR and 95% CI: 1.0260(1.0180, 1.0350) per IQR increase], PM10 [RR and 95% CI: 1.0070(1.000, 1.0140) per IQR increase] and IS [21]. Another multicity study reported that SO2 [RR and 95% CI: 1.0135 (1.0043, 1.0229) per IQR increase], NO2 [RR and 95% CI: 1.0294 (1.0178, 1.0412) per IQR increase] and PM10 [RR and 95% CI: 1.0103 (1.0004, 1.0204) per IQR increase] were significantly associated with IS onset [22]. In addition, meta-analyses also found significantly positive associations between SO2, NO2, PM10 and stroke [9, 10]. Previous studies used different study designs and model specifications, which made us unable to compare the non-linear and lag influence across cities. However, these findings of associations between high concentration and IS onset were similar as a whole. Thus, the associations between the three kinds of air pollutants and IS were unlikely to be mendacious because of potential confounding, defects of the design or statistical analysis, or publication and reporting bias. Conflicting evidences were found on the modification by gender or age in the associations between air pollution and IS [16, 21, 23, 24], and the underlying mechanisms were ambiguous.

Biologic mechanisms for these associations above have not been fully established, and majority of previous studies focused on particulate matter (PM). Several potential mechanisms have been proposed including systemic inflammation [25, 26], thrombosis [27,28,29], artery calcification [30], and vascular endothelial dysfunction [31], induced by exposure air pollution. The acute systemic inflammatory response with increased plasma fibrinogen, C-reactive protein and white blood cell [25, 32, 33], could be a trigger for inflammation and thrombokinesis. These air-pollution-related pathophysiologic changes may be associated with the occurrence of IS.

We found no association between CO and IS onset, although it was controversial in earlier studies. A study of 9 cities suggested a positive association between CO (RR and 95% CI: 1.0283 (1.0123, 1.0446) per 0.3 ppm increase) and IS onset in single-pollutant model [22], while the association may not persist when adjusted by other pollutants. Another study conducted in Copenhagen, Denmark found ambient CO was associated with an increased risk of IS in single-pollutant model and the association was attenuated to null after adjusting for PM [34]. In Taiwan, significant positive association was found between CO and IS admission in the single-pollutant model, but it became insignificant when controlled for other pollutants [35]. The lack of multi-pollutant models may distort the association. A time series study in Hong Kong found a protective influence between CO [RR and 95% CI was 0.980 (0.967, 0.993) per 0.3 ppm increase] and stroke admission [36]. However, the study did not differentiate between ischemic and hemorrhagic stroke, so it may be not sure that whether the association exists when we limited the outcome to IS. Also, Also, there are several other multicity-based studies in China found no association between CO and IS [21, 37], which supported our results.

The cumulative dose-response relationship between O3 and IS onset was V-shaped, which meant that higher or lower levels of O3 concentrations may both increase the risk of IS. However, evidences from previous studies were incongruous. Several studies reported significantly positive association between O3 and IS onset [38,39,40], or stroke emergency hospital visits [41]. Nevertheless, some other studies found no association between O3 and stroke incidence [11, 41,42,43]. One possible explanation for the inconsistent findings was that O3 may increase the risk of IS onset among susceptible or vulnerable populations [44]. Several other studies found insignificant negative association between O3 and stroke onset, with estimation of OR with 95% CI as 0.97(0.87, 1.07) per 10 ppb increase in Nueces County, Texas [42], and 0.98(0.96, 1.01) per 10 ppb increase in South Carolina [43], respectively. A recent study in Changzhou, China reported a protective influence of ambient O3 on stroke onset, and the RR with 95% CI was 0.9966(0.9944, 0.9988) per IQR increase [45]. Major ozonated auto-hemo-therapy has been used in the treatment in ischemic disorders, including acute cerebral infarction [46, 47]. In addition, the neuroprotective dose-response curve for O3 after a stroke was shaped as U with effective level range of 80–120 μg/mL in rat models [48]. Therefore, the consistency of this biological evidence with our epidemiological results indicates necessity of conducting more investigations to validate complicate association of O3 with IS and then to determine range of low-risk O3 concentration.

There are some strengths of the present study. First, the data on IS were derived from the majority of general hospitals with the ability to diagnose and treat stoke in Shenzhen, and the sample size was relatively large. Therefore, the present results may be representative of the authentic associations between air pollution and IS in the study area. Second, we only included the new IS cases, as exclusion of the recurrent cases can better reflect the influence of air pollution on IS onset. Third, as far as we know, it was the first time that a time-stratified case-crossover design combining with DLNM has been used to study the between air pollution and IS onset. The case-crossover design is a reliable method to estimate the short-term health influence of exposure to air pollution [49], and the DLNM can better depict the non-linear and delayed influence of air pollution on IS onset. In addition, we found the V-shaped relationship between O3 and IS for the first time, which may important implications for the prevention and treatment of IS. However, several limitations should be addressed. Firstly, the present results were derived from data of only one city. Although Shenzhen is a representative of cities with advanced economy and low levels of air pollution, it should be cautious to extrapolate our results to other areas. Secondly, the date of diagnosis, rather than the time of stroke symptom onset, was used in the analysis, which may result in temporal misalignment between air pollution exposure and IS incidence and the underestimate of exposure effects [50]. However, because stoke is a kind of disease requiring prompt hospitalization and treatment, temporal misalignment is expected to play a minor role in our estimation. As the admission or diagnosis may occur in the days after onset of symptoms, it should be gingerly to explain the lag influence of air pollution on IS onset. Finally, the data on air pollution was evaluated using the arithmetic mean from 9 fixed stations, but not individual exposure, which may cause exposure measurement error, leading to underestimation of the influence of air pollution [51].

Although the increased risk of IS triggered by air pollution is relatively small for each individual, the public health implications are very important because a large number of populations are at risk for IS and exposed to unavoided and modifiable air pollution. Further studies with more accurate time of stroke symptom onset are needed to further validate our findings. Additional studies should be conducted to investigate the specific components of air pollutants that play a role in the association between air pollution and IS onset.

Conclusions

In conclusion, our study suggested that the short-term exposure to SO2, NO2, and PM10 was significantly associated with increased IS risk, except for CO. The cumulative exposure-response curve between O3 and IS onset was V-shaped, which may provide an epidemiological evidence for the neuroprotective function. These findings may have significant public health implications for the prevention of IS. Further studies on this topic are warranted to validate our research.

Availability of data and materials

The datasets generated and analyzed during the current study are not publicly available due to no permission but are available from the corresponding author on reasonable request.

Abbreviations

CI:

Confidence interval

DLNM:

Distributed lag nonlinear model

ICD-10:

International Classification of Diseases, 10th Revision

IQR:

Interquartile range

IS:

Ischemic stroke

RR :

Relative risk

References

  1. GBD 2013 Mortality and Causes of Death Collaborators. Global, regional, and national age-sex specific all-cause and cause-specific mortality for 240 causes of death, 1990–2013: A systematic analysis for the global burden of disease study 2013. Lancet. 2015;385(9963):117–71.

    Article  Google Scholar 

  2. Murray CJ, Barber RM, Foreman KJ, et al. Global, regional, and national disability-adjusted life years (DALYs) for 306 diseases and injuries and healthy life expectancy (HALE) for 188 countries, 1990-2013: Quantifying the epidemiological transition. Lancet. 2015;386(10009):2145–91.

    Article  Google Scholar 

  3. Feigin VL, Krishnamurthi RV, Parmar P, et al. Update on the global burden of ischemic and hemorrhagic stroke in 1990-2013: the GBD 2013 study. Neuroepidemiology. 2015;45(3):161–76.

    Article  Google Scholar 

  4. Wang W, Jiang B, Sun H, et al. Prevalence, incidence, and mortality of stroke in China: results from a Nationwide population-based survey of 480 687 adults. Circulation. 2017;135(8):759–71.

    Article  Google Scholar 

  5. Liu L, Wang D, Wong KS, Wang Y. Stroke and stroke care in China: huge burden, significant workload, and a national priority. Stroke. 2011;42(12):3651–4.

    Article  Google Scholar 

  6. Lelieveld J, Evans JS, Fnais M, Giannadaki D, Pozzer A. The contribution of outdoor air pollution sources to premature mortality on a global scale. Nature. 2015;525(7569):367–71.

    Article  CAS  Google Scholar 

  7. Newby DE, Mannucci PM, Tell GS, et al. Expert position paper on air pollution and cardiovascular disease. Eur Heart J. 2015;36(2):83–93b.

    Article  CAS  Google Scholar 

  8. Brook RD, Newby DE, Rajagopalan S. The global threat of outdoor ambient air pollution to cardiovascular health: time for intervention. JAMA Cardiol. 2017;2(4):353–4.

    Article  Google Scholar 

  9. Shah AS, Lee KK, McAllister DA, et al. Short term exposure to air pollution and stroke: systematic review and meta-analysis. BMJ (Clinical research ed). 2015;350:h1295.

    Google Scholar 

  10. Yang WS, Wang X, Deng Q, Fan WY, Wang WY. An evidence-based appraisal of global association between air pollution and risk of stroke. Int J Cardiol. 2014;175(2):307–13.

    Article  Google Scholar 

  11. Wellenius GA, Burger MR, Coull BA, et al. Ambient air pollution and the risk of acute ischemic stroke. Arch Intern Med. 2012;172(3):229–34.

    Article  Google Scholar 

  12. Butland BK, Atkinson RW, Crichton S, et al. Air pollution and the incidence of ischaemic and haemorrhagic stroke in the South London stroke register: a case-cross-over analysis. J Epidemiol Community Health. 2017;71(7):707–12.

    Article  CAS  Google Scholar 

  13. Vivanco-Hidalgo RM, Wellenius GA, Basagana X, et al. Short-term exposure to traffic-related air pollution and ischemic stroke onset in Barcelona, Spain. Environ Res. 2018;162:160–5.

    Article  CAS  Google Scholar 

  14. Dong H, Yu Y, Yao S, et al. Acute effects of air pollution on ischaemic stroke onset and deaths: a time-series study in Changzhou. China. 2018;8(7):e020425.

    Google Scholar 

  15. Guo P, Wang Y, Feng W, et al. Ambient air pollution and risk for ischemic stroke: a short-term exposure assessment in South China. Int J Environ Res Public Health. 2017;14(9):1091.

    Article  CAS  Google Scholar 

  16. Huang F, Luo Y, Guo Y, et al. Particulate matter and hospital admissions for stroke in Beijing, China: modification effects by ambient temperature. J Am Heart Assoc. 2016;5(7):e003437.

    Google Scholar 

  17. Zanobetti A, Schwartz J. The effect of particulate air pollution on emergency admissions for myocardial infarction: a multicity case-crossover analysis. Environ Health Perspect. 2005;113(8):978–82.

    Article  Google Scholar 

  18. Janes H, Sheppard L, Lumley T. Case-crossover analyses of air pollution exposure data: referent selection strategies and their implications for bias. Epidemiology. 2005;16(6):717–26.

    Article  Google Scholar 

  19. Gasparrini A. Modeling exposure-lag-response associations with distributed lag non-linear models. Stat Med. 2014;33(5):881–99.

    Article  Google Scholar 

  20. Gasparrini A, Armstrong B, Kenward MG. Distributed lag non-linear models. Stat Med. 2010;29(21):2224–34.

    Article  CAS  Google Scholar 

  21. Liu H, Tian Y, Xu Y, et al. Association between ambient air pollution and hospitalization for ischemic and hemorrhagic stroke in China: a multicity case-crossover study. Environ Pollut. 2017;230:234–41.

    Article  CAS  Google Scholar 

  22. Wellenius GA, Schwartz J, Mittleman MA. Air pollution and hospital admissions for ischemic and hemorrhagic stroke among medicare beneficiaries. Stroke. 2005;36(12):2549–53.

    Article  CAS  Google Scholar 

  23. Zheng S, Wang M, Wang S, Tao Y, Shang K. Short-term effects of gaseous pollutants and particulate matter on daily hospital admissions for cardio-cerebrovascular disease in Lanzhou: evidence from a heavily polluted city in China. Int J Environ Res Public Health. 2013;10(2):462–77.

    Article  CAS  Google Scholar 

  24. Villeneuve PJ, Chen L, Stieb D, Rowe BH. Associations between outdoor air pollution and emergency department visits for stroke in Edmonton, Canada. Eur J Epidemiol. 2006;21(9):689–700.

    Article  Google Scholar 

  25. Rich DQ, Kipen HM, Huang W, et al. Association between changes in air pollution levels during the Beijing Olympics and biomarkers of inflammation and thrombosis in healthy young adults. Jama. 2012;307(19):2068–78.

    Article  CAS  Google Scholar 

  26. Mills NL, Tornqvist H, Robinson SD, et al. Diesel exhaust inhalation causes vascular dysfunction and impaired endogenous fibrinolysis. Circulation. 2005;112(25):3930–6.

    Article  CAS  Google Scholar 

  27. Bind MA, Baccarelli A, Zanobetti A, et al. Air pollution and markers of coagulation, inflammation, and endothelial function: associations and epigene-environment interactions in an elderly cohort. Epidemiology. 2012;23(2):332–40.

    Article  Google Scholar 

  28. Franchini M, Mannucci PM. Thrombogenicity and cardiovascular effects of ambient air pollution. Blood. 2011;118(9):2405–12.

    Article  CAS  Google Scholar 

  29. Lucking AJ, Lundback M, Mills NL, et al. Diesel exhaust inhalation increases thrombus formation in man. Eur Heart J. 2008;29(24):3043–51.

    Article  CAS  Google Scholar 

  30. Kaufman JD, Adar SD, Barr RG, et al. Association between air pollution and coronary artery calcification within six metropolitan areas in the USA (the Multi-Ethnic Study of Atherosclerosis and Air Pollution): a longitudinal cohort study. Lancet. 2016;388(10045):696–704.

    Article  CAS  Google Scholar 

  31. Tornqvist H, Mills NL, Gonzalez M, et al. Persistent endothelial dysfunction in humans after diesel exhaust inhalation. Am J Respir Crit Care Med. 2007;176(4):395–400.

    Article  CAS  Google Scholar 

  32. Pekkanen J, Brunner EJ, Anderson HR, Tiittanen P, Atkinson RW. Daily concentrations of air pollution and plasma fibrinogen in London. Occup Environ Med. 2000;57(12):818–22.

    Article  CAS  Google Scholar 

  33. Schwartz J. Air pollution and blood markers of cardiovascular risk. Environ Health Perspect. 2001;109(Suppl 3):405–9.

    Article  CAS  Google Scholar 

  34. Andersen ZJ, Olsen TS, Andersen KK, Loft S, Ketzel M, Raaschou-Nielsen O. Association between short-term exposure to ultrafine particles and hospital admissions for stroke in Copenhagen, Denmark. Eur Heart J. 2010;31(16):2034–40.

    Article  CAS  Google Scholar 

  35. Tsai SS, Goggins WB, Chiu HF, Yang CY. Evidence for an association between air pollution and daily stroke admissions in Kaohsiung, Taiwan. Stroke. 2003;34(11):2612–6.

    Article  CAS  Google Scholar 

  36. Tian L, Qiu H, Pun VC, Ho KF, Chan CS, Yu IT. Carbon monoxide and stroke: a time series study of ambient air pollution and emergency hospitalizations. Int J Cardiol. 2015;201:4–9.

    Article  Google Scholar 

  37. Tian Y, Liu H, Zhao Z, Xiang X. Association between ambient air pollution and daily hospital admissions for ischemic stroke: A nationwide time-series analysis. PLoS Med. 2018;15(10):e1002668.

    Article  CAS  Google Scholar 

  38. Xu X, Sun Y, Ha S, Talbott EO, Lissaker CT. Association between ozone exposure and onset of stroke in Allegheny County, Pennsylvania, USA, 1994-2000. Neuroepidemiology. 2013;41(1):2–6.

    Article  CAS  Google Scholar 

  39. Henrotin JB, Besancenot JP, Bejot Y, Giroud M. Short-term effects of ozone air pollution on ischaemic stroke occurrence: a case-crossover analysis from a 10-year population-based study in Dijon, France. Occup Environ Med. 2007;64(7):439–45.

    Article  CAS  Google Scholar 

  40. Wing JJ, Sanchez BN, Adar SD, et al. Synergism of short-term air pollution exposures and neighborhood disadvantage on initial stroke severity. Stroke. 2017;48(11):3126–9.

    Article  CAS  Google Scholar 

  41. Carlsen HK, Forsberg B, Meister K, Gislason T, Oudin A. Ozone is associated with cardiopulmonary and stroke emergency hospital visits in Reykjavik, Iceland 2003-2009. Environ Health. 2013;12(28):1–8.

    Google Scholar 

  42. Wing JJ, Adar SD, Sanchez BN, Morgenstern LB, Smith MA, Lisabeth LD. Short-term exposures to ambient air pollution and risk of recurrent ischemic stroke. Environ Res. 2017;152:304–7.

    Article  CAS  Google Scholar 

  43. Montresor-Lopez JA, Yanosky JD, Mittleman MA, et al. Short-term exposure to ambient ozone and stroke hospital admission: a case-crossover analysis. J Exp Sci Environ Epidemiol. 2016;26(2):162–6.

    Article  CAS  Google Scholar 

  44. Bell ML, Zanobetti A, Dominici F. Who is more affected by ozone pollution? A systematic review and meta-analysis. Am J Epidemiol. 2014;180(1):15–28.

    Article  Google Scholar 

  45. Yu Y, Dong H, Yao S, Ji M, Yao X, Zhang Z. Protective effects of ambient ozone on incidence and outcomes of ischemic stroke in Changzhou, China: a time-series study. Int J Environ Res Public Health. 2017;14(12):1610–20.

    Article  CAS  Google Scholar 

  46. Sancak EB, Turkon H, Cukur S, et al. Major Ozonated Autohemotherapy preconditioning ameliorates kidney ischemia-reperfusion injury. Inflammation. 2016;39(1):209–17.

    Article  CAS  Google Scholar 

  47. Wu X, Li Z, Liu X, et al. Major ozonated autohemotherapy promotes the recovery of upper limb motor function in patients with acute cerebral infarction. Neural Regen Res. 2013;8(5):461–8.

    CAS  Google Scholar 

  48. Frosini M, Contartese A, Zanardi I, Travagli V, Bocci V. Selective ozone concentrations may reduce the ischemic damage after a stroke. Free Radic Res. 2012;46(5):612–8.

    Article  CAS  Google Scholar 

  49. Jaakkola JJ. Case-crossover design in air pollution epidemiology. Eur Respir J Suppl. 2003;40:81s–5s.

    Article  CAS  Google Scholar 

  50. Lokken RP, Wellenius GA, Coull BA, et al. Air pollution and risk of stroke: underestimation of effect due to misclassification of time of event onset. Epidemiology. 2009;20(1):137–42.

    Article  Google Scholar 

  51. Goldman GT, Mulholland JA, Russell AG, et al. Impact of exposure measurement error in air pollution epidemiology: effect of error type in time-series studies. Environ Health. 2011;10:61–72.

    Article  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by the National Natural Science Foundation of China (No. 81573262) from National Institutes of Health.

Author information

Author notes

  1. Zhinghui Wang, Ji Peng, and Peiyi Liu contribute equally to this paper.

  2. Jinquan Cheng and Ping Yin contribute equally to this work and should be considered as co-corresponding authors.

Authors and Affiliations

  1. Department of Epidemiology and Biostatistics, School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, 13 Hangkong Rd, Wuhan, 430030, Hubei, China

    Zhinghui Wang, Yanran Duan, Hongyan Li, Siyu Yan & Ping Yin

  2. Shenzhen Center for Chronic Disease Control, Shenzhen, China

    Ji Peng

  3. Department of Molecular Epidemiology, Shenzhen Center for Disease Control and Prevention, Shenzhen, China

    Peiyi Liu, Suli Huang & Ying Wen

  4. Department of Occupational and Environment Health, Key Laboratory of Environment and Health, School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China

    Peiyi Liu

  5. Department of Public Health Promotion, Shenzhen Center for Disease Control and Prevention, Shenzhen, China

    Yi Liao

  6. Shenzhen Center for Disease Control and Prevention, 8 Longyuan Rd, Shenzhen, 518055, Guangdong, China

    Jinquan Cheng & Ping Yin

Authors
  1. Zhinghui Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Ji Peng
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Peiyi Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Yanran Duan
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Suli Huang
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Ying Wen
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Yi Liao
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Hongyan Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Siyu Yan
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Jinquan Cheng
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Ping Yin
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

ZW, JP, PL, JC, and PY contributed to the conception and design of the study. ZW led the data analysis with support from JP, PL. SH, and JC contributed to the ischemic stroke data collection and data quality control. YW and YL contributed to the exposure data collection. YD, HL, and SY contributed to literature review. SY and JP supervised the analysis and generating of results. ZW drafted and finalised the paper with input from JP, PL, JC, and PY. All authors contributed to data interpretation, and reviewed and approved the final version manuscript. PY and JC had full access to all the data in the study and had final responsibility for the decision to submit for publication.

Corresponding authors

Correspondence to Jinquan Cheng or Ping Yin.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Additional file

Additional file 1.

Supplemental materials.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Wang, Z., Peng, J., Liu, P. et al. Association between short-term exposure to air pollution and ischemic stroke onset: a time-stratified case-crossover analysis using a distributed lag nonlinear model in Shenzhen, China. Environ Health 19, 1 (2020). https://doi.org/10.1186/s12940-019-0557-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12940-019-0557-4

Keywords

  • Ischemic stroke
  • Air pollution
  • Case-crossover design
  • Distributed lag nonlinear model
  • Quasi-Poisson regression

Environmental Health

ISSN: 1476-069X