CCBYNC Open access
Research

Long term exposure to ambient air pollution and incidence of acute coronary events: prospective cohort study and meta-analysis in 11 European cohorts from the ESCAPE Project

BMJ 2014; 348 doi: http://dx.doi.org/10.1136/bmj.f7412 (Published 21 January 2014) Cite this as: BMJ 2014;348:f7412
  1. Giulia Cesaroni, senior researcher1,
  2. Francesco Forastiere, research director1,
  3. Massimo Stafoggia, senior researcher1,
  4. Zorana J Andersen, associate professor in epidemiology23,
  5. Chiara Badaloni, research fellow1,
  6. Rob Beelen, senior researcher4,
  7. Barbara Caracciolo, researcher56,
  8. Ulf de Faire, senior professor of cardiovascular epidemiology7,
  9. Raimund Erbel, professor8,
  10. Kirsten T Eriksen, researcher2,
  11. Laura Fratiglioni, professor in geriatric epidemiology5910,
  12. Claudia Galassi, medical epidemiologist11,
  13. Regina Hampel, research fellow12,
  14. Margit Heier, research fellow1213,
  15. Frauke Hennig, research fellow14,
  16. Agneta Hilding, researcher15,
  17. Barbara Hoffmann, professor1416,
  18. Danny Houthuijs, senior researcher17,
  19. Karl-Heinz Jöckel, professor18,
  20. Michal Korek, doctoral student7,
  21. Timo Lanki, chief researcher19,
  22. Karin Leander, researcher7,
  23. Patrik K E Magnusson, professor20,
  24. Enrica Migliore, epidemiologist11,
  25. Caes-Göran Ostenson, professor15,
  26. Kim Overvad, professor2122,
  27. Nancy L Pedersen, professor of genetic epidemiology20,
  28. Juha Pekkanen J, professor19,
  29. Johanna Penell, researcher7,
  30. Göran Pershagen, professor7,
  31. Andrei Pyko, research fellow7,
  32. Ole Raaschou-Nielsen, head of research group2,
  33. Andrea Ranzi, project manager in environmental epidemiology23,
  34. Fulvio Ricceri, research fellow24,
  35. Carlotta Sacerdote, medical epidemiologist11,
  36. Veikko Salomaa, research professor25,
  37. Wim Swart, researcher17,
  38. Anu W Turunen, researcher19,
  39. Paolo Vineis, professor in epidemiology2426,
  40. Gudrun Weinmayr, research associate1427,
  41. Kathrin Wolf, research fellow12,
  42. Kees de Hoogh, senior research officer26,
  43. Gerard Hoek, associate professor4,
  44. Bert Brunekreef, professor428,
  45. Annette Peters, professor12
  1. 1Department of Epidemiology, Lazio Regional Health Service, 00198 Rome, Italy
  2. 2Danish Cancer Society Research Center, 2100 Copenhagen, Denmark
  3. 3Center for Epidemiology and Screening, Department of Public Health, University of Copenhagen, 1014 Copenhagen, Denmark
  4. 4Institute for Risk Assessment Sciences, Utrecht University, PO Box 80178, 3508 TD Utrecht, Netherlands
  5. 5Aging Research Center, Department of Neurobiology, Care Sciences and Society, Karolinska Institutet, S-113 30 Stockholm, Sweden
  6. 6Stress Research Institute, Stockholm University, SE-106 91 Stockholm, Sweden
  7. 7Institute of Environmental Medicine, Karolinska Institutet, Box 210 SE-171 77 Stockholm, Sweden
  8. 8West German Heart Center, University Hospital of Essen, 45122 Essen, Germany
  9. 9Stockholm Gerontology Research Center, 113 30 Stockholm, Sweden
  10. 10Division of Clinical Geriatrics, Karolinska University Hospital, SE-141 86 Stockholm, Sweden
  11. 11Unit of Cancer Epidemiology, “Città della Salute e della Scienza” Hospital, University of Turin, and Center for Cancer Prevention Piemonte, 10126 Turin, Italy
  12. 12Helmholtz Zentrum München, German Research Center for Environmental Health, Institute of Epidemiology II, 85764 Neuherberg, Germany
  13. 13Central Hospital of Augsburg, MONICA/KORA Myocardial Infarction Registry, D-86156 Augsburg, Germany
  14. 14IUF-Leibniz Research Institute for Environmental Medicine, 40225 Düsseldorf, Germany
  15. 15Department of Molecular Medicine and Surgery, Karolinska Institutet, Karolinska University Hospital, SE-171 77 Stockholm, Sweden
  16. 16Medical Faculty, University of Düsseldorf, 40225 Düsseldorf, Germany
  17. 17National Insititute for Public Health and the Environment, PO Box 1, 3720 BA Bilthoven, Netherlands
  18. 18Institute for Medical Informatics, Biometry and Epidemiology, University Hospital of Essen, D-45122 Essen, Germany
  19. 19National Institute for Health and Welfare “THL”, Department of Environmental Health, PO Box 95, FI-70701 Kuopio, Finland
  20. 20Department of Medical Epidemiology and Biostatistics, Karolinska Institutet, SE-171 77 Stockholm, Sweden
  21. 21Section for Epidemiology, Department of Public Health, Aarhus University, DK-8000 Aarhus, Denmark
  22. 22Department of Cardiology, Cardiovascular Research Center, Aalborg University Hospital, DK-9000 Aalborg, Denmark
  23. 23Environmental Health Reference Centre-Regional Agency for Environmental Prevention of Emilia-Romagna, 41121 Modena, Italy
  24. 24Molecular and Genetic Epidemiology Unit, HuGeF-Human Genetics Foundation-Turin, 10126 Turin, Italy
  25. 25THL-National Institute for Health and Welfare, Department of Chronic Disease Prevention, POB 30, FI-00271 Helsinki, Finland
  26. 26MRC-HPA Centre for Environment and Health, Department of Epidemiology and Biostatistics, Imperial College London, St Mary’s Campus, London W2 1PG, UK
  27. 27Institute of Epidemiology and Medical Biometry, Ulm University, 89069 Ulm, Germany
  28. 28Julius Center for Health Sciences and Primary Care, University Medical Center Utrecht, PO Box 85500, 3508 GA Utrecht, Netherlands
  1. Correspondence to: G Cesaroni g.cesaroni{at}deplazio.it
  • Accepted 4 December 2013

