Selected Publications


GWAS meta-analysis of psoriasis identifies new susceptibility alleles impacting disease mechanisms and therapeutic targets
abstract
Psoriasis is a common, debilitating immune-mediated skin disease. Genetic studies have identified biological mechanisms of psoriasis risk, including those targeted by effective therapies. However, the genetic liability to psoriasis is not fully explained by variation at robustly identified risk loci. To move towards a saturation map of psoriasis susceptibility we meta-analysed 18 GWAS comprising 36,466 cases and 458,078 controls and identified 109 distinct psoriasis susceptibility loci, including 45 that have not been previously reported. These include susceptibility variants at loci in which the therapeutic targets IL17RA and AHR are encoded, and deleterious coding variants supporting potential new drug targets (including in STAP2, CPVL and POU2F3). We conducted a transcriptome-wide association study to identify regulatory effects of psoriasis susceptibility variants and cross-referenced these against single cell expression profiles in psoriasis-affected skin, highlighting roles for the transcriptional regulation of haematopoietic cell development and epigenetic modulation of interferon signalling in psoriasis pathobiology.
by: Nick Dand, Philip E Stuart, John Bowes, David Ellinghaus, Joanne Nititham, Jake R Saklatvala, Maris Teder-Laving, Laurent F Thomas, Tanel Traks, Steffen Uebe, Gunter Assmann, David Baudry, Frank Behrens, Allison C Billi, Matthew A Brown, Harald Burkhardt, Francesca Capon, Raymond Chung, Charles J Curtis, Michael Duckworth, Eva Ellinghaus, Oliver FitzGerald, Sascha Gerdes, Christopher E M Griffiths, Susanne Gulliver, Philip Helliwell, Pauline Ho, Per Hoffmann, Oddgeir L Holmen, Zhi-Ming Huang, Kristian Hveem, Deepak Jadon, Michaela Köhm, Cornelia Kraus, Céline Lamacchia, Sang Hyuck Lee, Feiyang Ma, Satveer K Mahil, Neil McHugh, Ross McManus, Ellen H Modalsli, Michael J Nissen, Markus Nöthen, Vinzenz Oji, Jorge R Oksenberg, Matthew T Patrick, Bethany E Perez-White, Andreas Ramming
medRxiv [Preprint]. Oct 5, 2023
PMID: 37873414


MHConstructor: A high-throughput, haplotype-informed solution to the MHC assembly challenge
abstract
The extremely high levels of genetic polymorphism within the human major histocompatibility complex (MHC) limit the usefulness of reference-based alignment methods for sequence assembly. We incorporate a short read de novo assembly algorithm into a workflow for novel application to the MHC. MHConstructor is a containerized pipeline designed for high-throughput, haplotype-informed, reproducible assembly of both whole genome sequencing and target-capture short read data in large, population cohorts. To-date, no other self-contained tool exists for the generation of de novo MHC assemblies from short read data. MHConstructor facilitates wide-spread access to high quality, alignment-free MHC sequence analysis. Keywords: Major histocompatibility complex; de novo assembly; haplotype; human leukocyte antigen genes; short read sequencing.
by: Wade KJ, Suseno R, Kizer K, Williams J, Boquett J, Caillier S, Pollock NR, Renschen A, Santaniello A, Oksenberg JR, Norman PJ, Augusto DG, Hollenbach JA.
bioRxiv [Preprint]. May 21, 2024
PMID: 38826378


A Predictive Autoantibody Signature in Multiple Sclerosis
abstract
Although B cells are implicated in multiple sclerosis (MS) pathophysiology, a predictive or diagnostic autoantibody remains elusive. Here, the Department of Defense Serum Repository (DoDSR), a cohort of over 10 million individuals, was used to generate whole-proteome autoantibody profiles of hundreds of patients with MS (PwMS) years before and subsequently after MS onset. This analysis defines a unique cluster of PwMS that share an autoantibody signature against a common motif that has similarity with many human pathogens. These patients exhibit antibody reactivity years before developing MS symptoms and have higher levels of serum neurofilament light (sNfL) compared to other PwMS. Furthermore, this profile is preserved over time, providing molecular evidence for an immunologically active prodromal period years before clinical onset. This autoantibody reactivity was validated in samples from a separate incident MS cohort in both cerebrospinal fluid (CSF) and serum, where it is highly specific for patients eventually diagnosed with MS. This signature is a starting point for further immunological characterization of this MS patient subset and may be clinically useful as an antigen-specific biomarker for high-risk patients with clinically- or radiologically-isolated neuroinflammatory syndromes. Keywords: Multiple sclerosis; PhIP-Seq; antigen; neurofilament light chain; phage display; serology.
by: Colin R Zamecnik, Gavin M Sowa, Ahmed Abdelhak, Ravi Dandekar, Rebecca D Bair, Kristen J Wade, Christopher M Bartley, Asritha Tubati, Refujia Gomez, Camille Fouassier, Chloe Gerungan, Jessica Alexander, Anne E Wapniarski, Rita P Loudermilk, Erica L Eggers, Kelsey C Zorn, Kirtana Ananth, Nora Jabassini, Sabrina A Mann, Nicholas R Ragan, Adam Santaniello, Roland G Henry, Sergio E Baranzini, Scott S Zamvil, Riley M Bove, Chu-Yueh Guo, Jeffrey M Gelfand, Richard Cuneo, H-Christian von Büdingen, Jorge R Oksenberg, Bruce Ac Cree, Jill A Hollenbach, Ari J Green, Stephen L Hauser, Mitchell T Wallin, Joseph L DeRisi, Michael R Wilson
medRxiv [Preprint]. May 15, 2023
PMID: 37205595


Integration of epigenetic and genetic profiles identifies multiple sclerosis disease-critical cell types and genes
abstract
Genome-wide association studies (GWAS) successfully identified multiple sclerosis (MS) susceptibility variants. Despite this notable progress, understanding the biological context of these associations remains challenging, due in part to the complexity of linking GWAS results to causative genes and cell types. Here, we aimed to address this gap by integrating GWAS data with single-cell and bulk chromatin accessibility data and histone modification profiles from immune and nervous systems. MS-GWAS associations are significantly enriched in regulatory regions of microglia and peripheral immune cell subtypes, especially B cells and monocytes. Cell-specific polygenic risk scores were developed to examine the cumulative impact of the susceptibility genes on MS risk and clinical phenotypes, showing significant associations with risk and brain white matter volume. The findings reveal enrichment of GWAS signals in B cell and monocyte/microglial cell-types, consistent with the known pathology and presumed targets of effective MS therapeutics.
by: Qin Ma, Hengameh Shams, Alessandro Didonna, Sergio E Baranzini, Bruce A C Cree, Stephen L Hauser, Roland G Henry, Jorge R Oksenberg
Commun Biol. Mar 30, 2023
PMID: 36997638


Predicting disease severity in multiple sclerosis using multimodal data and machine learning
abstract
Background: Multiple sclerosis patients would benefit from machine learning algorithms that integrates clinical, imaging and multimodal biomarkers to define the risk of disease activity. Methods: We have analysed a prospective multi-centric cohort of 322 MS patients and 98 healthy controls from four MS centres, collecting disability scales at baseline and 2 years later. Imaging data included brain MRI and optical coherence tomography, and omics included genotyping, cytomics and phosphoproteomic data from peripheral blood mononuclear cells. Predictors of clinical outcomes were searched using Random Forest algorithms. Assessment of the algorithm performance was conducted in an independent prospective cohort of 271 MS patients from a single centre. Results: We found algorithms for predicting confirmed disability accumulation for the different scales, no evidence of disease activity (NEDA), onset of immunotherapy and the escalation from low- to high-efficacy therapy with intermediate to high-accuracy. This accuracy was achieved for most of the predictors using clinical data alone or in combination with imaging data. Still, in some cases, the addition of omics data slightly increased algorithm performance. Accuracies were comparable in both cohorts. Conclusion: Combining clinical, imaging and omics data with machine learning helps identify MS patients at risk of disability worsening. Keywords: Imaging; Machine learning; Multiple sclerosis; Omics; Precision medicine.
by: Magi Andorra, Ana Freire, Irati Zubizarreta, Nicole Kerlero de Rosbo, Steffan D Bos, Melanie Rinas, Einar A Høgestøl, Sigrid A de Rodez Benavent, Tone Berge, Synne Brune-Ingebretse, Federico Ivaldi, Maria Cellerino, Matteo Pardini, Gemma Vila, Irene Pulido-Valdeolivas, Elena H Martinez-Lapiscina, Sara Llufriu, Albert Saiz, Yolanda Blanco, Eloy Martinez-Heras, Elisabeth Solana, Priscilla Bäcker-Koduah, Janina Behrens, Joseph Kuchling, Susanna Asseyer, Michael Scheel, Claudia Chien, Hanna Zimmermann, Seyedamirhosein Motamedi, Josef Kauer-Bonin, Alex Brandt, Julio Saez-Rodriguez, Leonidas G Alexopoulos, Friedemann Paul, Hanne F Harbo, Hengameh Shams, Jorge Oksenberg, Antonio Uccelli, Ricardo Baeza-Yates, Pablo Villoslada
J Neurol. Mar 2024
PMID: 38133801


Multiple sclerosis: two decades of progress
abstract
None provided
by: Bruce A C Cree, Jorge R Oksenberg, Stephen L Hauser
Lancet Neurol. Mar 2022
PMID: 35182500


Plasma neurofilament light chain levels suggest neuroaxonal stability following therapeutic remyelination in people with multiple sclerosis
abstract
Background: Chronic demyelination is a major contributor to axonal vulnerability in multiple sclerosis (MS). Therefore, remyelination could provide a potent neuroprotective strategy. The ReBUILD trial was the first study showing evidence for successful remyelination following treatment with clemastine in people with MS (pwMS) with no evidence of disease activity or progression (NEDAP). Whether remyelination was associated with neuroprotection remains unexplored. Methods: Plasma neurofilament light chain (NfL) levels were measured from ReBUILD trial's participants. Mixed linear effect models were fit for individual patients, epoch and longitudinal measurements to compare NfL concentrations between samples collected during the active and placebo treatment period. Results: NfL concentrations were 9.6% lower in samples collected during the active treatment with clemastine (n=53, geometric mean=6.33 pg/mL) compared to samples collected during treatment with placebo (n=73, 7.00 pg/mL) (B=-0.035 [-0.068 to -0.001], p=0.041). Applying age- and body mass index-standardised NfL Z-scores and percentiles revealed similar results (0.04 vs 0.35, and 27.5 vs 33.3, p=0.023 and 0.042, respectively). Higher NfL concentrations were associated with more delayed P100 latencies (B=1.33 [0.26 to 2.41], p=0.015). In addition, improvement of P100 latencies between visits was associated with a trend for lower NfL values (B=0.003 [-0.0004 to 0.007], p=0.081). Based on a Cohen's d of 0.248, a future 1:1 parallel-arm placebo-controlled study using a remyelinating agent with comparable effect as clemastine would need 202 subjects per group to achieve 80% power. Conclusions: In pwMS, treatment with the remyelinating agent clemastine was associated with a reduction of blood NfL, suggesting that neuroprotection is achievable and measurable with therapeutic remyelination.
by: Ahmed Abdelhak, Christian Cordano, W John Boscardin, Eduardo Caverzasi, Jens Kuhle, Brandon Chan, Jeffrey M Gelfand, Hao H Yiu, Frederike C Oertel, Alexandra Beaudry-Richard, Shivany Condor Montes, Jorge R Oksenberg, Argentina Lario Lago, Adam Boxer, Julio C Rojas-Martinez, Fanny M Elahi, Jonah R Chan, Ari J Green
J Neurol Neurosurg Psychiatry. Jun 16, 2022
PMID: 35710320


Locus for severity implicates CNS resilience in progression of multiple sclerosis
abstract
Multiple sclerosis (MS) is an autoimmune disease of the central nervous system (CNS) that results in significant neurodegeneration in the majority of those affected and is a common cause of chronic neurological disability in young adults1,2. Here, to provide insight into the potential mechanisms involved in progression, we conducted a genome-wide association study of the age-related MS severity score in 12,584 cases and replicated our findings in a further 9,805 cases. We identified a significant association with rs10191329 in the DYSF-ZNF638 locus, the risk allele of which is associated with a shortening in the median time to requiring a walking aid of a median of 3.7 years in homozygous carriers and with increased brainstem and cortical pathology in brain tissue. We also identified suggestive association with rs149097173 in the DNM3-PIGC locus and significant heritability enrichment in CNS tissues. Mendelian randomization analyses suggested a potential protective role for higher educational attainment. In contrast to immune-driven susceptibility3, these findings suggest a key role for CNS resilience and potentially neurocognitive reserve in determining outcome in MS.
by: International Multiple Sclerosis Genetics Consortium; MultipleMS Consortium
Nature. Jul 2023
PMID: 37380766


Spinal Cord Atrophy Predicts Progressive Disease in Relapsing Multiple Sclerosis
abstract
Objective: A major challenge in multiple sclerosis (MS) research is the understanding of silent progression and Progressive MS. Using a novel method to accurately capture upper cervical cord area from legacy brain MRI scans we aimed to study the role of spinal cord and brain atrophy for silent progression and conversion to secondary progressive disease (SPMS). Methods: From a single-center observational study, all RRMS (n = 360) and SPMS (n = 47) patients and 80 matched controls were evaluated. RRMS patient subsets who converted to SPMS (n = 54) or silently progressed (n = 159), respectively, during the 12-year observation period were compared to clinically matched RRMS patients remaining RRMS (n = 54) or stable (n = 147), respectively. From brain MRI, we assessed the value of brain and spinal cord measures to predict silent progression and SPMS conversion. Results: Patients who developed SPMS showed faster cord atrophy rates (-2.19%/yr) at least 4 years before conversion compared to their RRMS matches (-0.88%/yr, p < 0.001). Spinal cord atrophy rates decelerated after conversion (-1.63%/yr, p = 0.010) towards those of SPMS patients from study entry (-1.04%). Each 1% faster spinal cord atrophy rate was associated with 69% (p < 0.0001) and 53% (p < 0.0001) shorter time to silent progression and SPMS conversion, respectively. Interpretation: Silent progression and conversion to secondary progressive disease are predominantly related to cervical cord atrophy. This atrophy is often present from the earliest disease stages and predicts the speed of silent progression and conversion to Progressive MS. Diagnosis of SPMS is rather a late recognition of this neurodegenerative process than a distinct disease phase.
by: Antje Bischof, Nico Papinutto, Anisha Keshavan, Anand Rajesh, Gina Kirkish, Xinheng Zhang, Jacob M Mallott, Carlo Asteggiano, Simone Sacco, Tristan J Gundel, Chao Zhao, William A Stern, Eduardo Caverzasi, Yifan Zhou, Refujia Gomez, Nicholas R Ragan, Adam Santaniello, Alyssa H Zhu, Jeremy Juwono, Carolyn J Bevan, Riley M Bove, Elizabeth Crabtree, Jeffrey M Gelfand, Douglas S Goodin, Jennifer S Graves, Ari J Green, Jorge R Oksenberg, Emmanuelle Waubant, Michael R Wilson, Scott S Zamvil; University of California, San Francisco MS-EPIC Team; Bruce A C Cree, Stephen L Hauser, Roland G Henry
Ann Neurol. Feb 2022
PMID: 34878197


Polygenic risk score association with multiple sclerosis susceptibility and phenotype in Europeans
abstract
Polygenic inheritance plays a pivotal role in driving multiple sclerosis susceptibility, an inflammatory demyelinating disease of the CNS. We developed polygenic risk scores (PRS) of multiple sclerosis and assessed associations with both disease status and severity in cohorts of European descent. The largest genome-wide association dataset for multiple sclerosis to date (n = 41 505) was leveraged to generate PRS scores, serving as an informative susceptibility marker, tested in two independent datasets, UK Biobank [area under the curve (AUC) = 0.73, 95% confidence interval (CI): 0.72-0.74, P = 6.41 × 10-146] and Kaiser Permanente in Northern California (KPNC, AUC = 0.8, 95% CI: 0.76-0.82, P = 1.5 × 10-53). Individuals within the top 10% of PRS were at higher than 5-fold increased risk in UK Biobank (95% CI: 4.7-6, P = 2.8 × 10-45) and 15-fold higher risk in KPNC (95% CI: 10.4-24, P = 3.7 × 10-11), relative to the median decile. The cumulative absolute risk of developing multiple sclerosis from age 20 onwards was significantly higher in genetically predisposed individuals according to PRS. Furthermore, inclusion of PRS in clinical risk models increased the risk discrimination by 13% to 26% over models based only on conventional risk factors in UK Biobank and KPNC, respectively. Stratifying disease risk by gene sets representative of curated cellular signalling cascades, nominated promising genetic candidate programmes for functional characterization. These pathways include inflammatory signalling mediation, response to viral infection, oxidative damage, RNA polymerase transcription, and epigenetic regulation of gene expression to be among significant contributors to multiple sclerosis susceptibility. This study also indicates that PRS is a useful measure for estimating susceptibility within related individuals in multicase families. We show a significant association of genetic predisposition with thalamic atrophy within 10 years of disease progression in the UCSF-EPIC cohort (P < 0.001), consistent with a partial overlap between the genetics of susceptibility and end-organ tissue injury. Mendelian randomization analysis suggested an effect of multiple sclerosis susceptibility on thalamic volume, which was further indicated to be through horizontal pleiotropy rather than a causal effect. In summary, this study indicates important, replicable associations of PRS with enhanced risk assessment and radiographic outcomes of tissue injury, potentially informing targeted screening and prevention strategies. Keywords: multiple sclerosis; pathway-specific risk score; phenotype association; polygenic risk score.
by: Hengameh Shams, Xiaorong Shao, Adam Santaniello, Gina Kirkish, Adil Harroud, Qin Ma, Noriko Isobe; University of California San Francisco MS-EPIC Team; Catherine A Schaefer, Jacob L McCauley, Bruce A C Cree, Alessandro Didonna, Sergio E Baranzini, Nikolaos A Patsopoulos, Stephen L Hauser, Lisa F Barcellos, Roland G Henry, Jorge R Oksenberg
Brain. Feb 13, 2023
PMID: 35253861


High-resolution DNA methylation screening of the major histocompatibility complex in multiple sclerosis
abstract
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. Keywords: DNA methylation; DNA methylation quantitative trait loci; differentially methylated regions; major histocompatibility complex; multiple sclerosis.
by: Qin Ma, Danillo G Augusto, Gonzalo Montero-Martin, Stacy J Caillier, Kazutoyo Osoegawa, Bruce A C Cree, Stephen L Hauser, Alessandro Didonna, Jill A Hollenbach, Paul J Norman, Marcelo Fernandez-Vina, Jorge R Oksenberg
Front Neurol. Dec 8, 2023
PMID: 38145128


Ancestral risk modification for multiple sclerosis susceptibility detected across the Major Histocompatibility Complex in a multi-ethnic population
abstract
The Major Histocompatibility Complex (MHC) makes the largest genetic contribution to multiple sclerosis (MS) susceptibility, with 32 independent effects across the region explaining 20% of the heritability in European populations. Variation is high across populations with allele frequency differences and population-specific risk alleles identified. We sought to identify MHC-specific MS susceptibility variants and assess the effect of ancestral risk modification within 2652 Latinx and Hispanic individuals as well as 2435 Black and African American individuals. We have identified several novel susceptibility alleles which are rare in European populations including HLA-B*53:01, and we have utilized the differing linkage disequilibrium patterns inherent to these populations to identify an independent role for HLA-DRB1*15:01 and HLA-DQB1*06:02 on MS risk. We found a decrease in Native American ancestry in MS cases vs controls across the MHC, peaking near the previously identified MICB locus with a decrease of ~5.5% in Hispanics and ~0.4% in African Americans. We have identified several susceptibility variants, including within the MICB gene region, which show global ancestry risk modification and indicate ancestral differences which may be due in part to correlated environmental factors. We have also identified several susceptibility variants for which MS risk is modified by local ancestry and indicate true ancestral genetic differences; including HLA-DQB1*06:02 for which MS risk for European allele carriers is almost two times the risk for African allele carriers. These results validate the importance of investigating MS susceptibility at an ancestral level and offer insight into the epidemiology of MS phenotypic diversity.
by: Ashley H Beecham, Lilyana Amezcua, Angel Chinea, Clara P Manrique, Lissette Gomez, Andrea Martinez, Gary W Beecham, Nikolaos A Patsopoulos, Tanuja Chitnis, Howard L Weiner, Philip L De Jager, Esteban G Burchard, Brett T Lund, Kathryn C Fitzgerald, Peter A Calabresi, Silvia R Delgado, Jorge R Oksenberg, Jacob L McCauley
PLoS One. Dec 22, 2022
PMID: 36548255


Specific hypomethylation programs underpin B cell activation in early multiple sclerosis
abstract
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. Keywords: B cell; hypomethylation; multiple sclerosis.
by: Qin Ma, Stacy J Caillier, Shaun Muzic; University of California San Francisco MS-EPIC Team; Michael R Wilson, Roland G Henry, Bruce A C Cree, Stephen L Hauser, Alessandro Didonna, Jorge R Oksenberg
Proc Natl Acad Sci U S A. Dec 21, 2021
PMID: 34911760


Neurofilament Light Chain Elevation and Disability Progression in Multiple Sclerosis
abstract
Importance: Mechanisms contributing to disability accumulation in multiple sclerosis (MS) are poorly understood. Blood neurofilament light chain (NfL) level, a marker of neuroaxonal injury, correlates robustly with disease activity in people with MS (MS); however, data on the association between NfL level and disability accumulation have been conflicting. Objective: To determine whether and when NfL levels are elevated in the context of confirmed disability worsening (CDW). Design, setting, and participants: This study included 2 observational cohorts: results from the Expression, Proteomics, Imaging, Clinical (EPIC) study at the University of California San Francisco (since 2004) were confirmed in the Swiss Multiple Sclerosis Cohort (SMSC), a multicenter study in 8 centers since 2012. Data were extracted from EPIC in April 2022 (sampling July 1, 2004, to December 20, 2016) and SMSC in December 2022 (sampling June 6, 2012, to September 2, 2021). The study included 2 observational cohorts in tertiary MS centers. All participants of both cohorts with available NfL results were included in the study, and no eligible participants were excluded or declined to participate. Exposure: Association between NfL z scores and CDW. Main outcome measures: CDW was defined as Expanded Disability Status Scale (EDSS) worsening that was confirmed after 6 or more months and classified into CDW associated with clinical relapses (CDW-R) or independent of clinical relapses (CDW-NR). Visits were classified in relation to the disability worsening events into CDW(-2) for 2 visits preceding event, CDW(-1) for directly preceding event, CDW(event) for first diagnosis of EDSS increase, and the confirmation visit. Mixed linear and Cox regression models were used to evaluate NfL dynamics and to assess the association of NfL with future CDW, respectively. Results: A total of 3906 EPIC visits (609 participants; median [IQR] age, 42.0 [35.0-50.0] years; 424 female [69.6%]) and 8901 SMSC visits (1290 participants; median [IQR] age, 41.2 [32.5-49.9] years; 850 female [65.9%]) were included. In CDW-R (EPIC, 36 events; SMSC, 93 events), NfL z scores were 0.71 (95% CI, 0.35-1.07; P < .001) units higher at CDW-R(-1) in EPIC and 0.32 (95% CI, 0.14-0.49; P < .001) in SMSC compared with stable MS samples. NfL elevation could be detected preceding CDW-NR (EPIC, 191 events; SMSC, 342 events) at CDW-NR(-2) (EPIC: 0.23; 95% CI, 0.01-0.45; P = .04; SMSC: 0.28; 95% CI, 0.18-0.37; P < .001) and at CDW-NR(-1) (EPIC: 0.27; 95% CI, 0.11-0.44; P < .001; SMSC: 0.09; 95% CI, 0-0.18; P = .06). Those findings were replicated in the subgroup with relapsing-remitting MS. Time-to-event analysis confirmed the association between NfL levels and future CDW-R within approximately 1 year and CDW-NR (in approximately 1-2 years). Conclusions and relevance: This cohort study documents the occurrence of NfL elevation in advance of clinical worsening and may hint to a potential window of ongoing dynamic central nervous system pathology that precedes the diagnosis of CDW.
by: Ahmed Abdelhak, Pascal Benkert, Sabine Schaedelin, W John Boscardin, Christian Cordano, Johanna Oechtering, Kirtana Ananth, Cristina Granziera, Lester Melie-Garcia, Shivany Condor Montes, Alexandra Beaudry-Richard, Lutz Achtnichts, Frederike C Oertel, Patrice H Lalive, David Leppert, Stefanie Müller, Roland G Henry, Caroline Pot, Amandine Matthias, Anke Salmen, Jorge R Oksenberg, Giulio Disanto, Chiara Zecca, Marcus D'Souza, Renaud Du Pasquier, Claire Bridel, Claudio Gobbi, Ludwig Kappos, Stephen L Hauser, Bruce A C Cree, Jens Kuhle, Ari J Green, UCSF, MS EPIC, and the SMSC Study Teams
JAMA Neurol. Dec 1, 2023
PMID: 37930670


Genetic susceptibility to multiple sclerosis in African Americans
abstract
Objective: To explore the nature of genetic-susceptibility to multiple sclerosis (MS) in African-Americans. Background: Recently, the number of genetic-associations with MS has exploded although the MS-associations of specific haplotypes within the major histocompatibility complex (MHC) have been known for decades. For example, the haplotypes HLA-DRB1*15:01~HLA-DQB1*06:02, and HLA-DRB1*03:01~ HLA-DQB1*02:01 have odds ratios (ORs) for an MS-association orders of magnitude stronger than many of these newly-discovered associations. Nevertheless, all these haplotypes are part of much larger conserved extended haplotypes (CEHs), which span both the Class I and Class II MHC regions. African-Americans are at greater risk of developing MS compared to a native Africans but at lesser risk compared to Europeans. It is the purpose of this manuscript to explore the relationship between MS-susceptibility and the CEH make-up of our African-American cohort. Design/methods: The African-American (AA) cohort consisted of 1,305 patients with MS and 1,155 controls, who self-identified as being African-American. For comparison, we used the 18,492 controls and 11,144 MS-cases from the predominantly European Wellcome Trust Case Control Consortium (WTCCC) and the 28,557 phased native Africans from the multinational "Be the Match" registry. The WTCCC and the African-Americans were phased at each of five HLA loci (HLA-A, HLA-C, HLA-B, HLA-DRB1 and HLA-DQB1) and the at 11 SNPs (10 of which were in non-coding regions) surrounding the Class II region of the DRB1 gene using previously-published probabilistic phasing algorithms. Results: Of the 32 most frequent CEHs, 18 (56%) occurred either more frequently or exclusively in Africans) whereas 9 (28%) occurred more frequently or exclusively in Europeans. The remaining 5 CEHs occurred in neither control group although, likely, these were African in origin. Eight of these CEHs carried the DRB1*15:03~DQB1*06:02~a36 haplotype and three carried the DRB1*15:01~DQB1*06:02~a1 haplotype. In African Americans, a single-copy of the European CEH (03:01_07:02_07:02_15:01_06:02_a1) was associated with considerable MS-risk (OR = 3.30; p = 0.0001)-similar to that observed in the WTCCC (OR = 3.25; p<10-168). By contrast, the MS-risk for the European CEH (02:01_07:02_07:02_15:01_06:02_a1) was less (OR = 1.49; ns)-again, similar to the WTCCC (OR = 2.2; p<10-38). Moreover, four African haplotypes were "protective" relative to a neutral reference, to three European CEHs, and also to the five other African CEHs. Conclusions: The common CEHs in African Americans are divisible into those that are either African or European in origin, which are derived without modification from their source population. European CEHs, linked to MS-risk, in general, had similar impacts in African-Americans as they did in Europeans. By contrast, African CEHs had mixed MS-risks. For a few, the MS-risk exceeded that in a neutral-reference group whereas, for many others, these CEHs were "protective"-perhaps providing a partial rationale for the lower MS-risk in African-Americans compared to European-Americans.
by: Douglas S Goodin, Jorge R Oksenberg, Venceslas Douillard, Pierre-Antoine Gourraud, Nicolas Vince
PLoS One. Aug 9, 2021
PMID: 34370753


Shared Genetic Risk Factors for Multiple Sclerosis/Psoriasis Suggest Involvement of Interleukin-17 and Janus Kinase-Signal Transducers and Activators of Transcription Signaling
abstract
Objective: Psoriasis and multiple sclerosis (MS) are complex immune diseases that are mediated by T cells and share multiple comorbidities. Previous studies have suggested psoriatic patients are at higher risk of MS; however, causal relationships between the two conditions remain unclear. Through epidemiology and genetics, we provide a comprehensive understanding of the relationship, and share molecular factors between psoriasis and MS. Methods: We used logistic regression, trans-disease meta-analysis and Mendelian randomization. Medical claims data were included from 30 million patients, including 141,544 with MS and 742,919 with psoriasis. We used genome-wide association study summary statistics from 11,024 psoriatic, 14,802 MS cases, and 43,039 controls for trans-disease meta-analysis, with additional summary statistics from 5 million individuals for Mendelian randomization. Results: Psoriatic patients have a significantly higher risk of MS (4,637 patients with both diseases; odds ratio [OR] 1.07, p = 1.2 × 10-5 ) after controlling for potential confounders. Using inverse variance and equally weighted trans-disease meta-analysis, we revealed >20 shared and opposing (direction of effect) genetic loci outside the major histocompatibility complex that showed significant genetic colocalization (in COLOC and COLOC-SuSiE v5.1.0). Co-expression analysis of genes from these loci further identified distinct clusters that were enriched among pathways for interleukin-17/tumor necrosis factor-α (OR >39, p < 1.6 × 10-3 ) and Janus kinase-signal transducers and activators of transcription (OR 35, p = 1.1 × 10-5 ), including genes, such as TNFAIP3, TYK2, and TNFRSF1A. Mendelian randomization found psoriasis as an exposure has a significant causal effect on MS (OR 1.04, p = 5.8 × 10-3 ), independent of type 1 diabetes (OR 1.05, p = 4.3 × 10-7 ), type 2 diabetes (OR 1.08, p = 2.3 × 10-3 ), inflammatory bowel disease (OR 1.11, p = 1.6 × 10-11 ), and vitamin D level (OR 0.75, p = 9.4 × 10-3 ). Interpretation: By investigating the shared genetics of psoriasis and MS, along with their modifiable risk factors, our findings will advance innovations in treatment for patients suffering from comorbidities.
by: Matthew T Patrick, Rajan P Nair, Kevin He, Philip E Stuart, Allison C Billi, Xiang Zhou, Johann E Gudjonsson, Jorge R Oksenberg, James T Elder, Lam C Tsoi
Ann Neurol. Aug 2023
PMID: 37127916


Epigenetic control of ataxin-1 in multiple sclerosis
abstract
Objective: ATXN1 encodes the polyglutamine protein ataxin-1, which we have demonstrated exerting an immunomodulatory function in the context of central nervous system (CNS) autoimmunity, in addition to its classical role in the neurodegenerative disorder spinocerebellar ataxia type 1 (SCA1). In this study, we dissected the contribution of DNA methylation to the regulation of ATXN1 in multiple sclerosis (MS). Methods: We interrogated a DNA methylation dataset previously generated via bisulfate DNA sequencing (BS-seq) in sorted peripheral immune cytotypes (CD4+ and CD8+ T cells, CD19+ B cells, and CD14+ monocytes) isolated from untreated MS patients at symptoms onset. Results: Here, we report that ATXN1 undergoes hypo-methylation at four distinct regions upon MS, exclusively in B cells. We also highlight how these differentially methylated sites overlap with other regulatory epigenetic marks and MS risk variants. Lastly, we employ luciferase assays to assess the functionality of these regions, showing that the loss of methylation leads to an increase in ATXN1 expression. Interpretation: Altogether, these findings provide biological insights into ataxin-1 regulation in the immune system as well as into the molecular mechanisms underlying MS risk.
by: Qin Ma, Jorge R Oksenberg, Alessandro Didonna
Ann Clin Transl Neurol. Aug 2022
PMID: 35903875


An autoantibody signature predictive for multiple sclerosis
abstract
Although B cells are implicated in multiple sclerosis (MS) pathophysiology, a predictive or diagnostic autoantibody remains elusive. In this study, the Department of Defense Serum Repository (DoDSR), a cohort of over 10 million individuals, was used to generate whole-proteome autoantibody profiles of hundreds of patients with MS (PwMS) years before and subsequently after MS onset. This analysis defines a unique cluster in approximately 10% of PwMS who share an autoantibody signature against a common motif that has similarity with many human pathogens. These patients exhibit antibody reactivity years before developing MS symptoms and have higher levels of serum neurofilament light (sNfL) compared to other PwMS. Furthermore, this profile is preserved over time, providing molecular evidence for an immunologically active preclinical period years before clinical onset. This autoantibody reactivity was validated in samples from a separate incident MS cohort in both cerebrospinal fluid and serum, where it is highly specific for patients eventually diagnosed with MS. This signature is a starting point for further immunological characterization of this MS patient subset and may be clinically useful as an antigen-specific biomarker for high-risk patients with clinically or radiologically isolated neuroinflammatory syndromes.
by: Colin R Zamecnik, Gavin M Sowa, Ahmed Abdelhak, Ravi Dandekar, Rebecca D Bair, Kristen J Wade, Christopher M Bartley, Kerry Kizer, Danillo G Augusto, Asritha Tubati, Refujia Gomez, Camille Fouassier, Chloe Gerungan, Colette M Caspar, Jessica Alexander, Anne E Wapniarski, Rita P Loudermilk, Erica L Eggers, Kelsey C Zorn, Kirtana Ananth, Nora Jabassini, Sabrina A Mann, Nicholas R Ragan, Adam Santaniello, Roland G Henry, Sergio E Baranzini, Scott S Zamvil, Joseph J Sabatino Jr, Riley M Bove, Chu-Yueh Guo, Jeffrey M Gelfand, Richard Cuneo, H-Christian von Büdingen, Jorge R Oksenberg, Bruce A C Cree, Jill A Hollenbach, Ari J Green, Stephen L Hauser, Mitchell T Wallin, Joseph L DeRisi, Michael R Wilson
Nat Med. Apr 19, 2024
PMID: 38641750


Synaptic injury in the inner plexiform layer of the retina is associated with progression in multiple sclerosis
abstract
While neurodegeneration underlies the pathological basis for permanent disability in multiple sclerosis (MS), predictive biomarkers for progression are lacking. Using an animal model of chronic MS, we find that synaptic injury precedes neuronal loss and identify thinning of the inner plexiform layer (IPL) as an early feature of inflammatory demyelination-prior to symptom onset. As neuronal domains are anatomically segregated in the retina and can be monitored longitudinally, we hypothesize that thinning of the IPL could represent a biomarker for progression in MS. Leveraging our dataset with over 800 participants enrolled for more than 12 years, we find that IPL atrophy directly precedes progression and propose that synaptic loss is predictive of functional decline. Using a blood proteome-wide analysis, we demonstrate a strong correlation between demyelination, glial activation, and synapse loss independent of neuroaxonal injury. In summary, monitoring synaptic injury is a biologically relevant approach that reflects a potential driver of progression. Keywords: EAE; blood biomarkers; demyelination; experimental autoimmune encephalomyelitis; optical coherence tomography; retina; serum proteomics; synaptic injury.
by: Christian Cordano, Sebastian Werneburg, Ahmed Abdelhak, Daniel J Bennett, Alexandra Beaudry-Richard, Greg J Duncan, Frederike C Oertel, W John Boscardin, Hao H Yiu, Nora Jabassini, Lauren Merritt, Sonia Nocera, Jung H Sin, Isaac P Samana, Shivany Y Condor Montes, Kirtana Ananth, Antje Bischof, UCSF MS EPIC Team, Bardia Nourbakhsh, Stephen L Hauser, Bruce A C Cree, Ben Emery, Dorothy P Schafer, Jonah R Chan, Ari J Green
Cell Rep Med. Apr 16, 2024
PMID: 38574736


A short HLA-DRA isoform binds the HLA-DR2 heterodimer on the outer domain of the peptide-binding site
abstract
The human leukocyte antigen (HLA) locus encodes a large group of proteins governing adaptive and innate immune responses. Among them, HLA class II proteins form α/β heterodimers on the membrane of professional antigen-presenting cells (APCs), where they display both, self and pathogen-derived exogenous antigens to CD4+ T lymphocytes. We have previously shown that a shorter HLA-DRA isoform (sHLA-DRA) lacking 25 amino acids can be presented onto the cell membrane via binding to canonical HLA-DR2 heterodimers. Here, we employed atomistic molecular dynamics simulations to decipher the binding position of sHLA-DRA and its structural impact on functional regions of the HLA-DR2 molecule. We show that a loop region exposed only in the short isoform (residues R69 to G83) is responsible for binding to the outer domain of the HLA-DR2 peptide-binding site, and experimentally validated the critical role of F76 in mediating such interaction. Additionally, sHLA-DRA allosterically modifies the peptide-binding pocket conformation. In summary, this study unravels key molecular mechanisms underlying sHLA-DRA function, providing important insights into the role of full-length proteins in structural modulation of HLA class II receptors. Keywords: Antigen presentation; Molecular dynamics; Protein-protein binding; Structural modulation.
by: Hengameh Shams, Jill A Hollenbach, Atsuko Matsunaga, Mohammad R K Mofrad, Jorge R Oksenberg, Alessandro Didonna
Arch Biochem Biophys. Apr 15, 2022
PMID: 35218721


