Visualizing molecular profiles of glioblastoma with GBM-BioDP

3D Brain Vascular Niche Model Captures Invasive Behavior and Gene Signatures of Glioblastoma

Glioblastoma (GBM) is a lethal brain cancer with no effective treatment; understanding how GBM cells respond to tumor microenvironment remains challenging as conventional cell cultures lack proper cytoarchitecture while in vivo animal models present complexity all at once. Developing a culture system to bridge the gap is thus crucial. Here, we employed a multicellular approach using human glia and vascular cells to optimize a 3-dimensional (3D) brain vascular niche model that enabled not only long-term culture of patient derived GBM cells but also recapitulation of key features of GBM heterogeneity, in particular invasion behavior and vascular association. Comparative transcriptomics of identical patient derived GBM cells in 3D and in vivo xenotransplants models revealed that glia-vascular contact induced genes concerning neural/glia development, synaptic regulation, as well as immune suppression. This gene signature displayed region specific enrichment in the leading edge and microvascular proliferation zones in human GBM and predicted poor prognosis. Gene variance analysis also uncovered histone demethylation and xylosyltransferase activity as main themes for gene adaption of GBM cells in vivo . Furthermore, our 3D model also demonstrated the capacity to provide a quiescence and a protective niche against chemotherapy. In summary, an advanced 3D brain vascular model can bridge the gap between 2D cultures and in vivo models in capturing key features of GBM heterogeneity and unveil previously unrecognized influence of glia-vascular contact for transcriptional adaption in GBM cells featuring neural/synaptic interaction and immunosuppression.

C-Myc/H19/miR-29b axis downregulates nerve/glial (NG)2 expression in glioblastoma multiforme

Nerve/glial antigen (NG)2 is highly expressed in glioblastoma multiforme (GBM). However, the underlying mechanisms of its upregulated expression are largely unknown. In silico analyses reveal that the tumor-suppressive miR-29b targets NG2. We used GBM-based data from The Cancer Genome Atals databases to analyze the expression pattern of miR-29b and different target genes, including NG2. Moreover, we investigated the regulatory function of miR-29b on NG2 expression and NG2-related signaling pathways. We further studied upstream mechanisms affecting miR-29b-dependent NG2 expression. We found that miR-29b downregulates NG2 expression directly and indirectly via the transcription factor Sp1. Furthermore, we identified the NG2 coreceptor platelet-derived growth factor receptor (PDGFR)α as an additional miR-29b target. As shown by a panel of functional cell assays, a reduced miR-29b-dependent NG2 expression suppresses tumor cell proliferation and migration. Signaling pathway analyses revealed that this is associated with a decreased ERK1/2 activity. In addition, we found that the long noncoding RNA H19 and c-Myc act as upstream repressors of miR-29b in GBM cells, resulting in an increased NG2 expression. These findings indicate that the c-Myc/H19/miR-29b axis crucially regulates NG2 expression in GBM and, thus, represents a target for the development of future GBM therapies.

CA9, CYFIP2 and LGALS3BP-A Novel Biomarker Panel to Aid Prognostication in Glioma

Brain cancer is a devastating and life-changing disease. Biomarkers are becoming increasingly important in addressing clinical issues, including in monitoring tumour progression and assessing survival and treatment response. The goal of this study was to identify prognostic biomarkers associated with glioma progression. Discovery proteomic analysis was performed on a small cohort of astrocytomas that were diagnosed as low-grade and recurred at a higher grade. Six proteins were chosen to be validated further in a larger cohort. Three proteins, CA9, CYFIP2, and LGALS3BP, were found to be associated with glioma progression and, in univariate analysis, could be used as prognostic markers. However, according to the results of multivariate analysis, these did not remain significant. These three proteins were then combined into a three-protein panel. This panel had a specificity and sensitivity of 0.7459 for distinguishing between long and short survival. In silico data confirmed the prognostic significance of this panel.

Hypoxic niches attract and sequester tumor-associated macrophages and cytotoxic T cells and reprogram them for immunosuppression

Glioblastoma (GBM), a highly lethal brain cancer, is notorious for immunosuppression, but the mechanisms remain unclear. Here, we documented a temporospatial patterning of tumor-associated myeloid cells (TAMs) corresponding to vascular changes during GBM progression. As tumor vessels transitioned from the initial dense regular network to later scant and engorged vasculature, TAMs shifted away from perivascular regions and trafficked to vascular-poor areas. This process was heavily influenced by the immunocompetence state of the host. Utilizing a sensitive fluorescent UnaG reporter to track tumor hypoxia, coupled with single-cell transcriptomics, we revealed that hypoxic niches attracted and sequestered TAMs and cytotoxic T lymphocytes (CTLs), where they were reprogrammed toward an immunosuppressive state. Mechanistically, we identified chemokine CCL8 and cytokine IL-1β as two hypoxic-niche factors critical for TAM trafficking and co-evolution of hypoxic zones into pseudopalisading patterns. Therefore, perturbation of TAM patterning in hypoxic zones may improve tumor control.

Glioma-BioDP: database for visualization of molecular profiles to improve prognosis of brain cancer

Cancer researchers often seek user-friendly interactive tools for validation, exploration, analysis, and visualization of molecular profiles in cancer patient samples. To aid researchers working on the both low- and high-grade gliomas, we developed Glioma-BioDP, a web tool for exploration and visualization of RNA and protein expression profiles of interest in these tumor types. Glioma-BioDP is user friendly application that include expression data from both the low- and high-grade glioma patient samples from The Cancer Genome Atlas and enabled querying by mRNA, microRNA, and protein level expression data from Illumina HiSeq and RPPA platforms respectively. Glioma-BioDP provides advance query interface and enables users to explore the association of genes, proteins, and miRNA expression with molecular and/or histological subtypes of gliomas, surgical resection status and survival. The prognostic significance and visualization of the selected expression profiles can be explored using interactive utilities provided. This tool may also enable validation and generation of new hypotheses of novel therapies impacting gliomas that aid in personalization of treatment for optimum outcomes.

Glioblastoma cells have increased capacity to repair radiation-induced DNA damage after migration to the olfactory bulb

BACKGROUND

The invasive nature of GBM combined with the diversity of brain microenvironments creates the potential for a topographic heterogeneity in GBM radioresponse. Investigating the mechanisms responsible for a microenvironment-induced differential GBM response to radiation may provide insights into the molecules and processes mediating GBM radioresistance.

METHODS

Using a model system in which human GBM stem-like cells implanted into the right striatum of nude mice migrate throughout the right hemisphere (RH) to the olfactory bulb (OB), the radiation-induced DNA damage response was evaluated in each location according to γH2AX and 53BP1 foci and cell cycle phase distribution as determined by flow cytometry and immunohistochemistry. RNAseq was used to compare transcriptomes of tumor cells growing in the OB and the RH. Protein expression and neuron-tumor interaction were defined by immunohistochemistry and confocal microscopy.

RESULTS

After irradiation, there was a more rapid dispersal of γH2AX and 53BP1 foci in the OB versus in the RH, indicative of increased double strand break repair capacity in the OB and consistent with the OB providing a radioprotective niche. With respect to the cell cycle, by 6 h after irradiation there was a significant loss of mitotic tumor cells in both locations suggesting a similar activation of the G2/M checkpoint. However, by 24 h post-irradiation there was an accumulation of G2 phase cells in the OB, which continued out to at least 96 h. Transcriptome analysis showed that tumor cells in the OB had higher expression levels of DNA repair genes involved in non-homologous end joining and genes related to the spindle assembly checkpoint. Tumor cells in the OB were also found to have an increased frequency of soma-soma contact with neurons.

CONCLUSION

GBM cells that have migrated to the OB have an increased capacity to repair radiation-induced double strand breaks and altered cell cycle regulation. These results correspond to an upregulation of genes involved in DNA damage repair and cell cycle control. Because the murine OB provides a source of radioresistant tumor cells not evident in other experimental systems, it may serve as a model for investigating the mechanisms mediating GBM radioresistance.

