The Spatial Landscape of Progression and Immunoediting in Primary Melanoma at Single-Cell Resolution

ICG-labeled PD-L1-antagonistic affibody dimer for tumor imaging and enhancement of tumor photothermal-immunotherapy

In clinical practice, tumor-targeting diagnosis and immunotherapy against programmed death ligand 1 (PD-L1) have a significant impact. In this research, a PD-L1-antagonistic affibody dimer (ZPD-L1) was successfully prepared through Escherichia coli expression system, and conjugated with the photosensitizer of ICG via N-hydroxysuccinimide (NHS) ester to develop a novel tumor-targeting agent (ICG-ZPD-L1) for both tumor imaging diagnosis and photothermal-immunotherapy simultaneously. In vitro, ZPD-L1 could specifically bind to PD-L1-positive LLC and MC38 tumor cells, and ICG-ZPD-L1-mediated photothermal therapy (PTT) also showed excellent phototoxicity to these tumor cells. In vivo, ICG-ZPD-L1 selectively enriched into the PD-L1-positive MC38 tumor tissues, and the high-contrast optical imaging of tumors was obtained. ICG-ZPD-L1-mediated PTT exhibited a potent anti-tumor effect in vivo due to its remarkable photothermal properties. Furthermore, ICG-ZPD-L1-mediated PTT significantly induced the immunogenic cell death (ICD) of primary tumors, promoted maturation of dendritic cells (DCs), up-regulated anti-tumor immune response, enhanced immunotherapy, and superiorly inhibited the growth of metastatic tumors. In addition, ICG-ZPD-L1 showed favorable biosafety throughout the brief duration of treatment. In summary, these results suggest that ICG-ZPD-L1 is a multifunctional tumor-targeting drug integrating tumor imaging diagnosis and photothermal-immunotherapy, and has great guiding significance for the diagnosis and treatment of clinical PD-L1-positive tumor patients.

Integrating spatial and single-cell transcriptomics reveals tumor heterogeneity and intercellular networks in colorectal cancer

Single cell RNA sequencing (scRNA-seq), a powerful tool for studying the tumor microenvironment (TME), does not preserve/provide spatial information on tissue morphology and cellular interactions. To understand the crosstalk between diverse cellular components in proximity in the TME, we performed scRNA-seq coupled with spatial transcriptomic (ST) assay to profile 41,700 cells from three colorectal cancer (CRC) tumor-normal-blood pairs. Standalone scRNA-seq analyses revealed eight major cell populations, including B cells, T cells, Monocytes, NK cells, Epithelial cells, Fibroblasts, Mast cells, Endothelial cells. After the identification of malignant cells from epithelial cells, we observed seven subtypes of malignant cells that reflect heterogeneous status in tumor, including tumor_CAV1, tumor_ATF3_JUN | FOS, tumor_ZEB2, tumor_VIM, tumor_WSB1, tumor_LXN, and tumor_PGM1. By transferring the cellular annotations obtained by scRNA-seq to ST spots, we annotated four regions in a cryosection from CRC patients, including tumor, stroma, immune infiltration, and colon epithelium regions. Furthermore, we observed intensive intercellular interactions between stroma and tumor regions which were extremely proximal in the cryosection. In particular, one pair of ligands and receptors (C5AR1 and RPS19) was inferred to play key roles in the crosstalk of stroma and tumor regions. For the tumor region, a typical feature of TMSB4X-high expression was identified, which could be a potential marker of CRC. The stroma region was found to be characterized by VIM-high expression, suggesting it fostered a stromal niche in the TME. Collectively, single cell and spatial analysis in our study reveal the tumor heterogeneity and molecular interactions in CRC TME, which provides insights into the mechanisms underlying CRC progression and may contribute to the development of anticancer therapies targeting on non-tumor components, such as the extracellular matrix (ECM) in CRC. The typical genes we identified may facilitate to new molecular subtypes of CRC.

A general algorithm for consensus 3D cell segmentation from 2D segmented stacks

Cell segmentation is the fundamental task. Only by segmenting, can we define the quantitative spatial unit for collecting measurements to draw biological conclusions. Deep learning has revolutionized 2D cell segmentation, enabling generalized solutions across cell types and imaging modalities. This has been driven by the ease of scaling up image acquisition, annotation and computation. However 3D cell segmentation, which requires dense annotation of 2D slices still poses significant challenges. Labelling every cell in every 2D slice is prohibitive. Moreover it is ambiguous, necessitating cross-referencing with other orthoviews. Lastly, there is limited ability to unambiguously record and visualize 1000's of annotated cells. Here we develop a theory and toolbox, u-Segment3D for 2D-to-3D segmentation, compatible with any 2D segmentation method. Given optimal 2D segmentations, u-Segment3D generates the optimal 3D segmentation without data training, as demonstrated on 11 real life datasets, >70,000 cells, spanning single cells, cell aggregates and tissue.

Early Immune Remodeling Steers Clinical Response to First-Line Chemoimmunotherapy in Advanced Gastric Cancer

Adding anti-programmed cell death protein 1 (anti-PD-1) to 5-fluorouracil (5-FU)/platinum improves survival in some advanced gastroesophageal adenocarcinomas (GEA). To understand the effects of chemotherapy and immunotherapy, we conducted a phase II first-line trial (n = 47) sequentially adding pembrolizumab to 5-FU/platinum in advanced GEA. Using serial biopsy of the primary tumor at baseline, after one cycle of 5-FU/platinum, and after the addition of pembrolizumab, we transcriptionally profiled 358,067 single cells to identify evolving multicellular tumor microenvironment (TME) networks. Chemotherapy induced early on-treatment multicellular hubs with tumor-reactive T-cell and M1-like macrophage interactions in slow progressors. Faster progression featured increased MUC5A and MSLN containing treatment resistance programs in tumor cells and M2-like macrophages with immunosuppressive stromal interactions. After pembrolizumab, we observed increased CD8 T-cell infiltration and development of an immunity hub involving tumor-reactive CXCL13 T-cell program and epithelial interferon-stimulated gene programs. Strategies to drive increases in antitumor immune hub formation could expand the portion of patients benefiting from anti-PD-1 approaches.

SIGNIFICANCE

The benefit of 5-FU/platinum with anti-PD-1 in first-line advanced gastric cancer is limited to patient subgroups. Using a trial with sequential anti-PD-1, we show coordinated induction of multicellular TME hubs informs the ability of anti-PD-1 to potentiate T cell-driven responses. Differential TME hub development highlights features that underlie clinical outcomes. This article is featured in Selected Articles from This Issue, p. 695.

psudo: Exploring Multi-Channel Biomedical Image Data with Spatially and Perceptually Optimized Pseudocoloring

Over the past century, multichannel fluorescence imaging has been pivotal in myriad scientific breakthroughs by enabling the spatial visualization of proteins within a biological sample. With the shift to digital methods and visualization software, experts can now flexibly pseudocolor and combine image channels, each corresponding to a different protein, to explore their spatial relationships. We thus propose psudo, an interactive system that allows users to create optimal color palettes for multichannel spatial data. In psudo, a novel optimization method generates palettes that maximize the perceptual differences between channels while mitigating confusing color blending in overlapping channels. We integrate this method into a system that allows users to explore multi-channel image data and compare and evaluate color palettes for their data. An interactive lensing approach provides on-demand feedback on channel overlap and a color confusion metric while giving context to the underlying channel values. Color palettes can be applied globally or, using the lens, to local regions of interest. We evaluate our palette optimization approach using three graphical perception tasks in a crowdsourced user study with 150 participants, showing that users are more accurate at discerning and comparing the underlying data using our approach. Additionally, we showcase psudo in a case study exploring the complex immune responses in cancer tissue data with a biologist.

