Article Text

Original research
Ambient air pollution exposure and chronic bronchitis in the Lifelines cohort
  1. Dany Doiron1,
  2. Jean Bourbeau1,
  3. Kees de Hoogh2,3,
  4. Anna L Hansell4,5
  1. 1Respiratory Epidemiology and Clinical Research Unit, Research Institute of the McGill University, Centre for Outcomes Research and Evaluation, Research Institute of the McGill University Health Centre, Montréal, Québec, Canada
  2. 2Swiss Tropical and Public Health Institute, Basel, Basel-Stadt, Switzerland
  3. 3University of Basel, Basel, Basel-Stadt, Switzerland
  4. 4Centre for Environmental Health and Sustainability, University of Leicester, Leicester, UK
  5. 5MRC Centre for Environment and Health, Department of Epidemiology and Biostatistics, School of Public Health, Imperial College London, London, UK
  1. Correspondence to Dr Dany Doiron, Respiratory Epidemiology and Clinical Research Unit, Research Institute of the McGill University, Centre for Outcomes Research and Evaluation, Research Institute of the McGill University Health Centre, Montreal, Canada; dany.doiron{at}


Background Few large studies have assessed the relationship of long-term ambient air pollution exposure with the prevalence and incidence of symptoms of chronic bronchitis and cough.

Methods We leveraged Lifelines cohort data on 132 595 (baseline) and 65 009 (second assessment) participants linked to ambient air pollution estimates. Logistic regression models adjusted for sex, age, educational attainment, body mass index, smoking status, pack-years smoking and environmental tobacco smoke at home were used to assess associations of air pollution with prevalence and incidence of chronic bronchitis (winter cough and sputum almost daily for ≥3 months/year), chronic cough (winter cough almost daily for ≥3 months/year) and prevalence of cough and sputum symptoms, irrespective of duration.

Results Associations were seen for all pollutants for prevalent cough or sputum symptoms. However, for prevalent and incident chronic bronchitis, statistically significant associations were seen for nitrogen dioxide (NO2) and black carbon (BC) but not for fine particulate matter (PM2.5). For prevalent chronic bronchitis, associations with NO2 showed OR: 1.05 (95% CI: 1.02 to 1.08) and with BC OR: 1.06 (95% CI: 1.03 to 1.09) expressed per IQR; corresponding results for incident chronic bronchitis were NO2 OR: 1.07 (95% CI: 1.02 to 1.13) and BC OR: 1.07 (95% CI: 1.02 to 1.13). In subgroup analyses, slightly stronger associations were observed among women, never smokers and younger individuals.

Conclusion This is the largest analysis to date to examine cross-sectional and longitudinal associations between ambient air pollution and chronic bronchitis. NO2 and BC air pollution was associated with increased odds of prevalent and incident chronic bronchitis.

  • COPD epidemiology

This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See:

Statistics from

Key messages

What is the key question?

  • Is long-term ambient air pollution exposure associated with prevalence and incidence of chronic bronchitis and/or cough or sputum symptoms?

What is the bottom line?

  • Nitrogen dioxide and black carbon air pollutant concentrations in the northern Netherlands were clearly associated with increased odds of prevalent and incident chronic bronchitis, with higher effect sizes seen among women, never smokers and younger individuals.

Why read on?

  • This is the largest study so far on the relationship between ambient air pollution exposure and chronic bronchitis or related symptoms; even at relatively low concentrations compared with urban locations within Europe and internationally, statistically significant associations of long-term air pollution exposure on prevalent and incident chronic bronchitis were observed.


Chronic bronchitis (CB) is characterised by persistent cough and sputum production and is a common phenotype of chronic obstructive pulmonary disease (COPD). Reported prevalence of CB ranges from 3.4% to 22% in the adult population, and from 14% to 74% among COPD patients.1 In the presence of chronic airflow obstruction, CB has been related to excess decline in lung function2 3 and to a greater prevalence of other respiratory symptoms and exacerbations.2 4–8 Although there is still debate on the relevance of CB in individuals with normal spirometry, recent studies have suggested it to be of worse prognosis. In individuals with normal lung function, CB can carry a risk of excess lung function decline and COPD development.9–11 Chronic cough and sputum symptoms have been demonstrated to be associated with worse health-related quality of life not only in adults with COPD but also in those without airflow obstruction,4 7 8 12 and a recent population-based study showed that CB symptoms have a larger burden on quality of life than chronic airflow obstruction or asthma.8 CB symptoms may not be as benign as previously thought when also considering their association with a higher risk of respiratory-related mortality.2

