Time-resolved single-cell transcriptomics defines immune trajectories in glioblastoma.

Deciphering the cell-state transitions underlying immune adaptation across time is fundamental for advancing biology. Empirical in vivo genomic technologies that capture cellular dynamics are currently lacking. We present Zman-seq, a single-cell technology recording transcriptomic dynamics across time by introducing time stamps into circulating immune cells, tracking them in tissues for days. Applying Zman-seq resolved cell-state and molecular trajectories of the dysfunctional immune microenvironment in glioblastoma. Within 24 hours of tumor infiltration, cytotoxic natural killer cells transitioned to a dysfunctional program regulated by TGFB1 signaling. Infiltrating monocytes differentiated into immunosuppressive macrophages, characterized by the upregulation of suppressive myeloid checkpoints Trem2, Il18bp, and Arg1, over 36 to 48 hours. Treatment with an antagonistic anti-TREM2 antibody reshaped the tumor microenvironment by redirecting the monocyte trajectory toward pro-inflammatory macrophages. Zman-seq is a broadly applicable technology, enabling empirical measurements of differentiation trajectories, which can enhance the development of more efficacious immunotherapies.

Glioblastoma evolution and heterogeneity from a 3D whole-tumor perspective.

Treatment failure for the lethal brain tumor glioblastoma (GBM) is attributed to intratumoral heterogeneity and tumor evolution. We utilized 3D neuronavigation during surgical resection to acquire samples representing the whole tumor mapped by 3D spatial coordinates. Integrative tissue and single-cell analysis revealed sources of genomic, epigenomic, and microenvironmental intratumoral heterogeneity and their spatial patterning. By distinguishing tumor-wide molecular features from those with regional specificity, we inferred GBM evolutionary trajectories from neurodevelopmental lineage origins and initiating events such as chromothripsis to emergence of genetic subclones and spatially restricted activation of differential tumor and microenvironmental programs in the core, periphery, and contrast-enhancing regions. Our work depicts GBM evolution and heterogeneity from a 3D whole-tumor perspective, highlights potential therapeutic targets that might circumvent heterogeneity-related failures, and establishes an interactive platform enabling 360° visualization and analysis of 3D spatial patterns for user-selected genes, programs, and other features across whole GBM tumors.

A holistic view of the malignant organism we call glioblastoma.

Tumors are not simply a chaotic mass of mutated cells but can follow complex organizational principles, including in space. In this issue of Cell, Mathur and colleagues reconstruct a 3D genomic, epigenomic, and transcriptomic spatial cartograph of glioblastoma, offering a "whole-tumor" perspective with patterns of clonal expansion that are embedded in neurodevelopmental hierarchy.

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.

Glioblastoma hijacks neuronal mechanisms for brain invasion.

Glioblastomas are incurable tumors infiltrating the brain. A subpopulation of glioblastoma cells forms a functional and therapy-resistant tumor cell network interconnected by tumor microtubes (TMs). Other subpopulations appear unconnected, and their biological role remains unclear. Here, we demonstrate that whole-brain colonization is fueled by glioblastoma cells that lack connections with other tumor cells and astrocytes yet receive synaptic input from neurons. This subpopulation corresponds to neuronal and neural-progenitor-like tumor cell states, as defined by single-cell transcriptomics, both in mouse models and in the human disease. Tumor cell invasion resembled neuronal migration mechanisms and adopted a Lévy-like movement pattern of probing the environment. Neuronal activity induced complex calcium signals in glioblastoma cells followed by the de novo formation of TMs and increased invasion speed. Collectively, superimposing molecular and functional single-cell data revealed that neuronal mechanisms govern glioblastoma cell invasion on multiple levels. This explains how glioblastoma's dissemination and cellular heterogeneity are closely interlinked.

Invasive glioma cells: The malignant pioneers that follow the current.

Glioblastoma is a lethal, diffusely invasive brain cancer that is robustly regulated by the activity of the brain itself, in part through neuron-to-glioma synaptic communication. Venkataramani et al. have conceptually advanced understanding of glioblastoma interactions with neural circuits, demonstrating that conduction of electrochemical signals via neuron-to-glioma synapses drives glioma invasion.

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.

A Patient-Derived Glioblastoma Organoid Model and Biobank Recapitulates Inter- and Intra-tumoral Heterogeneity.

