Mapping single-cell developmental potential in health and disease with interpretable deep learning.

Single-cell RNA sequencing (scRNA-seq) has transformed our understanding of cell fate in developmental systems. However, identifying the molecular hallmarks of potency - the capacity of a cell to differentiate into other cell types - has remained challenging. Here, we introduce CytoTRACE 2, an interpretable deep learning framework for characterizing potency and differentiation states on an absolute scale from scRNA-seq data. Across 31 human and mouse scRNA-seq datasets encompassing 28 tissue types, CytoTRACE 2 outperformed existing methods for recovering experimentally determined potency levels and differentiation states covering the entire range of cellular ontogeny. Moreover, it reconstructed the temporal hierarchy of mouse embryogenesis across 62 timepoints; identified pan-tissue expression programs that discriminate major potency levels; and facilitated discovery of cellular phenotypes in cancer linked to survival and immunotherapy resistance. Our results illuminate a fundamental feature of cell biology and provide a broadly applicable platform for delineating single-cell differentiation landscapes in health and disease.

Spatially Segregated Macrophage Populations Predict Distinct Outcomes In Colon Cancer.

Tumor-associated macrophages are transcriptionally heterogeneous, but the spatial distribution and cell interactions that shape macrophage tissue roles remain poorly characterized. Here, we spatially resolve five distinct human macrophage populations in normal and malignant human breast and colon tissue and reveal their cellular associations. This spatial map reveals that distinct macrophage populations reside in spatially segregated micro-environmental niches with conserved cellular compositions that are repeated across healthy and diseased tissue. We show that IL4I1+ macrophages phagocytose dying cells in areas with high cell turnover and predict good outcome in colon cancer. In contrast, SPP1+ macrophages are enriched in hypoxic and necrotic tumor regions and portend worse outcome in colon cancer. A subset of FOLR2+ macrophages is embedded in plasma cell niches. NLRP3+ macrophages co-localize with neutrophils and activate an inflammasome in tumors. Our findings indicate that a limited number of unique human macrophage niches function as fundamental building blocks in tissue.

Sarcoma microenvironment cell states and ecosystems are associated with prognosis and predict response to immunotherapy.

Characterization of the diverse malignant and stromal cell states that make up soft tissue sarcomas and their correlation with patient outcomes has proven difficult using fixed clinical specimens. Here, we employed EcoTyper, a machine-learning framework, to identify the fundamental cell states and cellular ecosystems that make up sarcomas on a large scale using bulk transcriptomes with clinical annotations. We identified and validated 23 sarcoma-specific, transcriptionally defined cell states, many of which were highly prognostic of patient outcomes across independent datasets. We discovered three conserved cellular communities or ecotypes associated with underlying genomic alterations and distinct clinical outcomes. We show that one ecotype defined by tumor-associated macrophages and epithelial-like malignant cells predicts response to immune-checkpoint inhibition but not chemotherapy and validate our findings in an independent cohort. Our results may enable identification of patients with soft tissue sarcomas who could benefit from immunotherapy and help develop new therapeutic strategies.

Single-cell mapping identifies MSI+ cells as a common origin for diverse subtypes of pancreatic cancer.

Identifying the cells from which cancers arise is critical for understanding the molecular underpinnings of tumor evolution. To determine whether stem/progenitor cells can serve as cells of origin, we created a Msi2-CreERT2 knock-in mouse. When crossed to CAG-LSL-MycT58A mice, Msi2-CreERT2 mice developed multiple pancreatic cancer subtypes: ductal, acinar, adenosquamous, and rare anaplastic tumors. Combining single-cell genomics with computational analysis of developmental states and lineage trajectories, we demonstrate that MYC preferentially triggers transformation of the most immature MSI2+ pancreas cells into multi-lineage pre-cancer cells. These pre-cancer cells subsequently diverge to establish pancreatic cancer subtypes by activating distinct transcriptional programs and large-scale genomic changes, and enforced expression of specific signals like Ras can redirect subtype specification. This study shows that multiple pancreatic cancer subtypes can arise from a common pool of MSI2+ cells and provides a powerful model to understand and control the programs that shape divergent fates in pancreatic cancer.

