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.

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.

Genome-wide association studies (GWAS) in migraine are providing the molecular basis of this heterogeneous disease, but the understanding of its aetiology is still incomplete. Although some biomarkers have currently been accepted for migraine, large amount of studies for identifying new ones is needed. The migraine-associated variant rs12355831:A>G (P=2 × 10−6), described in a GWAS of the International Headache Genetic Consortium, is localized in a non-coding sequence with unknown function. We sought to identify the causal variant and the genetic mechanism involved in the migraine risk. To this end, we integrated data of RNA sequences from the Genetic European Variation in Health and Disease (GEUVADIS) and genotypes from 1000 GENOMES of 344 lymphoblastoid cell lines (LCLs), to determine the expression quantitative trait loci (eQTLs) in the region. We found that the migraine-associated variant belongs to a linkage disequilibrium block associated with the expression of an acyl-coenzyme A synthetase 5 (ACSL5) transcript lacking exon 20 (ACSL5-Δ20). We showed by exon-skipping assay a direct causality of rs2256368-G in the exon 20 skipping of approximately 20 to 40% of ACSL5 RNA molecules. In conclusion, we identified the functional variant (rs2256368:A>G) affecting ACSL5 exon 20 skipping, as a causal factor linked to the migraine-associated rs12355831:A>G, suggesting that the activation of long-chain fatty acids by the spliced ACSL5-Δ20 molecules, a mitochondrial located enzyme, is involved in migraine pathology.