RNA splicing analysis deciphers developmental hierarchies and reveals therapeutic targets in adult glioma.

Widespread alterations in RNA alternative splicing (AS) have been identified in adult gliomas. However, their regulatory mechanism, biological significance, and therapeutic potential remain largely elusive. Here, using a computational approach with both bulk and single cell RNA-sequencing, we uncover a prognostic AS signature linked with neural developmental hierarchies. Using advanced iPSC glioma models driven by glioma driver mutations, we show that this AS signature could be enhanced by EGFRvIII and inhibited by in situ IDH1 mutation. Functional validation of two isoform switching events in CERS5 and MPZL1 shows regulations of sphingolipid metabolism and SHP2 signaling, respectively. Analysis of upstream RNA binding proteins reveals PTBP1 as a key regulator of the AS signature where targeting of PTBP1 suppresses tumor growth and promotes the expression of a neuron marker TUJ1 in glioma stem-like cells. Overall, our data highlights the role of AS in impacting glioma malignance and heterogeneity and its potential as a therapeutic vulnerability for treating adult gliomas.

Homeostatic iron regulatory protein drives glioblastoma growth via tumor cell-intrinsic and sex-specific responses.

Background: Glioblastoma (GBM) displays alterations in iron that drive proliferation and tumor growth. Iron regulation is complex and involves many regulatory mechanisms, including the homeostatic iron regulator (HFE) gene, which encodes the homeostatic iron regulatory protein. While HFE is upregulated in GBM and correlates with poor survival outcomes, the function of HFE in GBM remains unclear. Methods: We interrogated the impact of cell-intrinsic Hfe expression on proliferation and survival of intracranially implanted animals through genetic gain- and loss-of-function approaches in syngeneic mouse glioma models, along with in vivo immune assessments. We also determined the expression of iron-associated genes and their relationship to survival in GBM using public data sets and used transcriptional profiling to identify differentially expressed pathways in control compared to Hfe-knockdown cells. Results: Overexpression of Hfe accelerated GBM proliferation and reduced animal survival, whereas suppression of Hfe induced apoptotic cell death and extended survival, which was more pronounced in females and associated with attenuation of natural killer cells and CD8+ T cell activity. Analysis of iron gene signatures in Hfe-knockdown cells revealed alterations in the expression of several iron-associated genes, suggesting global disruption of intracellular iron homeostasis. Further analysis of differentially expressed pathways revealed oxidative stress as the top pathway upregulated following Hfe loss. Hfe knockdown indeed resulted in enhanced 55Fe uptake and generation of reactive oxygen species. Conclusions: These findings reveal an essential function for HFE in GBM cell growth and survival, as well as a sex-specific interaction with the immune response.

Patient-derived glioblastoma organoids reflect tumor heterogeneity and treatment sensitivity.

Background: Treatment resistance and tumor relapse are the primary causes of mortality in glioblastoma (GBM), with intratumoral heterogeneity playing a significant role. Patient-derived cancer organoids have emerged as a promising model capable of recapitulating tumor heterogeneity. Our objective was to develop patient-derived GBM organoids (PGO) to investigate treatment response and resistance. Methods: GBM samples were used to generate PGOs and analyzed using whole-exome sequencing (WES) and single-cell karyotype sequencing. PGOs were subjected to temozolomide (TMZ) to assess viability. Bulk RNA sequencing was performed before and after TMZ. Results: WES analysis on individual PGOs cultured for 3 time points (1-3 months) showed a high inter-organoid correlation and retention of genetic variants (range 92.3%-97.7%). Most variants were retained in the PGO compared to the tumor (range 58%-90%) and exhibited similar copy number variations. Single-cell karyotype sequencing demonstrated preservation of genetic heterogeneity. Single-cell multiplex immunofluorescence showed maintenance of cellular states. TMZ treatment of PGOs showed a differential response, which largely corresponded with MGMT promoter methylation. Differentially expressed genes before and after TMZ revealed an upregulation of the JNK kinase pathway. Notably, the combination treatment of a JNK kinase inhibitor and TMZ demonstrated a synergistic effect. Conclusions: Overall, these findings demonstrate the robustness of PGOs in retaining the genetic and phenotypic heterogeneity in culture and the application of measuring clinically relevant drug responses. These data show that PGOs have the potential to be further developed into avatars for personalized adaptive treatment selection and actionable drug target discovery and as a platform to study GBM biology.

