Integrative spatial analysis reveals a multi-layered organization of glioblastoma.

Glioma contains malignant cells in diverse states. Here, we combine spatial transcriptomics, spatial proteomics, and computational approaches to define glioma cellular states and uncover their organization. We find three prominent modes of organization. First, gliomas are composed of small local environments, each typically enriched with one major cellular state. Second, specific pairs of states preferentially reside in proximity across multiple scales. This pairing of states is consistent across tumors. Third, these pairwise interactions collectively define a global architecture composed of five layers. Hypoxia appears to drive the layers, as it is associated with a long-range organization that includes all cancer cell states. Accordingly, tumor regions distant from any hypoxic/necrotic foci and tumors that lack hypoxia such as low-grade IDH-mutant glioma are less organized. In summary, we provide a conceptual framework for the organization of cellular states in glioma, highlighting hypoxia as a long-range tissue organizer.

Discovering dominant tumor immune archetypes in a pan-cancer census.

Cancers display significant heterogeneity with respect to tissue of origin, driver mutations, and other features of the surrounding tissue. It is likely that individual tumors engage common patterns of the immune system-here "archetypes"-creating prototypical non-destructive tumor immune microenvironments (TMEs) and modulating tumor-targeting. To discover the dominant immune system archetypes, the University of California, San Francisco (UCSF) Immunoprofiler Initiative (IPI) processed 364 individual tumors across 12 cancer types using standardized protocols. Computational clustering of flow cytometry and transcriptomic data obtained from cell sub-compartments uncovered dominant patterns of immune composition across cancers. These archetypes were profound insofar as they also differentiated tumors based upon unique immune and tumor gene-expression patterns. They also partitioned well-established classifications of tumor biology. The IPI resource provides a template for understanding cancer immunity as a collection of dominant patterns of immune organization and provides a rational path forward to learn how to modulate these to improve therapy.

Pan-cancer profiling of tumor-infiltrating natural killer cells through transcriptional reference mapping.

The functional diversity of natural killer (NK) cell repertoires stems from differentiation, homeostatic, receptor-ligand interactions and adaptive-like responses to viral infections. In the present study, we generated a single-cell transcriptional reference map of healthy human blood- and tissue-derived NK cells, with temporal resolution and fate-specific expression of gene-regulatory networks defining NK cell differentiation. Transfer learning facilitated incorporation of tumor-infiltrating NK cell transcriptomes (39 datasets, 7 solid tumors, 427 patients) into the reference map to analyze tumor microenvironment (TME)-induced perturbations. Of the six functionally distinct NK cell states identified, a dysfunctional stressed CD56bright state susceptible to TME-induced immunosuppression and a cytotoxic TME-resistant effector CD56dim state were commonly enriched across tumor types, the ratio of which was predictive of patient outcome in malignant melanoma and osteosarcoma. This resource may inform the design of new NK cell therapies and can be extended through transfer learning to interrogate new datasets from experimental perturbations or disease conditions.

ZBTB46 coordinates angiogenesis and immunity to control tumor outcome.

Tumor angiogenesis and immunity show an inverse correlation in cancer progression and outcome1. Here, we report that ZBTB46, a repressive transcription factor and a widely accepted marker for classical dendritic cells (DCs)2,3, controls both tumor angiogenesis and immunity. Zbtb46 was downregulated in both DCs and endothelial cells by tumor-derived factors to facilitate robust tumor growth. Zbtb46 downregulation led to a hallmark pro-tumor microenvironment (TME), including dysfunctional vasculature and immunosuppressive conditions. Analysis of human cancer data revealed a similar association of low ZBTB46 expression with an immunosuppressive TME and a worse prognosis. In contrast, enforced Zbtb46 expression led to TME changes to restrict tumor growth. Mechanistically, Zbtb46-deficient endothelial cells were highly angiogenic, and Zbtb46-deficient bone marrow progenitors upregulated Cebpb and diverted the DC program to immunosuppressive myeloid lineage output, potentially explaining the myeloid lineage skewing phenomenon in cancer4. Conversely, enforced Zbtb46 expression normalized tumor vessels and, by suppressing Cebpb, skewed bone marrow precursors toward immunostimulatory myeloid lineage output, leading to an immune-hot TME. Remarkably, Zbtb46 mRNA treatment synergized with anti-PD1 immunotherapy to improve tumor management in preclinical models. These findings identify ZBTB46 as a critical factor for angiogenesis and for myeloid lineage skewing in cancer and suggest that maintaining its expression could have therapeutic benefits.

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.

