Comparative genomic and genetic analysis of glioblastoma-derived brain tumor-initiating cells and their parent tumors

A novel disulfide death-related genes prognostic signature identifies the role of IPO4 in glioma progression

BACKGROUND

"Disulfide death," a form of cellular demise, is triggered by the abnormal accumulation of intracellular disulfides under conditions of glucose deprivation. However, its role in the prognosis of glioma remains undetermined. Therefore, the main objective of this study is to establish prognostic signature based on disulfide death-related genes (DDRGs) and to provide new solutions in choosing the effective treatment of glioma.

METHODS

The RNA transcriptome, clinical information, and mutation data of glioma samples were sourced from The Cancer Genome Atlas (TCGA) and the Chinese Glioma Genome Atlas (CGGA), while normal samples were obtained from the Genotype-Tissue Expression (GTEx). DDRGs were compiled from previous studies and selected through differential analysis and univariate Cox regression analysis. The molecular subtypes were determined through consensus clustering analysis. Further, LASSO analysis was employed to select characteristic genes, and subsequently, a risk model comprising seven DDRGs was constructed based on multivariable Cox analysis. Kaplan-Meier survival curves were employed to assess survival differences between high and low-risk groups. Additionally, functional analyses (GO, KEGG, GSEA) were conducted to explore the potential biological functions and signaling pathways of genes associated with the model. The study also explored immune checkpoint (ICP) genes, immune cell infiltration levels, and immune stromal scores. Finally, the effect of Importin-4(IPO4) on glioma has been further confirmed through RT-qPCR, Western blot, and cell functional experiments.

RESULTS

7 genes associated with disulfide death were obtained and two subgroups of patients with different prognosis and clinical characteristics were identified. Risk signature was subsequently developed and proved to serve as an prognostic predictor. Notably, the high-risk group exhibited an immunosuppressive microenvironment characterized by a high concentration of M2 macrophages and regulatory T cells (Tregs). In contrast, the low-risk group showed lower half-maximal inhibitory concentration (IC50) values. Therefore, patients in the high-risk group may benefit more from immunotherapy, while patients in the low-risk group may benefit more from chemotherapy. In addition, in vitro experiments have shown that inhibition of the expression of IPO4 leads to a significant reduction in the proliferation, migration, and invasion of glioma cells.

CONCLUSION

This study identified two glioma subtypes and constructed a prognostic signature based on DDRGs. The signature has the potential to optimize the selection of patients for immune- and chemotherapy and provided a potential therapeutic target for glioma.

Enhancing temozolomide antiglioma response by inhibiting O6-methylguanine-DNA methyltransferase with selected phytochemicals: in silico and in vitro approach

The aim of our study was to investigate the potential of rutin, catechin, dehydrozingerone, naringenin, and quercetin, both alone and in combination with temozolomide, to inhibit the expression of O6-methylguanine-DNA methyltransferase (MGMT) in glioma cells. MGMT has been shown to be a major cause of temozolomide resistance in glioma. Our study used both in silico and in vitro methods to assess the inhibitory activity of these phytochemicals on MGMT, with the goal of identifying the most effective combination of compounds for reducing temozolomide resistance. After conducting an initial in silico screening of natural compounds against MGMT protein, five phytochemicals were chosen based on their high docking scores and favorable binding energies. From the molecular docking and simulation studies, we found that quercetin showed a good inhibitory effect of MGMT with its high binding affinity. C6 glioma cells showed increased cytotoxicity when treated with the temozolomide and quercetin combination. It was understood from the isobologram and combination index plot that the drug combination showed a synergistic effect at the lowest dose. Quercetin when combined with temozolomide significantly decreased the MGMT levels in C6 cells in comparison with the other drugs as estimated by ELISA. The percentage of apoptotic cells increased significantly in the temozolomide-quercetin group indicating the potency of quercetin in decreasing the resistance of temozolomide as confirmed by acridine orange/ethidium bromide staining. Our experiment hence suggests that temozolomide resistance can be reduced by combining the drug with quercetin which will serve as an effective therapeutic target for glioblastoma treatment.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03821-7.

Identification of M5c regulator-medicated methylation modification patterns for prognosis and immune microenvironment in glioma

Glioma is a common intracranial tumor and is generally associated with poor prognosis. Recently, numerous studies illustrated the importance of 5-methylcytosine (m5C) RNA modification to tumorigenesis. However, the prognostic value and immune correlation of m5C in glioma remain unclear. We obtained RNA expression and clinical information from The Cancer Genome Atlas (TCGA) and The Chinese Glioma Genome Atlas (CGGA) datasets to analyze. Nonnegative matrix factorization (NMF) was used to classify patients into two subgroups and compare these patients in survival and clinicopathological characteristics. CIBERSORT and single-sample gene-set algorithm (ssGSEA) methods were used to investigate the relationship between m5C and the immune environment. The Weighted correlation network analysis (WGCNA) and univariate Cox proportional hazard model (CoxPH) were used to construct a m5C-related signature. Most of m5C RNA methylation regulators presented differential expression and prognostic values. There were obvious relationships between immune infiltration cells and m5C regulators, especially NSUN7. In the m5C-related module from WGCNA, we found SEPT3, CHI3L1, PLBD1, PHYHIPL, SAMD8, RAP1B, B3GNT5, RER1, PTPN7, SLC39A1, and MXI1 were prognostic factors for glioma, and they were used to construct the signature. The great significance of m5C-related signature in predicting the survival of patients with glioma was confirmed in the validation sets and CGGA cohort.

Establishment and Comprehensive Molecular Characterization of an Immortalized Glioblastoma Cell Line from a Brazilian Patient

