No evidence for association between rs10191329 severity locus and longitudinal disease severity in 1813 relapse-onset multiple sclerosis patients from the MSBase registry

BACKGROUND

The International Multiple Sclerosis Genetics Consortium and MultipleMS Consortium recently reported a genetic variant associated with multiple sclerosis (MS) severity. However, it remains unclear if these variants remain associated with more robust, longitudinal measures of disease severity.

METHODS

We examined the top variant, rs10191329, from Harroud et al.'s study in 1813 relapse-onset MS patients from the MSBase Registry to assess association with longitudinal disease severity.

RESULTS

Our analysis revealed no significant association between rs10191329 genotype and longitudinal binary disease severity (p > 0.05).

CONCLUSION

These findings highlight the complexity of genetic factors mediating long-term MS outcomes and the need for further research.

Multifaceted Analysis of Cerebrospinal Fluid and Serum from Progressive Multiple Sclerosis Patients: Potential Role of Vitamin C and Metal Ion Imbalance in the Divergence of Primary Progressive Multiple Sclerosis and Secondary Progressive Multiple Sclerosis

The progressive forms of multiple sclerosis (MS) primary progressive MS (PPMS) and secondary progressive MS (SPMS) are clinically distinguished by the rate at which symptoms worsen. Little is however known about the pathological mechanisms underlying the differential rate of accumulation of pathological changes. In this study, ¹H NMR spectroscopy was used to measure low-molecular-weight metabolites in paired cerebrospinal fluid (CSF) and serum of PPMS, SPMS, and control patients, as well as to determine lipoproteins and glycoproteins in serum samples. Additionally, neurodegenerative and inflammatory markers, neurofilament light (NFL) and chitinase-3-like protein 1 (CHI3L1), and the concentration of seven metal elements, Mg, Mn, Cu, Fe, Pb, Zn, and Ca, were also determined in both CSF and serum. The results indicate that the pathological changes associated with progressive MS are mainly localized in the central nervous system (CNS). More so, PPMS and SPMS patients with comparable disability status are pathologically similar in relation to neurodegeneration, neuroinflammation, and some metabolites that distinguish them from controls. However, the rapid progression of PPMS from the onset may be driven by a combination of neurotoxicity induced by heavy metals coupled with diminished CNS antioxidative capacity associated with differential intrathecal ascorbate retention and imbalance of Mg and Cu.

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 1,813 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−08), the other in males (rs698805; βmale = -1.5395, P = 4.35 × 10−08), providing some evidence for sex dimorphism in multiple sclerosis severity. Tissue enrichment and pathway analyses identified an overrepresentation of genes expressed in central nervous system 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.

Identification of the genetic mechanism that associates L3MBTL3 to multiple sclerosis

Multiple sclerosis (MS) is a complex and demyelinating disease of the central nervous system. One of the challenges of the post-GWAS era is to understand the molecular basis of statistical associations to reveal gene networks and potential therapeutic targets. The L3MBTL3 locus has been associated with MS risk by GWAS. To identify the causal variant of the locus, we performed fine mapping in a cohort of 3440 MS patients and 1688 healthy controls. The variant that best explained the association was rs6569648 (P = 4.13E-10, OR = 0.71, 95% CI = 0.64-0.79), which tagged rs7740107, located in intron 7 of L3MBTL3. The rs7740107 (A/T) variant has been reported to be the best expression and splice quantitative trait locus (eQTL and sQTL) of the region in up to 35 human GTEx tissues. By sequencing RNA from blood of 17 MS patients and quantification by digital qPCR, we determined that this eQTL/sQTL originated from the expression of a novel short transcript starting in intron 7 near rs7740107. The short transcript was translated into three proteins starting at different translation initiation codons. These N-terminal truncated proteins lacked the region where L3MBTL3 interacts with the transcriptional regulator RBPJ (Recombination Signal Binding Protein for Immunoglobulin Kappa J Region) which, in turn, regulates the Notch signaling pathway. Our data and other functional studies suggest that the genetic mechanism underlying the MS association of rs7740107 affects not only the expression of L3MBTL3 isoforms, but might also involve the Notch signaling pathway.

NLRP3 inflammasome as prognostic factor and therapeutic target in primary progressive multiple sclerosis patients

Primary progressive multiple sclerosis is a poorly understood disease entity with no specific prognostic biomarkers and scarce therapeutic options. We aimed to identify disease activity biomarkers in multiple sclerosis by performing an RNA sequencing approach in peripheral blood mononuclear cells from a discovery cohort of 44 untreated patients with multiple sclerosis belonging to different clinical forms and activity phases of the disease, and 12 healthy control subjects. A validation cohort of 58 patients with multiple sclerosis and 26 healthy control subjects was included in the study to replicate the RNA sequencing findings. The RNA sequencing revealed an interleukin 1 beta (IL1B) signature in patients with primary progressive multiple sclerosis. Subsequent immunophenotyping pointed to blood monocytes as responsible for the IL1B signature observed in this group of patients. Functional experiments at baseline measuring apoptosis-associated speck-like protein containing a CARD (ASC) speck formation showed that the NOD-leucine rich repeat and pyrin containing protein 3 (NLRP3) inflammasome was overactive in monocytes from patients with primary progressive multiple sclerosis, and canonical NLRP3 inflammasome activation with a combination of ATP plus lipopolysaccharide was associated with increased IL1B production in this group of patients. Primary progressive multiple sclerosis patients with high IL1B gene expression levels in peripheral blood mononuclear cells progressed significantly faster compared to patients with low IL1B levels based on the time to reach an EDSS of 6.0 and the Multiple Sclerosis Severity Score. In agreement with peripheral blood findings, both NLRP3 and IL1B expression in brain tissue from patients with primary progressive multiple sclerosis was mainly restricted to cells of myeloid lineage. Treatment of mice with a specific NLRP3 inflammasome inhibitor attenuated established experimental autoimmune encephalomyelitis disease severity and improved CNS histopathology. NLRP3 inflammasome-specific inhibition was also effective in reducing axonal damage in a model of lipopolysaccharide-neuroinflammation using organotypic cerebellar cultures. Altogether, these results point to a role of IL1B and the NLRP3 inflammasome as prognostic biomarker and potential therapeutic target, respectively, in patients with primary progressive multiple sclerosis.

