Not all roads lead to the immune system: the genetic basis of multiple sclerosis severity.

Multiple sclerosis is a leading cause of neurological disability in adults. Heterogeneity in multiple sclerosis clinical presentation has posed a major challenge for identifying genetic variants associated with disease outcomes. To overcome this challenge, we used prospectively ascertained clinical outcomes data from the largest international multiple sclerosis registry, MSBase. We assembled a cohort of deeply phenotyped individuals of European ancestry with relapse-onset multiple sclerosis. We used unbiased genome-wide association study and machine learning approaches to assess the genetic contribution to longitudinally defined multiple sclerosis severity phenotypes in 1813 individuals. Our primary analyses did not identify any genetic variants of moderate to large effect sizes that met genome-wide significance thresholds. The strongest signal was associated with rs7289446 (ß = -0.4882, P = 2.73 × 10-7), intronic to SEZ6L on chromosome 22. However, we demonstrate that clinical outcomes in relapse-onset multiple sclerosis are associated with multiple genetic loci of small effect sizes. Using a machine learning approach incorporating over 62 000 variants together with clinical and demographic variables available at multiple sclerosis disease onset, we could predict severity with an area under the receiver operator curve of 0.84 (95% CI 0.79-0.88). Our machine learning algorithm achieved positive predictive value for outcome assignation of 80% and negative predictive value of 88%. This outperformed our machine learning algorithm that contained clinical and demographic variables alone (area under the receiver operator curve 0.54, 95% CI 0.48-0.60). Secondary, sex-stratified analyses identified two genetic loci that met genome-wide significance thresholds. One in females (rs10967273; ßfemale = 0.8289, P = 3.52 × 10-8), the other in males (rs698805; ßmale = -1.5395, P = 4.35 × 10-8), providing some evidence for sex dimorphism in multiple sclerosis severity. Tissue enrichment and pathway analyses identified an overrepresentation of genes expressed in CNS compartments generally, and specifically in the cerebellum (P = 0.023). These involved mitochondrial function, synaptic plasticity, oligodendroglial biology, cellular senescence, calcium and G-protein receptor signalling pathways. We further identified six variants with strong evidence for regulating clinical outcomes, the strongest signal again intronic to SEZ6L (adjusted hazard ratio 0.72, P = 4.85 × 10-4). Here we report a milestone in our progress towards understanding the clinical heterogeneity of multiple sclerosis outcomes, implicating functionally distinct mechanisms to multiple sclerosis risk. Importantly, we demonstrate that machine learning using common single nucleotide variant clusters, together with clinical variables readily available at diagnosis can improve prognostic capabilities at diagnosis, and with further validation has the potential to translate to meaningful clinical practice change.

Predictive factors and early biomarkers of response in multiple sclerosis patients treated with natalizumab.

There are an increasing number of treatments available for multiple sclerosis (MS). The early identification of optimal responders to individual treatments is important to achieve individualized therapy. With this aim, we performed a multicenter retrospective longitudinal study including 186 MS patients treated with natalizumab who were followed for 2 years. We analyzed the following variables at recruitment: sex, current age, age at disease onset, disease duration, EDSS, number of T2 and Gd + lesions, IgG and IgM oligoclonal bands, HLA class II (DR, DRB, DQA, DQB, and DRB1*15:01), IgG and IgM antibody titers against human herpesvirus 6 (HHV-6) and the antibody response to Epstein-Barr virus (EBV) through the measurement of the anti-EBNA-1 and anti-VCA IgG titers, in relation to clinical response (no relapses or disability progression), and to NEDA-3 (no evidence of disease activity in terms of clinical response and no changes in MRI scans either) after 2-years follow-up. Baseline EDSS score, baseline EBNA-1 IgG titers and percentage change of HHV6 IgG titers between baseline and 6 month visits were significantly different in clinical responders and in NEDA-3 status (all of them remained significant in the multivariate analysis). We identified three variables for the early identification of natalizumab optimal responders in a rapid and cost-effective approach.

