Article Text

Original research
Multiple metal exposures associate with higher amyotrophic lateral sclerosis risk and mortality independent of genetic risk and correlate to self-reported exposures: a case-control study
Free
  1. Dae-Gyu Jang1,2,
  2. John F Dou3,
  3. Emily J Koubek1,2,
  4. Samuel Teener1,2,
  5. Lili Zhou4,
  6. Kelly M Bakulski3,
  7. Bhramar Mukherjee5,
  8. Stuart A Batterman6,
  9. Eva L Feldman1,2,
  10. Stephen A Goutman1,2
  1. 1 Department of Neurology, University of Michigan, Ann Arbor, Michigan, USA
  2. 2 NeuroNetwork for Emerging Therapies, University of Michigan, Ann Arbor, Michigan, USA
  3. 3 Department of Epidemiology, University of Michigan, Ann Arbor, Michigan, USA
  4. 4 Department of Biostatistics, Corewell Health, Royal Oak, Michigan, USA
  5. 5 Department of Biostatistics, University of Michigan, Ann Arbor, Michigan, USA
  6. 6 Department of Environmental Health Sciences, University of Michigan, Ann Arbor, Michigan, USA
  1. Correspondence to Dr Stephen A Goutman; sgoutman{at}med.umich.edu

Abstract

Background The pathogenesis of amyotrophic lateral sclerosis (ALS) involves both genetic and environmental factors. This study investigates associations between metal measures in plasma and urine, ALS risk and survival and exposure sources.

Methods Participants with and without ALS from Michigan provided plasma and urine samples for metal measurement via inductively coupled plasma mass spectrometry. ORs and HRs for each metal were computed using risk and survival models. Environmental risk scores (ERS) were created to evaluate the association between exposure mixtures and ALS risk and survival and exposure source. ALS (ALS-PGS) and metal (metal-PGS) polygenic risk scores were constructed from an independent genome-wide association study and relevant literature-selected single-nucleotide polymorphisms.

Results Plasma and urine samples from 454 ALS and 294 control participants were analysed. Elevated levels of individual metals, including copper, selenium and zinc, significantly associated with ALS risk and survival. ERS representing metal mixtures strongly associated with ALS risk (plasma, OR=2.95, CI=2.38–3.62, p<0.001; urine, OR=3.10, CI=2.43–3.97, p<0.001) and poorer ALS survival (plasma, HR=1.37, CI=1.20–1.58, p<0.001; urine, HR=1.44, CI=1.23–1.67, p<0.001). Addition of the ALS-PGS or metal-PGS did not alter the significance of metals with ALS risk and survival. Occupations with high potential of metal exposure associated with elevated ERS. Additionally, occupational and non-occupational metal exposures were associated with measured plasma and urine metals.

Conclusion Metals in plasma and urine associated with increased ALS risk and reduced survival, independent of genetic risk, and correlated with occupational and non-occupational metal exposures. These data underscore the significance of metal exposure in ALS risk and progression.

  • EPIDEMIOLOGY
  • ALS

Data availability statement

Data are available upon reasonable request.

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 IS ALREADY KNOWN ON THIS TOPIC

  • Metal exposure is linked to amyotrophic lateral sclerosis (ALS) risk and survival; however, information regarding the impact of metal mixtures and the influence of underlying genetic heterogeneity on the relationship between metals and ALS is lacking.

WHAT THIS STUDY ADDS

  • This study demonstrates that mixtures of metals are significantly associated with increased risk of ALS and reduced survival, independent of genetic risk. Additionally, occupations linked to increased metal exposure correlate with higher environmental risk scores, supporting the relationship between occupational metal exposure and increased ALS risk.

HOW THIS STUDY MIGHT AFFECT RESEARCH, PRACTICE OR POLICY

  • This study reinforces the relationship between metal exposure and ALS pathophysiology and underscores the importance of further research into the underlying mechanisms driving this association to inform targeted risk reduction strategies.

Introduction

Amyotrophic lateral sclerosis (ALS) causes cranial, trunk, limb and respiratory weakess.1 Pathogenesis is multifaceted, and the gene-environment hypothesis proposes that ALS results from environmental exposures acting on a genetically susceptible background.1 We previously reported that polygenic risk contributes to the development of ALS2 and environmental risk scores (ERS), representing persistent organic pollutants (POPs) mixtures, associated with ALS risk and survival.3–5 Thus, both genetic and environmental factors contribute to ALS in our Michigan ALS cohort.

Several studies using diverse biosamples, such as blood,6 cerebrospinal fluid (CSF)7 8 and teeth,9 are linked to metal exposures to ALS.10 11 While significant associations between metal exposures and ALS vary across studies, selenium, lead, copper, aluminium, cadmium and zinc have been identified as significant risk factors for ALS in multiple studies, enhancing the statistical evidence behind these associations.11 However, these studies are limited by relatively small sample sizes,10 restricting the generalisability of findings and the ability to detect subtle associations. While recent studies suggest the potential for gene-metal interactions to influence ALS outcomes,12 13 no studies have yet to comprehensively examine this interaction.

