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.

Molecular identity of human outer radial glia during cortical development.

Radial glia, the neural stem cells of the neocortex, are located in two niches: the ventricular zone and outer subventricular zone. Although outer subventricular zone radial glia may generate the majority of human cortical neurons, their molecular features remain elusive. By analyzing gene expression across single cells, we find that outer radial glia preferentially express genes related to extracellular matrix formation, migration, and stemness, including TNC, PTPRZ1, FAM107A, HOPX, and LIFR. Using dynamic imaging, immunostaining, and clonal analysis, we relate these molecular features to distinctive behaviors of outer radial glia, demonstrate the necessity of STAT3 signaling for their cell cycle progression, and establish their extensive proliferative potential. These results suggest that outer radial glia directly support the subventricular niche through local production of growth factors, potentiation of growth factor signals by extracellular matrix proteins, and activation of self-renewal pathways, thereby enabling the developmental and evolutionary expansion of the human neocortex.

Systematic identification of barriers to human iPSC generation.

Reprogramming of somatic cells to induced pluripotent stem cells (iPSCs) holds enormous promise for regenerative medicine. To elucidate endogenous barriers limiting this process, we systematically dissected human cellular reprogramming by combining a genome-wide RNAi screen, innovative computational methods, extensive single-hit validation, and mechanistic investigation of relevant pathways and networks. We identify reprogramming barriers, including genes involved in transcription, chromatin regulation, ubiquitination, dephosphorylation, vesicular transport, and cell adhesion. Specific a disintegrin and metalloproteinase (ADAM) proteins inhibit reprogramming, and the disintegrin domain of ADAM29 is necessary and sufficient for this function. Clathrin-mediated endocytosis can be targeted with small molecules and opposes reprogramming by positively regulating TGF-ß signaling. Genetic interaction studies of endocytosis or ubiquitination reveal that barrier pathways can act in linear, parallel, or feedforward loop architectures to antagonize reprogramming. These results provide a global view of barriers to human cellular reprogramming.

miRNA-independent function of long noncoding pri-miRNA loci.

Among the large, diverse set of mammalian long noncoding RNAs (lncRNAs), long noncoding primary microRNAs (lnc-pri-miRNAs) are those that host miRNAs. Whether lnc-pri-miRNA loci have important biological function independent of their cognate miRNAs is poorly understood. From a genome-scale lncRNA screen, lnc-pri-miRNA loci were enriched for function in cell proliferation, and in glioblastoma (i.e., GBM) cells with DGCR8 or DROSHA knockdown, lnc-pri-miRNA screen hits still regulated cell growth. To molecularly dissect the function of a lnc-pri-miRNA locus, we studied LOC646329 (also known as MIR29HG), which hosts the miR-29a/b1 cluster. In GBM cells, LOC646329 knockdown reduced miR-29a/b1 levels, and these cells exhibited decreased growth. However, genetic deletion of the miR-29a/b1 cluster (LOC646329-miR29Δ) did not decrease cell growth, while knockdown of LOC646329-miR29Δ transcripts reduced cell proliferation. The miR-29a/b1-independent activity of LOC646329 corresponded to enhancer-like activation of a neighboring oncogene (MKLN1), regulating cell propagation. The LOC646329 locus interacts with the MKLN1 promoter, and antisense oligonucleotide knockdown of the lncRNA disrupts these interactions and reduces the enhancer-like activity. More broadly, analysis of genome-wide data from multiple human cell types showed that lnc-pri-miRNA loci are significantly enriched for DNA looping interactions with gene promoters as well as genomic and epigenetic characteristics of transcriptional enhancers. Functional studies of additional lnc-pri-miRNA loci demonstrated cognate miRNA-independent enhancer-like activity. Together, these data demonstrate that lnc-pri-miRNA loci can regulate cell biology via both miRNA-dependent and miRNA-independent mechanisms.

Ensemble learning for classifying single-cell data and projection across reference atlases.

