Skip to main content

Investigation of the causal relationship between ALS and autoimmune disorders: a Mendelian randomization study



Epidemiological studies have reported an association between amyotrophic lateral sclerosis (ALS) and different autoimmune disorders. This study aims to explore the causal relationship between autoimmune disorders and ALS using Mendelian randomization (MR).


To test the genetically predicted effects of liability towards immune-related outcomes on ALS risk, we used summary statistics from the largest European genome-wide association studies (GWAS) for these disorders in a two-sample MR setting. To do this, we extracted single nucleotide polymorphisms (SNPs) from the GWAS, which strongly associated with the 12 traits, and queried their effects in a large European ALS GWAS (27,265 cases and 110,881 controls). To avoid bias in our MR instruments related to the complex linkage disequilibrium structure of the human leukocyte antigen (HLA) region, we excluded SNPs within this region from the analyses. We computed inverse-variance weighted (IVW) MR estimates and undertook sensitivity analyses using MR methods robust to horizontal pleiotropy. We also performed a reverse MR analysis testing the causal effects of ALS on the above autoimmune traits.


After applying Bonferroni correction for multiple testing, our MR analyses showed that the liability to autoimmune disorders does not affect ALS risk. Our reverse MR analysis also did not support the effects of liability to ALS on other autoimmune disorders. The results of the main IVW MR analyses were generally supported by our sensitivity MR analyses. The variance in the exposures explained by the sets of SNPs used as MR instruments ranged from 8.1 × 10−4 to 0.31. Our MR study was well-powered to detect effects as small as an odds ratio (OR) of 1.045 for ALS in the main MR and as small as an OR of 1.32 in the reverse MR.


Our MR study does not support a relationship between liability to autoimmune disorders and ALS risk in the European population. The associations observed in epidemiological studies could be partly attributed to shared biology or environmental confounders.

Peer Review reports


Amyotrophic lateral sclerosis (ALS) is a rapidly progressive neurodegenerative disorder that leads to the paralysis of almost all skeletal muscles, with a high mortality rate typically due to respiratory paralysis in the 3–5 years following diagnosis. The lifetime risk of ALS is estimated to be 1 in 400 [1, 2], and less than 10% of patients survive beyond 10 years [3, 4]. ALS is caused by genetic and environmental factors, with both familial and sporadic cases described. The genetic and environmental interplay in ALS varies depending on the innate penetrance of pathogenic variants, polygenic predisposition, and neurotoxicity of specific environmental modifiers [5,6,7]. Consequently, the identification of underlying pathologies that could trigger ALS might point to avenues for new treatments.

Neuroinflammation and autoimmunity have been reported to be associated with ALS in previous studies [8,9,10,11,12,13,14,15,16]. The presence of autoimmune disorders has been reported to increase the risk of ALS [17] which could suggest a shared genetic architecture between these diseases. Although epidemiological studies suggest associations between ALS and autoimmune disorders, it remains unknown if these associations are causal. Because of unmeasured confounding and reverse causation, these epidemiological studies are prone to bias hampering the potential for causal inference. Nonetheless, determining if there is a causal association between ALS and autoimmune disorders could point to specific biological pathways and inform prevention strategies.

Mendelian randomization (MR) is a study design that allows to explore the causality between an exposure of interest (here an autoimmune disorder) and an outcome (here ALS) using an instrumental variables (IVs) approach [18]. In MR, genetic variants strongly associated with an exposure, and satisfying specific assumptions are used as IVs to study the causal association with an outcome. Since these variants are randomly assigned at conception, this could reduce bias due to environmental confounders if MR is conducted properly. As such, the MR design is analogous to a randomized controlled trial. In the current study, we applied bidirectional MR to seek evidence of a causal association between autoimmune disorders (celiac disease (CeD), multiple sclerosis (MS), rheumatoid arthritis (RA), Crohn’s disease (CD), psoriasis (PsO), primary sclerosing cholangitis (PSC), asthma, irritable bowel syndrome (IBS), primary biliary cirrhosis (PBC), ulcerative colitis (UC), type 1 diabetes (T1D), systemic lupus erythematosus (SLE)) and ALS using summary statistics from the largest available genome-wide association studies (GWAS) in European populations for the above traits. The term “autoimmune” will be used to describe both autoimmune and autoinflammatory disorders (such as asthma).