The role of outdoor air pollution in respiratory mortality, COPD and asthma has been studied widely,13–16 but its role in the prevalence and incidence of CB needs to be elucidated. Ambient air pollution has been shown to increase hospital admissions for respiratory disease17 18 and has been associated with chronic respiratory outcomes such as the development of asthma19 and reduced lung function growth in children,20 as well as lung function decline and increased COPD prevalence and incidence in adults.21 22 In 2015, ambient fine particulate matter and ozone air pollution caused an estimated 1.1 million deaths from COPD worldwide.23 In the Western world, smoking is still in itself the most important indicator of risk of CB, but other exposures including ambient air pollution must be investigated. Few recent studies have examined the long-term effects of outdoor air pollution exposure on CB prevalence and incidence. The largest to date, involving 47 357 US women followed between 2003 and 2014, found statistically significant associations between particulate matter exposure and prevalence, but not incidence, of CB.24 In a 2014 study of five heterogeneous European cohorts (n~10 000), no associations were found for CB, but small increases in prevalence of sputum symptoms were reported among never smokers25—it did not assess incidence. The objective of this study was to assess the relationship of long-term ambient air pollution exposure with the prevalence and incidence of CB as well as cough or sputum symptoms in a large population-based study in the Netherlands with linked air pollution estimates.

Materials and methods

Study population

We analysed baseline (2006–2013) and second assessment (2014–2017) data from adults aged 18 to 93 years from the three northernmost, largely rural, provinces of Netherlands (Friesland, Groningen and Drenthe) enrolled in the Lifelines cohort study. General practitioners invited all their patients aged 25–50 years unless they met one or more exclusion criteria: severe psychiatric/physical illness, limited life expectancy (<5 years) or insufficient knowledge of the Dutch language to complete a study questionnaire. Recruited individuals were asked to invite family members to participate in the study. Inhabitants of the Friesland, Groningen and Drenthe provinces could also register themselves via the Lifelines website.26 Details on Lifelines population sampling and data collection methodology have been outlined elsewhere.26 27 Participants provided informed consent at recruitment for broad use of their data by local and international investigators. All necessary approvals were obtained from the Lifelines ethics and scientific review boards for this project.

Outcome assessment

The gold standard definition of CB, as established by the British Medical Research Council (MRC) in 1965, is cough and sputum production for 3 months of the year for at least two consecutive years.28 CB assessment items included in the Lifelines respiratory health questionnaire deviated slightly from this classical definition in that it: (a) specified ‘winter’ as the time period for recall of respiratory symptoms at baseline and, (b) did not allow assessment of whether symptoms had been present for ≥2 consecutive years. For our study, prevalent (baseline) CB cases were therefore defined as individuals who reported ‘coughing up phlegm when getting up, during the daytime, or at night in the winter almost daily for at least 3 months a year’. Incident CB was defined as individuals ‘coughing up mucus every day for a period of 3 months a year’, among those not reporting daily cough or daily cough with phlegm ≥3 months a year at baseline (table 1). As evidence suggests that chronic cough and chronic sputum production may represent different pathophysiological processes,3 29 we considered ‘cough without phlegm production almost daily for at least 3 months a year’ as a separate outcome at baseline and second assessment. To assess the effect of air pollution exposure with cough and phlegm symptoms, we also considered ‘usual cough’ and ‘usual sputum’ during the winter, but not specified as occurring 3 months a year, as prevalent outcomes in separate analyses (table 1). Lifelines CB questionnaire assessment items and recoding procedures used to define outcomes are included as online supplemental table S1.

Table 1

Chronic bronchitis and bronchitis-like symptom definitions

Air pollution exposure

