Multiplexed ion beam imaging of human breast tumors

Computational methods for alignment and integration of spatially resolved transcriptomics data

Most of the complex biological regulatory activities occur in three dimensions (3D). To better analyze biological processes, it is essential not only to decipher the molecular information of numerous cells but also to understand how their spatial contexts influence their behavior. With the development of spatially resolved transcriptomics (SRT) technologies, SRT datasets are being generated to simultaneously characterize gene expression and spatial arrangement information within tissues, organs or organisms. To fully leverage spatial information, the focus extends beyond individual two-dimensional (2D) slices. Two tasks known as slices alignment and data integration have been introduced to establish correlations between multiple slices, enhancing the effectiveness of downstream tasks. Currently, numerous related methods have been developed. In this review, we first elucidate the details and principles behind several representative methods. Then we report the testing results of these methods on various SRT datasets, and assess their performance in representative downstream tasks. Insights into the strengths and weaknesses of each method and the reasons behind their performance are discussed. Finally, we provide an outlook on future developments. The codes and details of experiments are now publicly available at https://github.com/YangLabHKUST/SRT_alignment_and_integration.

The tumor microenvironment of benign and malignant salivary gland tumors

BACKGROUND

Treatment of salivary gland tumors (SGTs) remains challenging. Little is known about the immune landscape of SGTs. We aimed to characterize the tumor microenvironment in benign and malignant SGTs.

METHODS

Eleven benign and nine malignant tumors were collected from patients undergoing curative intent surgery. Specimens were analyzed using mass cytometry by time-of-flight. Immune cell populations were manually gated, and T cells were clustered using the FlowSOM algorithm. Population frequencies were compared between high-grade and low-grade malignancies, corrected for multiple hypothesis testing.

RESULTS

There were trends towards increased CD4+ and CD8+ T cells among malignant tumors. High-grade malignancies exhibited trends towards higher frequencies of CD8+ PD-1+ CD39+ CD103+ exhausted T cells, CD4+ FoxP3+ TCF-1+ CD127- Tregs, and CD69+ CD25- CD4+ T cells compared to low-grade malignancies.

CONCLUSION

SGTs exhibit significant immunologic diversity. High-grade malignancies tended to have greater infiltration of exhausted CD8+ T cells and Tregs, which may guide future studies for immunotherapy strategies.

Probe-based mass spectrometry approaches for single-cell and single-organelle measurements

Exploring the chemical content of individual cells not only reveals underlying cell-to-cell chemical heterogeneity but is also a key component in understanding how cells combine to form emergent properties of cellular networks and tissues. Recent technological advances in many analytical techniques including mass spectrometry (MS) have improved instrumental limits of detection and laser/ion probe dimensions, allowing the analysis of micron and submicron sized areas. In the case of MS, these improvements combined with MS's broad analyte detection capabilities have enabled the rise of single-cell and single-organelle chemical characterization. As the chemical coverage and throughput of single-cell measurements increase, more advanced statistical and data analysis methods have aided in data visualization and interpretation. This review focuses on secondary ion MS and matrix-assisted laser desorption/ionization MS approaches for single-cell and single-organelle characterization, which is followed by advances in mass spectral data visualization and analysis.

Increasing Signal Intensity of Fluorescent Oligo-Labeled Antibodies to Enable Combination Multiplexing

Full-spectrum flow cytometry has increased antibody-based multiplexing, yet further increases remain potentially impactful. We recently proposed how fluorescence multiplexing using spectral imaging and combinatorics (MuSIC) could do so using tandem dyes and an oligo-based antibody labeling method. In this work, we found that such labeled antibodies had significantly lower signal intensities than conventionally labeled antibodies in human cell experiments. To improve signal intensity, we tested moving the fluorophores from the original external (ext.) 5' or 3' end-labeled orientation to internal (int.) fluorophore modifications. Cell-free spectrophotometer measurements showed a ∼6-fold signal intensity increase of the new int. configuration compared to the previous ext. configuration. Time-resolved fluorescence and fluorescence correlation spectroscopy showed that the ∼3-fold brightness difference is due to static quenching most likely by the oligo or solution in the ext. configuration. Spectral flow cytometry experiments using peripheral blood mononuclear cells show int. MuSIC probe-labeled antibodies (i) retained increased signal intensity while having no significant difference in the estimated % of CD8+ lymphocytes and (ii) labeled with Atto488, Atto647, and Atto488/647 combinations can be demultiplexed in triple-stained samples. The antibody labeling approach is general and can be broadly applied to many biological and diagnostic applications where spectral detection is available.

Using random forests to uncover the predictive power of distance-varying cell interactions in tumor microenvironments

Tumor microenvironments (TMEs) contain vast amounts of information on patient's cancer through their cellular composition and the spatial distribution of tumor cells and immune cell populations. Exploring variations in TMEs between patient groups, as well as determining the extent to which this information can predict outcomes such as patient survival or treatment success with emerging immunotherapies, is of great interest. Moreover, in the face of a large number of cell interactions to consider, we often wish to identify specific interactions that are useful in making such predictions. We present an approach to achieve these goals based on summarizing spatial relationships in the TME using spatial K functions, and then applying functional data analysis and random forest models to both predict outcomes of interest and identify important spatial relationships. This approach is shown to be effective in simulation experiments at both identifying important spatial interactions while also controlling the false discovery rate. We further used the proposed approach to interrogate two real data sets of Multiplexed Ion Beam Images of TMEs in triple negative breast cancer and lung cancer patients. The methods proposed are publicly available in a companion R package funkycells.

Machine learning integrative approaches to advance computational immunology

The study of immunology, traditionally reliant on proteomics to evaluate individual immune cells, has been revolutionized by single-cell RNA sequencing. Computational immunologists play a crucial role in analysing these datasets, moving beyond traditional protein marker identification to encompass a more detailed view of cellular phenotypes and their functional roles. Recent technological advancements allow the simultaneous measurements of multiple cellular components-transcriptome, proteome, chromatin, epigenetic modifications and metabolites-within single cells, including in spatial contexts within tissues. This has led to the generation of complex multiscale datasets that can include multimodal measurements from the same cells or a mix of paired and unpaired modalities. Modern machine learning (ML) techniques allow for the integration of multiple "omics" data without the need for extensive independent modelling of each modality. This review focuses on recent advancements in ML integrative approaches applied to immunological studies. We highlight the importance of these methods in creating a unified representation of multiscale data collections, particularly for single-cell and spatial profiling technologies. Finally, we discuss the challenges of these holistic approaches and how they will be instrumental in the development of a common coordinate framework for multiscale studies, thereby accelerating research and enabling discoveries in the computational immunology field.

spillR: spillover compensation in mass cytometry data

MOTIVATION

