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.

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.

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.

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.

Clinical and pathological findings in neurolymphomatosis: Preliminary association with gene expression profiles in sural nerves.

Although inflammation appears to play a role in neurolymphomatosis (NL), the mechanisms leading to degeneration in the peripheral nervous system are poorly understood. The purpose of this exploratory study was to identify molecular pathways underlying NL pathogenesis, combining clinical and neuropathological investigation with gene expression (GE) studies. We characterized the clinical and pathological features of eight patients with NL. We further analysed GE changes in sural nerve biopsies obtained from a subgroup of NL patients (n=3) and thirteen patients with inflammatory neuropathies as neuropathic controls. Based on the neuropathic symptoms and signs, NL patients were classified into three forms of neuropathy: chronic symmetrical sensorimotor polyneuropathy (SMPN, n=3), multiple mononeuropathy (MN, n=4) and acute motor-sensory axonal neuropathy (AMSAN, n=1). Predominantly diffuse malignant cells infiltration of epineurium was present in chronic SMPN, whereas endoneurial perivascular cells invasion was observed in MN. In contrast, diffuse endoneurium malignant cells localization occurred in AMSAN. We identified alterations in the expression of 1266 genes, with 115 up-regulated and 1151 down-regulated genes, which were mainly associated with ribosomal proteins (RP) and olfactory receptors (OR) signaling pathways, respectively. Among the top up-regulated genes were actin alpha 1 skeletal muscle (ACTA1) and desmin (DES). Similarly, in NL nerves ACTA1, DES and several RPs were highly expressed, associated with endothelial cells and pericytes abnormalities. Peripheral nerve involvement may be due to conversion towards a more aggressive phenotype, potentially explaining the poor prognosis. The candidate genes reported in this study may be a source of clinical biomarkers for NL.

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.

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.

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.

Genomic and functional evaluation of TNFSF14 in multiple sclerosis susceptibility.

Among multiple sclerosis (MS) susceptibility genes, the strongest non-human leukocyte antigen (HLA) signal in the Italian population maps to the TNFSF14 gene encoding LIGHT, a glycoprotein involved in dendritic cell (DC) maturation. Through fine-mapping in a large Italian dataset (4,198 patients with MS and 3,903 controls), we show that the TNFSF14 intronic SNP rs1077667 is the primarily MS-associated variant in the region. Expression quantitative trait locus (eQTL) analysis indicates that the MS risk allele is significantly associated with reduced TNFSF14 messenger RNA levels in blood cells, which is consistent with the allelic imbalance in RNA-Seq reads (P < 0.0001). The MS risk allele is associated with reduced levels of TNFSF14 gene expression (P < 0.01) in blood cells from 84 Italian patients with MS and 80 healthy controls (HCs). Interestingly, patients with MS are lower expressors of TNFSF14 compared to HC (P < 0.007). Individuals homozygous for the MS risk allele display an increased percentage of LIGHT-positive peripheral blood myeloid DCs (CD11c+, P = 0.035) in 37 HCs, as well as in in vitro monocyte-derived DCs from 22 HCs (P = 0.04). Our findings suggest that the intronic variant rs1077667 alters the expression of TNFSF14 in immune cells, which may play a role in MS pathogenesis.

Assessment of the genetic contribution to brain magnetic resonance imaging lesion load and atrophy measures in multiple sclerosis patients.

BACKGROUND AND PURPOSE: Multiple sclerosis (MS) susceptibility is influenced by genetics; however, little is known about genetic determinants of disease expression. We aimed at assessing genetic factors influencing quantitative neuroimaging measures in two cohorts of progressive MS (PMS) and relapsing-remitting MS (RRMS) patients. METHODS: Ninety-nine PMS and 214 RRMS patients underwent a 3-T brain magnetic resonance imaging (MRI) scan, with the measurement of five MRI metrics including T2 lesion volumes and measures of white matter, grey matter, deep grey matter, and hippocampal volumes. A candidate pathway strategy was adopted; gene set analysis was carried out to estimate cumulative contribution of genes to MRI phenotypes, adjusting for relevant confounders, followed by single nucleotide polymorphism (SNP) regression analysis. RESULTS: Seventeen Kyoto Encyclopedia of Genes and Genomes pathways and 42 Gene Ontology (GO) terms were tested. We additionally included in the analysis genes with enriched expression in brain cells. Gene set analysis revealed a differential pattern of association across the two cohorts, with processes related to sodium homeostasis being associated with grey matter volume in PMS (p = 0.002), whereas inflammatory-related GO terms such as adaptive immune response and regulation of inflammatory response appeared to be associated with T2 lesion volume in RRMS (p = 0.004 and p = 0.008, respectively). As for SNPs, the rs7104613T mapping to SPON1 gene was associated with reduced deep grey matter volume (ß = -0.731, p = 3.2*10-7 ) in PMS, whereas we found evidence of association between white matter volume and rs740948A mapping to SEMA3A gene (ß = 22.04, p = 5.5*10-6 ) in RRMS. CONCLUSIONS: Our data suggest a different pattern of associations between MRI metrics and functional processes across MS disease courses, suggesting different phenomena implicated in MS.

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.