Understanding the Tumor Microenvironment in Melanoma Patients with In-Transit Metastases and Its Impacts on Immune Checkpoint Immunotherapy Responses

Melanoma is the leading cause of global skin cancer-related death and currently ranks as the third most commonly diagnosed cancer in Australia. Melanoma patients with in-transit metastases (ITM), a type of locoregional metastasis located close to the primary tumor site, exhibit a high likelihood of further disease progression and poor survival outcomes. Immunotherapies, particularly immune checkpoint inhibitors (ICI), have demonstrated remarkable efficacy in ITM patients with reduced occurrence of further metastases and prolonged survival. The major challenge of immunotherapeutic efficacy lies in the limited understanding of melanoma and ITM biology, hindering our ability to identify patients who likely respond to ICIs effectively. In this review, we provided an overview of melanoma and ITM disease. We outlined the key ICI therapies and the critical immune features associated with therapy response or resistance. Lastly, we dissected the underlying biological components, including the cellular compositions and their communication networks within the tumor compartment, to enhance our understanding of the interactions between immunotherapy and melanoma, providing insights for future investigation and the development of drug targets and predictive biomarkers.

Tumor-infiltrating lymphocytes as a prognostic and predictive factor for Melanoma

INTRODUCTION

Tumor-infiltrating lymphocytes (TILs) have been investigated as prognostic factors in melanoma. Recent advancements in assessing the tumor microenvironment in the setting of more widespread use of immune checkpoint blockade have reignited interest in identifying predictive biomarkers. This review examines the function and significance of TILs in cutaneous melanoma, evaluating their potential as prognostic and predictive markers.

AREAS COVERED

A literature search was conducted on papers covering tumor infiltrating lymphocytes in cutaneous melanoma available online in PubMed and Web of Science from inception to 1 December 2023, supplemented by citation searching. This article encompasses the assessment of TILs, the role of TILs in the immune microenvironment, TILs as a prognostic factor, TILs as a predictive factor for immunotherapy response, and clinical applications of TILs in the treatment of cutaneous melanoma.

EXPERT OPINION

Tumor-infiltrating lymphocytes play a heterogeneous role in cutaneous melanoma. While they have historically been associated with improved survival, their status as independent prognostic or predictive factors remains uncertain. Novel methods of TIL assessment, such as determination of TIL subtypes and molecular signaling, demonstrate potential for predicting therapeutic response. Further, while their clinical utility in risk-stratification in melanoma treatment shows promise, a lack of consensus data hinders standardized application.

Integrative analysis of bulk and single-cell RNA-seq reveals the molecular characterization of the immune microenvironment and oxidative stress signature in melanoma

Background

The immune microenvironment and oxidative stress of melanoma show significant heterogeneity, which affects tumor growth, invasion and treatment response. Single-cell and bulk RNA-seq data were used to explore the heterogeneity of the immune microenvironment and oxidative stress of melanoma.

Methods

The R package Seurat facilitated the analysis of the single-cell dataset, while Harmony, another R package, was employed for batch effect correction. Cell types were classified using Uniform Manifold Approximation and Projection (UMAP). The Secreted Signaling algorithm from CellChatDB.human was applied to elucidate cell-to-cell communication patterns within the single-cell data. Consensus clustering analysis for the skin cutaneous melanoma (SKCM) samples was executed with the R package ConsensusClusterPlus. To quantify immune infiltrating cells, we utilized CIBERSORT, ESTIMATE, and TIMERxCell algorithms provided by the R package Immuno-Oncology Biological Research (IOBR). Single nucleotide variant (SNV) analysis was conducted using Maftools, an R package specifically designed for this purpose. Subsequently, the expression levels of PXDN and PAPSS2 genes were assessed in melanoma tissues compared to adjacent normal tissues. Furthermore, in vitro experiments were conducted to evaluate the proliferation and reactive oxygen species expression in melanoma cells following transfection with siRNA targeting PXDN and PAPSS2.

Results

Malignant tumor cell populations were reclassified based on a comprehensive single-cell dataset analysis, which yielded six distinct tumor subsets. The specific marker genes identified for these subgroups were then used to interrogate the Cancer Genome Atlas Skin Cutaneous Melanoma (TCGA-SKCM) cohort, derived from bulk RNA sequencing data, resulting in the delineation of two immune molecular subtypes. Notably, patients within the cluster2 (C2) subtype exhibited a significantly more favorable prognosis compared to those in the cluster1 (C1) subtype. An alignment of immune characteristics was observed between the C2 subtype and unique immune functional tumor cell subsets. Genes differentially expressed across these subtypes were subsequently leveraged to construct a predictive risk model. In vitro investigations further revealed elevated expression levels of PXDN and PAPSS2 in melanoma tissue samples. Functional assays indicated that modulation of PXDN and PAPSS2 expression could influence the production of reactive oxygen species (ROS) and the proliferative capacity of melanoma cells.

Conclusion

The constructed six-gene signature can be used as an immune response and an oxidative stress marker to guide the clinical diagnosis and treatment of melanoma.

Multiplexed 3D Analysis of Immune States and Niches in Human Tissue

Tissue homeostasis and the emergence of disease are controlled by changes in the proportions of resident and recruited cells, their organization into cellular neighbourhoods, and their interactions with acellular tissue components. Highly multiplexed tissue profiling (spatial omics)1 makes it possible to study this microenvironment in situ, usually in 4-5 micron thick sections (the standard histopathology format)2. Microscopy-based tissue profiling is commonly performed at a resolution sufficient to determine cell types but not to detect subtle morphological features associated with cytoskeletal reorganisation, juxtracrine signalling, or membrane trafficking3. Here we describe a high-resolution 3D imaging approach able to characterize a wide variety of organelles and structures at sub-micron scale while simultaneously quantifying millimetre-scale spatial features. This approach combines cyclic immunofluorescence (CyCIF) imaging4 of over 50 markers with confocal microscopy of archival human tissue thick enough (30-40 microns) to fully encompass two or more layers of intact cells. 3D imaging of entire cell volumes substantially improves the accuracy of cell phenotyping and allows cell proximity to be scored using plasma membrane apposition, not just nuclear position. In pre-invasive melanoma in situ5, precise phenotyping shows that adjacent melanocytic cells are plastic in state and participate in tightly localised niches of interferon signalling near sites of initial invasion into the underlying dermis. In this and metastatic melanoma, mature and precursor T cells engage in an unexpectedly diverse array of juxtracrine and membrane-membrane interactions as well as looser "neighbourhood" associations6 whose morphologies reveal functional states. These data provide new insight into the transitions occurring during early tumour formation and immunoediting and demonstrate the potential for phenotyping of tissues at a level of detail previously restricted to cultured cells and organoids.

SpatialCells: automated profiling of tumor microenvironments with spatially resolved multiplexed single-cell data

Cancer is a complex cellular ecosystem where malignant cells coexist and interact with immune, stromal and other cells within the tumor microenvironment (TME). Recent technological advancements in spatially resolved multiplexed imaging at single-cell resolution have led to the generation of large-scale and high-dimensional datasets from biological specimens. This underscores the necessity for automated methodologies that can effectively characterize molecular, cellular and spatial properties of TMEs for various malignancies. This study introduces SpatialCells, an open-source software package designed for region-based exploratory analysis and comprehensive characterization of TMEs using multiplexed single-cell data. The source code and tutorials are available at https://semenovlab.github.io/SpatialCells. SpatialCells efficiently streamlines the automated extraction of features from multiplexed single-cell data and can process samples containing millions of cells. Thus, SpatialCells facilitates subsequent association analyses and machine learning predictions, making it an essential tool in advancing our understanding of tumor growth, invasion and metastasis.

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.