Abstract

Objectives To study the effect of long term exposure to airborne pollutants on the incidence of acute coronary events in 11 cohorts participating in the European Study of Cohorts for Air Pollution Effects (ESCAPE).

Design Prospective cohort studies and meta-analysis of the results.

Setting Cohorts in Finland, Sweden, Denmark, Germany, and Italy.

Participants 100 166 people were enrolled from 1997 to 2007 and followed for an average of 11.5 years. Participants were free from previous coronary events at baseline.

Main outcome measures Modelled concentrations of particulate matter <2.5 μm (PM2.5), 2.5-10 μm (PMcoarse), and <10 μm (PM10) in aerodynamic diameter, soot (PM2.5 absorbance), nitrogen oxides, and traffic exposure at the home address based on measurements of air pollution conducted in 2008-12. Cohort specific hazard ratios for incidence of acute coronary events (myocardial infarction and unstable angina) per fixed increments of the pollutants with adjustment for sociodemographic and lifestyle risk factors, and pooled random effects meta-analytic hazard ratios.

Results 5157 participants experienced incident events. A 5 μg/m3 increase in estimated annual mean PM2.5 was associated with a 13% increased risk of coronary events (hazard ratio 1.13, 95% confidence interval 0.98 to 1.30), and a 10 μg/m3 increase in estimated annual mean PM10 was associated with a 12% increased risk of coronary events (1.12, 1.01 to 1.25) with no evidence of heterogeneity between cohorts. Positive associations were detected below the current annual European limit value of 25 μg/m3 for PM2.5 (1.18, 1.01 to 1.39, for 5 μg/m3 increase in PM2.5) and below 40 μg/m3 for PM10 (1.12, 1.00 to 1.27, for 10 μg/m3 increase in PM10). Positive but non-significant associations were found with other pollutants.

Conclusions Long term exposure to particulate matter is associated with incidence of coronary events, and this association persists at levels of exposure below the current European limit values.

Introduction

According to the recent report on the Global Burden of Disease, throughout the world particulate air pollution is estimated to cause 3.1 million deaths a year and 22% of disability adjusted life years (DALY) due to ischaemic heart disease.1 Several cohort studies have reported that long term exposure to air pollution is associated with mortality, in particular cardiovascular mortality.2 3 4 5 6 7 8 9 10 11 12 13 The evidence of an effect on incidence of cardiovascular events—that is, acute myocardial infarction and unstable angina—is less consistent and requires further investigations.12 14 15

