Summary Gliomas synaptically integrate into neural circuits. Prior work has demonstrated bidirectional interactions between neurons and glioma cells, with neuronal activity driving glioma growth and gliomas increasing neuronal excitability. In this study we wanted to know how glioma induced neuronal changes influence neural circuits underlying cognition and whether these interactions influence patient survival. We use intracranial brain recordings during lexical retrieval language tasks in awake humans in addition to site specific tumor tissue biopsies and cell biology experiments. We find that gliomas remodel functional neural circuitry such that task-relevant neural responses activate tumor-infiltrated cortex, beyond cortical excitation normally recruited in the healthy brain. Site-directed biopsies from functionally connected regions within the tumor are enriched for a glioblastoma subpopulation that exhibits a distinct synaptogenic and neuronotrophic phenotype. Tumor cells from functionally connected regions secrete the synaptogenic factor thrombospondin-1, which contributes to the differential neuron-glioma interactions observed in functionally connected tumor regions compared to tumor regions with less functional connectivity. The degree of functional connectivity between glioblastoma and the normal brain negatively impacts both patient survival and language task performance. These data demonstrate that high-grade gliomas functionally remodel neural circuits in the human brain, which both promotes tumor proliferation and impairs cognition.

Abstract Genetic perturbation screens with single cell readouts have enabled rich phenotyping of gene function and regulatory networks. These approaches have been challenging in vivo, especially in adult disease models such as cancer, which include mixtures of malignant and microenvironment cells. Glioblastoma (GBM) is a fatal cancer, and methods of systematically interrogating gene function and therapeutic targets in vivo , especially in combination with standard of care treatment such as radiotherapy, are lacking. Here, we iteratively develop a multiplex in vivo perturb-seq CRISPRi platform for single cell genetic screens in cancer and tumor microenvironment cells that leverages intracranial convection enhanced delivery (CED) of sgRNA libraries into models of GBM. Our platform enables potent silencing of drivers of in vivo growth and tumor maintenance, as well as genes that sensitize GBM to radiotherapy. We find radiotherapy rewires transcriptional responses to genetic perturbations in an in vivo dependent manner, revealing heterogenous patterns of treatment sensitization or resistance in GBM. Furthermore, we demonstrate targeting of genes that function in the tumor microenvironment, enabling alterations of ligand-receptor interactions between immune/stromal cells following in vivo CRISPRi perturbations. In sum, we demonstrate the utility of multiplexed perturb-seq for in vivo single cell dissection of adult cancer and normal tissue biology across multiple cell types in the context of therapeutic intervention, a platform with potential for broad application.