SUMMARY: Single-cell data are being generated at an accelerating pace. How best to project data across single-cell atlases is an open problem. We developed a boosted learner that overcomes the greatest challenge with status quo classifiers: low sensitivity, especially when dealing with rare cell types. By comparing novel and published data from distinct scRNA-seq modalities that were acquired from the same tissues, we show that this approach preserves cell-type labels when mapping across diverse platforms. AVAILABILITY AND IMPLEMENTATION: https://github.com/diazlab/ELSA. CONTACT: aaron.diaz@ucsf.edu. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

CONICS integrates scRNA-seq with DNA sequencing to map gene expression to tumor sub-clones.

Motivation: Single-cell RNA-sequencing (scRNA-seq) has enabled studies of tissue composition at unprecedented resolution. However, the application of scRNA-seq to clinical cancer samples has been limited, partly due to a lack of scRNA-seq algorithms that integrate genomic mutation data. Results: To address this, we present. CONICS: COpy-Number analysis In single-Cell RNA-Sequencing. CONICS is a software tool for mapping gene expression from scRNA-seq to tumor clones and phylogenies, with routines enabling: the quantitation of copy-number alterations in scRNA-seq, robust separation of neoplastic cells from tumor-infiltrating stroma, inter-clone differential-expression analysis and intra-clone co-expression analysis. Availability and implementation: CONICS is written in Python and R, and is available from https://github.com/diazlab/CONICS. Supplementary information: Supplementary data are available at Bioinformatics online.

Recurrent epimutations activate gene body promoters in primary glioblastoma.

Aberrant DNA hypomethylation may play an important role in the growth rate of glioblastoma (GBM), but the functional impact on transcription remains poorly understood. We assayed the GBM methylome with MeDIP-seq and MRE-seq, adjusting for copy number differences, in a small set of non-glioma CpG island methylator phenotype (non-G-CIMP) primary tumors. Recurrent hypomethylated loci were enriched within a region of chromosome 5p15 that is specified as a cancer amplicon and also encompasses TERT, encoding telomerase reverse transcriptase, which plays a critical role in tumorigenesis. Overall, 76 gene body promoters were recurrently hypomethylated, including TERT and the oncogenes GLI3 and TP73. Recurring hypomethylation also affected previously unannotated alternative promoters, and luciferase reporter assays for three of four of these promoters confirmed strong promoter activity in GBM cells. Histone H3 lysine 4 trimethylation (H3K4me3) ChIP-seq on tissue from the GBMs uncovered peaks that coincide precisely with tumor-specific decrease of DNA methylation at 200 loci, 133 of which are in gene bodies. Detailed investigation of TP73 and TERT gene body hypomethylation demonstrated increased expression of corresponding alternate transcripts, which in TP73 encodes a truncated p73 protein with oncogenic function and in TERT encodes a putative reverse transcriptase-null protein. Our findings suggest that recurring gene body promoter hypomethylation events, along with histone H3K4 trimethylation, alter the transcriptional landscape of GBM through the activation of a limited number of normally silenced promoters within gene bodies, in at least one case leading to expression of an oncogenic protein.

SCell: integrated analysis of single-cell RNA-seq data.

UNLABELLED: Analysis of the composition of heterogeneous tissue has been greatly enabled by recent developments in single-cell transcriptomics. We present SCell, an integrated software tool for quality filtering, normalization, feature selection, iterative dimensionality reduction, clustering and the estimation of gene-expression gradients from large ensembles of single-cell RNA-seq datasets. SCell is open source, and implemented with an intuitive graphical interface. Scripts and protocols for the high-throughput pre-processing of large ensembles of single-cell, RNA-seq datasets are provided as an additional resource. AVAILABILITY AND IMPLEMENTATION: Binary executables for Windows, MacOS and Linux are available at http://sourceforge.net/projects/scell, source code and pre-processing scripts are available from https://github.com/diazlab/SCellSupplementary information: Supplementary data are available at Bioinformatics online. CONTACT: aaron.diaz@ucsf.edu.

Single-cell sequencing maps gene expression to mutational phylogenies in PDGF- and EGF-driven gliomas.

Glioblastoma multiforme (GBM) is the most common and aggressive type of primary brain tumor. Epidermal growth factor (EGF) and platelet-derived growth factor (PDGF) receptors are frequently amplified and/or possess gain-of-function mutations in GBM However, clinical trials of tyrosine-kinase inhibitors have shown disappointing efficacy, in part due to intra-tumor heterogeneity. To assess the effect of clonal heterogeneity on gene expression, we derived an approach to map single-cell expression profiles to sequentially acquired mutations identified from exome sequencing. Using 288 single cells, we constructed high-resolution phylogenies of EGF-driven and PDGF-driven GBMs, modeling transcriptional kinetics during tumor evolution. Descending the phylogenetic tree of a PDGF-driven tumor corresponded to a progressive induction of an oligodendrocyte progenitor-like cell type, expressing pro-angiogenic factors. In contrast, phylogenetic analysis of an EGFR-amplified tumor showed an up-regulation of pro-invasive genes. An in-frame deletion in a specific dimerization domain of PDGF receptor correlates with an up-regulation of growth pathways in a proneural GBM and enhances proliferation when ectopically expressed in glioma cell lines. In-frame deletions in this domain are frequent in public GBM data.

HiTSelect: a comprehensive tool for high-complexity-pooled screen analysis.

Genetic screens of an unprecedented scale have recently been made possible by the availability of high-complexity libraries of synthetic oligonucleotides designed to mediate either gene knockdown or gene knockout, coupled with next-generation sequencing. However, several sources of random noise and statistical biases complicate the interpretation of the resulting high-throughput data. We developed HiTSelect, a comprehensive analysis pipeline for rigorously selecting screen hits and identifying functionally relevant genes and pathways by addressing off-target effects, controlling for variance in both gene silencing efficiency and sequencing depth of coverage and integrating relevant metadata. We document the superior performance of HiTSelect using data from both genome-wide RNAi and CRISPR/Cas9 screens. HiTSelect is implemented as an open-source package, with a user-friendly interface for data visualization and pathway exploration. Binary executables are available at http://sourceforge.net/projects/hitselect/, and the source code is available at https://github.com/diazlab/HiTSelect.

Polycomb-like 3 promotes polycomb repressive complex 2 binding to CpG islands and embryonic stem cell self-renewal.

Polycomb repressive complex 2 (PRC2) trimethylates lysine 27 of histone H3 (H3K27me3) to regulate gene expression during diverse biological transitions in development, embryonic stem cell (ESC) differentiation, and cancer. Here, we show that Polycomb-like 3 (Pcl3) is a component of PRC2 that promotes ESC self-renewal. Using mass spectrometry, we identified Pcl3 as a Suz12 binding partner and confirmed Pcl3 interactions with core PRC2 components by co-immunoprecipitation. Knockdown of Pcl3 in ESCs increases spontaneous differentiation, yet does not affect early differentiation decisions as assessed in teratomas and embryoid bodies, indicating that Pcl3 has a specific role in regulating ESC self-renewal. Consistent with Pcl3 promoting PRC2 function, decreasing Pcl3 levels reduces H3K27me3 levels while overexpressing Pcl3 increases H3K27me3 levels. Furthermore, chromatin immunoprecipitation and sequencing (ChIP-seq) reveal that Pcl3 co-localizes with PRC2 core component, Suz12, and depletion of Pcl3 decreases Suz12 binding at over 60% of PRC2 targets. Mutation of conserved residues within the Pcl3 Tudor domain, a domain implicated in recognizing methylated histones, compromises H3K27me3 formation, suggesting that the Tudor domain of Pcl3 is essential for function. We also show that Pcl3 and its paralog, Pcl2, exist in different PRC2 complexes but bind many of the same PRC2 targets, particularly CpG islands regulated by Pcl3. Thus, Pcl3 is a component of PRC2 critical for ESC self-renewal, histone methylation, and recruitment of PRC2 to a subset of its genomic sites.

Challenges in the discovery of tumor-specific alternative splicing-derived cell-surface antigens in glioma.

