Association of Plasma IL-6 and Hsp70 with HRV at Different Levels of PAHs Metabolites

Abstract

Background

Exposure to polycyclic aromatic hydrocarbons (PAHs) is associated with reduced heart rate variability (HRV), a strong predictor of cardiovascular diseases, but the mechanism is not well understood.

Objectives

We hypothesized that PAHs might induce systemic inflammation and stress response, contributing to altered cardiac autonomic function.

Methods

HRV indices were measured using a 3-channel digital Holter monitor in 800 coke oven workers. Plasma levels of interleukin-6 (IL-6) and heat shock protein 70 (Hsp70) were determined using ELISA. Twelve urinary PAHs metabolites (OH-PAHs) were measured by gas chromatography-mass spectrometry.

Results

We found that significant dose-dependent relationships between four urinary OH-PAHs and IL-6 (all Ptrend<0.05); and an increase in quartiles of IL-6 was significantly associated with a decrease in total power (TP) and low frequency (LF) (Ptrend = 0.014 and 0.006, respectively). In particular, elevated IL-6 was associated in a dose-dependent manner with decreased TP and LF in the high-PAHs metabolites groups (all Ptrend<0.05), but not in the low-PAHs metabolites groups. No significant association between Hsp70 and HRV in total population was found after multivariate adjustment. However, increased Hsp70 was significantly associated with elevated standard deviation of NN intervals (SDNN), TP and LF in the low-PAHs metabolites groups (all Ptrend<0.05). We also observed that both IL-6 and Hsp70 significantly interacted with multiple PAHs metabolites in relation to HRV.

Conclusions

In coke oven workers, increased IL-6 was associated with a dose-response decreased HRV in the high-PAHs metabolites groups, whereas increase of Hsp70 can result in significant dose-related increase in HRV in the low-PAHs metabolites groups.

About the Authors:

Jian Ye

Contributed equally to this work with: Jian Ye, Rui Zhu

Affiliation: Department of Occupational and Environmental Health and Ministry of Education Key Lab for Environment and Health, School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China

Rui Zhu

Contributed equally to this work with: Jian Ye, Rui Zhu

Affiliation: Department of Traditional Chinese Medicine, Union Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China

Xiaosheng He

Affiliation: Department of Occupational and Environmental Health and Ministry of Education Key Lab for Environment and Health, School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China

Yingying Feng

Affiliation: Department of Occupational and Environmental Health and Ministry of Education Key Lab for Environment and Health, School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China

Liangle Yang

Affiliation: Department of Occupational and Environmental Health and Ministry of Education Key Lab for Environment and Health, School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China

Xiaoyan Zhu

Affiliation: Department of Occupational and Environmental Health and Ministry of Education Key Lab for Environment and Health, School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China

Qifei Deng

Affiliation: Department of Occupational and Environmental Health and Ministry of Education Key Lab for Environment and Health, School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China

Tangchun Wu

Affiliation: Department of Occupational and Environmental Health and Ministry of Education Key Lab for Environment and Health, School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China

Xiaomin Zhang

* E-mail: [email protected]

Affiliation: Department of Occupational and Environmental Health and Ministry of Education Key Lab for Environment and Health, School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China

Introduction

Polycyclic aromatic hydrocarbons (PAHs) are the major component of coke oven emissions produced during incomplete combustion of natural or synthetic fuels [1], [2]. Recent studies, including those from our lab [3], [4], found that increased exposure to PAHs was associated with reduced cardiac autonomic function, assessed by heart rate variability (HRV), which is considered one of the main pathophysiologic pathways for air pollution-mediated adverse cardiac events [5]–[8]. However, the mechanisms through which the cardiac autonomic system responds to PAHs exposure have not been sufficiently understood.

Systemic inflammation may be one of the potential mechanisms linking air pollution including PAHs to cardiovascular autonomic decline. It has been reported that elevated levels of inflammation biomarkers, especially that of C-reactive protein (CRP) and interleukin-6 (IL-6), are inversely associated with various HRV indices in apparently normal adults [9]–[11] and in patients with coronary heart disease (CHD) [12], [13]. However, to our knowledge, no studies have evaluated the relationship between inflammation and HRV in occupational population exposed to different levels of PAHs.