Microenvironment drives cell state, plasticity, and drug response in pancreatic cancer.

Prognostically relevant RNA expression states exist in pancreatic ductal adenocarcinoma (PDAC), but our understanding of their drivers, stability, and relationship to therapeutic response is limited. To examine these attributes systematically, we profiled metastatic biopsies and matched organoid models at single-cell resolution. In vivo, we identify a new intermediate PDAC transcriptional cell state and uncover distinct site- and state-specific tumor microenvironments (TMEs). Benchmarking models against this reference map, we reveal strong culture-specific biases in cancer cell transcriptional state representation driven by altered TME signals. We restore expression state heterogeneity by adding back in vivo-relevant factors and show plasticity in culture models. Further, we prove that non-genetic modulation of cell state can strongly influence drug responses, uncovering state-specific vulnerabilities. This work provides a broadly applicable framework for aligning cell states across in vivo and ex vivo settings, identifying drivers of transcriptional plasticity and manipulating cell state to target associated vulnerabilities.

A modular master regulator landscape controls cancer transcriptional identity.

Despite considerable efforts, the mechanisms linking genomic alterations to the transcriptional identity of cancer cells remain elusive. Integrative genomic analysis, using a network-based approach, identified 407 master regulator (MR) proteins responsible for canalizing the genetics of individual samples from 20 cohorts in The Cancer Genome Atlas (TCGA) into 112 transcriptionally distinct tumor subtypes. MR proteins could be further organized into 24 pan-cancer, master regulator block modules (MRBs), each regulating key cancer hallmarks and predictive of patient outcome in multiple cohorts. Of all somatic alterations detected in each individual sample, >50% were predicted to induce aberrant MR activity, yielding insight into mechanisms linking tumor genetics and transcriptional identity and establishing non-oncogene dependencies. Genetic and pharmacological validation assays confirmed the predicted effect of upstream mutations and MR activity on downstream cellular identity and phenotype. Thus, co-analysis of mutational and gene expression profiles identified elusive subtypes and provided testable hypothesis for mechanisms mediating the effect of genetic alterations.

Single-cell landscape of the ecosystem in early-relapse hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) has high relapse and low 5-year survival rates. Single-cell profiling in relapsed HCC may aid in the design of effective anticancer therapies, including immunotherapies. We profiled the transcriptomes of ∼17,000 cells from 18 primary or early-relapse HCC cases. Early-relapse tumors have reduced levels of regulatory T cells, increased dendritic cells (DCs), and increased infiltrated CD8+ T cells, compared with primary tumors, in two independent cohorts. Remarkably, CD8+ T cells in recurrent tumors overexpressed KLRB1 (CD161) and displayed an innate-like low cytotoxic state, with low clonal expansion, unlike the classical exhausted state observed in primary HCC. The enrichment of these cells was associated with a worse prognosis. Differential gene expression and interaction analyses revealed potential immune evasion mechanisms in recurrent tumor cells that dampen DC antigen presentation and recruit innate-like CD8+ T cells. Our comprehensive picture of the HCC ecosystem provides deeper insights into immune evasion mechanisms associated with tumor relapse.

A proteogenomic portrait of lung squamous cell carcinoma.

Lung squamous cell carcinoma (LSCC) remains a leading cause of cancer death with few therapeutic options. We characterized the proteogenomic landscape of LSCC, providing a deeper exposition of LSCC biology with potential therapeutic implications. We identify NSD3 as an alternative driver in FGFR1-amplified tumors and low-p63 tumors overexpressing the therapeutic target survivin. SOX2 is considered undruggable, but our analyses provide rationale for exploring chromatin modifiers such as LSD1 and EZH2 to target SOX2-overexpressing tumors. Our data support complex regulation of metabolic pathways by crosstalk between post-translational modifications including ubiquitylation. Numerous immune-related proteogenomic observations suggest directions for further investigation. Proteogenomic dissection of CDKN2A mutations argue for more nuanced assessment of RB1 protein expression and phosphorylation before declaring CDK4/6 inhibition unsuccessful. Finally, triangulation between LSCC, LUAD, and HNSCC identified both unique and common therapeutic vulnerabilities. These observations and proteogenomics data resources may guide research into the biology and treatment of LSCC.

