DNA Methylation in the Anti-Mullerian Hormone Gene and the Risk of Disease Activity in Multiple Sclerosis.

OBJECTIVE: Multiple sclerosis (MS) has a complex pathobiology, with genetic and environmental factors being crucial players. Understanding the mechanisms underlying heterogeneity in disease activity is crucial for tailored treatment. We explored the impact of DNA methylation, a key mechanism in the genetics-environment interplay, on disease activity in MS. METHODS: Peripheral immune methylome profiling using Illumina Infinium MethylationEPIC BeadChips was conducted on 249 untreated relapsing-remitting MS patients, sampled at the start of disease-modifying treatment (DMT). A differential methylation analysis compared patients with evidence of disease activity (EDA) to those with no evidence of disease activity (NEDA) over 2 years from DMT start. Utilizing causal inference testing (CIT) and Mendelian randomization (MR), we sought to elucidate the relationships between DNA methylation, gene expression, genetic variation, and disease activity. RESULTS: Four differentially methylated regions (DMRs) were identified between EDA and NEDA. Examining the influence of single nucleotide polymorphisms (SNPs), 923 variants were found to account for the observed differences in the 4 DMRs. Importantly, 3 out of the 923 SNPs, affecting DNA methylation in a DMR linked to the anti-Mullerian hormone (AMH) gene, were associated with disease activity risk in an independent cohort of 1,408 MS patients. CIT and MR demonstrated that DNA methylation in AMH acts as a mediator for the genetic risk of disease activity. INTERPRETATION: This study uncovered a novel molecular pathway implicating the interaction between DNA methylation and genetic variation in the risk of disease activity in MS, emphasizing the role of sex hormones, particularly the AMH, in MS pathobiology. ANN NEUROL 2024;96:289-301.

Vitamin D affects the risk of disease activity in multiple sclerosis.

BACKGROUND: Vitamin D (VitD) affects the risk of multiple sclerosis (MS), but the impact on disease activity is controversial. We assessed whether VitD is associated with the No-Evidence of Disease Activity-3 (NEDA-3) status at 2 years from disease-modifying treatment (DMT) start, and whether this association is causal or the result of confounding factors. Furthermore, we explored if a genetic predisposition to higher VitD levels affects the risk of disease activity. METHODS: 230 untreated relapsing-remitting MS patients underwent serum 25-OH-vitamin-D measurement, and the association between seasonally adjusted VitD and disease activity was tested. Modelling a Polygenic Risk Score from a Genome-Wide Association Study on ~400 000 individuals, we studied the impact of genetic predisposition to higher VitD on the NEDA-3 status in 1408 independent MS patients. Two-sample Mendelian randomisation (MR) was used to assess causality. RESULTS: Lower baseline VitD was associated with decreased probability of NEDA-3 at 2 years (p=0.019). Particularly, VitD levels <20 ng/mL conferred an over twofold risk of disease activity (OR 2.36, 95% CI 1.30 to 3.88, p=0.0037). Genetic predisposition to higher VitD levels was associated with delayed age at MS onset (p=0.018) and with a higher probability of NEDA-3 status (p=0.034). MR confirmed causality between VitD and the risk of disease activity (p=0.041). CONCLUSIONS: VitD levels before DMT start affect the risk of disease activity in MS. Genetic predisposition to higher VitD levels confers a lower risk of disease activity and is associated with delayed MS onset. Our work prompts future prospective studies regarding VitD supplementation and lifestyle interventions to hamper disease activity in MS.

Broadening the Genetic Spectrum of Painful Small-Fiber Neuropathy through Whole-Exome Study in Early-Onset Cases.

Small-Fiber Neuropathy (SFN) is a disorder of the peripheral nervous system, characterised by neuropathic pain; approximately 11% of cases are linked to variants in Voltage-Gated Sodium Channels (VGSCs). This study aims to broaden the genetic knowledge on painful SFN by applying Whole-Exome Sequencing (WES) in Early-Onset (EO) cases. A total of 88 patients from Italy (n = 52) and the Netherlands (n = 36), with a disease onset at age ≤ 45 years old and a Pain Numerical Rating Score ≥ 4, were recruited. After variant filtering and classification, WES analysis identified 142 potentially causative variants in 93 genes; 8 are Pathogenic, 15 are Likely Pathogenic, and 119 are Variants of Uncertain Significance. Notably, an enrichment of variants in transient receptor potential genes was observed, suggesting their role in pain modulation alongside VGSCs. A pathway analysis performed by comparing EO cases with 40 Italian healthy controls found enriched mutated genes in the "Nicotinic acetylcholine receptor signaling pathway". Targeting this pathway with non-opioid drugs could offer novel therapeutic avenues for painful SFN. Additionally, with this study we demonstrated that employing a gene panel of reported mutated genes could serve as an initial screening tool for SFN in genetic studies, enhancing clinical diagnostics.

