Search

Patients with glioblastoma (GBM) need bold new approaches to their treatment, yet progress has been hindered by a relative inability to dynamically track treatment response, mechanisms of resistance, evolution of targetable mutations, and changes in mutational burden. We are writing on behalf of a multidisciplinary group of academic neuro-oncology professionals who met at the collaborative Christopher Davidson Forum at Washington University in St Louis in the fall of 2019. We propose a dramatic but necessary change to the routine management of patients with GBM to advance the field: to routinely biopsy recurrent GBM at the time of presumed recurrence. Data derived from these samples will identify true recurrence vs treatment effect, avoid treatments with little chance of success, enable clinical trial access, and aid in the scientific advancement of our understanding of GBM.

Human tumors are composed of diverse malignant and nonmalignant cells, generating a complex ecosystem that governs tumor biology and response to treatments. Recent technological advances have enabled the characterization of tumors at single-cell resolution, providing a compelling strategy to dissect their intricate biology. Here we describe recent developments in single-cell expression profiling and the studies applying them in clinical settings. We highlight some of the powerful insights gleaned from these studies for tumor classification, stem cell programs, tumor microenvironment, metastasis, and response to targeted and immune therapies.Significance:Intratumor heterogeneity (ITH) has been a major barrier to our understanding of cancer. Single-cell genomics is leading a revolution in our ability to systematically dissect ITH. In this review, we focus on single-cell expression profiling and lessons learned in key aspects of human tumor biology.

Glioma-associated oncogene homolog 3 (GLI3), whose main function is to inhibit GLI1, has been associated with neuronal differentiation in medulloblastoma. However, it is not clear what molecular subtype(s) show increased GLI3 expression. GLI3 levels were assessed by immunohistochemistry in 2 independent cohorts, including a total of 88 cases, and found to be high in both WNT- and SHH-activated medulloblastoma. Analysis of bulk mRNA expression data and single cell RNA sequencing studies confirmed that GLI1 and GLI3 are highly expressed in SHH-activated medulloblastoma, whereas GLI3 but not GLI1 is highly expressed in WNT-activated medulloblastoma. Immunohistochemical analysis has shown that GLI3 is expressed inside the neuronal differentiated nodules of SHH-activated medulloblastoma, whereas GLI1/2 are expressed in desmoplastic areas. In contrast, GLI3 is diffusely expressed in WNT-activated medulloblastoma, whereas GLI1 is suppressed. Our data suggest that GLI3 may be a master regulator of neuronal differentiation and morphology in these subgroups.

Ewing’s Sarcoma Ewing’s sarcoma, an aggressive cancer of bone and soft tissue, primarily affects children and young adults. A t(11;22) translocation is noted in 85 to 90% of cases. Management inclu...

SummaryRecent data showed promising signs of objective tumor responses in subsets of patients with low grade glioma treated with inhibitors of mutant IDH (IDHi) (1). However, the molecular and cellular underpinnings of such responses are not known. Here, we profiled 6,039 transcriptomes by single-cell or single-nucleus RNA-sequencing isolated from three IDH-mutant oligodendroglioma patients with clinical response to IDHi. Importantly, the tissues were sampled on-drug, four weeks from treatment initiation and our dataset includes a matched pre- and on-treatment sample pair. We integrate our findings with analysis of 8,241 transcriptomes from seven untreated samples, 134 bulk samples from the TCGA and experimental models.(2,3) We find that IDHi treatment induces a robust differentiation towards glial lineages, accompanied by a depletion of stem-like cells and a reduction of cell proliferation. Our study provides the first evidence in patients of the differentiating potential of IDHi on the cellular hierarchies that drive oligodendrogliomas.

Tumor cell extensions called tumor microtubes (TMs) in glioma resemble neurites during neurodevelopment and connect glioma cells to a network that has considerable relevance for tumor progression and therapy resistance. The determination of interconnectivity in individual tumors has been challenging and the impact of tumor cell connectivity on patient survival remained unresolved so far. Here, a connectivity signature from single-cell RNA-sequenced (scRNA-Seq) xenografted primary glioblastoma (GB) cells was established and clinically validated. Thirty-four of 40 connectivity genes were related to neurogenesis, neural tube development or glioma progression, including the TM-network-relevant GAP43 gene. Astrocytic-like and mesenchymal-like GB cells had the highest connectivity signature scores in scRNA-Seq data of patient-derived xenografts and patient samples. In 230 human GBs, high connectivity correlated with the mesenchymal expression subtype, TP53 wildtype, and with dismal patient survival. CHI3L1 was identified as a robust molecular marker of connectivity. Thus, the connectivity signature allows novel insights into brain tumor biology, provides a proof-of-principle that tumor cell connectivity is relevant for patients’ prognosis, and serves as a robust biomarker that can be used for future clinical trials.Statement of significanceIntegration of GB cells into functional networks drives tumor progression and resistance. Here, we established and validated a novel connectivity gene expression signature of single GB cells and whole tumors that can be easily applied to clinical and preclinical samples. It is shown that connectivity is determining prognosis combining molecular, functional and clinical insights into the disease.

Single-cell RNA sequencing has revealed extensive transcriptional cell state diversity in cancer, often observed independently of genetic heterogeneity, raising the central question of how malignant cell states are encoded epigenetically. To address this, here we performed multiomics single-cell profiling-integrating DNA methylation, transcriptome and genotype within the same cells-of diffuse gliomas, tumors characterized by defined transcriptional cell state diversity. Direct comparison of the epigenetic profiles of distinct cell states revealed key switches for state transitions recapitulating neurodevelopmental trajectories and highlighted dysregulated epigenetic mechanisms underlying gliomagenesis. We further developed a quantitative framework to directly measure cell state heritability and transition dynamics based on high-resolution lineage trees in human samples. We demonstrated heritability of malignant cell states, with key differences in hierarchal and plastic cell state architectures in IDH-mutant glioma versus IDH-wild-type glioblastoma, respectively. This work provides a framework anchoring transcriptional cancer cell states in their epigenetic encoding, inheritance and transition dynamics.

High-grade gliomas are lethal brain cancers whose progression is robustly regulated by neuronal activity. Activity-regulated release of growth factors promotes glioma growth, but this alone is insufficient to explain the effect that neuronal activity exerts on glioma progression. Here we show that neuron and glioma interactions include electrochemical communication through bona fide AMPA receptor-dependent neuron–glioma synapses. Neuronal activity also evokes non-synaptic activity-dependent potassium currents that are amplified by gap junction-mediated tumour interconnections, forming an electrically coupled network. Depolarization of glioma membranes assessed by in vivo optogenetics promotes proliferation, whereas pharmacologically or genetically blocking electrochemical signalling inhibits the growth of glioma xenografts and extends mouse survival. Emphasizing the positive feedback mechanisms by which gliomas increase neuronal excitability and thus activity-regulated glioma growth, human intraoperative electrocorticography demonstrates increased cortical excitability in the glioma-infiltrated brain. Together, these findings indicate that synaptic and electrical integration into neural circuits promotes glioma progression. Neurons form synapses onto glioma cells, and depolarization of glioma membranes promotes glioma growth in vivo, whereas blocking electrochemical signalling blocks tumour growth.

