Expression profiling of cerebrospinal fluid identifies dysregulated antiviral mechanisms in multiple sclerosis.

Despite the overwhelming evidence that multiple sclerosis is an autoimmune disease, relatively little is known about the precise nature of the immune dysregulation underlying the development of the disease. Reasoning that the CSF from patients might be enriched for cells relevant in pathogenesis, we have completed a high-resolution single-cell analysis of 96 732 CSF cells collected from 33 patients with multiple sclerosis (n = 48 675) and 48 patients with other neurological diseases (n = 48 057). Completing comprehensive cell type annotation, we identified a rare population of CD8+ T cells, characterized by the upregulation of inhibitory receptors, increased in patients with multiple sclerosis. Applying a Multi-Omics Factor Analysis to these single-cell data further revealed that activity in pathways responsible for controlling inflammatory and type 1 interferon responses are altered in multiple sclerosis in both T cells and myeloid cells. We also undertook a systematic search for expression quantitative trait loci in the CSF cells. Of particular interest were two expression quantitative trait loci in CD8+ T cells that were fine mapped to multiple sclerosis susceptibility variants in the viral control genes ZC3HAV1 (rs10271373) and IFITM2 (rs1059091). Further analysis suggests that these associations likely reflect genetic effects on RNA splicing and cell-type specific gene expression respectively. Collectively, our study suggests that alterations in viral control mechanisms might be important in the development of multiple sclerosis.

ADAMS project: a genetic Association study in individuals from Diverse Ancestral backgrounds with Multiple Sclerosis based in the UK.

