Identifying Gene-Treatment Interactions and Targetable Radiation Vulnerabilities in Glioblastoma through Coupling of In Vivo CRISPR Perturbation and Single Cell Transcriptomics

Abstract

Glioblastoma (GBM) is an incurable brain tumor comprised of dynamic malignant cell states and microenvironment components that underlie treatment resistance. Here we use genome-wide CRISPR/Cas9 functional genomics to define biological drivers and therapeutic vulnerabilities across human and mouse GBM models. To interrogate these mechanisms in the context of the tumor microenvironment and in vivo physiology, we established in vivo Perturb-seq intracranially, a technique coupling functional genomics with single cell transcriptomics, where each cell is an individual experiment. Orthotopic intracranial tumor models were established using human (GBM6, GBM43) or mouse (GL261, SB28) GBM cells stably expressing CRISPR interference (CRISPRi) machinery. Perturb-seq target selection for phenotyping of gene-treatment interactions was performed using genome-wide CRISPRi screens ± radiotherapy in cell cultures. Dual sgRNA lentivirus libraries were transduced either ex vivo prior to intracranial GBM cell transplantation or in vivo using intratumor convection enhanced delivery (CED). Transcriptional phenotyping was performed using single-cell RNA-seq with CRISPR direct capture following focal brain radiotherapy (2 Gy x 5) or mock treatment. GBM cell states were validated using single-nucleus RNA-seq data from 86 primary-recurrent patient-matched GBMs. Mechanistic and functional validation was performed using small molecule inhibitors, immunohistochemistry, clonogenic assays, and in vivo survival experiments. In vivo Perturb-seq ± radiotherapy of 48 genes underlying GBM radiotherapy responses, which were enriched for DNA damage response and metabolic pathways, was performed in > 425,000 single cells. Radiotherapy induced 16 distinct GBM cell states, and genetic perturbations reprogrammed these cell states in a treatment-dependent fashion. Quantitative modeling of gene/radiotherapy interactions using high dimensional manifolds revealed in vivo-specific genetic dependencies. We revealed a critical role for Prkdc, the catalytic subunit of DNA-dependent protein kinase (DNA-PK), as a radiotherapy sensitizer through regulation of cell intrinsic growth and oxidative stress pathways, and cell extrinsic interferon and signaling pathways that altered cell-cell interactions in vivo. These pathways were also disrupted in single-nucleus RNA-seq analysis of post-radiotherapy human GBM tumors. Inhibition of Prkdc using a Food and Drug Administration approved small molecule sensitized GBM cells to radiotherapy and extended survival in mice harboring human intracranial xenografts. We establish in vivo Perturb-seq in orthotopic GBM models as a platform for simultaneous functional genomic discovery and characterization of therapeutic targets, revealing an underappreciated role for Prkdc in GBM tumors in vivo that is targetable using small molecules. These tools are adaptable for a wide range of disease models and treatment modalities.