Ambient Air Pollution and Chronic Bronchitis in a Cohort of U.S. Women
Limited evidence links air pollution exposure to chronic cough and sputum production. Few reports have investigated the association between long-term exposure to air pollution and classically defined chronic bronchitis.
Our objective was to estimate the association between long-term exposure to particulate matter (diameter less than 10 micrometers, PM sub 10; less than 2.5 micrometers, PM sub 2.5), nitrogen dioxide open parenthesis N O sub 2 close parenthesis, and both incident and prevalent chronic bronchitis.
We estimated annual average PM sub 2.5, PM sub 10, and N O sub 2 concentrations using a national land-use regression model with spatial smoothing at home addresses of participants in a prospective nationwide U.S. cohort study of sisters of women with breast cancer. Incident chronic bronchitis and prevalent chronic bronchitis, cough and phlegm, were assessed by questionnaires.
Among 47,357 individuals with complete data, 1,383 had prevalent chronic bronchitis at baseline, and 647 incident cases occurred over 5.7-y average follow-up. No associations with incident chronic bronchitis were observed. Prevalent chronic bronchitis was associated with PM sub 10 [adjusted odds ratio (aOR) per interquartile range (IQR) difference open parenthesis 5.8 micrograms per cubic meter close parenthesis equals 1.07; 95% confidence interval (CI): 1.01, 1.13]. In never-smokers, PM sub 2.5 was associated with prevalent chronic bronchitis (a O R equals 1.18 per IQR difference; 95% CI: 1.04, 1.34), and N O sub 2 was associated with prevalent chronic bronchitis (a O R equals 1.10; 95 percent CI equals 1.01 to 1.20), cough (a O R equals 1.10; 95% CI: 1.05, 1.16), and phlegm (a O R equals 1.07; 95% CI: 1.01, 1.14); interaction p-values (nonsmokers vs. smokers) less than 0.05.
PM sub 10 exposure was related to chronic bronchitis prevalence. Among never-smokers, PM sub 2.5 and N O sub 2 exposure was associated with chronic bronchitis and component symptoms. Results may have policy ramifications for PM sub 10 regulation by providing evidence for respiratory health effects related to long-term PM sub 10 exposure. https://doi.org/10.1289/EHP2199
Chronic bronchitis is a common clinical condition defined by chronic cough and sputum production for at least 3 mo in 2 or more consecutive years (American Thoracic Society 1995). Prevalence estimates in the general population of adults range from 3.5 to 27% (Kim et al. 2011; Martinez et al. 2014; Montes De Oca et al. 2012). This wide range may reflect, in part, variability in case definitions. Chronic bronchitis is a phenotype of chronic obstructive pulmonary disease (COPD) (Kim and Criner 2013). Among persons with COPD, chronic bronchitis portends increased frequency and severity of exacerbations (Burgel et al. 2009; Kim et al. 2011). Among persons without COPD, chronic bronchitis symptoms predict an increased risk of developing COPD, lower health-related quality-of-life scores, and increased risk for all-cause mortality (de Marco et al. 2007; Guerra et al. 2009; Lindberg et al. 2005; Martinez et al. 2014).
Smoking is the primary risk factor for chronic bronchitis, but exposure to ambient air pollution may also contribute (Kim and Criner 2013). The relationship between short-term air pollution exposure and acute respiratory symptoms and hospitalizations is well established (Peacock et al. 2011; Peel et al. 2005; Sunyer 2001), but limited data suggest a relationship between long-term ambient pollution exposure and COPD (Schikowski et al. 2014). There is a paucity of data on the possible relationship between classically defined chronic bronchitis and long-term exposure to the criteria pollutants PM sub 2.5, PM sub 10 (particulate matter less than 2.5 micrometers and less than 10 micrometers in diameter, respectively), and nitrogen dioxide open parenthesis N O sub 2 close parenthesis. The sparse existing data provide inconsistent support for an association between PM sub 10 and chronic cough and phlegm, and between N O sub 2 and chronic cough (Bentayeb et al. 2010b; Cai et al. 2014; Schikowski et al. 2005; Zemp et al. 1999).
To address these relationships in a larger study, using specific outcome definitions and advanced exposure assessments, we investigated the association between residential exposure to PM sub 2.5, PM sub 10, and N O sub 2 and both incident and prevalent chronic bronchitis in a prospective nationwide cohort of more than 50,000 U.S. women participating in the National Institute of Environmental Health Sciences (NIEHS) Sister Study. We estimated exposure at individuals’ residential addresses. Taking advantage of the comprehensive survey, we uniformly classified cases of chronic bronchitis using the classical clinical definition.
The NIEHS Sister Study is a longitudinal cohort study of U.S. women with a sister diagnosed with breast cancer, but no personal breast cancer diagnosis at time of baseline interview (n equals 50,884). Women were enrolled between August 2003 and March 2009, and completed a baseline computer-assisted telephone survey. Follow-up telephone surveys were performed every 2 to 3 y. We analyzed data through the second follow-up survey (data release 4, data available through August 2014). Baseline and follow-up surveys queried participants on a wide range of health diagnoses and symptoms.
Of the 50,884 women participating in the NIEHS Sister Study, 1,234 (2.4%) were excluded for missing exposure data due to residential locations outside the modeling region or addresses that could not be geocoded (Figure 1). After excluding those missing baseline data on cough and phlegm, 47,357 individuals remained for analysis of prevalent outcomes. Of the 45,955 participants without chronic bronchitis symptoms at baseline, 6,111 (12.3%) were missing data on cough or phlegm for at least one of the two follow-up questionnaires, leaving 39,844 individuals for analysis of incident outcomes.
The Institutional Review Boards of the University of Washington and the NIEHS approved this study; all participants provided written informed consent.
Chronic bronchitis was defined according to the classical symptom-based definition of chronic cough productive of phlegm for at least 3 mo out of a year for a minimum of 2 consecutive years (American Thoracic Society 1995). Participants were asked about the presence of cough and phlegm independently, and the duration of each symptom was specified using questions derived from the British Medical Research Council adult respiratory symptom standardized questionnaire. Women with cough and phlegm symptoms, both present for at least 3 mo per year out of the previous 2 y, were considered to have chronic bronchitis. Prevalent chronic bronchitis was determined by meeting symptom-based criteria at the baseline questionnaire. In a sensitivity analysis, we included history of physician diagnosis of chronic bronchitis in the case definition. Incident chronic bronchitis was defined by satisfying the case definition at either the second follow-up survey, or both the first and second follow-ups among participants who did not have chronic bronchitis at baseline. Participants whose symptoms did not persist from first through second follow-up were not considered cases.
Secondary outcomes were chronic cough (3 or more months of cough for at least 2 consecutive years, regardless of phlegm production), chronic phlegm (3 or more months of phlegm production for at least 2 consecutive years, regardless of cough), and chronic cough or phlegm. Both prevalent chronic cough and chronic phlegm were defined by being present at baseline.
Ambient Air Pollution Exposure Assessment
Air pollution exposure was estimated using annual average PM sub 2.5, PM sub 10, and N O sub 2 levels at each participant’s current primary residence. Home addresses of participants were geocoded using ArcGIS (version 10; Esri). We estimated long-term exposure using year 2000 annual mean concentration levels for all pollutants. Measurements of PM sub 2.5, PM sub 10, and N O sub 2 concentrations from monitors using federal reference methods were obtained from the U.S. Environmental Protection Agency (EPA) Air Quality System database. After excluding locations with only seasonal coverage or large amounts of missing data, the observations were aggregated into annual averages. The annual averages were used to fit a universal kriging regression model for predicting at points within the contiguous United States. The models for PM sub 2.5 (Sampson et al. 2013) and N O sub 2 (Young et al. 2016) have been previously described in detail, and the model for PM sub 10 was fit in the same manner as the PM sub 2.5 model. Partial least squares, a dimension reduction technique, was used to select linear combinations of land use, roadway proximity, and other geographic covariates. The N O sub 2 prediction model additionally incorporated satellite data (Young et al. 2016). Spatial smoothing was included via an exponential covariance function. This model therefore incorporated land-use regression and spatial smoothing of values observed in the monitoring network. Model performance was evaluated using 10-fold cross-validation and for the year 2000. The cross-validated R super 2 was 0.85 for N O sub 2, 0.53 for PM sub 10, and 0.77 for PM sub 2.5. Exposure modeling was limited to the continental United States; participants from Alaska, Hawaii, and Puerto Rico were excluded (n equals 1,234).
To estimate the association between outcomes and pollutant exposures, we used multivariable logistic regression. Covariates were selected a priori based on plausible relationships and review of existing literature. Potential confounders, measured at baseline, were age (continuous), ethnicity (non-Hispanic white, non-Hispanic black, Hispanic, other), body mass index (continuous), education (high school or less, some college, associate or technical degree, bachelor’s degree, graduate degree), household income (continuous), occupational exposure to dust (ever/never) or vapors/fumes (ever/never), smoking status (never, former, current), tobacco pack-years (continuous), and years of secondhand smoke exposure since age 19 (continuous). After exclusion of individuals missing any of these covariates, 44,158 individuals were available for the analysis of prevalent outcomes. For incident outcomes, among those without chronic bronchitis at baseline, 38,006 individuals were available after excluding missing covariates.
A model adjusted for age alone was first performed followed by a fully adjusted model including all a priori identified covariates. PM sub 2.5, PM sub 10, and N O sub 2 were modeled separately. Given potential bias by length of follow-up time for the incident chronic bronchitis outcome, adjustment was made for duration of follow-up time (from baseline to second follow-up survey) using restricted cubic splines with four knots (Dinse and Lagakos 1983). We performed analyses stratified by baseline smoking status (ever/never) and tested for interaction using product terms for smoking status and pollutant exposure. A two-pollutant model was performed by including remaining copollutants in the fully adjusted model. In all instances, a p-value of less than 0.05 was considered significant.
Sensitivity analyses were prespecified and performed on the following subgroups in independent analyses: a) prevalent chronic bronchitis that included either symptom-based criteria or report of physician diagnosis; b) excluding baseline asthmatics [defined by history of physician diagnosis, recent (within 12 mo) asthma medication use, and self-reported current asthma], given clinical overlap between asthma and chronic bronchitis; and c) participants who lived at least 10 y at their primary residence. In the asthma sensitivity analysis, asthmatics who reported either current smoking or history of greater than or equal to 10 pack-years at baseline were not excluded given possibility of smoking-related symptoms leading to asthma misdiagnosis. Given concern that seasonal variation may affect results, an additional sensitivity analysis was performed by adjusting for season at the time of baseline and follow-up questionnaires.
Statistical analysis was performed using the statistical program Stata (version 11; StataCorp).
Participants were, on average, 55.4 y old at baseline open parenthesis SD equals 8.9 close parenthesis, 84.8% were white, 52.6% had a bachelor’s degree or higher level of education, 56.4% had never smoked cigarettes, and only 8% were current smokers. The proportion of black and Hispanic participants increased across tertiles of PM sub 2.5, PM sub 10, and N O sub 2 (Table 1). The distributions of occupational exposures, smoking history, and cumulative tobacco smoke exposure did not vary materially by ambient air pollution exposure levels (Table 1).
|Characteristic||PM sub 2.5 (micrograms per cubic meter)||N O sub 2 (ppb)||PM sub 10 (micrograms per cubic meter)|
|Age (years)||55.9 plus or minus 8.9||55.4 plus or minus 9||54.9 plus or minus 8.9||55.7 plus or minus 8.8||55.3 plus or minus 8.9||55.3 plus or minus 9||55.7 plus or minus 8.7||55.3 plus or minus 8.9||55.3 plus or minus 9|
|BMI||27.5 plus or minus 6||27.9 plus or minus 6.3||28.3 plus or minus 6.5||27.8 plus or minus 6.1||27.8 plus or minus 6.3||27.9 plus or minus 6.5||27.5 plus or minus 6.1||27.9 plus or minus 6.3||28.2 plus or minus 6.5|
|HS or less||14.7||15.1||14.8||17.6||14.8||12.3||15.1||15.1||14.5|
|Household income (USD)||41,218 plus or minus 25,849||43,570 plus or minus 26,914||43,743 plus or minus 28,293||39,995 plus or minus 24,817||43,064 plus or minus 26,339||45,472 plus or minus 29,533||42,473 plus or minus 25,999||43,104 plus or minus 26,892||42,953 plus or minus 28,242|
|Occupational Exposures (%)|
|Smoking status (%)|
|Pack-years among ever-smokers||14.5 plus or minus 15.3||14.7 plus or minus 15.2||14.8 plus or minus 15.5||15.1 plus or minus 15.6||14.3 plus or minus 15.1||14.5 plus or minus 15.2||14.8 plus or minus 15.5||14.5 plus or minus 15.2||14.6 plus or minus 15.2|
|Packs per day among current smokers||0.7 plus or minus 0.4||0.7 plus or minus 0.4||0.7 plus or minus 0.5||0.7 plus or minus 0.5||0.7 plus or minus 0.4||0.6 plus or minus 0.4||0.7 plus or minus 0.5||0.7 plus or minus 0.4||0.7 plus or minus 0.4|
|Adult secondhand smoke (years)||10.9 plus or minus 12.9||11.3 plus or minus 13.3||11.4 plus or minus 13.2||11.5 plus or minus 13.2||11.1 plus or minus 13.1||11.1 plus or minus 13.1||11.3 plus or minus 13.1||11.1 plus or minus 13.1||11.3 plus or minus 13.3|
|Number of years at primary residence||12.8 plus or minus 10.9||13.5 plus or minus 11.3||13.6 plus or minus 11.5||12.7 plus or minus 11||12.7 plus or minus 10.7||14.6 plus or minus 11.8||13.4 plus or minus 11||13.2 plus or minus 11.2||13.3 plus or minus 11.5|
|Lived at primary residence greater than or equal to 10 years (%)||50.7||53.5||53.9||50.1||50.8||57.1||53.6||52.2||52.4|
|Asthma at baseline (%)||5.7||5.6||5.9||5||6||6.3||5.4||5.9||6|
|Physician diagnosis of COPD (%)||1.5||1.5||1.5||1.6||1.6||1.4||1.5||1.5||1.5|
|Physician diagnosis of chronic bronchitis (%)||7.4||8||8.4||8||8.1||7.8||7.5||8||8.4|
The mean follow-up time was 5.7 y from enrollment to the second follow-up survey. During the follow-up period, there were 638 incident cases of chronic bronchitis, giving an estimated incidence rate of 2.8 cases per 1,000 person-years. At baseline, 1,351 (3.1%) women met symptom-based criteria for chronic bronchitis, whereas 4,698 (10.6%) participants reported ever having had a physician diagnosis of chronic bronchitis. Prevalent chronic cough was reported by 3,749 (8.5%) and chronic phlegm by 2,776 (6.3%) participants at baseline.
The median estimated exposure concentrations were 12.4 micrograms per cubic meteropen parenthesis interquartile range open square bracket IQR close square bracket close parenthesis equals 4.4 micrograms per cubic meter for PM sub 2.5, 2.16 micrograms per cubic meter (5.8 micrograms per cubic meter) for PM sub 10, and 11.7 parts per billion (7.3 parts per billion) for N O sub 2. The results of the age-adjusted and fully adjusted regression analyses are presented in Table 2. No statistically significant associations were found between incident chronic bronchitis and any of the air pollution exposures. Limiting the incidence analysis to long-term residents (greater than 10 years) did not appreciably alter the effect estimates (Table S1).
|Exposure and outcome||Cases||Age adjusted||Fully adjusted|
|OR (95% CI)||p-Value||OR (95% CI)||p-Value|
|PM sub 2.5|
|Incident chronic bronchitis||638||0.94 (0.84, 1.05)||0.256||0.94 (0.83, 1.06)||0.289|
|Prevalent (at baseline)|
|Chronic bronchitis||1,351||1.04 (0.97, 1.13)||0.276||1.04 (0.96, 1.13)||0.318|
|Chronic cough||3,749||1.03 (0.98, 1.08)||0.213||1.04 (0.99, 1.10)||0.103|
|Chronic phlegm||2,776||1.07 (1.02, 1.13)||0.010||1.04 (0.98, 1.10)||0.213|
|Chronic cough or phlegm||5,271||1.05 (1.01, 1.10)||0.015||1.04 (1.00, 1.09)||0.067|
|N O sub 2|
|Incident chronic bronchitis||638||0.95 (0.87, 1.03)||0.198||1.00 (0.92, 1.09)||0.974|
|Prevalent (at baseline)|
|Chronic bronchitis||1,351||1.00 (0.95, 1.06)||0.923||1.05 (0.99, 1.11)||0.136|
|Chronic cough||3,749||1.02 (0.99, 1.06)||0.215||1.06 (1.02, 1.10)||0.002|
|Chronic phlegm||2,776||1.01 (0.97, 1.05)||0.730||1.02 (0.98, 1.07)||0.266|
|Chronic cough or phlegm||5,271||1.02 (0.99, 1.05)||0.199||1.04 (1.01, 1.08)||0.008|
|PM sub 10|
|Incident chronic bronchitis||638||0.92 (0.85, 1.01)||0.066||0.98 (0.90, 1.08)||0.745|
|Prevalent (at baseline)|
|Chronic bronchitis||1,351||1.06 (1.01, 1.12)||0.027||1.07 (1.01, 1.13)||0.019|
|Chronic cough||3,749||1.04 (1.00, 1.07)||0.045||1.04 (1.00, 1.08)||0.030|
|Chronic phlegm||2,776||1.07 (1.03, 1.12)||less than 0.001||1.07 (1.02, 1.11)||0.002|
|Chronic cough or phlegm||5,271||1.05 (1.01, 1.08)||0.001||1.05 (1.02, 1.08)||0.003|
For prevalent chronic bronchitis, a statistically significant positive association was seen with PM sub 10 [odds ratio (OR) per IQR increase in PM sub 10 equals 1.07; 95% confidence interval (CI): 1.01, 1.13] (Table 2). Similar magnitudes of association with prevalent chronic bronchitis were seen for N O sub 2 (O R equals 1.05; 95% CI: 0.99, 1.11) and PM sub 2.5 (O R equals 1.04; 95% CI: 0.96, 1.13), but were not statistically significant. PM sub 10 was also statistically significantly associated with chronic cough (O R equals 1.04; 95% CI: 1.00, 1.08), chronic phlegm (O R equals 1.07; 95% CI: 1.02, 1.11), and chronic cough or phlegm (O R equals 1.05; 95% CI: 1.02, 1.08); coadjustment for PM sub 2.5 did not alter these effect estimates (Table S2). Adjustment of the PM sub 10 model for N O sub 2 resulted in a general attenuation of associations between prevalent symptoms and PM sub 10. This attenuation is likely due in part to the strong correlation between N O sub 2 and PM sub 10 (Pearson’s r: 0.59).
N O sub 2 showed a significant positive association with chronic cough (OR = 1.06; 95% CI: 1.02, 1.10) and chronic cough or phlegm (O R equals 1.04; 95% CI: 1.01, 1.08). In the N O sub 2 model, ORs were robust to coadjustment for PM sub 2.5 (Table S2). Coadjustment for PM sub 10 in the N O sub 2 model resulted in a loss of precision for the association between N O sub 2 and chronic cough or phlegm, and an overall decrease in size of effect estimates across all prevalent outcomes. The significant association with chronic cough was preserved (O R equals 1.05; 95% CI: 1.01, 1.10). ORs for all pollutants and outcomes were generally very similar between age-adjusted and fully adjusted models.
For prevalent chronic bronchitis, sensitivity analyses incorporating additional case requirements into the classical symptom-based definition showed comparable effect estimates in association with PM sub 10 (Table 3). For example, PM sub 10 was significantly associated with prevalent chronic bronchitis defined either by symptoms or including participant-reported physician diagnosis (OR per IQR increase equals 1.06; 95% CI: 1.02, 1.09). In an analysis of prevalent chronic bronchitis excluding baseline asthmatics, the effect estimate was similar but less precise, commensurate with the smaller sample size (OR for IQR increase in PM sub 10 equals 1.06; 95% CI: 0.99, 1.13). Similarly, exclusion of the 47% of participants who lived at their residence less than 10 y largely preserved the estimated association, but with loss of precision reflecting the smaller sample size (OR per IQR increase in PM sub 10 equals 1.07; 95% CI: 0.99, 1.16). Comparable sensitivity analyses involving PM sub 2.5 or N O sub 2 and prevalent chronic bronchitis yielded ORs that were similar in magnitude and direction to the primary models (Table 3). Effect estimates for all three pollutants were essentially unchanged by seasonal adjustment (Table S3).
|Exposure and case definitions (primary and sensitivity analyses)||na||Cases||Adjusted OR (95% CI)||p-Value|
|PM sub 2.5|
|Prevalent chronic bronchitis (primary case definition)||44,158||1,351||1.04 (0.96, 1.13)||0.318|
|Including physician diagnosis||44,099||4,698||1.04 (0.99, 1.09)||0.104|
|Excluding asthma at baseline||41,488||1,104||0.99 (0.90, 1.08)||0.798|
|Excluding those living at residence less than 10 years||23,273||720||0.97 (0.87, 1.09)||0.643|
|N O sub 2|
|Prevalent chronic bronchitis (primary case definition)||44,158||1,351||1.05 (0.99, 1.11)||0.136|
|Including physician diagnosis||44,099||4,698||1.02 (0.99, 1.06)||0.191|
|Excluding asthma at baseline||41,488||1,104||1.02 (0.96, 1.09)||0.492|
|Excluding those living at residence less than 10 years||23,273||720||1.03 (0.95, 1.11)||0.444|
|PM sub 10|
|Prevalent chronic bronchitis (primary case definition)||44,158||1,351||1.07 (1.01, 1.13)||0.019|
|Including physician diagnosis||44,099||4,698||1.06 (1.02, 1.09)||0.001|
|Excluding asthma at baseline||41,488||1,104||1.06 (0.99, 1.13)||0.077|
|Excluding those living at residence less than 10 years||23,273||720||1.07 (0.99, 1.16)||0.093|
In smoking-stratified analyses, we found evidence for stronger associations between all three air pollutants and prevalent outcomes in never-smokers (Table 4). PM sub 2.5 was strongly associated with prevalent chronic bronchitis among never-smokers (OR per IQR difference equals 1.18; 95% CI: 1.04, 1.34), and the difference by smoking status was statistically significant (p sub interaction equals 0.013). Similarly with N O sub 2, in never-smokers, significant associations were seen for all four prevalent outcomes: chronic bronchitis (O R equals 1.10; 95% CI: 1.01, 1.20), chronic cough (O R equals 1.10; 95% CI: 1.05, 1.16), chronic phlegm (O R equals 1.07; 95% CI: 1.01, 1.14), and chronic cough or phlegm (O R equals 1.09; 95% CI: 1.04, 1.13), and the differences by smoking status were statistically significant for both cough (p sub interaction equals 0.020), phlegm (p sub interaction equals 0.017), and cough or phlegm (p sub interaction equals 0.004). Corresponding ORs for PM sub 2.5 and N O sub 2 were close to null for ever-smokers. For PM sub 10, results did not differ significantly by smoking status, although the same pattern of stronger associations in never-smokers was seen (Table 4).
|Exposure and outcome||Never-smoker||Ever-smoker|
|Cases||OR (95% CI)||p-Value||Cases||OR (95% CI)||p-Value||P sub Interaction|
|PM sub 2.5|
|Incident chronic bronchitis||271||0.92 (0.77, 1.10)||0.382||367||0.95 (0.81, 1.11)||0.516||0.815|
|Prevalent (at baseline)|
|Chronic bronchitis||580||1.18 (1.04, 1.34)||0.011||771||0.96 (0.86, 1.07)||0.427||0.013|
|Chronic cough||1,802||1.07 (1.00, 1.15)||0.053||1,947||1.02 (0.95, 1.10)||0.545||0.345|
|Chronic phlegm||1,362||1.08 (1.00, 1.17)||0.059||1,414||1.00 (0.92, 1.09)||0.956||0.189|
|Chronic cough or phlegm||2,632||1.06 (0.99, 1.12)||0.077||2,639||1.03 (0.97, 1.10)||0.292||0.638|
|N O sub 2|
|Incident chronic bronchitis||271||1.03 (0.91, 1.17)||0.609||367||0.97 (0.86, 1.10)||0.660||0.498|
|Prevalent (at baseline)|
|Chronic bronchitis||580||1.10 (1.01, 1.20)||0.029||771||1.00 (0.92, 1.08)||0.955||0.097|
|Chronic cough||1,802||1.10 (1.05, 1.16)||less than 0.001||1,947||1.01 (0.96, 1.06)||0.642||0.020|
|Chronic phlegm||1,362||1.07 (1.01, 1.14)||0.014||1,414||0.97 (0.92, 1.03)||0.359||0.017|
|Chronic cough or phlegm||2,632||1.09 (1.04, 1.13)||less than 0.001||2,639||0.99 (0.95, 1.04)||0.806||0.004|
|PM sub 10|
|Incident chronic bronchitis||271||1.04 (0.91, 1.18)||0.587||367||0.95 (0.83, 1.08)||0.458||0.359|
|Prevalent (at baseline)|
|Chronic bronchitis||580||1.09 (1.00, 1.20)||0.055||771||1.06 (0.98, 1.14)||0.131||0.597|
|Chronic cough||1,802||1.07 (1.01, 1.12)||0.015||1,947||1.02 (0.97, 1.08)||0.391||0.260|
|Chronic phlegm||1,362||1.10 (1.04, 1.17)||0.002||1,414||1.04 (0.98, 1.10)||0.179||0.189|
|Chronic cough or phlegm||2,632||1.08 (1.03, 1.13)||less than 0.001||2,639||1.02 (0.98, 1.07)||0.352||0.081|
To our knowledge, this is the largest study to investigate the association between classically defined chronic bronchitis and long-term ambient air pollution exposure using a validated national exposure model. We did not find an association between incident chronic bronchitis and any of the three air pollution measures. However, exposure to higher concentrations of PM sub 10 was significantly associated with all prevalent outcomes: chronic bronchitis, chronic cough, chronic phlegm, and chronic cough or phlegm. These findings were robust to coadjustment for PM sub 2.5 in a two-pollutant model (Table S2). We also found N O sub 2 exposure was significantly associated with chronic cough and chronic cough or phlegm. To the best of our knowledge, no other study has shown an association between PM sub 10 and classically defined chronic bronchitis. These findings provide evidence that long-term ambient air pollution exposure, particularly PM sub 10, is a risk factor for chronic bronchitis and the chronic respiratory symptoms of cough and phlegm that define it.
Incident chronic bronchitis should be superior to prevalent chronic bronchitis for making causal inference regarding observed associations with air pollution. However, the relatively short follow-up duration (mean: 5.7 y) limited our power to detect an association between ambient air pollutants and incident chronic bronchitis. With the much larger number of cases of prevalent conditions, we had substantially higher power than for the incident analyses. One smaller study of nonsmoking Seventh Day Adventists in California has shown an association between incident chronic bronchitis and long-term exposure to PM sub 2.5; however, levels were in excess of 20 micrograms per cubic meter, a concentration almost double that observed in our study (Abbey et al. 1995).
Comparison to previous studies is limited due to substantial variability in defining chronic bronchitis and exposure estimation methods. The observed incidence rate of 2.5 cases per 1,000 person-years and prevalence of 2.9% are at the low end of the range reported in the literature (Cai et al. 2014; Cerveri et al. 2001; Huchon et al. 2002; Kim et al. 2011; Sobradillo et al. 1999). However, our study population was more than half nonsmoking women, and our estimates are in agreement with study populations with similar demographics (Montes De Oca et al. 2012; Sunyer et al. 2006). National prevalence and incidence figures for chronic bronchitis are lacking because they rely on physician diagnosis rather than the classical symptom-based diagnostic criteria (American Lung Association 2013). Including participant-reported physician diagnosis greatly increases the prevalence of chronic bronchitis in this study and likely elsewhere (Schikowski et al. 2005).
Our study provides evidence that PM sub 10 exposure is a risk factor for chronic bronchitis, while the existing literature suggests associations between PM sub 10 and various respiratory symptoms. A large cross-sectional study in Switzerland found an association between increased prevalence of chronic cough and phlegm with PM sub 10 exposure among never-smokers (Zemp et al. 1999). The European Study of Cohorts for Air Pollution Effects (ESCAPE) meta-analysis of five European cohorts similarly showed an association between PM sub 10 and prevalent chronic phlegm, but not chronic bronchitis, in never-smokers (Cai et al. 2014). A French study of elderly adults demonstrated increased prevalence of chronic cough associated with PM sub 10 exposure (Bentayeb et al. 2010a). Furthermore, in the Swiss Cohort Study on Air Pollution and Lung and Heart Diseases in Adults (SAPALDIA) cohort, decline in PM sub 10 over time was associated with a reduction in chronic cough and phlegm (Schindler et al. 2009). Our study suggests long-term PM sub 10 exposure is associated with prevalent chronic bronchitis, the distinct clinical entity, as well as the associated symptoms that define it.
In contrast to PM sub 2.5, which deposits within the distal alveoli, the preferential deposition of coarse particles within the conducting airways of the tracheobronchial tree provides biologic plausibility for the association between PM sub 10 and chronic bronchitis (Carvalho et al. 2011). Chronic airway epithelial inflammation and mucus metaplasia are the pathologic bases of chronic bronchitis (Kim and Criner 2013). Chronic bronchitis is associated with narrowing and mucus plugging of the nonalveolated conducting airways (Matsuba and Thurlbeck 1973). PM has been frequently implicated in triggering pro-inflammatory cascades within airway epithelial cells (Øvrevik et al. 2015). Certain PM sub 10 components, including transition metals and endotoxins, have been shown to drive airway inflammation and mucus hypersecretion via upregulation of transcription factors and generation of reactive oxygen species and oxidative stress (Longphre et al. 2000; Øvrevik et al. 2015).
For N O sub 2, we saw no associations with prevalent chronic bronchitis or chronic phlegm in the main analyses; however, there were positive associations with both chronic cough and chronic cough or phlegm. These associations were robust to coadjustment for PM sub 2.5. Adjustment for PM sub 10, in an attempt to isolate the effect of N O sub 2, was limited by the strong correlation between these copollutants. Previous studies have shown inconsistent associations between respiratory symptoms and N O sub 2. The ESCAPE meta-analysis showed no significant association between N O sub 2 and chronic bronchitis, cough, or phlegm (Cai et al. 2014). The German Study on the Influence of Air Pollution on Lung, Inflammation and Aging (SALIA) cohort analysis showed association with cough, but not chronic bronchitis, while the Swiss SAPALDIA study showed no association overall (Schikowski et al. 2005; Zemp et al. 1999). In both the SALIA and SAPALDIA studies, conducted in the 1990s, the annual mean N O sub 2 levels were more than double those in our study.
Among never-smokers, associations between prevalent symptoms and exposures were stronger for all pollutants. If we limit our interpretation of the stratified analyses to those with statistically significant interactions, we find stronger evidence for associations between both N O sub 2 and PM sub 2.5 and our prevalent outcomes among never-smokers. The finding of stronger associations with these two pollutants and outcomes in never-smokers is consistent with some previous literature findings. In the Swiss SAPALDIA study, N O sub 2 was related to chronic bronchitis only among nonsmokers. For PM sub 2.5, the ESCAPE meta-analysis found a positive association between chronic cough and PM sub 2.5 only among never-smokers (Cai et al. 2014). N O sub 2 was associated with chronic cough in a study of about 4,700 women in Germany, 74% of whom were never-smokers (Schikowski et al. 2005).
The reason for stronger associations with N O sub 2 in never-smokers in our study and others is unclear. Perhaps the airways of smokers are more tolerant to the irritant effects of ambient N O sub 2 than nonsmokers who, as a result, are dissuaded from smoking because of greater sensitivity to these effects. Alternatively, the effects of compounds in cigarette smoke might swamp the effects of ambient N O sub 2 exposure. Stronger associations with PM sub 2.5 in never-smokers could reflect, in part, the high dose of fine particles inhaled by smokers. The biologic effects of cigarette smoke may overwhelm the effects from long-term, low-level ambient air pollution and thus mask any association. For PM sub 10, we saw clear associations with chronic bronchitis in all subjects, possibly related to PM sub 10 favoring deposition in the conducting airways that are responsible for producing bronchitic symptoms. This association is apparent in the aggregate sample. In contrast, the associations between chronic bronchitis and PM sub 2.5 and N O sub 2 were only apparent in never-smokers. Theses pollutants have distribution patterns that tend to bypass the conducting airways for deposition and adsorption in the distal alveoli.
Outcome misclassification was reduced by using the symptom-defined definition of chronic bronchitis. However, the symptom-based questionnaire still has limitations; due to recall bias, it likely results in inclusion of cases with recent, but not necessarily chronic, symptoms. Overlap with asthma remains possible given the clinical similarities of these conditions. Sensitivity analyses excluding subjects with a physician diagnosis of asthma and active asthma symptoms at baseline attenuated the association with all pollutants. However, the exclusion of these participants with overlapping chronic bronchitis symptoms may have eliminated true cases of chronic bronchitis and thus reduced power.
The sensitivity analysis including individuals reporting doctor diagnosis of chronic bronchitis showed preserved associations with increased precision. Self-report of doctor diagnosis is likely to include more individuals with symptom duration shorter than the 2-y minimum required for the chronic bronchitis definition (i.e., may include participants who have received a diagnosis of acute bronchitis in the past). Given that chronic bronchitis is defined by duration of symptoms, directly asking subjects questions on cough and phlegm is preferable to asking about physician diagnosis, which requires accurate reporting by both parties. However, the association between PM sub 10 and prevalent chronic bronchitis remains robust, even with the less strict definition, suggesting that presence of symptoms is driving the relationship, rather than the specific duration of symptoms.
The objection has been raised that chronic bronchitis prevalence as reported on questionnaires may reflect recent symptoms and that prevalence and/or severity might then vary by season. Therefore, we undertook a sensitivity analysis adjusting for season of both baseline and follow-up questionnaire (Table S3). No change was observed in the effect estimates for outcomes associated with PM sub 2.5 or N O sub 2. The associations between prevalent chronic bronchitis and chronic cough and PM sub 10 were no longer statistically significant after adjusting for season, but the effect estimates remained largely unchanged and in the anticipated direction. Season of questionnaire administration does not seem to contribute significant bias in reporting of chronic bronchitis.
Our air pollution exposure estimates are based on a validated national model using land-use regression and spatial smoothing to capture within- and between-region air pollution variability and minimize exposure misclassification. This model is a considerable improvement over road proximity, regional fixed-monitor averages, and simple land-use regression models employed in prior chronic bronchitis investigations (Keller et al. 2014; Young et al. 2016). In addition, seasonal bias in the exposure should be mitigated by using annual averages and a chronic outcome whose case definition dictates that symptoms must span a minimum of 2 consecutive years. Air pollution estimates used year 2000 annual averages, predating baseline enrollment for all participants. While concentrations of criteria pollutants are declining nationally, spatial differences of annual average pollution concentrations account for the majority of variability in PM sub 2.5 measurement and were relatively stable across the study period (Kim et al. 2017). However, it is acknowledged that variability in the decline in pollution levels may contribute to exposure misclassification, and the resulting biases are difficult to predict. It is plausible that our observed lack of association for incident outcomes and positive effects for prevalent symptoms could be related to variable change in pollution levels; e.g., if the most polluted regions experienced more dramatic declines in levels than the cleaner areas.
Additional limitations exist. This study is limited to women, and the findings may not be broadly applicable to men. Outdoor ambient pollutant concentrations may not reflect the indoor exposures. Exposure measurement error owing to our use of residential addresses to characterize exposure owing both to subjects’ residential mobility, time spent away from home or indoors, and spatiotemporal trends is a limitation both of this study and epidemiological studies of air pollution health effects in general. Exposure measurement of this nature error is generally expected to bias associations toward the null rather than producing false positive associations.
In 2006, the EPA revoked the National Ambient Air Quality Standard for annual PM sub 10 due to insufficient data on health risks associated with long-term exposure to PM sub 10 as opposed to the finer PM sub 2.5 fraction (U.S. EPA 2006). The preceding long-term PM sub 10 exposure standard was an annual average of 50 micrograms per cubic meter, roughly double the mean concentration experienced by participants in this study. This study provides evidence that chronic respiratory health effects occur with long-term exposure to PM sub 10 at levels below the previous national standards. These results add to a limited body of evidence relating morbidity to long-term PM sub 10 exposure and consequently may have policy implications both nationally and globally.
The authors thank A. Gassett and C. Sack for data management and support with this analysis, as well as the staff at Social and Scientific Systems, Inc. and Westat, Inc. for overall study support and management of Sister Study data acquisition. The authors also thank the University of Washington Center for Clean Air Research (UW CCAR) and the Multi-Ethnic Study of Atherosclerosis and Air Pollution (MESA Air) for providing the infrastructure and development of the exposure models. MESA Air is conducted and supported by the U.S. EPA through grant RD831697 to the UW. This publication has not been formally reviewed by the U.S. EPA. The views expressed in this document are solely those of the authors and do not necessarily reflect those of the agency. U.S. EPA does not endorse any products or commercial services mentioned in this publication. The authors also thank the thousands of women who participate in the Sister Study.
This work was supported by the University of Washington Pulmonary and Critical Care Medicine Training Grant (T32 HL007287), the Intramural Research Program of the National Institutes of Health (NIH), and the National Institute of Environmental Health Sciences (NIEHS) (Z01 ES044005 and ES043012). Exposure modeling was supported by the U.S. EPA (RD831697 and R834796), and Biostatistics, Epidemiologic, and Bioinformatic Training in Environmental Health Training Grant was from NIEHS (T32ES015459).
Abbey DE, Ostro BE, Petersen F, Burchette RJ. 1995. Chronic respiratory symptoms associated with estimated long-term ambient concentrations of fine particulates less than 2.5 microns in aerodynamic diameter (PM2.5) and other air pollutants. J Expo Anal Environ Epidemiol 5(2):137–159, PMID: 7492903. Medline, Google Scholar
- American Lung Association.2013. Trends in COPD (Chronic Bronchitis and Emphysema): Morbidity and Mortality. http://www.lung.org/assets/documents/research/copd-trend-report.pdf [accessed
9 May 2017]. Google Scholar
- American Thoracic Society. 1995. Standards for the diagnosis and care of patients with chronic obstructive pulmonary disease. Am J Respir Crit Care Med
152(5 Part 2):S77–121, PMID: 7582322. Medline, Google Scholar Bentayeb M, Helmer C, Raherison C, Dartigues JF, Tessier JF, Annesi-Maesano I. 2010. Bronchitis-like symptoms and proximity air pollution in French elderly. Respir Med 104(6):880–888, PMID: 20129767, doi: 10.1016/j.rmed.2010.01.004. Crossref, Medline, Google Scholar Burgel PR, Nesme-Meyer P, Chanez P, Caillaud D, Carré P, Perez T, 2009. Cough and sputum production are associated with frequent exacerbations and hospitalizations in COPD subjects. Chest 135(4):975–982, PMID: 19017866, doi: 10.1378/chest.08-2062. Crossref, Medline, Google Scholar Cai Y, Schikowski T, Adam M, Buschka A, Carsin AE, Jacquemin B, 2014. Cross-sectional associations between air pollution and chronic bronchitis: an ESCAPE meta-analysis across five cohorts. Thorax 69(11):1005–1014, PMID: 25112730, doi: 10.1136/thoraxjnl-2013-204352. Crossref, Medline, Google Scholar Carvalho TC, Peters JI, Williams RO. 2011. Influence of particle size on regional lung deposition–what evidence is there?Int J Pharm 406(1–2):1–10, PMID: 21232585, doi: 10.1016/j.ijpharm.2010.12.040. Crossref, Medline, Google Scholar Cerveri I, Accordini S, Verlato G, Corsico A, Zoia MC, Casali L, 2001. Variations in the prevalence across countries of chronic bronchitis and smoking habits in young adults. Eur Respir J 18(1):85–92, PMID: 11510810. Crossref, Medline, Google Scholar de Marco R, Accordini S, Cerveri I, Corsico A, Antó JM, Künzli N, 2007. Incidence of chronic obstructive pulmonary disease in a cohort of young adults according to the presence of chronic cough and phlegm. Am J Respir Crit Care Med 175(1):32–39, PMID: 17008642, doi: 10.1164/rccm.200603-381OC. Crossref, Medline, Google Scholar Dinse GE, Lagakos SW. 1983. Regression analysis of tumour prevalence data. J R Stat Soc Ser C Appl Stat 32(3):236, doi: 10.2307/2347946. Crossref, Google Scholar Guerra S, Sherrill DL, Venker C, Ceccato CM, Halonen M, Martinez FD. 2009. Chronic bronchitis before age 50 years predicts incident airflow limitation and mortality risk. Thorax 64(10):894–900, PMID: 19581277, doi: 10.1136/thx.2008.110619. Crossref, Medline, Google Scholar Huchon GJ, Vergnenègre A, Neukirch F, Brami G, Roche N, Preux PM. 2002. Chronic bronchitis among French adults: high prevalence and underdiagnosis. Eur Respir J 20(4):806–812, PMID: 12412668, doi: 10.1183/09031936.02.00042002. Crossref, Medline, Google Scholar Keller JP, Olives C, Kim SY, Sheppard L, Sampson PD, Szpiro AA, 2014. A unified spatiotemporal modeling approach for predicting concentrations of multiple air pollutants in the Multi-Ethnic Study of Atherosclerosis and Air Pollution. Environ. Health Perspect 123(4):301–309, PMID: 25398188, doi: 10.1289/ehp.1408145. Link, Google Scholar Kim SY, Olives C, Sheppard L, Sampson PD, Larson TV, Keller JP, 2017. Historical prediction modeling approach for estimating long-term concentrations of PM sub 2.5 in cohort studies before the 1999 implementation of widespread monitoring. Environ Health Perspect 125(1):38–46, PMID: 27340825, doi: 10.1289/EHP131. Link, Google Scholar Kim V, Criner GJ. 2013. Chronic bronchitis and chronic obstructive pulmonary disease. Am J Respir Crit Care Med 187(3):228–237, PMID: 23204254, doi: 10.1164/rccm.201210-1843CI. Crossref, Medline, Google Scholar Kim V, Han MK, Vance GB, Make BJ, Newell JD, Hokanson JE, 2011. The chronic bronchitic phenotype of COPD: an analysis of the COPDGene Study. Chest 140(3):626–633, PMID: 21474571, doi: 10.1378/chest.10-2948. Crossref, Medline, Google Scholar Lindberg A, Jonsson AC, Rönmark E, Lundgren R, Larsson LG, Lundbäck B. 2005. Ten-year cumulative incidence of COPD and risk factors for incident disease in a symptomatic cohort. Chest 127(5):1544–1552, PMID: 15888826, doi: 10.1378/chest.127.5.1544. Crossref, Medline, Google Scholar Longphre M, Li D, Li J, Matovinovic E, Gallup M, Samet JM, 2000. Lung mucin production is stimulated by the air pollutant residual oil fly ash. Toxicol Appl Pharmacol 162(2):86–92, PMID: 10637131, doi: 10.1006/taap.1999.8838. Crossref, Medline, Google Scholar Martinez CH, Kim V, Chen Y, Kazerooni EA, Murray S, Criner GJ, 2014. The clinical impact of non-obstructive chronic bronchitis in current and former smokers. Respir Med 108(3):491–499, PMID: 24280543, doi: 10.1016/j.rmed.2013.11.003. Crossref, Medline, Google Scholar Matsuba K, Thurlbeck WM. 1973. Disease of the small airways in chronic bronchitis. Am Rev Respir Dis 107(4):552–558, PMID: 4697664, doi: 10.1164/arrd.1922.214.171.1242. Crossref, Medline, Google Scholar Montes de Oca M, Halbert RJ, Lopez MV, Perez-Padilla R, Tálamo C, Moreno D, 2012. The chronic bronchitis phenotype in subjects with and without COPD: the PLATINO study. Eur Respir J 40(1):28–36, PMID: 22282547, doi: 10.1183/09031936.00141611. Crossref, Medline, Google Scholar Øvrevik J, Refsnes M, Låg M, Holme J, Schwarze PE. 2015. Activation of proinflammatory responses in cells of the airway mucosa by particulate matter: oxidant- and non-oxidant-mediated triggering mechanisms. Biomolecules 5(3):1399–1440, PMID: 26147224, doi: 10.3390/biom5031399. Crossref, Medline, Google Scholar Peacock JL, Anderson HR, Bremner SA, Marston L, Seemungal TA, Strachan DP, 2011. Outdoor air pollution and respiratory health in patients with COPD. Thorax 66(7):591–596, PMID: 21459856, doi: 10.1136/thx.2010.155358. Crossref, Medline, Google Scholar Peel JL, Tolbert PE, Klein M, Metzger KB, Flanders WD, Todd K, 2005. Ambient air pollution and respiratory emergency department visits. Epidemiology 16(2):164–174, PMID: 15703530, doi: 10.1097/01.ede.0000152905.42113.db. Crossref, Medline, Google Scholar Sampson PD, Richards M, Szpiro AA, Bergen S, Sheppard L, Larson TV, 2013. A regionalized national universal kriging model using Partial Least Squares regression for estimating annual PM sub 2.5 concentrations in epidemiology. Atmos Environ (1994) 75:383–392, PMID: 24015108, doi: 10.1016/j.atmosenv.2013.04.015. Crossref, Medline, Google Scholar Schikowski T, Mills IC, Anderson HR, Cohen A, Hansell A, Kauffmann F, 2014. Ambient air pollution: a cause of COPD?Eur Respir J 43(1):250–263, PMID: 23471349, doi: 10.1183/09031936.00100112. Crossref, Medline, Google Scholar Schikowski T, Sugiri D, Ranft U, Gehring U, Heinrich J, Wichmann HE, 2005. Long-term air pollution exposure and living close to busy roads are associated with COPD in women. Respir Res 6:152, PMID: 16372913, doi: 10.1186/1465-9921-6-152. Crossref, Medline, Google Scholar Schindler C, Keidel D, Gerbase MW, Zemp E, Bettschart R, Brändli O, 2009. Improvements in PM10 exposure and reduced rates of respiratory symptoms in a cohort of Swiss adults (SAPALDIA). Am J Respir Crit Care Med 179(7):579–587, PMID: 19151198, doi: 10.1164/rccm.200803-388OC. Crossref, Medline, Google Scholar Sobradillo V, Miravitlles M, Jiménez CA, Gabriel R, Viejo JL, Masa JF, 1999. [Epidemiological study of chronic obstructive pulmonary disease in Spain (IBERPOC): prevalence of chronic respiratory symptoms and airflow limitation]. Arch. Bronconeumol 35(4):159–166, PMID: 10330536. Crossref, Medline, Google Scholar Sunyer J. 2001. Urban air pollution and chronic obstructive pulmonary disease: a review. Eur Respir J 17(5):1024–1033, PMID: 11488305, doi: 10.1183/09031936.01.17510240. Crossref, Medline, Google Scholar Sunyer J, Jarvis D, Gotschi T, Garcia-Esteban R, Jacquemin B, Aguilera I, 2006. Chronic bronchitis and urban air pollution in an international study. Occup Environ Med 63(12):836–843, PMID: 16847030, doi: 10.1136/oem.2006.027995. Crossref, Medline, Google Scholar
- U.S. EPA (Environmental Protection Agency). 2006. National Ambient Air Quality Standards for Particulate Matter. Final Rule. https://www.gpo.gov/fdsys/pkg/FR-2006-10-17/pdf/06-8477.pdf [accessed
29 December 2017]. Google Scholar Young MT, Bechle MJ, Sampson PD, Szpiro AA, Marshall JD, Sheppard L, 2016. Satellite-based N O sub 2 and model validation in a national prediction model based on universal kriging and land-use regression. Environ Sci Technol 50(7):3686–3694, PMID: 26927327, doi: 10.1021/acs.est.5b05099. Crossref, Medline, Google Scholar Zemp E, Elsasser S, Schindler C, Künzli N, Perruchoud AP, Domenighetti G, 1999. Long-term ambient air pollution and respiratory symptoms in adults (SAPALDIA study). The SAPALDIA Team. Am J Respir Crit Care Med 159(4 Pt 1):1257–1266, PMID: 10194174, doi: 10.1164/ajrccm.159.4.9807052. Crossref, Medline, Google Scholar