Immunology in the Era of Single-Cell Technologies.

Immune cells are characterized by diversity, specificity, plasticity, and adaptability-properties that enable them to contribute to homeostasis and respond specifically and dynamically to the many threats encountered by the body. Single-cell technologies, including the assessment of transcriptomics, genomics, and proteomics at the level of individual cells, are ideally suited to studying these properties of immune cells. In this review we discuss the benefits of adopting single-cell approaches in studying underappreciated qualities of immune cells and highlight examples where these technologies have been critical to advancing our understanding of the immune system in health and disease.

Integrative spatial analysis reveals a multi-layered organization of glioblastoma.

Glioma contains malignant cells in diverse states. Here, we combine spatial transcriptomics, spatial proteomics, and computational approaches to define glioma cellular states and uncover their organization. We find three prominent modes of organization. First, gliomas are composed of small local environments, each typically enriched with one major cellular state. Second, specific pairs of states preferentially reside in proximity across multiple scales. This pairing of states is consistent across tumors. Third, these pairwise interactions collectively define a global architecture composed of five layers. Hypoxia appears to drive the layers, as it is associated with a long-range organization that includes all cancer cell states. Accordingly, tumor regions distant from any hypoxic/necrotic foci and tumors that lack hypoxia such as low-grade IDH-mutant glioma are less organized. In summary, we provide a conceptual framework for the organization of cellular states in glioma, highlighting hypoxia as a long-range tissue organizer.

The cycling and aging mouse female reproductive tract at single-cell resolution.

The female reproductive tract (FRT) undergoes extensive remodeling during reproductive cycling. This recurrent remodeling and how it shapes organ-specific aging remains poorly explored. Using single-cell and spatial transcriptomics, we systematically characterized morphological and gene expression changes occurring in ovary, oviduct, uterus, cervix, and vagina at each phase of the mouse estrous cycle, during decidualization, and into aging. These analyses reveal that fibroblasts play central-and highly organ-specific-roles in FRT remodeling by orchestrating extracellular matrix (ECM) reorganization and inflammation. Our results suggest a model wherein recurrent FRT remodeling over reproductive lifespan drives the gradual, age-related development of fibrosis and chronic inflammation. This hypothesis was directly tested using chemical ablation of cycling, which reduced fibrotic accumulation during aging. Our atlas provides extensive detail into how estrus, pregnancy, and aging shape the organs of the female reproductive tract and reveals the unexpected cost of the recurrent remodeling required for reproduction.

Open-ST: High-resolution spatial transcriptomics in 3D.

Spatial transcriptomics (ST) methods unlock molecular mechanisms underlying tissue development, homeostasis, or disease. However, there is a need for easy-to-use, high-resolution, cost-efficient, and 3D-scalable methods. Here, we report Open-ST, a sequencing-based, open-source experimental and computational resource to address these challenges and to study the molecular organization of tissues in 2D and 3D. In mouse brain, Open-ST captured transcripts at subcellular resolution and reconstructed cell types. In primary head-and-neck tumors and patient-matched healthy/metastatic lymph nodes, Open-ST captured the diversity of immune, stromal, and tumor populations in space, validated by imaging-based ST. Distinct cell states were organized around cell-cell communication hotspots in the tumor but not the metastasis. Strikingly, the 3D reconstruction and multimodal analysis of the metastatic lymph node revealed spatially contiguous structures not visible in 2D and potential biomarkers precisely at the 3D tumor/lymph node boundary. All protocols and software are available at https://rajewsky-lab.github.io/openst.

Mapping the cellular biogeography of human bone marrow niches using single-cell transcriptomics and proteomic imaging.

Non-hematopoietic cells are essential contributors to hematopoiesis. However, heterogeneity and spatial organization of these cells in human bone marrow remain largely uncharacterized. We used single-cell RNA sequencing (scRNA-seq) to profile 29,325 non-hematopoietic cells and discovered nine transcriptionally distinct subtypes. We simultaneously profiled 53,417 hematopoietic cells and predicted their interactions with non-hematopoietic subsets. We employed co-detection by indexing (CODEX) to spatially profile over 1.2 million cells. We integrated scRNA-seq and CODEX data to link predicted cellular signaling with spatial proximity. Our analysis revealed a hyperoxygenated arterio-endosteal neighborhood for early myelopoiesis, and an adipocytic localization for early hematopoietic stem and progenitor cells (HSPCs). We used our CODEX atlas to annotate new images and uncovered mesenchymal stromal cell (MSC) expansion and spatial neighborhoods co-enriched for leukemic blasts and MSCs in acute myeloid leukemia (AML) patient samples. This spatially resolved, multiomic atlas of human bone marrow provides a reference for investigation of cellular interactions that drive hematopoiesis.