Chemotherapy induces myeloid-driven spatial T-cell exhaustion in ovarian cancer

To uncover the intricate, chemotherapy-induced spatiotemporal remodeling of the tumor microenvironment, we conducted integrative spatial and molecular characterization of 97 high-grade serous ovarian cancer (HGSC) samples collected before and after chemotherapy. Using single-cell and spatial analyses, we identify increasingly versatile immune cell states, which form spatiotemporally dynamic microcommunities at the tumor-stroma interface. We demonstrate that chemotherapy triggers spatial redistribution and exhaustion of CD8+ T cells due to prolonged antigen presentation by macrophages, both within interconnected myeloid networks termed "Myelonets" and at the tumor stroma interface. Single-cell and spatial transcriptomics identifies prominent TIGIT-NECTIN2 ligand-receptor interactions induced by chemotherapy. Using a functional patient-derived immuno-oncology platform, we show that CD8+T-cell activity can be boosted by combining immune checkpoint blockade with chemotherapy. Our discovery of chemotherapy-induced myeloid-driven spatial T-cell exhaustion paves the way for novel immunotherapeutic strategies to unleash CD8+ T-cell-mediated anti-tumor immunity in HGSC.

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.

Macrophage barrier in the tumor microenvironment and potential clinical applications

The tumor microenvironment (TME) constitutes a complex microenvironment comprising a diverse array of immune cells and stromal components. Within this intricate context, tumor-associated macrophages (TAMs) exhibit notable spatial heterogeneity. This heterogeneity contributes to various facets of tumor behavior, including immune response modulation, angiogenesis, tissue remodeling, and metastatic potential. This review summarizes the spatial distribution of macrophages in both the physiological environment and the TME. Moreover, this paper explores the intricate interactions between TAMs and diverse immune cell populations (T cells, dendritic cells, neutrophils, natural killer cells, and other immune cells) within the TME. These bidirectional exchanges form a complex network of immune interactions that influence tumor immune surveillance and evasion strategies. Investigating TAM heterogeneity and its intricate interactions with different immune cell populations offers potential avenues for therapeutic interventions. Additionally, this paper discusses therapeutic strategies targeting macrophages, aiming to uncover novel approaches for immunotherapy. Video Abstract.

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.

Residency Octree: A Hybrid Approach for Scalable Web-Based Multi-Volume Rendering

We present a hybrid multi-volume rendering approach based on a novel Residency Octree that combines the advantages of out-of-core volume rendering using page tables with those of standard octrees. Octree approaches work by performing hierarchical tree traversal. However, in octree volume rendering, tree traversal and the selection of data resolution are intrinsically coupled. This makes fine-grained empty-space skipping costly. Page tables, on the other hand, allow access to any cached brick from any resolution. However, they do not offer a clear and efficient strategy for substituting missing high-resolution data with lower-resolution data. We enable flexible mixed-resolution out-of-core multi-volume rendering by decoupling the cache residency of multi-resolution data from a resolution-independent spatial subdivision determined by the tree. Instead of one-to-one node-to-brick correspondences, each residency octree node is mapped to a set of bricks from different resolution levels. This makes it possible to efficiently and adaptively choose and mix resolutions, adapt sampling rates, and compensate for cache misses. At the same time, residency octrees support fine-grained empty-space skipping, independent of the data subdivision used for caching. Finally, to facilitate collaboration and outreach, and to eliminate local data storage, our implementation is a web-based, pure client-side renderer using WebGPU and WebAssembly. Our method is faster than prior approaches and efficient for many data channels with a flexible and adaptive choice of data resolution.

Laser capture microdissection provides a novel molecular profile of human primary cutaneous melanoma

Melanoma accounts for the majority of skin cancer-related mortality, highlighting the need to better understand melanoma initiation and progression. In-depth molecular analysis of neoplastic melanocytes in whole tissue biopsies may be diluted by inflammatory infiltration, which may obscure gene signatures specific to neoplastic cells. Thus, a method is needed to precisely uncover molecular changes specific to tumor cells from a limited sample of primary melanomas. Here, we performed laser capture microdissection (LCM) and gene expression profiling of patient-derived frozen sections of pigmented lesions and primary cutaneous melanoma. Compared to bulk tissue analysis, analysis of LCM-derived samples identified 9528 additional differentially expressed genes (DEGs) including melanocyte-specific genes like PMEL and TYR, with enriched of pathways related to cell proliferation. LCM methodology also identified potentially targetable kinases specific to melanoma cells that were not detected by bulk tissue analysis. Taken together, our data demonstrate that there are marked differences in gene expression profiles depending on the method of sample isolation. We found that LCM captured higher expression of melanoma-related genes while whole tissue biopsy identified a wider range of inflammatory markers. Taken together, our data demonstrate that LCM is a valid approach to identify melanoma-specific changes using a relatively small amount of primary patient-derived melanoma sample.

Spatial Patterning Analysis of Cellular Ensembles (SPACE) enables statistically robust discovery of complex spatial organization at the cell and tissue level

Spatial patterns of cells and other entities drive both physiologic and pathologic processes within tissues. While many imaging and transcriptomic methods document tissue organization, discerning these patterns is challenging, especially when they involve multiple entities in complex arrangements. To address this challenge, we present Spatial Patterning Analysis of Cellular Ensembles (SPACE), an R package for analysis of high-plex tissue images generated using any collection modality. Unlike existing platforms, SPACE detects context-dependent associations, quantitative gradients and orientations, and other organizational complexities. Via a robust information theoretic framework, SPACE explores all possible ensembles - single entities, pairs, triplets, and so on - and ranks the strongest patterns of tissue organization. Using lymph node images for which ground truth has been defined, we validate SPACE and demonstrate its advantages. We then use SPACE to reanalyze a public dataset of human tuberculosis granulomas, verifying known patterns and discovering new patterns with possible insights into disease progression.

Cell Spotter (CSPOT): A machine-learning approach to automated cell spotting and quantification of highly multiplexed tissue images

Highly multiplexed tissue imaging and in situ spatial profiling aim to extract single-cell data from specimens containing closely packed cells of diverse morphology. This is challenging due to the difficulty of accurately assigning boundaries between cells (segmentation) and then generating per-cell staining intensities. Existing methods use gating to convert per-cell intensity data to positive and negative scores; this is a common approach in flow cytometry, but one that is problematic in imaging. In contrast, human experts identify cells in crowded environments using morphological, neighborhood, and intensity information. Here we describe a computational approach (Cell Spotter or CSPOT) that uses supervised machine learning in combination with classical segmentation to perform automated cell type calling. CSPOT is robust to artifacts that commonly afflict tissue imaging and can replace conventional gating. The end-to-end Python implementation of CSPOT can be integrated into cloud-based image processing pipelines to substantially improve the speed, accuracy, and reproducibility of single-cell spatial data.

SpatialCells: Automated Profiling of Tumor Microenvironments with Spatially Resolved Multiplexed Single-Cell Data

Background

Cancer is a complex cellular ecosystem where malignant cells coexist and interact with immune, stromal, and other cells within the tumor microenvironment. Recent technological advancements in spatially resolved multiplexed imaging at single-cell resolution have led to the generation of large-scale and high-dimensional datasets from biological specimens. This underscores the necessity for automated methodologies that can effectively characterize the molecular, cellular, and spatial properties of tumor microenvironments for various malignancies.

Results

This study introduces SpatialCells, an open-source software package designed for region-based exploratory analysis and comprehensive characterization of tumor microenvironments using multiplexed single-cell data.

Conclusions

