Inflammatory and neurodegenerative serum protein biomarkers increase sensitivity to detect clinical and radiographic disease activity in multiple sclerosis.

The multifaceted nature of multiple sclerosis requires quantitative biomarkers that can provide insights related to diverse physiological pathways. To this end, proteomic analysis of deeply-phenotyped serum samples, biological pathway modeling, and network analysis were performed to elucidate inflammatory and neurodegenerative processes, identifying sensitive biomarkers of multiple sclerosis disease activity. Here, we evaluated the concentrations of > 1400 serum proteins in 630 samples from three multiple sclerosis cohorts for association with clinical and radiographic new disease activity. Twenty proteins were associated with increased clinical and radiographic multiple sclerosis disease activity for inclusion in a custom assay panel. Serum neurofilament light chain showed the strongest univariate correlation with gadolinium lesion activity, clinical relapse status, and annualized relapse rate. Multivariate modeling outperformed univariate for all endpoints. A comprehensive biomarker panel including the twenty proteins identified in this study could serve to characterize disease activity for a patient with multiple sclerosis.

High-resolution DNA methylation screening of the major histocompatibility complex in multiple sclerosis.

Background: The HLA-DRB1 gene in the major histocompatibility complex (MHC) region in chromosome 6p21 is the strongest genetic factor identified as influencing multiple sclerosis (MS) susceptibility. DNA methylation changes associated with MS have been consistently detected at the MHC region. However, understanding the full scope of epigenetic regulations of the MHC remains incomplete, due in part to the limited coverage of this region by standard whole genome bisulfite sequencing or array-based methods. Methods: We developed and validated an MHC capture protocol coupled with bisulfite sequencing and conducted a comprehensive analysis of the MHC methylation landscape in blood samples from 147 treatment naïve MS study participants and 129 healthy controls. Results: We identified 132 differentially methylated region (DMRs) within MHC region associated with disease status. The DMRs overlapped with established MS risk loci. Integration of the MHC methylome with human leukocyte antigen (HLA) genetic data indicate that the methylation changes are significantly associated with HLA genotypes. Using DNA methylation quantitative trait loci (mQTL) mapping and the causal inference test (CIT), we identified 643 cis-mQTL-DMRs paired associations, including 71 DMRs possibly mediating causal relationships between 55 single nucleotide polymorphisms (SNPs) and MS risk. Results: The results describe MS-associated methylation changes in MHC region and highlight the association between HLA genotypes and methylation changes. Results from the mQTL and CIT analyses provide evidence linking MHC region variations, methylation changes, and disease risk for MS.

Inflammatory and neurodegenerative serum protein biomarkers increase sensitivity to detect disease activity in multiple sclerosis.

Background/Objectives: Serum proteomic analysis of deeply-phenotyped samples, biological pathway modeling and network analysis were performed to elucidate the inflammatory and neurodegenerative processes of multiple sclerosis (MS) and identify sensitive biomarkers of MS disease activity (DA). Methods: Over 1100 serum proteins were evaluated in >600 samples from three MS cohorts to identify biomarkers of clinical and radiographic (gadolinium-enhancing lesions) new MS DA. Protein levels were analyzed and associated with presence of gadolinium-enhancing lesions, clinical relapse status (CRS), and annualized relapse rate (ARR) to create a custom assay panel. Results: Twenty proteins were associated with increased clinical and radiographic MS DA. Serum neurofilament light chain (NfL) showed the strongest univariate correlation with radiographic and clinical DA measures. Multivariate modeling significantly outperformed univariate NfL to predict gadolinium lesion activity, CRS and ARR. Discussion: These findings provide insight regarding correlations between inflammatory and neurodegenerative biomarkers and clinical and radiographic MS DA. Funding: Octave Bioscience, Inc (Menlo Park, CA).

Specific hypomethylation programs underpin B cell activation in early multiple sclerosis.

