DOI: https://doi.org/10.4414/smw.2012.13681
Abbreviations
95% CI Ninety-five percent confidence interval
AP Atmospheric pressure
AUC Area under the curve
CLR Conditional logistic regression
COPD Chronic obstructive pulmonary disease
FEV1/FVC Forced expiratory volume in 1 second over forced vital capacity
FEV1 Forced expiratory volume in one second
µg/m3 Microgram per cubic meter
NO2 Nitrogen dioxide
OR Odds ratio
PM2.5 Particulate matter of 2.5 micrometers in diameter
PM10 Particulate matter of 10 micrometers in diameter
RH Relative humidity
SAPALDIA Swiss Cohort Study on Air Pollution and Lung Diseases in Adults
T Temperature
TSP Total suspended particulates
The effect of air pollution on respiratory morbidity and mortality in adults is well documented in epidemiologic studies which examine severe outcomes including emergency room visits or hospitalisations. There is also evidence to suggest that individuals with underlying chronic respiratory disease including asthma or chronic obstructive pulmonary disease (COPD) are more susceptible to the acute effects of air pollutants including aggravation of respiratory symptoms, hospitalisations, and mortality [1–4].
Less is known about the effect of air pollution on respiratory outcomes in the context of a primary healthcare setting, such as visits to the doctor for general consultation or doctors’ house calls in relation to respiratory symptoms or disease [5–8]. Symptomatic respiratory episodes leading to subsequent primary care are likely to be less severe than those that lead to emergency care, but they may also be more common, resulting in time off from work for the doctor visit or for home recuperation, suggesting a large impact in terms of public health and societal costs. Doctor visits may serve as a useful indicator for the health effects of air pollution in the acute time domain. Whether adults with underlying chronic respiratory diseases are more susceptible to respiratory symptom-related doctor visits in association with air pollutant exposure is unknown.
The purpose of this study was to investigate the effects of current levels of ambient nitrogen dioxide (NO2) and total suspended particulates (TSP), as markers of traffic-related particulate matter, on doctor visits for respiratory problems in a large diary panel study of adults that was conducted in 1992 and 1993 as part of the Swiss Cohort Study on Air Pollution and Lung and Heart Diseases in Adults (SAPALDIA). While the data from the SAPALDIA diary study was collected nearly 20 years ago, the size of the diary panel study is unprecedented, and the time patterns of exposure-response over a long duration have not been addressed in the literature. Furthermore, while air quality in Switzerland may have changed over time, air quality in Switzerland from the early 1990’s is still a reasonable proxy of air quality in many other places in the world. Many of the diary study participants were also present with symptoms consistent with asthma, chronic bronchitis, and COPD, thus providing a unique opportunity to evaluate association lag structures in short-term effects of air pollutant levels related to primary care in adults with different patterns of underlying chronic disease.
SAPALDIA is a multicenter and population based prospective cohort study consisting of a random sample of 9,561 adults between the ages of 18–60 recruited from eight regions (Aarau, Basel, Davos, Geneva, Lugano, Montana, Payerne, and Wald) in Switzerland. Baseline examination occurred in 1991 where participants were administered a detailed respiratory health questionnaire and underwent allergy, lung function, and methacholine challenge testing [9, 10]. A subset of the study population at baseline (n = 5,383) was invited to enrol in a longitudinal diary study spanning 1992 and 1993 [11, 12]. Participants with a higher probability of reporting respiratory symptoms in the baseline survey were over-sampled for inclusion in the diary study to enhance the study's statistical power to detect associations between air pollution and lung diseases. Accordingly, meeting one of the following criteria assessed in the baseline survey was sufficient for a subject to be eligible for participation in the diary study: (1.) an affirmative answer to at least one of a set of 10 questions on respiratory symptoms; (2.) a ratio between forced expiratory volume in 1 second (FEV1) and forced vital capacity that was less than 80% of the value predicted by European Community Coal and Steel norms; (3.) a decline in FEV1, of at least 20% from baseline during the bronchochallenge test; and (4.) satisfying none of the above criteria and being a non-smoker. All participants fulfilling at least one of the first three criteria at the baseline survey were invited to enroll in the diary study. Participants fulfilling the fourth criterion were restricted to a random sample of 150 non-symptomatic non-smokers per area [11]. Of the 5,383 individuals invited to enrol, 3,279 participated (60.9%) [11, 12]. Ethical approval for this study was given by the ethics committees of the participating cantons in each study region.
The diary study participants completed up to 6 diary periods of 4 weeks in duration, in 1992 and 1993, during which they recorded daily occurrence of respiratory symptoms and visits to the doctor for respiratory symptoms. Of the 3,279 diary study participants, we included the 459 individuals in this analysis who reported at least one doctor visit for a total of 1,048 visits. A time-stratified case-crossover analysis was developed, where for each participant a day with a reported doctor visit (case-day) was matched with all days absent of a doctor visit (control-day) that occurred in the diary period by weekday [13]. Case and control-days were excluded from the analysis if the participant did not report being present in the study area on that day. For this analysis, there was complete information on NO2 exposure and other covariates in 914 case-days (i.e., reported doctor visits) and 2,362 control-days; for TSP, there were 892 case-days and 2,230 control-days with complete information.
Of the 459 participants, 374 were divided into subgroups of having asthma (n = 78), chronic bronchitis (n = 189) or COPD (n = 54), or who were healthy (n = 116). Presence of asthma was derived from positive responses to the following questions, “Have you ever had asthma?’’ and ‘‘was this confirmed by a doctor?’’. Chronic bronchitis was defined as self-report of cough or phlegm during the day or at night on most days for as much as 3 months each year for ≥2 years. COPD was defined as a ratio of forced expiratory volume in 1 second over forced vital capacity (FEV1/FVC) less than 70% and no previous history of asthma or current asthma medication use. Thus, while participants overlapped between the asthma and chronic bronchitis subgroups, the asthma and COPD subgroups were mutually exclusive. For this analysis, participants were defined as healthy if they did not meet the criteria for inclusion for asthma, chronic bronchitis, or COPD and if they were not present with bronchial hyper-responsiveness at the baseline survey; 14 of the 116 participants in the healthy subgroup were current smokers.
Doctor visits for respiratory symptoms were divided into two subcategories: lower respiratory symptoms (wheezing or tightness of chest ortrouble breathing orcough or phlegm during the day or the preceding night), and upper respiratory symptoms (sore throat or hoarseness or runny nose during the day or the preceding night and no lower respiratory symptoms).
Additional details on exposure measurements are provided in Monn et al. [14]. NO2 and TSP were monitored at fixed sites by local authorities in the eight SAPALDIA regions, and 24-hour means of NO2 (µg/m3) and half-hourly means of TSP (µg/m3) were used in this analysis. Particulate matter of 2.5 nor 10 micrometers in diameter (PM2.5and PM10, respectively) were not monitored on a daily basis at the time of the diary study, but a previous exposure assessment study demonstrated a high correlation between weekly averages of TSP and PM10 in the SAPALDIA regions [15]. In the same study, a high correlation was also observed between PM10 and NO2 in the urban regions. At the time of the diary study, emission sources for NO2 in Switzerland were predominately automobile traffic [15, 16].
Meteorological parameters including 24-hour means of temperature (T), relative humidity (RH), and atmospheric pressure (AP), were also collected for each SAPALDIA region. For chronological gaps of ≤2 consecutive days in the diary period with missing data for NO2, TSP, T, RH, or AP, linear interpolation between the day before and the day after was used to impute the missing data. Of the 3,276 person-days of observation with complete information specific to the NO2 analysis, 465 person-days consisted of imputed pollutant or meteorological data; 356 of the 3,122 person-days specific to the TSP analysis consisted of imputed pollutant or meteorological data.
Conditional logistic regression (CLR), performed with the PROC PHREG procedure in SAS V 9.2 (SAS Institute, Cary, NC, 2008), was used to estimate odds ratios (OR) and 95% confidence intervals for doctor visits for respiratory symptoms (combination of upper and lower respiratory symptoms) as a function of NO2 or TSP concentration (per 10 µg/m3) in a polynomial distributed lag model. The polynomial distributed lag model was derived from NO2 and TSP concentrations on the respective case or control-day of observation and the previous 13 days to account for substantial correlation between air pollution concentrations on days close together [17, 18], including terms up to the 5th degree order constrained in such a way that the resulting function had linear tails between lags 0 and 1 and between lags 12 and 13. Longer lag structures up to two weeks were also developed in this analysis to enhance stability of the distributed lag estimates, and because the previous literature demonstrates that lags beyond one week have some influence on respiratory mortality and morbidity [19, 20].
All models were adjusted for time-varying covariates including quadratic polynomial distributed lag function of T, RH, AP, T2, RH2,and RH*T on the respective day and the 6 preceding days (i.e., lag 0 to 6) up to ≤2 degrees for each term, and average T, RH, and AP from lag 7 to lag 13 days. Finally, there was additional adjustment for weekly average influenza incidence at the regional level, which was provided by the Federal Office of Public Health. All models were stratified by gender and smoking status (current smoker, former or never smoker); there was insufficient power to distinguish the effect of exposure between former and never smokers, particularly in the COPD subgroup, so we chose to combine former and never smokers. Doctor visits were also limited to lower respiratory symptoms in a secondary analysis.
For each subgroup, the average percent change in doctor visits for respiratory symptoms was estimated over 1 or 2 weeks following a 10 µg/m3 increase in NO2 and TSP concentration on a given day by computing the area under the curve (AUC), which was derived from summing logits lag 0 to 6 or lag 0 to 13 days from the CLR models and then exponentiating their sum. The difference in AUC between mutually exclusive subgroups was estimated to evaluate whether time patterns were distinguishable between the subgroups.
Different patterns of characteristics emerged between the subgroups as shown in table 1. Participants with COPD were more likely to be male, considerably older, report more than one doctor visit, and be a current smoker compared with participants with asthma and healthy participants. The majority of doctor visits were for lower rather than upper respiratory symptoms, but healthy participants had the highest proportion of doctor visits for upper respiratory symptoms.
Comparison of characteristics at the baseline SAPALDIA survey between participants who reported or did not report doctor visits was also evaluated within the asthma and COPD subgroups (supplementary file table 1). Among participants with COPD, the distributions of age, current smoking status, symptoms consistent with chronic bronchitis, and percent predicted FEV1 were not considerably different between those who reported and did not report doctor visits. Asthmatics who reported doctor visits had significantly lower percent predicted forced expiratory volume in one second (FEV1) and a significantly higher prevalence of self-reported asthma medication use.
The mean NO2 concentration observed on case and control-days, and the differences between them from the same risk set are summarised in table 2. On average, higher NO2concentrations were observed on case-days relative to control-days (i.e., at an exposure lag of 0 days). From an exposure lag of 1 day onward, the difference in NO2 concentration, relative to control-days, continued to rise and peaked at an exposure lag of 4 days, before trending downward into negative territory by a lag of 8 days. A high correlation was observed between NO2 and TSP concentrations (Pearson correlation r2 = 0.72) at exposure lag of 0 days.
The maximum OR for a doctor visit related to respiratory symptoms per 10 µg/m3 increase in NO2 in all participants occurred at an exposure lag of 4 days (OR: 1.038, 95%CI: 0.999–1.077), indicating a small increase in doctor visits. A bimodal lag distribution of peak effects of NO2 occurred in all subgroups, but the timing and strength of the peak effects varied by condition (fig. 1, table 3). In the asthma subgroup, the peak effects occurred at exposure lags of 5 and 13 days, but at shorter lags for participants with chronic bronchitis and COPD, all of which were statistically significant (p <0.05). In the healthy subgroup, the peak effects occurred at exposure lags of 5 and 10 days. Statistically significant inverse associations were also observed in the healthy and chronic bronchitis subgroups at lags of 2 and 9 days, respectively. Exposure lag patterns were broadly similar when doctor visits were limited to lower respiratory symptoms (data not shown).
Transformation of the AUC’s from the CLR models into average percent changes in doctor visits following an increase NO2 exposure is illustrated in figure 2. The largest increase in doctor visits was observed in the COPD subgroup; a 10 µg/m3 increase in NO2 concentration on a given day was associated with a 9.1% increase (95%CI: 3.2–15.4%) in doctor visits for respiratory symptoms on average over the first week following exposure. The increase in doctor visits for the COPD subgroup was reduced to 4.4% over two weeks following exposure, but remained statistically significant (95%CI: 0.6–8.3%). Average increases in doctor visits in both time periods were also observed in the asthma and chronic bronchitis subgroups, but were only statistically significant for those with chronic bronchitis over one week following an increase in exposure (4.2%, 95%CI: 1.2–7.2%). Average decreases in doctor visits were observed in the healthy subgroup for both time periods, but these changes were not statistically significant. The differences in average change in doctor visits were statistically significant in both time periods between the COPD and healthy subgroups, and over the first week between the COPD and asthma subgroups.
The peak effects of NO2 occurring with shorter exposure lags in the COPD subgroup, as shown in figures 1 and 2, were largely observed in current smokers (supplementary file table 2). There was little difference in effect lag patterns observed between current and non-smokers with chronic bronchitis. For participants with asthma, the peak effects of NO2 in both groups occurred at exposure lags of approximately 5–6 days and 13 days.
Differences in exposure lag patterns by gender also varied between subgroups (supplementary file table 3). In the COPD subgroup, the peak effects of NO2 occurred at exposure lags between 0 and 3 days in females and between 2 and 4 days in males. For the chronic bronchitis subgroup, elevated ORs were observed mainly in males at exposure lags between 0 and 1 days and between 4 and 6 days. Exposure lag patterns were similar between males and females with asthma.
The results for TSP and doctor visits (see supplementary file) were broadly similar to those observed for NO2. The adjusted OR for doctor visits in all subjects was highest at an exposure lag of 4 days (OR: 1.038, 95%CI: 1.010–1.065). In comparison with the asthma and healthy subgroups, the peak effect of TSP occurred at a shorter exposure lag in the COPD subgroup (supplementary file fig. 1A–1D, table 5), and the effect was largely apparent in smokers and females of this subgroup (supplementary file tables 6 and 7).
Table 1: Characteristics of participants who reported doctor visits for respiratory symptoms in the SAPALDIA diary study. | ||||
Participants with at least one visit to the doctor | ||||
Characteristics | Asthma | Bronchitis | COPD | Healthy |
Number of subjects | 78 | 189 | 54 | 116 |
n (%)* | n (%) | n (%) | n (%) | |
Sex | ||||
Male | 31 (39.7) | 81 (42.9) | 29 (53.7) | 34 (29.3) |
Female | 47 (60.3) | 108 (57.1) | 25 (46.3) | 82 (70.7) |
Age | ||||
<35 | 17 (21.8) | 39 (20.6) | 7 (13.0) | 33 (28.5) |
35–45 | 24 (30.8) | 50 (26.5) | 10 (18.5) | 41 (35.3) |
>45 | 37 (47.4) | 100 (52.9) | 37 (68.5) | 42 (36.2) |
mean (SD) | 43.6 (10.1) | 44.3 (9.9) | 48.2 (9.9) | 41.5 (10.2) |
Smoking status | ||||
Current smoker | 14 (17.3) | 80 (42.3) | 24 (46.0) | 14 (10.4) |
Former smoker | 24 (30.8) | 37 (19.6) | 13 (24.1) | 29 (25.0) |
Never smoker | 40 (51.3) | 72 (38.1) | 17 (31.5) | 73 (63.0) |
Medication** | ||||
Current use of asthma medication | 35 (44.9) | 22 (11.6) | ||
Use of inhalation steroids in last 12 months | 18 (23.1) | 12 (6.4) | ||
Lung function | ||||
Percent predicted FEV1, mean (SD) | 88.9 (21.9) | 100.2 (20.3) | 91.9 (16.2) | 112.1 (13.7) |
FEV1/FVC, %, mean (SD) | 71.7 (11.8) | 75.5 (10.3) | 63.4 (6.2) | 81.3 (5.6) |
PEF***, %, mean (SD) | 6.77 (2.33) | 7.04 (2.46) | 6.47 (2.55) | 7.63 (2.40) |
Visits to the doctor | ||||
1 | 43 (55.1) | 95 (50.3) | 24 (44.4) | 65 (56.0) |
>1 | 35 (44.9) | 94 (49.7) | 30 (55.6) | 51 (44.0) |
Visits to the doctor | ||||
Type | Asthma | Bronchitis | COPD | Healthy |
All visits | 201 | 478 | 166 | 231 |
Visits for upper respiratory symptoms | 20 | 63 | 42 | 86 |
Visits for lower respiratory symptoms | 181 | 415 | 124 | 145 |
* Summarised in n (%) unless noted otherwise. ** Self-reported medication use ascertained from SAPALDIA 1 survey (1991). *** Peak expiratory flow variability expressed as amplitude as a percent of mean. |
Table 2: Distribution of 24-hr mean concentrations of NO2 (µg/m3) for case and control-days from all participants (n = 457). | ||||||||
Lag | Case/Control-days* | Mean | SD | 10th | 25th | 50th | 75th | 90th |
Case | 34.4 | 22.2 | 8.5 | 15.4 | 31.1 | 49.7 | 63.3 | |
Control | 33.4 | 20.1 | 8.4 | 15.8 | 31.8 | 48.5 | 59.1 | |
Difference | 0.9 | 13.9 | –14.2 | –5.6 | 0.5 | 7.2 | 17.1 | |
1 | Case | 33.0 | 22.3 | 7.9 | 14.8 | 29.3 | 47.3 | 62.2 |
Control | 32.4 | 19.8 | 8.1 | 14.7 | 30.0 | 46.3 | 58.3 | |
Difference | 0.6 | 14 | –14.6 | –6.1 | 0.6 | 7.2 | 17 | |
2 | Case | 32.8 | 22.1 | 7.5 | 14.1 | 30.1 | 46.8 | 63.4 |
Control | 32.2 | 19.7 | 7.6 | 14.7 | 30.4 | 46.4 | 57.3 | |
Difference | 0.7 | 13.9 | –14.2 | –5.4 | 0.6 | 7.6 | 17.1 | |
3 | Case | 33.3 | 22.0 | 8.2 | 15.4 | 30.3 | 47.1 | 62.3 |
Control | 32.7 | 19.5 | 7.9 | 15.2 | 31.8 | 46.9 | 58.2 | |
Difference | 0.8 | 13.5 | –14.9 | –5.5 | 0.4 | 7.2 | 15.5 | |
4 | Case | 33.6 | 22.5 | 8.2 | 15.4 | 31.1 | 47.2 | 61.8 |
Control | 32.7 | 19.3 | 8.3 | 15.9 | 32.5 | 46.1 | 59.3 | |
Difference | 1.0 | 13.6 | –14.7 | –4.9 | 0.8 | 6.7 | 16.6 | |
5 | Case | 33.7 | 22.5 | 7.9 | 15.3 | 30.7 | 47.1 | 61.9 |
Control | 33.0 | 19.8 | 8.2 | 15.1 | 31.9 | 47.3 | 59.3 | |
Difference | 0.7 | 13.6 | –14.1 | –5.5 | 0.2 | 6.8 | 16.3 | |
6 | Case | 34.5 | 23.1 | 8.1 | 15.3 | 31.6 | 48.4 | 63.0 |
Control | 34.0 | 20.3 | 8.4 | 16.0 | 33.8 | 48.8 | 60.6 | |
Difference | 0.6 | 13.5 | –14.5 | –6.0 | -0.2 | 6.8 | 15.4 | |
7 | Case | 34.1 | 22.4 | 8.3 | 15.2 | 31.5 | 47.8 | 62.1 |
Control | 34.0 | 20.2 | 8.5 | 16.3 | 33.0 | 49.0 | 60.1 | |
Difference | 0.0 | 13.3 | –15.0 | –6.3 | –0.3 | 6.6 | 15.8 | |
8 | Case | 32.5 | 21.7 | 8.0 | 14.4 | 30.0 | 46.2 | 59.9 |
Control | 32.7 | 19.9 | 7.9 | 15.1 | 31.0 | 47.6 | 60.3 | |
Difference | –0.2 | 14.1 | –16.2 | –6.8 | –0.1 | 6.2 | 15.5 | |
9 | Case | 32.3 | 21.6 | 7.4 | 14.5 | 29.2 | 45.8 | 59.4 |
Control | 32.7 | 19.7 | 7.6 | 15.1 | 32.1 | 47.6 | 57.2 | |
Difference | –0.3 | 14.2 | –16.3 | –6.5 | –0.2 | 6.1 | 15.1 | |
10 | Case | 33.0 | 21.7 | 7.8 | 14.7 | 30.3 | 47.6 | 62.0 |
Control | 33.2 | 19.8 | 8.2 | 15.9 | 32.5 | 47.1 | 58.8 | |
Difference | –0.3 | 13.9 | –16.7 | –6.3 | 0.0 | 6.0 | 16.0 | |
11 | Case | 33.0 | 21.2 | 8.2 | 15.1 | 31.3 | 46.9 | 60.8 |
Control | 33.4 | 20.2 | 8.5 | 15.9 | 33.2 | 47.2 | 59.2 | |
Difference | –0.5 | 13.6 | –17.2 | –6.6 | –0.3 | 6.2 | 15.8 | |
12 | Case | 33.0 | 21.3 | 8.0 | 14.7 | 31.2 | 47.1 | 59.9 |
Control | 33.6 | 20.4 | 8.1 | 15.7 | 33.2 | 47.4 | 60.7 | |
Difference | –0.5 | 13.1 | –16.3 | –6.6 | –0.3 | 5.6 | 14.6 | |
13 | Case | 34.3 | 22.1 | 8.0 | 15.9 | 32.1 | 49.9 | 62.2 |
Control | 34.3 | 20.5 | 8.5 | 16.1 | 34.0 | 48.5 | 60.7 | |
Difference | 0.0 | 13.0 | –15.9 | –6.9 | –0.2 | 6.1 | 16.0 | |
* There were 914 case-days and 2,362 control days with complete information on NO2 concentration and other covariates. The difference in NO2 concentration was estimated from concentration measured on case-day minus the mean concentration measured on control-days occurring in the diary period. |
Table 3: Adjusted OR's* for doctor’s visit for respiratory symptoms per 10 µg/m3 increase in NO2. | ||||||||||||
Asthma | Bronchitis | COPD | Healthy | |||||||||
N doctor visits | 167 | 407 | 150 | 207 | ||||||||
OR** | 95%CI** | OR | 95%CI | OR | 95%CI | OR | 95%CI | |||||
Lag 0 | 1.019 | 0.849 | 1.223 | 1.116 | 1.001 | 1.244 | 1.098 | 0.908 | 1.328 | 0.978 | 0.832 | 1.148 |
Lag 1 | 0.953 | 0.878 | 1.036 | 1.047 | 0.996 | 1.101 | 1.134 | 1.029 | 1.250 | 0.924 | 0.852 | 1.001 |
Lag 2 | 0.892 | 0.789 | 1.009 | 0.983 | 0.913 | 1.059 | 1.171 | 1.017 | 1.348 | 0.873 | 0.769 | 0.990 |
Lag 3 | 0.937 | 0.828 | 1.060 | 0.981 | 0.914 | 1.053 | 1.146 | 1.002 | 1.311 | 0.919 | 0.822 | 1.027 |
Lag 4 | 1.067 | 0.966 | 1.178 | 1.033 | 0.970 | 1.099 | 1.073 | 0.969 | 1.188 | 1.027 | 0.941 | 1.122 |
Lag 5 | 1.152 | 1.024 | 1.295 | 1.071 | 0.992 | 1.156 | 1.020 | 0.896 | 1.160 | 1.075 | 0.959 | 1.205 |
Lag 6 | 1.135 | 1.023 | 1.258 | 1.067 | 1.004 | 1.134 | 1.006 | 0.896 | 1.131 | 1.035 | 0.945 | 1.134 |
Lag 7 | 1.048 | 0.945 | 1.162 | 1.027 | 0.967 | 1.091 | 1.028 | 0.929 | 1.138 | 0.965 | 0.883 | 1.055 |
Lag 8 | 0.962 | 0.858 | 1.079 | 0.975 | 0.905 | 1.051 | 1.064 | 0.954 | 1.188 | 0.933 | 0.835 | 1.043 |
Lag 9 | 0.930 | 0.853 | 1.013 | 0.939 | 0.884 | 0.997 | 1.087 | 0.988 | 1.196 | 0.978 | 0.898 | 1.066 |
Lag 10 | 0.970 | 0.868 | 1.085 | 0.934 | 0.871 | 1.002 | 1.070 | 0.953 | 1.200 | 1.085 | 0.975 | 1.208 |
Lag 11 | 1.051 | 0.932 | 1.185 | 0.959 | 0.889 | 1.035 | 1.004 | 0.893 | 1.128 | 1.107 | 0.983 | 1.248 |
Lag 12 | 1.119 | 1.036 | 1.208 | 0.994 | 0.948 | 1.042 | 0.921 | 0.838 | 1.012 | 0.996 | 0.928 | 1.069 |
Lag 13 | 1.191 | 1.007 | 1.409 | 1.029 | 0.930 | 1.139 | 0.845 | 0.700 | 1.021 | 0.896 | 0.768 | 1.046 |
* As estimated in conditional logistic regression at different lags and adjusted for flu prevalence, polynomial distributed lags of T, RH, AP, T2, RH2, and T*RH from lag 0–6 days up to 2 degrees of freedom, and average T, RH, AP, and T2 from lag 7–13 days. ** Bold font indicates p-values <0.05. |
We observed positive associations between ambient exposures to NO2 and visits to the doctor for respiratory symptoms on days following exposure, but the timing and magnitude in the effect of NO2varied by the underlying respiratory condition. Participants in the COPD subgroup were more vulnerable during the first week following an increase in exposure to NO2 compared with those in the asthma and healthy subgroups. The results for TSP were broadly similar to those observed for NO2.
Few studies have examined effects of ambient air pollutants on respiratory morbidity outcomes in adults in primary care rather than emergency care settings. Previous studies of air pollution and general practice consultations for asthma or lower respiratory tract diseases in London [5], and doctors’ house calls for upper or lower respiratory tract diseases in the greater Paris area [4] observed no significant associations with NO2. The former study did not consider exposure lags more than 3 days, but the latter study evaluated up to 15 days and also found that the effects of PM10 on house calls for lower respiratory diseases persisted until 4 days after exposure followed by a steady decrease.
The studies cited in the previous paragraph did not evaluate the susceptibility to respiratory health effects of air pollutants by pre-existing chronic respiratory disease. There is evidence from both epidemiologic and laboratory studies to suggest that adults with asthma and COPD are particularly sensitive to air pollutants using a wide range of outcomes that include aggravation of respiratory symptoms, emergency room visits, and hospital admissions [1, 2, 21, 22]. Evaluation of susceptibility factors may also help redefine exposure limits. The findings from our study indicate that the number of doctor visits for respiratory symptoms among individuals with COPD would increase by 9.1% over the course of one week following an increase of 10 µg/m3 of exposure to NO2. COPD is a common disease in most countries, with a reported prevalence in Europe and United States falling between 4 and 10% [23]. Thus, the present findings re-emphasise that air quality improvements would represent considerable medical and other societal savings, and should be included in general strategies for disease prevention [24].
The variation in lag distribution of effects between the disease subgroups is also of particular interest. The observed findings suggest that the time to doctor visit following increase in exposure to NO2 is more immediate for individuals who have at least mild COPD than for those with asthma. Both diseases are characterised by chronic inflammatory processes, but the general features are distinct with different pathophysiology. Airflow obstruction is poorly reversible for COPD in comparison with asthma, whose exacerbations can be managed more effectively with the use of bronchodilators and inhaled steroids. Thus, the variation in lag distribution of effects may reflect a comparison of disease severity between COPD and asthma. However, there may be alternative explanations for the observed findings. First, differences in lag time and associations with doctor visits observed between the subgroups may be attributable to differences in smoking and gender. Lagged effect patterns of NO2 exposure were similar for asthma and COPD when the analysis was limited to former and never smokers, but the effects were more immediate and higher for COPD than asthma in current smokers. Although it is known that there are gender differences in treatment seeking behaviour, the lagged effects of NO2 exposure on doctor visits were more immediate for COPD than for asthma in both males and females. Secondly, peak effects may have occurred at shorter exposure lags for participants with COPD and chronic bronchitis because they are less likely to take medication compared with participants with asthma. It is possible that the use of bronchodilators or inhaled steroids in asthma exacerbations may delay (or postpone) the decision to visit the doctor. However, effects of NO2 exposure did not differ considerably between asthmatics who reported use of asthma medication in the 1991 baseline survey and those who did not (data not shown). Molecular and genetic epidemiology may help to elucidate the biological mechanisms of susceptibility to respiratory health effects in air pollution, which has important links with systemic inflammation and oxidative stress, in individuals with COPD and asthma [25].
Furthermore, we observed a bimodal lag distribution of peak effects of NO2 and TSP in all subgroups, although there was variation in timing and magnitude by condition. For the COPD subgroup, the development of respiratory symptoms leading to a doctor visit on the second or third day after exposure to NO2 may reflect a more rapid mechanism of airway irritation in individuals with more severe obstruction, or with presence of other susceptibility factors including older age and co-morbidities. The second wave of doctor visits on 9–10 days following exposure may be due to delayed or slower effects of a similar mechanism in milder cases of obstruction but might also be an artefact. The presence of such a second peak would need to be confirmed with other data.
This study has multiple strengths and limitations that are inherent to its design. The subjects of the SAPALDIA diary study participated one to two years earlier in the baseline survey of SAPALDIA, which included information on spirometry and methacholine challenge testing, and respiratory symptoms. Many of the diary study participants were selected on the basis of presenting with respiratory symptoms or poor lung function at the baseline survey, which enabled us to make comparisons between subgroups [11, 12]. However, misclassification of subgroup is a possibility, particularly for COPD. While the prevalence of COPD in the younger age demographic of the general population is rare, 9 of the 54 participants (16.7%) in the COPD subgroup were less than 40 years of age, although the higher prevalence may reflect the selection of those who reported doctor visits. Bronchodilator tests were also not performed in the SAPALDIA baseline survey, so it was not possible to adhere to the more strict definition of COPD. Instead, validated questions for history of physician diagnosed asthma, current asthma, and asthma medication use in combination with lung function testing were used to minimise the inclusion of participants with potential asthma in the COPD subgroup. Furthermore, lag structures were similar between asthmatics who reported or didn’t report asthma medication use in the 1991 baseline survey (data not shown); if unrecognised and untreated asthmatics were excluded from the COPD subgroup, the lagged effect patterns between the COPD and asthma subgroups would be likely to differ even more. Finally, to ensure accurate classification of healthy participants, normal bronchial reactivity tests were required in addition to absence of asthma, bronchitis symptoms, and airflow obstruction. The distinct patterns of time-related effects of NO2 effects observed between the asthma, COPD, and healthy subgroups suggest that misclassification is likely to be minimal.
The case-crossover analysis, in which the same participant was both a case and a control, prevented confounding by gender, age, smoking status and other individual factors that do not vary within time strata. However, there may be unmeasured confounding by time-varying factors, including daily medication use, which was not available in this analysis. Moderate prevalence of asthma medication or inhalation steroids, according to data collected at the baseline survey of SAPALDIA, was only observed in the asthma subgroup. As described in the previous paragraph, lag structures were also similar between treated and untreated participants in the asthma subgroup. Thus, unmeasured confounding by daily medication use in asthmatic participants may be rather minimal.
While referent time periods are individually matched to case time periods in case-crossover studies, there are several different strategies to select referent periods, and the risk estimates may vary by selection strategy [26, 27]. Based on review of this literature, we chose our referent periods using a time-stratified design, in which control and case-days were matched by weekday within the diary period, thus providing sufficient control for time trend biases including seasonal trends [27].
The observed associations between daily levels of NO2 and TSP and reported doctor visits were statistically significant at specific lags, but were mostly of smaller magnitude with confidence intervals indicating wide variation. The CLR models were adjusted for time-varying factors including weekly flu prevalence at the regional level and complex combinations of meteorological variables. Effect estimates of similar magnitude were also observed in CLR models that were unadjusted for meteorological covariates (data not shown). Furthermore, these models were stratified by presence of chronic respiratory disease, gender, and smoking to account for the hypothesised heterogeneity of effects between subgroups.
The assignment of daily exposure was also based on fixed site monitors that have rather limited precision to indicate exposure of a particular participant on a given day, so there is potential for exposure misclassification of NO2 and TSP. However, relevant to TSP, the spatial homogeneity of particulate concentrations was demonstrated to be quite high in one of the regions for this study [28]. Finally, other potential causal components, including ultrafine particles, were not addressed and may explain the observed effects in this study. Ambient NO2 was and still is a proxy for combustion, and nearly 20 years ago, it was strongly related to traffic; while ultrafine particle data was not collected in this study, we may consider NO2 a proxy of traffic related PM.
In summary, an increase in exposure to NO2 and TSP was associated with a visit to the doctor for a few days following exposure. Respiratory problems leading to a doctor visit may have a faster dynamic in individuals with underlying COPD.
Supplementary file: Table 1 Characteristics* of subjects with asthma and COPD who report or did not doctor visits for respiratory symptoms in the SAPALDIA diary study. | ||||
Asthma | COPD | |||
Did not report doctor visit | Reported doctor visit | Did not report doctor visit | Reported doctor visit | |
Number of subjects | 251 | 78 | 296 | 54 |
n (%)** | n (%) | n (%) | n (%) | |
Sex | ||||
Male | 124 (49.4) | 31 (39.7) | 192 (64.9) | 29 (53.7) |
Female | 127 (50.6) | 47 (60.3) | 104 (35.1) | 25 (46.3) |
Age | ||||
<35 | 72 (28.7) | 17 (21.8) | 19 (6.4) | 7 (13.0) |
35–45 | 62 (24.7) | 24 (30.8) | 77 (26.0) | 10 (18.5) |
>45 | 117 (46.6) | 37 (47.4) | 200 (67.6) | 37 (68.5) |
mean (SD) | 42.3 (11.7) | 43.6 (10.1) | 49.0 (9.0) | 48.2 (9.9) |
Bronchitis symptoms | ||||
Yes | 92 (36.7) | 40 (51.3) | 107 (36.2) | 23 (42.6) |
No | 159 (63.4) | 38 (48.7) | 189 (63.9) | 31 (57.4) |
Smoking status | ||||
Current smoker | 60 (23.9) | 14 (17.3) | 104 (35.1) | 24 (46.0) |
Former or never smoker | 191 (76.1) | 64 (82.7) | 192 (64.9) | 30 (54.0) |
Current asthma medication use | ||||
Yes | 72 (28.7) | 35 (44.9) | ||
No | 179 (71.3) | 43 (55.1) | ||
Lung function | ||||
Percent predicted FEV1, mean (SD) | 97.8 (18.3) | 88.7 (21.9) | 95.2 (16.1) | 91.9 (16.2) |
* Characteristics summarised from the baseline survey of SAPALDIA (1991) ** Expressed as n (%) unless otherwise noted |
Supplementary file: Table 2 Adjusted OR's* for doctor’s visit for respiratory symptoms per 10 µg/m3 increase in NO2 after stratification by smoking status. | ||||||||||||
Current smokers | ||||||||||||
Asthma | Bronchitis | COPD | Healthy | |||||||||
N doctor visits | 44 | 196 | 75 | 25 | ||||||||
OR** | 95%CI** | OR | 95%CI | OR | 95%CI | OR | 95%CI | |||||
Lag 0 | 1.207 | 0.743 | 1.962 | 1.110 | 0.943 | 1.307 | 1.284 | 0.881 | 1.872 | N/A | ||
Lag 1 | 1.003 | 0.828 | 1.214 | 1.066 | 0.985 | 1.154 | 1.371 | 1.064 | 1.767 | |||
Lag 2 | 0.833 | 0.635 | 1.093 | 1.024 | 0.909 | 1.154 | 1.464 | 1.064 | 2.013 | |||
Lag 3 | 0.834 | 0.615 | 1.130 | 1.026 | 0.917 | 1.148 | 1.248 | 0.947 | 1.646 | |||
Lag 4 | 1.010 | 0.716 | 1.424 | 1.057 | 0.957 | 1.167 | 0.964 | 0.776 | 1.197 | |||
Lag 5 | 1.218 | 0.827 | 1.793 | 1.065 | 0.943 | 1.203 | 0.892 | 0.689 | 1.155 | |||
Lag 6 | 1.321 | 0.978 | 1.784 | 1.038 | 0.943 | 1.143 | 1.022 | 0.806 | 1.295 | |||
Lag 7 | 1.258 | 0.930 | 1.701 | 0.991 | 0.906 | 1.085 | 1.283 | 0.993 | 1.658 | |||
Lag 8 | 1.086 | 0.766 | 1.539 | 0.951 | 0.848 | 1.067 | 1.482 | 1.101 | 1.995 | |||
Lag 9 | 0.914 | 0.727 | 1.149 | 0.938 | 0.851 | 1.033 | 1.379 | 1.082 | 1.759 | |||
Lag 10 | 0.824 | 0.614 | 1.105 | 0.956 | 0.859 | 1.065 | 1.040 | 0.811 | 1.334 | |||
Lag 11 | 0.877 | 0.607 | 1.267 | 0.987 | 0.881 | 1.105 | 0.810 | 0.610 | 1.076 | |||
Lag 12 | 1.046 | 0.882 | 1.240 | 1.006 | 0.934 | 1.083 | 0.753 | 0.583 | 0.971 | |||
Lag 13 | 1.247 | 0.906 | 1.717 | 1.025 | 0.880 | 1.194 | 0.699 | 0.460 | 1.062 | |||
Former or never smokers | ||||||||||||
Asthma | Bronchitis | COPD | Healthy | |||||||||
N doctor visits | 123 | 210 | 75 | 182 | ||||||||
OR | 95%CI | OR | 95%CI | OR | 95%CI | OR | 95%CI | |||||
Lag 0 | 1.060 | 0.846 | 1.327 | 1.146 | 0.975 | 1.346 | 1.219 | 0.865 | 1.717 | 1.005 | 0.836 | 1.209 |
Lag 1 | 0.954 | 0.855 | 1.064 | 1.024 | 0.949 | 1.106 | 1.012 | 0.861 | 1.190 | 0.937 | 0.855 | 1.026 |
Lag 2 | 0.859 | 0.730 | 1.011 | 0.916 | 0.821 | 1.021 | 0.840 | 0.642 | 1.099 | 0.873 | 0.758 | 1.004 |
Lag 3 | 0.901 | 0.768 | 1.056 | 0.917 | 0.827 | 1.017 | 0.855 | 0.659 | 1.109 | 0.920 | 0.812 | 1.043 |
Lag 4 | 1.055 | 0.939 | 1.185 | 1.017 | 0.932 | 1.109 | 1.011 | 0.852 | 1.199 | 1.046 | 0.949 | 1.154 |
Lag 5 | 1.159 | 1.009 | 1.333 | 1.106 | 0.993 | 1.232 | 1.112 | 0.898 | 1.378 | 1.115 | 0.980 | 1.269 |
Lag 6 | 1.142 | 0.999 | 1.304 | 1.124 | 1.027 | 1.229 | 1.081 | 0.893 | 1.309 | 1.087 | 0.976 | 1.210 |
Lag 7 | 1.044 | 0.906 | 1.202 | 1.070 | 0.980 | 1.169 | 0.979 | 0.832 | 1.151 | 1.015 | 0.917 | 1.124 |
Lag 8 | 0.955 | 0.822 | 1.109 | 0.989 | 0.891 | 1.098 | 0.908 | 0.758 | 1.088 | 0.974 | 0.861 | 1.102 |
Lag 9 | 0.937 | 0.841 | 1.044 | 0.928 | 0.856 | 1.006 | 0.930 | 0.797 | 1.084 | 1.006 | 0.913 | 1.110 |
Lag 10 | 1.006 | 0.873 | 1.160 | 0.915 | 0.828 | 1.010 | 1.039 | 0.857 | 1.259 | 1.095 | 0.968 | 1.238 |
Lag 11 | 1.098 | 0.944 | 1.277 | 0.951 | 0.853 | 1.060 | 1.082 | 0.882 | 1.327 | 1.091 | 0.954 | 1.247 |
Lag 12 | 1.143 | 1.031 | 1.266 | 1.005 | 0.939 | 1.075 | 0.982 | 0.828 | 1.164 | 0.956 | 0.880 | 1.039 |
Lag 13 | 1.189 | 0.948 | 1.490 | 1.062 | 0.918 | 1.229 | 0.890 | 0.636 | 1.247 | 0.838 | 0.701 | 1.003 |
* As estimated in conditional logistic regression per 10 µg/m3 increase in NO2 at different lags and adjusted for flu prevalence, polynomial distributed lags of T, RH, AP, T2, RH2, and T*RH from lag 0-6 days up to 2 degrees of freedom, and average T, RH, AP, and T2 from lag 7–13 days. ** Bold font indicates p-values <0.05. |
Supplementary file: Table 3 Adjusted OR's* for doctor’s visit for respiratory symptoms per 10 µg/m3 increase in NO2 after stratification by gender. | ||||||||||||
Males | ||||||||||||
Asthma | Bronchitis | COPD | Healthy | |||||||||
N doctor visits | 81 | 168 | 93 | 59 | ||||||||
OR** | 95%CI** | OR | 95%CI | OR | 95%CI | OR | 95%CI | |||||
Lag 0 | 1.106 | 0.803 | 1.524 | 1.272 | 1.055 | 1.535 | 1.053 | 0.824 | 1.347 | 1.051 | 0.646 | 1.710 |
Lag 1 | 0.985 | 0.860 | 1.128 | 1.099 | 1.008 | 1.197 | 1.129 | 0.987 | 1.292 | 0.792 | 0.589 | 1.064 |
Lag 2 | 0.878 | 0.714 | 1.079 | 0.949 | 0.826 | 1.090 | 1.210 | 1.003 | 1.459 | 0.597 | 0.374 | 0.954 |
Lag 3 | 0.957 | 0.771 | 1.187 | 0.986 | 0.867 | 1.122 | 1.247 | 1.041 | 1.494 | 0.740 | 0.517 | 1.057 |
Lag 4 | 1.166 | 0.972 | 1.399 | 1.164 | 1.041 | 1.302 | 1.209 | 1.037 | 1.410 | 1.197 | 0.936 | 1.530 |
Lag 5 | 1.252 | 1.016 | 1.541 | 1.263 | 1.096 | 1.456 | 1.135 | 0.942 | 1.368 | 1.410 | 0.998 | 1.990 |
Lag 6 | 1.139 | 0.964 | 1.347 | 1.192 | 1.070 | 1.329 | 1.065 | 0.903 | 1.256 | 1.124 | 0.903 | 1.398 |
Lag 7 | 0.960 | 0.821 | 1.123 | 1.030 | 0.934 | 1.135 | 1.025 | 0.894 | 1.176 | 0.765 | 0.593 | 0.986 |
Lag 8 | 0.855 | 0.710 | 1.029 | 0.894 | 0.783 | 1.021 | 1.026 | 0.885 | 1.188 | 0.620 | 0.425 | 0.906 |
Lag 9 | 0.891 | 0.774 | 1.027 | 0.852 | 0.760 | 0.954 | 1.058 | 0.925 | 1.210 | 0.754 | 0.582 | 0.979 |
Lag 10 | 1.079 | 0.896 | 1.299 | 0.909 | 0.805 | 1.028 | 1.086 | 0.928 | 1.270 | 1.270 | 0.879 | 1.833 |
Lag 11 | 1.247 | 1.003 | 1.551 | 0.985 | 0.863 | 1.123 | 1.032 | 0.887 | 1.200 | 1.633 | 1.053 | 2.532 |
Lag 12 | 1.241 | 1.061 | 1.450 | 0.996 | 0.914 | 1.086 | 0.912 | 0.799 | 1.040 | 1.272 | 0.995 | 1.625 |
Lag 13 | 1.234 | 0.933 | 1.632 | 1.008 | 0.845 | 1.204 | 0.806 | 0.616 | 1.054 | 0.990 | 0.563 | 1.742 |
Females | ||||||||||||
Asthma | Bronchitis | COPD | Healthy | |||||||||
N doctor visits | 86 | 239 | 57 | 148 | ||||||||
OR | 95%CI | OR | 95%CI | OR | 95%CI | OR | 95%CI | |||||
Lag 0 | 0.913 | 0.699 | 1.191 | 1.038 | 0.894 | 1.206 | 2.330 | 0.989 | 5.486 | 0.966 | 0.795 | 1.172 |
Lag 1 | 0.902 | 0.788 | 1.032 | 1.021 | 0.952 | 1.096 | 1.958 | 1.215 | 3.156 | 0.933 | 0.845 | 1.031 |
Lag 2 | 0.891 | 0.738 | 1.075 | 1.005 | 0.911 | 1.108 | 1.646 | 0.947 | 2.864 | 0.902 | 0.775 | 1.049 |
Lag 3 | 0.957 | 0.797 | 1.149 | 0.986 | 0.896 | 1.086 | 1.307 | 0.814 | 2.098 | 0.935 | 0.819 | 1.068 |
Lag 4 | 1.091 | 0.941 | 1.264 | 0.977 | 0.899 | 1.061 | 1.051 | 0.806 | 1.370 | 1.018 | 0.916 | 1.130 |
Lag 5 | 1.191 | 1.002 | 1.415 | 0.982 | 0.890 | 1.083 | 0.948 | 0.638 | 1.408 | 1.071 | 0.932 | 1.231 |
Lag 6 | 1.198 | 1.023 | 1.403 | 0.995 | 0.917 | 1.080 | 0.970 | 0.689 | 1.366 | 1.069 | 0.954 | 1.199 |
Lag 7 | 1.122 | 0.949 | 1.326 | 1.004 | 0.923 | 1.091 | 1.074 | 0.801 | 1.441 | 1.034 | 0.929 | 1.150 |
Lag 8 | 1.018 | 0.851 | 1.217 | 0.998 | 0.904 | 1.101 | 1.194 | 0.823 | 1.734 | 1.002 | 0.878 | 1.143 |
Lag 9 | 0.941 | 0.826 | 1.072 | 0.974 | 0.901 | 1.053 | 1.242 | 0.917 | 1.682 | 1.000 | 0.899 | 1.112 |
Lag 10 | 0.928 | 0.775 | 1.110 | 0.945 | 0.859 | 1.040 | 1.168 | 0.832 | 1.639 | 1.024 | 0.902 | 1.161 |
Lag 11 | 0.994 | 0.825 | 1.198 | 0.942 | 0.852 | 1.043 | 1.028 | 0.694 | 1.522 | 1.012 | 0.881 | 1.164 |
Lag 12 | 1.107 | 0.995 | 1.231 | 0.970 | 0.909 | 1.036 | 0.907 | 0.687 | 1.197 | 0.947 | 0.869 | 1.033 |
Lag 13 | 1.233 | 0.952 | 1.598 | 0.999 | 0.868 | 1.149 | 0.800 | 0.466 | 1.374 | 0.886 | 0.742 | 1.059 |
* As estimated in conditional logistic regression per 10 µg/m3 increase in NO2 at different lags and adjusted for flu prevalence, polynomial distributed lags of T, RH, AP, T2, RH2, and T*RH from lag 0-6 days up to 2 degrees of freedom, and average T, RH, AP, and T2 from lag 7–13 days. ** Bold font indicates p-values <0.05. |
Supplementary file: Table 4 Distribution of 24-mean concentrations of TSP (µg/m3) for case and control-days. | ||||||||
Lag | Case/Control* | Mean | SD | 10th | 25th | 50th | 75th | 90th |
Case-days | 35.7 | 27.2 | 10.9 | 17.0 | 27.3 | 45.3 | 72.0 | |
Control-days | 33.9 | 21.5 | 12.7 | 17.7 | 12.7 | 44.5 | 62.0 | |
1 | Case-days | 34.3 | 26.0 | 10.0 | 15.6 | 26.4 | 45.0 | 70.0 |
Control-days | 33.4 | 20.3 | 12.3 | 18.5 | 28.5 | 43.8 | 62.7 | |
2 | Case-days | 34.2 | 26.0 | 9.8 | 14.8 | 26.0 | 47.0 | 73.0 |
Control-days | 32.7 | 20.0 | 12.3 | 17.5 | 28.0 | 42.4 | 63.0 | |
3 | Case-days | 34.9 | 26.7 | 10.4 | 15.8 | 25.9 | 48.0 | 72.0 |
Control-days | 33.1 | 20.0 | 12.0 | 17.3 | 28.7 | 44.1 | 61.5 | |
4 | Case-days | 35.4 | 27.1 | 10.0 | 16.2 | 27.0 | 48.0 | 72.0 |
Control-days | 33.1 | 20.1 | 12.0 | 17.3 | 28.5 | 44.3 | 62.7 | |
5 | Case-days | 35.7 | 27.6 | 9.8 | 15.6 | 28.0 | 46.2 | 74.1 |
Control-days | 33.9 | 21.0 | 12.2 | 17.3 | 29.0 | 46.0 | 65.6 | |
6 | Case-days | 36.2 | 28.4 | 10.4 | 16.0 | 28.0 | 47.0 | 75.0 |
Control-days | 34.7 | 21.2 | 12.5 | 18.7 | 30.4 | 45.0 | 63.8 | |
7 | Case-days | 35.3 | 28.1 | 10.0 | 16.0 | 27.0 | 45.0 | 72.0 |
Control-days | 34.9 | 21.5 | 12.5 | 18.7 | 31.0 | 45.3 | 64.0 | |
8 | Case-days | 33.9 | 25.9 | 9.4 | 16.0 | 26.5 | 45.0 | 68.0 |
Control-days | 33.8 | 20.5 | 12.5 | 18.1 | 29.4 | 44.8 | 62.3 | |
9 | Case-days | 33.4 | 25.2 | 9.8 | 15.6 | 26.0 | 44.0 | 66.0 |
Control-days | 33.7 | 20.7 | 12.5 | 17.9 | 28.5 | 46.0 | 61.0 | |
10 | Case-days | 33.9 | 25.4 | 9.7 | 15.4 | 27.2 | 45.8 | 67.0 |
Control-days | 34.0 | 21.1 | 12.5 | 17.9 | 28.6 | 46.0 | 64.5 | |
11 | Case-days | 34.0 | 25.1 | 10.0 | 15.0 | 27.0 | 46.0 | 69.2 |
Control-days | 34.5 | 21.0 | 13.0 | 18.1 | 28.3 | 46.1 | 64.3 | |
12 | Case-days | 33.9 | 26.8 | 10.0 | 15.0 | 26.2 | 42.9 | 70.0 |
Control-days | 35.2 | 21.7 | 12.8 | 18.7 | 29.3 | 47.3 | 65.3 | |
13 | Case-days | 34.8 | 27.1 | 9.8 | 15.5 | 27.0 | 45.0 | 69.7 |
Control-days | 35.7 | 22.0 | 12.7 | 19.1 | 30.3 | 47.7 | 66.0 | |
* There were 892 case-days and 2,230 control days with complete information on TSP concentration and other covariates. |
Supplementary file: Table 5Adjusted OR's* for doctor's visit for respiratory symptoms per 10 µg/m3 increase in TSP at separate lags in major subcategories of participants as estimated in conditional logistic regression. | ||||||||||||
All subjects | ||||||||||||
Asthma diagnosis | Bronchitis | COPD | Healthy | |||||||||
n visits | 170 | 414 | 148 | 181 | ||||||||
OR | 95%CI | OR | 95%CI | OR | 95%CI | OR | 95%CI | |||||
Lag 0 | 0.977 | 0.875 | 1.091 | 1.016 | 0.957 | 1.079 | 1.001 | 0.894 | 1.120 | 0.966 | 0.885 | 1.054 |
Lag 1 | 0.979 | 0.930 | 1.031 | 1.013 | 0.984 | 1.042 | 1.032 | 0.976 | 1.092 | 0.966 | 0.930 | 1.003 |
Lag 2 | 0.981 | 0.912 | 1.057 | 1.009 | 0.965 | 1.055 | 1.065 | 0.984 | 1.152 | 0.966 | 0.901 | 1.035 |
Lag 3 | 1.020 | 0.947 | 1.099 | 1.011 | 0.970 | 1.054 | 1.051 | 0.978 | 1.131 | 0.991 | 0.928 | 1.059 |
Lag 4 | 1.068 | 1.003 | 1.137 | 1.013 | 0.976 | 1.052 | 1.006 | 0.944 | 1.072 | 1.025 | 0.974 | 1.079 |
Lag 5 | 1.073 | 1.002 | 1.150 | 1.009 | 0.964 | 1.055 | 0.981 | 0.905 | 1.063 | 1.033 | 0.970 | 1.101 |
Lag 6 | 1.036 | 0.981 | 1.093 | 0.999 | 0.965 | 1.035 | 0.986 | 0.924 | 1.052 | 1.015 | 0.964 | 1.068 |
Lag 7 | 0.988 | 0.932 | 1.047 | 0.991 | 0.956 | 1.027 | 1.013 | 0.959 | 1.070 | 0.990 | 0.944 | 1.039 |
Lag 8 | 0.965 | 0.899 | 1.035 | 0.989 | 0.945 | 1.035 | 1.040 | 0.972 | 1.111 | 0.983 | 0.925 | 1.045 |
Lag 9 | 0.985 | 0.932 | 1.041 | 0.995 | 0.960 | 1.032 | 1.045 | 0.985 | 1.108 | 1.007 | 0.954 | 1.062 |
Lag 10 | 1.037 | 0.967 | 1.113 | 1.004 | 0.962 | 1.048 | 1.022 | 0.947 | 1.102 | 1.051 | 0.989 | 1.117 |
Lag 11 | 1.049 | 0.974 | 1.130 | 0.997 | 0.952 | 1.045 | 0.990 | 0.917 | 1.068 | 1.065 | 1.001 | 1.133 |
Lag 12 | 0.996 | 0.952 | 1.042 | 0.974 | 0.946 | 1.002 | 0.969 | 0.920 | 1.021 | 1.031 | 0.990 | 1.074 |
Lag 13 | 0.946 | 0.855 | 1.046 | 0.950 | 0.895 | 1.009 | 0.949 | 0.843 | 1.068 | 0.999 | 0.914 | 1.091 |
* All models were adjusted for flu prevalence, polynomial distributed lags of T, RH, AP, T2, RH2, and T*RH from lag 0–6 days up to 2 degrees of freedom, and average T, RH, AP, and T2 from lag 7–13 days. ** Bold font indicates p-values <0.05. |
Supplementary file: Table 6 Adjusted OR's* for doctor's visit per 10 µg/m3 increase in TSP at separate lags in major subcategories of participants after stratification by smoking status as estimated in conditional logistic regression. | |||||||||||||||||||||
Current smokers | Former or never smokers | ||||||||||||||||||||
Asthma | Bronchitis | COPD | Asthma | Bronchitis | COPD | Healthy | |||||||||||||||
n visits | 48 | 198 | 78 | 122 | 215 | 70 | 153 | ||||||||||||||
OR** | 95%CI** | OR | 95%CI | OR | 95%CI | OR | 95%CI | OR | 95%CI | OR | 95%CI | OR | 95%CI | ||||||||
Lag 0 | 1.062 | 0.749 | 1.506 | 1.030 | 0.941 | 1.126 | 1.043 | 0.871 | 1.250 | 0.974 | 0.853 | 1.111 | 1.003 | 0.916 | 1.098 | 1.002 | 0.808 | 1.242 | 0.910 | 0.808 | 1.024 |
Lag 1 | 0.934 | 0.787 | 1.109 | 1.025 | 0.980 | 1.073 | 1.089 | 0.984 | 1.205 | 0.998 | 0.937 | 1.062 | 1.003 | 0.959 | 1.048 | 0.981 | 0.880 | 1.093 | 0.944 | 0.899 | 0.991 |
Lag 2 | 0.821 | 0.658 | 1.025 | 1.021 | 0.953 | 1.093 | 1.137 | 0.996 | 1.298 | 1.022 | 0.931 | 1.122 | 1.002 | 0.937 | 1.073 | 0.960 | 0.813 | 1.133 | 0.979 | 0.899 | 1.066 |
Lag 3 | 0.847 | 0.680 | 1.056 | 1.017 | 0.955 | 1.082 | 1.071 | 0.956 | 1.200 | 1.050 | 0.955 | 1.154 | 1.012 | 0.950 | 1.078 | 0.967 | 0.826 | 1.133 | 1.022 | 0.940 | 1.110 |
Lag 4 | 0.993 | 0.799 | 1.234 | 1.013 | 0.956 | 1.074 | 0.953 | 0.860 | 1.056 | 1.058 | 0.981 | 1.141 | 1.021 | 0.967 | 1.077 | 1.004 | 0.880 | 1.147 | 1.051 | 0.986 | 1.119 |
Lag 5 | 1.113 | 0.876 | 1.416 | 1.009 | 0.940 | 1.082 | 0.906 | 0.789 | 1.040 | 1.036 | 0.954 | 1.125 | 1.016 | 0.951 | 1.084 | 1.043 | 0.889 | 1.224 | 1.049 | 0.970 | 1.135 |
Lag 6 | 1.127 | 0.946 | 1.343 | 1.001 | 0.950 | 1.054 | 0.939 | 0.838 | 1.052 | 0.999 | 0.933 | 1.069 | 1.002 | 0.951 | 1.056 | 1.069 | 0.939 | 1.217 | 1.026 | 0.964 | 1.092 |
Lag 7 | 1.053 | 0.870 | 1.275 | 0.989 | 0.938 | 1.042 | 1.018 | 0.930 | 1.115 | 0.973 | 0.902 | 1.049 | 0.995 | 0.942 | 1.050 | 1.079 | 0.967 | 1.203 | 1.003 | 0.945 | 1.064 |
Lag 8 | 0.965 | 0.753 | 1.238 | 0.973 | 0.909 | 1.042 | 1.081 | 0.970 | 1.205 | 0.978 | 0.894 | 1.070 | 1.004 | 0.940 | 1.073 | 1.078 | 0.943 | 1.232 | 0.998 | 0.924 | 1.079 |
Lag 9 | 0.922 | 0.759 | 1.121 | 0.958 | 0.906 | 1.012 | 1.068 | 0.970 | 1.177 | 1.019 | 0.951 | 1.092 | 1.032 | 0.979 | 1.088 | 1.075 | 0.944 | 1.223 | 1.022 | 0.955 | 1.094 |
Lag 10 | 0.940 | 0.762 | 1.159 | 0.951 | 0.889 | 1.017 | 0.985 | 0.875 | 1.110 | 1.073 | 0.989 | 1.165 | 1.059 | 0.994 | 1.128 | 1.076 | 0.915 | 1.266 | 1.059 | 0.981 | 1.143 |
Lag 11 | 0.967 | 0.769 | 1.217 | 0.964 | 0.897 | 1.036 | 0.934 | 0.831 | 1.051 | 1.062 | 0.973 | 1.160 | 1.032 | 0.964 | 1.106 | 1.084 | 0.923 | 1.273 | 1.059 | 0.980 | 1.145 |
Lag 12 | 0.957 | 0.836 | 1.095 | 0.994 | 0.955 | 1.036 | 0.956 | 0.876 | 1.044 | 0.982 | 0.929 | 1.039 | 0.955 | 0.912 | 1.001 | 1.092 | 0.959 | 1.244 | 1.012 | 0.966 | 1.061 |
Lag 13 | 0.946 | 0.702 | 1.275 | 1.026 | 0.938 | 1.122 | 0.979 | 0.805 | 1.190 | 0.908 | 0.806 | 1.023 | 0.884 | 0.802 | 0.974 | 1.101 | 0.853 | 1.421 | 0.968 | 0.870 | 1.076 |
* All models were adjusted for flu prevalence, polynomial distributed lags of T, RH, AP, T2, RH2, and T*RH from lag 0–6 days up to 2 degrees of freedom, and average T, RH, AP, and T2 from lag 7–13 days. ** Bold font indicates p-values <0.05. |
Supplementary file: Table 7Adjusted OR's* for doctor's visit per 10 µg/m3 increase in TSP at separate lags in major subcategories of participants after stratification by gender as estimated in conditional logistic regression. | ||||||||||||
Males | ||||||||||||
Asthma | Bronchitis | COPD | Healthy | |||||||||
n visits | 81 | 173 | 92 | 42 | ||||||||
OR** | 95%CI** | OR | 95%CI | OR | 95%CI | OR | 95%CI | |||||
Lag 0 | 1.140 | 0.929 | 1.399 | 1.014 | 0.896 | 1.147 | 1.016 | 0.882 | 1.170 | 0.764 | 0.486 | 1.202 |
Lag 1 | 1.052 | 0.956 | 1.159 | 1.020 | 0.961 | 1.082 | 1.032 | 0.960 | 1.109 | 0.896 | 0.753 | 1.067 |
Lag 2 | 0.971 | 0.850 | 1.111 | 1.025 | 0.951 | 1.105 | 1.048 | 0.949 | 1.158 | 1.051 | 0.715 | 1.545 |
Lag 3 | 0.986 | 0.864 | 1.126 | 1.044 | 0.971 | 1.123 | 1.057 | 0.964 | 1.158 | 1.331 | 0.885 | 2.002 |
Lag 4 | 1.061 | 0.944 | 1.193 | 1.056 | 0.988 | 1.129 | 1.048 | 0.965 | 1.138 | 1.505 | 1.089 | 2.080 |
Lag 5 | 1.087 | 0.946 | 1.250 | 1.041 | 0.962 | 1.126 | 1.028 | 0.927 | 1.139 | 1.303 | 0.932 | 1.824 |
Lag 6 | 1.043 | 0.936 | 1.161 | 1.005 | 0.946 | 1.066 | 1.008 | 0.929 | 1.094 | 0.942 | 0.728 | 1.219 |
Lag 7 | 0.972 | 0.883 | 1.070 | 0.969 | 0.916 | 1.025 | 1.001 | 0.934 | 1.073 | 0.698 | 0.541 | 0.901 |
Lag 8 | 0.931 | 0.822 | 1.056 | 0.953 | 0.885 | 1.026 | 1.013 | 0.932 | 1.102 | 0.653 | 0.477 | 0.894 |
Lag 9 | 0.954 | 0.859 | 1.060 | 0.966 | 0.909 | 1.027 | 1.041 | 0.966 | 1.121 | 0.847 | 0.678 | 1.059 |
Lag 10 | 1.030 | 0.926 | 1.145 | 0.996 | 0.930 | 1.068 | 1.060 | 0.960 | 1.172 | 1.328 | 1.031 | 1.712 |
Lag 11 | 1.054 | 0.940 | 1.183 | 0.995 | 0.926 | 1.070 | 1.025 | 0.931 | 1.128 | 1.544 | 1.139 | 2.093 |
Lag 12 | 0.986 | 0.905 | 1.074 | 0.952 | 0.904 | 1.003 | 0.943 | 0.882 | 1.009 | 1.175 | 0.945 | 1.461 |
Lag 13 | 0.922 | 0.773 | 1.099 | 0.910 | 0.814 | 1.019 | 0.869 | 0.739 | 1.020 | 0.894 | 0.572 | 1.400 |
Females | ||||||||||||
Asthma | Bronchitis | COPD | Healthy | |||||||||
n visits | 89 | 241 | 56 | 139 | ||||||||
OR | 95%CI | OR | 95%CI | OR | 95%CI | OR | 95%CI | |||||
Lag 0 | 0.947 | 0.803 | 1.116 | 1.041 | 0.959 | 1.130 | 1.038 | 0.588 | 1.832 | 0.983 | 0.868 | 1.114 |
Lag 1 | 0.971 | 0.900 | 1.047 | 1.017 | 0.976 | 1.061 | 1.315 | 0.943 | 1.833 | 0.973 | 0.919 | 1.030 |
Lag 2 | 0.995 | 0.887 | 1.116 | 0.994 | 0.935 | 1.058 | 1.667 | 1.031 | 2.694 | 0.962 | 0.879 | 1.054 |
Lag 3 | 1.029 | 0.918 | 1.154 | 0.985 | 0.930 | 1.044 | 1.350 | 0.948 | 1.923 | 0.975 | 0.893 | 1.066 |
Lag 4 | 1.065 | 0.966 | 1.175 | 0.989 | 0.939 | 1.043 | 0.870 | 0.668 | 1.133 | 1.012 | 0.949 | 1.080 |
Lag 5 | 1.084 | 0.971 | 1.210 | 0.996 | 0.937 | 1.059 | 0.761 | 0.519 | 1.116 | 1.047 | 0.967 | 1.133 |
Lag 6 | 1.078 | 0.981 | 1.185 | 1.001 | 0.954 | 1.050 | 0.968 | 0.754 | 1.242 | 1.063 | 0.993 | 1.139 |
Lag 7 | 1.052 | 0.948 | 1.169 | 1.004 | 0.954 | 1.056 | 1.443 | 1.027 | 2.027 | 1.057 | 0.993 | 1.125 |
Lag 8 | 1.018 | 0.906 | 1.143 | 1.008 | 0.946 | 1.073 | 1.850 | 1.081 | 3.167 | 1.032 | 0.958 | 1.111 |
Lag 9 | 0.987 | 0.908 | 1.073 | 1.012 | 0.964 | 1.063 | 1.640 | 1.051 | 2.558 | 1.002 | 0.938 | 1.070 |
Lag 10 | 0.972 | 0.858 | 1.101 | 1.015 | 0.956 | 1.077 | 1.067 | 0.815 | 1.397 | 0.982 | 0.907 | 1.064 |
Lag 11 | 0.982 | 0.860 | 1.122 | 1.005 | 0.940 | 1.073 | 0.861 | 0.632 | 1.172 | 0.992 | 0.916 | 1.074 |
Lag 12 | 1.007 | 0.941 | 1.078 | 0.982 | 0.944 | 1.022 | 1.073 | 0.871 | 1.321 | 1.023 | 0.976 | 1.073 |
Lag 13 | 1.032 | 0.877 | 1.215 | 0.960 | 0.887 | 1.040 | 1.337 | 0.848 | 2.107 | 1.056 | 0.949 | 1.174 |
* All models were adjusted for flu prevalence, polynomial distributed lags of T, RH, AP, T2, RH2, and T*RH from lag 0–6 days up to 2 degrees of freedom, and average T, RH, AP, and T2 from lag 7–13 days. ** Bold font indicates p-values <0.05. |
The SAPALDIA team:
Study directorate: T. Rochat (p), J. M. Gaspoz (c), N. Künzli (e/exp), N. M. Probst Hensch (e/g), C. Schindler (s).
Scientific team: U. Ackermann-Liebrich (e), J. C. Barthélémy (c), W. Berger (g), R. Bettschart (p), A. Bircher (a), G. Bolognini (p), O. Brändli (p), C. Brombach (n), M. Brutsche (p), L. Burdet (p), M. Frey (p), U. Frey (pd), M. W. Gerbase (p), D. Gold (e/c/p), E. de Groot (c), W. Karrer (p), R. Keller (p), B. Knöpfli (p), B. Martin (pa), D. Miedinger (o), U. Neu (exp), L. Nicod (p), M. Pons (p), F. Roche (c), T. Rothe (p), E. Russi (p), P. Schmid-Grendelmeyer (a), M. Tamm (P), A. Schmidt-Trucksäss (pa), A. Turk (p), J. Schwartz (e), D. Stolz (p), P. Straehl (exp), J.-M. Tschopp (p), A. von Eckardstein (cc), J. P. Zellweger (p), E Zemp Stutz (e).
Scientific team at coordinating centres: M. Adam (e/g), E. Boes (g), P. O. Bridevaux (p), D. Carballo (c), E. Corradi (e), I. Curjuric (e), J. Dratva (e), A. Di Pasquale (s), L. Grize (s), D. Keidel (s), S. Kriemler (pa), A. Kumar (g), M. Imboden (g), N. Maire (s), A. J. Mehta (e), F. Meier (e), H. Phuleria (exp), E. Schaffner (s), G. A. Thun (g) A. Ineichen (exp), M. Ragettli (e), M. Ritter (exp), T. Schikowski (e), G. Stern (pd), M. Tarantino (s), M. Tsai (e), M. Wanner (pa).
(a) allergology; (c) cardiology; (cc) clinical chemistry; (e) epidemiology; (exp) exposure; (g) genetic and molecular biology; (m) meteorology; (n) nutrition; (o) occupational health; (p) pneumology; (pa) physical activity; (pd) paediatrics; (s) statistics.
Local fieldworkers: Aarau: S. Brun, G. Giger, M. Sperisen, M. Stahel; Basel: C. Bürli, C. Dahler, N. Oertli, I. Harreh, F. Karrer, G. Novicic, N. Wyttenbacher; Davos: A. Saner, P. Senn, R. Winzeler; Geneva: F. Bonfils, B. Blicharz, C. Landolt, J. Rochat; Lugano: S. Boccia, E. Gehrig, M.-T. Mandia, G. Solari, B. Viscardi, Montana: A. P. Bieri, C. Darioly, M. Maire, Payerne: F. Ding, P. Danieli A. Vonnez; Wald: D. Bodmer, E. Hochstrasser, R. Kunz, C. Meier, J. Rakic, U. Schafroth, A. Walder.
Administrative staff: C. Gabriel, R. Gutknecht.
Acknowledgement: We would like to thank the whole SAPALDIA team for their contribution to the study. Additionally, the study could not have been done without the help of the study participants, technical and administrative support and the medical teams and field workers at the local study sites.
1 Dominici F, Peng RD, Bell ML, Pham L, McDermott A, Zeger SL, Samet JM. Fine particulate air pollution and hospital admission for cardiovascular and respiratory diseases. JAMA. 2006;295:1127–34.
2 Koenig JQ, Air pollution and asthma. J Allergy Clin Immunol. 1999;104:717–22.
3 Sunyer J, Schwartz J, Tobías A, Macfarlane D, Garcia J, Antó JM, Patients with chronic obstructive pulmonary disease are at increased risk of death associated with urban particle air pollution: a case-crossover analysis. Am J Epidemiol. 2000;151:50–6.
4 Perez L, Rapp R, Künzli N. The Year of the Lung: outdoor air pollution and lung health. Swiss Med Wkly. 2010;140:w13129.
5 Chardon B, Lefranc A, Granados D, Grémy I. Air pollution and doctors’ house calls for respiratory diseases in the Greater Paris area (2000-3). Occup Environ Med. 2007;64:320–4.
6 Hajat S, Haines A, Goubet SA, Anderson H. Association of air pollution with daily GP consultations for asthma and other lower respiratory conditions in London. Thorax. 1999;54:597–605.
7 Hwang JS, Chan CC. Effects of air pollution on daily clinic visits for lower respiratory tract illness. Am J Epidemiol. 2002;155:1–10.
8 Wong TW, Tam W, Tak Sun Yu I, Wun YT, Wong AH, Wong CM. Association between air pollution and general practitioner visits for respiratory diseases in Hong Kong. Thorax. 2006;61:585–91.
9 Martin BW, Ackermann-Liebrich U, Leuenberger P, Künzli N, Stutz EZ, Keller R, et al. SAPALDIA: methods and participation in the cross-sectional part of the Swiss Study on Air Pollution and Lung Diseases in Adults. Soz Praventivmed. 1997;42:67–84.
10 Ackermann-Liebrich U, Kuna-Dibbert B, Probst-Hensch NM, Schindler C, Felber Dietrich D, Stutz EZ, et al. Follow-up of the Swiss Cohort Study on Air Pollution and Lung Diseases in Adults (SAPALDIA 2) 1991–2003: methods and characterization of participants. Soz Praventivmed. 2005;50:245–63.
11 Kaiser R, Schindler C, Künzli N, Ackermann-Liebrich U, Heeb D, Medici TC, et al. Use of transition probabilities to estimate the effect of smoking on the duration of episodes of respiratory symptoms in diary data: the Swiss Study on Air Pollution and Lung Diseases in Adults (SAPALDIA). Am J Epidemiol. 1998;148:600–8.
12 Künzli N, Stutz EZ, Perruchoud AP, Brändli O, Tschopp JM, Bolognini G, et al. Peak flow variability in the SAPALDIA study and its validity in screening for asthma-related conditions. Am J Respir Crit Care Med. 1999;160:427–34.
13 Lumley T, Levy D. Bias in the case-crossover design: implications for studies of air pollution. Environmetrics. 2000;11:689–704.
14 Monn C, Brändli O, Schindler C, Ackermann-Liebrich U, Leuenberger P. Personal exposure to nitrogen dioxide in Switzerland. Sci Total Environ. 1998;215:243–51.
15 Monn C, Braendli O, Schaeppi G, Schindler C, Ackermann-Liebrich U, Leuenberger P. Particulate matter <10 μm (PM10) and total suspended particulates (TSP) in urban, rural and alpine air in Switzerland. Atmos Environ. 1995;29:2565–73.
16 Hueglin C, Buchmann B, Weber RO. Long-term observation of real-world road traffic emission factors on a motorway in Switzerland. Atmos Environ. 2006;40:3696–709.
17 Pope CA 3rd, Schwartz J. Time series for the analysis of pulmonary health data. Am J Respir Crit Care Med. 1996;154:S229–33.
18 Schwartz J. The distributed lag between air pollution and daily deaths. Epidemiol. 2000;11:320–6.
19 Zanobetti A, Schwartz J, Samoli E, Gryparis A, Touloumi G, Peacock J, et al. The temporal pattern of respiratory and heart disease mortality in response to air pollution. Environ Health Perspect. 2003;111:1188–93.
20 Peel JL, Tolbert PE, Klein M, Metzger KB, Flanders WD, Todd K, et al. Ambient air pollution and respiratory emergency department visits. Epidemiology. 2005;16:164–74.
21 Mar TF, Larson TV, Stier RA, Claiborn C, Koenig JQ. An analysis of the association between respiratory symptoms in subjects with asthma and daily air pollution in Spokane, Washington. Inhal Toxicol. 2004;16:809–15.
22 Sint T, Donohue JF, Ghio AJ. Ambient air pollution particles and the acute exacerbation of chronic obstructive pulmonary disease. Inhal Toxicol. 2008;20:25–9.
23 Halbert RJ, Isonaka S, George D, Iqbal A. Interpreting COPD prevalence estimates: what is the true burden of disease? Chest. 2003;123:1684–92.
24 Probst-Hensch N, Tanner M, Kessler C, Burri C, Künzli N. Prevention – a cost-effective way to fight the non-communicable disease epidemic: an academic perspective of the United Nations High-level NCD Meeting. Swiss Med Wkly. 2011;141:w13266.
25 Probst-Hensch NM. Chronic age-related diseases share risk factors: do they share pathophysiological mechanisms and why does that matter? Swiss Med Wkly. 2010;140:w13072.
26 Carracedo-Martínez E, Taracido M, Tobias A, Saez M, Figueiras A. Case-crossover analysis of air pollution health effects: a systematic review of methodology and application. Environ Health Perspect. 2010;118:1173–82.
27 Janes H, Sheppard L, Lumley T. Case-crossover analyses of air pollution exposure data: referent selection strategies and their implications for bias. Epidemiol. 2005;16:717–26.
28 Fischer A, Gehrig R, Hofer P. Russmessungen in der Aussenluft-Methodik und Reultate. Bern: Bundesamt fur Umwelt, Wld und Landschaft (BUWAL); 1997.
Authors’ contribution: C. Schindler, O. Brändli, W. Karrer, J.-M. Tschopp, and N. Künzli made substantial contributions to conception and design, and acquisition of data. A. J. Mehta, C. Schindler, and N. Künzli contributed to the analysis and interpretation of data. A. J. Mehta, C. Schindler, L. Perez, N. Probst-Hensch, J. Schwartz, O. Brändli, W. Karrer, J.-M. Tschopp, T. Rochat and N. Künzli were involved in drafting the manuscript or revising it critically for important intellectual content. A. J. Mehta, C. Schindler, L. Perez, N. Probst-Hensch, J. Schwartz, O. Brändli, W. Karrer, J.-M. Tschopp, T. Rochat and N. Künzli gave final approval of the version of this manuscript to be published.
Funding / potential competing interests: Research support: the Swiss National Science Foundation (grants no 33CSCO-108796, 3247BO-104283, 3247BO-104288, 3247BO-104284, 3247-065896, 3100-059302, 3200-052720, 3200-042532, 4026-028099), the Federal Office for Forest, Environment and Landscape, the Federal Office of Public Health, the Federal Office of Roads and Transport, the canton's government of Aargau, Basel-Stadt, Basel-Land, Geneva, Luzern, Ticino, Valais, and Zurich, the Swiss Lung League, the canton's Lung League of Basel Stadt/ Basel Landschaft, Geneva, Ticino, Valais, and Zurich.