Association of Sulfur, Transition Metals, and the Oxidative Potential of Outdoor PM2.5 with Acute Cardiovascular Events: A Case-Crossover Study of Canadian Adults.

BACKGROUND: We do not currently understand how spatiotemporal variations in the composition of fine particulate air pollution [fine particulate matter with aerodynamic diameter ≤2.5μm (PM2.5)] affects population health risks. However, recent evidence suggests that joint concentrations of transition metals and sulfate may influence the oxidative potential (OP) of PM2.5 and associated health impacts.
OBJECTIVES: The purpose of the study was to evaluate how combinations of transition metals/OP and sulfur content in outdoor PM2.5 influence associations with acute cardiovascular events.
METHODS: We conducted a national case-crossover study of outdoor PM2.5 and acute cardiovascular events in Canada between 2016 and 2017 (93,344 adult cases). Monthly mean transition metal and sulfur (S) concentrations in PM2.5 were determined prospectively along with estimates of OP using acellular assays for glutathione (OPGSH), ascorbate (OPAA), and dithiothreitol depletion (OPDTT). Conditional logistic regression models were used to estimate odds ratios (OR) [95% confidence intervals (CI)] for PM2.5 across strata of transition metals/OP and sulfur.
RESULTS: Among men, the magnitudes of observed associations were strongest when both transition metal and sulfur content were elevated. For example, an OR of 1.078 (95% CI: 1.049, 1.108) (per 10μg/m3) was observed for cardiovascular events in men when both copper and S were above the median, whereas a weaker association was observed when both elements were below median values (OR=1.019, 95% CI: 1.007, 1.031). A similar pattern was observed for OP metrics. PM2.5 was not associated with acute cardiovascular events in women. DISCUSSION: The combined transition metal and sulfur content of outdoor PM2.5 influences the strength of association with acute cardiovascular events in men. Regions with elevated concentrations of both sulfur and transition metals in PM2.5 should be examined as priority areas for regulatory interventions. https://doi.org/10.1289/EHP9449.

Introduction

Numerous epidemiological studies have documented the acute cardio-respiratory health impacts of outdoor fine particulate air pollution [fine particulate matter with aerodynamic diameter ()] (Achilleos et al. 2017; Atkinson et al. 2014; Yang et al. 2019). To date, most studies have relied on mass concentrations as the primary exposure of interest without considering spatial or temporal differences in particle composition that may affect overall health risks. Because is a complex and dynamic mixture of organic and inorganic components, treating all particle mass concentrations as equally harmful has clear limitations, and doing so may contribute to spatial heterogeneity in estimated health risks (Liu et al. 2019). Recently, an increasing number of studies have incorporated measures of particle oxidative potential (OP) to complement traditional mass-based measurements because oxidative stress is an important mechanism contributing to air pollution health effects (Bates et al. 2019; Rao et al. 2018). Moreover, recent analysis of atmospheric particles suggests that soluble metals and OP peak at the intersection of sulfate and transition metal concentrations in (Fang et al. 2017). Therefore, if soluble metals and OP are important determinants of the magnitude of health risks associated with , it may be necessary to consider both of these components simultaneously in epidemiological analyses.

In Canada, we previously reported evidence of effect modification by particle OP in a panel study of personal exposures and airway inflammation in children who have asthma (Maikawa et al. 2016), in a cohort study of outdoor and preterm birth (Lavigne et al. 2018), and in case-crossover studies of daily and emergency room visits for myocardial infarction (Weichenthal et al. 2016a) and asthma (Weichenthal et al. 2016b). In all of these studies, the short-term health impacts of were greater when glutathione-related oxidative potential () was higher as determined using a cell-free assay based on a synthetic respiratory tract lining fluid. However, other studies have not found evidence of effect modification by OP in evaluating the acute health impacts of (see review by Atkinson et al. 2016). To our knowledge no epidemiological studies to date have examined how the acute cardiovascular health risks of may vary based on combinations of transition metal and sulfur content or how spatiotemporal variations in combinations of these components may be related to particle OP. To address this need, we conducted a national-scale case-crossover study of outdoor mass concentrations and acute cardiovascular events including prospective monthly measurements of sulfur (as a marker for sulfate and indirectly metal solubility) and transition metals in as well as three different metrics of OP (as monthly means).