Glioblastoma (GBM) is the most common and aggressive primary brain tumor in adults, with few effective treatment strategies. The research on the development of new treatments is often constrained by the limitations of preclinical models, which fail to accurately replicate the disease's essential characteristics. Herein, we describe the obtention, molecular, and functional characterization of the GBM33 cell line. This cell line belongs to the GBM class according to the World Health Organization 2021 Classification of Central Nervous System Tumors, identified by methylation profiling. GBM33 expresses the astrocytic marker GFAP, as well as markers of neuronal origin commonly expressed in GBM cells, such as βIII-tubulin and neurofilament. Functional assays demonstrated an increased growth rate when compared to the U87 commercial cell line and a similar sensitivity to temozolamide. GBM33 cells retained response to serum starvation, with reduced growth and diminished activation of the Akt signaling pathway. Unlike LN-18 and LN-229 commercial cell lines, GBM33 is able to produce primary cilia upon serum starvation. In summary, the successful establishment and comprehensive characterization of this GBM cell line provide researchers with invaluable tools for studying GBM biology, identifying novel therapeutic targets, and evaluating the efficacy of potential treatments.

Ex vivo drug sensitivity screening predicts response to temozolomide in glioblastoma patients and identifies candidate biomarkers

BACKGROUND

Patient-derived glioma stem-like cells (GSCs) have become the gold-standard in neuro-oncological research; however, it remains to be established whether loss of in situ microenvironment affects the clinically-predictive value of this model. We implemented a GSC monolayer system to investigate in situ-in vitro molecular correspondence and the relationship between in vitro and patient response to temozolomide (TMZ).

METHODS

DNA/RNA-sequencing was performed on 56 glioblastoma tissues and 19 derived GSC cultures. Sensitivity to TMZ was screened across 66 GSC cultures. Viability readouts were related to clinical parameters of corresponding patients and whole-transcriptome data.

RESULTS

Tumour DNA and RNA sequences revealed strong similarity to corresponding GSCs despite loss of neuronal and immune interactions. In vitro TMZ screening yielded three response categories which significantly correlated with patient survival, therewith providing more specific prediction than the binary MGMT marker. Transcriptome analysis identified 121 genes related to TMZ sensitivity of which 21were validated in external datasets.

CONCLUSION

GSCs retain patient-unique hallmark gene expressions despite loss of their natural environment. Drug screening using GSCs predicted patient response to TMZ more specifically than MGMT status, while transcriptome analysis identified potential biomarkers for this response. GSC drug screening therefore provides a tool to improve drug development and precision medicine for glioblastoma.

Oncolytic herpes simplex viruses for the treatment of glioma and targeting glioblastoma stem-like cells

Glioblastoma (GBM) is one of the most lethal cancers, having a poor prognosis and a median survival of only about 15 months with standard treatment (surgery, radiation, and chemotherapy), which has not been significantly extended in decades. GBM demonstrates remarkable cellular heterogeneity, with glioblastoma stem-like cells (GSCs) at the apex. GSCs are a subpopulation of GBM cells that possess the ability to self-renew, differentiate, initiate tumor formation, and manipulate the tumor microenvironment (TME). GSCs are no longer considered a static population of cells with specific markers but are quite flexible phenotypically and in driving tumor heterogeneity and therapeutic resistance. In light of these features, they are a critical target for successful GBM therapy. Oncolytic viruses, in particular oncolytic herpes simplex viruses (oHSVs), have many attributes for therapy and are promising agents to target GSCs. oHSVs are genetically-engineered to selectively replicate in and kill cancer cells, including GSCs, but not normal cells. Moreover, oHSV can induce anti-tumor immune responses and synergize with other therapies, such as chemotherapy, DNA repair inhibitors, and immune checkpoint inhibitors, to potentiate treatment effects and reduce GSC populations that are partly responsible for chemo- and radio-resistance. Herein, we present an overview of GSCs, activity of different oHSVs, clinical trial results, and combination strategies to enhance efficacy, including therapeutic arming of oHSV. Throughout, the therapeutic focus will be on GSCs and studies specifically targeting these cells. Recent clinical trials and approval of oHSV G47Δ in Japan for patients with recurrent glioma demonstrate the efficacy and promise of oHSV therapy.

CD95 gene deletion may reduce clonogenic growth and invasiveness of human glioblastoma cells in a CD95 ligand-independent manner

CD95 (Fas/APO-1) is a multifunctional cell surface receptor with antithetic roles. First described to mediate cell death, interactions of CD95 with its natural ligand, CD95L, have also been described to induce tumor-promoting signaling leading to proliferation, invasion and stem cell maintenance, mainly in cancer cells that are resistant to CD95-mediated apoptosis. While activation of CD95-mediated apoptosis in cancer cells may not be clinically practicable due to toxicity, inhibition of tumor-promoting CD95 signaling holds therapeutic potential. In the present study, we characterized CD95 and CD95L expression in human glioma-initiating cells (GIC), a glioblastoma cell population with stem cell features, and investigated the consequences of CRISPR-Cas9-mediated CD95 or CD95L gene deletion. In vitro, GIC expressed CD95 but not CD95L and were sensitive to CD95-mediated apoptosis. Upon genetic deletion of CD95, GIC acquired resistance to CD95L-induced apoptosis but exhibited inferior clonogenic growth, sphere-forming capacity, and invasiveness compared with control cells, suggesting the existence of CD95L-independent constitutive CD95 signaling with tumor-promoting properties in GIC. In vivo, GIC expressed CD95 and a non-canonical form of CD95L lacking the CD95-binding region. CD95 genetic deletion did not prolong survival in immunocompromised GIC-bearing mice. Altogether, these data indicate that canonical CD95L may not be expressed in human GIC and suggest the existence of a CD95L-independent CD95-signaling pathway that maintains some malignancy traits of GIC. The lack of altered survival of tumor-bearing mice after genetic deletion of CD95 suggests that CD95 signaling is not essential to maintain the growth of human GIC xenografted into the brains of nude mice. The ligand-independent tumor-promoting role of constitutive CD95 in our GIC models in vitro highlights the complexity and challenges associated with targeting CD95 with therapeutic intent.

Tumor-targeting cell-penetrating peptide, p28, for glioblastoma imaging and therapy