Exome sequencing in multiple sclerosis families identifies 12 candidate genes and nominates biological pathways for the genesis of disease.

Multiple sclerosis (MS) is an inflammatory disease of the central nervous system characterized by myelin loss and neuronal dysfunction. Although the majority of patients do not present familial aggregation, Mendelian forms have been described. We performed whole-exome sequencing analysis in 132 patients from 34 multi-incident families, which nominated likely pathogenic variants for MS in 12 genes of the innate immune system that regulate the transcription and activation of inflammatory mediators. Rare missense or nonsense variants were identified in genes of the fibrinolysis and complement pathways (PLAU, MASP1, C2), inflammasome assembly (NLRP12), Wnt signaling (UBR2, CTNNA3, NFATC2, RNF213), nuclear receptor complexes (NCOA3), and cation channels and exchangers (KCNG4, SLC24A6, SLC8B1). These genes suggest a disruption of interconnected immunological and pro-inflammatory pathways as the initial event in the pathophysiology of familial MS, and provide the molecular and biological rationale for the chronic inflammation, demyelination and neurodegeneration observed in MS patients.

Inverse correlation of genetic risk score with age at onset in bout-onset and progressive-onset multiple sclerosis.

We correlated the weighted genetic risk score measured using 107 established susceptibility variants for multiple sclerosis (MS) with the age at onset in bout-onset (BOMS, n=906) and progressive-onset MS Italian patients (PrMS) (n=544). We observed an opposite relationship in the two disease courses: a higher weighted genetic risk score was associated with an earlier age at onset in BOMS (rho= -0.1; p=5 × 10(-3)) and a later age at onset in PrMS cases (rho=0.07; p=0.15) (p of difference of regression=1.4 × 10(-2)). These findings suggest that established MS risk variants anticipate the onset of the inflammatory phase, while they have no impact on, or even delay, the onset of the progressive phase.

Genetic Contribution to Medium-Term Disease Activity in Multiple Sclerosis.

Multiple sclerosis (MS) is a complex disorder characterized by high heterogeneity in terms of phenotypic expression, prognosis and treatment response. In the present study, we aimed to explore the genetic contribution to MS disease activity at different levels: genes, pathways and tissue-specific networks. Two cohorts of relapsing-remitting MS patients who started a first-line treatment (n = 1294) were enrolled to evaluate the genetic association with disease activity after 4 years of follow-up. The analyses were performed at whole-genome SNP and gene level, followed by the construction of gene-gene interaction networks specific for brain and lymphocytes. The resulting gene modules were evaluated to highlight key players from a topological and functional perspective. We identified 23 variants and 223 genes with suggestive association to 4-years disease activity, highlighting genes like PON2 involved in oxidative stress and in mitochondria functions and other genes, like ILRUN, involved in the modulation of the immune system. Network analyses led to the identification of a brain module composed of 228 genes and a lymphocytes module composed of 287 genes. The network analysis allowed us to prioritize genes relevant for their topological properties; among them, there are MPHOSPH9 (connector hub in both brain and lymphocyte module) and OPA1 (in brain module), two genes already implicated in MS. Modules showed the enrichment of both shared and tissue-specific pathways, mainly implicated in inflammation. In conclusion, our results suggest that the processes underlying disease activity act on shared mechanisms across brain and lymphocyte tissues.

Risk HLA Variants Affect the T-Cell Repertoire in Multiple Sclerosis.