Despite advancements in cancer immunotherapy, solid tumors remain formidable challenges. In glioma, profound inter- and intra-tumoral heterogeneity of antigen landscape hampers therapeutic development. Therefore, it is critical to consider alternative sources to expand the repertoire of targetable (neo-)antigens and improve therapeutic outcomes. Accumulating evidence suggests that tumor-specific alternative splicing (AS) could be an untapped reservoir of antigens. In this study, we investigated tumor-specific AS events in glioma, focusing on those predicted to generate major histocompatibility complex (MHC)-presentation-independent, cell-surface antigens that could be targeted by antibodies and chimeric antigen receptor-T cells. We systematically analyzed bulk RNA-sequencing datasets comparing 429 tumor samples (from The Cancer Genome Atlas) and 9166 normal tissue samples (from the Genotype-Tissue Expression project), and identified 13 AS events in 7 genes predicted to be expressed in more than 10% of the patients, including PTPRZ1 and BCAN, which were corroborated by an external RNA-sequencing dataset. Subsequently, we validated our predictions and elucidated the complexity of the isoforms using full-length transcript amplicon sequencing on patient-derived glioblastoma cells. However, analyses of the RNA-sequencing datasets of spatially mapped and longitudinally collected clinical tumor samples unveiled remarkable spatiotemporal heterogeneity of the candidate AS events. Furthermore, proteomics analysis did not reveal any peptide spectra matching the putative antigens. Our investigation illustrated the diverse characteristics of the tumor-specific AS events and the challenges of antigen exploration due to their notable spatiotemporal heterogeneity and elusive nature at the protein levels. Redirecting future efforts toward intracellular, MHC-presented antigens could offer a more viable avenue.

VC-resist glioblastoma cell state: vessel co-option as a key driver of chemoradiation resistance.

Glioblastoma (GBM) is a highly lethal type of cancer. GBM recurrence following chemoradiation is typically attributed to the regrowth of invasive and resistant cells. Therefore, there is a pressing need to gain a deeper understanding of the mechanisms underlying GBM resistance to chemoradiation and its ability to infiltrate. Using a combination of transcriptomic, proteomic, and phosphoproteomic analyses, longitudinal imaging, organotypic cultures, functional assays, animal studies, and clinical data analyses, we demonstrate that chemoradiation and brain vasculature induce cell transition to a functional state named VC-Resist (vessel co-opting and resistant cell state). This cell state is midway along the transcriptomic axis between proneural and mesenchymal GBM cells and is closer to the AC/MES1-like state. VC-Resist GBM cells are highly vessel co-opting, allowing significant infiltration into the surrounding brain tissue and homing to the perivascular niche, which in turn induces even more VC-Resist transition. The molecular and functional characteristics of this FGFR1-YAP1-dependent GBM cell state, including resistance to DNA damage, enrichment in the G2M phase, and induction of senescence/stemness pathways, contribute to its enhanced resistance to chemoradiation. These findings demonstrate how vessel co-option, perivascular niche, and GBM cell plasticity jointly drive resistance to therapy during GBM recurrence.

Normalization, bias correction, and peak calling for ChIP-seq.

Next-generation sequencing is rapidly transforming our ability to profile the transcriptional, genetic, and epigenetic states of a cell. In particular, sequencing DNA from the immunoprecipitation of protein-DNA complexes (ChIP-seq) and methylated DNA (MeDIP-seq) can reveal the locations of protein binding sites and epigenetic modifications. These approaches contain numerous biases which may significantly influence the interpretation of the resulting data. Rigorous computational methods for detecting and removing such biases are still lacking. Also, multi-sample normalization still remains an important open problem. This theoretical paper systematically characterizes the biases and properties of ChIP-seq data by comparing 62 separate publicly available datasets, using rigorous statistical models and signal processing techniques. Statistical methods for separating ChIP-seq signal from background noise, as well as correcting enrichment test statistics for sequence-dependent and sonication biases, are presented. Our method effectively separates reads into signal and background components prior to normalization, improving the signal-to-noise ratio. Moreover, most peak callers currently use a generic null model which suffers from low specificity at the sensitivity level requisite for detecting subtle, but true, ChIP enrichment. The proposed method of determining a cell type-specific null model, which accounts for cell type-specific biases, is shown to be capable of achieving a lower false discovery rate at a given significance threshold than current methods.

PP2Ac Deficiency Enhances Tumor Immunogenicity by Activating STING-Type I Interferon Signaling in Glioblastoma.