A TOX-ic axis of epigenetic stem cell maintenance and chemoresistance in colon cancer.

Cancer stem cells drive tumor growth and survival via self-renewal and therapeutic resistance, but the upstream mechanisms are not well defined. In this issue of PLOS Biology, a study in colon cancer reveals a new signalling network that links epigenetic regulation to these phenotypes.

Tissue factor is a critical regulator of radiation therapy-induced glioblastoma remodeling.

Radiation therapy (RT) provides therapeutic benefits for patients with glioblastoma (GBM), but inevitably induces poorly understood global changes in GBM and its microenvironment (TME) that promote radio-resistance and recurrence. Through a cell surface marker screen, we identified that CD142 (tissue factor or F3) is robustly induced in the senescence-associated ß-galactosidase (SA-ßGal)-positive GBM cells after irradiation. F3 promotes clonal expansion of irradiated SA-ßGal+ GBM cells and orchestrates oncogenic TME remodeling by activating both tumor-autonomous signaling and extrinsic coagulation pathways. Intratumoral F3 signaling induces a mesenchymal-like cell state transition and elevated chemokine secretion. Simultaneously, F3-mediated focal hypercoagulation states lead to activation of tumor-associated macrophages (TAMs) and extracellular matrix (ECM) remodeling. A newly developed F3-targeting agent potently inhibits the aforementioned oncogenic events and impedes tumor relapse in vivo. These findings support F3 as a critical regulator for therapeutic resistance and oncogenic senescence in GBM, opening potential therapeutic avenues.

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.

VRK1 Is a Synthetic-Lethal Target in VRK2-Deficient Glioblastoma.

Synthetic lethality is a genetic interaction that results in cell death when two genetic deficiencies co-occur but not when either deficiency occurs alone, which can be co-opted for cancer therapeutics. Pairs of paralog genes are among the most straightforward potential synthetic-lethal interactions by virtue of their redundant functions. Here, we demonstrate a paralog-based synthetic lethality by targeting vaccinia-related kinase 1 (VRK1) in glioblastoma (GBM) deficient of VRK2, which is silenced by promoter methylation in approximately two thirds of GBM. Genetic knockdown of VRK1 in VRK2-null or VRK2-methylated cells resulted in decreased activity of the downstream substrate barrier to autointegration factor (BAF), a regulator of post-mitotic nuclear envelope formation. Reduced BAF activity following VRK1 knockdown caused nuclear lobulation, blebbing, and micronucleation, which subsequently resulted in G2-M arrest and DNA damage. The VRK1-VRK2 synthetic-lethal interaction was dependent on VRK1 kinase activity and was rescued by ectopic expression of VRK2. In VRK2-methylated GBM cell line-derived xenograft and patient-derived xenograft models, knockdown of VRK1 led to robust tumor growth inhibition. These results indicate that inhibiting VRK1 kinase activity could be a viable therapeutic strategy in VRK2-methylated GBM. SIGNIFICANCE: A paralog synthetic-lethal interaction between VRK1 and VRK2 sensitizes VRK2-methylated glioblastoma to perturbation of VRK1 kinase activity, supporting VRK1 as a drug discovery target in this disease.

Maintaining Human Glioblastoma Cellular Diversity Ex vivo using Three-Dimensional Organoid Culture.

Glioblastoma (GBM) is the most commonly occurring primary malignant brain cancer with an extremely poor prognosis. Intra-tumoral cellular and molecular diversity, as well as complex interactions between tumor microenvironments, can make finding effective treatments a challenge. Traditional adherent or sphere culture methods can mask such complexities, whereas three-dimensional organoid culture can recapitulate regional microenvironmental gradients. Organoids are a method of three-dimensional GBM culture that better mimics patient tumor architecture, contains phenotypically diverse cell populations, and can be used for medium-throughput experiments. Although three-dimensional organoid culture is more laborious and time-consuming compared to traditional culture, it offers unique benefits and can serve to bridge the gap between current in vitro and in vivo systems. Organoids have established themselves as invaluable tools in the arsenal of cancer biologists to better understand tumor behavior and mechanisms of resistance, and their applications only continue to grow. Here, details are provided about methods for generating and maintaining GBM organoids. Instructions of how to perform organoid sample embedding and sectioning using both frozen and paraffin-embedding techniques, as well as recommendations for immunohistochemistry and immunofluorescence protocols on organoid sections, and measurement of total organoid cell viability, are all also described.

