The effect of smoking on multiple sclerosis: a mendelian randomization study
============================================================================

* Ruth E Mitchell
* Kirsty Bates
* Robyn E Wootton
* Adil Harroud
* J. Brent Richards
* George Davey Smith
* Marcus R Munafò

## Abstract

The causes of multiple sclerosis (MS) remain unknown. Smoking has been associated with MS in observational studies and is often thought of as an environmental risk factor. We used two-sample Mendelian Randomization (MR) to examined whether this association is causal using genetic variants identified in genome-wide association studies (GWAS) as associated with smoking. We assessed both smoking initiation and lifetime smoking behaviour (which captures smoking duration, heaviness and cessation). There was very limited evidence for a meaningful effect of smoking on MS susceptibility was measured using summary statistics from the International Multiple Sclerosis Genetics Consortium (IMSGC) meta-analysis, including 14,802 cases and 26,703 controls. There was no clear evidence for an effect of smoking on the risk of developing MS (smoking initiation: odds ratio [OR] 1.03, 95% confidence interval [CI] 0.92-1.61; lifetime smoking: OR 1.10, 95% CI 0.87-1.40). These findings suggest that smoking does not have a detrimental consequence on MS susceptibility. Further work is needed to determine the causal effect of smoking on MS progression.

## Background

Smoking is an avoidable environmental cause to many life-threatening diseases such as lung cancer, heart and respiratory disorders (1,2). There is emerging evidence linking cigarette smoke to conditions negatively affecting the central nervous system (CNS), like multiple sclerosis (MS) (3,4). MS is a chronic neurological disorder causing autoimmune breakdown of the myelin sheath surrounding axons in the CNS (5). The disease is characterised by periods of disease activity followed by remission and/or progressive neurological decline, resulting in increasing disability (6). Like most autoimmune conditions, there is no known specific cause; however, we know there is an interaction between genetic and environmental factors in susceptible individuals, that go on to develop the disorder (7). Unfortunately, there is no cure for MS(8), and people diagnosed often live with extreme disability (9). There are emerging treatments aimed at modifying the disease course (10), but they are not universally effective particularly with regards to the progressive form of the disease. Therefore, it is important to continue targeting prevention by means of establishing causal links. Evidence from observational epidemiological studies suggests that smoking increases MS risk (11). It is hypothesised from experimental studies that exposure to chemicals in cigarette smoke alter the immune cell balance in the lung (12,13) which in turn can lead to generalised pro-inflammatory effects that trigger autoimmunity (14,15), in genetically susceptible individuals (16,17). In addition, cigarette chemicals contribute mechanistically to MS pathobiology. Specifically, nicotine increases the permeability of the blood-brain barrier (18); cyanide contributes to demyelination (19); and nitric oxide causes degeneration of axons (20). There is evidence for an association between smoking and worsening symptoms, number of relapses, lesion load on MRI, brain atrophy rate (15) and the rapidity of disability progression in MS patients (4,21,22).

However, it is hard to make causal inferences from observational studies which can be biased by issues of reverse causation and residual confounding. One method which can be used to reduce these sources of bias is Mendelian Randomization (MR) (23). MR can be implemented through instrumental variable analysis that uses genetic variants to proxy the exposure (e.g, smoking) and estimate a causal effect of that exposure on the outcome (e.g, MS). The MR method makes three important assumptions: 1) the genetic variants must robustly predict the exposure, 2) the genetic variants must not be associated with any confounders and 3) the genetic variants must only affect the outcome through the exposure (24). To satisfy the first assumption we selected the most recently available genetic instruments from previously conducted genome-wide association studies (GWAS) associated with smoking behaviour (smoking initiation (25) and lifetime smoking (26)) that can be implemented in a two-sample MR context (Fig 1). The latter two assumptions can be violated by horizontal pleiotropy which occurs when the genetic variants affect the outcome other than through the exposure. We test for this possibility using multiple sensitivity analyses.

![Fig 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/06/24/2020.06.24.20138834/F1.medium.gif)