circSMARCA5 Is an Upstream Regulator of the Expression of miR-126-3p, miR-515-5p, and Their mRNA Targets, Insulin-like Growth Factor Binding Protein 2 (IGFBP2) and NRAS Proto-Oncogene, GTPase (NRAS) in Glioblastoma

The involvement of non-coding RNAs (ncRNAs) in glioblastoma multiforme (GBM) pathogenesis and progression has been ascertained but their cross-talk within GBM cells remains elusive. We previously demonstrated the role of circSMARCA5 as a tumor suppressor (TS) in GBM. In this paper, we explore the involvement of circSMARCA5 in the control of microRNA (miRNA) expression in GBM. By using TaqMan^® low-density arrays, the expression of 748 miRNAs was assayed in U87MG overexpressing circSMARCA5. Differentially expressed (DE) miRNAs were validated through single TaqMan^® assays in: (i) U87MG overexpressing circSMARCA5; (ii) four additional GBM cell lines (A172; CAS-1; SNB-19; U251MG); (iii) thirty-eight GBM biopsies; (iv) twenty biopsies of unaffected brain parenchyma (UC). Validated targets of DE miRNAs were selected from the databases TarBase and miRTarbase, and the literature; their expression was inferred from the GBM TCGA dataset. Expression was assayed in U87MG overexpressing circSMARCA5, GBM cell lines, and biopsies through real-time PCR. TS miRNAs 126-3p and 515-5p were upregulated following circSMARCA5 overexpression in U87MG and their expression was positively correlated with that of circSMARCA5 (r-values = 0.49 and 0.50, p-values = 9 × 10^-5 and 7 × 10^-5, respectively) in GBM biopsies. Among targets, IGFBP2 (target of miR-126-3p) and NRAS (target of miR-515-5p) mRNAs were positively correlated (r-value = 0.46, p-value = 0.00027), while their expression was negatively correlated with that of circSMARCA5 (r-values = -0.58 and -0.30, p-values = 0 and 0.019, respectively), miR-126-3p (r-value = -0.36, p-value = 0.0066), and miR-515-5p (r-value = -0.34, p-value = 0.010), respectively. Our data identified a new GBM subnetwork controlled by circSMARCA5, which regulates downstream miRNAs 126-3p and 515-5p, and their mRNA targets IGFBP2 and NRAS.

MicroRNA and mRNA Expression Changes in Glioblastoma Cells Cultivated under Conditions of Neurosphere Formation

Glioblastoma multiforme (GBM) is one of the most highly metastatic cancers. The study of the pathogenesis of GBM, as well as the development of targeted oncolytic drugs, require the use of actual cell models, in particular, the use of 3D cultures or neurospheres (NS). During the formation of NS, the adaptive molecular landscape of the transcriptome, which includes various regulatory RNAs, changes. The aim of this study was to reveal changes in the expression of microRNAs (miRNAs) and their target mRNAs in GBM cells under conditions of NS formation. Neurospheres were obtained from both immortalized U87 MG and patient-derived BR3 GBM cell cultures. Next generation sequencing analysis of small and long RNAs of adherent and NS cultures of GBM cells was carried out. It was found that the formation of NS proceeds with an increase in the level of seven and a decrease in the level of 11 miRNAs common to U87 MG and BR3, as well as an increase in the level of 38 and a decrease in the level of 12 mRNA/lncRNA. Upregulation of miRNAs hsa-miR: -139-5p; -148a-3p; -192-5p; -218-5p; -34a-5p; and -381-3p are accompanied by decreased levels of their target mRNAs: RTN4, FLNA, SH3BP4, DNPEP, ETS2, MICALL1, and GREM1. Downregulation of hsa-miR: -130b-5p, -25-5p, -335-3p and -339-5p occurs with increased levels of mRNA-targets BDKRB2, SPRY4, ERRFI1 and TGM2. The involvement of SPRY4, ERRFI1, and MICALL1 mRNAs in the regulation of EGFR/FGFR signaling highlights the role of hsa-miR: -130b-5p, -25-5p, -335-3p, and -34a-5p not only in the formation of NS, but also in the regulation of malignant growth and invasion of GBM. Our data provide the basis for the development of new approaches to the diagnosis and treatment of GBM.

Extracellular Vesicles Secreted by Glioma Stem Cells Are Involved in Radiation Resistance and Glioma Progression

Glioblastoma is the most aggressive brain tumour with short survival, partly due to resistance to conventional therapy. Glioma stem cells (GSC) are likely to be involved in treatment resistance, by releasing extracellular vesicles (EVs) containing specific molecular cargoes. Here, we studied the EVs secreted by glioma stem cells (GSC-EVs) and their effects on radiation resistance and glioma progression. EVs were isolated from 3 GSCs by serial centrifugation. NanoSight measurement, cryo-electron microscopy and live imaging were used to study the EVs size, morphology and uptake, respectively. The non-GSC glioma cell lines LN229 and U118 were utilised as a recipient cell model. Wound healing assays were performed to detect cell migration. Colony formation, cell viability and invadopodium assays were conducted to detect cell survival of irradiated recipient cells and cell invasion post GSC-EV treatment. NanoString miRNA global profiling was used to select for the GSC-EVs’ specific miRNAs. All three GSC cell lines secreted different amounts of EVs, and all expressed consistent levels of CD9 but different level of Alix, TSG101 and CD81. EVs were taken up by both LN229 and U118 recipient cells. In the presence of GSC-EVs, these recipient cells survived radiation exposure and initiated colony formation. After GSC-EVs exposure, LN229 and U118 cells exhibited an invasive phenotype, as indicated by an increase in cell migration. We also identified 25 highly expressed miRNAs in the GSC-EVs examined, and 8 of these miRNAs can target PTEN. It is likely that GSC-EVs and their specific miRNAs induced the phenotypic changes in the recipient cells due to the activation of the PTEN/Akt pathway. This study demonstrated that GSC-EVs have the potential to induce radiation resistance and modulate the tumour microenvironment to promote glioma progression. Future therapeutic studies should be designed to interfere with these GSC-EVs and their specific miRNAs.

GFAP Alternative Splicing and the Relevance for Disease – A Focus on Diffuse Gliomas

Glial fibrillary acidic protein (GFAP) is an intermediate filament protein that is characteristic for astrocytes and neural stem cells, and their malignant analogues in glioma. Since the discovery of the protein 50 years ago, multiple alternative splice variants of the GFAP gene have been discovered, leading to different GFAP isoforms. In this review, we will describe GFAP isoform expression from gene to protein to network, taking the canonical isoforms GFAPα and the main alternative variant GFAPδ as the starting point. We will discuss the relevance of studying GFAP and its isoforms in disease, with a specific focus on diffuse gliomas.

Glioblastoma multiforme: a multi-omics analysis of driver genes and tumour heterogeneity

Glioblastoma (GBM) is the most aggressive and common brain cancer in adults with the lowest life expectancy. The current neuro-oncology practice has incorporated genes involved in key molecular events that drive GBM tumorigenesis as biomarkers to guide diagnosis and design treatment. This study summarizes findings describing the significant heterogeneity of GBM at the transcriptional and genomic levels, emphasizing 18 driver genes with clinical relevance. A pattern was identified fitting the stem cell model for GBM ontogenesis, with an upregulation profile for MGMT and downregulation for ATRX, H3F3A, TP53 and EGFR in the mesenchymal subtype. We also detected overexpression of EGFR, NES, VIM and TP53 in the classical subtype and of MKi67 and OLIG2 genes in the proneural subtype. Furthermore, we found a combination of the four biomarkers EGFR, NES, OLIG2 and VIM with a remarkable differential expression pattern which confers them a strong potential to determine the GBM molecular subtype. A unique distribution of somatic mutations was found for the young and adult population, particularly for genes related to DNA repair and chromatin remodelling, highlighting ATRX, MGMT and IDH1. Our results also revealed that highly lesioned genes undergo differential regulation with particular biological pathways for young patients. This multi-omic analysis will help delineate future strategies related to the use of these molecular markers for clinical decision-making in the medical routine.

MicroRNA Expression Profile Distinguishes Glioblastoma Stem Cells from Differentiated Tumor Cells

