ORIGINAL RESEARCH article

Front. Environ. Sci., 25 January 2018

Sec. Air Pollution

Volume 6 - 2018 | https://doi.org/10.3389/fenvs.2018.00002

Exposure to Fine Particulate Matter Leads to Rapid Heart Rate Variability Changes

  • 1. Institute of Occupational Medicine Singapore, Singapore, Singapore

  • 2. Institut Universitaire Romand de Santé au Travail, University of Lausanne, Lausanne, Switzerland

  • 3. School of Materials Science and Engineering, Nanyang Technological University, Singapore, Singapore

  • 4. Institut Universitaire de Médecine Sociale et Préventive, University of Lausanne, Lausanne, Switzerland

  • 5. Federal Office for the Environment, Bern, Switzerland

Abstract

Introduction: Heart Rate Variability (HRV) reflects the adaptability of the heart to internal and external stimuli. Reduced HRV is a predictor of post-infarction mortality. We previously found in road maintenance workers HRV-increases several hours after exposure to fine particulate matter (PM2.5). This seemed to conflict with studies where PM-exposure acutely reduced HRV. We therefore assessed whether time from exposure to HRV-assessment could explain the differences observed.

Methods: On five non-consecutive days, workers carried nephelometers providing 1-min-interval PM2.5-exposure. Five-min HRV-intervals of SDNN (Standard Deviation of Normal to Normal beat intervals) and pNN50 (Percentage of the interval differences exceeding 50 ms) were extracted from 24-h electrocardiograms (ECGs). Following 60 min PM2.5-exposure, changes in HRV-parameters were assessed during 120-min visually and by regression analysis with control for time at work, at home, and during the night using autoregressive integrating moving average (ARIMA) models to account for autocorrelation of the time-series. Additional controls included changing the time windows and including body mass index (BMI) and age in the models.

Result: Pattern analysis of 12,669 data points showed high modulation of mean, standard deviation (SD), and time trend of HRV (SDNN and pNN50) at low, and much reduced modulation at high PM2.5-exposures. The time trend following exposure was highly symmetrical, resembling a funnel plot. Regression analysis showed significant associations of decreasing SDNN and pNN50 (average, SD, and absolute value of time trend) with increasing PM2.5-exposure, which remained significant when controlling for activity phases. Changing time windows did not change the pattern of response. Including BMI and age did not change the results.

Conclusions: The reduced modulation of HRV following PM2.5-exposure is striking. It suggests strong interference with homeostatic controls. Such an interference would represent a serious bodily burden, and could help explain acute cardiac events. In this model, the increase of HRV several hours later would reflect a recovery response.

Introduction

Heart Rate Variability (HRV) describes the variation of the timing between individual heart beats. It is a consequence of the heart's response to a series of internal and external stimuli such as respiration, baroreceptor feedback loop, and parasympathetic and sympathetic nervous system activity (Billman, ; Shaffer et al., ). Elevated HRV is a sign for a healthy cardiovascular system, while reduced HRV is a predictor of post-infarct mortality (Billman, ), cardiac events in healthy individuals (Tsuji et al., ), and incident heart failure in older at-risk adults (Patel et al., ). Most commonly assessed are the time- or the frequency-domain of HRV, though a large number of mathematical and statistical approaches are available (Electrophysiology Taskforce Guidelines, ; Sassi et al., ). Amongst the frequently used time domain parameters are SDNN (Standard Deviation of Normal to Normal beat intervals) and pNN50 (Percentage of the interval differences exceeding 50 ms): SDNN was found to reflect overall autonomic activity and cardiovascular fitness (Electrophysiology Taskforce Guidelines, ; Tsuji et al., ; Da Silva et al., ), While PNN50 is frequently used to assess parasympathetic activity (Ewing et al., ; Electrophysiology Taskforce Guidelines, ; Tsuji et al., ; Mietus et al., ).