In order to perform any MR study, the valid genetic variants (SNPs) used as instruments should satisfy three criteria: the relevance assumption—they should be strongly associated with the exposure of interest; the independence assumption—there is no shared common cause with the outcome; and the exclusion restriction assumption—SNPs only affect the outcome through the path of the exposure [19].

Study cohorts and GWAS

To perform our MR analyses, we used summary-level data from the largest publicly available GWAS for each trait (Table 1). Specifically, the summary statistics of the 12 autoimmune disorders GWAS are available from the GWAS catalog [20, 21] and the IEU OpenGWAS project [22] and the ALS GWAS from the Project MinE [23, 24]. Information on recruitment procedures and diagnostic criteria is detailed in the original publications. All cases and controls in these studies were of European ancestry. Also, there is no significant overlap between GWAS populations.

Table 1 Characteristics of the ALS and autoimmune disease GWAS cohorts and results of the power analysis

Mendelian randomization

IV selection

For the selection of genetic instruments from each of the eleven exposure GWASs, we used the default settings in the R package TwoSampleMR [36, 37]. Specifically, genome-wide significant (p-value < 5.0E−08) SNPs were extracted, and SNPs harboring the HLA region (chr6:27,477,797–34,448,354, hg19/GRCh37) were excluded [38]. For analyses in which there were insufficient IVs (due to SNPs not available in outcome GWAS), we used the LDproxyR tool to replace the IVs with proxy SNPs in high linkage disequilibrium (LD r2 > 0.8). Standard clumping parameters were used to select independent SNPs (clumping window of 10,000 kb, LD r2 cutoff 0.001). The proportion of variance of the exposures explained by the SNPs (R2) and F-statistics were calculated to estimate the strength of IVs [39] to satisfy the first MR assumption. The same approach was taken for the reverse MR; as such, we performed 24 bidirectional MR studies, where autoimmune disorders and ALS were regarded either as exposure or as the outcome.

Power calculation

Using an online Mendelian randomization power calculation, we set the power to 80% and reported the minimum odds ratio (OR) for each MR analysis [40, 41] (Table 1).

Mendelian randomization analyses

Bidirectional two-sample MR was employed in the TwoSample MR R package [36, 37]. After clumping IVs, we performed Steiger filtering to exclude SNPs explaining more variance in the outcome than in the exposure [36]. Then, we applied the inverse variance weighted (IVW) method to combine the effect of different IVs. The Wald ratio was calculated for each SNP, and the individual effect of each SNP was meta-analyzed using IVW to generate the concluding beta estimate, which was transformed to an OR [42,43,44]. To test the third MR assumption, we applied MR Egger to detect possible violations of instrumental variable assumptions due to directional horizontal pleiotropy [42]. Additionally, we used weighted median (WM) which is a median of the weighted estimates and provides a consistent effect even if 50% of IVs are pleiotropic [45]. Heterogeneity was tested using Cochran’s Q test in the IVW and MR-Egger methods [46]. To illustrate the results of the different MR methods, we constructed scatter plots using the TwoSampleMR package. Additionally, we used the MR-PRESSO test to detect outlier SNPs which may be biasing estimates through horizontal pleiotropy and adjust for these [47]. Finally, a leave-one-out analysis (LOO) was performed to detect if there is any single SNP disproportionately responsible for the result of each MR study.

Multivariable MR

UC and CD have been reported to have a genetic association [48]. To control for the pleiotropic effect in this study, we run multivariable MR [49] for the studied traits (UC, and CD). The significant SNPs were extracted, then combined and clumped. The analysis was performed using the TwoSampleMR R package.