Differential pre-malignant programs and microenvironment chart distinct paths to malignancy in human colorectal polyps.

Colorectal cancers (CRCs) arise from precursor polyps whose cellular origins, molecular heterogeneity, and immunogenic potential may reveal diagnostic and therapeutic insights when analyzed at high resolution. We present a single-cell transcriptomic and imaging atlas of the two most common human colorectal polyps, conventional adenomas and serrated polyps, and their resulting CRC counterparts. Integrative analysis of 128 datasets from 62 participants reveals adenomas arise from WNT-driven expansion of stem cells, while serrated polyps derive from differentiated cells through gastric metaplasia. Metaplasia-associated damage is coupled to a cytotoxic immune microenvironment preceding hypermutation, driven partly by antigen-presentation differences associated with tumor cell-differentiation status. Microsatellite unstable CRCs contain distinct non-metaplastic regions where tumor cells acquire stem cell properties and cytotoxic immune cells are depleted. Our multi-omic atlas provides insights into malignant progression of colorectal polyps and their microenvironment, serving as a framework for precision surveillance and prevention of CRC.

Spatially organized multicellular immune hubs in human colorectal cancer.

Immune responses to cancer are highly variable, with mismatch repair-deficient (MMRd) tumors exhibiting more anti-tumor immunity than mismatch repair-proficient (MMRp) tumors. To understand the rules governing these varied responses, we transcriptionally profiled 371,223 cells from colorectal tumors and adjacent normal tissues of 28 MMRp and 34 MMRd individuals. Analysis of 88 cell subsets and their 204 associated gene expression programs revealed extensive transcriptional and spatial remodeling across tumors. To discover hubs of interacting malignant and immune cells, we identified expression programs in different cell types that co-varied across tumors from affected individuals and used spatial profiling to localize coordinated programs. We discovered a myeloid cell-attracting hub at the tumor-luminal interface associated with tissue damage and an MMRd-enriched immune hub within the tumor, with activated T cells together with malignant and myeloid cells expressing T cell-attracting chemokines. By identifying interacting cellular programs, we reveal the logic underlying spatially organized immune-malignant cell networks.

Glycogen accumulation and phase separation drives liver tumor initiation.

Glucose consumption is generally increased in tumor cells to support tumor growth. Interestingly, we report that glycogen accumulation is a key initiating oncogenic event during liver malignant transformation. We found that glucose-6-phosphatase (G6PC) catalyzing the last step of glycogenolysis is frequently downregulated to augment glucose storage in pre-malignant cells. Accumulated glycogen undergoes liquid-liquid phase separation, which results in the assembly of the Laforin-Mst1/2 complex and consequently sequesters Hippo kinases Mst1/2 in glycogen liquid droplets to relieve their inhibition on Yap. Moreover, G6PC or another glycogenolysis enzyme-liver glycogen phosphorylase (PYGL) deficiency in both human and mice results in glycogen storage disease along with liver enlargement and tumorigenesis in a Yap-dependent manner. Consistently, elimination of glycogen accumulation abrogates liver growth and cancer incidence, whereas increasing glycogen storage accelerates tumorigenesis. Thus, we concluded that cancer-initiating cells adapt a glycogen storing mode, which blocks Hippo signaling through glycogen phase separation to augment tumor incidence.

FBXO44 promotes DNA replication-coupled repetitive element silencing in cancer cells.