Patient-derived explants as tumor models.

Tumors contain heterogeneous neoplastic cells and diverse stromal elements that collectively function as dynamic ecosystems, and this complicates predictive modeling ex vivo. LeBlanc et al. utilize single-cell analysis to demonstrate that patient-derived explants replicate tumor cell diversity and transient stromal cell types in patient surgical specimens. This suggests that patient-derived explants can be valuable as tumor models.

Three-dimensional organoid culture unveils resistance to clinical therapies in adult and pediatric glioblastoma.

BACKGROUND: Glioblastoma (GBM) is the most common primary brain tumor with a dismal prognosis. The inherent cellular diversity and interactions within tumor microenvironments represent significant challenges to effective treatment. Traditional culture methods such as adherent or sphere cultures may mask such complexities whereas three-dimensional (3D) organoid culture systems derived from patient cancer stem cells (CSCs) can preserve cellular complexity and microenvironments. The objective of this study was to determine if GBM organoids may offer a platform, complimentary to traditional sphere culture methods, to recapitulate patterns of clinical drug resistance arising from 3D growth. METHODS: Adult and pediatric surgical specimens were collected and established as organoids. We created organoid microarrays and visualized bulk and spatial differences in cell proliferation using immunohistochemistry (IHC) staining, and cell cycle analysis by flow cytometry paired with 3D regional labeling. We tested the response of CSCs grown in each culture method to temozolomide, ibrutinib, lomustine, ruxolitinib, and radiotherapy. RESULTS: GBM organoids showed diverse and spatially distinct proliferative cell niches and include heterogeneous populations of CSCs/non-CSCs (marked by SOX2) and cycling/senescent cells. Organoid cultures display a comparatively blunted response to current standard-of-care therapy (combination temozolomide and radiotherapy) that reflects what is seen in practice. Treatment of organoids with clinically relevant drugs showed general therapeutic resistance with drug- and patient-specific antiproliferative, apoptotic, and senescent effects, differing from those of matched sphere cultures. CONCLUSIONS: Therapeutic resistance in organoids appears to be driven by altered biological mechanisms rather than physical limitations of therapeutic access. GBM organoids may therefore offer a key technological approach to discover and understand resistance mechanisms of human cancer cells.

Altered lipid metabolism marks glioblastoma stem and non-stem cells in separate tumor niches.

Glioblastoma (GBM) displays marked cellular and metabolic heterogeneity that varies among cellular microenvironments within a tumor. Metabolic targeting has long been advocated as a therapy against many tumors including GBM, but how lipid metabolism is altered to suit different microenvironmental conditions and whether cancer stem cells (CSCs) have altered lipid metabolism are outstanding questions in the field. We interrogated gene expression in separate microenvironments of GBM organoid models that mimic the transition between nutrient-rich and nutrient-poor pseudopalisading/perinecrotic tumor zones using spatial-capture RNA-sequencing. We revealed a striking difference in lipid processing gene expression and total lipid content between diverse cell populations from the same patient, with lipid enrichment in hypoxic organoid cores and also in perinecrotic and pseudopalisading regions of primary patient tumors. This was accompanied by regionally restricted upregulation of hypoxia-inducible lipid droplet-associated (HILPDA) gene expression in organoid cores and pseudopalisading regions of clinical GBM specimens, but not lower-grade brain tumors. CSCs have low lipid droplet accumulation compared to non-CSCs in organoid models and xenograft tumors, and prospectively sorted lipid-low GBM cells are functionally enriched for stem cell activity. Targeted lipidomic analysis of multiple patient-derived models revealed a significant shift in lipid metabolism between GBM CSCs and non-CSCs, suggesting that lipid levels may not be simply a product of the microenvironment but also may be a reflection of cellular state. CSCs had decreased levels of major classes of neutral lipids compared to non-CSCs, but had significantly increased polyunsaturated fatty acid production due to high fatty acid desaturase (FADS1/2) expression which was essential to maintain CSC viability and self-renewal. Our data demonstrate spatially and hierarchically distinct lipid metabolism phenotypes occur clinically in the majority of patients, can be recapitulated in laboratory models, and may represent therapeutic targets for GBM.