Glioblastomas exhibit vast inter- and intra-tumoral heterogeneity, complicating the development of effective therapeutic strategies. Current in vitro models are limited in preserving the cellular and mutational diversity of parental tumors and require a prolonged generation time. Here, we report methods for generating and biobanking patient-derived glioblastoma organoids (GBOs) that recapitulate the histological features, cellular diversity, gene expression, and mutational profiles of their corresponding parental tumors. GBOs can be generated quickly with high reliability and exhibit rapid, aggressive infiltration when transplanted into adult rodent brains. We further demonstrate the utility of GBOs to test personalized therapies by correlating GBO mutational profiles with responses to specific drugs and by modeling chimeric antigen receptor T cell immunotherapy. Our studies show that GBOs maintain many key features of glioblastomas and can be rapidly deployed to investigate patient-specific treatment strategies. Additionally, our live biobank establishes a rich resource for basic and translational glioblastoma research.

An Integrative Model of Cellular States, Plasticity, and Genetics for Glioblastoma.

Diverse genetic, epigenetic, and developmental programs drive glioblastoma, an incurable and poorly understood tumor, but their precise characterization remains challenging. Here, we use an integrative approach spanning single-cell RNA-sequencing of 28 tumors, bulk genetic and expression analysis of 401 specimens from the The Cancer Genome Atlas (TCGA), functional approaches, and single-cell lineage tracing to derive a unified model of cellular states and genetic diversity in glioblastoma. We find that malignant cells in glioblastoma exist in four main cellular states that recapitulate distinct neural cell types, are influenced by the tumor microenvironment, and exhibit plasticity. The relative frequency of cells in each state varies between glioblastoma samples and is influenced by copy number amplifications of the CDK4, EGFR, and PDGFRA loci and by mutations in the NF1 locus, which each favor a defined state. Our work provides a blueprint for glioblastoma, integrating the malignant cell programs, their plasticity, and their modulation by genetic drivers.

Functional Enhancers Shape Extrachromosomal Oncogene Amplifications.

Non-coding regions amplified beyond oncogene borders have largely been ignored. Using a computational approach, we find signatures of significant co-amplification of non-coding DNA beyond the boundaries of amplified oncogenes across five cancer types. In glioblastoma, EGFR is preferentially co-amplified with its two endogenous enhancer elements active in the cell type of origin. These regulatory elements, their contacts, and their contribution to cell fitness are preserved on high-level circular extrachromosomal DNA amplifications. Interrogating the locus with a CRISPR interference screening approach reveals a diversity of additional elements that impact cell fitness. The pattern of fitness dependencies mirrors the rearrangement of regulatory elements and accompanying rewiring of the chromatin topology on the extrachromosomal amplicon. Our studies indicate that oncogene amplifications are shaped by regulatory dependencies in the non-coding genome.

Tumor hypoxia represses γδ T cell-mediated antitumor immunity against brain tumors.

The anatomic location and immunologic characteristics of brain tumors result in strong lymphocyte suppression. Consequently, conventional immunotherapies targeting CD8 T cells are ineffective against brain tumors. Tumor cells escape immunosurveillance by various mechanisms and tumor cell metabolism can affect the metabolic states and functions of tumor-infiltrating lymphocytes. Here, we discovered that brain tumor cells had a particularly high demand for oxygen, which affected γδ T cell-mediated antitumor immune responses but not those of conventional T cells. Specifically, tumor hypoxia activated the γδ T cell protein kinase A pathway at a transcriptional level, resulting in repression of the activatory receptor NKG2D. Alleviating tumor hypoxia reinvigorated NKG2D expression and the antitumor function of γδ T cells. These results reveal a hypoxia-mediated mechanism through which brain tumors and γδ T cells interact and emphasize the importance of γδ T cells for antitumor immunity against brain tumors.

Glucose-driven histone lactylation promotes the immunosuppressive activity of monocyte-derived macrophages in glioblastoma.

Immunosuppressive macrophages restrict anti-cancer immunity in glioblastoma (GBM). Here, we studied the contribution of microglia (MGs) and monocyte-derived macrophages (MDMs) to immunosuppression and mechanisms underlying their regulatory function. MDMs outnumbered MGs at late tumor stages and suppressed T cell activity. Molecular and functional analysis identified a population of glycolytic MDM expressing GLUT1 with potent immunosuppressive activity. GBM-derived factors promoted high glycolysis, lactate, and interleukin-10 (IL-10) production in MDMs. Inhibition of glycolysis or lactate production in MDMs impaired IL-10 expression and T cell suppression. Mechanistically, intracellular lactate-driven histone lactylation promoted IL-10 expression, which was required to suppress T cell activity. GLUT1 expression on MDMs was induced downstream of tumor-derived factors that activated the PERK-ATF4 axis. PERK deletion in MDM abrogated histone lactylation, led to the accumulation of intratumoral T cells and tumor growth delay, and, in combination with immunotherapy, blocked GBM progression. Thus, PERK-driven glucose metabolism promotes MDM immunosuppressive activity via histone lactylation.