SpatialCells efficiently streamlines the automated extraction of features from multiplexed single-cell data and can process samples containing millions of cells. Thus, SpatialCells facilitates subsequent association analyses and machine learning predictions, making it an essential tool in advancing our understanding of tumor growth, invasion, and metastasis.

2D and 3D multiplexed subcellular profiling of nuclear instability in human cancer

Nuclear atypia, including altered nuclear size, contour, and chromatin organization, is ubiquitous in cancer cells. Atypical primary nuclei and micronuclei can rupture during interphase; however, the frequency, causes, and consequences of nuclear rupture are unknown in most cancers. We demonstrate that nuclear envelope rupture is surprisingly common in many human cancers, particularly glioblastoma. Using highly-multiplexed 2D and super-resolution 3D-imaging of glioblastoma tissues and patient-derived xenografts and cells, we link primary nuclear rupture with reduced lamin A/C and micronuclear rupture with reduced lamin B1. Moreover, ruptured glioblastoma cells activate cGAS-STING-signaling involved in innate immunity. We observe that local patterning of cell states influences tumor spatial organization and is linked to both lamin expression and rupture frequency, with neural-progenitor-cell-like states exhibiting the lowest lamin A/C levels and greatest susceptibility to primary nuclear rupture. Our study reveals that nuclear instability is a core feature of cancer, and links nuclear integrity, cell state, and immune signaling.

Concepts of multi-level dynamical modelling: understanding mechanisms of squamous cell carcinoma development in Fanconi anemia

Fanconi anemia (FA) is a rare disease (incidence of 1:300,000) primarily based on the inheritance of pathogenic variants in genes of the FA/BRCA (breast cancer) pathway. These variants ultimately reduce the functionality of different proteins involved in the repair of DNA interstrand crosslinks and DNA double-strand breaks. At birth, individuals with FA might present with typical malformations, particularly radial axis and renal malformations, as well as other physical abnormalities like skin pigmentation anomalies. During the first decade of life, FA mostly causes bone marrow failure due to reduced capacity and loss of the hematopoietic stem and progenitor cells. This often makes hematopoietic stem cell transplantation necessary, but this therapy increases the already intrinsic risk of developing squamous cell carcinoma (SCC) in early adult age. Due to the underlying genetic defect in FA, classical chemo-radiation-based treatment protocols cannot be applied. Therefore, detecting and treating the multi-step tumorigenesis process of SCC in an early stage, or even its progenitors, is the best option for prolonging the life of adult FA individuals. However, the small number of FA individuals makes classical evidence-based medicine approaches based on results from randomized clinical trials impossible. As an alternative, we introduce here the concept of multi-level dynamical modelling using large, longitudinally collected genome, proteome- and transcriptome-wide data sets from a small number of FA individuals. This mechanistic modelling approach is based on the "hallmarks of cancer in FA", which we derive from our unique database of the clinical history of over 750 FA individuals. Multi-omic data from healthy and diseased tissue samples of FA individuals are to be used for training constituent models of a multi-level tumorigenesis model, which will then be used to make experimentally testable predictions. In this way, mechanistic models facilitate not only a descriptive but also a functional understanding of SCC in FA. This approach will provide the basis for detecting signatures of SCCs at early stages and their precursors so they can be efficiently treated or even prevented, leading to a better prognosis and quality of life for the FA individual.

An end-to-end workflow for multiplexed image processing and analysis

Multiplexed imaging enables the simultaneous spatial profiling of dozens of biological molecules in tissues at single-cell resolution. Extracting biologically relevant information, such as the spatial distribution of cell phenotypes from multiplexed tissue imaging data, involves a number of computational tasks, including image segmentation, feature extraction and spatially resolved single-cell analysis. Here, we present an end-to-end workflow for multiplexed tissue image processing and analysis that integrates previously developed computational tools to enable these tasks in a user-friendly and customizable fashion. For data quality assessment, we highlight the utility of napari-imc for interactively inspecting raw imaging data and the cytomapper R/Bioconductor package for image visualization in R. Raw data preprocessing, image segmentation and feature extraction are performed using the steinbock toolkit. We showcase two alternative approaches for segmenting cells on the basis of supervised pixel classification and pretrained deep learning models. The extracted single-cell data are then read, processed and analyzed in R. The protocol describes the use of community-established data containers, facilitating the application of R/Bioconductor packages for dimensionality reduction, single-cell visualization and phenotyping. We provide instructions for performing spatially resolved single-cell analysis, including community analysis, cellular neighborhood detection and cell-cell interaction testing using the imcRtools R/Bioconductor package. The workflow has been previously applied to imaging mass cytometry data, but can be easily adapted to other highly multiplexed imaging technologies. This protocol can be implemented by researchers with basic bioinformatics training, and the analysis of the provided dataset can be completed within 5-6 h. An extended version is available at https://bodenmillergroup.github.io/IMCDataAnalysis/ .

Single-cell and spatial analyses reveal a tradeoff between murine mammary proliferation and lineage programs associated with endocrine cues

Although distinct epithelial cell types have been distinguished in glandular tissues such as the mammary gland, the extent of heterogeneity within each cell type and the degree of endocrine control of this diversity across development are incompletely understood. By combining mass cytometry and cyclic immunofluorescence, we define a rich array of murine mammary epithelial cell subtypes associated with puberty, the estrous cycle, and sex. These subtypes are differentially proliferative and spatially segregate distinctly in adult versus pubescent glands. Further, we identify systematic suppression of lineage programs at the protein and RNA levels as a common feature of mammary epithelial expansion during puberty, the estrous cycle, and gestation and uncover a pervasive enrichment of ribosomal protein genes in luminal cells elicited specifically during progesterone-dominant expansionary periods. Collectively, these data expand our knowledge of murine mammary epithelial heterogeneity and connect endocrine-driven epithelial expansion with lineage suppression.

An update on methods for detection of prognostic and predictive biomarkers in melanoma

The approval of immunotherapy for stage II-IV melanoma has underscored the need for improved immune-based predictive and prognostic biomarkers. For resectable stage II-III patients, adjuvant immunotherapy has proven clinical benefit, yet many patients experience significant adverse events and may not require therapy. In the metastatic setting, single agent immunotherapy cures many patients but, in some cases, more intensive combination therapies against specific molecular targets are required. Therefore, the establishment of additional biomarkers to determine a patient's disease outcome (i.e., prognostic) or response to treatment (i.e., predictive) is of utmost importance. Multiple methods ranging from gene expression profiling of bulk tissue, to spatial transcriptomics of single cells and artificial intelligence-based image analysis have been utilized to better characterize the immune microenvironment in melanoma to provide novel predictive and prognostic biomarkers. In this review, we will highlight the different techniques currently under investigation for the detection of prognostic and predictive immune biomarkers in melanoma.

ARF6-dependent endocytic trafficking of the Interferon-γ receptor drives adaptive immune resistance in cancer

Adaptive immune resistance (AIR) is a protective process used by cancer to escape elimination by CD8+ T cells. Inhibition of immune checkpoints PD-1 and CTLA-4 specifically target Interferon-gamma (IFNγ)-driven AIR. AIR begins at the plasma membrane where tumor cell-intrinsic cytokine signaling is initiated. Thus, plasma membrane remodeling by endomembrane trafficking could regulate AIR. Herein we report that the trafficking protein ADP-Ribosylation Factor 6 (ARF6) is critical for IFNγ-driven AIR. ARF6 prevents transport of the receptor to the lysosome, augmenting IFNγR expression, tumor intrinsic IFNγ signaling and downstream expression of immunosuppressive genes. In murine melanoma, loss of ARF6 causes resistance to immune checkpoint blockade (ICB). Likewise, low expression of ARF6 in patient tumors correlates with inferior outcomes with ICB. Our data provide new mechanistic insights into tumor immune escape, defined by ARF6-dependent AIR, and support that ARF6-dependent endomembrane trafficking of the IFNγ receptor influences outcomes of ICB.