NF1 is recurrently mutated in glioblastoma yet the molecular landscape and efficacy of targeted therapies remain unclear. Here, we combine bulk and single cell genomics of human somatic NF1 mutant, IDH-wildtype glioblastomas with functional genomic analysis of cell lines and mouse intracranial tumor models to identify molecular subgroups within NF1 mutant glioblastomas and mechanisms underlying MEK inhibitor response. Targeted DNA sequencing showed homozygous deletion of the cell cycle regulator CDKN2A/B is a poor prognostic marker in somatic NF1 mutant, but not NF1 wildtype, tumors. DNA methylation array profiling revealed three epigenetic groups highlighted by distinct clinical features, co-mutation patterns, and reference methylation classifier identities. Genome-wide CRISPRi screens in glioblastoma cells revealed cell cycle regulators are conserved mediators of cell growth while response to the MEK inhibitor selumetinib converges on Ras/RAF/MEK pathway activation. Repression of the RAF regulator SHOC2 sensitizes glioblastomas to selumetinib in vitro and in vivo in mouse intracranial glioblastoma models. Single cell RNA-sequencing of mouse intracranial glioblastomas treated with the MEK inhibitor selumetinib reveals distinct responses between mesenchymal-like (MES-like) and non MES-like subpopulations suggesting non-MES like cells are intrinsically resistant to MEK inhibition. Finally, single nuclear RNA-sequencing (snRNA-seq) of human NF1 mutant, CDKN2A/B deleted glioblastomas reveals MES-like tumor cells are associated with selumetinib sensitivity signatures while non MES-like cells exhibit increased cell cycle progression and lack selumetinib sensitivity, further supporting the notion that MEK inhibition specifically targets MES-like tumor cell subpopulations. Taken together, our data underscores the heterogeneity between and within somatic NF1 mutant glioblastomas and delineates mechanisms of MEK inhibitor response across distinct tumor subpopulations, guiding the development of future therapeutic strategies that may synergize with MEK inhibition for NF1 mutant tumors. The tumor suppressor NF1 is mutated in 15% of glioblastomas, 1–3 the most common malignant brain tumor with poor outcomes and few effective treatments. 4 NF1 is a GTPase activating protein (GAP) that negatively regulates Ras, and thus, NF1 loss leads to induction of Ras/RAF/MEK/ERK signaling, driving tumorigenesis and comprising a targetable molecular cascade. 5,6 Genomic analysis of glioblastoma demonstrates NF1 mutation is associated with a mesenchymal-like (MES-like) transcriptomic tumor cell subpopulation and altered tumor microenvironment. 7,8 More broadly, DNA methylation analysis reveals multiple epigenetic subgroups with overlapping relationships to transcriptomic subtype and DNA alterations, underscoring the complex relationship between genetic drivers and molecular signatures. 9 While the updated 2021 Central Nervous System WHO tumor classification incorporates an ever increasing amount of molecular criteria for diffuse astrocytic tumors, 10 the existence and clinical significance of molecular subgroups within somatic NF1 mutant, IDH-wildtype glioblastomas based on genetic co-mutations, epigenetic profile, or transcriptomic signatures remain unclear. Preclinical data support the utility of MEK inhibition in NF1 mutant gliomas, 11,12 and the MEK inhibitor selumetinib is FDA approved for tumors arising in patients with syndromic neurofibromatosis type 1 (NF-1) harboring a germline NF1 mutation. 13,14 In NF-1 associated gliomas, MEK inhibition demonstrates efficacy in a limited case series, 15 and combined BRAF/MEK inhibition shows efficacy in BRAF p.V600E mutant gliomas, 16 further supporting the translational potential of Ras/Raf/MEK/ERK blockade within genetically defined glioma subtypes. Nevertheless, treatment resistance to molecular monotherapy remains a challenge, 17–20 and the mechanisms underlying MEK inhibitor resistance in NF1 mutant glioma are unknown. Here, we integrate targeted DNA sequencing, DNA methylation profiling, and single nuclear RNA-sequencing (snRNA-seq) of human patient somatic NF1 mutant, IDH-wildtype glioblastomas with single cell RNA-sequencing (scRNA-seq), genome-wide clustered regularly interspaced short palindromic repeats interference (CRISPRi) screens, and pharmacologic studies in cell lines and mouse intracranial glioblastoma models to define molecular subgroups and functional mediators of MEK inhibitor response. Targeted DNA sequencing of NF1 mutant, IDH-wildtype glioblastomas (n=186 tumors) revealed CDKN2A/B deletion was associated with poor outcomes in NF1 mutant, but not NF1 wildtype, glioblastomas. DNA methylation profiling (n=129 tumors) demonstrated three epigenetic subgroups with distinct clinical features, co- mutation patterns across cell cycle genes, and reference methylation classifier identities. Genome-wide CRISPRi screens in mouse SB28 and human GBM43 glioblastoma cells identified a conserved cell cycle gene network mediating cell growth, consistent with the clinical importance of additional hits affecting the cell cycle in human somatic NF1 mutant glioblastomas. Moreover, genome-wide mediators of selumetinib response converged upon two Ras pathway effectors mediating selumetinib sensitivity: BRAF and SHOC2. SHOC2 repression in glioblastoma cells significantly improved selumetinib response both in vitro and in intracranial allografts in vivo . Single cell RNA-sequencing (scRNA-seq) of mouse intracranial glioblastomas treated with the MEK inhibitor selumetinib revealed MES-like tumor cells correlated with CDKN2A retention and the CRISPRi screen selumetinib sensitivity signature, with selumetinib resistant cells displaying Ras pathway induction. In contrast, non-MES like tumor cells were CDKN2A deficient and lacked expression of the CRISPRi screen selumetinib sensitivity signature, with selumetinib resistant cells inducing a glial de-differentiation program. Finally, snRNA-seq of NF1 mutant, CDKN2A/B deleted, IDH-wildtype glioblastomas (n=9) showed non MES-like tumor cells exhibit increased cell cycle progression and were not associated with the CRISPRi screen selumetinib sensitivity signature. MES-like tumor cells within newly diagnosed, but not recurrent, tumors retained expression of the CRISPRi screen selumetinib sensitivity signature, suggesting resistance can arise both between and within specific transcriptomic glioblastoma cell tumor cell subpopulations. Taken together, our data identifies clinically important subgroups of NF1 mutant, IDH-wildtype glioblastomas and supports a model in which heterogeneity between tumors and within tumor cell subpopulations underlies MEK inhibitor response, supporting the need for additional synergistic therapeutic approaches beyond maximal Ras pathway blockade for NF1 mutant glioblastomas.

Summary Loss of function (LoF) of Tar-binding protein 43 (TDP-43) and mislocalization, together with TDP-43-positive and hyperphosphorylated inclusions, are found in postmortem tissue of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) patients, including those carrying LoF variants in the progranulin gene ( GRN ). Modelling TDP-43 pathology has been challenging in vivo and in vitro . We present a 3D-induced pluripotent stem cell (iPSC)-derived paradigm - mature brain organoids (mbOrg) - composed of cortical-like-astrocytes (iA) and neurons (iN). When devoid of GRN , mbOrgs spontaneously recapitulate TDP-43 mislocalization, hyperphosphorylation and LoF phenotypes. Mixing-and-matching genotypes in mbOrgs showed that GRN −/− iA are drivers for TDP-43 pathology. Finally, we rescued TDP-43 LoF by adding exogenous progranulin, demonstrating a link between TDP-43 LoF and progranulin expression. In conclusion, we present an iPSC-derived platform that shows striking features of human TDP-43 proteinopathy and provides a tool for mechanistic modelling of TDP-43 pathology and patient-tailored therapeutic screening for FTD and ALS. Highlights GRN −/− iPSC-derived 3D paradigm (mbOrg) composed of mature cortical-like astrocytes and neurons recapitulates features of TDP-43 proteinopathy GRN −/− cortical-like astrocytes are necessary and sufficient for the development of the TDP-43 loss of function phenotype in mbOrg. A TDP-43 phenotype can be rescued in neurons by treating neuron and astrocyte co-cultures with progranulin full length protein. eTOC blurb In this article, de Majo and colleagues present a novel 3D iPSC-derived model to study neurodegenerative disorders such as ALS and FTD. When devoid of GRN expression, these cultures present features of ALS and FTD associated pathology hardly ever observed in vitro . These phenotypes are shown to be primarily driven by diseased astrocytes and can be rescued by progranulin supplementation.