TissUUmaps 3: Improvements in interactive visualization, exploration, and quality assessment of large-scale spatial omics data

High-Resolution Spatial Map of the Human Facial Sebaceous Gland Reveals Marker Genes and Decodes Sebocyte Differentiation

The sebaceous gland is essential for skin homeostasis by producing sebum to lubricate and protect the skin. Dysfunctions in sebaceous gland activity are associated with skin disorders such as acne, seborrheic dermatitis, and alopecia. However, its cellular and molecular mechanisms in humans remain poorly understood as most studies have been conducted in mouse models. This study provides a comprehensive molecular analysis of the human sebaceous gland, focusing on cellular interactions, sebocyte differentiation, and, to our knowledge, previously unreported gene markers. By integrating Stereo-seq spatial transcriptomics, single-cell RNA sequencing, and validation by MERFISH, we identified four distinct stages of sebocyte differentiation, each characterized by unique gene signatures. These results reveal that sebocyte differentiation is a dynamic and complex process. Our findings enhance the understanding of sebaceous gland biology and provide a valuable reference for future research and the development of therapies for sebaceous gland-related disorders, including acne.

High-resolution spatial transcriptomics and cell lineage analysis reveal spatiotemporal cell fate determination during craniofacial development

The differentiation of post-migratory cranial neural crest cells (CNCCs) into distinct mesenchymal lineages is crucial for craniofacial development. Here we report a high-resolution spatiotemporal transcriptomic and cell-type atlas of CNCC-derived mesenchymal lineage diversification during mouse palatogenesis. We systematically defined each mesenchymal cell type by mapping their transcriptomic profiles to spatial identities. Integrative analysis of spatial transcriptomic data from E12.5 to E15.5 further revealed mesenchymal lineage establishment at or prior to initiation of palatogenesis. We also identified a heterogeneous Sox9+ mesenchymal progenitor population at the onset of palatal development, with subpopulations already activating early lineage-specific markers. In vivo lineage tracing using these early lineage-specific markers demonstrated that distinct mesenchymal populations are established as early as E10.5 to E11.5, preceding palatal development, and contribute to their respective lineages. Together, our findings reveal the comprehensive, dynamic molecular and cellular landscape of palate development and shed light on cell fate regulation during embryogenesis.

Optimizing Xenium In Situ data utility by quality assessment and best-practice analysis workflows

The Xenium In Situ platform is a new spatial transcriptomics product commercialized by 10x Genomics, capable of mapping hundreds of genes in situ at subcellular resolution. Given the multitude of commercially available spatial transcriptomics technologies, recommendations in choice of platform and analysis guidelines are increasingly important. Herein, we explore 25 Xenium datasets generated from multiple tissues and species, comparing scalability, resolution, data quality, capacities and limitations with eight other spatially resolved transcriptomics technologies and commercial platforms. In addition, we benchmark the performance of multiple open-source computational tools, when applied to Xenium datasets, in tasks including preprocessing, cell segmentation, selection of spatially variable features and domain identification. This study serves as an independent analysis of the performance of Xenium, and provides best practices and recommendations for analysis of such datasets.

A clinical prostate biopsy dataset with undetected cancer

Prostate cancer is a heterogeneous disease showing variability both among individuals and within a patient. While most cases are indolent, aggressive tumors require early intervention. Accurately predicting tumor behavior is challenging, contributing to overdiagnosis but also undertreatment. Current imaging methods may miss the most malignant areas, leading to biopsies often capturing non-malignant prostate tissue even if cancer is present elsewhere in the organ. This non-malignant tissue, however, holds potential as a source for novel diagnostic and prognostic markers. Our clinical dataset comprises men with raised prostate-specific antigen but whose initial prostate needle biopsies only contained benign tissue. Half of the paired patients remained cancer-free for over eight years, while the others were diagnosed with prostate cancer within 30 months of follow-up. We share these initial benign biopsies to enable the exploration of morphological changes in non-malignant tissue and the potential for improved diagnostic accuracy in the early identification of patients with prostate cancer.

DaedalusData: Exploration, Knowledge Externalization and Labeling of Particles in Medical Manufacturing - A Design Study