Heat shock protein 70 (Hsp70), one of the main HSP family members, has a dual role as chaperone and cytokine by inducing pro-inflammatory cytokines such as IL-6 production [14]. A number of studies have shown that circulating Hsp70 play a crucial role in the pathogenesis and/or prognosis of cardiovascular diseases [15]–[18]. We previously found that exposure to PAHs, as reflected by urinary 1-hydroxypyrene (1-OHP), resulted in a dose-dependent increase in levels of plasma Hsp70, suggesting that high Hsp70 levels may serve as a danger marker among workers exposure to PAHs [19]. We thus hypothesized that circulating Hsp70 may be associated with cardiac autonomic function in response to occupational PAHs exposure.

To test the above hypothesis, we examined the association of urinary monohydroxylated PAHs (OH-PAHs) with plasma IL-6 and Hsp70, and the effects of IL-6 and Hsp70 on HRV in 800 workers exposed to different levels of PAHs.

Materials and Methods

Ethics Statement

The study was approved by the Ethics and Human Subject Committee of Tongji Medical College and informed written consent was obtained from each subject.

Study Subjects

The cross-sectional survey was conducted throughout the workday (Monday to Friday) from Oct to Nov 2009 in Wuhan (Hubei, China). Eight hundred healthy workers who had worked for at least one year at different sites, including at the top, side, bottom, or adjacent workplaces of coke ovens, and at the offices, in the same coke oven plant, were included in this study. Subjects who had suffered from cardiovascular diseases, serious medical conditions affecting HRV or chronic inflammation in the previous 3 months were excluded from the study. In the morning, workers visited the occupational health examination center. 5 mL of venous blood and 20 mL of urine samples were collected at the beginning (pre-shift) of workshift from the workers who were on the day shift, and the end (post-shift) of workshift from those who were on the night shift. All blood and urine samples were divided into aliquots and frozen at −80°C until laboratory analysis. Questionnaires, physical examination and Holter monitoring followed to be completed. Information on demographic variables, occupational history, medications, history of disease, and lifestyle habits including smoking, passive smoking, alcohol consumption, physical activity, and diet were collected using standardized occupational questionnaires.

Determination of Urinary PAHs Metabolites

We measured 12 urinary PAHs metabolites [pyrene metabolite: 1-hydroxypyrene (1-OHP); naphthalene metabolites: 1-hydroxynaphthalene (1-OHNa), 2-OHNa; fluorene metabolites: 2-hydroxyfluorene (2-OHFlu), 9-OHFlu; phenanthrene metabolites: 1-hydroxyphenanthrene (1-OHPh), 2-OHPh, 3-OHPh, 4-OHPh, 9-OHPh; chrysene metabolite: 6-hydroxychrysene (6-OHChr); and benzo[a]pyrene metabolite: 3-hydroxybenzo[a]pyrene (3-OHBaP)] by gas chromatography-mass spectrometry (GC/MS, Agilent 6890N+5975B, Agilent Technologies Inc., Santa Clara, CA, USA) as previously described [4]. Briefly, each 3.0 mL urine sample was extracted three times to elevate the detection rate. The set of the standard curve was operated about every 100 samples and re-run about 10% of the total samples were used as quality control. The identification and quantification of urinary PAHs metabolites were based on retention time, mass-to-charge ratio, and peak area using a linear regression curve obtained from separate internal standard solutions. The limits of detection (LOD) for the urinary PAHs metabolites were in the range 0.1–1.4 μg/L; default values were replaced with 50% of the LOD. Valid urinary PAHs metabolite concentrations were calibrated by levels of urinary creatinine and calculated as μg/mmol creatinine. Because 6-OHChr and 3-OHBaP were below the limits of quantification, we only analyzed 10 metabolites of PAHs.

Detection of IL-6 and Hsp70 in Plasma

IL-6 and Hsp70 levels were measured in EDTA-plasma using IL-6 and Hsp70 enzyme-linked immunosorbent assay (ELISA) kits (NeoBioscience Technology Company, China and Stressgen Bioreagents Company, Victoria, BC, Canada). The Hsp70 assay only detects inducible Hsp70 and does not detect other Hsp70 family members such as constitutive Hsp70, Grp78, DnaK (Escherichia coli), or Hsp71 (Mycobacterium tuberculosis). The inter-assay and intra-assay coefficients of variation of the assays were <10%.

Measurements of HRV