Epigenetic changes have been consistently detected in different cell types in multiple sclerosis (MS). However, their contribution to MS pathogenesis remains poorly understood partly because of sample heterogeneity and limited coverage of array-based methods. To fill this gap, we conducted a comprehensive analysis of genome-wide DNA methylation patterns in four peripheral immune cell populations isolated from 29 MS patients at clinical disease onset and 24 healthy controls. We show that B cells from new-onset untreated MS cases display more significant methylation changes than other disease-implicated immune cell types, consisting of a global DNA hypomethylation signature. Importantly, 4,933 MS-associated differentially methylated regions in B cells were identified, and this epigenetic signature underlies specific genetic programs involved in B cell differentiation and activation. Integration of the methylome to changes in gene expression and susceptibility-associated regions further indicates that hypomethylated regions are significantly associated with the up-regulation of cell activation transcriptional programs. Altogether, these findings implicate aberrant B cell function in MS etiology.

High Resolution Haplotype Analyses of Classical HLA Genes in Families With Multiple Sclerosis Highlights the Role of HLA-DP Alleles in Disease Susceptibility.

Multiple sclerosis (MS) susceptibility shows strong genetic associations with HLA alleles and haplotypes. We genotyped 11 HLA genes in 477 non-Hispanic European MS patients and their 954 unaffected parents using a validated next-generation sequencing (NGS) methodology. HLA haplotypes were assigned unequivocally by tracing HLA allele transmissions. We explored HLA haplotype/allele associations with MS using the genotypic transmission disequilibrium test (gTDT) and multiallelic TDT (mTDT). We also conducted a case-control (CC) study with all patients and 2029 healthy unrelated ethnically matched controls. We performed separate analyses of 54 extended multi-case families by reviewing transmission of haplotype blocks. The haplotype fragment including DRB5*01:01:01~DRB1*15:01:01:01 was significantly associated with predisposition (gTDT: p < 2.20e-16; mTDT: p =1.61e-07; CC: p < 2.22e-16) as reported previously. A second risk allele, DPB1*104:01 (gTDT: p = 3.69e-03; mTDT: p = 2.99e-03; CC: p = 1.00e-02), independent from the haplotype bearing DRB1*15:01 was newly identified. The allele DRB1*01:01:01 showed significant protection (gTDT: p = 8.68e-06; mTDT: p = 4.50e-03; CC: p = 1.96e-06). Two DQB1 alleles, DQB1*03:01 (gTDT: p = 2.86e-03; mTDT: p = 5.56e-02; CC: p = 4.08e-05) and DQB1*03:03 (gTDT: p = 1.17e-02; mTDT: p = 1.16e-02; CC: p = 1.21e-02), defined at two-field level also showed protective effects. The HLA class I block, A*02:01:01:01~C*03:04:01:01~B*40:01:02 (gTDT: p = 5.86e-03; mTDT: p = 3.65e-02; CC: p = 9.69e-03) and the alleles B*27:05 (gTDT: p = 6.28e-04; mTDT: p = 2.15e-03; CC: p = 1.47e-02) and B*38:01 (gTDT: p = 3.20e-03; mTDT: p = 6.14e-03; CC: p = 1.70e-02) showed moderately protective effects independently from each other and from the class II associated factors. By comparing statistical significance of 11 HLA loci and 19 haplotype segments with both untruncated and two-field allele names, we precisely mapped MS candidate alleles/haplotypes while eliminating false signals resulting from 'hitchhiking' alleles. We assessed genetic burden for the HLA allele/haplotype identified in this study. This family-based study including the highest-resolution of HLA alleles proved to be powerful and efficient for precise identification of HLA genotypes associated with both, susceptibility and protection to development of MS.

Next-Generation Sequencing Identifies Extended HLA Class I and II Haplotypes Associated With Early-Onset and Late-Onset Myasthenia Gravis in Italian, Norwegian, and Swedish Populations.