In medical diagnostics of both early disease detection and routine patient care, particle-based contamination of in-vitro diagnostics consumables poses a significant threat to patients. Objective data-driven decision-making on the severity of contamination is key for reducing patient risk, while saving time and cost in quality assessment. Our collaborators introduced us to their quality control process, including particle data acquisition through image recognition, feature extraction, and attributes reflecting the production context of particles. Shortcomings in the current process are limitations in exploring thousands of images, data-driven decision making, and ineffective knowledge externalization. Following the design study methodology, our contributions are a characterization of the problem space and requirements, the development and validation of DaedalusData, a comprehensive discussion of our study's learnings, and a generalizable framework for knowledge externalization. DaedalusData is a visual analytics system that enables domain experts to explore particle contamination patterns, label particles in label alphabets, and externalize knowledge through semi-supervised label-informed data projections. The results of our case study and user study show high usability of DaedalusData and its efficient support of experts in generating comprehensive overviews of thousands of particles, labeling of large quantities of particles, and externalizing knowledge to augment the dataset further. Reflecting on our approach, we discuss insights on dataset augmentation via human knowledge externalization, and on the scalability and trade-offs that come with the adoption of this approach in practice.

Vitessce: integrative visualization of multimodal and spatially resolved single-cell data

Multiomics technologies with single-cell and spatial resolution make it possible to measure thousands of features across millions of cells. However, visual analysis of high-dimensional transcriptomic, proteomic, genome-mapped and imaging data types simultaneously remains a challenge. Here we describe Vitessce, an interactive web-based visualization framework for exploration of multimodal and spatially resolved single-cell data. We demonstrate integrative visualization of millions of data points, including cell-type annotations, gene expression quantities, spatially resolved transcripts and cell segmentations, across multiple coordinated views. The open-source software is available at http://vitessce.io .

Advanced single-cell and spatial analysis with high-multiplex characterization of circulating tumor cells and tumor tissue in prostate cancer: Unveiling resistance mechanisms with the CoDuCo in situ assay

BACKGROUND

Metastatic prostate cancer is a highly heterogeneous and dynamic disease and practicable tools for patient stratification and resistance monitoring are urgently needed. Liquid biopsy analysis of circulating tumor cells (CTCs) and circulating tumor DNA are promising, however, comprehensive testing is essential due to diverse mechanisms of resistance. Previously, we demonstrated the utility of mRNA-based in situ padlock probe hybridization for characterizing CTCs.

METHODS

We have developed a novel combinatorial dual-color (CoDuCo) assay for in situ mRNA detection, with enhanced multiplexing capacity, enabling the simultaneous analysis of up to 15 distinct markers. This approach was applied to CTCs, corresponding tumor tissue, cancer cell lines, and peripheral blood mononuclear cells for single-cell and spatial gene expression analysis. Using supervised machine learning, we trained a random forest classifier to identify CTCs. Image analysis and visualization of results was performed using open-source Python libraries, CellProfiler, and TissUUmaps.

RESULTS

Our study presents data from multiple prostate cancer patients, demonstrating the CoDuCo assay's ability to visualize diverse resistance mechanisms, such as neuroendocrine differentiation markers (SYP, CHGA, NCAM1) and AR-V7 expression. In addition, druggable targets and predictive markers (PSMA, DLL3, SLFN11) were detected in CTCs and formalin-fixed, paraffin-embedded tissue. The machine learning-based CTC classification achieved high performance, with a recall of 0.76 and a specificity of 0.99.

CONCLUSIONS

The combination of high multiplex capacity and microscopy-based single-cell analysis is a unique and powerful feature of the CoDuCo in situ assay. This synergy enables the simultaneous identification and characterization of CTCs with epithelial, epithelial-mesenchymal, and neuroendocrine phenotypes, the detection of CTC clusters, the visualization of CTC heterogeneity, as well as the spatial investigation of tumor tissue. This assay holds significant potential as a tool for monitoring dynamic molecular changes associated with drug response and resistance in prostate cancer.

High-parametric protein maps reveal the spatial organization in early-developing human lung