The measurement of HRV was conducted between 8∶00 am and 11∶00 am and measurement methods have been described previously [4]. After at least 5-min rest, each participant was seated comfortably on a chair and was fitted with a 3-channel digital Holter monitor (Lifecard CF; Del Mar Reynolds Medical, Inc., Irvine, USA) with a 1024 samples/second sampling rate for 10 minutes. We cleaned the participant’s skin with an alcohol wipe and abraded slightly to keep good lead contacts, and placed the separate electrode in the position according to the instruction and technical manual of Lifecard CF. All of the HRV indices were calculated on 5-min epoch in the entire recording. Only heart rates between 40 and 100 beats per minute were submitted to analyses [20]. We selected five consecutive minutes of ECG reading in the statistical analysis without atria and ventricular premature beats and flutter. The HRV spectrum was computed with a fast Fourier transform method and analyzed in both time and frequency domains. The measured time domain parameters include: a) SDNN (standard deviation of NN intervals in milliseconds), is an estimation of total HRV power; b) RMSSD (the root mean of square of successive differences between adjacent normal NN intervals in milliseconds), reflects the activities of parasympathetic nervous system. The frequency-domain variables include: a) low-frequency (LF msec2; 0.04–0.15Hz) may represent the combination of both parasympathetic and sympathetic activity of heart rate; b) high-frequency (HF msec2; 0.15–0.4Hz) shows the actions of the parasympathetic modulation of heart rate; c) total power (TP msec2; approximately≤0.4Hz) is a mixture of total variability of heart rate.

Statistical Analyses

We assessed the normality of all variables with the one-sample K-S test. Normal distributions of HRV indices and values of all urinary PAHs metabolites and plasma IL-6 and Hsp70 were obtained by natural logarithmic transformation. Because there is no well-established cutoff point for defining high versus low total OH-PAHs exposure, while 33rd and 67th percentiles [1] or the interquartile range (75th vs. 25th percentile) [21] was used as cutoffs in previous studies, we dichotomized total OH-PAHs at the 33th percentile as either low or high. Plasma IL-6 and Hsp70 were considered as categorical variables at different cut-off points. Baseline characteristics between urinary PAHs metabolites categories were compared using t-tests for continuous variables and chi-square tests for categorical variables. The trend for urinary PAHs metabolites levels in the four environmental exposure groups were tested by using simple linear regression. Pearson partial correlation coefficients (r) were computed to examine the associations among ln-transformed individual PAHs metabolites and ΣOH-PAHs. The trend for IL-6 and Hsp70 among quartiles of PAHs metabolites and the trend for HRV among quartiles of the IL-6 and Hsp70 were tested by multivariate linear regression models. Moreover, effects of IL-6 and Hsp70 on HRV stratified by PAHs metabolites were estimated using general linear models. Finally, to analyze potential interactions between IL-6 and PAHs metabolites as well as Hsp70 and PAHs metabolites, interaction product terms were included in the model. Potential confounders, including age, sex, length of work, smoking status, alcohol use, BMI, physical activity, working sites (top-oven, side and bottom-oven, adjunct-oven, and at the offices), and workshift (pre-shift and post-shift) for acute exposure – over the course of a day and weekday (Monday to Friday) for intermediate-term exposure – over the course of a week were used in the analysis. All statistical tests were derived from two-sided analyses and P<0.05 was considered as statistical significance. All analyses were conducted using SPSS (version 12.0).

Results

Characteristics of Study Population in Relation to Urinary Total OH-PAHs Levels

General characteristics of the workers classified by urinary total OH-PAHs levels are presented in Table 1. The distribution according to age, length of work, current smoker, current alcohol drinker, BMI and physical activity in the two groups was similar (all P>0.05), but differed in sex and workshift (pre-shift) (P<0.001). The median level of IL-6 was lower in the low total OH-PAHs group compared with the high total OH-PAHs group (1.86 pg/ml vs. 2.22 pg/ml), whereas Hsp70 level was similar. Averages of the five HRV indices were lower among those with the high total OH-PAHs levels, although SDNN and TP were not statistically significant.

[Figure omitted. See PDF.]

Table 1. General characteristics of the workers by urinary total OH-PAHs levels*.

https://doi.org/10.1371/journal.pone.0092964.t001

Urinary Levels of PAHs Metabolites at Different Working Sites

As our previous report [4], we defined office workers as the control group, and workers at adjunct-oven, and at side and bottom-oven, and top-oven as low, intermediate and high exposure groups, respectively. The concentrations of ΣOH-PAHs and nine urinary PAHs metabolites with the exception of 4-OHPh were significantly elevated with increasing environmental PAHs exposure levels (Table S1 in File S1, all Ptrend<0.01).

Correlations of Urinary PAHs Metabolites