Our IVW results showed that liability to autoimmune disorders does not affect ALS risk which is consistent with results from the other MR methods including MR Egger and weighted median (Table 2, Fig. 1). The results of the MR analyses investigating the causal relationship between ALS and 12 autoimmune disorders are shown in Table 2. The estimates (ORs) represent the effects on ALS risk of genetically predicted liability to each disease exposure. Our MR analyses had 80% power to detect the effect sizes of liability to exposures on ALS corresponding to OR from 1.045 to 2 in the main MR analyses and the effects of liability to ALS on the risk of autoimmune diseases corresponding to OR from 1.32 to 2.37 in the reverse MR analyses (Table 1). The variance in the exposures explained by their respective set of SNPs ranged from 0.08 to 31% (Tables 2 and 3). All instruments had an F-statistics of > 39, which is above the standard cutoff (> 10) indicating sufficient instrumental strength [50, 51] (Tables 2 and 3). Applying Bonferroni correction for multiple testing, a p-value below 2.1E−03 was considered as significant.

Table 2 Results of the MR analyses between liability to autoimmune disorders and the ALS risk
Fig. 1
figure 1

Scatter plots showing the MR effect of each exposure on ALS. The x-axis represents the genetic association with autoimmune disease risk; the y-axis represents the genetic association with the risk of ALS. Each line represents a different MR method. CD, Crohn’s disease; CeD, celiac disease; IBS, irritable bowel syndrome; MS, multiple sclerosis; PBC, primary biliary cirrhosis; PSC, primary sclerosing cholangitis; PsO, psoriasis; RA, rheumatoid arthritis; T1D, type 1 diabetes; UC, ulcerative colitis; SLE, systemic lupus erythematosus; MR, Mendelian randomization. No significant associations were detected

Table 3 Results of reverse MR analysis between liability to ALS and risk of autoimmune disorders

Sensitivity analyses were performed to detect the presence of horizontal pleiotropy. Significant heterogeneity was apparent in our IVs for RA (MR Egger, Q p-value = 6.8E−05; IVW, Q p-value = 3.2E−05) which is illustrated in our funnel plot and LOO plot (Additional file 1: Table S1, Additional file 1: Fig. S1). The MR-Egger intercept did not provide evidence for unbalanced horizontal pleiotropy as it was centered around zero for all MR analyses. However, MR-PRESSO identified outlier SNPs for RA (rs9275183), CD (rs42861), and CeD (rs13198474) (Additional file 1: Table S1, Additional file 1: Fig. S1). The distortion test did not suggest significant changes in the effect estimates after removing these outlier SNPs (Additional file 1: Table S1).

Additionally, performing reverse MR studies with liability to ALS as the exposure and risk of the 12 autoimmune diseases as outcomes, we found no significant association. In these studies, we used 2–10 ALS-associated SNPs including proxy SNPs as IVs, which explained 0.004–0.007 of the variances in ALS risk and had an F-statistic of more than 63. Significant heterogeneity was apparent in our IVs for ALS (when CeD is the outcome) (IVW, Q p-value = 4E−07). However, because of not enough instrumental variables for the MR of ALS and CeD, MR-PRESSO analyses and LOO plots were not possible to be provided. The results of all three MR methods were consistent (Table 3, Additional file 1: Fig. S2).

Using 68 and 44 SNPs robustly and independently associated with CD and UC, respectively, multivariable MR provided evidence that liability to these autoimmune traits does not affect ALS risk.


In this study, we performed a comprehensive bidirectional two-sample MR study to investigate the causal relationships between liability to various autoimmune disorders and ALS risk. Using this approach, we found no relationship between liability to these disorders and ALS, which implies that the reported epidemiological associations could be the result of unmeasured confounding or shared genetic architecture. To distinguish between a true negative result and a lack of validity of the MR studies, multiple sensitivity analyses were applied to ensure that the three MR assumptions were satisfied. Given the consistency of our MR findings across these different methods, we are confident about the validity of our MR analyses to exclude moderate to large causal effects of the exposures on the outcomes.

Our results oppose with an observational study by Turner et al. supporting an association between liability to autoimmune disorders and increased risk of ALS [17]. In another comprehensive study investigating the genetic correlation between ALS and 10 autoimmune disorders, Li et al. reported a positive genetic correlation between ALS and CeD, MS, RA, and SLE, which is in the same direction as those found in a study by Turner et al. [52]. One possible explanation for the described associations between autoimmune disorders and ALS in the absence of causal effect in our MR studies is pleiotropy. In the study by Li et al., using conjunctional FDR statistics, the authors identified shared genetic loci between ALS and autoimmune disorders; these loci encompassed membrane trafficking, vesicle-mediated transport, endoplasmic reticulum (ER) to Golgi anterograde transport, and transport to the Golgi and subsequent modification [52]. These results support the hypothesis that the pathology of ALS might be mediated by a dysfunction of the immune system. However, these findings do not imply that liability to these autoimmune disorders predisposes to ALS and cannot resolve the previous debate on whether dysregulated immunity is the cause of ALS or a consequence of it. Overall, these findings support a shared genetic architecture between ALS and autoimmune disorders, but it is likely that there is no causality between the autoimmune disorders and ALS based on our MR study.

Our study has several limitations. This is a European-based study, and our findings cannot be generalized to other populations. Additionally, some of our MR analyses did not have enough power to detect small effects, due to the limited variance of the exposures explained by the SNP instruments or the limited sample sizes of the outcome GWAS. In this direction, excluding ambiguous or palindromic SNPs from our MR instruments might have further affected the power of our MR studies. Larger GWAS on these autoimmune traits will boost the power of future MR studies to detect associations.


Our results do not support that liability to various autoimmune disorders affects ALS risk in Europeans and suggest that the observed associations could be a result of shared genetic effects or environmental confounders.

Availability of data and materials

All data generated and codes used in the current study are available in this published article and Additional file 1 associated with it.



Amyotrophic lateral sclerosis


Crohn’s disease


Celiac disease


Endoplasmic reticulum


Genome-wide association studies


Human leukocyte antigen


Irritable bowel syndrome


Instrumental variables


Inverse-variance weighted


Leave-one-out analyses


Mendelian randomization


Multiple sclerosis


Odds ratio


Primary biliary cirrhosis


Primary sclerosing cholangitis




Rheumatoid arthritis


Systemic lupus erythematosus


Single nucleotide polymorphisms


Type 1 diabetes


Ulcerative colitis