Natural Killer Cells Control Tumor Growth by Sensing a Growth Factor.

Many tumors produce platelet-derived growth factor (PDGF)-DD, which promotes cellular proliferation, epithelial-mesenchymal transition, stromal reaction, and angiogenesis through autocrine and paracrine PDGFRß signaling. By screening a secretome library, we found that the human immunoreceptor NKp44, encoded by NCR2 and expressed on natural killer (NK) cells and innate lymphoid cells, recognizes PDGF-DD. PDGF-DD engagement of NKp44 triggered NK cell secretion of interferon gamma (IFN)-γ and tumor necrosis factor alpha (TNF-α) that induced tumor cell growth arrest. A distinctive transcriptional signature of PDGF-DD-induced cytokines and the downregulation of tumor cell-cycle genes correlated with NCR2 expression and greater survival in glioblastoma. NKp44 expression in mouse NK cells controlled the dissemination of tumors expressing PDGF-DD more effectively than control mice, an effect enhanced by blockade of the inhibitory receptor CD96 or CpG-oligonucleotide treatment. Thus, while cancer cell production of PDGF-DD supports tumor growth and stromal reaction, it concomitantly activates innate immune responses to tumor expansion.

N6-methyladenine DNA Modification in Glioblastoma.

Genetic drivers of cancer can be dysregulated through epigenetic modifications of DNA. Although the critical role of DNA 5-methylcytosine (5mC) in the regulation of transcription is recognized, the functions of other non-canonical DNA modifications remain obscure. Here, we report the identification of novel N6-methyladenine (N6-mA) DNA modifications in human tissues and implicate this epigenetic mark in human disease, specifically the highly malignant brain cancer glioblastoma. Glioblastoma markedly upregulated N6-mA levels, which co-localized with heterochromatic histone modifications, predominantly H3K9me3. N6-mA levels were dynamically regulated by the DNA demethylase ALKBH1, depletion of which led to transcriptional silencing of oncogenic pathways through decreasing chromatin accessibility. Targeting the N6-mA regulator ALKBH1 in patient-derived human glioblastoma models inhibited tumor cell proliferation and extended the survival of tumor-bearing mice, supporting this novel DNA modification as a potential therapeutic target for glioblastoma. Collectively, our results uncover a novel epigenetic node in cancer through the DNA modification N6-mA.

Hypoxic niches attract and sequester tumor-associated macrophages and cytotoxic T cells and reprogram them for immunosuppression.

Glioblastoma (GBM), a highly lethal brain cancer, is notorious for immunosuppression, but the mechanisms remain unclear. Here, we documented a temporospatial patterning of tumor-associated myeloid cells (TAMs) corresponding to vascular changes during GBM progression. As tumor vessels transitioned from the initial dense regular network to later scant and engorged vasculature, TAMs shifted away from perivascular regions and trafficked to vascular-poor areas. This process was heavily influenced by the immunocompetence state of the host. Utilizing a sensitive fluorescent UnaG reporter to track tumor hypoxia, coupled with single-cell transcriptomics, we revealed that hypoxic niches attracted and sequestered TAMs and cytotoxic T lymphocytes (CTLs), where they were reprogrammed toward an immunosuppressive state. Mechanistically, we identified chemokine CCL8 and cytokine IL-1ß as two hypoxic-niche factors critical for TAM trafficking and co-evolution of hypoxic zones into pseudopalisading patterns. Therefore, perturbation of TAM patterning in hypoxic zones may improve tumor control.

Plexin-B2 Mediates Physiologic and Pathologic Functions of Angiogenin.