Almost all individual PAHs metabolites were significantly correlated with each other, except 3-OHPh and 9-OHFlu as well as 2-OHPh, 3-OHPh and 4-OHPh. The correlation coefficients between individual PAHs metabolites and ΣOH-PAHs ranged from 0.262 to 0.851. The strongest correlations between individual PAHs metabolites and ΣOH-PAHs was 1-OHNa (r = 0.851), followed by 1-OHP and 2-OHNa (r = 0.836 and 0.772), the lowest correlations was 4-OHPh (r = 0.262) (Table S2 in File S1).

Effects of PAHs Metabolites on IL-6 and Hsp70 as Well as HRV

After controlling for potential confounders included age, sex, length of work, smoking status, alcohol use, BMI, physical activity, working sites, workshift and weekday, increasing quartiles of four urinary OH-PAHs (1-OHNa, 1-OHPh, 9-OHPh and 1-OHP) were significantly associated, in a dose-responsive manner, with increased IL-6 (Ptrend = 0.004, 0.008, 0.0001 and 0.001, respectively) (Figure 1). However, no significant dose-response relationship between PAHs metabolites and Hsp70 was found after multivariate adjustment.

[Figure omitted. See PDF.]

Figure 1. Association of urinary PAHs metabolites quartiles and ln-transformed IL-6.

Ptrend was calculated from the regression models, adjusting for age, sex, length of work, smoking status, alcohol use, BMI, physical activity, working sites, workshift and weekday. 1-OHNa: Q1 (≤0.92 μg/mmol creatinine), Q2 (0.92–1.57 μg/mmol creatinine), Q3 (1.57–2.80 μg/mmol creatinine), Q4 (>2.80 μg/mmol creatinine); 1-OHPh: Q1 (≤0.40 μg/mmol creatinine), Q2 (0.40–0.84 μg/mmol creatinine), Q3 (0.84–1.55 μg/mmol creatinine), Q4 (>1.55 μg/mmol creatinine); 9-OHPh: Q1 (≤0.39 μg/mmol creatinine), Q2 (0.39–0.70 μg/mmol creatinine), Q3 (0.70–1.35 μg/mmol creatinine), Q4 (>1.35 μg/mmol creatinine); 1-OHP: Q1 (≤1.86 μg/mmol creatinine), Q2 (1.86–3.29 μg/mmol creatinine), Q3 (3.29–6.11 μg/mmol creatinine), Q4 (>6.11 μg/mmol creatinine).

https://doi.org/10.1371/journal.pone.0092964.g001

In addition, increasing quartiles of 2-OHNa and 1-OHPh were significantly associated with decreased HF (Ptrend = 0.028 and 0.012, respectively), but no significant association between other metabolites and HRV indices was observed (Table S3 in File S1).

Dose-dependent Decrease in HRV Associated with Increased IL-6

As shown in Figure 2, after adjustment for covariates, increasing quartiles of IL-6 levels were significantly associated with decreases in TP and LF (Ptrend = 0.014 and 0.006). Compared to the lowest quartile, individuals with the highest quartile of IL-6 were significantly associated with 2.8% and 4.8% reduction in TP and LF.

[Figure omitted. See PDF.]

Figure 2. Values for natural log transformed TP (A) and LF (B) by quartiles of IL-6.

Ptrend was calculated from the regression models, adjusting for age, sex, length of work, smoking status, alcohol use, BMI, physical activity, working sites, workshift and weekday. Box plots indicate the 10th quartile, 25th quartile, median, 75th quartile, and 90th quartile for each group. IL-6: Q1 (≤1.25 pg/ml), Q2 (1.25–2.13 pg/ml), Q3 (2.13–3.71 pg/ml), Q4 (>3.71 pg/ml).

https://doi.org/10.1371/journal.pone.0092964.g002

Association of IL-6 with HRV by PAHs Metabolites

Given that HRV had significant association with IL-6 levels, we next examined whether the relationship of IL-6 with HRV differed in low- versus high-PAHs metabolites. As shown in Table 2, after adjustment for confounders, we found that elevated IL-6 levels were significantly associated with decreases in TP and LF in the high-PAHs metabolites groups (ΣOH-PAHs, 1-OHP and 2-OHPh, all Ptrend<0.05), but not in the low-PAHs metabolites groups. Significant interactions between IL-6 and ΣOH-PAHs on TP (Pinteraction = 0.017), and between IL-6 and 1-OHP, 2-OHPh and ΣOH-PAHs on LF (Pinteraction = 0.027, 0.044 and 0.001, respectively) were also observed.

[Figure omitted. See PDF.]