Pilocytic astrocytoma (PA), the most common childhood brain tumor, is a low-grade glioma with a single driver BRAF rearrangement. Here, we perform scRNAseq in six PAs using methods that enabled detection of the rearrangement. When compared to higher-grade gliomas, a strikingly higher proportion of the PA cancer cells exhibit a differentiated, astrocyte-like phenotype. A smaller proportion of cells exhibit a progenitor-like phenotype with evidence of proliferation. These express a mitogen-activated protein kinase (MAPK) programme that was absent from higher-grade gliomas. Immune cells, especially microglia, comprise 40% of all cells in the PAs and account for differences in bulk expression profiles between tumor locations and subtypes. These data indicate that MAPK signaling is restricted to relatively undifferentiated cancer cells in PA, with implications for investigational therapies directed at this pathway., Pilocytic astrocytoma is a low-grade pediatric glioma, characterized by a single BRAF rearrangement. Here, Reitman and colleagues use single-cell RNA sequencing to reveal molecular hallmarks of the disease that might be targeted therapeutically.

Bulk genomic analyses and expression profiling of clinical specimens have shaped much of our understanding of cancer in patients. However, human tumors are intricate ecosystems composed of diverse cells, including malignant, immune, and stromal subsets, whose precise characterization is masked by bulk genomic methods. Single-cell genomic techniques have emerged as powerful approaches to dissect human tumors at the resolution of individual cells, providing a compelling approach to deciphering cancer biology. Here, we discuss some of the common themes emerging from initial studies of single-cell RNA sequencing in cancer and then highlight challenges in cancer biology for which emerging single-cell genomics methods may provide a compelling approach.

Human tumors are complex ecosystems where diverse cancer and noncancer cells interact to determine tumor biology and response to therapies. Genomic and transcriptomic methods have traditionally profiled these intricate ecosystems as bulk samples, thereby masking individual cellular programs and the variability among them. Recent advances in single-cell profiling have paved the way for studying tumors at the resolution of individual cells, providing a compelling strategy to bridge gaps in our understanding of human tumors. Here, we review methodologies for single-cell expression profiling of tumors and the initial studies deploying them in clinical contexts. We highlight how these studies uncover new biology and provide insights into drug resistance, stem cell programs, metastasis, and tumor classifications. We also discuss areas of technology development in single-cell genomics that provide new tools to address key questions in cancer biology. These emerging studies and technologies have the potential to revolutionize our understanding and management of human malignancies.

The transcriptomic classification of glioblastoma (GBM) has failed to predict survival and therapeutic vulnerabilities. A computational approach for unbiased identification of core biological traits of single cells and bulk tumors uncovered four tumor cell states and GBM subtypes distributed along neurodevelopmental and metabolic axes, classified as proliferative/progenitor, neuronal, mitochondrial and glycolytic/plurimetabolic. Each subtype was enriched with biologically coherent multiomic features. Mitochondrial GBM was associated with the most favorable clinical outcome. It relied exclusively on oxidative phosphorylation for energy production, whereas the glycolytic/plurimetabolic subtype was sustained by aerobic glycolysis and amino acid and lipid metabolism. Deletion of the glucose-proton symporter SLC45A1 was the truncal alteration most significantly associated with mitochondrial GBM, and the reintroduction of SLC45A1 in mitochondrial glioma cells induced acidification and loss of fitness. Mitochondrial, but not glycolytic/plurimetabolic, GBM exhibited marked vulnerability to inhibitors of oxidative phosphorylation. The pathway-based classification of GBM informs survival and enables precision targeting of cancer metabolism.

Communication between cancer cells and immune cells is a key determinant of the glioblastoma ecosystem and its response to therapies, but remains poorly understood. Here we leveraged single-cell RNA-sequencing (scRNA-seq) of human samples and mouse models, deconvolution analysis of bulk specimen from The Cancer Genome Atlas (TCGA) and functional approaches, to dissect cellular cross-talks in glioblastoma. We demonstrate that macrophages induce a transition of glioblastoma cells into mesenchymal-like (MES-like) states. This effect is mediated, both in vitro and in vivo, by macrophage-derived Oncostatin M (OSM) and its cognate receptor OSMR on glioblastoma cells. We show that MES-like glioblastoma states are also associated with increased expression of a mesenchymal program in macrophages and with increased cytotoxicity of T cells, highlighting extensive alterations of the immune microenvironment with potential therapeutic implications.

Targeting of the most aggressive tumor cell subpopulations is key for effective management of most solid malignancies. However, the metastable nature of tumor heterogeneity, which allows cells to transition between strong and weak tumorigenic phenotypes, and the lack of reliable markers of tumor-promoting properties hamper identification of the most relevant cells. To overcome these obstacles, we designed a functional microRNA (miR)-based live-cell reporter assay to identify highly tumorigenic cells in xenotransplants of primary Ewing sarcoma (EwS) 3D cultures. Leveraging the inverse relationship between cell pluripotency and miR-145 expression, we successfully separated highly tumorigenic, metastasis-prone (miR-145 low ) cells from poorly tumorigenic, nonmetastatic (miR-145 high ) counterparts. Gene expression and functional studies of the two cell populations identified the EPHB2 receptor as a prognostic biomarker in patients with EwS and a major promoter of metastasis. Our study provides a simple and powerful means to identify and isolate tumor cells that display aggressive behavior.

Histone mutations occur in approximately 4% of different cancer types. In 2012, mutations were found in the gene encoding histone variant H3.3 (H3F3A gene) in pediatric diffuse intrinsic pontine gliomas and pediatric hemispheric gliomas. Tumors with mutations in the H3F3A gene are generally characterized as histone mutated gliomas (HMGs) or diffuse midline gliomas. HMGs are a rare subtype of glial tumor that is malignant and fast growing, carrying a poor prognosis. In 2017, the Beau Biden Cancer Moonshot Program appropriated $1.7 billion toward cancer care in 10 select areas. The National Cancer Institute (NCI) was granted support to focus specifically on rare central nervous system (CNS) tumors through NCI-CONNECT. Its mission is to address the challenges and unmet needs in CNS cancer research and treatment by connecting patients, providers, researchers, and advocacy organizations to work in partnership. On September 27, 2018, NCI-CONNECT convened a workshop on histone mutated midline glioma, one of the 12 CNS cancers included in its initial portfolio. Three leaders in the field provided an overview of advances in histone mutated midline glioma research. These experts shared observations and experiences related to common scientific and clinical challenges in studying these tumors. Although the clinical focus of this workshop was on adult patients, one important objective was to start a collaborative dialogue between pediatric and adult clinicians and researchers. Meeting participants identified needs for diagnostic and treatment standards, disease biology and biological targets for this cancer, disease-specific trial designs, and developed a list of action items and future direction.

Ewing sarcoma (EwS) is associated with poor prognosis despite current multimodal therapy. Targeting of EWS-FLI1, the fusion protein responsible for its pathogenesis, and its principal downstream targets has not yet produced satisfactory therapeutic options, fueling the search for alternative approaches. Here, we show that the oncofetal RNA-binding protein LIN28B regulates the stability of EWS-FLI1 mRNA in ~10% of EwSs. LIN28B depletion in these tumors leads to a decrease in the expression of EWS-FLI1 and its direct transcriptional network, abrogating EwS cell self-renewal and tumorigenicity. Moreover, pharmacological inhibition of LIN28B mimics the effect of LIN28B depletion, suggesting that LIN28B sustains the emergence of a subset of EwS in which it also serves as an effective therapeutic target.