Splice-site variant in ACSL5: a marker promoting opposing effect on cell viability and protein expression.

Long-chain Acyl-CoA synthetases (ACSLs) activate fatty acids (FAs) by thioesterification with Coenzyme A (CoA), generating FA-CoAs. These products are essential for lipid metabolism and carcinogenesis. In previous study, we identified an intronic variant rs2256368:A>G, whose G allele promotes exon 20 skipping in up to 43% of ACSL5 transcripts but its functional relevance is unclear. Here, we compared the expression of splice (Spl) and nonsplice (NSpl) ACSL5 variants and the effect on cell viability under culture conditions that force cells to metabolize fatty acids. We found that lymphoblastoid cell lines from 1000 Genomes Project, bearing Spl genotypes, showed a reduced expression of total ACSL5 protein due to an inefficient translation of the Spl RNA. These cells impaired growth in cultures with phorbol myristate acetate-ionomycin (PMA-Io) or medium deprived of glucose, while production of reactive oxygen species increased in PMA-Io. Specific ACSL5-isoform transfection in HEK239T (kidney), U87 (astroglioma), and HOG (oligodendrocyte) cells showed the Spl protein to be the causal factor of cell-growth inhibition, despite its reduced protein expression. Our findings indicate that the variant rs2256368:A>G can predict a growth inhibitory activity, caused by the Spl isoform of ACSL5 protein, opposed to the activity of the NSpl. Deep understanding of its functioning might have application in metabolic diseases and cancer.

SP140 regulates the expression of immune-related genes associated with multiple sclerosis and other autoimmune diseases by NF-κB inhibition.

SP140 locus has been associated with multiple sclerosis (MS) as well as other autoimmune diseases by genome-wide association studies (GWAS). The causal variant of these associations (rs28445040-T) alters the splicing of the SP140 gene transcripts reducing the protein expression. We aimed to understand why the reduction of SP140 expression produced by the risk variant can increase the susceptibility to MS. To this end, we determined by RNA sequencing (RNA-seq) analysis the differentially expressed genes after SP140 silencing in lymphoblastoid cell lines (LCLs). We analyzed these genes by gene ontology (GO), comparative transcriptome profiles, enrichment of transcription factors (TFs) in the promoters of these genes and colocalization with GWAS risk variants. We also monitored the activity of the nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB) in SP140-silenced cells by luciferase reporter system. We identified 100 genes that were up-regulated and 22 genes down-regulated in SP140-silenced LCLs. GO analysis revealed that genes affected by SP140 were involved in regulation of cytokine production, inflammatory response and cell-cell adhesion. We observed enrichment of NF-κB TF in the promoter of up-regulated genes and NF-κB-increased activity in SP140-silenced cell lines. We showed enrichment of genes regulated by SP140 in GWAS-detected risk loci for MS (14.63 folds), Crohn's disease (4.82 folds) and inflammatory bowel disease (4.47 folds), not observed in other unrelated immune diseases. Our findings showed that SP140 is an important repressor of genes implicated in inflammation, suggesting that decreased expression of SP140, promoted by the rs28445040-T risk variant, may lead to up-regulation of these genes by means of NF-κB inhibition in B cells.

A functional variant that affects exon-skipping and protein expression of SP140 as genetic mechanism predisposing to multiple sclerosis.