Genetic susceptibility to myasthenia gravis (MG) associates with specific HLA alleles and haplotypes at the class I and II regions in various populations. Previous studies have only examined alleles at a limited number of HLA loci that defined only broad serotypes or alleles defined at the protein sequence level. Consequently, genetic variants in noncoding and untranslated HLA gene segments have not been fully explored but could also be important determinants for MG. To gain further insight into the role of HLA in MG, we applied next-generation sequencing to analyze sequence variation at eleven HLA genes in early-onset (EO) and late-onset (LO) non-thymomatous MG patients positive for the acetylcholine receptor (AChR) antibodies and ethnically matched controls from Italy, Norway, and Sweden. For all three populations, alleles and haplotype blocks present on the ancestral haplotype AH8.1 were associated with risk in AChR-EOMG patients. HLA-B*08:01:01:01 was the dominant risk allele in Italians (OR = 3.28, P = 1.83E-05), Norwegians (OR = 3.52, P = 4.41E-16), and in Swedes HLA-B*08:01 was the primary risk allele (OR = 4.24, P <2.2E-16). Protective alleles and haplotype blocks were identified on the HLA-DRB7, and HLA-DRB13.1 class II haplotypes in Italians and Norwegians, whereas in Swedes HLA-DRB7 exhibited the main protective effect. For AChR-LOMG patients, the HLA-DRB15.1 haplotype and associated alleles were significantly associated with susceptibility in all groups. The HLA-DR13-HLA-DR-HLA-DQ haplotype was associated with protection in all AChR-LOMG groups. This study has confirmed and extended previous findings that the immunogenetic predisposition profiles for EOMG and LOMG are distinct. In addition, the results are consistent with a role for non-coding HLA genetic variants in the pathogenesis of MG.

A splice acceptor variant in HLA-DRA affects the conformation and cellular localization of the class II DR alpha-chain.

Class II human leucocyte antigen (HLA) proteins are involved in the immune response by presenting pathogen-derived peptides to CD4+ T lymphocytes. At the molecular level, they are constituted by α/ß-heterodimers on the surface of professional antigen-presenting cells. Here, we report that the acceptor variant (rs8084) in the HLA-DRA gene mediates the transcription of an alternative version of the α-chain lacking 25 amino acids in its extracellular domain. Molecular dynamics simulations suggest this isoform undergoes structural refolding which in turn affects its stability and cellular trafficking. The short HLA-DRA isoform cannot reach the cell surface, although it is still able to bind the corresponding ß-chain. Conversely, it remains entrapped within the endoplasmic reticulum where it is targeted for degradation. Furthermore, we demonstrate that the short isoform can be transported to the cell membrane via interactions with the peptide-binding site of canonical HLA heterodimers. Altogether, our findings indicate that short HLA-DRA functions as a novel intact antigen for class II HLA molecules.

Cell type-specific transcriptomics identifies neddylation as a novel therapeutic target in multiple sclerosis.

Multiple sclerosis is an autoimmune disease of the CNS in which both genetic and environmental factors are involved. Genome-wide association studies revealed more than 200 risk loci, most of which harbour genes primarily expressed in immune cells. However, whether genetic differences are translated into cell-specific gene expression profiles and to what extent these are altered in patients with multiple sclerosis are still open questions in the field. To assess cell type-specific gene expression in a large cohort of patients with multiple sclerosis, we sequenced the whole transcriptome of fluorescence-activated cell sorted T cells (CD4+ and CD8+) and CD14+ monocytes from treatment-naive patients with multiple sclerosis (n = 106) and healthy subjects (n = 22). We identified 479 differentially expressed genes in CD4+ T cells, 435 in monocytes, and 54 in CD8+ T cells. Importantly, in CD4+ T cells, we discovered upregulated transcripts from the NAE1 gene, a critical subunit of the NEDD8 activating enzyme, which activates the neddylation pathway, a post-translational modification analogous to ubiquitination. Finally, we demonstrated that inhibition of NEDD8 activating enzyme using the specific inhibitor pevonedistat (MLN4924) significantly ameliorated disease severity in murine experimental autoimmune encephalomyelitis. Our findings provide novel insights into multiple sclerosis-associated gene regulation unravelling neddylation as a crucial pathway in multiple sclerosis pathogenesis with implications for the development of tailored disease-modifying agents.

Ataxin-1 regulates B cell function and the severity of autoimmune experimental encephalomyelitis.

