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.

RNA aggregates harness the danger response for potent cancer immunotherapy.

Cancer immunotherapy remains limited by poor antigenicity and a regulatory tumor microenvironment (TME). Here, we create "onion-like" multi-lamellar RNA lipid particle aggregates (LPAs) to substantially enhance the payload packaging and immunogenicity of tumor mRNA antigens. Unlike current mRNA vaccine designs that rely on payload packaging into nanoparticle cores for Toll-like receptor engagement in immune cells, systemically administered RNA-LPAs activate RIG-I in stromal cells, eliciting massive cytokine/chemokine response and dendritic cell/lymphocyte trafficking that provokes cancer immunogenicity and mediates rejection of both early- and late-stage murine tumor models. In client-owned canines with terminal gliomas, RNA-LPAs improved survivorship and reprogrammed the TME, which became "hot" within days of a single infusion. In a first-in-human trial, RNA-LPAs elicited rapid cytokine/chemokine release, immune activation/trafficking, tissue-confirmed pseudoprogression, and glioma-specific immune responses in glioblastoma patients. These data support RNA-LPAs as a new technology that simultaneously reprograms the TME while eliciting rapid and enduring cancer immunotherapy.

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.

The local microenvironment drives activation of neutrophils in human brain tumors.

Neutrophils are abundant immune cells in the circulation and frequently infiltrate tumors in substantial numbers. However, their precise functions in different cancer types remain incompletely understood, including in the brain microenvironment. We therefore investigated neutrophils in tumor tissue of glioma and brain metastasis patients, with matched peripheral blood, and herein describe the first in-depth analysis of neutrophil phenotypes and functions in these tissues. Orthogonal profiling strategies in humans and mice revealed that brain tumor-associated neutrophils (TANs) differ significantly from blood neutrophils and have a prolonged lifespan and immune-suppressive and pro-angiogenic capacity. TANs exhibit a distinct inflammatory signature, driven by a combination of soluble inflammatory mediators including tumor necrosis factor alpha (TNF-ɑ) and Ceruloplasmin, which is more pronounced in TANs from brain metastasis versus glioma. Myeloid cells, including tumor-associated macrophages, emerge at the core of this network of pro-inflammatory mediators, supporting the concept of a critical myeloid niche regulating overall immune suppression in human brain tumors.

Glioma progression is shaped by genetic evolution and microenvironment interactions.

The factors driving therapy resistance in diffuse glioma remain poorly understood. To identify treatment-associated cellular and genetic changes, we analyzed RNA and/or DNA sequencing data from the temporally separated tumor pairs of 304 adult patients with isocitrate dehydrogenase (IDH)-wild-type and IDH-mutant glioma. Tumors recurred in distinct manners that were dependent on IDH mutation status and attributable to changes in histological feature composition, somatic alterations, and microenvironment interactions. Hypermutation and acquired CDKN2A deletions were associated with an increase in proliferating neoplastic cells at recurrence in both glioma subtypes, reflecting active tumor growth. IDH-wild-type tumors were more invasive at recurrence, and their neoplastic cells exhibited increased expression of neuronal signaling programs that reflected a possible role for neuronal interactions in promoting glioma progression. Mesenchymal transition was associated with the presence of a myeloid cell state defined by specific ligand-receptor interactions with neoplastic cells. Collectively, these recurrence-associated phenotypes represent potential targets to alter disease progression.

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.

Inhibitory CD161 receptor identified in glioma-infiltrating T cells by single-cell analysis.

T cells are critical effectors of cancer immunotherapies, but little is known about their gene expression programs in diffuse gliomas. Here, we leverage single-cell RNA sequencing (RNA-seq) to chart the gene expression and clonal landscape of tumor-infiltrating T cells across 31 patients with isocitrate dehydrogenase (IDH) wild-type glioblastoma and IDH mutant glioma. We identify potential effectors of anti-tumor immunity in subsets of T cells that co-express cytotoxic programs and several natural killer (NK) cell genes. Analysis of clonally expanded tumor-infiltrating T cells further identifies the NK gene KLRB1 (encoding CD161) as a candidate inhibitory receptor. Accordingly, genetic inactivation of KLRB1 or antibody-mediated CD161 blockade enhances T cell-mediated killing of glioma cells in vitro and their anti-tumor function in vivo. KLRB1 and its associated transcriptional program are also expressed by substantial T cell populations in other human cancers. Our work provides an atlas of T cells in gliomas and highlights CD161 and other NK cell receptors as immunotherapy targets.

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.

Interrogation of the Microenvironmental Landscape in Brain Tumors Reveals Disease-Specific Alterations of Immune Cells.

Brain malignancies encompass a range of primary and metastatic cancers, including low-grade and high-grade gliomas and brain metastases (BrMs) originating from diverse extracranial tumors. Our understanding of the brain tumor microenvironment (TME) remains limited, and it is unknown whether it is sculpted differentially by primary versus metastatic disease. We therefore comprehensively analyzed the brain TME landscape via flow cytometry, RNA sequencing, protein arrays, culture assays, and spatial tissue characterization. This revealed disease-specific enrichment of immune cells with pronounced differences in proportional abundance of tissue-resident microglia, infiltrating monocyte-derived macrophages, neutrophils, and T cells. These integrated analyses also uncovered multifaceted immune cell activation within brain malignancies entailing converging transcriptional trajectories while maintaining disease- and cell-type-specific programs. Given the interest in developing TME-targeted therapies for brain malignancies, this comprehensive resource of the immune landscape offers insights into possible strategies to overcome tumor-supporting TME properties and instead harness the TME to fight cancer.

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.

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.

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.

