Figure 1. Odds ratio of ischemic stroke onset for US Environmental Protection Agency categories (good and moderate) of mean ambient fine particulate matter air pollution (PM2.5) levels in the 24 hours preceding stroke onset. Error bars indicate 95% CIs.
Figure 2. Odds ratio of ischemic stroke onset per interquartile range increase in concentration of ambient fine particulate matter air pollution (6.4 μg/m3) in the hours preceding stroke onset. Error bars indicate 95% CIs.
Wellenius GA, Burger MR, Coull BA, Schwartz J, Suh HH, Koutrakis P, Schlaug G, Gold DR, Mittleman MA. Ambient Air Pollution and the Risk of Acute Ischemic Stroke. Arch Intern Med. 2012;172(3):229–234. doi:10.1001/archinternmed.2011.732
Background The link between daily changes in level of ambient fine particulate matter (PM) air pollution (PM <2.5 μm in diameter [PM2.5]) and cardiovascular morbidity and mortality is well established. Whether PM2.5 levels below current US National Ambient Air Quality Standards also increase the risk of ischemic stroke remains uncertain.
Methods We reviewed the medical records of 1705 Boston area patients hospitalized with neurologist-confirmed ischemic stroke and abstracted data on the time of symptom onset and clinical characteristics. The PM2.5 concentrations were measured at a central monitoring station. We used the time-stratified case-crossover study design to assess the association between the risk of ischemic stroke onset and PM2.5 levels in the hours and days preceding each event. We examined whether the association with PM2.5 levels differed by presumed ischemic stroke pathophysiologic mechanism and patient characteristics.
Results The estimated odds ratio (OR) of ischemic stroke onset was 1.34 (95% CI, 1.13-1.58) (P < .001) following a 24-hour period classified as moderate (PM2.5 15-40 μg/m3) by the US Environmental Protection Agency's (EPA) Air Quality Index compared with a 24-hour period classified as good (≤15 μg/m3). Considering PM2.5 levels as a continuous variable, we found the estimated odds ratio of ischemic stroke onset to be 1.11 (95% CI, 1.03-1.20) (P = .006) per interquartile range increase in PM2.5 levels (6.4 μg/m3). The increase in risk was greatest within 12 to 14 hours of exposure to PM2.5 and was most strongly associated with markers of traffic-related pollution.
Conclusion These results suggest that exposure to PM2.5 levels considered generally safe by the US EPA increase the risk of ischemic stroke onset within hours of exposure.
Daily changes in levels of ambient fine particulate matter (PM) air pollution (PM <2.5 μm in diameter [PM2.5]) have been associated with higher risk of acute cardiovascular events, excess hospitalizations, and deaths.1 These cardiovascular effects of PM2.5 appear to be mediated through a combination of autonomic, hemostatic, inflammatory, and vascular endothelial disturbances with consequent changes in cardiac and vascular function.2- 7 Relying on the current evidence, the US Environmental Protection Agency (EPA) regulates mean daily and annual PM2.5 levels. Whether the current regulatory standards are sufficient to protect public health remains controversial.8
Although a number of studies have examined the association between other air pollutants and the risk of stroke,9- 12 fewer studies have evaluated the effects of PM2.5 on stroke risk, and the results remain equivocal. Specifically, some,13- 15 but not all,16- 18 studies suggest that PM2.5 exposure may increase the incidence of the combined end point of acute cerebrovascular diseases, an etiologically diverse group with multiple underlying pathophysiologic mechanisms. Prior work has suggested that ambient air pollution is more strongly associated with ischemic stroke than intracranial hemorrhage.19 However, few studies have specifically evaluated the link between PM2.5 exposure and ischemic stroke risk,14,20- 22 and only 1 of these studies assessed whether the associations differ by presumed ischemic stroke cause.22
We therefore evaluated the association between daily and hourly changes in PM2.5 levels and the risk of ischemic stroke onset among patients residing in the greater Boston area and admitted between 1999 and 2008 to the Beth Israel Deaconess Medical Center (BIDMC). During the study period, PM2.5 levels in the Boston area did not exceed current EPA standards.
This study was approved by the committee on clinical investigations at BIDMC. We identified 1763 consecutive patients 21 years or older admitted to the BIDMC between April 1, 1999, and October 31, 2008, with neurologist-confirmed ischemic stroke, excluding patients with in-hospital strokes or transient ischemic attacks. The BIDMC is a 650-bed teaching hospital of Harvard Medical School and designated as a primary stroke service hospital by the state. The stroke service consists of 5 vascular neurologists who see about 550 patients with stroke annually. As in previous studies,23 we excluded patients residing farther than 40 km from the Harvard ambient monitoring station to reduce exposure misclassification. We identified patients potentially eligible for this study by reviewing daily emergency department admission logs, stroke service admission logs, stroke service consult logs, and hospital electronic discharge records. For patients meeting eligibility criteria, we abstracted data on demographics, presenting symptoms, medical history (including history of stroke), and imaging results from each patient's medical record. We classified presumed stroke pathophysiologic mechanisms as (1) large-artery atherosclerosis, (2) small-vessel occlusion, (3) cardioembolism, (4) other determined cause or (5) undetermined cause, using the approach developed for the Trial of ORG 10172 in Acute Stroke Treatment (TOAST).24 Time of stroke symptom onset or time last seen normal as documented by the attending stroke neurologist at the time of hospital presentation was classified as exact, estimated, or unknown.22 Given the documented morning peak in ischemic stroke incidence,25 we assumed that stroke onset occurred at 9:00 AM for 221 patients for whom the date of stroke symptom onset was documented but not the time (13%). We excluded 58 patients for whom neither the date nor the time of stroke onset was documented (3%), leaving 1705 patients available for analysis.
The concentrations of PM2.5 and black carbon were measured continuously, and levels of sulfate particles (SO42-) were measured daily (9:00 AM-9:00 AM) at the Harvard ambient monitoring station, as previously described.23 Hourly measures of nitrogen dioxide (NO2), carbon monoxide (CO), and ozone (O3) were obtained from local monitoring sites operated by the Massachusetts Department of Environmental Protection and averaged. We obtained hourly meteorologic data from the National Weather Service station at Boston's Logan Airport and calculated apparent temperature, an index of thermal comfort, as previously described.26 In a secondary analysis, we estimated exposure to black carbon at each subject's home as a marker of traffic pollution using a validated temporal-spatial model.27 Black carbon predictions are based on over 6000 black carbon measurements at 82 locations in the greater Boston area, meteorologic and other characteristics of a given day, and measures of land use (eg, traffic and population density) at a given location.
We used the time-stratified case-crossover study design28 to assess the association between the risk of ischemic stroke onset and PM2.5 concentrations in the hours and days preceding each event. In this design, each subject's exposure prior to a case-defining event (case period) is compared with his or her own exposure experience during 1 or more control periods when the subject did not become a case (control period). Control periods were chosen such that exposures during the case period were compared with exposures occurring on other days of the same month falling on the same day of the week and time of day as the case period. The use of control periods from both before and after the index event is appropriate in this setting because individual events do not affect the distribution of future exposure in the overall study population.29 This approach effectively controls for seasonality, time trends, and chronic and slowly varying potential confounders because the exposure information for cases and controls within the same stratum come from the same calendar month.30
We performed conditional logistic regression, stratifying on each hospitalization, to obtain estimates of odds ratios (ORs) associated with PM2.5 exposure and corresponding 95% CIs. In all analyses, we controlled for ambient temperature and dew point temperature using natural cubic splines (3 degrees of freedom each) and barometric pressure modeled as a linear continuous variable. We first considered 2 categories of PM2.5 levels defined a priori by the EPA's Air Quality Index31 as either good (≤15 μg/m3) or moderate (15-40 μg/m3), excluding 11 days where PM2.5 levels exceeded 40 μg/m3 (0.3%). Next we considered 5 a priori categories of PM2.5 levels (break points at 5, 10, 15, and 20 μg/m3). Finally, we considered PM2.5 exposure as a continuous variable.
In all analyses, PM2.5 exposure was assessed relative to the time of stroke symptom onset. We separately evaluated the association between the risk of ischemic stroke onset and PM2.5 levels averaged at 4 different periods prior to stroke symptom onset (0 to <24 hours, 24 to <48 hours, 48 to <72 hours, and 72 to <96 hours). Results were subsequently confirmed using unconstrained distributed lag models such that pollutant levels at each time point were considered jointly in a single model. To explore associations with shorter exposures in more detail, we additionally evaluated the association between risk of stroke onset and PM2.5 levels averaged in 2-hour increments prior to stroke symptom onset (0 to <2 hours through 36 to <38 hours). We repeated these analyses for black carbon, NO2, CO, O3, and SO42-.
We evaluated whether the associations with PM2.5 concentrations differed by presumed ischemic stroke cause or according to the presence of major stroke risk factors including diabetes mellitus, atrial fibrillation, hypertension, and history of stroke or transient ischemic attack using fully stratified models. We used the χ2 test for homogeneity32 to evaluate whether associations differed significantly across subgroups. A 2-sided P < .05 was considered statistically significant. Analyses were performed using SAS version 9.2 (SAS Institute Inc) and the R statistical package, version 2.8.1.
The 1705 patients were predominantly white women, mean (SD) age 73.1 (14.5) years (Table 1). Small-vessel strokes were the most common determined cause of stroke (26%), followed by strokes due to cardioembolism (25%) and large-artery atherosclerosis (20%). The median delay time from stroke symptom onset to hospital presentation was 10 hours (25th percentile, 3 hours; 75th percentile, 26 hours). In-hospital mortality was 5.8%. Most patients resided less than 20 km from the Harvard ambient monitoring site in Boston (97%) (eTable 1).
For 2888 days during the study period (83% of study days), the air quality was classified as good, according to the EPA's Air Quality Index for PM2.5 (“air quality is satisfactory and poses little or no risk”33(p3)); for 572 days (16% of study days), it was classified as moderate (“air quality is acceptable; however, there may be a moderate health concern for a very small number of people”33(p3)).
The OR of ischemic stroke onset was 1.34 (95% CI, 1.13-1.58) (P < .001) following a 24-hour period during which the air quality was moderate compared with a period during which it was good. The association between PM2.5 level and risk of ischemic stroke onset was approximately linear (Figure 1). Considering PM2.5 concentration as a continuous variable, the OR of stroke onset was 1.11 (95% CI, 1.03-1.20) (P = .006) comparing the 75th to the 25th percentile of PM2.5 levels (6.4 μg/m3) over the previous 24 hours. The PM2.5 levels preceding stroke onset by more than 24 hours were not associated with higher risk.
The results were not materially different when we excluded from analysis patients living more than 20 km from the monitoring site, those presenting more than 48 hours after symptom onset, or those patients for whom time of symptom onset was not documented. Results also did not differ materially when we adjusted for apparent temperature or ozone levels (eTable 2). The results from an unconstrained distributed lag model were also similar.
Figure 2 shows that the OR of stroke onset was elevated immediately, peaked in association with mean PM2.5 levels 12 to 14 hours earlier (OR comparing the 75th to 25th percentile of PM2.5, 1.10 [95% CI, 1.04-1.17] [P = .001]), and decreased thereafter.
We considered the association between the risk of stroke onset and other pollutants (Table 2). Results for black carbon concentrations (either measured at the Harvard ambient monitoring site or estimated at each patient's home address) and NO2 levels were similar to those for PM2.5 levels, while results for CO, O3, and SO42- concentrations were not statistically significant.
We examined whether the association between PM2.5 levels and stroke onset varied across subgroups of patients with specific clinical characteristics (Table 3). Fine particulate matter air pollution levels were associated with stroke onset among patients with ischemic stroke classified as due to large-artery atherosclerosis and small-vessel disease but not cardioembolism. There was no evidence suggesting that the association varied according to the presence of comorbid diabetes mellitus, atrial fibrillation or hypertension, history of stroke or transient ischemic attack, or age.
In this Boston area study conducted while the air quality in the region was in attainment of EPA regulatory standards, we found that ischemic stroke risk was 34% higher (95% CI, 13%-58%) (P < .001) on days with moderate PM2.5 levels compared with days with good levels, according to the EPA's Air Quality Index.31 The relationship between higher PM2.5 levels and increased risk of stroke onset was linear, strongest within 12 hours of PM2.5 exposure, and observed among patients presenting with strokes classified as due to large-artery atherosclerosis or small-vessel occlusion but not cardioembolism. Stroke risk was more strongly associated with concentrations of black carbon and NO2, markers of traffic pollution, than with components linked to nontraffic sources.
From a public health perspective, we observed an association between PM2.5 levels and ischemic stroke onset in an area in attainment of the US National Ambient Air Quality Standards.33 Although the observed relative risk is modest, the number of strokes attributable to PM2.5 may be high given the high incidence of ischemic stroke and the fact that nearly everyone is exposed to ambient fine particulate matter. If the association observed in this study is causal, and if a linear dose-response function is assumed, a 2 μg/m3 reduction in mean PM2.5 levels (approximately 20%) during this time period might have averted approximately 6100 of the 184 000 stroke hospitalizations observed in the US Northeast region in 2007 alone.34
An analysis of Medicare beneficiaries in 204 US counties found a 0.4% (95% CI, 0.0%-0.9%) higher risk of admission for the combined end point of cerebrovascular disease per 10-μg/m3 increase in same-day PM2.5 levels.13 Smaller studies in Taiwan14 and Southern California15 have reported excess relative risks of approximately 2% per 10-μg/m3 increase in PM2.5 levels, while others have found no evidence of an association.16- 18 Of 4 prior studies that have specifically evaluated the association between PM2.5 concentrations and the risk of ischemic stroke, 2 Canadian studies found no association,20,22 while studies from Taiwan14 and Nueces County, Texas,21 found a 3% (95% CI, −1% to 7%) and 6% (95% CI, −1% to 13%) excess risk per 10-μg/m3 increase in PM2.5 levels, respectively.
The estimate from the current study scaled to a 10-μg/m3 increase in PM2.5 levels would be 18% (95% CI, 5%-34%), substantially larger than that of previous reports. We believe that this difference is attributable, at least in part, to the use in the current study of data on the timing of ischemic stroke onset. Most prior studies have assessed exposure to PM2.5 based on the calendar day of hospital admission, an exposure assessment strategy that can bias health effects estimates toward the null by as much as 60%.35 Our group's estimates in a previous study in a quality of care registry in Ontario, Canada,22 may have been biased toward the null partly by misclassification of stroke onset time for patients presenting later than the limited time when they might have benefited from thrombolytic therapy. Differences in outcome assessment methods, population and pollutant characteristics, and other aspects of the exposure assessment strategy very likely also contribute to heterogeneity across studies. For example, health effects of PM2.5 concentrations in North America are known to vary geographically depending on the local pollutant sources and components36 as well as community characteristics.37
In the current study we were able to estimate the time course of the association between PM2.5 levels and stroke onset in greater detail than has been previously possible. We observed an immediate increase in the OR for stroke onset that peaked 12 hours after PM2.5 exposure and decreased thereafter (Figure 2). Experimental studies in humans and animals have shown that exposure to concentrated ambient PM2.5 can induce increases in blood pressure and heart rate and reductions in heart rate variability within this time frame,7,38 suggesting that altered hemodynamics could play an important role. Other potential mechanisms include alterations in hemostatic factors, systemic inflammation, endothelial cell injury, and vascular dysfunction.4- 6 Although these physiologic intermediate factors have typically been investigated in association with PM2.5 exposures lasting a day or longer, there is some evidence suggesting that these effects may also follow exposures shorter than 24 hours.7,39- 41
The observation that PM2.5 exposure was more strongly associated with stroke onset in patients with strokes due to large-artery atherosclerosis is consistent with a mechanism involving altered hemodynamics and/or vascular dysfunction that results in disruption of a vulnerable atherosclerotic plaque with subsequent thrombosis and/or downstream embolism, as well as with results from a prior study.22 This result is also consistent with a Boston area study showing a higher risk of acute myocardial infarction within 2 hours of exposure to elevated levels of PM2.5.42 The mechanisms underlying the observed association between PM2.5 exposure and small-vessel stroke are less clear because the pathologic characteristics of these strokes remains poorly understood. However, evidence suggests that endothelial dysfunction and injury, potentially triggered by exposure to PM2.5 or its components, may contribute to the distinct nonatherosclerotic arteriopathy that likely underlies many small-vessel strokes.43 We did not find evidence to suggest that the presence of comorbid diabetes, hypertension, atrial fibrillation, or a history of stroke increased patients' vulnerability to PM2.5-related stroke.
Identifying the components or sources of PM2.5 responsible for the observed associations is of public health and regulatory interest. We found that the risk of stroke onset was most strongly associated with PM2.5 exposure, but also significantly associated with exposure to black carbon and NO2, markers of traffic pollution. This finding is in agreement with past studies suggesting that traffic pollution may trigger ischemic strokes.20,44,45 We did not find any association between risk of stroke onset and ozone levels, a secondary pollutant formed from the reaction of oxides of nitrogen and volatile organic compounds in the presence of sunlight. This is in agreement with most,14,16,17,20 but not all,46 prior studies. We also did not find any association between risk of stroke onset and levels of SO42-, consistent with prior studies using administrative data.16,47 In the Northeast, SO42- generally represents regional pollution from coal-fired power plants.48
Our study has some limitations. First, the use of air quality measures from a single monitoring site is expected to lead to exposure misclassification, increasing the width of our CIs but not otherwise biasing our results.49 However, PM2.5 levels measured at this monitoring site have been shown to be strong proxies for personal exposure to particles of ambient origin in communities surrounding Boston.50 In support of this, our results were not materially different when we restricted the analyses to patients living less than 20 km from the monitoring site, and the results were similar when we evaluated the association between black carbon estimated at patients' homes and ischemic stroke onset. Second, we did not study the association between PM2.5 exposure and stroke resulting in death prior to coming to medical attention. Third, in 13% of patients, we were able to estimate the day but not the time of stroke symptom onset, likely resulting in some exposure misclassification and biasing our results toward the null hypothesis of no association. Consistent with this notion, the association between PM2.5 exposure and ischemic stroke was more pronounced after exclusion of these patients from analysis. Fourth, this study involved patients from a single tertiary care center in Boston. Since the effects of PM2.5 exposure likely vary depending on population characteristics, pollution sources, and particle constituents, our results are not necessarily generalizable to other populations or geographic locations.
As in previous studies, important strengths of this study include detailed data on the time of stroke symptom onset22 and patient clinical characteristics.22,45 In particular, our use of data on the time of stroke symptom onset provides novel insights into the mechanisms by which PM2.5 exposure may increase the risk of ischemic stroke onset.
In conclusion, these results suggest that PM2.5 exposure increases the risk of ischemic stroke at levels below those currently considered safe under US regulations. These associations can be observed within hours of exposure and are most strongly associated with pollution from local or transported traffic emissions. If pollution levels decline with regulation, data on timing of stroke onset, patient clinical characteristics, and stroke mechanisms will be essential for proper evaluation of the clinical benefits of pollution control on stroke risk.
Correspondence: Gregory A. Wellenius, ScD, Center for Environmental Health and Technology, Brown University, 121 S Main St, PO Box G-121S, Room 213, Providence, RI 02912 (firstname.lastname@example.org).
Accepted for Publication: October 1, 2011.
Author Contributions: Drs Wellenius and Mittleman had full access to all of the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis. Study concept and design: Wellenius, Schwartz, Schlaug, Gold, and Mittleman. Acquisition of data: Wellenius, Burger, Schwartz, Suh, Koutrakis, Schlaug, and Mittleman. Analysis and interpretation of data: Wellenius, Coull, Schwartz, Koutrakis, Schlaug, and Mittleman. Drafting of the manuscript: Wellenius. Critical revision of the manuscript for important intellectual content: Wellenius, Burger, Coull, Schwartz, Suh, Koutrakis, Schlaug, Gold, and Mittleman. Statistical analysis: Wellenius, Coull, Schwartz, and Mittleman. Obtained funding: Wellenius, Suh, Koutrakis, Gold, and Mittleman. Administrative, technical, and material support: Burger, Suh, Schlaug, Gold, and Mittleman. Study supervision: Burger, Schlaug, and Mittleman.
Financial Disclosure: One or more of the authors have previously received or currently receive funding from the following sources: Health Effects Institute, Boston; Electric Power Research Institute, Palo Alto, California; US EPA; and National Institutes of Health (NIH).
Funding/Support: This project was supported by grants ES015774, ES009825, ES017125, and ES000002 from the National Institute of Environmental Health Sciences (NIEHS), NIH, and grants R832416 and RD83479801 from the US EPA.
Role of the Sponsors: The funding agencies were not involved in the design and conduct of the study; collection, management, analysis and interpretation of the data; and preparation, review or approval of the manuscript.
Disclaimer: The contents of this report are solely the responsibility of the authors and do not necessarily represent the official views of the NIEHS, NIH, or US EPA.
Additional Contributions: Elissa Wilker, ScD, provided helpful suggestions regarding the manuscript.