Weighted median


  1. Brown RH, Al-Chalabi A. Amyotrophic lateral sclerosis. N Engl J Med. 2017;377(2):162–72.

    Article  CAS  PubMed  Google Scholar 

  2. Taylor JP, Brown RH, Cleveland DW. Decoding ALS: from genes to mechanism. Nature. 2016;539(7628):197–206.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Chio A, Logroscino G, Hardiman O, Swingler R, Mitchell D, Beghi E, et al. Prognostic factors in ALS: a critical review. Amyotroph Lateral Scler. 2009;10(5-6):310–23.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Gladman M, Zinman L. The economic impact of amyotrophic lateral sclerosis: a systematic review. Expert Rev Pharmacoecon Outcomes Res. 2015;15(3):439–50.

    Article  PubMed  Google Scholar 

  5. Cooper-Knock J, Jenkins T, Shaw PJ. editors. Clinical and molecular aspects of motor neuron disease. Colloq Ser Genomic Mol Med. 2013;2(2):1–60.

    Article  Google Scholar 

  6. Ryan M, Heverin M, McLaughlin RL, Hardiman O. Lifetime risk and heritability of amyotrophic lateral sclerosis. JAMA Neurol. 2019;76(11):1367–74.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Trabjerg BB, Garton FC, van Rheenen W, Fang F, Henderson RD, Mortensen PB, et al. ALS in Danish Registries: Heritability and links to psychiatric and cardiovascular disorders. Neurol Genet. 2020;6(2).

  8. Zhang R, Hadlock KG, Do H, Yu S, Honrada R, Champion S, et al. Gene expression profiling in peripheral blood mononuclear cells from patients with sporadic amyotrophic lateral sclerosis (sALS). J Neuroimmunol. 2011;230(1-2):114–23.

    Article  CAS  PubMed  Google Scholar 

  9. Zhang R, Gascon R, Miller RG, Gelinas DF, Mass J, Hadlock K, et al. Evidence for systemic immune system alterations in sporadic amyotrophic lateral sclerosis (sALS). J Neuroimmunol. 2005;159(1-2):215–24.

    Article  CAS  PubMed  Google Scholar 

  10. Prinz M, Priller J. The role of peripheral immune cells in the CNS in steady state and disease. Nat Neurosci. 2017;20(2):136–44.

    Article  CAS  PubMed  Google Scholar 

  11. Murdock BJ, Zhou T, Kashlan SR, Little RJ, Goutman SA, Feldman EL. Correlation of peripheral immunity with rapid amyotrophic lateral sclerosis progression. JAMA Neurol. 2017;74(12):1446–54.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Mccombe PA, Henderson RD. The role of immune and inflammatory mechanisms in ALS. Curr Mol Med. 2011;11(3):246–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Kawamata T, Akiyama H, Yamada T, McGeer P. Immunologic reactions in amyotrophic lateral sclerosis brain and spinal cord tissue. Am J Clin Pathol. 1992;140(3):691.

    CAS  Google Scholar 

  14. Henkel JS, Engelhardt JI, Siklós L, Simpson EP, Kim SH, Pan T, et al. Presence of dendritic cells, MCP-1, and activated microglia/macrophages in amyotrophic lateral sclerosis spinal cord tissue. Ann Neurol. 2004;55(2):221–35.

    Article  CAS  PubMed  Google Scholar 

  15. Henkel JS, Beers DR, Wen S, Rivera AL, Toennis KM, Appel JE, et al. Regulatory T-lymphocytes mediate amyotrophic lateral sclerosis progression and survival. EMBO Mol Med. 2013;5(1):64–79.

    Article  CAS  PubMed  Google Scholar 

  16. Duarte F, Binet S, Lacomblez L, Bouche P, Preud’homme J-L, Meininger V. Quantitative analysis of monoclonal immunoglobulins in serum of patients with amyotrophic lateral sclerosis. J Neurol Sci. 1991;104(1):88–91.

    Article  CAS  PubMed  Google Scholar 

  17. Turner MR, Goldacre R, Ramagopalan S, Talbot K, Goldacre MJ. Autoimmune disease preceding amyotrophic lateral sclerosis: an epidemiologic study. Neurology. 2013;81(14):1222–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Burgess S, Small DS, Thompson SG. A review of instrumental variable estimators for Mendelian randomization. Stat Methodol. 2017;26(5):2333–55.

    Google Scholar 

  19. Davies NM, Holmes MV, Smith GD. Reading Mendelian randomisation studies: a guide, glossary, and checklist for clinicians. BMJ. 2018;362.

  20. Buniello A, MacArthur JAL, Cerezo M, Harris LW, Hayhurst J, Malangone C, et al. The NHGRI-EBI GWAS Catalog of published genome-wide association studies, targeted arrays and summary statistics 2019. Nucleic Acids Res. 2019;47(D1):D1005–D12.

    Article  CAS  PubMed  Google Scholar 

  21. GWAS Catalog. Accessed 20 March 2022.

  22. IEU OpenGWAS project. Accessed 20 March 2022.

  23. van Rheenen W, van der Spek R, Bakker M, van den Berg L, Veldink J, van Vugt J, et al. Common and rare variant association analyses in amyotrophic lateral sclerosis identify 15 risk loci with distinct genetic architectures and neuron-specific biology. Nat Genet. 2021;53(12):1636–48.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Project MinE. Accessed 20 March 2022.

  25. Valette K, Li Z, Bon-Baret V, Chignon A, Bérubé J-C, Eslami A, et al. Prioritization of candidate causal genes for asthma in susceptibility loci derived from UK Biobank. Commun Biol. 2021;4(1):1–15.

    Article  Google Scholar 

  26. De Lange KM, Moutsianas L, Lee JC, Lamb CA, Luo Y, Kennedy NA, et al. Genome-wide association study implicates immune activation of multiple integrin genes in inflammatory bowel disease. Nat Genet. 2017;49(2):256–61.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Trynka G, Hunt KA, Bockett NA, Romanos J, Mistry V, Szperl A, et al. Dense genotyping identifies and localizes multiple common and rare variant association signals in celiac disease. Nat Genet. 2011;43(12):1193–201.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Eijsbouts C, Zheng T, Kennedy NA, Bonfiglio F, Anderson CA, Moutsianas L, et al. Genome-wide analysis of 53,400 people with irritable bowel syndrome highlights shared genetic pathways with mood and anxiety disorders. Nat Genet. 2021;53(11):1543–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Consortium IMSG, ANZgene, IIBDGC, WTCCC2. Multiple sclerosis genomic map implicates peripheral immune cells and microglia in susceptibility. Science. 2019;365(6460):eaav7188.

    Article  Google Scholar 

  30. Cordell HJ, Fryett JJ, Ueno K, Darlay R, Aiba Y, Hitomi Y, et al. An international genome-wide meta-analysis of primary biliary cholangitis: novel risk loci and candidate drugs. J Hepatol. 2021;75(3):572–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Ji S-G, Juran BD, Mucha S, Folseraas T, Jostins L, Melum E, et al. Genome-wide association study of primary sclerosing cholangitis identifies new risk loci and quantifies the genetic relationship with inflammatory bowel disease. Nat Genet. 2017;49(2):269–73.

    Article  CAS  PubMed  Google Scholar 

  32. Tsoi LC, Spain SL, Knight J, Ellinghaus E, Stuart PE, Capon F, et al. Identification of 15 new psoriasis susceptibility loci highlights the role of innate immunity. Nat Genet. 2012;44(12):1341–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Okada Y, Wu D, Trynka G, Raj T, Terao C, Ikari K, et al. Genetics of rheumatoid arthritis contributes to biology and drug discovery. Nature. 2014;506(7488):376–81.

    Article  CAS  PubMed  Google Scholar 

  34. Forgetta V, Manousaki D, Istomine R, Ross S, Tessier M-C, Marchand L, et al. Rare genetic variants of large effect influence risk of type 1 diabetes. Diabetes. 2020;69(4):784–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Bentham J, Morris DL, Graham DSC, Pinder CL, Tombleson P, Behrens TW, et al. Genetic association analyses implicate aberrant regulation of innate and adaptive immunity genes in the pathogenesis of systemic lupus erythematosus. Nat Genet. 2015;47(12):1457–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Hemani G, Tilling K, Davey SG. Orienting the causal relationship between imprecisely measured traits using GWAS summary data. PLoS Genet. 2017;13(11):e1007081.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, et al. The MR-Base platform supports systematic causal inference across the human phenome. Elife. 2018;7:e34408.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Matzaraki V, Kumar V, Wijmenga C, Zhernakova A. The MHC locus and genetic susceptibility to autoimmune and infectious diseases. Genome Biol. 2017;18(1):1–21.

    Article  Google Scholar 

  39. Burgess S, Thompson SG, Collaboration CCG. Avoiding bias from weak instruments in Mendelian randomization studies. Int J Epidemiol. 2011;40(3):755–64.

    Article  PubMed  Google Scholar 

  40. Burgess S. Sample size and power calculations in Mendelian randomization with a single instrumental variable and a binary outcome. Int J Epidemiol. 2014;43(3):922–9.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Online sample size and power calculator for Mendelian randomization with a binary outcome. Accessed 20 March 2022.

  42. Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. 2015;44(2):512–25.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Burgess S, Butterworth A, Thompson SG. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol. 2013;37(7):658–65.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Burgess S, Bowden J, Fall T, Ingelsson E, Thompson SG. Sensitivity analyses for robust causal inference from Mendelian randomization analyses with multiple genetic variants. Epidemiology. 2017;28(1):30.

    Article  PubMed  Google Scholar 

  45. Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent estimation in Mendelian randomization with some invalid instruments using a weighted median estimator. Genet Epidemiol. 2016;40(4):304–14.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Bowden J, Del Greco MF, Minelli C, Zhao Q, Lawlor DA, Sheehan NA, et al. Improving the accuracy of two-sample summary-data Mendelian randomization: moving beyond the NOME assumption. Int J Epidemiol. 2019;48(3):728–42.

    Article  PubMed  Google Scholar 

  47. Verbanck M, Chen C-y, Neale B, Do R. Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet. 2018;50(5):693–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. McGovern DP, Gardet A, Törkvist L, Goyette P, Essers J, Taylor KD, et al. Genome-wide association identifies multiple ulcerative colitis susceptibility loci. Nat Genet. 2010;42(4):332–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Burgess S, Thompson SG. Multivariable Mendelian randomization: the use of pleiotropic genetic variants to estimate causal effects. Am J Epidemiol. 2015;181(4):251–60.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Burgess S, Smith GD, Davies NM, Dudbridge F, Gill D, Glymour MM, et al. Guidelines for performing Mendelian randomization investigations. Wellcome Open Res. 2019;4.

  51. Bowden J, Hemani G, Davey SG. Invited commentary: detecting individual and global horizontal pleiotropy in Mendelian randomization—a job for the humble heterogeneity statistic? Am J Epidemiol. 2018;187(12):2681–5.

    PubMed  PubMed Central  Google Scholar 

  52. Li CY, Yang TM, Ou RW, Wei QQ, Shang HF. Genome-wide genetic links between amyotrophic lateral sclerosis and autoimmune diseases. BMC Med. 2021;19(1):1–11.

    Article  Google Scholar 