Table 2. Association of plasma IL-6 with HRV by urinary PAHs metabolites in coke oven workers*.

https://doi.org/10.1371/journal.pone.0092964.t002

Association of Hsp70 with HRV by PAHs Metabolites

No significant associations between Hsp70 and HRV were found in total population after adjusting for confounders (data not shown). When we examined the relationship of Hsp70 and HRV stratified by PAHs metabolites, we found that there was a significant increase in SDNN, TP and LF associated with Hsp70 when multiple PAHs metabolites were low (including 1-OHP, 1-OHNa, 2-OHPh, 9-OHPh and ΣOH-PAHs, all Ptrend <0.05), but not when they were high. Tests for interaction between Hsp70 and multiple PAHs metabolites were significant on SDNN, TP and LF (all Pinteraction<0.05) (Table 3).

[Figure omitted. See PDF.]

Table 3. Association of plasma Hsp70 with HRV by urinary PAHs metabolites in coke oven workers*.

https://doi.org/10.1371/journal.pone.0092964.t003

Discussion

Increasing epidemiological evidences have indicated significant relationship between inflammatory markers, most commonly CRP and IL-6, and various HRV indices in CHD patients [12], [13], type I diabetes [22], and in apparently normal adults [9], [11].However, until now, the association between inflammatory markers and HRV in workers exposed to different levels of PAHs has been less clear. We hypothesized that inflammation and stress response might play an important role in cardiac autonomic dysfunction induced by PAHs and potential dose-response relationships among them. In the present study, we used IL-6 as a marker to test this hypothesis, because that IL-6 is a good indicator of cytokine cascade activation, accurately reflects the inflammatory status, and is stable as its half-life is longer than other proinflammatory cytokines [23]. We found a significant dose-dependent relationship between four urinary OH-PAHs and IL-6; and IL-6 was significantly associated, in a dose-responsive manner, with decreased HRV, suggesting that certain OH-PAHs rather than others may be a significant contributor to the autonomic dysfunction through an inflammation mechanism. What is new and intriguing about our results is there were significant interactions between some PAHs metabolites and IL-6 on HRV. The link between IL-6 and HRV was only significant in workers with high-PAHs metabolites, but not in those with low-PAHs metabolites, indicating that high-PAHs exposure amplified the association of inflammation and HRV.

PAHs have been documented as one of the most important components of coke oven emission. Measurement of the urinary OH-PAHs metabolites, as internal markers, is an important way of assessing exposure to PAHs, since it takes into account all absorption routes of exogenous compounds [24], [25]. Due to relatively short half-lives of PAHs, the information provided by biomonitoring of urinary OH-PAH is limited to recent exposure [26], [27]. Most potential health effects are likely associated with exposure over time, thus considering the exposure time is essential for data interpretation, especially for chemicals with short half-lives. Li et al [26] evaluated the variability of nine urinary OH-PAHs metabolites over time in non-occupationally exposed subjects and found that the within-day variance exceeded the between-day variance, and the between-subject variance out-weighted the between-day variance for certain metabolites. However, our results did not substantially change even after adjusting for indicator variables of exposure time such as workshift for potential acute exposure (over the course of a day) and weekday intermediate-term exposure (over the course of a week), which strengthened the robustness of the observed association.

Studies have shown that exposure to high levels PAHs was positively associated with the prevalence, mortality or morbidity of cardiovascular disease[1], [28], and one mechanism by which PAHs may play a role in the development of cardiovascular disease is through inflammation [28]–[30]. Thus, high-PAHs exposure has been speculated to contribute to increasing systemic inflammation and subsequent reduced autonomic function, which might explain the significant relationship between IL-6 and HRV in high-PAHs exposure. The hypothesis is supported by studies that have identified that elevated PAHs levels were associated with increased IL-6 or CRP in non-occupational exposure population [29], [31], and systemic inflammatory activity was associated with a decrease in HRV [11], [32]. However, previous studies have only focused on one aspect of the relationships among PAHs, inflammatory biomarkers or HRV; no studies so far have examined all of them together, thus we cannot compare directly with that of other published studies. Luttmann-Gibson et al [33] showed that systemic inflammation modified susceptibility to air pollution-associated autonomic dysfunction in elderly subjects, but the cross-sectional nature of our study precludes any inference about the causal direction of the association between inflammation and HRV related to PAHs exposure. Together with the results of prior studies, these findings make it clear that significant decreased HRV is indeed associated with inflammation among workers with high-PAHs exposure, implicating that pollution can lead to autonomic dysfunction through inflammation.