Glioblastoma (GBM) is an immunologically "cold" tumor that does not respond to current immunotherapy. Here, we demonstrate a fundamental role for the α-isoform of the catalytic subunit of protein phosphatase-2A (PP2Ac) in regulating glioma immunogenicity. Genetic ablation of PP2Ac in glioma cells enhanced double-stranded DNA (dsDNA) production and cGAS-type I IFN signaling, MHC-I expression, and tumor mutational burden. In coculture experiments, PP2Ac deficiency in glioma cells promoted dendritic cell (DC) cross-presentation and clonal expansion of CD8+ T cells. In vivo, PP2Ac depletion sensitized tumors to immune-checkpoint blockade and radiotherapy treatment. Single-cell analysis demonstrated that PP2Ac deficiency increased CD8+ T-cell, natural killer cell, and DC accumulation and reduced immunosuppressive tumor-associated macrophages. Furthermore, loss of PP2Ac increased IFN signaling in myeloid and tumor cells and reduced expression of a tumor gene signature associated with worse patient survival in The Cancer Genome Atlas. Collectively, this study establishes a novel role for PP2Ac in inhibiting dsDNA-cGAS-STING signaling to suppress antitumor immunity in glioma. SIGNIFICANCE: PP2Ac deficiency promotes cGAS-STING signaling in glioma to induce a tumor-suppressive immune microenvironment, highlighting PP2Ac as a potential therapeutic target to enhance tumor immunogenicity and improve response to immunotherapy.

Glioblastoma induces the recruitment and differentiation of hybrid neutrophils from skull bone marrow.

Tumor-associated neutrophil (TAN) effects on glioblastoma biology remain under-characterized. We show here that 'hybrid' neutrophils with dendritic features - including morphological complexity, expression of antigen presentation genes, and the ability to process exogenous peptide and stimulate MHCII-dependent T cell activation - accumulate intratumorally and suppress tumor growth in vivo . Trajectory analysis of patient TAN scRNA-seq identifies this phenotype as a polarization state which is distinct from canonical cytotoxic TANs and differentiates intratumorally from immature precursors absent in circulation. Rather, these hybrid-inducible immature neutrophils - which we identified in patient and murine glioblastomas - arise from local skull marrow. Through labeled skull flap transplantation and targeted ablation, we characterize calvarial marrow as a potent contributor of antitumoral myeloid APCs, including hybrid TANs and dendritic cells, which elicit T cell cytotoxicity and memory. As such, agents augmenting neutrophil egress from skull marrow - such as intracalvarial AMD3100 whose survival prolonging-effect in GBM we demonstrate - present therapeutic potential.

Single-cell RNA sequencing and spatial transcriptomics reveal cancer-associated fibroblasts in glioblastoma with protumoral effects.

Cancer-associated fibroblasts (CAFs) were presumed absent in glioblastoma given the lack of brain fibroblasts. Serial trypsinization of glioblastoma specimens yielded cells with CAF morphology and single-cell transcriptomic profiles based on their lack of copy number variations (CNVs) and elevated individual cell CAF probability scores derived from the expression of 9 CAF markers and absence of 5 markers from non-CAF stromal cells sharing features with CAFs. Cells without CNVs and with high CAF probability scores were identified in single-cell RNA-Seq of 12 patient glioblastomas. Pseudotime reconstruction revealed that immature CAFs evolved into subtypes, with mature CAFs expressing actin alpha 2, smooth muscle (ACTA2). Spatial transcriptomics from 16 patient glioblastomas confirmed CAF proximity to mesenchymal glioblastoma stem cells (GSCs), endothelial cells, and M2 macrophages. CAFs were chemotactically attracted to GSCs, and CAFs enriched GSCs. We created a resource of inferred crosstalk by mapping expression of receptors to their cognate ligands, identifying PDGF and TGF-ß as mediators of GSC effects on CAFs and osteopontin and HGF as mediators of CAF-induced GSC enrichment. CAFs induced M2 macrophage polarization by producing the extra domain A (EDA) fibronectin variant that binds macrophage TLR4. Supplementing GSC-derived xenografts with CAFs enhanced in vivo tumor growth. These findings are among the first to identify glioblastoma CAFs and their GSC interactions, making them an intriguing target.