Channel interference in mass cytometry can cause spillover and may result in miscounting of protein markers. Chevrier et al. introduce an experimental and computational procedure to estimate and compensate for spillover implemented in their R package CATALYST. They assume spillover can be described by a spillover matrix that encodes the ratio between the signal in the unstained spillover receiving and stained spillover emitting channel. They estimate the spillover matrix from experiments with beads. We propose to skip the matrix estimation step and work directly with the full bead distributions. We develop a nonparametric finite mixture model and use the mixture components to estimate the probability of spillover. Spillover correction is often a pre-processing step followed by downstream analyses, and choosing a flexible model reduces the chance of introducing biases that can propagate downstream.

RESULTS

We implement our method in an R package spillR using expectation-maximization to fit the mixture model. We test our method on simulated, semi-simulated, and real data from CATALYST. We find that our method compensates low counts accurately, does not introduce negative counts, avoids overcompensating high counts, and preserves correlations between markers that may be biologically meaningful.

AVAILABILITY AND IMPLEMENTATION

Our new R package spillR is on bioconductor at bioconductor.org/packages/spillR. All experiments and plots can be reproduced by compiling the R markdown file spillR_paper.Rmd at github.com/ChristofSeiler/spillR_paper.

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.

Development of a high dimensional imaging mass cytometry panel to investigate spatial organization of tissue microenvironment in formalin-fixed archival clinical tissues

To decipher the interactions between various components of the tumor microenvironment (TME) and tumor cells in a preserved spatial context, a multiparametric approach is essential. In this pursuit, imaging mass cytometry (IMC) emerges as a valuable tool, capable of concurrently analyzing up to 40 parameters at subcellular resolution. In this study, a set of antibodies was selected to spatially resolve multiple cell types and TME elements, including a comprehensive panel targeted at dissecting the heterogeneity of cancer-associated fibroblasts (CAF), a pivotal TME component. This antibody panel was standardized and optimized using formalin-fixed paraffin-embedded tissue (FFPE) samples from different organs/lesions known to express the markers of interest. The final composition of the antibody panel was determined based on the performance of conjugated antibodies in both immunohistochemistry (IHC) and IMC. Tissue images were segmented employing the Steinbock framework. Unsupervised clustering of single-cell data was carried out using a bioinformatics pipeline developed in R program. This paper provides a detailed description of the staining procedure and analysis workflow. Subsequently, the panel underwent validation on clinical FFPE samples from head and neck squamous cell carcinoma (HNSCC). The panel and bioinformatics pipeline established here proved to be robust in characterizing different TME components of HNSCC while maintaining a high degree of spatial detail. The platform we describe shows promise for understanding the clinical implications of TMA heterogeneity in large patient cohorts with FFPE tissues available in diagnostic biobanks worldwide.

Increasing Signal Intensity of Fluorescent Oligo-Labeled Antibodies to Enable Combination Multiplexing

Full-spectrum flow cytometry has increased antibody-based multiplexing, yet further increases remain potentially impactful. We recently proposed how fluorescence Multiplexing using Spectral Imaging and Combinatorics (MuSIC) could do so using tandem dyes and an oligo-based antibody labeling method. In this work, we found that such labeled antibodies had significantly lower signal intensity than conventionally-labeled antibodies in human cell experiments. To improve signal intensity, we tested moving the fluorophores from the original external (ext.) 5' or 3' end-labeled orientation to internal (int.) fluorophore modifications. Cell-free spectrophotometer measurements showed a ~6-fold signal intensity increase of the new int. configuration compared to the previous ext. configuration. Time-resolved fluorescence and fluorescence correlation spectroscopy showed that ~3-fold brightness difference is due to static quenching most likely by the oligo or solution in the ext. configuration. Spectral flow cytometry experiments using peripheral blood mononuclear cells show int. MuSIC probe-labeled antibodies (i) retained increased signal intensity while having no significant difference in the estimated % of CD8+ lymphocytes and (ii) labeled with Atto488, Atto647, and Atto488/647 combinations can be demultiplexed in triple-stained samples. The antibody labeling approach is general and can be broadly applied to many biological and diagnostic applications where spectral detection is available.

SCIMAP: A Python Toolkit for Integrated Spatial Analysis of Multiplexed Imaging Data

Multiplexed imaging data are revolutionizing our understanding of the composition and organization of tissues and tumors ("Catching up with Multiplexed Tissue Imaging," 2022). A critical aspect of such "tissue profiling" is quantifying the spatial relationships among cells at different scales from the interaction of neighboring cells to recurrent communities of cells of multiple types. This often involves statistical analysis of 107 or more cells in which up to 100 biomolecules (commonly proteins) have been measured. While software tools currently cater to the analysis of spatial transcriptomics data (Liu et al., 2022), there remains a need for toolkits explicitly tailored to the complexities of multiplexed imaging data including the need to seamlessly integrate image visualization with data analysis and exploration. We introduce SCIMAP, a Python package specifically crafted to address these challenges. With SCIMAP, users can efficiently preprocess, analyze, and visualize large datasets, facilitating the exploration of spatial relationships and their statistical significance. SCIMAP's modular design enables the integration of new algorithms, enhancing its capabilities for spatial analysis.

Spatiotemporal multi-omics: exploring molecular landscapes in aging and regenerative medicine

Aging and regeneration represent complex biological phenomena that have long captivated the scientific community. To fully comprehend these processes, it is essential to investigate molecular dynamics through a lens that encompasses both spatial and temporal dimensions. Conventional omics methodologies, such as genomics and transcriptomics, have been instrumental in identifying critical molecular facets of aging and regeneration. However, these methods are somewhat limited, constrained by their spatial resolution and their lack of capacity to dynamically represent tissue alterations. The advent of emerging spatiotemporal multi-omics approaches, encompassing transcriptomics, proteomics, metabolomics, and epigenomics, furnishes comprehensive insights into these intricate molecular dynamics. These sophisticated techniques facilitate accurate delineation of molecular patterns across an array of cells, tissues, and organs, thereby offering an in-depth understanding of the fundamental mechanisms at play. This review meticulously examines the significance of spatiotemporal multi-omics in the realms of aging and regeneration research. It underscores how these methodologies augment our comprehension of molecular dynamics, cellular interactions, and signaling pathways. Initially, the review delineates the foundational principles underpinning these methods, followed by an evaluation of their recent applications within the field. The review ultimately concludes by addressing the prevailing challenges and projecting future advancements in the field. Indubitably, spatiotemporal multi-omics are instrumental in deciphering the complexities inherent in aging and regeneration, thus charting a course toward potential therapeutic innovations.

Cross-domain information fusion for enhanced cell population delineation in single-cell spatial-omics data