Download references


We would like to thank the relevant consortia for making their data available. The NHGRI-EBI GWAS Catalog is funded by NHGRI Grant Number 2U41HG007823 and delivered by collaboration between the NHGRI, EMBL-EBI, and NCBI. KS is supported by a post-doctoral fellowship from the Canada First Research Excellence Fund (CFREF), awarded to McGill University for the Healthy Brains for Healthy Lives Initiative (HBHL), and FRQC post-doctoral fellowship. DM is a Fonds de Recherche du Quebec- Sante (FRQS) Junior 1 Scholar and has received a career development award from the Canadian Child Health Clinician Scientist Program (CCHCSP).


This project is funded by the foundation grant from the Canadian Institute of Health Research (332971).

Author information

Authors and Affiliations



PA, GAR, and PAD conceived the idea for the study. PA, KS, and JPR obtained the GWAS datasets. PA, KS, DM, and DS performed the data analyses. PA, DM, KS, and JPR interpreted the results of the data analyses. All authors wrote the manuscript. The authors read and approved the final manuscript.

Corresponding author

Correspondence to Guy A. Rouleau.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1.

Funnel plots and LOO plots which detect outlier SNPs in RA and CD. Figure S2. Scatter plots showing the effect of liability to ALS on risk of autoimmune disorders. Table S1. Sensitivity analysis, heterogeneity, and pleiotropy, investigating MR assumption violation. Table S2. Sensitivity analysis, heterogeneity, and pleiotropy, investigating MR reverse assumption violation. Code used in this study includes all the code used to produce all the results in this paper.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Alipour, P., Senkevich, K., Ross, J.P. et al. Investigation of the causal relationship between ALS and autoimmune disorders: a Mendelian randomization study. BMC Med 20, 382 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Mendelian randomization
  • Amyotrophic lateral sclerosis
  • Autoimmune disorders
  • Causal relationship