Glioblastoma (GBM) represents the most common and aggressive tumor of the brain. Despite the fact that several studies have recently addressed the molecular mechanisms underlying the disease, its etiology and pathogenesis are still poorly understood. GBM displays poor prognosis and its resistance to common therapeutic approaches makes it a highly recurrent tumor. Several studies have identified a subpopulation of tumor cells, known as GBM cancer stem cells (CSCs) characterized by the ability of self-renewal, tumor initiation and propagation. GBM CSCs have been shown to survive GBM chemotherapy and radiotherapy. Thus, targeting CSCs represents a promising approach to treat GBM. Recent evidence has shown that GBM is characterized by a dysregulated expression of microRNA (miRNAs). In this study we have investigated the difference between human GBM CSCs and their paired autologous differentiated tumor cells. Array-based profiling and quantitative Real-Time PCR (qRT-PCR) were performed to identify miRNAs differentially expressed in CSCs. The Cancer Genome Atlas (TCGA) data were also interrogated, and functional interpretation analysis was performed. We have identified 14 miRNAs significantly differentially expressed in GBM CSCs (p < 0.005). MiR-21 and miR-95 were among the most significantly deregulated miRNAs, and their expression was also associated to patient survival. We believe that the data provided here carry important implications for future studies aiming at elucidating the molecular mechanisms underlying GBM.

Oncolytic Virus Therapy Alters the Secretome of Targeted Glioblastoma Cells

Simple Summary

Proteins secreted by cancer cells in response to oncolytic virus anti-tumor therapy constitute the instructions for the immune cells. Yet as there are hundreds of proteins, including those encapsulated in vesicles, whose message drives the mobilization of immune cells, we aimed to decipher the instruction sent by cancer cells in response to therapy. Searching the cataloged vesicle and vesicle-free secreted proteins, we found that the proteins associated with the favorable survival of brain cancer patients were those that have the power to mobilize the immune cells. Thus, this approach established cancer-secreted contributors to the immune–therapeutic effect of the oncolytic virus.

Abstract

Oncolytic virus (OV) therapy, which is being tested in clinical trials for glioblastoma, targets cancer cells, while triggering immune cells. Yet OV sensitivity varies from patient to patient. As OV therapy is regarded as an anti-tumor vaccine, by making OV-infected cancer cells secrete immunogenic proteins, linking these proteins to transcriptome would provide a measuring tool to predict their sensitivity. A set of six patient-derived glioblastoma cells treated ex-vivo with herpes simplex virus type 1 (HSV1) modeled a clinical setting of OV infection. The cellular transcriptome and secreted proteome (separated into extracellular vesicles (EV) and EV-depleted fractions) were analyzed by gene microarray and mass-spectroscopy, respectively. Data validation and in silico analysis measured and correlated the secretome content with the response to infection and patient survival. Glioblastoma cells reacted to the OV infection in a seemingly dissimilar fashion, but their transcriptomes changed in the same direction. Therefore, the upregulation of transcripts encoding for secreted proteins implies a common thread in the response of cancer cells to infection. Indeed, the OV-driven secretome is linked to the immune response. While these proteins have distinct membership in either EV or EV-depleted fractions, it is their co-secretion that augments the immune response and associates with favorable patient outcomes.

The Role of Network Science in Glioblastoma

Network science has long been recognized as a well-established discipline across many biological domains. In the particular case of cancer genomics, network discovery is challenged by the multitude of available high-dimensional heterogeneous views of data. Glioblastoma (GBM) is an example of such a complex and heterogeneous disease that can be tackled by network science. Identifying the architecture of molecular GBM networks is essential to understanding the information flow and better informing drug development and pre-clinical studies. Here, we review network-based strategies that have been used in the study of GBM, along with the available software implementations for reproducibility and further testing on newly coming datasets. Promising results have been obtained from both bulk and single-cell GBM data, placing network discovery at the forefront of developing a molecularly-informed-based personalized medicine.

Combined PI3Kα-mTOR Targeting of Glioma Stem Cells

Glioblastoma (GBM) is the most common and lethal primary intrinsic tumour of the adult brain and evidence indicates disease progression is driven by glioma stem cells (GSCs). Extensive advances in the molecular characterization of GBM allowed classification into proneural, mesenchymal and classical subtypes, and have raised expectations these insights may predict response to targeted therapies. We utilized GBM neurospheres that display GSC characteristics and found activation of the PI3K/AKT pathway in sphere-forming cells. The PI3Kα selective inhibitor alpelisib blocked PI3K/AKT activation and inhibited spheroid growth, suggesting an essential role for the PI3Kα catalytic isoform. p110α expression was highest in the proneural subtype and this was associated with increased phosphorylation of AKT. Further, employing the GBM BioDP, we found co-expression of PIK3CA with the neuronal stem/progenitor marker NES was associated with poor prognosis in PN GBM patients, indicating a unique role for PI3Kα in PN GSCs. Alpelisib inhibited GSC neurosphere growth and these effects were more pronounced in GSCs of the PN subtype. The antineoplastic effects of alpelisib were substantially enhanced when combined with pharmacologic mTOR inhibition. These findings identify the alpha catalytic PI3K isoform as a unique therapeutic target in proneural GBM and suggest that pharmacological mTOR inhibition may sensitize GSCs to selective PI3Kα inhibition.

The nuclear DICER-circular RNA complex drives the deregulation of the glioblastoma cell microRNAome

The assortment of cellular microRNAs ("microRNAome") is a vital readout of cellular homeostasis, but the mechanisms that regulate the microRNAome are poorly understood. The microRNAome of glioblastoma is substantially down-regulated in comparison to the normal brain. Here, we find malfunction of the posttranscriptional maturation of the glioblastoma microRNAome and link it to aberrant nuclear localization of DICER, the major enzymatic complex responsible for microRNA maturation. Analysis of DICER's nuclear interactome reveals the presence of an RNA binding protein, RBM3, and of a circular RNA, circ2082, within the complex. Targeting of this complex by knockdown of circ2082 results in the restoration of cytosolic localization of DICER and widespread derepression of the microRNAome, leading to transcriptome-wide rearrangements that mitigate the tumorigenicity of glioblastoma cells in vitro and in vivo with correlation to favorable outcomes in patients with glioblastoma. These findings uncover the mechanistic foundation of microRNAome deregulation in malignant cells.

Present and Future of Anti-Glioblastoma Therapies: A Deep Look into Molecular Dependencies/Features

Glioblastoma (GBM) is aggressive malignant tumor residing within the central nervous system. Although the standard treatment options, consisting of surgical resection followed by combined radiochemotherapy, have long been established for patients with GBM, the prognosis is still poor. Despite recent advances in diagnosis, surgical techniques, and therapeutic approaches, the increased patient survival after such interventions is still sub-optimal. The unique characteristics of GBM, including highly infiltrative nature, hard-to-access location (mainly due to the existence of the blood brain barrier), frequent and rapid recurrence, and multiple drug resistance mechanisms, pose challenges to the development of an effective treatment. To overcome current limitations on GBM therapy and devise ideal therapeutic strategies, efforts should focus on an improved molecular understanding of GBM pathogenesis. In this review, we summarize the molecular basis for the development and progression of GBM as well as some emerging therapeutic approaches.

microRNA-10 and -221 modulate differential expression of Hippo signaling pathway in human astroglial tumors

Gliomas represent over 70% of all brain tumors, they are highly invasive and structurally vascular neoplasms. Despite the latest technological advance in neuro-surgery the survival of patients with high-grade glioma remains poor. The lack of robust treatment options has propelled the search for new markers that may able allow the identification of patients who can benefit from molecularly targeted therapies. The Hippo signaling pathway is considered as a key regulator of tissue homeostasis, cell proliferation and apoptosis, and alterations of this pathway seem to contribute to tumorigenesis. Yes-associated protein (YAP1) is a downstream target of the Hippo pathway which acts as a transcription co-activator. In cancer, YAP1 has been reported to function either as an oncogene or tumor suppressor, depending on the cell context. The aim of this study was to examine the expression of YAP1, Survivin and LATS1 kinase activity in human astroglial tumors with different grades of malignancy. Moreover, we also investigated the expression of miR-221 and miR-10b and their relationship with core molecules of the Hippo pathway. Our results showed the overexpression of YAP1 and Survivin as well as a decreased activity of large tumor suppressor 1 (LATS1) in high-grade glioblastoma versus anaplastic astrocytoma and low-grade glioma. Furthermore, we also demonstrated that miR-221 and miR-10b are specifically involved in Hippo signaling via LATS1 regulation and that their knockdown significantly decreased glioma cell proliferation. This preliminary data confirmed the crucial role of the Hippo pathway in cancer and suggested that miR-221 and miR-10b could be potential therapeutic targets for glioma treatment.