Fine particulate matter (PM) on roads originates to an important part from traffic and combustion engines which release particles and precursors of secondary aerosols into the air (Riediker et al., ; Putaud et al., ; Meier et al., ). High concentrations occur in particular along busy roads such as highways, which are depending on factors such as traffic counts but also weather and wind (Imhof et al., ; Zwack et al., ). Workers maintaining the highways will experience exposure to PM2.5 that originates from the general background, from the passing traffic and also the different work tasks (Meier et al., ). Exposure to airborne PM is a well-established risk for human health (Brook et al., ; Landrigan et al., ), in particular fine (PM2.5) and ultrafine PM increases both short- and long-term cardiovascular risk. PM negatively affects cardiovascular risk factors (reviewed by Brook et al., ) such as inflammation (Tsai et al., ), sodium regulation by the kidney (Tsai et al., ), pulse pressure (Tsai et al., ), and metabolic syndrome (Brook et al., ; Devlin et al., ). PM is associated with cardiac arrhythmias (Kowalska and Kocot, ; Mordukhovich et al., ; Song et al., ), arrhythmic events are more frequently detected and treated by implanted cardioverter-defibrillators (Dockery et al., ; Metzger et al., ; Link et al., ), and people spending time in traffic are more likely to suffer a heart attack within the hour (Peters et al., ). In line with these ill-health effects of PM was the observation that PM2.5 exposure is associated with reduced HRV: in controlled exposure of rats to PM2.5, HRV was reduced during ongoing exposure (Wagner et al., ). In a series of epidemiological (Pope et al., ; Wu et al., ; Huang et al., ; Nyhan et al., ; Mordukhovich et al., ; Lee et al., ) as well as controlled PM-exposure (Devlin et al., ; Weichenthal et al., ; Cole-Hunter et al., ; Hemmingsen et al., ) studies, small reductions in HRV within an hour after PM exposure were observed. However, a recent meta-analysis of panel studies investigating HRV effects of PM2.5 (Buteau and Goldberg, ) questioned the association between PM2.5 and frequently used indices of HRV. Indeed, in our own two studies in healthy workers, we had found increased HRV (SDNN and PNN50) in the morning after exposure to traffic particles (Riediker et al., ,; Riediker, ; Meier et al., ). In the study with highway patrol troopers, we had also assessed an intermediary point at the end of the working day, which was well after peak rush hour exposure (Riediker et al., ). At that intermediary point, no associations between PM2.5 and HRV were observed. The time from exposure to HRV analysis varies strongly between studies, an aspect that was not further assessed in the aforementioned review (Buteau and Goldberg, ). In this study, we aimed to understand if differences in time from exposure until the HRV-response assessment could explain the differences in HRV-association with PM2.5. For this purpose, we assessed the immediate HRV response to particles by analyzing the time-series of electrocardiograms (ECGs) and PM2.5 obtained in the highway maintenance workers study (Meier et al., ) where we had found increases of HRV in the morning after exposure. In addition, we aimed to better understand the way HRV was acutely changed and thus assessed (in addition to changes in average HRV) also the modulation of HRV over a 2-h time window following exposure and the time-trend of the series of 5-min averages of HRV during this time window.

Methods

We analyzed PM2.5-exposure and ECG data collected from a panel of 18 healthy male Swiss highway maintenance workers from 10 different maintenance centers in western Switzerland. The study was approved by the Ethics Committee of the University of Lausanne. All subjects gave written informed consent prior to participation in the study. The here presented analysis focused on the effects of continuously measured PM2.5 on HRV. The methodology to collect exposure (Meier et al., ) and health (Meier et al., ) data was previously described in detail. Briefly, subjects were equipped in the morning before the start of their work with small devices to record ECG and personal PM2.5-exposure for 24-h. Twelve ECG leads were attached by a trained researcher, then connected to a. digital ECG Holter monitor recorder (H12+ Digital, Mortara Instrument, Inc., Milwaukee, WI, USA).The high-resolution ECG records were afterwards processed with H-Scribe+ software, and inspected manually by an experienced cardiologist (Wayne E. Cascio, U.S. EPA, Chapel Hill, NC, USA) to remove artifacts and non-normal beats. Processed ECG-derived HRV data was recorded at 5-min resolution. PM2.5 was measured at 1-min resolution by light scattering nephelometers (DataRam pDR1000, Thermo Scientific, Waltham, MA, USA). This averaging time gives a precision of ±0.2% of reading or ±0.005 mg/m3, whichever is larger (manufacturer specification). We further improved the readout by correcting if for humidity (Richards, ). The monitor was attached to the body near the breathing zone while the worker was awake; and placed on the bed side table during sleep. The 1-min recordings were condensed into 5-min average intervals so that they matched the time intervals of the ECGs. Logbooks were used to identify activity phases at work, at home and at night, defined as the time between 22:00 and 06:00. Subjects participated up to five times on non-consecutive days.

