Pascale Haddad1*
Katherine Ogurtsova1
Sarah Lucht1,2Lina Glaubitz1Peter Höppe3
Dennis Nowak3Peter Angerer1Barbara Hoffmann1
- 1Institute for Occupational Social and Environmental Medicine, Centre for Health and Society, Medical Faculty and University Hospital Düsseldorf, Heinrich Heine University Düsseldorf, Düsseldorf, Germany
- 2Real-World Evidence & Insights, Cardinal Health, Dublin, OH, United States
- 3Institute and Outpatient Clinic for Occupational, Social and Environmental Medicine, Ludwig-Maximilians-University Munich, Munich, Germany
Background: Short-term exposure particulate matter with a diameter of 10 µm or less (PM10) and fine particulate matter (PM2.5) has been associated with heart rate variability (HRV), but exposure to ultrafine particles (UFP) has been less well examined. We investigated the associations between the HRV outcomes and short-term exposure to UFP, PM10 and PM2.5 among school-aged children and seniors.
Methods: CorPuScula (Coronary, Pulmonary and Sanguis) is a longitudinal, repeated-measure panel study conducted in 2000–2002 in Munich, Germany including 52 seniors (58–94 years old) with 899 observations and 50 children (6–10 years old) with 925 observations. A 10-min resting electrocardiogram was performed to assess resting HRV outcomes [Standard Deviation of Normal to Normal Intervals (SDNN), Root Mean Square of Successive Differences between Normal Heartbeats (RMSSD), Low Frequency power (LF), High Frequency power (HF), ration between low and high frequency (LF/HF)]. UFP and PM exposures were measured near the care home and school yard for seniors and children, respectively. Mean exposures during the day of examination (9–21 h) as well as 3-h, 12-h, 24-h, one-day, and two-day lags were assessed. Linear mixed-effect models were used to investigate the associations between short-term air pollution and HRV outcomes separately in children and seniors. The models were adjusted for sex, age, weather conditions (temperature, precipitation, and water vapor pressure), BMI, lifestyle and medical information. Two and multipollutant models adjusted for NO2 and O3 were performed.
Results: Among seniors, we observed increases in SDNN, LF, HF and LF/HF ratio after short-term exposure to UFP (hourly and daily lags) in contrast to decreases in SDNN and RMSSD after exposure to PM10. Associations were generally robust to two- and multipollutant adjustment. Among children, we observed increases of the LF/HF ratio after short-term exposures to UFP at lags 12 and 24 h. In contrast, we observed decreases of the ratio after exposure to PM2.5 and PM10. Results were largely unchanged for multipollutant modelling, however we found a more pronounced increase in SDNN and LF/HF (UFP lag 12 and 24 h) after adjusting for NO2.
Conclusions: Overall, among seniors, we observed associations of UFP and PM10 exposure with sympathetic responses of the ANS, which play an important role in sudden heart attacks or arrhythmia. Among children we found more inconsistent associations between UFP and a delayed increase in HRV. Adjusting for co-pollutants including NO2 and O3 yielded robust results.
1. Introduction
Heart rate (HR) is the number of heart beats per given time period, conventionally given per minute. Heart rate variability (HRV) is the fluctuation in the time intervals between successive heart beats (1). HRV is a non-invasive dynamic metric and assesses the activity of the autonomic nervous system (ANS), including both the sympathetic and parasympathetic nervous systems (SNS and PNS, respectively). Importantly, the ANS has a well-established role as an indicator of cardiovascular disease (CVD) risk among adults (2–5). The ANS controls involuntary actions like HR, body temperature, digestion, perspiration, and the widening or narrowing of blood vessels (1). In the past decades, several studies have attributed imbalance in the ANS to some CVD clinical conditions such as sudden death, coronary artery disease, or heart failure (6). The ANS regulates the HR and blood pressure in the short-term to cope with everyday situations (7). A higher HRV indicates better general health (8), whereas a reduced HRV has been reported among children with Attention Deficit/Hyperactivity Disorder (9) or overweight and obesity (10). Moreover, among adults, a reduced HRV is associated with the development of cardio-metabolic diseases such as diabetes (11).
The literature shows that ambient air pollution, and more specifically the exposure to fine particulate matter (PM), is a major threat to older people with a CVD history and causally contributes to CVD morbidity and mortality in general (12–17). However, fewer studies have assessed the effects of ultrafine particles (UFP) with an aerodynamic diameter of ≤100 nm on CVD. UFPs are of great interest since they have the ability—due to their small size—to reach the pulmonary alveoli, from where they can pass into the circulatory system and reach all organs (18). Specifically, the potential systemic translocation of UFP can lead to endothelial and vascular dysfunction, hypertension, thrombosis, atherosclerosis and eventually might promote vasoconstriction in the coronary arteries, increasing the risk of stroke, coronary heart disease, heart failure and ultimately mortality (18–23). Moreover, UFPs are thought to be more toxic than larger PM and to have independent adverse health effects since their high reactivity surface area enables them to adsorb a large amount of toxic metals and organic compounds. These compounds generate oxidative stress that leads to inflammation and eventually increases risk for cardiovascular and respiratory diseases (24, 25).
It has been hypothesized that the inhaled UFP, PM10 and PM2.5, can activate the pulmonary reflexes that modify the autonomic control of cardiovascular function and potentially mediate an immediate response (within a few hours) leading to an autonomic imbalance that can be assessed by HRV outcomes (26). The latest review and meta-analysis on short-term UFP exposure, including 12 studies published up to June 2022, found decreases in some HRV outcomes within hours and up to days of exposure (22). Specific studies on short-term UFP exposure (hours to days) have found consistent decreases in HRV among adults but very little information is available for children (21, 27–29). Moreover, recent studies showed consistent decreases in HRV outcomes with short-term exposures to fine PM (20, 21). To our knowledge, only one study reported the same inverse association between PM and HRV in children (30). Moreover, the use of multipollutant modelling in this area is also scare. Previous work has investigated the effects of PM2.5 and ozone (31) on HRV, and some studies have looked at the adverse effects of two-pollutant models (UFP and PM2.5) on HRV (32–34). However, it is also important to look at the effect of NO2 since UFP and NO2 share traffic as their main source and might be correlated. Without adjusting for co-pollutants, it is challenging to attribute any association with HRV to UFP alone.
We therefore investigated the association of short-term exposure to UFP, PM10 and PM2.5 with HRV among elementary school children and seniors in a care-home in single, two and multi-pollutant models taking into account other co-pollutants (i.e., NO2, O3 and PM2.5).
2. Materials and methods
2.1. Study design and participants
CorPuScula (Cor for coronary, Pu for pulmonary and S for sanguis or blood) is a longitudinal repeated-measures panel study conducted in 2000–2002 by the Institute and Outpatient Clinic for Occupational, Social and Environmental Medicine, Ludwig-Maximilians-University (LMU) in Munich, Germany. The study included 52 seniors and 50 children.
2.1.1. Seniors
The seniors were recruited via letter and an informational event at their care home in Munich (Wohnstift Augustinum Munich-North). Participation was voluntary and all those who consented to participate were invited to an initial one-hour long personal interview in the measurement van located outside their care home. Moreover, the participants filled in a standardized questionnaire and individual characteristics were collected including: age, height, weight, smoking habits, alcohol consumption, chronic illnesses, medication intake, and physician-diagnosed allergies. At the end of the initial interview, the participants gave written informed consent. Inclusion criteria for the seniors were as follows: (1) older than 55 years of age, (2) resident of the care home Wohnstift Augustinum Munich-North, and (3) time availability to complete the study over the upcoming year. The exclusion criteria included: (1) poor physical health that would prevent obtainment of clinical measurements, (2) current tobacco use, (3) presence of a pacemaker, and (4) intake of blood thinners (e.g., Marcumar). The participants were followed between May 2000 and July 2001 except during acute sickness. If possible, the appointments were scheduled at the same time of the day and the same day of the week to eliminate potential confounding from circadian rhythms or weekly effects. Routine examinations included the following: (1) a standardized questionnaire, (2) a bi-weekly 10-min electrocardiogram, (3) a bi-weekly blood sample, (4) blood-pressure measurement, and (5) lung-function measurements.
2.1.2. Children
In the same way, the children were recruited at the Strehleranger elementary school in Munich-Neuperlach by sending a letter describing the study along with a response form to the parents/guardians. Participation was voluntary, and all parents/guardians who gave their consent were invited alongside their children to an initial one-hour long personal interview and filled in the initial questionnaire. For all participants, at least one parent/guardian was required to be present at the first interview. The only inclusion criterion was attendance at the Strehleranger elementary school. The only exclusion criterion was exposure to tobacco at home or in private life. The children were followed between September 2001 and July 2002 except during vacation, school holidays, or acute sickness. Similarly to the seniors and following the same strategy, the appointments were rescheduled if possible. Routine examinations included the same measurements as in the senior cohort except for the withdrawal of a blood sample. In this paper and for both groups, we will focus on the HRV outcomes only. A flow diagram reflecting all the targeted number of observations vs. the ones available and eventually used in the analyses can be found in Supplementary Figures S1A,B. This study was approved by the Ethics Committee of the Medical Faculty at LMU (No. 130/00).
2.2. Exposure assessment: air pollutants and meteorological parameters
UFPs (≤100 nm in diameter) were measured as particle number concentration (PNC, particles/cm3). The measurements for the seniors were carried out with a condensation particle counter (TSI 3022A, lower particle size detection limit >7 nm) which was located approximately 3 km away from the care home. For children, a particle counter (TSI 3025A, lower particle size detection limit >3 nm) was placed on the school grounds, which is located on the Strehleranger road. The calendar daytime mean concentrations and individual hourly and lagged exposures prior to the examinations were calculated from hourly means.
Particulate matter PM10 and PM2.5 concentrations were measured using a Low Volume Sampler (LVS3; Leckel brand) located in the park of the Wohnstift Augustinum and the courtyard of the Strehleranger elementary school for both the seniors and the children's cohorts, respectively. Measures of fine PM were collected for 24 h but only the concentrations between 9:30 am until 9:30 pm were used as these were the main times in which the participants were active outdoors. The daytime mean concentrations were then calculated and used in the analyses as lags 0, 1 and 2.
Exposure to other air pollutants and meteorological parameters (i.e., temperature, water vapor pressure, precipitation and wind speed) were collected at the measuring station Johanneskirchen (L8.12). The station was located 9 km north and 7 km south of the senior home and the school, respectively. We obtained air pollution data from the station for nitrogen dioxide (NO2, μg/m3) and ozone (O3, μg/m3). Meteorological variables for the senior cohort were obtained from the measuring station of the Meteorological Institute of the LMU in Munich (MIM), which is approximately 8 km from the care home. For the children's cohort, we obtained the meteorological data from the MIM's measuring station Theresienstrasse (approximately 7 km away). The meteorological data included hourly mean values for air temperature (°C), water vapor pressure (hPa), precipitation (mm) and wind speed (m/s). We calculated the corresponding daytime mean values for the measured air pollution intervals (9:30 am to 9:30 pm).
2.3. Outcomes: HRV measures
Every other week, HRV was measured using a long-term ECG monitoring system from the company Customed (Custo Flash 200). This monitor used 2 channels alongside 5 electrodes. Raw ECG data was processed and cleaned prior to HRV analyses by an experienced doctor. Quality control was performed for the 10-min measurement interval, and extra systoles were identified and removed from the dataset.
For data collection, participants lay relaxed on the bed for at least 5 min before and during the entire 10 min of the ECG examinations. They were instructed not to speak or move and to breathe calmly and evenly during the examination. We extracted the time-domain measures standard deviation of normal-to-normal intervals (SDNN) representing both the SNS and PNS activity, and the root mean square of successive differences between normal heartbeats (RMSSD) that is more strongly influenced by the PNS than SDNN. The frequency-domain measures yielded information about the total power in the heart's rhythm explained by the variability of the normal-to-normal intervals. These include the low-frequency power (LF, 0.04–0.15 Hz) influenced by both SNS and the PNS, the high-frequency power (HF, 0.15–0.4 Hz) reflecting the PNS, and the frequency ratio of LF/HF to estimate the balance between SNS and PNS activity.
2.4. Statistical analysis
We analyzed associations between UFP, PM2.5/PM10 and HRV outcomes using linear mixed-effects regression models separately in adults and children. The time variable in the models was the date of measurement. The models also included the personal identification number as random effect. The autocorrelation structure included the time covariate and the personal identification number as a grouping factor within the observations that were assumed correlated. We used both unweighted (homoscedastic) as well as weighted (heteroscedastic) models where we allowed variances to differ by person, with a control for time trend (35). The weights were calculated as a ratio of the stratum standard errors to the within-group standard error. We assessed the adequacy of the heteroscedastic fit by re-examining plots of the standardized residuals vs. the fitted values by stratum to check if the standardized residuals had the same variability in each stratum. Also, the fit was tested with the ANOVA method. 95% confidence intervals (CIs) and p-values were computed using a Wald t-distribution approximation.
For each pollutant (UFP, PM10 and PM2.5), we calculated lagged exposures including both single-day lags (i.e., lag 0, lag 1 and lag 2) as well as individual lags for PM10 (3, 12 and 24 h) and UFP (1, 3, 12 and 24 h).
The 5 outcomes (SDNN, RSMMD, LF, HF and LF/HF) were analyzed separately. In regression analyses, we first estimated the crude association between the lag periods for single pollutants with each outcome. Second, we estimated the same associations adjusting for possible confounders as found in previous studies (21, 36). To identify potential confounding and minimal sufficient adjustment sets, we constructed a directed acyclic graph (DAG) using the DAGitty program (37) separately for seniors and children. For seniors, the second model included age (years), sex (male/female), alcohol [did they drink on the day of the examination (yes/no), repeated measure], BMI (kg/m2), day of the week, precipitation, temperature, and water vapor pressure (see Supplementary Figure S2A). For children, the second model included age (years), sex (male/female), BMI (kg/m2), day of the week, precipitation, temperature, and water vapor pressure (see Supplementary Figure S2B). We did not adjust for multiple comparisons because the tests were highly correlated, thus violating the assumption of independence (see Supplementary Table S1 for correlations of outcome variables and the results section for correlations of exposure variables).
Meteorological variables (i.e., temperature, precipitation, and water vapor) were assessed for linearity using models without exposure variables, which resulted in modelling temperature and water vapor pressure with natural splines using 4 knots and precipitation with no splines. The best fitted lags for meteorological variables, which varied by outcome and group, were used (see Supplementary Table S1).
We then conducted two and multipollutant models by simultaneously adding the best fitted concentration (lag 0–2) of NO2 and O3 to the main models to examine the robustness of the effects after further controlling for the exposure (see Supplementary Table S2). Also, for models where the main exposure was UFP, we further adjusted for PM2.5 and vice versa.
We also ran further sensitivity analyses to assess the robustness of our results. For the seniors, we investigated the inclusion of the intake of heart medication (yes/no if they were prescribed a heart meds at the time) before every examination. We also repeated analyses including season as an extra meteorological factor (4 seasons: winter, spring, summer, and autumn). For the children, we investigated the inclusion of having hay fever symptoms (yes/no) at the day of the examination as well as season.
For seniors and as examined in prior studies (38), we evaluated potential effect modification by inflammation markers (i.e., C-reactive protein as a continuous variable) and for blood pressure (hypertensive vs. not hypertensive). We conducted effect modification analyses through addition of a multiplicative interaction term between the exposure and the covariate of interest into the main model.
The models were estimated using restricted maximum likelihood (REML) and “nlminb” optimizer [R-package “state” (39)]. We employed the R functions varIdent() to define the variance function structure and lme() to fit the models [R-package “nlme” (35, 40)]. All statistical analyses were conducted in R version 4.3 (R Core Team) (39).
3. Results
In total, 47 seniors (37 females, mean age = 76.7) and 47 children (27 females, mean age = 8.1) were included in the analyses (Table 1). Six seniors were obese and a third were former smokers. The majority were on heart medication (61.7%) and half had hypertension. Among children, only one was obese, none reported asthma problems, and none were on regular medication.
Table 1. Characteristics of the study population at the initial examination.
Seniors were followed up between May 2000 and June 2001 with a total of 899 available observations on HRV and an average of 19.1 observations per senior. Of the 899 observations, 603 and 578 had complete UFP and fine PM measures, respectively, with no missing values (Supplementary Figure S1A). The mean values for SDNN and RMSSD were 30.7 ms (SD = 14.3) and 22.7 ms (SD = 12.2) respectively (Table 2). For the seniors, the mean UFP exposure was 27,471 (SD = 15,100) n/ml. Very high positive correlations were found between PM2.5 and PM10 with medium to high positive correlations between PM2.5, PM10, NO2 and UFP and medium to high negative correlations between UFP and O3 (Table 3).
Table 2. HRV outcome characteristics averaged over examination days.
Table 3. Air pollution characteristics averaged over examination days.
Children were followed up between September 2001 and July 2002 with a total of 925 available observations on HRV and an average of 19.7 observations per child. Of the 925 observations, only 325 had complete UFP measures but 710 had complete fine PM measures (Supplementary Figure S1B). The mean values for SDNN and RMSSD were 49.3 ms (SD = 16.3) and 46.3 ms (SD = 18.8) respectively (Table 2). The mean UFP exposure was 19,872 (SD = 5,482) n/ml. Similar to the seniors, a very high positive correlation was observed between PM2.5 and PM10. However very low negative correlations were present between PM2.5, PM10, NO2 and UFP (Table 3).
3.1. Seniors
In UFP single pollutant models among seniors (Figure 1 and Supplementary Table S3), we observed an increase in SDNN at lags 0 and 1 day. In contrast, at short term exposures to PM10, we observed decreases in SDNN (lags 12 and 24 h, and lag 1 day). No associations were found between UFP exposure and RMSSD. PM10 (12 h lag) was associated with a clear decrease in RMSSD. In relation to the frequency-domain outcomes, UFP exposure was associated with consistent increases in both LF and HF powers at all lags, whereas for PM10, only very acute exposure was associated with a transient increase in LF and HF. UFP exposures (lags 12 and 24 h and lags 1 and 2 days) were positively associated with the LF/HF ratio, but not PM10 was not. No associations were observed between exposure to PM2.5 and any outcomes.
Figure 1. Changes in outcomes per IQR exposure increase to UFP, PM10 and PM2.5 among seniors. The circular estimates represent the hours (1, 3, 12 or 24 hours) after exposure and the triangular estimates represent the daily lags (0, 1, 2).
In UFP models, we observed robust results after adjusting for NO2 and O3 throughout the 5 different outcomes. Results were also robust to adjustment for PM2.5 with the exception of a more pronounced increase in SDNN and LF at lag 0 and lag 1, thus also affecting the estimates at lag 0 and 1 in multipollutant modelling (Figure 2 and Supplementary Table S4). In PM10 models, results were robust to adjustment for NO2 and O3 in two and multipollutant modelling (Supplementary Figure S3).
Figure 2. Changes in outcomes per IQR exposure increase to UFP, PM10 and PM2.5 among seniors - two and multi-pollutant models. The circular estimates represent the hours (1, 3, 12 or 24 hours) after exposure and the triangular estimates represent the daily lags (0, 1, 2).
3.2. Children
Among children no clear associations were found between UFP and fine PM and the time-domain outcomes (Figure 3 and Supplementary Table S5). While point estimates for UFP were consistently below 0 for LF, all CIs were large and overlapped the null. In contrast, no patterns could be observed for UFP and HF. Similarly for PM10, consistent but weak decreases were seen for LF but no pattern for HF. Short-term exposures to UFP (12 and 24 h) were positively associated with the LF/HF ratio. PM2.5 and PM10 were both negatively associated with the ratio at lag 2.
Figure 3. Changes in outcomes per IQR exposure increase to UFP, PM10 and PM2.5 among children. The circular estimates represent the hours (1, 3, 12 or 24 hours) after exposure and the triangular estimates represent the daily lags (0, 1, 2).
In UFP models, results were robust to adjustment for PM2.5 and O3 throughout the 5 different outcomes. Models were also robust to adjustment for NO2 with the exception of a more pronounced increase in SDNN and LF/HF ratio at lag 12 and 24 h thus also affecting the estimates in multipollutant modelling (Figure 4 and Supplementary Table S6). In PM10 models, results were robust to adjustment for NO2 and O3 in two and multipollutant modelling (Supplementary Figure S4).
Figure 4. Changes in outcomes per IQR exposure increase to UFP, PM10 and PM2.5 among children - two and multi-pollutant models. The circular estimates represent the hours (1, 3, 12 or 24 hours) after exposure and the triangular estimates represent the daily lags (0, 1, 2).
3.3. Sensitivity and effect modification analyses
Sensitivity analyses showed that UFP and PM10 results were robust to adjustment for heart medication and season among seniors (Supplementary Figures S5, S6). Among children, model results for UFP and PM10 were largely unchanged after adjusting for season and hay fever (Supplementary Figures S7, S8). There were no clear differences in the effects when it came to analyses by hypertension (Supplementary Figures S9, S10) or CRP (results not shown) among seniors in any of the exposures.
4. Discussion
4.1. Our findings
Among seniors, we observed positive associations between SDNN, LF, HF and LF/HF ratio and short-term exposure to UFP (hourly and daily lags) in contrast to decreases in SDNN and RMSSD after exposure to PM10. Generally, the association between fine PM and all 5 HRV outcomes was less clear. For two and multipollutant modelling, we generally found UFP associations were robust to adjustment. Among children, we observed inconsistent associations between UFP and all HRV outcomes and mostly null associations between PM2.5 and PM10 and HRV outcomes. Results were robust for multipollutant modelling in general, however, adjustment for NO2 enhanced some patterns for the association of UFP with SDNN and LF/HF (UFP lag 12 and 24 h).
4.2. Comparing with the literature
In terms of short-term exposure to UFP, and in contrast to our findings, a recent meta-analysis by Zhang et al. (22) showed decreases in SDNN, RMSSD, LF and LF/HF within 6 h after exposure. Moreover, they found that daily lags on the same day or preceding days were not associated with HRV which is in contrast with our results of a positive association with SDNN, LF, HF and LF/HF. In terms of short-term exposure to fine PM, previous studies on the general population show inconsistent associations with HRV outcomes. Some observed decreases in time-domains (21, 41) and frequency-domain outcomes (20, 27, 42, 43), but also others observed increases in both domains (43–45) or no associations at all (46–48) as seen in our study.
Further, a very limited number of studies have looked at the exposure of ambient short-term UFP and fine PM on HRV outcomes among children. A study in California demonstrated that adolescents residing in neighborhoods characterized by higher concentrations of PM2.5 had a decrease in HRV specifically in HF band in response to social stress (49). Saenen et al. (30), examined the association between fine PM exposure at home and school on SDNN, RMSSD, HF and LF among school children aged between 9 and 12 years old and found a decrease in LF with an increase in outdoor PM10 exposure as well as decreases in RMSSD, HF and LF with an increase in indoor PM10 exposure. In contrast, we observed no consistent and clear associations between fine PM exposure and HRV. Chen et al. (29) investigated the short-term effects of indoor UFP on children aged between 11 and 14 years old in schools in Beijing and found decreases in SDNN at lag 1 and decreases in SDNN, LF and HF at lag 2. Similarly, we found decreases in LF at most hourly and daily lags.
4.3. Multipollutant modelling and sensitivity analyses
Our results for two and multipollutant models showed mostly robust results. A review on ambient air pollution and HRV conducted by Buteau and Goldberg (50) observed some decreases in HRV after the exposure to NO2 ranging between 5 min to 5 days prior measurements. In our study and among seniors, the correlations between pollutants were high. Even though UFP was correlated with the co-pollutants, including NO2, the robustness in our results indicates that there was no major collinearity influencing the reliability of the two and multi-pollutant models and that the associations seen between UFP and HRV were independent of fine PM, NO2 and O3. Moreover, among children, we observed a weaker correlation between UFP and the co-pollutants compared to the seniors. We adjusted our analyses for NO2 as surrogate for traffic related air pollution (TRAP), thus reducing confounding of the air pollutant and HRV associations by sources of traffic. However, the potential for residual confounding by traffic remains due to the imperfect correlation between NO2 and TRAP.. In terms of adjusting for O3 and as observed by a recent systematic review, it is important to consider adjustment for O3 since an association has been observed between short-term exposure to O3 and a decrease in HRV (51). Moreover, in a study conducted on HRV in a two-pollutant model including PM2.5 and O3, it was concluded that PM2.5 influenced HRV rather than O3 (31). We obtained comparable conclusions where, after the adjustment of O3 in two and multipollutant modelling, the results were robust leading to the affirmation of the effects of UFP and fine PM.
It is important to note that half of the seniors were on heart medication and the previous literature observed a decrease in HRV activity among people with heart disease (31, 52). This might explain the low SDNN in comparison to other studies in the general population. Upon adjusting for heart medication, the estimates were robust.
4.4. Risk of exposure misclassification
It is important to note that the placement of the pollutant measurement stations may lead to pollutant-specific exposure misclassification. In particular, among seniors, the fine PM measurement sampler was located in the park of the care home whereas for UFP, the counter was placed 3 km. away from the care home. For NO2 and O3, the measuring station was located 9 km north of the care home. We therefore expect a higher degree of exposure misclassification for UFP and NO2. However, simultaneous measurements of UFP at the care home showed that the temporal concentration changes at both locations were very similar. Among children, the exposure measurement does not account for exposures of children at home but only at school, potentially contributing to the inconsistent results in children.
4.5. Biological mechanisms
Previous studies have demonstrated that exposure to UFP and fine PM causes both inflammation and oxidative stress that can contribute to cardiovascular disease and eventually cardiovascular mortality (53). A number of experiments have shown the biological mechanisms behind the association of UFP and fine PM with cardiovascular diseases at the molecular level. The very first contact of the particles is with the nasal and pulmonary receptors in the respiratory epithelium that will lead to a disturbance of the autonomic dysfunction and perturbations in the heart rate. Secondly, oxidative stress initiates an inflammatory response by the release of biologically active intermediates. This localized imbalance might also lead to an oxidative stress response and eventually lung inflammation. HRV is an indicator of the ANS that depends on the SNS and PNS for transferring information and where both systems interplay to regulate HRV (53). An imbalance in the function of the SNS and PNS could have opposite effects on HRV leading either to an increase or a decrease in HRV. More specifically, the effects of an exposure to UFP and fine PM could be mediated either through the SNS leading to decreases in HRV or through the PNS leading to increases in HRV (33). The SDNN and LF band reflect both SNS and PNS activities while the RMSSD and the HF band reflect PNS activity. We therefore consider the LF/HF as a ratio between SNS/PNS.
Among seniors we observed a general decrease in SDNN, RMSSD and HF after the exposure to PM10 which would suggest a sympathetic response. In line with this, we also observed increases in the LF/HF ratio upon the exposure to UFP. This suggests greater SNS than PNS responses, which coincides with previous studies where greater SNS responses played an important role in sudden heart attacks or arrhythmia (22). The same response was observed among children with increase in the LF/HF ratio upon the exposure to UFP leading to the risks of arrythmia and heart disease (22).
4.6. Strengths and limitations
This study has several strengths, firstly the ability to consider a range of exposure windows for all air pollutants and thus the assessment of temporal exposure patterns, both hourly and daily lags. Second the availability of data on co-pollutants allowed us to adjust for those in two and multipollutant modelling. Thirdly, the design of Corpuscula allowed for repeated-measures. While the number of participants was small, we were able to analyze a large number of observations per participant, which was enough to detect some associations among the seniors.
One main limitation is the small number of observations we had among children in terms of UFP measures and when adjusting for multi-pollutants. Because of the large number of missing values, the weighted models did not converge, and we had to report unweighted results, thus the confidence intervals were very wide and we were not able to detect clear associations. Another limitation was that we did not have personal exposure measures for NO2 and O3, since these co-pollutants were not measured at the school or the care-home but instead 7–8 km away from the actual research area. Thus, our ability to detect accurate associations for these pollutants in multipollutant modelling might have been reduced. A further limitation is that the study was carried among seniors and children, and the results may not be transferable to other age groups or to the general population. Nevertheless, these are populations traditionally thought to be particularly vulnerable to the effects of air pollution.
5. Conclusions
Overall, among seniors, we observed associations of UFP and PM10 exposure with sympathetic responses of the ANS, which play an important role in sudden heart attacks or arrhythmia. Among children we found more inconsistent associations between UFP and a delayed increase in HRV. Adjusting for co-pollutants including NO2 and O3 yielded robust results.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving humans were approved by the Ethics Committee of the Medical Faculty at LMU (No. 130/00). The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants’ legal guardians/next of kin.
Author contributions
PH: Methodology, Writing – original draft, Formal Analysis. KO: Data curation, Formal Analysis, Writing – review & editing. SL: Data curation, Methodology, Writing – review & editing. LG: Data curation, Methodology, Writing – review & editing. PH: Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Writing – review & editing. DN: Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Writing – review & editing. PA: Project administration, Supervision, Writing – review & editing. BH: Formal Analysis, Methodology, Supervision, Writing – review & editing.
Funding
The authors declare financial support was received for the research, authorship, and/or publication of this article.
This study was funded by the “Bayerisches Staatsministerium für Umweltschutz und Konsumentenschutz” (Bavarian State Ministry for the Environment and Consumer Protection).
Acknowledgments
We thank the Bavarian State Office for Environmental Protection, and the Meteorological Institute of the of the Ludwig-Maximilians-University Munich who provided us with data on gaseous pollutants, and meteorological data, respectively. We would specifically thank Wolfgang Guth and Rudolf Schierl as main investigators in Corpuscula. In addition we also thank the director of the Augustinum residential home at the time, Benne, and the principal of the elementary school at Strehleranger, Fränkl, for their great support in recruiting the subjects and in the technical conduct of our investigations. Finally we would like to thank Haeran Jeong for her work on the referencing and organization of the manuscript.
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.
The authors DN, KO declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Publisher's note
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fepid.2023.1278506/full#supplementary-material
References
1. Shaffer F, Ginsberg JP. An overview of heart rate variability metrics and norms. Front Public Health. (2017) 5:258. doi: 10.3389/fpubh.2017.00258
2. Cygankiewicz I, Zareba W. Heart rate variability. In: Buijs RM, Swaab DF, editors. Autonomic nervous system. Science Direct (2013). p. 379–93.
3. Xhyheri B, Manfrini O, Mazzolini M, Pizzi C, Bugiardini R. Heart rate variability today. Prog Cardiovasc Dis. (2012) 55:321–31. doi: 10.1016/j.pcad.2012.09.001
4. Hägglund H, Uusitalo A, Peltonen JE, Koponen AS, Aho J, Tiinanen S, et al. Cardiovascular autonomic nervous system function and aerobic capacity in type 1 diabetes. Front Physiol. (2012) 3:356. doi: 10.3389/fphys.2012.00356
5. Wang Y, Zhou L, Jiang H, Yu L. Editorial: autonomic nervous system and cardiovascular diseases: from brain to heart. Front Physiol. (2022) 13:884832. doi: 10.3389/fphys.2022.884832
6. Dayer M, MacIver DH, Rosen SD. The central nervous system and heart failure. Future Cardiol. (2021) 17:363–81. doi: 10.2217/fca-2020-0059
7. Armstrong R, Wheen P, Brandon L, Maree A, Kenny RA. Heart rate: control mechanisms, pathophysiology and assessment of the neurocardiac system in health and disease. QJM. (2022) 115:806–12. doi: 10.1093/qjmed/hcab016
8. Jarczok MN, Kleber ME, Koenig J, Loerbroks A, Herr RM, Hoffmann K, et al. Investigating the associations of self-rated health: heart rate variability is more strongly associated than inflammatory and other frequently used biomarkers in a cross sectional occupational sample. PLoS One. (2015) 10:e0117196. doi: 10.1371/journal.pone.0117196
9. Rukmani MR, Seshadri SP, Thennarasu K, Raju TR, Sathyaprabha TN. Heart rate variability in children with attention-deficit/hyperactivity disorder: a pilot study. Ann Neurosci. (2016) 23:81–8. doi: 10.1159/000443574
10. Birch SL, Duncan MJ, Franklin C. Overweight and reduced heart rate variability in British children: an exploratory study. Prev Med. (2012) 55:430–2. doi: 10.1016/j.ypmed.2012.09.015
11. Young HA, Benton D. Heart-rate variability: a biomarker to study the influence of nutrition on physiological and psychological health? Behav Pharmacol. (2018) 29:140–51. doi: 10.1097/FBP.0000000000000383
12. Kaufman JD, Elkind MSV, Bhatnagar A, Koehler K, Balmes JR, Sidney S, et al. Guidance to reduce the cardiovascular burden of ambient air pollutants: a policy statement from the American heart association. Circulation. (2020) 142:e432–47. doi: 10.1161/CIR.0000000000000930
13. Franklin BA, Brook R, Pope CA. Air pollution and cardiovascular disease. Curr Probl Cardiol. (2015) 40:207–38. doi: 10.1016/j.cpcardiol.2015.01.003
14. Newman JD, Bhatt DL, Rajagopalan S, Balmes JR, Brauer M, Breysse PN, et al. Cardiopulmonary impact of particulate air pollution in high-risk populations: JACC state-of-the-art review. J Am Coll Cardiol. (2020) 76:2878–94. doi: 10.1016/j.jacc.2020.10.020
15. Wang F, Liang Q, Sun M, Ma Y, Lin L, Li T, et al. The relationship between exposure to PM2.5 and heart rate variability in older adults: a systematic review and meta-analysis. Chemosphere. (2020) 261:127635. doi: 10.1016/j.chemosphere.2020.127635
16. Strak M, Weinmayr G, Rodopoulou S, Chen J, de Hoogh K, Andersen ZJ, et al. Long term exposure to low level air pollution and mortality in eight European cohorts within the ELAPSE project: pooled analysis. Br Med J. (2021) 374:n1904. doi: 10.1136/bmj.n1904
17. Wolf K, Hoffmann B, Andersen ZJ, Atkinson RW, Bauwelinck M, Bellander T, et al. Long-term exposure to low-level ambient air pollution and incidence of stroke and coronary heart disease: a pooled analysis of six European cohorts within the ELAPSE project. Lancet Planet Health. (2021) 5:e620–32. doi: 10.1016/S2542-5196(21)00195-9
18. La Rovere MT, Pinna GD, Maestri R, Mortara A, Capomolla S, Febo O, et al. Short-term heart rate variability strongly predicts sudden cardiac death in chronic heart failure patients. Circulation. (2003) 107:565–70. doi: 10.1161/01.cir.0000047275.25795.17
19. Kinlay S, Ganz P. Role of endothelial dysfunction in coronary artery disease and implications for therapy. Am J Cardiol. (1997) 80:11I–6I. doi: 10.1016/s0002-9149(97)00793-5
20. Brook RD, Rajagopalan S, Pope CA III, Brook JR, Bhatnagar A, Diez-Roux AV, et al. Particulate matter air pollution and cardiovascular disease: an update to the scientific statement from the American heart association. Circulation. (2010) 121:2331–78. doi: 10.1161/CIR.0b013e3181dbece1
21. Breitner S, Peters A, Zareba W, Hampel R, Oakes D, Wiltshire J, et al. Ambient and controlled exposures to particulate air pollution and acute changes in heart rate variability and repolarization. Sci Rep. (2019) 9:1946. doi: 10.1038/s41598-019-38531-9
22. Zhang S, Breitner S, Pickford R, Lanki T, Okokon E, Morawska L, et al. Short-term effects of ultrafine particles on heart rate variability: a systematic review and meta-analysis. Environ Pollut. (2022) 314:120245. doi: 10.1016/j.envpol.2022.120245
23. Stone PH, Godleski JJ. First steps toward understanding the pathophysiologic link between air pollution and cardiac mortality. Am Heart J. (1999) 138:804–7. doi: 10.1016/S0002-8703(99)70002-5
24. Kwon H-S, Ryu MH, Carlsten C. Ultrafine particles: unique physicochemical properties relevant to health and disease. Exp Mol Med. (2020) 52:318–28. doi: 10.1038/s12276-020-0405-1
25. Moreno-Ríos AL, Tejeda-Benítez LP, Bustillo-Lecompte CF. Sources, characteristics, toxicity, and control of ultrafine particles: an overview. Geosci Front. (2022) 13:101147. doi: 10.1016/j.gsf.2021.101147
26. Andersen MHG, Frederiksen M, Saber AT, Wils RS, Fonseca AS, Koponen IK, et al. Health effects of exposure to diesel exhaust in diesel-powered trains. Part Fibre Toxicol. (2019) 16:21. doi: 10.1186/s12989-019-0306-4
27. Biel R, Danieli C, Shekarrizfard M, Minet L, Abrahamowicz M, Baumgartner J, et al. Acute cardiovascular health effects in a panel study of personal exposure to traffic-related air pollutants and noise in Toronto, Canada. Sci Rep. (2020) 10:16703. doi: 10.1038/s41598-020-73412-6
28. Hu J, Fan H, Li Y, Li H, Tang M, Wen J, et al. Fine particulate matter constituents and heart rate variability: a panel study in Shanghai, China. Sci Total Environ. (2020) 747:141199. doi: 10.1016/j.scitotenv.2020.141199
29. Chen C, Liu S, Dong W, Song Y, Chu M, Xu J, et al. Increasing cardiopulmonary effects of ultrafine particles at relatively low fine particle concentrations. Sci Total Environ. (2021) 751:141726. doi: 10.1016/j.scitotenv.2020.141726
30. Saenen ND, Provost EB, Cuypers A, Kicinski M, Pieters N, Plusquin M, et al. Child’s buccal cell mitochondrial DNA content modifies the association between heart rate variability and recent air pollution exposure at school. Environ Int. (2019) 123:39–49. doi: 10.1016/j.envint.2018.11.028
31. Mirowsky JE, Carraway MS, Dhingra R, Tong H, Neas L, Diaz-Sanchez D, et al. Exposures to low-levels of fine particulate matter are associated with acute changes in heart rate variability, cardiac repolarization, and circulating blood lipids in coronary artery disease patients. Environ Res. (2022) 214:113768. doi: 10.1016/j.envres.2022.113768
32. Peters A, Hampel R, Cyrys J, Breitner S, Geruschkat U, Kraus U, et al. Elevated particle number concentrations induce immediate changes in heart rate variability: a panel study in individuals with impaired glucose metabolism or diabetes. Part Fibre Toxicol. (2015) 12:7. doi: 10.1186/s12989-015-0083-7
33. Timonen KL, Vanninen E, de Hartog J, Ibald-Mulli A, Brunekreef B, Gold DR, et al. Effects of ultrafine and fine particulate and gaseous air pollution on cardiac autonomic control in subjects with coronary artery disease: the ULTRA study. J Expo Sci Environ Epidemiol. (2006) 16:332–41. doi: 10.1038/sj.jea.7500460
34. Weichenthal S, Hatzopoulou M, Goldberg MS. Exposure to traffic-related air pollution during physical activity and acute changes in blood pressure, autonomic and micro-vascular function in women: a cross-over study. Part Fibre Toxicol. (2014) 11:70. doi: 10.1186/s12989-014-0070-4
35. Pinheiro JC, Bates DM. Mixed-effects models in S and S-PLUS. New York, Berlin, Heidelberg: Springer (2000).
36. Huang C, Tang M, Li H, Wen J, Wang C, Gao Y, et al. Particulate matter air pollution and reduced heart rate variability: how the associations vary by particle size in Shanghai, China. Ecotoxicol Environ Saf. (2021) 208:111726. doi: 10.1016/j.ecoenv.2020.111726
37. Textor J, Hardt J, Knüppel S. DAGitty: a graphical tool for analyzing causal diagrams. Epidemiology. (2011) 22:745. doi: 10.1097/EDE.0b013e318225c2be
38. 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–9. doi: 10.1289/ehp.7447
41. Graff DW, Cascio WE, Rappold A, Hou H, Huang YC, Devlin RB. Exposure to concentrated coarse air pollution particles causes mild cardiopulmonary effects in healthy young adults. Environ Health Perspect. (2009) 117:1089–94. doi: 10.1289/ehp0900558
42. Vora R, Zareba W, Utell MJ, Pietropaoli AP, Chalupa D, Little EL, et al. Inhalation of ultrafine carbon particles alters heart rate and heart rate variability in people with type 2 diabetes. Part Fibre Toxicol. (2014) 11:31. doi: 10.1186/s12989-014-0031-y
43. Devlin RB, Smith CB, Schmitt MT, Rappold AG, Hinderliter A, Graff D, et al. Controlled exposure of humans with metabolic syndrome to concentrated ultrafine ambient particulate matter causes cardiovascular effects. Toxicol Sci. (2014) 140:61–72. doi: 10.1093/toxsci/kfu063
44. Schneider A, Hampel R, Ibald-Mulli A, Zareba W, Schmidt G, Schneider R, et al. Changes in deceleration capacity of heart rate and heart rate variability induced by ambient air pollution in individuals with coronary artery disease. Part Fibre Toxicol. (2010) 7:29. doi: 10.1186/1743-8977-7-29
45. Samet JM, Rappold A, Graff D, Cascio WE, Berntsen JH, Huang YC, et al. Concentrated ambient ultrafine particle exposure induces cardiac changes in young healthy volunteers. Am J Respir Crit Care Med. (2009) 179:1034–42. doi: 10.1164/rccm.200807-1043OC
46. Heusser K, Tank J, Holz O, May M, Brinkmann J, Engeli S, et al. Ultrafine particles and ozone perturb norepinephrine clearance rather than centrally generated sympathetic activity in humans. Sci Rep. (2019) 9:3641. doi: 10.1038/s41598-019-40343-w
47. Tobaldini E, Bollati V, Prado M, Fiorelli EM, Pecis M, Bissolotti G, et al. Acute particulate matter affects cardiovascular autonomic modulation and IFN-γ methylation in healthy volunteers. Environ Res. (2018) 161:97–103. doi: 10.1016/j.envres.2017.10.036
48. Byrd JB, Morishita M, Bard RL, Das R, Wang L, Sun Z, et al. Acute increase in blood pressure during inhalation of coarse particulate matter air pollution from an urban location. J Am Soc Hypertens. (2016) 10:133–9.e4. doi: 10.1016/j.jash.2015.11.015
49. Miller JG, Gillette JS, Manczak EM, Kircanski K, Gotlib IH. Fine particle air pollution and physiological reactivity to social stress in adolescence: the moderating role of anxiety and depression. Psychosom Med. (2019) 81:641–8. doi: 10.1097/PSY.0000000000000714
50. Buteau S, Goldberg MS. A structured review of panel studies used to investigate associations between ambient air pollution and heart rate variability. Environ Res. (2016) 148:207–47. doi: 10.1016/j.envres.2016.03.013
51. Zong Z, Zhang M, Xu K, Zhang Y, Hu C. Association between short-term exposure to ozone and heart rate variability: a systematic review and meta-analysis. Int J Environ Res Public Health. (2022) 19:11186. doi: 10.3390/ijerph191811186
52. Cowan MJ, Pike K, Burr RL. Effects of gender and age on heart rate variability in healthy individuals and in persons after sudden cardiac arrest. J Electrocardiol. (1994) 27(Suppl):1–9. doi: 10.1016/s0022-0736(94)80037-5
Keywords: ultrafine particles, heart rate variability, co-pollutants, children, seniors, sympathetic and para-sympathetic nervous systems
Citation: Haddad P, Ogurtsova K, Lucht S, Glaubitz L, Höppe P, Nowak D, Angerer P and Hoffmann B (2023) Short-term exposure to ultrafine and fine particulate matter with multipollutant modelling on heart rate variability among seniors and children from the CorPuScula (coronary, pulmonary, sanguis) longitudinal study in Germany. Front. Epidemiol. 3:1278506. doi: 10.3389/fepid.2023.1278506
Received: 16 August 2023; Accepted: 19 October 2023;
Published: 8 November 2023.
Edited by:
Kai Chen, Yale University, United StatesReviewed by:
Jason Wong, National Cancer Institute (NIH), United StatesSiqi Zhang, Helmholtz Association of German Research Centres (HZ), Germany
© 2023 Haddad, Ogurtsova, Lucht, Glaubitz, Höppe, Nowak, Angerer and Hoffmann. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Pascale Haddad pascale.haddad@med.uni-duesseldorf.de
Abbreviations ANS, autonomic nervous system; CO, carbon monoxide; CVD, cardiovascular disease; HF, high frequency power; HRV, heart rate variability; LF, low frequency power; MIM, Meteorological Institute of the LMU in Munich; NO2, nitrogen dioxide; O3, ozone; PM, particulate matter; PNS, parasympathetic nervous system; RMSSD, root mean square of successive differences between normal heartbeats; SDNN, standard deviation of normal to normal intervals; SNS, sympathetic nervous system; UFP, ultrafine particles.