Diversity of intratumoral regulatory T cells in B-cell non-Hodgkin lymphoma.

Tumor-infiltrating regulatory T cells (Tregs) contribute to an immunosuppressive tumor microenvironment. Despite extensive studies, the prognostic impact of tumor-infiltrating Tregs in B-cell non-Hodgkin lymphomas (B-NHLs) remains unclear. Emerging studies suggest substantial heterogeneity in the phenotypes and suppressive capacities of Tregs, emphasizing the importance of understanding Treg diversity and the need for additional markers to identify highly suppressive Tregs. Here, we applied single-cell RNA sequencing and T-cell receptor sequencing combined with high-dimensional cytometry to decipher the heterogeneity of intratumoral Tregs in diffuse large B-cell lymphoma and follicular lymphoma (FL), compared with that in nonmalignant tonsillar tissue. We identified 3 distinct transcriptional states of Tregs: resting, activated, and unconventional LAG3+FOXP3- Tregs. Activated Tregs were enriched in B-NHL tumors, coexpressed several checkpoint receptors, and had stronger immunosuppressive activity compared with resting Tregs. In FL, activated Tregs were found in closer proximity to CD4+ and CD8+ T cells than other cell types. Furthermore, we used a computational approach to develop unique gene signature matrices, which were used to enumerate each Treg subset in cohorts with bulk gene expression data. In 2 independent FL cohorts, activated Tregs was the major subset, and high abundance was associated with adverse outcome. This study demonstrates that Tregs infiltrating B-NHL tumors are transcriptionally and functionally diverse. Highly immunosuppressive activated Tregs were enriched in tumor tissue but absent in the peripheral blood. Our data suggest that a deeper understanding of Treg heterogeneity in B-NHL could open new paths for rational drug design, facilitating selective targeting to improve antitumor immunity.

CD4 T cells and toxicity from immune checkpoint blockade.

Immune-related toxicities, otherwise known as immune-related adverse events (irAEs), occur in a substantial fraction of cancer patients treated with immune checkpoint inhibitors (ICIs). Ranging from asymptomatic to life-threatening, ICI-induced irAEs can result in hospital admission, high-dose corticosteroid treatment, ICI discontinuation, and in some cases, death. A deeper understanding of the factors underpinning severe irAE development will be essential for improved irAE prediction and prevention, toward maximizing the benefits and safety profiles of ICIs. In recent work, we applied mass cytometry, single-cell RNA sequencing, single-cell V(D)J sequencing, bulk RNA sequencing, and bulk T-cell receptor (TCR) sequencing to identify pretreatment determinants of severe irAE development in patients with advanced melanoma. Across 71 patients separated into three cohorts, we found that two baseline features in circulation-elevated activated CD4 effector memory T-cell abundance and TCR diversity-are associated with severe irAE development, independent of the affected organ system within 3 months of ICI treatment initiation. Here, we provide an extended perspective on this work, synthesize and discuss related literature, and summarize practical considerations for clinical translation. Collectively, these findings lay a foundation for data-driven and mechanistic insights into irAE development, with the potential to reduce ICI morbidity and mortality in the future.

p53 governs an AT1 differentiation programme in lung cancer suppression.

Lung cancer is the leading cause of cancer deaths worldwide1. Mutations in the tumour suppressor gene TP53 occur in 50% of lung adenocarcinomas (LUADs) and are linked to poor prognosis1-4, but how p53 suppresses LUAD development remains enigmatic. We show here that p53 suppresses LUAD by governing cell state, specifically by promoting alveolar type 1 (AT1) differentiation. Using mice that express oncogenic Kras and null, wild-type or hypermorphic Trp53 alleles in alveolar type 2 (AT2) cells, we observed graded effects of p53 on LUAD initiation and progression. RNA sequencing and ATAC sequencing of LUAD cells uncovered a p53-induced AT1 differentiation programme during tumour suppression in vivo through direct DNA binding, chromatin remodelling and induction of genes characteristic of AT1 cells. Single-cell transcriptomics analyses revealed that during LUAD evolution, p53 promotes AT1 differentiation through action in a transitional cell state analogous to a transient intermediary seen during AT2-to-AT1 cell differentiation in alveolar injury repair. Notably, p53 inactivation results in the inappropriate persistence of these transitional cancer cells accompanied by upregulated growth signalling and divergence from lung lineage identity, characteristics associated with LUAD progression. Analysis of Trp53 wild-type and Trp53-null mice showed that p53 also directs alveolar regeneration after injury by regulating AT2 cell self-renewal and promoting transitional cell differentiation into AT1 cells. Collectively, these findings illuminate mechanisms of p53-mediated LUAD suppression, in which p53 governs alveolar differentiation, and suggest that tumour suppression reflects a fundamental role of p53 in orchestrating tissue repair after injury.