In the European Union the current annual limit for particulate matter <2.5 μm (PM2.5) is 25 µg/m3, which is far above that implemented in the United States (12 µg/m3). One obstacle in the European standard setting process in the past had been that the available estimates of the exposure-response associations of particulate matter were primarily based on studies conducted in North America.3 4 5 6 7 8 9 10 11 12 In response, the ESCAPE Study (European Study of Cohorts for Air Pollution Effects) was conducted between 2008 and 2012 to quantify the associations between exposures and health outcomes by using standardised methods for assessment of exposure and data from existing cohort studies.16 17

We estimated the association between long term exposure to particulate matter <2.5 μm (PM2.5), 2.5-10 μm (coarse), <10 μm (PM10) in aerodynamic diameter, soot (PM2.5 absorbance), nitrogen oxides (NOx and NO2), and traffic indicators and the incidence of coronary events. We evaluated effect modification by several individual characteristics, and we investigated the exposure-response relations of the pollutants below selected thresholds.

Video abstract

Methods

Design and population

This study is an analysis of cohort data obtained by ESCAPE to investigate the long term effects of exposure to air pollution on human health in Europe and a meta-analysis of the cohort specific results. The present study included 11 European cohorts from five countries with information about incident cases of acute coronary events and the most important potential confounders. The cohorts were in Finland (FINRISK)18; Sweden (the Swedish National Study on Aging and Care in Kungsholmen (SNAC-K), the Screening Across the Lifespan Twin Study, (SALT), the 60 year olds study, and the Stockholm Diabetes Prevention Program study (SDPP))19 20 21 22; Denmark (the Danish Diet, Cancer and Health cohort study [DCH])23; Germany (the Heinz Nixdorf Recall Study (HNR), the Cooperative Health Research in the Augsburg Region (KORA))24 25; and Italy (the European Prospective Investigation into Cancer and Nutrition in Turin (Epic-Turin), the International Study on Asthma and Allergies in Childhood in Turin and Rome (parents of representative samples of children participating in SIDRIA-Turin and SIDRIA-Rome)).26 27 These were all centres included in the ESCAPE intensive monitoring programme for particulate matter and nitrogen oxides.28 29 Table 1 shows specific settings, periods of recruitment, and person years at risk. With the exception of the Italian cohorts, all the cohorts enrolled representative samples of the adult population for the age groups considered in the single studies. Additional information on each cohort is presented in the appendix.

Table 1

Summary description of 11 European cohorts from five countries contributing data to study long term exposure to ambient air pollution and incidence of acute coronary events

View this table:

Outcome definition

For nine out of 11 cohorts, information on coronary events was derived by record linkage procedures with hospital discharge and mortality registries. We selected international classification of diseases (ICD) codes for “acute myocardial infarction” or “other acute and sub-acute forms of ischemic heart disease” (ICD-9-CM codes: 410, 411; ICD-10 codes: I21, I23, I20.0, I24) in principal diagnoses of hospital discharges. We also considered as incident cases those people who died outside hospital from ischaemic heart diseases, according to the death certificates (ICD-9: 410-414; ICD-10: I20-I25) and had no evidence of hospital admission for ischaemic heart disease in the 28 days before death and no evidence of admission for any cause in the two days before death. To identify incident cases, we excluded those who had had an acute coronary event or cerebrovascular event previously (previous admission to hospital with principal or secondary diagnoses with ICD-9 codes 410, 411, 430, 431, 434, 436 and ICD-10 codes I21, I23, I20.0, I24, I61, I63, I64). The reference period to evaluate previous events varied across the cohorts, with the minimum of three years in the two SIDRIA cohorts and all lifelong history for the KORA cohort. For the HNR study, incident cases were adjudicated by an independent end point committee based on medical records.24 For the KORA cohort, incident cases were determined by interview and medical history and validated through linkage with the myocardial infarction register.25

Exposure to air pollution

Exposure to air pollution in each area was estimated following a standard methods developed for the ESCAPE project and described elsewhere.28 29 Briefly, for each area under study, particulate matter of varying sizes measured in μm (shown as subscript)—that is, PM10, coarse PM, PM2.5, and PM2.5 absorbance (blackness of the PM2.5 exposed filter, determined by measurement of light reflectance as a marker for soot and black carbon)—was measured in 20 sites, and nitrogen oxides were measured in 40 sites in three separate two week periods (to cover different seasons) over one year (between 2008 and 2011). For each site, results from the three measurements were averaged to estimate the annual average, with adjustment for temporal variation by using a centrally located background reference site, which was operated for a whole year.30 31 By using several traffic and land use variables, we developed area specific land use regression (LUR) models to explain the spatial variation of each measured pollutant. These models were then used to estimate concentrations of air pollution at each participant’s residential address. Geographical variables typically evaluated include altitude, population density, industrial land use, green space, and traffic flows variables.30 31 In addition to concentrations of pollutants, we considered two traffic variables at the participant’s residence: traffic intensity on the nearest road (vehicles/day) and traffic load on major roads in a 100 m buffer (vehicles×m/day), defined as the sum of traffic intensity on roads with >5000 vehicles/day multiplied by the length of those roads in a 100 m buffer. To validate the models, we used the leave one out cross validation method— that is, systematically subtracting each of the monitoring points from the model one by one, and then comparing the predicted value for each monitoring location with the measured level at the location without using this measurement in the development of the model.30 31