Targeted resequencing reveals rare variants enrichment in multiple sclerosis susceptibility genes

Although genome-wide association studies have identified a number of common variants associated with multiple sclerosis (MS) susceptibility, little is known about the relevance of rare variants. Here, we aimed to explore the role of rare variants in 14 MS risk genes (FCRL1, RGS1, TIMMDC1, HHEX, CXCR5, LTBR, TSFM, GALC, TRAF3, STAT3, TNFSF14, IFI30, CD40, and CYP24A1) by targeted resequencing in an Iberian population of 524 MS cases and 546 healthy controls. Four rare variants-enriched regions within CYP24A1, FCRL1, RGS1, and TRAF3 were identified as significantly associated with MS. Functional studies revealed significantly decreased regulator of G protein signaling 1 (RGS1) gene expression levels in peripheral blood mononuclear cells from MS patients with RGS1 rare variants compared to noncarriers, whereas no significant differences in gene expression were observed for CYP24A1, FCRL1, and TRAF3 between rare variants carriers and noncarriers. Immunophenotyping showed significant decrease in RGS1 expression in peripheral blood B lymphocytes from MS patients with RGS1 rare variants relative to noncarriers. Lastly, peripheral blood mononuclear cell from MS patients carrying RGS1 rare variants showed significantly lower induction of RGS1 gene expression by interferon-β compared to MS patients lacking RGS1 variants. The presence of rare variants in RGS1 reinforce the ideas of high genetic heterogeneity and a role of rare variants in MS pathogenesis.

A New Risk Variant for Multiple Sclerosis at 11q23.3 Locus Is Associated with Expansion of CXCR5+ Circulating Regulatory T Cells

Genome-wide association studies and meta-analysis have contributed to the identification of more than 200 loci associated with multiple sclerosis (MS). However, a proportion of MS heritability remains unknown. We aimed to uncover new genetic variants associated with MS and determine their functional effects. For this, we resequenced the exons and regulatory sequences of 14 MS risk genes in a cohort of MS patients and healthy individuals (n = 1070) and attempted to validate a selection of signals through genotyping in an independent cohort (n = 5138). We identified three new MS-associated variants at C-X-C motif chemokine receptor 5 (CXCR5), Ts translation elongation factor, mitochondrial (TSFM) and cytochrome P450 family 24 subfamily A member 1 (CYP24A1). Rs10892307 resulted in a new signal at the CXCR5 region that explains one of the associations with MS within the locus. This polymorphism and three others in high linkage disequilibrium mapped within regulatory regions. Of them, rs11602393 showed allele-dependent enhancer activity in the forward orientation as determined by luciferase reporter assays. Immunophenotyping using peripheral blood mononuclear cells from MS patients associated the minor allele of rs10892307 with increased percentage of regulatory T cells expressing CXCR5. This work reports a new signal for the CXCR5 MS risk locus and points to rs11602393 as the causal variant. The expansion of CXCR5+ circulating regulatory T cells induced by this variant could cause its MS association.

The Rare IL22RA2 Signal Peptide Coding Variant rs28385692 Decreases Secretion of IL-22BP Isoform-1, -2 and -3 and Is Associated with Risk for Multiple Sclerosis

The IL22RA2 locus is associated with risk for multiple sclerosis (MS) but causative variants are yet to be determined. In a single nucleotide polymorphism (SNP) screen of this locus in a Basque population, rs28385692, a rare coding variant substituting Leu for Pro at position 16 emerged significantly (p = 0.02). This variant is located in the signal peptide (SP) shared by the three secreted protein isoforms produced by IL22RA2 (IL-22 binding protein-1(IL-22BPi1), IL-22BPi2 and IL-22BPi3). Genotyping was extended to a Europe-wide case-control dataset and yielded high significance in the full dataset (p = 3.17 × 10^-4). Importantly, logistic regression analyses conditioning on the main known MS-associated SNP at this locus, rs17066096, revealed that this association was independent from the primary association signal in the full case-control dataset. In silico analysis predicted both disruption of the alpha helix of the H-region of the SP and decreased hydrophobicity of this region, ultimately affecting the SP cleavage site. We tested the effect of the p.Leu16Pro variant on the secretion of IL-22BPi1, IL-22BPi2 and IL-22BPi3 and observed that the Pro16 risk allele significantly lowers secretion levels of each of the isoforms to around 50%-60% in comparison to the Leu16 reference allele. Thus, our study suggests that genetically coded decreased levels of IL-22BP isoforms are associated with augmented risk for MS.

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.

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.

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.

Exome sequencing study in patients with multiple sclerosis reveals variants associated with disease course

Background

It remains unclear whether disease course in multiple sclerosis (MS) is influenced by genetic polymorphisms. Here, we aimed to identify genetic variants associated with benign and aggressive disease courses in MS patients.

Methods

MS patients were classified into benign and aggressive phenotypes according to clinical criteria. We performed exome sequencing in a discovery cohort, which included 20 MS patients, 10 with benign and 10 with aggressive disease course, and genotyping in 2 independent validation cohorts. The first validation cohort encompassed 194 MS patients, 107 with benign and 87 with aggressive phenotypes. The second validation cohort comprised 257 patients, of whom 224 patients had benign phenotypes and 33 aggressive disease courses. Brain immunohistochemistries were performed using disease course associated genes antibodies.

Results

By means of single-nucleotide polymorphism (SNP) detection and comparison of allele frequencies between patients with benign and aggressive phenotypes, a total of 16 SNPs were selected for validation from the exome sequencing data in the discovery cohort. Meta-analysis of genotyping results in two validation cohorts revealed two polymorphisms, rs28469012 and rs10894768, significantly associated with disease course. SNP rs28469012 is located in CPXM2 (carboxypeptidase X, M14 family, member 2) and was associated with aggressive disease course (uncorrected p value < 0.05). SNP rs10894768, which is positioned in IGSF9B (immunoglobulin superfamily member 9B) was associated with benign phenotype (uncorrected p value < 0.05). In addition, a trend for association with benign phenotype was observed for a third SNP, rs10423927, in NLRP9 (NLR family pyrin domain containing 9). Brain immunohistochemistries in chronic active lesions from MS patients revealed expression of IGSF9B in astrocytes and macrophages/microglial cells, and expression of CPXM2 and NLRP9 restricted to brain macrophages/microglia.