Several variants in strong linkage disequilibrium (LD) at the SP140 locus have been associated with multiple sclerosis (MS), Crohn's disease (CD) and chronic lymphocytic leukemia (CLL). To determine the causal polymorphism, we have integrated high-density data sets of expression quantitative trait loci (eQTL), using GEUVADIS RNA sequences and 1000 Genomes genotypes, with MS-risk variants of the high-density Immunochip array performed by the International Multiple Sclerosis Genetic Consortium (IMSGC). The variants most associated with MS were also correlated with a decreased expression of the full-length RNA isoform of SP140 and an increase of an isoform lacking exon 7. By exon splicing assay, we have demonstrated that the rs28445040 variant was the causal factor for skipping of exon 7. Western blots of peripheral blood mononuclear cells from MS patients showed a significant allele-dependent reduction of the SP140 protein expression. To confirm the association of this functional variant with MS and to compare it with the best-associated variant previously reported by GWAS (rs10201872), a case-control study including 4384 MS patients and 3197 controls was performed. Both variants, in strong LD (r(2) = 0.93), were found similarly associated with MS [P-values, odds ratios: 1.9E-9, OR = 1.35 (1.22-1.49) and 4.9E-10, OR = 1.37 (1.24-1.51), respectively]. In conclusion, our data uncover the causal variant for the SP140 locus and the molecular mechanism associated with MS risk. In addition, this study and others previously reported strongly suggest that this functional variant may be shared with other immune-mediated diseases as CD and CLL.

TNFRSF1A polymorphisms rs1800693 and rs4149584 in patients with multiple sclerosis.

OBJECTIVES: To investigate the roles of 2 polymorphisms of the tumor necrosis factor (TNF) receptor superfamily member 1A (TNFRSF1A) gene, rs1800693 (a common variant) and rs4149584 (a coding polymorphism that results in an amino acid substitution-R92Q), as genetic modifiers of multiple sclerosis (MS), and to evaluate their potential functional implications in the disease. METHODS: The effects of rs1800693 and rs4149584 on 2 measures of disease severity, age at disease onset and Multiple Sclerosis Severity Score, were analyzed in 2,032 patients with MS. In a subgroup of patients, serum levels of the soluble form of TNF-R1 (sTNF-R1) were measured by ELISA; mRNA expression levels of the full-length TNF-R1 and Δ6-TNF-R1 isoform were investigated in peripheral blood mononuclear cells (PBMC) by real-time PCR; cell surface expression of the TNF-R1 was determined in T cells by flow cytometry. RESULTS: For rs4149584, R92Q carriers were younger at disease onset and progressed slower compared to noncarriers. However, no association with disease severity was observed for rs1800693. Serum levels of sTNF-R1 and mRNA expression levels of the full-length receptor were significantly increased in patients with MS carrying the R92Q mutation (p = 0.003 and p = 0.011, respectively), but similarly distributed among rs1800693 genotypes; cell surface TNF-R1 expression in T cells did not differ between rs4149584 and rs1800693 genotypes. The truncated soluble Δ6-TNF-R1 isoform was identified in PBMC from patients carrying the risk allele for rs1800693. CONCLUSIONS: These findings suggest that both rs1800693 and rs4149584 TNFRSF1A polymorphisms have functional consequences in the TNF-R1.

Genome-wide association study of multiple sclerosis confirms a novel locus at 5p13.1.

Multiple Sclerosis (MS) is the most common progressive and disabling neurological condition affecting young adults in the world today. From a genetic point of view, MS is a complex disorder resulting from the combination of genetic and non-genetic factors. We aimed to identify previously unidentified loci conducting a new GWAS of Multiple Sclerosis (MS) in a sample of 296 MS cases and 801 controls from the Spanish population. Meta-analysis of our data in combination with previous GWAS was done. A total of 17 GWAS-significant SNPs, corresponding to three different loci were identified:HLA, IL2RA, and 5p13.1. All three have been previously reported as GWAS-significant. We confirmed our observation in 5p13.1 for rs9292777 using two additional independent Spanish samples to make a total of 4912 MS cases and 7498 controls (ORpooled = 0.84; 95%CI: 0.80-0.89; p = 1.36 × 10-9). This SNP differs from the one reported within this locus in a recent GWAS. Although it is unclear whether both signals are tapping the same genetic association, it seems clear that this locus plays an important role in the pathogenesis of MS.

Sample reproducibility of genetic association using different multimarker TDTs in genome-wide association studies: characterization and a new approach.

