Starfysh integrates spatial transcriptomic and histologic data to reveal heterogeneous tumor-immune hubs.

Spatially resolved gene expression profiling provides insight into tissue organization and cell-cell crosstalk; however, sequencing-based spatial transcriptomics (ST) lacks single-cell resolution. Current ST analysis methods require single-cell RNA sequencing data as a reference for rigorous interpretation of cell states, mostly do not use associated histology images and are not capable of inferring shared neighborhoods across multiple tissues. Here we present Starfysh, a computational toolbox using a deep generative model that incorporates archetypal analysis and any known cell type markers to characterize known or new tissue-specific cell states without a single-cell reference. Starfysh improves the characterization of spatial dynamics in complex tissues using histology images and enables the comparison of niches as spatial hubs across tissues. Integrative analysis of primary estrogen receptor (ER)-positive breast cancer, triple-negative breast cancer (TNBC) and metaplastic breast cancer (MBC) tissues led to the identification of spatial hubs with patient- and disease-specific cell type compositions and revealed metabolic reprogramming shaping immunosuppressive hubs in aggressive MBC.

Coordinated Immune Cell Networks in the Bone Marrow Microenvironment Define the Graft versus Leukemia Response with Adoptive Cellular Therapy.

Understanding how intra-tumoral immune populations coordinate to generate anti-tumor responses following therapy can guide precise treatment prioritization. We performed systematic dissection of an established adoptive cellular therapy, donor lymphocyte infusion (DLI), by analyzing 348,905 single-cell transcriptomes from 74 longitudinal bone-marrow samples of 25 patients with relapsed myeloid leukemia; a subset was evaluated by protein-based spatial analysis. In acute myelogenous leukemia (AML) responders, diverse immune cell types within the bone-marrow microenvironment (BME) were predicted to interact with a clonally expanded population of ZNF683 + GZMB + CD8+ cytotoxic T lymphocytes (CTLs) which demonstrated in vitro specificity for autologous leukemia. This population, originating predominantly from the DLI product, expanded concurrently with NK and B cells. AML nonresponder BME revealed a paucity of crosstalk and elevated TIGIT expression in CD8+ CTLs. Our study highlights recipient BME differences as a key determinant of effective anti-leukemia response and opens new opportunities to modulate cell-based leukemia-directed therapy.

Mapping the evolution of T cell states during response and resistance to adoptive cellular therapy.

To elucidate mechanisms by which T cells eliminate leukemia, we study donor lymphocyte infusion (DLI), an established immunotherapy for relapsed leukemia. We model T cell dynamics by integrating longitudinal, multimodal data from 94,517 bone marrow-derived single T cell transcriptomes in addition to chromatin accessibility and single T cell receptor sequencing from patients undergoing DLI. We find that responsive tumors are defined by enrichment of late-differentiated T cells before DLI and rapid, durable expansion of early differentiated T cells after treatment, highly similar to "terminal" and "precursor" exhausted subsets, respectively. Resistance, in contrast, is defined by heterogeneous T cell dysfunction. Surprisingly, early differentiated T cells in responders mainly originate from pre-existing and novel clonotypes recruited to the leukemic microenvironment, rather than the infusion. Our work provides a paradigm for analyzing longitudinal single-cell profiling of scenarios beyond adoptive cell therapy and introduces Symphony, a Bayesian approach to infer regulatory circuitry underlying T cell subsets, with broad relevance to exhaustion antagonists across cancers.