Spatial transcriptomics reveal heterogeneous cell‒cell interactions among brain regions in cuprizone model consistent with multiple sclerosis lesions.

The cuprizone (CPZ) model is widely used for modeling demyelination in multiple sclerosis (MS) and for testing potential remyelination therapies. To better understand the underlying pathology of the CPZ model and evaluate its translatability, we integrated single-cell and spatial transcriptomics (ST) to investigate spatial cellular and molecular interactions during de- and remyelination in multiple brain regions. ST revealed global demyelination and neuroinflammation in the brain beyond the corpus callosum (CC), with region-specific differences. We identified oligodendroglia and microglia as two major cell types with significant transcriptomic changes in the model. CPZ-associated subclusters of oligodendroglia (marker genes Arap2, Dock10, Tenm4, Pex5l and Dock1) and microglia (marker genes ApoE, Axl, Cd9 and Lpl) were mapped to the CC by ST. During remyelination, while mature oligodendrocytes (MOL) nearly reversed their phenotype back to the control state, microglia remained associated with the demyelination phenotype. Ligand‒receptor (LR) pairing analyses predicted growth factor and phagocytic pathway enrichment during demyelination, which is consistent with changes in MS lesions, and microglia were predicted to be the major sender cells. LR pairing also predicted a high likelihood of interaction between oligodendroglia and microglia, and a novel interaction between MOL and oligodendrocyte precursor cells (OPC), underscoring their roles during de- and remyelination. Finally, astrocytes in the CPZ model had the greatest preservation of disease-associated modules in MS lesions, while MOL, OPC, and microglia showed moderate to low preservation, which overall suggests that the CPZ model has moderate translatability to chronically active MS lesions.