Ataxin-1 (ATXN1) is a ubiquitous polyglutamine protein expressed primarily in the nucleus where it binds chromatin and functions as a transcriptional repressor. Mutant forms of ataxin-1 containing expanded glutamine stretches cause the movement disorder spinocerebellar ataxia type 1 (SCA1) through a toxic gain-of-function mechanism in the cerebellum. Conversely, ATXN1 loss-of-function is implicated in cancer development and Alzheimer's disease (AD) pathogenesis. ATXN1 was recently nominated as a susceptibility locus for multiple sclerosis (MS). Here, we show that Atxn1-null mice develop a more severe experimental autoimmune encephalomyelitis (EAE) course compared to wildtype mice. The aggravated phenotype is mediated by increased T helper type 1 (Th1) cell polarization, which in turn results from the dysregulation of B cell activity. Ataxin-1 ablation in B cells leads to aberrant expression of key costimulatory molecules involved in proinflammatory T cell differentiation, including cluster of differentiation (CD)44 and CD80. In addition, comprehensive phosphoflow cytometry and transcriptional profiling link the exaggerated proliferation of ataxin-1 deficient B cells to the activation of extracellular signal-regulated kinase (ERK) and signal transducer and activator of transcription (STAT) pathways. Lastly, selective deletion of the physiological binding partner capicua (CIC) demonstrates the importance of ATXN1 native interactions for correct B cell functioning. Altogether, we report a immunomodulatory role for ataxin-1 and provide a functional description of the ATXN1 locus genetic association with MS risk.

Association Between Serum Neurofilament Light Chain Levels and Long-term Disease Course Among Patients With Multiple Sclerosis Followed up for 12 Years.

IMPORTANCE: Blood sample-based biomarkers that are associated with clinically meaningful outcomes for patients with multiple sclerosis (MS) have not been developed. OBJECTIVE: To evaluate the potential of serum neurofilament light chain (sNFL) measurements as a biomarker of disease activity and progression in a longitudinal MS data set. DESIGN, SETTING, AND PARTICIPANTS: Single-center, ongoing, prospective observational cohort study of 607 patients with MS from the longitudinal EPIC (Expression, Proteomics, Imaging, Clinical) study at the University of California, San Francisco from July 1, 2004, through August 31, 2017. Clinical evaluations and sample collection were performed annually for 5 years, then at different time points for up to 12 years, with a median follow-up duration of 10 (interquartile range, 7-11) years. Serum NFL levels were measured using a sensitive single molecule array platform and compared with clinical and magnetic resonance imaging variables with the use of univariable and multivariable analyses. MAIN OUTCOMES AND MEASURES: The main outcomes were disability progression defined as clinically significant worsening on the Expanded Disability Status Scale (EDSS) score and brain fraction atrophy. RESULTS: Mean (SD) age of the 607 study participants at study entry was 42.5 (9.8) years; 423 (69.7%) were women; and all participants were of non-Hispanic European descent. Of 3911 samples sequentially collected, 3904 passed quality control for quantification of sNFL. Baseline sNFL levels showed significant associations with EDSS score (ß, 1.080; 95% CI, 1.047-1.114; P < .001), MS subtype (ß, 1.478; 95% CI, 1.279-1.707; P < .001), and treatment status (ß, 1.120; 95% CI, 1.007-1.245; P = .04). A significant interaction between EDSS worsening and change in levels of sNFL over time was found (ß, 1.015; 95% CI, 1.007-1.023; P < .001). Baseline sNFL levels alone were associated with approximately 11.6% of the variance in brain fraction atrophy at year 10. In a multivariable analysis that considered sex, age, and disease duration, baseline sNFL levels were associated with 18.0% of the variance in brain fraction atrophy at year 10. After 5 years' follow-up, active treatment was associated with lower levels of sNFL, with high-potency treatments associated with the greater decreases in sNFL levels compared with platform therapies (high-potency vs untreated: ß, 0.946; 95% CI, 0.915-0.976; P < .001; high-potency vs platform: ß, 0.972; 95% CI, 0.948-0.998; P = .04). CONCLUSIONS AND RELEVANCE: This study found that statistically significant associations of sNFL with relevant clinical and neuroimaging outcomes in MS were confirmed and extended, supporting the potential of sNFL as an objective surrogate of ongoing MS disease activity. In this data set of patients with MS who received early treatment, the prognostic power of sNFL for relapse activity and long-term disability progression was limited. Further prospective studies are necessary to assess the assay's utility for decision-making in individual patients.

Next-generation sequencing reveals new information about HLA allele and haplotype diversity in a large European American population.