Background Single-cell RNA sequencing is the reference technique for characterizing the heterogeneity of the tumor microenvironment. The composition of the various cell types making up the microenvironment can significantly affect the way in which the immune system activates cancer rejection mechanisms. Understanding the cross-talk signals between immune cells and cancer cells is of fundamental importance for the identification of immuno-oncology therapeutic targets. Results We present a novel method, single-cell Tumor–Host Interaction tool (scTHI), to identify significantly activated ligand–receptor interactions across clusters of cells from single-cell RNA sequencing data. We apply our approach to uncover the ligand–receptor interactions in glioma using 6 publicly available human glioma datasets encompassing 57,060 gene expression profiles from 71 patients. By leveraging this large-scale collection we show that unexpected cross-talk partners are highly conserved across different datasets in the majority of the tumor samples. This suggests that shared cross-talk mechanisms exist in glioma. Conclusions Our results provide a complete map of the active tumor–host interaction pairs in glioma that can be therapeutically exploited to reduce the immunosuppressive action of the microenvironment in brain tumor.

Glioblastoma is a devastating form of brain cancer. To identify aspects of tumor heterogeneity that may illuminate drivers of tumor invasion, we created a glioblastoma tumor cell atlas with single-cell transcriptomics of cancer cells mapped onto a reference framework of the developing and adult human brain. We find that multiple GSC subtypes exist within a single tumor. Within these GSCs, we identify an invasive cell population similar to outer radial glia (oRG), a fetal cell type that expands the stem cell niche in normal human cortex. Using live time-lapse imaging of primary resected tumors, we discover that tumor derived oRG-like cells undergo characteristic mitotic somal translocation behavior previously only observed in human development, suggesting a reactivation of developmental programs. In addition, we show that PTPRZ1 mediates both mitotic somal translocation and glioblastoma tumor invasion. These data suggest that the presence of heterogeneous GSCs may underlie glioblastoma’s rapid progression and invasion.

Diffuse intrinsic pontine glioma (DIPG) is a universally fatal malignancy of the childhood central nervous system, with a median overall survival of 9–11 months. We have previously shown that primary DIPG tissue contains numerous tumor-associated macrophages, and substantial work has demonstrated a significant pathological role for adult glioma-associated macrophages. However, work over the past decade has highlighted many molecular and genomic differences between pediatric and adult high-grade gliomas. Thus, we directly compared inflammatory characteristics of DIPG and adult glioblastoma (GBM). We found that the leukocyte (CD45+) compartment in primary DIPG tissue samples is predominantly composed of CD11b + macrophages, with very few CD3+ T-lymphocytes. In contrast, T-lymphocytes are more abundant in adult GBM tissue samples. RNA sequencing of macrophages isolated from primary tumor samples revealed that DIPG- and adult GBM-associated macrophages both express gene programs related to ECM remodeling and angiogenesis, but DIPG-associated macrophages express substantially fewer inflammatory factors than their adult GBM counterparts. Examining the secretome of glioma cells, we found that patient-derived DIPG cell cultures secrete markedly fewer cytokines and chemokines than patient-derived adult GBM cultures. Concordantly, bulk and single-cell RNA sequencing data indicates low to absent expression of chemokines and cytokines in DIPG. Together, these observations suggest that the inflammatory milieu of the DIPG tumor microenvironment is fundamentally different than adult GBM. The low intrinsic inflammatory signature of DIPG cells may contribute to the lack of lymphocytes and non-inflammatory phenotype of DIPG-associated microglia/macrophages. Understanding the glioma subtype-specific inflammatory milieu may inform the design and application of immunotherapy-based treatments. Electronic supplementary material The online version of this article (10.1186/s40478-018-0553-x) contains supplementary material, which is available to authorized users.

Communication between cancer cells and immune cells is a key determinant of the glioblastoma ecosystem and its response to therapies, but remains poorly understood. Here we leveraged single-cell RNA-sequencing (scRNA-seq) of human samples and mouse models, deconvolution analysis of bulk specimens from The Cancer Genome Atlas (TCGA) and functional approaches to dissect cellular states and cross-talk in glioblastoma. We demonstrate that macrophages induce a transition of glioblastoma cells into mesenchymal-like (MES-like) states. This effect is mediated, both in vitro and in vivo, by macrophage-derived Oncostatin M (OSM) and its cognate receptor OSMR on glioblastoma cells. We show that MES-like glioblastoma states are associated with increased T cells cytotoxicity and potentially with better clinical response to immunotherapies. Overall, our work dissects the cellular interactions within the glioblastoma microenvironment, with potential implications for therapies.