Despite recent advances in cancer research, glioblastoma multiforme (GBM) remains a highly aggressive brain tumor as its treatment options are limited. The current standard treatment includes surgery followed by radiotherapy and adjuvant chemotherapy. However, surgery without image guidance is often challenging to achieve maximal safe resection as it is difficult to precisely discern the lesion to be removed from surrounding brain tissue. In addition, the efficacy of adjuvant chemotherapy is limited by poor penetration of therapeutics through the blood-brain barrier (BBB) into brain tissues, and the lack of tumor targeting. In this regard, we utilized a tumor-targeting cell-penetration peptide, p28, as a therapeutic agent to improve the efficacy of a current chemotherapeutic agent for GBM, and as a carrier for a fluorescence imaging agent for a clear identification of GBM. Here, we show that a near-infrared (NIR) imaging agent, ICG-p28 (a chemical conjugate of an FDA-approved NIR dye, indocyanine green ICG, and tumor-targeting p28 peptide) can preferentially localize tumors in multiple GBM animal models. Moreover, xenograft studies show that p28, as a therapeutic agent, can enhance the cytotoxic activity of temozolomide (TMZ), one of the few effective drugs for brain tumors. Collectively, our findings highlight the important role of the tumor-targeting peptide, which has great potential for intraoperative image-guided surgery and the development of new therapeutic strategies for GBM.

A Novel Defined Endoplasmic Reticulum Stress-Related lncRNA Signature for Prognosis Prediction and Immune Therapy in Glioma

Gliomas are a group of the most aggressive primary central nervous system tumors with limited treatment options. The abnormal expression of long non-coding RNA (lncRNA) is related to the prognosis of glioma. However, the role of endoplasmic reticulum (ER) stress-associated lncRNAs in glioma prognosis has not been reported. In this paper, we obtained ER stress-related lncRNAs by co-expression analysis, and then a risk signature composed of 6 ER stress-related lncRNAs was constructed using Cox regression analysis. Glioma samples in The Cancer Genome Atlas (TCGA) were separated into high- and low-risk groups based on the median risk score. Compared with the low-risk group, patients in the high-risk group had shorter survival times. Additionally, we verified the predictive ability of these candidate lncRNAs in the testing set. Three glioma patient subgroups (cluster 1/2/3) were identified by consensus clustering. We further analysed the abundance of immune-infiltrating cells and the expression levels of immune checkpoint molecules in both three subgroups and two risk groups, respectively. Immunotherapy and anticancer drug response prediction showed that ER stress-related lncRNA risk signature positively correlates with responding to immune checkpoints and chemosensitivity. Functional analysis showed that these gene sets are enriched in the malignant process of tumors. Finally, LINC00519 was chosen for functional experiments. The silence of LINC00519 restrained the migration and invasion of glioma cells. Hence, those results indicated that ER stress-related lncRNA risk signature could be a potential treatment target and a prognosis biomarker for glioma patients.

Glioblastoma, from disease understanding towards optimal cell-based in vitro models

Background

Glioblastoma (GBM) patients are notoriously difficult to treat and ultimately all succumb to disease. This unfortunate scenario motivates research into better characterizing and understanding this disease, and into developing novel research tools by which potential novel therapeutics and treatment options initially can be evaluated pre-clinically. Here, we provide a concise overview of glioblastoma epidemiology, disease classification, the challenges faced in the treatment of glioblastoma and current novel treatment strategies. From this, we lead into a description and assessment of advanced cell-based models that aim to narrow the gap between pre-clinical and clinical studies. Such invitro models are required to deliver reliable and meaningful data for the development and pre-validation of novel therapeutics and treatments.

Conclusions

The toolbox for GBM cell-based models has expanded substantially, with the possibility of 3D printing tumour tissues and thereby replicating invivo tissue architectures now looming on the horizon. A comparison of experimental cell-based model systems and techniques highlights advantages and drawbacks of the various tools available, based on which cell-based models and experimental approaches best suited to address a diversity of research questions in the glioblastoma research field can be selected.

Anti-glioblastoma effects of phenolic variants of benzoylphenoxyacetamide (BPA) with high potential for blood brain barrier penetration

Glioblastomas are the most aggressive brain tumors for which therapeutic options are limited. Current therapies against glioblastoma include surgical resection, followed by radiotherapy plus concomitant treatment and maintenance with temozolomide (TMZ), however, these standard therapies are often ineffective, and average survival time for glioblastoma patients is between 12 and 18 months. We have previously reported a strong anti-glioblastoma activity of several metabolic compounds, which were synthetized based compounds, which were synthetized based on the chemical structure of a common lipid-lowering drug, fenofibrate, and share a general molecular skeleton of benzoylphenoxyacetamide (BPA). Extensive computational analyses of phenol and naphthol moieties added to the BPA skeleton were performed in this study with the objective of selecting new BPA variants for subsequent compound preparation and anti-glioblastoma testing. Initially, 81 structural variations were considered and their physical properties such as solubility (logS), blood–brain partitioning (logBB), and probability of entering the CNS calculated by the Central Nervous System—Multiparameter Optimization (MPO-CNS) algorithm were evaluated. From this initial list, 18 compounds were further evaluated for anti-glioblastoma activity in vitro. Nine compounds demonstrated desirable glioblastoma cell toxicity in cell culture, and two of them, HR51, and HR59 demonstrated significantly improved capability of crossing the model blood–brain-barrier (BBB) composed of endothelial cells, astrocytes and pericytes.

ATX-101, a Peptide Targeting PCNA, Has Antitumor Efficacy Alone or in Combination with Radiotherapy in Murine Models of Human Glioblastoma