The human leukocyte antigen (HLA) genes are extremely polymorphic and are useful molecular markers to make inferences about human population history. However, the accuracy of the estimation of genetic diversity at HLA loci very much depends on the technology used to characterize HLA alleles; high-resolution genotyping of long-range HLA gene products improves the assessment of HLA population diversity as well as other population parameters compared to lower resolution typing methods. In this study we examined allelic and haplotype HLA diversity in a large healthy European American population sourced from the UCSF-DNA bank. A high-resolution next-generation sequencing method was applied to define non-ambiguous 3- and 4-field alleles at the HLA-A, HLA-C, HLA-B, HLA-DRB1, HLA-DRB3/4/5, HLA-DQA1, HLA-DQB1, HLA-DPA1, and HLA-DPB1 loci in samples provided by 2248 unrelated individuals. A number of population parameters were examined including balancing selection and various measurements of linkage disequilibrium were calculated. There were no detectable deviations from Hardy-Weinberg proportions at HLA-A, HLA-DRB1, HLA-DQA1 and HLA-DQB1. For the remaining loci moderate and significant deviations were detected at HLA-C, HLA-B, HLA-DRB3/4/5, HLA-DPA1 and HLA-DPB1 loci mostly from population substructures. Unique 4-field associations were observed among alleles at 2 loci and haplotypes extending large intervals that were not apparent in results obtained using testing methodologies with limited sequence coverage and phasing. The high diversity at HLA-DPA1 results from detection of intron variants of otherwise well conserved protein sequences. It may be speculated that divergence in exon sequences may be negatively selected. Our data provides a valuable reference source for future population studies that may allow for precise fine mapping of coding and non-coding sequences determining disease susceptibility and allo-immunogenicity.

Sex-specific Tau methylation patterns and synaptic transcriptional alterations are associated with neural vulnerability during chronic neuroinflammation.

The molecular events underlying the transition from initial inflammatory flares to the progressive phase of multiple sclerosis (MS) remain poorly understood. Here, we report that the microtubule-associated protein (MAP) Tau exerts a gender-specific protective function on disease progression in the MS model experimental autoimmune encephalomyelitis (EAE). A detailed investigation of the autoimmune response in Tau-deficient mice excluded a strong immunoregulatory role for Tau, suggesting that its beneficial effects are presumably exerted within the central nervous system (CNS). Spinal cord transcriptomic data show increased synaptic dysfunctions and alterations in the NF-kB activation pathway upon EAE in Tau-deficient mice as compared to wildtype animals. We also performed the first comprehensive characterization of Tau post-translational modifications (PTMs) in the nervous system upon EAE. We report that the methylation levels of the conserved lysine residue K306 are significantly decreased in the chronic phase of the disease. By combining biochemical assays and molecular dynamics (MD) simulations, we demonstrate that methylation at K306 decreases the affinity of Tau for the microtubule network. Thus, the down-regulation of this PTM might represent a homeostatic response to enhance axonal stability against an autoimmune CNS insult. The results, altogether, position Tau as key mediator between the inflammatory processes and neurodegeneration that seems to unify many CNS diseases.

Deconstruction of HLA-DRB1*04:01:01 and HLA-DRB1*15:01:01 class II haplotypes using next-generation sequencing in European-Americans with multiple sclerosis.

BACKGROUND: The association between HLA-DRB1*15:01 with multiple sclerosis (MS) susceptibility is well established, but the contribution of the tightly associated HLA-DRB5*01:01 allele has not yet been completely ascertained. Similarly, the effects of HLA-DRB1*04:01 alleles and haplotypes, defined at the full-gene resolution level with MS risk remains to be elucidated. OBJECTIVES: To characterize the molecular architecture of class II HLA-DR15 and HLA-DR4 haplotypes associated with MS. METHODS: Next-generation sequencing was used to determine HLA-DQB1, HLA-DQA1, and HLA-DRB1/4/5 alleles in 1403 unrelated European-American patients and 1425 healthy unrelated controls. Effect sizes of HLA alleles and haplotypes on MS risk were measured by odds ratio (OR) with 95% confidence intervals. RESULTS: HLA-DRB1*15:01:01:01SG (OR = 3.20, p < 2.2E-16), HLA-DRB5*01:01:01 (OR = 2.96, p < 2.2E-16), and HLA-DRB5*01:01:01v1_STR1 (OR = 8.18, p = 4.3E-05) alleles all occurred at significantly higher frequencies in MS patients compared to controls. The most significant predis-posing haplotypes were HLA-DQB1*06:02:01~ HLA-DQA1*01:02:01:01SG~HLA-DRB1*15:01:01:01SG~HLA-DRB5*01:01:01 and HLA-DQB1*06:02:01~HLA-DQA1*01:02:01:01SG~HLA-DRB1*15:01:01:01SG~HLA-DRB5*01:01:01v1_STR1 (OR = 3.19, p < 2.2E-16; OR = 9.30, p = 9.7E-05, respectively). Analyses of the HLA-DRB1*04 cohort in the absence of HLA-DRB1*15:01 haplotypes revealed that the HLA-DQB1*03:01:01:01~HLA-DQA1*03:03:01:01~HLA-DRB1*04:01:01:01SG~HLA-DRB4*01:03:01:01 haplotype was protective (OR = 0.64, p = 0.028), whereas the HLA-DQB1*03:02:01~HLA-DQA1*03:01:01~HLA-DRB1*04:01:01:01SG~HLA-DRB4*01:03:01:01 haplotype was associated with MS susceptibility (OR = 1.66, p = 4.9E-03). CONCLUSION: HLA-DR15 haplotypes, including genomic variants of HLA-DRB5, and HLA-DR4 haplotypes affect MS risk.

