You are here

  • Home
  • Archive
  • Volume 74, Issue 11
  • Application of linear mixed-effects model with LASSO to identify metal components associated with cardiac autonomic responses among welders: a repeated measures study

Article Text

Article menu

Download PDFPDF

XML
Workplace
Original Article
Application of linear mixed-effects model with LASSO to identify metal components associated with cardiac autonomic responses among welders: a repeated measures study
  1. Jinming Zhang1,
  2. Jennifer M Cavallari2,
  3. Shona C Fang3,
  4. Marc G Weisskopf4,
  5. Xihong Lin5,
  6. Murray A Mittleman4,6,
  7. David C Christiani1,7
  1. 1 Department of Environmental Health, Harvard T H Chan School of Public Health, Boston, USA
  2. 2 Department of Community Medicine and Health Care, University of Connecticut Health Center, Farmington, USA
  3. 3 Department of Epidemiology, New England Research Institute, Watertown, USA
  4. 4 Department of Epidemiology, Harvard T H Chan School of Public Health, Boston, USA
  5. 5 Department of Biostatistics, Harvard T H Chan School of Public Health, Boston, USA
  6. 6 Cardiovascular Epidemiology Research Unit, Beth Israel Deaconess Medical Center/Harvard Medical School, Boston, USA
  7. 7 Pulmonary and Critical Care Division, Massachusetts General Hospital/Harvard Medical School, Boston, USA
  1. Correspondence to Dr David C Christiani, 665 Huntington Avenue, Building I Room 1401, Boston, Massachusetts 02115, USA; dchris{at}hsph.harvard.edu

Abstract

Background Environmental and occupational exposure to metals is ubiquitous worldwide, and understanding the hazardous metal components in this complex mixture is essential for environmental and occupational regulations.

Objective To identify hazardous components from metal mixtures that are associated with alterations in cardiac autonomic responses.

Methods Urinary concentrations of 16 types of metals were examined and ‘acceleration capacity’ (AC) and ‘deceleration capacity’ (DC), indicators of cardiac autonomic effects, were quantified from ECG recordings among 54 welders. We fitted linear mixed-effects models with least absolute shrinkage and selection operator (LASSO) to identify metal components that are associated with AC and DC. The Bayesian Information Criterion was used as the criterion for model selection procedures.

Results Mercury and chromium were selected for DC analysis, whereas mercury, chromium and manganese were selected for AC analysis through the LASSO approach. When we fitted the linear mixed-effects models with ‘selected’ metal components only, the effect of mercury remained significant. Every 1 µg/L increase in urinary mercury was associated with −0.58 ms (−1.03, –0.13) changes in DC and 0.67 ms (0.25, 1.10) changes in AC.

Conclusion Our study suggests that exposure to several metals is associated with impaired cardiac autonomic functions. Our findings should be replicated in future studies with larger sample sizes.

  • metals
  • cardiovascular
  • epidemiology

This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/

https://doi.org/10.1136/oemed-2016-104067