The respiratory system, including the lungs, is essential for terrestrial life. While recent research has advanced our understanding of lung development, much still relies on animal models and transcriptome analyses. In this study conducted within the Human Developmental Cell Atlas (HDCA) initiative, we describe the protein-level spatiotemporal organization of the lung during the first trimester of human gestation. Using high-parametric tissue imaging with a 30-plex antibody panel, we analyzed human lung samples from 6 to 13 post-conception weeks, generating data from over 2 million cells across five developmental timepoints. We present a resource detailing spatially resolved cell type composition of the developing human lung, including proliferative states, immune cell patterns, spatial arrangement traits, and their temporal evolution. This represents an extensive single-cell resolved protein-level examination of the developing human lung and provides a valuable resource for further research into the developmental roots of human respiratory health and disease.

Spatial Transcriptome Mapping of the Desmoplastic Growth Pattern of Colorectal Liver Metastases by In Situ Sequencing Reveals a Biologically Relevant Zonation of the Desmoplastic Rim

PURPOSE

We describe the fibrotic rim formed in the desmoplastic histopathologic growth pattern (DHGP) of colorectal cancer liver metastasis (CLM) using in situ sequencing (ISS). The origin of the desmoplastic rim is still a matter of debate, and the detailed cellular organization has not yet been fully elucidated. Understanding the biology of the DHGP in CLM can lead to targeted treatment and improve survival.

EXPERIMENTAL DESIGN

We used ISS, targeting 150 genes, to characterize the desmoplastic rim by unsupervised clustering of gene coexpression patterns. The cohort comprised 10 chemo-naïve liver metastasis resection samples with a DHGP.

RESULTS

Unsupervised clustering of spatially mapped genes revealed molecular and cellular diversity within the desmoplastic rim. We confirmed the presence of the ductular reaction and cancer-associated fibroblasts. Importantly, we discovered angiogenesis and outer and inner zonation in the rim, characterized by nerve growth factor receptor and periostin expression.

CONCLUSIONS

ISS enabled the analysis of the cellular organization of the fibrous rim surrounding CLM with a DHGP and suggests a transition from the outer part of the rim, with nonspecific liver injury response, into the inner part, with gene expression indicating collagen synthesis and extracellular matrix remodeling influenced by the interaction with cancer cells, creating a cancer cell-supportive environment. Moreover, we found angiogenic processes in the rim. Our results provide a potential explanation of the origin of the rim in DHGP and lead to exploring novel targeted treatments for patients with CLM to improve survival.

SpatialOne: end-to-end analysis of visium data at scale

MOTIVATION

Spatial transcriptomics allow to quantify mRNA expression within the spatial context. Nonetheless, in-depth analysis of spatial transcriptomics data remains challenging and difficult to scale due to the number of methods and libraries required for that purpose.

RESULTS

Here we present SpatialOne, an end-to-end pipeline designed to simplify the analysis of 10x Visium data by combining multiple state-of-the-art computational methods to segment, deconvolve, and quantify spatial information; this approach streamlines the analysis of reproducible spatial-data at scale.

AVAILABILITY AND IMPLEMENTATION

SpatialOne source code and execution examples are available at https://github.com/Sanofi-Public/spatialone-pipeline, experimental data is available at https://zenodo.org/records/12605154. SpatialOne is distributed as a docker container image.

Points2Regions: Fast, interactive clustering of imaging-based spatial transcriptomics data