To address these gaps and build on existing research, we conducted a large cohort study to investigate the association between metal exposures and ALS, as well as the potential influence of gene-metal interactions on ALS. We investigated the impact of metal levels in plasma and urine, individually and as a mixture modelled by ERS, on ALS risk and survival. This approach captured a comprehensive picture of metal exposures from diverse sources, as well as the biotransformation within and elimination processes from the body.14 We also examined the relationship between genetic background and metal exposures in ALS risk and determined if measured metal mixtures are associated with self-reported occupational exposures. These findings provide valuable insights into the role of metals in underlying disease mechanisms and identify potentially modifiable ALS risk factors.

Methods

Participants

Participant recruitment details are previously published.2–5 11 15–19 Briefly, all ALS patients meeting the Gold Coast criteria seen in the University of Michigan Pranger ALS clinic are approached for enrollment and depositing of blood and urine samples into the University of Michigan ALS Biorepository, provided they are able to consent in English. Control participants are identified via population outreach using a University of Michigan research database, random address mailing and Facebook advertisements. These population outreach efforts target individuals falling with the same sex and age distribution of the ALS participants. All participants provided informed consent for this University of Michigan IRB approved the study (HUM28826) and demographic information, including age, race, education, military service history, medical history and supplements/vitamins usage. Blood and/or urine samples collected at the time of recruitment were stored in the Michigan ALS Biorepository. Plasma samples were collected between 26 February 2012 and 7 December 2022, and urine samples were collected between 23 September 2015 and 9 September 022. Additional information abstracted from medical records included ALS disease characteristics, including onset segment, age of symptom onset, age of diagnosis, ALS family history, El Escorial Criteria and revised ALS functional rating scale (ALSFRS-R) score. The El Escorial Criteria and ALSFRS-R obtained on initial presentation to the Pranger ALS clinic are used for analysis. C9orf72 expansion status was obtained from medical records and previously published data.2

Biosamples

Detailed methods are reported in online supplemental file 1. Briefly, plasma and urine samples were analysed by inductively coupled mass spectrometry (ICP-MS) at the Dartmouth Trace Element Analysis Core. The Dartmouth Trace Element Analysis Core did not provide normal reference ranges.

Supplemental material

Polygenic scores