Statistics from Altmetric.com

    Request Permissions

    If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

    What this paper adds

    • Welder are exposed to a mixture of variety types of metals existed in the welding fume during the welding process.

    • Previous studies have indicated that welders have higher risks of cardiovascular autonomic dysfunction due to occupational welding fumes exposure compared with the general population.

    • This paper suggests that occupational mercury, chromium and manganese exposures may be associated with higher risk of cardiac autonomic dysfunction among welders.

    • Interventions should aim at reducing metal-rich welding fumes exposure to mitigate cardiac autonomic alterations in the occupational setting.

    Introduction

    In recent years, toxic metal pollution has become one of the major environmental problems across the globe due to increase in industrialisation and anthropogenic activities. There has been increasing concern over the impact of metals exposure on cardiovascular diseases (CVD).1–3 Studies have shown that alterations in cardiac autonomic responses might be one of the potential mechanisms of cardiovascular effects of metals exposure.4 5 The imbalance between sympathetic and parasympathetic activities is associated with increased cardiovascular morbidity and mortality.6 7

    Conventional heart rate variability (HRV) parameters, including SD of NN intervals (SDNN), the root mean square of heartbeat interval differences (RMSSDs), high frequency (HF) and low frequency (LF), have been widely used as indicators of cardiac autonomic responses. However, recent studies suggested that these HRV parameters might not be accurate as they failed to account for heart rate, a potential confounder.8 One possible solution might be the phase-rectified signal averaging (PRSA) approach introduced by Bauer et al.9 Using this method, the deceleration and acceleration were separately characterised through two novel markers: acceleration capacity (AC) and deceleration capacity (DC). DC provides a measure of parasympathetic modulation of heart rate. Decreased DC, an indicator for the impaired cardiac autonomic functions, has been a strong predictor of cardiovascular mortality among patients after myocardial infarction.9 10

    Welders are exposed to various types of metals from welding fumes. Welding materials often consist of mixtures of metals, such as iron, manganese, chromium and nickel.11 Metal-rich welding fumes are small particles consisting of mostly vapourised metals that are suspended in air.12 Previous studies have reported a higher risk of cardiac autonomic dysfunction associated with fine particulate matters (PM2.5) exposure from welding fumes.13 14 However, few studies have examined the health effects of specific metal components. Understanding the hazardous metal species responsible for the adverse cardiac effects will have important implications for regulations in occupational settings. Therefore, in the current study, we examined the association between 16 metal species and AC and DC values. To address the challenge of disentangling effects of correlated metals in a mixture in the welding fume, we conducted model selection procedures via least absolute shrinkage and selection operator (LASSO) to identify metal components that are significantly associated with alterations in cardiac autonomic functions.

    Methods

    Study population

    The study population included a total of 54 welders from a local boilermaker union in Quincy, Massachusetts. Participants were recruited and monitored over four sampling occasions in the summer and winter months — June 2010, January 2011, June 2011 and June 2012, periods of boilermaker ‘off-season’ at the union welding school — and were invited to participate in multiple sampling occasions. All participants were male workers who (1) were 18 years of age or older; (2) were union apprentice or journeyman; (3) were free of CVD at first entry; and (4) provided urine and ECG samples at least once over the four sampling occasions.

    We monitored each participant for approximately 6 hours during a ‘work-shift’ at each sampling occasion. The major tasks were welding, cutting and grinding activities at a union welding school. Participants performed electric arc welding, using base metals of mild steel and stainless steel. Each participant provided a urine sample and resting ECG recordings at both prior (baseline) and post work. We collected questionnaires on demographics including age, height, weight, smoking status, medical history and medication use from participants at baseline. We obtained informed written consent from each participant prior to participation. The Harvard TH Chan School of Public Health Institutional Review Board approved our study protocol.

    ECG recordings and sample analysis

    During each sampling occasion, participants were fitted with a 7-lead ambulatory ECG Holter monitor. We collected 12 min resting ECG recordings from each participant at both prior and post work-shift. During the 12 min resting period, participants were asked to remain seated and quiet; walking, talking or eating was not allowed, as AC and DC are sensitive to these activities. The ECG samples were then sent to the Cardiovascular Epidemiology Research Unit of Beth Israel Deaconess Medical Center (Boston, Massachusetts, USA), where trained technicians blinded to exposure status processed and analysed these samples. The methods of processing ECG samples for AC and DC analysis have been discussed in a prior study.15 AC and DC quantities were computed and summarised every 5 min through the PRSA method, as described by Bauer et al. In brief, to compute DC values, RR intervals longer than the immediate preceding interval were defined as anchors; to compute AC values, RR intervals shorter than the immediate preceding interval were defined as anchors. Segments of interval data that had the same size around these anchors were identified. On alignment of all segments, RR intervals at all defined anchors (X0), immediately preceding (X1) and following the anchors (X-1) were averaged separately. The quantities of AC or DC were obtained by computing the difference between the sum of X0 and X1 and the sum of (X-1) and (X-2).9 To account for acclimation, the first 2 min of the 12 min resting ECG recordings was discarded and the remaining 10 min recordings were summarised as two non-overlapping 5 min AC or DC values (3rd–7th and 8th–12th minute intervals). We collected a total of 280 five-minute AC or DC values.

    Measurement of urinary metals and creatinine

    During each sampling occasion, participants provided urine samples at both baseline and postexposure. A total of 171 samples were collected during the four sampling occasions. Each urine sample was collected with a 120 mL sterile urine collection cup and temporarily kept on ice. An approximately 12 mL aliquot of urine was then transferred to 15 mL Falcon tubes for storage in a freezer at −20°C. The urine samples were sent to Brooks Rand Labs, where the urinary concentrations of 16 types of metals were determined. The inductively coupled plasma mass spectrometry (ICP-MS) was used to analyse aluminium (Al), copper (Cu), manganese (Mn), lead (Pb), zinc (Zn), cobalt (Co), vanadium (V), magnesium (Mg), molybdenum (Mo) and cadmium (Cd) concentrations, while the ICP-MS with dynamic reaction cell technology was used for iron (Fe), chromium (Cr), arsenic (As), selenium (Se) and nickel (Ni) analysis. In addition, the urinary mercury (Hg) concentration was determined using cold-vapour atomic fluorescence spectrometry. Measurements of both baseline and postexposure urinary metal species were obtained and used. For quality control purpose, duplicate, matrix spike and matrix spike duplicate analyses, as well as NIST SRM 1643e (certified reference material with trace elements in natural water), were used. The method detection limit for urinary metals measurements ranged from 0.01 to 3.30 µg/L. In addition, urinary creatinine levels were also determined.

    Statistical analysis

    In this study, we examined the urinary concentrations of 16 types of metals: Cr, Fe, Cu, Mn, Ni, Pb, V, Zn, Cd, As, Mg, Hg, Co, Mo, Se and Al, including metals that are commonly found in the welding fumes, as well as trace metals that have been associated with CVDs from previous research. A group of covariates were controlled for, including baseline age (a continuous variable in years), body mass index (BMI; calculated as weight in kilograms divided by height squared in metres), current smoking status (categorised into two groups: non-smoker, and former smoker or current smoker), time (baseline or postwork) and season (summer or winter when each sampling occasion occurred). These variables have been reported as predictors of cardiac autonomic responses. A variable ‘time in day’ was also adjusted for as AC and DC might exhibit a circadian pattern throughout the day.16 In addition, urinary creatinine level was included in the model as a covariate as studies suggested the potential bias might be introduced when using the creatinine-standardised metal concentration (the ratio of metal concentrations over creatinine levels).17

    Tibshirani18 introduced the idea of ‘least absolute shrinkage and selection operator’ (LASSO), which applies an L1 penalty to ‘minimize the residual sum of squares subject to the sum of the absolute value of the coefficients being less than a constant’.18 LASSO can be used when the study sample size is very small or when the number of parameters is larger than the number of samples. We conducted linear mixed-effects models with LASSO to further account for correlations of repeated measurements from the same subject. The linear mixed-effects models with random subject-specific intercepts assessing metals exposure on AC or DC changes (Yij) were fitted as follows:

    Yij = Xi β+γ1 Ageij+γ2 BMIij+γ3 Time in dayij+γ4 smoking statusij+γ5 Timeij+γ6 Seasonij+bi+εij

    where Yij indicates AC or DC measurements of subject i at jth time, εij~N(0, σ2 and bi ~N(0,σb 2), β=(βi1,…, βij)T indicates a vector of regression coefficients of metals species vector Xi=(Xi1,…, Xij) that is ‘selected’ via LASSO with L1 penalty. To identify and select metal component species associated with AC and DC, a penalty term was applied to all urinary metal species only, but not the other covariates (ie, they were ‘fixed’ covariates in the model).

    When L1 penalty is applied, a non-negative tuning parameter λ is used to determine the amount of shrinkage or penalisation. For example, when λ is small, the model coefficients of all metal species are weakly penalised and are similar to those estimated from a regular linear mixed-effects model; when λ is large enough, the model coefficients of all metal species are strongly penalised, the shrinkage tends towards the null (0), and the final model will only include the ‘fixed’ covariates; when λ takes a value in between the extreme values, the model coefficients of some metal species are 0, while those of the others are non-zero in the penalised model. In this study, we ran all models across the range of λs. The Bayesian Information Criterion (BIC)19 was used as the criterion for model selection, and the λ in a model with the smallest BIC was identified and selected. Since LASSO did not provide SEs for the selected metal species, we then fitted a final linear mixed-effects model including selected metal species and aforementioned ‘fixed’ covariates. In addition, we also conducted linear mixed-effects models to estimate pairwise correlations between metal species. All statistical analyses were performed with SAS V.9.3 and R V.3.2.2 (R Core Team 2015).

    Results

    Demographic descriptions of the study population are summarised in table 1. Among 54 male welders, there were 22 (40.7%) current smokers and 32 (59.3%) former smokers or non-smokers. The average age was 41.6 (ranged from 21.7 to 71.2 years). The average BMI was 28 kg/m2. The median urinary metals concentrations ranged from 0.18 μg/L to 77.55 mg/L at baseline and 0.19 μg/L to 82.30 mg/L at post work. The distribution of urinary metals concentrations was skewed (table 2).

    View this table:
    Table 1

    Demographics of study populations at first entry (n=54)

    View this table:
    Table 2

    Summaries of urinary metals concentrations at preshift and postshift

    The relationship between λ and model BIC, the criterion for model selection procedures, is shown in figure 1. For DC analysis, there was maximum shrinkage when λ=20; the model with the smallest BIC (λ=9) included Hg and Cr only, and the regression coefficients of all the other metal components were 0. For AC analysis, the model with the smallest BIC (λ=6) included Mn, in addition to Hg and Cr, and there was maximum shrinkage when λ=21.

    Figure 1
    Figure 1

    The relationship between λ, the penalisation factor, and model BIC, the criterion for model selection. AC, acceleration capacity; BIC, Bayesian Information Criterion; DC, deceleration capacity.

    When we fitted the linear mixed-effects models with ‘fixed’ covariates and ‘selected’ metal components only, there were negative associations of urinary Hg and Cr levels with DC, whereas there were positive associations of urinary Hg, Cr and Mn levels with AC. However, only the effects of Hg exposure exhibited statistical significance, and the effects of Cr and Mn exposure were marginally significant. A 1 µg/L increase in urinary Hg was associated with −0.58 ms (−1.03, –0.13) change in DC and 0.67 ms (0.25, 1.10) change in AC (table 3). The correlation coefficients of urinary Hg and Cr, Cr and Mn, as well as Hg and Mn that were estimated through linear mixed-effects models were 0.46, 0.37 and 0.29, respectively.

    View this table:
    Table 3

    Effects of copollutants exposure on AC or DC changes

    Discussion

    Overall our findings suggested that exposure to metals was associated with adverse cardiac autonomic responses. A total of 16 urinary metal species were examined, and their half-lives in urine vary from a few hours to a few months. We used linear mixed-effects models with LASSO shrinkage method to identify metal components that might be associated with AC or DC. We observed that urinary Hg and Cr concentrations were negatively associated with DC, whereas urinary Hg, Cr and Mn concentrations were positively associated with AC. The significant effects of Hg exposure on AC and DC persisted in the copollutant models, while the effects of Cr and Mn yielded marginal significance. Our findings suggested that metals exposure might induce imbalance between sympathetic and parasympathetic activities.

    Mercury is classified in three major groups: elemental mercury, inorganic mercury compounds and organic mercury. To date, several biomarkers have been used to assess Hg exposure, such as hair, blood and urine samples.20 21 Hg levels in blood and hair are indicators of organic Hg exposures, whereas urinary Hg concentration mainly reflects exposure to elemental or inorganic forms of Hg. In the current study, we used urinary Hg concentrations to estimate medium-term Hg exposure, as it has a half-life of approximately 1–3 months in urine. The baseline median value (0.52 µg/L) was similar to the US population medians from National Health22 and Nutrition Examination Survey (0.48 µg/L)22 and 17.5% lower than the values (0.63 µg/L) from a study in a dental professionals cohort.23

    While Hg has long been considered as an environmental and occupational neurotoxicant, there has been increasing concern about the cardiotoxic effects of Hg exposure. Several studies have documented the role of organic Hg exposure in triggering cardiac autonomic responses. For example, Lim et al 24 examined hair Hg levels and HRV in a cross-sectional study among 1589 residents living near industrial complexes in South Korea, and it was found that a 1 parts per million (ppm) increase in hair Hg concentration was associated with an 8.4% decline in HF parameter. In a clinical trial where 27 participants in the intervention group were instructed to consume bigeye tuna and swordfish once a week for 14 weeks, while 27 participants in the control group continued usual diets, participants in the intervention group had significantly higher ratio of LF to HF at the end of the study over the control group.25 On the other hand, epidemiological research investigating the cardiovascular effects of elemental or inorganic Hg exposure has mainly focused on blood pressure.26 27 To our knowledge, this is the first study reporting the association of element or inorganic Hg exposure and cardiac autonomic responses in humans.

    Furthermore, it is notable that these welders might be exposed to elemental or inorganic mercury from ambient environment rather than occupational settings, as Hg is not a common welding fume component. Studies have reported ubiquitous exposure to inorganic forms of Hg among the population in USA and Canada.22 28 Elemental mercury is often used in products like thermometers and dental amalgam.29 30 In the general population, elemental mercury from dental amalgam is one of the major sources of Hg exposure. Inorganic mercury compounds have been widely used in batteries or organic chemicals. In addition, urinary mercury levels could also change due to high fish consumption.31

    On the other hand, Cr and Mn are both common components that exist in welding fumes. These vapourised metal components become small particles that could suspend and react with oxygen in the air throughout the welding process. Cr can be found in two forms: trivalent chromium (Cr(III)) or hexavalent chromium (Cr(VI)). Cr(VI) is a well-established environmental pollutant and is much more toxic than Cr(III).32 Inhalation, digestion and dermal absorption were the main pathway of occupational Cr exposure. In the current study, we measured short-term Cr(VI) exposure through urinary Cr(VI) concentrations, which reflect exposure over approximately 1–2 days.33 We observed a significant association between occupational Cr exposures and impaired cardiac autonomic effects. Similarly, Cavallari et al 34 found that airborne Cr(VI) exposure was associated with reduced nocturnal RMSSD among this welder population. However, evidence from epidemiological studies was very limited and the cardiotoxicity of Cr(VI) exposure warrants further research.

    In addition, urinary Mn concentrations were used to estimate recent airborne manganese from welding fumes among these welders.35 Previous studies have documented the association of Mn exposure and annual cardiovascular mortality, as well as Mn being a modifier of the association between PM2.5 exposure and hospital admissions for CVD.36 Consistent with research on airborne Mn among manganese alloy workers37 and boilermaker construction workers,34 our studies suggested that urinary Mn levels were associated with impaired cardiac autonomic functions.

    Although the underlying mechanisms of metals exposure on cardiovascular effects are still under investigation, studies suggested that oxidative stress, inflammatory responses and enzymatic inhibition38–40 might play a role. Transition metals, such as Hg, Cr and Mn, may undergo redox cycling reactions and induce oxidative stress through enhanced generation of reactive oxygen species.38 Low concentrations of Hg exposure may induce decreased nitrogen monoxide (NO) bioavailability due to oxidative stress and further promote endothelial dysfunction in both resistance and conductance arteries in male Wistar rats.41 Furthermore, decreased HRV parameters (SDNN, SD1 (standard deviation of the instantaneous beat-to-beat variability) and SD2 (long term variability of continuous RR intervals)) following peripheral injection of interleukin 6, a proinflammatory cytokine, were documented in another toxicological study.40 Taken together, future studies examining oxidative stress and inflammation as the mediators linking metals exposure and cardiovascular effects may be essential for illustrating the underlying causal pathway.

    To our knowledge, this is the first repeated measurements study investigating metal components exposure and cardiac autonomic effects. We used linear mixed-effects models with LASSO for model selection procedures, and this approach has several advantages over conventional methods. Typically, researchers examined the effects of metal components by fitting a model including all components simultaneously. However, collinearity issues due to high correlations among metal components may induce inflated variances of the estimated regression coefficients and sap the statistical power.42 Another commonly used approach, stepwise regression, may have similar problems in the presence of collinearity, as well as yield biased coefficients that need shrinkage.18

    When we fitted ‘selected’ metals only in the copollutant model, the effects of Cr and Mn were only marginally significant. This may be due to the relatively small sample size and the modest correlations among these metal components. It is also notable that metals may often induce antagonistic, additive or synergistic responses in the presence of other metal components. For example, Papp et al 43 observed synergistic neurotoxicity of combined exposure of Hg, Mn and Pb among male Wistar rats.43 While the exact relationships between Hg, Cr and Mn in triggering cardiac autonomic effects were unclear, it is possible that model misspecification may also occur in this study as we assumed independent relationships among these components in the copollutant model. Future studies may examine the overall effect of the mixture by accounting for the possible interplay among metal species.

    Several other limitations to our study should be considered. First, we acknowledge that fitting a linear mixed-effects model including selected metal species by LASSO may not be ideal for estimating the SEs or p values and may yield narrow CIs, since they fail to include a penalty term used in LASSO. Second, urinary concentrations of some metals44 45 may not be ideal for assessing metals exposures, and non-differential exposure misclassification may be possible, and measurement errors may also occur due to individual variability in the excretion of metal and creatinine in urine. Third, there may be healthy worker effect if these welders tend to be healthier compared with non-working population. Finally, our study had a relatively small sample size and the results from our study need to be replicated in additional studies of metal-exposed workers.

    Acknowledgments

    We thank the participants and leadership from Local 29 of International Brotherhood of Boilermaker (Quincy, MA).

    References

    1. ↵
      2. Solenkova NV ,
      3. Newman JD ,
      4. Berger JS ,
      5. Thurston G , et al
      . Metal pollutants and cardiovascular disease: mechanisms and consequences of exposure. Am Heart J 2014;168:812–22.doi:10.1016/j.ahj.2014.07.007
      OpenUrlCrossRefPubMed
    2. ↵
      2. Navas-Acien A ,
      3. Selvin E ,
      4. Sharrett AR , et al
      . Lead, cadmium, smoking, and increased risk of peripheral arterial disease. Circulation 2004;109:3196–201.doi:10.1161/01.CIR.0000130848.18636.B2
      OpenUrlAbstract/FREE Full Text
    3. ↵
      2. Alissa EM ,
      3. Ferns GA
      . Heavy metal poisoning and cardiovascular disease. J Toxicol 2011;2011:1–21.doi:10.1155/2011/870125
      OpenUrl
    4. ↵
      2. Feng W ,
      3. He X ,
      4. Chen M , et al
      . Urinary metals and heart rate variability: a cross-sectional study of urban adults in Wuhan, China. Environ Health Perspect 2015;123:217–22.doi:10.1289/ehp.1307563
      OpenUrl
    5. ↵
      2. Magari SR ,
      3. Schwartz J ,
      4. Williams PL , et al
      . The association of particulate air metal concentrations with heart rate variability. Environ Health Perspect 2002;110:875–80.doi:10.1289/ehp.02110875
      OpenUrlCrossRefPubMedWeb of Science
    6. ↵
      2. Nolan J ,
      3. Batin PD ,
      4. Andrews R , et al
      . Prospective study of heart rate variability and mortality in chronic heart failure: results of the United Kingdom heart failure evaluation and assessment of risk trial (UK-heart). Circulation 1998;98:1510–6.doi:10.1161/01.CIR.98.15.1510
      OpenUrlAbstract/FREE Full Text
    7. ↵
      2. Chandola T ,
      3. Britton A ,
      4. Brunner E , et al
      . Work stress and coronary heart disease: what are the mechanisms? Eur Heart J 2008;29:640–8.doi:10.1093/eurheartj/ehm584
      OpenUrlCrossRefPubMed
    8. ↵
      2. Monfredi O ,
      3. Lyashkov AE ,
      4. Johnsen AB , et al
      . Biophysical characterization of the underappreciated and important relationship between heart rate variability and heart rate. Hypertension 2014;64:1334–43.doi:10.1161/HYPERTENSIONAHA.114.03782
      OpenUrlCrossRef
    9. ↵
      2. Bauer A ,
      3. Kantelhardt JW ,
      4. Barthel P , et al
      . Deceleration capacity of heart rate as a predictor of mortality after myocardial infarction: cohort study. Lancet 2006;367:1674–81.doi:10.1016/S0140-6736(06)68735-7
      OpenUrlCrossRefPubMedWeb of Science
    10. ↵
      2. Schneider A ,
      3. Hampel R ,
      4. Ibald-Mulli A , 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
      OpenUrlCrossRefPubMed
    11. ↵
      2. Antonini JM
      . Health effects of welding. Crit Rev Toxicol 2003;33:61–103.doi:10.1080/713611032
      OpenUrlCrossRefPubMedWeb of Science
    12. ↵
      2. Wong JY ,
      3. Fang SC ,
      4. Grashow R , et al
      . The relationship between occupational metal exposure and arterial compliance. J Occup Environ Med 2015;57:355–60.doi:10.1097/JOM.0000000000000427
      OpenUrl
    13. ↵
      2. Fan T ,
      3. Fang SC ,
      4. Cavallari JM , et al
      . Heart rate variability and DNA methylation levels are altered after short-term metal fume exposure among occupational welders: a repeated-measures panel study. BMC Public Health 2014;14:1279.doi:10.1186/1471-2458-14-1279
    14. ↵
      2. Scharrer E ,
      3. Hessel H ,
      4. Kronseder A , et al
      . Heart rate variability, hemostatic and acute inflammatory blood parameters in healthy adults after short-term exposure to welding fume. Int Arch Occup Environ Health 2007;80:265–72.doi:10.1007/s00420-006-0127-2
      OpenUrlCrossRefPubMedWeb of Science
    15. ↵
      2. Umukoro PE ,
      3. Cavallari JM ,
      4. Fang SC , et al
      . Short-term metal particulate exposures decrease cardiac acceleration and deceleration capacities in welders: a repeated-measures panel study. Occup Environ Med 2016;73:91–6.doi:10.1136/oemed-2015-103052
      OpenUrlAbstract/FREE Full Text
    16. ↵
      2. Cavallari JM ,
      3. Fang SC ,
      4. Mittleman MA , et al
      . Circadian variation of heart rate variability among welders. Occup Environ Med 2010;67:717–9.doi:10.1136/oem.2010.055210
      OpenUrlAbstract/FREE Full Text
    17. ↵
      2. Barr DB ,
      3. Wilder LC ,
      4. Caudill SP , et al
      . Urinary creatinine concentrations in the U.S. population: implications for urinary biologic monitoring measurements. Environ Health Perspect 2005;113:192–200.doi:10.1289/ehp.7337
      OpenUrlPubMedWeb of Science
    18. ↵
      2. Tibshirani R
      . Regression shrinkage and selection via the lasso. J R Stat Soc Series B Stat Methodol 1996;58:267–88.
      OpenUrl
    19. ↵
      2. Schwarz G
      . Estimating the dimension of a Model. The Annals of Statistics 1978;6:461–4.doi:10.1214/aos/1176344136
      OpenUrlCrossRefWeb of Science
    20. ↵
      2. McDowell MA ,
      3. Dillon CF ,
      4. Osterloh J , et al
      . Hair mercury levels in U.S. children and women of childbearing age: reference range data from NHANES 1999-2000. Environ Health Perspect 2004;112:1165–71.doi:10.1289/ehp.7046
      OpenUrlCrossRefPubMedWeb of Science
    21. ↵
      2. Somers EC ,
      3. Ganser MA ,
      4. Warren JS , et al
      . Mercury exposure and antinuclear antibodies among females of Reproductive Age in the United States: nhanes. Environ Health Perspect 2015;123:792–8.doi:10.1289/ehp.1408751
      OpenUrl
    22. ↵
      2. Crinnion WJ
      . The CDC fourth national report on human exposure to environmental chemicals: what it tells Us about our toxic burden and how it assist environmental medicine physicians. Altern Med Rev 2010;15:101–9.
      OpenUrlPubMed
    23. ↵
      2. Goodrich JM ,
      3. Wang Y ,
      4. Gillespie B , et al
      . Methylmercury and elemental mercury differentially associate with blood pressure among dental professionals. Int J Hyg Environ Health 2013;216:195–201.doi:10.1016/j.ijheh.2012.03.001
      OpenUrl
    24. ↵
      2. Lim S ,
      3. Chung HU ,
      4. Paek D
      . Low dose mercury and heart rate variability among community residents nearby to an industrial complex in Korea. Neurotoxicology 2010;31:10–16.doi:10.1016/j.neuro.2009.10.001
      OpenUrlPubMed
    25. ↵
      2. Yaginuma-Sakurai K ,
      3. Murata K ,
      4. Shimada M , et al
      . Intervention study on cardiac autonomic nervous effects of methylmercury from seafood. Neurotoxicol Teratol 2010;32:240–5.doi:10.1016/j.ntt.2009.08.009
      OpenUrlCrossRefPubMed
    26. ↵
      2. Valera B ,
      3. Dewailly E ,
      4. Poirier P
      . Environmental mercury exposure and blood pressure among Nunavik Inuit adults. Hypertension 2009;54:981–6.doi:10.1161/HYPERTENSIONAHA.109.135046
      OpenUrlCrossRef
    27. ↵
      2. Kobal AB ,
      3. Horvat M ,
      4. Prezelj M , et al
      . The impact of long-term past exposure to elemental mercury on antioxidative capacity and lipid peroxidation in mercury miners. J Trace Elem Med Biol 2004;17:261–74.doi:10.1016/S0946-672X(04)80028-2
      OpenUrlCrossRefPubMedWeb of Science
    28. ↵
      2. Nicolae A ,
      3. Ames H ,
      4. Quiñonez C
      . Dental amalgam and urinary mercury concentrations: a descriptive study. BMC Oral Health 2013;13:44.doi:10.1186/1472-6831-13-44
    29. ↵
      2. Kingman A ,
      3. Albertini T ,
      4. Brown LJ
      . Mercury concentrations in urine and whole blood associated with amalgam exposure in a US military population. J Dent Res 1998;77:461–71.doi:10.1177/00220345980770030501
      OpenUrlCrossRefPubMedWeb of Science
    30. ↵
      2. Woods JS ,
      3. Martin MD ,
      4. Leroux BG , et al
      . The contribution of dental amalgam to urinary mercury excretion in children. Environ Health Perspect 2007;115:1527–31.doi:10.1289/ehp.10249
      OpenUrlPubMedWeb of Science
    31. ↵
      2. Passos CJ ,
      3. Mergler D ,
      4. Lemire M , et al
      . Fish consumption and bioindicators of inorganic mercury exposure. Sci Total Environ 2007;373:68–76.doi:10.1016/j.scitotenv.2006.11.015
      OpenUrlCrossRefPubMedWeb of Science
    32. ↵
      2. Berardi R ,
      3. Pellei C ,
      4. Valeri G , et al
      . Chromium exposure and germinal embryonal carcinoma: first two cases and review of the literature. J Toxicol Environ Health A 2015;78:1–6.doi:10.1080/15287394.2015.958416
      OpenUrl
    33. ↵
      2. Wilbur S ,
      3. Abadin H ,
      4. Fay M , et al
      . Toxicological Profile for Chromium. Atlanta, 2012.
    34. ↵
      2. Cavallari JM ,
      3. Eisen EA ,
      4. Fang SC , et al
      . PM2.5 metal exposures and nocturnal heart rate variability: a panel study of boilermaker construction workers. Environ Health 2008;7:36.doi:10.1186/1476-069X-7-36
      OpenUrlCrossRefPubMed
    35. ↵
      2. Williams M ,
      3. Todd GD ,
      4. Roney N , et al
      . Toxicological Profile for Manganese. Atlanta (GA) 2012.
    36. ↵
      2. Jiang Y ,
      3. Zheng W
      . Cardiovascular toxicities upon manganese exposure. Cardiovasc Toxicol 2005;5:345–54.doi:10.1385/CT:5:4:345
      OpenUrlCrossRefPubMed
    37. ↵
      2. Barrington WW ,
      3. Angle CR ,
      4. Willcockson NK , et al
      . Autonomic function in manganese alloy workers. Environ Res 1998;78:50–8.doi:10.1006/enrs.1997.3826
      OpenUrlPubMed
    38. ↵
      2. Jomova K ,
      3. Valko M
      . Advances in metal-induced oxidative stress and human disease. Toxicology 2011;283:65–87.doi:10.1016/j.tox.2011.03.001
      OpenUrlCrossRefPubMedWeb of Science
    39. ↵
      2. Chuang HC ,
      3. Hsueh TW ,
      4. Chang CC , et al
      . Nickel-regulated heart rate variability: the roles of oxidative stress and inflammation. Toxicol Appl Pharmacol 2013;266:298–306.doi:10.1016/j.taap.2012.11.006
      OpenUrl
    40. ↵
      2. Hajiasgharzadeh K ,
      3. Mirnajafi-Zadeh J ,
      4. Mani AR
      . Interleukin-6 impairs chronotropic responsiveness to cholinergic stimulation and decreases heart rate variability in mice. Eur J Pharmacol 2011;673:70–7.doi:10.1016/j.ejphar.2011.10.013
      OpenUrlCrossRefPubMed
    41. ↵
      2. Wiggers GA ,
      3. Peçanha FM ,
      4. Briones AM , et al
      . Low mercury concentrations cause oxidative stress and endothelial dysfunction in conductance and resistance arteries. Am J Physiol Heart Circ Physiol 2008;295:H1033–H1043.doi:10.1152/ajpheart.00430.2008
      OpenUrlAbstract/FREE Full Text
    42. ↵
      2. Tormod N ,
      3. Bjørn-Helge M
      . Understanding the collinearity problem in regression and discriminant analysis. Journal of Chemometrics 2001;15:413–26.
      OpenUrlCrossRefWeb of Science
    43. ↵
      2. Papp A ,
      3. Pecze L ,
      4. Szabó A , et al
      . Effects on the central and peripheral nervous activity in rats elicited by acute administration of lead, mercury and manganese, and their combinations. J Appl Toxicol 2006;26:374–80.doi:10.1002/jat.1152
      OpenUrlCrossRefPubMed
    44. ↵
      2. Paustenbach DJ ,
      3. Panko JM ,
      4. Fredrick MM , et al
      . Urinary chromium as a biological marker of environmental exposure: what are the limitations? Regul Toxicol Pharmacol 1997;26:S23–S34.doi:10.1006/rtph.1997.1135
      OpenUrlPubMed
    45. ↵
      2. Bader M ,
      3. Dietz MC ,
      4. Ihrig A , et al
      . Biomonitoring of manganese in blood, urine and axillary hair following low-dose exposure during the manufacture of dry cell batteries. Int Arch Occup Environ Health 1999;72:521–7.doi:10.1007/s004200050410
      OpenUrlCrossRefPubMed

    Footnotes

    • Contributors JZ contributed to all aspects of data acquisition, data analysis and interpretation, and preparation of the manuscript. DCC contributed to the study conception and design, data collection, analysis and interpretation. JMC and SCF contributed to the data collection and interpretation of the results. MAM, MGW and XL contributed to the data analysis and interpretation of the results. All authors reviewed and approved the manuscript.

    • Competing interests None declared.

    • Patient consent Obtained.

    • Ethics approval The Harvard T H Chan School of Public Health Institutional Review Board.

    • Provenance and peer review Not commissioned; externally peer reviewed.