We used annual mean concentrations of particulate matter with a diameter of less than 2.5 µm (PM2.5), nitrogen dioxide (NO2) and black carbon (BC) developed in the context of the Effects of Low-Level Air Pollution: A Study in Europe (ELAPSE) project.30 31 Land use regression models in ELAPSE were based on 2010 monitoring data from the AirBase dataset for PM2.5 and NO2 and from the European Study of Cohorts for Air Pollution Effects (ESCAPE) project for BC. Predictor variables in the ELAPSE PM2.5 model included satellite-derived PM2.5 data, long range chemical transport model data for PM2.5, altitude, sum of all roads length, natural land areas, ports and residential areas. NO2 predictors included NO2 chemical transport model data, sum of all roads and major roads lengths, natural land area, ports and residential areas. Because BC is a major component of PM2.5, PM2.5 satellite and chemical transport model data were included in the BC models in addition to major roads, all roads, residential areas, urban green spaces and east–west trend. Across western Europe, final models explained 72%, 59% and 54% of the spatial variation in measured PM2.5, NO2 and BC concentrations, respectively.30 When compared with ESCAPE project measured PM2.5 and NO2 concentrations for the Netherlands, the ELAPSE models performed very well for NO2 (R2=75.9) but less so for PM2.5 (R2=12.6).30 Residential addresses of Lifelines participants at recruitment were geocoded, which allowed linkage of ELAPSE air pollution estimates.

Statistical analyses

We conducted descriptive statistics for sociodemographic factors, health outcomes and air pollution. Multivariable logistic regression models were used to explore associations of exposure to PM2.5, NO2 and BC with prevalent CB at baseline. Logistic regression models were also used to estimate associations between air pollution exposure and incident CB at second assessment. Sociodemographic and behavioural risk factors to include as potential confounders were selected a priori through literature search. Prevalent CB analyses were adjusted for sex, age (continuous), educational attainment (low: junior secondary/lower vocational or less, medium: senior secondary/secondary vocational, high: higher vocational/university), body mass index (continuous), smoking status (never, past, current), pack-years smoking and exposure to environmental tobacco smoke at home (yes/no). We also explored prevalent (baseline) and incident (second assessment) chronic cough, as well as prevalent usual cough and usual sputum reported at baseline, using logistics regression models adjusted for the same covariates as CB models. Incident analyses were additionally adjusted for follow-up time (continuous) due to potential bias by length of follow-up. In order to explore the temporal stability of CB, we examined the proportion of individuals with persistent CB—that is, present at baseline and at second assessment—and that of individuals with ‘resolved’ CB—that is, present at baseline but not at second assessment.

Sensitivity analyses were conducted for the following subgroups: women, individuals aged 50 years and older, individuals without post-secondary education, never smokers and asthmatics. For both prevalent and incident outcomes, we conducted two-pollutant models which included NO2 and PM2.5 in order to differentiate the roles of co-occurring air pollutants on respiratory outcomes. To assess the impact of a change of residence and potential for exposure misclassification on incident CB results, we excluded individuals who reported moving houses between baseline and follow-up assessments in sensitivity analyses. Finally, all analyses were conducted on participants with complete data for respiratory outcomes, air pollution exposure and confounders. Statistical analyses were conducted using the R statistical package, V.3.5.2 (20 December 2018).32


Study population

Figure 1 shows Lifelines study population, exclusions and missing data. Of the 152 141 participants recruited at baseline, 4029 had missing data for either air pollution or CB and 15 243 had missing data for covariates, resulting in a sample of 132 595 individuals for prevalent CB analyses, that is, the ‘prevalent cohort’. After excluding 13 761 individuals reporting either CB or daily cough at least 3 months a year at baseline, and 53 825 with either no follow-up data or missing data for CB at second asessement, 65 009 participants were included in incident CB analyses, that is, the ‘incident cohort’. Descriptive analyses (table 2) showed that participants on average were in their mid-40s, three out of five were women and one-third had higher education. About half of individuals were never smokers, while one-third and one-fifth were former and current smokers, respectively. The proportion of never smokers at baseline was slightly higher for women (49.3%) relative to men (44.8%). CB prevalence was 6.1% at baseline while CB incidence was 4.5% at second assessment. Chronic cough prevalence and incidence were 7.4% and 6.1%, respectively. Mean follow-up between baseline and second assessment was 46.5 months (SD=14.02 months), with a range of 13–122 months. Significant differences were found for all variables between the baseline population used in our analyses (n=132 595) and individuals with missing data for respiratory outcomes, air pollution exposure or confounders (n=19 546) (online supplemental table S2). Mean±SD average ambient concentrations of PM2.5, NO2 and BC at baseline residential locations of Lifelines participants were 14.8±1.02 µg/m3, 21.5±4.99 µg/m3 and 1.22±0.23 10−5 m−1, respectively. PM2.5 was moderately correlated with NO2 (r=0.57) and BC (r=0.48), while a high correlation was seen between NO2 and BC (0.89) (table 3).