PURPOSE: Genetic studies of multiple sclerosis (MS) susceptibility and severity have focused on populations of European ancestry. Studying MS genetics in other ancestral groups is necessary to determine the generalisability of these findings. The genetic Association study in individuals from Diverse Ancestral backgrounds with Multiple Sclerosis (ADAMS) project aims to gather genetic and phenotypic data on a large cohort of ancestrally-diverse individuals with MS living in the UK. PARTICIPANTS: Adults with self-reported MS from diverse ancestral backgrounds. Recruitment is via clinical sites, online (https://app.mantal.co.uk/adams) or the UK MS Register. We are collecting demographic and phenotypic data using a baseline questionnaire and subsequent healthcare record linkage. We are collecting DNA from participants using saliva kits (Oragene-600) and genotyping using the Illumina Global Screening Array V.3. FINDINGS TO DATE: As of 3 January 2023, we have recruited 682 participants (n=446 online, n=55 via sites, n=181 via the UK MS Register). Of this initial cohort, 71.2% of participants are female, with a median age of 44.9 years at recruitment. Over 60% of the cohort are non-white British, with 23.5% identifying as Asian or Asian British, 16.2% as Black, African, Caribbean or Black British and 20.9% identifying as having mixed or other backgrounds. The median age at first symptom is 28 years, and median age at diagnosis is 32 years. 76.8% have relapsing-remitting MS, and 13.5% have secondary progressive MS. FUTURE PLANS: Recruitment will continue over the next 10 years. Genotyping and genetic data quality control are ongoing. Within the next 3 years, we aim to perform initial genetic analyses of susceptibility and severity with a view to replicating the findings from European-ancestry studies. In the long term, genetic data will be combined with other datasets to further cross-ancestry genetic discoveries.

Multivariate profile and acute-phase correlates of cognitive deficits in a COVID-19 hospitalised cohort.

Background: Preliminary evidence has highlighted a possible association between severe COVID-19 and persistent cognitive deficits. Further research is required to confirm this association, determine whether cognitive deficits relate to clinical features from the acute phase or to mental health status at the point of assessment, and quantify rate of recovery. Methods: 46 individuals who received critical care for COVID-19 at Addenbrooke's hospital between 10th March 2020 and 31st July 2020 (16 mechanically ventilated) underwent detailed computerised cognitive assessment alongside scales measuring anxiety, depression and post-traumatic stress disorder under supervised conditions at a mean follow up of 6.0 (± 2.1) months following acute illness. Patient and matched control (N = 460) performances were transformed into standard deviation from expected scores, accounting for age and demographic factors using N = 66,008 normative datasets. Global accuracy and response time composites were calculated (G_SScore & G_RT). Linear modelling predicted composite score deficits from acute severity, mental-health status at assessment, and time from hospital admission. The pattern of deficits across tasks was qualitatively compared with normal age-related decline, and early-stage dementia. Findings: COVID-19 survivors were less accurate (G_SScore=-0.53SDs) and slower (G_RT=+0.89SDs) in their responses than expected compared to their matched controls. Acute illness, but not chronic mental health, significantly predicted cognitive deviation from expected scores (G_SScore (p=​​0.0037) and G_RT (p = 0.0366)). The most prominent task associations with COVID-19 were for higher cognition and processing speed, which was qualitatively distinct from the profiles of normal ageing and dementia and similar in magnitude to the effects of ageing between 50 and 70 years of age. A trend towards reduced deficits with time from illness (r∼=0.15) did not reach statistical significance. Interpretation: Cognitive deficits after severe COVID-19 relate most strongly to acute illness severity, persist long into the chronic phase, and recover slowly if at all, with a characteristic profile highlighting higher cognitive functions and processing speed. Funding: This work was funded by the National Institute for Health Research (NIHR) Cambridge Biomedical Research Centre (BRC), NIHR Cambridge Clinical Research Facility (BRC-1215-20014), the Addenbrooke's Charities Trust and NIHR COVID-19 BioResource RG9402. AH is funded by the UK Dementia Research Institute Care Research and Technology Centre and Imperial College London Biomedical Research Centre. ETB and DKM are supported by NIHR Senior Investigator awards. JBR is supported by the Wellcome Trust (220258) and Medical Research Council (SUAG/051 G101400). VFJN is funded by an Academy of Medical Sciences/ The Health Foundation Clinician Scientist Fellowship. The views expressed are those of the authors and not necessarily those of the NHS, the NIHR or the Department of Health and Social Care.

Genetic associations at regulatory phenotypes improve fine-mapping of causal variants for 12 immune-mediated diseases.

The resolution of causal genetic variants informs understanding of disease biology. We used regulatory quantitative trait loci (QTLs) from the BLUEPRINT, GTEx and eQTLGen projects to fine-map putative causal variants for 12 immune-mediated diseases. We identify 340 unique loci that colocalize with high posterior probability (≥98%) with regulatory QTLs and apply Bayesian frameworks to fine-map associations at each locus. We show that fine-mapping credible sets derived from regulatory QTLs are smaller compared to disease summary statistics. Further, they are enriched for more functionally interpretable candidate causal variants and for putatively causal insertion/deletion (INDEL) polymorphisms. Finally, we use massively parallel reporter assays to evaluate candidate causal variants at the ITGA4 locus associated with inflammatory bowel disease. Overall, our findings suggest that fine-mapping applied to disease-colocalizing regulatory QTLs can enhance the discovery of putative causal disease variants and enhance insights into the underlying causal genes and molecular mechanisms.

Mononeuritis multiplex: an unexpectedly frequent feature of severe COVID-19.

The prolonged mechanical ventilation that is often required by patients with severe COVID-19 is expected to result in significant intensive care unit-acquired weakness (ICUAW) in many of the survivors. However, in our post-COVID-19 follow-up clinic we have found that, as well as the anticipated global weakness related to loss of muscle mass, a significant proportion of these patients also have disabling focal neurological deficits relating to multiple axonal mononeuropathies. Amongst the 69 patients with severe COVID-19 that have been discharged from the intensive care units in our hospital, we have seen 11 individuals (16%) with such a mononeuritis multiplex. In many instances, the multi-focal nature of the weakness in these patients was initially unrecognised as symptoms were wrongly assumed to relate simply to "critical illness neuromyopathy". While mononeuropathy is well recognised as an occasional complication of intensive care, our experience suggests that such deficits are surprisingly frequent and often disabling in patients recovering from severe COVID-19.

Complex genetic signatures in immune cells underlie autoimmunity and inform therapy.

We report on the influence of ~22 million variants on 731 immune cell traits in a cohort of 3,757 Sardinians. We detected 122 significant (P < 1.28 × 10-11) independent association signals for 459 cell traits at 70 loci (53 of them novel) identifying several molecules and mechanisms involved in cell regulation. Furthermore, 53 signals at 36 loci overlapped with previously reported disease-associated signals, predominantly for autoimmune disorders, highlighting intermediate phenotypes in pathogenesis. Collectively, our findings illustrate complex genetic regulation of immune cells with highly selective effects on autoimmune disease risk at the cell-subtype level. These results identify drug-targetable pathways informing the design of more specific treatments for autoimmune diseases.

Author Correction: Complex genetic signatures in immune cells underlie autoimmunity and inform therapy.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

The glycosyltransferase EXTL2 promotes proteoglycan deposition and injurious neuroinflammation following demyelination.

BACKGROUND: Chondroitin sulfate proteoglycans (CSPGs) are potent inhibitors of axonal regrowth and remyelination. More recently, they have also been highlighted as a modulator of macrophage infiltration into the central nervous system in experimental autoimmune encephalomyelitis, an inflammatory model of multiple sclerosis. METHODS: We interrogated results from single nucleotide polymorphisms (SNPs) lying in or close to genes regulating CSPG metabolism in the summary results from two publicly available systematic studies of multiple sclerosis (MS) genetics. A demyelinating injury model in the spinal cord of exostosin-like 2 deficient  (EXTL2-/-) mice was used to investigate the effects of dysregulation of EXTL2 on remyelination. Cell cultures of bone marrow-derived macrophages and primary oligodendrocyte precursor cells and neurons were supplemented with purified CSPGs or conditioned media to assess potential mechanisms of action. RESULTS: The strongest evidence for genetic association was seen for SNPs mapping to the region containing the glycosyltransferase exostosin-like 2 (EXTL2), an enzyme that normally suppresses CSPG biosynthesis. Six of these SNPs showed genome-wide significant evidence for association in one of the studies with concordant and nominally significant effects in the second study. We then went on to show that a demyelinating injury to the spinal cord of EXTL2-/- mice resulted in excessive deposition of CSPGs in the lesion area. EXTL2-/- mice had exacerbated axonal damage and myelin disruption relative to wild-type mice, and increased representation of microglia/macrophages within lesions. In tissue culture, activated bone marrow-derived macrophages from EXTL2-/- mice overproduce tumor necrosis factor α (TNFα) and matrix metalloproteinases (MMPs). CONCLUSIONS: These results emphasize CSPGs as a prominent modulator of neuroinflammation and they highlight CSPGs accumulating in lesions in promoting axonal injury.

Transcript specific regulation of expression influences susceptibility to multiple sclerosis.

Genome-wide association studies (GWAS) have identified over 100 loci containing single nucleotide variants (SNVs) that influence the risk of developing multiple sclerosis (MS). Most of these loci lie in non-coding regulatory regions of the genome that are active in immune cells and are therefore thought to modify risk by altering the expression of key immune genes. To explore this hypothesis we screened genes flanking MS-associated variants for evidence of allele specific expression (ASE) by quantifying the transcription of coding variants in linkage disequilibrium with MS-associated SNVs. In total, we were able to identify and successfully analyse 200 such coding variants (from 112 genes) in both CD4+ and CD8+ T cells from 106 MS patients and 105 controls. Fifty-six of these coding variants (from 43 genes) showed statistically significant evidence of ASE in one or both cell types. In the Lck interacting transmembrane adaptor 1 gene (LIME1), for example, we were able to show that in both cell types, the MS-associated variant rs2256814 increased the expression of some transcripts while simultaneously reducing the expression of other transcripts. In CD4+ cells from an additional independent set of 96 cases and 93 controls we were able to replicate the effect of this SNV on the balance of alternate LIME1 transcripts using qPCR (p = 5 × 10-24). Our data thus indicate that some of the MS-associated SNVs identified by GWAS likely exert their effects on risk by distorting the balance of alternate transcripts rather than by changing the overall level of gene expression.

Do Cerebral Small Vessel Disease and Multiple Sclerosis Share Common Mechanisms of White Matter Injury?

Background and Purpose- The role of inflammation in ischemic white matter disease is increasingly recognized, and further understanding of the pathophysiology might inform future treatment strategies. Multiple sclerosis (MS) is a chronic autoimmune condition in which inflammation plays a central role that also affects the white matter. We hypothesized that white matter injury might share common mechanisms and used statistical genetics techniques to assess whether having genetically elevated white matter hyperintensity (WMH) volume was associated with increased MS risk. Methods- We investigated the genetic association in 2 cohorts with magnetic resonance imaging-quantified ischemic white matter lesion volume (WMH in stroke; n=2797 and UK Biobank; n=8353) and 14 802 cases of MS and 26 703 controls from the International Multiple Sclerosis Genetics Consortium. We further performed individual-level polygenic risk score calculations for MS and measures of structural white matter disease in UK Biobank. Finally, we looked for evidence of overlapping risk across the whole genome. Results- There was no association of genetic variants influencing MS with WMH volume using summary statistics in the WMH in stroke cohort (relative risk score =1.014; 95% CI, 0.936-1.110) or in the UK Biobank cohort (relative risk score =1.030; 95% CI, 0.932-1.117). Conversely, assessing the contribution of single nucleotide polymorphisms significantly associated with WMH on the risk of MS there was no significant association (relative risk score =0.930; 95% CI, 0.736-1.191). There were no significant associations between polygenic risk scores calculations; these results were robust to the selection of single nucleotide polymorphisms at a range of significance thresholds. Whole genome analysis did not reveal any overlap of risk between the traits. Conclusions- Our results do not provide evidence to suggest a shared mechanism of white matter damage in ischemia and MS. We propose that inflammation acts in distinct pathways because of the differing nature of the primary insult.

Effect of age at puberty on risk of multiple sclerosis: A mendelian randomization study.

OBJECTIVE: To investigate the potential for a causal effect of age at puberty on multiple sclerosis (MS) susceptibility using a mendelian randomization (MR) approach. METHODS: We used 372 genetic variants strongly associated with age at menarche in a genome-wide association study (GWAS) involving 329,245 women. The genetic architecture of pubertal timing across both sexes is highly correlated (genetic correlation [r g] = 0.75, p = 1.2 × 10-79), allowing these variants to provide reliable insight into pubertal timing in males as well. The effect of pubertal timing on risk of MS was measured with summary statistics from a GWAS of 14,802 cases with MS and 26,703 controls from the International Multiple Sclerosis Genetics Consortium. Multivariable MR controlling for effects of body mass index (BMI) using genetic data from additional consortia investigated whether pubertal effects on MS were dependent on weight status. RESULTS: A 1-year increase in genetically predicted age at puberty decreased odds of MS by 8% (odds ratio [OR] 0.92, 95% confidence interval [CI] 0.86-0.99, p = 0.03). However, multivariable MR analysis showed that after accounting for effects on adult BMI, the association of age at puberty with MS susceptibility attenuated (OR 0.96, 95% CI 0.88-1.04, p = 0.36). Similar results were obtained when childhood BMI was incorporated. Sensitivity analyses provided no evidence of major bias from genetic pleiotropy. CONCLUSIONS: We found support for an association between higher age at puberty and decreased risk of MS with a magnitude comparable to that reported in observational studies. This effect appears to be largely mediated by the strong association between age at puberty and obesity. A large causal effect of pubertal timing independent of BMI is unlikely.

Exploring rare and low-frequency variants in the Saguenay-Lac-Saint-Jean population identified genes associated with asthma and allergy traits.

The Saguenay-Lac-Saint-Jean (SLSJ) region is located in northeastern Quebec and is known for its unique demographic history and founder effect. As founder populations are enriched with population-specific variants, we characterized the variants distribution in SLSJ and compared it with four European populations (Finnish, Sweden, United Kingdom and France), of which the Finnish population is another founder population. Targeted sequencing of the coding and non-coding immune regulatory regions of the SLSJ asthma familial cohort and the four European populations were performed. Rare and low-frequency coding and non-coding regulatory variants identified in the SLSJ population were then investigated for variant- and gene-level associations with asthma and allergy-related traits (eosinophil percentage, immunoglobulin (Ig) E levels and lung function). Our data showed that (1) rare or deleterious variants were not enriched in the two founder populations as compared with the three non-founder European populations; (2) a larger proportion of founder population-specific variants occurred with higher frequencies; and (3) low-frequency variants appeared to be more deleterious. Furthermore, a rare variant, rs1386931, located in the 3'-UTR of CXCR6 and intron of FYCO1 was found to be associated with eosinophil percentage. Gene-based analyses identified NRP2, MRPL44 and SERPINE2 to be associated with various asthma and allergy-related traits. Our study demonstrated the usefulness of using a founder population to identify new genes associated with asthma and allergy-related traits; thus better understand the genes and pathways implicated in pathophysiology.

Multiple sclerosis risk variants alter expression of co-stimulatory genes in B cells.

The increasing evidence supporting a role for B cells in the pathogenesis of multiple sclerosis prompted us to investigate the influence of known susceptibility variants on the surface expression of co-stimulatory molecules in these cells. Using flow cytometry we measured surface expression of CD40 and CD86 in B cells from 68 patients and 162 healthy controls that were genotyped for the multiple sclerosis associated single nucleotide polymorphisms (SNPs) rs4810485, which maps within the CD40 gene, and rs9282641, which maps within the CD86 gene. We found that carrying the risk allele rs4810485*T lowered the cell-surface expression of CD40 in all tested B cell subtypes (in total B cells P ≤ 5.10 × 10-5 in patients and ≤4.09 × 10-6 in controls), while carrying the risk allele rs9282641*G increased the expression of CD86, with this effect primarily seen in the naïve B cell subset (P = 0.048 in patients and 5.38 × 10-5 in controls). In concordance with these results, analysis of RNA expression demonstrated that the risk allele rs4810485*T resulted in lower total CD40 expression (P = 0.057) but with an increased proportion of alternative splice-forms leading to decoy receptors (P = 4.00 × 10-7). Finally, we also observed that the risk allele rs4810485*T was associated with decreased levels of interleukin-10 (P = 0.020), which is considered to have an immunoregulatory function downstream of CD40. Given the importance of these co-stimulatory molecules in determining the immune reaction that appears in response to antigen our data suggest that B cells might have an important antigen presentation and immunoregulatory role in the pathogenesis of multiple sclerosis.

Population- and individual-specific regulatory variation in Sardinia.

Genetic studies of complex traits have mainly identified associations with noncoding variants. To further determine the contribution of regulatory variation, we combined whole-genome and transcriptome data for 624 individuals from Sardinia to identify common and rare variants that influence gene expression and splicing. We identified 21,183 expression quantitative trait loci (eQTLs) and 6,768 splicing quantitative trait loci (sQTLs), including 619 new QTLs. We identified high-frequency QTLs and found evidence of selection near genes involved in malarial resistance and increased multiple sclerosis risk, reflecting the epidemiological history of Sardinia. Using family relationships, we identified 809 segregating expression outliers (median z score of 2.97), averaging 13.3 genes per individual. Outlier genes were enriched for proximal rare variants, providing a new approach to study large-effect regulatory variants and their relevance to traits. Our results provide insight into the effects of regulatory variants and their relationship to population history and individual genetic risk.

Age Related Multiple Sclerosis Severity Score: Disability ranked by age.

BACKGROUND: The Multiple Sclerosis Severity Score (MSSS) is obtained by normalising the Expanded Disability Status Scale (EDSS) score for disease duration and has been a valuable tool in cross-sectional studies. OBJECTIVE: To assess whether use of age rather than the inherently ambiguous disease duration was a feasible approach. METHOD: We pooled disability data from three population-based cohorts and developed an Age Related Multiple Sclerosis Severity (ARMSS) score by ranking EDSS scores based on the patient's age at the time of assessment. We established the power to detect a difference between groups afforded by the ARMSS score and assessed its relative consistency over time. RESULTS: The study population included 26058 patients from Sweden ( n = 11846), Canada ( n = 6179) and the United Kingdom ( n = 8033). There was a moderate correlation between EDSS and disease duration ( r = 0.46, 95% confidence interval (CI): 0.45-0.47) and between EDSS and age ( r = 0.44, 95% CI: 0.43-0.45). The ARMSS scores showed comparable power to detect disability differences between groups to the updated and original MSSS. CONCLUSION: Since age is typically unbiased and readily obtained, and the ARMSS and MSSS were comparable, the ARMSS may provide a more versatile tool and could minimise study biases and loss of statistical power caused by inaccurate or missing onset dates.

Large differences in adiponectin levels have no clear effect on multiple sclerosis risk: A Mendelian randomization study.

BACKGROUND: Mendelian randomization (MR) studies have demonstrated strong support for an association between genetically increased body mass index and risk of multiple sclerosis (MS). The adipokine adiponectin may be a potential mechanism linking body mass to risk of MS. OBJECTIVE: To evaluate whether genetically increased adiponectin levels influence risk of MS. METHODS: Using genome-wide significant single nucleotide polymorphisms (SNPs) for adiponectin, we undertook an MR study to estimate the effect of adiponectin on MS. This method prevents bias due to reverse causation and minimizes bias due to confounding. Sensitivity analyses were performed to evaluate the assumptions of MR. RESULTS: MR analyses did not support a role for genetically elevated adiponectin in risk of MS (odds ratio (OR) = 0.93 per unit increase in natural-log-transformed adiponectin, equivalent to a two-standard deviation increase in adiponectin on the absolute scale; 95% confidence interval (CI) = 0.66-1.33; p = 0.61). Further MR analysis suggested that genetic variation at the adiponectin gene, which influences adiponectin level, does not impact MS risk. Sensitivity analyses, including MR-Egger regression, suggested no bias due to pleiotropy. CONCLUSION: Lifelong genetically increased adiponectin levels in humans have no clear effect on risk of MS. Other biological factors driving the association between body mass and MS should be investigated.

Resolving TYK2 locus genotype-to-phenotype differences in autoimmunity.

Thousands of genetic variants have been identified, which contribute to the development of complex diseases, but determining how to elucidate their biological consequences for translation into clinical benefit is challenging. Conflicting evidence regarding the functional impact of genetic variants in the tyrosine kinase 2 (TYK2) gene, which is differentially associated with common autoimmune diseases, currently obscures the potential of TYK2 as a therapeutic target. We aimed to resolve this conflict by performing genetic meta-analysis across disorders; subsequent molecular, cellular, in vivo, and structural functional follow-up; and epidemiological studies. Our data revealed a protective homozygous effect that defined a signaling optimum between autoimmunity and immunodeficiency and identified TYK2 as a potential drug target for certain common autoimmune disorders.

Immunoseq: the identification of functionally relevant variants through targeted capture and sequencing of active regulatory regions in human immune cells.

BACKGROUND: The observation that the genetic variants identified in genome-wide association studies (GWAS) frequently lie in non-coding regions of the genome that contain cis-regulatory elements suggests that altered gene expression underlies the development of many complex traits. In order to efficiently make a comprehensive assessment of the impact of non-coding genetic variation in immune related diseases we emulated the whole-exome sequencing paradigm and developed a custom capture panel for the known DNase I hypersensitive site (DHS) in immune cells - "Immunoseq". RESULTS: We performed Immunoseq in 30 healthy individuals where we had existing transcriptome data from T cells. We identified a large number of novel non-coding variants in these samples. Relying on allele specific expression measurements, we also showed that our selected capture regions are enriched for functional variants that have an impact on differential allelic gene expression. The results from a replication set with 180 samples confirmed our observations. CONCLUSIONS: We show that Immunoseq is a powerful approach to detect novel rare variants in regulatory regions. We also demonstrate that these novel variants have a potential functional role in immune cells.

Increased THEMIS First Exon Usage in CD4+ T-Cells Is Associated with a Genotype that Is Protective against Multiple Sclerosis.

Multiple sclerosis is an autoimmune disease of the central nervous system. Genome wide association studies have identified over 100 common variants associated with multiple sclerosis, the majority of which implicate immunologically relevant genes, particularly those involved in T-cell development. SNP rs13204742 at the THEMIS/PTPRK locus is one such variant. Here, we have demonstrated mutually exclusive use of exon 1 and 2 amongst 16 novel THEMIS isoforms. We also show inverse correlation between THEMIS expression in human CD4+ T-cells and dosage of the multiple sclerosis risk allele at rs13204742, driven by reduced expression of exon 1- containing isoforms. In silico analysis suggests that this may be due to cell-specific, allele-dependent binding of the transcription factors FoxP3 and/or E47. Research exploring the functional implications of GWAS variants is important for gaining an understanding of disease pathogenesis, with the ultimate aim of identifying new therapeutic targets.

Obesity and Multiple Sclerosis: A Mendelian Randomization Study.

BACKGROUND: Observational studies have reported an association between obesity, as measured by elevated body mass index (BMI), in early adulthood and risk of multiple sclerosis (MS). However, bias potentially introduced by confounding and reverse causation may have influenced these findings. Therefore, we elected to perform Mendelian randomization (MR) analyses to evaluate whether genetically increased BMI is associated with an increased risk of MS. METHODS AND FINDINGS: Employing a two-sample MR approach, we used summary statistics from the Genetic Investigation of Anthropometric Traits (GIANT) consortium and the International MS Genetics Consortium (IMSGC), the largest genome-wide association studies for BMI and MS, respectively (GIANT: n = 322,105; IMSGC: n = 14,498 cases and 24,091 controls). Seventy single nucleotide polymorphisms (SNPs) were genome-wide significant (p < 5 x 10-8) for BMI in GIANT (n = 322,105) and were investigated for their association with MS risk in the IMSGC. The effect of each SNP on MS was weighted by its effect on BMI, and estimates were pooled to provide a summary measure for the effect of increased BMI upon risk of MS. Our results suggest that increased BMI influences MS susceptibility, where a 1 standard deviation increase in genetically determined BMI (kg/m2) increased odds of MS by 41% (odds ratio [OR]: 1.41, 95% CI 1.20-1.66, p = 2.7 x 10-5, I2 = 0%, 95% CI 0-29). Sensitivity analyses, including MR-Egger regression, and the weighted median approach provided no evidence of pleiotropic effects. The main study limitations are that, while these sensitivity analyses reduce the possibility that pleiotropy influenced our results, residual pleiotropy is difficult to exclude entirely. CONCLUSION: Genetically elevated BMI is associated with risk of MS, providing evidence for a causal role for obesity in MS etiology. While obesity has been associated with many late-life outcomes, these findings suggest an important consequence of childhood and/or early adulthood obesity.