BACKGROUND AND OBJECTIVES: The major histocompatibility complex (MHC) locus has a predominant role in the genetic predisposition to multiple sclerosis (MS), with 32 associations found to be involved. We aimed to investigate the impact of MHC MS-risk alleles on T-cell repertoire in patients with MS. METHODS: We studied 161 untreated patients with relapsing-remitting MS for whom Class I and II human leukocyte antigen (HLA) alleles were inferred from whole-genome genotyping data, and T-cell receptor (TCR) CDR3 sequences were obtained through next-generation sequencing. T-cell repertoire features including diversity, public clones, and architecture were evaluated. RESULTS: We identified 5 MS-risk loci associated with TCR diversity: HLA-DRB1*15:01 (7.65 × 10-3), rs9271366 (1.96 × 10-3), rs766848979 A (1.89 × 10-2), rs9277626 (2.95 × 10-2), and rs11751659 (1.92 × 10-2), with evidence of expanded clonotypes in carriers of risk alleles. Moreover, HLA-DRB1*15:01 (4.99 × 10-3), rs9271366 (6.54 × 10-3), rs1049079 C (4.37 × 10-2), AA DQΒ1 position -5 L (1.05 × 10-3), and AA DQΒ1 position 221 Q (9.39 × 10-4) showed an association with the CDR3 aminoacidic sequence architecture, suggesting an impact on the antigen recognition breadth as well. Evaluating the sharing of clones across MS-risk allele carrier individuals revealed the presence of highly shared clonotypes predicted to target viral antigens, including Epstein-Barr virus. DISCUSSION: Our study supports the association between MHC-risk alleles and macrofeatures of the T-cell repertoire in the context of MS. Further studies are needed to understand the underlying molecular mechanisms.

Combining Clinical and Genetic Data to Predict Response to Fingolimod Treatment in Relapsing Remitting Multiple Sclerosis Patients: A Precision Medicine Approach.

A personalized approach is strongly advocated for treatment selection in Multiple Sclerosis patients due to the high number of available drugs. Machine learning methods proved to be valuable tools in the context of precision medicine. In the present work, we applied machine learning methods to identify a combined clinical and genetic signature of response to fingolimod that could support the prediction of drug response. Two cohorts of fingolimod-treated patients from Italy and France were enrolled and divided into training, validation, and test set. Random forest training and robust feature selection were performed in the first two sets respectively, and the independent test set was used to evaluate model performance. A genetic-only model and a combined clinical-genetic model were obtained. Overall, 381 patients were classified according to the NEDA-3 criterion at 2 years; we identified a genetic model, including 123 SNPs, that was able to predict fingolimod response with an AUROC= 0.65 in the independent test set. When combining clinical data, the model accuracy increased to an AUROC= 0.71. Integrating clinical and genetic data by means of machine learning methods can help in the prediction of response to fingolimod, even though further studies are required to definitely extend this approach to clinical applications.

BDNF Val66Met Polymorphism Is Associated With Motor Recovery After Rehabilitation in Progressive Multiple Sclerosis Patients.

Background: Rehabilitation is fundamental for progressive multiple sclerosis (MS), but predictive biomarkers of motor recovery are lacking, making patient selection difficult. Motor recovery depends on synaptic plasticity, in which the Brain-Derived Neurotrophic Factor (BDNF) is a key player, through its binding to the Neurotrophic-Tyrosine Kinase-2 (NTRK2) receptor. Therefore, genetic polymorphisms in the BDNF pathway may impact motor recovery. The most well-known polymorphism in BDNF gene (rs6265) causes valine to methionine substitution (Val66Met) and it influences memory and motor learning in healthy individuals and neurodegenerative diseases. To date, no studies have explored whether polymorphisms in BDNF or NTRK2 genes may impact motor recovery in MS. Objectives: To assess whether genetic variants in BDNF and NTRK2 genes affect motor recovery after rehabilitation in progressive MS. Methods: The association between motor recovery after intensive neurorehabilitation and polymorphisms in BDNF (rs6265) and NTKR2 receptor (rs2289656 and rs1212171) was assessed using Six-Minutes-Walking-Test (6MWT), 10-Metres-Test (10MT) and Nine-Hole-Peg-Test (9HPT) in 100 progressive MS patients. Results: We observed greater improvement at 6MWT after rehabilitation in carriers of the BDNF Val66Met substitution, compared to BDNF Val homozygotes (p = 0.024). No significant association was found for 10MT and 9HPT. NTRK2 polymorphisms did not affect the results of motor function tests. Conclusion: BDNF Val66Met was associated with walking function improvement after rehabilitation in progressive MS patients. This result is in line with previous evidence showing a protective effect of Val66Met substitution on brain atrophy in MS. Larger studies are needed to explore its potential as a predictive biomarker of rehabilitation outcome.

Transcriptional effects of fingolimod treatment on peripheral T cells in relapsing remitting multiple sclerosis patients.