Figure 1

Lifelines study population, missing data and exclusions.

Table 2

Lifelines study population characteristics

Table 3

Baseline air pollution concentrations descriptive statistics

Prevalence of CB and respiratory symptoms

Logistic regression results for prevalent outcomes are presented in table 4. Higher odds of CB were seen for exposure to NO2 and BC air pollution but not for PM2.5, after adjusting for confounding factors. Statistically significant positive associations were also seen for prevalent chronic cough and exposure to ambient BC, as well as for odds of usual cough and usual sputum and exposure to each ambient BC, NO2 and PM2.5.

Table 4

Prevalent chronic bronchitis and cough and sputum symptoms at baseline per IQR increase in each pollutant

In subgroup analyses (table 5), women and never smokers showed similar patterns of associations with slightly higher effect sizes compared with full analyses. Among never smokers, higher odds of CB were seen for NO2 and marginally higher for BC exposure. Associations between each pollutant and prevalent usual cough and usual sputum were also stronger among asthmatics. Restricting association to individuals with no post-secondary education yielded similar associations for each pollutant and each outcome as seen in the overall population. Finally, individuals aged 50 years and older generally showed reduced effect sizes compared with overall population. No associations were found between any air pollutants and prevalent respiratory outcomes in this subgroup except for the following: prevalent usual cough was associated with exposure to NO2 and BC and chronic cough was associated with BC.

Table 5

Subgroup analyses for prevalent chronic bronchitis and cough and sputum symptoms at baseline, per IQR increase in each pollutant

Incidence of CB and respiratory symptoms

Table 6 shows results of incident CB and chronic cough at second assessment. Consistent with results for CB prevalence, statistically significant positive associations were seen between incident CB and exposure to ambient NO2 and BC, but not for PM2.5. Incident chronic cough was not associated with air pollution exposure. Per IQR increase in NO2 and BC, ORs for incident CB were higher among women and never smokers relative to full analyses. When limiting analyses to never smokers, associations between NO2 and BC exposure and incident chronic cough became larger and statistically significant. Approximately the same proportion of individuals reported persistent CB (ie, present at baseline and second assessment; 1.4%) as did individuals reporting ‘resolved’ CB (ie, present at baseline but not at second assessment; 1.6%) (online supplemental table S3).

Table 6

Incident chronic bronchitis and chronic cough at second assessment, per IQR increase in each pollutant

Sensitivity analyses

Of the 65 009 participants with complete data for incident analyses, 9130 individuals (14.1%) reported moving residences between baseline and follow-up assessments. When limiting incident CB analyses to non-movers (n=55 602), positive, though not statistically significant, associations were found for NO2 air pollution while associations for BC remained statistically significant (online supplemental table S4). In two-pollutant models, associations between PM2.5 exposure and prevalent respiratory symptoms were attenuated to the null after adjustment for NO2 exposure (online supplemental table S5). Conversely, associations between NO2 and CB prevalence and incidence were slightly stronger when adjusting for PM2.5, while NO2 associations with prevalent and incident respiratory symptoms were robust to PM2.5 adjustment.


To our knowledge, this study is the largest to date to explore the associations of long-term exposure to outdoor air pollution with CB. Exposure to ambient NO2 and BC air pollution, but not PM2.5, were associated with increased odds of both prevalent and incident CB. Prevalent but not incident chronic cough was also associated with BC air pollution exposure. All three pollutants explored also showed positive associations with prevalence of usual cough and usual sputum. In both prevalence and incidence analyses, we demonstrated slightly stronger associations of respiratory outcomes with exposure to NO2 and BC air pollution among women, never smokers and younger individuals. Stronger associations were also seen in asthmatics for exposure to air pollution and prevalent cough and sputum symptoms.