Challenges in the discovery of tumor-specific alternative splicing-derived cell-surface antigens in glioma.

Background: Despite advancements in cancer immunotherapy, solid tumors remain formidable challenges. In glioma, profound inter-and intra-tumoral heterogeneity of antigen landscape hampers therapeutic development. Therefore, it is critical to consider alternative sources to expand the repertoire of targetable (neo-)antigens and improve therapeutic outcomes. Accumulating evidence suggests that tumor-specific alternative splicing (AS) could be an untapped reservoir of neoantigens. Results: In this study, we investigated tumor-specific AS events in glioma, focusing on those predicted to generate major histocompatibility complex (MHC)-presentation-independent, cell-surface neoantigens that could be targeted by antibodies and chimeric antigen receptor (CAR)-T cells. We systematically analyzed bulk RNA-sequencing datasets comparing 429 tumor samples (from The Cancer Genome Atlas [TCGA]) and 9,166 normal tissue samples (from the Genotype-Tissue Expression project [GTEx]), and identified 13 AS events in 7 genes predicted to be expressed in more than 10% of the patients, including PTPRZ1 and BCAN , which were corroborated by an external RNA-sequencing dataset. Subsequently, we validated our predictions and elucidated the complexity of the isoforms using full-length transcript amplicon sequencing on patient-derived glioblastoma cells. However, analyses of the RNA-sequencing datasets of spatially mapped and longitudinally collected clinical tumor samples unveiled remarkable spatiotemporal heterogeneity of the candidate AS events. Furthermore, proteomics analysis did not reveal any peptide spectra matching the putative neoantigens. Conclusions: Our investigation illustrated the diverse characteristics of the tumor-specific AS events and the challenges of antigen exploration due to their notable spatiotemporal heterogeneity and elusive nature at the protein levels. Redirecting future efforts toward intracellular, MHC-presented antigens could offer a more viable avenue.

Tumor-wide RNA splicing aberrations generate immunogenic public neoantigens.

T-cell-mediated immunotherapies are limited by the extent to which cancer-specific antigens are homogenously expressed throughout a tumor. We reasoned that recurrent splicing aberrations in cancer represent a potential source of tumor-wide and public neoantigens, and to test this possibility, we developed a novel pipeline for identifying neojunctions expressed uniformly within a tumor across diverse cancer types. Our analyses revealed multiple neojunctions that recur across patients and either exhibited intratumor heterogeneity or, in some cases, were tumor-wide. We identified CD8+ T-cell clones specific for neoantigens derived from tumor-wide and conserved neojunctions in GNAS and RPL22 , respectively. TCR-engineered CD8 + T-cells targeting these mutations conferred neoantigen-specific tumor cell eradication. Furthermore, we revealed that cancer-specific dysregulation in splicing factor expression leads to recurrent neojunction expression. Together, these data reveal that a subset of neojunctions are both intratumorally conserved and public, providing the molecular basis for novel T-cell-based immunotherapies that address intratumoral heterogeneity.

A single-cell atlas of glioblastoma evolution under therapy reveals cell-intrinsic and cell-extrinsic therapeutic targets.

Recent longitudinal studies of glioblastoma (GBM) have demonstrated a lack of apparent selection pressure for specific DNA mutations in recurrent disease. Single-cell lineage tracing has shown that GBM cells possess a high degree of plasticity. Together this suggests that phenotype switching, as opposed to genetic evolution, may be the escape mechanism that explains the failure of precision therapies to date. We profiled 86 primary-recurrent patient-matched paired GBM specimens with single-nucleus RNA, single-cell open-chromatin, DNA and spatial transcriptomic/proteomic assays. We found that recurrent GBMs are characterized by a shift to a mesenchymal phenotype. We show that the mesenchymal state is mediated by activator protein 1. Increased T-cell abundance at recurrence was prognostic and correlated with hypermutation status. We identified tumor-supportive networks of paracrine and autocrine signals between GBM cells, nonmalignant neuroglia and immune cells. We present cell-intrinsic and cell-extrinsic targets and a single-cell multiomics atlas of GBM under therapy.