Covariates

All cohorts had a common set of potential confounders and effect modifiers at baseline: marital status (coded as single, married/living with partner, divorced/separated, widowed; for the SDPP cohort only the binary variable “living with partner” was available; for SIDRIA-Rome all participants were living as couples at baseline), education (primary school or less, up to secondary school or equivalent, university degree and more), occupation (employed, unemployed, homemaker/housewife, retired), smoking status (current, former, never), duration of smoking (years), smoking intensity (cigarettes/day), hypertension, and diabetes. Most cohorts had information on additional cardiovascular risk factors such as body mass index (BMI; coded as <25, 25-29, ≥30), physical activity (<1 hour/week, about 1 hour/week, >2 hours/week), alcohol consumption (never, 1-3 drinks/week, 3-6 drinks/week, >6 drinks/week), and three cohorts had information on use of drugs for hormone replacement therapy. Four of 11 cohorts had data on cholesterol concentrations. Ten of 11 cohorts had data on noise exposure, estimated at the residential addresses from European noise exposure assessment in 2007, and we used an eight class categorical variable of 5 dB of exposure from 45 dB to 75 dB and more (see appendix). Each cohort also used an area based socioeconomic status indicator because socioeconomic indicators at the area level are predictors of morbidity, access to care, and lifestyle risk factors in addition to individual socioeconomic characteristics.32

Statistical analyses

To evaluate the association between exposure to air pollution and incidence of coronary events we performed the analyses in two stages. Firstly, we analysed each cohort using a common protocol for confounders, outcomes, and statistical modelling. We used Cox proportional hazards regression models (hazard ratios) with age as the time scale. A script with statistical code was provided to all cohort specific analysts. The results were evaluated centrally at the Department of Epidemiology in Rome. Secondly, we carried out a random effects meta-analysis to pool results.33 We calculated I2 statistics and P values for χ2 test from Cochran’s Q to quantify the heterogeneity among studies.34

We first estimated hazard ratios adjusted for age, sex, and year of enrolment only (model 1), then adjusted for a common set of individual covariates (marital status, education, occupation, smoking status, smoking duration, and smoking intensity; model 2), and then adjusted for area based socioeconomic status indicators (model 3). When investigating the effects of traffic variables we also adjusted all models for background NO2 levels. We estimated the effects of exposure to air pollution using fixed increments of pollutants (20 µg/m3 for NOx, 10 µg/m3 for PM10 and NO2, 5 µg/m3 for PM2.5 and PMcoarse,10−5/m for PM2.5 absorbance).

We performed several sensitivity analyses to deal with potential sources of bias and to consider the confounding role of additional cardiovascular risk factors. Firstly, we adjusted for possible intermediate variables (diabetes, hypertension) available in all the cohorts; then for physical activity, alcohol consumption, and BMI (available in eight cohorts); finally we added to previously mentioned factors the cholesterol concentration (available in four cohorts). We then analysed the role of living in low/high urbanised settings, of noise exposure, and of residential stability (restricting the analysis to those living at the same baseline address during the follow-up). We stratified the Cox model 3 for predictors that did not meet the proportionality hazard assumption. We analysed the influential role of the largest cohort by excluding the DCH study from the analysis. Finally, we evaluated the role of the performance of land use regression models, stratifying cohorts by cross validation R2> or ≤60%.

To explore the effect of clustering by area of residence in the association between exposure and mortality—that is, residents in the same area usually share similar characteristics (socioeconomic status, health, access to services) and have similar environmental and air conditions—we performed a frailty model to measure the role each area played.4

We evaluated potential effect modification by adding to model 3 an interaction term of exposure and one effect modifier at a time (sex, educational level, smoking status, BMI, hypertension, and residence in low or high urbanised settings). We used the likelihood ratio test to compare the models with and without interaction terms. For age during follow-up, we estimated effect modification by adding an interaction term between exposure and a time dependent categorical variable indicating age group (<60, 60-74, >75).