Transfer of MicroRNA via Macrophage-Derived Extracellular Vesicles Promotes Proneural-to-Mesenchymal Transition in Glioma Stem Cells

Proneural-to-mesenchymal transition (PMT) is a common process in glioblastoma (GBM) progression that leads to increased radiotherapy resistance. However, the mechanism underlying PMT is poorly understood. Here, we found that tumor-associated macrophages triggered PMT in glioma stem cells (GSC) via small extracellular vesicles (sEV). sEVs from monocyte-derived macrophages transferred miR-27a-3p, miR-22-3p, and miR-221-3p to GSCs, and these miRNAs promoted several mesenchymal phenotypes in proneural (PN) GSCs by simultaneously targeting CHD7 We found that CHD7 played a critical role in the maintenance of the PN phenotype, and CHD7 knockdown significantly promoted PMT in GSCs via the RelB/P50 and p-STAT3 pathways. The induction of PMT by sEVs containing miR-27a-3p, miR-22-3p, and miR-221-3p in a xenograft nude mouse model exacerbated radiotherapy resistance and thus decreased the benefits of radiotherapy. Collectively, these findings identified macrophage-derived sEVs as key regulators of PMT in GSCs and demonstrated that CHD7 is a novel inhibitor of PMT.

A Serum Proteomic Signature Predicting Survival in Patients with Glioblastoma

Purpose

Glioblastoma (GBM) is the most common form of brain tumor and has a uniformly poor prognosis. Development of prognostic biomarkers in easily accessible serum samples have the potential to improve the outcomes of patients with GBM through personalized therapy planning.

Material/Methods

In this study pre-treatment serum samples from 30 patients newly diagnosed with GBM were evaluated using a 40-protein multiplex ELISA platform. Analysis of potentially relevant gene targets using The Cancer Genome Atlas database was done using the Glioblastoma Bio Discovery Portal (GBM-BioDP). A ten-biomarker subgroup of clinically relevant molecules was selected using a functional grouping analysis of the 40 plex genes with two genes selected from each group on the basis of degree of variance, lack of co-linearity with other biomarkers and clinical interest. A Multivariate Cox proportional hazard approach was used to analyze the relationship between overall survival (OS), gene expression, and resection status as covariates.

Results

Thirty of 40 of the MSD molecules mapped to known genes within TCGA and separated the patient cohort into two main clusters centered predominantly around a grouping of classical and proneural versus the mesenchymal subtype as classified by Verhaak. Using the values for the 30 proteins in a prognostic index (PI) demonstrated that patients in the entire cohort with a PI below the median lived longer than those patients with a PI above the median (HR 1.8, p=0.001) even when stratified by both age and MGMT status. This finding was also consistent within each Verhaak subclass and highly significant (range p=0.0001-0.011). Additionally, a subset of ten proteins including, CRP, SAA, VCAM1, VEGF, MDC, TNFA, IL7, IL8, IL10, IL16 were found to have prognostic value within the TCGA database and a positive correlation with overall survival in GBM patients who had received gross tumor resection followed by conventional radiation therapy and temozolomide treatment concurrent with the addition of valproic acid.

Conclusion

These findings demonstrate that proteomic approaches to the development of prognostic assays for treatment of GBM may hold potential clinical value.

MUDENG Expression Profiling in Cohorts and Brain Tumor Biospecimens to Evaluate Its Role in Cancer

Mu-2-related death-inducing gene (MUDENG, MuD) has been reported to be involved in the tumor necrosis factor-related apoptosis-inducing ligand (TRAIL)-associated apoptotic pathway of glioblastoma multiforme (GBM) cells; however, its expression level, interactors, and role in tumors are yet to be discovered. To investigate whether MuD expression correlates with cancer progression, we analyzed The Cancer Genome Atlas (TCGA) database using UALCAN and Gene Expression Profiling Interactive Analysis (GEPIA). Differential expression of MuD was detected in 6 and 10 cancer types, respectively. Validation performed using data from the Gene Expression Omnibus database showed that MuD expression is downregulated in KIRC tumor and correlate with higher chance of survival. Upregulation of MuD expression in GBM tumors was detected through GEPIA and high MuD expression correlated with higher survival in proneural GBM, whereas the opposite was observed in classical GBM subtype. GBM biospecimens analysis shows that MuD protein level was upregulated in three of six specimens, whereas mRNA level remained relatively unaltered. Therefore, MuD may exert differential effects according to subtypes, and/or be subjected to post-translational regulation in GBM. Correlation analysis between GBM cohort database and experiments using GBM cell lines revealed its positive effect on regulation of protein phosphatase 2 regulatory subunit B’Epsilon (PPP2R5E) and son of sevenless homolog 2 (SOS2). STRING database analysis indicated that the components of adaptor protein complexes putatively interacted with MuD but showed no correlation in terms of survival of patients with different GBM subtypes. In summary, we analyzed the expression of MuD in publicly available cancer patient data sets, GBM cell lines, and biospecimens to demonstrate its potential role as a biomarker for cancer prognosis and identified its candidate interacting molecules.

Identification of miRNA signature associated with BMP2 and chemosensitivity of TMZ in glioblastoma stem-like cells

Glioblastoma multiform (GBM) is the most lethal intracranial tumor in adults. Glioblastoma stem-like cells (GSCs) are responsible for tumorigenesis and chemotherapy resistance. BMPs are known to increase temozolomide (TMZ) response in GSCs, however, the intracellular molecular mechanism remains largely unknown. In this study, we built a GSC cell model called U87S, and performed RNA sequencing to identify differentially expressed (DE) miRNA profiles in U87S cells treated with BMP2, TMZ or combined BMP2 and TMZ respectively. Bioinformatics analysis revealed that most DE miRNAs were involved in the cancer pathways, suggesting their crucial roles in gliomagenesis. Eight miRNAs from RNA-seq were validated. Four out of these miRNAs (has-miR-199a-3p, hsa-miR-374b-5p, hsa-miR-320d, and hsa-miR-339-5p) were found significantly up-regulated in GBM tumor tissues. One of them, hsa-miR-199a-3p, was significantly correlated with the survival of GBM patients, and differentially expressed in U87S cells. Expression of hsa-miR-199a-3p was up-regulated by BMP. Overexpression of hsa-miR-199a-3p in U87S cells inhibited cell viability and enhanced the cytotoxicity of TMZ. And activation of BMP boosted the effect of hsa-miR-199a-3p on cell viability and TMZ-mediated cytotoxicity. Besides, expressions of five predicted targets of hsa-miR-199a-3p were evaluated. Four of them were differentially expressed in GBM tumors. And one of them, SLC22A18, was associated with the survival of GBM patients. In the end, a hsa-miR-199a-3p-mediated ceRNA network was constructed for the convenience of future study. Together, our data provided DE miRNA expression profiles associated with BMP2 and TMZ in GSCs, which might lead to finding out miRNA-based target therapies that specially target GSCs.

GLUT1 and TUBB4 in Glioblastoma Could be Efficacious Targets