Aim: To investigate the transcriptional changes induced by Fingolimod (FTY) in T cells of relapsing remitting multiple sclerosis patients. Patients & methods: Transcriptomic changes after 6 months of FTY therapy were evaluated on T cells from 24 relapsing remitting multiple sclerosis patients through RNA-sequencing, followed by technical validation and pathway analysis. Results: Among differentially expressed genes, CX3CR1 and CCR7 resulted strongly up- and downregulated, respectively. Two relevant genes were validated with quantitative PCR and we largely confirmed findings from two previous microarray-based studies with similar design. Pathway analysis pointed to an involvement of processes related to immune function and cell migration. Conclusion: Our data support the evidence that FTY induces major transcriptional changes in genes involved in immune response and cell trafficking in T lymphocytes.

Exploring the Association of HLA Genetic Risk Burden on Thalamic and Hippocampal Atrophy in Multiple Sclerosis Patients.

Multiple sclerosis (MS) is a complex disease of the central nervous system for which human leukocyte antigen (HLA) alleles are major contributors to susceptibility. Several investigations have focused on the relationship between HLA and clinical parameters, while few studies have evaluated its correlation with brain magnetic resonance imaging (MRI) measures. We investigated the association between the HLA genetic burden (HLAGB), originating from the most updated HLA alleles associated with MS, and neuroimaging endophenotypes, with a specific focus on brain atrophy metrics. A monocentric Italian cohort of 334 MS patients with imputed HLA alleles and cross-sectional volumetric measures of white matter (WM), gray matter (GM), hippocampus, thalamus and T2-hyperintense lesions was investigated. Linear regression models with covariate adjustment were fitted for each metric. We detected no effect of HLAGB on WM and GM volumes. Interestingly, we found a marginal correlation between higher HLAGB and lower hippocampal volume (ß = -0.142, p = 0.063) and a nominal association between higher HLAGB and lower thalamic volume (ß = -0.299, p = 0.047). No association was found with T2 lesion volumes. The putative impact of higher HLAGB on hippocampus and thalamus suggests, if replicated in independent cohorts, a possible cumulative contribution of HLA risk loci on brain volumetric traits linked to clinical deficits in MS.

Involvement of NINJ2 Protein in Inflammation and Blood-Brain Barrier Transmigration of Monocytes in Multiple Sclerosis.

Multiple sclerosis (MS) is an inflammatory neurodegenerative disorder of the central nervous system (CNS). The migration of immune cells into the CNS is essential for its development, and plasma membrane molecules play an important role in triggering and maintaining the inflammation. We previously identified ninjurin2, a plasma membrane protein encoded by NINJ2 gene, as involved in the occurrence of relapse under Interferon-ß treatment in MS patients. The aim of the present study was to investigate the involvement of NINJ2 in inflammatory conditions and in the migration of monocytes through the blood-brain barrier (BBB). We observed that NINJ2 is downregulated in monocytes and in THP-1 cells after stimulation with the pro-inflammatory cytokine LPS, while in hCMEC/D3 cells, which represent a surrogate of the BBB, LPS stimulation increases its expression. We set up a transmigration assay using an hCMEC/D3 transwell-based model, finding a higher transmigration rate of monocytes from MS subjects compared to healthy controls (HCs) in the case of an activated hCMEC/D3 monolayer. Moreover, a positive correlation between NINJ2 expression in monocytes and monocyte migration rate was observed. Overall, our results suggest that ninjurin2 could be involved in the transmigration of immune cells into the CNS in pro-inflammatory conditions. Further experiments are needed to elucidate the exact molecular mechanisms.

A multi-step genomic approach prioritized TBKBP1 gene as relevant for multiple sclerosis susceptibility.

BACKGROUND: Over 200 genetic loci have been associated with multiple sclerosis (MS) explaining ~ 50% of its heritability, suggesting that additional mechanisms may account for the "missing heritability" phenomenon. OBJECTIVE: To analyze a large cohort of Italian individuals to identify markers associated with MS with potential functional impact in the disease. METHODS: We studied 2571 MS and 3234 healthy controls (HC) of continental Italian origin. Discovery phase included a genome wide association study (1727 MS, 2258 HC), with SNPs selected according to their association in the Italian cohort only or in a meta-analysis of signals with a cohort of European ancestry (4088 MS, 7144 HC). Top associated loci were then tested in two Italian cohorts through array-based genotyping (903 MS, 884 HC) and pool-based target sequencing (588 MS, 408 HC). Finally, functional prioritization through conditional eQTL and mQTL has been performed. RESULTS: Top associated signals overlap with already known MS loci on chromosomes 3 and 17. Three SNPs (rs4267364, rs8070463, rs67919208), all involved in the regulation of TBKBP1, were prioritized to be functionally relevant. CONCLUSIONS: No evidence of novel signal of association with MS specific for the Italian continental population has been found; nevertheless, two MS loci seems to play a relevant role, raising the interest to further investigations for TBKBP1 gene.

A Whole-Genome Sequencing Study Implicates GRAMD1B in Multiple Sclerosis Susceptibility.

While the role of common genetic variants in multiple sclerosis (MS) has been elucidated in large genome-wide association studies, the contribution of rare variants to the disease remains unclear. Herein, a whole-genome sequencing study in four affected and four healthy relatives of a consanguineous Italian family identified a novel missense c.1801T > C (p.S601P) variant in the GRAMD1B gene that is shared within MS cases and resides under a linkage peak (LOD: 2.194). Sequencing GRAMD1B in 91 familial MS cases revealed two additional rare missense and two splice-site variants, two of which (rs755488531 and rs769527838) were not found in 1000 Italian healthy controls. Functional studies demonstrated that GRAMD1B, a gene with unknown function in the central nervous system (CNS), is expressed by several cell types, including astrocytes, microglia and neurons as well as by peripheral monocytes and macrophages. Notably, GRAMD1B was downregulated in vessel-associated astrocytes of active MS lesions in autopsied brains and by inflammatory stimuli in peripheral monocytes, suggesting a possible role in the modulation of inflammatory response and disease pathophysiology.

An Investigation of the Role of Common and Rare Variants in a Large Italian Multiplex Family of Multiple Sclerosis Patients.

Known multiple sclerosis (MS) susceptibility variants can only explain half of the disease's estimated heritability, whereas low-frequency and rare variants may partly account for the missing heritability. Thus, here we sought to determine the occurrence of rare functional variants in a large Italian MS multiplex family with five affected members. For this purpose, we combined linkage analysis and next-generation sequencing (NGS)-based whole exome and whole genome sequencing (WES and WGS, respectively). The genetic burden attributable to known common MS variants was also assessed by weighted genetic risk score (wGRS). We found a significantly higher burden of common variants in the affected family members compared to that observed among sporadic MS patients and healthy controls (HCs). We also identified 34 genes containing at least one low-frequency functional variant shared among all affected family members, showing a significant enrichment in genes involved in specific biological processes-particularly mRNA transport-or neurodegenerative diseases. Altogether, our findings point to a possible pathogenic role of different low-frequency functional MS variants belonging to shared pathways. We propose that these rare variants, together with other known common MS variants, may account for the high number of affected family members within this MS multiplex family.

Contribution of Rare and Low-Frequency Variants to Multiple Sclerosis Susceptibility in the Italian Continental Population.

Genome-wide association studies identified over 200 risk loci for multiple sclerosis (MS) focusing on common variants, which account for about 50% of disease heritability. The goal of this study was to investigate whether low-frequency and rare functional variants, located in MS-established associated loci, may contribute to disease risk in a relatively homogeneous population, testing their cumulative effect (burden) with gene-wise tests. We sequenced 98 genes in 588 Italian patients with MS and 408 matched healthy controls (HCs). Variants were selected using different filtering criteria based on allelic frequency and in silico functional impacts. Genes showing a significant burden (n = 17) were sequenced in an independent cohort of 504 MS and 504 HC. The highest signal in both cohorts was observed for the disruptive variants (stop-gain, stop-loss, or splicing variants) located in EFCAB13, a gene coding for a protein of an unknown function (p < 10-4). Among these variants, the minor allele of a stop-gain variant showed a significantly higher frequency in MS versus HC in both sequenced cohorts (p = 0.0093 and p = 0.025), confirmed by a meta-analysis on a third independent cohort of 1298 MS and 1430 HC (p = 0.001) assayed with an SNP array. Real-time PCR on 14 heterozygous individuals for this variant did not evidence the presence of the stop-gain allele, suggesting a transcript degradation by non-sense mediated decay, supported by the evidence that the carriers of the stop-gain variant had a lower expression of this gene (p = 0.0184). In conclusion, we identified a novel low-frequency functional variant associated with MS susceptibility, suggesting the possible role of rare/low-frequency variants in MS as reported for other complex diseases.

