- Research
- Open access
- Published:
Causal relationship between multiple sclerosis and cortical structure: a Mendelian randomization study
Journal of Translational Medicine volume 22, Article number: 83 (2024)
Abstract
Background
Observational studies have suggested an association between multiple sclerosis (MS) and cortical structure, but the results have been inconsistent.
Objective
We used two-sample Mendelian randomization (MR) to assess the causal relationship between MS and cortical structure.
Methods
MS data as the exposure trait, including 14,498 cases and 24,091 controls, were obtained from the International Multiple Sclerosis Genetics Consortium. Genome-wide association study (GWAS) data for cortical surface area (SAw/nw) and thickness (THw/nw) in 51,665 individuals of European ancestry were obtained from the ENIGMA Consortium. The inverse-variance weighted (IVW) method was used as the primary analysis for MR. Sensitivity analyses were conducted to evaluate heterogeneity and pleiotropy. Enrichment analysis was performed on MR analyses filtered by sensitivity analysis.
Results
After IVW and sensitivity analysis filtering, only six surviving MR results provided suggestive evidence supporting a causal relationship between MS and cortical structure, including lingual SAw (p = .0342, beta (se) = 5.7127 (2.6969)), parahippocampal SAw (p = .0224, beta (se) = 1.5577 (0.6822)), rostral middle frontal SAw (p = .0154, beta (se) = − 9.0301 (3.7281)), cuneus THw (p = .0418, beta (se) = − 0.0020 (0.0010)), lateral orbitofrontal THw (p = .0281, beta (se) = 0.0025 (0.0010)), and lateral orbitofrontal THnw (p = .0417, beta (se) = 0.0029 (0.0014)). Enrichment analysis suggested that leukocyte cell-related pathways, JAK-STAT signaling pathway, NF-kappa B signaling pathway, cytokine-cytokine receptor interaction, and prolactin signaling pathway may be involved in the effect of MS on cortical morphology.
Conclusion
Our results provide evidence supporting a causal relationship between MS and cortical structure. Enrichment analysis suggests that the pathways mediating brain morphology abnormalities in MS patients are mainly related to immune and inflammation-driven pathways.
Introduction
Multiple sclerosis (MS) is a multifactorial central nervous system autoimmune disease that leads to demyelination and axonal loss, with a male-to-female ratio of 1:3 [1]. MS affects over 2 million people, imposing a significant global health burden [2]. Some studies have reported that approximately 40–70% of MS patients experience cognitive dysfunction [3, 4], and up to 60% have neuropsychiatric symptoms [5, 6]. These neuropsychiatric abnormalities may be related to changes in brain connectivity. Accumulating evidence suggests that gray matter structures are widely affected in MS patients [7,8,9]. A 1-year follow-up study of 75 relapsing–remitting MS (RRMS) and 11 progressive MS (PMS) patients showed that RRMS patients had thinning of the frontal and temporal cortices, while PMS patients had a general decrease in cortical thickness (TH) [10]. Several other studies have confirmed that RRMS patients exhibit a decrease in average TH across the entire brain [11,12,13]. Nygaard et al. found that RRMS patients had similar cortical surface area (SA) to healthy controls [14], but another study provided evidence supporting a decrease in central anterior SA in MS patients [15].
However, these reports of the association between MS and brain structure are observational, and the conclusions are not consistent. Moreover, these studies are vulnerable to confounding factors and reverse causality, making causal inference difficult using traditional observational epidemiology. More importantly, the extent to which MS affects cortical structures and the underlying mechanisms remain unclear. Therefore, a tailored approach to exploring the causal relationship between MS and the cortical brain structure is crucial.
Mendelian randomization (MR) is a widely used approach in the study of neurological disorders, which employs single nucleotide polymorphisms (SNPs) as instrumental variables (IVs) to explore causal relationships between exposures and outcomes [16,17,18,19]. MR minimizes confounding and avoids reverse causation [17]. The fundamental principle in mitigating confounding factors within MR studies is rooted in the random allocation of genetic variations during meiosis, ensuring that one trait is typically unrelated to others. This method serves to avoid reverse causation, given that the genetic variations utilized to represent the impact of exposure remain unaltered by the occurrence and progression of outcomes [17]. There are three key assumptions in MR analysis. First, the genetic variation used as IVs should be associated with the risk factor of interest. Second, the genetic variation used should not be associated with potential confounders. Third, the selected genetic variation should only affect the outcome through the risk factor [20]. Here, we conducted a two-sample MR study using publicly available genome-wide association study (GWAS) data to determine the causal relationship between MS and brain structure.
Methods
Data sources and genetic instruments
Participants
We obtained MS GWAS data from the International Multiple Sclerosis Genetics Consortium (IMSGC). The dataset comprises 38,589 individuals of European ancestry, including 14,498 cases and 24,091 healthy controls. Disease diagnosis was conducted by neurologists familiar with MS, following widely recognized diagnostic criteria [21,22,23]. Disease severity was assessed using the Expanded Disability Status Score (EDSS) [24] and the Multiple Sclerosis Severity Score (MSSS) [25]. The overall age of onset for MS is 33.1 years. The study identified 48 new susceptibility loci for MS, highlighting the role of NF-kappa B in the disease's pathobiology. In addition, the overlap between MS and other autoimmune diseases was calculated using an immune chip analysis. Approximately 22% of the MS signals overlapped with at least one other autoimmune disease, primarily including inflammatory bowel disease (~ 9.1%), primary biliary cirrhosis (~ 9.1%), Crohn’s disease (~ 9.1%), ulcerative colitis (~ 7.3%), celiac disease (~ 4.5%), rheumatoid arthritis (~ 4.5%), and autoimmune thyroid disease (~ 2.7%) [26]. For further details on this study, please refer to the original publication [26].
Brain cortical structure
The summary-level data on the cortical structure is derived from the work of the Enhancing NeuroImaging Genetics through Meta-Analysis Consortium (ENIGMA)-Genetics working group. Grasby et al. utilized magnetic resonance imaging data from 51,665 individuals to investigate the cortical SA and TH [27]. The GWAS meta-analysis identified loci that affect regional SA near genes involved in the Wnt signaling pathway. Specifically, the study divided cortical structures into 34 functionally specialized regions based on the Desikan-Killiany atlas, and regional SA and TH were analyzed with and without whole-brain weighting (w indicates weighted regions, while nw indicates unweighted regions) [27]. Additional details can be found in the original research study [27]. Therefore, we utilized this GWAS data to analyze the causal effects of MS on both whole-brain and the 34 functional regions' cortical SA and TH, resulting in a total of 138 analyses.
Selection of instrumental variables
Based on the basic assumptions of MR, we first extracted MS-related IVs with a threshold of P < 5E−8. Second, we clumped the IVs based on the European 1000 Genomes Project (clumped R2 < 0.001, window size = 10 Mb). Third, proxy genetic variants were used when no corresponding SNPs were available (LD R2 threshold was set at 0.8). Fourth, we homogenized the SNPs and removed palindromic SNPs. To ensure the IVs had strong statistical power, we selected SNPs with an F-statistic greater than 10. We used the formula R2 * (N−k−1)/[ (1−R2) * k] to calculate the F-statistic, where N is the sample size of MS, k is the number of SNPs, and R2 is the proportion of MS variation explained by each SNP. R2 was calculated by the formula: 2* beta2* (1-eaf) * eaf, where eaf is the allele frequency of the effect, and beta is the estimate of the genetic effect of the variant on MS [28, 29]. Our MR study followed the STROBE-MR Statement guidelines [30].
MR analysis and sensitivity analysis
The inverse-variance weighted (IVW) method with random effects is used as the main approach for MR analysis. The IVW approach combines the Wald ratio estimates of each SNP to obtain an overall estimate of the causal effect. Although IVW allows for heterogeneity [20, 31], it is susceptible to pleiotropic bias [32, 33]. MR-Egger, weighted median, simple mode, and weighted mode methods rely on assumptions different from those of IVW and are relatively robust to horizontal pleiotropy, which can complement IVW to make MR estimates more reliable [34]. Cochran's Q test is used to detect heterogeneity. The MR-Egger intercept test generates a non-zero intercept, indicating the presence of directional pleiotropy [35]. The MR Pleiotropy Residual Sum and Outlier (MR-PRESSO) test is used to detect potential outliers and provides corrected estimates of MR to account for horizontal pleiotropy after removing such outliers [36]. Leave-one-out analysis checks whether MR estimates are driven by individual SNPs. Additionally, we employed the Steiger test to mitigate the impact of reverse causation [37, 38]. For significant results determined by IVW (p < 0.05), we use the PhenoScanner online tool to search for the second phenotypes of genetic variants, including body mass index, obesity, smoking, drinking, neuropsychiatric disease, hypertension, and hyperlipemia, to evaluate whether these MR estimates are overturned by potential confounding factors [39, 40]. These tests have different assumptions but are useful for explaining other causal pathways apart from the hypothesized pathway [34]. Multiple comparisons are corrected by the Bonferroni method. The meaningful and nominally significant thresholds were set at p < 0.05/138 = 0.0004 and p < 0.05, respectively, out of 138 tests conducted. The analysis was performed using the TwoSampleMR package (version 0.5.6) in R software (version 4.2.1).
Causal genomic loci and enrichment analysis
To further explore the potential mechanisms underlying cortical structural changes mediated by MS, we examined gene sets with distinct causal effects (increasing or decreasing TH and SA) on cortical structure. These sets were identified through a series of sensitivity tests and filtered MR analysis results. We then performed Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis on these gene sets. The analysis was completed using the clusterProfiler package (version v4.4.4) and the biomaRt package (version 2.52.0) in R software (version 4.2.1).
Results
MR analysis and sensitivity analysis
Out of a total of 138 MR analyses, we used the IVW method to preliminarily establish nominal causal associations between MS and cortical structure in 9 brain regions. These regions include lingual SAw, parahippocampal SAw, postcentral SAw, rostral middle frontal SAw, lingual SAnw, cuneus THw, lateral orbitofrontal THw, superior temporal THw, and lateral orbitofrontal THnw (Figs. 1, 2).
Sensitivity analysis showed that the MR-Egger intercept was removed in the lingual SAnw due to the detection of directional pleiotropy (intercept = − 2.51, se = 1.19, p = 0.04). Directional pleiotropy was not detected in the remaining MR estimates (MR-Egger intercept pvalue > 0.05, Additional file 1: Table S1). Although heterogeneity was detected in the superior temporal THw (Cochran's Q = 68.52, p = 0.01), the random effects IVW method accepted the presence of heterogeneity. The MR-PRESSO analysis confirmed the presence of potential outliers in the superior temporal THw region, but their identification was not feasible (p = 0.004). However, no outliers were detected in the other MR results (MR-PRESSO pvalue > 0.05, Additional file 1: Table S1). The remaining four supplemental methods for MR estimation are shown in Table S2. Leave-one-out analysis showed that none of the MR estimates were driven by a single SNP (Additional file 1: Figs. S1–S9). All MR estimates passed the Steiger test, maintaining consistency with the previous results.
In addition, we further conducted sensitivity analyses using the PhenoScanner tool for the remaining MR estimates. We found that rs7923837 and rs11554159 were associated with body mass index, rs8070345 with frequency of alcohol consumption, rs11154801 with schizophrenia, and rs2857700 with Parkinson's disease, total cholesterol, and self-reported hypertension. After removing the SNPs associated with these risk factors, the MR estimates for postcentral SAw and superior temporal THw were no longer consistent with the previous ones (pvalue > 0.05, Additional file 1: Table S3).
Therefore, after a series of sensitivity analyses, only six surviving MR results suggested a causal relationship between MS and cortical structure, including lingual SAw (p = 0.0342, beta (se) = 5.7127 (2.6969)), parahippocampal SAw (p = 0.0224, beta (se) = 1.5577 (0.6822)), rostral middle frontal SAw (p = 0.0154, beta (se) = − 9.0301 (3.7281)), cuneus THw (p = 0.0418, beta (se) = − 0.0020 (0.0010)), lateral orbitofrontal THw (p = 0.0281, beta (se) = 0.0025 (0.0010)), and lateral orbitofrontal THnw (p = 0.0417, beta (se) = 0.0029 (0.0014)) (Fig. 2). It is noteworthy that these MR estimates were limited to regional functional levels, and there was no evidence to support a potential causal relationship between MS and global cortical TH or SA (Fig. 1).
Causal loci and enrichment analysis
Although MS has different effects on the TH and SA of the cerebral cortex, we found that the gene sets underlying the causal relationship between MS and brain morphology were completely overlapping. A total of 41 genes were implicated in the causal link between MS and cortical structure (Additional file 1: Table S4). Notably, 13 key genes including ELMO1, SOCS1, STAT3, STAT4, IL7R, IL22RA2, TNFRSF1A, TNFSF14, IL2RA, TNFAIP3, BCL10, MAPK3, CD86, CXCR5, and LTBR were found to regulate the changes in cortical structure caused by MS (Additional file 1: Fig. S10). GO analysis revealed that these pathways were mainly enriched in leukocyte proliferation, positive regulation of T cell activation, positive regulation of leukocyte cell–cell adhesion, and receptor signaling pathways via STAT (Fig. 3A). KEGG analysis showed that these pathways were mainly involved in JAK-STAT signaling, NF-kappa B signaling, cytokine-cytokine receptor interaction, and prolactin signaling pathways (Fig. 3B).
Discussion
To our knowledge, this is the first MR study investigating the causal relationship between MS and brain structure. We implemented a large-scale GWAS summary data analysis on MS and brain structure, ultimately discovering nominal associations between MS and 6 cortical structures, including decreased rostral middle frontal SAw, cuneus THw, and increased lingual SAw, parahippocampal SAw, lateral orbitofrontal THw, and lateral orbitofrontal THnw. Sensitivity analysis strengthened the robustness of our MR estimates.
Our MR analysis suggested that MS primarily affects cortical structures concentrated in the frontal and temporal lobes, partially overlapping with previous observational studies, although the effects were not entirely consistent. Previous studies have demonstrated atrophy in the precuneus and cuneus regions in MS patients [41, 42], and our findings indicated comparable declines in the TH of the cuneus region. Furthermore, we identified a causal association between MS and increased lingual SA. While most visual impairments in MS are due to optic neuritis [43], these regions are located in the occipital lobe and are closely related to vision, which may explain some of the mechanisms underlying visual impairments in MS patients. Unfortunately, the brain structure GWAS by Grasby et al. did not include the optic nerve [27], so we were unable to determine the causal effect of MS on the optic nerve in the current study.
Furthermore, we presented compelling evidence supporting a causal association between MS and increased SA in the parahippocampal gyrus, elevated TH in the lateral orbitofrontal cortex, and decreased SA in the rostral middle frontal region. Accumulating evidence has suggested a close association between the parahippocampal gyrus and olfactory and cognitive impairments [44,45,46], with reports of 20 to 45% of MS patients experiencing olfactory dysfunction [47]. Moreover, the parahippocampal gyrus is directly connected to multiple regions of the prefrontal cortex [44], and olfactory dysfunction may be associated with damage to the subfrontal and temporal lobes [48]. Steenwijk et al. observed that the pattern of cortical thickness changes in individuals with MS primarily involves bilateral temporal poles and insular cortex [9]. Several observational studies have reported cortical thinning in the frontal and temporal lobes of MS patients [10, 49,50,51]. Our MR estimates similarly suggest damage to the temporal and frontal lobes, especially a decrease in rostral middle frontal SA. Therefore, it is reasonable to speculate that damage to the olfactory pathway from the temporal lobe to the frontal cortex, especially the orbitofrontal cortex, may contribute to the olfactory decline in MS patients [46, 47]. It is worth noting that the parahippocampal cortex is associated with various cognitive processes, including visual spatial processing and episodic memory [44, 45]. Although our results suggest an association between MS and the parahippocampal cortex, further research is needed to determine whether MS patients experience neurobehavioral symptoms through the impact on the parahippocampal gyrus.
Our findings of increased cortical SA or TH in some regions, though unexpected, are supported by limited observational studies. A magnetic resonance imaging study of patients with clinically isolated syndrome suggestive of MS within 2 years of onset showed increased volume in frontal, parietal, temporal, and cerebellar gray matter regions at 3 months [52]. Another study focusing on mildly disabled patients with RRMS showed increased TH in two visual cortical areas, the left hemisphere's inferior occipital gyrus and the right hemisphere's cuneus [53]. The researchers speculated that cortical reorganization in adjacent brain regions and compensatory mechanisms recruited by increased cortical networks could reduce disease activity and even result in compensatory hypertrophy or brain edema, leading to increased cortical volume [40, 53]. Furthermore, unconventional cortical changes have been observed in other neurological diseases, such as larger cortical SA in early-stage Parkinson's disease patients compared to controls [54]. An ENIGMA study showed a correlation between smoking and increased brain structural changes [55]. Thus, MS may not be limited to the pattern of gray matter atrophy alone. However, the underlying mechanisms are not fully understood and require further clarification. Additionally, we noted inconsistent trends in cortical SA and TH morphology, which may be attributed to their different anatomical developmental patterns [40].
Research has confirmed that alterations in the function of regulatory T cells and/or effector B cells and T cells disrupt peripheral tolerance mechanisms in patients with MS, thereby promoting the generation of pro-inflammatory mediators such as cytokines. The communication between the peripheral immune system and the central nervous system, facilitated by messenger molecules like cytokines, leads to neuroinflammation and immune responses in the brain and spinal cord regions [56]. Through causal gene enrichment analysis, we identified pathways mediating cortical changes in MS that predominantly cluster in leukocyte cell-related pathways, cytokine-cytokine receptor interaction, JAK-STAT signaling pathway, and NF-kappa B signaling pathway. These pathways are intricately linked to various immune and inflammatory processes [56,57,58,59]. Simultaneously, gray matter damage is closely associated with peripheral-induced neuroinflammation and immune responses [60,61,62]. Given these findings, we speculate that the structural abnormalities in the brains of MS patients are primarily driven by immune and inflammatory pathways. Our findings suggest that the current disease-modifying treatment strategy [56], primarily focused on anti-inflammatory approaches, is well-founded. This is further supported by certain drugs, such as selective sphingosine 1-phosphate receptor subtype 1 [61]. Additionally, considering the successful application of Janus kinase inhibitors (Jakinibs) in the treatment of rheumatic diseases [59], we propose that Jakinibs could serve as candidate drugs for treating cortical changes in MS, pending further research to confirm this hypothesis.
The present study should be considered in light of its strengths and limitations. Our study has several strengths. Firstly, we employed a two-sample MR design, which minimizes confounding factors and reverses causation that is inevitable in observational studies. Secondly, we strengthened the statistical power by using a stringent P-value threshold and F-statistics. Thirdly, sensitivity analyses provided no evidence of horizontal pleiotropy, reinforcing the robustness of our MR estimates. However, our study has several limitations. Firstly, our study population was limited to individuals of European ancestry, which restricts the generalizability of our results to other populations. Future studies in other populations are warranted. Secondly, a large body of research has focused on RRMS and PMS patients, whereas our MR estimates were based on summary-level data and do not provide more specific subtype information. This implies that stratified analyses may lead to inconsistent conclusions. However, researchers have increasingly recognized that different phenotypes of MS may represent different courses of the same disease [63], making our MR analysis appropriate from that perspective. Thirdly, although the causality between MS and cortical structure is suggestive, Bonferroni correction is considered conservative [64].
In summary, we provide robust evidence supporting a causal relationship between genetic proxies of MS and cortical structure. Enrichment analysis suggests that the pathways mediating brain morphology abnormalities in MS patients are mainly related to immune and inflammation-driven pathways. Further research is needed to validate these findings.
Availability of data and materials
Data supporting the findings of this study are available from the article/additional material.
Abbreviations
- MS:
-
Multiple sclerosis
- RRMS:
-
Relapsing-remitting multiple sclerosis
- PMS:
-
Progressive multiple sclerosis
- TH:
-
Thickness
- SA:
-
Surface area
- MR:
-
Mendelian randomization
- SNPs:
-
Single nucleotide polymorphisms
- IVs:
-
Instrumental variables
- GWAS:
-
Genome-wide association study
- IMSGC:
-
International Multiple Sclerosis Genetics Consortium
- ENIGMA:
-
Enhancing NeuroImaging Genetics through Meta-Analysis Consortium
- IVW:
-
The inverse-variance weighted method with random effects
- MR-PRESSO:
-
MR Pleiotropy Residual Sum and Outlier test
- GO:
-
Gene Ontology
- KEGG:
-
Kyoto Encyclopedia of Genes and Genomes
References
Sparaco M, Bonavita S. The role of sex hormones in women with multiple sclerosis: from puberty to assisted reproductive techniques. Front Neuroendocrinol. 2021;60: 100889.
Feigin VL, et al. Global, regional, and national burden of neurological disorders during 1990–2015: a systematic analysis for the Global Burden of Disease Study 2015. Lancet Neurol. 2017;16(11):877–97.
Manca R, Sharrack B, Paling D, Wilkinson ID, Venneri A. Brain connectivity and cognitive processing speed in multiple sclerosis: a systematic review. J Neurol Sci. 2018;388:115–27.
DeLuca GC, Yates RL, Beale H, Morrow SA. Cognitive impairment in multiple sclerosis: clinical, radiologic and pathologic insights. Brain Pathol (Zurich, Switzerland). 2015;25(1):79–98.
Murphy R, O’Donoghue S, Counihan T, McDonald C, Calabresi PA, Ahmed MA, et al. Neuropsychiatric syndromes of multiple sclerosis. J Neurol Neurosurg Psychiatry. 2017;88(8):697–708.
Marrie RA, Horwitz R, Cutter G, Tyry T, Campagnolo D, Vollmer T. The burden of mental comorbidity in multiple sclerosis: frequent, underdiagnosed, and undertreated. Mult Scler (Houndmills, Basingstoke, England). 2009;15(3):385–92.
Colato E, Stutters J, Tur C, Narayanan S, Arnold DL, Gandini Wheeler-Kingshott CAM, et al. Predicting disability progression and cognitive worsening in multiple sclerosis using patterns of grey matter volumes. J Neurol Neurosurg Psychiatry. 2021;92(9):995–1006.
Cagol A, Schaedelin S, Barakovic M, Benkert P, Todea RA, Rahmanzadeh R, et al. Association of brain atrophy with disease progression independent of relapse activity in patients with relapsing multiple sclerosis. JAMA Neurol. 2022;79(7):682–92.
Steenwijk MD, Geurts JJ, Daams M, Tijms BM, Wink AM, Balk LJ, et al. Cortical atrophy patterns in multiple sclerosis are non-random and clinically relevant. Brain J Neurol. 2016;139(Pt 1):115–26.
Hidalgo de la Cruz M, Valsasina P, Gobbi C, Gallo A, Zecca C, Bisecco A, et al. Longitudinal cortical thinning progression differs across multiple sclerosis phenotypes and is clinically relevant: a multicentre study. Mult Scler (Houndmills, Basingstoke, England). 2021;27(6):827–40.
Tillema JM, Hulst HE, Rocca MA, Vrenken H, Steenwijk MD, Damjanovic D, et al. Regional cortical thinning in multiple sclerosis and its relation with cognitive impairment: a multicenter study. Mult Scler (Houndmills, Basingstoke, England). 2016;22(7):901–9.
Orbach L, Menascu S, Hoffmann C, Miron S, Achiron A. Focal cortical thinning in patients with stable relapsing-remitting multiple sclerosis: cross-sectional-based novel estimation of gray matter kinetics. Neuroradiology. 2018;60(2):179–87.
Geisseler O, Pflugshaupt T, Bezzola L, Reuter K, Weller D, Schuknecht B, et al. Cortical thinning in the anterior cingulate cortex predicts multiple sclerosis patients’ fluency performance in a lateralised manner. NeuroImage Clin. 2016;10:89–95.
Nygaard GO, Walhovd KB, Sowa P, Chepkoech JL, Bjørnerud A, Due-Tønnessen P, et al. Cortical thickness and surface area relate to specific symptoms in early relapsing-remitting multiple sclerosis. Mult Scler (Houndmills, Basingstoke, England). 2015;21(4):402–14.
Gorgoraptis N, Wheeler-Kingshott CA, Jenkins TM, Altmann DR, Miller DH, Thompson AJ, et al. Combining tractography and cortical measures to test system-specific hypotheses in multiple sclerosis. Mult Scler (Houndmills, Basingstoke, England). 2010;16(5):555–65.
Zhu G, Zhou S, Xu Y, Gao R, Zhang M, Zeng Q, et al. Chickenpox and multiple sclerosis: a Mendelian randomization study. J Med Virol. 2023;95(1): e28315.
Yuan S, Xiong Y, Larsson SC. An atlas on risk factors for multiple sclerosis: a Mendelian randomization study. J Neurol. 2021;268(1):114–24.
Sun D, Wang R, Shi Z, Chen H, Zhou H. Causal relationship of COVID-19 on Alzheimer’s disease: an updated Mendelian randomization study. Neurol Sci. 2023. https://doi.org/10.1007/s10072-023-06807-8.
Wang X, Shi Z, Zhao Z, Chen H, Lang Y, Kong L, et al. The causal relationship between neuromyelitis optica spectrum disorder and other autoimmune diseases. Front Immunol. 2022;13: 959469.
Burgess S, Butterworth A, Thompson SG. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol. 2013;37(7):658–65.
Poser CM, Paty DW, Scheinberg L, McDonald WI, Davis FA, Ebers GC, et al. New diagnostic criteria for multiple sclerosis: guidelines for research protocols. Ann Neurol. 1983;13(3):227–31.
McDonald WI, Compston A, Edan G, Goodkin D, Hartung HP, Lublin FD, et al. Recommended diagnostic criteria for multiple sclerosis: guidelines from the International Panel on the diagnosis of multiple sclerosis. Ann Neurol. 2001;50(1):121–7.
Polman CH, Reingold SC, Edan G, Filippi M, Hartung HP, Kappos L, et al. Diagnostic criteria for multiple sclerosis: 2005 revisions to the “McDonald Criteria.” Ann Neurol. 2005;58(6):840–6.
Kurtzke JF. Rating neurologic impairment in multiple sclerosis: an expanded disability status scale (EDSS). Neurology. 1983;33(11):1444–52.
Roxburgh RH, Seaman SR, Masterman T, Hensiek AE, Sawcer SJ, Vukusic S, et al. Multiple sclerosis severity score: using disability and disease duration to rate disease severity. Neurology. 2005;64(7):1144–51.
Beecham AH, Patsopoulos NA, Xifara DK, Davis MF, Kemppinen A, Cotsapas C, et al. Analysis of immune-related loci identifies 48 new susceptibility variants for multiple sclerosis. Nat Genet. 2013;45(11):1353–60.
Grasby KL, Jahanshad N, Painter JN, Colodro-Conde L, Bralten J, Hibar DP, et al. The genetic architecture of the human cerebral cortex. Science (New York NY). 2020. https://doi.org/10.1126/science.aay6690.
Burgess S, Thompson SG. Avoiding bias from weak instruments in Mendelian randomization studies. Int J Epidemiol. 2011;40(3):755–64.
Pierce BL, Ahsan H, Vanderweele TJ. Power and instrument strength requirements for Mendelian randomization studies using multiple genetic variants. Int J Epidemiol. 2011;40(3):740–52.
Skrivankova VW, Richmond RC, Woolf BAR, Yarmolinsky J, Davies NM, Swanson SA, et al. Strengthening the reporting of observational studies in epidemiology using Mendelian randomization: The STROBE-MR Statement. JAMA. 2021;326(16):1614–21.
Burgess S, Davey Smith G, Davies NM, Dudbridge F, Gill D, Glymour MM, et al. Guidelines for performing Mendelian randomization investigations. Wellcome open research. 2019;4:186.
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–42.
Cai J, Li X, Wu S, Tian Y, Zhang Y, Wei Z, et al. Assessing the causal association between human blood metabolites and the risk of epilepsy. J Transl Med. 2022;20(1):437.
Choi KW, Chen CY, Stein MB, Klimentidis YC, Wang MJ, Koenen KC, et al. Assessment of bidirectional relationships between physical activity and depression among adults: a 2-sample Mendelian randomization study. JAMA Psychiat. 2019;76(4):399–408.
Williams JA, Burgess S, Suckling J, Lalousis PA, Batool F, Griffiths SL, et al. Inflammation and brain structure in schizophrenia and other neuropsychiatric disorders: a Mendelian randomization study. JAMA Psychiat. 2022;79(5):498–507.
Verbanck M, Chen CY, 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.
Hemani G, Tilling K, Davey SG. Orienting the causal relationship between imprecisely measured traits using GWAS summary data. PLoS Genet. 2017;13(11): e1007081.
Lin J, Zhou J, Xu Y. Potential drug targets for multiple sclerosis identified through Mendelian randomization analysis. Brain J Neurol. 2023;146(8):3364–72.
Kamat MA, Blackshaw JA, Young R, Surendran P, Burgess S, Danesh J, et al. PhenoScanner V2: an expanded tool for searching human genotype-phenotype associations. Bioinformatics (Oxford, England). 2019;35(22):4851–3.
Chen X, Kong J, Pan J, Huang K, Zhou W, Diao X, et al. Kidney damage causally affects the brain cortical structure: A Mendelian randomization study. EBioMedicine. 2021;72: 103592.
Eshaghi A, Marinescu RV, Young AL, Firth NC, Prados F, Jorge Cardoso M, et al. Progression of regional grey matter atrophy in multiple sclerosis. Brain J Neurol. 2018;141(6):1665–77.
Solana E, Martinez-Heras E, Montal V, Vilaplana E, Lopez-Soley E, Radua J, et al. Regional grey matter microstructural changes and volume loss according to disease duration in multiple sclerosis patients. Sci Rep. 2021;11(1):16805.
Dhanapalaratnam R, Markoulli M, Krishnan AV. Disorders of vision in multiple sclerosis. Clin Exp Optom. 2022;105(1):3–12.
Peng X, Burwell RD. Beyond the hippocampus: the role of parahippocampal-prefrontal communication in context-modulated behavior. Neurobiol Learn Mem. 2021;185: 107520.
Aminoff EM, Kveraga K, Bar M. The role of the parahippocampal cortex in cognition. Trends Cogn Sci. 2013;17(8):379–90.
Douaud G, Lee S, Alfaro-Almagro F, Arthofer C, Wang C, McCarthy P, et al. SARS-CoV-2 is associated with changes in brain structure in UK Biobank. Nature. 2022;604(7907):697–707.
Lucassen EB, Turel A, Knehans A, Huang X, Eslinger P. Olfactory dysfunction in multiple sclerosis: a scoping review of the literature. Mult Scler Relat Disord. 2016;6:1–9.
Zorzon M, Ukmar M, Bragadin LM, Zanier F, Antonello RM, Cazzato G, et al. Olfactory dysfunction and extent of white matter abnormalities in multiple sclerosis: a clinical and MR study. Mult Scler (Houndmills, Basingstoke, England). 2000;6(6):386–90.
Cruz-Gomez ÁJ, Forero L, Lozano-Soto E, Cano-Cano F, Sanmartino F, Rashid-López R, et al. Cortical thickness and serum NfL explain cognitive dysfunction in newly diagnosed patients with multiple sclerosis. Neurol Neuroimmunol Neuroinflamm. 2021. https://doi.org/10.1212/NXI.0000000000001074.
Pareto D, Sastre-Garriga J, Auger C, Vives-Gilabert Y, Delgado J, Tintoré M, et al. Juxtacortical lesions and cortical thinning in multiple sclerosis. AJNR Am J Neuroradiol. 2015;36(12):2270–6.
Treaba CA, Herranz E, Barletta VT, Mehndiratta A, Ouellette R, Sloane JA, et al. The relevance of multiple sclerosis cortical lesions on cortical thinning and their clinical impact as assessed by 70-T MRI. J Neurol. 2021;268(7):2473–81.
Rocca MA, Preziosa P, Mesaros S, Pagani E, Dackovic J, Stosic-Opincal T, et al. Clinically isolated syndrome suggestive of multiple sclerosis: dynamic patterns of gray and white matter changes—a 2-year MR imaging study. Radiology. 2016;278(3):841–53.
Achiron A, Chapman J, Tal S, Bercovich E, Gil H, Achiron A. Superior temporal gyrus thickness correlates with cognitive performance in multiple sclerosis. Brain Struct Funct. 2013;218(4):943–50.
Jubault T, Gagnon JF, Karama S, Ptito A, Lafontaine AL, Evans AC, et al. Patterns of cortical thickness and surface area in early Parkinson’s disease. Neuroimage. 2011;55(2):462–7.
Chye Y, Mackey S, Gutman BA, Ching CRK, Batalla A, Blaine S, et al. Subcortical surface morphometry in substance dependence: an ENIGMA addiction working group study. Addict Biol. 2020;25(6): e12830.
Dendrou CA, Fugger L, Friese MA. Immunopathology of multiple sclerosis. Nat Rev Immunol. 2015;15(9):545–58.
Barnabei L, Laplantine E, Mbongo W, Rieux-Laucat F, Weil R. NF-κB: at the borders of autoimmunity and inflammation. Front Immunol. 2021;12: 716469.
Hu X, Li J, Fu M, Zhao X, Wang W. The JAK/STAT signaling pathway: from bench to clinic. Signal Transduct Target Ther. 2021;6(1):402.
Banerjee S, Biehl A, Gadina M, Hasni S, Schwartz DM. JAK-STAT signaling as a target for inflammatory and autoimmune diseases: current and future prospects. Drugs. 2017;77(5):521–46.
Mittli D. Inflammatory processes in the prefrontal cortex induced by systemic immune challenge: focusing on neurons. Brain Behav Immunity Health. 2023;34: 100703.
Cohan SL, Benedict RHB, Cree BAC, DeLuca J, Hua LH, Chun J. The two sides of siponimod: evidence for brain and immune mechanisms in multiple sclerosis. CNS Drugs. 2022;36(7):703–19.
Zuroff LR, Benjamins JA, Bar-Or A, Lisak RPJN. Inflammatory mechanisms underlying cortical injury in progressive multiple sclerosis. Neuroimmunol Neuroinflamm. 2021;8:111.
Kuhlmann T, Moccia M, Coetzee T, Cohen JA, Correale J, Graves J, et al. Multiple sclerosis progression: time for a new mechanism-driven framework. Lancet Neurol. 2023;22(1):78–88.
Zeng R, Wang J, Jiang R, Yang J, Zheng C, Huihuan W, et al. Investigating causality and shared genetic architecture between neurodegenerative disorders and inflammatory bowel disease. Aging Dis. 2022. https://doi.org/10.14336/AD.2022.12209.
Acknowledgements
The authors express their gratitude for the generous sharing of GWAS summary statistics provided by the International Multiple Sclerosis Genetics Consortium (IMSGC) and the ENIGMA Consortium. We would also like to extend our appreciation to the ChipPlot tool (https://www.chiplot.online/) and the online analysis and drawing tool (http://bioinformatics.com.cn/) for their invaluable contributions to image processing. Furthermore, we would like to acknowledge the assistance of Dr. Qiuxing Lin from the Department of Neurology at West China Hospital of Sichuan University for his expertise in visualizing the images.
Funding
This study was supported by grants from the Department of Science and Technology of Sichuan Province (2022YFS0315 to ZHY), the 1·3·5 project for disciplines of excellence – Clinical Research Incubation Project, West China Hospital, Sichuan University (21HXFH041 to ZHY), the National Natural Science Foundation of China (82201494 to SZY), and the Natural Science Foundation of Sichuan Province (2022NSFSC1432 to WXF).
Author information
Authors and Affiliations
Contributions
Conceptualization and design: HZ, XW; Data collection and integration: DS, RW, QD, YZ, HC, ZS; Data analysis and interpretation: DS, RW, QD, YZ, HC, ZS, XW, HZ; Writing of the first draft of the manuscript: DS; Manuscript review and editing: DS, RW, QD, YZ, HC, ZS, XW, HZ. All authors reviewed and agreed to the final version of the manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
Not applicable. The present MR analysis was based on summary data from previous genome-wide association studies that had achieved written informed consent and ethics approval.
Competing interests
The authors declare 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: Table S1.
Results of sensitivity analysis. Table S2. Results of MR-Egger,weighted median,simple mode,weighted mode methods. Table S3. IVW results after removing confounders. Table S4. Nearest genes from causal SNPs. Figure S1. Leave-one-out analysis of lingual SAw. Figure S2. Leave-one-out analysis of parahippocampal SAw. Figure S3. Leave-one-out analysis of postcentral SAw. Figure S4. Leave-one-out analysis of rostral middle frontal SAw. Figure S5. Leave-one-out analysis of lingual SAnw. Figure S6. Leave-one-out analysis of cuneus THw. Figure S7. Leave-one-out analysis of lateral orbitofrontal THw. Figure S8. Leave-one-out analysis of superior temporal THw. Figure S9. Leave-one-out analysis of lateral orbitofrontal THnw. Figure S10. Network diagram for pathway analysis.
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 http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Sun, D., Wang, R., Du, Q. et al. Causal relationship between multiple sclerosis and cortical structure: a Mendelian randomization study. J Transl Med 22, 83 (2024). https://doi.org/10.1186/s12967-024-04892-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12967-024-04892-7