A Standardized and Reproducible Workflow for Membrane Glass Slides in Routine Histology and Spatial Proteomics

Defining the molecular phenotype of single cells in situ is key for understanding tissue architecture in health and disease. Advanced imaging platforms have recently been joined by spatial omics technologies, promising unparalleled insights into the molecular landscape of biological samples. Furthermore, high-precision laser microdissection (LMD) of tissue on membrane glass slides is a powerful method for spatial omics technologies and single-cell type spatial proteomics in particular. However, current histology protocols have not been compatible with glass membrane slides and LMD for automated staining platforms and routine histology procedures. This has prevented the combination of advanced staining procedures with LMD. In this study, we describe a novel method for handling glass membrane slides that enables automated eight-color multiplexed immunofluorescence staining and high-quality imaging followed by precise laser-guided extraction of single cells. The key advance is the glycerol-based modification of heat-induced epitope retrieval protocols, termed "G-HIER." We find that this altered antigen-retrieval solution prevents membrane distortion. Importantly, G-HIER is fully compatible with current antigen retrieval workflows and mass spectrometry-based proteomics and does not affect proteome depth or quality. To demonstrate the versatility of G-HIER for spatial proteomics, we apply the recently introduced deep visual proteomics technology to perform single-cell type analysis of adjacent suprabasal and basal keratinocytes of human skin. G-HIER overcomes previous incompatibility of standard and advanced staining protocols with membrane glass slides and enables robust integration with routine histology procedures, high-throughput multiplexed imaging, and sophisticated downstream spatial omics technologies.

Multiplex imaging in immuno-oncology

Multiplex imaging has emerged as an invaluable tool for immune-oncologists and translational researchers, enabling them to examine intricate interactions among immune cells, stroma, matrix, and malignant cells within the tumor microenvironment (TME). It holds significant promise in the quest to discover improved biomarkers for treatment stratification and identify novel therapeutic targets. Nonetheless, several challenges exist in the realms of study design, experiment optimization, and data analysis. In this review, our aim is to present an overview of the utilization of multiplex imaging in immuno-oncology studies and inform novice researchers about the fundamental principles at each stage of the imaging and analysis process.

SpaRx: elucidate single-cell spatial heterogeneity of drug responses for personalized treatment

Spatial cellular authors heterogeneity contributes to differential drug responses in a tumor lesion and potential therapeutic resistance. Recent emerging spatial technologies such as CosMx, MERSCOPE and Xenium delineate the spatial gene expression patterns at the single cell resolution. This provides unprecedented opportunities to identify spatially localized cellular resistance and to optimize the treatment for individual patients. In this work, we present a graph-based domain adaptation model, SpaRx, to reveal the heterogeneity of spatial cellular response to drugs. SpaRx transfers the knowledge from pharmacogenomics profiles to single-cell spatial transcriptomics data, through hybrid learning with dynamic adversarial adaption. Comprehensive benchmarking demonstrates the superior and robust performance of SpaRx at different dropout rates, noise levels and transcriptomics coverage. Further application of SpaRx to the state-of-the-art single-cell spatial transcriptomics data reveals that tumor cells in different locations of a tumor lesion present heterogenous sensitivity or resistance to drugs. Moreover, resistant tumor cells interact with themselves or the surrounding constituents to form an ecosystem for drug resistance. Collectively, SpaRx characterizes the spatial therapeutic variability, unveils the molecular mechanisms underpinning drug resistance and identifies personalized drug targets and effective drug combinations.

Progress in research on tumor microenvironment-based spatial omics technologies

Spatial omics technology integrates the concept of space into omics research and retains the spatial information of tissues or organs while obtaining molecular information. It is characterized by the ability to visualize changes in molecular information and yields intuitive and vivid visual results. Spatial omics technologies include spatial transcriptomics, spatial proteomics, spatial metabolomics, and other technologies, the most widely used of which are spatial transcriptomics and spatial proteomics. The tumor microenvironment refers to the surrounding microenvironment in which tumor cells exist, including the surrounding blood vessels, immune cells, fibroblasts, bone marrow-derived inflammatory cells, various signaling molecules, and extracellular matrix. A key issue in modern tumor biology is the application of spatial omics to the study of the tumor microenvironment, which can reveal problems that conventional research techniques cannot, potentially leading to the development of novel therapeutic agents for cancer. This paper summarizes the progress of research on spatial transcriptomics and spatial proteomics technologies for characterizing the tumor immune microenvironment.

The spread of interferon-γ in melanomas is highly spatially confined, driving nongenetic variability in tumor cells

Interferon-γ (IFNγ) is a critical antitumor cytokine that has varied effects on different cell types. The global effect of IFNγ in the tumor depends on which cells it acts upon and the spatial extent of its spread. Reported measurements of IFNγ spread vary dramatically in different contexts, ranging from nearest-neighbor signaling to perfusion throughout the entire tumor. Here, we apply theoretical considerations to experiments both in vitro and in vivo to study the spread of IFNγ in melanomas. We observe spatially confined niches of IFNγ signaling in 3-D mouse melanoma cultures and human tumors that generate cellular heterogeneity in gene expression and alter the susceptibility of affected cells to T cell killing. Widespread IFNγ signaling only occurs when niches overlap due to high local densities of IFNγ-producing T cells. We measured length scales of ~30 to 40 μm for IFNγ spread in B16 mouse melanoma cultures and human primary cutaneous melanoma. Our results are consistent with IFNγ spread being governed by a simple diffusion-consumption model and offer insight into how the spatial organization of T cells contributes to intratumor heterogeneity in inflammatory signaling, gene expression, and immune-mediated clearance. Solid tumors are often viewed as collections of diverse cellular "neighborhoods": Our work provides a general explanation for such nongenetic cellular variability due to confinement in the spread of immune mediators.

An introduction to representation learning for single-cell data analysis

Single-cell-resolved systems biology methods, including omics- and imaging-based measurement modalities, generate a wealth of high-dimensional data characterizing the heterogeneity of cell populations. Representation learning methods are routinely used to analyze these complex, high-dimensional data by projecting them into lower-dimensional embeddings. This facilitates the interpretation and interrogation of the structures, dynamics, and regulation of cell heterogeneity. Reflecting their central role in analyzing diverse single-cell data types, a myriad of representation learning methods exist, with new approaches continually emerging. Here, we contrast general features of representation learning methods spanning statistical, manifold learning, and neural network approaches. We consider key steps involved in representation learning with single-cell data, including data pre-processing, hyperparameter optimization, downstream analysis, and biological validation. Interdependencies and contingencies linking these steps are also highlighted. This overview is intended to guide researchers in the selection, application, and optimization of representation learning strategies for current and future single-cell research applications.

SpaRx: Elucidate single-cell spatial heterogeneity of drug responses for personalized treatment

Spatial cellular heterogeneity contributes to differential drug responses in a tumor lesion and potential therapeutic resistance. Recent emerging spatial technologies such as CosMx SMI, MERSCOPE, and Xenium delineate the spatial gene expression patterns at the single cell resolution. This provides unprecedented opportunities to identify spatially localized cellular resistance and to optimize the treatment for individual patients. In this work, we present a graph-based domain adaptation model, SpaRx, to reveal the heterogeneity of spatial cellular response to drugs. SpaRx transfers the knowledge from pharmacogenomics profiles to single-cell spatial transcriptomics data, through hybrid learning with dynamic adversarial adaption. Comprehensive benchmarking demonstrates the superior and robust performance of SpaRx at different dropout rates, noise levels, and transcriptomics coverage. Further application of SpaRx to the state-of-art single-cell spatial transcriptomics data reveals that tumor cells in different locations of a tumor lesion present heterogenous sensitivity or resistance to drugs. Moreover, resistant tumor cells interact with themselves or the surrounding constituents to form an ecosystem for drug resistance. Collectively, SpaRx characterizes the spatial therapeutic variability, unveils the molecular mechanisms underpinning drug resistance, and identifies personalized drug targets and effective drug combinations.

Key Points

We have developed a novel graph-based domain adaption model named SpaRx, to reveal the heterogeneity of spatial cellular response to different types of drugs, which bridges the gap between pharmacogenomics knowledgebase and single-cell spatial transcriptomics data.SpaRx is developed tailored for single-cell spatial transcriptomics data and is provided available as a ready-to-use open-source software, which demonstrates high accuracy and robust performance.SpaRx uncovers that tumor cells located in different areas within tumor lesion exhibit varying levels of sensitivity or resistance to drugs. Moreover, SpaRx reveals that tumor cells interact with themselves and the surrounding microenvironment to form an ecosystem capable of drug resistance.

Advances in Mapping Tumor Progression from Precancer Atlases

Tissue profiling technologies present opportunities for understanding transition from precancerous lesions to malignancy, which may impact risk stratification, prevention, and even cancer treatment. A human precancer atlas building effort is ongoing to tackle the significant challenge of decoding the heterogeneity among cells, specimens, and patients. Here, we discuss the findings resulting from atlases built across precancer types, including those found in colon, breast, lung, stomach, cervix, and skin, using bulk, single-cell, and spatial profiling strategies. We highlight two main themes that emerge across precancer types: the ordering of molecular events that occur during tumor progression and the fluctuation of microenvironmental response during precancer progression. We further highlight the key challenges of data integration across large cohorts of patients, and the need for computational tools to reliably annotate and quality control high-volume, high-dimensional data.

MDM2 Inhibition Enhances Immune Checkpoint Inhibitor Efficacy by Increasing IL15 and MHC Class II Production

The treatment of patients with metastatic melanoma with immune checkpoint inhibitors (ICI) leads to impressive response rates but primary and secondary resistance to ICI reduces progression-free survival. Novel strategies that interfere with resistance mechanisms are key to further improve patient outcome during ICI therapy. P53 is often inactivated by mouse-double-minute-2 (MDM2), which may decrease immunogenicity of melanoma cells. We analyzed primary patient-derived melanoma cell lines, performed bulk sequencing analysis of patient-derived melanoma samples, and used melanoma mouse models to investigate the role of MDM2-inhibition for enhanced ICI therapy. We found increased expression of IL15 and MHC-II in murine melanoma cells upon p53 induction by MDM2-inhibition. MDM2-inhibitor induced MHC-II and IL15-production, which was p53 dependent as Tp53 knockdown blocked the effect. Lack of IL15-receptor in hematopoietic cells or IL15 neutralization reduced the MDM2-inhibition/p53-induction-mediated antitumor immunity. P53 induction by MDM2-inhibition caused anti-melanoma immune memory as T cells isolated from MDM2-inhibitor-treated melanoma-bearing mice exhibited anti-melanoma activity in secondary melanoma-bearing mice. In patient-derived melanoma cells p53 induction by MDM2-inhibition increased IL15 and MHC-II. IL15 and CIITA expressions were associated with a more favorable prognosis in patients bearing WT but not TP53-mutated melanoma.

IMPLICATIONS

MDM2-inhibition represents a novel strategy to enhance IL15 and MHC-II-production, which disrupts the immunosuppressive tumor microenvironment. On the basis of our findings, a clinical trial combining MDM2-inhibition with anti-PD-1 immunotherapy for metastatic melanoma is planned.

A review on deep learning applications in highly multiplexed tissue imaging data analysis

Since its introduction into the field of oncology, deep learning (DL) has impacted clinical discoveries and biomarker predictions. DL-driven discoveries and predictions in oncology are based on a variety of biological data such as genomics, proteomics, and imaging data. DL-based computational frameworks can predict genetic variant effects on gene expression, as well as protein structures based on amino acid sequences. Furthermore, DL algorithms can capture valuable mechanistic biological information from several spatial "omics" technologies, such as spatial transcriptomics and spatial proteomics. Here, we review the impact that the combination of artificial intelligence (AI) with spatial omics technologies has had on oncology, focusing on DL and its applications in biomedical image analysis, encompassing cell segmentation, cell phenotype identification, cancer prognostication, and therapy prediction. We highlight the advantages of using highly multiplexed images (spatial proteomics data) compared to single-stained, conventional histopathological ("simple") images, as the former can provide deep mechanistic insights that cannot be obtained by the latter, even with the aid of explainable AI. Furthermore, we provide the reader with the advantages/disadvantages of DL-based pipelines used in preprocessing highly multiplexed images (cell segmentation, cell type annotation). Therefore, this review also guides the reader to choose the DL-based pipeline that best fits their data. In conclusion, DL continues to be established as an essential tool in discovering novel biological mechanisms when combined with technologies such as highly multiplexed tissue imaging data. In balance with conventional medical data, its role in clinical routine will become more important, supporting diagnosis and prognosis in oncology, enhancing clinical decision-making, and improving the quality of care for patients. Since its introduction into the field of oncology, deep learning (DL) has impacted clinical discoveries and biomarker predictions. DL-driven discoveries and predictions in oncology are based on a variety of biological data such as genomics, proteomics, and imaging data. DL-based computational frameworks can predict genetic variant effects on gene expression, as well as protein structures based on amino acid sequences. Furthermore, DL algorithms can capture valuable mechanistic biological information from several spatial "omics" technologies, such as spatial transcriptomics and spatial proteomics. Here, we review the impact that the combination of artificial intelligence (AI) with spatial omics technologies has had on oncology, focusing on DL and its applications in biomedical image analysis, encompassing cell segmentation, cell phenotype identification, cancer prognostication, and therapy prediction. We highlight the advantages of using highly multiplexed images (spatial proteomics data) compared to single-stained, conventional histopathological ("simple") images, as the former can provide deep mechanistic insights that cannot be obtained by the latter, even with the aid of explainable AI. Furthermore, we provide the reader with the advantages/disadvantages of the DL-based pipelines used in preprocessing the highly multiplexed images (cell segmentation, cell type annotation). Therefore, this review also guides the reader to choose the DL-based pipeline that best fits their data. In conclusion, DL continues to be established as an essential tool in discovering novel biological mechanisms when combined with technologies such as highly multiplexed tissue imaging data. In balance with conventional medical data, its role in clinical routine will become more important, supporting diagnosis and prognosis in oncology, enhancing clinical decision-making, and improving the quality of care for patients.

Expanded vacuum-stable gels for multiplexed high-resolution spatial histopathology

Cellular organization and functions encompass multiple scales in vivo. Emerging high-plex imaging technologies are limited in resolving subcellular biomolecular features. Expansion Microscopy (ExM) and related techniques physically expand samples for enhanced spatial resolution, but are challenging to be combined with high-plex imaging technologies to enable integrative multiscaled tissue biology insights. Here, we introduce Expand and comPRESS hydrOgels (ExPRESSO), an ExM framework that allows high-plex protein staining, physical expansion, and removal of water, while retaining the lateral tissue expansion. We demonstrate ExPRESSO imaging of archival clinical tissue samples on Multiplexed Ion Beam Imaging and Imaging Mass Cytometry platforms, with detection capabilities of > 40 markers. Application of ExPRESSO on archival human lymphoid and brain tissues resolved tissue architecture at the subcellular level, particularly that of the blood-brain barrier. ExPRESSO hence provides a platform for extending the analysis compatibility of hydrogel-expanded biospecimens to mass spectrometry, with minimal modifications to protocols and instrumentation.

High-plex immunofluorescence imaging and traditional histology of the same tissue section for discovering image-based biomarkers

Precision medicine is critically dependent on better methods for diagnosing and staging disease and predicting drug response. Histopathology using hematoxylin and eosin (H&E)-stained tissue (not genomics) remains the primary diagnostic method in cancer. Recently developed highly multiplexed tissue imaging methods promise to enhance research studies and clinical practice with precise, spatially resolved single-cell data. Here, we describe the 'Orion' platform for collecting H&E and high-plex immunofluorescence images from the same cells in a whole-slide format suitable for diagnosis. Using a retrospective cohort of 74 colorectal cancer resections, we show that immunofluorescence and H&E images provide human experts and machine learning algorithms with complementary information that can be used to generate interpretable, multiplexed image-based models predictive of progression-free survival. Combining models of immune infiltration and tumor-intrinsic features achieves a 10- to 20-fold discrimination between rapid and slow (or no) progression, demonstrating the ability of multimodal tissue imaging to generate high-performance biomarkers.

Functionally and Metabolically Divergent Melanoma-Associated Macrophages Originate from Common Bone-Marrow Precursors

Tumor-associated macrophages (TAMs) can be widely heterogeneous, based on their ontogeny and function, and driven by the tissue-specific niche. TAMs are highly abundant in the melanoma tumor microenvironment (TME), usually correlating with worse prognoses. However, the understanding of their diversity may be harnessed for therapeutic purposes. Here, we used the clinically relevant YUMM1.7 model to study melanoma TAM origin and dynamics during tumor progression. In i.d. YUMM1.7 tumors, we identified distinct TAM subsets based on F4/80 expression, with the F4/80^high fraction increasing over time and displaying a tissue-resident-like phenotype. While skin-resident macrophages showed mixed ontogeny, F4/80⁺ TAM subsets in the melanoma TME originated almost exclusively from bone-marrow precursors. A multiparametric analysis of the macrophage phenotype showed a temporal divergence of the F4/80⁺ TAM subpopulations, which also differed from the skin-resident subsets and their monocytic precursors. Overall, the F4/80⁺ TAMs displayed co-expressions of M1- and M2-like canonical markers, while RNA sequencing showed differential immunosuppressive and metabolic profiles. Gene-set enrichment analysis (GSEA) revealed F4/80^high TAMs to rely on oxidative phosphorylation, with increased proliferation and protein secretion, while F4/80^low cells had high pro-inflammatory and intracellular signaling pathways, with lipid and polyamine metabolism. Overall, we provide an in-depth characterization of and compelling evidence for the BM-dependency of melanoma TAMs. Interestingly, the transcriptomic analysis of these BM-derived TAMs matched macrophage subsets with mixed ontogeny, which have been observed in other tumor models. Our findings may serve as a guide for identifying potential ways of targeting specific immunosuppressive TAMs in melanoma.

Tumor cell heterogeneity drives spatial organization of the intratumoral immune response in squamous cell skin carcinoma

Intratumoral heterogeneity (ITH)-defined as genetic and cellular diversity within a tumor-is linked to failure of immunotherapy and an inferior anti-tumor immune response. The underlying mechanism of this association is unknown. To address this question, we modeled heterogeneous tumors comprised of a pro-inflammatory ("hot") and an immunosuppressive ("cold") tumor population, labeled with YFP and RFP tags respectively to enable precise spatial tracking. The resulting mixed-population tumors exhibited distinct regions comprised of YFP+ (hot) cells, RFP+ (cold) cells, or a mixture. We found that tumor regions occupied by hot tumor cells (YFP+) harbored more total T cells and a higher frequency of Th1 cells and IFNγ+ CD8 T cells compared to regions occupied by cold tumor cells (RFP+), whereas immunosuppressive macrophages showed the opposite spatial pattern. We identified the chemokine CX3CL1, produced at higher levels by our cold tumors, as a mediator of intratumoral macrophage accumulation, particularly immunosuppressive CD206Hi macrophages. Furthermore, we examined the response of heterogeneous tumors to a therapeutic combination of PD-1 blockade and CD40 agonist on a region-by-region basis. While the combination successfully increases Th1 abundance in "cold" tumor regions, it fails to bring overall T cell activity to the same level as seen in "hot" regions. The presence of the "cold" cells thus ultimately leads to a failure of the therapy to induce tumor rejection. Collectively, our results demonstrate that the organization of heterogeneous tumor cells has a profound impact on directing the spatial organization and function of tumor-infiltrating immune cells as well as on responses to immunotherapy.

Current perspectives on Vaxinia virus: an immuno-oncolytic vector in cancer therapy

Viruses are being researched as cutting-edge therapeutic agents in cancer due to their selective oncolytic action against malignancies. Immuno-oncolytic viruses are a potential category of anticancer treatments because they have natural features that allow viruses to efficiently infect, replicate, and destroy cancer cells. Oncolytic viruses may be genetically modified; engineers can use them as a platform to develop additional therapy modalities that overcome the limitations of current treatment approaches. In recent years, researchers have made great strides in the understanding relationship between cancer and the immune system. An increasing corpus of research is functioning on the immunomodulatory functions of oncolytic virus (OVs). Several clinical studies are currently underway to determine the efficacy of these immuno-oncolytic viruses. These studies are exploring the design of these platforms to elicit the desired immune response and to supplement the available immunotherapeutic modalities to render immune-resistant malignancies amenable to treatment. This review will discuss current research and clinical developments on Vaxinia immuno-oncolytic virus.

Recent advancements in tumour microenvironment landscaping for target selection and response prediction in immune checkpoint therapies achieved through spatial protein multiplexing analysis

Immune checkpoint therapies have significantly advanced cancer treatment. Nevertheless, the high costs and potential adverse effects associated with these therapies highlight the need for better predictive biomarkers to identify patients who are most likely to benefit from treatment. Unfortunately, the existing biomarkers are insufficient to identify such patients. New high-dimensional spatial technologies have emerged as a valuable tool for discovering novel biomarkers by analysing multiple protein markers at a single-cell resolution in tissue samples. These technologies provide a more comprehensive map of tissue composition, cell functionality, and interactions between different cell types in the tumour microenvironment. In this review, we provide an overview of how spatial protein-based multiplexing technologies have fuelled biomarker discovery and advanced the field of immunotherapy. In particular, we will focus on how these technologies contributed to (i) characterise the tumour microenvironment, (ii) understand the role of tumour heterogeneity, (iii) study the interplay of the immune microenvironment and tumour progression, (iv) discover biomarkers for immune checkpoint therapies (v) suggest novel therapeutic strategies.

Functionally and metabolically divergent melanoma-associated macrophages originate from common bone-marrow precursors

Melanomas display high numbers of tumor-associated macrophages (TAMs), which correlate with worse prognosis. Harnessing macrophages for therapeutic purposes has been particularly challenging due to their heterogeneity, based on their ontogeny and function and driven by the tissue-specific niche. In the present study, we used the YUMM1.7 model to better understand melanoma TAM origin and dynamics during tumor progression, with potential therapeutic implications. We identified distinct TAM subsets based on F4/80 expression, with the F4/80 ^high fraction increasing over time and displaying tissue-resident-like phenotype. While skin-resident macrophages showed mixed on-togeny, F4/80 ⁺ TAM subsets in i.d. YUMM1.7 tumors originated almost exclusively from bone-marrow precursors. Mul-tiparametric analysis of macrophage phenotype showed a temporal divergence of F4/80 ⁺ TAM subpopulations, which also differed from skin-resident subsets, and from their monocytic precursors. Overall, F4/80 ⁺ TAMs displayed co-ex-pression of M1- and M2-like canonical markers, while RNA-seq and pathway analysis showed differential immunosup-pressive and metabolic profiles. GSEA showed F4/80 ^high TAMs to rely on oxidative phosphorylation, with increased proliferation and protein secretion while F4/80 ^low cells had high pro-inflammatory and intracellular signaling pathways, with lipid and polyamine metabolism. Overall, the present in-depth characterization provides further evidence of the ontogeny of the evolving melanoma TAMs, whose gene expression profiles matched recently-identified TAM clusters in other tumor models and human cancers. These findings provide evidence for potentially targeting specific immunosup-pressive TAMs in advanced tumor stages.