Intracellular Hsp70 has been found to provide a spectrum of protection against any of a variety of stresses, preventing damage of cells and tissues [34], [35], while plasma extracellular Hsp70 (eHsp70) has been shown to have protective or deleterious effects in cardiovascular events [15], [16], [36]–[39]. Some studies have shown that increased Hsp70 level was associated with low risk of cardiovascular diseases [15], [16], [36], [37]. In contrast, several studies including some from our lab suggested that eHsp70 might possibly play a harmful role in the pathogenesis and progression of cardiovascular diseases [38], [39]. Although autonomic dysfunction have been associated with and may even precede various forms of cardiovascular diseases [40], no study has reported the effect of plasma Hsp70 on HRV. In this study, we did not find significant association between Hsp70 and various HRV indices in total population. However, we observed significant and consistent interactions between multiple PAHs metabolites and Hsp70 on HRV. There was a dose-response increase in SDNN, LF and TP in relation to elevated Hsp70 levels when PAHs metabolites was low but not when it was high, suggesting that increased plasma Hsp70 may have played a protective role in autonomic function when workers were exposed to relatively low levels of PAHs. Our previous studies showed that exposure to coke oven emission resulted in a dose-dependent increase in levels of plasma Hsp70 [19], but we did not explore whether there was a reverse relationship between PAHs levels and Hsp70 when workers were exposed to low coke oven emission as the relatively small sample size (n = 101) [19]. We speculate that circulating Hsp70 levels can be in the range which has been shown to exert protective effect against light damage when exposure to low-PAHs levels, and subsequent increase of HRV.

Several limitations of the present study merit consideration. First, the cross-sectional design is not able to prove the existence of a causal relationship whether inflammation leads to decreased HRV or vice versa. Second, monitoring of the 5-min HRV represents only short-term cardiac autonomic nerve modulation and the definitive clinical significance should be elucidated. Third, we did not adjust for multiple testing because of a high correlation both in various HRV indices and in PAHs metabolites. Over adjustment for multiple comparisons may increase the type II error and render truly important associations insignificant. Lastly, we could not rule out residual and unmeasured confounders, although we adjusted for a wide range of confounding factors.

Conclusions

In this study, we observed the dose-response relationships between urinary OH-PAHs and IL-6 as well as IL-6 and HRV indices. In particular, a prominent effect of IL-6 on HRV was found in the high-PAHs metabolites groups, whereas high Hsp70 levels may have a protective effect on HRV in low-PAHs metabolites groups. However, further researches are needed to confirm these findings in prospective studies and elucidate the possible mechanisms of these associations.

Supporting Information

[Figure omitted. See PDF.]

File S1.

This file contains Table S1 through Table S3. Table S1, Urinary levels of PAHs metabolites at different working sites. Table S2, The correlation coefficients (r) among individual PAHs metabolites and ΣOH-PAHs. Table S3, Effects of PAHs metabolites on HRV.

https://doi.org/10.1371/journal.pone.0092964.s001

(DOC)

Acknowledgments

We are particularly grateful to all volunteers for participating in the present study and staffs in the Institute of Industrial Health, Wuhan Iron & Steel Corporation.

Author Contributions

Conceived and designed the experiments: JY RZ X. Zhang. Performed the experiments: JY RZ XH YF LY X. Zhu QD. Analyzed the data: JY RZ X. Zhang. Contributed reagents/materials/analysis tools: TW X. Zhang. Wrote the paper: JY RZ X. Zhang. Revised the manuscript: JY RZ TW X. Zhang.

Citation: Ye J, Zhu R, He X, Feng Y, Yang L, Zhu X, et al. (2014) Association of Plasma IL-6 and Hsp70 with HRV at Different Levels of PAHs Metabolites. PLoS ONE 9(4): e92964. https://doi.org/10.1371/journal.pone.0092964

References

1. Xu X, Cook RL, Ilacqua VA, Kan H, Talbott EO, et al. (2010) Studying associations between urinary metabolites of polycyclic aromatic hydrocarbons (PAHs) and cardiovascular diseases in the United States. Sci Total Environ 408: 4943–4948.

2. International Programme on Chemical Safety (IPCS) (1998) Environmental health criteria 202. Selected non-heterocyclic polycyclic aromatic hydrocarbons. WHO, Geneva.

3. Lee MS, Magari S, Christiani DC (2011) Cardiac autonomic dysfunction from occupational exposure to polycyclic aromatic hydrocarbons. Occup Environ Med 68: 474–478.