Cell population delineation and identification is an essential step in single-cell and spatial-omics studies. Spatial-omics technologies can simultaneously measure information from three complementary domains related to this task: expression levels of a panel of molecular biomarkers at single-cell resolution, relative positions of cells, and images of tissue sections, but existing computational methods for performing this task on single-cell spatial-omics datasets often relinquish information from one or more domains. The additional reliance on the availability of "atlas" training or reference datasets limits cell type discovery to well-defined but limited cell population labels, thus posing major challenges for using these methods in practice. Successful integration of all three domains presents an opportunity for uncovering cell populations that are functionally stratified by their spatial contexts at cellular and tissue levels: the key motivation for employing spatial-omics technologies in the first place. In this work, we introduce Cell Spatio- and Neighborhood-informed Annotation and Patterning (CellSNAP), a self-supervised computational method that learns a representation vector for each cell in tissue samples measured by spatial-omics technologies at the single-cell or finer resolution. The learned representation vector fuses information about the corresponding cell across all three aforementioned domains. By applying CellSNAP to datasets spanning both spatial proteomic and spatial transcriptomic modalities, and across different tissue types and disease settings, we show that CellSNAP markedly enhances de novo discovery of biologically relevant cell populations at fine granularity, beyond current approaches, by fully integrating cells' molecular profiles with cellular neighborhood and tissue image information.

Ultrasensitive and Multiplexed Protein Imaging with Clickable and Cleavable Fluorophores

Single-cell spatial proteomic analysis holds great promise to advance our understanding of the composition, organization, interaction, and function of the various cell types in complex biological systems. However, the current multiplexed protein imaging technologies suffer from low detection sensitivity, limited multiplexing capacity, or are technically demanding. To tackle these issues, here, we report the development of a highly sensitive and multiplexed in situ protein profiling method using off-the-shelf antibodies. In this approach, the protein targets are stained with horseradish peroxidase (HRP) conjugated antibodies and cleavable fluorophores via click chemistry. Through repeated cycles of target staining, fluorescence imaging, and fluorophore cleavage, many proteins can be profiled in single cells in situ. Applying this approach, we successfully quantified 28 different proteins in human formalin-fixed paraffin-embedded (FFPE) tonsil tissue, which represents the highest multiplexing capacity among the tyramide signal amplification (TSA) methods. Based on their unique protein expression patterns and their microenvironment, ∼820,000 cells in the tissue are classified into distinct cell clusters. We also explored the cell-cell interactions between these varied cell clusters and observed that different subregions of the tissue are composed of cells from specific clusters.

Spatial Biology of Breast Cancer

Spatial findings have shaped on our understanding of breast cancer. In this review, we discuss how spatial methods, including spatial transcriptomics and proteomics and the resultant understanding of spatial relationships, have contributed to concepts regarding cancer progression and treatment. In addition to discussing traditional approaches, we examine how emerging multiplex imaging technologies have contributed to the field and how they might influence future research.

Clonal hematopoiesis and its impact on the aging osteo-hematopoietic niche

Clonal hematopoiesis (CH) defines a premalignant state predominantly found in older persons that increases the risk of developing hematologic malignancies and age-related inflammatory diseases. However, the risk for malignant transformation or non-malignant disorders is variable and difficult to predict, and defining the clinical relevance of specific candidate driver mutations in individual carriers has proved to be challenging. In addition to the cell-intrinsic mechanisms, mutant cells rely on and alter cell-extrinsic factors from the bone marrow (BM) niche, which complicates the prediction of a mutant cell's fate in a shifting pre-malignant microenvironment. Therefore, identifying the insidious and potentially broad impact of driver mutations on supportive niches and immune function in CH aims to understand the subtle differences that enable driver mutations to yield different clinical outcomes. Here, we review the changes in the aging BM niche and the emerging evidence supporting the concept that CH can progressively alter components of the local BM microenvironment. These alterations may have profound implications for the functionality of the osteo-hematopoietic niche and overall bone health, consequently fostering a conducive environment for the continued development and progression of CH. We also provide an overview of the latest technology developments to study the spatiotemporal dependencies in the CH BM niche, ideally in the context of longitudinal studies following CH over time. Finally, we discuss aspects of CH carrier management in clinical practice, based on work from our group and others.

In situ single-cell therapeutic response imaging facilitated by the TRIPODD fluorescence imaging platform

Purpose: Small molecule drugs such as tyrosine kinase inhibitors (TKIs) targeting tumoral molecular dependencies have become standard of care for numerous cancer types. Notably, epidermal growth factor receptor (EGFR) TKIs (e.g., erlotinib, afatinib, osimertinib) are the current first-line treatment for non-small cell lung cancer (NSCLC) due to their improved therapeutic outcomes for EGFR mutated and overexpressing disease over traditional platinum-based chemotherapy. However, many NSCLC tumors develop resistance to EGFR TKI therapy causing disease progression. Currently, the relationship between in situ drug target availability (DTA), local protein expression and therapeutic response cannot be accurately assessed using existing analytical tools despite being crucial to understanding the mechanism of therapeutic efficacy. Procedure: We have previously reported development of our fluorescence imaging platform termed TRIPODD (Therapeutic Response Imaging through Proteomic and Optical Drug Distribution) that is capable of simultaneous quantification of single-cell DTA and protein expression with preserved spatial context within a tumor. TRIPODD combines two complementary fluorescence imaging techniques: intracellular paired agent imaging (iPAI) to measure DTA and cyclic immunofluorescence (cyCIF), which utilizes oligonucleotide conjugated antibodies (Ab-oligos) for spatial proteomic expression profiling on tissue samples. Herein, TRIPODD was modified and optimized to provide a downstream analysis of therapeutic response through single-cell DTA and proteomic response imaging. Results: We successfully performed sequential imaging of iPAI and cyCIF resulting in high dimensional imaging and biomarker assessment to quantify single-cell DTA and local protein expression on erlotinib treated NSCLC models. Pharmacodynamic and pharmacokinetic studies of the erlotinib iPAI probes revealed that administration of 2.5 mg/kg each of the targeted and untargeted probe 4 h prior to tumor collection enabled calculation of DTA values with high Pearson correlation to EGFR, the erlotinib molecular target, expression in the tumors. Analysis of single-cell biomarker expression revealed that a single erlotinib dose was insufficient to enact a measurable decrease in the EGFR signaling cascade protein expression, where only the DTA metric detected the presence of bound erlotinib. Conclusion: We demonstrated the capability of TRIPODD to evaluate therapeutic response imaging to erlotinib treatment as it relates to signaling inhibition, DTA, proliferation, and apoptosis with preserved spatial context.

Spatiotemporal architecture of immune cells and cancer-associated fibroblasts in high-grade serous ovarian carcinoma

High-grade serous ovarian carcinoma (HGSOC), the deadliest form of ovarian cancer, is typically diagnosed after it has metastasized and often relapses after standard-of-care platinum-based chemotherapy, likely due to advanced tumor stage, heterogeneity, and immune evasion and tumor-promoting signaling from the tumor microenvironment. To understand how spatial heterogeneity contributes to HGSOC progression and early relapse, we profiled an HGSOC tissue microarray of patient-matched longitudinal samples from 42 patients. We found spatial patterns associated with early relapse, including changes in T cell localization, malformed tertiary lymphoid structure (TLS)-like aggregates, and increased podoplanin-positive cancer-associated fibroblasts (CAFs). Using spatial features to compartmentalize the tissue, we found that plasma cells distribute in two different compartments associated with TLS-like aggregates and CAFs, and these distinct microenvironments may account for the conflicting reports about the role of plasma cells in HGSOC prognosis.