Cell proliferation requires the orchestrated actions of a myriad of proteins regulating DNA replication, DNA repair and damage tolerance, and cell cycle. Proliferating cell nuclear antigen (PCNA) is a master regulator which interacts with multiple proteins functioning in these processes, and this makes PCNA an attractive target in anticancer therapies. Here, we show that a cell-penetrating peptide containing the AlkB homolog 2 PCNA-interacting motif (APIM), ATX-101, has antitumor activity in a panel of human glioblastoma multiforme (GBM) cell lines and patient-derived glioma-initiating cells (GICs). Their sensitivity to ATX-101 was not related to cellular levels of PCNA, or p53, PTEN, or MGMT status. However, ATX-101 reduced Akt/mTOR and DNA-PKcs signaling, and a correlation between high Akt activation and sensitivity for ATX-101 was found. ATX-101 increased the levels of γH2AX, DNA fragmentation, and apoptosis when combined with radiotherapy (RT). In line with the in vitro results, ATX-101 strongly reduced tumor growth in two subcutaneous xenografts and two orthotopic GBM models, both as a single agent and in combination with RT. The ability of ATX-101 to sensitize cells to RT is promising for further development of this compound for use in GBM.

Tumor heterogeneity assessed by sequencing and fluorescence in situ hybridization (FISH) data

MOTIVATION

Computational reconstruction of clonal evolution in cancers has become a crucial tool for understanding how tumors initiate and progress and how this process varies across patients. The field still struggles, however, with special challenges of applying phylogenetic methods to cancers, such as the prevalence and importance of copy number alteration (CNA) and structural variation (SV) events in tumor evolution, which are difficult to profile accurately by prevailing sequencing methods in such a way that subsequent reconstruction by phylogenetic inference algorithms is accurate.

RESULTS

In the present work, we develop computational methods to combine sequencing with multiplex interphase fluorescence in situ hybridization (miFISH) to exploit the complementary advantages of each technology in inferring accurate models of clonal CNA evolution accounting for both focal changes and aneuploidy at whole-genome scales. By integrating such information in an integer linear programming (ILP) framework, we demonstrate on simulated data that incorporation of FISH data substantially improves accurate inference of focal CNA and ploidy changes in clonal evolution from deconvolving bulk sequence data. Analysis of real glioblastoma data for which FISH, bulk sequence, and single cell sequence are all available confirms the power of FISH to enhance accurate reconstruction of clonal copy number evolution in conjunction with bulk and optionally single-cell sequence data.

AVAILABILITY

Source code is available on Github at https://github.com/CMUSchwartzLab/FISH_deconvolution.

Control of brain tumor growth by reactivating myeloid cells with niacin

Glioblastomas are generally incurable partly because monocytes, macrophages, and microglia in afflicted patients do not function in an antitumor capacity. Medications that reactivate these macrophages/microglia, as well as circulating monocytes that become macrophages, could thus be useful to treat glioblastoma. We have discovered that niacin (vitamin B3) is a potential stimulator of these inefficient myeloid cells. Niacin-exposed monocytes attenuated the growth of brain tumor-initiating cells (BTICs) derived from glioblastoma patients by producing anti-proliferative interferon-α14. Niacin treatment of mice bearing intracranial BTICs increased macrophage/microglia representation within the tumor, reduced tumor size, and prolonged survival. These therapeutic outcomes were negated in mice depleted of circulating monocytes or harboring interferon-α receptor-deleted BTICs. Combination treatment with temozolomide enhanced niacin-promoted survival. Monocytes from glioblastoma patients had increased interferon-α14 upon niacin exposure and were reactivated to reduce BTIC growth in culture. We highlight niacin, a common vitamin that can be quickly translated into clinical application, as an immune stimulator against glioblastomas.

Hallmarks of tumor-associated microglia response to experimental U87 human glioblastoma xenograft

Glioblastoma (GBM) is one of the deadliest primary brain neoplasm, heavily infiltrated with tumor-associated microglia/macrophages (TAM), which has received a great deal of interest. Bearing in mind that the number of peripheral macrophages by the 14^th day is negligible, in our study TAM were referred to as microglia. Here we evaluated histopathological characterization of TAM and kinetics of their infiltration in U87 orthotopic GBM, a commonly used model in preclinical research. To mimic different stages of GBM growth, we evaluated three-time points. Our data showed that the highest areal density of TAM was 7 days after GBM inoculation, with ability to proliferate early after initiation of GBM growth. The areal density of TAM within the tumor correlated with GBM growth and proliferation. Moreover, microglia underwent substantial morphological changes upon exposure to GBM cells. A transition from ramified morphology in peritumoral area to ameboid shape with larger soma and shortened, thick branches in the tumor core was observed. Higher areal fraction of blood vessels also correlated with the areal density of TAM. Given these pro-invasive features of microglia, this GBM model represents a good basis for further testing microglia as a target and new strategy to fight with.

In silico identification of the prognostic biomarkers and therapeutic targets associated with cancer stem cell characteristics of glioma

Glioma is the common histological subtype of malignancy in central nervous system, with a high morbidity and mortality. Cancer stem cells (CSCs) play an important role in regulating the tumorigenesis and progression of glioma; however, the prognostic biomarkers and therapeutic targets associated with CSC characteristics have not been fully acknowledged in glioma. In order to identify the prognostic stemness-related genes (SRGs) of glioma in silico, the RNA sequencing data of patients with glioma were retrieved from The Cancer Genome Atlas (TCGA) databases. The mRNA expression-based stemness index (mRNAsi) was significantly associated with the glioma histologic grade, isocitrate dehydrogenase 1 (IDH1) mutation and overall survival of glioma patients by the nonparametric test and Kaplan–Meier survival analysis. A total of 340 SRGs were identified as the overlapped stemness-related differential expressed genes (DEGs) of different histologic grade screened by the univariate Cox analysis. Based on 11 prognostic SRGs, the predict nomogram was constructed with the AUC of 0.832. Moreover, the risk score of the nomogram was an independent prognostic factor, indicating its significant applicability. Besides other eight reported biomarkers of glioma, we found that F2RL2, CLCNKA and LOXL4 were first identified as prognostic biomarkers for glioma. In conclusion, this bioinformatics study demonstrates the mRNAsi as a reliable index for the IDH1 mutation, histologic grade and OS of glioma patients and provides a well-applied model for predicting the OS for patients with glioma based on prognostic SRGs. Additionally, this in silico study also identifies three novel prognostic biomarkers (F2RL2, CLCNKA and LOXL4) for glioma patients.