Methods

Study Design

A time-stratified case-crossover design (Janes et al. 2005) was used to estimate associations between short-term changes in outdoor mass concentrations and hospital admissions for cardiovascular outcomes across strata of OP, mass proportions (percent) of transition metals [i.e., copper (Cu), iron (Fe), nickel (Ni), manganese (Mn), zinc (Zn)], and S (sulfur) in monthly samples (described below). The time-stratified case-crossover design selects reference periods matched on the same day of the week, month, and year as event days (i.e., if a case occurs on the first Friday of June 2021, the references days are all other Fridays in June 2021). Our analysis included the following 34 locations across Canada (Figure 1) (all cities had a single monitor unless otherwise indicated): Red Deer, Edmonton (two sites), St. Albert, Regina, Hamilton, Victoria, Ottawa (two sites), Montreal, Winnipeg, St. John’s, Courtenay, Windsor, Duncan, Saskatoon, Quebec City, Kelowna, Prince George, Kamloops, London, Mt. Pearl, Prince Albert, Nanaimo, Fort Mackay, Whitehorse, Swift Current, Fort McMurray, Brandon, Calgary, Yellowknife, Athabasca Valley, Quesnel, Halifax, Saint John, and Fredericton.

Figure 1.

Locations of monitoring sites for components and oxidative potential across Canada (2016–2017).

All adult cases (ages 18–100 y) of acute cardiovascular events were identified between 1 June 2016 and 31 December 2017 to overlap with the timing of monthly OP/component measurements. Specifically, hospital admissions for cardiovascular outcomes [International Classification of Diseases (ICD), 10th revision, Code I00-I99] were identified using the Discharge Abstract Database (DAD) maintained by the Canadian Institute for Health Information (CIHI) along with information on age, sex, and residential postal code. The DAD captures hospital admissions across Canada with the exception of Quebec (Gibson et al. 2008). Quebec data were obtained from the Quebec Ministry of Health and Social Services through MED-ÉCHO (Ministère de la Santé et des Services Sociaux 2021). Cases were included in case-crossover analyses if at the time of admission their residential 6-digit postal code (about the size of a city block face) was within of a monitoring site for daily (which in most cases was the same location where OP and components were measured). Fatal cases were included in the analysis, but readmissions were excluded (i.e., if the same person was admitted twice, only the first visit was included). Ethics approval for this study was granted through data sharing agreement between Health Canada and CIHI.

Monthly Sample Collection

Integrated samples were collected on Teflon™ filters at each study location for 2 wk each month using cascade impactors operating at a flow rate of per minute. In most cities, monthly samplers were located at provincial monitoring sites (i.e., the same locations that provided daily data) except for Ottawa ( from the location of daily measurements) and Montreal ( from the location of daily measurements) where OP/component measurements were collected outside private residences (located at ground level in Ottawa and approximately 6 meters off the ground in Montreal). Following gravimetric analyses, monthly filters were analyzed for sulfur and transition metal content using X-ray fluorescence (U.S. Environmental Protection Agency Method IO-3.3) and oxidative potential as described below. The following transition metals were selected for inclusion in our analyses based on previous evidence suggesting an association with particle OP: Cu, Fe, Zn, Ni, and Mn (Charrier and Anastasio 2012; Fang et al. 2017; Gao et al. 2020; Verma et al. 2014). It is important to note that although analyses were conducted across strata of specific metals, we cannot conclusively attribute the results to specific transition metals, because trace metals tended to be highly correlated. If more than one sample was collected from a site in a given month, samples were averaged to obtain a single monthly estimate.

OP Measurements