To provide information about the health effects below specific threshold values (20, 30, and 40 µg/m3 for PM10, and 15, 20, and 25 µg/m3 for PM2.5), we studied the effect of fixed increments of PM2.5 and PM10 among people with air pollution concentrations at residences only below these thresholds.

We used STATA software (versions 10, 11, and 12) for all the analyses, with the exception of frailty models for which we used R (www.r-project.org).

Results

Table 1 shows a summary description of the 11 European cohorts. Table 2 shows individual characteristics of participants and table 3 their exposures to air pollution. Additional characteristics of the participants are provided in table 4. The enrolment period covered 15 years, and the average follow-up was 11.5 years. The participating cohorts varied in characteristics, availability of data on covariates, and levels of exposure. The proportion of participants included in the study ranged from 82.4% to 99.3% of the original cohorts (overall 93.1%), after exclusion of missing values on any of the covariates in model 3. There were no differences in exposure levels between included and excluded participants. Among 100 166 participants included in the study and followed for 1 154 386 person years, there were 5157 incident cases. The average PM2.5 level at residence ranged from 7.3 µg/m3 in Sweden (SD=1.3) to 31.0 µg/m3 (SD=1.7) in northern Italy, and all other pollutants had similar patterns. All exposure models had good performance (R2 ≥0.61), and the prediction ability of models for nitrogen dioxides (based on a larger number of measurements sites) was slightly higher than those of particulate models.

Table 2

Individual baseline characteristics from 11 European cohorts from five countries contributing data to study long term exposure to ambient air pollution and incidence of acute coronary events. Values are numbers (percentages) unless stated otherwise

View this table:
Table 3

Air pollution exposure at residential address of individuals from 11 European cohorts from five countries contributing data to study long term exposure to ambient air pollution and incidence of acute coronary events. Figures are means and 5-95th centile ranges, correlation coefficients between PM2.5 and other pollutants, and R2 of models and their cross validation

View this table:
Table 4

Additional baseline characteristics of study populations in analysis of long term exposure to ambient air pollution and incidence of acute coronary events. Values are numbers (percentages) unless stated otherwise

View this table:

Table 5 shows the pooled hazard ratios with 95% confidence intervals for incidence of coronary events for fixed increments of all pollutants. We observed the strongest association for PM10. While in model 1 (adjusted for age, sex, and calendar period) all particulate matter indicators were strongly associated with incidence of coronary events, when we adjusted for marital status, education, occupation, smoking status, smoking duration, smoking intensity, and socioeconomic area indicator (model 3), only PM10 showed a significant association (hazard ratio 1.12, 95% confidence interval 1.01 to 1.25, for each 10 µg/m3 increase). There was also an association for PM2.5 (1.13, 0.98 to 1.30, for each 5 µg/m3 increase) and for coarse particles, whereas only small positive associations were found for nitrogen oxides. There was no evidence of an effect of traffic variables. There was no evidence of statistical heterogeneity between the cohort specific effect estimates (I2<5%). The figure shows the forest plots of 10 µg/m3 PM10 and 5 µg/m3 PM2.5 increments (model 3).

Figure1

Hazard ratios of incident coronary events per 10 µg/m3 PM10 and 5 µg/m3 PM2.5. Cohort specific and meta-analytic results

Table 5

Association between exposure to pollutants and incidence of coronary events. Results expressed per fixed increments

View this table:

Table 6 shows the results of the sensitivity analyses for PM10 and PM2.5 from the base model (model 3), reporting the number of cohorts and participants included in each analysis. Overall, there were only a few marginal changes in the effect estimates across the various models, and there was no evidence of heterogeneity among cohorts. When we considered diabetes and hypertension in the adjustment, the effect estimates for PM10 and for PM2.5 were slightly lower (hazard ratio 1.11 (95% confidence interval 1.00 to 1.24) and 1.11 (0.96 to 1.28), respectively). When we considered physical activity, alcohol consumption, and BMI in the eight cohorts with available information, the association between particulate matter and incidence of coronary events remained stable. Adjustment for cholesterol in the four cohorts with this variable did not alter the results. Adjustment for location of residence (urban, suburban/rural) slightly increased the effect estimates. In 63 121 participants who did not change their address during the follow-up, the effect estimate of air pollution was higher than in the whole population (hazard ratio 1.16 (1.01 to 1.32) for 10 μg/m3 PM10, and 1.18 (0.98 to 1.42) for 5 μg/m3 PM2.5), possibly because of more accurate exposure assignment. Exclusion of the large DCH study did not influence the results. When we restricted the meta-analysis to nine cohorts with exposure models with high cross validation (R2 >60) for PM10 and six cohorts for PM2.5, we found stronger associations (hazard ratio 1.18 (1.05 to 1.33) for 10 μg/m3 PM10, and 1.35 (1.04 to 1.74) for 5 μg/m3 PM2.5). When we took individual and area level covariates into account, we found no evidence of clustering in the neighbourhoods (data not shown).

