Article Text


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. 1Department of Environmental Health, Harvard T H Chan School of Public Health, Boston, USA
  2. 2Department of Community Medicine and Health Care, University of Connecticut Health Center, Farmington, USA
  3. 3Department of Epidemiology, New England Research Institute, Watertown, USA
  4. 4Department of Epidemiology, Harvard T H Chan School of Public Health, Boston, USA
  5. 5Department of Biostatistics, Harvard T H Chan School of Public Health, Boston, USA
  6. 6Cardiovascular Epidemiology Research Unit, Beth Israel Deaconess Medical Center/Harvard Medical School, Boston, USA
  7. 7Pulmonary 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}


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:

Statistics from

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.


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.


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 Ageij2 BMIij3 Time in dayij4 smoking statusij5 Timeij6 Seasonij+biij

where Yij indicates AC or DC measurements of subject i at jth time, εij~N(0, σ2 and bi ~N(0,σb2), β=(βi1,…, βij)T indicates a vector of regression coefficients of metals species vector Xi=(Xi1,…, Xij) that is ‘selected’ via LASSO with Lpenalty. 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).


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).

Table 1

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

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

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.

Table 3

Effects of copollutants exposure on AC or DC changes

AC, acceleration capacity; DC, deceleration capacity.


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 al24 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 al34 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 al43 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.


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


View Abstract


  • 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.

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.