Imaging-based spatial transcriptomics techniques generate data in the form of spatial points belonging to different mRNA classes. A crucial part of analyzing the data involves the identification of regions with similar composition of mRNA classes. These biologically interesting regions can manifest at different spatial scales. For example, the composition of mRNA classes on a cellular scale corresponds to cell types, whereas compositions on a millimeter scale correspond to tissue-level structures. Traditional techniques for identifying such regions often rely on complementary data, such as pre-segmented cells, or lengthy optimization. This limits their applicability to tasks on a particular scale, restricting their capabilities in exploratory analysis. This article introduces "Points2Regions," a computational tool for identifying regions with similar mRNA compositions. The tool's novelty lies in its rapid feature extraction by rasterizing points (representing mRNAs) onto a pyramidal grid and its efficient clustering using a combination of hierarchical and k -means clustering. This enables fast and efficient region discovery across multiple scales without relying on additional data, making it a valuable resource for exploratory analysis. Points2Regions has demonstrated performance similar to state-of-the-art methods on two simulated datasets, without relying on segmented cells, while being several times faster. Experiments on real-world datasets show that regions identified by Points2Regions are similar to those identified in other studies, confirming that Points2Regions can be used to extract biologically relevant regions. The tool is shared as a Python package integrated into TissUUmaps and a Napari plugin, offering interactive clustering and visualization, significantly enhancing user experience in data exploration.

Deciphering brain cellular and behavioral mechanisms: Insights from single-cell and spatial RNA sequencing

The brain is a complex computing system composed of a multitude of interacting neurons. The computational outputs of this system determine the behavior and perception of every individual. Each brain cell expresses thousands of genes that dictate the cell's function and physiological properties. Therefore, deciphering the molecular expression of each cell is of great significance for understanding its characteristics and role in brain function. Additionally, the positional information of each cell can provide crucial insights into their involvement in local brain circuits. In this review, we briefly overview the principles of single-cell RNA sequencing and spatial transcriptomics, the potential issues and challenges in their data processing, and their applications in brain research. We further outline several promising directions in neuroscience that could be integrated with single-cell RNA sequencing, including neurodevelopment, the identification of novel brain microstructures, cognition and behavior, neuronal cell positioning, molecules and cells related to advanced brain functions, sleep-wake cycles/circadian rhythms, and computational modeling of brain function. We believe that the deep integration of these directions with single-cell and spatial RNA sequencing can contribute significantly to understanding the roles of individual cells or cell types in these specific functions, thereby making important contributions to addressing critical questions in those fields. This article is categorized under: RNA Evolution and Genomics > Computational Analyses of RNA RNA in Disease and Development > RNA in Development RNA in Disease and Development > RNA in Disease.

Spatial omics techniques and data analysis for cancer immunotherapy applications

In-depth profiling of cancer cells/tissues is expanding our understanding of the genomic, epigenomic, transcriptomic, and proteomic landscape of cancer. However, the complexity of the cancer microenvironment, particularly its immune regulation, has made it difficult to exploit the potential of cancer immunotherapy. High-throughput spatial omics technologies and analysis pipelines have emerged as powerful tools for tackling this challenge. As a result, a potential revolution in cancer diagnosis, prognosis, and treatment is on the horizon. In this review, we discuss the technological advances in spatial profiling of cancer around and beyond the central dogma to harness the full benefits of immunotherapy. We also discuss the promise and challenges of spatial data analysis and interpretation and provide an outlook for the future.

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.

DEPICTER: Deep representation clustering for histology annotation

Automatic segmentation of histopathology whole-slide images (WSI) usually involves supervised training of deep learning models with pixel-level labels to classify each pixel of the WSI into tissue regions such as benign or cancerous. However, fully supervised segmentation requires large-scale data manually annotated by experts, which can be expensive and time-consuming to obtain. Non-fully supervised methods, ranging from semi-supervised to unsupervised, have been proposed to address this issue and have been successful in WSI segmentation tasks. But these methods have mainly been focused on technical advancements in algorithmic performance rather than on the development of practical tools that could be used by pathologists or researchers in real-world scenarios. In contrast, we present DEPICTER (Deep rEPresentatIon ClusTERing), an interactive segmentation tool for histopathology annotation that produces a patch-wise dense segmentation map at WSI level. The interactive nature of DEPICTER leverages self- and semi-supervised learning approaches to allow the user to participate in the segmentation producing reliable results while reducing the workload. DEPICTER consists of three steps: first, a pretrained model is used to compute embeddings from image patches. Next, the user selects a number of benign and cancerous patches from the multi-resolution image. Finally, guided by the deep representations, label propagation is achieved using our novel seeded iterative clustering method or by directly interacting with the embedding space via feature space gating. We report both real-time interaction results with three pathologists and evaluate the performance on three public cancer classification dataset benchmarks through simulations. The code and demos of DEPICTER are publicly available at https://github.com/eduardchelebian/depicter.