Overlapping variants in the blood, tissues and cell lines for patients with intracranial meningiomas are predominant in stem cell-related genes

OBJECTIVE

Bulk tissue genomic analysis of meningiomas identified common somatic mutations, however, it often excluded blood-related variants. In contrast, genomic characterisation of primary cell lines that can provide critical information regarding growth and proliferation, have been rare. In our work, we identified the variants that are present in the blood, tissues and corresponding cell lines that are likely to be predictive, tumorigenic and progressive.

METHOD

Whole-exome sequencing was used to identify variants and distinguish related pathways that exist in 42 blood, tissues and corresponding cell lines (BTCs) samples for patients with intracranial meningiomas. Conventional sequencing was used for the confirmation of variants. Integrative analysis of the gene expression for the corresponding samples was utilised for further interpretations.

RESULTS

In total, 926 BTC variants were detected, implicating 845 genes. A pathway analysis of all BTC genes with damaging variants indicated the 'cell morphogenesis involved in differentiation' stem cell-related pathway to be the most frequently affected pathway. Concordantly, five stem cell-related genes, GPRIN2, ALDH3B2, ASPN, THSD7A and SIGLEC6, showed BTC variants in at least five of the patients. Variants that were heterozygous in the blood and homozygous in the tissues or the corresponding cell lines were rare (average: 1.3 ± 0.3%), and included variants in the RUNX2 and CCDC114 genes. An analysis comparing the variants detected only in tumours with aggressive features indicated a total of 240 BTC genes, implicating the 'homophilic cell adhesion via plasma membrane adhesion molecules' pathway, and identifying the stem cell-related transcription coactivator NCOA3/AIB1/SRC3 as the most frequent BTC gene. Further analysis of the possible impact of the poly-Q mutation present in the NCOA3 gene indicated associated deregulation of 15 genes, including the up-regulation of the stem cell related SEMA3D gene and the angiogenesis related VEGFA gene.

CONCLUSION

Stem cell-related pathways and genes showed high prevalence in the BTC variants, and novel variants in stem cell-related genes were identified for meningioma. These variants can potentially be used as predictive, tumorigenic and progressive biomarkers for meningioma.

Glioma-derived IL-33 orchestrates an inflammatory brain tumor microenvironment that accelerates glioma progression

Despite a deeper molecular understanding, human glioblastoma remains one of the most treatment refractory and fatal cancers. It is known that the presence of macrophages and microglia impact glioblastoma tumorigenesis and prevent durable response. Herein we identify the dual function cytokine IL-33 as an orchestrator of the glioblastoma microenvironment that contributes to tumorigenesis. We find that IL-33 expression in a large subset of human glioma specimens and murine models correlates with increased tumor-associated macrophages/monocytes/microglia. In addition, nuclear and secreted functions of IL-33 regulate chemokines that collectively recruit and activate circulating and resident innate immune cells creating a pro-tumorigenic environment. Conversely, loss of nuclear IL-33 cripples recruitment, dramatically suppresses glioma growth, and increases survival. Our data supports the paradigm that recruitment and activation of immune cells, when instructed appropriately, offer a therapeutic strategy that switches the focus from the cancer cell alone to one that includes the normal host environment.

Prognosis Analysis and Validation of m6A Signature and Tumor Immune Microenvironment in Glioma

Glioma is one of the most typical intracranial tumors, comprising about 80% of all brain malignancies. Several key molecular signatures have emerged as prognostic biomarkers, which indicate room for improvement in the current approach to glioma classification. In order to construct a more veracious prediction model and identify the potential prognosis-biomarker, we explore the differential expressed m⁶A RNA methylation regulators in 665 gliomas from TCGA-GBM and TCGA-LGG. Consensus clustering was applied to the m6A RNA methylation regulators, and two glioma subgroups were identified with a poorer prognosis and a higher grade of WHO classification in cluster 1. The further chi-squared test indicated that the immune infiltration was significantly enriched in cluster 1, indicating a close relation between m⁶A regulators and immune infiltration. In order to explore the potential biomarkers, the weighted gene co-expression network analysis (WGCNA), along with Least absolute shrinkage and selection operator (LASSO), between high/low immune infiltration and m⁶A cluster 1/2 groups were utilized for the hub genes, and four genes (TAGLN2, PDPN, TIMP1, EMP3) were identified as prognostic biomarkers. Besides, a prognostic model was constructed based on the four genes with a good prediction and applicability for the overall survival (OS) of glioma patients (the area under the curve of ROC achieved 0.80 (0.76-0.83) and 0.72 (0.68-0.76) in TCGA and Chinese Glioma Genome Atlas (CGGA), respectively). Moreover, we also found PDPN and TIMP1 were highly expressed in high-grade glioma from The Human Protein Atlas database and both of them were correlated with m6A and immune cell marker in glioma tissue samples. In conclusion, we construct a novel prognostic model which provides new insights into glioma prognosis. The PDPN and TIMP1 may serve as potential biomarkers for prognosis of glioma.

Tumor Copy Number Deconvolution Integrating Bulk and Single-Cell Sequencing Data

Characterizing intratumor heterogeneity (ITH) is crucial to understanding cancer development, but it is hampered by limits of available data sources. Bulk DNA sequencing is the most common technology to assess ITH, but involves the analysis of a mixture of many genetically distinct cells in each sample, which must then be computationally deconvolved. Single-cell sequencing is a promising alternative, but its limitations-for example, high noise, difficulty scaling to large populations, technical artifacts, and large data sets-have so far made it impractical for studying cohorts of sufficient size to identify statistically robust features of tumor evolution. We have developed strategies for deconvolution and tumor phylogenetics combining limited amounts of bulk and single-cell data to gain some advantages of single-cell resolution with much lower cost, with specific focus on deconvolving genomic copy number data. We developed a mixed membership model for clonal deconvolution via non-negative matrix factorization balancing deconvolution quality with similarity to single-cell samples via an associated efficient coordinate descent algorithm. We then improve on that algorithm by integrating deconvolution with clonal phylogeny inference, using a mixed integer linear programming model to incorporate a minimum evolution phylogenetic tree cost in the problem objective. We demonstrate the effectiveness of these methods on semisimulated data of known ground truth, showing improved deconvolution accuracy relative to bulk data alone.