Statistical analyses

From the time-series of PM2.5 recordings, moving average concentrations were calculated for 60-min time-windows. From the HRV-recordings (5-min interval SDNN and pNN50), the following summary variables were calculated for moving time-windows of 120 min: mean HRV-value, Standard Deviation (SD) of the HRV, the time-trend (expressed as the slope over time) of HRV during this time window and the absolute value of the time trend. These HRV summary-variables were compared to the preceding PM2.5-exposure (example shown in Figure 1). The time intervals were defined based on the hypothesis that effects would occur already within the first 2-h after exposure, and the attempt to have a sufficiently large number of data-points per interval to allow for a robust assessment. The consequences of changing the averaging intervals were tested in a sensitivity analysis. The influence of activity phase at the moment of exposure (at work, home, or night) was assessed both, by stratification and calculation of the effect estimates for each strata, and by including activity phase as a co-factor in the models. The ARIMA methodology was used to account for second-order autocorrelation of the time-series. The models were computed individually for each person and the resulting coefficients were combined using a meta-analysis approach to obtain the population coefficients and variances. To account for the skewed distribution, PM2.5 was log-transformed (in basis 2) when fitting the models. To enable a comparison among results involving different cardiac parameters, all these parameters were standardized for each individual by subtracting the mean from the individual observations and dividing the difference by the parameter's SD. The potential influence of age and body mass index (BMI) on the individuals' effect estimates was assessed by using Spearman rank test. Calculations were performed using the basic functions of the R statistical software package Forecast version 7.3.

Figure 1

Results

Data were collected on 50 days from 18 workers, resulting in 88 personal observation-days. On 77 of these days both, valid PM2.5 and HRV data were available to calculate HRV and PM2.5 parameters. The resulting 12,669 ECG data points and matching exposure data were used for analysis (summarized in Table 1). Mean PM2.5-exposure at work was about double than that at home and about four times the night values. SDNN during wake phases was slightly smaller than during the night, pNN50 was lowest at work and highest during the night.

Table 1

ParameterAll activity phases mean (range)Work only mean (range)Home only mean (range)Night only mean (range)
n1,19413,8873,6054,449
PM2.5 (60 min mean) (μg/m3)51.35 (0.026–3333)88.48 (0.368–3333)48.17 (0.026–2129)21.5 (0.026–148.1)
SDNN (120 min mean) (ms)72.44 (14.91–240.2)63.28 (22.52–158.30)63.4 (17.74–205.3)87.75 (14.91–240.2)
pNN50 (120 min mean) (%)13.32 (0–86.67)5.63 (0–56.38)11.04 (0–77)21.89 (0–86.67)

Summary statistics of assessed data.

Plotting mean and SD of 120 min of SDNN and pNN50 against 60 min of preceding PM2.5-exposure (Figure 2) showed a wide range of values at low concentrations and a smaller range at higher concentrations, while visually, there was no clear difference in mean value for different exposure levels. The SD of 5-min HRV over the 2-h following exposure showed a similar pattern as the mean (Figure 2). However, plotting the time trend of SDNN and pNN50 following the different exposure levels (Figure 2) revealed a very symmetrical pattern of time trend-values being highly variable at low exposure and close to zero at elevated exposures, yielding a typical appearance of a funnel plot. This pattern was highly significant, as revealed in the regression analysis below. A similar pattern was observed for all activity phases (work, home, night).

Figure 2

Regression analysis of SDNN and pNN50 parameters (mean, SD and absolute value of the time trend) against exposure (Table 2) showed significant negative PM2.5 effect estimates for all three parameters of SDNN and pNN50. Control for activity phases at work, at home and during the night reduced the effect estimates but all remained significant at the 5% significance level. Neither age nor BMI were statistically significantly related to the individual effect estimates. Stratified analyses also show significant effect estimates in most strata (Table 3).