Monthly samples were extracted into high-performance liquid chromatography high-performance liquid chromatography (HPLC) grade methanol by vortexing at 1,800 rpm for 20 min and sonicating for 10 min. Decanted methanol was evaporated under a gentle flow of nitrogen. PM samples were resuspended in ultrapure water containing 5% HPLC methanol to a final storage concentration of PM/mL. Resuspended samples were analyzed in triplicate using three OP metrics: the ascorbate (AA), glutathione (GSH), and dithiothreitol (DTT) assays. Ascorbate and glutathione oxidative potential ( and ) were assessed using the acellular respiratory tract lining fluid (RTLF) OP assay (Maikawa et al. 2016; Weichenthal et al. 2019). Briefly, samples were incubated at a concentration of in a 96-well plate for 4 h at 37°C with synthetic respiratory tract lining fluid (RTLF) containing of each AA, GSH, and uric acid in an ultraviolet-visible plate reader (SpectraMax 190; Molecular Devices) alongside positive controls ( , 0.02% ) and blanks. Ascorbate depletion was calculated over the 4-h incubation period, and GSH depletion was measured using the glutathione-reductase enzyme recycling assay (Baker et al. 1990; Godri et al. 2010; Griffith 1980; Mudway et al. 2004).

Dithiothreitol OP () was assessed using an adapted version of the DTT assay (Cho et al. 2005). Briefly, resuspended samples were incubated with DTT in a 96-well plate alongside positive controls [ ], blanks, and DTT standards (containing DTT) for 35 min at 37°C, with constant shaking. After 5, 15, 25, and 35 min, the remaining DTT was measured by adding 5,5′-dithiobis(2-nitrobenzoic acid) (DTNB) to each well and measuring absorbance at . Samples were initially analyzed at a concentration of ; however, if DTT depletion exceeded 25% after 35 min, the sample was reanalyzed at a lower concentration. All OP values are expressed in units of picomoles per minute per microgram; values below detection were replaced with half the limit of detection.

Daily Air Pollution Data

Daily mean outdoor concentrations of , and were obtained from fixed-site monitoring stations operated by the National Air Pollution Surveillance network maintained by Environment and Climate Change Canada. was calculated as a weighted average of and using redox potentials as the weights {i.e., , where V refers to volts; Weichenthal et al. 2016a}. The redox-weighted measure accounts for the fact that is a stronger oxidant than . Case and control periods were assigned outdoor concentrations based on these fixed-site monitoring data. Daily mean temperature and relative humidity data were also provided by Environment and Climate Change Canada using mean 24-h values from the closest weather station.

Statistical Analyses

Conditional logistic regression models were used to estimate associations [odds ratios (OR) and 95% confidence intervals (CI)] between lag-0 (i.e., same day) concentrations (including the case day) and hospital admissions for cardiovascular outcomes adjusting for continuous measures of lag-0 mean temperature and concentrations. This lag period was selected based on past analyses of the acute cardiovascular health impacts of that suggested that same-day exposures were most strongly associated with the outcome (Weichenthal et al. 2016a). Natural cubic splines with 3 knots were also examined for temperature and but gave results similar to those of models including linear terms (Table S1), likely owing to the restricted range of values across within-subject comparisons (i.e., the exposure variance of interest is case-crossover studies relates to the distribution of pollutant differences between case and control periods, which is generally much less than day-to-day variations across the study period) (Künzli and Schindler 2005). Relative humidity was not correlated with outdoor () and was not included as a confounder in the analyses. Lag-0 concentrations were not correlated with lag-0 concentrations () but were included a priori as a covariate in our models because we previously observed an association between daily and acute myocardial infarction (Weichenthal et al. 2016a). Lag-0 concentrations were not correlated with monthly average OP () or mass proportions of transition metals (i.e., Cu, Fe, Zn, Ni, Mn) in (). Other time-varying confounding factors were not examined in the models. Models for 3-d mean concentrations were evaluated, but results were similar (but slightly attenuated) to those for lag-0, and thus we focus only on results for lag-0 concentrations (Table S2). A cluster variance estimator was used to account for potential within-city clustering of observations.

All analyses were conducted separately for men and women, because initial analyses looking only at mass concentrations suggested stronger associations among men (Table S2; interaction -value for lag-0 ). The interaction between and sex was evaluated by including a product term in conditional logistic regression models; a -value less than 0.05 was considered statistically significant. Among men and women, we also examined models stratified by age (i.e., above/below the median age of 70 y) but did not perform formal tests for effect modification by age.