Conclusions

Genetic variants located in CPXM2, IGSF9B, and NLRP9 have the potential to modulate disease course in MS patients and may be used as disease activity biomarkers to identify patients with divergent disease courses. Altogether, the reported results from this study support the influence of genetic factors in MS disease course and may help to better understand the complex molecular mechanisms underlying disease pathogenesis.

Electronic supplementary material

The online version of this article (10.1186/s12974-018-1307-1) contains supplementary material, which is available to authorized users.

A comparison of genomic profiles of complex diseases under different models

Background

Various approaches are being used to predict individual risk to polygenic diseases from data provided by genome-wide association studies. As there are substantial differences between the diseases investigated, the data sets used and the way they are tested, it is difficult to assess which models are more suitable for this task.

Results

We compared different approaches for seven complex diseases provided by the Wellcome Trust Case Control Consortium (WTCCC) under a within-study validation approach. Risk models were inferred using a variety of learning machines and assumptions about the underlying genetic model, including a haplotype-based approach with different haplotype lengths and different thresholds in association levels to choose loci as part of the predictive model. In accordance with previous work, our results generally showed low accuracy considering disease heritability and population prevalence. However, the boosting algorithm returned a predictive area under the ROC curve (AUC) of 0.8805 for Type 1 diabetes (T1D) and 0.8087 for rheumatoid arthritis, both clearly over the AUC obtained by other approaches and over 0.75, which is the minimum required for a disease to be successfully tested on a sample at risk, which means that boosting is a promising approach. Its good performance seems to be related to its robustness to redundant data, as in the case of genome-wide data sets due to linkage disequilibrium.

Conclusions

In view of our results, the boosting approach may be suitable for modeling individual predisposition to Type 1 diabetes and rheumatoid arthritis based on genome-wide data and should be considered for more in-depth research.

Electronic supplementary material

The online version of this article (doi:10.1186/s12920-015-0157-2) contains supplementary material, which is available to authorized users.

The multiple sclerosis-associated regulatory variant rs10877013 affects expression of CYP27B1 and VDR under inflammatory or vitamin D stimuli

Background:

Vitamin D deficit is considered an important risk factor for many inflammatory and autoimmune diseases.

Objective:

To investigate the influence of the multiple sclerosis (MS)-associated regulatory variant rs10877013 on the expression of genes involved in vitamin D activation ( CYP27B1), vitamin D receptor ( VDR), and vitamin D degradation ( CYP24A1) under inflammatory environment or vitamin D.

Methods:

We used lipopolysaccharide and interferon-gamma (LPS+IFNγ) activated monocytes from 119 individuals and vitamin D-stimulated lymphoblastoid cell lines (LCLs, n = 109) of 1000 genomes to quantify the mRNA expression of vitamin D genes by quantitative reverse transcription polymerase chain reaction (RT-qPCR).

Results:

We found that CYP27B1 mRNA expression level was associated with the rs10877013 genotypes ( p = 5.0E-6) in LPS+IFNγ treated monocytes, but not in vitamin D-stimulated LCLs. Inversely, rs10877013 genotypes were associated with VDR expression in LCLs ( p = 6.0E-4) but not in monocytes. Finally, CYP24A1 was highly induced by the active form of vitamin D and its expression correlated with the expression of VDR in LCLs but neither the MS-associated variant in the region (rs2248359) nor any other variant located in 1 Mb around CYP24A1 was associated with its expression.

Conclusions:

The MS-associated variant rs10877013 is a genetic determinant that affects the functioning of the vitamin D system linking environmental and genetic factors.

Influence of the LILRA3 Deletion on Multiple Sclerosis Risk: Original Data and Meta-Analysis

Background

Multiple sclerosis (MS) is a neurodegenerative, autoimmune disease of the central nervous system. Genome-wide association studies (GWAS) have identified over hundred polymorphisms with modest individual effects in MS susceptibility and they have confirmed the main individual effect of the Major Histocompatibility Complex. Additional risk loci with immunologically relevant genes were found significantly overrepresented. Nonetheless, it is accepted that most of the genetic architecture underlying susceptibility to the disease remains to be defined. Candidate association studies of the leukocyte immunoglobulin-like receptor LILRA3 gene in MS have been repeatedly reported with inconsistent results.

Objectives

In an attempt to shed some light on these controversial findings, a combined analysis was performed including the previously published datasets and three newly genotyped cohorts. Both wild-type and deleted LILRA3 alleles were discriminated in a single-tube PCR amplification and the resulting products were visualized by their different electrophoretic mobilities.

Results and Conclusion

Overall, this meta-analysis involved 3200 MS patients and 3069 matched healthy controls and it did not evidence significant association of the LILRA3 deletion [carriers of LILRA3 deletion: p = 0.25, OR (95% CI) = 1.07 (0.95–1.19)], even after stratification by gender and the HLA-DRB1*15:01 risk allele.

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.

A new risk variant for multiple sclerosis at the immunoglobulin heavy chain locus associates with intrathecal IgG, IgM index and oligoclonal bands

Background:

Recent findings have shown a correlation between the intrathecal IgG index and variants at the immunoglobulin heavy chain (IGHC) locus in patients with multiple sclerosis (MS).

Objectives:

The objective of this paper is to analyse the association of the locus with MS susceptibility and its relationship with intrathecal immunoglobulin (Ig) parameters.

Methods:

We genotyped the rs11621145 variant, located at the IGHC locus, in 2726 patients with MS and 2133 healthy controls. Associations of intrathecal IgG and IgM indexes with rs11621145 were analysed by linear regression analysis in 538 MS patients.

Results:

We found that rs11621145 showed statistically significant evidence for association with susceptibility to MS (odds ratio = 0.69, p = 1.053E-09), though validation of this result in additional cohorts would be desirable. We confirmed the association between the IgG index and the rs11621145 ( p = 6.85E-07, Beta = 0.207). Furthermore, rs11621145 was inversely correlated with IgM index ( p = 7.24E-04, Beta = -0.277), and therefore marks a decreased likelihood of presenting IgM oligoclonal bands (odds ratio = 0.38, p = 2.35E-06).