Glioblastoma multiforme (GBM) is the most aggressive and deadly brain tumor, portending a median 13-month survival even following gross total resection with adjuvant chemotherapy and radiotherapy. This prognosis necessitates improved therapies for the disease. A target of interest for novel chemotherapies is the Warburg Effect, which describes the tumor's shift away from oxidative phosphorylation towards glycolysis. Here, we elucidate GLUT1 (Glucose transporter 1) and one of its associated binding partners, TUBB4 (Tubulin 4), as potentially druggable targets in GBM. Using data mining approach, we demonstrate that GLUT1 is overexpressed as a function of tumor grade in astrocytoma's and that its overexpression is associated with poorer prognosis. Using both mass spectrometry performed on hGBM (human glioblastoma patient specimen) and in silico modeling, we show that GLUT1 interacts with TUBB4, and more accurately demonstrates GLUT1's binding with fasentin. Proximity ligation assay (PLA) and immunoprecipitation studies confirm GLUT1 interaction with TUBB4. Treatment of GSC33 and GSC28 cells with TUBB4 inhibitor, CR-42-24, reduces the expression of GLUT1 however, TUBB4 expression is unaltered upon fasentin treatment. Using human pluripotent stem cell antibody array, we demonstrate reduced levels of Oct3/4, Nanog, Sox2, Sox17, Snail and VEGFR2 (Vascular endothelial growth factor receptor 2) upon CR-42-24 treatment. Overall, our data confirm that silencing TUBB4 or GLUT1 reduce GSC tumorsphere formation, self-renewal and proliferation in vitro. These findings suggest GLUT1 and its binding partner TUBB4 as druggable targets that warrant further investigation in GBM.

The Quiescent Metabolic Phenotype of Glioma Stem Cells

Introduction

Glioblastoma (GBM) is the most common primary malignant brain tumor in humans and, even with aggressive treatment that includes surgical resection, radiation (IR), and chemotherapy administration, prognosis is poor due to tumor recurrence. There is evidence that within GBMs a small number of glioma stem-like cells (GSLCs) exist, which are thought to be therapy resistant and are thus capable of repopulating a tumor after treatment. Like most cancers, GBMs largely employ aerobic glycolysis to create ATP, a phenomenon known as the Warburg Effect. There is no consensus on the metabolic characteristics of cancer stem cells. GSLCs have been shown to rely more heavily on oxidative phosphorylation, but there is also evidence that cancer stem cells can adapt their metabolism by fluctuating between energy pathways or acquiring intermediate metabolic phenotypes. We hypothesized that the metabolism of GSLCs differs from that of differentiated GBM tumor cell lines, and that the steady state metabolism would be differentially altered following radiation treatment.

Materials and Methods

We evaluated the oxygen consumption rate, extracellular acidification rate, and metabolic enzyme levels of GBM cell lines and GSLCs before and after irradiation using extracellular flux assays. We also measured absolute metabolite levels in these cells via mass spectroscopy with and without radiation treatment.

Results

GSLCs were found to be significantly more quiescent in comparison to adherent GBM cell lines, highlighted by lower glycolytic and maximal respiratory capacities as well as lower oxygen consumption and extracellular acidification rates. Analysis of individual metabolite concentrations revealed lower total metabolite concentrations overall but also elevated levels of metabolites in different energy pathways for GSLCs compared to GBM cell lines. Additionally, the metabolism of both GSLCs and GBM cell lines were found to be altered by IR.

Conclusions

While there is not one metabolic alteration that distinguishes irradiated GSLC metabolism from that of GBM cell lines, therapies targeting more metabolically quiescent tumor cells and thus the resistant GSLC population may increase a cancer's sensitivity to radiotherapy.

Targeting Hyaluronan Interactions for Glioblastoma Stem Cell Therapy

Even with rigorous treatments, glioblastoma multiforme (GBM) has an abysmal median survival rate, greatly due to the drug-resistant glioblastoma stem cell (GSC) population. GSCs are known to remodel their microenvironment, but the precise role of extracellular matrix components hyaluronic acid (HA) and hyaluronidases (HAases) on the GSC population is still largely unknown. Our objective was to determine how HAase can sensitize GSCs to chemotherapy drugs by disrupting the HA-CD44 signaling. GBM cell line U87-MG and patient-derived D456 cells were grown in GSC-enriching media and treated with HA or HAase. Expressions of GSC markers, HA-related genes, and drug resistance genes were measured via flow cytometry, confocal microscopy, and qRT-PCR. Proliferation after combined HAase and temozolomide (TMZ) treatment was measured via WST-8. HA supplementation promoted the expression of GSC markers and CD44 in GBM cells cultured in serum-free media. Conversely, HAase addition inhibited GSC gene expression while promoting CD44 expression. Finally, HAase sensitized GBM cells to TMZ. We propose a combined treatment of HAase and chemotherapy drugs by disrupting the stemness-promoting HA to target GSCs. This combination therapy shows promise even when temozolomide treatment alone causes resistance.

MicroRNA-451 Inhibits Migration of Glioblastoma while Making It More Susceptible to Conventional Therapy

Malignant glioblastoma (GBM, glioma) is the most common and aggressive primary adult brain tumor. The prognosis of GBM patients remains poor, despite surgery, radiation and chemotherapy. The major obstacles for successful remedy are invasiveness and therapy resistance of GBM cells. Invasive glioma cells leave primary tumor core and infiltrate surrounding normal brain leading to inevitable recurrence, even after surgical resection, radiation and chemotherapy. Therapy resistance allowing for selection of more aggressive and resistant sub-populations including GBM stem-like cells (GSCs) upon treatment is another serious impediment to successful treatment. Through their regulation of multiple genes, microRNAs can orchestrate complex programs of gene expression and act as master regulators of cellular processes. MicroRNA-based therapeutics could thus impact broad cellular programs, leading to inhibition of invasion and sensitization to radio/chemotherapy. Our data show that miR-451 attenuates glioma cell migration in vitro and invasion in vivo. In addition, we have found that miR-451 sensitizes glioma cells to conventional chemo- and radio-therapy. Our data also show that miR-451 is regulated in vivo by AMPK pathway and that AMPK/miR-451 loop has the ability to switch between proliferative and migratory pattern of glioma cells behavior. We therefore postulate that AMPK/miR-451 negative reciprocal feedback loop allows GBM cells/GSCs to adapt to tumor “ecosystem” by metabolic and behavioral flexibility, and that disruption of such a loop reduces invasiveness and diminishes therapy resistance.

Cholesterol Metabolism: A Potential Therapeutic Target in Glioblastoma

Glioblastoma is a highly lethal adult brain tumor with no effective treatments. In this review, we discuss the potential to target cholesterol metabolism as a new strategy for treating glioblastomas. Twenty percent of cholesterol in the body is in the brain, yet the brain is unique among organs in that it has no access to dietary cholesterol and must synthesize it de novo. This suggests that therapies targeting cholesterol synthesis in brain tumors might render their effects without compromising cell viability in other organs. We will describe cholesterol synthesis and homeostatic feedback pathways in normal brain and brain tumors, as well as various strategies for targeting these pathways for therapeutic intervention.

Immune evasion mediated by PD-L1 on glioblastoma-derived extracellular vesicles

Binding of programmed death ligand-1 (PD-L1) to programmed cell death protein-1 (PD1) leads to cancer immune evasion via inhibition of T cell function. One of the defining characteristics of glioblastoma, a universally fatal brain cancer, is its profound local and systemic immunosuppression. Glioblastoma has also been shown to generate extracellular vesicles (EVs), which may play an important role in tumor progression. We thus hypothesized that glioblastoma EVs may be important mediators of immunosuppression and that PD-L1 could play a role. We show that glioblastoma EVs block T cell activation and proliferation in response to T cell receptor stimulation. PD-L1 was expressed on the surface of some, but not of all, glioblastoma-derived EVs, with the potential to directly bind to PD1. An anti-PD1 receptor blocking antibody significantly reversed the EV-mediated blockade of T cell activation but only when PD-L1 was present on EVs. When glioblastoma PD-L1 was up-regulated by IFN-γ, EVs also showed some PD-L1-dependent inhibition of T cell activation. PD-L1 expression correlated with the mesenchymal transcriptome profile and was anatomically localized in the perinecrotic and pseudopalisading niche of human glioblastoma specimens. PD-L1 DNA was present in circulating EVs from glioblastoma patients where it correlated with tumor volumes of up to 60 cm³. These results suggest that PD-L1 on EVs may be another mechanism for glioblastoma to suppress antitumor immunity and support the potential of EVs as biomarkers in tumor patients.

Modulation of miR-21 signaling by MPS1 in human glioblastoma