Table 6

Incidence of coronary events, results from sensitivity analyses. Figures are pooled hazard ratios (and 95% CI) for 10 µg/m3 PM10 and 5 µg/m3 PM2.5, I2 statistics, and P value for heterogeneity

View this table:

Effect modification by age showed stronger PM2.5 effects for those aged 60-74 (hazard ratio 1.25, 95% confidence interval 1.03 to 1.51) and for those >75 (1.18, 0.85 to 1.64) than among those aged under 60 (0.91, 0.71 to 1.15), with P=0.11 for effect modification. For all other effect modifiers, the P value of effect modification was far from significant (≥0.22).

Table 7 shows the results of the threshold analysis for PM2.5 and PM10. When we restricted the analysis to participants from nine cohorts exposed to PM2.5 concentrations below the current European limit value (<25 µg/m3), there was a 18% increased risk per 5 µg/m3 increase in PM2.5 (hazard ratio 1.18, 95% confidence interval 1.01 to 1.39). Similarly, for PM10 below the current European limit value (40 µg/m3) there was a 12% increased risk per 10 µg/m3 increase in PM10 (1.12, 1.00 to 1.27). The effect was present even at exposure levels <15 µg/m3 for PM2.5 and <20 µg/m3 for PM10. This result was not an artefact caused by cohort selection at different thresholds as it persisted when we the restricted analysis to the seven cohorts with data available in each threshold (table 7, right columns).

Table 7

Results from random-effects meta-analyses for adjusted association between incidence of coronary events and exposure to PM10 and PM2.5 below various threshold values. Figures are pooled hazard ratios (95% CI) for 10 µg/m3 PM10 and 5 µg/m3 PM2.5, I2 statistic, and P value for heterogeneity

View this table:

Discussion

This European multicentre study found that long term exposure to particulate air pollution is associated with an increased risk of first coronary events. For both particulate matter with diameter <10 μm and <2.5 μm (PM10 and PM2.5), we found associations at levels below the current European limits. The results were substantially robust to adjustment for confounding and model specification.

Comparison with other studies

Our effect estimates are similar to or higher than those from two cohort studies from the US. For an increase of 10 µg/m3 in PM2.5, Miller and colleagues found a 21% (hazard ratio 1.21, 95% confidence interval 1.04 to 1.42) increased risk of first coronary heart disease in postmenopausal women, and a (non-significant) 6% higher risk of myocardial infarction.9 In the US Nurses’ Health Study,11 there was an increased incidence of coronary heart disease for a 10 µg/m3 increase in PM2.5 (hazard ratio 1.11, 0.79 to 1.55) and for a 10 µg/m3 increase in coarse particles (1.04, 0.82 to 1.32). In the UK, Atkinson and colleagues found a non-significant 1% (−2% to 5%) higher risk in incident myocardial infarction per 3 µg/m3 increase in PM10.14 In the California Teachers Study, Lipsett and colleagues found no evidence of increased myocardial infarction incidence per 10µg/m3 increments in PM2.5 and PM10.12 The results are also in line with the evidence produced so far by studies on the short term effects of particulate matter on ischaemic heart disease.35 36

We found a suggestion of a higher risk associated with particulate exposure in participants aged over 60. Evidence on which age range is the most susceptible is mixed. Miller and colleagues found a higher risk in those aged >70 compared with <70,9 while Gan and colleagues found higher risks in first admissions to hospital for coronary heart disease in both those aged <60 and 60-69 than those aged ≥70.5

We found an effect of PM2.5 below 25, 20, and 15 μg/m3, and our results highlight the possible harmful health effects of fine particles well below the current EU annual limit for PM2.5. It is clear that generalisability of these findings to the entire European population is not straightforward because the cohorts are not representative of the European population. Besides the differences among the cohorts in age range and underlying risk profiles, however, we detected no heterogeneity in the effect estimates among the cohorts.

There are several possible mechanisms through which exposure to air pollution can affect the cardiovascular system, including systemic inflammation, systemic oxidative stress, thrombosis and coagulation, changes in blood pressure, progression of atherosclerosis, and reduced heart rate variability.15 Most of this evidence comes from studies of short term exposures and short term responses of cardiovascular function. The pathways implicated provide plausible biological mechanisms based on repeated exacerbations of cardiovascular risk factors potentially leading to long term progression of coronary artery disease.15 37