Stratified analyses were conducted for lag-0 across strata defined by monthly measurements of , , , metals, and S in . To maximize the number of cases in each cell (between approximately 10,000–18,000 for men and 7,000–12,000 for women), strata for transition metals, OP, and S were based on median values in monthly samples [this provided a sample size 2–3 times greater than a recent case-crossover study of the acute cardiovascular health impacts of in Canada at similar concentrations (Weichenthal et al. 2017)]. These analyses provided estimates of associations between lag-0 concentrations and acute cardiovascular outcomes across various combinations of monthly average S and metals/OP. Using these models, we also calculated -values for interactions between lag-0 and monthly S (by including an interaction term between lag-0 and an indicator variable for S) separately for strata defined by median metals/OP concentrations. This calculation was done to determine whether the pattern of effect modification by S also depended on metal/OP concentrations. A third strata for S was also evaluated ( percentile, including between 1,200 and 3,500 cases) to examine patterns on the upper end of the S distribution but was not included in formal tests of interactions. Similarly, we calculated -values for interactions between lag-0 and monthly metals/OP (by including an interaction term between lag-0 and indicator variables for metals/OP) separately for strata defined by median S concentrations. This calculation was done to determine whether the pattern of effect modification by metals/OP also depended on S concentrations. Finally, three-way interactions were also examined including interaction terms between lag-0 and indicators variables (above/below monthly median) for metals/OP (separately) and S. All ORs reflect a change in outdoor to facilitate comparisons with other studies of outdoor . All statistical analyses were conducted using R (version 4.0; R Development Core Team). Three-dimensional plots were generated using generalized additive models in the mgcv package in R with three knots (for each dependent variable) (Wood 2021).

Results

In total, 93,344 cases of acute cardiovascular events were included in the analysis; 3,765 eligible cases (3.9%) were excluded because of missing data for daily . Cases were predominantly male (54,338 men and 39,006 women), with a mean age of 70 y []. Daily mean concentrations were low (; ), but monthly OP, transition metal, and S concentrations varied substantially across the study period (Table 1; descriptive data for element mass proportions are shown in Table S3). Distributions of monthly components (i.e., transition metals and S) were based on 685 samples, whereas 636 samples were available for , 637 for , and 634 for .

Table 1

Descriptive statistics for daily data used in case-crossover analyses and monthly estimates of (micrograms per cubic meter), oxidative potential [glutathione (), ascorbate (), and dithiothreitol-related oxidative potential ()], and components (Canada, 2016–2017).

Pollutant/component	Mean	SD	5th	25th	50th	75th	95th
Daily air pollutants and temperature
PM_2.5 (μg/m^3)	7.25	6.86	1.79	3.83	5.87	8.83	16.5
O_3 (ppb)	23.1	8.37	9.70	17.4	22.8	28.9	36.9
NO_2 (ppb)	8.54	5.85	2.00	4.46	7.17	11.0	20.0
O_x (ppb)	18.2	5.03	10.3	14.7	17.9	21.6	26.6
Temperature (°C)	7.68	11.2	−13.4	0.0404	9.14	16.9	22.4
Monthly average components
PM_2.5 (μg/m^3)	7.67	5.47	3.27	5.17	6.82	9.08	14.3
OP (pmol/min/μg)
OP^GSH	3.45	1.98	1.03	2.09	3.09	4.35	7.62
OP^AA	2.91	0.91	1.46	2.29	2.79	3.50	4.67
OP^DTT	11.9	7.46	2.12	6.05	10.2	16.6	26.5
Elements (ng/m^3)
S	261	139	97.6	162	242	314	520
Cu	3.18	3.08	0.501	1.06	1.95	4.28	9.00
Fe	93.7	61.1	20.6	43.3	83.4	130	205
Ni	0.460	0.81	0.0554	0.121	0.207	0.419	2.12
Mn	3.26	2.67	0.555	1.34	2.82	4.19	7.52
Zn	12.4	13.3	1.99	3.93	6.83	15.6	39.4

Note: Cu, copper; Fe, iron; Mn, manganese; Ni, nickel; OP, oxidative potential; , ascorbate; , dithiothreitol; , glutathione; ppb, parts per billion; S, sulfur; SD, standard deviation; Zn, zinc.