Unraveling cellular complexity with transient adapters in highly multiplexed super-resolution imaging.

Mapping the intricate spatial relationships between the many different molecules inside a cell is essential to understanding cellular functions in all their complexity. Super-resolution fluorescence microscopy offers the required spatial resolution but struggles to reveal more than four different targets simultaneously. Exchanging labels in subsequent imaging rounds for multiplexed imaging extends this number but is limited by its low throughput. Here, we present a method for rapid multiplexed super-resolution microscopy that can, in principle, be applied to a nearly unlimited number of molecular targets by leveraging fluorogenic labeling in conjunction with transient adapter-mediated switching for high-throughput DNA-PAINT (FLASH-PAINT). We demonstrate the versatility of FLASH-PAINT with four applications: mapping nine proteins in a single mammalian cell, elucidating the functional organization of primary cilia by nine-target imaging, revealing the changes in proximity of thirteen different targets in unperturbed and dissociated Golgi stacks, and investigating and quantifying inter-organelle contacts at 3D super-resolution.

Spatial proteomics in neurons at single-protein resolution.

To understand biological processes, it is necessary to reveal the molecular heterogeneity of cells by gaining access to the location and interaction of all biomolecules. Significant advances were achieved by super-resolution microscopy, but such methods are still far from reaching the multiplexing capacity of proteomics. Here, we introduce secondary label-based unlimited multiplexed DNA-PAINT (SUM-PAINT), a high-throughput imaging method that is capable of achieving virtually unlimited multiplexing at better than 15 nm resolution. Using SUM-PAINT, we generated 30-plex single-molecule resolved datasets in neurons and adapted omics-inspired analysis for data exploration. This allowed us to reveal the complexity of synaptic heterogeneity, leading to the discovery of a distinct synapse type. We not only provide a resource for researchers, but also an integrated acquisition and analysis workflow for comprehensive spatial proteomics at single-protein resolution.

Glioblastoma evolution and heterogeneity from a 3D whole-tumor perspective.

Treatment failure for the lethal brain tumor glioblastoma (GBM) is attributed to intratumoral heterogeneity and tumor evolution. We utilized 3D neuronavigation during surgical resection to acquire samples representing the whole tumor mapped by 3D spatial coordinates. Integrative tissue and single-cell analysis revealed sources of genomic, epigenomic, and microenvironmental intratumoral heterogeneity and their spatial patterning. By distinguishing tumor-wide molecular features from those with regional specificity, we inferred GBM evolutionary trajectories from neurodevelopmental lineage origins and initiating events such as chromothripsis to emergence of genetic subclones and spatially restricted activation of differential tumor and microenvironmental programs in the core, periphery, and contrast-enhancing regions. Our work depicts GBM evolution and heterogeneity from a 3D whole-tumor perspective, highlights potential therapeutic targets that might circumvent heterogeneity-related failures, and establishes an interactive platform enabling 360° visualization and analysis of 3D spatial patterns for user-selected genes, programs, and other features across whole GBM tumors.

3D reconstruction of a gastrulating human embryo.

Progress in understanding early human development has been impeded by the scarcity of reference datasets from natural embryos, particularly those with spatial information during crucial stages like gastrulation. We conducted high-resolution spatial transcriptomics profiling on 38,562 spots from 62 transverse sections of an intact Carnegie stage (CS) 8 human embryo. From this spatial transcriptomic dataset, we constructed a 3D model of the CS8 embryo, in which a range of cell subtypes are identified, based on gene expression patterns and positional register, along the anterior-posterior, medial-lateral, and dorsal-ventral axis in the embryo. We further characterized the lineage trajectories of embryonic and extra-embryonic tissues and associated regulons and the regionalization of signaling centers and signaling activities that underpin lineage progression and tissue patterning during gastrulation. Collectively, the findings of this study provide insights into gastrulation and post-gastrulation development of the human embryo.

Cellular architecture of evolving neuroinflammatory lesions and multiple sclerosis pathology.

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.

The FXR1 network acts as a signaling scaffold for actomyosin remodeling.

It is currently not known whether mRNAs fulfill structural roles in the cytoplasm. Here, we report the fragile X-related protein 1 (FXR1) network, an mRNA-protein (mRNP) network present throughout the cytoplasm, formed by FXR1-mediated packaging of exceptionally long mRNAs. These mRNAs serve as an underlying condensate scaffold and concentrate FXR1 molecules. The FXR1 network contains multiple protein binding sites and functions as a signaling scaffold for interacting proteins. We show that it is necessary for RhoA signaling-induced actomyosin reorganization to provide spatial proximity between kinases and their substrates. Point mutations in FXR1, found in its homolog FMR1, where they cause fragile X syndrome, disrupt the network. FXR1 network disruption prevents actomyosin remodeling-an essential and ubiquitous process for the regulation of cell shape, migration, and synaptic function. Our findings uncover a structural role for cytoplasmic mRNA and show how the FXR1 RNA-binding protein as part of the FXR1 network acts as an organizer of signaling reactions.