Computational immune synapse analysis reveals T-cell interactions in distinct tumor microenvironments

The tumor microenvironment (TME) and the cellular interactions within it can be critical to tumor progression and treatment response. Although technologies to generate multiplex images of the TME are advancing, the many ways in which TME imaging data can be mined to elucidate cellular interactions are only beginning to be realized. Here, we present a novel approach for multipronged computational immune synapse analysis (CISA) that reveals T-cell synaptic interactions from multiplex images. CISA enables automated discovery and quantification of immune synapse interactions based on the localization of proteins on cell membranes. We first demonstrate the ability of CISA to detect T-cell:APC (antigen presenting cell) synaptic interactions in two independent human melanoma imaging mass cytometry (IMC) tissue microarray datasets. We then generate melanoma histocytometry whole slide images and verify that CISA can detect similar interactions across data modalities. Interestingly, CISA histoctyometry analysis also reveals that T-cell:macrophage synapse formation is associated with T-cell proliferation. We next show the generality of CISA by extending it to breast cancer IMC images, finding that CISA quantifications of T-cell:B-cell synapses are predictive of improved patient survival. Our work demonstrates the biological and clinical significance of spatially resolving cell-cell synaptic interactions in the TME and provides a robust method to do so across imaging modalities and cancer types.

Spatial distribution of tumor-infiltrating T cells indicated immune response status under chemoradiotherapy plus PD-1 blockade in esophageal cancer

Background

The spatial distribution of tumor-infiltrating T cells and its dynamics during chemoradiotherapy combined with PD-1 blockade is little known in esophageal squamous cell carcinoma (ESCC).

Methods

We applied the multiplex immunofluorescence method to identify T cells (CD4⁺, CD8⁺ T cells, and their PD-1^- or PD-1⁺ subsets) and myeloid-derived cells (CD11c⁺ dendritic cells, CD68⁺ macrophages, and their PD-L1⁺ subpopulations) in paired tumor biopsies (n = 36) collected at baseline and during combination (40 Gy of radiation) from a phase Ib trial (NCT03671265) of ESCC patients treated with first-line chemoradiotherapy plus anti-PD-1 antibody camrelizumab. We used the FoundationOne CDx assay to evaluate tumor mutational burden (TMB) in baseline tumor biopsies (n = 14). We dynamically assessed the nearest distance and proximity of T-cell subsets to tumor cells under combination and estimated the association between T-cell spatial distribution and combination outcome, myeloid-derived subsets, TMB, and patient baseline characteristics.

Findings

We found that the tumor compartment had lower T-cell subsets than the stromal compartment but maintained a comparable level under combination. Both before and under combination, PD-1^- T cells were located closer than PD-1⁺ T cells to tumor cells; T cells, dendritic cells, and macrophages showed the highest accumulation in the 5-10-μm distance. Higher CD4⁺ T cells in the tumor compartment and a shorter nearest distance of T-cell subsets at baseline predicted poor OS. Higher baseline CD4⁺ T cells, dendritic cells, and macrophages were associated with worse OS in less than 10-μm distance to tumor cells, but related with better OS in the farther distance. Higher on-treatment PD-1-positive-expressed CD4⁺ and CD8⁺ T cells within the 100-μm distance to tumor cells predicted longer OS. T cells, dendritic cells, and macrophages showed a positive spatial correlation. Both high TMB and smoking history were associated with a closer location of T cells to tumor cells at baseline.

Conclusions

We firstly illustrated the T-cell spatial distribution in ESCC. Combining chemoradiotherapy with PD-1 blockade could improve the antitumor immune microenvironment, which benefits the treatment outcome. Further understanding the precision spatiality of tumor-infiltrating T cells would provide new evidence for the tumor immune microenvironment and for the combination treatment with immunotherapy.

Robust biomarker discovery through multiplatform multiplex image analysis of breast cancer clinical cohorts

Spatial profiling of tissues promises to elucidate tumor-microenvironment interactions and enable development of spatial biomarkers to predict patient response to immunotherapy and other therapeutics. However, spatial biomarker discovery is often carried out on a single patient cohort or imaging technology, limiting statistical power and increasing the likelihood of technical artifacts. In order to analyze multiple patient cohorts profiled on different platforms, we developed methods for comparative data analysis from three disparate multiplex imaging technologies: 1) cyclic immunofluorescence data we generated from 102 breast cancer patients with clinical follow-up, in addition to publicly available 2) imaging mass cytometry and 3) multiplex ion-beam imaging data. We demonstrate similar single-cell phenotyping results across breast cancer patient cohorts imaged with these three technologies and identify cellular abundance and proximity-based biomarkers with prognostic value across platforms. In multiple platforms, we identified lymphocyte infiltration as independently associated with longer survival in triple negative and high-proliferation breast tumors. Then, a comparison of nine spatial analysis methods revealed robust spatial biomarkers. In estrogen receptor-positive disease, quiescent stromal cells close to tumor were more abundant in good prognosis tumors while tumor neighborhoods of mixed fibroblast phenotypes were enriched in poor prognosis tumors. In triple-negative breast cancer (TNBC), macrophage proximity to tumor and B cell proximity to T cells were greater in good prognosis tumors, while tumor neighborhoods of vimentin-positive fibroblasts were enriched in poor prognosis tumors. We also tested previously published spatial biomarkers in our ensemble cohort, reproducing the positive prognostic value of isolated lymphocytes and lymphocyte occupancy and failing to reproduce the prognostic value of tumor-immune mixing score in TNBC. In conclusion, we demonstrate assembly of larger clinical cohorts from diverse platforms to aid in prognostic spatial biomarker identification and validation.

Spatial transcriptomics in human skin research

Spatial transcriptomics is a revolutionary technique that enables researchers to characterise tissue architecture and localisation of gene expression. A plethora of technologies that map gene expression are currently being developed, aiming to facilitate spatially resolved, high-dimensional assessment of gene transcription in the context of human skin research. Knowing which gene is expressed by which cell and in which location within skin, facilitates understanding of skin function and dysfunction in both health and disease. In this review, we summarise the available spatial transcriptomic methods and we describe their application to a broad spectrum of dermatological diseases.

Localized immune surveillance of primary melanoma in the skin deciphered through executable modeling

While skin is a site of active immune surveillance, primary melanomas often escape detection. Here, we have developed an in silico model to determine the local cross-talk between melanomas and Langerhans cells (LCs), the primary antigen-presenting cells at the site of melanoma development. The model predicts that melanomas fail to activate LC migration to lymph nodes until tumors reach a critical size, which is determined by a positive TNF-α feedback loop within melanomas, in line with our observations of murine tumors. In silico drug screening, supported by subsequent experimental testing, shows that treatment of primary tumors with MAPK pathway inhibitors may further prevent LC migration. In addition, our in silico model predicts treatment combinations that bypass LC dysfunction. In conclusion, our combined approach of in silico and in vivo studies suggests a molecular mechanism that explains how early melanomas develop under the radar of immune surveillance by LC.