Quantitative performance assessment of Ultivue multiplex panels in formalin-fixed, paraffin-embedded human and murine tumor specimens

We present a rigorous validation strategy to evaluate the performance of Ultivue multiplex immunofluorescence panels. We have quantified the accuracy and precision of four different multiplex panels (three human and one mouse) in tumor specimens with varying levels of T cell density. Our results show that Ultivue panels are typically accurate wherein the relative difference in cell proportion between a multiplex image and a 1-plex image is less than 20% for a given biomarker. Ultivue panels exhibited relatively high intra-run precision (CV ≤ 25%) and relatively low inter-run precision (CV >> 25%) which can be remedied by using local intensity thresholding to gate biomarker positivity. We also evaluated the reproducibility of cell-cell distance estimates measured from multiplex images which show high intra- and inter-run precision. We introduce a new metric, multiplex labeling efficiency, which can be used to benchmark the overall fidelity of the multiplex data across multiple batch runs. Taken together our results provide a comprehensive characterization of Ultivue panels and offer practical guidelines for analyzing multiplex images.

SpaceANOVA: Spatial Co-occurrence Analysis of Cell Types in Multiplex Imaging Data Using Point Process and Functional ANOVA

Multiplex imaging platforms have enabled the identification of the spatial organization of different types of cells in complex tissue or the tumor microenvironment. Exploring the potential variations in the spatial co-occurrence or colocalization of different cell types across distinct tissue or disease classes can provide significant pathological insights, paving the way for intervention strategies. However, the existing methods in this context either rely on stringent statistical assumptions or suffer from a lack of generalizability. We present a highly powerful method to study differential spatial co-occurrence of cell types across multiple tissue or disease groups, based on the theories of the Poisson point process and functional analysis of variance. Notably, the method accommodates multiple images per subject and addresses the problem of missing tissue regions, commonly encountered due to data-collection complexities. We demonstrate the superior statistical power and robustness of the method in comparison with existing approaches through realistic simulation studies. Furthermore, we apply the method to three real data sets on different diseases collected using different imaging platforms. In particular, one of these data sets reveals novel insights into the spatial characteristics of various types of colorectal adenoma.

MIM-CyCIF: masked imaging modeling for enhancing cyclic immunofluorescence (CyCIF) with panel reduction and imputation

Cyclic Immunofluorescence (CyCIF) can quantify multiple biomarkers, but panel capacity is limited by technical challenges. We propose a computational panel reduction approach that can impute the information content from 25 markers using only 9 markers, learning co-expression and morphological patterns while concurrently increasing speed and panel content and decreasing cost. We demonstrate strong correlations in predictions and generalizability across breast and colorectal cancer, illustrating applicability of our approach to diverse tissue types.

The intricate dance of tumor evolution: Exploring immune escape, tumor migration, drug resistance, and treatment strategies

Recent research has unveiled fascinating insights into the intricate mechanisms governing tumor evolution. These studies have illuminated how tumors adapt and proliferate by exploiting various factors, including immune evasion, resistance to therapeutic drugs, genetic mutations, and their ability to adapt to different environments. Furthermore, investigations into tumor heterogeneity and chromosomal aberrations have revealed the profound complexity that underlies the evolution of cancer. Emerging findings have also underscored the role of viral influences in the development and progression of cancer, introducing an additional layer of complexity to the field of oncology. Tumor evolution is a dynamic and complex process influenced by various factors, including immune evasion, drug resistance, tumor heterogeneity, and viral influences. Understanding these elements is indispensable for developing more effective treatments and advancing cancer therapies. A holistic approach to studying and addressing tumor evolution is crucial in the ongoing battle against cancer. The main goal of this comprehensive review is to explore the intricate relationship between tumor evolution and critical aspects of cancer biology. By delving into this complex interplay, we aim to provide a profound understanding of how tumors evolve, adapt, and respond to treatment strategies. This review underscores the pivotal importance of comprehending tumor evolution in shaping effective approaches to cancer treatment.

Applicability of Spatial Technology in Cancer Research

This review explores spatial mapping technologies in cancer research, highlighting their crucial role in understanding the complexities of the tumor microenvironment (TME). The TME, which is an intricate ecosystem of diverse cell types, has a significant impact on tumor dynamics and treatment outcomes. This review closely examines cutting-edge spatial mapping technologies, categorizing them into capture-, imaging-, and antibody-based approaches. Each technology was scrutinized for its advantages and disadvantages, factoring in aspects such as spatial profiling area, multiplexing capabilities, and resolution. Additionally, we draw attention to the nuanced choices researchers face, with capture-based methods lending themselves to hypothesis generation, and imaging/antibody-based methods that fit neatly into hypothesis testing. Looking ahead, we anticipate a scenario in which multi-omics data are seamlessly integrated, artificial intelligence enhances data analysis, and spatiotemporal profiling opens up new dimensions.

Understanding Spondyloarthritis Pathogenesis: The Promise of Single-Cell Profiling

PURPOSE OF REVIEW

Single-cell profiling, either in suspension or within the tissue context, is a rapidly evolving field. The purpose of this review is to outline recent advancements and emerging trends with a specific focus on studies in spondyloarthritis.

RECENT FINDINGS

The introduction of sequencing-based approaches for the quantification of RNA, protein, or epigenetic modifications at single-cell resolution has provided a major boost to discovery-driven research. Fluorescent flow cytometry, mass cytometry, and image-based cytometry continue to evolve. Spatial transcriptomics and imaging mass cytometry have extended high-dimensional analysis to cells in tissues. Applications in spondyloarthritis include the indexing and functional characterization of cells, discovery of disease-associated cell states, and identification of signatures associated with therapeutic responses. Single-cell TCR-seq has provided evidence for clonal expansion of CD8+ T cells in spondyloarthritis. The use of single-cell profiling approaches in spondyloarthritis research is still in its early stages. Challenges include high cost and limited availability of diseased tissue samples. To harness the full potential of the rapidly expanding technical capabilities, large-scale collaborative efforts are imperative.

DELVE: feature selection for preserving biological trajectories in single-cell data

Single-cell technologies can measure the expression of thousands of molecular features in individual cells undergoing dynamic biological processes. While examining cells along a computationally-ordered pseudotime trajectory can reveal how changes in gene or protein expression impact cell fate, identifying such dynamic features is challenging due to the inherent noise in single-cell data. Here, we present DELVE, an unsupervised feature selection method for identifying a representative subset of molecular features which robustly recapitulate cellular trajectories. In contrast to previous work, DELVE uses a bottom-up approach to mitigate the effects of confounding sources of variation, and instead models cell states from dynamic gene or protein modules based on core regulatory complexes. Using simulations, single-cell RNA sequencing, and iterative immunofluorescence imaging data in the context of cell cycle and cellular differentiation, we demonstrate how DELVE selects features that better define cell-types and cell-type transitions. DELVE is available as an open-source python package: https://github.com/jranek/delve .

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.

Three-dimensional assessments are necessary to determine the true, spatially-resolved composition of tissues

Methods for spatially resolved cellular profiling using thinly cut sections have enabled in-depth quantitative tissue mapping to study inter-sample and intra-sample differences in normal human anatomy and disease onset and progression. These methods often profile extremely limited regions, which may impact the evaluation of heterogeneity due to tissue sub-sampling. Here, we applied CODA, a deep learning-based tissue mapping platform, to reconstruct the three-dimensional (3D) microanatomy of grossly normal and cancer-containing human pancreas biospecimens obtained from individuals who underwent pancreatic resection. To compare inter- and intra-sample heterogeneity, we assessed bulk and spatially resolved tissue composition in a cohort of two-dimensional (2D) whole slide images (WSIs) and a cohort of thick slabs of pancreas tissue that were digitally reconstructed in 3D from serial sections. To demonstrate the marked under sampling of 2D assessments, we simulated the number of WSIs and tissue microarrays (TMAs) necessary to represent the compositional heterogeneity of 3D data within 10% error to reveal that tens of WSIs and hundreds of TMA cores are sometimes needed. We show that spatial correlation of different pancreatic structures decay significantly within a span of microns, demonstrating that 2D histological sections may not be representative of their neighboring tissues. In sum, we demonstrate that 3D assessments are necessary to accurately assess tissue composition in normal and abnormal specimens and in order to accurately determine neoplastic content. These results emphasize the importance of intra-sample heterogeneity in tissue mapping efforts.

Quality Control for Single Cell Analysis of High-plex Tissue Profiles using CyLinter

Tumors are complex assemblies of cellular and acellular structures patterned on spatial scales from microns to centimeters. Study of these assemblies has advanced dramatically with the introduction of high-plex spatial profiling. Image-based profiling methods reveal the intensities and spatial distributions of 20-100 proteins at subcellular resolution in 103-107 cells per specimen. Despite extensive work on methods for extracting single-cell data from these images, all tissue images contain artefacts such as folds, debris, antibody aggregates, optical aberrations and image processing errors that arise from imperfections in specimen preparation, data acquisition, image assembly, and feature extraction. We show that these artefacts dramatically impact single-cell data analysis, obscuring meaningful biological interpretation. We describe an interactive quality control software tool, CyLinter, that identifies and removes data associated with imaging artefacts. CyLinter greatly improves single-cell analysis, especially for archival specimens sectioned many years prior to data collection, such as those from clinical trials.

Unveiling spatial complexity in solid tumor immune microenvironments through multiplexed imaging

Deciphering cellular components and the spatial interaction network of the tumor immune microenvironment (TIME) of solid tumors is pivotal for understanding biologically relevant cross-talks and, ultimately, advancing therapies. Multiplexed tissue imaging provides a powerful tool to elucidate spatial complexity in a holistic manner. We established and cross-validated a comprehensive immunophenotyping panel comprising over 121 markers for multiplexed tissue imaging using MACSima™ imaging cyclic staining (MICS) alongside an end-to-end analysis workflow. Applying this panel and workflow to primary cancer tissues, we characterized tumor heterogeneity, investigated potential therapeutical targets, conducted in-depth profiling of cell types and states, sub-phenotyped T cells within the TIME, and scrutinized cellular neighborhoods of diverse T cell subsets. Our findings highlight the advantage of spatial profiling, revealing immunosuppressive molecular signatures of tumor-associated myeloid cells interacting with neighboring exhausted, PD1high T cells in the TIME of hepatocellular carcinoma (HCC). This study establishes a robust framework for spatial exploration of TIMEs in solid tumors and underscores the potency of multiplexed tissue imaging and ultra-deep cell phenotyping in unraveling clinically relevant tumor components.

A Hitchhiker's guide to high-dimensional tissue imaging with multiplexed ion beam imaging

Advancements in multiplexed tissue imaging technologies are vital in shaping our understanding of tissue microenvironmental influences in disease contexts. These technologies now allow us to relate the phenotype of individual cells to their higher-order roles in tissue organization and function. Multiplexed Ion Beam Imaging (MIBI) is one of such technologies, which uses metal isotope-labeled antibodies and secondary ion mass spectrometry (SIMS) to image more than 40 protein markers simultaneously within a single tissue section. Here, we describe an optimized MIBI workflow for high-plex analysis of Formalin-Fixed Paraffin-Embedded (FFPE) tissues following antigen retrieval, metal isotope-conjugated antibody staining, imaging using the MIBI instrument, and subsequent data processing and analysis. While this workflow is focused on imaging human FFPE samples using the MIBI, this workflow can be easily extended to model systems, biological questions, and multiplexed imaging modalities.

Epstein-Barr Virus Orchestrates Spatial Reorganization and Immunomodulation within the Classic Hodgkin Lymphoma Tumor Microenvironment

Classic Hodgkin Lymphoma (cHL) is a tumor composed of rare malignant Hodgkin and Reed-Sternberg (HRS) cells nested within a T-cell rich inflammatory immune infiltrate. cHL is associated with Epstein-Barr Virus (EBV) in 25% of cases. The specific contributions of EBV to the pathogenesis of cHL remain largely unknown, in part due to technical barriers in dissecting the tumor microenvironment (TME) in high detail. Herein, we applied multiplexed ion beam imaging (MIBI) spatial pro-teomics on 6 EBV-positive and 14 EBV-negative cHL samples. We identify key TME features that distinguish between EBV-positive and EBV-negative cHL, including the relative predominance of memory CD8 T cells and increased T-cell dysfunction as a function of spatial proximity to HRS cells. Building upon a larger multi-institutional cohort of 22 EBV-positive and 24 EBV-negative cHL samples, we orthogonally validated our findings through a spatial multi-omics approach, coupling whole transcriptome capture with antibody-defined cell types for tu-mor and T-cell populations within the cHL TME. We delineate contrasting transcriptomic immunological signatures between EBV-positive and EBV-negative cases that differently impact HRS cell proliferation, tumor-immune interactions, and mecha-nisms of T-cell dysregulation and dysfunction. Our multi-modal framework enabled a comprehensive dissection of EBV-linked reorganization and immune evasion within the cHL TME, and highlighted the need to elucidate the cellular and molecular fac-tors of virus-associated tumors, with potential for targeted therapeutic strategies.

A Spatial Multi-Modal Dissection of Host-Microbiome Interactions within the Colitis Tissue Microenvironment

The intricate and dynamic interactions between the host immune system and its microbiome constituents undergo dynamic shifts in response to perturbations to the intestinal tissue environment. Our ability to study these events on the systems level is significantly limited by in situ approaches capable of generating simultaneous insights from both host and microbial communities. Here, we introduce Microbiome Cartography (MicroCart), a framework for simultaneous in situ probing of host features and its microbiome across multiple spatial modalities. We demonstrate MicroCart by comprehensively investigating the alterations in both gut host and microbiome components in a murine model of colitis by coupling MicroCart with spatial proteomics, transcriptomics, and glycomics platforms. Our findings reveal a global but systematic transformation in tissue immune responses, encompassing tissue-level remodeling in response to host immune and epithelial cell state perturbations, and bacterial population shifts, localized inflammatory responses, and metabolic process alterations during colitis. MicroCart enables a deep investigation of the intricate interplay between the host tissue and its microbiome with spatial multiomics.

Multiplex protein imaging in tumour biology

Tissue imaging has become much more colourful in the past decade. Advances in both experimental and analytical methods now make it possible to image protein markers in tissue samples in high multiplex. The ability to routinely image 40-50 markers simultaneously, at single-cell or subcellular resolution, has opened up new vistas in the study of tumour biology. Cellular phenotypes, interaction, communication and spatial organization have become amenable to molecular-level analysis, and application to patient cohorts has identified clinically relevant cellular and tissue features in several cancer types. Here, we review the use of multiplex protein imaging methods to study tumour biology, discuss ongoing attempts to combine these approaches with other forms of spatial omics, and highlight challenges in the field.

Molecular analysis of the extracellular microenvironment: from form to function

The extracellular matrix (ECM) proteome represents an important component of the tissue microenvironment that controls chemical flux and induces cell signaling through encoded structure. The analysis of the ECM represents an analytical challenge through high levels of post-translational modifications, protease-resistant structures, and crosslinked, insoluble proteins. This review provides a comprehensive overview of the analytical challenges involved in addressing the complexities of spatially profiling the extracellular matrix proteome. A synopsis of the process of synthesizing the ECM structure, detailing inherent chemical complexity, is included to present the scope of the analytical challenge. Current chromatographic and spatial techniques addressing these challenges are detailed. Capabilities for multimodal multiplexing with cellular populations are discussed with a perspective on developing a holistic view of disease processes that includes both the cellular and extracellular microenvironment.

High-Plex Assessment of Biomarkers in Tumors

The assessment of biomarkers plays a critical role in the diagnosis and treatment of many cancers. Biomarkers not only provide diagnostic, prognostic, or predictive information but also can act as effective targets for new pharmaceutical therapies. As the utility of biomarkers increases, it becomes more important to utilize accurate and efficient methods for biomarker discovery and, ultimately, clinical assessment. High-plex imaging studies, defined here as assessment of 8 or more biomarkers on a single slide, have become the method of choice for biomarker discovery and assessment of biomarker spatial context. In this review, we discuss methods of measuring biomarkers in slide-mounted tissue samples, detail the various high-plex methods that allow for the simultaneous assessment of multiple biomarkers in situ, and describe the impact of high-plex biomarker assessment on the future of anatomic pathology.

Current updates in sarcoma biomarker discovery: emphasis on next-generation sequencing-based methods

Soft tissue sarcomas comprise a heterogeneous group of neoplasms. Although soft tissue malignancies make up only 2% of adult cancers, classification based on histomorphology presents a diagnostic challenge. Characterisation of soft tissue sarcomas by molecular analysis is rapidly evolving to improve diagnostic accuracy and develop targeted therapies. This review highlights the advances in molecular techniques, including current next-generation sequencing-based assays (fusion detection by RNA sequencing, targeted/whole exome sequencing, microRNA profiling), as well as emerging methods (liquid biopsies, DNA methylation profiling, single-cell molecular profiling and next-generation immunohistochemistry) for future clinical applications.

Semi-automated approaches for interrogating spatial heterogeneity of tissue samples

Tissues are spatially orchestrated ecosystems composed of heterogeneous cell populations and non-cellular elements. Tissue components' interactions shape the biological processes that govern homeostasis and disease, thus comprehensive insights into tissues' composition are crucial for understanding their biology. Recently, advancements in the spatial biology field enabled the in-depth analyses of tissue architecture at single-cell resolution, while preserving the structural context. The increasing number of biomarkers analyzed, together with whole tissue imaging, generate datasets approaching several hundreds of gigabytes in size, which are rich sources of valuable knowledge but require investments in infrastructure and resources for extracting quantitative information. The analysis of multiplex whole-tissue images requires extensive training and experience in data analysis. Here, we showcase how a set of open-source tools can allow semi-automated image data extraction to study the spatial composition of tissues with a focus on tumor microenvironment (TME). With the use of Lunaphore COMET platform, we interrogated lung cancer specimens where we examined the expression of 20 biomarkers. Subsequently, the tissue composition was interrogated using an in-house optimized nuclei detection algorithm followed by a newly developed image artifact exclusion approach. Thereafter, the data was processed using several publicly available tools, highlighting the compatibility of COMET-derived data with currently available image analysis frameworks. In summary, we showcased an innovative semi-automated workflow that highlights the ease of adoption of multiplex imaging to explore TME composition at single-cell resolution using a simple slide in, data out approach. Our workflow is easily transferrable to various cohorts of specimens to provide a toolset for spatial cellular dissection of the tissue composition.

Harmonizing the Generation and Pre-publication Stewardship of FAIR Image data