High-resolution alignment of single-cell and spatial transcriptomes with CytoSPACE.

Recent studies have emphasized the importance of single-cell spatial biology, yet available assays for spatial transcriptomics have limited gene recovery or low spatial resolution. Here we introduce CytoSPACE, an optimization method for mapping individual cells from a single-cell RNA sequencing atlas to spatial expression profiles. Across diverse platforms and tissue types, we show that CytoSPACE outperforms previous methods with respect to noise tolerance and accuracy, enabling tissue cartography at single-cell resolution.

Profiling Cellular Ecosystems at Single-Cell Resolution and at Scale with EcoTyper.

Tissues are composed of diverse cell types and cellular states that organize into distinct ecosystems with specialized functions. EcoTyper is a collection of machine learning tools for the large-scale delineation of cellular ecosystems and their constituent cell states from bulk, single-cell, and spatially resolved gene expression data. In this chapter, we provide a primer on EcoTyper and demonstrate its use for the discovery and recovery of cell states and ecosystems from healthy and diseased tissue specimens.

A spatial map of human macrophage niches reveals context-dependent macrophage functions in colon and breast cancer.

Tumor-associated macrophages (TAMs) display heterogeneous phenotypes. Yet the exact tissue cues that shape macrophage functional diversity are incompletely understood. Here we discriminate, spatially resolve and reveal the function of five distinct macrophage niches within malignant and benign breast and colon tissue. We found that SPP1 TAMs reside in hypoxic and necrotic tumor regions, and a novel subset of FOLR2 tissue resident macrophages (TRMs) supports the plasma cell tissue niche. We discover that IL4I1 macrophages populate niches with high cell turnover where they phagocytose dying cells. Significantly, IL4I1 TAMs abundance correlates with anti-PD1 treatment response in breast cancer. Furthermore, NLRP3 inflammasome activation in NLRP3 TAMs correlates with neutrophil infiltration in the tumors and is associated with poor outcome in breast cancer patients. This suggests the NLRP3 inflammasome as a novel cancer immunetherapy target. Our work uncovers context-dependent roles of macrophage subsets, and suggests novel predictive markers and macrophage subset-specific therapy targets.

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.

Multiomic analysis for optimization of combined focal and immunotherapy protocols in murine pancreatic cancer.

Background: Although combination immunotherapies incorporating local and systemic components have shown promising results in treating solid tumors, varied tumor microenvironments (TMEs) can impact immunotherapeutic efficacy. Method: We designed and evaluated treatment strategies for breast and pancreatic cancer combining magnetic resonance-guided focused ultrasound (MRgFUS) ablation and antibody therapies. With a combination of single-cell sequencing, spectral flow cytometry, and histological analyses, we profiled an immune-suppressed KPC (Kras+/LSL-G12D; Trp53+/LSL-R172H; Pdx1-Cre) pancreatic adenocarcinoma (MT4) model and a dense epithelial neu deletion (NDL) HER2+ mammary adenocarcinoma model with a greater fraction of lymphocytes, natural killer cells and activated dendritic cells. We then performed gene ontology analysis, spectral and digital cytometry to assess the immune response to combination immunotherapies and correlation with survival studies. Result: Based on gene ontology analysis, adding ablation to immunotherapy enriched immune cell migration pathways in the pancreatic cancer model and extensively enriched wound healing pathways in the breast cancer model. With CIBERSORTx digital cytometry, aCD40 + aPD-1 immunotherapy combinations enhanced dendritic cell activation in both models. In the MT4 TME, adding the combination of aCD40 antibody and checkpoint inhibitors (aPD-1 and aCTLA-4) with ablation was synergistic, increasing activated natural killer cells and T cells in distant tumors. Furthermore, ablation with immunotherapy upregulated critical Ly6c myeloid remodeling phenotypes that enhance T-cell effector function and increased granzyme and protease encoding genes by as much as 100-fold. Ablation combined with immunotherapy then extended survival in the MT4 model to a greater extent than immunotherapy alone. Conclusion: In summary, TME profiling informed a successful multicomponent treatment protocol incorporating ablation and facilitated differentiation of TMEs in which ablation is most effective.