It is unlikely that our findings are driven by an excess in coronary disease mortality because in these and some other cohorts participating in the ESCAPE project we could not find much evidence of an association between exposure to air pollution and mortality for ischaemic heart disease or myocardial infarction.38

Our main positive results for coronary events are for an association with particulate matter, and the evidence for other pollutants such as nitrogen oxides could be considered as suggestive based on the small positive estimates. Several predictor variables explained the spatial variation of PM10 and PM2.5 in ESCAPE study areas,30 including not only traffic variables but also population density, industrial sources, urban green, and altitude. Estimated PM10 at place of residence therefore represents not only traffic emissions. In contrast, the spatial variation in PM2.5 absorbance was more exclusively explained by traffic variables. Our findings suggest that sources other than vehicular traffic can have an important role on the risk of acute ischaemic heart disease.

Strengths and weaknesses

In this first multicentre European study to investigate long term effects of air pollution on incidence of coronary events, we used a standardised approach for assessment of exposure with respect to measurements and modelling.28 29 30 31 Both the measurements and the land use regression models were centrally validated, and data analyses were centrally planned and conducted in the same way in each cohort. The European locations captured in this project provide a wide range of exposures and mixtures that vary from north to south. Although the design of the study did not allow us to benefit from the large variability of pollutant exposure between cohorts, we took advantage of the exposure variation within cohorts. The use of existing cohorts, not designed for air pollution epidemiology, might be seen as a limitation of the ESCAPE project. Indeed, we were able to perform a planned meta-analysis and make the most of existing studies in an efficient way. Therefore, we consider these data the first systematic assessment of the impact of ambient particles on incident coronary artery disease in Europe, highlighting that the burden of disease might be underestimated when estimates of mortality are considered alone.15 39

The assessment of exposure was conducted in 2008-11, while the cohorts were enrolled over 15 years, starting in 1992. The use of modelled concentrations based on measurements taken long after the beginning of the study assumes that spatial contrasts were stable over time, an assumption that has been supported from observations in different settings.40 41 42 In the ESCAPE project considerable effort was made to take this problem into account; when possible, we back-extrapolated air pollution concentrations.17 Analyses of the association between back-extrapolated concentrations and incidence of coronary events made no important difference.

We collected individual data on many cardiovascular risk factors to deal with confounding. When we compared the results of the “base model” with a model containing additional cardiovascular risk factors, we obtained similar results. It should be noted that variables used to adjust for area level socioeconomic status were different from cohort to cohort, ranging from the percentage of low income inhabitants in a 5 km grid in KORA to a census block deprivation index in the Italian cohorts. We could not determine whether this might have affected the results, but general adjustment for area level socioeconomic indicators did not have much influence in cohort specific analyses. Data on the most relevant cardiovascular risk factors (smoking, diabetes or hypertension, BMI, physical activity) were available for almost all cohorts, therefore a strong bias in the effect estimates caused by confounding is unlikely. Nevertheless, we cannot rule out residual confounding by, for example, differences in diet. In addition, we observed evidence for effect modification, but even within this large multicentre study, the power to reliably detect effect modification is limited and therefore we selected a limited set of interactions to evaluate a priori.

All cohorts, with the exception of the DCH study, enrolled fewer than 10 000 participants, making the statistical power to detect effects of air pollution in single studies limited. The characteristics of participants were heterogeneous, with two cohorts having younger participants than the others. Younger age could be a reason for the smaller effects found in some of the cohorts. Finally, in nine out of 11 cohorts outcomes were ascertained from routine administrative databases, which could have led to a less precise effect estimate from misclassification of the outcome.

When we included only the studies with the best exposure assessment models (performance of the land use regression leave-one-out cross validation R2 >60%), we found increased effect estimates. This could be because of reduced misclassification of exposure. The cohorts with validation R2 ≤60%, however, had young participants and (in particular the Italian ones) had higher rates of smoking and a better exposure assessment might be not the only explanation.

Conclusions

Our study suggests an association between long term exposure to inhalable particulate matter and incidence of coronary events. These associations remained for exposure concentrations below the current European limits. The results of this study, together with other ESCAPE findings, support lowering of European limits for particulate air pollution to adequately protect public health.

What is already known on this topic

  • Ambient particulate matter air pollution is estimated to cause 3.1 million deaths worldwide per year, and 22% of DALYs (disability adjusted life years) from ischaemic heart disease

  • The association between long term exposure to air pollution and incidence of coronary events remains controversial