Table 2

ModelStandardized SDNN effect estimates per log2(PM2.5)p-valueStandardized pNN50 effect estimates per log2(PM2.5)p-value
Crude: mean ~ PM2.5−1.449<0.001−0.9970.008
Controlled: mean ~ PM2.5−1.026<0.001−0.8290.012
Crude: SD ~ PM2.5−1.158<0.001−0.977<0.001
Controlled: SD ~ PM2.5−0.685<0.001−0.800<0.001
Crude: abs(time trend) ~ PM2.5−0.2790.003−0.3470.001
Controlled: abs(time trend) ~ PM2.5−0.1490.047−0.2840.002

Association of 1-h exposure to PM2.5 with mean, SD, and absolute value of the time trend of SDNN and pNN50 (controlled for activity phase) in the following 2-h.

Coefficients are standardized effect estimates per log2(PM2.5).

Table 3

ModelSDNN coefficientVariancep-valuepNN50 coefficientVariancep-value
WORK ONLY
Mean—PM2.5−0.5920.4400.0280.1390.0150.356
SD—PM2.5−0.4380.0830.013−0.4780.0100.002
Absolute slope—PM2.5−0.0890.0010.370−0.1670.0000.007
HOME ONLY
Mean—PM2.5−0.4100.9140.095−0.9400.5110.009
SD—PM2.5−0.6250.1200.000−0.6010.0370.003
Absolute slope—PM2.5−0.0600.0020.606−0.2670.0010.046
NIGHT ONLY
Mean—PM2.5−0.7951.8580.0320.3470.2280.344
SD—PM2.5−1.1430.5470.001−1.0950.0810.002
Absolute slope—PM2.5−0.1710.0050.163−0.6990.0010.000

Stratification by activity phase of the analysis of the association of short-term exposure to PM2.5 for mean, SD, and absolute value of the slope of SDNN and pNN50.

Coefficients are standardized effect estimates per log2(PM2.5).

In a sensitivity analysis, the influence of the averaging times was tested by shortening and increasing the intervals. The visual pattern remained similar at all tested interval lengths and combinations (examples of tested variants shown in Figure 3). The effect estimates (all tested combinations shown in Table 4) for mean and SD of SDNN and pNN50 remained significant in most tested variants, while the absolute value of the time trends of SDNN and pNN50 were significant only for the longer time windows after control for activity (at least 40 min of exposure, at least 80 min of HRV).

Figure 3

Table 4