Identification of a minority population of LMO2+ breast cancer cells that integrate into the vasculature and initiate metastasis.

Metastasis is responsible for most breast cancer-related deaths; however, identifying the cellular determinants of metastasis has remained challenging. Here, we identified a minority population of immature THY1+/VEGFA+ tumor epithelial cells in human breast tumor biopsies that display angiogenic features and are marked by the expression of the oncogene, LMO2. Higher abundance of LMO2+ basal cells correlated with tumor endothelial content and predicted poor distant recurrence-free survival in patients. Using MMTV-PyMT/Lmo2CreERT2 mice, we demonstrated that Lmo2 lineage-traced cells integrate into the vasculature and have a higher propensity to metastasize. LMO2 knockdown in human breast tumors reduced lung metastasis by impairing intravasation, leading to a reduced frequency of circulating tumor cells. Mechanistically, we find that LMO2 binds to STAT3 and is required for STAT3 activation by tumor necrosis factor-α and interleukin-6. Collectively, our study identifies a population of metastasis-initiating cells with angiogenic features and establishes the LMO2-STAT3 signaling axis as a therapeutic target in breast cancer metastasis.

Inferring gene expression from cell-free DNA fragmentation profiles.

Profiling of circulating tumor DNA (ctDNA) in the bloodstream shows promise for noninvasive cancer detection. Chromatin fragmentation features have previously been explored to infer gene expression profiles from cell-free DNA (cfDNA), but current fragmentomic methods require high concentrations of tumor-derived DNA and provide limited resolution. Here we describe promoter fragmentation entropy as an epigenomic cfDNA feature that predicts RNA expression levels at individual genes. We developed 'epigenetic expression inference from cell-free DNA-sequencing' (EPIC-seq), a method that uses targeted sequencing of promoters of genes of interest. Profiling 329 blood samples from 201 patients with cancer and 87 healthy adults, we demonstrate classification of subtypes of lung carcinoma and diffuse large B cell lymphoma. Applying EPIC-seq to serial blood samples from patients treated with PD-(L)1 immune-checkpoint inhibitors, we show that gene expression profiles inferred by EPIC-seq are correlated with clinical response. Our results indicate that EPIC-seq could enable noninvasive, high-throughput tissue-of-origin characterization with diagnostic, prognostic and therapeutic potential.

Integrative molecular and clinical profiling of acral melanoma links focal amplification of 22q11.21 to metastasis.

Acral melanoma, the most common melanoma subtype among non-White individuals, is associated with poor prognosis. However, its key molecular drivers remain obscure. Here, we perform integrative genomic and clinical profiling of acral melanomas from 104 patients treated in North America (n = 37) or China (n = 67). We find that recurrent, late-arising focal amplifications of cytoband 22q11.21 are a leading determinant of inferior survival, strongly associated with metastasis, and linked to downregulation of immunomodulatory genes associated with response to immune checkpoint blockade. Unexpectedly, LZTR1 - a known tumor suppressor in other cancers - is a key candidate oncogene in this cytoband. Silencing of LZTR1 in melanoma cell lines causes apoptotic cell death independent of major hotspot mutations or melanoma subtypes. Conversely, overexpression of LZTR1 in normal human melanocytes initiates processes associated with metastasis, including anchorage-independent growth, formation of spheroids, and an increase in MAPK and SRC activities. Our results provide insights into the etiology of acral melanoma and implicate LZTR1 as a key tumor promoter and therapeutic target.

T cell characteristics associated with toxicity to immune checkpoint blockade in patients with melanoma.

Severe immune-related adverse events (irAEs) occur in up to 60% of patients with melanoma treated with immune checkpoint inhibitors (ICIs). However, it is unknown whether a common baseline immunological state precedes irAE development. Here we applied mass cytometry by time of flight, single-cell RNA sequencing, single-cell V(D)J sequencing, bulk RNA sequencing and bulk T cell receptor (TCR) sequencing to study peripheral blood samples from patients with melanoma treated with anti-PD-1 monotherapy or anti-PD-1 and anti-CTLA-4 combination ICIs. By analyzing 93 pre- and early on-ICI blood samples and 3 patient cohorts (n = 27, 26 and 18), we found that 2 pretreatment factors in circulation-activated CD4 memory T cell abundance and TCR diversity-are associated with severe irAE development regardless of organ system involvement. We also explored on-treatment changes in TCR clonality among patients receiving combination therapy and linked our findings to the severity and timing of irAE onset. These results demonstrate circulating T cell characteristics associated with ICI-induced toxicity, with implications for improved diagnostics and clinical management.

Atlas of clinically distinct cell states and ecosystems across human solid tumors.

Determining how cells vary with their local signaling environment and organize into distinct cellular communities is critical for understanding processes as diverse as development, aging, and cancer. Here we introduce EcoTyper, a machine learning framework for large-scale identification and validation of cell states and multicellular communities from bulk, single-cell, and spatially resolved gene expression data. When applied to 12 major cell lineages across 16 types of human carcinoma, EcoTyper identified 69 transcriptionally defined cell states. Most states were specific to neoplastic tissue, ubiquitous across tumor types, and significantly prognostic. By analyzing cell-state co-occurrence patterns, we discovered ten clinically distinct multicellular communities with unexpectedly strong conservation, including three with myeloid and stromal elements linked to adverse survival, one enriched in normal tissue, and two associated with early cancer development. This study elucidates fundamental units of cellular organization in human carcinoma and provides a framework for large-scale profiling of cellular ecosystems in any tissue.

The landscape of tumor cell states and ecosystems in diffuse large B cell lymphoma.

Biological heterogeneity in diffuse large B cell lymphoma (DLBCL) is partly driven by cell-of-origin subtypes and associated genomic lesions, but also by diverse cell types and cell states in the tumor microenvironment (TME). However, dissecting these cell states and their clinical relevance at scale remains challenging. Here, we implemented EcoTyper, a machine-learning framework integrating transcriptome deconvolution and single-cell RNA sequencing, to characterize clinically relevant DLBCL cell states and ecosystems. Using this approach, we identified five cell states of malignant B cells that vary in prognostic associations and differentiation status. We also identified striking variation in cell states for 12 other lineages comprising the TME and forming cell state interactions in stereotyped ecosystems. While cell-of-origin subtypes have distinct TME composition, DLBCL ecosystems capture clinical heterogeneity within existing subtypes and extend beyond cell-of-origin and genotypic classes. These results resolve the DLBCL microenvironment at systems-level resolution and identify opportunities for therapeutic targeting (https://ecotyper.stanford.edu/lymphoma).

Integrated spatial multiomics reveals fibroblast fate during tissue repair.

In the skin, tissue injury results in fibrosis in the form of scars composed of dense extracellular matrix deposited by fibroblasts. The therapeutic goal of regenerative wound healing has remained elusive, in part because principles of fibroblast programming and adaptive response to injury remain incompletely understood. Here, we present a multimodal -omics platform for the comprehensive study of cell populations in complex tissue, which has allowed us to characterize the cells involved in wound healing across both time and space. We employ a stented wound model that recapitulates human tissue repair kinetics and multiple Rainbow transgenic lines to precisely track fibroblast fate during the physiologic response to skin injury. Through integrated analysis of single cell chromatin landscapes and gene expression states, coupled with spatial transcriptomic profiling, we are able to impute fibroblast epigenomes with temporospatial resolution. This has allowed us to reveal potential mechanisms controlling fibroblast fate during migration, proliferation, and differentiation following skin injury, and thereby reexamine the canonical phases of wound healing. These findings have broad implications for the study of tissue repair in complex organ systems.