As shown in Figure 2, correlations between metals and OP were stronger when S concentrations were higher. For example, when S was below the median, was only weakly correlated with Fe (), Cu (), Mn(), Zn (), and Ni (). When S concentrations were above the 90th percentile, these correlations increased substantially: Fe (), Cu (), Mn (), Zn (), and Ni (). A similar pattern was observed for , with stronger correlations between and metals at higher S concentrations. For , correlations with metals did not change as dramatically across levels of S content, but was more strongly correlated with when S was higher ( vs. ) (Figure S1). Figure 3 highlights how and values varied across values of both metal (i.e., Fe, Cu, Zn) and S content in . These relationships were generally nonlinear with concave shapes for both and except for the plot for , Cu, and S, which was nearly linear with a steeper slope observed between and Cu than between and S.

Figure 2.

Spearman correlations between monthly mean oxidative potential [glutathione (), dithiothreitol (), and ascorbate () depletion (pmol/min/μg)] and transition metal concentrations (Canada, 2016–2017) [copper (Cu), iron (Fe), manganese (Mn), nickel (Ni), and zinc (Zn) in nanograms per cubic meter] at high () (A) and low () (B) sulfur (S) concentrations.

Figure 3.

Relationships between glutathione-related oxidative potential [ (picomoles per minute per microgram)] and mass proportions (percent) of iron (Fe) (A), copper (Cu) (B), zinc (Zn) (C), and sulfur (S) in when the mass proportion of S is above the median (Canada, 2016–2017). Plots for ascorbate-related oxidative potential () are shown in panels (D), (E), and (F).

Overall, daily mean mass concentrations were weakly associated with acute cardiovascular events (, 95% CI: 1.005, 1.025) with positive associations limited to men (, 95% CI: 1.014, 1.036) and no association observed among women (, 95% CI: 0.979, 1.019) (interaction -value for lag-0 and ) (Table S2). Among men, the magnitude of this risk was similar among younger () (, 95% CI: 1.004, 1.042) and older men (, 95% CI: 1.004, 1.052) (Table S2). When analyses were conducted across strata of transition metal and S content in , a consistent pattern of stronger associations was observed among men when both transition metal and S content were elevated. For example, an OR of 1.078 (95% CI: 1.049, 1.108) (per ) was observed for cardiovascular events in men when both Cu and S were above the median, whereas a weaker association was observed when both elements were below median values (, 95% CI: 1.007, 1.031) (Table S4). This trend is shown in Figure 4 for Fe, Cu, and Zn, and complete results are shown in Table S4. When metal concentrations were above median values, -values for interaction terms between lag-0 and an indicator variable for S content (i.e., above/below median values) were 0.076 for Mn, 0.045 for Cu, and 0.00037 for Zn (Table S4). S content did not modify associations between and acute cardiovascular events in men when metal concentrations were below median values (Table S4). A similar pattern was not observed in women, and in several cases was inversely associated with hospital admissions for acute cardiovascular outcomes among women (Table S5). Within strata of S (i.e., above/below median values), interaction -values between lag-0 and metals among men were as follows: Cu (high S: 0.12; low S: 0.67); Fe (high S: 0.31; low S: 0.85); Mn (high S: 0.33; low S: 0.50); Ni (high S: 0.24; low S: 0.59); Zn (high S: 0.39; low S: 0.024). In models examining three-way interactions between lag-0 , metals (above/below median), and S (above/below median) in men, -values for the three-way interaction terms were as follows: Cu (0.22), Fe (0.61), Ni (0.93), Mn (0.37), and Zn (0.054).

Figure 4.

Variations in the strength of associations between lag-0 (per ) and acute cardiovascular events among men across categories of mass proportions of sulfur (S) and copper (Cu) (A) and (B), iron (Fe) (C) and (D), and zinc (Zn) (E) and (F) in . All conditional logistic regression models are adjusted for lag-0 mean temperature and . High and Low refer to mass proportions above () and below () median values. See Table S4 for corresponding numeric data.

When analyses were conducted across strata of OP and S, a similar pattern was observed for all three assays, with the strongest associations among men occurring when both OP and S were elevated (there was no pattern among women) (Figure 5; Tables S4 and S5). When OP was above the median, S content modified associations between and acute cardiovascular events among men for all three assays with interaction -values of 0.010, 0.0052, and 0.0067 observed for , , and , respectively. The strongest association was observed in the 90th percentile of S content when was high (; 95% CI: 1.095, 1.237); was weakly associated with acute cardiovascular events when both S content and were below median values (; 95% CI: 1.003, 1.027). Similar patterns were also observed for and (Table S4). Interaction -values for and OP (above/below median) within strata of S mass proportions (above/below median) were as follows: (high S: 0.066; low S: 0.74); (high S: 0.59; low S: 0.12); (high S: 0.18; low S: 0.32). Little evidence of effect modification by , , or was observed among men in analyses not considering S content, although stronger associations were observed between and acute cardiovascular events among men when ( above median: ; 95% CI: 1.013, 1.065; below median: ; 95% CI: 1.006, 1.029) and ( above median: ; 95% CI: 1.015, 1.042; below median: ; 95% CI: 0.984, 1.046) were elevated (Table S6). In models examining three-way interactions between lag-0 , OP (above/below median), and S (above/below median) in men, -values for the three-way interaction terms were as follows: (0.16), (0.58), (0.13).

Figure 5.

Variations in the strength of associations between lag-0 (per ) and acute cardiovascular events among men across categories of mass proportions of sulfur (S) in and glutathione () (A) and (B), dithiothreitol () (C) and (D), and ascorbate-related oxidative potential () (E) and (F) (Canada, 2016–2017). All conditional logistic regression models are adjusted for lag-0 mean temperature and . High and Low refer to values above () and below () median values. See Table S4 for corresponding numeric data.

Finally, given the observed pattern of stronger associations between outdoor and acute cardiovascular events in men when both transition metal and S content were elevated, we also examined how the proportions of these components related to overall mass concentrations. This question was of interest because if these proportions are not constant across mass distributions, they could play a role determining the shape of concentration–response curves. These results are shown in Figure 6 and indicate that for a given mass proportion of metals in (Fe and Cu in this example) the mass proportion of S increases as mass concentrations decrease. To verify that this pattern was not unique to Canada, we also downloaded data for mass concentrations, Fe, Cu, and S for the year 2017 from the Interagency Monitoring of Protected Visual Environments (IMPROVE) network in the United States (1,836 monthly averages values calculated using daily data) (IMPROVE 2021). As shown in Figure S2, for a given mass proportion of Cu/Fe, S mass proportions in also increased as mass concentrations decreased in the U.S. data.

Figure 6.

Generalized additive model plots for the mass proportion (percent) of sulfur (S), iron (Fe) (A), copper (Cu) (B), and (Canada, 2016–2017). mass concentrations (micrograms per cubic meter) are plotted up to the 99th percentile of monthly measurements.

Discussion

is a complex mixture that varies in composition across both space and time. In this study, we examined the acute cardiovascular health impacts of mass concentrations across strata of metals, OP, and S. Among men, outdoor mass concentrations were most strongly associated with acute cardiovascular outcomes when mass proportions of both metals and S were elevated. A similar pattern was observed for OP and S. Moreover, OP appears to capture heterogeneity in health effects among men because OP is more strongly associated with metals when the S content of is higher. was not associated with acute cardiovascular outcomes among women.

Recently, an increasing number of epidemiological studies have examined OP as a complementary metric to traditional particle mass concentrations (Bates et al. 2019; Gao et al. 2020; Molina et al. 2020). To date, few studies have examined the role of oxidative potential in modifying the acute health impacts of . Specifically, we previously reported that modified associations between outdoor mass concentrations and acute myocardial infarction (Weichenthal et al. 2016a), asthma (Weichenthal et al. 2016a), and airway inflammation in children with asthma (Maikawa et al. 2016). Similarly, two time-series studies noted stronger associations between and acute respiratory outcomes than for (Bates et al. 2015; Abrams et al. 2017), but other studies have failed to replicate this finding (Atkinson et al. 2016). With respect to composition, epidemiological evidence to date has been somewhat inconsistent (Schlesinger 2007; Lippmann 2014; Adams et al. 2015), and we have yet to conclusively identify specific components/sources of that are most relevant to health. Nevertheless, a recent analysis of atmospheric particles suggests that combined transition metal and S content may influence the OP of because sulfate increases particle acidity, which in turn increases metal dissolution and solubility, thus allowing metals to participate in redox reactions contributing to oxidative stress (Fang et al. 2017). Together, these mechanisms may help to explain spatiotemporal variations in the adverse health effects of , but epidemiological evidence in this area remains limited. In this study, we examined the combined impacts of transition metal, OP, and S content (as a marker for sulfate and indirectly metal solubility) on associations between outdoor mass concentrations and the risk of acute cardiovascular events and noted several interesting findings.

First, the magnitudes of associations between daily mean outdoor and acute cardiovascular outcomes in men were stronger when both transition metal and S content were elevated. In addition, evidence for effect modification was clearest across categories of S content when the transition metal content of was elevated (particularly for Cu, Mn, and Zn). A similar pattern was apparent for OP, with stronger associations observed among men when both OP and S content were higher. Moreover, analysis of correlations between OP and transition metals at high and low S content suggested that OP and metals were more strongly correlated when S content was elevated. Therefore, our findings suggest that OP metrics may capture heterogeneity in the acute cardiovascular health impacts of at least in part because they are more strongly associated with metals when S content is higher (i.e., when metal content appears to be most relevant to health). This finding was particularly true for and and is consistent with previous evidence suggesting that these metrics are more strongly correlated with metals in (e.g., Cu, Fe, Ni, Zn), whereas is more strongly related to organic components (e.g., organic carbon, water-soluble organic carbon) (Gao et al. 2020). In generalized additive model plots, nonlinear relationships were observed between transition metals, S, and OP, often with a convex shape (i.e., with OP peaking at certain mass proportions of metals and S). Although a detailed examination of particle chemistry was beyond the scope of this investigation, these plots suggest that component interactions influence overall OP in a complex manner, and further investigation is needed to better understand these relationships because they may be useful in guiding future regulatory interventions.

Our analysis of transition metal and S content across the distribution of mass concentrations indicates that for a given mass proportion of transition metals the mass proportion of S increases as mass concentrations decrease. This finding could have important implications for future analyses of the shapes of concentration–response relationships for outdoor mass concentrations. Specifically, if combined transition metal and S content are important determinants of -associated health impacts, one would expect to observe steeper increases in concentration–response curves at lower concentrations. Indeed, patterns of steeper increases in concentrations–response curves at lower concentrations have been reported around the world in time-series analyses of outdoor and mortality (Liu et al. 2019) as well as in cohort studies of long-term exposures to outdoor (Burnett et al. 2018; Christidis et al. 2019; Crouse et al. 2012, 2015; Pappin et al. 2019; Pinault et al. 2016, 2017), including pooled cohort study estimates used to estimate the global burden of disease attributable to outdoor air pollution (GBD 2019 Risk Factor Collaborators 2020). Although we could not specifically address this question in the current study, our findings combined with existing patterns of steeper concentration–response relationships at lower mass concentrations in existing studies of outdoor is intriguing. If confirmed, this hypothesis could shed new light on why we continue to see adverse health impacts of outdoor at low mass concentrations and why these relationships tend to be stronger at the lower end of exposure distributions.

The combined importance of transition metals and sulfur in is also interesting from the perspective of public health interventions. For example, it may be helpful to systematically identify regions (and time periods within regions) with elevated mass proportions of both components as a first step in prioritizing areas for future risk management activities. Moreover, it is possible that some threshold may exist for acute cardiovascular health impacts with respect to the relative proportions of S and transition metals in , but we did not have the power to examine this question in the present study. Sources of metals in are diverse and include emissions from vehicles, industry, shipping, and brake wear (Bates et al. 2019; Birmili et al. 2006; Lough et al. 2005; Quiterio et al. 2004), whereas S (i.e., sulfate) primarily comes from point sources of fossil fuel combustion (e.g., coal-fired power plants for electricity generation) (Manoli et al. 2002; Moldanová et al. 2009; Oeder et al. 2015). The efficiency of targeting sources of metals or sulfate (or both) in a given region will depend on various logistical and economic considerations. However, in many cases the diverse nature of metal sources may make reductions in S concentrations more feasible. Importantly, if replicated in future analyses, our findings suggest that regulatory actions aimed at reducing sulfate may reduce the acute cardiovascular health impacts of even if metal concentrations remain unchanged. Indeed, numerous previous regulations have targeted sources of S in , including limits on sulfur in gasoline (U.S. EPA 2021) and marine fuel oil (International Maritime Organization 2019), and our findings add strength to the justification for these regulatory actions from an environmental health perspective. More broadly, widespread integration of composition/OP data into the risk assessment/risk management process will require new exposure models capable of estimating these parameters over broad spatiotemporal scales. Although models are currently available to predict large-scale variations in components using remote sensing/chemical transport models (van Donkelaar et al. 2019; Xu et al. 2019), large-scale models for oxidative potential are only beginning to emerge (Daellenbach et al. 2020). Nevertheless, an increasing number of large-scale exposure models for OP are likely to be developed in the coming years as ground-level data support continues to increase over broad geographic areas.