Although comparisons with previous studies are limited given different definitions and questionnaires used for CB and symptoms assessment, as well as variations in air pollution exposure estimation methods, evidence reported in past literature are in line with associations seen in the current study. Our results confirm stronger associations previously reported among never smokers. A large study involving 47 357 US women followed from 2003 to 2014, found significant associations between particulate matter with a diameter <10 µm (PM10) and prevalence (OR: 1.07, 95% CI: 1.01 to 1.13, per 5.8 µg/m3 increase), but not incidence, of CB.24 When limiting analyses to never smokers, this study showed that prevalent (ie, baseline) CB was also associated with exposure to PM2.5 (OR: 1.18, 95% CI: 1.04 to 1.34, per 4.4 µg/m3 increase) and NO2 (OR: 1.10, 95% CI: 1.01 to 1.20, per 7.3 ppb increase). While no associations were found between CB and ambient air pollution in a 2014 cross-sectional analysis of five European cohorts (n>10 000) in the ESCAPE project, increases in sputum symptoms were reported among never smokers for PM10 (OR: 1.32, 95% CI: 1.02 to 1.71, per 10 µg/m3 increase) and particulate matter with a diameter between 2.5 µm and 10 µm or PMcoarse (OR: 1.31, 95% CI: 1.05 to 1.64, per 5 µg/m3 increase). A study of 9651 individuals enrolled in the Swiss Study on Air Pollution and Lung Disease in Adults (SAPALDIA) found positive associations between a 10 µg/m3 increase in annual mean concentrations of PM10 and reported prevalence of chronic sputum (OR: 1.35, 95% CI: 1.11 to 1.65) and chronic cough or sputum (OR: 1.27, 95% CI: 1.08 to 1.50) among never smokers, but not among current or former smokers.33 Stronger associations between air pollution and other respiratory outcomes such as COPD and lung function among never smokers have also been reported.22 34 Smoking could saturate the effects of other exposures on respiratory outcomes; additional impacts of air pollution might therefore be smaller and more difficult to detect among past and current tobacco smokers.

As in our study, higher susceptibility of women to the impacts of ambient air pollution exposure on bronchitis symptoms has been reported. In a European Community Respiratory Health Survey analysis involving 6924 subjects in 10 European countries, associations between chronic sputum and ambient NO2 concentrations were seen in women (OR: 1.76, 95% CI: 1.04 to 2.98, per 30 µg/m3 increase), but not in men.35 In line with air pollution research, evidence shows that compared with men with similar levels of chronic airflow obstruction, women may be more susceptible to smoking effects on lung function and respiratory symptoms than men.36 37 Biological factors such as differences in airway reactivity and metabolism of inhaled particles between men and women have been proposed for such sex differences.38 39 A better estimation of air pollution exposure due to more time spent at home could also explain larger associations among women. Although overall indoor air quality was not directly measured, this potential source of bias is largely mitigated by the use of environmental tobacco smoke exposure as an adjustment factor in all regression models.

Evidence for effect modification by age on the relationship between air pollution exposure and CB is limited. In our study, subgroup analyses found few significant associations in individuals 50 years and older. This might partly explain the unexpected null associations found by Cai et al25; three of the five cohorts included in the ESCAPE Study had a mean age >50 years at baseline assessment. A longitudinal study of 1421 individuals aged 21–80 at recruitment also found CB to be an independent risk factors for incident COPD and mortality risk among subjects younger than 50 years, but not among subjects 50 years and older.40 Additional research exploring age-stratified effects of air pollution exposure on CB is therefore warranted.