Mapping the transcriptome: Realizing the full potential of spatial data analysis.

RNA sequencing in situ allows for whole-transcriptome characterization at high resolution, while retaining spatial information. These data present an analytical challenge for bioinformatics-how to leverage spatial information effectively? Properties of data with a spatial dimension require special handling, which necessitate a different set of statistical and inferential considerations when compared to non-spatial data. The geographical sciences primarily use spatial data and have developed methods to analye them. Here we discuss the challenges associated with spatial analysis and examine how we can take advantage of practice from the geographical sciences to realize the full potential of spatial information in transcriptomic datasets.

Multiplexed 3D atlas of state transitions and immune interaction in colorectal cancer.

Advanced solid cancers are complex assemblies of tumor, immune, and stromal cells characterized by high intratumoral variation. We use highly multiplexed tissue imaging, 3D reconstruction, spatial statistics, and machine learning to identify cell types and states underlying morphological features of known diagnostic and prognostic significance in colorectal cancer. Quantitation of these features in high-plex marker space reveals recurrent transitions from one tumor morphology to the next, some of which are coincident with long-range gradients in the expression of oncogenes and epigenetic regulators. At the tumor invasive margin, where tumor, normal, and immune cells compete, T cell suppression involves multiple cell types and 3D imaging shows that seemingly localized 2D features such as tertiary lymphoid structures are commonly interconnected and have graded molecular properties. Thus, while cancer genetics emphasizes the importance of discrete changes in tumor state, whole-specimen imaging reveals large-scale morphological and molecular gradients analogous to those in developing tissues.

Single-cell spatial transcriptome reveals cell-type organization in the macaque cortex.

Elucidating the cellular organization of the cerebral cortex is critical for understanding brain structure and function. Using large-scale single-nucleus RNA sequencing and spatial transcriptomic analysis of 143 macaque cortical regions, we obtained a comprehensive atlas of 264 transcriptome-defined cortical cell types and mapped their spatial distribution across the entire cortex. We characterized the cortical layer and region preferences of glutamatergic, GABAergic, and non-neuronal cell types, as well as regional differences in cell-type composition and neighborhood complexity. Notably, we discovered a relationship between the regional distribution of various cell types and the region's hierarchical level in the visual and somatosensory systems. Cross-species comparison of transcriptomic data from human, macaque, and mouse cortices further revealed primate-specific cell types that are enriched in layer 4, with their marker genes expressed in a region-dependent manner. Our data provide a cellular and molecular basis for understanding the evolution, development, aging, and pathogenesis of the primate brain.

Molecular and spatial signatures of mouse brain aging at single-cell resolution.

The diversity and complex organization of cells in the brain have hindered systematic characterization of age-related changes in its cellular and molecular architecture, limiting our ability to understand the mechanisms underlying its functional decline during aging. Here, we generated a high-resolution cell atlas of brain aging within the frontal cortex and striatum using spatially resolved single-cell transcriptomics and quantified changes in gene expression and spatial organization of major cell types in these regions over the mouse lifespan. We observed substantially more pronounced changes in cell state, gene expression, and spatial organization of non-neuronal cells over neurons. Our data revealed molecular and spatial signatures of glial and immune cell activation during aging, particularly enriched in the subcortical white matter, and identified both similarities and notable differences in cell-activation patterns induced by aging and systemic inflammatory challenge. These results provide critical insights into age-related decline and inflammation in the brain.

Multimodal charting of molecular and functional cell states via in situ electro-sequencing.

Paired mapping of single-cell gene expression and electrophysiology is essential to understand gene-to-function relationships in electrogenic tissues. Here, we developed in situ electro-sequencing (electro-seq) that combines flexible bioelectronics with in situ RNA sequencing to stably map millisecond-timescale electrical activity and profile single-cell gene expression from the same cells across intact biological networks, including cardiac and neural patches. When applied to human-induced pluripotent stem-cell-derived cardiomyocyte patches, in situ electro-seq enabled multimodal in situ analysis of cardiomyocyte electrophysiology and gene expression at the cellular level, jointly defining cell states and developmental trajectories. Using machine-learning-based cross-modal analysis, in situ electro-seq identified gene-to-electrophysiology relationships throughout cardiomyocyte development and accurately reconstructed the evolution of gene expression profiles based on long-term stable electrical measurements. In situ electro-seq could be applicable to create spatiotemporal multimodal maps in electrogenic tissues, potentiating the discovery of cell types and gene programs responsible for electrophysiological function and dysfunction.

