Different states of stemness of glioblastoma stem cells sustain glioblastoma subtypes indicating novel clinical biomarkers and high-efficacy customized therapies
Journal of Experimental & Clinical Cancer Research volume 42, Article number: 244 (2023)
Glioblastoma (GBM) is the most malignant among gliomas with an inevitable lethal outcome. The elucidation of the physiology and regulation of this tumor is mandatory to unravel novel target and effective therapeutics. Emerging concepts show that the minor subset of glioblastoma stem cells (GSCs) accounts for tumorigenicity, representing the true target for innovative therapies in GBM.
Here, we isolated and established functionally stable and steadily expanding GSCs lines from a large cohort of GBM patients. The molecular, functional and antigenic landscape of GBM tissues and their derivative GSCs was highlited in a side-by-side comprehensive genomic and transcriptomic characterization by ANOVA and Fisher’s exact tests. GSCs’ physio-pathological hallmarks were delineated by comparing over time in vitro and in vivo their expansion, self-renewal and tumorigenic ability with hierarchical linear models for repeated measurements and Kaplan–Meier method. Candidate biomarkers performance in discriminating GBM patients’ classification emerged by classification tree and patients’ survival analysis.
Here, distinct biomarker signatures together with aberrant functional programs were shown to stratify GBM patients as well as their sibling GSCs population into TCGA clusters. Of importance, GSCs cells were demonstrated to fully resemble over time the molecular features of their patient of origin. Furthermore, we pointed out the existence of distinct GSCs subsets within GBM classification, inherently endowed with different self-renewal and tumorigenic potential. Particularly, classical GSCs were identified by more undifferentiated biological hallmarks, enhanced expansion and clonal capacity as compared to the more mature, relatively slow-propagating mesenchymal and proneural cells, likely endowed with a higher potential for infiltration either ex vivo or in vivo. Importantly, the combination of DCX and EGFR markers, selectively enriched among GSCs pools, almost exactly predicted GBM patients’ clusters together with their survival and drug response.
In this study we report that an inherent enrichment of distinct GSCs pools underpin the functional inter-cluster variances displayed by GBM patients. We uncover two selectively represented novel functional biomarkers capable of discriminating GBM patients’ stratification, survival and drug response, setting the stage for the determination of patient-tailored diagnostic and prognostic strategies and, mostly, for the design of appropriate, patient-selective treatment protocols.
Glioblastoma (GBM), which is classified as Isocitrate dehydrogenase (IDH1) wild-type (IDH1wt) according to 2021 WHO , is the most common and aggressive form among primary malignant brain tumors. Despite aggressive therapies and advances in genomic and molecular classifications, the five-year overall rate of survival after developing this cancer is still only 6,9% with a median survival of ≤ 15 months [2, 3]. The high mortality is mainly due to the rapidly recurs of this tumor and several factors contribute to it, such as a high intratumoral phenotypic heterogeneity and plasticity, as well as the feature of the hGBM cells to rapidly migrate and infiltrate within brain tissue, preventing a complete removal and inducing resistance to chemotherapy and radiotherapy. It follows that the residual cells develop alternative evolutionary paths that drive the growth of recurrent tumors and contribute to the treatment failure [4,5,6].
In order to identify the molecular processes underlying heterogeneity and plasticity of GBM cells, a bioinformatic analysis of the gene expression profile was initially performed by the Cancer Genome Atlas Consortium (TCGA), allowing to define the existence of four molecular clusters: proneural (PN), classical (CL), mesenchymal (MS) and neural (NE). The latter was excluded as it was contaminated with normal neural tissue [7, 8]. These subtypes of GBM identified by distinct molecular profiles are described to differ in their clinical courses and drug responses, being responsible for the failure of multimodal therapies, including RT, chemotherapy and other targeted therapies [7, 9,10,11,12,13,14,15,16,17].
Subsequently, Single-cell RNA-sequencing (scRNA-seq) studies have shown that several subtypes may exist simultaneously in different region of the same tumor and that molecular clusters can change over time and through therapy [8, 18, 19]. The combination of Single-cell seq and TCGA analysis, demonstrates that in GBM exist 4 cellular states that are reminiscent of canonical neurodevelopmental cell types: 1) neural-progenitor-like (NPC-like); 2) oligodendrocyte-progenitor-like (OPC-like); 3) astrocyte-like (AC-like), and 4) mesenchymal-like (MES-like). Each state or the combination of two of them is consistent with the three molecular subtypes previously described by TGCA. Actually, TGCA-CL and TCGA-MS correspond to AC-like and MES-like state, while TCGA-PN to OPC-like and NPC-like ones. As a result, each tumor turns out to be composed of cells that are in multiple cellular states and the tumor microenvironment could affect the transition among these different conditions or proliferation . Although the identification of molecular subtypes has been of extreme importance by itself in understanding the molecular heterogeneity of the GBM, the source of a functional heterogeneity remains unclear and has had almost no translational impact on the clinical context and development of specific therapies.
The search for innovative and more effective therapies for GBM has benefited from the discovery that a relatively rare type of cancer cell, which possesses the cardinal features of an aberrant neural stem cell (NSCs), lie at the root of GBM insurgency propagation and perpetuation in humans, the glioblastoma stem cells (GSCs). This subpopulation of cells, endowed with tumor-propagating ability, express embryonic or tissue stem cell genes and, consistent with its role in supporting the relapse after therapy, is reported to be resistant to existing standard therapies [18, 21,22,23,24,25,26,27]. Importantly, GSCs cells have also been recently described to play a key role in underlying the phenotypic diversity and plasticity of GBM, which is strongly influenced also by the microenvironment [27,28,29]. Within the complexity of intratumoral heterogeneity, Richards and colleagues outlined that, down to the transcriptional and genetic level, GSCs can be dynamically retrieved in a combination of two phenotypic states (conditions), reproducing both neurodevelopmental and inflammatory functional programs. The switch from one program to another can occur based on early somatic alteration and copy number variation (CNV) state of GSCs and, in part, by cytokine signaling . Patient-derived GSCs have also been described by different transcriptional profiles [8, 30, 31] and, importantly, to inherently display distinct level of stem cell markers depending on their GBM cluster of origin [32, 33]. Nevertheless, the intrinsic GSCs different subsets within the tumor itself and whether these differences might account for the functional inter-cluster variances in GBM subgroups are still under-investigated. Thus, GSCs remain a crucial elusive and difficult cellular target in glioblastoma treatment [6, 34].
By means of a standardized approach faithfully modeling GBM based on a cell system of stable, fully characterized GSCs lines, here we report that each transcriptional cluster of GBM patients contain different GSCs subpopulations endowed with distinctive molecular, functional and antigenic phenotypes, which are closely related to inherent different states of stemness. This allow the identification of new subtype-associated functional biomarkers predicting GBM patient stratification and survival, thus opening new and exciting possibilities for patient-tailored diagnostic and prognostic purposes, and, importantly, for the definition of individual patient-specific, drug-responsiveness therapeutic protocols by tackling different subset of GSCs.
Sample cohort features and population analysis
The sample cohort analyzed in this study was composed of 93 glioblastoma post-surgery tissues from patients who underwent neurosurgical resection at IRCCS National Neurologic Institute “C. Besta” and 34 GSCs lines established from them together with patients’ clinical data and matched peripheral blood (available from 27 patients). Clinical and sequencing information (GBM subtype, overall survival, IDH1, TERT promoter and EGFR status) are provided in Table 1. Tumor samples were collected in accordance to the ethical guidelines of the 2013 Declaration of Helsinki after approval of the institutional ethic board (protocol n°02 and protocol n°61) with signed informed consent. Material was anonymized at the time of collection. All samples were from patients with a confirmed diagnosis of GBM by a pathologist and classified according to the WHO guidelines. GSCs cells population, clonogenic and differentiation analyses were performed as in [22, 32]. The authenticity of each cell line was last checked in January 2021 by CNV profiling.
In vivo studies
All animal analyses were performed according to the Guidelines for the Care and Use of Laboratory Animals and experimental protocols approved by the Italian Ministry of Health (805/2016-PR and 119/2019-PR). Anesthesia and analgesics were used in order to minimize any suffering of the animals. Tumorigenic capacity of GSCs lines and their in vivo phenotypic hallmarks were examined by stereotactic injection of 3 × 10^5 firefly luciferase-transduced (F-luc) GSCs cells from different subtypes into the right striatum of SCID mice (Charles River Lab) as previously described [22, 32, 33, 35]. Mice were checked daily for any signs of distress and monitored weekly by ventral and dorsal views with In vivo Lumina (Xenogen, PerkinElmer Inc) and tumor growth was indirectly quantified. Mice were then sacrificed at different endpoints as soon as they became symptomatic according to the subtype of the GSCs injected, and their brains were collected after transcardiac perfusion and processed as previously described [22, 32, 33, 35, 36].
To evaluate the relationship between the level of EGFR and DCX and patients’ outcome, 117 IDH1 wild-type GBM patients were selected in the TCGA dataset  and mRNA expression data with corresponding clinical information downloaded from https://xenabrowser.net/datapages/. Optimal cutoff between high and low mRNA expression groups were determined through the R package “survminer”. 83 high-grade glioma patients of our cohort with available clinical follow up data were stratified into TCGA subtypes (n = 13 TCGA-PN, n = 28 TCGA-MS and n = 42 TCGA-CL) and their survival data analyzed with Kaplan–Meier plot by GraphPad Prism v.7.0 software. Overall comparisons were performed by Log-rank test. P-values < 0.05 were considered significant.
In vitro data were analyzed using R and GraphPad Prism v7.0 software with statistical test selected according to the variance and distribution of data. Gene expression analysis and differences among groups originated from microarray sequencing was performed by ANOVA test with a q-value < 0.05 considered significant. Sequencing and SNPs arrays data were examined by comparing each cluster to the others with a two-tailed Fisher’s exact test. Canonical pathway enrichment analysis on somatic variants was performed with Pathscore , which is based on a collection of ‘canonical pathways’ from the Molecular Signatures Database  including pathways from the KEGG, Biocarta, Reactome and Nature-NCI databases. Analysis of the biological functions of differently altered gene copy number and differently expressed mRNAs was performed by Ingenuity Pathway Analysis (IPA; Qiagen, http://www.ingenuity.com/) and R software [32, 33, 35, 40] with a right-tailed Fischer’s exact test. P-values were further adjusted and cutoff for significance was set as q-values < 0.05 and activating z-score threshold ≥ 2 or inhibiting z-score threshold ≤ 2.
To assess cell stability, CNV mean values at different time points were compared. According to their difference, we classified as amplification a difference larger than 0.5 and as deletion a difference smaller than 0.5. Therefore, we cross-tabulated the amplification/deletion status at two different time points in order to assess cell stability. As observed for each cell line, changes in copy number status (i.e. from amplification to deletion or from deletion to amplification) were observed for a very small and negligible number of genes (< 5%). To assess and compared trend over time, in vitro GSCs growth curves and in vivo GSCs-tumors growth rate across GBM clusters was analyzed with hierarchical linear models for repeated measurements by SAS Statistical Package Release 9.4 (SAS Institute) [41, 42]. As outcome, the log-transformed number of cell number. Unequally spaced time occasions within the experiments were counted by a spatial power correlation type . A 5% cutoff was applied to validate data significance. Overall survival of GSCs-implanted mice among GBM subtypes was examined with GraphPad Prism v7.0 software by Kaplan-Meyer curves based on survival algorithms determined by the log rank Mantel-Cox and Gehan-Breslow Wilcoxon tests. Groups were compared by respective median survival of number of days taken to reach 50% morbidity. P < 0.05 as significant.
Molecular hallmarks of GBM patients and their sibling GSCs lines across TCGA clusters
To get insight into the distinct aberrant molecular signature differentiating GBM subtypes [7, 20, 37, 43], high-grade glioma tissues, matched blood DNA and their derivative GSCs lines were sequenced and then analyzed for their expression program. As shown by the integrated matrix in Fig. 1 and Supplementary Fig. 1, we found that each GBM patient together with the cognate GSCs line displayed a distinctive landscape of somatic mutations and copy number alterations across GBM clusters. Mutation-calling analysis revealed a comparable number of somatic variants between GBM patients and their GSCs, whereas significant differences emerged when the three molecular clusters were compared to each other (**P = 0.005, Kruskal–Wallis test) (Supplementary Fig. 2A and Supplementary Tables 1 and 2). Analysis of molecular abnormalities’ distribution confirmed that, among other genes, TCGA-CL samples predominantly harbored somatic mutations in EGFR gene (***P < 0.001, two-sided Fisher’s exact test), also reporting its constitutive variant EGFRvIII, known to confer enhanced tumorigenic behavior and pharmacological-resistance [11, 27] (***P < 0.001) (Fig. 1, Supplementary Fig. 2B, C and Supplementary Table 2). Chromosome (chr) 7 amplification, with focal amplification at 7p11.2 (EGFR), paired with chr10 loss, which mostly co-occurred (46%) with the constitutively active variant EGFRvIII, was also highlighted together with significant gains at 3p12.1 and 4p15.33 (q value = 0.02) (RHOA, EPHB1 and BMPR1B, CAMK2D) (Fig. 2A and Supplementary Table 3). Consistently, classical cases were mostly highlighted by the highest level of EGFR, NES, AQP4, CDKN2A, BCAN, DCLK1, VCAN and DLL1 (***P < 0.001, **P < 0.01, *P < 0.05, Kruskal–Wallis test) (Fig. 2B and Supplementary Fig. 2D-F). PTEN and NF1 genes, usually reported as commuted, have been observed harboring the highest rate of mutation in mesenchymal cases, together with RB1 and PIK3CA (***P < 0.001, **P < 0.01). In the same cluster, focal amplifications at 21q21.3 (q value = 0.03) (TIAM1, ABCG1, OLIG1) and 2p24.3 (q value = 0.07) (MYCN) were mainly reported together with an enrichment for genes associated with immune/complement responses and inflammation (CD109, C1R, PAPPA, CTSB) and underlining an invasive “mesenchymal” signature (LOX, CHI3L1, MET, CD44). The highest rate of PDGFRA (*P < 0.05), EPHA2 and BMPR2 variations was retrieved in proneural cases, also mostly characterized by focal amplification at 4q12 (q value = 3.2881E-23) (PDGFRA), associated in the 29% of case with the amplification at 12q13.3, harboring CDK4 gene , and 7p14.3 (q value = 0.02) (DLX5, EPHB6, HOXA10). Proneural patients and their GSCs were also found to predominantly express developmental genes (SLITRK2, DLX1, PDGFRA, OLIG2, NKX2) and the highest level of the known stemness-associated genes (SOX11, DCX, HOXA7, EPHA2 and PROM1). Remarkably, considerable overlap in genetic aberrations between pairs of GBM patient’s tissue and GSCs line was observed (Supplementary Fig. 3A), being the latter also reported genetically stable over time (Supplementary Table 4).
Gene enrichment analysis on SNVs and CNVs [44, 45] and functional analysis performed on the global gene expression profiles of GBM patients and their GSCs outlined the existence of distinct, cluster-related, biological damaged profiles and enrichment for gene sets (Fig. 1 and Supplementary Fig. 3B). Particularly, genetic alterations of immune system, DNA replication and TNFR1 pathways occurred exclusively in mesenchymal cases, which also harbored aberrations involved in key biological processes as immune response, migration and invasion of cells and inflammatory response (P < 0.001, two-sided Fisher’s exact test, corrected for multiple testing using a Familywise Error Rate) (Fig. 1 and Supplementary Table 4). The same cluster was described by an expression pattern, which included genes regulating invasion, metastasis, angiogenesis and activation of immune system (Fig. 1, Supplementary Fig. 3B and Supplementary Tables 5 and 6). Conversely, calcium signaling, adherens junction and axon guidance pathways (P < 0.05), ERBB network signaling, necrosis and apoptosis functions were predominantly disrupted in classical patients, who also reported significant enhancements in the expression of selective signaling including proliferation and cycling and cell viability. Cellular homeostasis, differentiation of oligodendrocytes and neuroglia, activation of dendritic cells, proliferation of phagocytes and antigen presenting cells and permeability of vascular system functions were mainly targeted in proneural cases, identified by genes involved in CNS development and function, neuritogenesis, tissue and cellular development.
These findings underlined that peculiar pathogenic features discriminate each GBM patient and its sibling GSCs line across GBM clusters.
Distinct functional pools of GSCs distinguish GBM patients across molecular subtypes
We next tested as to whether the differential molecular profiles identified across GBM patients might be associated to an inherent different functional state of GSCs. As shown in Fig. 3A, the specific enhancement in the expression of proliferative signaling distinctive for classical patients was reinforced by the typical rounded morphology displayed by neurospheres obtained through the enrichment of cells isolated from the post-surgery specimen [22, 32, 33, 35]. When the tumor proliferation rate was measured at the protein level in either GBM tissues or in their derivative GSCs by counting the cells positive for the proliferation markers Ki67 and MELK, TCGA-CL samples were confirmed to propagate more extensively than tumors and cells from either TCGA-MS or TCGA-PN subtype (Supplementary Fig. 4A, B). Irregularly shaped clones, characterized by the presence of many protruding and elongated cells, suggestive of increased cell adhesion, were instead generated in GSCs cultures from TCGA-MS cases, enriched in the invasive/infiltrative and EMT signaling pathways. Tumors from these patients were also shown to be highly vascularized, with a total vascular area > threefold higher than in classical tumors, as measured by CD31 + cells (***P < 0.001, **P < 0.01, *P < 0.05, one-way ANOVA, n = 5) (Supplementary Fig. 4A, B). Meanwhile, GSCs from the TCGA-PN counterpart, pinpointed by high levels of signaling involved in neurogenesis and the highest expression of the markers of stemness EphA2 and SSEA1 [35, 46], displayed the typical regular morphology, being much similar to human normal neural stem cells (NSCs) (Fig. 3A and Supplementary Fig. 4A, B).
Next, we assessed whether the increased migratory and invasive phenotype of mesenchymal population might well represent a cell-autonomous trait of these cells by in vitro migration assays [32, 33, 47]. As clearly reported in Fig. 3B, mesenchymal GSCs were able to invade much more efficiently displaying a significant higher infiltrative phenotype as compared to that one of classical GSCs. The pattern displayed by proneural GSCs was similar to that one of mesenchymal GSCs. When GSCs’ global growth trend was compared across subtypes, the inherently expansion rate (Fig. 3C) and self-renewing capacity (Fig. 3D) of mesenchymal and proneural GSCs was shown to remain strikingly lower as compared to their classical counterpart, suggesting that GBM clusters might contain GSCs’ subpopulations intrinsically endowed with differential self-renewal ability (P < 0.0001 TCGA-CL vs. TCGA-MS and P = 0.028, TCGA-CL vs. TCGA-PN, hierarchical linear model for repeated measurements and ***P < 0.0001, ANOVA multiple comparison test, respectively) (Fig. 3C). Strikingly, confidence intervals emerged by cut-offs analysis suggested that growth slope between 0.2532 and 0.2887 should likely classify GSCs cells as classical, between 0.2214 and 0.2609 as proneural and between 0.1888 and 0.2304 as mesenchymal (Fig. 3C). Consistently, differentiated [22, 32, 33, 35] fast-growing GSCs from the TCGA-CL cluster were shown to comprise the highest fraction of abnormal cells which promiscuously co-expressed the astroglial marker GFAP with the neuronal marker Tuj1 [22, 48] as compared to the terminally differentiated mesenchymal and proneural cultures, which displayed a more “mature”, GFAP-positive astroglial phenotype (***P < 0.0001,**P < 0.001,*P < 0.05, Kruskal–Wallis ANOVA, dunn’s test) (Supplementary Fig. 4C).
Findings so far lent to the conclusion that distinct subsets of GSCs endowed with differential state of stemness inherently exist within GBM molecular clusters and might distinguish their classification.
A distinct tumorigenic potential identifies GSCs subpopulation across subtypes
To finally assess a functional association between the different biological hallmarks of GBM subtypes and an inherent enrichment of a specific GSCs subpopulation, we generated GBM xenografts by transplanting GSCs into the brain of immunocompromised SCID mice. As expected, consistent with the in vitro observations, upon orthotopic injection of distinct sets of GSCs from different subtypes, interesting differential tumorigenic and lethal potentials were observed (Fig. 4). More specifically, classical GSCs established tumors which were fast growing (P < 0.001 vs. TCGA-PN and TCGA-MS, hierarchical linear models for repeated measurements) but less invasive as compared to proneural (P = 0.009 vs. TCGA-MS) and mesenchymal ones (Fig. 4A, B). Furthermore, confidence intervals emerged from the comparisons of growth rates indicated that tumors displaying a slope between 0.079 and 0.106 should be likely classified as classical, between 0.040 and 0.063 as proneural and between 0.023 and 0.041 as mesenchymal (Fig. 4A).
Remarkably, as confirmed by serial reconstruction performed at the histological and stereological level in Fig. 4C, the xenografting of classical GSCs established more extended and delimitated tumor masses as compared to those of mesenchymal GSCs, which gave rise to intracranial tumors with a strikingly faster and broader infiltration pattern. Proneural GBM generated upon intracranial transplantation, stood in the middle of this spectrum. Consistently, analysis of overall survival reported that mice receiving classical GSCs died significantly earlier (n = 20 mice; median = 64 days) than those xenografted with mesenchymal (n = 20 mice; median = 132 days) and proneural (n = 17 mice; median = 101 days) cells (P < 0.001, Log-rank test) (Fig. 4D).
All of these data confirmed that distinct GSCs populations might discriminate GBM patients’ phenotype across molecular clusters.
Classification and prediction of GBM patients by means of stemness-related markers
Having perceived in GBM molecular subtypes an association between a distinct biological profile and an inherent enrichment of different GSCs pools, we finally tested the combination of two biomarkers emerged selectively represented for their ability to discriminate and predict GBM clusters. Among the candidates described (Fig. 2B and Supplementary Fig. 2D-F), classification tree analysis revealed that three heterogeneous expression of doublecortin (DCX) and EGFR genes, the former associated with neurodevelopmental and progenitor cell signatures and the latter with stemness [20, 49,50,51,52], almost perfectly segregated GBM patients. In details, patients showing an EGFR expression larger or equal to 9.43 were very likely (93.8%) allocated to the classical cluster. On the other hand, cases described by an EGFR expression lower than 9.43 and, concomitantly, a DCX expression lower than 6.14 belonged were included without error (100%) within the mesenchymal subtype class. Finally, tissues showing EGFR level lower than 9.43 and, concomitantly, a DCX expression larger or equal to 6.14 were surely (100%) characterized by a proneural fingerprint (Fig. 5A). Remarkably, plots in Fig. 5B reinforced the correlation between the inherent complementary regulation of these biomarkers and patient’s TCGA cluster of origin, revealing tissues from the classical subtype identified by an EGFRHigh-DCXLow profile, whereas proneural GBM by a complementary EGFRLow-DCXHigh percentage. Mesenchymal cases accounted for the EGFRLow-DCXLow subpopulation. These results were further confirmed when analyzing EGFR and DCX protein expression in the public dataset . Among GBM subtypes, the highest EGFR level was shown to be distinctive of the classical cluster, whereas DCX protein of the proneural group (Fig. 5C). As shown by Kaplan–Meier plots in Fig. 5D, when analyzing and comparing the outcome of patients defined by the three different patterns of DCX and EGFR expression according to our dataset, cases described by low EGFR and, concomitantly, high DCX level died earlier (n = 74; median = 12 months) than those with an EGFRLow-DCXLow profile (n = 26; median = 13 months) and those depicting high EGFR and low DCX expression, who was the one with the longest survival (n = 17; median = 20 months)(P = 0.0334, Log-rank test and P = 0.0221, Gehan-Breslow-Wilcoxon tests). Remarkably, this was well in agreement with survival curves of GBM patients across subtypes (n = 83 IDH1 wild-type GBM patients), in which proneural cases was shown to exhibit more lethal tumors (n = 13; median = 12 months), than the mesenchymal (n = 28; median = 13 months) and classical ones (n = 42; median = 19 months)(P = 0.0047, Log-rank test and P = 0.0186, Gehan-Breslow-Wilcoxon test) (Fig. 5E).
Data here demonstrate that stemness-related markers differently enriched among GSCs subsets can discriminate and predict the molecular subtype of GBM patients, representing new opportunities for patient-tailored diagnostic and prognostic purposes.
In the last decade, high-throughput profiling analysis have helped decipher glioblastoma heterogeneous molecular profile, outcome and therapeutic response (7–9; 17–20). This lent to the concept of developing molecularly targeted therapies, to provide personalized treatments that are also more effective and less toxic than conventional chemotherapy [54, 55]. Thus, new therapeutics approaches have recently been focused on tackling a number of signaling pathways implicated in GBM development, like those related to intrinsic survival, cell cycle and apoptosis or metabolism, and intrinsically altered molecules [54, 56, 57]. While this approach promises to significantly improve the therapeutic scenario, the development of new markers, the identification of specific molecular targets and the overall process of developing therapeutics for GBM has been hampered by the lack of information as to the actual identity and nature of the normal cell type(s) that were hit by transformation and more importantly by the absence in a detailed understanding of the heterogeneity that exists within this tumor and a fundamental understanding of the cell types that are involved in tumor promotion and relapse. The discovery that GBMs embody cells endowed with tumor-initiating ability and all of the functional features that define stem cells of the CNS, the so-called glioblastoma stem cells (GSCs), has opened new frontiers for the development of potential new therapeutic approaches for this tumor. Major efforts are now directed at targeting GSCs for new therapeutic purposes, endowed with high and selective toxicity towards the tumor but innocuous towards normal cells [58,59,60].
In this context, here, we isolated, established and characterized [22, 32, 33, 35] a large enough culture collection of functionally stable GSCs, collected and studied their surgery tissues of origin (grade IV glioma). In considering GBM heterogeneity and its recent stratification in different subtypes and cellular programs [7, 20, 37, 43], each single GSCs primary line has been subjected to a comprehensive differential analysis for their molecular and antigenic profile, in comparison with the cognate tissue of origin, to each-other and to normal neural stem cells (NSCs) and fully characterized for their stem-like and tumor-initiating characteristics. Here we described that distinct and unique molecular, functional and antigenic profiles of GSCs cells, particularly with regard to stem cells properties and expression of novel critical regulators, were related not only to the transcriptional cluster of their GBM patient of origin but, strikingly, to inherently different states of stemness. In agreement with our previous studies [32, 33], here we found that classical GSCs pool was more undifferentiated and endowed with enhanced expansion capacity as compared to the more mature, relatively slow-propagating proneural and mesenchymal cells, likely endowed with a higher potential for infiltration.
We first identified critical aberrant pathways and networks specific to the distinct GSCs molecular phenotypes. We found that TCGA-MS GBM tissues and their sibling GSCs were defined by highly invasive genetic and mutational programs, which included mediators of tissue invasiveness and angiogenesis, cytoskeleton rearrangements, markers involved in the extracellular communication and cell migration and invasion with a typical invasive signature. TCGA-CL specimens were regulated by effectors included in genomic abnormalities and biological pathways controlling cell growth, proliferation and cycling, while TCGA-PN tissues and their GSCs were defined by CNS development and function, tissue and cellular development mutational and biological patterns and the highest levels of the known “stemness” markers (Fig. 2 and Supplementary Figs. 1–4). Notably, in agreement with our previous studies [22, 32, 35] when comparing side-by-side GBM specimens and sibling GSCs lines, we found that the latter faithfully fully resembled the main aberrant physio-pathological and phenotypic features of their tissue and, thus, of their patient of origin (Supplementary Fig. 3A), suggesting once more the usefulness of GSCs-tailored anti-GBM strategies instead of the “one fits all” GBM therapeutic approaches.
We confirmed the whole molecular features when functionally fully characterizing in vitro GSCs’ fingerprint and physio-pathological hallmarks. Although all of them was shown to possess the full complement of neural stem cells (NSCs) characteristics [22, 32, 33, 61], we show how each GSCs line was endowed with distinct ability to perpetuate and expand across glioblastoma sub-classes. Integrating and confirming their molecular features, we retrieved that GSCs derived from the classical subtype displayed a higher clonal efficiency and were fast growing ex-vivo and much less invasive (Fig. 3B-D). Instead, mesenchymal and proneural GSCs inherently presented the lowest expansion rate and self-renewal capacity, while invading much more efficiently. Focusing on the in vivo functional analyses, despite being all of GSCs lines tested able to produce tumors that closely resemble the main architectural features of glioblastoma [22, 32, 33, 35], we found that these cells differed quite significantly across subtype in their inherent ability to spread throughout the brain parenchyma. The xenografting of mesenchymal GSCs established intracranial human GBMs with a strikingly faster and broader infiltration pattern than those of classical GSCs, which gave rise to larger, more defined neoplastic masses in the mouse brain (Fig. 4C). The proneural GSCs seemed to display a pattern similar to that of mesenchymal GSCs. This not only provided us a better understanding of the intrinsic heterogeneity of GSCs’ basic physiology but also the evidence that GBM subgroups contain GSCs inherently endowed with differential self-renewal ability, thus underpinning their functional inter-cluster variances. Furthermore, the lower self-renewal capacity and expansion rate and the higher invasion efficiency retrieved in GSCs from mesenchymal and proneural subtype, as compared to that one of classical cluster, might also explain the variation in response to therapy of their patients of origin.
The search for innovative, more effective therapies for GBM is now focusing on targeting the relatively rare GSCs cells, that are responsible for tumor growth, resistance to treatments and recurrence. From this perspective, having observed the existence of distinct GSCs’ stemness-related phenotypes across GBM cluster, on their primary specimen and clinical data of the donor patient, we also found that the inherently, stemness-related, different abundance of DCX and EGFR perfectly discriminated GBM patients predicting their stratification into TCGA clusters and their survival (Fig. 5), thus identifying potential subtype-associated biomarkers for patient-tailored diagnostic and prognostic purposes.
In considering the heterogeneous sensitivity of GBM patients at large to standard treatments, which may mostly come from belonging to functionally and molecularly distinct clusters, this also allow to recognize those tumors that may be the most appropriate target for different therapeutic protocols, opening new and exciting possibilities for the definition of novel patient-specific drug-responsiveness protocols that may help to devise appropriate patient-selective therapies, such as patient-tailored chemotherapy regimen, differentiation therapies and combined protocols.
In this study, we identified distinct GSCs cells subsets within GBM patients’ clusters, which are inherently endowed with different self-renewal, invasive and tumorigenic ability underpinning inter-cluster variance together with diverse clinical outcome and response to therapy. Two biomarkers emerged selectively enriched among these different states of stemness, DCX and EGFR, are capable of predicting GBM patients’ clusters, their overall survival and drug sensitivity. These findings may not only promote our understanding of the intrinsic intra-heterogeneity of GBM cases and their GSCs but also provide new potential critical targets opening new and exciting possibilities for GBM prognosis and more effective, individual patient-specific treatment.
Availability of data and materials
Transcriptome, Cytoscan array and Targeted sequencing raw data are available in the Arrayexpress repository under the accession codes E-MTAB-12962, E-MTAB-12963 and E-MTAB-12964, respectively.
Glioblastoma stem cells
Isocitrate dehydrogenase 1
Classical Cancer Genome Atlas
Proneural Cancer Genome Atlas
Mesenchymal Cancer Genome Atlas
Neural stem cells
Copy number variation
Single nucleotide variants
False discovery rate
Telomerase Reverse Transcriptase
Epidermal Growth Factor Receptor
Days post transplantation
Komori T. The molecular framework of pediatric-type diffuse gliomas: shifting toward the revision of the WHO classification of tumors of the central nervous system. Brain Tumor Pathol. 2021;38(1):1–3. PubMed PMID: 33398472. eng.
Stupp R, Roila F, Group EGW. Malignant glioma: ESMO clinical recommendations for diagnosis, treatment and follow-up. Ann Oncol. 2009;20(Suppl 4):126–8 PubMed PMID: 19454432. eng.
Ostrom QT, Gittleman H, Truitt G, Boscia A, Kruchko C, Barnholtz-Sloan JS. CBTRUS statistical report: primary brain and other central nervous system tumors diagnosed in the United States in 2011-2015. Neuro Oncol. 2018;20(suppl_4):iv1–86 PubMed PMID: 30445539. PMCID: PMC6129949. eng.
Kim H, Zheng S, Amini SS, Virk SM, Mikkelsen T, Brat DJ, et al. Whole-genome and multisector exome sequencing of primary and post-treatment glioblastoma reveals patterns of tumor evolution. Genome Res. 2015;25(3):316–27 PubMed PMID: 25650244. PMCID: PMC4352879. Epub 20150203. eng.
Wang J, Cazzato E, Ladewig E, Frattini V, Rosenbloom DI, Zairis S, et al. Clonal evolution of glioblastoma under therapy. Nat Genet. 2016;48(7):768–76 PubMed PMID: 27270107. PMCID: PMC5627776. Epub 20160606. eng.
Binda E, Reynolds BA, Vescovi AL. Glioma stem cells: turpis omen in nomen? (The evil in the name?). J Intern Med. 2014;276(1):25–40 PubMed PMID: 24708237. eng.
Verhaak RG, Hoadley KA, Purdom E, Wang V, Qi Y, Wilkerson MD, et al. Integrated genomic analysis identifies clinically relevant subtypes of glioblastoma characterized by abnormalities in PDGFRA, IDH1, EGFR, and NF1. Cancer Cell. 2010;17(1):98–110 PubMed PMID: 20129251. PMCID: PMC2818769. eng.
Wang Q, Hu B, Hu X, Kim H, Squatrito M, Scarpace L, et al. Tumor Evolution of glioma-intrinsic gene expression subtypes associates with immunological changes in the microenvironment. Cancer Cell. 2017;32(1):42–56.e6. PubMed PMID: 28697342. PMCID: PMC5599156. eng.
Malmström A, Grønberg BH, Marosi C, Stupp R, Frappaz D, Schultz H, et al. Temozolomide versus standard 6-week radiotherapy versus hypofractionated radiotherapy in patients older than 60 years with glioblastoma: the Nordic randomised, phase 3 trial. Lancet Oncol. 2012;13(9):916–26 PubMed PMID: 22877848. Epub 20120808. eng.
Yang P, Zhang W, Wang Y, Peng X, Chen B, Qiu X, et al. IDH mutation and MGMT promoter methylation in glioblastoma: results of a prospective registry. Oncotarget. 2015;6(38):40896–906 PubMed PMID: 26503470. PMCID: PMC4747376. eng.
Brennan CW, Verhaak RG, McKenna A, Campos B, Noushmehr H, Salama SR, et al. The somatic genomic landscape of glioblastoma. Cell. 2013;155(2):462–77 PubMed PMID: 24120142. PMCID: PMC3910500. eng.
Noushmehr H, Weisenberger DJ, Diefes K, Phillips HS, Pujara K, Berman BP, et al. Identification of a CpG island methylator phenotype that defines a distinct subgroup of glioma. Cancer Cell. 2010;17(5):510–22 PubMed PMID: 20399149. PMCID: PMC2872684. Epub 20100415. eng.
Hovinga KE, McCrea HJ, Brennan C, Huse J, Zheng J, Esquenazi Y, et al. EGFR amplification and classical subtype are associated with a poor response to bevacizumab in recurrent glioblastoma. J Neurooncol. 2019;142(2):337–45 PubMed PMID: 30680510. PMCID: PMC6752204. Epub 20190124. eng.
Quail DF, Joyce JA. The microenvironmental landscape of brain tumors. Cancer Cell. 2017;31(3):326–41 PubMed PMID: 28292436. PMCID: PMC5424263. eng.
Heiland DH, Haaker G, Delev D, Mercas B, Masalha W, Heynckes S, et al. Comprehensive analysis of PD-L1 expression in glioblastoma multiforme. Oncotarget. 2017;8(26):42214–25 PubMed PMID: 28178682. PMCID: PMC5522061. eng.
Martinez-Lage M, Lynch TM, Bi Y, Cocito C, Way GP, Pal S, et al. Immune landscapes associated with different glioblastoma molecular subtypes. Acta Neuropathol Commun. 2019;7(1):203 PubMed PMID: 31815646. PMCID: PMC6902522. Epub 20191129. eng.
Kaffes I, Szulzewsky F, Chen Z, Herting CJ, Gabanic B, Velázquez Vega JE, et al. Human Mesenchymal glioblastomas are characterized by an increased immune cell presence compared to Proneural and Classical tumors. Oncoimmunology. 2019;8(11):e1655360 PubMed PMID: 31646100. PMCID: PMC6791439. Epub 20190822. eng.
Patel AP, Tirosh I, Trombetta JJ, Shalek AK, Gillespie SM, Wakimoto H, et al. Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma. Science. 2014;344(6190):1396–401 PMCID: PMC4123637. Epub 2014/06/12. eng.
Sottoriva A, Spiteri I, Piccirillo SG, Touloumis A, Collins VP, Marioni JC, et al. Intratumor heterogeneity in human glioblastoma reflects cancer evolutionary dynamics. Proc Natl Acad Sci U S A. 2013;110(10):4009–14 PubMed PMID: 23412337. PMCID: PMC3593922. Epub 2013/02/14. eng.
Neftel C, Laffy J, Filbin MG, Hara T, Shore ME, Rahme GJ, et al. An integrative model of cellular states, plasticity, and genetics for glioblastoma. Cell. 2019;178(4):835–49.e21. PubMed PMID: 31327527. PMCID: PMC6703186. Epub 2019/07/18. eng.
Singh SK, Clarke ID, Hide T, Dirks PB. Cancer stem cells in nervous system tumors. Oncogene. 2004;23(43):7267–73 PubMed PMID: 15378086. eng.
Galli R, Binda E, Orfanelli U, Cipelletti B, Gritti A, De Vitis S, et al. Isolation and characterization of tumorigenic, stem-like neural precursors from human glioblastoma. Cancer Res. 2004;64(19):7011–21 PubMed PMID: 15466194. eng.
Hemmati HD, Nakano I, Lazareff JA, Masterman-Smith M, Geschwind DH, Bronner-Fraser M, et al. Cancerous stem cells can arise from pediatric brain tumors. Proc Natl Acad Sci U S A. 2003;100(25):15178–83 PubMed PMID: 14645703. PMCID: PMC299944. eng.
Reya T, Morrison SJ, Clarke MF, Weissman IL. Stem cells, cancer, and cancer stem cells. Nature. 2001;414(6859):105–11 PubMed PMID: 11689955. Epub 2001/11/02. eng.
Bao S, Wu Q, McLendon RE, Hao Y, Shi Q, Hjelmeland AB, et al. Glioma stem cells promote radioresistance by preferential activation of the DNA damage response. Nature. 2006;444(7120):756–60 PubMed PMID: 17051156. eng.
Chen J, Li Y, Yu TS, McKay RM, Burns DK, Kernie SG, et al. A restricted cell population propagates glioblastoma growth after chemotherapy. Nature. 2012;488(7412):522–6 PubMed PMID: 22854781. PMCID: Pmc3427400. Epub 2012/08/03. eng.
Körber V, Yang J, Barah P, Wu Y, Stichel D, Gu Z, et al. Evolutionary Trajectories of IDH. Cancer Cell. 2019;35(4):692–704.e12. PubMed PMID: 30905762. Epub 2019/03/21. eng.
Dirkse A, Golebiewska A, Buder T, Nazarov PV, Muller A, Poovathingal S, et al. Stem cell-associated heterogeneity in Glioblastoma results from intrinsic tumor plasticity shaped by the microenvironment. Nat Commun. 2019;10(1):1787 PubMed PMID: 30992437. PMCID: PMC6467886. Epub 20190416. eng.
Richards LM, Whitley OKN, MacLeod G, Cavalli FMG, Coutinho FJ, Jaramillo JE, et al. Gradient of Developmental and Injury Response transcriptional states defines functional vulnerabilities underpinning glioblastoma heterogeneity. Nat Cancer. 2021;2(2):157–73 PubMed PMID: 35122077. Epub 20210104. eng.
Shen Y, Grisdale CJ, Islam SA, Bose P, Lever J, Zhao EY, et al. Comprehensive genomic profiling of glioblastoma tumors, BTICs, and xenografts reveals stability and adaptation to growth environments. Proc Natl Acad Sci U S A. 2019;116(38):19098–108 PubMed PMID: 31471491. PMCID: PMC6754609. Epub 20190830. eng.
Stringer BW, Day BW, D’Souza RCJ, Jamieson PR, Ensbey KS, Bruce ZC, et al. A reference collection of patient-derived cell line and xenograft models of proneural, classical and mesenchymal glioblastoma. Sci Rep. 2019;9(1):4902 PubMed PMID: 30894629. PMCID: PMC6427001. Epub 20190320. eng.
Binda E, Visioli A, Giani F, Trivieri N, Palumbo O, Restelli S, et al. Wnt5a drives an invasive phenotype in human glioblastoma stem-like cells. Cancer Res. 2017;77(4):996–1007 PubMed PMID: 28011620. Epub 2016/12/23. eng.
Trivieri N, Visioli A, Mencarelli G, Cariglia MG, Marongiu L, Pracella R, et al. Growth factor independence underpins a paroxysmal, aggressive Wnt5a(High)/EphA2(Low) phenotype in glioblastoma stem cells, conducive to experimental combinatorial therapy. J Exp Clin Cancer Res. 2022;41(1):139 PubMed PMID: 35414102. PMCID: PMC9004109. Epub 20220412. eng.
Parada LF, Dirks PB, Wechsler-Reya RJ. Brain tumor stem cells remain in play. J Clin Oncol. 2017;35(21):2428–31 PubMed PMID: 28640710. PMCID: PMC5516484. Epub 2017/06/22. eng.
Binda E, Visioli A, Giani F, Lamorte G, Copetti M, Pitter KL, et al. The EphA2 receptor drives self-renewal and tumorigenicity in stem-like tumor-propagating cells from human glioblastomas. Cancer Cell. 2012;22(6):765–80 PubMed PMID: 23238013. eng..
Visioli A, Giani F, Trivieri N, Pracella R, Miccinilli E, Cariglia MG, et al. Stemness underpinning all steps of human colorectal cancer defines the core of effective therapeutic strategies. EBioMedicine. 2019;44:346–60 PubMed PMID: 31056474. Epub 2019/05/02. eng.
Network CGAR. Comprehensive genomic characterization defines human glioblastoma genes and core pathways. Nature. 2008;455(7216):1061–8 PubMed PMID: 18772890. PMCID: PMC2671642. eng.
Gaffney SG, Townsend JP. PathScore: a web tool for identifying altered pathways in cancer data. Bioinformatics. 2016;32(23):3688–90 PubMed PMID: 27503224. Epub 2016/08/08. eng.
Liberzon A, Subramanian A, Pinchback R, Thorvaldsdóttir H, Tamayo P, Mesirov JP. Molecular signatures database (MSigDB) 3.0. Bioinformatics. 2011;27(12):1739–40 PubMed PMID: 21546393. PMCID: PMC3106198. Epub 20110505. eng.
Trivieri N, Pracella R, Cariglia MG, Panebianco C, Parrella P, Visioli A, et al. BRAF(V600E) mutation impinges on gut microbial markers defining novel biomarkers for serrated colorectal cancer effective therapies. J Exp Clin Cancer Res. 2020;39(1):285 PubMed PMID: 33317591. PMCID: PMC7737386. Epub 20201214. eng.
Singer JD, Willett JB. Applied longitudinal data analysis: modeling change and event occurrence. Oxford: Oxford University Press; 2003. xx, 644.
Diggle P, Liang K-Y, Zeger SL. Analysis of longitudinal data. Oxford New York: Clarendon Press; Oxford University Press; 1994. xi, p. 253.
Phillips HS, Kharbanda S, Chen R, Forrest WF, Soriano RH, Wu TD, et al. Molecular subclasses of high-grade glioma predict prognosis, delineate a pattern of disease progression, and resemble stages in neurogenesis. Cancer Cell. 2006;9(3):157–73 PubMed PMID: 16530701. eng.
D’Angelo F, Ceccarelli M, Tala, Garofano L, Zhang J, Frattini V, et al. The molecular landscape of glioma in patients with Neurofibromatosis 1. Nat Med. 2019;25(1):176–87 PubMed PMID: 30531922. PMCID: PMC6857804. Epub 2018/12/10. eng.
Tokheim CJ, Papadopoulos N, Kinzler KW, Vogelstein B, Karchin R. Evaluating the evaluation of cancer driver genes. Proc Natl Acad Sci U S A. 2016;113(50):14330–5 PubMed PMID: 27911828. PMCID: PMC5167163. Epub 20161122. eng.
Dell’Albani P. Stem cell markers in gliomas. Neurochem Res. 2008;33(12):2407–15 PubMed PMID: 18493853. eng.
Pennacchietti S, Michieli P, Galluzzo M, Mazzone M, Giordano S, Comoglio PM. Hypoxia promotes invasive growth by transcriptional activation of the met protooncogene. Cancer Cell. 2003;3(4):347–61 PubMed PMID: 12726861. eng.
Gritti A, Parati EA, Cova L, Frolichsthal P, Galli R, Wanke E, et al. Multipotential stem cells from the adult mouse brain proliferate and self-renew in response to basic fibroblast growth factor. J Neurosci. 1996;16(3):1091–100 PubMed PMID: 8558238. eng.
Völker J, Engert J, Völker C, Bieniussa L, Schendzielorz P, Hagen R, et al. Isolation and characterization of neural stem cells from the rat inferior colliculus. Stem Cells Int. 2019;2019:5831240 PubMed PMID: 31781242. PMCID: PMC6875198. Epub 20191029. eng.
Odrzywolski A, Jarosz B, Kiełbus M, Telejko I, Ziemianek D, Knaga S, et al. Profiling glioblastoma cases with an expression of DCX, OLIG2 and NES. Int J Mol Sci. 2021;22(24):13217 PubMed PMID: 34948016. PMCID: PMC8708973. Epub 20211208. eng.
Brown JP, Couillard-Després S, Cooper-Kuhn CM, Winkler J, Aigner L, Kuhn HG. Transient expression of doublecortin during adult neurogenesis. J Comp Neurol. 2003;467(1):1–10 PubMed PMID: 14574675. eng.
Liffers K, Lamszus K, Schulte A. EGFR amplification and glioblastoma stem-like cells. Stem Cells Int. 2015;2015:427518 PubMed PMID: 26136784. PMCID: PMC4468289. Epub 20150602. eng.
LB W, A K, MA G, JE K, S C, Y L, et al. Proteogenomic and metabolomic characterization of human glioblastoma. Cancer cell. 2021 02/11/2021. PubMed PMID: 33577785.
El Atat O, Naser R, Abdelkhalek M, Habib RA, El Sibai M. Molecular targeted therapy: a new avenue in glioblastoma treatment. Oncol Lett. 2023;25(2):46 PubMed PMID: 36644133. PMCID: PMC9811647. Epub 20221215. eng.
Garrett AM, Lastakchi S, McConville C. The personalisation of glioblastoma treatment using whole exome sequencing: a pilot study. Genes (Basel). 2020;11(2):173 PubMed PMID: 32041307. PMCID: PMC7074406. Epub 20200206. eng.
Alifieris C, Trafalis DT. Glioblastoma multiforme: Pathogenesis and treatment. Pharmacol Ther. 2015;152:63–82 PubMed PMID: 25944528. Epub 20150502. eng.
Han B, Wang R, Chen Y, Meng X, Wu P, Li Z, et al. QKI deficiency maintains glioma stem cell stemness by activating the SHH/GLI1 signaling pathway. Cell Oncol (Dordr). 2019;42(6):801–13 PubMed PMID: 31292920. Epub 20190710. eng.
Mitchell K, Troike K, Silver DJ, Lathia JD. The evolution of the cancer stem cell state in glioblastoma: emerging insights into the next generation of functional interactions. Neuro Oncol. 2021;23(2):199–213 PubMed PMID: 33173943. PMCID: PMC7906055. eng.
Lin Q, Ba T, Ho J, Chen D, Cheng Y, Wang L, et al. First-in-human trial of EphA2-Redirected CAR T-Cells in patients with recurrent glioblastoma: a preliminary report of three cases at the starting dose. Front Oncol. 2021;11:694941 PubMed PMID: 34235085. PMCID: PMC8256846. Epub 20210621. eng.
Prager BC, Bhargava S, Mahadev V, Hubert CG, Rich JN. Glioblastoma stem cells: driving resilience through chaos. Trends Cancer. 2020;6(3):223–35 PubMed PMID: 32101725. PMCID: PMC8779821. Epub 20200203. eng.
Vescovi AL, Parati EA, Gritti A, Poulin P, Ferrario M, Wanke E, et al. Isolation and cloning of multipotential stem cells from the embryonic human CNS and establishment of transplantable human neural stem cell lines by epigenetic stimulation. Exp Neurol. 1999;156(1):71–83 PubMed PMID: 10192778. eng.
We are thankful to Gilberto S. Pazienza for his loving cooperation and to Lucia Sergisergi for kindly providing the luciferase lentivirus.
This work was supported by grants from the “Italian Ministery of Health” (GR-2011–02351534 and Progetto Ricerca Corrente 2018–20; 2020–2023) to EB and from “Associazione Italiana Cancro” (IG-22027) to ALV.
Ethics approval and consent to participate
The cohort of glioma patients analyzed in this study was enrolled at IRCCS National Neurologic Institute “C. Besta” after approval of the institutional Ethic board number protocol n°02 and protocol n°61. All the subject with a confirmed diagnosis of GBM agreed to participate according to the ethical guidelines of the 2013 Declaration of Helsinki with signed informed consent. Material was anonymized at the time of collection. All animal procedures were performed according to the Guidelines for the Care and Use of Laboratory Animals and experimental protocols approved by the Italian Ministry of Health (805/2016-PR and 119/2019-PR).
Consent for publication
All authors approved the manuscript and consent publication.
The authors declare no competing interests. AL Vescovi has ownership interest in Hyperstem SA.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Visioli, A., Trivieri, N., Mencarelli, G. et al. Different states of stemness of glioblastoma stem cells sustain glioblastoma subtypes indicating novel clinical biomarkers and high-efficacy customized therapies. J Exp Clin Cancer Res 42, 244 (2023). https://doi.org/10.1186/s13046-023-02811-0