Control for activityTime-interval:Exposure_HRVCoefficientvariancep-valuesCoefficientVariancep-values
Mean SDNNMean pNN50
Crude30_60−0.2770.1060.001−0.3010.0170.004
Crude40_40−0.2450.1370.000−0.1950.0200.031
Crude40_80−0.6520.2850.000−0.4790.0330.010
Crude50_100−0.9070.5820.001−0.5510.0560.037
Crude60_1201.4490.8000.0000.9970.1060.008
Crude60_60−0.5550.3620.000−0.4370.0500.020
Crude80_80−1.1150.8870.000−0.7320.1020.029
Crude120_120−2.6641.5050.000−2.4350.2090.000
Controlled30_60−0.1090.0740.128−0.1550.0080.042
Controlled40_40−0.1490.1140.020−0.1010.0130.201
Controlled40_80−0.3850.1770.006−0.2450.0150.066
Controlled50_100−0.5750.3640.013−0.2450.0150.066
Controlled60_1201.0260.4710.0000.8290.0570.012
Controlled60_60−0.2090.2960.150−0.1750.0320.270
Controlled80_80−0.7140.7530.014−0.3860.0640.171
Controlled120_120−2.0251.1400.000−1.9150.1280.000
SD SDNNSD pNN50
Crude30_60−0.2590.0350.000−0.2630.0070.001
Crude40_40−0.2100.0320.000−0.1700.0050.001
Crude40_80−0.4920.1070.000−0.4940.0110.000
Crude50_100−0.7580.1400.000−0.6990.0180.000
Crude60_1201.1580.1740.0000.9770.0280.000
Crude60_60−0.3780.1590.002−0.3770.0180.004
Crude80_80−0.7390.2080.000−0.7980.0320.000
Crude120_120−2.1110.3650.000−2.1680.0650.000
Controlled30_60−0.1360.0300.011−0.1420.0040.018
Controlled40_40−0.0990.0300.007−0.0640.0020.055
Controlled40_80−0.2600.0680.013−0.3070.0050.000
Controlled50_100−0.4340.0720.001−0.5330.0120.000
Controlled60_1200.6850.0950.0000.8000.0170.000
Controlled60_60−0.1830.1290.097−0.2500.0120.023
Controlled80_80−0.4170.1470.007−0.6110.0230.001
Controlled120_120−1.3650.2640.000−1.7380.0500.000
abs(slope SDNN)abs(slope pNN50)
Crude30_60−0.1030.0020.000−0.0710.0000.027
Crude40_40−0.0900.0030.000−0.0880.0000.000
Crude40_80−0.1590.0020.001−0.1330.0000.011
Crude50_100−0.2400.0020.000−0.2580.0000.001
Crude60_1200.2790.0020.0030.3470.0000.001
Crude60_60−0.1290.0030.000−0.0970.0000.017
Crude80_80−0.2050.0020.000−0.2370.0010.005
Crude120_120−0.3290.0030.003−0.6930.0010.000
Controlled30_60−0.0360.0010.059−0.0110.0000.635
Controlled40_40−0.0180.0020.259−0.0230.0000.173
Controlled40_80−0.0670.0010.070−0.0590.0000.149
Controlled50_100−0.1470.0010.004−0.1970.0000.004
Controlled60_1200.1490.0010.0470.2840.0000.002
Controlled60_60−0.0360.0010.140−0.0290.0000.354
Controlled80_80−0.0840.0010.032−0.1220.0000.054
Controlled120_120−0.1920.0010.025−0.5530.0010.001

Analysis of the sensitivity of the effect estimates to changing the averaging times of PM2.5-exposure and HRV-values, with and without control for activity phases.

The main model is highlighted in bold-italic. Coefficients are standardized effect estimates per log2(PM2.5).

Discussion

Acute exposure to PM2.5 was followed by reduced HRV, which is in agreement with previous studies (Pope et al., ; Wu et al., ; Huang et al., ; Nyhan et al., ; Mordukhovich et al., ; Lee et al., ). Different from those studies, we also assessed whether HRV shows a time-trend to increase or decrease after high exposures. The scatterplots showing the time trends of SDNN and pNN50 vs. PM2.5 have a remarkable appearance of a “funnel plot”: whereas we had approximately as many increasing and decreasing time trends, compensating each other and averaging zero, whatever the value of PM2.5, the values of the time trends (whether increasing or decreasing) were consistently distributed over a wide range for low exposures and distributed in a very narrow range for high exposures (Figure 2).

The pattern of reduced modulation of HRV at elevated PM2.5-levels exists for all three activity phases assessed (work, home, at night), though the pattern is easiest to be recognized visually for the working hours. This is remarkable because the exposure at home and at night is representing what the normal (Swiss) population experiences, which is quite low. In contrast, the working exposure is reflective of the workplace along the highways with occasionally very high exposure levels (Riediker et al., ; Meier et al., ). In many major cities in low-income countries around the globe similar. Concentrations are at similar or even higher levels than those found at the workplace of these Swiss highway maintenance teams (Landrigan et al., ). The here observed striking work-time effects can be expected in those city dwellers. Our findings raise serious questions about the long-term health of the populations of such cities, who are continuously exposed to such levels. However, note that our study was limited to a small sample of healthy male Swiss adults of good socioeconomic status. It would seem necessary to assess whether other population groups such as females, children, other ethnicities, unhealthy, and lower socioeconomic status may respond differently to the same level of exposure.