[Fig 1.](http://medrxiv.org/content/early/2020/06/24/2020.06.24.20138834/F1)

Fig 1. Directed acyclic graph of the Mendelian randomization framework investigating the causal relationship between smoking and multiple sclerosis.
Instrumental variable assumptions: IV1: the instruments must be associated with the exposure; IV2: the instruments must influence MS only through smoking; IV3: the instruments must not associate with measured or unmeasured confounders in the smoking to MS relationship.

## Results

### Smoking initiation

The inverse-variance weighted MR estimate (OR 1.03, 95% confidence interval (CI) 0.92 − 1.16) revealed no strong evidence for a causal effect of the genetic risk of smoking initiation on incidence of MS (Fig 3). This was consistent across all MR methods employed, providing further support for the result as each MR method has different assumptions and therefore tests for different violations of those assumptions. Indeed, the weighted median and weighted mode only allow SNPs in the largest homogeneous cluster to contribute to the overall estimate and provide estimates with confidence intervals overlapping the null (Fig 3 and Supplementary Fig 1). The 371 SNPs used as genetic proxies for smoking initiation (Fig 2a and Supplementary Table 1) had an F statistic of 44.90 indicating a strong instrument and that weak instrument bias was unlikely to be influencing the effect estimates. There was evidence of heterogeneity with a large Cochran’s Q statistic of 559.48, p=6.65×10−6 and the MR-PRESSO global test value of 562.12, p<0.000125. However, this not indicative of directional horizontal pleiotropy given the consistent MR Egger estimate (OR 1.13, 95% CI 0.67 to 1.91), small intercept (0.0017, p=0.73) and symmetrical funnel plot (Supplementary Fig 2). Similarly, MR-RAPS is robust to systematic and idiosyncratic pleiotropy, accounting for weak instruments, pleiotropy and extreme outliers, and gave a similar causal estimate (OR 1.05, 95% CI 0.93 to 1.17). Furthermore, MR-PRESSO removes individual SNPs that contribute to heterogeneity disproportionately more than expected in order to reduce heterogeneity. The MR-PRESSO outlier corrected causal estimate was 1.040 (95% CI 1.040 to 1.041). Therefore, the second IV assumption (known as the exclusion restriction assumption) of MR has not been violated and directional pleiotropy is unlikely to be biasing the estimates, even though the outlier removal automatically leads to over precise estimates. Leave-one-out and single SNP analyses (Supplementary Fig 3 and 4) were conducted as sensitivity tests sequentially omitting one SNP at a time and performing MR using a single SNP respectively to assess the sensitivity of the results to individual variants. These indicated that there is not a single SNP driving the association whose effect is being masked in the overall analysis. The exclusion of exposure variants located within the MHC did not alter the null association between smoking initiation and incidence of MS (Supplementary Table 2).

![Fig 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/06/24/2020.06.24.20138834/F2.medium.gif)

[Fig 2.](http://medrxiv.org/content/early/2020/06/24/2020.06.24.20138834/F2)

Fig 2. Flowchart for selection of genetic variants associated with smoking initiation (A) and lifetime smoking (B).
Abbreviations: GSCAN: GWAS & Sequencing Consortium of Alcohol and Nicotine use; GWAS: Genome-Wide Association Study; LD: Linkage Disequilibrium; MS: multiple sclerosis; IMSGC: International Multiple Sclerosis Genetics Consortium; MR: Mendelian Randomization; SNP: Single-Nucleotide Polymorphism.

![Fig 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/06/24/2020.06.24.20138834/F3.medium.gif)

[Fig 3:](http://medrxiv.org/content/early/2020/06/24/2020.06.24.20138834/F3)

Fig 3: Two-sample Mendelian Randomization estimates of the association between smoking initiation and incidence of multiple sclerosis.
Odds ratios are expressed per unit increase in log odds of ever smoking regularly (smoking initiation). MR: Mendelian Randomization; OR=Odds Ratio; CI=Confidence Intervals; p.val: p value.

### Lifetime smoking

There was no clear evidence for a causal effect of the genetic risk of lifetime smoking on incidence of MS (Fig 4). The 125 SNPs used as genetic proxies for lifetime smoking (Fig 2b and Supplementary Table 3) had an F statistic of 44.05 indicating a strong instrument that is unlikely to cause the effect estimates to be affected by weak instrument bias. The inverse-variance weighted MR analysis estimate (OR 1.10, 95% CI 0.87 to 1.40) revealed no strong evidence for a causal effect of the genetic risk of lifetime smoking on incidence of MS and was consistent across all MR methods employed (Fig 4 and Supplementary Fig 5). There was evidence of heterogeneity among the individual SNP effect estimates for lifetime smoking with a large Cochran Q statistic (156.18, p=0.02) and MR-PRESSO global test estimate of 158.4895, p=0.03. However, this was not supported by the symmetrical funnel plot (Supplementary Fig 6) nor by any outliers detected in the MR-PRESSO test. Furthermore, the small MR Egger intercept (−0.003, p=0.69) and consistent MR Egger estimate (OR 1.34, 95% CI 0.49 to 3.65) suggests that the magnitude of potential bias from directional pleiotropy is low. Furthermore, there was no single SNP driving the association whose effect is being masked in the overall estimate as demonstrated by the leave-one-out and single SNP sensitivity analyses (Supplementary Fig 7 and 8). MR excluding the lifetime smoking associated variant located within the MHC region yielded consistent results overlapping the null (Supplementary Table 4).

![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/06/24/2020.06.24.20138834/F4.medium.gif)

[Figure 4:](http://medrxiv.org/content/early/2020/06/24/2020.06.24.20138834/F4)

Figure 4: Two-sample Mendelian Randomization estimates of the association between lifetime smoking and incidence of multiple sclerosis.
Odds ratios are expressed per 1 standard deviation increase of the lifetime smoking index. MR: Mendelian Randomization; OR=Odds Ratio; CI=Confidence Intervals; p.val: p value.

A bidirectional analysis shows that there was no clear evidence that a genetic predisposition to MS is associated with either smoking initiation or lifetime smoking (Supplementary Table 5 and 7). MR of MS associated variants located within the MHC region yielded consistent results overlapping the null (Supplementary Table 6 and 8).

## Discussion

This study uses the MR method to estimate the causal effect of smoking on risk for MS. Using a two-sample MR design in 14,802 MS cases and 26,703 controls, we found little evidence that both genetically predicted smoking initiation and lifetime smoking are associated with MS risk. These findings suggest that smoking is not a clear environmental risk factor for MS susceptibility and are in line with a recent independent study (27). Although a small effect cannot be entirely excluded, the relatively narrow confidence intervals, particularly for smoking initiation, make a clinically relevant effect less likely.

This contradicts previously reported observational studies that show an association with MS risk among smokers, compared to non-smokers, of a meta-analysed effect estimate odds ratio of 1.5 (4,11). The studies included limitations such as self-report MS diagnosis (28), participation rate less than 80% (29–31) and loss to follow up (32). Additionally, observational studies may have heterogenous results due to how smoking status was defined (11). The strength of association and causality between smoking and MS risk has been suggested due to a dose-dependent relationship in duration and intensity of smoking (4,33) as well as from the interaction between compounds present in cigarettes and specific genetic HLA variants, which include the presence of HLA-DRB1*15, the absence of HLA-A*0201 (34) and specific N-acetyltranferase 1 (NAT1) polymorphisms (35). These genetic variants facilitated epitope cross-reactivity and activation of T cells and smoking may strongly influence the risk of MS observed with these HLA genotypes. However, other studies have failed to replicate this interaction (36,37). In order to test this interaction in an MR casual inference context, a factorial MR design in MS patients with and without those alleles would be required. This was not possible in the present study due to the use of GWAS summary statistics. Observational estimates may have also been biased by residual or unmeasured confounding from factors influencing both smoking status and MS. For example, co-morbidities and socioeconomic status may influence the likelihood of being a smoker and having MS (15,38).

Reverse causation could also partly explain the discrepancy between our MR results and observational studies especially as MS onset may occur long before the first clinical symptoms (39). For instance, this prodromal phase is characterized in part by a higher risk of depression and anxiety up to 10 years prior to MS diagnosis (40), and these in turn are associated with a higher rate of smoking. This study sought to reduce bias from confounding and reverse causation by using a MR design given genetic variants are much less associated with confounders than directly measured environmental exposures (41) (here smoking) and genetic variants are fixed over our lifetime ensuring directionality of effect. This is a major strength of this study in establishing causality in the relationship between smoking and MS risk. Additionally, MR reverse direction MR was performed and shows that reverse causation is unlikely to be playing a role. A further strength of this study is the use of robust genetic instruments which are strong predictors of smoking behaviour. Finally, we used multiple MR methods and sensitivity analyses to test for bias from directional horizontal pleiotropy. Our estimates were consistent across these multiple methods, strengthening our conclusions.

The current study cannot inform us about the effects of smoking on MS symptom severity, disability or progression of disease. Indeed, smoking shows an association with disease progression, disease activity (new lesions on magnetic resonance imaging (MRI); clinical relapse rates) and brain atropy (15). Observational studies have shown an association between smoking and progression from relapsing remitting MS to secondary progressive MS with a dose-response relationship (42–46) as well as a faster rate increasing Expanded Disability Status Scale (EDSS) (22). However, more research in this area is needed to for a definitive conclusion of an effect and specific mechanisms of action. As new methods are being developed to assess disease progression using MR (47), when a GWAS of MS progression becomes available, future studies should explore the association between smoking and the different measures of MS progression in a MR framework.

The instrument predicting smoking initiation and lifetime smoking were broadly distinct (only 9 SNPs overlapping). The measure of lifetime smoking exposure takes into account smoking status and, among ever smokers, duration, heaviness and cessation. Although our lifetime smoking instrument captured smoking heaviness in part, however we were unable to explore whether there was a dose-response relationship between the number of cigarettes smoked and the likelihood of developing MS given we were unable to stratify the MS GWAS by smoking status. Most, but not all (30,48,49), evidence to date seems to suggest that there is a positive correlation between the amount smoked and the severity of illness (4,31,37,43,50–53). It might be that rather than a causal relationship between smoking and MS risk, that smoking instead accelerates the disease process in those that would have already developed MS.

Limitations of this study are, firstly, that although we assessed pleiotropy using MR methods that account for pleiotropic effects, pleiotropy can only be addressed indirectly, and some SNPs may relate to MS risk through pathways other than smoking. We did not find evidence for bias for horizontal pleiotropy using the MR Egger intercept test nor the funnel plots which did not reveal evidence of directional, or unbalanced, pleiotropy. Secondly, this study was a stwo-sample MR using MS meta-analysis summary statistics and therefore this does not allow for gene-environment interaction or sex stratified analysis.

In conclusion, we find no clear evidence for a causal effect of smoking on the risk of developing MS. Previous observational results may have been due to confounding factors, which we have avoided through our analysis. Future research should focus on the effect of smoking on the disease course of MS and its effect on progression.

## Methods

### Genetic instruments for smoking

#### Smoking initiation

We used the most recent GWAS of smoking initiation from the GWAS & Sequencing Consortium of Alcohol and Nicotine use (GSCAN) consortium which identified 378 conditionally independent genome-wide significant SNPs in a sample of 1,232,091 individuals of European ancestry. These genetic variants explain 2% of the variance in smoking initiation(25).

#### Lifetime smoking

In order to incorporate measures of smoking heaviness without having to stratify on smoking status (which is not possible in the two-sample MR context without a stratified GWAS of MS), we used the GWAS of lifetime smoking conducted in 462,690 individuals of European ancestry from the UK Biobank (26). Lifetime smoking is a combination of smoking initiation, duration, heaviness and cessation described in detail elsewhere (26). This GWAS identified 126 independent genome-wide significant SNPs that explain 0.36% of the variance(26).

### Genetic variants associated with multiple sclerosis

Effect estimates and standard errors for smoking associated SNPs on MS susceptibility were obtained from the summary statistics of the discovery cohorts of the latest International Multiple Sclerosis Genetics Consortium (IMSGC) meta-analysis, including 14,802 cases and 26,703 controls (54). All details relating to demographic characteristics, MS case ascertainment and eligibility criteria for the meta-analysis can be found in the original publication (54). For SNPs not available in the IMSGC dataset, we identified proxy SNPs in high linkage disequilibrium (r2 > 0.8) using an online tool LDlink [[https://ldlink.nci.nih.gov/?tab=ldproxy](https://ldlink.nci.nih.gov/?tab=ldproxy)], giving a total of 371 SNPs for smoking initiation instrument and 125 SNPs for lifetime smoking (Figure 2 and Supplementary Table 1 and 2).

### Mendelian Randomization analyses

A two sample MR was undertaken to obtain effect estimates of genetically predicted smoking on MS susceptibility, using both initiation and lifetime proxy measures. MR and sensitivity analyses were performed in R (version 3.5.1) using the TwoSampleMR R package ([https://mrcieu.github.io/TwoSampleMR/](https://mrcieu.github.io/TwoSampleMR/)) (55) with effect estimates compared across five different methods: inverse variance weighted (IVW); MR Egger (56); weighted median (57); weighted mode (58); robust adjusted profile score (RAPS) (59); pleiotropy residual sum and outlier (PRESSO) (60). Given the different assumptions that each of these methods make about the nature of pleiotropy, consistency in the point estimate across the methods strengthens causal evidence (61). The IWV method is the main analysis and the other methods provide sensitivity analyses. Instrumental variable analysis of MR is based on a ratio of the regressions of the genetic instrument-outcome association (weighted smoking associated SNPs with MS from IMSGC) on the genetic instrument-exposure association (smoking associated SNPs with smoking initiation or lifetime smoking in the independent smoking GWASs). For smoking initiation, the odds ratios are expressed per unit increase in log odds of ever smoking regularly (smoking initiation); for lifetime smoking, the odds ratios are expressed per 1 standard deviation increase of the lifetime smoking index.

Additional sensitivity analyses were performed in order to formally test for potential violations of MR assumptions. The mean F statistic was calculated as an indicator of instrument strength (a value of >10 indicates a strong instrument) and the Cochran’s Q statistic was assessed as a measure of heterogeneity for the IVW method to estimate whether the individual SNP effects of smoking on MS were inconsistent. The MR Egger intercept was assessed to detect directional pleiotropy where the genetic instruments would be influencing MS through another pathway other than smoking. To identify potentially influential SNPs, which could be driven for example by horizontal pleiotropy, we used leave-one-out and single-SNP MR analyses. Additionally, due to the strong genetic signal for MS within the MHC region and high potential for pleiotropy, MR analysis excluding exposure variants located within the extended major histocompatibility (MHC) region was performed (defined as base positions 24,000,000 to 35,000,000 on chromosome 6 [GRCh37].

## Data Availability

The data underlying the results presented in the study are available from the GWAS studies cited.

*   Received June 24, 2020.
*   Revision received June 24, 2020.
*   Accepted June 24, 2020.


*   © 2020, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/)

## References

1.  1.O’Keeffe LM, Taylor G, Huxley RR, Mitchell P, Woodward M, Peters SAE. Smoking as a risk factor for lung cancer in women and men: A systematic review and metaanalysis. Vol. 8, BMJ Open. BMJ Publishing Group; 2018.
    
    

2.  2.Aune D, Schlesinger S, Norat T, Riboli E. Tobacco smoking and the risk of sudden cardiac death: a systematic review and meta-analysis of prospective studies. Eur J Epidemiol [Internet]. 2018 Jun 7 [cited 2020 Jan 20];33(6):509–21. Available from: [http://link.springer.com/10.1007/s10654-017-0351-y](http://link.springer.com/10.1007/s10654-017-0351-y)
    
    

3.  3.George MF, Briggs FBS, Shao X, Gianfrancesco MA, Kockum I, Harbo HF, et al. Multiple sclerosis risk loci and disease severity in 7,125 individuals from 10 studies. Neurol Genet [Internet]. 2016 Aug [cited 2017 Mar 28];2(4):e87. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/27540591](http://www.ncbi.nlm.nih.gov/pubmed/27540591)
    
    

4.  4.Poorolajal J, Bahrami M, Karami M, Hooshmand E. Effect of smoking on multiple sclerosis: A meta-analysis. J Public Heal (United Kingdom). 2017 Jun 1;39(2):312–20.
    
    

5.  5.Dobson R, Giovannoni G. Multiple sclerosis – a review [Internet]. Vol. 26, European Journal of Neurology. Blackwell Publishing Ltd; 2019 [cited 2020 Mar 31]. p. 27–40. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/30300457](http://www.ncbi.nlm.nih.gov/pubmed/30300457)
    
    

6.  6.Filippi M, Bar-Or A, Piehl F, Preziosa P, Solari A, Vukusic S, et al. Multiple sclerosis. Nat Rev Dis Prim. 2018 Dec 1;4(1):1–27.
    
    

7.  7.Olsson T, Barcellos LF, Alfredsson L. Interactions between genetic, lifestyle and environmental risk factors for multiple sclerosis. Vol. 13, Nature Reviews Neurology. Nature Publishing Group; 2016. p. 26–36.
    
    

8.  8.Gohil K. Multiple sclerosis: Progress, but no cure. Vol. 40, P and T. Medi Media USA Inc; 2015. p. 604–5.
    
    

9.  9.Kister I, Chamot E, Salter AR, Cutter GR, Bacon TE, Herbert J. Disability in multiple sclerosis: A reference for patients and clinicians. Neurology. 2013 Mar 12;80(11):1018–24.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1212/WNL.0b013e3182872855&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23427319&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F24%2F2020.06.24.20138834.atom) 

10. 10.Scolding N, Barnes D, Cader S, Chataway J, Chaudhuri A, Coles A, et al. Association of British Neurologists: Revised (2015) guidelines for prescribing disease-modifying treatments in multiple sclerosis. Pract Neurol. 2015 Jan 8;15(4):273–9.
    
    [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MTE6InByYWN0bmV1cm9sIjtzOjU6InJlc2lkIjtzOjg6IjE1LzQvMjczIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDYvMjQvMjAyMC4wNi4yNC4yMDEzODgzNC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

11. 11.Degelman ML, Herman KM. Smoking and multiple sclerosis: A systematic review and meta-analysis using the Bradford Hill criteria for causation. Vol. 17, Multiple Sclerosis and Related Disorders. Elsevier B.V.; 2017. p. 207–16.
    
    

12. 12.Ammitzbøll C, Börnsen L, Romme Christensen J, Ratzer R, Romme Nielsen B, Søndergaard HB, et al. Smoking reduces circulating CD26 hi CD161 hi MAIT cells in healthy individuals and patients with multiple sclerosis. J Leukoc Biol. 2017 May;101(5):1211–20.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1189/jlb.3A0616-267R&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28179539&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F24%2F2020.06.24.20138834.atom) 

13. 13.Makrygiannakis D, Hermansson M, Ulfgren AK, Nicholas AP, Zendman AJW, Eklund A, et al. Smoking increases peptidylarginine deiminase 2 enzyme expression in human lungs and increases citrullination in BAL cells. Ann Rheum Dis [Internet]. 2008 Oct [cited 2020 Apr 2];67(10):1488–92. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/18413445](http://www.ncbi.nlm.nih.gov/pubmed/18413445)
    
    

14. 14.Odoardi F, Sie C, Streyl K, Ulaganathan VK, Schläger C, Lodygin D, et al. T cells become licensed in the lung to enter the central nervous system. Nature [Internet]. 2012 Aug 30 [cited 2020 Apr 2];488(7413):675–9. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/22914092](http://www.ncbi.nlm.nih.gov/pubmed/22914092)
    
    

15. 15.Rosso M, Chitnis T. Association between Cigarette Smoking and Multiple Sclerosis: A Review [Internet]. Vol. 77, JAMA Neurology. American Medical Association; 2020 [cited 2020 Apr 2]. p. 245–53. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/31841592](http://www.ncbi.nlm.nih.gov/pubmed/31841592)
    
    

16. 16.Hedström AK, Katsoulis M, Hössjer O, Bomfim IL, Oturai A, Sondergaard HB, et al. The interaction between smoking and HLA genes in multiple sclerosis: replication and refinement. Eur J Epidemiol. 2017 Oct 1;32(10):909–19.
    
    

17. 17.Sawcer S, Hellenthal G. The major histocompatibility complex and multiple sclerosis: a smoking gun? Brain [Internet]. 2011 Mar [cited 2020 Jan 20];134(3):638–40. Available from: [https://academic.oup.com/brain/article-lookup/doi/10.1093/brain/awq384](https://academic.oup.com/brain/article-lookup/doi/10.1093/brain/awq384)
    
    

18. 18.Chen JL, Wei L, Bereczki D, Hans FJ, Otsuka T, Acuff V, et al. Nicotine raises the influx of permeable solutes across the rat blood-brain barrier with little or no capillary recruitment. J Cereb Blood Flow Metab [Internet]. 1995 Jul 29 [cited 2020 Apr 2];15(4):687–98. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/7790419](http://www.ncbi.nlm.nih.gov/pubmed/7790419)
    
    

19. 19.Philbrick DJ, Hopkins JB, Hill DC, Alexander JC, Thomson RG. Effects of prolonged cyanide and thiocyanate feeding in rats. J Toxicol Environ Health. 1979;5(4):579–92.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=490674&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F24%2F2020.06.24.20138834.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1979HF30600001&link_type=ISI) 

20. 20.Smith KJ, Kapoor R, Hall SM, Davies M. Electrically active axons degenerate when exposed to nitric oxide. Ann Neurol [Internet]. 2001 Apr 1 [cited 2020 Apr 2];49(4):470–6. Available from: [http://doi.wiley.com/10.1002/ana.96](http://doi.wiley.com/10.1002/ana.96)
    
    

21. 21.Wingerchuk DM. Smoking: Effects on multiple sclerosis susceptibility and disease progression. Vol. 5, Therapeutic Advances in Neurological Disorders. 2012. p. 13–22.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/1756285611425694&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22276073&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F24%2F2020.06.24.20138834.atom) 

22. 22.Heydarpour P, Manouchehrinia A, Beiki O, Mousavi SE, Abdolalizadeh A, Lakeh MM, et al. Smoking and worsening disability in multiple sclerosis: A meta-analysis. Acta Neurol Scand. 2018 Jul 1;138(1):62–9.
    
    

23. 23.Davey Smith G, Ebrahim S. ‘Mendelian randomization’: can genetic epidemiology contribute to understanding environmental determinants of disease?*. Int J Epidemiol [Internet]. 2003 Feb [cited 2019 Jan 28];32(1):1–22. Available from: [https://academic.oup.com/ije/article-lookup/doi/10.1093/ije/dyg070](https://academic.oup.com/ije/article-lookup/doi/10.1093/ije/dyg070)
    
    

24. 24.Davey Smith G, Hemani G. Mendelian randomization: genetic anchors for causal inference in epidemiological studies. Hum Mol Genet [Internet]. 2014 Sep 15 [cited 2019 Jan 28];23(R1):R89–98. Available from: [https://academic.oup.com/hmg/article-lookup/doi/10.1093/hmg/ddu328](https://academic.oup.com/hmg/article-lookup/doi/10.1093/hmg/ddu328)
    
    

25. 25.Liu M, Jiang Y, Wedow R, Li Y, Brazel DM, Chen F, et al. Association studies of up to 1.2 million individuals yield new insights into the genetic etiology of tobacco and alcohol use. Vol. 51, Nature Genetics. Nature Publishing Group; 2019. p. 237–44.
    
    

26. 26.Wootton RE, Richmond RC, Stuijfzand BG, Lawn RB, Sallis HM, Taylor GMJ, et al. Evidence for causal effects of lifetime smoking on risk for depression and schizophrenia: a Mendelian randomisation study. Psychol Med [Internet]. 2019 Nov 6 [cited 2020 Jan 20];1–9. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/31689377](http://www.ncbi.nlm.nih.gov/pubmed/31689377)
    
    

27. 27.Vandebergh M, Goris A. Smoking and multiple sclerosis risk: a Mendelian randomization study. J Neurol [Internet]. 2020 [cited 2020 Jun 22]; Available from: [https://pubmed.ncbi.nlm.nih.gov/32529581/](https://pubmed.ncbi.nlm.nih.gov/32529581/)
    
    

28. 28.Riise T, Nortvedt MW, Ascherio A. Smoking is a risk factor for multiple sclerosis. Neurology. 2003 Oct 28;61(8):1122–4.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1212/01.WNL.0000081305.66687.D2&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14581676&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F24%2F2020.06.24.20138834.atom) 

29. 29.O’Gorman C, Broadley SA. Smoking and multiple sclerosis: evidence for latitudinal and temporal variation. J Neurol. 2014;261(9):1677–83.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00415-014-7397-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24923244&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F24%2F2020.06.24.20138834.atom) 

30. 30.Salzer J, Hallmans G, Nyström M, Stenlund H, Wadell G, Sundström P. Smoking as a risk factor for multiple sclerosis. Mult Scler J. 2013;19(8):1022–7.
    
    

31. 31.Hedström AK, Bäärnhielm M, Olsson T, Alfredsson L. Tobacco smoking, but not Swedish snuff use, increases the risk of multiple sclerosis. Neurology. 2009 Sep 1;73(9):696–701.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1212/WNL.0b013e3181b59c40&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19720976&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F24%2F2020.06.24.20138834.atom) 

32. 32.Thorogood M, Hannaford PC. The influence of oral contraceptives on the risk of multiple sclerosis. BJOG An Int J Obstet Gynaecol. 1998;105(12):1296–9.
    
    

33. 33.Hedström AK, Olsson T, Alfredsson L. Smoking is a major preventable risk factor for multiple sclerosis. Mult Scler. 2016 Jul 1;22(8):1021–6.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/1352458515609794&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26459151&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F24%2F2020.06.24.20138834.atom) 

34. 34.Hedström AK, Sundqvist E, Bä M, Nordin N, Hillert J, Kockum I, et al. Smoking and two human leukocyte antigen genes interact to increase the risk for multiple sclerosis. A J Neurol [Internet]. [cited 2020 Apr 28]; Available from: [https://academic.oup.com/brain/article-abstract/134/3/653/446931](https://academic.oup.com/brain/article-abstract/134/3/653/446931)
    
    

35. 35.Briggs FBS, Acuna B, Shen L, Ramsay P, Quach H, Bernstein A, et al. Smoking and risk of multiple sclerosis: Evidence of modification by NAT1 variants. Epidemiology. 2014;25(4):605–14.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/EDE.0000000000000089&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24625537&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F24%2F2020.06.24.20138834.atom) 

36. 36.Petersen ER, Oturai AB, Koch-Henriksen N, Magyari M, Sørensen PS, Sellebjerg F, et al. Smoking affects the interferon beta treatment response in multiple sclerosis. Neurology. 2018 Feb 13;90(7):e593–600.
    
    

37. 37.Simon KC, Van Der Mei IAF, Munger KL, Ponsonby A, Dickinson J, Dwyer T, et al. Combined effects of smoking, anti-EBNA antibodies, and HLA-DRB1*1501 on multiple sclerosis risk. Neurology. 2010;74(17):1365–71.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1212/WNL.0b013e3181dad57e&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20375311&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F24%2F2020.06.24.20138834.atom) 

38. 38.Wingerchuk DM. Smoking: Effects on multiple sclerosis susceptibility and disease progression. Vol. 5, Therapeutic Advances in Neurological Disorders. SAGE Publications; 2012. p. 13–22.
    
    

39. 39.Disanto G, Zecca C, MacLachlan S, Sacco R, Handunnetthi L, Meier UC, et al. Prodromal symptoms of multiple sclerosis in primary care. Ann Neurol [Internet]. 2018 Jun 1 [cited 2020 Apr 3];83(6):1162–73. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/29740872](http://www.ncbi.nlm.nih.gov/pubmed/29740872)
    
    

40. 40.M F, Ae TM G, Mr M. The Association of Cigarette Smoking With Depression and Anxiety: A Systematic Review. Nicotine Tob Res. 2017;19(1).
    
    

41. 41.Smith GD, Ebrahim S. “Mendelian randomization”: Can genetic epidemiology contribute to understanding environmental determinants of disease? Vol. 32, International Journal of Epidemiology. 2003. p. 1–22.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dyg070&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12689998&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F24%2F2020.06.24.20138834.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000182341300001&link_type=ISI) 

42. 42.Koch M, Van Harten A, Uyttenboogaart M, De Keyser J. Cigarette smoking and progression in multiple sclerosis. Neurology. 2007 Oct;69(15):1515–20.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1212/01.wnl.0000277658.78381.db&link_type=DOI) 

43. 43.Sundström P, Nyström L. Smoking worsens the prognosis in multiple sclerosis. Mult Scler. 2008;14(8):1031–5.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/1352458508093615&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18632778&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F24%2F2020.06.24.20138834.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000258715000004&link_type=ISI) 

44. 44.Healy BC, Ali EN, Guttmann CRG, Chitnis T, Glanz BI, Buckle G, et al. Smoking and disease progression in multiple sclerosis. Arch Neurol. 2009 Jul;66(7):858–64.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/archneurol.2009.122&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19597087&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F24%2F2020.06.24.20138834.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000268189200007&link_type=ISI) 

45. 45.Hernán MA, Jick SS, Logroscino G, Olek MJ, Ascherio A, Jick H. Cigarette smoking and the progression of multiple sclerosis. Brain [Internet]. 2005 Jun [cited 2020 Apr 3];128(Pt 6):1461–5. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/15758034](http://www.ncbi.nlm.nih.gov/pubmed/15758034)
    
    

46. 46.Roudbari SA, Ansar MM, Yousefzad A. Smoking as a risk factor for development of Secondary Progressive Multiple Sclerosis: A study in IRAN, Guilan. J Neurol Sci. 2013 Jul 15;330(1–2):52–5.
    
    

47. 47.Paternoster L, Tilling K, Davey Smith G. Genetic epidemiology and Mendelian randomization for informing disease therapeutics: Conceptual and methodological challenges. Barsh GS, editor. PLOS Genet [Internet]. 2017 Oct 5 [cited 2017 Dec 12];13(10):e1006944. Available from: [http://www.ncbi.nlm.nih.gov/pubmed/28981501](http://www.ncbi.nlm.nih.gov/pubmed/28981501)
    
    

48. 48.Maghzi AH, Etemadifar M, Heshmat-Ghahdarijani K, Moradi V, Nonahal S, Ghorbani A, et al. Cigarette smoking and the risk of multiple sclerosis: A sibling case-control study in Isfahan, Iran. Neuroepidemiology. 2011 Dec;37(3–4):238–42.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1159/000332765&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22156601&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F24%2F2020.06.24.20138834.atom) 

49. 49.Jafari N, Hoppenbrouwers IA, Hop WCJ, Breteler MMB, Hintzen RQ. Cigarette smoking and risk of MS in multiplex families. Mult Scler. 2009;15(11):1363–7.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/1352458509345907&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19825892&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F24%2F2020.06.24.20138834.atom) 

50. 50.Asadollahi S, Fakhri M, Heidari K, Zandieh A, Vafaee R, Mansouri B. Cigarette smoking and associated risk of multiple sclerosis in the Iranian population. J Clin Neurosci. 2013 Dec 1;20(12):1747–50.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jocn.2013.01.018&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23972559&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F24%2F2020.06.24.20138834.atom) 

