Multiple sclerosis (MS) is a neurological disease characterized by multifocal lesions and smoldering pathology. Although single-cell analyses provided insights into cytopathology, evolving cellular processes underlying MS remain poorly understood. We investigated the cellular dynamics of MS by modeling temporal and regional rates of disease progression in mouse experimental autoimmune encephalomyelitis (EAE). By performing single-cell spatial expression profiling using in situ sequencing (ISS), we annotated disease neighborhoods and found centrifugal evolution of active lesions. We demonstrated that disease-associated (DA)-glia arise independently of lesions and are dynamically induced and resolved over the disease course. Single-cell spatial mapping of human archival MS spinal cords confirmed the differential distribution of homeostatic and DA-glia, enabled deconvolution of active and inactive lesions into sub-compartments, and identified new lesion areas. By establishing a spatial resource of mouse and human MS neuropathology at a single-cell resolution, our study unveils the intricate cellular dynamics underlying MS.

Tuberculosis (TB) and sarcoidosis are both granulomatous diseases. Here, we compared the immunological microenvironments of granulomas from TB and sarcoidosis patients using in situ sequencing (ISS) transcriptomic analysis and multiplexed immunolabeling of tissue sections. TB lesions consisted of large necrotic and cellular granulomas, whereas “multifocal” granulomas with macrophages or epitheloid cell core and a T-cell rim were observed in sarcoidosis samples. The necrotic core in TB lesions was surrounded by macrophages and encircled by a dense T-cell layer. Within the T-cell layer, compact B-cell aggregates were observed in most TB samples. These B-cell clusters were vascularized and could contain defined B-/T-cell and macrophage-rich areas. The ISS of 40–60 immune transcripts revealed the enriched expression of transcripts involved in homing or migration to lymph nodes, which formed networks at single-cell distances in lymphoid areas of the TB lesions. Instead, myeloid-annotated regions were enriched in CD68, CD14, ITGAM, ITGAX, and CD4 mRNA. CXCL8 and IL1B mRNA were observed in granulocytic areas in which M. tuberculosis was also detected. In line with ISS data indicating tertiary lymphoid structures, immune labeling of TB sections expressed markers of high endothelial venules, follicular dendritic cells, follicular helper T cells, and lymph-node homing receptors on T cells. Neither ISS nor immunolabeling showed evidence of tertiary lymphoid aggregates in sarcoidosis samples. Together, our finding suggests that despite their heterogeneity, the formation of tertiary immune structures is a common feature in granulomas from TB patients.

Background: Oligodendrocytes are glial cells that support and insulate axons in the central nervous system through the production of myelin. Oligodendrocytes arise throughout embryonic and early postnatal development from oligodendrocyte precursor cells (OPCs), and recent work demonstrated that they are a transcriptional heterogeneous cell population, but the regional and functional implications of this heterogeneity are less clear. Here, we apply in situ sequencing (ISS) to simultaneously probe the expression of 124 marker genes of distinct oligodendrocyte populations, providing comprehensive maps of the corpus callosum, cingulate, motor, and somatosensory cortex in the brain, as well as gray matter (GM) and white matter (WM) regions in the spinal cord, at postnatal (P10), juvenile (P20), and young adult (P60) stages. We systematically compare the abundances of these populations and investigate the neighboring preference of distinct oligodendrocyte populations.

Results: We observed that oligodendrocyte lineage progression is more advanced in the juvenile spinal cord compared to the brain, corroborating with previous studies. We found myelination still ongoing in the adult corpus callosum while it was more advanced in the cortex. Interestingly, we also observed a lateral-to-medial gradient of oligodendrocyte lineage progression in the juvenile cortex, which could be linked to arealization, as well as a deep-to-superficial gradient with mature oligodendrocytes preferentially accumulating in the deeper layers of the cortex. The ISS experiments also exposed differences in abundances and population dynamics over time between GM and WM regions in the brain and spinal cord, indicating regional differences within GM and WM, and we found that neighboring preferences of some oligodendroglia populations are altered from the juvenile to the adult CNS.

Conclusions: Overall, our ISS experiments reveal spatial heterogeneity of oligodendrocyte lineage progression in the brain and spinal cord and uncover differences in the timing of oligodendrocyte differentiation and myelination, which could be relevant to further investigate functional heterogeneity of oligodendroglia, especially in the context of injury or disease.

Mycobacterium tuberculosis (Mtb) bacilli are the causative agent of tuberculosis (TB), a major killer of mankind. Although it is widely accepted that local interactions between Mtb and the immune system in the tuberculous granuloma determine whether the outcome of infection is controlled or disseminated, these have been poorly studied due to methodological constraints. We have recently used a spatial transcriptomic technique, in situ sequencing (ISS), to define the spatial distribution of immune transcripts in TB mouse lungs. To further contribute to the understanding of the immune microenvironments of Mtb and their local diversity, we here present two complementary automated bacteria-guided analysis pipelines. These position 33 ISS-identified immune transcripts in relation to single bacteria and bacteria clusters. The analysis was applied on new ISS data from lung sections of Mtb-infected C57BL/6 and C3HeB/FeJ mice. In lungs from C57BL/6 mice early and late post infection, transcripts that define inflammatory macrophages were enriched at subcellular distances to bacteria, indicating the activation of infected macrophages. In contrast, expression patterns associated to antigen presentation were enriched in non-infected cells at 12 weeks post infection. T-cell transcripts were evenly distributed in the tissue. In Mtb-infected C3HeB/FeJ mice, transcripts characterizing activated macrophages localized in apposition to small bacteria clusters, but not in organized granulomas. Despite differences in the susceptibility to Mtb, the transcript patterns found around small bacteria clusters of C3HeB/FeJ and C57BL/6 mice were similar. Altogether, the presented tools allow us to characterize in depth the immune cell populations and their activation that interact with Mtb in the infected lung.