MIF and D-DT are potential disease severity modifiers in male MS subjects.

Little is known about mechanisms that drive the development of progressive multiple sclerosis (MS), although inflammatory factors, such as macrophage migration inhibitory factor (MIF), its homolog D-dopachrome tautomerase (D-DT), and their common receptor CD74 may contribute to disease worsening. Our findings demonstrate elevated MIF and D-DT levels in males with progressive disease compared with relapsing-remitting males (RRMS) and female MS subjects, with increased levels of CD74 in females vs. males with high MS disease severity. Furthermore, increased MIF and D-DT levels in males with progressive disease were significantly correlated with the presence of two high-expression promoter polymorphisms located in the MIF gene, a -794CATT5-8 microsatellite repeat and a -173 G/C SNP. Conversely, mice lacking MIF or D-DT developed less-severe signs of experimental autoimmune encephalomyelitis, a murine model of MS, thus implicating both homologs as copathogenic contributors. These findings indicate that genetically controlled high MIF expression (and D-DT) promotes MS progression in males, suggesting that these two factors are sex-specific disease modifiers and raising the possibility that aggressive anti-MIF treatment of clinically isolated syndrome or RRMS males with a high-expresser genotype might slow or prevent the onset of progressive MS. Additionally, selective targeting of MIF:CD74 signaling might provide an effective, trackable therapeutic approach for MS subjects of both sexes.

Multiple sclerosis risk loci and disease severity in 7,125 individuals from 10 studies.

OBJECTIVE: We investigated the association between 52 risk variants identified through genome-wide association studies and disease severity in multiple sclerosis (MS). METHODS: Ten unique MS case data sets were analyzed. The Multiple Sclerosis Severity Score (MSSS) was calculated using the Expanded Disability Status Scale at study entry and disease duration. MSSS was considered as a continuous variable and as 2 dichotomous variables (median and extreme ends; MSSS of ≤5 vs >5 and MSSS of <2.5 vs ≥7.5, respectively). Single nucleotide polymorphisms (SNPs) were examined individually and as both combined weighted genetic risk score (wGRS) and unweighted genetic risk score (GRS) for association with disease severity. Random-effects meta-analyses were conducted and adjusted for cohort, sex, age at onset, and HLA-DRB1*15:01. RESULTS: A total of 7,125 MS cases were analyzed. The wGRS and GRS were not strongly associated with disease severity after accounting for cohort, sex, age at onset, and HLA-DRB1*15:01. After restricting analyses to cases with disease duration ≥10 years, associations were null (p value ≥0.05). No SNP was associated with disease severity after adjusting for multiple testing. CONCLUSIONS: The largest meta-analysis of established MS genetic risk variants and disease severity, to date, was performed. Results suggest that the investigated MS genetic risk variants are not associated with MSSS, even after controlling for potential confounders. Further research in large cohorts is needed to identify genetic determinants of disease severity using sensitive clinical and MRI measures, which are critical to understanding disease mechanisms and guiding development of effective treatments.

Association of HLA Genetic Risk Burden With Disease Phenotypes in Multiple Sclerosis.