Monopolar spindle 1 (MPS1) is an essential spindle assembly checkpoint (SAC) kinase involved in determining spindle integrity. Beyond its mitotic functions, it has been implicated in several other signaling pathways. Our earlier studies have elaborated on role of MPS1 in glioblastoma (GBM) radiosensitization. In this study using reverse phase protein arrays (RPPAs), we assessed MPS1 mediated cell signaling pathways and demonstrated that inhibiting MPS1 could upregulate the expression of the tumor suppressor PDCD4 and MSH2 genes, by down regulating micro RNA-21 (miR-21). In GBMs miR-21 expression is significantly elevated and is associated with chemo and radioresistance. Both MPS1 and miR-21 depletion suppressed GBM cell proliferation, whereas, ectopic expression of miR-21 rescued GBM cell growth from MPS1 inhibition. Further, we demonstrate that MPS1 mediates phosphorylation of SMAD3 but not SMAD2 in GBM cells; A possible mechanism behind miR-21 modulation by MPS1. Collectively, our results shed light onto an important role of MPS1 in TGF-β/SMAD signaling via miR-21 regulation. We also, show the prognostic effect of miR-21, PDCD4 and MSH2 levels to patient survival across different GBM molecular subtypes. This scenario in which miR-21 is modulated by MPS1 inhibition may be exploited as a potential target for effective GBM therapy.

Survival kinase genes present prognostic significance in glioblastoma

Cancer biomarkers with a strong predictive power for diagnosis/prognosis and a potential to be therapeutic targets have not yet been fully established. Here we employed a loss-of-function screen in glioblastoma (GBM), an infiltrative brain tumor with a dismal prognosis, and identified 20 survival kinase genes (SKGs). Survival analyses using The Cancer Genome Atlas (TCGA) datasets revealed that the expression of CDCP1, CDKL5, CSNK1E, IRAK3, LATS2, PRKAA1, STK3, TBRG4, and ULK4 stratified GBM prognosis with or without temozolomide (TMZ) treatment as a covariate. For the first time, we found that GBM patients with a high level of NEK9 and PIK3CB had a greater chance of having recurrent tumors. The expression of CDCP1, IGF2R, IRAK3, LATS2, PIK3CB, ULK4, or VRK1 in primary GBM tumors was associated with recurrence-related prognosis. Notably, the level of PIK3CB in recurrent tumors was much higher than that in newly diagnosed ones. Congruent with these results, genes in the PI3K/AKT pathway showed a significantly strong correlation with recurrence rate, further highlighting the pivotal role of PIK3CB in the disease progression. Importantly, 17 SKGs together presented a novel GBM prognostic signature. SKGs identified herein are associated with recurrence rate and present prognostic significance in GBM, thereby becoming attractive therapeutic targets.

Discovery and validation of a glioblastoma co-expressed gene module

Tumors exhibit complex patterns of aberrant gene expression. Using a knowledge-independent, noise-reducing gene co-expression network construction software called KINC, we created multiple RNAseq-based gene co-expression networks relevant to brain and glioblastoma biology. In this report, we describe the discovery and validation of a glioblastoma-specific gene module that contains 22 co-expressed genes. The genes are upregulated in glioblastoma relative to normal brain and lower grade glioma samples; they are also hypo-methylated in glioblastoma relative to lower grade glioma tumors. Among the proneural, neural, mesenchymal, and classical glioblastoma subtypes, these genes are most-highly expressed in the mesenchymal subtype. Furthermore, high expression of these genes is associated with decreased survival across each glioblastoma subtype. These genes are of interest to glioblastoma biology and our gene interaction discovery and validation workflow can be used to discover and validate co-expressed gene modules derived from any co-expression network.

A Dexamethasone-regulated Gene Signature Is Prognostic for Poor Survival in Glioblastoma Patients

BACKGROUND

Dexamethasone is reported to induce both tumor-suppressive and tumor-promoting effects. The purpose of this study was to identify the genomic impact of dexamethasone in glioblastoma stem cell (GSC) lines and its prognostic value; furthermore, to identify drugs that can counter these side effects of dexamethasone exposure.

METHODS

We utilized 3 independent GSC lines with tumorigenic potential for this study. Whole-genome expression profiling and pathway analyses were done with dexamethasone-exposed and control cells. GSCs were also co-exposed to dexamethasone and temozolomide. Risk scores were calculated for most affected genes, and their associations with survival in The Cancer Genome Atlas and Repository of Molecular Brain Neoplasia Data databases. In silico Connectivity Map analysis identified camptothecin as antagonist to dexamethasone-induced negative effects.

RESULTS

Pathway analyses predicted an activation of dexamethasone network (z-score: 2.908). Top activated canonical pathways included "role of breast cancer 1 in DNA damage response" (P=1.07E-04). GSCs were protected against temozolomide-induced apoptosis when coincubated with dexamethasone. Altered cellular functions included cell movement, cell survival, and apoptosis with z-scores of 2.815, 5.137, and -3.122, respectively. CCAAT/enhancer binding protein beta (CEBPB) was activated in a dose dependent manner specifically in slow-dividing "stem-like" cells. CEBPB was activated in dexamethasone-treated orthotopic tumors. Patients with high risk scores had significantly shorter survival. Camptothecin was validated as potential partial neutralizer of dexamethasone-induced oncogenic effects.

CONCLUSIONS

Dexamethasone exposure induces a genetic program and CEBPB expression in GSCs that adversely affects key cellular functions and response to therapeutics. High risk scores associated with these genes have negative prognostic value in patients. Our findings further suggest camptothecin as a potential neutralizer of adverse dexamethasone-mediated effects.

A three-microRNA signature identifies two subtypes of glioblastoma patients with different clinical outcomes

Glioblastoma multiforme (GBM) is the most common and malignant primary brain tumor in adults, characterized by aggressive growth, limited response to therapy, and inexorable recurrence. Because of the extremely unfavorable prognosis of GBM, it is important to develop more effective diagnostic and therapeutic strategies based on biologically and clinically relevant patient stratification systems. Analyzing a collection of patient‐derived GBM stem‐like cells (GSCs) by gene expression profiling, nuclear magnetic resonance spectroscopy, and signal transduction pathway activation, we identified two GSC clusters characterized by different clinical features. Due to the widely documented role played by microRNAs (miRNAs) in the tumorigenesis process, in this study we explored whether these two GBM patient subtypes could also be discriminated by different miRNA signatures. Global miRNA expression pattern was analyzed by oblique principal component analysis and principal component analysis. By a combined inferential strategy on PCA results, we identified a reduced set of three miRNAs – miR‐23a, miR‐27a, and miR‐9* (miR‐9‐3p) – able to discriminate the proneural‐ and mesenchymal‐like GSC phenotypes as well as mesenchymal and proneural subtypes of primary GBM included in The Cancer Genome Atlas (TCGA) data set. Kaplan–Meier analysis showed a significant correlation between the selected miRNAs and overall survival in 429 GBM specimens from TCGA‐identifying patients who had an unfavorable outcome. The survival prognostic capability of the three‐miRNA signatures could have important implications for the understanding of the biology of GBM subtypes and could be useful in patient stratification to facilitate interpretation of results from clinical trials.

MicroRNA Signatures and Molecular Subtypes of Glioblastoma: The Role of Extracellular Transfer

Despite the importance of molecular subtype classification of glioblastoma (GBM), the extent of extracellular vesicle (EV)-driven molecular and phenotypic reprogramming remains poorly understood. To reveal complex subpopulation dynamics within the heterogeneous intratumoral ecosystem, we characterized microRNA expression and secretion in phenotypically diverse subpopulations of patient-derived GBM stem-like cells (GSCs). As EVs and microRNAs convey information that rearranges the molecular landscape in a cell type-specific manner, we argue that intratumoral exchange of microRNA augments the heterogeneity of GSC that is reflected in highly heterogeneous profile of microRNA expression in GBM subtypes.

•

MicroRNA signatures reveal tissue heterogeneity in defined glioblastoma subtypes

•

GSC EV/microRNA acts via cell-dependent targeting, propagating intratumoral heterogeneity

•

EV/microRNAs modify molecular landscape, acting in tumor anatomic sites