Two polygenic risk scores (PGS) were used. The first, an ALS-PGS, was constructed as previously reported2 but updated with weights from a recent large-scale ALS genome-wide association study (GWAS, (online supplemental table 1).20 PRSice-2 software was used for clumping and pruning single-nucleotide polymorphisms (SNPs) using default settings, followed by weight application from the parent GWAS. No cut-off for SNPs (p value threshold of 1.0) was used. Notably, the ALS-PGS was exclusively constructed for participants not included in the large ALS GWAS. A metal PGS was constructed to represent metal metabolism (metal-PGS), using relevant SNPs from literature related to a higher body burden of metals or poorer outcomes from xenobiotic exposures (online supplemental table 2). As the selected SNPs originated from various studies with diverse outcomes, weights were not applied based on effect estimates; instead, the metal-PGS was constructed based on the count of risk alleles.

Supplemental material

Supplemental material

Genetic samples

Online supplemental file 1 shows genotyping analysis of sample batches. QC filtering strategies at the SNP level are detailed in online supplemental figure 1, with sample-level filtering strategies illustrated in online supplemental figure 2.

Supplemental material

Missing data and multiple imputation

Measurements of metals in plasma and urine are subject to missingness due to detection limits. Five metals detected below the limit of detection (LOD) in over 40% of samples were excluded from ERS construction. Certain metals were manually excluded from ERS construction because the biofluid source was unsuitable for accurate and representative measurement.14 For plasma ERS, arsenic,14 cadmium14 and silver were excluded.21 For urine ERS, manganese,14 silver21 and strontium were excluded.22 Metal concentrations below LOD were imputed by Embedded Image .

Certain adjustment covariates exhibited missing data. In the risk model, these covariates included military service (n=50 for plasma and n=37 for urine). In survival analyses, these covariates included family history of ALS (n=12 for plasma and n=8 for urine), onset segment (n=2 for both plasma and urine), El Escorial criteria (n=17 for plasma, and n=20 for urine), ALSFRS-R (n=2 for plasma and n=1 for urine) and time between symptom onset and diagnosis (n=5 for plasma and n=4 for urine). Multiple imputation with chained equations was used for missing adjustment covariates and stratified by case-control status for each sample source. The control imputation model included metals (log-transformed), age at sample collection, sex and military service. The case imputation model additionally included age at diagnosis, onset segment, ALS family history, El Escorial criteria, ALSFRS-R, smoking status, education and cumulative hazard rate (Nelson-Aalen estimator) from the diagnosis to the last follow-up. Twenty multiple imputation datasets were generated using a predictive mean matching model.

Risk

For ALS risk models, single metal logistic regression models adjusted for age at sample collection, sex and military service were constructed. Models were fit independently across 20 imputation datasets, and results across models were pooled using Rubin’s rule. P values were adjusted for multiple comparisons by Benjamini-Hochberg correction. To assess the potential confounding effects of education and smoking, sensitivity analyses were conducted by including education or smoking as additional covariates, respectively.

ERS were constructed to assess multiple metals concurrently. As metal levels differ between plasma and urine samples, separate ERS scores for plasma (ERS PR; plasma risk ERS) and urine (ERSUR; urine risk ERS) were developed. An adaptive elastic net approach applied to the combined 20 imputed datasets using the R-package miselect 23 determined ERS weights. To optimise performance of the adaptive elastic net, hyperparameters were carefully tuned through a five-fold cross-validation process to minimise deviance based on the log likelihood for a logistic regression model.

To test for effect modification between genetic predisposition to ALS and elevated metal levels, logistic regression models evaluated the interaction between plasma/urine levels and ALS-PGS. To test for the effect modification between altered metal metabolism genetics and elevated metal levels, logistic regression models evaluated the interaction between plasma/urine levels and metal-PGS. Adjustment covariates for both interaction models included age at sample collection, sex, military service and the first five genetic principal components.2 Of note, the risk model did not adjust for C9orf72 gene positivity status as there were no C9orf72-positive control cases.2

Survival

Cox proportional hazards models adjusted for age at diagnosis, sex, family history of ALS, bulbar versus non-bulbar onset, diagnostic El Escorial criteria, ALSFRS-R and time between symptom onset and diagnosis were fit for each metal model post-diagnosis survival with subsequent Benjamini-Hochberg correction. Sensitivity analyses, in which smoking was included as an additional covariate, were conducted to assess the potential confounding effect of smoking on the results. To ensure that the definition of survival time did not impact the results, a sensitivity analysis was conducted by defining survival time from the date of sample collection to the event of interest.

Survival ERS was constructed via the application of a Cox regression with ridge penalty to the combined 20 imputation datasets, following published methods.3 4 Again, separate ERS were developed for plasma (ERSPS; plasma survival ERS) and urine (ERSUS; urine survival ERS). Hyperparameters of the ridge penalty were tuned through a five-fold cross-validation process with the aim of minimising the deviance based on the log partial likelihood for a Cox model. Kaplan-Meier survival curves depicted the survival functions of quartiles based on ERS. Survival curves were adjusted for covariates via the inverse probability weight method. Propensity scores were computed using multinomial logistic regression models, with ERS quartiles as outcome variables and adjustment covariates (as stated above) as predictors.

Cox survival models assessed interactions between genetic factors (ALS-/metal-PGS) and plasma/urine metal levels. Models were adjusted for age at diagnosis, sex, family history of ALS, bulbar versus non-bulbar onset, initial El-Escorial criteria, ALSFRS-R, time between symptom onset and diagnosis, the first five genetic principal components and C9orf72 gene positivity.

Occupational and non-occupational exposures

Occupational11 and non-occupational24 data were generated as previously described. Briefly, participants provided a list of prior job titles and responsibilities that were sorted into standard occupational classification (SOC) codes (online supplemental table 3) and answered questions regarding known metal exposures in their most recent occupation and non-occupational settings. SOC codes were separated into high- and low-metal exposure groups based on previously published results.11 Distributions of ERSPR and ERSUR by high-/low-metal exposure groups were compared using Wilcoxon tests. Associations between occupational/non-occupational metal exposures and urine/plasma metals were assessed using the partial least squares method feature in the R package ‘mixOmics’25 and estimated associations were visualised in heatmap plots.

Associations between metal concentrations and metal-PGS

An analysis was conducted to assess the associations between the metal-PGS and log concentrations of metals in plasma and urine. The associations were examined using linear regression models adjusted for sex and two genetic ancestry principal components.

Results

Participants and metal analysis

The full population cohort included 454 ALS and 295 control participants, with plasma samples from 387 ALS and 281 control participants (table 1, figure 1) and urine samples from 322 ALS and 194 control participants (table 2, figure 1). In both the plasma (table 1) and urine (table 2) cohorts, ALS individuals were older, more frequently male and had lower educational attainment than control participants. ALS characteristics were consistent between plasma and urine cohorts with 93% and 91% of participants presenting with classical ALS, 31% and 27% meeting El Escorial Definite ALS criteria and 26% and 27% having bulbar onset, respectively (tables 1 and 2). Both plasma and urine samples were available for 255 ALS and 181 control participants (online supplemental table 4).

Figure 1

Study overview. (A) Participants with (n=454) or without (n=294) amyotrophic lateral sclerosis (ALS) provided blood and/or urine samples. Environmental risk scores (ERS) of multiple metals associated with ALS risk and survival were calculated for both blood and urine samples. Polygenic risk scores (PGS) for ALS risk or genes related to metal metabolism were available for a subset of participants. (B) Blood and urine samples were collected from study participants and analysed by inductively coupled plasma mass spectrometry (ICP-MS) to measure the concentrations of various metals. After removing metals in which over 40% of the samples fell under the limit of detection (LOD), the statistical weight of each metal was determined, and overall ERS were calculated. Individual ERS for plasma and urine samples were developed for risk (plasma, ERSPR; urine, ERSUR) and survival (plasma, ERSPS; urine, ERSUS). These scores represent the cumulative exposure to various metals and serve as a comprehensive indicator of environmental risk. Image created using BioRender.com. (C) DNA samples from ALS and control participants were genotyped to obtain single-nucleotide polymorphisms (SNP) data. ALS-PGS was computed based on statistical weights from an independent genome-wide association study. The metal-PGS was constructed using unweighted SNPs linked to higher body burden of metals or poorer outcomes from xenobiotic exposures based on a review of multiple publications from the literature.

Table 1

ALS and control participant demographics: plasma

Table 2

Amyotrophic lateral sclerosis and control participant demographics: urine

Among the 42 familial ALS cases in our plasma sample cohort, 9 patients had a C9orf72 mutation, 7 had other mutations, 15 had no mutations identified and 11 had no available test results (online supplemental table 5A). In the 30 familial ALS cases in our urine sample cohort, 7 patients had a C9orf72 mutation, 5 had other mutations, 12 had no mutations identified and 6 had no available test results (online supplemental table 5B).

Concentrations of the metals detected in plasma and urine samples are listed in online supplemental tables 6 and 7, respectively. Metals found below the LOD in over 40% of plasma samples included aluminium, beryllium, mercury, molybdenum and uranium, while in urine samples, these were beryllium, chromium, manganese, mercury and silver. Among plasma metals, levels of cadmium, copper, lead, selenium, vanadium and zinc were significantly (padjusted <0.05) higher in ALS versus control participants, while antimony, arsenic, barium, and strontium levels were lower (online supplemental table 6). In urine, concentrations of aluminium, barium, cadmium, copper, iron, molybdenum, selenium, strontium, tin, uranium, vanadium and zinc were significantly higher in ALS versus control participants (online supplemental table 7). Correlations among metals were more pronounced in urine samples than plasma, and few metals showed significant interclass correlation between plasma and urine measures (online supplemental figure 3).

ALS risk

Plasma

ORs for metals ranged from 0.58 for antimony to 1.58 for selenium (figure 2). Metals significantly associated with elevated ALS risk included copper, lead, selenium and zinc. Conversely, antimony, barium, nickel and strontium were linked with reduced ALS risk. Metals included in the ERSPR, to represent the metal mixture, are detailed in online supplemental table 8. Analysis revealed a significant association between the ERSPR and ALS risk (online supplemental figure 4A), with a one SD increase in ERSPR based on the control population corresponding to an OR=2.95 (p<0.001, figure 2).

Figure 2

Plasma single metal and mixture associations with amyotrophic lateral sclerosis (ALS) risk. Single metal logistic regression models based on plasma samples from ALS (n=387) and control (n=281). The outcome is case/control status, the predictors are log-transformed, standardised metal levels, and the covariates are continuous age at sample collection, sex and military service. Mixture metal ERSPR (plasma risk environmental risk score). %BDL, percentage of samples below detection limit; ALS, amyotrophic lateral sclerosis; BH, Benjamini-Hochberg, red font are significant correlations; OR, odds ratio corresponding to one SD increase in log-transformed metals.

The model was then adjusted for ALS-PGS to account for the cumulative burden of genetic predisposition and environmental exposure. Following adjustment for ALS-PGS (online supplemental figure 5; online supplemental table 9A), copper, selenium and zinc associated with increased ALS risk, while antimony and strontium were linked to decreased risk. Notably, the ERSPR for metal mixtures remained significantly elevated with OR=2.71 (p<0.0001). Next, models with metals and PGS interaction terms explored the interplay between ALS polygenetic risk, metal exposure and disease susceptibility. Only chromium exhibited a statistically significant interaction with the ALS-PGS (p<0.01).

We next incorporated a PGS constructed with genes related to metal metabolism (metal-PGS) into our risk model to account for metal metabolism variability. Adjustment for metal-PGS reaffirmed associations with copper, zinc and selenium to increased ALS risk, and antimony, nickel and strontium with decreased ALS risk (online supplemental figure 6 online supplemental table 9B). Again, the ERSPR remained significantly elevated with an OR=2.93 (p<0.0001). Copper and molybdenum showed significant associations with the metal-PGS.

Urine

Analysis of urine samples revealed associations between barium, cadmium, copper, molybdenum, selenium, tin, uranium, vanadium and zinc with increased ALS risk (figure 3); no metals were associated with decreased risk. Consistent with plasma models, the ERSUR exhibited a significant association with ALS risk (OR=3.10, p<0.001; online supplemental table 8; online supplemental figure 4B). These findings remained consistent after adjusting for ALS-PGS (online supplemental figure 7; online supplemental table 9C) and metal-PGS (online supplemental figure 8; online supplemental table 9D). Significant metal and gene interactions were observed between molybdenum and ALS-PGS, as well as arsenic, barium, nickel and vanadium with metal-PGS.

Figure 3

Urine single metal and mixture associations with amyotrophic lateral sclerosis (ALS) risk. Single metal logistic regression models based on urine samples from ALS (n=322) and controls (n=194). The outcome is case/control status, the predictors are log-transformed, standardised metal levels, and the covariates are continuous age at sample, sex and military service. Mixture metal ERSUR (urine risk environmental risk score). %BDL, the percentage of samples below detection limit; ALS, amyotrophic lateral sclerosis; BH, Benjamini-Hochberg, red font are significant correlations; ORcorresponding to one SD increase in log-transformed metals.

ALS survival

Plasma

No metals exhibited a significant impact on ALS survival (figure 4). As expected, the ERSPS representing the cumulative effects of multiple metals (online supplemental table 8) associated with poorer survival (HR=1.37, p<0.001; online supplemental figure 9A). ALS participants were stratified into quartiles according to their ERSPS, revealing a striking difference in survival times between quartile 4 (1.59 years; highest exposure) and quartile 1 (3.75 years; lowest exposure; figure 5A). On adjustment for ALS-PGS, cadmium, manganese, lead and tin negatively impacted ALS survival, with interactions detected between ALS-PGS and molybdenum, selenium and strontium (online supplemental figure 10). Following adjustment for metal-PGS, no significant associations were found between metal main effects and ALS survival. However, significant interactions were identified between metal-PGS and iron, manganese, thallium and the ERSPS (online supplemental figure 11).

Figure 4

Cohort survival analysis based on plasma samples. Cox proportional hazards survival models from plasma samples from amyotrophic lateral sclerosis (ALS) cases (n=387). The outcome is the survival time since diagnosis (in years), the predictors are log-transformed, standardised metal levels, and adjustment covariates are age at diagnosis, sex, family history of ALS, onset segment, El Escorial criteria, ALSFRS-R and time between symptom onset and diagnosis. %BDL, the percentage of samples below detection limit; BH, Benjamini-Hochberg, red font are significant correlations; CI, confidence interval; ERSPS, plasma survival environmental risk score; HR corresponding to one SD increase in log-transformed metals.

Figure 5

Adjusted survival curves stratified by environmental risk score (ERS) quartile. Kaplan-Meier survival curves of ERS quartiles adjusted for covariates with the inverse probability weights method for (A) plasma (ERSPS) and (B) urine (ERSUS) sample sets. Adjusted covariates are age at sample collection, sex, military service and family history of ALS. Dashed lines indicate the median survival in each ERS strata. The estimated adjusted median survival times for plasma (A) are 3.75 years for Quartile 1 (Q1), 2.48 years for Quartile 2 (Q2), 1.99 years for Quartile 3 (Q3) and 1.59 years for Quartile 4 (Q4). Estimated adjusted median survival times for urine (B) are 3.24 years for Q1, 2.48 years for Q2, 2.10 years for Q3 and 1.79 years for Q4. Adjusted survival curves and adjusted median survival time estimates are pooled across all 20 imputed datasets.

Urine

In contrast to plasma findings, multiple urine metals were associated with shorter ALS survival, including cadmium, copper, nickel, tin and zinc (figure 6). However, consistent with plasma results, the ERSUS (8online supplemental table 8) correlated with shorter ALS survival (HR=1.44, p<0.001) (online supplemental figure 9B). Survival curves based on ERSUS are depicted in figure 5B, with estimated median survival times of 3.24 years in quartile 1 and 1.79 years in quartile 4, translating to a 1.45-year survival difference. After adjusting for ALS-PGS, copper, nickel, tin and zinc continued to correlate with poorer survival, while cadmium did not (online supplemental figure 12). Significant interactions were detected between antimony, barium, beryllium, molybdenum, nickel, selenium, uranium, zinc and the ERSUS and the ALS-PGS (p<0.1). No associations were found between metals and ALS survival following adjustment for metal-PGS, although the ERSUS remained correlated with shorter survival (online supplemental figure 13). Interactions between metals and metal-PGS were observed for beryllium, copper, thallium, tin, vanadium, zinc and the ERSUS.

Figure 6

Cohort survival analysis based on urine samples. Cox proportional hazards survival models from urine samples from amyotrophic lateral sclerosis (ALS) cases (n=322). The outcome is survival time since diagnosis (in years), the predictors are log-transformed, standardised metal levels, and adjustment covariates are age at diagnosis, sex, family history of amyotrophic lateral sclerosis (ALS), onset segment, El Escorial criteria, ALSFRS-R and time between symptom onset and diagnosis. %BDL, the percentage of samples below detection limit; BH, Benjamini-Hochberg, red font are significant correlations; ERSUS, urine survival environmental risk score; HR corresponding to one SD increase in log-transformed metals.

Metal concentrations and metal-PGS

The associations between metal concentrations and metal-PGS are summarised in online supplemental table 10. Notably, an association was observed with beryllium in plasma, as well as inverse associations between aluminium and mercury in urine, with metal-PGS scores (p<0.05).

Occupational and non-occupational metal exposures

Plasma

ERSPR scores for individuals grouped into the high-metal exposure group (job titles of ‘Building and Grounds Cleaning and Maintenance’, ‘Construction and Extraction’, ‘Installation, Maintenance and Repair’ and ‘Production Occupations’) were significantly higher than ERSPR scores for those grouped into the low-metal exposure group (online supplemental figure 14A). Furthermore, occupational metal exposures were positively associated with antimony, iron and thallium levels (online supplemental figure 15). Finally, reported exposures to metals outside the occupational setting were positively associated with increased aluminium, cadmium, chromium, lead, manganese, uranium, vanadium and zinc levels (online supplemental figure 16).

Urine

Again, individuals in the high-metal exposure group based on SOC codes had higher ERSUR scores than those in the low-metal exposure group (online supplemental figure 14B). Self-reported occupational metal exposures positively associated with mercury and thallium levels (online supplemental figure 17). Additionally, self-reported metal exposure in non-occupational activities demonstrated positive associations with arsenic, cadmium, copper, iron, lead, mercury, molybdenum, nickel, selenium, strontium, thallium, tin, uranium and zinc (online supplemental figure 18).

Supplement/vitamin usage

Recognising the potential impact of supplement usage on metal levels, we compared metal levels characterised by participants’ supplement/vitamin usage. In plasma samples, significant differences in median metal concentrations between individuals using supplements/vitamins and those abstaining from them were only detected for selenium and vanadium (online supplemental table 11). No significant differences were observed for urine measures (online supplemental table 12). Overall, supplement/vitamin usage does not appear to substantially impact outcomes.

Sensitivity analyses

To address potential residual confounding related to education, a sensitivity analysis was conducted by including education as an additional covariate in the logistic regression models. The results were consistent with our original findings in both plasma samples (online supplemental table 13) and urine samples (online supplemental table 14). This suggests that the observed associations are robust and not significantly influenced by differences in educational levels.

Additionally, to address the potential impact of smoking behaviour on our findings, we included smoking status (Former Smoker, Smoker, Non-Smoker) as an additional covariate in both logistic regression and survival analyses. The findings continued to align closely with the original analyses for both plasma and urine samples, as detailed in online supplemental table 15 (risk) and online supplemental table 16 (survival) for plasma, as well as online supplemental table 17 (risk) and online supplemental table 18 (survival) for urine. Notably, the sensitivity analyses demonstrated significant associations between elevated levels of aluminium and antimony in urine with increased ALS risk and higher selenium in urine with shorter ALS survival, which were not observed in the original analyses (padjusted <0.05). Despite these changes, the ORs and HRs for these metals remained consistent between the original and sensitivity analyses, indicating that the overall strength and direction of the associations did not significantly change.

A sensitivity analysis using survival time defined as the time from the date of sample collection was conducted to address potential bias from the time gap between ALS diagnosis and sample collection. The results were consistent with our original findings in both plasma (online supplemental table 19) and urine (online supplemental table 20) samples, reinforcing the robustness of our initial analyses.

Discussion

ALS is influenced by genetic and environmental factors, including metals.10 By measuring plasma and urine metal levels, we confirmed that individual metals, such as copper, selenium and zinc, were associated with greater ALS risk and shorter survival (figure 7). Notably, the most robust associations with both risk and survival were observed with cumulative ERS, representing a mixture of multiple metals. Despite its influence on ALS susceptibility, genetic background did not significantly alter the relationship between metal exposure and ALS risk or survival. Finally, engagement in occupations with heightened metal exposure was associated with increased ERS. Overall, our report supports previous research and underscores the pivotal role of metals in ALS pathophysiology within this Michigan-based cohort.

Figure 7

Study summary. (A) Associations of individual metals and environmental risk scores (ERS) with amyotrophic lateral sclerosis (ALS) risk and survival. Variables with ORs of >1 or HRs of >1 are indicated as red upward arrows, while those with OR of <1 or HR of <1 are indicated as green downward arrows. (B) Associations between ERSs and occupational metal exposures. Occupations and occupational metal exposures are reported by questionnaire. Occupational groups are defined by Standard Occupational Classification (SOC) codes. Four occupational groups identified in the literature as exhibiting high metal exposure levels: ‘Building and Grounds Cleaning and Maintenance’ (SOC 37), ‘Construction and Extraction’ (SOC 47), ‘Installation, Maintenance and Repair’ (SOC 49) and ‘Production Occupations’ (SOC 51) had higher ERSPR and ERSUR. (C) We did not observe evidence of interaction between polygenic risk scores and ERS measures on ALS risk and survival models.

The associations between copper, selenium and zinc and ALS are consistent with previous clinical studies.6–9 26 Exposure to these metals can occur through diet, drinking water, supplements and environmental pollutants. In preclinical models, the accumulation of these metals induces oxidative stress, disrupts mitochondrial function and impairs protein folding,8 27 28 all mechanisms implicated in ALS.1 For selenium, excess levels promote SOD1 accumulation in the mitochondria of human neuroblastoma cells,8 a feature also observed in neurons during ALS pathogenesis.8 In regards to zinc, the superoxide dismutase 1 (SOD1) enzyme protects cells from oxidative stress, and zinc ions are essential for SOD1 stability and enzymatic activity. Mutations in SOD1 are associated with ALS,1 and alterations in zinc levels cause SOD1 misfolding and aggregation.29 Indeed, chronic oral administration of zinc sulfate to transgenic mice overexpressing the human-mutated form of SOD1 (G93A) decreases survival.30 While the exact mechanisms linking these metals to ALS are not fully understood, accumulating evidence suggests that their altered levels contribute to ALS pathophysiology.

In the current study, antimony, nickel and strontium are associated with decreased ALS risk. Research investigating the effects of antimony on human health are limited. However, antimony is associated with neurobehavioural changes and neurotoxicity in preclinical models.31 Regarding nickel, a population-based case-control study in Denmark did not detect an association between occupational exposure to nickel and ALS.32 It is surprising that nickel would be protective against ALS, considering it influences oxidative stress, increases DNA damage and induces apoptosis.32 Similarly, the protective role of strontium contradicts with our previous study, which analysed human permanent teeth and found that strontium during late childhood and adolescence is correlated with ALS.9 This discrepancy could stem from the biosample analysed, as strontium uptake in teeth occurs during tooth formation, reflecting exposure levels during that specific developmental stage.9 Thus, the potential protective effects of antimony, nickel and strontium against ALS warrant further confirmation and investigation.

The use of plasma and urine samples to measure metal exposures offers several advantages, for example, urine sampling is relatively easy and non-invasive, and both reflect integrated exposure. However, measurements must be interpreted cautiously. Levels in plasma and urine may fluctuate due to diet, hydration status and circadian rhythms, potentially leading to variability in measurements14; though, it is worth noting that supplement or vitamin use did not impact metal levels in the current study. Additionally, blood and urine often represent recent, rather than chronic, exposures. Moreover, some metals may not be well-represented in plasma or urine. Zinc is distributed intracellularly in the musculoskeletal system, and plasma zinc does not represent intracellular accumulation.14 Urinary zinc is an accurate biomarker but again represents recent exposures.14 Biomarkers such as hair, teeth and nails can provide a historical record of metal accumulation over time and have been used to study chronic metal exposure and ALS.9 33 34 However, challenges remain regarding biomarker suitability, variability in sample quality and growth rate.14 Thus, selecting the appropriate biomarker is crucial for accurately detecting metal levels and assessing exposure.

Most previous epidemiological investigations involving metals and ALS associations centred on individual metals as opposed to mixtures. However, a mixture approach using the ERS better reflects reality. Here, we found that metal mixtures, represented by ERS, had a greater effect on ALS risk and survival than single metal exposures. ERS derived from the quantification of metal levels in biofluids have been used to explore the relationship between metal exposure and disease across various health domains, including obesity and its comorbidities,35 glucose homeostasis36 and heart rate regulation.37 Thus, ERS are a valuable tool for studying multi-pollutant exposures and their impact on disease.

The gene-environment hypothesis suggests that interactions between genetic factors and environmental exposures influence an individual’s risk for developing a particular disease.1 Several studies suggest that environmental factors, such as metals, interact with genetic variants to influence ALS onset, progression and severity.12 13 Additionally, previous studies successfully integrated environmental and genetic factors into risk prediction models for cancer.38 39 However, we found that accounting for ALS risk genes or genes involved in metal metabolism did not influence the relationship between metal exposure and ALS. Relationships between genes, environmental exposures and disease risks are complex, and adjusting for underlying genetic heterogeneity may not fully capture the complicated nature of these interactions. Additionally, our current knowledge of the genetic basis of ALS and metal metabolism may be incomplete and other genetic factors or pathways involved in ALS risk or metal metabolism may be unaccounted for in our current PGS. Construction of our metal-PGS was limited by the lack of available relevant GWAS of metal levels and is, in turn, a weak genetic instrument for assessing any potential interaction with ALS-PGS. Its association with only a few metals further confirms the need for a more robust future metal-PGS before drawing meaningful conclusions. We also do not know how well the metal-PGS reflects metal metabolism. Factors such as efficiency of elimination and excretion of metals from the body may explain inverse associations observed with urine metals. Future studies require a newly constructed metal-PGS in the setting of well-characterised exposures and samples measured in a relevant period to determine true interaction with metal exposures.

Finally, we previously showed significant associations between workplace metal exposure and ALS risk and identified SOC codes that were likely to report occupational exposure to metals.11 Here, ALS participants categorised into a high-metal exposure group exhibited higher ERSs compared with those in the low-metal exposure groups. These results support the finding that occupations linked to metal exposures associate with elevated ALS risk. Furthermore, occupational and non-occupational metal exposures24 correlated with metal levels in plasma and urine. This association underscores the potential impact of both activities on overall metal exposure in individuals, emphasising the necessity of accounting for occupational and environmental sources when evaluating overall exposure risk.

Recognising these exposures as modifiable risk factors may open avenues for targeted prevention. By avoiding high-risk activities associated with occupational and non-occupational metal exposures, individuals might lower their overall exposure levels, potentially mitigating ALS risk. However, the strength of these associations varied, with some being weak or even negative. This variability highlights the complex nature of metal exposures and their measurements. Given this variability, further research is needed to better understand these associations and their implications for ALS pathogenesis.

This study benefits from a large sample size and biomarker-based exposure measurements, reducing recall bias. However, limitations exist. First, the study represents a cross-sectional analysis and may miss potential fluctuations in metal levels over time. Consequently, this snapshot may not accurately reflect lifetime exposures as metal levels are subject to sampling variation and metabolic pathways, potentially altering their concentrations over extended periods. Second, sampling after disease onset may confound results, blurring the distinction between disease progression effects and true baseline measures and future longitudinal prospective studies will be crucial. Third, while our findings suggest an association between metals and ALS, they do not establish causation. Fourth, unlike the ALS-PGS, the metal-PGS does not have a large GWAS from which to draw weights and its construction is more ad hoc and less robust than standard methods used for the ALS-PGS. Large genetic studies of xenobiotic metabolism would be critical to further understand gene-environment interactions. Fifth, there is a potential limitation related to the time gap between diagnosis and sample collection. This variability introduces a form of left truncation, as patients must survive until the time of sample collection to be included in the study, potentially biasing our findings by excluding patients with shorter survival times who may have different exposure profiles. Future studies with more uniform timing of sample collection relative to diagnosis are necessary to validate our findings and mitigate this potential source of bias. Our findings and mitigate this potential source of bias. Sixth, although we adjusted for age, sex and military service, other potential confounders, such as socioeconomic status and lifestyle factors, may not be fully accounted for. Sensitivity analyses indicated that education and smoking did not significantly impact our findings, but we recognise the possibility of unidentified confounding factors that could influence the results. Lastly, given the study used a Michigan-based cohort, validation across different cohorts is critical for broader generalisability.

To conclude, we demonstrate that metal mixtures, as represented by cumulative ERS, associate with increased ALS risk and reduced ALS survival, regardless of genetic predisposition. These associations correlate with self-reported occupational metal exposure. Using quantitative assessments of metal levels in biofluids like plasma or urine could offer a pathway to predict future ALS risk, providing a promising avenue for early identification and intervention. Overall, these results offer valuable insights into the complexities of ALS pathogenesis.

Data availability statement

Data are available upon reasonable request.

Ethics statements

Patient consent for publication

Ethics approval

This study involved human participants and was approved by the IRB at the University of Michigan (HUM28826). Participants gave informed consent to participate in the study before taking part.

Acknowledgments

We are indebted to the participants who provided samples. We thank Crystal Pacut, Stacey Jacoby, PhD, Madeleine Batra, Hasan Farid, MS, Caroline Piecuch, Sushant Obeja and Adam Patterson for study support. We thank Masha Savelieff, PhD for expert editorial support with figures.

References

Supplementary materials

Footnotes

  • Contributors D-GJ, JFD, LZ, KMB, BM, SB, ELF and SAG conceived and designed the study. ST and SAG contributed to the acquisition of the data. D-GJ, JD and KMB performed the statistical analyses. D-GJ, EK and SAG interpreted the data, drafted the text and prepared the figures. All authors critically reviewed and approved the final version of the manuscript. SAG is the guarantor of the study.

  • Funding Funding was provided by the National Institute of Neurological Disorders and Stroke (NINDS) (R01NS127188); National Institute of Environmental Health Sciences (NIEHS) (K23ES027221, R01ES030049); Centers for Disease Control and Prevention (R01TS000344); ALS Association (20-IIA-532, 20-PP-661); NeuroNetwork for Emerging Therapies; Peter R. Clark Fund for ALS Research; Robert and Katherine Jacobs Environmental Health Initiative; Richard Stravitz Foundation; Coleman Therapeutic Discovery Fund; Scott L. Pranger ALS Clinic Fund; the Dr. Randall W. Whitcomb Fund for ALS Genetics; University of Michigan. Metals analysis was carried out at the Dartmouth Trace Element Core Facility, which is supported by Dartmouth Cancer Center with NCI Cancer Center Support Grant 5P30 CA023108.

  • Competing interests ELF: Listed as inventors on a patent, Issue number US10660895, held by University of Michigan titled 'Methods for Treating Amyotrophic Lateral Sclerosis' that targets immune pathways for use in ALS therapeutics. SAG: Listed as inventors on a patent, Issue number US10660895, held by University of Michigan titled 'Methods for Treating Amyotrophic Lateral Sclerosis' that targets immune pathways for use in ALS therapeutics. Scientific consulting for Evidera.

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

  • Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.