Boxplots included baseline air pollution concentrations in 2000 or for the years of 2006-2008 (for black carbon) and follow-up air pollution exposures aggregated from the year of the baseline examination to that of the follow-up clinic examination from 2000 to 2018 for each participant at each scan. The boxes indicate the interquartile range (IQR) and the line in the center indicates the median concentration. Whiskers extend to 1.5 times the IQR to the most distant observation within that distance (ie, the largest or smallest observations within quartile 3 + 1.5 × IQR or quartile 1 − 1.5 × IQR). Outlier observations are shown as circles. Exposure data are shown by study areas, which included 1037 participants (2837 scans) in Winston-Salem, NC; 1182 (3227 scans) in New York, NY, and Rockland County, NY; 1054 (2842 scans) in Baltimore, MD; 1021 (2933 scans) in St Paul, MN; 1142 (3396 scans) in Chicago, IL; and 1424 (3339 scans) in the Los Angeles Basin area.
Metropolitan area temporal trends in air pollution exposures are shown based on annual mean concentrations for all participants in each area from continuous 2-week average concentrations from 2000 through follow-up. The exposure data included 1037 participants in Winston-Salem, NC; 1182 in New York, NY, and Rockland County, NY; 1054 in Baltimore, MD; 1021 in St Paul, MN; 1142 in Chicago, IL; and 1424 in the Los Angeles Basin area.
A, The linear longitudinal association of baseline exposure to O3, PM2.5, and NOx in 2000 or 2006 to 2008 for black carbon (18 902 scans) and mean follow-up air pollution exposures from the year of the baseline examination to that of the follow-up clinic examination from 2000 to 2018 (18 574 scans) for each scan. Progression of percent emphysema shown via computed tomographic (CT) imaging (percentage of lung pixels less than −950 Hounsfield units) over 10 years shown from linear mixed model. Associations were presented per rounded interquartile range increment of exposure to O3, PM2.5, NOx, and black carbon. Models are adjusted for age, sex, race/ethnicity, height, weight, temperature, smoking status, secondhand smoking, pack-years, cigarettes per day, body mass index, physical activity, income, employment status, education, neighborhood socioeconomic status (SES) index, study region, the interaction between SES and study region, and interactions of these variables with time as well as variables without interactions with time, such as CT scanners, pixel size, and baseline exposure. B, The concentration-response curve with 95% CI for the overall change of percent emphysema (progression rate associated with O3 concentrations over follow-up) with the same model adjustments as in panel A. The curve is based on a mixed model that includes a thin plate regression spline with 4 degrees of freedom to more flexibly assess the potentially nonlinear relationship. The highest and lowest 1% of overall concentrations have been trimmed for visualization because the relationship at the extremes is uncertain and might rely on concentrations from a single region. The histogram shows the distribution of O3 concentration over follow-up in this cohort, with the contribution of each area to the total distribution indicated by color; the height of each bar in the histogram represents the number of scans analyzed in each 0.5 parts per billion (ppb) bin of O3 concentration.
Associations are presented as rounded interquartile range (IQR) increment of exposure to O3, PM2.5, NOx, and black carbon (BC). A, Results of a single-pollutant model derived from main analyses for the associations between each of the air pollutants and progression of percent emphysema (percentage of lung pixels less than −950 Hounsfield units) are shown. The multipollutant model presents the associations of the fully adjusted model with all measured air pollutants modeled simultaneously. Linear combination models were implemented by combining the associations from any pairs of the air pollutants, based on associations from the multi-pollutant model.
A, The associations with baseline air pollution exposure in 2000 or 2006 to 2008 (for black carbon) in 18 902 scans. B, The associations with follow-up air pollution exposures aggregated from the year of the baseline examination to that of the follow-up clinic examination from 2000 to 2018 in 18 574 scans. Details of the exact numbers for the associations from the single-pollutant, multipollutant and linear combination models are shown in eTable 6 in the Supplement.
This Video Abstract reviews a cohort study examining associations of long-term exposure to air pollutants with progression of CT-assessed percent emphysema and with decline in lung function in adults in 6 US metropolitan regions.
eMethods. Statistical modeling
eTable 1 Unadjusted percent emphysema progression and lung function decline per year by quartile of O3 exposure concentrations in the MESA cohort
eTable 2 Descriptive characteristics [mean (SD) or %] of study participants at baseline and unadjusted longitudinal outcomes by study areas in the analysis of lung function
eTable 3 Pearson correlation coefficients between predicted air pollution exposures at the MESA participants’ homes within the MESA cities (N=6860)
eTable 4 Main and sensitivity analyses of effect estimates (95% CIs) of air pollution concentrations, assessed at baseline, with longitudinal changes of percent emphysema on CT, per 10 years from staged models
eTable 5 Main and sensitivity analyses of effect estimates (95% CIs) of long-term air pollution concentrations, assessed during follow-up, with longitudinal change of percent emphysema on CT over 10 years from staged models
eTable 6 Effect estimates (95% CIs) for the associations between air pollutant exposures at baseline or over follow-up and progression of percent emphysema in single-pollutant, multiple-pollutant and linear combination models, per IQR increment for O3 (3 ppb), PM2.5 (2 µg/m3), NOx (10 ppb) and black carbon (0.2 µg/m3).
eTable 7 Longitudinal changes in percent emphysema (95% CI) per over ten years per increase of O3, PM2.5, NOx and black carbon, assessed at baseline, stratified by potential effect modifiers
eTable 8 Longitudinal changes in percent emphysema (95% CI) over ten years per increase of long-term O3, PM2.5, NOx concentration, assessed during following, stratified by potential effect modifiers
eTable 9 Multiplicative interactions between air pollutant exposures for effect of air pollution on progression of percent emphysema over ten years (N=6860)
eTable 10 Main and sensitivity analyses of effect estimates (95% CIs) of long-term air pollution concentrations, assessed during following with longitudinal changes of lung function over 10 years from staged models
eTable 11 Longitudinal changes in lung function (95% CI) per 3 ppb increase of O3 assessed during following over 10 years stratified by personal factors
eTable 12 Main and sensitivity analyses of effect estimates (95% CIs) of exposure to air pollution concentrations assessed at baseline with longitudinal changes of lung function over 10 years from staged models
eFigure 1 Participant recruitment, retention, and flow of outcome testing in the analysis of percent emphysema and lung function in the Multi-Ethnic Study of Atherosclerosis.
eFigure 2 Concentration-response curves with 95%CI for the overall change of percent emphysema (progression rate associated with air pollution concentrations, assessed at baseline and over follow-up (N=6860)
eFigure 3 Concentration-response curve with 95%CI for the change of FEV1 (rate of change associated with air pollution concentrations, assessed at baseline and over follow-up, N=3636).
eFigure 4 Effect estimates for the associations between air pollutants and lung function decline (N=3636). Results from single-pollutant model, multiple-pollutant model and linear combination model for the effect estimates of multiple air pollutant exposures (A: at baseline for FEV1, B: over follow-up for FEV1, C: at baseline for FVC, and D: over follow-up for FVC), from multi-pollutant models.
Customize your JAMA Network experience by selecting one or more topics from the list below.
Wang M, Aaron CP, Madrigano J, et al. Association Between Long-term Exposure to Ambient Air Pollution and Change in Quantitatively Assessed Emphysema and Lung Function. JAMA. 2019;322(6):546–556. doi:10.1001/jama.2019.10255
Is there an association between ambient air pollutants and progression of emphysema and changes in lung function in the general population?
In this cohort study conducted between 2000 and 2018 that included 5780 participants in 6 US metropolitan regions followed up for a median of 10 years, there was a statistically significant association between baseline ambient concentrations of ambient ozone (O3), fine particulate matter (PM2.5), oxides of nitrogen (NOx), and black carbon with greater increases in emphysema assessed quantitatively using computed tomographic (CT) imaging. Concentrations of O3 and NOx, but not concentrations of PM2.5, over study follow-up were also associated with increases in emphysema. Baseline ambient O3 was significantly associated with a faster decline in forced expiratory volume in the first second (FEV1).
Long-term exposure to ambient air pollutants, especially O3, was significantly associated with increasing emphysema assessed quantitatively using CT imaging and with worsening lung function.
While air pollutants at historical levels have been associated with cardiovascular and respiratory diseases, it is not known whether exposure to contemporary air pollutant concentrations is associated with progression of emphysema.
To assess the longitudinal association of ambient ozone (O3), fine particulate matter (PM2.5), oxides of nitrogen (NOx), and black carbon exposure with change in percent emphysema assessed via computed tomographic (CT) imaging and lung function.
Design, Setting, and Participants
This cohort study included participants from the Multi-Ethnic Study of Atherosclerosis (MESA) Air and Lung Studies conducted in 6 metropolitan regions of the United States, which included 6814 adults aged 45 to 84 years recruited between July 2000 and August 2002, and an additional 257 participants recruited from February 2005 to May 2007, with follow-up through November 2018.
Residence-specific air pollutant concentrations (O3, PM2.5, NOx, and black carbon) were estimated by validated spatiotemporal models incorporating cohort-specific monitoring, determined from 1999 through the end of follow-up.
Main Outcomes and Measures
Percent emphysema, defined as the percent of lung pixels less than −950 Hounsfield units, was assessed up to 5 times per participant via cardiac CT scan (2000-2007) and equivalent regions on lung CT scans (2010-2018). Spirometry was performed up to 3 times per participant (2004-2018).
Among 7071 study participants (mean [range] age at recruitment, 60 [45-84] years; 3330 [47.1%] were men), 5780 were assigned outdoor residential air pollution concentrations in the year of their baseline examination and during the follow-up period and had at least 1 follow-up CT scan, and 2772 had at least 1 follow-up spirometric assessment, over a median of 10 years. Median percent emphysema was 3% at baseline and increased a mean of 0.58 percentage points per 10 years. Mean ambient concentrations of PM2.5 and NOx, but not O3, decreased substantially during follow-up. Ambient concentrations of O3, PM2.5, NOx, and black carbon at study baseline were significantly associated with greater increases in percent emphysema per 10 years (O3: 0.13 per 3 parts per billion [95% CI, 0.03-0.24]; PM2.5: 0.11 per 2 μg/m3 [95% CI, 0.03-0.19]; NOx: 0.06 per 10 parts per billion [95% CI, 0.01-0.12]; black carbon: 0.10 per 0.2 μg/m3 [95% CI, 0.01-0.18]). Ambient O3 and NOx concentrations, but not PM2.5 concentrations, during follow-up were also significantly associated with greater increases in percent emphysema. Ambient O3 concentrations, but not other pollutants, at baseline and during follow-up were significantly associated with a greater decline in forced expiratory volume in 1 second per 10 years (baseline: 13.41 mL per 3 parts per billion [95% CI, 0.7-26.1]; follow-up: 18.15 mL per 3 parts per billion [95% CI, 1.59-34.71]).
Conclusions and Relevance
In this cohort study conducted between 2000 and 2018 in 6 US metropolitan regions, long-term exposure to ambient air pollutants was significantly associated with increasing emphysema assessed quantitatively using CT imaging and lung function.
Quiz Ref IDChronic lower respiratory disease is the fourth leading cause of death in the United States and third leading cause of death worldwide.1,2 Most of this mortality is due to chronic obstructive pulmonary disease (COPD), defined by persistent airflow limitation, measured via spirometry, and pulmonary emphysema, defined by destruction of the lung parenchyma.3Quiz Ref ID Emphysema can be measured quantitatively on computed tomographic (CT) imaging and can be observed without spirometric COPD.4,5 The percentage of emphysema-like lung on CT scans (“percent emphysema”) is associated with the degree of lung function impairment and symptoms as well as mortality in patients with COPD.6 Percent emphysema is also associated with dyspnea and disease exacerbations independent of lung function among smokers7 and with incident COPD and all-cause mortality among individuals without COPD in the general population.4,8
Ambient air pollution is a major risk factor for poor health worldwide.9 The relationship of air pollution exposures with respiratory morbidity and mortality has been documented in the short-term10 and, by a more limited number of prospective cohorts, in the long-term.11,12 There is elevated risk of death associated with ozone (O3) and fine particulate matter less than 2.5 microns in aerodynamic diameter (PM2.5) exposure in patients with COPD11 and a faster decline in lung function associated with long-term exposure to traffic-related air pollutants in healthy populations.13 No longitudinal studies of air pollution and emphysema progression on CT scan, to our knowledge, have been reported.
This cohort study examined the hypothesis that long-term exposure to air pollutants is associated with progression of percent emphysema shown on CT imaging, as well as decline in lung function, in a well-characterized, multiethnic cohort of adults.
The Multi-Ethnic Study of Atherosclerosis (MESA) recruited participants aged 45 to 84 years without clinical cardiovascular disease from 2000 to 2002 in 6 US metropolitan areas. Investigators recruited white, African American, Hispanic, and Chinese American participants to achieve prespecified recruitment goals by clinical site; race/ethnicity was based on participant self-identification. The MESA study included 4 race/ethnic categories to create a more representative and diverse population than was represented in prior National Heart, Lung, and Blood Institute cohort studies. The MESA Air Pollution Study (MESA Air) recruited additional participants with comparable ages and without cardiovascular disease from 2005 to 2007.14 The MESA Lung Study enrolled a subsample of 3813 MESA participants from 2004 to 2006,15 additional MESA participants from 2010 to 2012, and selected all MESA participants in 2017 to 2018. No selection criteria were related to risk or presence of lung disease. The current analysis includes all MESA participants with valid air pollutant exposure estimates, outcomes, and covariates.
The protocols were approved by the institutional review boards of the collaborating institutions. Written informed consent was obtained from all participants.
Percent emphysema was calculated based on all usable CT scans, including cardiac and full-lung scans. All MESA participants underwent a cardiac CT scan at baseline from 2000 to 2002. Cardiac CT scans were acquired at suspended full inspiration.16 A second CT scan was conducted for 84% of participants in the cohort between 2002 and 2005 and for a subset of participants who did not have a prior scan in 2004 or 2005 and consented to a scan at this examination from 2005 to 2007. Participants in the MESA Lung Study underwent a full-lung CT scan at suspended full inspiration following a standardized protocol from 2010 to 201217 and from 2017 to 2018 with modification to account for scanner evolution.
Trained readers performed percent emphysema measurements at a central reading center, without knowledge of other participant information, using modified Pulmonary Analysis Software Suite for cardiac scans and the updated version, Apollo 1.2 (VIDA Diagnostics), for full-lung scans. For cardiac CT scans, which were performed in duplicate at each examination, the scan with the greater air volume was used except in cases in which 1 scan failed to meet quality control criteria.18
Percent emphysema was defined as the percentage of lung pixels below −950 Hounsfield units (HU), a threshold selected based on pathologic comparisons5 and prognostic significance in this cohort.4 The −950-HU threshold on all scans was adjusted for attenuation of air outside the chest to account for scanner variation.18 Because cardiac CT images approximate 66% of lung volume (from the carina to the lung bases), the upper-third portion of full-lung scans was excluded to follow the same lung region over time. Percent emphysema on cardiac scans was reproducible (scan-rescan intraclass correlation coefficient [ICC], 0.91), comparable between scanners (ICC, 0.94), and valid compared with full-lung CT scans (ICC, 0.93).18
CT technology changed over the 18 years of the study. Therefore, sensitivity analyses were performed using a hidden Markov measure field approach to the measurement of percent emphysema, which reduces interscanner and interprotocol variability by adapting the characteristics of each image under varying CT imaging protocols and differences in intensity distributions due to inspiration level.19 In additional sensitivity analyses, scans with inspiratory levels less than 80% of the participant’s maximum lung volume on scans were excluded and an alternative metric, the HU of the lowest 15th percentile of the lung frequency distribution (PD15), was examined in each participant. A higher PD15 value indicates less emphysema.
Spirometry was conducted in the MESA Lung Study from 2004 to 2006 and repeated twice from 2010 to 2012 and 2017 to 2018 in accordance with American Thoracic Society guidelines and with the same dry-rolling seal spirometers (Occupational Marketing Inc).15 All results were reviewed by 1 investigator. Airflow obstruction was defined at the first measurement of lung function as forced expiratory volume in the first second (FEV1)/forced vital capacity (FVC) ratio less than 0.70.
Quiz Ref IDThe method used to estimate long-term outdoor air pollution concentrations at each participant’s home address has been previously described in detail.20,21 In brief, spatiotemporal exposure models were developed for O3, PM2.5, and oxides of nitrogen (NOx) in each study region based on measurements (1999-2018) from the US Environmental Protection Agency (EPA) Air Quality System and the spatially dense cohort-specific monitoring performed for MESA Air.22 Mean outdoor concentrations of black carbon were estimated only for the years 2006 to 2008 because comparable air quality system data were not available at other times (the EPA’s measurement method for black carbon was not stable over follow-up). City-specific models incorporated a large number of variables covering a wide diversity of geographic features, such as traffic, industrial emissions, population density, and land use. The performance of these models ranged from good to excellent as assessed by the overall cross-validation R2 evaluated at participant residence locations (range, 0.59-0.92).
To account for strong temporal trends in some air pollutants and the knowledge that prior exposure can alter ventilatory function trajectories even after exposure stops (as seen with smoking23), 2 approaches were used to assign long-term mean exposure to each participant: (1) baseline exposure, assessed as pollutant concentration for the year of the baseline examination in 2000 (or 2006-2008 for black carbon) because this exposure measure more closely approximates prior epochs of this exposure and (2) exposure over follow-up, assessed as time-varying mean concentrations between the year of the baseline examination and follow-up time of the repeated outcome measure.
The principal interest of the study was to assess the longitudinal relationship between long-term exposure to air pollutants and progression of percent emphysema and decline in lung function over time (assessed as the interaction between air pollutant and time in the model). As described in more detail in the Supplement, these associations were estimated using linear mixed effects models with random intercept (and slope for the emphysema model).
The primary adjusted model for the percent emphysema and pulmonary function outcomes included cross-sectional terms, including baseline age, sex, race/ethnicity, study region, education, height, weight, body mass index (time-varying), smoking (status, pack-years, cigarettes per day for current smokers, and secondhand smoke exposure [both time-varying]), long-term mean city-specific temperature (time-varying), income, employment outside the home, physical activity, neighborhood socioeconomic status (SES) index,24 and the interaction between SES and study region; longitudinal terms for interactions between these variables and study time (years since baseline); and a transient effects term. While most key covariates were 99% to 100% complete, 5% of participants were missing complete information to calculate pack-years of smoking. We assigned those participants a value of zero and adjusted for an indicator variable reflecting missing pack-years data; we performed sensitivity analyses to confirm that results were consistent between this approach and alternative approaches (excluding missing subjects and assigning the median value for pack-years). Scanner parameters (eg, CT scanner model, milliampere-second adjustment, pixel size) were also included in the percent emphysema model, while the model of lung function additionally controlled for occupational exposure to fumes and gases. Each of the air pollution exposure metrics (at baseline or over follow-up as described above) was investigated separately and included as a main effect and in interaction with time. The associations of exposure to air pollutants with outcomes were compared with the association of pack-years in the same model. In secondary analyses, the regression models were extended by also adjusting for some self-reported transient respiratory conditions (eg, cold, flu, bronchitis, pneumonia) in the lung function models and for high-attenuation areas on CT imaging (percentage of imaged lung volume with attenuation values between −600 and −250 HU) in the emphysema model.
Effect modification of the associations with air pollutants was examined for age, sex, race/ethnicity, obesity, smoking status, SES level, baseline airflow obstruction, study region, and scanner type using 3-way interaction terms between study time, the effect modifier, and exposure concentration. A significant difference between subgroups was indicated by a P value less than .05 for the F test of the 3-way interaction.
To examine whether pollutants acted independently or jointly on lung outcomes, a model with linear combinations was implemented based on estimates from the fully adjusted models that included all the pollutants simultaneously. Concentration-response curves were generated using a thin plate regression spline with 4 degrees of freedom for long-term exposure to the air pollutants assessed at baseline and during follow-up.
The significance threshold was .05 and all tests were 2-sided. Additional sensitivity analyses were conducted by (1) including copollutants to assess independent, combined, and multiplicative associations on the study outcomes, (2) adding an additional random intercept for neighborhood clusters, (3) examining the relationship between pollutants and relative progression of percent emphysema using a log-transformed outcome, (4) using alternative definitions for emphysema, (5) excluding data with low lung volume on CT scan, (6) examining the effect of missing data by excluding participants who underwent only 1 examination, and (7) restricting to participants who had both CT and spirometry measurements. Analyses were performed using SAS version 9.4 (SAS Institute).
This study included 7071 participants (6814 MESA participants and 257 MESA Air study participants). Of 7069 participants with percent emphysema measures, 6860 were assigned outdoor residential air pollution concentrations in the year of their baseline examination and during the follow-up period (eFigure 1 in the Supplement). Of these participants, 5780 (84%) had at least 1 follow-up CT scan over a median (range) of 10 (1-18) years. Air pollution concentrations could be assigned to 3636 of 3813 MESA Lung participants with at least 1 spirometric lung function assessment, 2772 (76%) of whom had a second spirometric assessment at a median of 10 years.
Of the 6860 participants at baseline, 3126 (46%) were lifelong nonsmokers, 2595 (38%) were white, and 4344 (63%) reported at least some college education (Table). Baseline median (interquartile range [IQR]) percent emphysema was 3% (4.5%) and the mean rate of change was 0.58 percentage points (95% CI, 0.38-0.78) over 10 years. The rate of change varied by site and exposure levels (eTable 1 in the Supplement). Regarding lung function, 826 of 3636 participants (22%) had airflow obstruction at baseline (eTable 2 in the Supplement) and the cohort exhibited a mean decline in FEV1 of 309 mL (95% CI, 299-319) and in FVC of 331 mL (95% CI, 317-345) over 10 years.
Air pollution concentrations varied substantially across the study regions over 18 years (Figure 1 and Figure 2). While annual mean concentrations of PM2.5 and NOx declined over follow-up, O3 concentrations did not (Figure 2). Correlations of residential ambient O3 (2000-2018) with residential ambient PM2.5 and NOx within each study region were negative and of moderate degree (r = −0.32 for PM2.5 and −0.42 for NOx) (eTable 3 in the Supplement).
Among the 6860 included participants, greater exposures to O3, PM2.5, and NOx assessed at study baseline and to black carbon averaged through 2006 to 2008 were all significantly associated with an increased progression of percent emphysema over 10 years, per rounded IQR increment (O3: 0.13 per 3 parts per billion [ppb] [95% CI, 0.03-0.24]; PM2.5: 0.11 per 2 μg/m3 [95% CI, 0.03-0.19]; NOx: 0.06 per 10 ppb [95% CI, 0.01 to 0.12]; black carbon: 0.10 per 0.2 μg/m3 [95% CI, 0.01-0.18]).
In the analyses of air pollutant exposure over follow-up, a 3 ppb–higher long-term mean O3 exposure assessed over follow-up was significantly associated with an increased progression of 0.18 percentage points in percent emphysema over 10 years (95% CI, 0.08-0.28) (Figure 3). This increase is equal to the association of 29 pack-years of smoking (each 10 pack-years of smoking was significantly associated with an increased progression of 0.06 percentage points [95% CI, 0.02-0.11] of percent emphysema) or 3 years of aging in this cohort. The association between pollutant concentrations over follow-up and faster progression of percent emphysema was also statistically significant for NOx (increase, 0.12 percentage points per 10 ppb [95% CI, 0.04-0.19]), but not for PM2.5 (increase, −0.04 per 2 μg/m3 [95% CI, −0.15 to 0.08]).
Results for O3, PM2.5, and NOx were not sensitive to accounting for neighborhood clusters, use of log-transformed outcome or alternative definitions of quantitatively assessed emphysema (ie, hidden Markov measure field, PD15), restriction to participants having both percent emphysema and lung function measures, or exclusion of data with low scan lung volume or participants with a single examination (eTables 4 and 5 in the Supplement). Addition of PM2.5, NOx, or black carbon as a covariate similarly did not weaken the associations with exposure to O3 at baseline or over follow-up, nor did addition of O3 affect the association with baseline PM2.5.
Combining the associations of the multiple air pollutants on change in percent emphysema resulted in a greater magnitude of associations than when assessing associations of individual pollutants (Figure 4 and eTable 6 in the Supplement). The concentration-response function between long-term O3 exposure assessed over follow-up and progression of percent emphysema has a sigmoidal appearance, most strongly linear in the area with the most participants (where CIs are narrower) and with suggestion of attenuation of the association at both lower and higher concentrations (where CIs are wider) (Figure 3). While there was no clear-cut evidence for effect modification across the subgroups studied (eTables 7-9 in the Supplement), greater progression of emphysema with O3 assessed over follow-up was observed in the group with airflow limitation at baseline (increase of 0.35 percentage points [95% CI, 0.18-0.51]) compared with the group without airflow limitation (increase, 0.15 percentage points [95% CI, 0.03-0.27]).
Among 3636 participants, 3 ppb in long-term O3 exposure assessed over follow-up was significantly associated with a 18.15-mL greater decline in FEV1 (95% CI, 1.59-34.71) and a 40.19-mL greater decline in FVC (95% CI, 17.88-62.49) over 10 years (eTable 10 and eFigure 3 in the Supplement). These associations were robust when additionally adjusted for baseline respiratory symptoms. The associations were generally not sensitive to inclusion of copollutants. Associations between O3 and FEV1 and FVC were of greater magnitude among current smokers; associations for FEV1 were of greater magnitude among participants with airflow obstruction at baseline (eTable 11 in the Supplement).
There was no significant association between other air pollutants and lung function decline or between O3 and FEV1/FVC ratio (eTable 10 in the Supplement). No indication of combined association was observed for FEV1 and FVC with other pollutants (eFigure 4 in the Supplement). Associations with FEV1 and FVC were statistically significant for exposure to O3 at baseline, but not for the other pollutants (eTable 12 in the Supplement).
Quiz Ref IDHigher residential concentrations of O3, PM2.5, and NOx at study baseline, of black carbon averaged from 2006 to 2008, and of O3 and NOx assessed over study follow-up were significantly associated with greater increases in percent emphysema, assessed via CT imaging, over a median of 10 years. Findings were most robust and of greatest magnitude for O3. Because percent emphysema is related to respiratory symptoms, hospitalizations,7 and mortality even among individuals without airflow obstruction,4,6 these associations in a community-based population demonstrate novel evidence that air pollution contributes to worsening lung health.
Quiz Ref IDDespite existing regulations to prevent short-term excursions of O3 levels, long-term average concentrations of O3, which were associated with changes in lung structure and function, did not decline during the years of observation. Absent new control strategies, these levels are not expected to decline as climate change advances.25 This contrasts with the observed trends in PM2.5 and NOx that highlight the success of past regulatory initiatives to control these 2 pollutants. Because long-term concentrations of O3 at current levels were strongly and consistently associated with both progression of emphysema and decline in lung function in this study, more effective control strategies to reduce O3 concentrations may be needed to protect lung health.
Ground-level O3 is a powerful oxidizing agent and common air pollutant worldwide. Toxicologic studies have shown that O3 causes persistent pulmonary inflammation and structural alterations throughout the lung that may play a role in chronic lower respiratory disease, including emphysema.26 In a mouse model, exposure to O3 for 6 weeks induced a chronic inflammatory process, with increased protease expression, epithelial apoptosis and alveolar enlargement, and destruction, mimicking emphysema.27 In animals, PM2.5 exposure leads to neutrophilic pulmonary inflammation and oxidative stress,28 and in healthy adults, controlled exposure to PM2.5 leads to increased systemic inflammation and endothelial microparticles, including those of lung origin.29 This may be relevant because endothelial apoptosis causes emphysema in animals and endothelial microparticles have been linked to lower lung function and, specifically, percent emphysema in individuals with mild COPD and smokers without COPD.30
The current study found that faster progression of percent emphysema and decline in FEV1 were significantly associated with long-term O3 concentrations among the subset of participants with airflow limitation at baseline. This finding is consistent with prior studies demonstrating worse outcomes related to air pollution in patients with lung disease31 and supports the conclusion in the 2013 EPA Integrated Science Assessment that individuals with lung disease are a high-risk group for O3-related health effects.32 Findings for O3 and loss of lung function were greater among current smokers, perhaps because of additive effects of airway inflammation and oxidative stress.33
Higher exposures to PM2.5 and NOx at baseline and NOx (and no other pollutants) over follow-up were also significantly associated with faster emphysema progression. When all pollutants were considered together, there was evidence that the association for O3 was of the greatest magnitude, which is consistent with other literature,12 and for an additive effect of copollutants, which is consistent with animal studies.33
Unlike O3 and NOx, PM2.5 exposure over follow-up was not significantly associated with emphysema progression. Development of emphysema can be a slow, lifelong process and, once initiated, additional areas of emphysema in the adjacent lung can be produced from biomechanical forces, even in the absence of further exposures, such as cigarette smoke.34 These findings might explain the increased rate of emphysema progression significantly associated with PM2.5 at baseline but not during follow-up. Baseline measures likely reflect exposure at and prior to baseline, particularly for PM2.5 concentrations, which declined substantially over the 10 years of follow-up. In addition, the changes in scanner technology over follow-up as PM2.5 concentrations decreased may have made the follow-up analyses less sensitive for emphysema progression than the baseline analyses.
While O3 exposure was significantly associated with decline in lung function, the other pollutants were not; statistical power for spirometry was smaller than for percent emphysema on CT because the number of observations for spirometry were less than half those for percent emphysema and occurred during the latter period of follow-up when nonozone pollutant concentrations had declined. While air pollutants, including O3, PM2.5, NOx, and black carbon, have been associated consistently with respiratory events,10,12 relatively few studies have reported a longitudinal association between O3, PM2.5, and NOx and decline in lung function in older adults, and findings have not been completely consistent. Single-site studies have reported longitudinal associations between PM2.513 and black carbon35 and a faster decline in lung function, and a larger multi-center study in Europe found PM10 to be associated with a faster lung function decline.36 However, another large European study did not find any longitudinal associations between air pollutants and lung function decline.37
To our knowledge, this is the first longitudinal study to assess the association between long-term exposure to air pollutants and progression of percent emphysema on CT in a large, community-based multiethnic cohort. Other strengths of this study include the large sample size, the long period of follow-up, and the fine-scale assessment of residential-level outdoor air pollution concentrations from cohort-specific monitoring and advanced statistical modeling.
This study has several limitations. First, outdoor air pollution concentrations, especially in the case of O3, may not fully reflect individual air pollution exposures and dose in all microenvironments38; outdoor concentrations do not explain all variation in indoor concentrations, and most individuals spend the majority of their time indoors.39 Second, percent emphysema was measured only in the lower two-thirds of the lung. However, percent emphysema measured in the lower two-thirds of the lung correlates well with full-lung percent emphysema (ICC, 0.93) in this cohort and a cohort of smokers,18 and percent emphysema measured on cardiac scans was associated with dyspnea, spirometric obstruction, and mortality in this cohort.4,8 Third, CT scanners changed over the 18 years of data collection; however, analyses using an advanced image processing approach designed to account for scanner variation yielded similar results, as did sensitivity analyses using various approaches to adjust for stratification by scanner type (eTables 4 and 5 in the Supplement). There is debate about the optimal parameterization of the histogram of the lung attenuation; however, sensitivity analyses using PD15 yielded consistent results. Fourth, no volume correction of CT data was performed,40 given the debate over its utility and the partial lung imaging; however, results were consistent when participants with a difference greater than 20% in lung volume on CT scans were excluded.
In the United States between 2000 and 2018, long-term exposure to ambient air pollutants was significantly associated with increases in emphysema assessed quantitatively via CT imaging and lung function.
Corresponding Author: Joel D. Kaufman, MD, MPH, University of Washington, Box 354695, 4225 Roosevelt Way NE, Ste 100, Seattle, WA 98105 (firstname.lastname@example.org).
Accepted for Publication: June 24, 2019.
Author Contributions: Drs Kaufman and Barr 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. The first authorship roles were shared by Drs Wang, Aaron, and Madrigano and last authorship roles were shared by Drs Kaufman and Barr.
Concept and design: Wang, Hoffman, Kinney, Sheppard, Smith, Vedal, Kaufman, Barr.
Acquisition, analysis, or interpretation of data: Wang, Aaron, Madrigano, Hoffman, Angelini, Yang, Laine, Vetterli, Sampson, Sheppard, Szpiro, Adar, Smith, Kirwa, Lederer, Diez Roux, Vedal, Kaufman, Barr.
Drafting of the manuscript: Wang, Laine, Vedal, Kaufman, Barr.
Critical revision of the manuscript for important intellectual content: All authors.
Statistical analysis: Wang, Aaron, Madrigano, Angelini, Yang, Laine, Sampson, Szpiro, Adar, Kirwa, Vedal, Kaufman, Barr.
Obtained funding: Hoffman, Smith, Kaufman, Barr.
Administrative, technical, or material support: Wang, Hoffman, Angelini, Yang, Laine, Vetterli, Smith, Kirwa, Lederer, Kaufman.
Supervision: Hoffman, Laine, Smith, Lederer, Vedal, Kaufman.
Conflict of Interest Disclosures: Dr Wang reported receiving grants from University of Washington (US EPA RD831697, US EPA RD83479601-01, NIEHS K24ES013195, NIEHS P30ES07033) during the conduct of the study. Dr Aaron reported receiving grants from the National Institutes of Health (NIH) during the conduct of the study, grants from the Stony Wold-Herbert Fund and the Alpha1 Foundation, and personal fees from Lancet Respiratory Medicine outside the submitted work. Dr Madrigano reported receiving grants from National Oceanic and Atmospheric Administration during the conduct of the study and personal fees from the NIH outside the submitted work. Dr Hoffman reported receiving grants from the NIH; is a founder and shareholder of VIDA Diagnostics, a company commercializing lung image analysis software developed, in part, at the University of Iowa; and holds patents for an apparatus for analyzing CT images to determine the presence of pulmonary tissue pathology (US6466687B1), an apparatus for image display and analysis (WO1990016056A1), and a method for multiscale meshing of branching biological structures (US20110093243A1). Dr Yang reported receiving grants from the NIH (NIH R01-HL121270) during the conduct of the study. Dr Sampson reported receiving grants from the EPA (RD831697, RD83479601-01) during the conduct of the study. Dr Sheppard reported receiving grants from the EPA during the conduct of the study and grants from the NIH and personal fees from the Health Effects Institute outside the submitted work. Dr Szpiro reported receiving grants from the EPA and the National Institute of Environmental Health Sciences during the conduct of the study and personal fees from Health Effects Institute and the Electric Power Research Institute outside the submitted work. Dr Smith reported receiving grants from the NIH during the conduct of the study and grants from Quebec Health Research Fund, AstraZeneca, and McGill University Health Center Research Institute outside the submitted work. Dr Lederer reported receiving personal fees from Roche, Sanofi Genzyme, Philips Respironics, Fibrogen, Global Blood Therapeutics, Boehringer-Ingelheim, Veracyte, and Galapagos outside the submitted work; institutional grant support from Fibrogen, Global Blood Therapeutics, and Boehringer-Ingelhim; performing unpaid consulting work for Galecto, Pliant Therapeutics, and Bristol-Myers Squibb; and is now a full-time employee of Regeneron Pharmaceuticals (but was a was a full-time employee of Columbia University during the conduct of the study). Dr Diez-Roux reported receiving grants from the EPA and the NIH during the conduct of the study and Wellcome Trust outside the submitted work. Dr Kaufman reported receiving grants from the EPA and the NIH during the conduct of the study and the US National Institutes for Occupational Safety and Health, the Health Effects Institute, the Kresge Foundation, and the Global Alliance for Clean Cookstoves outside the submitted work. Dr Vedal reported receiving grants from the EPA during the conduct of the study and support for a research chair from AXA Research Fund outside the submitted work. Dr Barr reported receiving grants from the COPD Foundation and the Alpha1 Foundation outside the submitted work and grants from the EPA and NIH during the conduct of the study. No other disclosures were reported.
Funding/Support: This article was developed under a STAR research assistance agreement, No. RD831697 (MESA Air) and RD-83830001 (MESA Air Next Stage), awarded by the US Environmental Protection Agency (EPA) and the University of Washington Center for Clean Air Research (UW CCAR, EPA RD83479601-01). It has not been reviewed by the EPA. The MESA Lung Study was funded by the National Institutes of Health/National Heart, Lung, and Blood Institute grants R01-HL077612, RC1-HL100543, R01-HL093081, R01-HL121270, and R01-HL130605. MESA was funded by National Institutes of Health/National Heart, Lung, and Blood Institute contracts N01-HC-95159, N01-HC-95160, N01-HC-95161, N01-HC-95162, N01-HC-95163, N01-HC-95164, N01-HC-95165, N01-HC-95166, N01-HC-95167, N01-HC-95168, and N01-HC-95169 and by grants UL1-TR-000040 and UL1-TR-001079 from NCRR; the work was also supported by grants K24ES013195, P30ES07033, and R01ES023500 from the National Institute of Environmental Health Sciences.
Role of the Funder/Sponsor: The funders had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.
Disclaimer: The views expressed in this article are solely those of the authors and the EPA does not endorse any products or commercial services mentioned in this publication.
Additional Contributions: The authors thank the other investigators, the staff, and the participants of the MESA study for their valuable contributions. A full list of participating MESA investigators and institutions can be found at http://www.mesa-nhlbi.org.