Glioblastoma Stem Cells: Driving Resilience through Chaos.

Glioblastoma is an aggressive and heterogeneous tumor in which glioblastoma stem cells (GSCs) are at the apex of an entropic hierarchy and impart devastating therapy resistance. The high entropy of GSCs is driven by a permissive epigenetic landscape and a mutational landscape that revokes crucial cellular checkpoints. The GSC population encompasses a complex array of diverse microstates that are defined and maintained by a wide variety of attractors including the complex tumor ecosystem and therapeutic intervention. Constant dynamic transcriptional fluctuations result in a highly adaptable and heterogeneous entity primed for therapy evasion and survival. Analyzing the transcriptional, epigenetic, and metabolic landscapes of GSC dynamics in the context of a stochastically fluctuating tumor network will provide novel strategies to target resistant populations of GSCs in glioblastoma.

Luteolin inhibits Musashi1 binding to RNA and disrupts cancer phenotypes in glioblastoma cells.

RNA binding proteins have emerged as critical oncogenic factors and potential targets in cancer therapy. In this study, we evaluated Musashi1 (Msi1) targeting as a strategy to treat glioblastoma (GBM); the most aggressive brain tumor type. Msi1 expression levels are often high in GBMs and other tumor types and correlate with poor clinical outcome. Moreover, Msi1 has been implicated in chemo- and radio-resistance. Msi1 modulates a range of cancer relevant processes and pathways and regulates the expression of stem cell markers and oncogenic factors via mRNA translation/stability. To identify Msi1 inhibitors capable of blocking its RNA binding function, we performed a ~ 25,000 compound fluorescence polarization screen. NMR and LSPR were used to confirm direct interaction between Msi1 and luteolin, the leading compound. Luteolin displayed strong interaction with Msi1 RNA binding domain 1 (RBD1). As a likely consequence of this interaction, we observed via western and luciferase assays that luteolin treatment diminished Msi1 positive impact on the expression of pro-oncogenic target genes. We tested the effect of luteolin treatment on GBM cells and showed that it reduced proliferation, cell viability, colony formation, migration and invasion of U251 and U343 GBM cells. Luteolin also decreased the proliferation of patient-derived glioma initiating cells (GICs) and tumor-organoids but did not affect normal astrocytes. Finally, we demonstrated the value of combined treatments with luteolin and olaparib (PARP inhibitor) or ionizing radiation (IR). Our results show that luteolin functions as an inhibitor of Msi1 and demonstrates its potential use in GBM therapy.

Reciprocal Signaling between Glioblastoma Stem Cells and Differentiated Tumor Cells Promotes Malignant Progression.

Glioblastoma is the most lethal primary brain tumor; however, the crosstalk between glioblastoma stem cells (GSCs) and their supportive niche is not well understood. Here, we interrogated reciprocal signaling between GSCs and their differentiated glioblastoma cell (DGC) progeny. We found that DGCs accelerated GSC tumor growth. DGCs preferentially expressed brain-derived neurotrophic factor (BDNF), whereas GSCs expressed the BDNF receptor NTRK2. Forced BDNF expression in DGCs augmented GSC tumor growth. To determine molecular mediators of BDNF-NTRK2 paracrine signaling, we leveraged transcriptional and epigenetic profiles of matched GSCs and DGCs, revealing preferential VGF expression by GSCs, which patient-derived tumor models confirmed. VGF serves a dual role in the glioblastoma hierarchy by promoting GSC survival and stemness in vitro and in vivo while also supporting DGC survival and inducing DGC secretion of BDNF. Collectively, these data demonstrate that differentiated glioblastoma cells cooperate with stem-like tumor cells through BDNF-NTRK2-VGF paracrine signaling to promote tumor growth.

Therapeutic targeting of ependymoma as informed by oncogenic enhancer profiling.

Genomic sequencing has driven precision-based oncology therapy; however, the genetic drivers of many malignancies remain unknown or non-targetable, so alternative approaches to the identification of therapeutic leads are necessary. Ependymomas are chemotherapy-resistant brain tumours, which, despite genomic sequencing, lack effective molecular targets. Intracranial ependymomas are segregated on the basis of anatomical location (supratentorial region or posterior fossa) and further divided into distinct molecular subgroups that reflect differences in the age of onset, gender predominance and response to therapy. The most common and aggressive subgroup, posterior fossa ependymoma group A (PF-EPN-A), occurs in young children and appears to lack recurrent somatic mutations. Conversely, posterior fossa ependymoma group B (PF-EPN-B) tumours display frequent large-scale copy number gains and losses but have favourable clinical outcomes. More than 70% of supratentorial ependymomas are defined by highly recurrent gene fusions in the NF-κB subunit gene RELA (ST-EPN-RELA), and a smaller number involve fusion of the gene encoding the transcriptional activator YAP1 (ST-EPN-YAP1). Subependymomas, a distinct histologic variant, can also be found within the supratetorial and posterior fossa compartments, and account for the majority of tumours in the molecular subgroups ST-EPN-SE and PF-EPN-SE. Here we describe mapping of active chromatin landscapes in 42 primary ependymomas in two non-overlapping primary ependymoma cohorts, with the goal of identifying essential super-enhancer-associated genes on which tumour cells depend. Enhancer regions revealed putative oncogenes, molecular targets and pathways; inhibition of these targets with small molecule inhibitors or short hairpin RNA diminished the proliferation of patient-derived neurospheres and increased survival in mouse models of ependymomas. Through profiling of transcriptional enhancers, our study provides a framework for target and drug discovery in other cancers that lack known genetic drivers and are therefore difficult to treat.

Zika virus has oncolytic activity against glioblastoma stem cells.

Glioblastoma is a highly lethal brain cancer that frequently recurs in proximity to the original resection cavity. We explored the use of oncolytic virus therapy against glioblastoma with Zika virus (ZIKV), a flavivirus that induces cell death and differentiation of neural precursor cells in the developing fetus. ZIKV preferentially infected and killed glioblastoma stem cells (GSCs) relative to differentiated tumor progeny or normal neuronal cells. The effects against GSCs were not a general property of neurotropic flaviviruses, as West Nile virus indiscriminately killed both tumor and normal neural cells. ZIKV potently depleted patient-derived GSCs grown in culture and in organoids. Moreover, mice with glioblastoma survived substantially longer and at greater rates when the tumor was inoculated with a mouse-adapted strain of ZIKV. Our results suggest that ZIKV is an oncolytic virus that can preferentially target GSCs; thus, genetically modified strains that further optimize safety could have therapeutic efficacy for adult glioblastoma patients.

MYC-Regulated Mevalonate Metabolism Maintains Brain Tumor-Initiating Cells.

Metabolic dysregulation drives tumor initiation in a subset of glioblastomas harboring isocitrate dehydrogenase (IDH) mutations, but metabolic alterations in glioblastomas with wild-type IDH are poorly understood. MYC promotes metabolic reprogramming in cancer, but targeting MYC has proven notoriously challenging. Here, we link metabolic dysregulation in patient-derived brain tumor-initiating cells (BTIC) to a nexus between MYC and mevalonate signaling, which can be inhibited by statin or 6-fluoromevalonate treatment. BTICs preferentially express mevalonate pathway enzymes, which we find regulated by novel MYC-binding sites, validating an additional transcriptional activation role of MYC in cancer metabolism. Targeting mevalonate activity attenuated RAS-ERK-dependent BTIC growth and self-renewal. In turn, mevalonate created a positive feed-forward loop to activate MYC signaling via induction of miR-33b. Collectively, our results argue that MYC mediates its oncogenic effects in part by altering mevalonate metabolism in glioma cells, suggesting a therapeutic strategy in this setting. Cancer Res; 77(18); 4947-60. ©2017 AACR.