Cell state evolution underlies tumor development and response to therapy, but mechanisms specifying cancer cell states and intratumor heterogeneity are incompletely understood. Schwannomas are the most common tumors of the peripheral nervous system and are treated with surgery and ionizing radiation. Schwannomas can oscillate in size for many years after radiotherapy, suggesting treatment may reprogram schwannoma cells or the tumor microenvironment. Here we show epigenetic reprogramming shapes the cellular landscape of schwannomas. We find schwannomas are comprised of 2 molecular groups distinguished by reactivation of neural crest development pathways or misactivation of nerve injury mechanisms that specify cancer cell states and the architecture of the tumor immune microenvironment. Schwannoma molecular groups can arise independently, but ionizing radiation is sufficient for epigenetic reprogramming of neural crest to immune-enriched schwannoma by remodeling chromatin accessibility, gene expression, and metabolism to drive schwannoma cell state evolution and immune cell infiltration. To define functional genomic mechanisms underlying epigenetic reprograming of schwannomas, we develop a technique for simultaneous interrogation of chromatin accessibility and gene expression coupled with genetic and therapeutic perturbations in single-nuclei. Our results elucidate a framework for understanding epigenetic drivers of cancer evolution and establish a paradigm of epigenetic reprograming of cancer in response to radiotherapy.

Schwann cell derived tumors comprising schwannomas, neurofibromas, and malignant peripheral nerve sheath tumors (MPNSTs) are the most common cancers of the peripheral nervous system and often arise in patients with neurofibromatosis type-1 (NF-1) or type-2 (NF-2). NF-1 is caused by loss of NF1, a negative regulator of Ras signaling, and NF-2 is caused by loss of NF2, a pleiotropic tumor suppressor with numerous functions including inhibition of PAK signaling. However, whether functional interactions exist between the NF1 and NF2 tumor suppressors remain unclear. More broadly, there are currently no effective molecular therapies for patients with Schwann cell tumors beyond the MEK inhibitor selumetinib to treat neurofibromas in patients with NF-1. Here, we integrate DNA methylation profiling, whole exome sequencing, bulk and single-cell RNA sequencing, biochemistry, and pharmacology across human samples, cell lines, and mouse xenografts to identify cellular de-differentiation as a driver of malignant transformation and selumetinib resistance. Single nuclear RNA-sequencing of human neurofibromas (n = 3) or MPNSTs (n = 3) revealed a total of 13 cell types with increased proliferating, de-differentiated tumor cell populations in MPNST samples. Single cell RNA-sequencing of MPNST mouse xenografts revealed persistence of de-differentiated cell populations in selumetinib treated samples compared to vehicle control, suggesting cellular de-differentiation underlies treatment resistance. A genome-wide CRISPRi screen for mediators of selumetinib response in NF1 deficient neurofibroma cells revealed NF2 loss drives selumetinib resistance. Consistently, NF2 suppression in NF1 deficient neurofibroma cells caused Schwann cell de-differentiation and activation of PAK, a serine threonine kinase. Translationally, a small molecule PAK inhibitor in combination with selumetinib formed an effective therapy in mouse MPNST xenografts. In sum, we elucidate a paradigm of de-differentiation driving malignant transformation and treatment resistance in Schwann cell tumors, uncovering a functional link between the NF1 and NF2 tumor suppressors that sheds light on a novel therapeutic vulnerability.

DNA methylation profiling provides robust classification of nervous system tumors, but mechanisms driving epigenetic identity of individual tumor types are incompletely understood. Integrating DNA methylation profiling (n=76), RNA sequencing (n=24), single-cell RNA-sequencing (n=9), and mass cytometry (n=9), we discovered vestibular schwannomas are comprised of two epigenetic groups distinguished by neural crest development pathways or repair and regeneration pathways driving immune infiltration. Analyses of preoperative magnetic resonance imaging studies (n=66) or paired primary and recurrent schwannomas (n=13) suggested radiotherapy was sufficient but not necessary for epigenetic reprogramming of neural crest enriched schwannomas into immune enriched schwannomas. In support of this hypothesis, DNA methylation profiling, RNA sequencing, single-cell RNA sequencing, proteomic mass spectrometry, and lymphocyte migration assays demonstrated radiotherapy epigenetically reprogramed viable schwannoma cells to secrete immunomodulatory signals and recruit lymphocytes in vitro. Genome-wide CRISPRi screens identified histone acetyltransferases or DNA methyltransferases driving schwannoma radiotherapy responses, including the epigenetic regulators KDM5C or KDM1A. CRISPRi and lymphocyte migration assays ± radiotherapy confirmed KDM5C drives schwannoma immune infiltration whereas KDM1A inhibits schwannoma immune infiltration. To define genomic mechanisms underlying epigenetic group identity, we performed pooled CRISPRi screening coupled with single-cell RNA sequencing (Perturb-seq) of 44 schwannoma markers. In parallel, we developed single nuclei profiling of chromatin accessibility through paired ATAC sequencing and RNA sequencing coupled with pooled CRISPRi screening (snARC-seq) of 54 epigenetic regulators identified by our genome-wide CRISPRi screen. Functional genomic approaches revealed the tyrosine phosphatase PTPRG as a regulator of survival, and KDM5C and KDM1A as regulators of inflammation. In summary, we report two epigenetic groups of schwannomas and mechanisms underlying epigenetic group identity using a new functional genomic technique allowing for simultaneous interrogation of single-cell epigenetic and gene expression changes in the context of genetic and therapeutic perturbations. These data elucidate a novel epigenetic mechanism of action of radiotherapy.

Schwann cell tumors are the most common cancers of the peripheral nervous system and can arise sporadically or in patients with neurofibromatosis type-1 (NF-1) or type-2 (NF-2). NF-1 is caused by loss of NF1, a negative regulator of Ras signaling. NF-2 is caused by loss of NF2, a pleiotropic tumor suppressor that inhibits PAK signaling. Functional interactions between the NF1 and NF2 tumor suppressors and broader mechanisms underlying malignant transformation of the Schwann lineage are unclear. Here, we integrate DNA methylation profiling, whole exome sequencing, bulk and single-cell RNA sequencing, biochemistry, and pharmacology across human samples, patient-derived cell lines, and mouse xenografts to identify cellular de-differentiation mechanisms driving malignant transformation and treatment resistance in Schwann cell tumors. Our data show molecular groups of Schwann cell tumors are distinguished by de-differentiation trajectories that drive resistance to MEK inhibition, the only approved molecular therapy for patients with NF-1. Functional genomic screening for mediators of MEK inhibitor responses in NF1-deficient tumor cells reveals NF2 loss and PAK activation underlie Schwann cell tumor de-differentiation and MEK inhibitor resistance. In support of these findings, we identify a group of de-differentiated Schwann cell tumors with concurrent loss of NF1 and NF2, and find combination molecular therapy inhibiting MEK and PAK is an effective treatment for de-differentiated Schwann cell tumor xenografts. In sum, we elucidate a paradigm of de-differentiation driving malignant transformation and treatment resistance, uncovering a functional link between the NF1 and NF2 tumor suppressors that sheds light on a novel therapeutic vulnerability.