Angiogenin (ANG) is a secreted ribonuclease (RNase) with cell-type- and context-specific roles in growth, survival, and regeneration. Although these functions require receptor-mediated endocytosis and appropriate subcellular localization, the identity of the cell surface receptor remains undefined. Here, we show that plexin-B2 (PLXNB2) is the functional receptor for ANG in endothelial, cancer, neuronal, and normal hematopoietic and leukemic stem and progenitor cells. Mechanistically, PLXNB2 mediates intracellular RNA processing that contribute to cell growth, survival, and regenerative capabilities of ANG. Antibodies generated against the ANG-binding site on PLXNB2 restricts ANG activity in vitro and in vivo, resulting in inhibition of established xenograft tumors, ANG-induced neurogenesis and neuroprotection, levels of pro-self-renewal transcripts in hematopoietic and patient-derived leukemic stem and progenitor cells, and reduced progression of leukemia in vivo. PLXNB2 is therefore required for the physiological and pathological functions of ANG and has significant therapeutic potential in solid and hematopoietic cancers and neurodegenerative diseases.

Epigenetic Activation of WNT5A Drives Glioblastoma Stem Cell Differentiation and Invasive Growth.

Glioblastoma stem cells (GSCs) are implicated in tumor neovascularization, invasiveness, and therapeutic resistance. To illuminate mechanisms governing these hallmark features, we developed a de novo glioblastoma multiforme (GBM) model derived from immortalized human neural stem/progenitor cells (hNSCs) to enable precise system-level comparisons of pre-malignant and oncogene-induced malignant states of NSCs. Integrated transcriptomic and epigenomic analyses uncovered a PAX6/DLX5 transcriptional program driving WNT5A-mediated GSC differentiation into endothelial-like cells (GdECs). GdECs recruit existing endothelial cells to promote peritumoral satellite lesions, which serve as a niche supporting the growth of invasive glioma cells away from the primary tumor. Clinical data reveal higher WNT5A and GdECs expression in peritumoral and recurrent GBMs relative to matched intratumoral and primary GBMs, respectively, supporting WNT5A-mediated GSC differentiation and invasive growth in disease recurrence. Thus, the PAX6/DLX5-WNT5A axis governs the diffuse spread of glioma cells throughout the brain parenchyma, contributing to the lethality of GBM.

Lactate is an epigenetic metabolite that drives survival in model systems of glioblastoma.

Lactate accumulates to a significant amount in glioblastomas (GBMs), the most common primary malignant brain tumor with an unfavorable prognosis. However, it remains unclear whether lactate is metabolized by GBMs. Here, we demonstrated that lactate rescued patient-derived xenograft (PDX) GBM cells from nutrient-deprivation-mediated cell death. Transcriptome analysis, ATAC-seq, and ChIP-seq showed that lactate entertained a signature of oxidative energy metabolism. LC/MS analysis demonstrated that U-13C-lactate elicited substantial labeling of TCA-cycle metabolites, acetyl-CoA, and histone protein acetyl-residues in GBM cells. Lactate enhanced chromatin accessibility and histone acetylation in a manner dependent on oxidative energy metabolism and the ATP-citrate lyase (ACLY). Utilizing orthotopic PDX models of GBM, a combined tracer experiment unraveled that lactate carbons were substantially labeling the TCA-cycle metabolites. Finally, pharmacological blockage of oxidative energy metabolism extended overall survival in two orthotopic PDX models in mice. These results establish lactate metabolism as a novel druggable pathway for GBM.

Connect four with glioblastoma stem cell factors.

Hierarchical cell state models, wherein a few stem-like tumor-propagating cells repopulate the tumor after therapy, are often invoked in cancer. Suvà et al. demonstrate a plastic developmental hierarchy in glioma cell populations by characterizing the epigenetic states of phenotypically distinct cells and identifying four factors sufficient to reprogram differentiated cells into a tumorigenic stem-like state.

Mathematical modeling of PDGF-driven glioblastoma reveals optimized radiation dosing schedules.

Glioblastomas (GBMs) are the most common and malignant primary brain tumors and are aggressively treated with surgery, chemotherapy, and radiotherapy. Despite this treatment, recurrence is inevitable and survival has improved minimally over the last 50 years. Recent studies have suggested that GBMs exhibit both heterogeneity and instability of differentiation states and varying sensitivities of these states to radiation. Here, we employed an iterative combined theoretical and experimental strategy that takes into account tumor cellular heterogeneity and dynamically acquired radioresistance to predict the effectiveness of different radiation schedules. Using this model, we identified two delivery schedules predicted to significantly improve efficacy by taking advantage of the dynamic instability of radioresistance. These schedules led to superior survival in mice. Our interdisciplinary approach may also be applicable to other human cancer types treated with radiotherapy and, hence, may lay the foundation for significantly increasing the effectiveness of a mainstay of oncologic therapy. PAPERCLIP: