VC-resist glioblastoma cell state: vessel co-option as a key driver of chemoradiation resistance.

Glioblastoma (GBM) is a highly lethal type of cancer. GBM recurrence following chemoradiation is typically attributed to the regrowth of invasive and resistant cells. Therefore, there is a pressing need to gain a deeper understanding of the mechanisms underlying GBM resistance to chemoradiation and its ability to infiltrate. Using a combination of transcriptomic, proteomic, and phosphoproteomic analyses, longitudinal imaging, organotypic cultures, functional assays, animal studies, and clinical data analyses, we demonstrate that chemoradiation and brain vasculature induce cell transition to a functional state named VC-Resist (vessel co-opting and resistant cell state). This cell state is midway along the transcriptomic axis between proneural and mesenchymal GBM cells and is closer to the AC/MES1-like state. VC-Resist GBM cells are highly vessel co-opting, allowing significant infiltration into the surrounding brain tissue and homing to the perivascular niche, which in turn induces even more VC-Resist transition. The molecular and functional characteristics of this FGFR1-YAP1-dependent GBM cell state, including resistance to DNA damage, enrichment in the G2M phase, and induction of senescence/stemness pathways, contribute to its enhanced resistance to chemoradiation. These findings demonstrate how vessel co-option, perivascular niche, and GBM cell plasticity jointly drive resistance to therapy during GBM recurrence.

ABL1 and ABL2 promote medulloblastoma leptomeningeal dissemination.

Background: Medulloblastoma is the most common malignant pediatric brain tumor, and leptomeningeal dissemination (LMD) of medulloblastoma both portends a poorer prognosis at diagnosis and is incurable at recurrence. The biological mechanisms underlying LMD are unclear. The Abelson (ABL) tyrosine kinase family members, ABL1 and ABL2, have been implicated in cancer cell migration, invasion, adhesion, metastasis, and chemotherapy resistance, and are upstream mediators of the oncogene c-MYC in fibroblasts and lung cancer cells. However, their role in medulloblastoma has not yet been explored. The purpose of this work was to elucidate the role of ABL1/2 in medulloblastoma LMD. Methods: ABL1 and ABL2 mRNA expression of patient specimens was analyzed. shRNA knockdowns of ABL1/2 and pharmacologic inhibition of ABL1/2 were used for in vitro and in vivo analyses of medulloblastoma LMD. RNA sequencing of ABL1/2 genetic knockdown versus scrambled control medulloblastoma was completed. Results: ABL1/2 mRNA is highly expressed in human medulloblastoma and pharmacologic inhibition of ABL kinases resulted in cytotoxicity. Knockdown of ABL1/2 resulted in decreased adhesion of medulloblastoma cells to the extracellular matrix protein, vitronectin (P = .0013), and significantly decreased tumor burden in a mouse model of medulloblastoma LMD with improved overall survival (P = .0044). Furthermore, both pharmacologic inhibition of ABL1/2 and ABL1/2 knockdown resulted in decreased expression of c-MYC, identifying a putative signaling pathway, and genes/pathways related to oncogenesis and neurodevelopment were differentially expressed between ABL1/2 knockdown and control medulloblastoma cells. Conclusions: ABL1 and ABL2 have potential roles in medulloblastoma LMD upstream of c-MYC expression.

Inhibition of the pseudokinase MLKL alters extracellular vesicle release and reduces tumor growth in glioblastoma.

Extracellular vesicles (EVs) are lipid-based nanosized particles that convey biological material from donor to recipient cells. EVs play key roles in glioblastoma progression because glioblastoma stem-like cells (GSCs) release pro-oncogenic, pro-angiogenic, and pro-inflammatory EVs. However, the molecular basis of EV release remains poorly understood. Here, we report the identification of the pseudokinase MLKL, a crucial effector of cell death by necroptosis, as a regulator of the constitutive secretion of EVs in GSCs. We find that genetic, protein, and pharmacological targeting of MLKL alters intracellular trafficking and EV release, and reduces GSC expansion. Nevertheless, this function ascribed to MLKL appears independent of its role during necroptosis. In vivo, pharmacological inhibition of MLKL reduces the tumor burden and the level of plasmatic EVs. This work highlights the necroptosis-independent role of MLKL in vesicle release and suggests that interfering with EVs is a promising therapeutic option to sensitize glioblastoma cells.

The transcriptional landscape of Shh medulloblastoma.

Sonic hedgehog medulloblastoma encompasses a clinically and molecularly diverse group of cancers of the developing central nervous system. Here, we use unbiased sequencing of the transcriptome across a large cohort of 250 tumors to reveal differences among molecular subtypes of the disease, and demonstrate the previously unappreciated importance of non-coding RNA transcripts. We identify alterations within the cAMP dependent pathway (GNAS, PRKAR1A) which converge on GLI2 activity and show that 18% of tumors have a genetic event that directly targets the abundance and/or stability of MYCN. Furthermore, we discover an extensive network of fusions in focally amplified regions encompassing GLI2, and several loss-of-function fusions in tumor suppressor genes PTCH1, SUFU and NCOR1. Molecular convergence on a subset of genes by nucleotide variants, copy number aberrations, and gene fusions highlight the key roles of specific pathways in the pathogenesis of Sonic hedgehog medulloblastoma and open up opportunities for therapeutic intervention.

PRMT5 inhibition disrupts splicing and stemness in glioblastoma.

Glioblastoma (GBM) is a deadly cancer in which cancer stem cells (CSCs) sustain tumor growth and contribute to therapeutic resistance. Protein arginine methyltransferase 5 (PRMT5) has recently emerged as a promising target in GBM. Using two orthogonal-acting inhibitors of PRMT5 (GSK591 or LLY-283), we show that pharmacological inhibition of PRMT5 suppresses the growth of a cohort of 46 patient-derived GBM stem cell cultures, with the proneural subtype showing greater sensitivity. We show that PRMT5 inhibition causes widespread disruption of splicing across the transcriptome, particularly affecting cell cycle gene products. We identify a GBM splicing signature that correlates with the degree of response to PRMT5 inhibition. Importantly, we demonstrate that LLY-283 is brain-penetrant and significantly prolongs the survival of mice with orthotopic patient-derived xenografts. Collectively, our findings provide a rationale for the clinical development of brain penetrant PRMT5 inhibitors as treatment for GBM.

Gradient of Developmental and Injury Response transcriptional states defines functional vulnerabilities underpinning glioblastoma heterogeneity.

Glioblastomas harbor diverse cell populations, including rare glioblastoma stem cells (GSCs) that drive tumorigenesis. To characterize functional diversity within this population, we performed single-cell RNA sequencing on >69,000 GSCs cultured from the tumors of 26 patients. We observed a high degree of inter- and intra-GSC transcriptional heterogeneity that could not be fully explained by DNA somatic alterations. Instead, we found that GSCs mapped along a transcriptional gradient spanning two cellular states reminiscent of normal neural development and inflammatory wound response. Genome-wide CRISPR-Cas9 dropout screens independently recapitulated this observation, with each state characterized by unique essential genes. Further single-cell RNA sequencing of >56,000 malignant cells from primary tumors found that the majority organize along an orthogonal astrocyte maturation gradient yet retain expression of founder GSC transcriptional programs. We propose that glioblastomas grow out of a fundamental GSC-based neural wound response transcriptional program, which is a promising target for new therapy development.

Metabolic Regulation of the Epigenome Drives Lethal Infantile Ependymoma.

Posterior fossa A (PFA) ependymomas are lethal malignancies of the hindbrain in infants and toddlers. Lacking highly recurrent somatic mutations, PFA ependymomas are proposed to be epigenetically driven tumors for which model systems are lacking. Here we demonstrate that PFA ependymomas are maintained under hypoxia, associated with restricted availability of specific metabolites to diminish histone methylation, and increase histone demethylation and acetylation at histone 3 lysine 27 (H3K27). PFA ependymomas initiate from a cell lineage in the first trimester of human development that resides in restricted oxygen. Unlike other ependymomas, transient exposure of PFA cells to ambient oxygen induces irreversible cellular toxicity. PFA tumors exhibit a low basal level of H3K27me3, and, paradoxically, inhibition of H3K27 methylation specifically disrupts PFA tumor growth. Targeting metabolism and/or the epigenome presents a unique opportunity for rational therapy for infants with PFA ependymoma.

Tumor necrosis factor overcomes immune evasion in p53-mutant medulloblastoma.

Many immunotherapies act by enhancing the ability of cytotoxic T cells to kill tumor cells. Killing depends on T cell recognition of antigens presented by class I major histocompatibility complex (MHC-I) proteins on tumor cells. In this study, we showed that medulloblastomas lacking the p53 tumor suppressor do not express surface MHC-I and are therefore resistant to immune rejection. Mechanistically, this is because p53 regulates expression of the peptide transporter Tap1 and the aminopeptidase Erap1, which are required for MHC-I trafficking to the cell surface. In vitro, tumor necrosis factor (TNF) or lymphotoxin-ß receptor agonist can rescue expression of Erap1, Tap1 and MHC-I on p53-mutant tumor cells. In vivo, low doses of TNF prolong survival and synergize with immune checkpoint inhibitors to promote tumor rejection. These studies identified p53 as a key regulator of immune evasion and suggest that TNF could be used to enhance sensitivity of tumors to immunotherapy.

Locoregional delivery of CAR T cells to the cerebrospinal fluid for treatment of metastatic medulloblastoma and ependymoma.

Recurrent medulloblastoma and ependymoma are universally lethal, with no approved targeted therapies and few candidates presently under clinical evaluation. Nearly all recurrent medulloblastomas and posterior fossa group A (PFA) ependymomas are located adjacent to and bathed by the cerebrospinal fluid, presenting an opportunity for locoregional therapy, bypassing the blood-brain barrier. We identify three cell-surface targets, EPHA2, HER2 and interleukin 13 receptor α2, expressed on medulloblastomas and ependymomas, but not expressed in the normal developing brain. We validate intrathecal delivery of EPHA2, HER2 and interleukin 13 receptor α2 chimeric antigen receptor T cells as an effective treatment for primary, metastatic and recurrent group 3 medulloblastoma and PFA ependymoma xenografts in mouse models. Finally, we demonstrate that administration of these chimeric antigen receptor T cells into the cerebrospinal fluid, alone or in combination with azacytidine, is a highly effective therapy for multiple metastatic mouse models of group 3 medulloblastoma and PFA ependymoma, thereby providing a rationale for clinical trials of these approaches in humans.

Stalled developmental programs at the root of pediatric brain tumors.

Childhood brain tumors have suspected prenatal origins. To identify vulnerable developmental states, we generated a single-cell transcriptome atlas of >65,000 cells from embryonal pons and forebrain, two major tumor locations. We derived signatures for 191 distinct cell populations and defined the regional cellular diversity and differentiation dynamics. Projection of bulk tumor transcriptomes onto this dataset shows that WNT medulloblastomas match the rhombic lip-derived mossy fiber neuronal lineage and embryonal tumors with multilayered rosettes fully recapitulate a neuronal lineage, while group 2a/b atypical teratoid/rhabdoid tumors may originate outside the neuroectoderm. Importantly, single-cell tumor profiles reveal highly defined cell hierarchies that mirror transcriptional programs of the corresponding normal lineages. Our findings identify impaired differentiation of specific neural progenitors as a common mechanism underlying these pediatric cancers and provide a rational framework for future modeling and therapeutic interventions.

Recurrent noncoding U1 snRNA mutations drive cryptic splicing in SHH medulloblastoma.

In cancer, recurrent somatic single-nucleotide variants-which are rare in most paediatric cancers-are confined largely to protein-coding genes1-3. Here we report highly recurrent hotspot mutations (r.3A>G) of U1 spliceosomal small nuclear RNAs (snRNAs) in about 50% of Sonic hedgehog (SHH) medulloblastomas. These mutations were not present across other subgroups of medulloblastoma, and we identified these hotspot mutations in U1 snRNA in only <0.1% of 2,442 cancers, across 36 other tumour types. The mutations occur in 97% of adults (subtype SHHδ) and 25% of adolescents (subtype SHHα) with SHH medulloblastoma, but are largely absent from SHH medulloblastoma in infants. The U1 snRNA mutations occur in the 5' splice-site binding region, and snRNA-mutant tumours have significantly disrupted RNA splicing and an excess of 5' cryptic splicing events. Alternative splicing mediated by mutant U1 snRNA inactivates tumour-suppressor genes (PTCH1) and activates oncogenes (GLI2 and CCND2), and represents a target for therapy. These U1 snRNA mutations provide an example of highly recurrent and tissue-specific mutations of a non-protein-coding gene in cancer.

Engineering Genetic Predisposition in Human Neuroepithelial Stem Cells Recapitulates Medulloblastoma Tumorigenesis.

Human neural stem cell cultures provide progenitor cells that are potential cells of origin for brain cancers. However, the extent to which genetic predisposition to tumor formation can be faithfully captured in stem cell lines is uncertain. Here, we evaluated neuroepithelial stem (NES) cells, representative of cerebellar progenitors. We transduced NES cells with MYCN, observing medulloblastoma upon orthotopic implantation in mice. Significantly, transcriptomes and patterns of DNA methylation from xenograft tumors were globally more representative of human medulloblastoma compared to a MYCN-driven genetically engineered mouse model. Orthotopic transplantation of NES cells generated from Gorlin syndrome patients, who are predisposed to medulloblastoma due to germline-mutated PTCH1, also generated medulloblastoma. We engineered candidate cooperating mutations in Gorlin NES cells, with mutation of DDX3X or loss of GSE1 both accelerating tumorigenesis. These findings demonstrate that human NES cells provide a potent experimental resource for dissecting genetic causation in medulloblastoma.

High-resolution structural genomics reveals new therapeutic vulnerabilities in glioblastoma.

We investigated the role of 3D genome architecture in instructing functional properties of glioblastoma stem cells (GSCs) by generating sub-5-kb resolution 3D genome maps by in situ Hi-C. Contact maps at sub-5-kb resolution allow identification of individual DNA loops, domain organization, and large-scale genome compartmentalization. We observed differences in looping architectures among GSCs from different patients, suggesting that 3D genome architecture is a further layer of inter-patient heterogeneity for glioblastoma. Integration of DNA contact maps with chromatin and transcriptional profiles identified specific mechanisms of gene regulation, including the convergence of multiple super enhancers to individual stemness genes within individual cells. We show that the number of loops contacting a gene correlates with elevated transcription. These results indicate that stemness genes are hubs of interaction between multiple regulatory regions, likely to ensure their sustained expression. Regions of open chromatin common among the GSCs tested were poised for expression of immune-related genes, including CD276 We demonstrate that this gene is co-expressed with stemness genes in GSCs and that CD276 can be targeted with an antibody-drug conjugate to eliminate self-renewing cells. Our results demonstrate that integrated structural genomics data sets can be employed to rationally identify therapeutic vulnerabilities in self-renewing cells.

Childhood cerebellar tumours mirror conserved fetal transcriptional programs.

Study of the origin and development of cerebellar tumours has been hampered by the complexity and heterogeneity of cerebellar cells that change over the course of development. Here we use single-cell transcriptomics to study more than 60,000 cells from the developing mouse cerebellum and show that different molecular subgroups of childhood cerebellar tumours mirror the transcription of cells from distinct, temporally restricted cerebellar lineages. The Sonic Hedgehog medulloblastoma subgroup transcriptionally mirrors the granule cell hierarchy as expected, while group 3 medulloblastoma resembles Nestin+ stem cells, group 4 medulloblastoma resembles unipolar brush cells, and PFA/PFB ependymoma and cerebellar pilocytic astrocytoma resemble the prenatal gliogenic progenitor cells. Furthermore, single-cell transcriptomics of human childhood cerebellar tumours demonstrates that many bulk tumours contain a mixed population of cells with divergent differentiation. Our data highlight cerebellar tumours as a disorder of early brain development and provide a proximate explanation for the peak incidence of cerebellar tumours in early childhood.

Genome-Wide CRISPR-Cas9 Screens Expose Genetic Vulnerabilities and Mechanisms of Temozolomide Sensitivity in Glioblastoma Stem Cells.

Glioblastoma therapies have remained elusive due to limitations in understanding mechanisms of growth and survival of the tumorigenic population. Using CRISPR-Cas9 approaches in patient-derived GBM stem cells (GSCs) to interrogate function of the coding genome, we identify actionable pathways responsible for growth, which reveal the gene-essential circuitry of GBM stemness and proliferation. In particular, we characterize members of the SOX transcription factor family, SOCS3, USP8, and DOT1L, and protein ufmylation as important for GSC growth. Additionally, we reveal mechanisms of temozolomide resistance that could lead to combination strategies. By reaching beyond static genome analysis of bulk tumors, with a genome-wide functional approach, we reveal genetic dependencies within a broad range of biological processes to provide increased understanding of GBM growth and treatment resistance.

Intratumoral Genetic and Functional Heterogeneity in Pediatric Glioblastoma.

Pediatric glioblastoma (pGBM) is a lethal cancer with no effective therapies. To understand the mechanisms of tumor evolution in this cancer, we performed whole-genome sequencing with linked reads on longitudinally resected pGBM samples. Our analyses showed that all diagnostic and recurrent samples were collections of genetically diverse subclones. Clonal composition rapidly evolved at recurrence, with less than 8% of nonsynonymous single-nucleotide variants being shared in diagnostic-recurrent pairs. To track the origins of the mutational events observed in pGBM, we generated whole-genome datasets for two patients and their parents. These trios showed that genetic variants could be (i) somatic, (ii) inherited from a healthy parent, or (iii) de novo in the germlines of pGBM patients. Analysis of variant allele frequencies supported a model of tumor growth involving slow-cycling cancer stem cells that give rise to fast-proliferating progenitor-like cells and to nondividing cells. Interestingly, radiation and antimitotic chemotherapeutics did not increase overall tumor burden upon recurrence. These findings support an important role for slow-cycling stem cell populations in contributing to recurrences, because slow-cycling cell populations are expected to be less prone to genotoxic stress induced by these treatments and therefore would accumulate few mutations. Our results highlight the need for new targeted treatments that account for the complex functional hierarchies and genomic heterogeneity of pGBM. SIGNIFICANCE: This work challenges several assumptions regarding the genetic organization of pediatric GBM and highlights mutagenic programs that start during early prenatal development.Graphical Abstract: http://cancerres.aacrjournals.org/content/canres/79/9/2111/F1.large.jpg.

Heterogeneity within the PF-EPN-B ependymoma subgroup.

Posterior fossa ependymoma comprise three distinct molecular variants, termed PF-EPN-A (PFA), PF-EPN-B (PFB), and PF-EPN-SE (subependymoma). Clinically, they are very disparate and PFB tumors are currently being considered for a trial of radiation avoidance. However, to move forward, unraveling the heterogeneity within PFB would be highly desirable. To discern the molecular heterogeneity within PFB, we performed an integrated analysis consisting of DNA methylation profiling, copy-number profiling, gene expression profiling, and clinical correlation across a cohort of 212 primary posterior fossa PFB tumors. Unsupervised spectral clustering and t-SNE analysis of genome-wide methylation data revealed five distinct subtypes of PFB tumors, termed PFB1-5, with distinct demographics, copy-number alterations, and gene expression profiles. All PFB subtypes were distinct from PFA and posterior fossa subependymomas. Of the five subtypes, PFB4 and PFB5 are more discrete, consisting of younger and older patients, respectively, with a strong female-gender enrichment in PFB5 (age: p = 0.011, gender: p = 0.04). Broad copy-number aberrations were common; however, many events such as chromosome 2 loss, 5 gain, and 17 loss were enriched in specific subtypes and 1q gain was enriched in PFB1. Late relapses were common across all five subtypes, but deaths were uncommon and present in only two subtypes (PFB1 and PFB3). Unlike the case in PFA ependymoma, 1q gain was not a robust marker of poor progression-free survival; however, chromosome 13q loss may represent a novel marker for risk stratification across the spectrum of PFB subtypes. Similar to PFA ependymoma, there exists a significant intertumoral heterogeneity within PFB, with distinct molecular subtypes identified. Even when accounting for this heterogeneity, extent of resection remains the strongest predictor of poor outcome. However, this biological heterogeneity must be accounted for in future preclinical modeling and personalized therapies.

Fate mapping of human glioblastoma reveals an invariant stem cell hierarchy.

Human glioblastomas harbour a subpopulation of glioblastoma stem cells that drive tumorigenesis. However, the origin of intratumoural functional heterogeneity between glioblastoma cells remains poorly understood. Here we study the clonal evolution of barcoded glioblastoma cells in an unbiased way following serial xenotransplantation to define their individual fate behaviours. Independent of an evolving mutational signature, we show that the growth of glioblastoma clones in vivo is consistent with a remarkably neutral process involving a conserved proliferative hierarchy rooted in glioblastoma stem cells. In this model, slow-cycling stem-like cells give rise to a more rapidly cycling progenitor population with extensive self-maintenance capacity, which in turn generates non-proliferative cells. We also identify rare 'outlier' clones that deviate from these dynamics, and further show that chemotherapy facilitates the expansion of pre-existing drug-resistant glioblastoma stem cells. Finally, we show that functionally distinct glioblastoma stem cells can be separately targeted using epigenetic compounds, suggesting new avenues for glioblastoma-targeted therapy.