Together with the molecular knowledge of genes and proteins, biological images promise to significantly enhance the scientific understanding of complex cellular systems and to advance predictive and personalized therapeutic products for human health. For this potential to be realized, quality-assured image data must be shared among labs at a global scale to be compared, pooled, and reanalyzed, thus unleashing untold potential beyond the original purpose for which the data was generated. There are two broad sets of requirements to enable image data sharing in the life sciences. One set of requirements is articulated in the companion White Paper entitled "Enabling Global Image Data Sharing in the Life Sciences," which is published in parallel and addresses the need to build the cyberinfrastructure for sharing the digital array data (arXiv:2401.13023 [q-bio.OT], https://doi.org/10.48550/arXiv.2401.13023). In this White Paper, we detail a broad set of requirements, which involves collecting, managing, presenting, and propagating contextual information essential to assess the quality, understand the content, interpret the scientific implications, and reuse image data in the context of the experimental details. We start by providing an overview of the main lessons learned to date through international community activities, which have recently made considerable progress toward generating community standard practices for imaging Quality Control (QC) and metadata. We then provide a clear set of recommendations for amplifying this work. The driving goal is to address remaining challenges, and democratize access to common practices and tools for a spectrum of biomedical researchers, regardless of their expertise, access to resources, and geographical location.

Leveraging spatial omics for the development of precision sarcoma treatments

Sarcomas are rare and heterogeneous cancers that arise from bone or soft tissue, and are the second most prevalent solid cancer in children and adolescents. Owing to the complex nature of pediatric sarcomas, the development of therapeutics for pediatric sarcoma has seen little progress in the past decades. Existing treatments are largely limited to chemotherapy, radiation, and surgery. Limited knowledge of the sarcoma tumor microenvironment (TME) and of well-defined target antigens in the different subtypes necessitates an alternative investigative approach to improve treatments. Recent advances in spatial omics technologies have enabled a more comprehensive study of the TME in multiple cancers. In this opinion article we discuss advances in our understanding of the TME of some cancers enabled by spatial omics technologies, and we explore how these technologies might advance the development of precision treatments for sarcoma, especially pediatric sarcoma.

Mapping Cell Atlases at the Single-Cell Level

Recent advancements in single-cell technologies have led to rapid developments in the construction of cell atlases. These atlases have the potential to provide detailed information about every cell type in different organisms, enabling the characterization of cellular diversity at the single-cell level. Global efforts in developing comprehensive cell atlases have profound implications for both basic research and clinical applications. This review provides a broad overview of the cellular diversity and dynamics across various biological systems. In addition, the incorporation of machine learning techniques into cell atlas analyses opens up exciting prospects for the field of integrative biology.

A spatial cell atlas of neuroblastoma reveals developmental, epigenetic and spatial axis of tumor heterogeneity

Neuroblastoma is a pediatric cancer arising from the developing sympathoadrenal lineage with complex inter- and intra-tumoral heterogeneity. To chart this complexity, we generated a comprehensive cell atlas of 55 neuroblastoma patient tumors, collected from two pediatric cancer institutions, spanning a range of clinical, genetic, and histologic features. Our atlas combines single-cell/nucleus RNA-seq (sc/scRNA-seq), bulk RNA-seq, whole exome sequencing, DNA methylation profiling, spatial transcriptomics, and two spatial proteomic methods. Sc/snRNA-seq revealed three malignant cell states with features of sympathoadrenal lineage development. All of the neuroblastomas had malignant cells that resembled sympathoblasts and the more differentiated adrenergic cells. A subset of tumors had malignant cells in a mesenchymal cell state with molecular features of Schwann cell precursors. DNA methylation profiles defined four groupings of patients, which differ in the degree of malignant cell heterogeneity and clinical outcomes. Using spatial proteomics, we found that neuroblastomas are spatially compartmentalized, with malignant tumor cells sequestered away from immune cells. Finally, we identify spatially restricted signaling patterns in immune cells from spatial transcriptomics. To facilitate the visualization and analysis of our atlas as a resource for further research in neuroblastoma, single cell, and spatial-omics, all data are shared through the Human Tumor Atlas Network Data Commons at www.humantumoratlas.org.

Mapping the single cell spatial immune landscapes of the melanoma microenvironment

Melanoma is a highly immunogenic malignancy with an elevated mutational burden, diffuse lymphocytic infiltration, and one of the highest response rates to immune checkpoint inhibitors (ICIs). However, over half of all late-stage patients treated with ICIs will either not respond or develop progressive disease. Spatial imaging technologies are being increasingly used to study the melanoma tumor microenvironment (TME). The goal of such studies is to understand the complex interplay between the stroma, melanoma cells, and immune cell-types as well as their association with treatment response. Investigators seeking a better understanding of the role of cell location within the TME and the importance of spatial expression of biomarkers are increasingly turning to highly multiplexed imaging approaches to more accurately measure immune infiltration as well as to quantify receptor-ligand interactions (such as PD-1 and PD-L1) and cell-cell contacts. CyTOF-IMC (Cytometry by Time of Flight - Imaging Mass Cytometry) has enabled high-dimensional profiling of melanomas, allowing researchers to identify complex cellular subpopulations and immune cell interactions with unprecedented resolution. Other spatial imaging technologies, such as multiplexed immunofluorescence and spatial transcriptomics, have revealed distinct patterns of immune cell infiltration, highlighting the importance of spatial relationships, and their impact in modulating immunotherapy responses. Overall, spatial imaging technologies are just beginning to transform our understanding of melanoma biology, providing new avenues for biomarker discovery and therapeutic development. These technologies hold great promise for advancing personalized medicine to improve patient outcomes in melanoma and other solid malignancies.

cytoviewer: an R/Bioconductor package for interactive visualization and exploration of highly multiplexed imaging data

BACKGROUND

Highly multiplexed imaging enables single-cell-resolved detection of numerous biological molecules in their spatial tissue context. Interactive visualization of multiplexed imaging data is crucial at any step of data analysis to facilitate quality control and the spatial exploration of single cell features. However, tools for interactive visualization of multiplexed imaging data are not available in the statistical programming language R.

RESULTS

Here, we describe cytoviewer, an R/Bioconductor package for interactive visualization and exploration of multi-channel images and segmentation masks. The cytoviewer package supports flexible generation of image composites, allows side-by-side visualization of single channels, and facilitates the spatial visualization of single-cell data in the form of segmentation masks. As such, cytoviewer improves image and segmentation quality control, the visualization of cell phenotyping results and qualitative validation of hypothesis at any step of data analysis. The package operates on standard data classes of the Bioconductor project and therefore integrates with an extensive framework for single-cell and image analysis. The graphical user interface allows intuitive navigation and little coding experience is required to use the package. We showcase the functionality and biological application of cytoviewer by analysis of an imaging mass cytometry dataset acquired from cancer samples.

CONCLUSIONS

The cytoviewer package offers a rich set of features for highly multiplexed imaging data visualization in R that seamlessly integrates with the workflow for image and single-cell data analysis. It can be installed from Bioconductor via https://www.bioconductor.org/packages/release/bioc/html/cytoviewer.html . The development version and further instructions can be found on GitHub at https://github.com/BodenmillerGroup/cytoviewer .

ESQmodel: biologically informed evaluation of 2-D cell segmentation quality in multiplexed tissue images

MOTIVATION

Single cell segmentation is critical in the processing of spatial omics data to accurately perform cell type identification and analyze spatial expression patterns. Segmentation methods often rely on semi-supervised annotation or labeled training data which are highly dependent on user expertise. To ensure the quality of segmentation, current evaluation strategies quantify accuracy by assessing cellular masks or through iterative inspection by pathologists. While these strategies each address either the statistical or biological aspects of segmentation, there lacks a unified approach to evaluating segmentation accuracy.

RESULTS

In this article, we present ESQmodel, a Bayesian probabilistic method to evaluate single cell segmentation using expression data. By using the extracted cellular data from segmentation and a prior belief of cellular composition as input, ESQmodel computes per cell entropy to assess segmentation quality by how consistent cellular expression profiles match with cell type expectations.

AVAILABILITY AND IMPLEMENTATION

Source code is available on Github at: https://github.com/Roth-Lab/ESQmodel.

MAPS: pathologist-level cell type annotation from tissue images through machine learning

Highly multiplexed protein imaging is emerging as a potent technique for analyzing protein distribution within cells and tissues in their native context. However, existing cell annotation methods utilizing high-plex spatial proteomics data are resource intensive and necessitate iterative expert input, thereby constraining their scalability and practicality for extensive datasets. We introduce MAPS (Machine learning for Analysis of Proteomics in Spatial biology), a machine learning approach facilitating rapid and precise cell type identification with human-level accuracy from spatial proteomics data. Validated on multiple in-house and publicly available MIBI and CODEX datasets, MAPS outperforms current annotation techniques in terms of speed and accuracy, achieving pathologist-level precision even for typically challenging cell types, including tumor cells of immune origin. By democratizing rapidly deployable and scalable machine learning annotation, MAPS holds significant potential to expedite advances in tissue biology and disease comprehension.

Spatial omics: An innovative frontier in aging research

Disentangling the impact of aging on health and disease has become critical as population aging progresses rapidly. Studying aging at the molecular level is complicated by the diverse aging profiles and dynamics. However, the examination of cellular states within aging tissues in situ is hampered by the lack of high-resolution spatial data. Emerging spatial omics technologies facilitate molecular and spatial analysis of tissues, providing direct access to precise information on various functional regions and serving as a favorable tool for unraveling the heterogeneity of aging. In this review, we summarize the recent advances in spatial omics application in multi-organ aging research, which has enhanced the understanding of aging mechanisms from multiple standpoints. We also discuss the main challenges in spatial omics research to date, the opportunities for further developing the technology, and the potential applications of spatial omics in aging and aging-related diseases.

Strategies for studying immune and non-immune human and canine mammary gland cancer tumour infiltrate

The tumour microenvironment (TME) is usually defined as a cell environment associated with tumours or cancerous stem cells where conditions are established affecting tumour development and progression through malignant cell interaction with non-malignant cells. The TME is made up of endothelial, immune and non-immune cells, extracellular matrix (ECM) components and signalling molecules acting specifically on tumour and non-tumour cells. Breast cancer (BC) is the commonest malignant neoplasm worldwide and the main cause of mortality in women globally; advances regarding BC study and understanding it are relevant for acquiring novel, personalised therapeutic tools. Studying canine mammary gland tumours (CMGT) is one of the most relevant options for understanding BC using animal models as they share common epidemiological, clinical, pathological, biological, environmental, genetic and molecular characteristics with human BC. In-depth, detailed investigation regarding knowledge of human BC-related TME and in its canine model is considered extremely relevant for understanding changes in TME composition during tumour development. This review addresses important aspects concerned with different methods used for studying BC- and CMGT-related TME that are important for developing new and more effective therapeutic strategies for attacking a tumour during specific evolutionary stages.

OPTIMAL: An OPTimized Imaging Mass cytometry AnaLysis framework for benchmarking segmentation and data exploration

Analysis of imaging mass cytometry (IMC) data and other low-resolution multiplexed tissue imaging technologies is often confounded by poor single-cell segmentation and suboptimal approaches for data visualization and exploration. This can lead to inaccurate identification of cell phenotypes, states, or spatial relationships compared to reference data from single-cell suspension technologies. To this end we have developed the "OPTimized Imaging Mass cytometry AnaLysis (OPTIMAL)" framework to benchmark any approaches for cell segmentation, parameter transformation, batch effect correction, data visualization/clustering, and spatial neighborhood analysis. Using a panel of 27 metal-tagged antibodies recognizing well-characterized phenotypic and functional markers to stain the same Formalin-Fixed Paraffin Embedded (FFPE) human tonsil sample tissue microarray over 12 temporally distinct batches we tested several cell segmentation models, a range of different arcsinh cofactor parameter transformation values, 5 different dimensionality reduction algorithms, and 2 clustering methods. Finally, we assessed the optimal approach for performing neighborhood analysis. We found that single-cell segmentation was improved by the use of an Ilastik-derived probability map but that issues with poor segmentation were only really evident after clustering and cell type/state identification and not always evident when using "classical" bivariate data display techniques. The optimal arcsinh cofactor for parameter transformation was 1 as it maximized the statistical separation between negative and positive signal distributions and a simple Z-score normalization step after arcsinh transformation eliminated batch effects. Of the five different dimensionality reduction approaches tested, PacMap gave the best data structure with FLOWSOM clustering out-performing phenograph in terms of cell type identification. We also found that neighborhood analysis was influenced by the method used for finding neighboring cells with a "disc" pixel expansion outperforming a "bounding box" approach combined with the need for filtering objects based on size and image-edge location. Importantly, OPTIMAL can be used to assess and integrate with any existing approach to IMC data analysis and, as it creates .FCS files from the segmentation output and allows for single-cell exploration to be conducted using a wide variety of accessible software and algorithms familiar to conventional flow cytometrists.

Spatial modelling of the tumor microenvironment from multiplex immunofluorescence images: methods and applications

Spatial modelling methods have gained prominence with developments in high throughput imaging platforms. Multiplex immunofluorescence (mIF) provides the scope to examine interactions between tumor and immune compartment at single cell resolution using a panel of antibodies that can be chosen based on the cancer type or the clinical interest of the study. The markers can be used to identify the phenotypes and to examine cellular interactions at global and local scales. Several translational studies rely on key understanding of the tumor microenvironment (TME) to identify drivers of immune response in immunotherapy based clinical trials. To improve the success of ongoing trials, a number of retrospective approaches can be adopted to understand differences in response, recurrence and progression by examining the patient's TME from tissue samples obtained at baseline and at various time points along the treatment. The multiplex immunofluorescence (mIF) technique provides insight on patient specific cell populations and their relative spatial distribution as qualitative measures of a favorable treatment outcome. Spatial analysis of these images provides an understanding of the intratumoral heterogeneity and clustering among cell populations in the TME. A number of mathematical models, which establish clustering as a measure of deviation from complete spatial randomness, can be applied to the mIF images represented as spatial point patterns. These mathematical models, developed for landscape ecology and geographic information studies, can be applied to the TME after careful consideration of the tumor type (cold vs. hot) and the tumor immune landscape. The spatial modelling of mIF images can show observable engagement of T cells expressing immune checkpoint molecules and this can then be correlated with single-cell RNA sequencing data.

Spatial Proteomics toward Subcellular Resolution by Coupling Deep Ultraviolet Laser Ablation with Nanodroplet Sample Preparation

Multiplexed molecular profiling of tissue microenvironments, or spatial omics, can provide critical insights into cellular functions and disease pathology. The coupling of laser microdissection with mass spectrometry-based proteomics has enabled deep and unbiased mapping of >1000 proteins. However, the throughput of laser microdissection is often limited due to tedious two-step procedures, sequential laser cutting, and sample collection. The two-step procedure also hinders the further improvement of spatial resolution to <10 μm as needed for subcellular proteomics. Herein, we developed a high-throughput and high-resolution spatial proteomics platform by seamlessly coupling deep ultraviolet (DUV) laser ablation (LA) with nanoPOTS (Nanodroplet Processing in One pot for Trace Samples)-based sample preparation. We demonstrated the DUV-LA system can quickly isolate and collect tissue samples at a throughput of ∼30 spots/min and a spatial resolution down to 2 μm from a 10 μm thick human pancreas tissue section. To improve sample recovery, we developed a proximity aerosol collection approach by placing DMSO droplets close to LA spots. We demonstrated the DUV-LA-nanoPOTS platform can detect an average of 1312, 1533, and 1966 proteins from ablation spots with diameters of 7, 13, and 19 μm, respectively. In a proof-of-concept study, we isolated and profiled two distinct subcellular regions of the pancreas tissue revealed by hematoxylin and eosin (H&E) staining. Quantitative proteomics revealed proteins specifically enriched to subcellular compartments.