Repetitive elements (REs) compose ∼50% of the human genome and are normally transcriptionally silenced, although the mechanism has remained elusive. Through an RNAi screen, we identified FBXO44 as an essential repressor of REs in cancer cells. FBXO44 bound H3K9me3-modified nucleosomes at the replication fork and recruited SUV39H1, CRL4, and Mi-2/NuRD to transcriptionally silence REs post-DNA replication. FBXO44/SUV39H1 inhibition reactivated REs, leading to DNA replication stress and stimulation of MAVS/STING antiviral pathways and interferon (IFN) signaling in cancer cells to promote decreased tumorigenicity, increased immunogenicity, and enhanced immunotherapy response. FBXO44 expression inversely correlated with replication stress, antiviral pathways, IFN signaling, and cytotoxic T cell infiltration in human cancers, while a FBXO44-immune gene signature correlated with improved immunotherapy response in cancer patients. FBXO44/SUV39H1 were dispensable in normal cells. Collectively, FBXO44/SUV39H1 are crucial repressors of RE transcription, and their inhibition selectively induces DNA replication stress and viral mimicry in cancer cells.

Synthetic lethality-mediated precision oncology via the tumor transcriptome.

Precision oncology has made significant advances, mainly by targeting actionable mutations in cancer driver genes. Aiming to expand treatment opportunities, recent studies have begun to explore the utility of tumor transcriptome to guide patient treatment. Here, we introduce SELECT (synthetic lethality and rescue-mediated precision oncology via the transcriptome), a precision oncology framework harnessing genetic interactions to predict patient response to cancer therapy from the tumor transcriptome. SELECT is tested on a broad collection of 35 published targeted and immunotherapy clinical trials from 10 different cancer types. It is predictive of patients' response in 80% of these clinical trials and in the recent multi-arm WINTHER trial. The predictive signatures and the code are made publicly available for academic use, laying a basis for future prospective clinical studies.

The PP2A-Integrator-CDK9 axis fine-tunes transcription and can be targeted therapeutically in cancer.

Gene expression by RNA polymerase II (RNAPII) is tightly controlled by cyclin-dependent kinases (CDKs) at discrete checkpoints during the transcription cycle. The pausing checkpoint following transcription initiation is primarily controlled by CDK9. We discovered that CDK9-mediated, RNAPII-driven transcription is functionally opposed by a protein phosphatase 2A (PP2A) complex that is recruited to transcription sites by the Integrator complex subunit INTS6. PP2A dynamically antagonizes phosphorylation of key CDK9 substrates including DSIF and RNAPII-CTD. Loss of INTS6 results in resistance to tumor cell death mediated by CDK9 inhibition, decreased turnover of CDK9 phospho-substrates, and amplification of acute oncogenic transcriptional responses. Pharmacological PP2A activation synergizes with CDK9 inhibition to kill both leukemic and solid tumor cells, providing therapeutic benefit in vivo. These data demonstrate that fine control of gene expression relies on the balance between kinase and phosphatase activity throughout the transcription cycle, a process dysregulated in cancer that can be exploited therapeutically.

Spatially confined sub-tumor microenvironments in pancreatic cancer.

Intratumoral heterogeneity is a critical frontier in understanding how the tumor microenvironment (TME) propels malignant progression. Here, we deconvolute the human pancreatic TME through large-scale integration of histology-guided regional multiOMICs with clinical data and patient-derived preclinical models. We discover "subTMEs," histologically definable tissue states anchored in fibroblast plasticity, with regional relationships to tumor immunity, subtypes, differentiation, and treatment response. "Reactive" subTMEs rich in complex but functionally coordinated fibroblast communities were immune hot and inhabited by aggressive tumor cell phenotypes. The matrix-rich "deserted" subTMEs harbored fewer activated fibroblasts and tumor-suppressive features yet were markedly chemoprotective and enriched upon chemotherapy. SubTMEs originated in fibroblast differentiation trajectories, and transitory states were notable both in single-cell transcriptomics and in situ. The intratumoral co-occurrence of subTMEs produced patient-specific phenotypic and computationally predictable heterogeneity tightly linked to malignant biology. Therefore, heterogeneity within the plentiful, notorious pancreatic TME is not random but marks fundamental tissue organizational units.

The bone microenvironment invigorates metastatic seeds for further dissemination.

Metastasis has been considered as the terminal step of tumor progression. However, recent genomic studies suggest that many metastases are initiated by further spread of other metastases. Nevertheless, the corresponding pre-clinical models are lacking, and underlying mechanisms are elusive. Using several approaches, including parabiosis and an evolving barcode system, we demonstrated that the bone microenvironment facilitates breast and prostate cancer cells to further metastasize and establish multi-organ secondary metastases. We uncovered that this metastasis-promoting effect is driven by epigenetic reprogramming that confers stem cell-like properties on cancer cells disseminated from bone lesions. Furthermore, we discovered that enhanced EZH2 activity mediates the increased stemness and metastasis capacity. The same findings also apply to single cell-derived populations, indicating mechanisms distinct from clonal selection. Taken together, our work revealed an unappreciated role of the bone microenvironment in metastasis evolution and elucidated an epigenomic reprogramming process driving terminal-stage, multi-organ metastases.

Glioblastomas acquire myeloid-affiliated transcriptional programs via epigenetic immunoediting to elicit immune evasion.

Glioblastoma multiforme (GBM) is an aggressive brain tumor for which current immunotherapy approaches have been unsuccessful. Here, we explore the mechanisms underlying immune evasion in GBM. By serially transplanting GBM stem cells (GSCs) into immunocompetent hosts, we uncover an acquired capability of GSCs to escape immune clearance by establishing an enhanced immunosuppressive tumor microenvironment. Mechanistically, this is not elicited via genetic selection of tumor subclones, but through an epigenetic immunoediting process wherein stable transcriptional and epigenetic changes in GSCs are enforced following immune attack. These changes launch a myeloid-affiliated transcriptional program, which leads to increased recruitment of tumor-associated macrophages. Furthermore, we identify similar epigenetic and transcriptional signatures in human mesenchymal subtype GSCs. We conclude that epigenetic immunoediting may drive an acquired immune evasion program in the most aggressive mesenchymal GBM subtype by reshaping the tumor immune microenvironment.

Single-cell protein activity analysis identifies recurrence-associated renal tumor macrophages.

Clear cell renal carcinoma (ccRCC) is a heterogeneous disease with a variable post-surgical course. To assemble a comprehensive ccRCC tumor microenvironment (TME) atlas, we performed single-cell RNA sequencing (scRNA-seq) of hematopoietic and non-hematopoietic subpopulations from tumor and tumor-adjacent tissue of treatment-naive ccRCC resections. We leveraged the VIPER algorithm to quantitate single-cell protein activity and validated this approach by comparison to flow cytometry. The analysis identified key TME subpopulations, as well as their master regulators and candidate cell-cell interactions, revealing clinically relevant populations, undetectable by gene-expression analysis. Specifically, we uncovered a tumor-specific macrophage subpopulation characterized by upregulation of TREM2/APOE/C1Q, validated by spatially resolved, quantitative multispectral immunofluorescence. In a large clinical validation cohort, these markers were significantly enriched in tumors from patients who recurred following surgery. The study thus identifies TREM2/APOE/C1Q-positive macrophage infiltration as a potential prognostic biomarker for ccRCC recurrence, as well as a candidate therapeutic target.

ERα is an RNA-binding protein sustaining tumor cell survival and drug resistance.

Estrogen receptor α (ERα) is a hormone receptor and key driver for over 70% of breast cancers that has been studied for decades as a transcription factor. Unexpectedly, we discover that ERα is a potent non-canonical RNA-binding protein. We show that ERα RNA binding function is uncoupled from its activity to bind DNA and critical for breast cancer progression. Employing genome-wide cross-linking immunoprecipitation (CLIP) sequencing and a functional CRISPRi screen, we find that ERα-associated mRNAs sustain cancer cell fitness and elicit cellular responses to stress. Mechanistically, ERα controls different steps of RNA metabolism. In particular, we demonstrate that ERα RNA binding mediates alternative splicing of XBP1 and translation of the eIF4G2 and MCL1 mRNAs, which facilitates survival upon stress conditions and sustains tamoxifen resistance of cancer cells. ERα is therefore a multifaceted RNA-binding protein, and this activity transforms our knowledge of post-transcriptional regulation underlying cancer development and drug response.