The 18th International HLA & Immunogenetics workshop project report: Creating fully representative MHC reference haplotypes
abstract
A fundamental endeavor of the International Histocompatibility and Immunogenetics Workshop (IHIW) was assembling a collection of DNA samples homozygous through the MHC genomic region. This collection proved invaluable for assay development in the histocompatibility and immunogenetics field, for generating the human reference genome, and furthered our understanding of MHC diversity. Defined by their HLA-A, -B, -C and -DRB1 alleles, the combined frequency of the haplotypes from these individuals is ~20% in Europe. Thus, a significant proportion of MHC haplotypes, both common and rare throughout the world, and including many associated with disease, are not yet represented. In this workshop component, we are collecting the next generation of MHC-homozygous samples, to expand, diversify and modernize this critical community resource that has been foundational to the field. We asked laboratories worldwide to identify samples homozygous through all HLA class I and/or HLA class II genes, or through whole-genome SNP genotyping or sequencing, to have extensive homozygosity tracts within the MHC region. The focus is non-Europeans or those having HLA haplotypes less common in Europeans. Through this effort, we have obtained samples from 537 individuals representing 294 distinct haplotypes, as determined by their HLA class I and II alleles, and an additional 50 haplotypes distinct in HLA class I or II alleles. Although we have expanded the diversity, many populations remain underrepresented, particularly from Africa, and we encourage further participation. The data will serve as a resource for investigators seeking to characterize variation across the MHC genomic region for disease and population studies. Keywords: HLA; IHIW; MHC; ancestry; assembly; ethnicity; homozygous; long‐read; reference haplotype; sequencing; short‐read.
by: Pollock NR, Farias TDJ, Kichula KM, Sauter J, Scholz S, Nii-Trebi NI, Khor SS, Tokunaga K, Voorter CE, Groeneweg M, Augusto DG, Arrieta-Bolaños E, Mayor NP, Edinur HA, ElGhazali G, Issler HC, Petzl-Erler ML, Oksenberg JR, Marin WM, Hollenbach JA, Gendzekhadze K, Cita R, Stelet V, Rajalingam R, Koskela S, Clancy J, Chatzistamatiou T, Houwaart T, Kulski J, Guethlein LA, Parham P, Schmidt AH, Dilthey A, Norman PJ.
HLA. June 2024
PMID: 38923286


Progressive multifocal leukoencephalopathy genetic risk variants for pharmacovigilance of immunosuppressant therapies
abstract
Background: Progressive multifocal leukoencephalopathy (PML) is a rare and often lethal brain disorder caused by the common, typically benign polyomavirus 2, also known as JC virus (JCV). In a small percentage of immunosuppressed individuals, JCV is reactivated and infects the brain, causing devastating neurological defects. A wide range of immunosuppressed groups can develop PML, such as patients with: HIV/AIDS, hematological malignancies (e.g., leukemias, lymphomas, and multiple myeloma), autoimmune disorders (e.g., psoriasis, rheumatoid arthritis, and systemic lupus erythematosus), and organ transplants. In some patients, iatrogenic (i.e., drug-induced) PML occurs as a serious adverse event from exposure to immunosuppressant therapies used to treat their disease (e.g., hematological malignancies and multiple sclerosis). While JCV infection and immunosuppression are necessary, they are not sufficient to cause PML. Methods: We hypothesized that patients may also have a genetic susceptibility from the presence of rare deleterious genetic variants in immune-relevant genes (e.g., those that cause inborn errors of immunity). In our prior genetic study of 184 PML cases, we discovered 19 candidate PML risk variants. In the current study of another 152 cases, we validated 4 of 19 variants in both population controls (gnomAD 3.1) and matched controls (JCV+ multiple sclerosis patients on a PML-linked drug ≥ 2 years). Results: The four variants, found in immune system genes with strong biological links, are: C8B, 1-57409459-C-A, rs139498867; LY9 (alias SLAMF3), 1-160769595-AG-A, rs763811636; FCN2, 9-137779251-G-A, rs76267164; STXBP2, 19-7712287-G-C, rs35490401. Carriers of any one of these variants are shown to be at high risk of PML when drug-exposed PML cases are compared to drug-exposed matched controls: P value = 3.50E-06, OR = 8.7 [3.7-20.6]. Measures of clinical validity and utility compare favorably to other genetic risk tests, such as BRCA1 and BRCA2 screening for breast cancer risk and HLA-B*15:02 pharmacogenetic screening for pharmacovigilance of carbamazepine to prevent Stevens-Johnson Syndrome and Toxic Epidermal Necrolysis. Conclusion: For the first time, a PML genetic risk test can be implemented for screening patients taking or considering treatment with a PML-linked drug in order to decrease the incidence of PML and enable safer use of highly effective therapies used to treat their underlying disease. Keywords: JC virus; PML; immunodeficiency; multiple sclerosis; natalizumab; pharmacovigilance; progressive multifocal leukoencephalopathy; serious adverse event.
by: Eli Hatchwell, Edward B Smith 3rd, Shapour Jalilzadeh, Christopher D Bruno, Yassine Taoufik, Houria Hendel-Chavez, Roland Liblau, David Brassat, Guillaume Martin-Blondel, Heinz Wiendl, Nicholas Schwab, Irene Cortese, Maria Chiara Monaco, Luisa Imberti, Ruggero Capra, Jorge R Oksenberg, Jacques Gasnault, Bruno Stankoff, Todd A Richmond, David M Rancour, Igor J Koralnik, Barbara A Hanson, Eugene O Major, Christina R Chow, Peggy S Eis
Front Neurol. Dec 14, 2022
PMID: 36588876