Author Correction: Next Generation Sequencing of Pooled Samples: Guideline for Variants' Filtering.

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

A pharmacogenetic study implicates NINJ2 in the response to Interferon-ß in multiple sclerosis.

BACKGROUND: Multiple sclerosis (MS) is a disease in which biomarker identification is fundamental to predict response to treatments and to deliver the optimal drug to patients. We previously found an association between rs7298096, a polymorphism upstream to the NINJ2 gene, and the 4-year response to interferon-ß (IFNß) treatment in MS patients. OBJECTIVES: To analyse the association between rs7298096 and time to first relapse (TTFR) during IFNß therapy in MS patients and to better investigate its functional role. METHODS: Survival analysis was applied in three MS cohorts from different countries (n = 1004). We also studied the role of the polymorphism on gene expression using GTEx portal and a luciferase assay. We interrogated GEO datasets to explore the relationship between NINJ2 expression, IFNß and TTFR. RESULTS: Rs7298096AA patients show a shorter TTFR than rs7298096G-carriers (Pmeta-analysis = 3 × 10-4, hazard ratio = 1.41). Moreover, rs7298096AA is associated with a higher NINJ2 expression in blood (p = 7.0 × 10-6), which was confirmed in vitro (p = 0.009). Finally, NINJ2 expression is downregulated by IFNß treatment and related to TTFR. CONCLUSIONS: Rs7298096 could influence MS disease activity during IFNß treatment by modulating NINJ2 expression in blood. The gene encodes for an adhesion molecule involved in inflammation and endothelial cells activation, supporting its role in MS.

No C9orf72 repeat expansion in patients with primary progressive multiple sclerosis.

Pathological repeat expansion (RE) of the C9orf72 hexanucleotide sequence is associated to amyotrophic lateral sclerosis (ALS) and frontotemporal dementia disease continuum, although other heterogeneous clinical phenotypes have been documented. The occurrence of multiple sclerosis (MS) in some C9orf72 carriers with a more severe ALS disease course has suggested a possible modifying role for MS. However, C9orf72 RE seems not to play a role in MS pathogenesis. In this study, we screened C9orf72 in 189 Italian patients with primary progressive MS (PPMS), a rare clinical form characterized by less inflammation over neurodegenerative features. We failed to detect C9orf72 RE, but a significant representation of intermediate alleles (≥ 20 units) was observed in our PPMS cohort (2.1%) compared to healthy controls (0%, p < 0.05). In the normal range, allele distribution showed a trimodal pattern (2,5,8-repeat units) in PPMS and healthy controls with no significant difference. Our findings further demonstrate that C9orf72 RE is not genetically associated to MS spectrum, but suggest that intermediate alleles may represent risk factors as already reported for Parkinson disease.

Laser capture microdissection for transcriptomic profiles in human skin biopsies.

BACKGROUND: The acquisition of reliable tissue-specific RNA sequencing data from human skin biopsy represents a major advance in research. However, the complexity of the process of isolation of specific layers from fresh-frozen human specimen by laser capture microdissection, the abundant presence of skin nucleases and RNA instability remain relevant methodological challenges. We developed and optimized a protocol to extract RNA from layers of human skin biopsies and to provide satisfactory quality and amount of mRNA sequencing data. RESULTS: The protocol includes steps of collection, embedding, freezing, histological coloration and relative optimization to preserve RNA extracted from specific components of fresh-frozen human skin biopsy of 14 subjects. Optimization of the protocol includes a preservation step in RNALater® Solution, the control of specimen temperature, the use of RNase Inhibitors and the time reduction of the staining procedure. The quality of extracted RNA was measured using the percentage of fragments longer than 200 nucleotides (DV200), a more suitable measurement for successful library preparation than the RNA Integrity Number (RIN). RNA was then enriched using the TruSeq® RNA Access Library Prep Kit (Illumina®) and sequenced on HiSeq® 2500 platform (Illumina®). Quality control on RNA sequencing data was adequate to get reliable data for downstream analysis. CONCLUSIONS: The described implemented and optimized protocol can be used for generating transcriptomics data on skin tissues, and it is potentially applicable to other tissues. It can be extended to multicenter studies, due to the introduction of an initial step of preservation of the specimen that allowed the shipment of biological samples.