CTLA-4 blockade induces a microglia-Th1 cell partnership that stimulates microglia phagocytosis and anti-tumor function in glioblastoma.

The limited efficacy of immunotherapies against glioblastoma underscores the urgency of better understanding immunity in the central nervous system. We found that treatment with αCTLA-4, but not αPD-1, prolonged survival in a mouse model of mesenchymal-like glioblastoma. This effect was lost upon the depletion of CD4+ T cells but not CD8+ T cells. αCTLA-4 treatment increased frequencies of intratumoral IFNγ-producing CD4+ T cells, and IFNγ blockade negated the therapeutic impact of αCTLA-4. The anti-tumor activity of CD4+ T cells did not require tumor-intrinsic MHC-II expression but rather required conventional dendritic cells as well as MHC-II expression on microglia. CD4+ T cells interacted directly with microglia, promoting IFNγ-dependent microglia activation and phagocytosis via the AXL/MER tyrosine kinase receptors, which were necessary for tumor suppression. Thus, αCTLA-4 blockade in mesenchymal-like glioblastoma promotes a CD4+ T cell-microglia circuit wherein IFNγ triggers microglia activation and phagocytosis and microglia in turn act as antigen-presenting cells fueling the CD4+ T cell response.

Neural Precursor-Derived Pleiotrophin Mediates Subventricular Zone Invasion by Glioma.

The lateral ventricle subventricular zone (SVZ) is a frequent and consequential site of pediatric and adult glioma spread, but the cellular and molecular mechanisms mediating this are poorly understood. We demonstrate that neural precursor cell (NPC):glioma cell communication underpins this propensity of glioma to colonize the SVZ through secretion of chemoattractant signals toward which glioma cells home. Biochemical, proteomic, and functional analyses of SVZ NPC-secreted factors revealed the neurite outgrowth-promoting factor pleiotrophin, along with required binding partners SPARC/SPARCL1 and HSP90B, as key mediators of this chemoattractant effect. Pleiotrophin expression is strongly enriched in the SVZ, and pleiotrophin knock down starkly reduced glioma invasion of the SVZ in the murine brain. Pleiotrophin, in complex with the binding partners, activated glioma Rho/ROCK signaling, and ROCK inhibition decreased invasion toward SVZ NPC-secreted factors. These findings demonstrate a pathogenic role for NPC:glioma interactions and potential therapeutic targets to limit glioma invasion. PAPERCLIP.

Sequencing of N6-methyl-deoxyadenosine at single-base resolution across the mammalian genome.

Although DNA N6-methyl-deoxyadenosine (6mA) is abundant in bacteria and protists, its presence and function in mammalian genomes have been less clear. We present Direct-Read 6mA sequencing (DR-6mA-seq), an antibody-independent method, to measure 6mA at base resolution. DR-6mA-seq employs a unique mutation-based strategy to reveal 6mA sites as misincorporation signatures without any chemical or enzymatic modulation of 6mA. We validated DR-6mA-seq through the successful mapping of the well-characterized G(6mA)TC motif in the E. coli DNA. As expected, when applying DR-6mA-seq to mammalian systems, we found that genomic DNA (gDNA) 6mA abundance is generally low in most mammalian tissues and cells; however, we did observe distinct gDNA 6mA sites in mouse testis and glioblastoma cells. DR-6mA-seq provides an enabling tool to detect 6mA at single-base resolution for a comprehensive understanding of DNA 6mA in eukaryotes.

Haploinsufficiency of phosphodiesterase 10A activates PI3K/AKT signaling independent of PTEN to induce an aggressive glioma phenotype.

Glioblastoma is universally fatal and characterized by frequent chromosomal copy number alterations harboring oncogenes and tumor suppressors. In this study, we analyzed exome-wide human glioblastoma copy number data and found that cytoband 6q27 is an independent poor prognostic marker in multiple data sets. We then combined CRISPR-Cas9 data, human spatial transcriptomic data, and human and mouse RNA sequencing data to nominate PDE10A as a potential haploinsufficient tumor suppressor in the 6q27 region. Mouse glioblastoma modeling using the RCAS/tv-a system confirmed that Pde10a suppression induced an aggressive glioma phenotype in vivo and resistance to temozolomide and radiation therapy in vitro. Cell culture analysis showed that decreased Pde10a expression led to increased PI3K/AKT signaling in a Pten-independent manner, a response blocked by selective PI3K inhibitors. Single-nucleus RNA sequencing from our mouse gliomas in vivo, in combination with cell culture validation, further showed that Pde10a suppression was associated with a proneural-to-mesenchymal transition that exhibited increased cell adhesion and decreased cell migration. Our results indicate that glioblastoma patients harboring PDE10A loss have worse outcomes and potentially increased sensitivity to PI3K inhibition.

Glioblastoma microenvironment-from biology to therapy.

Glioblastoma (GBM) is the most aggressive primary brain cancer. These tumors exhibit high intertumoral and intratumoral heterogeneity in neoplastic and nonneoplastic compartments, low lymphocyte infiltration, and high abundance of myeloid subsets that together create a highly protumorigenic immunosuppressive microenvironment. Moreover, heterogeneous GBM cells infiltrate adjacent brain tissue, remodeling the neural microenvironment to foster tumor electrochemical coupling with neurons and metabolic coupling with nonneoplastic astrocytes, thereby driving growth. Here, we review heterogeneity in the GBM microenvironment and its role in low-to-high-grade glioma transition, concluding with a discussion of the challenges of therapeutically targeting the tumor microenvironment and outlining future research opportunities.