Short-term increases in daily mortality related to airborne concentrations of particulate matter (diameter ≤10 μm [PM10]; or ≤2.5 μm [PM2.5]) have been well established based on large multicity studies conducted in several parts of the world.1 However, some aspects of the PM-mortality association are still research priorities2 and identifying specific subpopulations most sensitive to air pollutants is one of the key questions.3–7 Several individual characteristics have been indicated as potential effect modifiers of the air pollution-mortality association, including age, sex, socioeconomic status,8–9 preexisting cardiovascular and respiratory diseases,4,10,11 diabetes,12–13 obesity,14 and use of statins.15 Recently, the modification of the PM10-mortality association was evaluated by several individual factors in a large case-crossover study of 20 US cities, including socio-demographic variables, location of death, and secondary diagnoses on the death certificate.16 Although the PM10-mortality association in Europe has been well-described, and some effect modifiers at the population level have been identified,17–18 no individual-level effect modifiers have been evaluated. Even so, clarification of individual susceptibility to air pollution might be useful to understand the physio-pathologic mechanisms leading to the clinical health effects and to reduce exposure of specific subgroups during periods of increased air pollution.
We aimed to identify various conditions that may render persons more likely to die because of particulate air pollution. In particular, we evaluated whether socio-demographic characteristics, place of death and individual clinical attributes are effect modifiers of the PM10-mortality association in 9 Italian cities. As in a previous study on vulnerability to high temperature,19–20 extensive record linkage procedures were used to characterize subjects for previous and concurrent morbidity.
The study population consisted of 321,024 subjects aged 35 years or older from 9 Italian cities (in north-south order: Turin, Milan, Mestre, Bologna, Florence, Pisa, Rome, Taranto, and Palermo) who died of natural causes (International Classification of Diseases, 9th revision—ICD IX: 1–799) in one of the above-listed cities between 1997 and 2004. Individual records were retrieved from Regional Registers of Causes of Death, which also record age and sex. Census blocks of residence were collected for 4 centers (Turin, Milan, Bologna, Rome) by record linkage with city-specific population registers. This was used to attribute the median income of the census block of residence to each subject (family income in 1998, provided by the Ministry of Finance), as an area-based indicator of socio-economic status. Income was divided into 4 categories from the city-specific distribution: up to the 20th percentile, 20th–50th percentile, 51st–80th percentile, above 80th percentile.
Individual fiscal codes were used to link records from the city-specific mortality datasets with the Regional Hospital Discharge files (which include hospitalizations in public and private hospitals nationwide of all resident citizens) to retrieve individual information about history of hospitalizations over the past 2 years.
Considering location of death, date of hospitalization and date of death, we could distinguish the following: out-of-hospital deaths (neither admission nor discharge within the last 4-week period), subjects discharged between 2 and 28 days before death, in-hospital deaths, and deaths in a nursing home (information on nursing homes available only for Florence, Milan and Turin). For those who died in-hospital, the number of days since they were hospitalized was considered (0–1, 2–28, more than 28 days). In addition, information about the specific hospital ward where the patient was 2 days before death (general medicine, medical or surgical wards with low-moderate care load, medical or surgical wards with high care load, intensive care units) was available for 4 (Bologna, Milan, Rome, and Turin) of the 9 cities.
Chronic and acute conditions for all subjects were defined based on the time elapsed between hospitalization and death. For each individual, all hospital admissions in the 2 years preceding death (excluding the last 28 days) were selected, considering both primary causes of admission and secondary contributing diagnoses. Each subject was classified according to a list of 17 groups of diagnoses chosen by adapting the Elixhauser list of comorbidities.21 The 28-day window was applied to distinguish chronic conditions from suddenly deteriorating health in the few weeks or days before death. Finally, for subjects who died in-hospital, we evaluated the acute clinical factors that contributed to death considering a group of conditions noted as primary diagnoses on the hospital records during the last 28 days.
Data on PM10 from fixed monitoring sites representing concentrations in residential areas were collected for each city according to standard procedures already employed in several European studies of air pollution.22 Three different methods have been used to measure PM10 in Italy: beta automatic, tapered element oscillating microbalance (TEOM) and gravimetric. The beta automatic method is well in agreement with the gravimetric method, whereas TEOM underestimates PM10, especially at high levels of concentrations. Therefore, we applied conversion factors when data from TEOM were available. We obtained hourly or daily values of PM10, according to the different methods of measure, through city-specific air monitoring networks managed by Regional Environmental Protection Agency or local authority. For valid 24-hour mean values, at least 75% of the hourly data for each monitor needed to be present. For each city we estimated a daily average by averaging monitor-specific daily means.23 Missing data on the aggregate level were replaced using a formula adapted from the APHEA (Air Pollution and Health—A European Approach) method.22 In particular, when the daily information for a station was missing, it was imputed as a weighted mean of values from other stations. When the daily data of a pollutant were missing for all available stations, the value was left missing. The average of the current and the preceding day means PM10 (lag 0–1) was used as the exposure variable on the basis of previous investigations.17,23
Daily information on temperature, humidity and barometric pressure were provided by the Italian Air Force Meteorological Service. Apparent temperature24 (a composite index that takes into account air temperature and humidity at the same time) was calculated. Information on city-specific influenza epidemics were collected from the national surveillance system. Influenza epidemics were defined on the basis of the weekly estimates of influenza incidence, as reported by the Italian National Institute of Health.
The case-crossover approach25 was applied to study the effect of PM10 (lag 0–1) on all noninjury mortality. This is a special case-control design in which every case serves as its own control, where control records are defined as days on which the subject did not experience the study outcome. Since each case-set is made up of a single subject, all individual time-invariant characteristics are perfectly matched, and this removes the need for adjustment for each potential known or unknown individual confounder. Control-days were selected according to the time-stratified approach, by dividing the study period into monthly strata and matching controls with cases on the same days of the week within each stratum.26
A conditional logistic regression was thus applied to each city, while controlling for influenza epidemics, population decreases during vacation periods, holidays, and barometric pressure, while time-trends and day of the week were controlled for by design. Controlling for the confounding effect of apparent temperature was more complex. Many studies have documented the J-shaped relationship between high summer temperatures and total mortality,19,27–29 while a few have depicted the relationship between cold temperatures and mortality as linear with negative slope.30 In the present study we combined the 2 components. The left side of the mortality-temperature curve was approximated by a linear term of the apparent temperature averaged from the previous 6 days (lag 1–6) (up to 9°C) to catch the latent effect of cold temperatures. The central and right side of the curve instead were fitted by superimposing a piecewise linear spline of apparent temperature (lag 0–1) with 2 city-specific inner knots chosen based on visual inspection, AIC, and segmented regression.31
For each city, the effect of 10 μg/m3 change in PM10 (lag 0–1) on daily mortality was studied. Since we were interested in effect modification by several individual conditions, no attempt was made to explore different lags. Because age is the strongest determinant of mortality, and is typically correlated with many of the modifiers examined, we standardized the stratum-specific estimates by age group in the effect-modification analyses. The weights were chosen as the relative frequencies of the age groups 65–74, 75–84 and 85+ years in the total population of people 65+ year old in the 9 cities together.
Finally, a random-effects meta-analysis of the city-specific estimates was performed, using the maximum likelihood method.32 All results are expressed as percent increase in risk of death and 95% confidence intervals (95% CI) associated with a 10 μg/m3 variation in PM10 (lag 0–1). Effect modification was tested, and P values of the comparisons between the putative effect modifiers and their reference categories are reported. The statistical analyses were performed using R software version 2.3.0,33 and Stata 8.0.34
Table 1 shows the study period in each city, the number of deaths included in the analysis, and the environmental characteristics, including the number of fixed monitoring sites. The period of study varied from a minimum of 3 years (Mestre, Palermo, and Taranto) to a maximum of 7 years (Rome and Turin). Annual PM10 levels ranged from an average of 35.1 μg/m3 in Pisa to 71.5 μg/m3 in Turin. Mean apparent temperature was lower in the northern cities (Milan and Turin), intermediate in central cities (Rome, Florence), and higher in southern cities (Palermo and Taranto).
The effect of PM10 increased with age (Table 2), with a null effect on 35–64 year olds and a large effect on 85+ year olds. Therefore, the overall effect on the 65+ year old group was slightly higher (0.75% increase [95%CI = 0.42% to 1.09%]) than the corresponding effect in the 35+ years age group (0.60% increase [0.31% to 0.89%]). Because age was found to be such a significant effect modifier, we restricted all subsequent analyses among the elderly (65+ years) and performed the age adjustment in all effect modification analyses as indicated in the methods. Figure 1 shows the city-specific and combined results of the association between PM10 and overall mortality in the 65+ years group. All city-specific point estimates were positive, although they were not homogenous, with the lowest value found in Taranto and the highest in Florence.
Table 2 shows the age-adjusted pooled results of the association between PM10 and mortality by sex, area-based-income and location of death; P values for the effect modification are also reported. The PM10-mortality association showed very similar effect estimates for men and women. There was a weaker PM10 effect (0.30% vs. 0.80%) among the most affluent people (the upper quintile of median income of the census block of residence) compared with the least affluent subjects (lowest quintile), although the interaction term was not statistically significant. There was no clear effect modification of the location of death. The effect of particulate matter on mortality remained strong for those who died in a hospital (0.65% [0.33% to 0.97%]) and was comparable to the effect detected in those who died out-of-the hospital (0.71% [0.32% to 1.11]). The highest risk was found for those who died at home who had been recently discharged from a hospital (1.34% [0.49% to 2.20%]); the lowest risk was for those who died in nursing homes (although the information was limited to 3 cities).
Table 3 displays the combined results of primary and secondary contributing causes of admissions between 29 days and 2 years before death for the 17 groups of diagnoses. The effect estimates were similar for those who were hospitalized in the 2-year period (175,625 subjects) and those who were not (99,458 subjects) (0.63% vs. 0.80% increase). When effect modification was evaluated with a statistical test, no chronic condition seemed to confer a specific susceptibility. However, increased risks were found for patients with diabetes (1.03% [0.28% to 1.79%]), patients with hypertension (0.54% [0.00% to 1.10%]), and with chronic pulmonary diseases (0.84% [0.17% to 1.52%). There were also a number of diseases that showed below average effect estimates (for malignant neoplasms, 0.17% [0.49% to 0.84%]; for coagulation defects: −3.46% [−6.42% to 0.40%]; for other ischemic heart diseases: 0.02% [−0.65% to 0.69%]).
Table 4 reports the results of the analyses restricted to in-hospital deaths (140,400 subjects). The age effect was confirmed, with the highest risk among those 85+years (0.82% [0.27% to 1.37%]). There were no sex differences. There were slight differences in the effect of PM10 according to the time elapsed since the hospitalization (0–1 day, 2–28 days, more than 28 days) with a tendency of a stronger effect among those admitted more recently. The effect among people last admitted 2–28 days before death, who were thus inside the hospital when exposed to the lagged PM10, was 0.63% (0.22% to 1.05%). The in-hospital PM10 effect was different by ward, ie, it was higher for patients staying in general medicine (0.48%) or “low-load” care units (0.75%), while it was null for those staying in high-load or intensive care units (−1.11% and 0.07%, respectively) (P values for effect modification = 0.14 and 0.74).
Table 5 shows the results of the analysis conducted among 111,511 subjects who died in-hospital according to the primary diagnosis that led to their last hospitalization (≤28 days before death). The conditions with the largest effect estimates were diseases of pulmonary circulation (acute and chronic pulmonary heart diseases) (4.56% [0.75% to 8.51%]) and heart failure (1.67% [0.30% to 3.04%]). When the analyses of acute conditions were restricted to 86,717 patients admitted to the hospital 2–28 days before death, and who were therefore indoors during the PM10 0–1 lag effect (Table 6), effect modification for pulmonary circulation disorders (4.83% [0.08 to 9.81], P value for effect modification = 0.081) and heart failure (2.44% [0.86 to 4.05], P value for effect modification = 0.017) was suggested.
The study indicates an overall PM10 effect on total mortality similar to what has been reported in the European APHEA study22 and in a recent meta-analysis.35 However, there are individual effect modifiers of the PM10 effect. Higher estimates were found for very elderly people, and social class differences were suggested, particularly lower effects among the most affluent subjects. The risk was higher in noninstitutionalized people, especially among those recently discharged from the hospital. However, a PM10 effect was detected also among hospitalized patients, with the exception of those in more intensive care. Admission to the hospital in the past 2 years was not an important risk modifier. Finally, among recently hospitalized patients, pulmonary heart diseases and acute heart failure clearly increased the PM10-related risk of death.
Several studies have found a stronger effect of PM10 on mortality in the elderly,16,17 and our results confirm that risk increases with age; subjects 85+ years of age had almost double the risk of 65–74 years olds. The cause of this greater susceptibility may be chronic medical conditions or the loss of compensation mechanisms in the cardio-respiratory system.11 In contrast to previous findings,36 but in agreement with a large American study,16 we did not find a sex effect. Previous research has found effect modification by social class when comparing area-specific analyses of low-SES and high-SES communities,9,37 but there was little evidence of effect modification in a large study in the US38 that considered national data. Another American study16 that used individual education level found that risk in low educational groups was almost twice as high as in highly educated people. Similarly, an analysis that classified residents of Rome based on an area-based index of socioeconomic status (SES) did not find a PM10 effect on daily mortality in high SES groups but found a strong effect in the other SES groups.39 Our study indicates that social class modifies the effect of air pollution 2-fold between those in the low (0.80%) and in the high (0.30%) income categories. Differential susceptibility may be the result of disparities in lifestyle factors, such as tobacco and alcohol use, physical activity, and occupational exposure to toxins, as wells as to reduced access to community and hospital care. Differential exposure to air pollution is an additional possibility, although this was not seen in 1 study.30
The present study provides new insights regarding place of death as an effect modifier. In North American studies, most of the effect of environmental factors such as air pollution and high temperatures has been observed among people who died outside the hospital16,40; a lower effect has been found among in-patients. The use of air filtration and air conditioning in US hospitals has been indicated as the primary explanation. In Italy, air filtration is used only in modern hospitals and the majority did not have air conditioning (until the 2003 heat wave), the exception being intensive care units. In-hospital mortality has been found associated with high temperatures in a study conducted in 4 Italian cities,19–20 and lack of air conditioning was the likely explanation. Similarly, hospital environments in Italy do not offer protection against air pollution because of the lack of air conditioning, with the exception of specialized and intensive care units where air filtration is common. Usually, windows tend to be open in summer and in the warmers parts of spring and fall. In our study people recently discharged from a hospital had a higher PM10 effect (1.34% increase). The reasons for such high vulnerability should be studied.
We were expecting our detailed search of chronic conditions using previous hospital admissions to allow us to single out the specific clinical factors associated with vulnerability to PM10. Some conditions were strong candidates. For instance, studies41–43 have reported hypertension as a modifier of the effects of air pollution. A similar phenomenon may occur in subjects with chronic obstructive pulmonary disease (COPD); Sunyer et al,44 De Leon et al,45 and Forastiere et al46 have reported higher risks of circulatory deaths in subjects with COPD. Finally, there are suggestions that the association between particulate matter and health effects is stronger in persons with increased baseline systemic inflammation and oxidative stress, such as patients with diabetes, obese individuals, and people not using statins.15,47 Contrary to those findings, we were not able to detect specific conditions with particularly high effect estimates, based on previous hospital admissions. In general, the strength of the PM10 effect was the same for those who had been hospitalized in the 2 preceding years and those who had not. It should be noted, however, that the PM10 effects were robust among people previously hospitalized for diabetes, hypertension and COPD. Overall, the results suggest that the associated effects of PM10 are widespread and not necessarily limited to persons with severe conditions requiring hospital care. On the other hand, cardiovascular, respiratory and metabolic diseases are highly prevalent in the general population and do not always require hospital admission. Although our work did not specifically address the issue of mortality displacement, our results may suggest that the association between particulate matter and mortality does not involve only an advance in the timing of death by few days in a selected group of frail individuals.48
The observation that people with a previous diagnosis of cancer, coagulation defects or other ischemic heart diseases were not affected by air pollution should be noted; there was statistical evidence of effect modification. We have no ready explanation for these findings, although one may speculate that the diseases themselves or their treatments may interfere with the effects of airborne particles. Patients with malignant cancer are at high risk of death and they usually receive chemotherapy and radio-therapy that produce very high oxidative stress, so that the additional strain produced by air pollution may be not important. Recent literature49 suggests that air pollution is associated with changes in global coagulation function with a tendency towards hyper-coagulability after short-term exposure to air pollution. People who already have coagulation problems may be less susceptible to this mechanism especially if chronic coagulation interfering drugs are used. Finally, people with chronic ischemic heart diseases may be taking antithrombotic drugs (eg, aspirin) and agents that lower cholesterol (eg, statins) that would modify their risk profile. These findings need confirmation.
We found evidence that the PM10 risk is doubled among patients with heart failure and quadrupled among those with disorders of pulmonary circulation. Heart failure due to left-ventricular systolic dysfunction is usually associated with pulmonary hypertension. Increased blood pressure in lung vessels is also secondary to pulmonary embolism and COPD, leading to right heart failure and cor pulmonale. The prevalence of these conditions in the general population is relatively high, especially among the elderly (5%–10% above 65 years).50 Very few studies have been able to consider clinical events in the last day or days before death as possible effect modifiers of the air pollution-mortality association, and suggestive indications are available for heart failure. Goldberg et al4,11 estimated the association between daily mortality and daily concentrations of particles in the ambient air of Montreal, Quebec among individuals with cardiovascular and respiratory conditions diagnosed before death. People who had acute lower respiratory diseases and congestive heart failure shortly before death had a particulate matter-associated mortality about twice that found over the entire study population. Kwon et al51 selected a cohort of congestive heart failure patients from the medical insurance files of a hospital among residents of Seoul, South Korea. Effects attributable to PM10 among the ill cohort appeared larger than among the general population and a stronger effect of urban particles was found in patients for whom a shorter period (<1 month) had elapsed since the index hospitalization than among patients who had survived the first month. Finally, in a large investigation of daily hospital admission rates from 1999 to 2002 in 204 US counties,52 the largest PM2.5 effect was found on heart failure.
The findings regarding heart failure and pulmonary heart conditions may help in understanding the mechanisms of the PM10 effect. Goldberg et al11 have already suggested that ambient particles may affect the heart indirectly by modifying endothelin homeostasis in the lungs. He speculated that elevated circulating endothelins, associated with oxidative stress linked to particulate matter exposure, exacerbate congestive heart failure, precipitate ischemic heart disease, promote infarct extension, and also possibly affect individuals with conditions involving endothelial dysfunctions such as diabetes and atherosclerosis. Endothelin-1, an endothelium-derived mediator, has fibrogenic and inflammatory properties, is a potent vasoconstrictive agent, has inotropic and mitogenic actions, modulates salt and water homeostasis and plays an important role in maintaining vascular tone and blood pressure. Consequently, endogenous endothelin-1 contributes to the pathophysiology of conditions associated with sustained vasoconstriction, such as heart failure53 and unregulated pulmonary hypertension.54
Several strengths of the present study deserve consideration. This study involved 9 cities and more than 300,000 deaths. The overall results are robust and the case-crossover analysis gives estimates consistent with those obtained in the APHEA study17 and in a previous Italian time-series analysis.23 Linking individual data records offered the opportunity to exploit individual information that is rarely available in other European countries. Adjusting for age in evaluating effect modification allowed us to assess the independent role of several factors that are linked to age.
Some limitations must also be taken into account. Chronic clinical conditions are based on hospital admissions and suffer from the limits of accuracy of the source used.55 In addition, simple hospital discharge summaries provide only an approximation of acute illnesses or exacerbations of chronic diseases. We tried to increase the sensitivity of the definition of chronic conditions by using all hospital admissions in the period of 29 days to 2 years before death and by considering both primary and contributory causes. Studying the clinical conditions that characterized the subjects in the few weeks before death is more complex. The information from the last hospitalization was used, which in most cases ended in death, and, to increase the specificity of the definition, the analysis was restricted to primary causes only. Additional data could have been useful to better define chronic susceptibility (individual habits, smoking status, obesity, etc.) but, unlike in the study conducted by Goldberg et al,11 we did not have access to other databases from healthcare services outside the hospital, or to filled drug prescriptions. Finally, despite the large size of the study, we had limited statistical power to detect effect modification.
In conclusion, we found that the greatest risk of death associated with high concentrations of PM10 was among the elderly, among people of lower socio-economic status, and among in-patients not in intensive care and with altered pulmonary circulation and heart failure.
SISTI Group: Italian Study on Susceptibility to Temperature and Air Pollution. Members: Bologna—Agostini D, De Lisio S, Miglio R, Pandolfi P, Scarnato C; Firenze—Biggeri A, Chellini E, Mallone S; Mestre—Simonato L, Tessari R; Milano—Bisanti L, Rognoni M, Russo A; Palermo—Cernigliaro A, Scondotto S; Pisa—Vigotti MA; Roma—Belleudi V, de'Donato F, Forastiere F, Michelozzi P, Perucci CA, Picciotto S, Stafoggia M; Taranto—Primerano R, Serinelli M; Torino—Berti G, Cadum E, Caranci N, Chiusolo M, Demaria M.
1. Pope CA III, Dockery DW. Health effects of fine particulate air pollution: lines that connect [Review]. J Air Waste Manage Assoc
2. National Research Council. Research Priorities for Airborne Particulate Matter. IV. Continuing Research Progress.
Washington, DC: The National Academies Press; 2004.
3. Schwartz J. What are people dying of on high air pollution days? Environ Res
4. Goldberg MS, Bailar JC III, Burnett RT, et al. Identifying subgroups of the general population that may be susceptible to short-term increases in particulate air pollution: a time-series study in Montreal, Quebec. Res Rep Health Eff Inst
5. Utell MJ, Frampton MW, Zareba W, et al. Cardiovascular effects associated with air pollution: potential mechanisms and methods of testing. Inhal Toxicol
6. Pope CA III. Epidemiology of fine particulate air pollution and human health: biologic mechanisms and who's at risk? Environ Health Perspect
. 2000;108(suppl 4):713–723.
7. Peters A. Susceptible subgroups: the challenge of studying interactions. Epidemiology
8. Zanobetti A, Schwartz J, Gold D. Are there sensitive subgroups for the effects of airborne particles? Environ Health Perspect
9. O'Neill MS, Jerrett M, Kawachi I, et al. Health, wealth, and air pollution: advancing theory and methods [Review]. Environ Health Perspect
10. Zanobetti A, Schwartz J. Particulate air pollution, progression, and survival after myocardial infarction. Environ Health Perspect
11. Goldberg MS, Burnett RT, Bailar JC III, et al. Identification of persons with cardiorespiratory conditions who are at risk of dying from the acute effects of ambient air particles. Environ Health Perspect
. 2001;109(suppl 4):487–494.
12. Zanobetti A, Schwartz J. Are diabetics more susceptible to the health effects of airborne particles? Am J Respir Crit Care Med
13. Bateson TF, Schwartz J. Who is sensitive to the effects of particulate air pollution on mortality? A case-crossover analysis of effect modifiers. Epidemiology
14. Dubowsky SD, Suh H, Schwartz J, et al. Diabetes, obesity, and hypertension may enhance associations between air pollution and markers of systemic inflammation. Environ Health Perspect
15. Schwartz J, Park SK, O'Neill MS, et al. Glutathione-S-transferase M1, obesity, statins, and autonomic effects of particles: gene-by-drug-by-environment interaction. Am J Respir Crit Care Med
16. Zeka A, Zanobetti A, Schwartz J. Individual-level modifiers of the effects of particulate matter on daily mortality. Am J Epidemiol
17. Katsouyanni K, Touloumi G, Samoli E, et al. Confounding and effect modification in the short-term effects of ambient particles on total mortality: results from 29 European cities within the APHEA2 project. Epidemiology
18. Analitis A, Katsouyanni K, Dimakopoulou K, et al. Short-term effects of ambient particles on cardiovascular and respiratory mortality. Epidemiology
19. Stafoggia M, Forastiere F, Agostini D, et al. Vulnerability to heat-related mortality: a multicity, population-based, case-crossover analysis. Epidemiology
20. Stafoggia M, Forastiere F, Agostini D, et al. Factors affecting in-hospital heat-related mortality: a multi-city case-crossover analysis. J Epidemiol Community Health
21. Elixhauser A, Steiner C, Harris DR, et al. Comorbidity measures for use with administrative data. Med Care
22. Katsouyanni K, Schwartz J, Spix C, et al. Short term effects of air pollution on health: a European approach using epidemiologic time series data: the APHEA protocol. J Epidemiol Community Health
. 1996;50(suppl 1):S12–S18.
23. Biggeri A, Bellini P, Terracini B. [Meta-analysis of the Italian studies on short-term effects of air pollution—MISA 1996–2002.] Epidemiol Prev
. 2004;28(4–5 suppl):4–100. Italian.
24. Kalkstein LS, Valimont KM. An evaluation of summer discomfort in the United States using a relative climatological index. Bull Am Meteorol Soc
25. Maclure M. The case-crossover design: a method for studying transient effects on the risk of acute events. Am J Epidemiology
26. Levy D, Lumley T, Sheppard L, et al. Referent selection in case-crossover analyses of acute health effects of air pollution. Epidemiology
27. Basu R, Samet JM. Relation between elevated ambient temperature and mortality: a review of the epidemiologic evidence. Epidemiol Rev
28. Braga ALF, Zanobetti A, Schwartz J. The time course of weather-related deaths. Epidemiology
29. Curriero FC, Heiner KS, Samet JM, et al. Temperature and mortality in eleven cities of the Eastern United States. Am J Epidemiol
30. Cold exposure and winter mortality from ischaemic heart disease, cerebrovascular disease, respiratory disease, and all causes in warm and cold regions of Europe. The Eurowinter Group. Lancet
31. Muggeo VM. Estimating regression models with unknown break-points. Stat Med
32. van Houwelingen HC, Arends LR, Stijnen T. Advanced methods in meta-analysis: multivariate approach and meta-regression. Stat Med
33. R Development Core Team. R: a language and environment for statistical computing. (ISBN 3-900051-00-3). Vienna, Austria: R Foundation for Statistical Computing, 2004. Available at: http://www.R-project.org
34. StataCorp. 2003. Stata Statistical Software: Release 8. College Station, TX: StataCorp LP.
35. Anderson HR, Atkinson RW, Peacock JL, et al. Meta-analysis of time-series studies and panel studies of particulate matter (PM) and ozone (O3). Report of a WHO task group. Copenhagen: WHO Regional Office for Europe; 2004. Available at: http://www.euro.who.int/document/e82792.pdf
36. Zanobetti A, Schwartz J. Race, gender, and social status as modifiers of the effects of PM10 on mortality. J Occup Environ Med
37. Jerrett M, Burnett RT, Brook J, et al. Do socioeconomic characteristics modify the short term association between air pollution and mortality? Evidence from a zonal time series in Hamilton, Canada. J Epidemiol Community Health
38. Health Effects Institute. Special report: revised analyses of time-series studies of air pollution and health. Part II. Revised analyses of the National Morbidity, Mortality, and Air Pollution Study. In: Revised Analyses of Selected Time-Series Studies
. Boston, MA: Health Effects Institute; 2003.
39. Forastiere F, Stafoggia M, Tasco C, et al. Socioeconomic status, particulate air pollution, and daily mortality: differential exposure or differential susceptibility. Am J Ind Med
40. Medina-Ramón M, Zanobetti A, Cavanagh DP, et al. Extreme temperatures and mortality: assessing effect modification by personal characteristics and specific cause of death in a multi-city case-only analysis. Environ Health Perspect
41. Holguin F, Tellez-Rojo MM, Hernandez M, et al. Air pollution and heart rate variability among the elderly in Mexico City. Epidemiology
42. Liao D, Duan Y, Whitsel EA, et al. Association of higher levels of ambient criteria pollutants with impaired cardiac autonomic control: a population-based study. Am J Epidemiol
43. Park SK, O'Neill MS, Vokonas PS, et al. Effects of air pollution on heart rate variability: the VA normative aging study. Environ Health Perspect
44. Sunyer J, Schwartz J, Tobias A, et al. Patients with chronic obstructive pulmonary disease are a susceptible population of dying due to urban particles. Am J Epidemiol
45. De Leon SF, Thurston GD, Ito K. Contribution of respiratory disease to nonrespiratory mortality associations with air pollution. Am J Respir Crit Care Med
46. Forastiere F, Stafoggia M, Picciotto S, et al. A case-crossover analysis of out-of-hospital coronary deaths and air pollution in Rome, Italy. Am J Respir Crit Care Med
47. O'Neill MS, Veves A, Zanobetti A, et al. Diabetes enhances vulnerability to particulate air pollution-associated impairment in vascular reactivity and endothelial function. Circulation
48. Dominici F, McDermott A, Zeger SL, et al. Airborne particulate matter and mortality: timescale effects in four US cities. Am J Epidemiol
49. Baccarelli A, Zanobetti A, Martinelli I, et al. Effects of exposure to air pollution on blood coagulation. J Thromb Haemost
50. Davies M, Hobbs F, Davis R, et al. Prevalence of left-ventricular systolic dysfunction and heart failure in the Echocardiographic Heart of England Screening study: a population based study. Lancet
51. Kwon HJ, Cho SH, Nyberg F, et al. Effects of ambient air pollution on daily mortality in a cohort of patients with congestive heart failure. Epidemiology
52. Dominici F, Peng RD, Bell ML, et al. Fine particulate air pollution and hospital admission for cardiovascular and respiratory diseases. JAMA
53. Giannessi D, Del Ry S, Vitale RL. The role of endothelins and their receptors in heart failure [Review]. Pharmacol Res
54. Keith IM. The role of endogenous lung neuropeptides in regulation of the pulmonary circulation. Physiol Res
55. Jollis JG, Ancukiewicz M, DeLong ER, et al. Discordance of databases designed for claims payment versus clinical information systems. Implications for outcomes research. Ann Intern Med