Spatial transcriptomics reveals segregation of tumor cell states in glioblastoma and marked immunosuppression within the perinecrotic niche.

Glioblastoma (GBM) remains an untreatable malignant tumor with poor patient outcomes, characterized by palisading necrosis and microvascular proliferation. While single-cell technology made it possible to characterize different lineage of glioma cells into neural progenitor-like (NPC-like), oligodendrocyte-progenitor-like (OPC-like), astrocyte-like (AC-like) and mesenchymal like (MES-like) states, it does not capture the spatial localization of these tumor cell states. Spatial transcriptomics empowers the study of the spatial organization of different cell types and tumor cell states and allows for the selection of regions of interest to investigate region-specific and cell-type-specific pathways. Here, we obtained paired 10x Chromium single-nuclei RNA-sequencing (snRNA-seq) and 10x Visium spatial transcriptomics data from three GBM patients to interrogate the GBM microenvironment. Integration of the snRNA-seq and spatial transcriptomics data reveals patterns of segregation of tumor cell states. For instance, OPC-like tumor and NPC-like tumor significantly segregate in two of the three samples. Our differentially expressed gene and pathway analyses uncovered significant pathways in functionally relevant niches. Specifically, perinecrotic regions were more immunosuppressive than the endogenous GBM microenvironment, and perivascular regions were more pro-inflammatory. Our gradient analysis suggests that OPC-like tumor cells tend to reside in areas closer to the tumor vasculature compared to tumor necrosis, which may reflect increased oxygen requirements for OPC-like cells. In summary, we characterized the localization of cell types and tumor cell states, the gene expression patterns, and pathways in different niches within the GBM microenvironment. Our results provide further evidence of the segregation of tumor cell states and highlight the immunosuppressive nature of the necrotic and perinecrotic niches in GBM.

An atlas of causal and mechanistic drivers of interpatient heterogeneity in glioma.

Grade IV glioma, formerly known as glioblastoma multiforme (GBM) is the most aggressive and lethal type of brain tumor, and its treatment remains challenging in part due to extensive interpatient heterogeneity in disease driving mechanisms and lack of prognostic and predictive biomarkers. Using mechanistic inference of node-edge relationship (MINER), we have analyzed multiomics profiles from 516 patients and constructed an atlas of causal and mechanistic drivers of interpatient heterogeneity in GBM (gbmMINER). The atlas has delineated how 30 driver mutations act in a combinatorial scheme to causally influence a network of regulators (306 transcription factors and 73 miRNAs) of 179 transcriptional "programs", influencing disease progression in patients across 23 disease states. Through extensive testing on independent patient cohorts, we share evidence that a machine learning model trained on activity profiles of programs within gbmMINER significantly augments risk stratification, identifying patients who are super-responders to standard of care and those that would benefit from 2 nd line treatments. In addition to providing mechanistic hypotheses regarding disease prognosis, the activity of programs containing targets of 2 nd line treatments accurately predicted efficacy of 28 drugs in killing glioma stem-like cells from 43 patients. Our findings demonstrate that interpatient heterogeneity manifests from differential activities of transcriptional programs, providing actionable strategies for mechanistically characterizing GBM from a systems perspective and developing better prognostic and predictive biomarkers for personalized medicine.

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

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

Gene regulatory network topology governs resistance and treatment escape in glioma stem-like cells.

Poor prognosis and drug resistance in glioblastoma (GBM) can result from cellular heterogeneity and treatment-induced shifts in phenotypic states of tumor cells, including dedifferentiation into glioma stem-like cells (GSCs). This rare tumorigenic cell subpopulation resists temozolomide, undergoes proneural-to-mesenchymal transition (PMT) to evade therapy, and drives recurrence. Through inference of transcriptional regulatory networks (TRNs) of patient-derived GSCs (PD-GSCs) at single-cell resolution, we demonstrate how the topology of transcription factor interaction networks drives distinct trajectories of cell state transitions in PD-GSCs resistant or susceptible to cytotoxic drug treatment. By experimentally testing predictions based on TRN simulations, we show that drug treatment drives surviving PD-GSCs along a trajectory of intermediate states, exposing vulnerability to potentiated killing by siRNA or a second drug targeting treatment-induced transcriptional programs governing non-genetic cell plasticity. Our findings demonstrate an approach to uncover TRN topology and use it to rationally predict combinatorial treatments that disrupts acquired resistance in GBM.

