Molecular pixelation: spatial proteomics of single cells by sequencing.

The spatial distribution of cell surface proteins governs vital processes of the immune system such as intercellular communication and mobility. However, fluorescence microscopy has limited scalability in the multiplexing and throughput needed to drive spatial proteomics discoveries at subcellular level. We present Molecular Pixelation (MPX), an optics-free, DNA sequence-based method for spatial proteomics of single cells using antibody-oligonucleotide conjugates (AOCs) and DNA-based, nanometer-sized molecular pixels. The relative locations of AOCs are inferred by sequentially associating them into local neighborhoods using the sequence-unique DNA pixels, forming >1,000 spatially connected zones per cell in 3D. For each single cell, DNA-sequencing reads are computationally arranged into spatial proteomics networks for 76 proteins. By studying immune cell dynamics using spatial statistics on graph representations of the data, we identify known and new patterns of spatial organization of proteins on chemokine-stimulated T cells, highlighting the potential of MPX in defining cell states by the spatial arrangement of proteins.

Semla: a versatile toolkit for spatially resolved transcriptomics analysis and visualization.

SUMMARY: Spatially resolved transcriptomics technologies generate gene expression data with retained positional information from a tissue section, often accompanied by a corresponding histological image. Computational tools should make it effortless to incorporate spatial information into data analyses and present analysis results in their histological context. Here, we present semla, an R package for processing, analysis, and visualization of spatially resolved transcriptomics data generated by the Visium platform, that includes interactive web applications for data exploration and tissue annotation. AVAILABILITY AND IMPLEMENTATION: The R package semla is available on GitHub (https://github.com/ludvigla/semla), under the MIT License, and deposited on Zenodo (https://doi.org/10.5281/zenodo.8321645). Documentation and tutorials with detailed descriptions of usage can be found at https://ludvigla.github.io/semla/.

Spatial multimodal analysis of transcriptomes and metabolomes in tissues.

We present a spatial omics approach that combines histology, mass spectrometry imaging and spatial transcriptomics to facilitate precise measurements of mRNA transcripts and low-molecular-weight metabolites across tissue regions. The workflow is compatible with commercially available Visium glass slides. We demonstrate the potential of our method using mouse and human brain samples in the context of dopamine and Parkinson's disease.

Expansion spatial transcriptomics.

Capture array-based spatial transcriptomics methods have been widely used to resolve gene expression in tissues; however, their spatial resolution is limited by the density of the array. Here we present expansion spatial transcriptomics to overcome this limitation by clearing and expanding tissue prior to capturing the entire polyadenylated transcriptome with an enhanced protocol. This approach enables us to achieve higher spatial resolution while retaining high library quality, which we demonstrate using mouse brain samples.

Profiling spatiotemporal gene expression of the developing human spinal cord and implications for ependymoma origin.

The spatiotemporal regulation of cell fate specification in the human developing spinal cord remains largely unknown. In this study, by performing integrated analysis of single-cell and spatial multi-omics data, we used 16 prenatal human samples to create a comprehensive developmental cell atlas of the spinal cord during post-conceptional weeks 5-12. This revealed how the cell fate commitment of neural progenitor cells and their spatial positioning are spatiotemporally regulated by specific gene sets. We identified unique events in human spinal cord development relative to rodents, including earlier quiescence of active neural stem cells, differential regulation of cell differentiation and distinct spatiotemporal genetic regulation of cell fate choices. In addition, by integrating our atlas with pediatric ependymomas data, we identified specific molecular signatures and lineage-specific genes of cancer stem cells during progression. Thus, we delineate spatiotemporal genetic regulation of human spinal cord development and leverage these data to gain disease insight.

High cardiomyocyte diversity in human early prenatal heart development.

Cardiomyocytes play key roles during cardiogenesis, but have poorly understood features, especially in prenatal stages. Here, we characterized human prenatal cardiomyocytes, 6.5-7 weeks post-conception, by integrating single-cell RNA sequencing, spatial transcriptomics, and ligand-receptor interaction information. Using a computational workflow developed to dissect cell type heterogeneity, localize cell types, and explore their molecular interactions, we identified eight types of developing cardiomyocyte, more than double compared to the ones identified in the Human Developmental Cell Atlas. These have high variability in cell cycle activity, mitochondrial content, and connexin gene expression, and are differentially distributed in the ventricles, including outflow tract, and atria, including sinoatrial node. Moreover, cardiomyocyte ligand-receptor crosstalk is mainly with non-cardiomyocyte cell types, encompassing cardiogenesis-related pathways. Thus, early prenatal human cardiomyocytes are highly heterogeneous and develop unique location-dependent properties, with complex ligand-receptor crosstalk. Further elucidation of their developmental dynamics may give rise to new therapies.

Spatially resolved transcriptomic profiling of degraded and challenging fresh frozen samples.

Spatially resolved transcriptomics has enabled precise genome-wide mRNA expression profiling within tissue sections. The performance of methods targeting the polyA tails of mRNA relies on the availability of specimens with high RNA quality. Moreover, the high cost of currently available spatial resolved transcriptomics assays requires a careful sample screening process to increase the chance of obtaining high-quality data. Indeed, the upfront analysis of RNA quality can show considerable variability due to sample handling, storage, and/or intrinsic factors. We present RNA-Rescue Spatial Transcriptomics (RRST), a workflow designed to improve mRNA recovery from fresh frozen specimens with moderate to low RNA quality. First, we provide a benchmark of RRST against the standard Visium spatial gene expression protocol on high RNA quality samples represented by mouse brain and prostate cancer samples. Then, we test the RRST protocol on tissue sections collected from five challenging tissue types, including human lung, colon, small intestine, pediatric brain tumor, and mouse bone/cartilage. In total, we analyze 52 tissue sections and demonstrate that RRST is a versatile, powerful, and reproducible protocol for fresh frozen specimens of different qualities and origins.

B cell expansion hinders the stroma-epithelium regenerative cross talk during mucosal healing.

Therapeutic promotion of intestinal regeneration holds great promise, but defining the cellular mechanisms that influence tissue regeneration remains an unmet challenge. To gain insight into the process of mucosal healing, we longitudinally examined the immune cell composition during intestinal damage and regeneration. B cells were the dominant cell type in the healing colon, and single-cell RNA sequencing (scRNA-seq) revealed expansion of an IFN-induced B cell subset during experimental mucosal healing that predominantly located in damaged areas and associated with colitis severity. B cell depletion accelerated recovery upon injury, decreased epithelial ulceration, and enhanced gene expression programs associated with tissue remodeling. scRNA-seq from the epithelial and stromal compartments combined with spatial transcriptomics and multiplex immunostaining showed that B cells decreased interactions between stromal and epithelial cells during mucosal healing. Activated B cells disrupted the epithelial-stromal cross talk required for organoid survival. Thus, B cell expansion during injury impairs epithelial-stromal cell interactions required for mucosal healing, with implications for the treatment of IBD.

Spatially resolved clonal copy number alterations in benign and malignant tissue.

Defining the transition from benign to malignant tissue is fundamental to improving early diagnosis of cancer1. Here we use a systematic approach to study spatial genome integrity in situ and describe previously unidentified clonal relationships. We used spatially resolved transcriptomics2 to infer spatial copy number variations in >120,000 regions across multiple organs, in benign and malignant tissues. We demonstrate that genome-wide copy number variation reveals distinct clonal patterns within tumours and in nearby benign tissue using an organ-wide approach focused on the prostate. Our results suggest a model for how genomic instability arises in histologically benign tissue that may represent early events in cancer evolution. We highlight the power of capturing the molecular and spatial continuums in a tissue context and challenge the rationale for treatment paradigms, including focal therapy.

SnapShot: Spatial transcriptomics.

Spatially resolved transcriptomics methodologies using RNA sequencing principles have and will continue to contribute to decode the molecular landscape of tissues. Linking quantitative sequencing data with tissue morphology empowers profiling of cellular morphology and transcription over time and space in health and disease. To view this SnapShot, open or download the PDF.

Integrated single cell and spatial transcriptomics reveal autoreactive differentiated B cells in joints of early rheumatoid arthritis.

B cells play a significant role in established Rheumatoid Arthritis (RA). However, it is unclear to what extent differentiated B cells are present in joint tissue already at the onset of disease. Here, we studied synovial biopsies (n = 8) captured from untreated patients at time of diagnosis. 3414 index-sorted B cells underwent RNA sequencing and paired tissue pieces were subjected to spatial transcriptomics (n = 4). We performed extensive bioinformatics analyses to dissect the local B cell composition. Select plasma cell immunoglobulin sequences were expressed as monoclonal antibodies and tested by ELISA. Memory and plasma cells were found irrespective of autoantibody status of the patients. Double negative memory B cells were prominent, but did not display a distinct transcriptional profile. The tissue architecture implicate both local B cell maturation via T cell help and plasma cell survival niches with a strong CXCL12-CXCR4 axis. The immunoglobulin sequence analyses revealed clonality between the memory B and plasma cell pools further supporting local maturation. One of the plasma cell-derived antibodies displayed citrulline autoreactivity, demonstrating local autoreactive plasma cell differentiation in joint biopsies captured from untreated early RA. Hence, plasma cell niches are not a consequence of chronic inflammation, but are already present at the time of diagnosis.

Clonal relations in the mouse brain revealed by single-cell and spatial transcriptomics.

The mammalian brain contains many specialized cells that develop from a thin sheet of neuroepithelial progenitor cells. Single-cell transcriptomics revealed hundreds of molecularly diverse cell types in the nervous system, but the lineage relationships between mature cell types and progenitor cells are not well understood. Here we show in vivo barcoding of early progenitors to simultaneously profile cell phenotypes and clonal relations in the mouse brain using single-cell and spatial transcriptomics. By reconstructing thousands of clones, we discovered fate-restricted progenitor cells in the mouse hippocampal neuroepithelium and show that microglia are derived from few primitive myeloid precursors that massively expand to generate widely dispersed progeny. We combined spatial transcriptomics with clonal barcoding and disentangled migration patterns of clonally related cells in densely labeled tissue sections. Our approach enables high-throughput dense reconstruction of cell phenotypes and clonal relations at the single-cell and tissue level in individual animals and provides an integrated approach for understanding tissue architecture.

The spatial transcriptomic landscape of the healing mouse intestine following damage.

The intestinal barrier is composed of a complex cell network defining highly compartmentalized and specialized structures. Here, we use spatial transcriptomics to define how the transcriptomic landscape is spatially organized in the steady state and healing murine colon. At steady state conditions, we demonstrate a previously unappreciated molecular regionalization of the colon, which dramatically changes during mucosal healing. Here, we identified spatially-organized transcriptional programs defining compartmentalized mucosal healing, and regions with dominant wired pathways. Furthermore, we showed that decreased p53 activation defined areas with increased presence of proliferating epithelial stem cells. Finally, we mapped transcriptomics modules associated with human diseases demonstrating the translational potential of our dataset. Overall, we provide a publicly available resource defining principles of transcriptomic regionalization of the colon during mucosal healing and a framework to develop and progress further hypotheses.

Three-dimensional spatial transcriptomics uncovers cell type localizations in the human rheumatoid arthritis synovium.

The inflamed rheumatic joint is a highly heterogeneous and complex tissue with dynamic recruitment and expansion of multiple cell types that interact in multifaceted ways within a localized area. Rheumatoid arthritis synovium has primarily been studied either by immunostaining or by molecular profiling after tissue homogenization. Here, we use Spatial Transcriptomics, where tissue-resident RNA is spatially labeled in situ with barcodes in a transcriptome-wide fashion, to study local tissue interactions at the site of chronic synovial inflammation. We report comprehensive spatial RNA-Seq data coupled to cell type-specific localization patterns at and around organized structures of infiltrating leukocyte cells in the synovium. Combining morphological features and high-throughput spatially resolved transcriptomics may be able to provide higher statistical power and more insights into monitoring disease severity and treatment-specific responses in seropositive and seronegative rheumatoid arthritis.

The discovAIR project: a roadmap towards the Human Lung Cell Atlas.

The Human Cell Atlas (HCA) consortium aims to establish an atlas of all organs in the healthy human body at single-cell resolution to increase our understanding of basic biological processes that govern development, physiology and anatomy, and to accelerate diagnosis and treatment of disease. The Lung Biological Network of the HCA aims to generate the Human Lung Cell Atlas as a reference for the cellular repertoire, molecular cell states and phenotypes, and cell-cell interactions that characterise normal lung homeostasis in healthy lung tissue. Such a reference atlas of the healthy human lung will facilitate mapping the changes in the cellular landscape in disease. The discovAIR project is one of six pilot actions for the HCA funded by the European Commission in the context of the H2020 framework programme. discovAIR aims to establish the first draft of an integrated Human Lung Cell Atlas, combining single-cell transcriptional and epigenetic profiling with spatially resolving techniques on matched tissue samples, as well as including a number of chronic and infectious diseases of the lung. The integrated Human Lung Cell Atlas will be available as a resource for the wider respiratory community, including basic and translational scientists, clinical medicine, and the private sector, as well as for patients with lung disease and the interested lay public. We anticipate that the Human Lung Cell Atlas will be the founding stone for a more detailed understanding of the pathogenesis of lung diseases, guiding the design of novel diagnostics and preventive or curative interventions.

Super-resolved spatial transcriptomics by deep data fusion.

Current methods for spatial transcriptomics are limited by low spatial resolution. Here we introduce a method that integrates spatial gene expression data with histological image data from the same tissue section to infer higher-resolution expression maps. Using a deep generative model, our method characterizes the transcriptome of micrometer-scale anatomical features and can predict spatial gene expression from histology images alone.

Spatial Transcriptomics to define transcriptional patterns of zonation and structural components in the mouse liver.

Reconstruction of heterogeneity through single cell transcriptional profiling has greatly advanced our understanding of the spatial liver transcriptome in recent years. However, global transcriptional differences across lobular units remain elusive in physical space. Here, we apply Spatial Transcriptomics to perform transcriptomic analysis across sectioned liver tissue. We confirm that the heterogeneity in this complex tissue is predominantly determined by lobular zonation. By introducing novel computational approaches, we enable transcriptional gradient measurements between tissue structures, including several lobules in a variety of orientations. Further, our data suggests the presence of previously transcriptionally uncharacterized structures within liver tissue, contributing to the overall spatial heterogeneity of the organ. This study demonstrates how comprehensive spatial transcriptomic technologies can be used to delineate extensive spatial gene expression patterns in the liver, indicating its future impact for studies of liver function, development and regeneration as well as its potential in pre-clinical and clinical pathology.

Spatial deconvolution of HER2-positive breast cancer delineates tumor-associated cell type interactions.

In the past decades, transcriptomic studies have revolutionized cancer treatment and diagnosis. However, tumor sequencing strategies typically result in loss of spatial information, critical to understand cell interactions and their functional relevance. To address this, we investigate spatial gene expression in HER2-positive breast tumors using Spatial Transcriptomics technology. We show that expression-based clustering enables data-driven tumor annotation and assessment of intra- and interpatient heterogeneity; from which we discover shared gene signatures for immune and tumor processes. By integration with single cell data, we spatially map tumor-associated cell types to find tertiary lymphoid-like structures, and a type I interferon response overlapping with regions of T-cell and macrophage subset colocalization. We construct a predictive model to infer presence of tertiary lymphoid-like structures, applicable across tissue types and technical platforms. Taken together, we combine different data modalities to define a high resolution map of cellular interactions in tumors and provide tools generalizing across tissues and diseases.

A single-cell and spatially resolved atlas of human breast cancers.

Breast cancers are complex cellular ecosystems where heterotypic interactions play central roles in disease progression and response to therapy. However, our knowledge of their cellular composition and organization is limited. Here we present a single-cell and spatially resolved transcriptomics analysis of human breast cancers. We developed a single-cell method of intrinsic subtype classification (SCSubtype) to reveal recurrent neoplastic cell heterogeneity. Immunophenotyping using cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) provides high-resolution immune profiles, including new PD-L1/PD-L2+ macrophage populations associated with clinical outcome. Mesenchymal cells displayed diverse functions and cell-surface protein expression through differentiation within three major lineages. Stromal-immune niches were spatially organized in tumors, offering insights into antitumor immune regulation. Using single-cell signatures, we deconvoluted large breast cancer cohorts to stratify them into nine clusters, termed 'ecotypes', with unique cellular compositions and clinical outcomes. This study provides a comprehensive transcriptional atlas of the cellular architecture of breast cancer.

The spatial RNA integrity number assay for in situ evaluation of transcriptome quality.

The RNA integrity number (RIN) is a frequently used quality metric to assess the completeness of rRNA, as a proxy for the corresponding mRNA in a tissue. Current methods operate at bulk resolution and provide a single average estimate for the whole sample. Spatial transcriptomics technologies have emerged and shown their value by placing gene expression into a tissue context, resulting in transcriptional information from all tissue regions. Thus, the ability to estimate RNA quality in situ has become of utmost importance to overcome the limitation with a bulk rRNA measurement. Here we show a new tool, the spatial RNA integrity number (sRIN) assay, to assess the rRNA completeness in a tissue wide manner at cellular resolution. We demonstrate the use of sRIN to identify spatial variation in tissue quality prior to more comprehensive spatial transcriptomics workflows.