Genomic and Phenotypic Characterization of a Broad Panel of Patient-Derived Xenografts Reflects the Diversity of Glioblastoma

PURPOSE

Glioblastoma is the most frequent and lethal primary brain tumor. Development of novel therapies relies on the availability of relevant preclinical models. We have established a panel of 96 glioblastoma patient-derived xenografts (PDX) and undertaken its genomic and phenotypic characterization.

EXPERIMENTAL DESIGN

PDXs were established from glioblastoma, IDH-wildtype (n = 93), glioblastoma, IDH-mutant (n = 2), diffuse midline glioma, H3 K27M-mutant (n = 1), and both primary (n = 60) and recurrent (n = 34) tumors. Tumor growth rates, histopathology, and treatment response were characterized. Integrated molecular profiling was performed by whole-exome sequencing (WES, n = 83), RNA-sequencing (n = 68), and genome-wide methylation profiling (n = 76). WES data from 24 patient tumors was compared with derivative models.

RESULTS

PDXs recapitulate many key phenotypic and molecular features of patient tumors. Orthotopic PDXs show characteristic tumor morphology and invasion patterns, but largely lack microvascular proliferation and necrosis. PDXs capture common and rare molecular drivers, including alterations of TERT, EGFR, PTEN, TP53, BRAF, and IDH1, most at frequencies comparable with human glioblastoma. However, PDGFRA amplification was absent. RNA-sequencing and genome-wide methylation profiling demonstrated broad representation of glioblastoma molecular subtypes. MGMT promoter methylation correlated with increased survival in response to temozolomide. WES of 24 matched patient tumors showed preservation of most genetic driver alterations, including EGFR amplification. However, in four patient-PDX pairs, driver alterations were gained or lost on engraftment, consistent with clonal selection.

CONCLUSIONS

Our PDX panel captures the molecular heterogeneity of glioblastoma and recapitulates many salient genetic and phenotypic features. All models and genomic data are openly available to investigators.

EGFR blockade in GBM brain tumor stem cells synergizes with JAK2/STAT3 pathway inhibition to abrogate compensatory mechanisms in vitro and in vivo

Background

The EGFR pathway is frequently mutated in glioblastoma (GBM). However, to date, EGFR therapies have not demonstrated efficacy in clinical trials. Poor brain penetration of conventional inhibitors, lack of patient stratification for EGFR status, and mechanisms of resistance are likely responsible for the failure of EGFR-targeted therapy. We aimed to address these elements in a large panel of molecularly diverse patient-derived GBM brain tumor stem cells (BTSCs).

Methods

In vitro growth inhibition and on-target efficacy of afatinib, pacritinib, or a combination were assessed by cell viability, neurosphere formation, cytotoxicity, limiting dilution assays, and western blotting. In vivo efficacy was assessed with mass spectrometry, immunohistochemistry, magnetic resonance imaging, and intracranial xenograft models.

Results

We show that afatinib and pacritinib decreased BTSC growth and sphere-forming capacity in vitro. Combinations of the 2 drugs were synergistic and abrogated the activation of STAT3 signaling observed upon EGFR inhibition in vitro and in vivo. We further demonstrate that the brain-penetrant EGFR inhibitor, afatinib, improved survival in EGFRvIII mt orthotopic xenograft models. However, upregulation of the oncogenic STAT3 signaling pathway was observed following afatinib treatment. Combined inhibition with 2 clinically relevant drugs, afatinib and pacritinib, synergistically decreased BTSC viability and abrogated this compensatory mechanism of resistance to EGFR inhibition. A significant decrease in tumor burden in vivo was observed with the combinatorial treatment.

Conclusions

These data demonstrate that brain-penetrant combinatorial therapies targeting the EGFR and STAT3 signaling pathways hold therapeutic promise for GBM.

Comprehensive genomic profiling of glioblastoma tumors, BTICs, and xenografts reveals stability and adaptation to growth environments

Glioblastoma multiforme (GBM) is the most deadly brain tumor, and currently lacks effective treatment options. Brain tumor-initiating cells (BTICs) and orthotopic xenografts are widely used in investigating GBM biology and new therapies for this aggressive disease. However, the genomic characteristics and molecular resemblance of these models to GBM tumors remain undetermined. We used massively parallel sequencing technology to decode the genomes and transcriptomes of BTICs and xenografts and their matched tumors in order to delineate the potential impacts of the distinct growth environments. Using data generated from whole-genome sequencing of 201 samples and RNA sequencing of 118 samples, we show that BTICs and xenografts resemble their parental tumor at the genomic level but differ at the mRNA expression and epigenomic levels, likely due to the different growth environment for each sample type. These findings suggest that a comprehensive genomic understanding of in vitro and in vivo GBM model systems is crucial for interpreting data from drug screens, and can help control for biases introduced by cell-culture conditions and the microenvironment in mouse models. We also found that lack of MGMT expression in pretreated GBM is linked to hypermutation, which in turn contributes to increased genomic heterogeneity and requires new strategies for GBM treatment.

Intertumoral heterogeneity in patient-specific drug sensitivities in treatment-naïve glioblastoma

Background

A major barrier to effective treatment of glioblastoma (GBM) is the large intertumoral heterogeneity at the genetic and cellular level. In early phase clinical trials, patient heterogeneity in response to therapy is commonly observed; however, how tumor heterogeneity is reflected in individual drug sensitivities in the treatment-naïve glioblastoma stem cells (GSC) is unclear.