Conclusions:

Our results suggest that the polymorphism of the IGHC locus could be altering the switching of the Ig isotype in B cells and it may be interfering with T-dependent and T-independent antibody responses.

HERV-W polymorphism in chromosome X is associated with multiple sclerosis risk and with differential expression of MSRV

Background

Multiple Sclerosis (MS) is an autoimmune demyelinating disease that occurs more frequently in women than in men. Multiple Sclerosis Associated Retrovirus (MSRV) is a member of HERV-W, a multicopy human endogenous retroviral family repeatedly implicated in MS pathogenesis. MSRV envelope protein is elevated in the serum of MS patients and induces inflammation and demyelination but, in spite of this pathogenic potential, its exact genomic origin and mechanism of generation are unknown. A possible link between the HERV-W copy on chromosome Xq22.3, that contains an almost complete open reading frame, and the gender differential prevalence in MS has been suggested.

Results

MSRV transcription levels were higher in MS patients than in controls (U-Mann–Whitney; p = 0.004). Also, they were associated with the clinical forms (Spearman; p = 0.0003) and with the Multiple Sclerosis Severity Score (MSSS) (Spearman; p = 0.016). By mapping a 3 kb region in Xq22.3, including the HERV-W locus, we identified three polymorphisms: rs6622139 (T/C), rs6622140 (G/A) and rs1290413 (G/A). After genotyping 3127 individuals (1669 patients and 1458 controls) from two different Spanish cohorts, we found that in women rs6622139 T/C was associated with MS susceptibility: [χ²; p = 0.004; OR (95% CI) = 0.50 (0.31-0.81)] and severity, since CC women presented lower MSSS scores than CT (U-Mann–Whitney; p = 0.039) or TT patients (U-Mann–Whitney; p = 0.031). Concordantly with the susceptibility conferred in women, rs6622139*T was associated with higher MSRV expression (U-Mann–Whitney; p = 0.003).

Conclusions

Our present work supports the hypothesis of a direct involvement of HERV-W/MSRV in MS pathogenesis, identifying a genetic marker on chromosome X that could be one of the causes underlying the gender differences in MS.

MANBA, CXCR5, SOX8, RPS6KB1 and ZBTB46 are genetic risk loci for multiple sclerosis

A recent genome-wide association study reported five loci for which there was strong, but sub-genome-wide significant evidence for association with multiple sclerosis risk. The aim of this study was to evaluate the role of these potential risk loci in a large and independent data set of ≈ 20,000 subjects. We tested five single nucleotide polymorphisms rs228614 (MANBA), rs630923 (CXCR5), rs2744148 (SOX8), rs180515 (RPS6KB1), and rs6062314 (ZBTB46) for association with multiple sclerosis risk in a total of 8499 cases with multiple sclerosis, 8765 unrelated control subjects and 958 trios of European descent. In addition, we assessed the overall evidence for association by combining these newly generated data with the results from the original genome-wide association study by meta-analysis. All five tested single nucleotide polymorphisms showed consistent and statistically significant evidence for association with multiple sclerosis in our validation data sets (rs228614: odds ratio = 0.91, P = 2.4 × 10(-6); rs630923: odds ratio = 0.89, P = 1.2 × 10(-4); rs2744148: odds ratio = 1.14, P = 1.8 × 10(-6); rs180515: odds ratio = 1.12, P = 5.2 × 10(-7); rs6062314: odds ratio = 0.90, P = 4.3 × 10(-3)). Combining our data with results from the previous genome-wide association study by meta-analysis, the evidence for association was strengthened further, surpassing the threshold for genome-wide significance (P < 5 × 10(-8)) in each case. Our study provides compelling evidence that these five loci are genuine multiple sclerosis susceptibility loci. These results may eventually lead to a better understanding of the underlying disease pathophysiology.

Closing the case of APOE in multiple sclerosis: no association with disease risk in over 29 000 subjects

BACKGROUND

Single nucleotide polymorphisms (SNPs) rs429358 (ε4) and rs7412 (ε2), both invoking changes in the amino-acid sequence of the apolipoprotein E (APOE) gene, have previously been tested for association with multiple sclerosis (MS) risk. However, none of these studies was sufficiently powered to detect modest effect sizes at acceptable type-I error rates. As both SNPs are only imperfectly captured on commonly used microarray genotyping platforms, their evaluation in the context of genome-wide association studies has been hindered until recently.

METHODS

We genotyped 12 740 subjects hitherto not studied for their APOE status, imputed raw genotype data from 8739 subjects from five independent genome-wide association studies datasets using the most recent high-resolution reference panels, and extracted genotype data for 8265 subjects from previous candidate gene assessments.

RESULTS

Despite sufficient power to detect associations at genome-wide significance thresholds across a range of ORs, our analyses did not support a role of rs429358 or rs7412 on MS susceptibility. This included meta-analyses of the combined data across 13 913 MS cases and 15 831 controls (OR=0.95, p=0.259, and OR 1.07, p=0.0569, for rs429358 and rs7412, respectively).

CONCLUSION

Given the large sample size of our analyses, it is unlikely that the two APOE missense SNPs studied here exert any relevant effects on MS susceptibility.

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.

Fine mapping and functional analysis of the multiple sclerosis risk gene CD6

CD6 has recently been identified and validated as risk gene for multiple sclerosis (MS), based on the association of a single nucleotide polymorphism (SNP), rs17824933, located in intron 1. CD6 is a cell surface scavenger receptor involved in T-cell activation and proliferation, as well as in thymocyte differentiation. In this study, we performed a haptag SNP screen of the CD6 gene locus using a total of thirteen tagging SNPs, of which three were non-synonymous SNPs, and replicated the recently reported GWAS SNP rs650258 in a Spanish-Basque collection of 814 controls and 823 cases. Validation of the six most strongly associated SNPs was performed in an independent collection of 2265 MS patients and 2600 healthy controls. We identified association of haplotypes composed of two non-synonymous SNPs [rs11230563 (R225W) and rs2074225 (A257V)] in the 2^nd SRCR domain with susceptibility to MS (P_{max(T) permutation} = 1×10⁻⁴). The effect of these haplotypes on CD6 surface expression and cytokine secretion was also tested. The analysis showed significantly different CD6 expression patterns in the distinct cell subsets, i.e. – CD4⁺ naïve cells, P = 0.0001; CD8⁺ naïve cells, P<0.0001; CD4⁺ and CD8⁺ central memory cells, P = 0.01 and 0.05, respectively; and natural killer T (NKT) cells, P = 0.02; with the protective haplotype (RA) showing higher expression of CD6. However, no significant changes were observed in natural killer (NK) cells, effector memory and terminally differentiated effector memory T cells. Our findings reveal that this new MS-associated CD6 risk haplotype significantly modifies expression of CD6 on CD4⁺ and CD8⁺ T cells.

