Research | Open | Open Peer Review | Published:
Traffic-related air pollution exposures and changes in heart rate variability in Mexico City: A panel study
Environmental Healthvolume 12, Article number: 7 (2013)
While air pollution exposures have been linked to cardiovascular outcomes, the contribution from acute gas and particle traffic-related pollutants remains unclear. Using a panel study design with repeated measures, we examined associations between personal exposures to traffic-related air pollutants in Mexico City and changes in heart rate variability (HRV) in a population of researchers aged 22 to 56 years.
Participants were monitored for approximately 9.5 hours for eight days while operating a mobile laboratory van designed to characterize traffic pollutants while driving in traffic and “chasing” diesel buses. We examined the association between HRV parameters (standard deviation of normal-to-normal intervals (SDNN), power in high frequency (HF) and low frequency (LF), and the LF/HF ratio) and the 5-minute maximum (or average in the case of PM2.5) and 30-, 60-, and 90-minute moving averages of air pollutants (PM2.5, O3, CO, CO2, NO2, NOx, and formaldehyde) using single- and two-pollutant linear mixed-effects models.
Short-term exposure to traffic-related emissions was associated with statistically significant acute changes in HRV. Gaseous pollutants – particularly ozone – were associated with reductions in time and frequency domain components (α = 0.05), while significant positive associations were observed between PM2.5 and SDNN, HF, and LF. For ozone and formaldehyde, negative associations typically increased in magnitude and significance with increasing averaging periods. The associations for CO, CO2, NO2, and NOx were similar with statistically significant associations observed for SDNN, but not HF or LF. In contrast, PM2.5 increased these HRV parameters.
Results revealed an association between traffic-related PM exposures and acute changes in HRV in a middle-aged population when PM exposures were relatively low (14 μ g/m3) and demonstrate heterogeneity in the effects of different pollutants, with declines in HRV – especially HF – with ozone and formaldehyde exposures, and increases in HRV with PM2.5 exposure. Given that exposure to traffic-related emissions is associated with increased risk of cardiovascular morbidity and mortality, understanding the mechanisms by which traffic-related emissions can cause cardiovascular disease has significant public health relevance.
Many studies have demonstrated the association between air pollution exposure, specifically fine particulate matter (PM2.5), and increased cardiovascular morbidity and mortality [1–9]. Short-term PM exposures have been linked to acute cardiovascular events including increased odds of having a myocardial infarction, cardiac arrhythmia, and venous thrombosis , while long-term exposure to PM has been associated with increased risk and progression of atherosclerosis .
While the exact biological mechanism linking exposure to PM and cardiovascular outcomes remains unknown [3, 12], alterations in heart rate variability (HRV) are thought to be one of the pathophysiologic pathways whereby PM affects the cardiovascular system . HRV is an indicator of the relative balance of parasympathetic and sympathetic autonomic control of the heart rate, and changes (both increases and decreases) in this metric have been associated with both ambient and traffic-related PM air pollution [14–17].
Typically, higher PM concentrations have been associated with decreased HRV in elderly populations and in patients with current or underlying cardiovascular disease [12, 18, 19]. Findings from studies on the association between HRV and PM in younger populations, however, have been inconsistent. In a panel study of 76 young college students, HRV indices declined in single-pollutant models with PM10, PM2.5, sulfate, nitrate, and ozone (O3) . In a controlled exposure study (mean age = 27 years), exposure to concentrated ambient particles had no consistent effect on HRV indices . Similarly, no consistent effect of diesel exhaust on HRV was observed in a separate double-blind, crossover, controlled-exposure study (mean age = 32 years for healthy subjects, mean age = 41 years for those with metabolic syndrome) . In this case, the controlled-exposure studies may not have accurately simulated actual environmental conditions. In an occupational panel of young boilermakers (mean age = 38 years), significant increases in an HRV index (standard deviation of normal-to-normal intervals or SDNN) were observed for every 1 μ g/m3 increase in lead and vanadium concentrations . Cardiovascular comorbidity, such as hypertension, which are more prevalent in older populations, has been shown to increase susceptibility to fine particulate matter-mediated reductions in HRV . This phenomenon may determine why older individuals have a different response to air pollution-mediated changes in cardiac autonomic regulation, when compared to younger persons.
Several studies have demonstrated a stronger association between cardiovascular endpoints and particles originating from traffic as compared to other sources [16, 18, 25–27]. The study of traffic-related air pollutants, however, is complicated due to the nature of traffic exposures, which may vary over short distances and thus limit the use of centralized exposure monitoring data in epidemiological studies.
To reduce the potential for exposure misclassification of traffic-related emissions, several studies have monitored in-vehicle pollutant exposure in young populations. Increases in HRV were observed in association with in-vehicle PM2.5 exposure in a group of young North Carolina State Highway Patrol troopers (mean age = 27 years) . In a population of young, highly-exposed taxi drivers in Beijing (mean age = 36 years), low PM2.5 exposures were associated with relatively high HRV, whereas higher PM2.5 exposures were associated with relatively low HRV .
In addition to the discrepancy in the relationship between PM exposure and HRV response, there is a lack of knowledge about the health effects due to the potential synergy between PM2.5 and ambient gaseous co-pollutants . Several studies have simultaneously evaluated the association of traffic-related PM and gaseous pollutants with HRV [12, 21, 22, 28, 30, 31]. This is particularly important as exposure to traffic-related emissions in close proximity is characterized by a rich mixture of fine PM and gaseous pollutants that are different from the mixture of the background air pollution exposure .
To address the aforementioned gaps in the literature, we examined the association between real-time, traffic-related air pollution exposures and acute sub-clinical cardiovascular outcomes in a middle-aged population. We used a panel study design with repeated measures to account for personal factors while enhancing the statistical power to detect associations through high and variable exposure levels. Our hypothesis was that there is an inverse exposure-response relationship between HRV parameters and traffic-related pollutants including PM2.5, carbon dioxide (CO2), carbon monoxide (CO), nitrogen dioxide (NO2), nitrogen oxides (NOx), O3, and formaldehyde.
Participants and study design
This repeated-measures panel study was conducted February 11–23, 2002, in the Mexico City Metropolitan Area as part of an effort is to contribute to the understanding of the air quality problem in megacities by conducting measurements and modeling studies of atmospheric pollutants . A convenience sample of the sixteen researchers, between the ages of 22 and 56 years, associated with the project participated in the study. A self-administered questionnaire was used to collect personal data, including sex, age, smoking status and hypertension history. The study design and methods were reviewed and approved by the human subjects committee at the National Institute of Public Health. All participants signed an informed consent form before participating in the study.
Real-time measurements of PM2.5, CO2, CO, NO2, NOx, O3, and formaldehyde were collected in a van-based mobile laboratory . The van pursued specific vehicles to measure their emissions, drove transects across the city to capture the spatial variation in pollutants, and parked at several locations in the city . The nine individual mobile episodes lasted 1–10 hours. The van drove past a variety of point sources throughout the city, including residential (e.g. biomass burning), industrial (e.g. metal welding, factories, oil burning), livestock, landfill, and sewage treatment sources.
During data collection, approximately six participants at a time wore an Aria Digital Holter Monitor and rode in the van or stood outside the van when it was in a stationary location. In-vehicle and ambient air exchange (and corresponding participant exposures) occurred through frequently open windows in the van’s cab and ambient air penetration through air conditioning vents and through open doors during the van’s frequent short stops during day and longer stops for meals.
Air pollutant measurements
All air pollution measurements were taken in a van-based mobile laboratory developed by Aerodyne Research Inc. [34, 36]. The van was designed to sample and characterize mobile and fixed-site emission plumes, as well as characterize gaseous and particulate emissions from selected classes of vehicles, including heavy-duty diesel trucks, buses, and colectivos (ubiquitous small gasoline or condensed natural gas powered microbuses). It was outfitted with state-of-the-art, fast-response instruments, including a non-dispersive infrared (NDIR) unit (Li-Cor LI 6262) for CO2, an Aerodyne tunable infrared laser differential absorption spectrometer (TILDAS) for NO2 and formaldehyde (HCHO), a NDIR analyzer for CO, a chemiluminscent analyzer (Thermo 42C) for nitrogen oxides (NOx), a UV monitor (Thermo Environmental 49–003) for O3, and an aerosol photometer for PM2.5 (TSI Dustrak 8520) .
Sampler inlets were positioned well above and forward of the vehicle engine and generator exhaust outlets . The Dustrak was mounted on a shelf in the van, and the stainless steel and Tygon sampling lines were designed to minimize particle deposition .
The intrinsic measurement period varied depending on the type of instrument used. The basic measurement interval was ~1 second for CO2, NO2, NOx, formaldehyde and PM2.5; the interval increased to ~20 seconds for O3 and CO. Except for a few power or computer loss periods, the instruments measured and recorded continuously. For each gaseous pollutants, i.e. O3, CO, CO2, NO2, NOx, and formaldehyde, the maximum value in a given 5-minute window was recorded. From these 5-minute maxes, 30-, 60-, and 90-minute average maximums were calculated (by averaging the 5-minute maxes). For PM2.5, 5-minute mean concentrations were recorded and 30-, 60-, and 90-minute means were calculated from these values.
The PM2.5 measurements reflected some uncertainty resulting from the calibration of the aerosol photometer . The Dustrak was calibrated against multiple 24-hour PM2.5 gravimetric samples throughout the field campaign and applied to the factory-calibrated readings (field calibration factor was 0.34±0.02). This method depends on scattering efficiencies, which, in turn, are a function of optical properties and particle size distributions. Though scattering efficiencies of ambient and diesel particles are similar, gasoline particles may differ in this efficiency, which may cause the calibration for individual vehicles to vary by a factor of two or more.
Additionally, distinct “spike events” in the UV O3 monitor on board the van were observed when the van was sampling the ambient diluted exhaust from on-road diesel vehicles . Associated fine particles were assumed to have caused the observed interference. This type of interference could lead to a mean measured O3 concentration that is at most 3% higher than actual concentrations.
Heart rate variability measurements
HRV was obtained from analysis of the ambulatory electrocardiogram recorded using an Aria Digital Holter Monitor (Del Mar Reynolds, US). Participants were allowed to participate on multiple occasions up to eight days for a total of 48 person-days. Electrodes were placed on the right parasternal and precordial areas, and continuous monitoring occurred between 7:05 and 20:05 with the average monitoring period occurring from 9:30 to 19:00. HRV parameters were calculated in 5-minute epochs and abnormal or ectopic beats were manually removed. Digitized Aria Digital Halter recordings were analyzed using a Marquette MARS Workstation that provided an algorithm for HRV analysis and interpolation for removed aberrant QRS complexes (Del Mar Reynolds, US).
HRV parameters included the standard deviation of normal R-R internals (SDNN), which is a time domain measure of overall HRV, and high- (0.18-0.40 Hz) and low-frequency power (0.03-0.15 Hz), which are respectively representative of predominantly parasympathetic and sympathetic autonomic cardiac regulation. (The LF/HF ratio represents the relative balance between sympathetic-vagal nervous activity .
All heart rate variability measures were log-transformed using the natural logarithm to help meet regression assumptions. No predictors were transformed, as model assumptions were reasonably met by transforming the outcome variables only.
The associations between each of the HRV parameters and individual pollutants were examined using linear mixed-effects models (α = 0.05). The gaseous pollutants O3, CO, and NO2 were also individually evaluated in two-pollutant models with PM2.5, and potential multicollinearity between the two pollutants was assessed using the condition number . A condition number greater than 30 suggests moderate multicollinearity could be present, whereas a value over 100 indicates severe multicollinearity is likely occurring. Multicollinarity is a concern because it may affect the stability of point estimates and the accuracy of their inference, leading to incorrect conclusions about associations between a set of predictors and an outcome.
A random subject intercept was incorporated into each model with an exponential covariance structure to account for the unevenly spaced 5-minute measurements over 12 hours, performed over eight unequally spaced days. The exponential covariance structure allows for the correlation within the same subject to decay over time . All analyses were conducted using one or two air pollutant models over periods of 5, 30, 60, and 90 minutes, as described previously.
Models were adjusted for fixed and time-varying factors that were potential confounders. Fixed factors included sex, age (linear), smoking status, and ethnic origin (Mexican or other). Time-varying factors included a categorical variable for time of day (06:00 to 11:59, 12:00 to 15:59, and 16:00 to 20:05), and study day (eight categories). To examine the influence of outlying exposure values, the smallest and largest 5% of pollutant values for a given averaging period (10% total) were excluded in a separate analysis .
Final results (β-values) are presented as the estimated percent change of a given HRV outcome per interquartile range (IQR) increase in the exposure to each air pollutant, controlling for sex, age, smoking status, ethnic origin, time of day, and study day. Estimates were calculated as β = [exp(β′ × IQR) – 1] × 100%, where β′ was the estimated effect of a pollutant from the mixed-model . Similarly, the 95% confidence interval was achieved from the following transformation: [exp(IQR × CI′)– 1] × 100%, where CI′ represents the estimated 95% confidence interval for β′ from the mixed-model. Statistical analyses were performed in R (version 2.13.2) and SAS (version 9.2).
Table 1 includes detailed information on the characteristics of the 16 study participants, including SDNN, high- and low-frequency spectral HRV domain (HF and LF, respectively), and the ratio of LF/HF. The majority (69%) were males. Fifty percent of participants were of Mexican ethnic origin; the remaining were Caucasian American. Participants’ ages ranged from 22 to 56 years, with a mean of 35 years (SD = 11.8) and a median of 31 (IQR = 16.5). The majority were non-smokers (87%), with only one individual reporting hypertension (6%).
Table 2 shows the exposure characteristics over different periods (5-, 30-, 60-, and 90-min intervals) for PM2.5, O3, CO, CO2, formaldehyde, NO2, and NOx. NO2 (mean: 130 ppb) exceeded the one-hour National Ambient Air Quality Standard (NAAQS) concentration (100 ppb) . CO mean one-hour exposure (6 ppm) was below the corresponding NAAQS (35 ppm).
Table 3 gives Spearman’s rank correlation coefficients for the exposure variables from the 5-min time period. Several of the pollutant measurements were strongly correlated (i.e. ρ > 0.70): CO2 and CO, formaldehyde and CO, formaldehyde and CO2, NO2 and CO, NO2 and CO2, NOx and CO, NOx and CO2, and NOx and NO2.
We estimated associations between the HRV measures and exposure to pollutants over different moving averages (from 5 to 90 min) after adjusting for potential confounders (Table 4 and Figure 1). Positive associations were observed between HRV parameters (SDNN, HF, and LF) and PM2.5 exposures, but the LF/HF ratio was negatively associated with PM2.5. The largest percent increases were observed for HF over the 90-min averaging period: a 7.74% (95% CI: 2.3 to 13.3) increase in HF was found per IQR 90-min PM2.5 (8.3 μ g/m3).
No positive associations were found between O3, CO, CO2, NO2, NOx, and formaldehyde exposures and SDNN, HF, and LF. For ozone and formaldehyde, negative associations increased in magnitude and significance with increasing averaging periods for SDNN, HF, and LF (ozone only). The largest effect for each pollutant was observed for HF, with a 16% (95% CI: 9.04 to 23.4) decline per IQR 90-min ozone (65 ppb) and a 12% (95% CI: 3.1 to 20.3) decline per IQR 90-min formaldehyde (34 ppb). The associations for CO, CO2, NO2, and NOx (which were all correlated) were similar with statistically significant associations observed for SDNN. The influence of averaging period differed for each pollutant, and we observed the largest declines in SDNN for each pollutant as follows: 4.2% (95% CI: 1.8 to 6.5) per IQR 30-min CO (10 ppm), 4.1% (95% CI: 1.8 to 6.3) per IQR 60-min CO2 (104 ppm), 3.9% (95% CI: 1.7 to 6.1) per IQR 60-min NO2 (155 ppb), and 4.4% (95% CI: 2.2 to 6.5) per IQR 30-min NOx (240 ppb). For this same group of four pollutants, no statistically significant associated were observed for HF, LF, or the LF/HF ratio, with the exception of significance between CO2 and HF and LF.
Across all HRV outcomes and pollutants, study day and time of day were significant predictors of the outcome. Ethnic origin and age were sometimes significant, but gender and smoking status were never significant predictors.
We investigated potential confounding by the gaseous pollutants O3, CO, and NO2 by including each of them individually in a two-pollutant model with PM2.5 (Table 5). The largest condition number (a measure of multicollinearity) for the fixed effects in the two-pollutant models was 18.05 (for PM2.5 and CO with a 30-min averaging period), and the values across all models were similar in magnitude. Condition numbers of this size indicate multicollinearity could be having a weak effect on the coefficient estimates . However, since a value greater than 100 is often used as a benchmark for significant multicollinearity, estimation and inference likely wasn’t affected in this instance.
In general, adjusting for O3 had little effect on the PM2.5 estimates across all HRV outcomes. The magnitudes of the PM2.5 estimates tended to fluctuate only slightly and the directions never changed. The significance of the results did not change with the exception of the 90-min window and SDNN (which became significant in the two-pollutant model).
The findings after adjusting for CO were similar. The PM2.5 magnitudes did not change notably and the directions of the estimates remained consistent. For SDNN 5-min and LF 60-min, the significance of the PM2.5 estimates in the two-pollutant models differed from those in the single-pollutant, while all others remained unchanged.
Lastly, adjusting for NO2, the PM2.5 estimate associated with the 5-min window and SDNN was not significant in the single-pollutant model, but significant in the two-pollutant model. No other results changed. As with the other gaseous confounders, adjusting for NO2 had little or no effect on the size or direction of the PM2.5 estimates.
To evaluate the influence of outlying exposure values on reported associations, we removed observations containing the highest and lowest 5% of pollutant concentrations (10% total). An average of 322 data points were removed for a given model, while an average of 2988 remained (across all individuals) (Table 6). Within a given individual, the smallest number of observations was 20 for this analysis.
The associations between PM2.5 and HF for the 30-, 60-, and 90-min averaging periods were smaller in magnitude and no longer significant after removal of the high/low values, suggesting these associations may be influenced by extreme exposure values. For other pollutants, the reported IQR systematically decreased, but the direction and significance of the estimates generally remained the same.
Short-term exposure to traffic-related emissions was associated with significant acute changes in HRV in this panel study of researchers that participated in the Mexico City Air Pollution Campaign . Gaseous pollutants – particularly ozone – were associated with reductions in time and frequency domain components. In contrast, PM2.5 increased these HRV parameters. Like Riediker , our results show a positive association between PM2.5 (mean PM2.5 mass concentration = 23.0 μg/m3) and two frequency domain HRV parameters (HF and LF), a result contrary to those observed in multiple studies of elderly populations exposed to PM air pollution related to traffic [13, 43, 44] (median PM2.5 mass concentration in each study, respectively = 10 μg/m3, 8.92 μg/m3 (2-hr), 7.7 μg/m3 (5 minute). Our results also demonstrate a negative association between the LF/HF ratio and PM2.5, CO, and CO2 that is significant for the 30-min averaging period. When assessed in normalized units (i.e. ratio), LF and HF provide quantitative indicators of neural control of the sinoatrial node and provide a synthetic index of the sympathovagal balance , which may be predictive of the development of ventricular arrhythmias . Our results, similar to other studies discussed below, suggest that the air pollution-associated acute changes in HRV parameters (increased or reduced) that have the highest potential for increasing the likelihood of a subsequent cardiac arrhythmia remain to be determined. Nonetheless, our results have significant implications to our understanding of how air pollution leads to an increase risk of cardiac arrhythmias. While having a reduced HRV is a risk factor for increased cardiovascular mortality; it is yet to be shown whether short-term environmental exposures associated with acute HRV changes can lead to cardiac arrhythmias in humans. This pathway, however, seems biologically plausible given that a) acute reductions or increases in HRV parameters have been associated with the onset of ventricular tachycardia [47–49], and b) short-term increases in ambient pollutants (both particles and gases) have been associated with increased likelihood of having a discharge from an implantable cardiac defibrillator .
Wu et al.  measured real-time, in-vehicle, traffic-related PM2.5 (56.6 μg/m3 (daily average)), and gaseous co-pollutants (CO, NO2, and NO) in a young population (n = 11, mean age = 35.5 years). They showed that IQR increases in PM2.5 mass concentrations (5–240 min moving averages) were associated with declines in three, 5-min HRV indices (SDNN, LF, and HF). Results, however, from their regression models for each subject showed heterogeneity among responses, i.e. several subjects in the study had positive associations with traffic-related PM exposures for the three, 5-min HRV indices. Further, their smoothed curves showing associations between PM exposure and 5-min HRV indices indicated that lower PM exposures were associated with increases in HRV, whereas higher PM exposures were associated with decreases in HRV. Similarly, to explore potential confounding by the gaseous co-pollutants, we separately included CO, NOx, and NO in a two-pollutant model with PM2.5. While adjusting for the co-pollutants decreased the precision of the estimates, our overall results were generally consistent with estimates not adjusted for co-pollutants. Overall, the results of Wu et al., like ours, suggest that differences in HRV response to traffic-related PM pollution may be related to differences in exposure levels, though factors impacting the heterogeneity of responses remain unclear.
Other studies in young populations illustrate the heterogeneity of responses between PM exposure and HRV. In a study of mail carriers between 25 and 46 years of age in Taiwan, no significant HRV effects from PM (mean sample time of approximately four hours using a personal cascade impactor sampler with a pump) were observed despite reported higher traffic-related air pollution exposure levels as compared to other studies (median PM2.5 exposure = 61.3 μ g/m3) . Similarly, PM2.5 exposure assessed by a light-scattering method (pDR real time instrument, 0–180 min averaging time) did not affect the SDNN measurements in 40 young, healthy residents of the Mexican metropolitan area . In a population of nine young highway patrol officers (mean age = 27.3 years), Riediker et al.  showed that in-vehicle exposure to PM2.5 (measured both by gravimetric and light-scattering methods) was associated with all time domain parameters (% difference between adjacent normal RR intervals that are greater than 50 msec or PNN50, SDNN, and mean cycle length, a phrase to emphasize that the interval between consecutive beats, rather than the heart rate, is being analyzed), as well as HF power and the power ratio LF/HF on the morning after the shift. In the occupational literature, Magari et al.  reported statistically significant associations using a two-hour lagged mean heart rate. Specifically, the PM2.5 average air concentration (measured using a light scattering instrument, mean 1160 μg/m3) showed an average increase of 1.67 msec (95% CI: 0.11 to 3.22) in the mean heart rate, for every 1 mg/m3 increase in the average PM2.5 concentration. The variety of exposure assessment techniques and corresponding averaging times may also contribute to the perceived heterogeneity of responses.
In addition to different PM exposure metrics, studies have used various methods to estimate the impact of gaseous co-pollutants. Some have used fixed-site monitors [13, 20, 30, 43] and/or time activity data  to estimate traffic-related ambient air pollution exposure, including gaseous co-pollutants. Our results showing acute reductions in HRV in association with ozone and other gaseous pollutants is in agreement with other recent studies .
In a repeated measures study of 46 subjects (43–75 years of age), Zanobetti et al.  suggested that pollutant mixtures may influence cardiac tone as both PM and O3 (30–120 min moving averages) had independent associations with reduced HRV in two-pollutant models. Using a five minute time resolution to evaluate the acute effects of residential outdoor ozone exposure and HRV changes, Jia et al.  showed in 20 elderly subjects that, after adjusting for other pollutants and subject characteristics, there was a reduction in the HF component of 4.87% (95% CI:0.97 to 8.62) per 10 ppb increment of O3. A similar result was observed in a study of patients recently discharged from the hospital for acute coronary disease-related complications. Here, the two-hour and five-day O3 moving averages were associated with reductions in time domain components of HRV indicative of parasympathetic function, whereas NO2 was associated with reductions in the HF spectral component .
A major limitation of these studies (as referenced in both) was the use of fixed-site monitors to estimate personal exposure measurements. Suh and Zanobetti  showed that changes in HRV – especially those associated with parasympathetic control – were significantly and negatively associated with elemental carbon, and, to a lesser degree, NOx when measurements of personal exposure (but not ambient, outdoor or indoor concentrations) were used to estimate their exposures. Interestingly and importantly, associations between personal exposure measurements and HRV were detectable only for these traffic-related pollutants. Non-significant findings with HRV were detected for 24-hour ambient concentrations and personal exposures to more spatially uniform regional pollutants, i.e. PM2.5 and O3.
Similar results were reported in a large controlled-exposure study examining the association between HRV and combined exposure to concentrated ambient particles (CAPs) and O3. In the participating group of healthy young adults (n = 50, mean age = 27.1 years), no consistent pattern of changes in HRV indices was detected among PM2.5 and O3 exposure categories . Despite the absence of a clear trend in the categorical exposure models, the dose–response analysis demonstrated a trend toward a negative linear association between CAPs mass concentration and change in several HRV indices, with a statically significant relationship for the LF HRV measure. No relationship existed without accounting for O3.
Given the similarities between these studies with our results, it is possible to speculate that, regardless of age or underlying cardiovascular disease, acute ozone exposures are associated with reductions in HRV. With respect to acute PM2.5 exposures, there are notable age-dependent and/or possibly underlying cardiovascular diseases that may influence the HRV responses to PM. In addition, acute changes in HRV (even those characterized by an increase in the HF domain) can be associated with the onset of a ventricular arrhythmia .
To the best of our knowledge, the literature contains varied to little information on the relationship between HRV parameters and potential associations between NOx, CO2, and formaldehyde. The selection of exposure measures reported here was not based on an a priori reason that some evidence existed to support association between exposure and HRV. Rather we evaluated potential associations within available exposures as part of the Mexico City Air Pollution Campaign. In 2004, the Multi-Ethnic Study of Atherosclerosis and Air Pollution was launched to investigate the relation between individual-level estimates of long-term air pollution exposure and the progression of subclinical atherosclerosis and the incidence of cardiovascular disease and includes exposure assessments of fine particulate matter, NOx, and black carbon; the majority of data collection will be completed in 2014 . Some have shown that certain nitrogen oxides, e.g. NO2, do not affect heart rate variability at concentrations high for urban background levels and in the absence of other pollutants . Others have shown that ambient NO2 concentrations were inversely associated with SDNN and positively associated with LF/HF. (β = 1.4; 95% CI, 0.35 to 2.5) 2 hr after the start of cycling . Our results demonstrate that both NO2 and NOx are significantly associated with SDNN. No studies were located that demonstrated a cardiovascular effects in humans after inhalation exposure to CO2 or formaldehyde . In our case, it’s likely that CO2 and/or formaldehyde exposures were acting as a surrogate for other vehicular gaseous pollutants.
One limitation of this work is that this analysis did not control for stress resulting from traffic. Psychological stress can influence both HRV and autonomic function . Since we did not measure participant stress while in traffic, we cannot disentangle the impact of both pollution and stress exposures resulting from being in traffic. In addition to traffic-related stress, participants may have observed the measured pollutant concentrations while being transported in the van-based mobile laboratory. Knowledge of their exposures may have stimulated the participants' sympathetic tone. Secondly, any exposure misclassification can lead to biased estimates of exposure-response estimates, particularly in cases with multiple correlated exposures where the direction of the bias is uncertain . While exposure error biases are often assumed to bias toward the null, a more complicated situation arises in cases like ours when two or more exposures are measured with error and are correlated with each other. This may lead to bias in both directions and with varying degree. A final caveat is that chemical exposure factors are known to affect HRV [61, 62], but these exposures are unlikely to be present in our study and thus, would be unlikely to affect our findings.
This study revealed an association between traffic-related PM2.5 exposure and acute changes in HRV in a population aged 22 to 56 years when PM2.5 exposures were relatively low (14 μ g/m3). Results also demonstrate heterogeneity in the effects of the different pollutants, with declines in HRV – especially HF – with ozone and formaldehyde exposures, and increases in HRV with PM2.5 exposure. Our findings support the need for additional research emphasizing the impacts of co-pollutants . Given that exposure to traffic-related emissions has been associated with increased risk of cardiovascular morbidity and mortality, understanding the mechanisms by which traffic-related emissions can cause cardiovascular disease has significant public health relevance.
Concentrated ambient particles
Power in high frequency
Heart rate variability
National Ambient Air Quality Standards
Fine particulate matter
Standard deviation of normal-to-normal intervals
Tunable infrared laser differential absorption spectrometer.
Lipsett MJ, Ostro BD, Reynolds P, Goldberg D, Hertz A, Jerrett M, Smith DF, Garcia C, Chang ET, Bernstein L: Long-term Exposure to Air Pollution and Cardiorespiratory Disease in the California Teachers Study Cohort. Am J Respir Crit Care Med. 2011, 184: 828-835. 10.1164/rccm.201012-2082OC.
Polichetti G, Cocco S, Spinali A, Trimarco V, Nunziata A: Effects of particulate matter (PM(10), PM(2.5) and PM(1)) on the cardiovascular system. Toxicology. 2009, 261: 1-8. 10.1016/j.tox.2009.04.035.
Pope CA, Burnett RT, Thurston GD, Thun MJ, Calle EE, Krewski D, Godleski JJ: Cardiovascular mortality and long-term exposure to particulate air pollution: epidemiological evidence of general pathophysiological pathways of disease. Circulation. 2004, 109: 71-77.
Chuang KJ, Coull BA, Zanobetti A, Suh H, Schwartz J, Stone PH, Litonjua A, Speizer FE, Gold DR: Particulate air pollution as a risk factor for ST-segment depression in patients with coronary artery disease. Circulation. 2008, 118: 1314-1320. 10.1161/CIRCULATIONAHA.108.765669.
Gan WQ, Koehoorn M, Davies HW, Demers PA, Tamburic L, Brauer M: Long-term exposure to traffic-related air pollution and the risk of coronary heart disease hospitalization and mortality. Environ Health Perspect. 2011, 119: 501-507.
Van Hee VC, Adar SD, Szpiro AA, Barr RG, Bluemke DA, Diez Roux AV, Gill EA, Sheppard L, Kaufman JD: Exposure to traffic and left ventricular mass and function: the Multi-Ethnic Study of Atherosclerosis. Am J Respir Crit Care Med. 2009, 179: 827-834. 10.1164/rccm.200808-1344OC.
Dockery DW: Epidemiologic evidence of cardiovascular effects of particulate air pollution. Environ Health Perspect. 2001, 109 (Suppl 4): 483-486.
Goldberg MS, Burnett RT, Yale JF, Valois MF, Brook JR: Associations between ambient air pollution and daily mortality among persons with diabetes and cardiovascular disease. Environ Res. 2006, 100: 255-267. 10.1016/j.envres.2005.04.007.
Larrieu S, Jusot JF, Blanchard M, Prouvost H, Declercq C, Fabre P, Pascal L, Tertre AL, Wagner V, Riviere S: Short term effects of air pollution on hospitalizations for cardiovascular diseases in eight French cities: the PSAS program. Sci Total Environ. 2007, 387: 105-112. 10.1016/j.scitotenv.2007.07.025.
Brook RD, Franklin B, Cascio W, Hong Y, Howard G, Lipsett M, Luepker R, Mittleman M, Samet J, Smith SC, Tager I: Air pollution and cardiovascular disease: a statement for healthcare professionals from the Expert Panel on Population and Prevention Science of the American Heart Association. Circulation. 2004, 109: 2655-2671. 10.1161/01.CIR.0000128587.30041.C8.
Kunzli N, Jerrett M, Garcia-Esteban R, Basagana X, Beckermann B, Gilliland F, Medina M, Peters J, Hodis HN, Mack WJ: Ambient air pollution and the progression of atherosclerosis in adults. PLoS One. 2010, 5: 1-10.
Wu S, Deng F, Niu J, Huang Q, Liu Y, Guo X: Association of heart rate variability in taxi drivers with marked changes in particulate air pollution in Beijing in 2008. Environ Health Perspect. 2010, 118: 87-91.
Schwartz J, Litonjua A, Suh H, Verrier M, Zanobetti A, Syring M, Nearing B, Verrier R, Stone P, MacCallum G: Traffic related pollution and heart rate variability in a panel of elderly subjects. Thorax. 2005, 60: 455-461. 10.1136/thx.2004.024836.
Anonymous: Heart rate variability: standards of measurement, physiological interpretation and clinical use. Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology. Circulation. 1996, 93: 1043-1065. 10.1161/01.CIR.93.5.1043.
Gold DR, Litonjua A, Schwartz J, Lovett E, Larson A, Nearing B, Allen G, Verrier M, Cherry R, Verrier R: Ambient pollution and heart rate variability. Circulation. 2000, 101: 1267-1273. 10.1161/01.CIR.101.11.1267.
Park SK, O'Neill MS, Vokonas PS, Sparrow D, Schwartz J: Effects of air pollution on heart rate variability: the VA normative aging study. Environ Health Perspect. 2005, 113: 304-309. 10.1289/ehp.113-a304.
Stein PK, Kleiger RE: Insights from the study of heart rate variability. Annu Rev Med. 1999, 50: 249-261. 10.1146/annurev.med.50.1.249.
Fan ZT, Meng Q, Weisel C, Laumbach R, Ohman-Strickland P, Shalat S, Hernandez MZ, Black K: Acute exposure to elevated PM2.5 generated by traffic and cardiopulmonary health effects in healthy older adults. J Expo Sci Environ Epidemiol. 2009, 19: 525-533. 10.1038/jes.2008.46.
Rowan WH, Campen MJ, Wichers LB, Watkinson WP: Heart rate variability in rodents: uses and caveats in toxicological studies. Cardiovasc Toxicol. 2007, 7: 28-51. 10.1007/s12012-007-0004-6.
Chuang KJ, Chan CC, Su TC, Lee CT, Tang CS: The effect of urban air pollution on inflammation, oxidative stress, coagulation, and autonomic dysfunction in young adults. Am J Respir Crit Care Med. 2007, 176: 370-376. 10.1164/rccm.200611-1627OC.
Fakhri AA, Ilic LM, Wellenius GA, Urch B, Silverman F, Gold DR, Mittleman MA: Autonomic effects of controlled fine particulate exposure in young healthy adults: effect modification by ozone. Environ Health Perspect. 2009, 117: 1287-1292.
Peretz A, Kaufman JD, Trenga CA, Allen J, Carlsten C, Aulet MR, Adar SD, Sullivan JH: Effects of diesel exhaust inhalation on heart rate variability in human volunteers. Environ Res. 2008, 107: 178-184. 10.1016/j.envres.2008.01.012.
Magari SR, Schwartz J, Williams PL, Hauser R, Smith TJ, Christiani DC: The association of particulate air metal concentrations with heart rate variability. Environ Health Perspect. 2002, 110: 875-880. 10.1289/ehp.02110875.
Holguin F, Tellez-Rojo MM, Hernandez M, Cortez M, Chow JC, Watson JG, Mannino D, Romieu I: Air pollution and heart rate variability among the elderly in Mexico City. Epidemiology. 2003, 14: 521-527. 10.1097/01.ede.0000081999.15060.ae.
Delfino RJ, Gillen DL, Tjoa T, Staimer N, Polidori A, Arhami M, Sioutas C, Longhurst J: Electrocardiographic ST-segment depression and exposure to traffic-related aerosols in elderly subjects with coronary artery disease. Environ Health Perspect. 2011, 119: 196-202. 10.1289/ehp.1103645.
Gold DR, Litonjua AA, Zanobetti A, Coull BA, Schwartz J, MacCallum G, Verrier RL, Nearing BD, Canner MJ, Suh H, Stone PH: Air pollution and ST-segment depression in elderly subjects. Environ Health Perspect. 2005, 113: 883-887. 10.1289/ehp.7737.
Laden F, Neas LM, Dockery DW, Schwartz J: Association of fine particulate matter from different sources with daily mortality in six U.S. cities. Environ Health Perspect. 2000, 108: 941-947. 10.1289/ehp.00108941.
Riediker M, Cascio WE, Griggs TR, Herbst MC, Bromberg PA, Neas L, Williams RW, Devlin RB: Particulate matter exposure in cars is associated with cardiovascular effects in healthy young men. Am J Respir Crit Care Med. 2004, 169: 934-940. 10.1164/rccm.200310-1463OC.
NRC: Research Priorities for Airborne Particulate Matter IV: Continuing Research Progress. 2004, Washington, DC: National Academies Press
Jia X, Song X, Shima M, Tamura K, Deng F, Guo X: Acute effect of ambient ozone on heart rate variability in healthy elderly subjects. J Expo Sci Environ Epidemiol. 2011, 21 (5): 541-547. 10.1038/jes.2011.18.
Laumbach RJ, Rich DQ, Gandhi S, Amorosa L, Schneider S, Zhang J, Ohman-Strickland P, Gong J, Lelyanov O, Kipen HM: Acute changes in heart rate variability in subjects with diabetes following a highway traffic exposure. J Occup Environ Med. 2010, 52: 324-331. 10.1097/JOM.0b013e3181d241fa.
HEI: Traffic-Related Air Pollution: A Critical Review of the Literature on Emissions, Exposure, and Health Effects. Book Traffic-Related Air Pollution: A Critical Review of the Literature on Emissions, Exposure, and Health Effects. 2010, City: Health Effects Institute
Report on the Mexico City Metropolitan Area (MCMA). 2003, [http://mce2.org/newsletter/nwsltr_3/english/field.html], Field Measurement Campaign
Kolb CEHS, McManus JB, Shorter JH, Zahniser MS, Nelson DD, Jayne JT, Canagaratha MR, Worsnop MDR: Mobile laboratory with rapid response instruments for real-time measurments of urban and regional trace gas and particulate distributions and emission source characteristics. Environ Sci Technol. 2004, 38: 5694-5703. 10.1021/es030718p.
Marr LC, Grogan LA, Wohrnschimmel H, Molina LT, Smith TJ, Garshick E: Vehicle traffic as a source of particulate polycyclic aromatic hydrocarbon exposure in the Mexico City metropolitan area. Environ Sci Technol. 2004, 38: 2584-2592. 10.1021/es034962s.
Jiang MML, Dunlea EJ: Vehicle flight emissions of black carbon, polycyclic aromatic hydrocarbons, and other pollutants measured by a mobile laboratory in Mexico City. Atmos Chem Phys. 2005, 5: 3377-3387. 10.5194/acp-5-3377-2005.
Dunlea EJ HS, Nelson DD, Volkamer RM, Lamb BK, Allwine EJ, Grutter M, Ramos Villegas CR, Marquex C, Blanco S, Cardenas B, Kolb CE, Molina LT, Molina MJ: Technical note: Evaluation of standard ultraviolet absorption ozone monitors in a polluted urban environment. Atmos Chem Phys. 2006, 6: 3163-3180. 10.5194/acp-6-3163-2006.
Jiang M, Marr LC, Dunlea EJ, Herndon SC, Jayne JT, Kolb CE, Knighton WB, Rogers TM, Zavala M, Molina LT, Molina MJ: Vehicle fleet emissions of black carbon, polycyclic aromatic hydrocarbons, and other pollutants measured by a mobile laboratory in Mexico City. Atmos Chem Phys. 2005, 5: 3377-3387. 10.5194/acp-5-3377-2005.
Stein PK, Bosner MS, Kleiger RE, Conger BM: Heart rate variability: a measure of cardiac autonomic tone. Am Heart J. 1994, 127: 1376-1381. 10.1016/0002-8703(94)90059-0.
Belsley DA, Kuh E, Welsch RE: Regression Diagnostics. 1980, New York: John Wiley & Sons
Fitzmaurice G, Laird N, Ware J: Applied longitudinal anlysis. 2004, New York: John Wiley & Sons
National Ambient Air Quality Standards (NAAQS). [http://www.epa.gov/air/criteria.html]
Zanobetti A, Gold DR, Stone PH, Suh HH, Schwartz J, Coull BA, Speizer FE: Reduction in heart rate variability with traffic and air pollution in patients with coronary artery disease. Environ Health Perspect. 2010, 118: 324-330.
Adar SD, Gold DR, Coull BA, Schwartz J, Stone PH, Suh H: Focused exposures to airborne traffic particles and heart rate variability in the elderly. Epidemiology. 2007, 18: 95-103. 10.1097/01.ede.0000249409.81050.46.
Malliani A, Lombardi F, Pagani M, Cerutti S: Power Spectral-Analysis of Cardiovascular Variability in Patients at Risk for Sudden Cardiac Death. J Cardiovasc Electr. 1994, 5: 274-286. 10.1111/j.1540-8167.1994.tb01164.x.
Vanoli E, Adamson PB, Foreman RD, Schwartz PJ: Prediction of unexpected sudden death among healthy dogs by a novel marker of autonomic neural activity. Heart Rhythm. 2008, 5: 300-305. 10.1016/j.hrthm.2007.10.021.
Lerma C, Wessel N, Schirdewan A, Kurths J, Glass L: Ventricular arrhythmias and changes in heart rate preceding ventricular tachycardia in patients with an implantable cardioverter defibrillator. Med Biol Eng Comput. 2008, 46: 715-727. 10.1007/s11517-008-0326-y.
Lombardi F, Porta A, Marzegalli M, Favale S, Santini M, Vincenti A, De Rosa A: Heart rate variability patterns before ventricular tachycardia onset in patients with an implantable cardioverter defibrillator. Participating Investigators of ICD-HRV Italian Study Group. Am J Cardiol. 2000, 86: 959-963. 10.1016/S0002-9149(00)01130-9.
Mani V, Wu X, Wood MA, Ellenbogen KA, Hsia PW: Variation of spectral power immediately prior to spontaneous onset of ventricular tachycardia/ventricular fibrillation in implantable cardioverter defibrillator patients. J Cardiovasc Electrophysiol. 1999, 10: 1586-1596. 10.1111/j.1540-8167.1999.tb00223.x.
Link MS, Dockery DW: Air pollution and the triggering of cardiac arrhythmias. Am J Epidemiol. 176 (9): 827-837.
Wu CF, Kuo IC, Su TC, Li YR, Lin LY, Chan CC, Hsu SC: Effects of personal exposure to particulate matter and ozone on arterial stiffness and heart rate variability in healthy adults. Am J Epidemiol. 2010, 171: 1299-1309. 10.1093/aje/kwq060.
Vallejo M, Ruiz S, Hermosillo AG, Borja-Aburto VH, Cardenas M: Ambient fine particles modify heart rate variability in young healthy adults. J Expo Sci Environ Epidemiol. 2006, 16: 125-130. 10.1038/sj.jea.7500447.
Wheeler A, Zanobetti A, Gold DR, Schwartz J, Stone P, Suh HH: The relationship between ambient air pollution and heart rate variability differs for individuals with heart and pulmonary disease. Environ Health Perspect. 2006, 114: 560-566.
Suh HH, Zanobetti A: Exposure error masks the relationship between traffic-related air pollution and heart rate variability. J Occup Environ Med. 2010, 52: 685-692. 10.1097/JOM.0b013e3181e8071f.
Kaufman JD, Adar SD, Allen RW, Barr RG, Budoff MJ, Burke GL, Casillas AM, Cohen MA, Curl CL, Daviglus ML: Prospective Study of Particulate Air Pollution Exposures, Subclinical Atherosclerosis, and Clinical Cardiovascular Disease: The Multi-Ethnic Study of Atherosclerosis and Air Pollution (MESA Air). Am J Epidemiol. 2012, epub ahead of print
Scaife A, Barclay J, Hillis GS, Srinivasan J, Macdonald DW, Ross JA, Ayres JG: Lack of effect of nitrogen dioxide exposure on heart rate variability in patients with stable coronary heart disease and impaired left ventricular systolic function. Occup Environ Med. 2012, 69: 587-591. 10.1136/oemed-2011-100126.
Weichenthal S, Kulka R, Dubeau A, Martin C, Wang D, Dales R: Traffic-related air pollution and acute changes in heart rate variability and respiratory function in urban cyclists. Environ Health Perspect. 2011, 119: 1373-1378. 10.1289/ehp.1003321.
Toxicological Profile for Formaldehyde. [http://www.atsdr.cdc.gov/toxprofiles/TP.asp?id=220&tid=39]
Anderson JL, Adams CD, Antman EM, Bridges CR, Califf RM, Casey DE, Chavey WE, Fesmire FM, Hochman JS, Levin TN: 2011 ACCF/AHA Focused Update Incorporated Into the ACC/AHA 2007 Guidelines for the Management of Patients With Unstable Angina/Non-ST-Elevation Myocardial Infarction: a report of the American College of Cardiology Foundation/American Heart Association Task Force on Practice Guidelines. Circulation. 2011, 123: e426-579. 10.1161/CIR.0b013e318212bb8b.
Zeka A, Schwartz J: Estimating the independent effects of multiple pollutants in the presence of measurement error: an application of a measurement-error-resistant technique. Environ Health Perspect. 2004, 112: 1686-1690. 10.1289/ehp.7286.
Bae S, Kim JH, Lim YH, Park HY, Hong YC: Associations of bisphenol A exposure with heart rate variability and blood pressure. Hypertension. 2012, 60: 786-793. 10.1161/HYPERTENSIONAHA.112.197715.
Togo F, Takahashi M: Heart rate variability in occupational health –a systematic review. Ind Health. 2009, 47: 589-602. 10.2486/indhealth.47.589.
Research Priorities for Airborne Particulate Matter IV: Continuing Research Progress. Edited by: Council NR. 2004, Washington DC: National Academies Press
Authors LM, MM and FH were funded through the Mexico City Project at the Massachusetts Institute of Technology.
The authors declare that they have no completing interests.
KNS drafted the manuscript. JC and MOH performed the statistical analysis. ML performed the field work. LM and MM provided logistical and financial support. FH conceived of the study and participated in its design and coordination. All authors read and approved the final manuscript.