Methods

We cultured 12 patient-derived primary GBMs as tumorspheres and validated tumor stem cell properties by functional assays. Using automated high-throughput screening (HTS), we evaluated sensitivity to 461 anticancer drugs in a collection covering most FDA-approved anticancer drugs and investigational compounds with a broad range of molecular targets. Statistical analyses were performed using one-way ANOVA and Spearman correlation.

Results

Although tumor stem cell properties were confirmed in GSC cultures, their in vitro and in vivo morphology and behavior displayed considerable tumor-to-tumor variability. Drug screening revealed significant differences in the sensitivity to anticancer drugs (p < 0.0001). The patient-specific vulnerabilities to anticancer drugs displayed a heterogeneous pattern. They represented a variety of mechanistic drug classes, including apoptotic modulators, conventional chemotherapies, and inhibitors of histone deacetylases, heat shock proteins, proteasomes and different kinases. However, the individual GSC cultures displayed high biological consistency in drug sensitivity patterns within a class of drugs. An independent laboratory confirmed individual drug responses.

Conclusions

This study demonstrates that patient-derived and treatment-naïve GSC cultures maintain patient-specific traits and display intertumoral heterogeneity in drug sensitivity to anticancer drugs. The heterogeneity in patient-specific drug responses highlights the difficulty in applying targeted treatment strategies at the population level to GBM patients. However, HTS can be applied to uncover patient-specific drug sensitivities for functional precision medicine.

Electronic supplementary material

The online version of this article (10.1186/s12885-019-5861-4) contains supplementary material, which is available to authorized users.

The efficacy of a coordinated pharmacological blockade in glioblastoma stem cells with nine repurposed drugs using the CUSP9 strategy

PURPOSE

Constructed from a theoretical framework, the coordinated undermining of survival paths in glioblastoma (GBM) is a combination of nine drugs approved for non-oncological indications (CUSP9; aprepitant, auranofin, captopril, celecoxib, disulfiram, itraconazole, minocycline, quetiapine, and sertraline) combined with temozolomide (TMZ). The availability of these drugs outside of specialized treatment centers has led patients to embark on combination treatments without systematic follow-up. However, no experimental data on efficacy using the CUSP9 strategy in GBM have been reported.

METHODS

Using patient-derived glioblastoma stem cell (GSC) cultures from 15 GBM patients, we described stem cell properties of individual cultures, determined the dose-response relationships of the drugs in the CUSP9, and assessed the efficacy the CUSP9 combination with TMZ in concentrations clinically achievable. The efficacy was evaluated by cell viability, cytotoxicity, and sphere-forming assays in both primary and recurrent GSC cultures.

RESULTS

We found that CUSP9 with TMZ induced a combination effect compared to the drugs individually (p < 0.0001). Evaluated by cell viability and cytotoxicity, 50% of the GSC cultures displayed a high sensitivity to the drug combination. In clinical plasma concentrations, the effect of the CUSP9 with TMZ was superior to TMZ monotherapy (p < 0.001). The Wnt-signaling pathway has been shown important in GSC, and CUSP9 significantly reduces Wnt-activity.

CONCLUSIONS

Adding experimental data to the theoretical rationale of CUSP9, our results demonstrate that the CUSP9 treatment strategy can induce a combination effect in both treatment-naïve and pretreated GSC cultures; however, predicting response in individual cultures will require further profiling of GSCs.

The Genetic Landscape of Human Glioblastoma and Matched Primary Cancer Stem Cells Reveals Intratumour Similarity and Intertumour Heterogeneity

Glioblastoma (GBM) is the most malignant human brain tumour, characterized by rapid progression, invasion, intense angiogenesis, high genomic instability, and resistance to therapies. Despite countless experimental researches for new therapeutic strategies and promising clinical trials, the prognosis remains extremely poor, with a mean survival of less than 14 months. GBM aggressive behaviour is due to a subpopulation of tumourigenic stem-like cells, GBM stem cells (GSCs), which hierarchically drive onset, proliferation, and tumour recurrence. The morbidity and mortality of this disease strongly encourage exploring genetic characteristics of GSCs. Here, using array-CGH platform, we investigated genetic and genomic aberration profiles of GBM parent tumour (n = 10) and their primarily derived GSCs. Statistical analysis was performed by using R software and complex heatmap and corrplot packages. Pearson correlation and K-means algorithm were exploited to compare genetic alterations and to group similar genetic profiles in matched pairs of GBM and derived GSCs. We identified, in both GBM and matched GSCs, recurrent copy number alterations, as chromosome 7 polysomy, chromosome 10 monosomy, and chromosome 9p21deletions, which are typical features of primary GBM, essential for gliomagenesis. These observations suggest a condition of strong genomic instability both in GBM as GSCs. Our findings showed the robust similarity between GBM mass and GSCs (Pearson corr.≥0.65) but also highlighted a marked variability among different patients. Indeed, the heatmap reporting Gain/Loss State for 21022 coding/noncoding genes demonstrated high interpatient divergence. Furthermore, K-means algorithm identified an impairment of pathways related to the development and progression of cancer, such as angiogenesis, as well as pathways related to the immune system regulation, such as T cell activation. Our data confirmed the preservation of the genomic landscape from tumour tissue to GSCs, supporting the relevance of this cellular model to test in vitro new target therapies for GBM.

Genetically distinct glioma stem-like cell xenografts established from paired glioblastoma samples harvested before and after molecularly targeted therapy