Mature oligodendrocytes (MOLs) show transcriptional heterogeneity, the functional consequences of which are unclear. MOL heterogeneity might correlate with the local environment or their interactions with different neuron types. Here, we show that distinct MOL populations have spatial preference in the mammalian central nervous system (CNS). We found that MOL type 2 (MOL2) is enriched in the spinal cord when compared to the brain, while MOL types 5 and 6 (MOL5/6) increase their contribution to the OL lineage with age in all analyzed regions. MOL2 and MOL5/6 also have distinct spatial preference in the spinal cord regions where motor and sensory tracts run. OL progenitor cells (OPCs) are not specified into distinct MOL populations during development, excluding a major contribution of OPC intrinsic mechanisms determining MOL heterogeneity. In disease, MOL2 and MOL5/6 present different susceptibility during the chronic phase following traumatic spinal cord injury. Our results demonstrate that the distinct MOL populations have different spatial preference and different responses to disease.

Objective- The Wnt/beta-catenin pathway orchestrates development of the blood-brain barrier, but the downstream mechanisms involved at different developmental windows and in different central nervous system (CNS) tissues have remained elusive. Approach and Results- Here, we create a new mouse model allowing spatiotemporal investigations of Wnt/beta-catenin signaling by induced overexpression of Axin1, an inhibitor of beta-catenin signaling, specifically in endothelial cells (Axin1(iEC)-(OE)). AOE (Axin1 overexpression) in Axin1(iEC)-(OE) mice at stages following the initial vascular invasion of the CNS did not impair angiogenesis but led to premature vascular regression followed by progressive dilation and inhibition of vascular maturation resulting in forebrain-specific hemorrhage 4 days post-AOE. Analysis of the temporal Wnt/beta-catenin driven CNS vascular development in zebrafish also suggested that Axin1(iEC)-(OE) led to CNS vascular regression and impaired maturation but not inhibition of ongoing angiogenesis within the CNS. Transcriptomic profiling of isolated, beta-catenin signaling-deficient endothelial cells during early blood-brain barrier-development (E11.5) revealed ECM (extracellular matrix) proteins as one of the most severely deregulated clusters. Among the 20 genes constituting the forebrain endothelial cell-specific response signature, 8 (Adamtsl2, Apod, Ctsw, Htra3, Pglyrp1, Spock2, Ttyh2, and Wfdc1) encoded bona fide ECM proteins. This specific beta-catenin-responsive ECM signature was also repressed in Axin1(iEC)-(OE) and endothelial cell-specific beta-catenin-knockout mice (Ctnnb1-KOiEC) during initial blood-brain barrier maturation (E14.5), consistent with an important role of Wnt/beta-catenin signaling in orchestrating the development of the forebrain vascular ECM. Conclusions- These results suggest a novel mechanism of establishing a CNS endothelium-specific ECM signature downstream of Wnt-beta-catenin that impact spatiotemporally on blood-brain barrier differentiation during forebrain vessel development.

Single-cell transcriptomics provides us with completely new insights into the molecular diversity of different cell types and the different states they can adopt. The technique generates inventories of cells that constitute the building blocks of multicellular organisms. However, since the method requires isolation of discrete cells, information about the original location within tissue is lost. Therefore, it is not possible to draw detailed cellular maps of tissue architecture and their positioning in relation to other cells. In order to better understand the cellular and tissue function of multicellular organisms, we need to map the cells within their physiological, morphological, and anatomical context and space. In this review, we will summarize and compare the different methods of in situ RNA analysis and the most recent developments leading to more comprehensive and highly multiplexed spatially resolved transcriptomic approaches. We will discuss their highlights and advantages as well as their limitations and challenges and give an outlook on promising future applications and directions both within basic research as well as clinical integration.

Midbrain dopamine (mDA) neurons constitute a heterogenous group of cells that have been intensely studied, not least because their degeneration causes major symptoms in Parkinson's disease. Understanding the diversity of mDA neurons - previously well characterized anatomically - requires a systematic molecular classification at the genome-wide gene expression level. Here, we use single cell RNA sequencing of isolated mouse neurons expressing the transcription factor Pitx3, a marker for mDA neurons. Analyses include cells isolated during development up until adulthood and the results are validated by histological characterization of newly identified markers. This identifies seven neuron subgroups divided in two major branches of developing Pitx3-expressing neurons. Five of them express dopaminergic markers, while two express glutamatergic and GABAergic markers, respectively. Analysis also indicate evolutionary conservation of diversity in humans. This comprehensive molecular characterization will provide a valuable resource for elucidating mDA neuron subgroup development and function in the mammalian brain.

Neural crest cells are embryonic progenitors that generate numerous cell types in vertebrates. With single-cell analysis, we show that mouse trunk neural crest cells become biased toward neuronal lineages when they delaminate from the neural tube, whereas cranial neural crest cells acquire ectomesenchyme potential dependent on activation of the transcription factor Twist1. The choices that neural crest cells make to become sensory, glial, autonomic, or mesenchymal cells can be formalized as a series of sequential binary decisions. Each branch of the decision tree involves initial coactivation of bipotential properties followed by gradual shifts toward commitment. Competing fate programs are coactivated before cells acquire fate-specific phenotypic traits. Determination of a specific fate is achieved by increased synchronization of relevant programs and concurrent repression of competing fate programs.