IMPORTANCE: Although multiple HLA alleles associated with multiple sclerosis (MS) risk have been identified, genotype-phenotype studies in the HLA region remain scarce and inconclusive. OBJECTIVES: To investigate whether MS risk-associated HLA alleles also affect disease phenotypes. DESIGN, SETTING, AND PARTICIPANTS: A cross-sectional, case-control study comprising 652 patients with MS who had comprehensive phenotypic information and 455 individuals of European origin serving as controls was conducted at a single academic research site. Patients evaluated at the Multiple Sclerosis Center at University of California, San Francisco between July 2004 and September 2005 were invited to participate. Spinal cord imaging in the data set was acquired between July 2013 and March 2014; analysis was performed between December 2014 and December 2015. MAIN OUTCOMES AND MEASURES: Cumulative HLA genetic burden (HLAGB) calculated using the most updated MS-associated HLA alleles vs clinical and magnetic resonance imaging outcomes, including age at onset, disease severity, conversion time from clinically isolated syndrome to clinically definite MS, fractions of cortical and subcortical gray matter and cerebral white matter, brain lesion volume, spinal cord gray and white matter areas, upper cervical cord area, and the ratio of gray matter to the upper cervical cord area. Multivariate modeling was applied separately for each sex data set. RESULTS: Of the 652 patients with MS, 586 had no missing genetic data and were included in the HLAGB analysis. In these 586 patients (404 women [68.9%]; mean [SD] age at disease onset, 33.6 [9.4] years), HLAGB was higher than in controls (median [IQR], 0.7 [0-1.4] and 0 [-0.3 to 0.5], respectively; P = 1.8 × 10-27). A total of 619 (95.8%) had relapsing-onset MS and 27 (4.2%) had progressive-onset MS. No significant difference was observed between relapsing-onset MS and primary progressive MS. A higher HLAGB was associated with younger age at onset and the atrophy of subcortical gray matter fraction in women with relapsing-onset MS (standard ß = -1.20 × 10-1; P = 1.7 × 10-2 and standard ß = -1.67 × 10-1; P = 2.3 × 10-4, respectively), which were driven mainly by the HLA-DRB1*15:01 haplotype. In addition, we observed the distinct role of the HLA-A*24:02-B*07:02-DRB1*15:01 haplotype among the other common DRB1*15:01 haplotypes and a nominally protective effect of HLA-B*44:02 to the subcortical gray atrophy (standard ß = -1.28 × 10-1; P = 5.1 × 10-3 and standard ß = 9.52 × 10-2; P = 3.6 × 10-2, respectively). CONCLUSIONS AND RELEVANCE: We confirm and extend previous observations linking HLA MS susceptibility alleles with disease progression and specific clinical and magnetic resonance imaging phenotypic traits.

A non-synonymous single-nucleotide polymorphism associated with multiple sclerosis risk affects the EVI5 interactome.

Despite recent progress in the characterization of genetic loci associated with multiple sclerosis (MS) risk, the ubiquitous linkage disequilibrium operating across the genome has stalled efforts to distinguish causative variants from proxy single-nucleotide polymorphisms (SNPs). Here, we have identified through fine mapping and meta-analysis EVI5 as the most plausible disease risk gene within the 1p22.1 locus. We further show that an exonic SNP associated with risk induces changes in superficial hydrophobicity patterns of the coiled-coil domain of EVI5, which, in turns, affects the EVI5 interactome. Immunoprecipitation of wild-type and mutated EVI5 followed by mass spectrometry generated a roster of disease-specific interactors functionally linked to lipid metabolism. Among the exclusive binding partners of the risk variant, we describe the novel interaction with sphingosine 1-phosphate lyase (SGPL1)-a key enzyme for the creation of the sphingosine-1 phosphate gradient, which is relevant to the pathogenic process and therapeutic management of MS.

Mitochondrial DNA sequence variation in multiple sclerosis.