4. Li X, Feng Y, Deng H, Zhang W, Kuang D, et al. (2012) The dose-response decrease in heart rate variability: any association with the metabolites of polycyclic aromatic hydrocarbons in coke oven workers? PLoS One 7: e44562.

5. Rich DQ, Zareba W, Beckett W, Hopke PK, Oakes D, et al. (2012) Are ambient ultrafine, accumulation mode, and fine particles associated with adverse cardiac responses in patients undergoing cardiac rehabilitation? Environ Health Perspect 120: 1162–1169.

6. Weichenthal S, Kulka R, Dubeau A, Martin C, Wang D, et al. (2011) Traffic-related air pollution and acute changes in heart rate variability and respiratory function in urban cyclists. Environ Health Perspect 119: 1373–1378.

7. Wu S, Deng F, Niu J, Huang Q, Liu Y, et al. (2010) Association of heart rate variability in taxi drivers with marked changes in particulate air pollution in Beijing in 2008. Environ Health Perspect 118: 87–91.

8. Zanobetti A, Gold DR, Stone PH, Suh HH, Schwartz J, et al. (2010) Reduction in heart rate variability with traffic and air pollution in patients with coronary artery disease. Environ Health Perspect 118: 324–330.

9. Lampert R, Bremner JD, Su S, Miller A, Lee F, et al. (2008) Decreased heart rate variability is associated with higher levels of inflammation in middle-aged men. Am Heart J 156: 759 e751–757.

10. Sloan RP, McCreath H, Tracey KJ, Sidney S, Liu K, et al. (2007) RR interval variability is inversely related to inflammatory markers: the CARDIA study. Mol Med 13: 178–184.

11. von Kanel R, Nelesen RA, Mills PJ, Ziegler MG, Dimsdale JE (2008) Relationship between heart rate variability, interleukin-6, and soluble tissue factor in healthy subjects. Brain Behav Immun 22: 461–468.

12. Haensel A, Mills PJ, Nelesen RA, Ziegler MG, Dimsdale JE (2008) The relationship between heart rate variability and inflammatory markers in cardiovascular diseases. Psychoneuroendocrinology 33: 1305–1312.

13. von Kanel R, Carney RM, Zhao S, Whooley MA (2011) Heart rate variability and biomarkers of systemic inflammation in patients with stable coronary heart disease: findings from the Heart and Soul Study. Clin Res Cardiol 100: 241–247.

14. Asea A, Kraeft SK, Kurt-Jones EA, Stevenson MA, Chen LB, et al. (2000) HSP70 stimulates cytokine production through a CD14-dependant pathway, demonstrating its dual role as a chaperone and cytokine. Nat Med 6: 435–442.

15. Pockley AG, Georgiades A, Thulin T, de Faire U, Frostegard J (2003) Serum heat shock protein 70 levels predict the development of atherosclerosis in subjects with established hypertension. Hypertension 42: 235–238.

16. Zhu J, Quyyumi AA, Wu H, Csako G, Rott D, et al. (2003) Increased serum levels of heat shock protein 70 are associated with low risk of coronary artery disease. Arterioscler Thromb Vasc Biol 23: 1055–1059.

17. Zhang X, Xu Z, Zhou L, Chen Y, He M, et al. (2010) Plasma levels of Hsp70 and anti-Hsp70 antibody predict risk of acute coronary syndrome Cell Stress and Chaperones. 15: 675–686.

18. Satoh M, Shimoda Y, Akatsu T, Ishikawa Y, Minami Y, et al. (2006) Elevated circulating levels of heat shock protein 70 are related to systemic inflammatory reaction through monocyte Toll signal in patients with heart failure after acute myocardial infarction. Eur J Heart Fail 8: 810–815.

19. Yang X, Zheng J, Bai Y, Tian F, Yuan J, et al. (2007) Using lymphocyte and plasma Hsp70 as biomarkers for assessing coke oven exposure among steel workers. Environ Health Perspect 115: 1573–1577.

20. Magari SR, Hauser R, Schwartz J, Williams PL, Smith TJ, et al. (2001) Association of heart rate variability with occupational and environmental exposure to particulate air pollution. Circulation 104: 986–991.

21. Clark JD 3rd, Serdar B, Lee DJ, Arheart K, Wilkinson JD, et al. (2012) Exposure to polycyclic aromatic hydrocarbons and serum inflammatory markers of cardiovascular disease. Environ Res 117: 132–137.