Identification of a functional variant in the KIF5A-CYP27B1-METTL1-FAM119B locus associated with multiple sclerosis

Background and aim

Several studies have highlighted the association of the 12q13.3–12q14.1 region with coeliac disease, type 1 diabetes, rheumatoid arthritis and multiple sclerosis (MS); however, the causal variants underlying diseases are still unclear. The authors sought to identify the functional variant of this region associated with MS.

Methods

Tag-single nucleotide polymorphism (SNP) analysis of the associated region encoding 15 genes was performed in 2876 MS patients and 2910 healthy Caucasian controls together with expression regulation analyses.

Results

rs6581155, which tagged 18 variants within a region where 9 genes map, was sufficient to model the association. This SNP was in total linkage disequilibrium (LD) with other polymorphisms that associated with the expression levels of FAM119B, AVIL, TSFM, TSPAN31 and CYP27B1 genes in different expression quantitative trait loci studies. Functional annotations from Encyclopedia of DNA Elements (ENCODE) showed that six out of these rs6581155-tagged-SNPs were located in regions with regulatory potential and only one of them, rs10877013, exhibited allele-dependent (ratio A/G=9.5-fold) and orientation-dependent (forward/reverse=2.7-fold) enhancer activity as determined by luciferase reporter assays. This enhancer is located in a region where a long-range chromatin interaction among the promoters and promoter-enhancer of several genes has been described, possibly affecting their expression simultaneously.

Conclusions

This study determines a functional variant which alters the enhancer activity of a regulatory element in the locus affecting the expression of several genes and explains the association of the 12q13.3–12q14.1 region with MS.

In-vivo and in-vitro analysis of retroviral vectors carrying the N-ras oncogene

We have analyzed, in vivo and in vitro, the behavior of two retroviral vectors carrying the genomic or cDNA N-ras oncogene to study the role of N-ras in the initiation and development of thymic lymphomas. The vector bearing the genomic gene produced an array of transcripts originating from the LTR and the oncogene promoter. The majority of the transcripts initiated at the LTR did not carry the packaging signal producing low titer clones. The cDNA vector produced two transcripts correctly spliced and the titers obtained were as high as 10(6) pfu/ml. Bone marrow cell infection and grafting of lethally irradiated mice was performed. The integrated vector in blood cells was followed at different times, observing that the provirus can disappear and reappear in peripheral blood cells during the course of the experiment. This observation fits with the hypothesis of clonal contribution of small number of stem cells in the renewal of blood cells. No tumors were detected in the infected animals, probably due to low expression of the integrated provirus. These experiments provide information on the advantages and disadvantages of genomic versus cDNA constructs in retroviral vectors.

Increasing power by using haplotype similarity in a multimarker transmission/disequilibrium test

It is already known that power in multimarker transmission/disequilibrium tests may improve with the number of markers as some associations may require several markers to be captured. However, a mechanism such as haplotype grouping must be used to avoid incremental complexity with the number of markers. 2G, a state-of-the-art transmission/disequilibrium test, implements this mechanism to its maximum extent by grouping haplotypes into only two groups, high and low-risk haplotypes, so that the test has only one degree of freedom regardless of the number of markers. The test checks whether those haplotypes more often transmitted from parents to offspring are truly high-risk haplotypes. In this paper we use haplotype similarity as prior knowledge to classify haplotypes as high or low risk ones and start with those haplotypes in which the prior will have lower impact i.e. those with the largest differences between transmission and non-transmission counts. If their counts are very different, the prior knowledge has little effect and haplotypes are classified as low or high risk as 2G does. We show a substantial gain in power achieved by this approach, in both simulation and real data sets.

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-Groups TDT (), a test which under the hypothesis of no linkage, asymptotically follows a distribution with 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 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, 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.

Multiple sclerosis risk variant HLA-DRB1*1501 associates with high expression of DRB1 gene in different human populations

The human leukocyte antigen (HLA) DRB1*1501 has been consistently associated with multiple sclerosis (MS) in nearly all populations tested. This points to a specific antigen presentation as the pathogenic mechanism though this does not fully explain the disease association. The identification of expression quantitative trait loci (eQTL) for genes in the HLA locus poses the question of the role of gene expression in MS susceptibility. We analyzed the eQTLs in the HLA region with respect to MS-associated HLA-variants obtained from genome-wide association studies (GWAS). We found that the Tag of DRB1*1501, rs3135388 A allele, correlated with high expression of DRB1, DRB5 and DQB1 genes in a Caucasian population. In quantitative terms, the MS-risk AA genotype carriers of rs3135388 were associated with 15.7-, 5.2- and 8.3-fold higher expression of DQB1, DRB5 and DRB1, respectively, than the non-risk GG carriers. The haplotype analysis of expression-associated variants in a Spanish MS cohort revealed that high expression of DRB1 and DQB1 alone did not contribute to the disease. However, in Caucasian, Asian and African American populations, the DRB1*1501 allele was always highly expressed. In other immune related diseases such as type 1 diabetes, inflammatory bowel disease, ulcerative colitis, asthma and IgA deficiency, the best GWAS-associated HLA SNPs were also eQTLs for different HLA Class II genes. Our data suggest that the DR/DQ expression levels, together with specific structural properties of alleles, seem to be the causal effect in MS and in other immunopathologies rather than specific antigen presentation alone.

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)].

Chitinase 3-like 1 plasma levels are increased in patients with progressive forms of multiple sclerosis

Background:

Chitinase 3-like 1 (CHI3L1) is upregulated in a wide variety of inflammatory conditions. Recent studies have pointed to a role of CHI3L1 in multiple sclerosis (MS) pathogenesis.