Multimarker Transmission/Disequilibrium Tests (TDTs) are very robust association tests to population admixture and structure which may be used to identify susceptibility loci in genome-wide association studies. Multimarker TDTs using several markers may increase power by capturing high-degree associations. However, there is also a risk of spurious associations and power reduction due to the increase in degrees of freedom. In this study we show that associations found by tests built on simple null hypotheses are highly reproducible in a second independent data set regardless the number of markers. As a test exhibiting this feature to its maximum, we introduce the multimarker 2-Groups TDT (mTDT(2G)), a test which under the hypothesis of no linkage, asymptotically follows a χ2 distribution with 1 degree of freedom regardless the number of markers. The statistic requires the division of parental haplotypes into two groups: disease susceptibility and disease protective haplotype groups. We assessed the test behavior by performing an extensive simulation study as well as a real-data study using several data sets of two complex diseases. We show that mTDT(2G) test is highly efficient and it achieves the highest power among all the tests used, even when the null hypothesis is tested in a second independent data set. Therefore, mTDT(2G) turns out to be a very promising multimarker TDT to perform genome-wide searches for disease susceptibility loci that may be used as a preprocessing step in the construction of more accurate genetic models to predict individual susceptibility to complex diseases.

Replication study of 10 genes showing evidence for association with multiple sclerosis: validation of TMEM39A, IL12B and CBLB [correction of CLBL] genes.

BACKGROUND AND OBJECTIVES: Ten genes previously showing different evidence of association with multiple sclerosis have been selected to validate. METHODS: Eleven polymorphisms were genotyped with the iPLEX™ Sequenom in a well-powered collection of Spanish origin including 2863 multiple sclerosis cases and 2930 controls. RESULTS: Replication extended to the following polymorphisms: PKN2 (rs305217), GTF2B (rs7538427), EPHA4 (rs1517440), YTHDF3 (rs12115114), ANKFN1 (rs17758761) and PTPRM (rs4798571), which did not reach the threshold of significance in a follow-up of the first genome-wide association study (GWAS) conducted in multiple sclerosis; TMEM39A (rs1132200), which appeared as a newly identified susceptibility gene in the same study; a gene previously reaching GWAS significance in Italy, CBLB (rs9657904); IL12B (rs6887695, rs10045431), a susceptibility gene shared by diverse autoimmune diseases and, finally, another gene showing inconclusive association with multiple sclerosis, CNR1 (rs1049353). CONCLUSIONS: Pooled analysis corroborated the effect on MS predisposition of three genes: TMEM39A [rs1132200: p(M-H)=0.001; OR(M-H) (95% CI)= 0.84 (0.75-0.93)], IL12B [rs6887695: p(M-H)=0.03; OR(M-H) (95% CI)= 1.09 (1.01-1.17)] and CBLB [rs9657904: p(M-H)=0.01; OR(M-H) (95% CI)= 0.89 (0.81-0.97)].

CD40: novel association with Crohn's disease and replication in multiple sclerosis susceptibility.

BACKGROUND: A functional polymorphism located at -1 from the start codon of the CD40 gene, rs1883832, was previously reported to disrupt a Kozak sequence essential for translation. It has been consistently associated with Graves' disease risk in populations of different ethnicity and genetic proxies of this variant evaluated in genome-wide association studies have shown evidence of an effect in rheumatoid arthritis and multiple sclerosis (MS) susceptibility. However, the protective allele associated with Graves' disease or rheumatoid arthritis has shown a risk role in MS, an effect that we aimed to replicate in the present work. We hypothesized that this functional polymorphism might also show an association with other complex autoimmune condition such as inflammatory bowel disease, given the CD40 overexpression previously observed in Crohn's disease (CD) lesions. METHODOLOGY: Genotyping of rs1883832C>T was performed in 1564 MS, 1102 CD and 969 ulcerative colitis (UC) Spanish patients and in 2948 ethnically matched controls by TaqMan chemistry. PRINCIPAL FINDINGS: The observed effect of the minor allele rs1883832T was replicated in our independent Spanish MS cohort [p = 0.025; OR (95% CI) = 1.12 (1.01-1.23)]. The frequency of the minor allele was also significantly higher in CD patients than in controls [p = 0.002; OR (95% CI) = 1.19 (1.06-1.33)]. This increased predisposition was not detected in UC patients [p = 0.5; OR (95% CI) = 1.04 (0.93-1.17)]. CONCLUSION: The impact of CD40 rs1883832 on MS and CD risk points to a common signaling shared by these autoimmune conditions.