Dynamic mapping of proteome trafficking within and between living cells by TransitID.

The ability to map trafficking for thousands of endogenous proteins at once in living cells would reveal biology currently invisible to both microscopy and mass spectrometry. Here, we report TransitID, a method for unbiased mapping of endogenous proteome trafficking with nanometer spatial resolution in living cells. Two proximity labeling (PL) enzymes, TurboID and APEX, are targeted to source and destination compartments, and PL with each enzyme is performed in tandem via sequential addition of their small-molecule substrates. Mass spectrometry identifies the proteins tagged by both enzymes. Using TransitID, we mapped proteome trafficking between cytosol and mitochondria, cytosol and nucleus, and nucleolus and stress granules (SGs), uncovering a role for SGs in protecting the transcription factor JUN from oxidative stress. TransitID also identifies proteins that signal intercellularly between macrophages and cancer cells. TransitID offers a powerful approach for distinguishing protein populations based on compartment or cell type of origin.

Spatiotemporal insight into early pregnancy governed by immune-featured stromal cells.

Endometrial decidualization connecting embryo implantation and placentation is transient but essential for successful pregnancy, which, however, is not systematically investigated. Here, we use a scStereo-seq technology to spatially visualize and define the dynamic functional decidual hubs assembled by distinct immune, endothelial, trophoblast, and decidual stromal cells (DSCs) in early pregnant mice. We unravel the DSC transdifferentiation trajectory and surprisingly discover a dual-featured type of immune-featured DSCs (iDSCs). We find that immature DSCs attract immune cells and induce decidual angiogenesis at the mesenchymal-epithelial transition hub during decidualization initiation. iDSCs enable immune cell recruitment and suppression, govern vascularization, and promote cytolysis at immune cell assembling and vascular hubs, respectively, to establish decidual homeostasis at a later stage. Interestingly, dysfunctional and spatially disordered iDSCs cause abnormal accumulation of immune cells in the vascular hub, which disrupts decidual hub specification and eventually leads to pregnancy complications in DBA/2-mated CBA/J mice.

Molecular cartography uncovers evolutionary and microenvironmental dynamics in sporadic colorectal tumors.

Colorectal cancer exhibits dynamic cellular and genetic heterogeneity during progression from precursor lesions toward malignancy. Analysis of spatial multi-omic data from 31 human colorectal specimens enabled phylogeographic mapping of tumor evolution that revealed individualized progression trajectories and accompanying microenvironmental and clonal alterations. Phylogeographic mapping ordered genetic events, classified tumors by their evolutionary dynamics, and placed clonal regions along global pseudotemporal progression trajectories encompassing the chromosomal instability (CIN+) and hypermutated (HM) pathways. Integrated single-cell and spatial transcriptomic data revealed recurring epithelial programs and infiltrating immune states along progression pseudotime. We discovered an immune exclusion signature (IEX), consisting of extracellular matrix regulators DDR1, TGFBI, PAK4, and DPEP1, that charts with CIN+ tumor progression, is associated with reduced cytotoxic cell infiltration, and shows prognostic value in independent cohorts. This spatial multi-omic atlas provides insights into colorectal tumor-microenvironment co-evolution, serving as a resource for stratification and targeted treatments.

Spatiotemporal transcriptome atlas reveals the regional specification of the developing human brain.

Different functional regions of brain are fundamental for basic neurophysiological activities. However, the regional specification remains largely unexplored during human brain development. Here, by combining spatial transcriptomics (scStereo-seq) and scRNA-seq, we built a spatiotemporal developmental atlas of multiple human brain regions from 6-23 gestational weeks (GWs). We discovered that, around GW8, radial glia (RG) cells have displayed regional heterogeneity and specific spatial distribution. Interestingly, we found that the regional heterogeneity of RG subtypes contributed to the subsequent neuronal specification. Specifically, two diencephalon-specific subtypes gave rise to glutamatergic and GABAergic neurons, whereas subtypes in ventral midbrain were associated with the dopaminergic neurons. Similar GABAergic neuronal subtypes were shared between neocortex and diencephalon. Additionally, we revealed that cell-cell interactions between oligodendrocyte precursor cells and GABAergic neurons influenced and promoted neuronal development coupled with regional specification. Altogether, this study provides comprehensive insights into the regional specification in the developing human brain.