Objective:

The objective of this study was to investigate the role of plasma CHI3L1 in MS clinical course and disease activity and to evaluate the effect of interferon-beta (IFNβ) treatment on protein levels.

Methods:

Plasma CHI3L1 levels were determined by ELISA in 57 healthy controls (HC), 220 untreated MS patients [66 primary progressive MS patients (PPMS), 30 secondary progressive MS patients (SPMS), and 124 relapsing–remitting MS patients (RRMS), 94 during clinical remission and 30 during relapse], and 32 MS patients receiving IFNβ treatment. A polymorphism of the CHI3L1 gene, rs4950928, was genotyped in 3274 MS patients and 3483 HC.

Results:

Plasma CHI3L1 levels were significantly increased in patients with progressive forms of MS compared with RRMS patients and HC. CHI3L1 levels were similar between RRMS patients in relapse and remission. A trend towards decreased CHI3L1 levels was observed in IFNβ-treated patients. Allele C of rs4950928 was significantly associated with PPMS patients and with higher plasma CHI3L1 levels.

Conclusions:

These findings point to a role of CHI3L1 in patients with progressive forms of MS, particularly in those with PPMS.

High ACSL5 transcript levels associate with systemic lupus erythematosus and apoptosis in Jurkat T lymphocytes and peripheral blood cells

BACKGROUND

Systemic lupus erythematosus (SLE) is a prototypical autoimmune disease in which increased apoptosis and decreased apoptotic cells removal has been described as most relevant in the pathogenesis. Long-chain acyl-coenzyme A synthetases (ACSLs) have been involved in the immunological dysfunction of mouse models of lupus-like autoimmunity and apoptosis in different in vitro cell systems. The aim of this work was to assess among the ACSL isoforms the involvement of ACSL2, ACSL4 and ACSL5 in SLE pathogenesis.

FINDINGS

With this end, we determined the ACSL2, ACSL4 and ACSL5 transcript levels in peripheral blood mononuclear cells (PBMCs) of 45 SLE patients and 49 healthy controls by quantitative real time-PCR (q-PCR). We found that patients with SLE had higher ACSL5 transcript levels than healthy controls [median (range), healthy controls = 16.5 (12.3-18.0) vs. SLE = 26.5 (17.8-41.7), P = 3.9×10 E-5] but no differences were found for ACSL2 and ACSL4. In in vitro experiments, ACSL5 mRNA expression was greatly increased when inducing apoptosis in Jurkat T cells and PBMCs by Phorbol-Myristate-Acetate plus Ionomycin (PMA+Io). On the other hand, short interference RNA (siRNA)-mediated silencing of ACSL5 decreased induced apoptosis in Jurkat T cells up to the control levels as well as decreased mRNA expression of FAS, FASLG and TNF.

CONCLUSIONS

These findings indicate that ACSL5 may play a role in the apoptosis that takes place in SLE. Our results point to ACSL5 as a potential novel functional marker of pathogenesis and a possible therapeutic target in SLE.

Interaction of PLP with GFP-MAL2 in the human oligodendroglial cell line HOG

The velocity of the nerve impulse conduction of vertebrates relies on the myelin sheath, an electrically insulating layer that surrounds axons in both the central and peripheral nervous systems, enabling saltatory conduction of the action potential. Oligodendrocytes are the myelin-producing glial cells in the central nervous system. A deeper understanding of the molecular basis of myelination and, specifically, of the transport of myelin proteins, will contribute to the search of the aetiology of many dysmyelinating and demyelinating diseases, including multiple sclerosis. Recent investigations suggest that proteolipid protein (PLP), the major myelin protein, could reach myelin sheath by an indirect transport pathway, that is, a transcytotic route via the plasma membrane of the cell body. If PLP transport relies on a transcytotic process, it is reasonable to consider that this myelin protein could be associated with MAL2, a raft protein essential for transcytosis. In this study, carried out with the human oligodendrocytic cell line HOG, we show that PLP colocalized with green fluorescent protein (GFP)-MAL2 after internalization from the plasma membrane. In addition, both immunoprecipitation and immunofluorescence assays, indicated the existence of an interaction between GFP-MAL2 and PLP. Finally, ultrastructural studies demonstrated colocalization of GFP-MAL2 and PLP in vesicles and tubulovesicular structures. Taken together, these results prove for the first time the interaction of PLP and MAL2 in oligodendrocytic cells, supporting the transcytotic model of PLP transport previously suggested.

Genome-wide association filtering using a highly locus-specific transmission/disequilibrium test

Multimarker transmission/disequilibrium tests (TDTs) are powerful association and linkage tests used to perform genome-wide filtering in the search for disease susceptibility loci. In contrast to case/control studies, they have a low rate of false positives for population stratification and admixture. However, the length of a region found in association with a disease is usually very large because of linkage disequilibrium (LD). Here, we define a multimarker proportional TDT (mTDT ( P )) designed to improve locus specificity in complex diseases that has good power compared to the most powerful multimarker TDTs. The test is a simple generalization of a multimarker TDT in which haplotype frequencies are used to weight the effect that each haplotype has on the whole measure. Two concepts underlie the features of the metric: the 'common disease, common variant' hypothesis and the decrease in LD with chromosomal distance. Because of this decrease, the frequency of haplotypes in strong LD with common disease variants decreases with increasing distance from the disease susceptibility locus. Thus, our haplotype proportional test has higher locus specificity than common multimarker TDTs that assume a uniform distribution of haplotype probabilities. Because of the common variant hypothesis, risk haplotypes at a given locus are relatively frequent and a metric that weights partial results for each haplotype by its frequency will be as powerful as the most powerful multimarker TDTs. Simulations and real data sets demonstrate that the test has good power compared with the best tests but has remarkably higher locus specificity, so that the association rate decreases at a higher rate with distance from a disease susceptibility or disease protective locus.

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.

STAT3 locus in inflammatory bowel disease and multiple sclerosis susceptibility