Although this study had a number of important strengths, including prospective data collection for multiple OP metrics and components across Canada, it is important to note several limitations. First, it was not possible to prospectively measure daily composition at multiple sites across Canada, and thus we relied on estimates of monthly mean composition and OP based on repeated 2-wk samples collected at each location. As a result, it is possible that a given month was misclassified with respect to composition and OP if the 2-wk measurement period was not representative of the true monthly average value. Moreover, this error could affect observed trends in risk estimates across strata of metals, S, and OP, making it more difficult to identify true patterns across strata (although findings were generally consistent across strata of all the components examined). Similarly, we could not identify which specific metals may be more or less harmful because all of the metals tended to be correlated and all demonstrated a similar pattern of increased risks at elevated metal and S content in . Likewise, we did not measure all possible time-varying factors potentially related to the outcome (i.e., stress), but it seems unlikely that these would be systematically related to daily consistently within strata of metals/OP and S. In addition, the nature of our hypothesis required many models to be examined, and chance findings could have affected our results. However, given the consistency of the results among men, this seems like an unlikely explanation for the overall patterns observed.

A second limitation is related to our use of all cardiovascular outcomes as opposed to specific types of cardiovascular events (e.g., myocardial infarction). The use of all cardiovascular outcomes was necessary, given the relatively short time period under investigation, the fact that cases were limited to those living within of a monitor (to reduce exposure misclassification, given high spatial variability for OP and metals) (Weichenthal et al. 2019), and the fact that monitoring was limited to many small to midsize cities that do not give rise to large numbers of cases. As a result, if only contributes to an increased risk of specific types of acute cardiovascular events (e.g., Weichenthal et al. 2017), our analyses may underestimate the true magnitude of association between and these specific outcomes. This aspect may be particularly relevant to our null findings among women because the use of all cardiovascular events may mask important associations with specific types of cardiovascular events among women. Indeed, although we observed clear and consistent patterns of associations between outdoor and acute cardiovascular outcomes among men, was not associated with acute cardiovascular events in women. We do not have an obvious explanation for this result, but it was consistent in all of the models we examined. Bell et al. (2015) also reported contrasting results between sexes (with a stronger association among women), but more generally our results highlight that sex-/gender-specific analyses should be more common in air pollution epidemiology because improved understanding of heterogeneity across sex/gender profiles may help to improve public health interventions and/or patient-level health information (Clougherty 2010).

Finally, our study was limited to the acute health impacts of , whereas long-term exposures are most important for overall burden of disease (GBD 2019 Risk Factor Collaborators 2020). However, our results highlight substantial spatiotemporal variations in outdoor oxidative potential and composition across Canada, and future cohort studies could leverage this information to explore possible effect modification for outcomes, including cancer incidence and cause-specific/nonaccidental mortality. We previously addressed this question in a cohort study conducted in Ontario, Canada, and noted a stronger association between lung cancer mortality and when was higher (Weichenthal et al. 2016c); however, future studies should aim to replicate these results on a larger scale.

In summary, our findings suggest that the mass proportions of transition metals and S play an important role in determining the strength of association between outdoor and the risk of acute cardiovascular events in men. Moreover, our results indicate that OP metrics capture this trend in part because they are more strongly correlated with transition metals in when S content is higher. These results provide new information on why we continue to see adverse health effects of at low mass concentrations. Specifically, at a given mass proportion of metals, the mass proportion of S increased as decreased; this finding may help to explain the repeated observation of steeper concentration–response relationships at the lower end of exposure distributions for . Identifying regions with elevated levels of both transition metal and S content in may be an efficient means of prioritizing areas and sources for future regulatory interventions.

Click here for additional data file.

Click here for additional data file.