OBJECTIVE: To assess the influence of common mitochondrial DNA (mtDNA) sequence variation on multiple sclerosis (MS) risk in cases and controls part of an international consortium. METHODS: We analyzed 115 high-quality mtDNA variants and common haplogroups from a previously published genome-wide association study among 7,391 cases from the International Multiple Sclerosis Genetics Consortium and 14,568 controls from the Wellcome Trust Case Control Consortium 2 project from 7 countries. Significant single nucleotide polymorphism and haplogroup associations were replicated in 3,720 cases and 879 controls from the University of California, San Francisco. RESULTS: An elevated risk of MS was detected among haplogroup JT carriers from 7 pooled clinic sites (odds ratio [OR] = 1.15, 95% confidence interval [CI] = 1.07-1.24, p = 0.0002) included in the discovery study. The increased risk of MS was observed for both haplogroup T (OR = 1.17, 95% CI = 1.06-1.29, p = 0.002) and haplogroup J carriers (OR = 1.11, 95% CI = 1.01-1.22, p = 0.03). These haplogroup associations with MS were not replicated in the independent sample set. An elevated risk of primary progressive (PP) MS was detected for haplogroup J participants from 3 European discovery populations (OR = 1.49, 95% CI = 1.10-2.01, p = 0.009). This elevated risk was borderline significant in the US replication population (OR = 1.43, 95% CI = 0.99-2.08, p = 0.058) and remained significant in pooled analysis of discovery and replication studies (OR = 1.43, 95% CI = 1.14-1.81, p = 0.002). No common individual mtDNA variants were associated with MS risk. CONCLUSIONS: Identification and validation of mitochondrial genetic variants associated with MS and PPMS may lead to new targets for treatment and diagnostic tests for identifying potential responders to interventions that target mitochondria.

An ImmunoChip study of multiple sclerosis risk in African Americans.

The aims of this study were: (i) to determine to what degree multiple sclerosis-associated loci discovered in European populations also influence susceptibility in African Americans; (ii) to assess the extent to which the unique linkage disequilibrium patterns in African Americans can contribute to localizing the functionally relevant regions or genes; and (iii) to search for novel African American multiple sclerosis-associated loci. Using the ImmunoChip custom array we genotyped 803 African American cases with multiple sclerosis and 1516 African American control subjects at 130 135 autosomal single nucleotide polymorphisms. We conducted association analysis with rigorous adjustments for population stratification and admixture. Of the 110 non-major histocompatibility complex multiple sclerosis-associated variants identified in Europeans, 96 passed stringent quality control in our African American data set and of these, >70% (69) showed over-representation of the same allele amongst cases, including 21 with nominally significant evidence for association (one-tailed test P < 0.05). At a further eight loci we found nominally significant association with an alternate correlated risk-tagging single nucleotide polymorphism from the same region. Outside the regions known to be associated in Europeans, we found seven potentially associated novel candidate multiple sclerosis variants (P < 10(-4)), one of which (rs2702180) also showed nominally significant evidence for association (one-tailed test P = 0.034) in an independent second cohort of 620 African American cases and 1565 control subjects. However, none of these novel associations reached genome-wide significance (combined P = 6.3 × 10(-5)). Our data demonstrate substantial overlap between African American and European multiple sclerosis variants, indicating common genetic contributions to multiple sclerosis risk.

SNP imputation bias reduces effect size determination.

Imputation is a commonly used technique that exploits linkage disequilibrium to infer missing genotypes in genetic datasets, using a well-characterized reference population. While there is agreement that the reference population has to match the ethnicity of the query dataset, it is common practice to use the same reference to impute genotypes for a wide variety of phenotypes. We hypothesized that using a reference composed of samples with a different phenotype than the query dataset would introduce imputation bias. To test this hypothesis we used GWAS datasets from Amyotrophic Lateral Sclerosis (ALS), Parkinson Disease (PD), and Crohn's Disease (CD). First, we masked and then performed imputation of 100 disease-associated markers and 100 non-associated markers from each study. Two references for imputation were used in parallel: one consisting of healthy controls and another consisting of patients with the same disease. We assessed the discordance (imprecision) and bias (inaccuracy) of imputation by comparing predicted genotypes to those assayed by SNP-chip. We also assessed the bias on the observed effect size when the predicted genotypes were used in a GWAS study. When healthy controls were used as reference for imputation, a significant bias was observed, particularly in the disease-associated markers. Using cases as reference significantly attenuated this bias. For nearly all markers, the direction of the bias favored the non-risk allele. In GWAS studies of the three diseases (with healthy reference controls from the 1000 genomes as reference), the mean OR for disease-associated markers obtained by imputation was lower than that obtained using original assayed genotypes. We found that the bias is inherent to imputation as using different methods did not alter the results. In conclusion, imputation is a powerful method to predict genotypes and estimate genetic risk for GWAS. However, a careful choice of reference population is needed to minimize biases inherent to this approach.