Surgical Management and Advances in the Treatment of Glioma.

The care of patients with both high-grade glioma and low-grade glioma necessitates an interdisciplinary collaboration between neurosurgeons, neuro-oncologists, neurologists and other practitioners. In this review, we aim to detail the considerations, approaches and advances in the neurosurgical care of gliomas. We describe the impact of extent-of-resection in high-grade and low-grade glioma, with particular focus on primary and recurrent glioblastoma. We address advances in surgical methods and adjunct technologies such as intraoperative imaging and fluorescence guided surgery that maximize extent-of-resection while minimizing the potential for iatrogenic neurological deficits. Finally, we review surgically-mediated therapies other than resection and discuss the role of neurosurgery in emerging paradigm-shifts in inter-disciplinary glioma management such as serial tissue sampling and "window of opportunity trials".

CD8+ T cells maintain killing of MHC-I-negative tumor cells through the NKG2D-NKG2DL axis.

The accepted paradigm for both cellular and anti-tumor immunity relies upon tumor cell killing by CD8+ T cells recognizing cognate antigens presented in the context of target cell major histocompatibility complex (MHC) class I (MHC-I) molecules. Likewise, a classically described mechanism of tumor immune escape is tumor MHC-I downregulation. Here, we report that CD8+ T cells maintain the capacity to kill tumor cells that are entirely devoid of MHC-I expression. This capacity proves to be dependent instead on interactions between T cell natural killer group 2D (NKG2D) and tumor NKG2D ligands (NKG2DLs), the latter of which are highly expressed on MHC-loss variants. Necessarily, tumor cell killing in these instances is antigen independent, although prior T cell antigen-specific activation is required and can be furnished by myeloid cells or even neighboring MHC-replete tumor cells. In this manner, adaptive priming can beget innate killing. These mechanisms are active in vivo in mice as well as in vitro in human tumor systems and are obviated by NKG2D knockout or blockade. These studies challenge the long-advanced notion that downregulation of MHC-I is a viable means of tumor immune escape and instead identify the NKG2D-NKG2DL axis as a therapeutic target for enhancing T cell-dependent anti-tumor immunity against MHC-loss variants.

WDR5 represents a therapeutically exploitable target for cancer stem cells in glioblastoma.

Glioblastomas (GBMs) are heterogeneous, treatment-resistant tumors driven by populations of cancer stem cells (CSCs). However, few molecular mechanisms critical for CSC population maintenance have been exploited for therapeutic development. We developed a spatially resolved loss-of-function screen in GBM patient-derived organoids to identify essential epigenetic regulators in the SOX2-enriched, therapy-resistant niche and identified WDR5 as indispensable for this population. WDR5 is a component of the WRAD complex, which promotes SET1 family-mediated Lys4 methylation of histone H3 (H3K4me), associated with positive regulation of transcription. In GBM CSCs, WDR5 inhibitors blocked WRAD complex assembly and reduced H3K4 trimethylation and expression of genes involved in CSC-relevant oncogenic pathways. H3K4me3 peaks lost with WDR5 inhibitor treatment occurred disproportionally on POU transcription factor motifs, including the POU5F1(OCT4)::SOX2 motif. Use of a SOX2/OCT4 reporter demonstrated that WDR5 inhibitor treatment diminished cells with high reporter activity. Furthermore, WDR5 inhibitor treatment and WDR5 knockdown altered the stem cell state, disrupting CSC in vitro growth and self-renewal, as well as in vivo tumor growth. These findings highlight the role of WDR5 and the WRAD complex in maintaining the CSC state and provide a rationale for therapeutic development of WDR5 inhibitors for GBM and other advanced cancers.

