Burden of rare coding variants in an Italian cohort of familial multiple sclerosis.

BACKGROUND: Multiple Sclerosis (MS) is a chronic inflammatory and neurodegenerative demyelinating disease of the central nervous system. It is a complex and heterogeneous disease caused by a combination of genetic and environmental factors, and it can cluster in families. OBJECTIVE: to evaluate at gene-level the aggregate contribution of predicted damaging low-frequency and rare variants to MS risk in multiplex families. METHODS: We performed whole exome sequencing (WES) in 28 multiplex MS families with at least 3 MS cases (81 affected and 42 unaffected relatives) and 38 unrelated healthy controls. A gene-based burden test was then performed, focusing on two sets of candidate genes: i) literature-driven selection and ii) data-driven selection. RESULTS: We identified 11 genes enriched with predicted damaging low-frequency and rare variants in MS compared to healthy individuals. Among them, UBR2 and DST were the two genes with the strongest enrichment (p = 5 × 10-4 and 3 × 10-4, respectively); interestingly enough the association signal in UBR2 is driven by rs62414610, which was present in 25% of analysed families. CONCLUSION: Despite limitations, this is one of the first studies evaluating the aggregate contribution of predicted damaging low-frequency and rare variants in MS families using WES data. A replication effort in independent cohorts is warranted to validate our findings and to evaluate the role of identified genes in MS pathogenesis.

Transcriptomic Analysis of Peripheral Monocytes upon Fingolimod Treatment in Relapsing Remitting Multiple Sclerosis Patients.

Fingolimod (FTY), a second-line oral drug approved for relapsing remitting Multiple Sclerosis (RRMS) acts in preventing lymphocyte migration outside lymph nodes; moreover, several lines of evidence suggest that it also inhibits myeloid cell activation. In this study, we investigated the transcriptional changes induced by FTY in monocytes in order to better elucidate its mechanism of action. CD14+ monocytes were collected from 24 RRMS patients sampled at baseline and after 6 months of treatment and RNA profiles were obtained through next-generation sequencing. We conducted pathway and sub-paths analysis, followed by centrality analysis of cell-specific interactomes on differentially expressed genes (DEGs). We investigated also the predictive role of baseline monocyte transcription profile in influencing the response to FTY therapy. We observed a marked down-regulation effect (60 down-regulated vs. 0 up-regulated genes). Most of the down-regulated DEGs resulted related with monocyte activation and migration like IL7R, CCR7 and the Wnt signaling mediators LEF1 and TCF7. The involvement of Wnt signaling was also confirmed by subpaths analyses. Furthermore, pathway and network analyses showed an involvement of processes related to immune function and cell migration. Baseline transcriptional profile of the HLA class II gene HLA-DQA1 and HLA-DPA1 were associated with evidence of disease activity after 2 years of treatment. Our data support the evidence that FTY induces major transcriptional changes in monocytes, mainly regarding genes involved in cell trafficking and immune cell activation. The baseline transcriptional levels of genes associated with antigen presenting function were associated with disease activity after 2 years of FTY treatment.

Identification of differential DNA methylation associated with multiple sclerosis: A family-based study.

Multiple Sclerosis (MS) is caused by a still unknown interplay between genetic and environmental factors. Epigenetics, including DNA methylation, represents a model for environmental factors to influence MS risk. Twenty-six affected and 26 unaffected relatives from 8 MS multiplex families were analysed in a multicentric Italian study using MeDIP-Seq, followed by technical validation and biological replication in two additional families of differentially methylated regions (DMRs) using SeqCap Epi Choice Enrichment kit (Roche®). Associations from MeDIP-Seq across families were combined with aggregation statistics, yielding 162 DMRs at FDR ≤ 0.1. Technical validation and biological replication led to 2 hypo-methylated regions, which point to NTM and BAI3 genes, and to 2 hyper-methylated regions in PIK3R1 and CAPN13. These 4 novel regions contain genes of potential interest that need to be tested in larger cohorts of patients.

Clinical response to Nabiximols correlates with the downregulation of immune pathways in multiple sclerosis.