The ELAPSE air pollution models used to estimate air pollution exposure in our study were robust explaining 58% and 66% of respectively NO2 and PM2.5 concentration in hold-out validation.30 When validated against ESCAPE monitoring data, the ELAPSE model explained 76% of the spatial variability in the measured NO2 at the ESCAPE sites, but only 12.6% for PM2.5. The low correlation between ELAPSE predicted and ESCAPE measured PM2.5 concentrations could be due to the small variation in PM2.5 concentrations (SD=1.8 µg/m3) in the Netherlands and the small number of PM2.5 ESCAPE sites (n=34) used for validation. Participants’ exposure was estimated based on 2010 measurement data although recruitment and follow-up CB and respiratory symptom measurements ranged from 2006 to 2017. A study by de Hoogh et al30 investigated the stability of the 2010 NO2 Western European LUR model by comparing the 2010 NO2 predictions with predictions from NO2 LUR models for 2005 and 2000 at 454 randomly selected points in the Netherlands.30 They found a good correlation between the 2010 predictions and the 2005 (88%) and the 2000 (82%) NO2 predictions. Eeftens et al developed NO2 LUR models for 2007 and 1999–2000 and applying these to NO2 measurements from respectively 1999–2000 and 2007.41 They found that the 2007 and 1999–2000 models were able to explain 77% and 81% of the spatial variability in the 1999–2000 and 2007 NO2 measurements, respectively. Both these studies give confidence that NO2 models for the Netherlands remain spatially stable over time during the study period.

This study has important public health implications. Annual mean NO2 exposure in our study was 21.5 µg/m3, about half the 40 µg/m3 threshold established by the WHO to protect public health. Even at relatively low concentrations compared with urban locations within Europe and internationally, statistically significant associations of long-term NO2 air pollution exposure on CB and cough and sputum symptoms were observed. The stronger associations of CB prevalence and incidence per IQR increase of BC compared with PM2.5 might also indicate that combustion-related particulate substances that dominate BC are more harmful to lung health compared with overall PM2.5. Fine particulate matter air pollution sources include combustion as well as mechanical (eg, brake and tire wear) traffic emissions, residential biomass burning, suspended dust and industry. BC is an important component of PM2.5 emitted largely from gas and diesel engines and power plants, whereas NO2 comes mainly from traffic-related fossil fuel combustion. While we did not have source-specific information for estimated air pollutants, positive associations seen for NO2 and BC reported in our study suggest that policies targeting reduction of air pollution from combustion emissions, as well as traffic abatement measures might help maximise health benefits.42 Positive associations were found for PM2.5 exposure and prevalent usual cough and usual sputum, but not for CB. Non-significant associations for PM2.5 exposure might in part be explained by the lower performance (ie, lower R2) of the ELAPSE PM2.5 model, which could be related to more important exposure misclassification for this pollutant. However, the differences may not all be due to exposure misclassification bias. As noted above, it may be that local traffic-related PM, which NO2 and BC provide a proxy measure of, are the most important exposure in relation to generation of respiratory symptoms. This hypothesis is supported by results from two-pollutant models in our study; PM2.5 exposure and prevalent respiratory symptoms were attenuated to the null after adjustment for NO2 exposure (online supplemental table S5), suggesting local traffic exposure component of PM (as captured by NO2) is an important driver of findings. Previous studies have reported associations between CB and cough and sputum symptoms with larger fractions of particulate matter (PM10 and PMcoarse), but not fine particulates.24 25 However, given coarse particulate matter was not modelled in the ELAPSE project, its impact on respiratory outcomes could not be investigated.

Varying case definitions for CB have been used in epidemiological studies, leading to differences in reported prevalence of the disease. Relative to the classical CB definition established by the MRC, our study used a less stringent definition by not specifying persistence of symptoms for 2 years or more. Prevalent CB in our study was also limited to winter symptoms. In the Latin American Project for Research in Pulmonary Obstruction Study, the prevalence of CB doubled when comparing the classical definition to a less stringent definition.7 Nonetheless, very similar clinical phenotypes have been found between individuals characterised with classically defined CB, compared with individuals reporting CB with less stringent definitions. For example, Kim and colleagues showed that regardless of the definitions used, individual reporting cough and sputum had a higher level of dyspnoea, more respiratory symptoms and a larger number of total and severe exacerbation-like events.43 The use of less stringent CB definitions is also helpful in detecting chronic respiratory symptoms related to long-term air pollution exposure that would not be identified using the gold standard CB definition.