A single-cell based precision medicine approach using glioblastoma patient-specific models.

Glioblastoma (GBM) is a heterogeneous tumor made up of cell states that evolve over time. Here, we modeled tumor evolutionary trajectories during standard-of-care treatment using multi-omic single-cell analysis of a primary tumor sample, corresponding mouse xenografts subjected to standard of care therapy, and recurrent tumor at autopsy. We mined the multi-omic data with single-cell SYstems Genetics Network AnaLysis (scSYGNAL) to identify a network of 52 regulators that mediate treatment-induced shifts in xenograft tumor-cell states that were also reflected in recurrence. By integrating scSYGNAL-derived regulatory network information with transcription factor accessibility deviations derived from single-cell ATAC-seq data, we developed consensus networks that modulate cell state transitions across subpopulations of primary and recurrent tumor cells. Finally, by matching targeted therapies to active regulatory networks underlying tumor evolutionary trajectories, we provide a framework for applying single-cell-based precision medicine approaches to an individual patient in a concurrent, adjuvant, or recurrent setting.

A Systems Approach to Brain Tumor Treatment.

Brain tumors are among the most lethal tumors. Glioblastoma, the most frequent primary brain tumor in adults, has a median survival time of approximately 15 months after diagnosis or a five-year survival rate of 10%; the recurrence rate is nearly 90%. Unfortunately, this prognosis has not improved for several decades. The lack of progress in the treatment of brain tumors has been attributed to their high rate of primary therapy resistance. Challenges such as pronounced inter-patient variability, intratumoral heterogeneity, and drug delivery across the blood-brain barrier hinder progress. A comprehensive, multiscale understanding of the disease, from the molecular to the whole tumor level, is needed to address the intratumor heterogeneity resulting from the coexistence of a diversity of neoplastic and non-neoplastic cell types in the tumor tissue. By contrast, inter-patient variability must be addressed by subtyping brain tumors to stratify patients and identify the best-matched drug(s) and therapies for a particular patient or cohort of patients. Accomplishing these diverse tasks will require a new framework, one involving a systems perspective in assessing the immense complexity of brain tumors. This would in turn entail a shift in how clinical medicine interfaces with the rapidly advancing high-throughput (HTP) technologies that have enabled the omics-scale profiling of molecular features of brain tumors from the single-cell to the tissue level. However, several gaps must be closed before such a framework can fulfill the promise of precision and personalized medicine for brain tumors. Ultimately, the goal is to integrate seamlessly multiscale systems analyses of patient tumors and clinical medicine. Accomplishing this goal would facilitate the rational design of therapeutic strategies matched to the characteristics of patients and their tumors. Here, we discuss some of the technologies, methodologies, and computational tools that will facilitate the realization of this vision to practice.

Single-cell CUT&Tag analysis of chromatin modifications in differentiation and tumor progression.

Methods for quantifying gene expression1 and chromatin accessibility2 in single cells are well established, but single-cell analysis of chromatin regions with specific histone modifications has been technically challenging. In this study, we adapted the CUT&Tag method3 to scalable nanowell and droplet-based single-cell platforms to profile chromatin landscapes in single cells (scCUT&Tag) from complex tissues and during the differentiation of human embryonic stem cells. We focused on profiling polycomb group (PcG) silenced regions marked by histone H3 Lys27 trimethylation (H3K27me3) in single cells as an orthogonal approach to chromatin accessibility for identifying cell states. We show that scCUT&Tag profiling of H3K27me3 distinguishes cell types in human blood and allows the generation of cell-type-specific PcG landscapes from heterogeneous tissues. Furthermore, we used scCUT&Tag to profile H3K27me3 in a patient with a brain tumor before and after treatment, identifying cell types in the tumor microenvironment and heterogeneity in PcG activity in the primary sample and after treatment.