BackgroundRecurrent high-grade glioma (HGG) represents a significant clinical unmet need with expected survival between 6 to 9 months. Oncolytic viruses are a new therapeutic approach for solid tumors that deploy oncolytic activity combined with local and systemic immune activation. CAN-3110 (rQNestin34.5v2) is an oncolytic herpes simplex virus (HSV), modified to encode the HSV1 ICP34.5 protein under control of the nestin promoter. Selective expression of nestin in brain tumors confers tumor-restricted replication of CAN-3110. We conducted an open-label dose-escalation phase 1 clinical trial in patients with recurrent HGG to evaluate safety, tolerability, and immunological changes after CAN-3110 treatment.MethodsThirty patients with biopsy-confirmed recurrent HGG were enrolled from September 2017 to February 2020. CAN-3110 was injected intratumorally starting at 1x106 plaque forming units (pfu) and dose-escalated by half log to 1x1010 pfu. Patients also received standard of care. Peripheral blood mononuclear cells (PBMCs), plasma and tumor samples were collected for analysis at different time-points post treatment. We evaluated HSV antigen expression in tumor tissue. RNA sequencing and T cell receptor (TCR) rearrangement analysis was performed in matched tissue and PBMCs. Cytokine profiling was completed in 29 patients at baseline, day 2, and day 28 post treatment.ResultsEighteen patients were recruited at their first recurrence and 12 at the second recurrence. Three patients presented with multifocal disease. Tumor volume ranged from 357.4 to and 54,036.1mm3 (median 7,733.9mm3, SDV 15,610.2). CAN-3110 was well-tolerated with no dose-limiting toxicity. Median overall survival was 11.7 months. We demonstrated persistence of HSV antigen and CD8+ T cell infiltrates at the site of injected tumor. Preliminary analysis revealed expansion of shared TCR clonotypes and upregulation of pro-inflammatory genes in post-treatment tumors and peripheral blood samples. Longitudinal modeling of cytokine profiling demonstrated increased levels of IL-6, VEGF alpha, CCL2 and IL1-RA and a decrease in GCP-2 levels at day 2 post-treatment (p ConclusionsIntratumoral administration of CAN-3110 appears well-tolerated in recurrent HGG. Histologic, molecular, and cytokine analyses demonstrate persistence of viral antigen as well as local and systemic immune activation after treatment.Ethics ApprovalThe study was approved by the Office for Human Research Studies at Dana-Farber Cancer Institute, Protocol Number 16–557.

Glioblastomas are malignant neoplasms composed of diverse cell populations. This intratumoral diversity has an underlying architecture, with a hierarchical relationship through clonal evolution from a common ancestor. Therapies are limited by emergence of resistant subclones from this phylogenetic reservoir. To characterize this clonal ancestral origin of recurrent tumors, we determined phylogenetic relationships using whole exome sequencing of pre-treatment IDH1/2 wild-type glioblastoma specimens, matched to post-treatment autopsy samples (n = 9) and metastatic extracranial post-treatment autopsy samples (n = 3). We identified “truncal” genetic events common to the evolutionary ancestry of the initial specimen and later recurrences, thereby inferring the identity of the precursor cell population. Mutations were identified in a subset of cases in known glioblastoma genes such as NF1(n = 3), TP53(n = 4) and EGFR(n = 5). However, by phylogenetic analysis, there were no protein-coding mutations as recurrent truncal events across the majority of cases. In contrast, whole copy-loss of chromosome 10 (12 of 12 cases), copy-loss of chromosome 9p21 (11 of 12 cases) and copy-gain in chromosome 7 (10 of 12 cases) were identified as shared events in the majority of cases. Strikingly, mutations in the TERT promoter were also identified as shared events in all evaluated pairs (9 of 9). Thus, we define four truncal non-coding genomic alterations that represent early genomic events in gliomagenesis, that identify the persistent cellular reservoir from which glioblastoma recurrences emerge. Therapies to target these key early genomic events are needed. These findings offer an evolutionary explanation for why precision therapies that target protein-coding mutations lack efficacy in GBM., Brain cancer: Non-coding genomic changes fuel glioblastoma growth Non-coding and structural alterations may be early drivers of brain cancer development. A team led by Priscilla Brastianos and Tracy Batchelor from Massachusetts General Hospital, Boston, USA, analyzed the genetic landscape of glioblastoma by comparing pre-treatment and autopsy tumor specimens from 12 patients who died of the aggressive brain cancer. They identified a common set of four genetic events that occurred early in the evolution of nearly every patient’s cancer: three losses or gains of chromosome regions or entire chromosomes, and mutations in the gene-activating promoter of TERT, which encodes an enzyme implicated in the cancer’s growth. The findings help explain why therapies that target protein-coding mutations don’t work in brain cancer when they do in other tumor types. They also point to new drug targets.

Single-cell RNA sequencing is the reference technique to characterize the heterogeneity of tumor microenvironment and can be efficiently used to discover cross-talk mechanisms between immune cells and cancer cells. We present a novel method, single cell Tumor-Host Interaction tool (scTHI), to identify significantly activated ligand-receptor interactions across clusters of cells from single-cell RNA sequencing data. We apply our approach to uncover the ligand-receptor interactions in glioma using six publicly available human glioma datasets encompassing 71 patients. We provide a comprehensive map of the signalling mechanisms between malignant cells and non-malignant cells in glioma uncovering potential novel therapeutic targets.

Medulloblastoma is a malignant childhood cerebellar tumor comprised of distinct molecular subgroups. Whereas genomic characteristics of these subgroups are well defined, the extent to which cellular diversity underlies their divergent biology and clinical behaviour remains largely unexplored. We used single-cell transcriptomics to investigate intra- and inter-tumoral heterogeneity in twenty-five medulloblastomas spanning all molecular subgroups. WNT, SHH, and Group 3 tumors comprised subgroup-specific undifferentiated and differentiated neuronal-like malignant populations, whereas Group 4 tumors were exclusively comprised of differentiated neuronal-like neoplastic cells. SHH tumors closely resembled granule neurons of varying differentiation states that correlated with patient age. Group 3 and Group 4 tumors exhibited a developmental trajectory from primitive progenitor-like to more mature neuronal-like cells, whose relative proportions distinguished these subgroups. Cross-species transcriptomics defined distinct glutamatergic populations as putative cells-of-origin for SHH and Group 4 subtypes. Collectively, these data provide novel insights into the cellular and developmental states underlying subtype-specific medulloblastoma biology.

Synovial sarcoma (SyS) is an aggressive neoplasm driven by the SS18–SSX fusion, and is characterized by low T cell infiltration. Here, we studied the cancer–immune interplay in SyS using an integrative approach that combines single-cell RNA sequencing (scRNA-seq), spatial profiling and genetic and pharmacological perturbations. scRNA-seq of 16,872 cells from 12 human SyS tumors uncovered a malignant subpopulation that marks immune-deprived niches in situ and is predictive of poor clinical outcomes in two independent cohorts. Functional analyses revealed that this malignant cell state is controlled by the SS18–SSX fusion, is repressed by cytokines secreted by macrophages and T cells, and can be synergistically targeted with a combination of HDAC and CDK4/CDK6 inhibitors. This drug combination enhanced malignant-cell immunogenicity in SyS models, leading to induced T cell reactivity and T cell–mediated killing. Our study provides a blueprint for investigating heterogeneity in fusion-driven malignancies and demonstrates an interplay between immune evasion and oncogenic processes that can be co-targeted in SyS and potentially in other malignancies. Single-cell transcriptional profiling of primary human synovial sarcoma tumors suggests that combinatorial treatment with HDAC and CDK4/CDK6 inhibitors could enhance tumor immunogenicity.

Single cell genomics is essential to chart the complex tumor ecosystem. While single cell RNA-Seq (scRNA-Seq) profiles RNA from cells dissociated from fresh tumor tissues, single nucleus RNA-Seq (snRNA-Seq) is needed to profile frozen or hard-to-dissociate tumors. Each strategy requires modifications to fit the unique characteristics of different tissue and tumor types, posing a barrier to adoption. Here, we developed a systematic toolbox for profiling fresh and frozen clinical tumor samples using scRNA-Seq and snRNA-Seq, respectively. We tested eight tumor types of varying tissue and sample characteristics (resection, biopsy, ascites, and orthotopic patient-derived xenograft): lung cancer, metastatic breast cancer, ovarian cancer, melanoma, neuroblastoma, pediatric sarcoma, glioblastoma, pediatric high-grade glioma, and chronic lymphocytic leukemia. Analyzing 212,498 cells and nuclei from 39 clinical samples, we evaluated protocols by cell quality, recovery rate, and cellular composition. We optimized protocols for fresh tissue dissociation for different tumor types using a decision tree to account for the technical and biological variation between clinical samples. We established methods for nucleus isolation from OCT embedded and fresh-frozen tissues, with an optimization matrix varying mechanical force, buffer, and detergent. scRNA-Seq and snRNA-Seq from matched samples recovered the same cell types and intrinsic expression profiles, but at different proportions. Our work provides direct guidance across a broad range of tumors, including criteria for testing and selecting methods from the toolbox for other tumors, thus paving the way for charting tumor atlases.

Summary 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.

T and B cells are the two known lineages of adaptive immune cells. Here, we describe a previously unknown lymphocyte that is a dual expresser (DE) of TCR and BCR and key lineage markers of both B and T cells. In type 1 diabetes (T1D), DEs are predominated by one clonotype that encodes a potent CD4 T cell autoantigen in its antigen binding site. Molecular dynamics simulations revealed that this peptide has an optimal binding register for diabetogenic HLA-DQ8. In concordance, a synthetic version of the peptide forms stable DQ8 complexes and potently stimulates autoreactive CD4 T cells from T1D patients, but not healthy controls. Moreover, mAbs bearing this clonotype are autoreactive against CD4 T cells and inhibit insulin tetramer binding to CD4 T cells. Thus, compartmentalization of adaptive immune cells into T and B cells is not absolute, and violators of this paradigm are likely key drivers of autoimmune diseases.

The identity and degree of heterogeneity of glial progenitors and their contributions to brain tumor malignancy remain elusive. By applying lineage-targeted single-cell transcriptomics, we uncover an unanticipated diversity of glial progenitor pools with unique molecular identities in developing brain. Our analysis identifies distinct transitional intermediate states and their divergent developmental trajectories in astroglial and oligodendroglial lineages. Moreover, intersectional analysis uncovers analogous intermediate progenitors during brain tumorigenesis, wherein oligodendrocyte-progenitor intermediates are abundant, hyper-proliferative, and progressively reprogrammed toward a stem-like state susceptible to further malignant transformation. Similar actively cycling intermediate progenitors are prominent components in human gliomas with distinct driver mutations. We further unveil lineage-driving networks underlying glial fate specification and identify Zfp36l1 as necessary for oligodendrocyte-astrocyte lineage transition and glioma growth. Together, our results resolve the dynamic repertoire of common and divergent glial progenitors during development and tumorigenesis and highlight Zfp36l1 as a molecular nexus for balancing glial cell-fate decision and controlling gliomagenesis.

There are ~1,400 oligodendroglioma diagnoses in the U.S. annually, with up to 15,000 persons living with this disease today. However, there have not been any new treatments approved for oligodendroglioma in the past 25 years. A lack of funding and advocacy has contributed to this, but the lack of reagents allowing researchers to study the disease may be the most critical factor. Oligodendroglioma tissue and disease models, including cell lines and xenograft models, are crucial resources needed for researchers to understand oligodendroglioma biology, identify new treatment strategies, and perform preclinical research. Oligo Nation has created a collection of oligo tissue large enough to enable adequately-sized studies and meet the needs of researchers interested in studying Oligodendroglioma. Oligo Nation is also investing in developing new models of the disease as well as acquiring existing models to make them easily accessible to the research community. A number of resources are now available to support academic researchers: An oligodendroglioma biobank with over 100 patient tissue samples; WGS and RNA sequencing data, including long reads for these samples; An oligodendroglioma cell line repository with 3 lines available (to date) to researchers; An oligodendroglioma PDX repository will be available in the coming year through the NCI Patient Derived Model Repository; GEMM models that include both IDH1 (R132) and CIC mutations, two of the most common mutations in oligodendroglioma, will be available by the end of 2020. Extensive analysis including high-throughput and CRISPR knockout screens will be performed on all existing and new cell lines. The data generated on these models will be stored in a public repository and made available to the research community. Making these reagents available to the research community will hopefully increase research in oligodendroglioma, and will set us on a path to discover improved and ultimately curative treatments for this disease.

Immune checkpoint inhibitors (ICI) have revolutionized oncologic treatment for metastatic melanoma. With improved systemic control, there has been increasing prevalence of patients with brain metastases. Recent evidence has demonstrated intracranial responses in a subset of these patients treated with ICI. We hypothesize that the response to ICI in melanoma brain metastases (MBM) is reflective of unique features within the tumor microenvironment of the brain. A cohort of 27 patients, encompassing 8 pre- and 19 post-immunotherapy MBM underwent single cell RNA sequencing (Smart-Seq2). The cohort includes patients with longitudinal cranial resections and simultaneously resected, spatially distinct tumors. Each tumor underwent unsupervised transcriptomic analysis, differential gene expression, inferred copy number variation, and T-cell receptor (TCR) clonotyping. Published extracranial melanoma single cell datasets were used to compare the tumor microenvironment of the brain and periphery in response to ICI. A total of 14,027 cells (6,189 malignant, 7,838 non-malignant) were sequenced. Brain metastases demonstrated a heterogeneous distribution of macrophage states. Intracranial macrophages were found to be more tumor-supportive than their extracranial counterparts. MBM also included a distribution of reactive neutrophils and astrocytes. Analysis across pre- and post-treatment MBM demonstrated an increase in clonally expanded T cells in patients responding to ICI. Across longitudinal brain metastases collected from the same patients, there was evidence of identical T cell clones across timepoints and locations. Single cell sequencing of MBM provides insights into the cellular composition of the tumor and microenvironment. Our data suggest the cellular heterogeneity within MBM is unique when compared to extracranial disease. Additionally, T cell clonal expansion is found following ICI and T cells of the same clonotype infiltrate spatially and temporally separated brain metastases. These findings raise potential therapeutic implications as we learn to target the differential features of the innate and adaptive immune system within brain metastases and their extracranial counterparts.

Glioblastoma is an incurable malignancy that remains poorly understood. By large-scale single cell RNA-seq and follow up experiments we found that malignant cells in glioblastoma exist in four main cellular states that recapitulate distinct steps of neural development, are influenced by the tumor microenvironment, and exhibit plasticity. Importantly, we show that while each glioblastoma contains similar cellular states, their relative frequency varies between tumors and is influenced by chromosomal aberrations as well as by the tumor microenvironment. In particular, we demonstrate that macrophages induce a transition of glioblastoma cells into mesenchymal-like (MES-like) states. This effect is mediated, both in vitro and in vivo, by macrophage-derived ligands that interact with receptors on glioblastoma cells. Finally, we show that MES-like glioblastoma states is associated with altered T cells cytotoxicity and response to immunotherapies. Overall, this work dissects the cellular diversity and interactions within the glioblastoma microenvironment, with potential implications for therapies. Citation Format: Toshi Hara, Rony Chanoch, Mario Suva, Itay Tirosh. Dissecting glioblastoma by single cell RNA-seq [abstract]. In: Proceedings of the AACR Virtual Special Conference on Tumor Heterogeneity: From Single Cells to Clinical Impact; 2020 Sep 17-18. Philadelphia (PA): AACR; Cancer Res 2020;80(21 Suppl):Abstract nr IA24.

Human diffuse gliomas are incurable brain tumors, where cellular state diversity fuels tumor progression and resistance to therapy. Single-cell RNA-sequencing (scRNAseq) studies recently charted the cellular states of the two major categories of human gliomas, IDH-mutant gliomas (IDH-MUT) and IDH-wildtype glioblastoma (GBM), showing that malignant cells partly recapitulate neurodevelopmental trajectories. This raises the central questions of how cell states are encoded epigenetically and whether unidirectional hierarchies or more plastic state transitions govern glioma cellular architectures. To address these questions, we generated multi-omics single-cell profiling, integrating DNA methylation (DNAme), transcriptome and genotyping of 1,728 cells from 11 GBM and IDH-MUT primary patient samples. The assessment of DNAme intra-tumoral heterogeneity of malignant cells revealed that single-cell DNAme profiles within tumors span multiple bulk subtypes, are associated with important biological features of malignant cells, and may be confounded by the tumor micro-environment. Such sources of intra-tumoral heterogeneity in bulk profiles are important to recognize, as DNAme profiling is increasingly being utilized for bulk clinical brain tumor classification. The direct comparison of the methylomes of distinct glioma cellular states revealed Polycomb repressive complex 2 (PRC2) targets DNAme as a key switch in the differentiation of malignant GBM cells. In contrast, dissecting aberrant circuits of hypermethylation and gene expression in IDH-MUT gliomas, we observed a decoupling of the promoter methylation-expression relationship, with disruption of CTCF insulation and enhancer vulnerability which increases with cellular differentiation. To define cell state transition dynamics directly in patient samples, we generated high-resolution lineage histories of glioma cells using heritable DNAme changes, and projected the scRNAseq-derived cell states onto the lineage trees. This analysis demonstrated that cell states are heritable across malignant gliomas and, while in IDH-MUT differentiation far outpaces de-differentiation, GBM harbors a higher degree of cellular state plasticity allowing reversion of GBM cells from a differentiated to a stem-like state. Overall, our work provides detailed insights into gliomagenesis, dissecting the epigenetic encoding, regulatory programs, and dynamics of the cellular states that drive human gliomas. Importantly, it also carries significant translational implication, as the high degree of de-differentiation in GBM challenges the paradigm of therapeutically targeting glioma stem-like cells to deprive tumors of their ability to regenerate. Citation Format: Federico Gaiti, Ronan Chaligne, Dana Silverbush, Joshua S. Schiffman, Hannah R. Weisman, Lloyd Kluegel, Simon Gritsch, Sunil D. Deochand, Alyssa R. Richman, Johanna Klughammer, Tommaso Biancalani, Christoph Muus, Caroline Sheridan, Alicia Alonso, Franco Izzo, Orit Rozenblatt-Rosen, Aviv Regev, Mario L. Suvà, Dan A. Landau. Deciphering differentiation hierarchies, heritability and plasticity in human gliomas via single-cell multi-omics [abstract]. In: Proceedings of the AACR Virtual Special Conference on Tumor Heterogeneity: From Single Cells to Clinical Impact; 2020 Sep 17-18. Philadelphia (PA): AACR; Cancer Res 2020;80(21 Suppl):Abstract nr PO-019.

In situ transgenesis methods such as virus and electroporation can create somatic transgenic mice quickly, but they lack the exquisite control over copy number, zygosity, and locus specificity. We have recently established mosaic analysis by dual recombinase-mediated cassette exchange (MADR), which permits stable labeling of mutant cells expressing transgenic elements from precisely-defined chromosomal loci. MADR provides a toolkit of elements for combinatorial labeling, inducible/reversible transgene manipulation, VCre recombinase expression, and genetic manipulation of human cells. Further, we have demonstrated the versatility of MADR by creating glioma models with mixed, reporter-identified zygosity or with “personalized” driver mutations from pediatric glioma. For example, introducing H3f3a mutation variants with MADR regulates the spatiotemporal profile of glioma, and single-cell RNA and ATAC sequencing analysis demonstrates a recapitulation of developmental hierarchy seen in K27M-mutant human glioma. Moreover, we have generated novel models of supratentorial ependymoma using patient-derived oncofusion transgenes. These models display a high degree of fidelity and we now compare these models on a single-cell level with our previous models and human tumor cell transcriptomes. In addition, we now demonstrate the ability to generalize the MADR technology to other non-CNS tissues using local plasmid delivery. Finally, we have engineered human cells to allow for MADR transgenesis and somatic transgenic organoids. These combined approaches will enable researchers to discovery disease mechanisms and test therapeutics in more physiologically relevant cancer models.MADR is extensible to thousands of existing mouse lines and can be adapted to human cells, providing a flexible platform to democratize the generation of somatic transgenic disease models. Citation Format: Gi Bum Kim, David Rincon Fernandez Pacheco, David Saxon, Amy Yang, Sara Sabat, Marina Dutra-Clarke, Rachelle Levy, Ashley Watkins, Hannah Park, Aslam Abbasi Akhtar, Paul W. Linesch, Naomi Kobritz, Swasty S. Shandra, Katie Grausam, Alberto Ayala-Sarmiento, Jessica Molina, Kristyna Sedivakova, Daniel S. Gareau, Mariella G. Filbin, Serguei Bannykh, Jie Tang, Mario L. Suva, Bin Chen, Moise Danielpour, Joshua J. Breunig. MADR: Rapid generation of somatic mosaics to model cancer in mice and human organoids [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr NG01.

Summary Ewing sarcoma (EwS) is associated with poor prognosis despite current multimodal therapy. Targeting of EWS-FLI1, the fusion protein responsible for its pathogenesis, and its principal downstream targets has not yet produced satisfactory therapeutic options, fueling the search for alternative approaches. Here, we show that the oncofetal RNA-binding protein LIN28B regulates the stability of EWS-FLI1 mRNA in ~10% of EwSs. LIN28B depletion in these tumors leads to a decrease in the expression of EWS-FLI1 and its direct transcriptional network, abrogating EwS cell self-renewal and tumorigenicity. Moreover, pharmacological inhibition of LIN28B mimics the effect of LIN28B depletion, suggesting that LIN28B sustains the emergence of a subset of EwS in which it also serves as an effective therapeutic target.

Summary Medulloblastoma is a malignant childhood cerebellar tumour comprised of distinct molecular subgroups. Whereas genomic characteristics of these subgroups are well defined, the extent to which cellular diversity underlies their divergent biology and clinical behaviour remains largely unexplored. We used single-cell transcriptomics to investigate intra- and inter-tumoural heterogeneity in twenty-five medulloblastomas spanning all molecular subgroups. WNT, SHH, and Group 3 tumours comprised subgroup-specific undifferentiated and differentiated neuronal-like malignant populations, whereas Group 4 tumours were exclusively comprised of differentiated neuronal-like neoplastic cells. SHH tumours closely resembled granule neurons of varying differentiation states that correlated with patient age. Group 3 and Group 4 tumours exhibited a developmental trajectory from primitive progenitor-like to more mature neuronal-like cells, whose relative proportions distinguished these subgroups. Cross-species transcriptomics defined distinct glutamatergic populations as putative cells-of-origin for SHH and Group 4 subtypes. Collectively, these data provide novel insights into the cellular and developmental states underlying subtype-specific medulloblastoma biology.

Diffuse gliomas are the most common human primary brain tumors and remain incurable. They are complex entities in which diverse genetic and nongenetic effects determine tumor biology and clinical course. Our current understanding of gliomas in patients is primarily based on genomic and transcriptomic methods that have profiled them as bulk, providing critical information yet masking the diversity of cells within each tumor. Recent advances in single-cell DNA and RNA profiling have paved the way to studying tumors at cellular resolution. Here, we review initial studies deploying single-cell analysis in clinical glioma samples, with a focus on RNA expression profiling. We highlight how these studies provide new insights into glioma biology, tumor heterogeneity, cancer cell lineages, cancer stem cell programs, the tumor microenvironment, and glioma classification.

Brain tumor-initiating cells (BTICs) have been identified as key contributors to therapy resistance, recurrence, and progression of diffuse gliomas, particularly glioblastoma (GBM). BTICs are elusive therapeutic targets that reside across the blood-brain barrier, underscoring the urgent need to develop novel therapeutic strategies. Additionally, intratumoral heterogeneity and adaptations to therapeutic pressure by BTICs impede the discovery of effective anti-BTIC therapies and limit the efficacy of individual gene targeting. Recent discoveries in the genetic and epigenetic determinants of BTIC tumorigenesis offer novel opportunities for RNAi-mediated targeting of BTICs. Here we show that BTIC growth arrest in vitro and in vivo is accomplished via concurrent siRNA knockdown of four transcription factors (SOX2, OLIG2, SALL2, and POU3F2) that drive the proneural BTIC phenotype delivered by multiplexed siRNA encapsulation in the lipopolymeric nanoparticle 7C1. Importantly, we demonstrate that 7C1 nano-encapsulation of multiplexed RNAi is a viable BTIC-targeting strategy when delivered directly in vivo in an established mouse brain tumor. Therapeutic potential was most evident via a convection-enhanced delivery method, which shows significant extension of median survival in two patient-derived BTIC xenograft mouse models of GBM. Our study suggests that there is potential advantage in multiplexed targeting strategies for BTICs and establishes a flexible nonviral gene therapy platform with the capacity to channel multiplexed RNAi schemes to address the challenges posed by tumor heterogeneity.

Tumor subclasses differ according to the genotypes and phenotypes of malignant cells as well as the composition of the tumor microenvironment (TME).We dissected these influences in isocitrate dehydrogenase (IDH)-mutant gliomas by combining 14,226 single-cell RNA sequencing (RNA-seq) profiles from 16 patient samples with bulk RNA-seq profiles from 165 patient samples. Differences in bulk profiles between IDH-mutant astrocytoma and oligodendroglioma can be primarily explained by distinct TME and signature genetic events, whereas both tumor types share similar developmental hierarchies and lineages of glial differentiation. As tumor grade increases, we find enhanced proliferation of malignant cells, larger pools of undifferentiated glioma cells, and an increase in macrophage over microglia expression programs in TME. Our work provides a unifying model for IDH-mutant gliomas and a general framework for dissecting the differences among human tumor subclasses., National Cancer Institute (U.S.) (Grant P30-CA14051)

Reversal of immune cell exhaustion through immune checkpoint blockade (ICB) has become a first-line approach for patients with metastatic melanoma due to success in controlling extracranial tumors. However, patients frequently experience discordant responses, with extracranial response and intracranial progression. We hypothesize that ICB exerts selective pressure leading to the clonal evolution of treatment resistant clones that ultimately culminate in disease progression. We collected a cohort of 97 patients, encompassing 312 pre- and post-immunotherapy melanoma tumors, for whole exome sequencing (WES) and included primary, extracranial, or intracranial samples. Each tumor was analyzed for somatic mutations, copy number alterations, neoantigen profile, and patient specific phylogenetic trees were constructed encompassing a tumor’s genetic subclones. Heterogeneity of the tumor microenvironment was evaluated using high multiplicity single-cell immunofluorescent staining (CycIF). Single cell sequencing was performed on fresh tissue from 4 pre-treatment and 14 post-immunotherapy melanoma brain metastases using the Smart-Seq2 protocol. WES of pre- and post-immunotherapy tumors yielded distinct patterns of clonal evolution and immunoediting within brain metastases compared to their extracranial counterparts, including mutations in B2M. In paired pre- and post-immunotherapy samples, CycIF demonstrated decreased in CD8 infiltration and increased CD45RO, FOXP3, and PD-L1 staining suggesting less cytotoxic, terminally differentiated T cells in resistant tumors. Single cell sequencing analysis of 3,974 tumor and immune cells demonstrated patient-specific tumor clustering and gene expression profiles mediating resistance to ICB. In conclusion, we document, for the first time, evidence of ongoing branched evolution during immunotherapy in brain metastases with divergence compared to systemic sites of disease. Next-generation sequencing provides novel insights into clonal evolution mediating discordant responses of intra- and extracranial sites and immunosuppressive features of the intracranial tumor microenvironment. Targeting these mechanisms of resistance provide potential therapeutic avenues for patients with progressive intracranial disease.

High-grade gliomas are a lethal group of cancers whose progression is robustly regulated by neuronal activity. Activity-regulated release of growth factors into the tumor microenvironment represents part of the mechanism by which neuronal activity influences glioma growth, but this alone is insufficient to explain the magnitude of the effect that activity exerts on glioma progression. Here, we report that neuron-glioma interactions include electrochemical communication through both bona fide synapses and activity-dependent potassium flux. Single cell transcriptomic analyses revealed unambiguous expression of synaptic genes by malignant glioma cells, and neuron to glioma synaptic structures were evident by electron microscopy. Whole cell patch clamp electrophysiology demonstrated AMPAR-mediated excitatory neurotransmission between pre-synaptic neurons and post-synaptic glioma cells. Millisecond timescale excitatory post-synaptic currents (EPSCs) were found in a subpopulation of glioma cells, reminiscent of the axon-glial synapses between neurons and normal oligodendrocyte precursor cells (OPCs). Neuronal activity also evokes a second, non-synaptic electrophysiological response characterized by a prolonged (>1 sec) depolarization in a subpopulation of glioma cells. These longer duration currents are blocked by tetrodotoxin or barium and induced by potassium, indicating neuronal activity-dependent potassium flux reminiscent of astrocyte currents. The amplitude of the prolonged currents is reduced by gap junction inhibitors, supporting the concept that gap junction-mediated tumor interconnections can function to amplify evoked potassium currents in an electrically coupled network. As membrane depolarization of normal neural precursor cells can regulate proliferation, differentiation and survival, and glioma cells exhibit two distinct mechanisms of neuronal activity-evoked membrane depolarization, we tested the hypothesis that membrane depolarization promotes glioma growth. Using in vivooptogenetic techniques to depolarize xenografted glioma cells, we found that glioma membrane depolarization robustly promoted proliferation, while pharmacologically or genetically blocking electrochemical signaling inhibited glioma xenograft growth and extended mouse survival. Together, these findings indicate that electrical circuit integration promotes glioma progression.

Single-cell RNA sequencing (scRNAseq) has been performed across a range of intermediate to high-grade gliomas that each harbor multiple driver mutations. Pilocytic astrocytoma (PA), the most common childhood brain tumor, is a low-grade glioma. PAs frequently harbor oncogenic KIAA1549-BRAF fusions, and exhibit low rates of other driver mutations. While PAs exhibit a favorable prognosis compared to the higher-grade gliomas, treatment morbidity and tumor recurrence can represent major challenges for some PA patients. We performed scRNAseq of both tumor and non-tumor cells in newly-diagnosed PAs that contained KIAA1549-BRAF rearrangements. To confidently distinguish tumor cells from nontumor cells, we sorted cells by glial progenitor marker A2B5 status and profiled KIAA1549-BRAF fusion status of cells using a sensitive quantitative PCR-based assay. Results were validated using RNA in situ hybridization. When compared to higher-grade gliomas, a higher proportion of the PA tumor cells exhibited a differentiated, astrocyte-like phenotype. A smaller proportion of cells exhibited a progenitor-like phenotype with evidence of proliferation. These progenitor-like tumor cells expressed a mitogen-activated protein kinase (MAPK) program that was absent from higher-grade gliomas. Similar patterns of expression of genes associated with the astrocyte-like and MAPK gene programs were also seen in formalin fixed, paraffin embedded PA tissues using RNA in situ hybridization. Immune cells, especially microglia, comprised almost half of all cells in the PAs and accounted for differences in bulk expression profiles between tumor locations and subtypes. These single cell transcriptional data reveal a transcriptional developmental hierarchy in a pediatric low grade glioma that is skewed towards more mature brain cells compared to higher-grade gliomas. The results indicate that future analyses of bulk PA tissues should attempt to account for considerable infiltration by nontumor cells. Finally, the finding that a MAPK gene program is not uniformly expressed in PA tumor cells has implications for ongoing clinical investigations of therapies directed at the MAPK pathway for the treatment of PA. Citation Format: Zachary J. Reitman, Brenton Paolella, Guillaume Bergthold, Kristine Pelton, Robert Jones, Sarah Becker, Claire E. Sinai, Haley Malkin, Ying Huang, Leslie Grimmett, Zachary T. Herbert, Yu Sun, Jessica Weatherbee, John Alberta, John Daley, Orit Rozenblatt-Rosen, Rosalind Segal, Daphne Haas-Kogan, Mariella G. Filbin, Mario L. Suva, Aviv Regev, Charles Stiles, Mark W. Kieran, Liliana Goumnerova, Keith L. Ligon, Alex K. Shalek, Pratiti Bandopadhayay, Rameen Beroukhim. Single cell RNA sequencing reveals mitogenic and progenitor gene programs inBRAF-rearranged pilocytic astrocytomas [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 3647.

We have developed a simple and generalizable in vivo method for brain tumor modeling, mosaic analysis by dual recombinase-mediated cassette exchange (MADR). MADR allows for stable labeling of mutant cells expressing transgenic elements from a precisely-defined chromosomal locus. We have demonstrated the power and versatility of MADR by creating novel glioma models with mixed, reporter-defined zygosity, or with “personalized” H3. 3-containing driver mutation signatures from pediatric glioma--each manipulation altering the spatiotemporal profile of resulting tumors. Further we have generated ependymoma models by employing patient-derived fusion driver mutations. Notably, each model displays divergent spatiotemporal tumor expansion profiles and cellular phenotypes. Now, we use single-cell RNA-seq to compare these models to their cells of mutation and to human datasets to elucidate the fundamental transcriptional programs and markers within and across these tumor types. This investigation will assess and credential these models against their clinical counterparts by scrutinizing the resulting datasets and validating their clinical relevance as pre-clinical models for therapeutic discovery and testing.

Pediatric high-grade gliomas (pHGGs) are a lethal group of cancers whose progression is robustly regulated by neuronal activity. Activity-regulated release of growth factors into the tumor microenvironment represents part of the mechanism by which neuronal activity influences pHGG growth, but this alone is insufficient to explain the magnitude of the effect that activity exerts on glioma progression. Here, we report that neuron-glioma interactions include bona fide synaptic communication. Single cell transcriptomic analyses of primary pediatric and adult glioma samples reveal unambiguous expression of synaptic genes by malignant glioma cells. Whole cell patch clamp recordings from xenografted, pediatric patient-derived glioma cells revealed the existence of AMPAR-mediated excitatory neurotransmission between pre-synaptic glutamatergic neurons and post-synaptic glioma cells. Millisecond timescale excitatory post-synaptic currents (EPSCs) that are depolarizing were observed in a subpopulation of pHGG cells and are associated with activity-induced glioma cell calcium transients. These excitatory axon-glioma synapses are reminiscent of the axon-glial synapses formed between neurons and oligodendrocyte precursor cells. A second electrophysiological response characterized by a prolonged (>1 sec) depolarization in response to neuronal activity was also observed. These longer duration currents are blocked by gap junction inhibitors, supporting the concept that gap junction-mediated tumor interconnections, such as observed with tumor microtubes, can function as an electrically coupled network. As neurotransmitter-mediated depolarization of normal neural precursor cells can profoundly affect precursor cell proliferation, differentiation and survival, we tested the hypothesis that depolarizing currents in pediatric glioma cells promote tumor growth. Using in vivooptogenetic techniques to depolarize xenografted pHGG cells expressing channelrhodopsin-2 (ChR2), we found that glioma depolarization robustly promoted proliferation, while expression of a dominant-negative AMPAR subunit (GluA2) that blocks neuron-glioma signaling inhibited pHGG xenograft growth and extended mouse survival. These findings suggest that integration of pediatric glioma into neural circuits promotes tumor progression.

Systematic use of next-generation sequencing technologies has transformed the field of cancer genomics, leading to the identification of genetic alterations in an unanticipated number of genes and regulatory elements. Here, we review recent advances in brain tumour genomics and highlight how these findings improve classification and diagnosis of brain tumours.The studies discussed in this review have shed light on different areas of neuro-oncology. In-depth analysis of paediatric low-grade gliomas as well as paediatric glioblastomas has clarified our molecular understanding of these diseases, clearly distinguishing them from their adult counterparts. Unexpected novel mutations have been discovered in adult low-grade astrocytomas and in glioblastomas. Novel studies also highlighted candidate tumour suppressor genes located on the chromosome arms frequently deleted in oligodendrogliomas. Finally, we review recent discoveries in the molecular landscapes of medulloblastomas and meningiomas.These recent studies begin to provide an in-depth view of the molecular routes leading to brain tumour development. The findings will be critical for refining classification systems and improving clinical management.

Brain tumor initiating cells (BTICs), also known as cancer stem cells, hijack high-affinity glucose uptake active normally in neurons to maintain energy demands. Here we link metabolic dysregulation in human BTICs to a nexus between MYC and de novo purine synthesis, mediating glucose-sustained anabolic metabolism. Inhibiting purine synthesis abrogated BTIC growth, self-renewal and in vivo tumor formation by depleting intracellular pools of purine nucleotides, supporting purine synthesis as a potential therapeutic point of fragility. In contrast, differentiated glioma cells were unaffected by the targeting of purine biosynthetic enzymes, suggesting selective dependence of BTICs. MYC coordinated the control of purine synthetic enzymes, supporting its role in metabolic reprogramming. Elevated expression of purine synthetic enzymes correlated with poor prognosis in glioblastoma patients. Collectively, our results suggest that stem-like glioma cells reprogram their metabolism to self-renew and fuel the tumor hierarchy, revealing potential BTIC cancer dependencies amenable to targeted therapy.

Although human tumours are shaped by the genetic evolution of cancer cells, evidence also suggests that they display hierarchies related to developmental pathways and epigenetic programs in which cancer stem cells (CSCs) can drive tumour growth and give rise to differentiated progeny1. Yet, unbiased evidence for CSCs in solid human malignancies remains elusive. Here we profile 4,347 single cells from six IDH1 or IDH2 mutant human oligodendrogliomas by RNA sequencing (RNA-seq) and reconstruct their developmental programs from genome-wide expression signatures. We infer that most cancer cells are differentiated along two specialized glial programs, whereas a rare subpopulation of cells is undifferentiated and associated with a neural stem cell expression program. Cells with expression signatures for proliferation are highly enriched in this rare subpopulation, consistent with a model in which CSCs are primarily responsible for fuelling the growth of oligodendroglioma in humans. Analysis of copy number variation (CNV) shows that distinct CNV sub-clones within tumours display similar cellular hierarchies, suggesting that the architecture of oligodendroglioma is primarily dictated by developmental programs. Subclonal point mutation analysis supports a similar model, although a full phylogenetic tree would be required to definitively determine the effect of genetic evolution on the inferred hierarchies. Our single-cell analyses provide insight into the cellular architecture of oligodendrogliomas at single-cell resolution and support the cancer stem cell model, with substantial implications for disease management.

Gliomas are the most common primary human brain tumors and occur in both adults and children. Over the past few years, systematic large-scale genomic and epigenomic profiling has provided unprecedented insight into their pathogenesis, uncovering alterations in an unanticipated number of genes and regulatory elements. In this review, we discuss recent discoveries about the genomics and epigenomics of adult and pediatric gliomas and highlight how some of the founding genetic mutations reshape the cancer epigenome. These studies provide an in-depth view of the molecular routes leading to glioma development, offer insight into the cancer stem cell model, help refine classifications, and should lay the foundation for improved clinical care.

BACKGROUND: The impact of individualized neoepitope vaccination targeting neoantigens arising from tumor-specific mutations for GBM, a low mutation burden tumor with an immunologically cold tumor microenvironment, as well as that of concurrently administered dexamethasone (dex), are unknown. METHODS: Individualized vaccination of up to 20 synthetic, long neoepitope peptides with high predicted HLA binding affinity admixed with poly-ICLC, were administered subcutaneously using a prime-boost schedule after RT to newly diagnosed, MGMT unmethylated, at least partially resected, GBM patients without progression after radiation (RT) in our phase 1b study. RESULTS: 9 of 10 screened patients had sufficient (10) identified neoepitope peptides. 8 patients without progression after RT received vaccine consisting of a median of 12 peptides (range, 7–20) beginning a median of 19.9 wks (range 17.1–24.7) after surgery. Adverse events were limited to infrequent grade 1/2 local reactions and fatigue. Median PFS and OS were 7.5 mths (90% CI: 6.2, 9.7) and 16.8 mths (90% CI: 9.6, 21.3). Evaluation of neoepitope-specific immune responses and tumor immune infiltrate analyses were performed on five patients with pre- and post-vaccination samples. Three patients on dex for post-RT edema during vaccine had no immune responses and no change in tumor infiltrating effector cells. In contrast 2 patients not on dex had robust, de novo immune responses against multiple predicted personal neoantigens including polyfunctional neoantigen-specific CD4+ and CD8+ T cell responses that were enriched for memory and activated phenotypes as well as increased numbers of tumor-infiltrating CD4+ and CD8+ T cells. T cell receptor analysis from one patient identified identical clonotypes isolated from post-vaccination tumor tissue and peripheral blood including two clonotypes specific for ARHGAP35, a neoantigen targeted by vaccination. CONCLUSIONS: Individualized, multi-neoepitope vaccines are feasible, safe and capable of generating systemic and intra-tumoral immune responses in GBM patients that appear to be abrogated by dex.