Craters on the melanoma surface facilitate tumor-immune interactions and demonstrate pathologic response to checkpoint blockade in humans.

Immunotherapy leads to cancer eradication despite the tumor's immunosuppressive environment. Here, we used extended long-term in-vivo imaging and high-resolution spatial transcriptomics of endogenous melanoma in zebrafish, and multiplex imaging of human melanoma, to identify domains that facilitate immune response during immunotherapy. We identified crater-shaped pockets at the margins of zebrafish and human melanoma, rich with beta-2 microglobulin (B2M) and antigen recognition molecules. The craters harbor the highest density of CD8+ T cells in the tumor. In zebrafish, CD8+ T cells formed prolonged interactions with melanoma cells within craters, characteristic of antigen recognition. Following immunostimulatory treatment, the craters enlarged and became the major site of activated CD8+ T cell accumulation and tumor killing that was B2M dependent. In humans, craters predicted immune response to ICB therapy, showing response better than high T cell infiltration. This marks craters as potential new diagnostic tool for immunotherapy success and targets to enhance ICB response.

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.

Duvelisib plus romidepsin in relapsed/refractory T cell lymphomas: a phase 1b/2a trial.

PI3K-δ inhibitors have shown impressive activity in lymphoid malignancies but have been hampered by autoimmune and infectious toxicities, leading to market withdrawals. We previously demonstrated activity of the PI3K-δγ inhibitor duvelisib in T cell lymphomas (TCLs) that was associated with inflammatory adverse events. As reported here, we conducted a phase 1b/2a study of duvelisib in combination with either romidepsin (n = 66) or bortezomib (n = 32) in patients with relapsed/refractory TCL and found that the addition of romidepsin, but not bortezomib, appeared to increase efficacy while attenuating PI3K inhibitor-driven toxicity. The primary endpoint of the study was to determine the safety and maximum tolerated dose of duvelisib, which was 75 mg twice daily when combined with romidepsin versus 25 mg twice daily when combined with bortezomib. The most common adverse events were neutropenia (42%, 25/59) and fatigue (37%, 22/59) in patients treated with duvelisib and romidepsin and diarrhea (48%, 11/23) and neutropenia (30%, 7/23) in patients treated with duvelisib and bortezomib. Duvelisib and romidepsin resulted in less grade 3/4 hepatotoxicity (14%, 8/59) compared to 40% (14/35) in our previous study with duvelisib monotherapy. This was associated with reductions in circulating inflammatory mediators and myeloid cell inflammatory gene expression. Secondary endpoints of overall and complete response rates were 55% (35/64) and 34% (22/64) for patients treated with duvelisib and romidepsin and 34% (11/32) and 13% (4/32) for patients treated with duvelisib and bortezomib. Among patients with peripheral T cell lymphomas (PTCLs), overall and complete response rates of duvelisib and romidepsin were 56% (27/48) and 44% (21/48), respectively, with exploratory analyses showing increased response rates in patients with a follicular helper T cell subtype. These findings support further development of combined PI3K and histone deacetylase (HDAC) inhibition in TCLs and suggest a unique strategy to enable PI3K inhibitor-based combinations for additional patient populations. ClinicalTrials.gov identifier: NCT02783625 .

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.

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. A critical aspect of such tissue profiling is quantifying the spatial relationship 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 10^7 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, 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.

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.

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 trafficking 3 . 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) imaging 4 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 situ 5 , 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" associations 6 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.

Development and validation of time-to-event models to predict metastatic recurrence of localized cutaneous melanoma.

BACKGROUND: The recent expansion of immunotherapy for stage IIB/IIC melanoma highlights a growing clinical need to identify patients at high risk of metastatic recurrence and, therefore, most likely to benefit from this therapeutic modality. OBJECTIVE: To develop time-to-event risk prediction models for melanoma metastatic recurrence. METHODS: Patients diagnosed with stage I/II primary cutaneous melanoma between 2000 and 2020 at Mass General Brigham and Dana-Farber Cancer Institute were included. Melanoma recurrence date and type were determined by chart review. Thirty clinicopathologic factors were extracted from electronic health records. Three types of time-to-event machine-learning models were evaluated internally and externally in the distant versus locoregional/nonrecurrence prediction. RESULTS: This study included 954 melanomas (155 distant, 163 locoregional, and 636 1:2 matched nonrecurrences). Distant recurrences were associated with worse survival compared to locoregional/nonrecurrences (HR: 6.21, P < .001) and to locoregional recurrences only (HR: 5.79, P < .001). The Gradient Boosting Survival model achieved the best performance (concordance index: 0.816; time-dependent AUC: 0.842; Brier score: 0.103) in the external validation. LIMITATIONS: Retrospective nature and cohort from one geography. CONCLUSIONS: These results suggest that time-to-event machine-learning models can reliably predict the metastatic recurrence from localized melanoma and help identify high-risk patients who are most likely to benefit from immunotherapy.

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.

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.