The Judicious Use of Stereotactic Radiosurgery and Hypofractionated Stereotactic Radiotherapy in the Management of Large Brain Metastases.

Brain metastases are the most common intracranial malignant tumor in adults and are a cause of significant morbidity and mortality for cancer patients. Large brain metastases, defined as tumors with a maximum dimension >2 cm, present a unique clinical challenge for the delivery of stereotactic radiosurgery (SRS) as patients often present with neurologic symptoms that require expeditious treatment that must also be balanced against the potential consequences of surgery and radiation therapy-namely, leptomeningeal disease (LMD) and radionecrosis (RN). Hypofractionated stereotactic radiotherapy (HSRT) and pre-operative SRS have emerged as novel treatment techniques to help improve local control rates and reduce rates of RN and LMD for this patient population commonly managed with post-operative SRS. Recent literature suggests that pre-operative SRS can potentially half the risk of LMD compared to post-operative SRS and that HSRT can improve risk of RN to less than 10% while improving local control when meeting the appropriate goals for biologically effective dose (BED) and dose-volume constraints. We recommend a 3- or 5-fraction regimen in lieu of SRS delivering 15 Gy or less for large metastases or resection cavities. We provide a table comparing the BED of commonly used SRS and HSRT regimens, and provide an algorithm to help guide the management of these challenging clinical scenarios.

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

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

Unravelling subclonal heterogeneity and aggressive disease states in TNBC through single-cell RNA-seq.

Triple-negative breast cancer (TNBC) is an aggressive subtype characterized by extensive intratumoral heterogeneity. To investigate the underlying biology, we conducted single-cell RNA-sequencing (scRNA-seq) of >1500 cells from six primary TNBC. Here, we show that intercellular heterogeneity of gene expression programs within each tumor is variable and largely correlates with clonality of inferred genomic copy number changes, suggesting that genotype drives the gene expression phenotype of individual subpopulations. Clustering of gene expression profiles identified distinct subgroups of malignant cells shared by multiple tumors, including a single subpopulation associated with multiple signatures of treatment resistance and metastasis, and characterized functionally by activation of glycosphingolipid metabolism and associated innate immunity pathways. A novel signature defining this subpopulation predicts long-term outcomes for TNBC patients in a large cohort. Collectively, this analysis reveals the functional heterogeneity and its association with genomic evolution in TNBC, and uncovers unanticipated biological principles dictating poor outcomes in this disease.

Surgical management of superior petrosal sinus dural arteriovenous fistulae with dominant internal carotid artery supply.

Background While technological advances have improved the efficacy of endovascular techniques for tentorial dural arteriovenous fistulae (DAVF), superior petrosal sinus (SPS) DAVF with dominant internal carotid artery (ICA) supply frequently require surgical intervention to achieve a definitive cure. Methods To compare the angiographic and clinical outcomes of endovascular and surgical interventions in patients with SPS DAVF, the records of all patients with tentorial DAVF from August 2010 to November 2015 were reviewed. Results Within this cohort, eight patients with nine SPS DAVF were eligible for evaluation. Five DAVF were initially treated with endovascular embolization, while four underwent surgical occlusion without embolization. Of the SPS DAVF treated with embolization, two (40%) remained occluded on follow-up, while the remaining three (60%) persisted/recurred and required surgical intervention for definitive closure. Of the four SPS DAVF treated with primary surgical occlusion, all four (100%) remained closed on follow-up. In addition, of the three SPS DAVF that persisted/recurred following embolization and required subsequent surgical closure, all three (100%) remained occluded on follow-up. Two (100%) SPS DAVF that were successfully treated with embolization had major or minor external carotid artery supply, while the three (100%) persistent lesions had major ICA supply via the meningohypophyseal trunk (MHT). Three (75%) of the four SPS DAVF treated with primary surgical occlusion had dominant MHT supply. Conclusion Complete endovascular closure of SPS DAVF with dominant ICA supply via the MHT may be difficult to achieve, while upfront surgical intervention is associated with a high rate of complete occlusion.

Rapid Intraoperative in Situ Synthetic Cranioplasty.

Craniectomy is a frequently performed neurosurgical procedure, and coverage of the cranial defect is necessary for protection of the underlying brain, cosmesis, and patient satisfaction. We report a new technique for intraoperative in situ synthetic cranioplasty that provides one-step resection of skull osteomas and reconstruction of cranial defects. Strategies of intraoperative cranioplasty are reviewed. A 48-year-old man who presented with a suspected benign osteoma over his forehead was offered surgical excision and primary cranioplasty in a one-step procedure using hydroxyapatite bone cement, a dural prosthetic, and a resorbable plate. Following craniectomy around the lesion, there was evidence of dural and bone involvement. The craniectomy was enlarged, and the involved dura was resected. SYNTHECEL dura repair was used to repair the dural defect and at the same time fashioned to form a receptacle for the cranioplasty by fixation of the dural substrate to the cut vertical bone edges. DirectInject hydroxyapatite bone cement was used to fill the receptacle and contoured to the curvature of the adjacent skull. A Delta resorbable plate was then placed over the bone cement and fixed to the skull. This technique provided a satisfactory cosmetic outcome following craniectomy for benign skull tumor excision. When possible, one-step surgery with primary cranioplasty should be considered.

Single-Cell Transcriptomic Analysis of Primary and Metastatic Tumor Ecosystems in Head and Neck Cancer.

The diverse malignant, stromal, and immune cells in tumors affect growth, metastasis, and response to therapy. We profiled transcriptomes of ∼6,000 single cells from 18 head and neck squamous cell carcinoma (HNSCC) patients, including five matched pairs of primary tumors and lymph node metastases. Stromal and immune cells had consistent expression programs across patients. Conversely, malignant cells varied within and between tumors in their expression of signatures related to cell cycle, stress, hypoxia, epithelial differentiation, and partial epithelial-to-mesenchymal transition (p-EMT). Cells expressing the p-EMT program spatially localized to the leading edge of primary tumors. By integrating single-cell transcriptomes with bulk expression profiles for hundreds of tumors, we refined HNSCC subtypes by their malignant and stromal composition and established p-EMT as an independent predictor of nodal metastasis, grade, and adverse pathologic features. Our results provide insight into the HNSCC ecosystem and define stromal interactions and a p-EMT program associated with metastasis.

Decoupling genetics, lineages, and microenvironment in IDH-mutant gliomas by single-cell RNA-seq.

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.

Adaptive Chromatin Remodeling Drives Glioblastoma Stem Cell Plasticity and Drug Tolerance.

Glioblastoma, the most common and aggressive malignant brain tumor, is propagated by stem-like cancer cells refractory to existing therapies. Understanding the molecular mechanisms that control glioblastoma stem cell (GSC) proliferation and drug resistance may reveal opportunities for therapeutic interventions. Here we show that GSCs can reversibly transition to a slow-cycling, persistent state in response to targeted kinase inhibitors. In this state, GSCs upregulate primitive developmental programs and are dependent upon Notch signaling. This transition is accompanied by widespread redistribution of repressive histone methylation. Accordingly, persister GSCs upregulate, and are dependent on, the histone demethylases KDM6A/B. Slow-cycling cells with high Notch activity and histone demethylase expression are present in primary glioblastomas before treatment, potentially contributing to relapse. Our findings illustrate how cancer cells may hijack aspects of native developmental programs for deranged proliferation, adaptation, and tolerance. They also suggest strategies for eliminating refractory tumor cells by targeting epigenetic and developmental pathways.

Single-cell RNA-seq supports a developmental hierarchy in human oligodendroglioma.

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