STAT3 (signal transducer and activator of transcription 3) signaling is a critical component of Th17-dependent autoimmune processes. Genome-wide association studies (GWAS) have revealed the role of the STAT3 gene in inflammatory bowel disease (IBD) susceptibility, although confirmation in clinical subphenotypes is warranted. Mice with targeted deletion of Stat3 in T cells are resistant to experimental autoimmune encephalomyelitis, which is a multiple sclerosis (MS) model. Moreover, increased phosphorylated STAT3 was reported in T cells of patients evolving from clinically isolated syndrome to defined MS and in relapsing patients. These evidences led us to analyze the role of STAT3 in Crohn's disease (CD), ulcerative colitis (UC) and MS risk. Polymorphisms in the STAT3 region (rs3809758/rs744166/rs1026916/rs12948909) were genotyped and the inferred haplotypes were subsequently analyzed in 860 IBD and 1540 MS Spanish patients and 1720 ethnically matched controls. The haplotype conformed by the risk alleles of each polymorphism was significantly associated with both clinical phenotypes of IBD (CD: P=0.005, odds ratio 1.25, 95% confidence interval 1.06-1.46; and UC: P=0.002, odds ratio 1.19, 95% confidence interval 1.02-1.38). No evidence of association was detected for MS. The originally described association of IBD with STAT3 polymorphisms is corroborated for the two clinical phenotypes, CD and UC, in an independent population. A major role of this gene in MS seems unlikely.

Polymorphisms in the IL2, IL2RA and IL2RB genes in multiple sclerosis risk

Interleukin (IL)-2/IL-2R signalling promotes proliferation and survival of activated T cells and has an essential non-redundant role in the production of regulatory T cells. Associations with different autoimmune diseases of polymorphisms in a linkage disequilibrium block in which the IL2/IL21 genes map (4q27), and also in genes encoding the IL2RA and IL2RB subunits (located in 10p15 and 22q13, respectively), were identified through genome-wide studies. Polymorphisms in these three genes were studied in 430 multiple sclerosis (MS) patients and in 550 ethnically matched controls from Madrid (Spain). Replication and meta-analysis with results from an independent cohort of 771 MS patients and 759 controls from Andalucía (Spain) confirmed the association of polymorphisms in the IL2RA gene (P(Mantel-Haenszel,) odds ratio (OR)(M-H) (95% confidence interval, CI) for rs2104286: 0.0001, 0.75 (0.65-0.87); for rs11594656/rs35285258: 0.004, 1.19 (1.06-1.34); for rs41295061: 0.03, 0.77 (0.60-0.98)); showed a trend for association of the IL2/IL21 rs6822844 (P(M-H)=0.07, OR(M-H) (95% CI)=0.86 (0.73-1.01)), but did not corroborate the association for IL2RB. Regression analyses of the combined Spanish cohort revealed the independence of two IL2RA association signals: rs2104286 and rs11594656/rs35285258. The relevant role of the IL2RA gene on MS susceptibility adds support to its common effect on autoimmune risk and the suggestive association of IL2/IL21 warrants further investigation.

Tag-SNP analysis of the GFI1-EVI5-RPL5-FAM69 risk locus for multiple sclerosis

A recent genome-wide association study conducted by the International Multiple Sclerosis Genetic Consortium (IMSGC) identified, among others, a number of putative multiple sclerosis (MS) susceptibility variants at position 1p22. Twenty-one SNPs positively associated with MS were located at the GFI-EVI5-RPL5-FAM69A locus. In this study, we performed an analysis and fine mapping of this locus, genotyping eight Tag-SNPs in 732 MS patients and 974 controls from Spain. We observed an association with MS in three of eight Tag-SNPs: rs11804321 (P=0.008, OR=1.29; 95% CI=1.08-1.54), rs11808092 (P=0.048, OR=1.19; 95% CI=1.03-1.39) and rs6680578 (P=0.0082, OR=1.23; 95% CI=1.07-1.41). After correcting for multiple comparisons and using logistic regression analysis to test the addition of each SNP to the most associated SNPs, we observed that rs11804321 alone was sufficient to model the association. This Tag-SNP captures two SNPs in complete linkage disequilibrium (r(2)=1), both located within the 17th intron of the EVI5 gene. Our findings agree with the corresponding data of the recent IMSGC study and present new genetic evidence that points to EVI5 as a factor of susceptibility to MS.

Hexose-6-phosphate dehydrogenase: a new risk gene for multiple sclerosis

A recent genome-wide association study (GWAS) performed by the The Wellcome Trust Case-Control Consortium based on 12 374 nonsynonymous single-nucleotide polymorphisms (SNPs) provided evidence for several genes involved in multiple sclerosis (MS) susceptibility. In this study, we aimed at verifying the association of 19 SNPs with MS, with P-values < or =0.005, in an independent cohort of 732 patients and 974 controls, all Caucasian from the South of Spain. We observed an association of the rs17368528 polymorphism with MS (P=0.04, odds ratio (OR)=0.801, 95% confidence interval (CI)=0.648-0.990). The association of this polymorphism with MS was further validated in an independent set of 1318 patients from the Canadian Collaborative Project (P=0.04, OR=0.838, 95% CI=0.716-0.964). This marker is located on chromosome 1p36.22, which is 1 Mb away from the MS-associated kinesin motor protein KIF1B, although linkage disequilibrium was not observed between these two markers. The rs17368528 SNP results in an amino-acid substitution (proline to leucine) in the fifth exon of the hexose-6-phosphate dehydrogenase (H6PD) gene, in which some variants have been reported to attenuate or abolish H6PD activity, in individuals with cortisone reductase deficiency. This study corroborates the association of one locus determined by GWAS and points to H6PD as a new candidate gene for MS.

Multiple sclerosis association study with the TENR-IL2-IL21 region in a Spanish population

Polymorphisms from the TENR-IL2-IL21 block in the 4q27 chromosome were recently associated with type 1 diabetes, celiac disease, rheumatoid arthritis and psoriasis. We undertook this study to investigate the potential role of polymorphisms rs3136534, rs6822844 and rs2069762 (-330 T/G IL2) in multiple sclerosis (MS) (805 patients of Spanish Caucasian origin and 952 health controls). We did not find evidence for association with any single nucleotide polymorphisms (SNPs) tested. Allele and genotype frequencies of the SNPs, which were studied, were similar in DRB1*15-positive or DRB1*15-negative patients. After stratification of MS patients by clinical course, a weak association was observed with rs2069762 G allele and haplotype bearing this allele with secondary progressive MS, although these cases represent 22% of the MS cases. Our results did not show major influence of TENR-IL2-IL21 locus on susceptibility or disease progression in MS. However, we could not exclude completely the effect in MS for this region. Additional studies, using much larger sample sizes and analysis of additional polymorphisms in the gene and its flanking region, will be required to ascertain their contributions to MS susceptibility.