Impact of multiple sclerosis risk loci in postinfectious neurological syndromes.

BACKGROUND: The genetic component of multiple sclerosis (MS) is now set to 200 autosomal common variants. However, it is unclear how genetic knowledge be clinically used in the differential diagnosis between MS and other inflammatory conditions like adult-onset postinfectious neurological syndromes (PINS). The aim of this study was to investigate whether PINS and MS have a shared genetic background using an updated polygenic risk scores. METHODS: Eighty-eight PINS patients have been consecutively recruited between 1996 and 2016 at Mondino Foundation of Pavia, diagnosed according to clinical, MRI and CSF findings and followed-up for several years. Patients were typed using Illumina array, and genotypes imputed using the 1000 Genomes Project reference panel. A weighted genetic risk score (wGRS) has been calculated based on autosomal MS risk loci derived from large-scale studies, and an HLA genetic burden (HLAGB) was also calculated on loci associated to MS. RESULTS: PINS occurred as an episode of myelitis in 44% of patients, encephalomyelitis in 44%, and encephalitis in remaining cases, with an involvement of peripheral nervous system in 41% of patients. Mean age of onset was 50.1 years, and female:male ratio was 1.4. Patients were followed-up for a mean of 7.2 years, and at last visit 55% had a low disability grade (mRS 0-1). Disease was monophasic in 67% of patients, relapsing in 18% and chronic-progressive in 15%. The wGRS of PINS cases was comparable to 370 healthy controls, while significantly lower compared to 907 bout-onset MS (BOMS) cases (wGRS= 20.9 vs 21.2; p<0.0001). The difference was even larger for PINS with peripheral nervous system involvement (wGRS=20.6) vs BOMS. CONCLUSION: The distinction between MS and PINS is not easy to make in clinical practice. However, our study shows that the new set of MS risk alleles does not confer increased susceptibility to PINS. These data support the importance to discriminate these cases from MS with pathophysiological and therapeutic implications.

CHRNA7 Gene and Response to Cholinesterase Inhibitors in an Italian Cohort of Alzheimer's Disease Patients.

Previous studies suggest that genetic variants in CHRNA7, which encodes for the major subunit of the acetylcholine receptor (α7-nAChR), are associated with the clinical response to cholinesterase inhibitors (ChEI) in Alzheimer's disease (AD) patients. We sought to replicate the association of two SNPs in the CHRNA7 gene, rs6494223 and rs8024987, with response to ChEI treatment in an Italian cohort of 169 AD patients, further extending the study to gene-level analysis. None of the tested variants was associated with clinical response. However, rs6494223 showed a consistent effect direction (OR = 1.4; p = 0.17), which after meta-analysis with previous study yielded a significant result (OR = 1.57, p = 0.02, I2 = 0%).

Impact of MS genetic loci on familial aggregation, clinical phenotype, and disease prediction.

OBJECTIVE: To investigate the role of known multiple sclerosis (MS)-associated genetic variants in MS familial aggregation, clinical expression, and accuracy of disease prediction in sporadic and familial cases. METHODS: A total of 1,443 consecutive patients were screened for MS and familial autoimmune history in a hospital-based Italian cohort. Among them, 461 sporadic and 93 familial probands were genotyped for 107 MS-associated polymorphisms. Their effect sizes were combined to calculate the weighted genetic risk score (wGRS). RESULTS: Family history of MS was reported by 17.2% of probands, and 33.8% reported a familial autoimmune disorder, with autoimmune thyroiditis and psoriasis being the most frequent. No difference in wGRS was observed between sporadic and familial MS cases. In contrast, a lower wGRS was observed in probands with greater familial aggregation (>1 first-degree relative or >2 relatives with MS) (p = 0.03). Also, female probands of familial cases with greater familial aggregation had a lower wGRS than sporadic cases (p = 0.0009) and male probands of familial cases (p = 0.04). An inverse correlation between wGRS and age at onset was observed (p = 0.05). The predictive performance of the genetic model including all known MS variants was modest but greater in sporadic vs familial cases (area under the curve = 0.63 and 0.57). CONCLUSIONS: Additional variants outside the known MS-associated loci, rare variants, and/or environmental factors may explain disease occurrence within families; in females, hormonal and epigenetic factors probably have a predominant role in explaining familial aggregation. The inclusion of these additional factors in future versions of aggregated genetic measures could improve their predictive ability.