What this study adds

  • There is an association between outdoor particulate matter and incidence of acute coronary events, even for exposure levels below the current European limits

  • The burden of disease attributable to outdoor particulate matter might be underestimated if only estimates of mortality are considered

  • The results of this study support lowering of the EU limits for particulate matter air pollution

Notes

Cite this as: BMJ 2014;348:f7412

Footnotes

  • We thank Marjan Tewis, Marieke Oldenwening, Marloes Eeftens, Ulrich Quass, Lee-Jane Sally Liu for their help with exposure assessment and data management within ESCAPE.

  • For the Finnish part, additional funding came from the Academy of Finland (project number 129317). Exposure assessment was performed by Tarja Yli-Tuomi, Pekka Taimisto, and Arto Pennanen from the Department of Environmental Health. Mortality, area-level SES, and building data were provided by Statistics Finland.

  • The four Swedish cohorts (SNAC-K, SALT, 60y, SDPP) were partially funded by the Swedish Environmental Protection Agency, the Swedish Council for Working Life and Social Research and the Swedish Heart-Lung Foundation. The SALT cohort was additionally supported by NIH grant AG-08724. The 60 year cohort was additionally funded by the Stockholm County Council and the Swedish Research Council (longitudinal research and 0593). The SDPP cohort was additionally funded by the Stockholm County Council; the Swedish Research Council; the Swedish Diabetes Association; and the Novo Nordisk Scandinavia. We thank the Heinz Nixdorf Foundation, the German Ministry of Education and Science (BMBF), the German Aerospace Center (Deutsches Zentrum für Luft- und Raumfahrt (DLR)), and the German Research Foundation (DFG HO 3314/2-1 and JO170/8-1), Bonn, Germany, for their generous support of Heinz Nixdorf Recall study. The KORA research platform and the MONICA Augsburg studies were initiated and financed by the Helmholtz Zentrum München, German Research Center for Environmental Health, which is funded by the German Federal Ministry of Education and Research and by the State of Bavaria. The SIDRIA cohort study was partially funded by the Italian Ministry of Health. Giuseppe Costa provided data from the Turin Longitudinal Study for the follow-up of the SIDRIA-Turin cohort. Exposure assessment in Turin was performed by Daniela Raffaele and Marco Gilardetti. Noise assessment in Turin was performed by Daniele Grasso and Jacopo Fogola (Regional Agency for the Protection of the Environment of Piedmont). We thank Simone Bucci, Patrizio Pasquinelli, and Eleonora Zirro for exposure assessment in Rome.

  • Contributors: GC contributed to the design, exposure assessment, data analyses, and drafted the manuscript; FF contributed to the design, provided local cohort data and drafted the manuscript; MS contributed to the design, statistical script, and data analyses; ZJA, RH contributed to the statistical script and data analysis; CB, FH, EM, JPen, APy, FR, and AT contributed to the data analyses; RB, GH contributed to the design, exposure assessment and statistical script; BC, UdF, RE, LF, MH, AH, KHJ, KL, PKEM, CGO, KO, NLP, JPek, CS, and VS provided local cohort data; CG contributed to exposure assessment and provided local cohort data; BH, ORN, GP, and PV contributed to the design and provided local cohort data; MK contributed to the exposure assessment and the data analyses; KTE, DH, TL, AR, WS, and KdH contributed to exposure assessment; GW and KW contributed to the statistical script; BB and APe contributed to the design and drafted the manuscript. All authors contributed to critical reading of and comments to the manuscript, interpretation of data and approved the final draft. GC is guarantor.

  • Funding: The research leading to these results has received funding from the European Community’s Seventh Framework Program (FP7/2007-2011) under grant agreement No 211250. The sources of funding had no role in study design, data collection, analyses, interpretation, and decision to submit the article for publication.

  • Competing interests: All authors have completed the ICMJE uniform disclosure form at www.icmje.org/coi_disclosure.pdf and declare: no support from any organisation for the submitted work; no financial relationships with any organisations that might have an interest in the submitted work in the previous three years; no other relationships or activities that could appear to have influenced the submitted work.

  • Ethical approval: The original cohort studies were approved by appropriate institutional review boards complying with all relevant national, state, and local regulations.

  • Transparency statement: The lead author affirms that this manuscript is an honest, accurate, and transparent account of the study being reported; that no important aspects of the study have been omitted; and that any discrepancies from the study as planned have been explained.

  • Data sharing: meta-analytic data and statistical code are available from the corresponding author.

This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 3.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 and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/3.0/.

References