The T244I variant of the interleukin-7 receptor-alpha gene and multiple sclerosis

Several but not all studies have provided evidence for the association between multiple sclerosis (MS) and the T244I variant of the interleukin-7 receptor-alpha gene (IL7RA), rs6897932. We performed a new replication case-control study in 599 MS patients and 594 healthy controls, all Caucasians from the south of Spain. The genotype and allele frequencies differed between MS cases and controls. The IL7RA rs6897932 C allele and the CC genotype were found to be factors for disease susceptibility [per allele odds ratio (OR) 1.32, 95% CI 1.1-1.6, P=0.0031; per CC genotype vs TT + TC genotypes, OR 1.5, 95% CI 1.18-1.87, P=0.0007]. The combined data analysis included 3324 cases and 5032 controls of Europeans and Americans of European origin resulting in stronger association with similar OR (P=1.9 x 10E-9). These findings in our sample support previous reported association studies between IL7RA rs6897932 and MS.

Interferon regulatory factor 5 (IRF5) gene variants are associated with multiple sclerosis in three distinct populations

Background:

IRF5 is a transcription factor involved both in the type I interferon and the toll-like receptor signalling pathways. Previously, IRF5 has been found to be associated with systemic lupus erythematosus, rheumatoid arthritis and inflammatory bowel diseases. Here we investigated whether polymorphisms in the IRF5 gene would be associated with yet another disease with features of autoimmunity, multiple sclerosis (MS).

Methods:

We genotyped nine single nucleotide polymorphisms and one insertion-deletion polymorphism in the IRF5 gene in a collection of 2337 patients with MS and 2813 controls from three populations: two case–control cohorts from Spain and Sweden, and a set of MS trio families from Finland.

Results:

Two single nucleotide polymorphism (SNPs) (rs4728142, rs3807306), and a 5 bp insertion-deletion polymorphism located in the promoter and first intron of the IRF5 gene, showed association signals with values of p<0.001 when the data from all cohorts were combined. The predisposing alleles were present on the same common haplotype in all populations. Using electrophoretic mobility shift assays we observed allele specific differences in protein binding for the SNP rs4728142 and the 5 bp indel, and by a proximity ligation assay we demonstrated increased binding of the transcription factor SP1 to the risk allele of the 5 bp indel.

Conclusion:

These findings add IRF5 to the short list of genes shown to be associated with MS in more than one population. Our study adds to the evidence that there might be genes or pathways that are common in multiple autoimmune diseases, and that the type I interferon system is likely to be involved in the development of these diseases.

OAS1 gene haplotype confers susceptibility to multiple sclerosis

Multiple sclerosis (MS) is associated with genetic susceptibility and unknown environmental triggers, possible viral infections, but the specific etiological mechanism that subsequently develops into an inflammatory/autoimmune cascade of events is poorly understood. Recently, genetic variants of 2',5'- oligoadenylate synthetase 1 (OAS1) gene, a critical enzyme involved in innate antivirus response, have been associated with differential enzyme activity and type 1 diabetes in both case-control and family studies. We hypothesized that polymorphisms in the OAS1 gene could influence the susceptibility to MS. To test this hypothesis, we conducted a case-control study of 333 patients with MS and 424 healthy controls and genotyped two OAS1 single nucleotide polymorphisms (SNPs) by restriction fragment length polymorphism method: rs 10774671, A/G SNP altering the splicing site at the seventh exon, and rs 3741981, a nonsynonymous (Ser162Gly) A/G SNP in the third exon. Haplotype but not single-marker analysis revealed an association of the haplotype created by the G allele at rs 10774671 and the A allele at rs 3741981 with the susceptibility to MS (P value = 8.8 x 10(-5)). Subjects carrying this haplotype had an increased risk of MS comparing with those not carrying it (odds ratio = 4.7, 95% confidence interval 2.1-10.9). Our findings indicate that the OAS1 gene polymorphisms may confer susceptibility to MS or serve as markers of functional variants and suggest that OAS1 activity is involved in the etiology of the disease. Future studies in a larger sample and association analysis with functional variants will clarify the role of the OAS1 gene in the susceptibility to MS.

The 1858T PTPN22 gene variant contributes to a genetic risk of type 1 diabetes in a Ukrainian population

The 1858T variant of the protein tyrosine phosphatase gene, PTPN22, is associated with an increased risk of several autoimmune diseases. The aim of this study has been to investigate the possible association of 1858C-->T PTPN22 polymorphism and type 1 diabetes (T1D) in Caucasians from Ukraine. Overall, the distribution of 1858 PTPN22 genotypes differed significantly between the T1D patient group (n = 296) and the control group (n = 242) (P = 0.0036). When both groups were classified according to sex, the TT genotype and T allele showed a statistically significant higher frequency in T1D female patients (5.9 and 22.8%, respectively) in comparison with the female controls (0 and 11.9%) (P = 0.008 for both analyses). The patients with the TT genotype were significantly younger at the onset of T1D compared with those with genotypes TC and CC (P = 0.035 and 0.019, respectively). In our Ukrainian Caucasian cohort, we confirmed the association between T1D and the PTPN22,1858T allele.

IFNAR1 and IFNAR2 polymorphisms confer susceptibility to multiple sclerosis but not to interferon-beta treatment response

We investigated the role of three polymorphisms in the IFNAR1 (SNPs 18417 and -408) and IFNAR2 (SNP 11876) genes in multiple sclerosis (MS) susceptibility and in the IFNbeta treatment response in a group of 147 patients and 210 controls undergoing interferon therapy during the last 2 years. Only the 18417 and the 11876 SNPs showed an association with disease susceptibility (p=0.001 and 0.035, respectively) although no differential genotype distribution were observed between interferon responders and non-responder MS patients. No alteration of the expression level of IFNAR-1 was observed with respect to the -408 genotypes or to interferon treatment response. These data suggest a role for the IFNAR pathway in susceptibility to MS.