The decrease of both, average SDNN and PNN50 can be interpreted as a reduced autonomic activity (Ewing et al., ; Electrophysiology Taskforce Guidelines, ; Mietus et al., ). This type of acute response was reported by other groups before (Devlin et al., ; Pope et al., ; Wu et al., ; Huang et al., ; Nyhan et al., ; Wagner et al., ; Weichenthal et al., ; Cole-Hunter et al., ; Hemmingsen et al., ; Mordukhovich et al., ; Lee et al., ). Interestingly, SDNN at high exposure levels is just slightly below the average SDNN, while PNN50 drops to near zero, which is suggestive of mostly parasympathetic withdrawal. Our analysis also shows that the random modulation of HRV (SD and time trend of SDNN and PNN50) is strongly reduced following elevated PM2.5-exposure, suggesting that this interference with regulatory autonomic mechanisms lasts for a while. PM2.5 was previously reported to interfere with salt regulation (Tsai et al., ), vascular function (Hemmingsen et al., ), and blood pressure (Tsai et al., ). It would be interesting to experimentally study if this could be a possible explanation for the observed reduced modulation of HRV over the 2-h window. Acutely altered HRV often represents a disturbance of normal homeostatic processes (Billman, ), which can pose a serious health concern especially for people having elevated cardiovascular risk factors such as those with previous infarction or existing arrhythmic disorders (Tsuji et al., ; Clyne et al., ; Patel et al., ).

The acute decrease of HRV following exposure identified in the regression analysis confirms our hypothesis that changing the time from exposure to HRV-response assessment explains the difference in direction observed in our previous analyses in the same road maintenance workers (Meier et al., ) and also other healthy individuals (Riediker et al., ), where we observed an increase in mean HRV associated with exposure to PM2.5 in the morning after exposure. Taken together, the previous analyses (Riediker et al., ; Meier et al., ) and this one are in support of the idea that immediately after exposure, air pollutants disturb cardiac function, which will fade after a few hours and be followed later on by a recovery response.

Another novel, and possibly troubling new insight of our analysis is the observed funnel plot-like reduced modulation of HRV following exposure. It supports the idea of an interference with homeostatic control mechanisms. It seems to pose a serious health risk if such a strong response was observed in people with already partly derailed control mechanisms, such as those having pre-existing cardiovascular disease, previous infarction, heart failure, or arrhythmic disorders. The acute impairment of homeostatic control would also explain the later occurring increase of HRV observed in the same workers (Meier et al., ) and also in an independent populations of healthy highway patrol men (Riediker et al., ), which would represent the period of homeostatic recovery.

Availability of data and material

The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.

Statements

Author contributions

All authors contributed to the writing of the article. MR: led the analysis and led the writing of the article; YF: conducted the statistical analysis and created figures and tables; VR: reviewed and contributed to the statistical planning and interpretation; RM: collected the health and exposure data, prepared the database, and helped with the statistical analysis; MB: contributed to the statistical assessment and provided homeostasis interpretation; WEC treated the ECG recordings, extracted 5-min HRV data, and contributed to the cardiac interpretation of the findings.

Funding

Data analyzed in this study was collected as part of a study sponsored by the Swiss National Science Foundation (PDFMP3_124954). The here presented data analysis was sponsored by the Institute for Work and Health, University of Lausanne, Lausanne, Switzerland.

Acknowledgments

We thank Wayne E. Cascio of the National Health and Environmental Effects Research Laboratory, U.S. Environmental Protection Agency for his important support in the initial data collection. He also gave valuable input to the interpretation of the findings of the here presented analysis.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

Summary

Keywords

heart rate variability, fine particulate matter, acute exposure, homeostatic interference, ARIMA methodology

Citation

Riediker M, Franc Y, Bochud M, Meier R and Rousson V (2018) Exposure to Fine Particulate Matter Leads to Rapid Heart Rate Variability Changes. Front. Environ. Sci. 6:2. doi: 10.3389/fenvs.2018.00002

Received

25 September 2017

Accepted

08 January 2018

Published

25 January 2018

Volume

6 - 2018

Edited by

Maurice Millet, Université de Strasbourg, France

Reviewed by

Eduardo Landulfo, Instituto de Pesquisas Energéticas e Nucleares (CNEN), Brazil; Satish Vutukuru, North Carolina State University, United States

Updates

Copyright

*Correspondence: Michael Riediker

This article was submitted to Air Pollution, a section of the journal Frontiers in Environmental Science

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics