Abstract
Multiple sclerosis (MS) is an autoimmune disease characterized by demyelination of the central nervous system (CNS). Autologous hematopoietic cell transplantation (HCT) shows promising benefits for relapsing–remitting MS in open-label clinical studies, but the cellular mechanisms underlying its therapeutic effects remain unclear. Using single-nucleus RNA sequencing, we identify a reactive myeloid cell state in chronic experimental autoimmune encephalitis (EAE) associated with neuroprotection and immune suppression. HCT in EAE mice results in an increase of the neuroprotective myeloid state, improvement of neurological deficits, reduced number of demyelinated lesions, decreased number of effector T cells and amelioration of reactive astrogliosis. Enhancing myeloid cell incorporation after a modified HCT further improved these neuroprotective effects. These data suggest that myeloid cell manipulation or replacement may be an effective therapeutic strategy for chronic inflammatory conditions of the CNS.
This is a preview of subscription content, access via your institution
Access options
Access Nature and 54 other Nature Portfolio journals
Get Nature+, our best-value online-access subscription
$29.99 / 30 days
cancel any time
Subscribe to this journal
Receive 12 print issues and online access
$209.00 per year
only $17.42 per issue
Buy this article
- Purchase on Springer Link
- Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout
Similar content being viewed by others
Data availability
The murine sequencing datasets generated and analyzed in the current study are available in the Gene Expression Omnibus repository under accession numbers GSE217529 and GSE242512. The human spinal cord dataset29 was accessed as six fastq files (three control and three MS samples) through Sequence Read Archive (accession number PRJNA726991). Data from previously published human brain NucSeq datasets ((three control and five MS samples; Cell Browser dataset ID, chronic-ms)30; (five control and four MS samples; Cell Browser dataset ID, oligodendrocyte-ms)32; and (nine control and twelve MS samples; Cell Browser dataset ID, ms)31) were obtained as expression matrices from the UCSC Cell Browser (https://cells.ucsc.edu)67. Reference genomes (refdata-gex-mm10-2020-A and refdata-gex-GRCh38-2020-A) were accessed at www.10xgenomics.com. Source data are provided with this paper.
Code availability
Computational analyses were performed using freely available software packages as described in the Methods. Central to the analysis pipeline of NucSeq data in this study are the software packages 10× Genomics Cell Ranger (v.6.1.2; https://www.10xgenomics.com/support/software/cell-ranger)70 and R packages Seurat (v.4.2.0)68, sctransform (v.0.3.5)61, fgsea (v.1.24.0)56, MAST (v.1.20.0)66, VISION (v.3.0.0)26 and biomaRt (v.2.50.3)69. Additional information regarding the computational methodology of this project is available from the corresponding author upon reasonable request.
References
Wallin, M. T. et al. Global, regional, and national burden of multiple sclerosis 1990–2016: a systematic analysis for the Global Burden of Disease Study 2016. Lancet Neurol. 18, 269–285 (2019).
Miller, A. E. et al. Autologous hematopoietic stem cell transplant in multiple sclerosis: recommendations of the national multiple sclerosis society. JAMA Neurol. 78, 241–246 (2021).
Noseworthy, J. H., Lucchinetti, C., Rodriguez, M. & Weinshenker, B. G. Multiple sclerosis. N. Engl. J. Med. 343, 938–952 (2000).
Yong, H. Y. F. & Yong, V. W. Mechanism-based criteria to improve therapeutic outcomes in progressive multiple sclerosis. Nat. Rev. Neurol. 18, 40–55 (2022).
Mancardi, G. L. et al. Autologous hematopoietic stem cell transplantation in multiple sclerosis: a phase II trial. Neurology 84, 981–988 (2015).
Nash, R. A. et al. High-dose immunosuppressive therapy and autologous hematopoietic cell transplantation for relapsing-remitting multiple sclerosis (HALT-MS): a 3-year interim report. JAMA Neurol. 72, 159–169 (2015).
Atkins, H. L. et al. Immunoablation and autologous haemopoietic stem-cell transplantation for aggressive multiple sclerosis: a multicentre single-group phase 2 trial. Lancet 388, 576–585 (2016).
Moore, J. J. et al. Prospective phase II clinical trial of autologous haematopoietic stem cell transplant for treatment refractory multiple sclerosis. J. Neurol. Neurosurg. Psychiatry 90, 514–521 (2019).
Burt, R. K. et al. Effect of nonmyeloablative hematopoietic stem cell transplantation vs continued disease-modifying therapy on disease progression in patients with relapsing–remitting multiple sclerosis: a randomized clinical trial. JAMA 321, 165–174 (2019).
Muraro, P. A. et al. Thymic output generates a new and diverse TCR repertoire after autologous stem cell transplantation in multiple sclerosis patients. J. Exp. Med. 201, 805–816 (2005).
Muraro, P. A. et al. T cell repertoire following autologous stem cell transplantation for multiple sclerosis. J. Clin. Invest. 124, 1168–1172 (2014).
Massey, J. et al. Haematopoietic stem cell transplantation results in extensive remodelling of the clonal T cell repertoire in multiple sclerosis. Front. Immunol. 13, 798300 (2022).
Ruder, J. et al. Dynamics of T cell repertoire renewal following autologous hematopoietic stem cell transplantation in multiple sclerosis. Sci. Transl. Med. 14, eabq1693 (2022).
Gold, R. Understanding pathogenesis and therapy of multiple sclerosis via animal models: 70 years of merits and culprits in experimental autoimmune encephalomyelitis research. Brain 129, 1953–1971 (2006).
Psenicka, M. W., Smith, B. C., Tinkey, R. A. & Williams, J. L. Connecting neuroinflammation and neurodegeneration in multiple sclerosis: Are oligodendrocyte precursor cells a nexus of disease? Front. Cell. Neurosci. 15, 654284 (2021).
Williams, J. L. et al. Astrocyte–T cell crosstalk regulates region‐specific neuroinflammation. Glia 68, 1361–1374 (2020).
Meijer, M. et al. Epigenomic priming of immune genes implicates oligodendroglia in multiple sclerosis susceptibility. Neuron 110, 1193–1210.e13 (2022).
Ajami, B. et al. Single-cell mass cytometry reveals distinct populations of brain myeloid cells in mouse neuroinflammation and neurodegeneration models. Nat. Neurosci. 21, 541–551 (2018).
Guerrero, B. L. & Sicotte, N. L. Microglia in multiple sclerosis: Friend or foe? Front. Immunol. 11, 374 (2020).
Miron, V. E. et al. M2 microglia and macrophages drive oligodendrocyte differentiation during CNS remyelination. Nat. Neurosci. 16, 1211–1218 (2013).
Lampron, A. et al. Inefficient clearance of myelin debris by microglia impairs remyelinating processes. J. Exp. Med. 212, 481–495 (2015).
Falcão, A. M. et al. Disease-specific oligodendrocyte lineage cells arise in multiple sclerosis. Nat. Med. 24, 1837–1844 (2018).
Jordão, M. J. C. et al. Single-cell profiling identifies myeloid cell subsets with distinct fates during neuroinflammation. Science 363, eaat7554 (2019).
Sailor, K. A. et al. Hematopoietic stem cell transplantation chemotherapy causes microglia senescence and peripheral macrophage engraftment in the brain. Nat. Med. 28, 517–527 (2022).
Shibuya, Y. et al. Treatment of a genetic brain disease by CNS-wide microglia replacement. Sci. Transl. Med. 14, eabl9945 (2022).
DeTomaso, D. et al. Functional interpretation of single cell similarity maps. Nat. Commun. 10, 4376 (2019).
Hahn, O. et al. Atlas of the aging mouse brain reveals white matter as vulnerable foci. Cell 186, 4117–4133.e22 (2023).
International Multiple Sclerosis Genetics Consortium et al. Multiple sclerosis genomic map implicates peripheral immune cells and microglia in susceptibility. Science 365, eaav7188 (2019).
Trobisch, T. et al. Cross-regional homeostatic and reactive glial signatures in multiple sclerosis. Acta Neuropathol. 144, 987–1003 (2022).
Absinta, M. et al. A lymphocyte–microglia–astrocyte axis in chronic active multiple sclerosis. Nature 597, 709–714 (2021).
Schirmer, L. et al. Neuronal vulnerability and multilineage diversity in multiple sclerosis. Nature 573, 75–82 (2019).
Jäkel, S. et al. Altered human oligodendrocyte heterogeneity in multiple sclerosis. Nature 566, 543–547 (2019).
Hohsfield, L. A. et al. Effects of long-term and brain-wide colonization of peripheral bone marrow-derived myeloid cells in the CNS. J. Neuroinflammation 17, 279 (2020).
Xu, Z. et al. Efficient strategies for microglia replacement in the central nervous system. Cell Rep. 32, 108041 (2020).
Cronk, J. C. et al. Peripherally derived macrophages can engraft the brain independent of irradiation and maintain an identity distinct from microglia. J. Exp. Med. 215, 1627–1647 (2018).
Lee, H. K. et al. Daam2-PIP5K Is a regulatory pathway for Wnt signaling and therapeutic target for remyelination in the CNS. Neuron 85, 1227–1243 (2015).
Liddelow, S. A. et al. Neurotoxic reactive astrocytes are induced by activated microglia. Nature 541, 481–487 (2017).
Escartin, C. et al. Reactive astrocyte nomenclature, definitions, and future directions. Nat. Neurosci. 24, 312–325 (2021).
Haimon, Z. et al. Cognate microglia–T cell interactions shape the functional regulatory T cell pool in experimental autoimmune encephalomyelitis pathology. Nat. Immunol. 23, 1749–1762 (2022).
Van Hove, H. et al. A single-cell atlas of mouse brain macrophages reveals unique transcriptional identities shaped by ontogeny and tissue environment. Nat. Neurosci. 22, 1021–1035 (2019).
Bennett, F. C. et al. A combination of ontogeny and CNS environment establishes microglial identity. Neuron 98, 1170–1183.e8 (2018).
Hammond, T. R. et al. Single-cell RNA sequencing of microglia throughout the mouse lifespan and in the injured brain reveals complex cell-state changes. Immunity 50, 253–271.e6 (2019).
Keren-Shaul, H. et al. A unique microglia type associated with restricting development of Alzheimer’s disease. Cell 169, 1276–1290.e17 (2017).
Shen, X., Qiu, Y., Wight, A. E., Kim, H.-J. & Cantor, H. Definition of a mouse microglial subset that regulates neuronal development and proinflammatory responses in the brain. Proc. Natl Acad. Sci. 119, e2116241119 (2022).
Schroeter, M., Zickler, P., Denhardt, D. T., Hartung, H.-P. & Jander, S. Increased thalamic neurodegeneration following ischaemic cortical stroke in osteopontin-deficient mice. Brain 129, 1426–1437 (2006).
Tagliabracci, V. S. et al. A single kinase generates the majority of the secreted phosphoproteome. Cell 161, 1619–1632 (2015).
Masuda, T. et al. Spatial and temporal heterogeneity of mouse and human microglia at single-cell resolution. Nature 566, 388–392 (2019).
Abel, A. M. et al. IQGAP1: insights into the function of a molecular puppeteer. Mol. Immunol. 65, 336–349 (2015).
Okuyama, Y. et al. IQGAP1 restrains T-cell cosignaling mediated by OX40. FASEB J. 34, 540–554 (2020).
Huang, Y. et al. Repopulated microglia are solely derived from the proliferation of residual microglia after acute depletion. Nat. Neurosci. 21, 530–540 (2018).
Elmore, M. R. P. et al. Colony-stimulating factor 1 receptor signaling is necessary for microglia viability, unmasking a microglia progenitor cell in the adult brain. Neuron 82, 380–397 (2014).
Hwang, D. et al. CSF-1 maintains pathogenic but not homeostatic myeloid cells in the central nervous system during autoimmune neuroinflammation. Proc. Natl Acad. Sci. 119, e2111804119 (2022).
Hagan, N. et al. CSF1R signaling is a regulator of pathogenesis in progressive MS. Cell Death Dis. 11, 904 (2020).
Nissen, J. C., Thompson, K. K., West, B. L. & Tsirka, S. E. Csf1R inhibition attenuates experimental autoimmune encephalomyelitis and promotes recovery. Exp. Neurol. 307, 24–36 (2018).
Wlodarczyk, A. et al. CSF1R stimulation promotes increased neuroprotection by CD11c+ microglia in EAE. Front. Cell. Neurosci. 12, 523 (2019).
Korotkevich, G. et al. Fast gene set enrichment analysis. Preprint at https://doi.org/10.1101/060012 (2021).
Cold Spring Harbor Protocols. Red blood cell lysis buffer. Cold Spring Harbor Protocols 2006:db.rec390 (2006).
Hahn, O. et al. CoolMPS for robust sequencing of single-nuclear RNAs captured by droplet-based method. Nucleic Acids Res. 49, e11 (2021).
Liu, S., Lockhart, J. R., Fontenard, S., Berlett, M. & Ryan, T. M. Mapping the chromosomal insertion site of the GFP transgene of UBC-GFP mice to the MHC locus. J. Immunol. 204, 1982–1987 (2020).
Stuart, T. et al. Comprehensive integration of single-cell data. Cell 177, 1888–1902.e21 (2019).
Hafemeister, C. & Satija, R. Normalization and variance stabilization of single-cell RNA-seq data using regularized negative binomial regression. Genome Biol. 20, 296 (2019).
Russ, D. E. et al. A harmonized atlas of mouse spinal cord cell types and their spatial organization. Nat. Commun. 12, 5722 (2021).
Blum, J. A. et al. Single-cell transcriptomic analysis of the adult mouse spinal cord reveals molecular diversity of autonomic and skeletal motor neurons. Nat. Neurosci. 24, 572–583 (2021).
Zeisel, A. et al. Molecular architecture of the mouse nervous system. Cell 174, 999–1014.e22 (2018).
Marques, S. et al. Oligodendrocyte heterogeneity in the mouse juvenile and adult central nervous system. Science 352, 1326–1329 (2016).
Finak, G. et al. MAST: a flexible statistical framework for assessing transcriptional changes and characterizing heterogeneity in single-cell RNA sequencing data. Genome Biol. 16, 278 (2015).
Speir, M. L. et al. UCSC Cell Browser: visualize your single-cell data. Bioinformatics 37, 4578–4580 (2021).
Hao, Y. et al. Integrated analysis of multimodal single-cell data. Cell 184, 3573–3587.e29 (2021).
Durinck, S., Spellman, P. T., Birney, E. & Huber, W. Mapping identifiers for the integration of genomic datasets with the R/Bioconductor package biomaRt. Nat. Protoc. 4, 1184–1191 (2009).
Zheng, G. X. Y. et al. Massively parallel digital transcriptional profiling of single cells. Nat. Commun. 8, 14049 (2017).
Acknowledgements
We would like to thank all members of the Wernig laboratory, G. Wang, C. Bennett, B. Ajami, J. Blum and N. Or-Geva for helpful discussions throughout the project; A. Lang and M. Vangipuram for administrative support; and the FACS core at the Institute for Stem Cell Biology and Regenerative Medicine and Stanford Neuroscience Microscopy Service (supported by the National Institutes of Health, NS069375) for technical support. We thank the scientists at Plexxikon Inc. for their generous gift of PLX5622. Funding was provided by a Howard Hughes Medical Institute Faculty Scholar Award (to M.W.); the German Research Foundation (Deutsche Forschungsgemeinschaft; MA 8492/1-1 to M.M.D.M.); the National Multiple Sclerosis Society and the American Brain Foundation (FAN-2207-39823 to D.W.); and the New York Stem Cell Foundation Druckenmiller Fellowship (NYSCF-D-F74 to Y.Y.).
Author information
Authors and Affiliations
Contributions
M.M.D.M., M.W., L.S., R.D. and C.T. designed and conceived the study. M.M.D.M., A.N., Y.S., A.S., D.W., Y.Y., R.D. and C.T. worked with the study animals. M.M.D.M., M.A., A.F., O.H. and T.W.C. performed the NucSeq experiments and analysis. D.W., A.S., M.M.D.M. and A.N. performed immunostaining, microscopy and image analysis. M.M.D.M., M.W., O.H., A.N., D.W. and L.S. performed data analysis and interpretation. M.M.D.M. and M.W. drafted and edited the original manuscript; all authors reviewed, revised and approved the final version of the paper.
Corresponding author
Ethics declarations
Competing interests
The authors declare that they have no conflict of interest related to this study. PLX5622 was provided by Plexxikon Inc. under a material transfer agreement between Stanford University and Plexxikon Inc.
Peer review
Peer review information
Nature Neuroscience thanks Jessica Williams and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data
Extended Data Fig. 1 Characterization of chronic EAE via spinal cord NucSeq.
a, Experimental design and groups of the NucSeq main cohort. b, Absolute number of nuclei per library. c, Number of features and counts per library. d, Annotated UMAP plot of 68,593 nuclei integrated from all 8 libraries. Clustering based on the first 20 principal components and a resolution of 0.6. e, Annotation and canonical marker genes of clusters. f, Significant DEGs (padj < 0.05, based on MAST) for main cell types between the Control and EAE condition. The top 15 intersections are illustrated.
Extended Data Fig. 2 Flowcytometric analysis of murine spinal cord nuclei reveals an increase in the NeuN- Olig2- fraction in EAE.
a, Gating strategy of fluorescence activated nuclei sorting and representative gates for each animal. b, Quantification of gated populations. Two-sided Welch t-test. Total n = 14 animals. Box plot elements represent median (center line), first and third quartiles (lower and upper hinges) and smallest/highest value with at most 1.5*IQR (inter-quartile range) from the hinge (whiskers).
Extended Data Fig. 3 EAE alters the transcriptional program in both astrocytes and myeloid cells.
a, GSEA-derived top 20 significantly enriched pathways (Control vs EAE) of the Gene Ontology biological process term for the myeloid and astrocyte clusters (p values based on GSEA). b, c, Representative genes of selected enriched pathways.
Extended Data Fig. 4 EAE induces a proinflammatory gene set enrichment in oligodendrocytes.
a, GSEA-derived top 16 significantly enriched pathways (Control vs EAE) of the Gene Ontology biological process term for oligodendrocyte lineage clusters (p values based on GSEA). b, Representative genes of selected enriched pathways. c, Immunofluorescent stain for Olig2 and complement component 4 (C4) in the spinal cord white matter. White arrows demonstrate spatial association between C4 and oligodendrocytes. Stain was evaluated in three independent EAE and three independent Control spinal cords. Scale bar is 20 µm.
Extended Data Fig. 5 EAE-specific transcriptional changes are not induced by CFA injection without MOG-immunization.
a, Experimental design of the control cohort. Animals with Complete Freund’s Adjuvant (CFA) were sacrificed and analyzed 25 days after injection. Animals injected with phosphate buffered saline (PBS) were age matched to the animals of the main cohort at the timepoint of analysis. b, Absolute number of nuclei per animal sample. c, Number of features and counts per animal sample. d, Annotation and canonical marker genes of cell clusters. Clustering based on the first 20 principal components and a resolution of 0.6. e, UMAP plot of 17,127 nuclei integrated from the murine control cohort. f, Distribution of main cell clusters of the PBS and CFA groups. g, Fraction of myeloid nuclei. Dots represent individual animals (n = 8 animals total). Two-sided Welch t-test. h, Absolute number and top five intersections of significant DEGs (padj < 0.05) between PBS and CFA group. i, DEGs of the myeloid cluster between PBS and CFA group. Significant (padj < 0.05) genes are color coded if >1.5 fold up (red) or down (blue) regulated. j, Expression of EAE Scores in Control and EAE animals of the main cohort (left panels), and PBS and CFA animals of the control cohort (right panels). Box plot elements represent median (center line), first and third quartiles (lower and upper hinges) and smallest/highest value with at most 1.5*IQR (inter-quartile range) from the hinge (whiskers).
Extended Data Fig. 6 Integration of NucSeq data derived from MS patients.
a, Absolute number of nuclei per original human dataset. b, Number of features and counts per original human dataset. c, UMAP plot of 132,425 nuclei integrated from the original human brain datasets. Lesion-dependent conditions of the original datasets are shown (A_, Absinta; J_, Jäkel; S_, Schirmer; A, active; A/CA, acute/chronic active; CA, chronic active; CI, chronic inactive; remyel, remyelinating; wm, white matter), which were combined into 4 groups: control, active/chronic active (A/CA), chronic inactive (CI), and MS white matter (WM). d, Annotation and canonical marker genes of brain clusters. Clustering based on the first 12 principal components and a resolution of 0.4. e, UMAP plot of 5419 nuclei integrated from the original individual libraries derived from human spinal cord samples. f, Annotation and canonical marker genes of spinal cord clusters. Clustering based on the first 12 principal components and a resolution of 0.4.
Extended Data Fig. 7 Active MS brain lesions demonstrate myeloid cell expansion and activation.
a, Relative distribution of main cluster groups between different original datasets and control/MS conditions. b, Relative proportion of the myeloid and lymphocyte cell clusters. Dots represent proportions per individual original sample (n = 54 total samples derived from brain and 6 from spinal cord). A/CA, active / chronic active lesion; CI, inactive lesion; WM, white matter. Two-sided Welch t-test. c, Number of significant DEGs (padj < 0.05, based on MAST) per main cell type between the Control and MS condition. Only cell types with significant DEGs are shown. d, Expression of the human orthologs of the cell type specific EAE Scores in the myeloid, oligodendrocyte, and astrocyte subclusters (p values based on two-sided Mann–Whitney U test with Holm method adjustment). e, UMAP plots of the subclusters of the immune cells of the human brain (left panel, 7822 nuclei) and spinal cord (right panel, 794 nuclei) datasets. f, Annotation and marker genes of brain (upper panel) and spinal cord (lower panel) immune subclusters. g, Distribution of cell populations between Control and MS samples for the human brain (left panel) and spinal cord (right panel). Notably, the category ‘Remyel’ consists of only 2 nuclei. h, Expression of the human ortholog of the EAE Myeloid Score in myeloid subclusters of the human brain (upper panel) and spinal cord (lower panel) datasets. i, Expression of IQGAP1 in the brain myeloid cluster. Box plot elements represent median (center line), first and third quartiles (lower and upper hinges) and smallest/highest value with at most 1.5*IQR (inter-quartile range) from the hinge (whiskers).
Extended Data Fig. 8 BMT-induced morphological changes in myeloid cells are present in both spinal cord and brain.
a, Representative FACS gating strategy to quantify GFP chimerism in peripheral blood. b, Representative images of GFP+ chimerism in brain myeloid cells. Scale bars are 100 µm. c, Quantification of GFP positive and negative Iba1+ cells. Dots represent the average of multiple measurements per animal. A total of 17 animals or 26 animals were analyzed for spinal cord or brain regions, respectively. d, Pearson correlation coefficient (r) and scatter plot demonstrate the correlation between brain and spinal cord GFP+ chimerism. Regression line based on linear model. Two-sided test. e, A gene signature for CDMCs was based on previously published 1296 DEGs between CDMCs and endogenous microglia25 and a score was calculated with the VISION pipeline. Expression in the myeloid cluster is demonstrated. f, The fractions of nuclei in the myeloid cluster with a CDMC score expression of >−0.05 are shown (total n = 8 sequencing libraries derived from 14 spinal cords). g, Quantification of maximum branch length of ramified myeloid cells of the spinal cord in different conditions. Each dot represents one cell (total n = 210), dot colors represent different animals (n = 12). Two-sided Mann–Whitney U test. h, Representative immunofluorescent images of ramified myeloid cells in the brain. Brightness/contrast was adjusted individually per image for morphological assessment due to differing Iba1 intensity between conditions. Scale bars are 20 µm. i, Morphological analysis of ramified myeloid cells of the cortex in different conditions. Each dot represents one cell (total n = 162), dot colors represent different animals (n = 15). Two-sided Mann–Whitney U test. Box plot elements represent median (center line), first and third quartiles (lower and upper hinges) and smallest/highest value with at most 1.5*IQR (inter-quartile range) from the hinge (whiskers).
Extended Data Fig. 9 Transcriptional changes in astrocytes and oligodendrocytes in the context of EAE and BMT.
a, Venn diagrams of DEGs (padj < 0.05, based on MAST with Bonferroni correction) between selected conditions for main oligodendrocyte and astrocyte clusters. b, Scatterplots show the relationship of significant DEGs (padj < 0.05 for both conditions, based on MAST with Bonferroni correction) of disease-associated (Control vs EAE) and treatment-associated (EAE vs EAE-BMT or EAE vs EAE-BMT-PLX) conditions of different oligodendrocyte and astrocyte clusters. Purple numbers represent the DEGs per quadrant. c, Gene set enrichment analysis based on differentially expressed genes between the EAE and EAE-BMT groups in selected clusters. Shown are the top six Gene Ontology (GO) Biological Process (BP) pathways based on GSEA-derived p values. Terms are abbreviated to fit the format. A complete list can be found in Supplementary Table 3. d, Correlation between white matter demyelination and clinical score. Regression lines based on linear models; r = Pearson correlation coefficient. Two-sided test. e, Expression of marker genes and EAE Astro Score in astrocyte subclusters. f, Expression profile of DEGs shared with genes of a multiple sclerosis genome wide association study (GWAS) for astrocyte subclusters. Expression values (average log2 fold-change) are normalized by columns. If no significant DEG was present for a cell type, the value 0 was applied.
Extended Data Fig. 10 Transcriptional changes in myeloid cells in the context of EAE and BMT.
a, Marker gene expression of different T cell populations. NK, natural killer cells. Teff, effector T cells. Treg, regulatory T cells. b, Expression of GFP in nuclei of the myeloid cluster. Ordered towards the front based on expression. c, Expression of selected microglia marker genes in clusters M.1, M.2, and M.5 in different conditions. Adjusted p values based on DEG analysis (MAST with Bonferroni correction), compare also Supplementary Table 16. d, Venn diagram showing significant DEGs (padj <0.05, based on MAST with Bonferroni correction) of the myeloid cluster between different conditions. e, Scatterplots show the relationship of significant DEGs (padj < 0.05 for both conditions, MAST with Bonferroni correction) of disease-associated (Control vs EAE) and treatment-associated (EAE vs EAE-BMT or EAE vs EAE-BMT-PLX) conditions of the myeloid cluster. Purple numbers represent the DEGs per quadrant. f, Expression profile of DEGs shared with genes of a multiple sclerosis genome wide association study (GWAS) for myeloid subclusters. Expression values (average log2 fold-change) are normalized by columns. If no significant DEG was present for a cell type, the value 0 was applied.
Supplementary information
Supplementary Information
Supplementary File 1
Supplementary Table 1
Marker genes of main cell clusters of the main EAE NucSeq dataset. Based on Seurat’s FindAllMarkers function (Wilcoxon Rank Sum test; adjusted P value based on Bonferroni correction). ‘pct.1’ represents the percentage of cells in the cluster of interest, in which the gene is expressed. ‘pct.2’ represents the percentage of cells in all other clusters, in which the gene is expressed.
Supplementary Table 2
Differentially expressed genes of main cell clusters of the main EAE NucSeq dataset. Calculation was based on the MAST algorithm with a log-fold-change threshold of 0.25, minimum detection fraction of 0.1 and Bonferroni correction of P values. ‘pct.1’ and ‘pct.2’ represent the percentages of cells in the cluster of interest (celltype) in condition 1 (cond1) and 2 (cond2), respectively, in which the gene is expressed.
Supplementary Table 3
Enriched pathways of main cell clusters of the main EAE NucSeq dataset (GSEA). Organized per pathway type (pathway.category), compared conditions (cond1 and cond2), and cell cluster (celltype).
Supplementary Table 4
Gene signatures and mouse orthologous GWAS genes used in this study (as indicated in the ‘category’ column).
Supplementary Table 5
Marker genes of main cell clusters of the control cohort NucSeq dataset. Based on Seurat’s FindAllMarkers function (Wilcoxon Rank Sum test; adjusted P value based on Bonferroni correction). ‘pct.1’ represents the percentage of cells in the cluster of interest, in which the gene is expressed. ‘pct.2’ represents the percentage of cells in all other clusters, in which the gene is expressed.
Supplementary Table 6
Differentially expressed genes of main cell clusters of the control cohort NucSeq dataset. Calculation was based on the MAST algorithm with a log-fold-change threshold of 0.25, minimum detection fraction of 0.1 and Bonferroni correction of P values. ‘pct.1’ and ‘pct.2’ represent the percentages of cells in the cluster of interest (celltype) in condition 1 (cond1) and 2 (cond2), respectively, in which the gene is expressed.
Supplementary Table 7
Marker genes of main cell clusters of the human brain NucSeq dataset. Based on Seurat’s FindAllMarkers function (Wilcoxon Rank Sum test; adjusted P value based on Bonferroni correction). ‘pct.1’ represents the percentage of cells in the cluster of interest, in which the gene is expressed. ‘pct.2’ represents the percentage of cells in all other clusters, in which the gene is expressed.
Supplementary Table 8
Marker genes of main cell clusters of the human spinal cord NucSeq dataset. Based on Seurat’s FindAllMarkers function (Wilcoxon Rank Sum test; adjusted P value based on Bonferroni correction). ‘pct.1’ represents the percentage of cells in the cluster of interest, in which the gene is expressed. ‘pct.2’ represents the percentage of cells in all other clusters, in which the gene is expressed.
Supplementary Table 9
Differentially expressed genes (control vs MS) of main cell clusters of the human brain NucSeq dataset. Calculation was based on the MAST algorithm with a log-fold-change threshold of 0.25, minimum detection fraction of 0.1 and Bonferroni correction of P values. ‘pct.1’ and ‘pct.2’ represent the percentages of cells in the cluster of interest (celltype) in condition 1 (cond1) and 2 (cond2), respectively, in which the gene is expressed.
Supplementary Table 10
Differentially expressed genes (control vs MS) of main cell clusters of the human spinal cord NucSeq dataset. Calculation was based on the MAST algorithm with a log-fold-change threshold of 0.25, minimum detection fraction of 0.1 and Bonferroni correction of P values. ‘pct.1’ and ‘pct.2’ represent the percentages of cells in the cluster of interest (celltype) in condition 1 (cond1) and 2 (cond2), respectively, in which the gene is expressed.
Supplementary Table 11
Marker genes of immune subclusters of the human brain NucSeq dataset. Based on Seurat’s FindAllMarkers function (Wilcoxon Rank Sum test; adjusted P value based on Bonferroni correction). ‘pct.1’ represents the percentage of cells in the cluster of interest, in which the gene is expressed. ‘pct.2’ represents the percentage of cells in all other clusters, in which the gene is expressed.
Supplementary Table 12
Marker genes of immune subclusters of the human spinal cord NucSeq dataset. Based on Seurat’s FindAllMarkers function (Wilcoxon Rank Sum test; adjusted P value based on Bonferroni correction). ‘pct.1’ represents the percentage of cells in the cluster of interest, in which the gene is expressed. ‘pct.2’ represents the percentage of cells in all other clusters, in which the gene is expressed.
Supplementary Table 13
Marker genes of astrocyte subcluster of the EAE NucSeq dataset. Based on Seurat’s FindAllMarkers function (Wilcoxon Rank Sum test; adjusted P value based on Bonferroni correction). ‘pct.1’ represents the percentage of cells in the cluster of interest, in which the gene is expressed. ‘pct.2’ represents the percentage of cells in all other clusters, in which the gene is expressed.
Supplementary Table 14
Marker genes of immune subclusters of the EAE NucSeq dataset. Based on Seurat’s FindAllMarkers function (Wilcoxon Rank Sum test; adjusted P value based on Bonferroni correction). ‘pct.1’ represents the percentage of cells in the cluster of interest, in which the gene is expressed. ‘pct.2’ represents the percentage of cells in all other clusters, in which the gene is expressed.
Supplementary Table 15
Marker genes of TC&NK subclusters of the EAE NucSeq dataset. Based on Seurat’s FindAllMarkers function (Wilcoxon Rank Sum test; adjusted P value based on Bonferroni correction). ‘pct.1’ represents the percentage of cells in the cluster of interest, in which the gene is expressed. ‘pct.2’ represents the percentage of cells in all other clusters, in which the gene is expressed.
Supplementary Table 16
Differentially expressed genes of immune cell subclusters of the EAE NucSeq dataset. Calculation was based on the MAST algorithm with a log-fold-change threshold of 0.25, minimum detection fraction of 0.1 and Bonferroni correction of P values. ‘pct.1’ and ‘pct.2’ represent the percentages of cells in the cluster of interest (celltype) in condition 1 (cond1) and 2 (cond2), respectively, in which the gene is expressed.
Supplementary Table 17
Enriched pathways of myeloid subcluster marker genes of the EAE NucSeq dataset (GSEA). Organized per pathway type (pathway.category) and cell subcluster (cluster).
Source data
Source Data Fig. 1
Statistical Source Data
Source Data Fig. 2
Statistical Source Data
Source Data Fig. 3
Statistical Source Data
Source Data Fig. 4
Statistical Source Data
Source Data Fig. 5
Statistical Source Data
Source Data Fig. 6
Statistical Source Data
Source Data Extended Data Fig./Table 1
Statistical Source Data
Source Data Extended Data Fig./Table 2
Statistical Source Data
Source Data Extended Data Fig./Table 3
Statistical Source Data
Source Data Extended Data Fig./Table 4
Statistical Source Data
Source Data Extended Data Fig./Table 5
Statistical Source Data
Source Data Extended Data Fig./Table 6
Statistical Source Data
Source Data Extended Data Fig./Table 7
Statistical Source Data
Source Data Extended Data Fig./Table 8
Statistical Source Data
Source Data Extended Data Fig./Table 9
Statistical Source Data
Source Data Extended Data Fig./Table 10
Statistical Source Data
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Mader, M.MD., Napole, A., Wu, D. et al. Myeloid cell replacement is neuroprotective in chronic experimental autoimmune encephalomyelitis. Nat Neurosci 27, 901–912 (2024). https://doi.org/10.1038/s41593-024-01609-3
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41593-024-01609-3