Single nuclei RNAseq stratifies multiple sclerosis patients into three distinct white matter glia responses

Abstract

AbstractThe lack of understanding as to the cellular and molecular basis of clinical and genetic heterogeneity in progressive multiple sclerosis (MS) has hindered the search for new effective therapies and biomarkers. Here, to address this gap, we analysed 740,000 single nuclei RNAseq profiles of 165 samples of white matter (WM) lesions, normal appearing WM, grey matter (GM) lesions and normal appearing GM from 55 MS patients and 28 controls. We find that gene expression changes in response to MS are highly cell-type specific in WM and GM lesions but are largely shared within an individual cell-type across lesions, following a continuum rather than discrete lesion-specific molecular programs. The major biological determinants of variability in gene expression in MS samples relate to individual patient effects, rather than to lesion types or other metadata. Using multi-omics factor analysis (MOFA+), we identify three subgroups of MS patients with distinct oligodendrocyte composition and WM glial gene expression signatures, suggestive of engagement of different pathological/regenerative processes. The discovery of these three patterns significantly advances our mechanistic understanding of progressive MS, provides a framework to use molecular biomarkers to stratify patients for best therapeutic approaches for progressive MS, and highlights the need for precision-medicine approaches to address heterogeneity among MS patients.