A perspective on FAIR quality control in multiplexed imaging data processing

Multiplexed imaging approaches are getting increasingly adopted for imaging of large tissue areas, yielding big imaging datasets both in terms of the number of samples and the size of image data per sample. The processing and analysis of these datasets is complex owing to frequent technical artifacts and heterogeneous profiles from a high number of stained targets To streamline the analysis of multiplexed images, automated pipelines making use of state-of-the-art algorithms have been developed. In these pipelines, the output quality of one processing step is typically dependent on the output of the previous step and errors from each step, even when they appear minor, can propagate and confound the results. Thus, rigorous quality control (QC) at each of these different steps of the image processing pipeline is of paramount importance both for the proper analysis and interpretation of the analysis results and for ensuring the reusability of the data. Ideally, QC should become an integral and easily retrievable part of the imaging datasets and the analysis process. Yet, limitations of the currently available frameworks make integration of interactive QC difficult for large multiplexed imaging data. Given the increasing size and complexity of multiplexed imaging datasets, we present the different challenges for integrating QC in image analysis pipelines as well as suggest possible solutions that build on top of recent advances in bioimage analysis.

Spatial-Live: A lightweight and versatile tool for single cell spatial-omics data visualization

Single cell spatial-omics data visualization plays a pivotal role in unraveling the intricate spatial organization and heterogeneity of cellular systems. Although various software tools and packages have been developed for this purpose, challenges persist in terms of user-friendly accessibility, data integration, and interactivity. In this study, we introduce Spatial-Live, a lightweight and versatile viewer tool designed for flexible single-cell spatial-omics data visualization. Spatial-Live overcomes the fundamental limitations of two-dimensional (2D) orthographic modes by employing a layer-stacking strategy, enabling efficient rendering of diverse data types with interactive features, and enhancing visualization with richer information in a unified three-dimensional (3D) space.

The spatial landscape of gene expression isoforms in tissue sections

In situ capturing technologies add tissue context to gene expression data, with the potential of providing a greater understanding of complex biological systems. However, splicing variants and full-length sequence heterogeneity cannot be characterized at spatial resolution with current transcriptome profiling methods. To that end, we introduce spatial isoform transcriptomics (SiT), an explorative method for characterizing spatial isoform variation and sequence heterogeneity using long-read sequencing. We show in mouse brain how SiT can be used to profile isoform expression and sequence heterogeneity in different areas of the tissue. SiT reveals regional isoform switching of Plp1 gene between different layers of the olfactory bulb, and the use of external single-cell data allows the nomination of cell types expressing each isoform. Furthermore, SiT identifies differential isoform usage for several major genes implicated in brain function (Snap25, Bin1, Gnas) that are independently validated by in situ sequencing. SiT also provides for the first time an in-depth A-to-I RNA editing map of the adult mouse brain. Data exploration can be performed through an online resource (https://www.isomics.eu), where isoform expression and RNA editing can be visualized in a spatial context.

Visualization and quality control tools for large-scale multiplex tissue analysis in TissUUmaps3

Large-scale multiplex tissue analysis aims to understand processes such as development and tumor formation by studying the occurrence and interaction of cells in local environments in, for example, tissue samples from patient cohorts. A typical procedure in the analysis is to delineate individual cells, classify them into cell types, and analyze their spatial relationships. All steps come with a number of challenges, and to address them and identify the bottlenecks of the analysis, it is necessary to include quality control tools in the analysis workflow. This makes it possible to optimize the steps and adjust settings in order to get better and more precise results. Additionally, the development of automated approaches for tissue analysis requires visual verification to reduce skepticism with regard to the accuracy of the results. Quality control tools could be used to build users' trust in automated approaches. In this paper, we present three plugins for visualization and quality control in large-scale multiplex tissue analysis of microscopy images. The first plugin focuses on the quality of cell staining, the second one was made for interactive evaluation and comparison of different cell classification results, and the third one serves for reviewing interactions of different cell types.