Tumor-derived fibulin-3 activates pro-invasive NF-κB signaling in glioblastoma cells and their microenvironment

Molecular profiling of glioblastomas has revealed the presence of key signaling hubs that contribute to tumor progression and acquisition of resistance. One of these main signaling mechanisms is the NF-κB pathway, which integrates multiple extracellular signals into transcriptional programs for tumor growth, invasion, and maintenance of the tumor-initiating population. We show here that an extracellular protein released by glioblastoma cells, fibulin-3, drives oncogenic NF-κB in the tumor and increases NF-κB activation in peritumoral astrocytes. Fibulin-3 expression correlates with a NF-κB-regulated “invasive signature” linked to poorer survival, being a possible tissue marker for regions of active tumor progression. Accordingly, fibulin-3 promotes glioblastoma invasion in a manner that requires NF-κB activation both in the tumor cells and their microenvironment. Mechanistically, we found that fibulin-3 activates the metalloprotease ADAM17 by competing with its endogenous inhibitor, TIMP3. This results in sustained release of soluble TNFα by ADAM17, which in turn activates TNF receptors and canonical NF-κB signaling. Taken together, our results underscore fibulin-3 as a novel extracellular signal with strong activating effect on NF-κB in malignant gliomas. Because fibulin-3 is produced de novo in these tumors and is absent from normal brain we propose that targeting the fibulin-3/NF-κB axis may provide a novel avenue to disrupt oncogenic NF-κB signaling in combination therapies for malignant brain tumors.

Computational analysis of the mesenchymal signature landscape in gliomas

BACKGROUND

Epithelial to mesenchymal transition, and mimicking processes, contribute to cancer invasion and metastasis, and are known to be responsible for resistance to various therapeutic agents in many cancers. While a number of studies have proposed molecular signatures that characterize the spectrum of such transition, more work is needed to understand how the mesenchymal signature (MS) is regulated in non-epithelial cancers like gliomas, to identify markers with the most prognostic significance, and potential for therapeutic targeting.

RESULTS

Computational analysis of 275 glioma samples from "The Cancer Genome Atlas" was used to identify the regulatory changes between low grade gliomas with little expression of MS, and high grade glioblastomas with high expression of MS. TF (transcription factor)-gene regulatory networks were constructed for each of the cohorts, and 5 major pathways and 118 transcription factors were identified as involved in the differential regulation of the networks. The most significant pathway - Extracellular matrix organization - was further analyzed for prognostic relevance. A 20-gene signature was identified as having prognostic significance (HR (hazard ratio) 3.2, 95% CI (confidence interval) = 1.53-8.33), after controlling for known prognostic factors (age, and glioma grade). The signature's significance was validated in an independent data set. The putative stem cell marker CD44 was biologically validated in glioma cell lines and brain tissue samples.

CONCLUSIONS

Our results suggest that the differences between low grade gliomas and high grade glioblastoma are associated with differential expression of the signature genes, raising the possibility that targeting these genes might prolong survival in glioma patients.

Disabled cell density sensing leads to dysregulated cholesterol synthesis in glioblastoma

A hallmark of cellular transformation is the evasion of contact-dependent inhibition of growth. To find new therapeutic targets for glioblastoma, we looked for pathways that are inhibited by high cell density in astrocytes but not in glioma cells. Here we report that glioma cells have disabled the normal controls on cholesterol synthesis. At high cell density, astrocytes turn off cholesterol synthesis genes and have low cholesterol levels, but glioma cells keep this pathway on and maintain high cholesterol. Correspondingly, cholesterol pathway upregulation is associated with poor prognosis in glioblastoma patients. Densely-plated glioma cells increase oxygen consumption, aerobic glycolysis, and the pentose phosphate pathway to synthesize cholesterol, resulting in a decrease in reactive oxygen species, TCA cycle intermediates, and ATP. This constitutive cholesterol synthesis is controlled by the cell cycle, as it can be turned off by cyclin-dependent kinase inhibitors and it correlates with disabled cell cycle control though loss of p53 and RB. Finally, glioma cells, but not astrocytes, are sensitive to cholesterol synthesis inhibition downstream of the mevalonate pathway, suggesting that specifically targeting cholesterol synthesis might be an effective treatment for glioblastoma.

Reverse phase protein arrays enable glioblastoma molecular subtyping

In the present study we investigated the phosphorylation status of the 12 most important signaling cascades in glioblastomas. More than 60 tumor and control biopsies from tumor center and periphery (based on neuronavigation) were subjected to selective protein expression analysis using reverse-phase protein arrays (RPPA) incubated with antibodies against posttranslationally modified cancer pathway proteins. The ratio between phosphorylated (or modified) and non-phosphorylated protein was assessed. All samples were histopathologically validated and proteomic profiles correlated with clinical and survival data. By RPPA, we identified three distinct activation patterns within glioblastoma defined by the ratios of pCREB1/CREB1, NOTCH-ICD/NOTCH1, and pGSK3β/GSK3β, respectively. These subclasses demonstrated distinct overall survival patterns in a cohort of patients from a single-institution and in an analysis of publicly available data. In particular, a high pGSK3β/GSK3β-ratio was associated with a poor survival. Wnt-activation/GSK3β-inhibition in U373 and U251 cell lines halted glioma cell proliferation and migration. Gene expression analysis was used as an internal quality control of baseline proteomic data. The protein expression and phosphorylation had a higher resolution, resulting in a better class-subdivision than mRNA based stratification data. Patients with different proteomic profiles from multiple biopsies showed a worse overall survival. The CREB1-, NOTCH1-, GSK3β-phosphorylation status correlated with glioma grades. RPPA represent a fast and reliable tool to supplement morphological diagnosis with pathway-specific information in individual tumors. These data can be exploited for molecular stratification and possible combinatorial treatment planning. Further, our results may optimize current glioma grading algorithms.

MNK Inhibition Disrupts Mesenchymal Glioma Stem Cells and Prolongs Survival in a Mouse Model of Glioblastoma

Glioblastoma multiforme remains the deadliest malignant brain tumor, with glioma stem cells (GSC) contributing to treatment resistance and tumor recurrence. We have identified MAPK-interacting kinases (MNK) as potential targets for the GSC population in glioblastoma multiforme. Isoform-level subtyping using The Cancer Genome Atlas revealed that both MNK genes (MKNK1 and MKNK2) are upregulated in mesenchymal glioblastoma multiforme as compared with other subtypes. Expression of MKNK1 is associated with increased glioma grade and correlated with the mesenchymal GSC marker, CD44, and coexpression of MKNK1 and CD44 predicts poor survival in glioblastoma multiforme. In established and patient-derived cell lines, pharmacologic MNK inhibition using LY2801653 (merestinib) inhibited phosphorylation of the eukaryotic translation initiation factor 4E, a crucial effector for MNK-induced mRNA translation in cancer cells and a marker of transformation. Importantly, merestinib inhibited growth of GSCs grown as neurospheres as determined by extreme limiting dilution analysis. When the effects of merestinib were assessed in vivo using an intracranial xenograft mouse model, improved overall survival was observed in merestinib-treated mice. Taken together, these data provide strong preclinical evidence that pharmacologic MNK inhibition targets mesenchymal glioblastoma multiforme and its GSC population.

IMPLICATIONS

These findings raise the possibility of MNK inhibition as a viable therapeutic approach to target the mesenchymal subtype of glioblastoma multiforme. Mol Cancer Res; 14(10); 984-93. ©2016 AACR.

The Long Non-coding RNA HIF1A-AS2 Facilitates the Maintenance of Mesenchymal Glioblastoma Stem-like Cells in Hypoxic Niches

Long non-coding RNAs (lncRNAs) have an undefined role in the pathobiology of glioblastoma multiforme (GBM). These tumors are genetically and phenotypically heterogeneous with transcriptome subtype-specific GBM stem-like cells (GSCs) that adapt to the brain tumor microenvironment, including hypoxic niches. We identified hypoxia-inducible factor 1 alpha-antisense RNA 2 (HIF1A-AS2) as a subtype-specific hypoxia-inducible lncRNA, upregulated in mesenchymal GSCs. Its deregulation affects GSC growth, self-renewal, and hypoxia-dependent molecular reprogramming. Among the HIF1A-AS2 interactome, IGF2BP2 and DHX9 were identified as direct partners. This association was needed for maintenance of expression of their target gene, HMGA1. Downregulation of HIF1A-AS2 led to delayed growth of mesenchymal GSC tumors, survival benefits, and impaired expression of HMGA1 in vivo. Our data demonstrate that HIF1A-AS2 contributes to GSCs' speciation and adaptation to hypoxia within the tumor microenvironment, acting directly through its interactome and targets and indirectly by modulating responses to hypoxic stress depending on the subtype-specific genetic context.

Extracellular Vesicles from High-Grade Glioma Exchange Diverse Pro-oncogenic Signals That Maintain Intratumoral Heterogeneity

A lack of experimental models of tumor heterogeneity limits our knowledge of the complex subpopulation dynamics within the tumor ecosystem. In high-grade gliomas (HGG), distinct hierarchical cell populations arise from different glioma stem-like cell (GSC) subpopulations. Extracellular vesicles (EV) shed by cells may serve as conduits of genetic and signaling communications; however, little is known about how HGG heterogeneity may impact EV content and activity. In this study, we performed a proteomic analysis of EVs isolated from patient-derived GSC of either proneural or mesenchymal subtypes. EV signatures were heterogeneous, but reflected the molecular make-up of the GSC and consistently clustered into the two subtypes. EV-borne protein cargos transferred between proneural and mesenchymal GSC increased protumorigenic behaviors in vitro and in vivo Clinically, analyses of HGG patient data from the The Cancer Genome Atlas database revealed that proneural tumors with mesenchymal EV signatures or mesenchymal tumors with proneural EV signatures were both associated with worse outcomes, suggesting influences by the proportion of tumor cells of varying subtypes in tumors. Collectively, our findings illuminate the heterogeneity among tumor EVs and the complexity of HGG heterogeneity, which these EVs help to maintain. Cancer Res; 76(10); 2876-81. ©2016 AACR.

Metabolic/Proteomic Signature Defines Two Glioblastoma Subtypes With Different Clinical Outcome

Glioblastoma (GBM) is one of the deadliest human cancers. Because of the extremely unfavorable prognosis of GBM, it is important to develop more effective diagnostic and therapeutic strategies based on biologically and clinically relevant subclassification systems. Analyzing a collection of seventeen patient-derived glioblastoma stem-like cells (GSCs) by gene expression profiling, NMR spectroscopy and signal transduction pathway activation, we identified two GSC clusters, one characterized by a pro-neural-like phenotype and the other showing a mesenchymal-like phenotype. Evaluating the levels of proteins differentially expressed by the two GSC clusters in the TCGA GBM sample collection, we found that SRC activation is associated with a GBM subgroup showing better prognosis whereas activation of RPS6, an effector of mTOR pathway, identifies a subgroup with a worse prognosis. The two clusters are also differentiated by NMR spectroscopy profiles suggesting a potential prognostic stratification based on metabolic evaluation. Our data show that the metabolic/proteomic profile of GSCs is informative of the genomic/proteomic GBM landscape, which differs among tumor subtypes and is associated with clinical outcome.

Quality control material for the detection of somatic mutations in fixed clinical specimens by next-generation sequencing

BACKGROUND

Targeted next generation sequencing (NGS) technology to assess the mutational status of multiple genes on formalin-fixed, paraffin embedded (FFPE) tumors is rapidly being adopted in clinical settings, where quality control (QC) practices are required. Establishing reliable FFPE QC materials for NGS can be challenging and/or expensive. Here, we established a reliable and cost-effective FFPE QC material for routine utilization in the Ion AmpliSeq™ Cancer Hotspot Panel v2 (CHP2) assay.

METHODS

The performance characteristics of the CHP2 assay were determined by sequencing various cell line mixtures and 55 different FFPE tumors on the Ion Torrent PGM platform. A FFPE QC material was prepared from a mixture of cell lines derived from different cancers, comprising single nucleotide variants and small deletions on actionable genes at different allelic frequencies.

RESULTS

The CHP2 assay performed with high precision and sensitivity when custom variant calling pipeline parameters where established. In addition, all expected somatic variants in the QC material were consistently called at variant frequencies ranging from 9.1 % (CV = 11.1 %) to 37.9 % (CV = 2.8 %).

CONCLUSIONS

The availability of a reliable and cost-effective QC material is instrumental in assessing the performance of this or any targeted NGS assay that detects somatic variants in fixed solid tumor specimens.

Kinome-wide shRNA screen identifies the receptor tyrosine kinase AXL as a key regulator for mesenchymal glioblastoma stem-like cells

Glioblastoma is a highly lethal cancer for which novel therapeutics are urgently needed. Two distinct subtypes of glioblastoma stem-like cells (GSCs) were recently identified: mesenchymal (MES) and proneural (PN). To identify mechanisms to target the more aggressive MES GSCs, we combined transcriptomic expression analysis and kinome-wide short hairpin RNA screening of MES and PN GSCs. In comparison to PN GSCs, we found significant upregulation and phosphorylation of the receptor tyrosine kinase AXL in MES GSCs. Knockdown of AXL significantly decreased MES GSC self-renewal capacity in vitro and inhibited the growth of glioblastoma patient-derived xenografts. Moreover, inhibition of AXL with shRNA or pharmacologic inhibitors also increased cell death significantly more in MES GSCs. Clinically, AXL expression was elevated in the MES GBM subtype and significantly correlated with poor prognosis in multiple cancers. In conclusion, we identified AXL as a potential molecular target for novel approaches to treat glioblastoma and other solid cancers.

•

shRNA screen identified kinases that alter GSC viability in a subtype-dependent manner

•

AXL is highly expressed in mesenchymal GSCs

•

Targeting AXL decreases mesenchymal GSC self-renewal, viability, and tumorigenicity

•

AXL expression predicts poor prognosis in several tumor types

miR-182 integrates apoptosis, growth, and differentiation programs in glioblastoma

Glioblastoma multiforme (GBM) is a lethal, therapy-resistant brain cancer consisting of numerous tumor cell subpopulations, including stem-like glioma-initiating cells (GICs), which contribute to tumor recurrence following initial response to therapy. Here, we identified miR-182 as a regulator of apoptosis, growth, and differentiation programs whose expression level is correlated with GBM patient survival. Repression of Bcl2-like12 (Bcl2L12), c-Met, and hypoxia-inducible factor 2α (HIF2A) is of central importance to miR-182 anti-tumor activity, as it results in enhanced therapy susceptibility, decreased GIC sphere size, expansion, and stemness in vitro. To evaluate the tumor-suppressive function of miR-182 in vivo, we synthesized miR-182-based spherical nucleic acids (182-SNAs); i.e., gold nanoparticles covalently functionalized with mature miR-182 duplexes. Intravenously administered 182-SNAs penetrated the blood-brain/blood-tumor barriers (BBB/BTB) in orthotopic GBM xenografts and selectively disseminated throughout extravascular glioma parenchyma, causing reduced tumor burden and increased animal survival. Our results indicate that harnessing the anti-tumor activities of miR-182 via safe and robust delivery of 182-SNAs represents a novel strategy for therapeutic intervention in GBM.

Biomarkers of glioblastoma multiforme

Glioblastoma multiforme (GBM) is the most common and lethal primary malignancy of the central nervous system. Modern treatments using surgery and/or chemotherapy and/or radiotherapy are improving survival of patients, but prognosis is still very poor, depending inter alia on the patients' individual genomic traits. Most GBMs are primary; however, secondary GBMs have a better prognosis. Aberrant gene expression and copy number alterations make it possible to identify four subtypes: classical, mesenchymal, proneural, and neural. More and more biomarkers continue to be identified in GBM patients. Such biomarkers are related with varying degrees of specificity to one or more of GBM's subtypes and, in many instances, may provide useful information about prognosis. Biomarkers fall into either the imaging or molecular category. Molecular biomarkers are identified by use of such platforms as genomics, proteomics, and metabolomics. In the future, biomarkers, either individually or in some combination, will more reliably identify the pathogenic type of GBM and determine choice of therapy.