TP63 fusions drive multicomplex enhancer rewiring, lymphomagenesis, and EZH2 dependence.

Gene fusions involving tumor protein p63 gene (TP63) occur in multiple T and B cell lymphomas and portend a dismal prognosis for patients. The function and mechanisms of TP63 fusions remain unclear, and there is no target therapy for patients with lymphoma harboring TP63 fusions. Here, we show that TP63 fusions act as bona fide oncogenes and are essential for fusion-positive lymphomas. Transgenic mice expressing TBL1XR1::TP63, the most common TP63 fusion, develop diverse lymphomas that recapitulate multiple human T and B cell lymphomas. Here, we identify that TP63 fusions coordinate the recruitment of two epigenetic modifying complexes, the nuclear receptor corepressor (NCoR)-histone deacetylase 3 (HDAC3) by the N-terminal TP63 fusion partner and the lysine methyltransferase 2D (KMT2D) by the C-terminal TP63 component, which are both required for fusion-dependent survival. TBL1XR1::TP63 localization at enhancers drives a unique cell state that involves up-regulation of MYC and the polycomb repressor complex 2 (PRC2) components EED and EZH2. Inhibiting EZH2 with the therapeutic agent valemetostat is highly effective at treating transgenic lymphoma murine models, xenografts, and patient-derived xenografts harboring TP63 fusions. One patient with TP63-rearranged lymphoma showed a rapid response to valemetostat treatment. In summary, TP63 fusions link partner components that, together, coordinate multiple epigenetic complexes, resulting in therapeutic vulnerability to EZH2 inhibition.

Immune Profiling of Dermatologic Adverse Events from Checkpoint Blockade using Tissue Cyclic Immunofluorescence.

In this study, we demonstrate the utility of whole-slide CyCIF (tissue-based cyclic immunofluorescence) imaging for characterizing immune cell infiltrates in immune checkpoint inhibitor (ICI)-induced dermatologic adverse events (dAEs). We analyzed six cases of ICI-induced dAEs, including lichenoid, bullous pemphigoid, psoriasis, and eczematous eruptions, comparing immune profiling results obtained using both standard immunohistochemistry (IHC) and CyCIF. Our findings indicate that CyCIF provides more detailed and precise single-cell characterization of immune cell infiltrates than IHC, which relies on semi-quantitative scoring by pathologists. This pilot study highlights the potential of CyCIF to advance our understanding of the immune environment in dAEs by revealing tissue-level spatial patterns of immune cell infiltrates, allowing for more precise phenotypic distinctions and deeper exploration of disease mechanisms. By demonstrating that CyCIF can be performed on friable tissues, such as bullous pemphigoid, we provide a foundation for future studies to examine the drivers of specific dAEs using larger cohorts of phenotyped toxicity and suggest a broader role for highly multiplexed tissue imaging in phenotyping the immune mediated disease that they resemble.

Lymphocyte networks are dynamic cellular communities in the immunoregulatory landscape of lung adenocarcinoma.

Lymphocytes are key for immune surveillance of tumors, but our understanding of the spatial organization and physical interactions that facilitate lymphocyte anti-cancer functions is limited. We used multiplexed imaging, quantitative spatial analysis, and machine learning to create high-definition maps of lung tumors from a Kras/Trp53-mutant mouse model and human resections. Networks of interacting lymphocytes ("lymphonets") emerged as a distinctive feature of the anti-cancer immune response. Lymphonets nucleated from small T cell clusters and incorporated B cells with increasing size. CXCR3-mediated trafficking modulated lymphonet size and number, but T cell antigen expression directed intratumoral localization. Lymphonets preferentially harbored TCF1+ PD-1+ progenitor CD8+ T cells involved in responses to immune checkpoint blockade (ICB) therapy. Upon treatment of mice with ICB or an antigen-targeted vaccine, lymphonets retained progenitor and gained cytotoxic CD8+ T cell populations, likely via progenitor differentiation. These data show that lymphonets create a spatial environment supportive of CD8+ T cell anti-tumor responses.

Single Cell Spatial Biology for Precision Cancer Medicine.

In cancer, complex ecosystems of interacting cell types play fundamental roles in tumor development, progression, and response to therapy. However, the cellular organization, community structure, and spatially defined microenvironments of human tumors remain poorly understood. With the emergence of new technologies for high-throughput spatial profiling of complex tissue specimens, it is now possible to identify clinically significant spatial features with high granularity. In this PSB workshop, we will highlight recent advances in this area and explore how single cell spatial profiling can advance precision cancer medicine.

Multiplexed and reproducible high content screening of live and fixed cells using Dye Drop.

High-throughput measurement of cells perturbed using libraries of small molecules, gene knockouts, or different microenvironmental factors is a key step in functional genomics and pre-clinical drug discovery. However, it remains difficult to perform accurate single-cell assays in 384-well plates, limiting many studies to well-average measurements (e.g., CellTiter-Glo®). Here we describe a public domain Dye Drop method that uses sequential density displacement and microscopy to perform multi-step assays on living cells. We use Dye Drop cell viability and DNA replication assays followed by immunofluorescence imaging to collect single-cell dose-response data for 67 investigational and clinical-grade small molecules in 58 breast cancer cell lines. By separating the cytostatic and cytotoxic effects of drugs computationally, we uncover unexpected relationships between the two. Dye Drop is rapid, reproducible, customizable, and compatible with manual or automated laboratory equipment. Dye Drop improves the tradeoff between data content and cost, enabling the collection of information-rich perturbagen-response datasets.

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

Cutaneous melanoma is a highly immunogenic malignancy that is surgically curable at early stages but life-threatening when metastatic. Here we integrate high-plex imaging, 3D high-resolution microscopy, and spatially resolved microregion transcriptomics to study immune evasion and immunoediting in primary melanoma. We find that recurrent cellular neighborhoods involving tumor, immune, and stromal cells change significantly along a progression axis involving precursor states, melanoma in situ, and invasive tumor. Hallmarks of immunosuppression are already detectable in precursor regions. When tumors become locally invasive, a consolidated and spatially restricted suppressive environment forms along the tumor-stromal boundary. This environment is established by cytokine gradients that promote expression of MHC-II and IDO1, and by PD1-PDL1-mediated cell contacts involving macrophages, dendritic cells, and T cells. A few millimeters away, cytotoxic T cells synapse with melanoma cells in fields of tumor regression. Thus, invasion and immunoediting can coexist within a few millimeters of each other in a single specimen. SIGNIFICANCE: The reorganization of the tumor ecosystem in primary melanoma is an excellent setting in which to study immunoediting and immune evasion. Guided by classic histopathology, spatial profiling of proteins and mRNA reveals recurrent morphologic and molecular features of tumor evolution that involve localized paracrine cytokine signaling and direct cell-cell contact. This article is highlighted in the In This Issue feature, p. 1397.

MCMICRO: a scalable, modular image-processing pipeline for multiplexed tissue imaging.

Highly multiplexed tissue imaging makes detailed molecular analysis of single cells possible in a preserved spatial context. However, reproducible analysis of large multichannel images poses a substantial computational challenge. Here, we describe a modular and open-source computational pipeline, MCMICRO, for performing the sequential steps needed to transform whole-slide images into single-cell data. We demonstrate the use of MCMICRO on tissue and tumor images acquired using multiple imaging platforms, thereby providing a solid foundation for the continued development of tissue imaging software.

Overcoming IMiD resistance in T-cell lymphomas through potent degradation of ZFP91 and IKZF1.

Immunomodulatory (IMiD) agents like lenalidomide and pomalidomide induce the recruitment of IKZF1 and other targets to the CRL4CRBN E3 ubiquitin ligase, resulting in their ubiquitination and degradation. These agents are highly active in B-cell lymphomas and a subset of myeloid diseases but have compromised effects in T-cell lymphomas (TCLs). Here, we show that 2 factors determine resistance to IMiDs among TCLs. First, limited CRBN expression reduces IMiD activity in TCLs but can be overcome by newer-generation degrader CC-92480. Using mass spectrometry, we show that CC-92480 selectively degrades IKZF1 and ZFP91 in TCL cells with greater potency than pomalidomide. As a result, CC-92480 is highly active against multiple TCL subtypes and showed greater efficacy than pomalidomide across 4 in vivo TCL models. Second, we demonstrate that ZFP91 functions as a bona fide transcription factor that coregulates cell survival with IKZF1 in IMiD-resistant TCLs. By activating keynote genes from WNT, NF-kB, and MAP kinase signaling, ZFP91 directly promotes resistance to IKZF1 loss. Moreover, lenalidomide-sensitive TCLs can acquire stable resistance via ZFP91 rewiring, which involves casein kinase 2-mediated c-Jun inactivation. Overall, these findings identify a critical transcription factor network within TCLs and provide clinical proof of concept for the novel therapy using next-generation degraders.