Intratumoural heterogeneity underlies tumour escape from molecularly targeted therapy in glioblastoma. A cell-based model preserving the evolving molecular profiles of a tumour during treatment is key to understanding the recurrence mechanisms and development of strategies to overcome resistance. In this study, we established a matched pair of glioblastoma stem-like cell (GSC) cultures from patient glioblastoma samples before and after epidermal growth factor receptor (EGFR)-targeted therapy. A patient with recurrent glioblastoma (MGG70R) harboring focal, high-level EGFR amplification received the irreversible EGFR tyrosine kinase inhibitor dacomitinib. The tumour that subsequently recurred (MGG70RR) showed diploid EGFR, suggesting inhibitor-mediated elimination of EGFR-amplified tumour cells and propagation of EGFR non-amplified cell subpopulations. The MGG70R-GSC line established from MGG70R formed xenografts retaining EGFR amplification and EGFR overexpression, while MGG70RR-GSC established from MGG70RR generated tumours that lacked EGFR amplification and EGFR overexpression. MGG70R-GSC-derived intracranial xenografts were more proliferative than MGG70RR-GSC xenografts, which had upregulated mesenchymal markers, mirroring the pathological observation in the corresponding patient tumours. In vitro MGG70R-GSC was more sensitive to EGFR inhibitors than MGG70RR-GSC. Thus, these molecularly distinct GSC lines recapitulated the subpopulation alteration that occurred during glioblastoma evasion of targeted therapy, and offer a valuable model facilitating therapeutic development for recurrent glioblastoma.

The first-in-class alkylating deacetylase inhibitor molecule tinostamustine shows antitumor effects and is synergistic with radiotherapy in preclinical models of glioblastoma

Background

The use of alkylating agents such as temozolomide in association with radiotherapy (RT) is the therapeutic standard of glioblastoma (GBM). This regimen modestly prolongs overall survival, also if, in light of the still dismal prognosis, further improvements are desperately needed, especially in the patients with O6-methylguanine-DNA-methyltransferase (MGMT) unmethylated tumors, in which the benefit of standard treatment is less. Tinostamustine (EDO-S101) is a first-in-class alkylating deacetylase inhibitor (AK-DACi) molecule that fuses the DNA damaging effect of bendamustine with the fully functional pan-histone deacetylase (HDAC) inhibitor, vorinostat, in a completely new chemical entity.

Methods

Tinostamustine has been tested in models of GBM by using 13 GBM cell lines and seven patient-derived GBM proliferating/stem cell lines in vitro. U87MG and U251MG (MGMT negative), as well as T98G (MGMT positive), were subcutaneously injected in nude mice, whereas luciferase positive U251MG cells and patient-derived GBM stem cell line (CSCs-5) were evaluated the orthotopic intra-brain in vivo experiments.

Results

We demonstrated that tinostamustine possesses stronger antiproliferative and pro-apoptotic effects than those observed for vorinostat and bendamustine alone and similar to their combination and irrespective of MGMT expression. In addition, we observed a stronger radio-sensitization of single treatment and temozolomide used as control due to reduced expression and increased time of disappearance of γH2AX indicative of reduced signal and DNA repair. This was associated with higher caspase-3 activation and reduction of RT-mediated autophagy. In vivo, tinostamustine increased time-to-progression (TTP) and this was additive/synergistic to RT. Tinostamustine had significant therapeutic activity with suppression of tumor growth and prolongation of DFS (disease-free survival) and OS (overall survival) in orthotopic intra-brain models that was superior to bendamustine, RT and temozolomide and showing stronger radio sensitivity.

Conclusions

Our data suggest that tinostamustine deserves further investigation in patients with glioblastoma.

Electronic supplementary material

The online version of this article (10.1186/s13045-018-0576-6) contains supplementary material, which is available to authorized users.

Multi-omics analysis of primary glioblastoma cell lines shows recapitulation of pivotal molecular features of parental tumors

Background

Glioblastoma (GBM) is the deadliest primary brain cancer in adults. Emerging innovative therapies hold promise for personalized cancer treatment. Improving therapeutic options depends on research relying on relevant preclinical models. In this line we have established in the setting of the GlioTex project (GBM and Experimental Therapeutics) a GBM patient-derived cell line (GBM-PDCL) library. A multi-omic approach was used to determine the molecular landscape of PDCL and the extent to which they represent GBM tumors.

Methods

Single nucleotide polymorphism array, expression arrays, exome sequencing, and RNA sequencing were used to measure and compare the molecular landscapes of 20 samples representing 10 human GBM tumors and paired GBM-PDCLs.

Results

Copy number variations were similar for a median of 85% of the genome and for 59% of the major focal events. Somatic point mutations were similar in a median of 41%. Mutations in GBM driver and "druggable" genes were maintained in 67% of events. Mutations that were not conserved in the PDCL were mainly low allelic fraction and/or non-driver mutations. Based on RNA expression profiling, PDCLs cluster closely to their parental tumor with overexpression of pathways associated with cancer progression in PDCL.

Conclusions

Overall, PDCLs recapitulate pivotal molecular alterations of paired-parental tumors supporting their use as a preclinical model of GBM. However, some driver aberrations are lost or gained in the passage from tumor to PDCL. Our results support using PDCL as a relevant preclinical model of GBM. Further investigations of changes between PDCLs and their parental tumor may provide insights into GBM biology.

Glioblastoma-associated microglia and macrophages: targets for therapies to improve prognosis

Glioblastoma is the most common and most malignant primary adult human brain tumour. Diagnosis of glioblastoma carries a dismal prognosis. Treatment resistance and tumour recurrence are the result of both cancer cell proliferation and their interaction with the tumour microenvironment. A large proportion of the tumour microenvironment consists of an inflammatory infiltrate predominated by microglia and macrophages, which are thought to be subverted by glioblastoma cells for tumour growth. Thus, glioblastoma-associated microglia and macrophages are logical therapeutic targets. Their emerging roles in glioblastoma progression are reflected in the burgeoning research into therapeutics directed at their modification or elimination. Here, we review the biology of glioblastoma-associated microglia and macrophages, and model systems used to study these cells in vitro and in vivo. We discuss translation of results using these model systems and review recent advances in immunotherapies targeting microglia and macrophages in glioblastoma. Significant challenges remain but medications that affect glioblastoma-associated microglia and macrophages hold considerable promise to improve the prognosis for patients with this disease.