Validation of the CD6 and TNFRSF1A loci as risk factors for multiple sclerosis in Spain.

A recent meta-analysis of genome-wide association screens coupled to a replication exercise in a combined US/UK collection led to the identification of 4 single nucleotide polymorphisms (SNPs) in three gene loci, i.e. TNFRSF1A, CD6 and IRF8, as novel risk factors for multiple sclerosis with genome-wide level of significance. In the present study, using a combined all-Spain collection of 2515 MS patients and 2942 healthy controls, we demonstrate significant association of rs17824933 in CD6 (P(CMH)=0.004; OR=1.14; 95% CI 1.04-1.24) and of rs1860545 in TNFRSF1A (P(CMH)=0.001; OR=1.15; 95% CI 1.06-1.25) with MS, while the low-frequency coding non-synonymous SNP rs4149584 in TNFRSF1A displayed a trend for association (P(CMH)=0.062; OR=1.27; 95% CI 0.99-1.63). This data reinforce a generic role for CD6 and TNFRSF1A in susceptibility to MS, extending to populations of southern European ancestry.

Characterization of the MAL2-positive compartment in oligodendrocytes.

Oligodendrocytes (OLs), the myelin-producing cells of the central nervous system, segregate different surface subdomains at the plasma membrane as do other differentiated cells such as polarized epithelia and neurons. To generate the complex membrane system that characterizes myelinating OLs, large amounts of membrane proteins and lipids need to be synthesized and correctly targeted. In polarized epithelia, a considerable fraction of apical proteins are transported by an indirect pathway involving a detour to the basolateral membrane before being internalized and transported across the cell to the apical membrane by a process known as transcytosis. The apical recycling endosome (ARE) or its equivalent, the subapical compartment (SAC), of hepatocytes is an intracellular trafficking station involved in the transcytotic pathway. MAL2, an essential component of the machinery for basolateral-to-apical transcytosis, is an ARE/SAC resident protein. Here, we show that, after differentiation, murine oligodendrocyte precursor and human oligodendroglioma derived cell lines, Oli-neu and HOG, respectively, up-regulate the expression of MAL2 and accumulate it in an intracellular compartment, exhibiting a peri-centrosomal localization. In these oligodendrocytic cell lines, this compartment shares some of the main features of the ARE/SAC, such as colocalization with Rab11a, sensitivity to disruption of the microtubule cytoskeleton with nocodazole, and lack of internalized transferrin. Therefore, we suggest that the MAL2-positive compartment in oligodendrocytic cells could be a structure analogous to the ARE/SAC and might have an important role in the sorting of proteins and lipids for myelin assembly during oligodendrocyte differentiation.

Effects of the multiple sclerosis associated -330 promoter polymorphism in IL2 allelic expression.

The -330 IL2 gene promoter polymorphism has been associated with multiple sclerosis (MS) [J. Neuroimmunol. 119 (2001) 101], but the basis underlying this association remains unknown to date. In the present work, we have found that IL2 promoter-luciferase constructs, transfected in Jurkat cell line, showed twofold higher levels of gene expression in the -330 G allele. However, the transcriptional effect of this polymorphism in lymphocytes showed that the G allele was related to lower expression of IL2. This difference increased in the patient group. Divergence between in vivo and in vitro influence of the -330 IL2 promoter polymorphic site suggests the existence of additional unknown polymorphisms affecting gene regulation. Our data show an increased IL2 expression among GT and TT genotypes previously associated with susceptibility to MS.