BACKGROUND AND PURPOSE: Nabiximols (Sativex® ) is a cannabinoid-based compound used for the treatment of moderate to severe spasticity in multiple sclerosis (MS). The aim of the study was to investigate the effect of the administration of Nabiximols on blood transcriptome profile of patients with MS and to interpret it in the context of pathways and networks. METHODS: Whole-genome expression profiling was performed in whole blood of 33 subjects with MS at baseline and after 4 weeks of drug treatment. Patients were classified as responders (n = 19) and non-responders (n = 14). Pathway and network analyses on genes modulated by the drug were performed, followed by in vitro stimulation of peripheral blood mononuclear cells with pro-inflammatory agents to support the immunomodulatory properties of the drug. RESULTS: Individual effect size was modest; however, we observed a downregulation of several immune-related pathways after 4 weeks of treatment, which was more pronounced when restricting analyses to responders. Interesting hub molecules functionally related to the immune system emerged from network analysis, including NFKB1, FYN, MAP14 and TP53. The immunomodulatory properties of the drug were confirmed through in vitro assays in peripheral blood mononuclear cells collected from patients with MS. CONCLUSIONS: Our findings support the immunomodulatory activity of cannabinoids in patients with MS. Further studies in more specific cell types are needed to refine these results.

Response to interferon-beta treatment in multiple sclerosis patients: a genome-wide association study.

Up to 50% of multiple sclerosis (MS) patients do not respond to interferon-beta (IFN-ß) treatment and determination of response requires lengthy clinical follow-up of up to 2 years. Response predictive genetic markers would significantly improve disease management. We aimed to identify IFN-ß treatment response genetic marker(s) by performing a two-stage genome-wide association study (GWAS). The GWAS was carried out using data from 151 Australian MS patients from the ANZgene/WTCCC2 MS susceptibility GWAS (responder (R)=51, intermediate responders=24 and non-responders (NR)=76). Of the single-nucleotide polymorphisms (SNP) that were validated in an independent group of 479 IFN-ß-treated MS patients from Australia, Spain and Italy (R=273 and NR=206), eight showed evidence of association with treatment response. Among the replicated associations, the strongest was observed for FHIT (Fragile Histidine Triad; combined P-value 6.74 × 10-6) and followed by variants in GAPVD1 (GTPase activating protein and VPS9 domains 1; combined P-value 5.83 × 10-5) and near ZNF697 (combined P-value 8.15 × 10-5).

Pharmacogenetic study of long-term response to interferon-ß treatment in multiple sclerosis.

The aim of the study is the identification of genetic factors that influence the long-term response to interferon-ß (IFNß) (4-year follow-up). We performed a genome-wide association study in 337 IFNß-treated Italian multiple sclerosis patients at the extreme of treatment response, and we meta-analyzed association effects, integrating results with pathway analysis, gene-expression profiling of IFNß-stimulated peripheral blood mononuclear cells from 20 healthy controls (HC) and expression quantitative locus (eQTL) analyses. From meta-analysis, 43 markers were associated at P<10-4, and two of them (rs7298096 and rs4726460) pointed to two genes, NINJ2 and TBXAS1, that were significantly downregulated after IFNß stimulation in HC (P=3.1 × 10-9 and 5.6 × 10-10). We also observed an eQTL effect for the allele associated with favorable treatment response (rs4726460A); moreover, TBXAS1 appeared downregulated upon IFNß administration (ß=-0.39; P=0.02). Finally, we found an enrichment of pathways related to inflammatory processes and presynaptic membrane, the latter with involvement of genes related to glutamatergic system (GRM3 and GRIK2), confirming its potential role in the response to IFNß.

Analysis of genes, pathways and networks involved in disease severity and age at onset in primary-progressive multiple sclerosis.

BACKGROUND: The role of genetic factors in influencing the clinical expression of multiple sclerosis (MS) is unclear. OBJECTIVE: The objective of this paper is to identify genes, pathways and networks implicated in age at onset (AAO) and severity, measured using the Multiple Sclerosis Severity Score (MSSS), of primary-progressive MS (PPMS). METHODS: We conducted a genome-wide association study (GWAS) of 470 PPMS patients of Italian origin:. Allelic association of 296,589 SNPs with AAO and MSSS was calculated. Pathway and network analyses were also conducted using different tools. RESULTS: No single association signal exceeded genome-wide significance in AAO and MSSS analyses. Nominally associated genes to AAO and MSSS were enriched in both traits for 10 pathways, including: "oxidative phosphorylation" (FDRAAO=9*10(-4); FDRMSSS=3.0*10(-2)), "citrate (TCA) cycle" (FDRAAO=1.6*10(-2); FDRMSSS=3.2*10(-3)), and "B cell receptor signaling" (FDRAAO=3.1*10(-2); FDRMSSS=2.2*10(-3)). In addition, an enrichment of "chemokine signaling pathway" (FDR=9*10(-4)) for AAO and of "leukocyte transendothelial migration" (FDR=2.4*10(-3)) for MSSS trait was observed, among others. Network analysis revealed that p53 and CREB1 were central hubs for AAO and MSSS traits, respectively. CONCLUSIONS: Despite the fact that no major effect signals emerged in the present GWAS, our data suggest that genetic variants acting in the context of oxidative stress and immune dysfunction could modulate the onset and severity of PPMS.