22. Gonzalez-Clemente JM, Vilardell C, Broch M, Megia A, Caixas A, et al. (2007) Lower heart rate variability is associated with higher plasma concentrations of IL-6 in type 1 diabetes. Eur J Endocrinol 157: 31–38.

23. Tateishi Y, Oda S, Nakamura M, Watanabe K, Kuwaki T, et al. (2007) Depressed heart rate variability is associated with high IL-6 blood level and decline in the blood pressure in septic patients. Shock 28: 549–553.

24. Huang W, Caudill SP, Grainger J, Needham LL, Patterson DG Jr (2006) Levels of 1-hydroxypyrene and other monohydroxy polycyclic aromatic hydrocarbons in children: a study based on U.S. reference range values. Toxicol Lett 163: 10–19.

25. Li Z, Sandau CD, Romanoff LC, Caudill SP, Sjodin A, et al. (2008) Concentration and profile of 22 urinary polycyclic aromatic hydrocarbon metabolites in the US population. Environ Res 107: 320–331.

26. Li Z, Romanoff LC, Lewin MD, Porter EN, Trinidad DA, et al. (2010) Variability of urinary concentrations of polycyclic aromatic hydrocarbon metabolite in general population and comparison of spot, first-morning, and 24-h void sampling. J Expo Sci Environ Epidemiol 20: 526–535.

27. Riojas-Rodriguez H, Schilmann A, Marron-Mares AT, Masera O, Li Z, et al. (2011) Impact of the improved patsari biomass stove on urinary polycyclic aromatic hydrocarbon biomarkers and carbon monoxide exposures in rural Mexican women. Environ Health Perspect 119: 1301–1307.

28. Burstyn I, Kromhout H, Partanen T, Svane O, Langard S, et al. (2005) Polycyclic aromatic hydrocarbons and fatal ischemic heart disease. Epidemiology 16: 744–750.

29. Everett CJ, King DE, Player MS, Matheson EM, Post RE, et al. (2010) Association of urinary polycyclic aromatic hydrocarbons and serum C-reactive protein. Environ Res 110: 79–82.

30. Chuang HC, Fan CW, Chen KY, Chang-Chien GP, Chan CC (2012) Vasoactive alteration and inflammation induced by polycyclic aromatic hydrocarbons and trace metals of vehicle exhaust particles. Toxicol Lett 214: 131–136.

31. Delfino RJ, Staimer N, Tjoa T, Arhami M, Polidori A, et al. (2010) Association of biomarkers of systemic inflammation with organic components and source tracers in quasi-ultrafine particles. Environ Health Perspect 118: 756–762.

32. Tracey KJ (2009) Reflex control of immunity. Nat Rev Immunol 9: 418–428.

33. Luttmann-Gibson H, Suh HH, Coull BA, Dockery DW, Sarnat SE, et al. (2010) Systemic inflammation, heart rate variability and air pollution in a cohort of senior adults. Occup Environ Med 67: 625–630.

34. Giffard RG, Han RQ, Emery JF, Duan M, Pittet JF (2008) Regulation of apoptotic and inflammatory cell signaling in cerebral ischemia - The complex roles of heat shock protein 70. Anesthesiology 109: 339–348.

35. Bielecka-Dabrowa A, Barylski M, Mikhailidis DP, Rysz J, Banach M (2009) HSP 70 and atherosclerosis protector - or activator? Expert Opinion on Therapeutic Targets 13: 307–317.

36. Dulin E, Garcia-Barreno P, Guisasola MC (2010) Extracellular heat shock protein 70 (HSPA1A) and classical vascular risk factors in a general population. Cell Stress Chaperones 15: 929–937.

37. Pockley AG, Calderwood SK, Multhoff G (2009) The atheroprotective properties of Hsp70: a role for Hsp70-endothelial interactions? Cell Stress Chaperones 14: 545–553.

38. Dybdahl B, Slordahl SA, Waage A, Kierulf P, Espevik T, et al. (2005) Myocardial ischaemia and the inflammatory response: release of heat shock protein 70 after myocardial infarction. Heart 91: 299–304.

39. Zhang X, Xu Z, Zhou L, Chen Y, He M, et al. (2010) Plasma levels of Hsp70 and anti-Hsp70 antibody predict risk of acute coronary syndrome. Cell Stress Chaperones 15: 675–686.

40. Gidron Y, Kupper N, Kwaijtaal M, Winter J, Denollet J (2007) Vagus-brain communication in atherosclerosis-related inflammation: a neuroimmunomodulation perspective of CAD. Atherosclerosis 195: e1–9.