51. 51.Hedström AK, Hillert J, Olsson T, Alfredsson L. Smoking and multiple sclerosis susceptibility. Eur J Epidemiol. 2013 Nov;28(11):867–74.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10654-013-9853-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24146047&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F24%2F2020.06.24.20138834.atom) 

52. 52.Ghadirian P, Dadgostar B, Azani R, Maisonneuve P. A case-control study of the association between socio-demographic, lifestyle and medical history factors and multiple sclerosis. Can J Public Heal. 2001;92(4):281–5.
    
    

53. 53.Mouhieddine TH, Darwish H, Fawaz L, Yamout B, Tamim H, Khoury SJ. Risk factors for multiple sclerosis and associations with anti-EBV antibody titers. Clin Immunol. 2015 May 1;158(1):59–66.
    
    

54. 54.Patsopoulos NA, Baranzini SE, Santaniello A, Shoostari P, Cotsapas C, Wong G, et al. Multiple sclerosis genomic map implicates peripheral immune cells and microglia in susceptibility. Science (80-). 2019 Sep 27;365(6460).
    
    

55. 55.Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, et al. The MRBase platform supports systematic causal inference across the human phenome. Elife [Internet]. 2018 May 30 [cited 2019 Jan 28];7. Available from: [https://elifesciences.org/articles/34408](https://elifesciences.org/articles/34408)
    
    

56. 56.Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol [Internet]. 2015 Apr 1 [cited 2019 Jan 28];44(2):512–25. Available from: [https://academic.oup.com/ije/article-lookup/doi/10.1093/ije/dyv080](https://academic.oup.com/ije/article-lookup/doi/10.1093/ije/dyv080)
    
    

57. 57.Bowden J, Smith GD, Haycock PC, Burgess S. Consistent Estimation in Mendelian Randomization with Some Invalid Instruments Using a Weighted Median Estimator. Genet Epidemiol [Internet]. 2016 [cited 2018 Nov 16];40:304–14. Available from: [https://onlinelibrary.wiley.com/doi/pdf/10.1002/gepi.21965](https://onlinelibrary.wiley.com/doi/pdf/10.1002/gepi.21965)
    
    

58. 58.Hartwig FP, Davey Smith G, Bowden J. Robust inference in summary data Mendelian randomization via the zero modal pleiotropy assumption. Int J Epidemiol [Internet]. 2017 Dec 1 [cited 2019 Jan 28];46(6):1985–98. Available from: [https://academic.oup.com/ije/article/46/6/1985/3957932](https://academic.oup.com/ije/article/46/6/1985/3957932)
    
    

59. 59.Zhao Q, Wang J, Hemani G, Bowden J, Small DS. Statistical inference in two-sample summary-data Mendelian randomization using robust adjusted profile score. 2018 Jan 29 [cited 2020 Jan 23]; Available from: [http://arxiv.org/abs/1801.09652](http://arxiv.org/abs/1801.09652)
    
    

60. 60.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 [Internet]. 2018 May 23 [cited 2019 Jan 28];50(5):693–8. Available from: [http://www.nature.com/articles/s41588-018-0099-7](http://www.nature.com/articles/s41588-018-0099-7)
    
    

61. 61.Lawlor DA, Tilling K, Smith GD. Triangulation in aetiological epidemiology. Int J Epidemiol. 2016 Dec 1;45(6):1866–86.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dyw314&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28108528&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F24%2F2020.06.24.20138834.atom)