Most studies exploring associations between CB and air pollution have been cross-sectional, in part because of lack of longitudinal data but also because of concerns that reporting of CB may vary over time. We found that a small proportion of individuals (1.6%) reported CB remittance at follow-up (online supplemental table S3). Part of this variation might be due to slightly different CB assessment items used for baseline and second assessment respiratory questionnaires. Longitudinal studies have shown that individuals with new or persistent CB symptoms are more likely to have continued smoking during follow-up, while individuals with ‘resolved’ CB (present at baseline but absent at follow-up) are more likely to have quit smoking.44 Given the temporal stability in the spatial patterns of air pollutant concentrations in the Netherlands,41 we expect participant-specific air pollution exposure to have been relatively constant between baseline and second assessment and do not expect differential exposure misclassification to be a major problem.

The use of ORs as an estimate of risk ratios (RR) may have overestimated risks of air pollution exposure in relation to respiratory health outcomes and we note that our effect sizes are small. However, ORs are accepted to provide a reasonable approximation of risk when disease prevalence/incidence is rare (ie, <10%)45 46; in our study CB prevalence and incidence were 6.1% and 4.5%, respectively, and chronic cough prevalence and incidence were 7.4% and 6.1%, respectively. Further, use of adjusted ORs from logistic regressions allowed direct comparison with the literature on air pollution and prevalent and incident CB and cough or phlegm symptoms, as most comparable previous studies report ORs. A limitation of the study is the large number of participants (n=19 546) with missing data for outcomes, confounders or exposures included in baseline regression models. Missingness did not appear to be at random (online supplemental table S2), which limited our ability to conduct imputation. While missing data does not invalidate findings from our study, it may impact their generalisability. A key strength of our study relative to past investigations was the large sample size, which allowed us to detect relatively small effect sizes and explore associations in potentially vulnerable subgroups. Our results are reinforced by recent cross-sectional analyses of the Lifelines cohort which showed positive associations between ambient air pollution exposure and wheeze and shortness of breath.47 Lastly, we used standardised air pollution models from the ELAPSE project, which are now being used in numerous epidemiological studies across Europe studying the association between air pollution and multiple health endpoints. The results of this study will therefore complement these ongoing studies.

In conclusion, this is the largest analysis to date to examine cross-sectional and longitudinal associations between ambient air pollution and CB. Higher NO2 and BC air pollutant concentrations, largely driven by local traffic-related pollution, were associated with increased odds of prevalent and incident CB. We observed higher effect sizes among women, never smokers and younger individuals. These results are of public health importance considering they are seen at relatively low air pollutant concentrations compared with urban locations within Europe and internationally.


The authors wish to acknowledge the services of the Lifelines Cohort Study, the contributing research centres delivering data to Lifelines and all the study participants. The Lifelines Cohort Study initiative has been made possible by subsidy from the Dutch Ministry of Health, Welfare and Sport, the Dutch Ministry of Economic Affairs, the University Medical Centre Groningen (UMCG the Netherlands), University Groningen and the Northern Provinces of the Netherlands.


Supplementary materials

  • Supplementary Data

    This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.


  • Contributors DD, JB, KdH and ALH contributed to the development of the study design; DD conducted the statistical analyses, wrote the first draft of the paper and acts guarantor of the manuscript. JB, KdH and ALH commented on the manuscript drafts and all authors approved the final version.

  • Funding This work was supported by the National Institute for Health Research (NIHR) Health Protection Research Unit in Environmental Exposures and Health, a partnership between Public Health England, the Health and Safety Executive and the University of Leicester, under grant number NIHR200901 (to ALH). The views expressed are those of the author(s) and not necessarily those of the NIHR, Public Health England, the Health and Safety Executive or the Department of Health and Social Care.

  • Competing interests JB reports grants from CIHR, grants from Canadian Respiratory Research Network (CRRN), personal fees from Canadian Thoracic Society, personal fees from CHEST, grants from Foundation of the MUHC, grants from Aerocrine, grants and personal fees from AstraZeneca, grants and personal fees from Boehringer Ingelheim, grants and personal fees from Grifols, grants and personal fees from GlaxoSmithKline, grants and personal fees from Novartis, grants and personal fees from Trudell, all outside the submitted work.

  • Patient consent for publication Not required.

  • Ethics approval Lifelines received ethics approval from the Medical Ethics Review Board of the University Medical Centre Groningen (approval letter ref: M07.052740).

  • Data availability statement Data may be obtained from a third party and are not publicly available. Lifelines provides access to data and samples for research on healthy ageing. Information on how to submit a data access application can be found at:

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.