Next Article in Journal
PD-1-Associated Gene Expression Signature of Neoadjuvant Trastuzumab-Treated Tumors Correlates with Patient Survival in HER2-Positive Breast Cancer
Next Article in Special Issue
Metformin as Potential Therapy for High-Grade Glioma
Previous Article in Journal
Bioactive Compounds: Multi-Targeting Silver Bullets for Preventing and Treating Breast Cancer
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 11
Issue 10
10.3390/cancers11101564
cancers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic and Lineage Classification of Glioma-Initiating Cells Identifies a Clinically Relevant Glioblastoma Model

by
Norihiko Saito
*,
Nozomi Hirai
,
Kazuya Aoki
,
Sho Sato
,
Ryo Suzuki
,
Yu Hiramoto
,
Satoshi Fujita
,
Haruo Nakayama
,
Morito Hayashi
,
Takatoshi Sakurai
and
Satoshi Iwabuchi
Department of Neurosurgery, Toho University Ohashi Medical Center, Tokyo 153-8515, Japan
*
Author to whom correspondence should be addressed.
Cancers 2019, 11(10), 1564; https://doi.org/10.3390/cancers11101564
Submission received: 25 August 2019 / Revised: 4 October 2019 / Accepted: 8 October 2019 / Published: 15 October 2019
(This article belongs to the Special Issue Recent Advances in Glioblastoma)
Browse Figures
Versions Notes

Abstract

:
The Cancer Genome Atlas (TCGA) project described a robust gene expression-based molecular classification of glioblastoma (GBM), but the functional and biological significance of the subclasses has not been determined. The present comprehensive analysis of 25 glioma-initiating cell (GIC) lines classifies GIC lines into four subtypes (classical, mesenchymal, proneural, and neural) that are closely related to the TCGA GBM subclasses and display distinct lineage characteristics and differentiation behavior that recapitulate neural development. More importantly, the GIC subtypes exhibit distinct biological phenotypes in relation to self-renewal capacity, proliferation, invasiveness, and angiogenic potential in vitro and in vivo. In addition, the GIC subtypes exhibit divergent patterns of signaling pathway activation and deactivation of the Wnt, Notch, and TGF-β pathways. These results will improve drug discovery targeting certain genetic mutation in glioblastoma and improve the development of precision medicine.

1. Introduction

Despite advances in our understanding of the molecular aspects of glioblastoma (GBM), the prognosis for glioblastoma remains dismal [1,2,3,4]. Recent studies indicate that some neoplastic cells within human glioblastoma are capable of self-renewal and multilineage differentiation—properties associated with normal neural stem cells [5,6,7]. These stem-like tumor cells, known as glioma-initiating cells (GICs), are responsible for tumor initiation and recurrence after therapy, which makes them attractive targets for novel GBM therapies [8,9,10]. Other studies indicate that GICs retain relevant molecular features of human GBM and therefore may enable development of better preclinical models for evaluating tumor biology and therapeutics [11].
GICs cultured under neural stem cell conditions exhibit heterogeneous biological characteristics and molecular profiles. Previous studies classified GICs based on CD133 expression and found that the CD133-negative subtype had a lower proliferation index [12]. Gunther et al. [13] proposed that two major GIC subtypes exist. One exhibits a full stem-like phenotype, with spherical growth in vitro, expression of CD133, capacity for broad neuroglial differentiation, and high tumorigenicity and invasive growth in vivo. The other displays a restricted stem-like phenotype that only partially satisfies the relevant criteria. More recent work from Lottaz et al. [14] divided 17 GICs into two groups, one expressing a “proneural” gene signature resembling normal fetal brain stem cells and the other expressing a mesenchymal gene signature, which is more similar to that of adult neural stem cells. These findings suggest that differences in GBM tumors may be related to variation in GIC ancestors.
Several signaling pathways are involved in maintaining stemness and tumorigenicity of GICs. The TGF-β pathway is active in high-grade gliomas and is associated with a poor prognosis [15,16]. It increases the self-renewal capacity of GICs through the induction of Leukemia inhibitory factor (LIF) and Sex determining region Y-box 2 (Sox2) and subsequent activation of the JAK-STAT pathway [7,17]. Notch signaling is required in order to prevent neuronal differentiation and promote neural stem cell maintenance for further commitment into a glial lineage [18]. In mice, blockade of Notch signaling reduced GIC growth and clonogenicity in vitro, as well as tumor formation [10,19]. Other pathways involved in GIC maintenance and function include the basic fibroblast growth factor (bFGF), Wnt, and the hedgehog signaling pathways [20,21,22]. Signal transduction pathways that modulate GIC growth and differentiation have been extensively studied; however, it remains unclear if these signals equally or differentially contribute to stem cell maintenance in GICs.
A distinguishing feature of GBM is marked genetic heterogeneity within individual tumors. Phillips and colleagues described a 35-gene signature that divided GBM into three subtypes (which resembled different stages of neurogenesis) that correlated with prognosis and established a pattern of disease progression [23]. Recently, The Cancer Genome Atlas (TCGA) included a robust gene expression-based molecular classification of GBM into proneural, neural, classical, and mesenchymal subtypes and highlighted strong relationships with the different neural lineages [24]. Further studies are needed in order to interpret the biological and functional significance of these GBM subclasses.
In this study, we analyzed the gene expression profiles of 25 GICs derived from fresh GBM tissues and identified four GIC subtypes, which closely resemble the recently identified TCGA GBM subclasses. We also observed diverse lineage characteristics in the GIC subtypes and differentiation behavior recapitulating neural development. In addition, the GIC subtypes exhibited distinct biological phenotypes in relation to self-renewal capacity, proliferation, invasiveness, and angiogenic potential in vitro and in vivo. More importantly, the GIC subtypes exhibited divergent patterns of signaling pathway activation and deactivation of the Wnt, Notch, and TGF-β pathways. This comprehensive genetics-based classification of GICs lays the groundwork for a better molecular understanding of GBM pathway signaling, which could assist in the development of personalized therapies for certain GBM patient subgroups.

2. Results

2.1. Identification of Four GIC Subtypes by Cluster Analysis of Gene Expression Profiling Data

To identify potential GIC subtypes, we performed cluster analysis of gene expression data (Affymetrix U133A2) from 25 GIC cell lines generated from fresh GBM tissues (Table S1). Two approaches were used for the analysis. First, we performed supervised clustering using 1461 probe sets matching 840 genes reported to be differentially expressed in GBM subclasses, which effectively separated GBMs into four previously described subclasses [24]. This separated 25 GICs into two major clusters-one comprising subclusters of GICs with high expressions of “neural” and “proneural” genes and another comprising subclusters of GICs with high expressions of “classical” and “mesenchymal” genes (Figure 1a). Second, we performed unsupervised clustering with the top 1600 varied probe sets. Surprisingly, this yielded four subclusters that exactly matched those identified in supervised clustering (Figure 1b). There were 275 probe sets (~20% of total) common to these two methods (Figure 1c). This second method was better able to separate “mesenchymal” genes, presumably because more “mesenchymal” gene probe sets were present in this clustering. Notably, the “proneural” subclass showed expression of “proneural” and “classical” subtypes. The fact that both methods yielded the same classification pattern confirmed our clustering results and further indicates that subtypes specific to GBM gene expressions might be determined by varied GIC subtypes that drive proliferation and progression of these tumors. Some genes highly enriched in GIC subtypes were bone morphogenetic protein 4 (BMP4), doublecortin (DCX), p16INK4a, and inhibitor of differentiation 2 (ID2) for the neural subtype; oligodendrocyte lineage transcription factor 2 (OLIG2), NK2 homeobox 2 (NKX2-2), Notch1, and Notch3 for the proneural and classical subtypes; and transforming growth factor beta receptor 2 (TGFBR2), CD44, and caveolin 1 (CAV1) for the mesenchymal subtype. Principal component analysis (PCA) of gene expression of the four GIC subtypes (Figure 1d) indicated that the proneural and classical GICs were most closely related and that neural and mesenchymal GICs deviated in a direction opposite to them, which supports the hypothesis that the “proneural” group has two GIC subtypes.

2.2. Distribution of Frequently Mutated Genes Across GIC Subtype Is Relevant to TCGA GBM Subclass

Mutation analysis was performed for a panel of genes known to be important in GBM. The p53 mutation was enriched in neural and mesenchymal subtypes, whereas the phosphatase and tensin homolog (PTEN) mutation was prevalent in classical and mesenchymal GICs. Mutation of p53 and IDH1 was present in neural GICs, whereas PI3K-kinase mutations (PI3CA) were present only in proneural GICs. In addition, co-mutation of p53, IDH1, and Akt 3 was observed in the neural subtype; however, co-mutation of PI3K genes was rare in all subtypes (Table 1). No EGFR mutation was present in GICs.

2.3. GIC Subtypes Exhibit Distinct Lineage Characteristics That Recapitulate Neural Development

To clarify the biological importance of the subtypes, we used probe set data presented by Cahoy et al. [25] to define gene sets associated with neurons, oligodendrocytes, astrocytes, and cultured astroglial cells. These four gene sets were used to calculate a single-sample gene set–enrichment analysis (ssGSEA) enrichment score for all samples, to indicate how closely expression in a sample reflects the expected expression pattern of the gene set. This analysis yielded a number of patterns that associated each subtype with expression patterns from purified murine neural cell types (Figure 2a). The proneural subtype was highly enriched with proneural, oligodendrocytic, and astrocytic signatures, whereas the classical group was exclusively associated with the murine astrocytic signature. Neural subtype was associated with proneural signature but also showed strong enrichment for genes differentially expressed by neurons. This suggests that the neural subtype is more committed than the proneural subtype in differentiated lineages. The mesenchymal subtype was exclusively associated with the in vitro cultured astroglial signature. These results are highly consistent with the lineage profiles of TCGA GBM subclasses [24].
Immunostaining was used to examine expressions of a panel of neural lineage markers in GIC subtypes (Figure 2b). GICs differentially expressed lineage markers related to specific classes of neural stem/progenitors and exhibited a bidirectional tendency toward neural or mesenchymal development. Notably, expressions of neural stem cell markers Sox2, Olig2, and Nestin were higher in proneural and classical subtypes, whereas expressions of glial fibrillary acidic protein (GFAP) and S100 calcium binding protein B (S100β), an astrocytic marker, were higher in neural GICs. Glial progenitor markers such as A2B5 and glutamate aspartate transporter (GLAST) were also expressed in proneural and classical subtypes. Mesenchymal GICs highly expressed the mesenchymal marker YKL-40, but positivity on neural lineage markers was limited. These findings suggest that classical and proneural GICs are more similar to the normal human fetal brain–derived stem cell line HFB2050. Further, classical and proneural subtypes are closely related and are a more primitive stem cell, whereas the neural subtype includes cells that are more differentiated. We also compared marker expression pattern in in vitro culture and xenografts in relation to GIC subtype (Figure 2c). Our results showed that GIC subtype lineage characteristics were retained in the xenograft: classical and proneural xenografts highly expressed Nestin, Olig2, and YKL40, and the mesenchymal xenograft had high YKL-40 expression but low expressions of Nestin, GFAP, and Olig2.

2.4. GIC Subtypes Exhibited Distinct Differentiation Behavior In Vitro

The differentiation potential of GIC subtypes was examined in culture under 1% (fetal bovine serum) FBS plus 1 μM retinoic acid. Classical and proneural GICs were able to differentiate into astrocytic, neuronal, and oligodendrocytic lineages. While neural GICs were mainly able to differentiate into astrocytic and neuronal lineages, mesenchymal GICs seemed resistant to FBS and retinoic acid (RA) treatment and failed to differentiate into neural lineage cells (Figure 3, Table S2).

2.5. Differential Growth Characteristics in GIC Subtypes

To identify growth characteristics of GIC subtypes and evaluate their capacity to form tumor spheres, we examined their growth properties in vitro. The GIC subtypes exhibited varied potential to form primary and secondary tumor spheres (Figure 4a,b). While classical GICs readily formed large tumor spheres, mesenchymal GICs did not. Instead, they formed loosely packed tumor spheres that were smaller than those of other subtypes.

2.6. GIC Subtypes Exhibit Varied Biological Behavior In Vitro and In Vivo

The in vivo biological behavior of GIC subtypes was studied by orthotopic injection of cells into mouse brain. The histological phenotypes of tumors formed differed in relation to GIC subtype (Table 2). Mesenchymal GICs exhibited high tumorgenicity in vivo despite the absence of CD133 expression. Mesenchymal xenografts were highly invasive, whereas neural xenografts were the least invasive, forming a well-defined tumor margin separating tumor from surrounding normal tissues (Figure 5a). In addition, the neural xenograft subtype was the least angiogenic, as indicated by von Willebrand factor(vWF) staining (Figure 5a). These results are consistent with those of in vitro culture, which showed that neural GICs produced significantly lower amounts of vascular endothelial growth factor (VEGF) and matrix metallopeptidase 9 (MMP-9) (Figure 5b)—key regulators of tumor angiogenesis and invasion, respectively. The proliferative potential of GICs was examined by ki67 staining in cultured GIC and xenografts (Figure 5a,b). Ki67 staining showed that classical GICs were most active in proliferation in vitro and in vivo; they exhibited intratumoral necrotic regions that were coextensive with extravasating blood cells. Thus, the in vitro and in vivo characteristics of GIC biological behaviors support the bioinformatic classification of GICs.

2.7. Signaling Pathway Activation and Deactivation Regulate GIC Subtype

GICs exhibited subtype-specific lineage markers and differentiation potentials, which is reminiscent of the hierarchical normal neural development process in which signaling pathways (e.g., the Notch and Wnt pathways) are differentially tuned. ANOVA analysis was used to identify genes that were significantly differentially expressed among GIC subtypes (p ≤ 0.05). The genes identified overlapped with a manually curated gene list comprising known genes involved in the Wnt, Notch, and TGF-β signaling pathways, which revealed differential activation and deactivation of signaling pathways within GIC subtypes. In mesenchymal GICs, the TGF-β pathway component was highly expressed, concurrent with deactivation of Notch and Wnt, as physiological inhibitors of Notch and Wnt are abundant in this subtype. Notch pathway components were highly enriched in classical and proneural GICs but deactivated in neural GICs (Figure 6).

3. Discussion

Accumulating evidence from TCGA has yielded a robust gene expression–based molecular classification of GBM into proneural, neural, classical, and mesenchymal subtypes [24]. Identification of valid GBM subtype counterparts in GICs represents an important aid to studying GBM subtypes, in particular for modeling and predicting therapeutic response. However, numerous important questions remain unanswered, such as the extent to which the GIC model preserves the genetic and biological features of GBMs during in vitro culture. Here, we report the establishment of GIC lines and identify four GIC subtypes by means of clustering analysis of gene expression profiles. Our results show that unsupervised clustering of GIC gene expression analysis matched perfectly with that from a supervised clustering using 840 genes previously used to separate GBM tumors into four subtypes. These findings confirm that, despite in vitro handling of GIC lines, the transcriptome of GIC lines resembles that of primary GBM tumors. Our study identified four subtypes (neural, proneural, classical, and mesenchymal) in GIC lines, though unsupervised clustering contained only ~20% of supervised clustering probe sets, which suggests that these subtypes are intrinsic within GIC lines and may be drivers of GBM subtypes.
Recent genome-wide profiling studies have reported genetic abnormalities associated with GBM subtypes [24,26]. These studies reported p53 gene mutations, mostly in proneural and mesenchymal GBM, and PTEN gene mutations in classical and mesenchymal GBM [11,24]. Our GIC study showed a similar distribution profile, although frequencies were higher, suggesting an intrinsic association between GBM and GIC. In addition, a higher frequency of mutations in PIK3CA and PIK3R1 was observed in GIC subtypes, but none coexisted in the same GIC line; thus; mutation in either PI3K gene may be sufficient to drive GIC growth. Although some studies found no AKT1 mutation in GBM [27,28], we identified the AKT3E17K mutation in GIC. AKT3 has a pivotal role in human GBM biology [29]; therefore, assessing the functional role of AKT3 activation by somatic mutations in GBM is relevant in identifying its role in this aggressive disease. In addition, IDH1 is mutated in >80% of secondary GBM, although <10% of primary GBM harbor these alterations [30,31]. Recently, TCGA revealed IDH mutation in proneural GBM and frequent co-mutation with p53 [24]. Our data confirm the low frequency of IDH1 mutation in a GIC line derived from primary GBM and revealed co-mutation of IDH1 and p53 in a proneural GIC line, which suggests that co-mutation of IDH1 and p53 is crucial in maintaining this GIC subtype. An important finding of TCGA analysis was that EGFR gene mutation was present in almost half of GBM tumors examined. However, no EGFR mutation was seen in our GIC lines, suggesting that cells with mutated EGFR might be lost or selected against during culture.
The developing and adult nervous system has distinct classes of neural stem/progenitors in the lineage hierarchy. Recent studies reported that glioma cells expressing lineage markers such as A2B5 [32], NG2 [33], CD44 [29], and even GFAP also meet the criteria for tumorigenic stem cells, suggesting that GIC originate from a broader spectrum of neural lineages. Our findings indicate that GICs are a heterogeneous population and that classical and proneural subtypes are more primitive, that neural subtype is more differentiated, and that the mesenchymal subtype seems to deviate from neural lineage through mesenchymal transition expressing the mesenchymal marker YKL40. The lineage hierarchy status of GIC subtypes may confer variation in the tendency to differentiate; the proneural and classical GIC subtypes, for example, highly express markers of the stem/early progenitors and can differentiate into trilineage cells. By contrast, the neural subtype exhibits higher GFAP expression, like adult subventricular zone (SVZ) astrocytes, and has less potential upon induction to differentiate into bilineage cells. The mesenchymal subtype of GIC is almost incapable of differentiating into neural lineage cells. Thus, the fact that gliomas have a large, continuous histological spectrum with regard to proportions of various differentiated and anaplastic cells may be attributable to the phenotype of the underlying tumor-initiating cells. If combined with mutation analysis, it is plausible that neural stem/progenitors or even differentiated cells can be converted to a stem cell state through genetic disruptions of specific sets of genes.
Tumor stem cells reside in a microenvironment known as the “stem cell niche”, which maintains them in a stem-like state [34]. Multiple recent studies have clearly demonstrated that—by delivering special signals to balance cell proliferation, self-renewal, and differentiation-the stem cell niche within glioma is a key regulator of stem cell fate [35,36,37]. The present data show that Wnt, Notch, and TGF-β signals are prominent in controlling cell deactivation of proliferation, self-renewal, and differentiation in GICs. Whether these signals equally or differentially contribute to stem cell maintenance in GICs remains unresolved [10,11]. We found that GICs exhibited subtype-specific activation and deactivation of these signaling pathways. Interestingly, in the mesenchymal subtype, TGF-β signaling was concurrent with attenuation of Wnt and Notch signaling. In contrast, in the proneural subtype, Notch signaling, rather than TGF-β, was highly activated. Disruption of these signaling pathways by specific inhibitors impaired growth of GIC subtypes. These results suggest that GIC subtypes might be maintained by varied activation of signaling pathways. Therefore, targeting these GIC subtypes by utilizing potential differences in stem cell–niche interaction might allow successful stem cell targeting and improve outcomes of cancer treatment.

4. Materials and Methods

4.1. Cell Lines and Cell Culture

The 25 glioma initiating cell lines were maintained in neurosphere medium by using a previously described method [38] to isolate neurosphere-forming cells from surgical specimens of human GBM. The study was approved by the Institutional Review Board of Toho University (H22-62). These GIC lines were cultured as GBM neurospheres in DMEM/F12 medium supplemented with B27 (Invitrogen, Grand Island, NY, USA), L-glutamine (GIBCO), penicillin/streptomycin, and growth factors (20 ng/mL EGF and 20 ng/mL FGF-2; Invitrogen). The human fetal brain neural stem cell line HFB2050 was kindly provided by Evan Snyder and cultured as described previously. For differentiation analysis, GICs were grown on polylysine-coated coverslips and treated with retinoic acid (1µM) (Sigma-Aldrich, St. Louis, MO, USA) plus 1% fetal bovine serum (Sigma-Aldrich, St. Louis, MO, USA) in DMEN/F12 for 5 days.

4.2. Gene Arrays and Bioinformatic

GIC expression profiles were generated on an Affymetrix GeneChip HG-U133A 2.0 microarray platform [37]. Clustering analysis was performed with Cluster and Tree View software (Cluster 3.0, Software copyright Stanford University, USA). For supervised clustering, we used 1461 probe sets that matched 840 previously identified subtype-specific genes to cluster GIC samples. The parameters used were correlation uncentred and average linkage. For unsupervised clustering, we used the top 1600 varied probes among all GIC samples and clustered both on samples and probe sets. The parameters used were Spearman ranking and average linkage. PCA analysis of the top 1600 varied probes was done with the R software package, after background removal, and the first and second components were scatter-plotted. Single-sample gene set–enrichment analysis (ssGSEA) scores were calculated by using previously reported methods [24]. Briefly, genes were ordered by their expression values. The empirical cumulative distribution functions (ECDF) of genes in the signature and in the remaining genes were calculated. The enrichment score was defined as the sum of the difference between the weighted ECDF of the genes in the signature and the ECDF of the remaining genes. This calculation was reported for all signatures and samples. A positive score indicates a high ranking of up-genes and a low ranking of down-genes in the signature. A negative score does not indicate the opposite, but rather lack of effect. Gene sets for four GBM expression subtypes were obtained from Verhaak et al. [24]. The gene signatures for six different neural lineages were obtained from Cahoy et al. [27]. Lineage cell-related gene sets were generated by using the transcriptome database presented in Cahoy et al. [25]. A positive enrichment score indicates a positive correlation between genes in the gene set and the expression of profile tumor samples; a negative enrichment score indicates an inverse correlation.

4.3. Somatic Mutation Analysis

Genomic DNA from all GIC samples was purified and subjected to phi29 polymerase multiple strand-displacement whole-genome amplification with the Sequeon OncoCarta Panel [39]. All analyses were performed in duplicate and repeated to measure assay reproducibility. Exon sequencing was done for the p53 and PTEN genes. We selected only nonsynonymous coding mutations that were previously reported as somatic mutations in human cancer, as described in the Cosmic, PubMed, and TCGA datasets.

4.4. Immunofluorescence Staining

Immunofluorescence staining was performed as described previously [40]. Cells were seeded at a concentration of 2 × 105 cells/well in six-well plates with polylysine-coated coverslips inside and left for 24 h in the incubator. The cells were washed once with phosphate-buffered saline (PBS) (Sigma-Aldrich, St. Louis, MO, USA) before being fixed with 4% formaldehyde in PBS for 20 min. After another PBS wash, the cells were permeabilized with 0.1% Triton X-100 in PBS for 5 minutes and then blocked with 5% goat serum in PBS (containing 0.1% Tween 20) for 1 h. Cells were then incubated with primary antibodies (anti-Nestin and anti-GFAP, anti–YKL-40 from Cell Signaling, anti-CNPase from Millipore, anti-TuJ-1 from Sigma-Aldrich, and anti-Olig2 from Abcam) overnight at 4 °C. After 2 washes with PBS (containing 0.1% Tween 20), the cells were incubated with the indicated fluorescently labeled secondary antibody in darkness at room temperature for 1 h. The cells were counterstained with Vectashield sealant containing 4′,6-diamidino-2-phenylindole (DAPI) (Vector Laboratories, CA, USA) and examined under a confocal laser scanning microscope (Carl Zeiss Microscopy, Germany). Images were scanned with a DFC340FX camera (Leica Microsystems, Germany). The percentage of positively stained cells in each section was determined at a magnification of 400× by counting 500 cells in a randomly selected field.

4.5. Immunohistochemical Staining

Sections (thickness 5 μm) of formalin-fixed, paraffin-embedded whole brain from animal specimens were stained with the primary antibodies listed in Table 3. The sections were visualized by using a diaminobenzidine substrate kit. The slides were examined under a bright-field microscope [10].

4.6. ELISA

We used a human VEGF enzyme-linked immunosorbent assay (ELISA) to quantify secretory VEGF and MMP-9 in conditioned media, in accordance with the manufacturer’s instructions (R&D Systems, Minneapolis, MN, USA). For analysis of MMP-9 and VEGF, 1 × 105 GICs were plated in a 12-well plate for 48 h, and VEGF was quantified in the conditioned media, as described previously. Data analysis was performed with Microplate Manager Software (Molecular Device, BioRad, Hercules, CA, USA). The VEGF/MMP-9 concentration in the samples was calculated by interpolation from the standard curve. Total protein was determined, to normalize ELISA.

4.7. Neurosphere Formation Assay

For the primary neurosphere formation assay, GICs were seeded in 96-well plates (100 cells/well) and incubated for 7 days. For the secondary neurosphere formation assay, primary neurospheres were dissociated into single cells, 50 of which were seeded in each well of 96-well plates in the absence of inhibitors. After 7 days, total neurosphere numbers were counted with a dissection microscope using 4 × 4 objective magnification.

4.8. Animal Models

The mice were housed and cared for at the animal care facility of Toho University School of Medicine in accordance with the institution’s guidelines for the care and use of laboratory animals. The experimental protocol was approved by the Animal Research Committee, Toho University School of Medicine (ARC/TUSM-R16-14). To create the intracranial xenograft model, we engrafted GICs into the caudate nucleus of each mouse, using a previously described guide-screw system [10]. Briefly, 5 × 105 dissociated tumor cells in DMEM/F-12 serum-free media (5 µM) were injected stereotactically into Nude (nu/nu) 6–8-week-old mice. Mice were monitored daily and euthanized when they became moribund. At necropsy, all organs were analyzed grossly and microscopically.

4.9. Statistical Analysis

Statistical analysis was performed by using the Student unpaired t-test. Results are presented as the mean of at least three independent experiments. Survival analysis was performed by using the log-rank analysis module in SPSS 10.0 (SPSS Inc, Chicago, IL, USA). Differences were considered significant at p < 0.05, for all comparisons.

5. Conclusions

In conclusion, the present comprehensive analysis revealed distinct GIC subtypes that are closely related to TCGA GBM. They appear to derive from different cells of the neural lineage through varied molecular processes, such as gene mutation and aberrant signaling activation. Importantly, the lineage profiles and biological phenotypes of GIC subtypes were retained by in vivo xenografts. These findings improve our understanding of GBM development and the heterogeneity of these tumors and may hasten development of more personalized and improved therapies targeting GICs. Thus, GIC is biologically and molecularly a more relevant model system for preclinical studies of therapeutic intervention and for improving our understanding of the molecular biology of human GBM.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-6694/11/10/1564/s1, Table S1: Background information on GICs, Table S2: Summary of GIC subtypes differentiation potential in vitro. GIC subtypes cultured under differentiation condition (1% FBS + 1 µM RA) exhibited differential potential of differentiation into neural lineage.

Author Contributions

N.S. designed and performed experiments, collected, analyzed, and interpreted data, and wrote the manuscript; N.H., K.A., S.S., R.S., Y.H., S.F., H.N., M.H., and T.S. performed experiments and interpreted data; S.I. supervised the study.

Funding

This work was supported by JSPS KAKENHI Grant Number 15K19980.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chang, S.M.; Theodosopoulos, P.; Lamborn, K.; Malec, M.; Rabbitt, J.; Page, M.; Prados, M.D. Temozolomide in the treatment of recurrent malignant glioma. Cancer 2004, 100, 605–611. [Google Scholar] [CrossRef] [PubMed]
  2. Wen, P.Y.; Kesari, S. Malignant gliomas in adults. N. Engl. J. Med. 2008, 359, 492–507. [Google Scholar] [CrossRef] [PubMed]
  3. Stupp, R.; Mason, W.P.; van den Bent, M.J.; Weller, M.; Fisher, B.; Taphoorn, M.J.; Belanger, K.; Brandes, A.A.; Marosi, C.; Bogdahn, U.; et al. Radiotherapy plus concomitant and adjuvant temozolomide for glioblastoma. N. Engl. J. Med. 2005, 352, 987–996. [Google Scholar] [CrossRef] [PubMed]
  4. Purow, B.; Schiff, D. Advances in the genetics of glioblastoma: Are we reaching critical mass? Nat. Rev. Neurol. 2009, 5, 419–426. [Google Scholar] [CrossRef]
  5. Singh, S.K.; Hawkins, C.; Clarke, I.D.; Squire, J.A.; Bayani, J.; Hide, T.; Henkelman, R.M.; Cusimano, M.D.; Dirks, P.B. Identification of human brain tumour initiating cells. Nature 2004, 432, 396–401. [Google Scholar] [CrossRef]
  6. Kondo, T.; Setoguchi, T.; Taga, T. Persistence of a small subpopulation of cancer stem-like cells in the C6 glioma cell line. Proc. Natl. Acad. Sci. USA 2004, 101, 781–786. [Google Scholar] [CrossRef] [Green Version]
  7. Ikushima, H.; Todo, T.; Ino, Y.; Takahashi, M.; Saito, N.; Miyazawa, K.; Miyazono, K. Glioma-initiating cells retain their tumorigenicity through integration of the Sox axis and Oct4 protein. J. Biol. Chem. 2011, 286, 41434–41441. [Google Scholar] [CrossRef]
  8. Bao, S.; Wu, Q.; McLendon, R.E.; Hao, Y.; Shi, Q.; Hjelmeland, A.B.; Dewhirst, M.W.; Bigner, D.D.; Rich, J.N. Glioma stem cells promote radioresistance by preferential activation of the DNA damage response. Nature 2006, 444, 756–760. [Google Scholar] [CrossRef]
  9. Liu, G.; Yuan, X.; Zeng, Z.; Tunici, P.; Ng, H.; Abdulkadir, I.R.; Lu, L.; Irvin, D.; Black, K.L.; Yu, J.S. Analysis of gene expression and chemoresistance of CD133+ cancer stem cells in glioblastoma. Mol. Cancer 2006, 5, 67. [Google Scholar] [CrossRef]
  10. Saito, N.; Fu, J.; Zheng, S.; Yao, J.; Wang, S.; Liu, D.D.; Yuan, Y.; Sulman, E.P.; Lang, F.F.; Colman, H.; et al. A high Notch pathway activation predicts response to gamma secretase inhibitors in proneural subtype of glioma tumor-initiating cells. Stem Cells 2014, 32, 301–312. [Google Scholar] [CrossRef]
  11. Saito, N.; Hirai, N.; Aoki, K.; Suzuki, R.; Fujita, S.; Nakayama, H.; Hayashi, M.; Ito, K.; Sakurai, T.; Iwabuchi, S. The oncogene addiction switch from NOTCH to PI3K requires simultaneous targeting of NOTCH and PI3K pathway inhibition in glioblastoma. Cancers 2019, 11, e121. [Google Scholar] [CrossRef] [PubMed]
  12. Beier, D.; Hau, P.; Proescholdt, M.; Lohmeier, A.; Wischhusen, J.; Oefner, P.J.; Aigner, L.; Brawanski, A.; Bogdahn, U.; Beier, C.P. CD133(+) and CD133(−) glioblastoma-derived cancer stem cells show differential growth characteristics and molecular profiles. Cancer Res. 2007, 67, 4010–4015. [Google Scholar] [CrossRef] [PubMed]
  13. Gunther, H.S.; Schmidt, N.O.; Phillips, H.S.; Kemming, D.; Kharbanda, S.; Soriano, R.; Modrusan, Z.; Meissner, H.; Westphal, M.; Lamszus, K. Glioblastoma-derived stem cell-enriched cultures form distinct subgroups according to molecular and phenotypic criteria. Oncogene 2008, 27, 2897–2909. [Google Scholar] [CrossRef] [PubMed]
  14. Lottaz, C.; Beier, D.; Meyer, K.; Kumar, P.; Hermann, A.; Schwarz, J.; Junker, M.; Oefner, P.J.; Bogdahn, U.; Wischhusen, J.; et al. Transcriptional profiles of CD133+ and CD133− glioblastoma-derived cancer stem cell lines suggest different cells of origin. Cancer Res. 2010, 70, 2030–2040. [Google Scholar] [CrossRef] [PubMed]
  15. Massague, J. TGFbeta in cancer. Cell 2008, 134, 215–230. [Google Scholar] [CrossRef] [PubMed]
  16. Bruna, A.; Darken, R.S.; Rojo, F.; Ocana, A.; Penuelas, S.; Arias, A.; Paris, R.; Tortosa, A.; Mora, J.; Baselga, J.; et al. High TGFbeta-Smad activity confers poor prognosis in glioma patients and promotes cell proliferation depending on the methylation of the PDGF-B gene. Cancer Cell 2007, 11, 147–160. [Google Scholar] [CrossRef] [PubMed]
  17. Penuelas, S.; Anido, J.; Prieto-Sanchez, R.M.; Folch, G.; Barba, I.; Cuartas, I.; Garcia-Dorado, D.; Poca, M.A.; Sahuquillo, J.; Baselga, J.; et al. TGF-beta increases glioma-initiating cell self-renewal through the induction of LIF in human glioblastoma. Cancer Cell 2009, 15, 315–327. [Google Scholar] [CrossRef]
  18. Lino, M.M.; Merlo, A.; Boulay, J.L. Notch signaling in glioblastoma: A developmental drug target? BMC Med. 2010, 8, e72. [Google Scholar] [CrossRef]
  19. Fan, X.; Khaki, L.; Zhu, T.S.; Soules, M.E.; Talsma, C.E.; Gul, N.; Koh, C.; Zhang, J.; Li, Y.M.; Maciaczyk, J.; et al. NOTCH pathway blockade depletes CD133-positive glioblastoma cells and inhibits growth of tumor neurospheres and xenografts. Stem Cells 2010, 28, 5–16. [Google Scholar]
  20. Lee, J.; Kotliarova, S.; Kotliarov, Y.; Li, A.; Su, Q.; Donin, N.M.; Pastorino, S.; Purow, B.W.; Christopher, N.; Zhang, W.; et al. Tumor stem cells derived from glioblastomas cultured in bFGF and EGF more closely mirror the phenotype and genotype of primary tumors than do serum-cultured cell lines. Cancer Cell 2006, 9, 391–403. [Google Scholar] [CrossRef] [Green Version]
  21. Zheng, H.; Ying, H.; Wiedemeyer, R.; Yan, H.; Quayle, S.N.; Ivanova, E.V.; Paik, J.H.; Zhang, H.; Xiao, Y.; Perry, S.R.; et al. PLAGL2 regulates Wnt signaling to impede differentiation in neural stem cells and gliomas. Cancer Cell 2010, 17, 497–509. [Google Scholar] [CrossRef]
  22. Clement, V.; Sanchez, P.; de Tribolet, N.; Radovanovic, I.; Ruiz i Altaba, A. HEDGEHOG-GLI1 signaling regulates human glioma growth, cancer stem cell self-renewal, and tumorigenicity. Curr. Biol. 2007, 17, 165–172. [Google Scholar] [CrossRef] [PubMed]
  23. Phillips, H.S.; Kharbanda, S.; Chen, R.; Forrest, W.F.; Soriano, R.H.; Wu, T.D.; Misra, A.; Nigro, J.M.; Colman, H.; Soroceanu, L.; 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, 157–173. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Verhaak, R.G.; Hoadley, K.A.; Purdom, E.; Wang, V.; Qi, Y.; Wilkerson, M.D.; Miller, C.R.; Ding, L.; Golub, T.; Mesirov, J.P.; et al. Integrated genomic analysis identifies clinically relevant subtypes of glioblastoma characterized by abnormalities in PDGFRA, IDH1, EGFR, and NF1. Cancer Cell 2010, 17, 98–110. [Google Scholar] [CrossRef]
  25. Cahoy, J.D.; Emery, B.; Kaushal, A.; Foo, L.C.; Zamanian, J.L.; Christopherson, K.S.; Xing, Y.; Lubischer, J.L.; Krieg, P.A.; Krupenko, S.A.; et al. A transcriptome database for astrocytes, neurons, and oligodendrocytes: A new resource for understanding brain development and function. J. Neurosci. 2008, 28, 264–278. [Google Scholar] [CrossRef]
  26. Chu, T.; Rice, E.J.; Booth, G.T.; Salamanca, H.H.; Wang, Z.; Core, L.J.; Longo, S.L.; Corona, R.J.; Chin, L.S.; Lis, J.T.; et al. Chromatin run-on and sequencing maps the transcriptional regulatory landscape of glioblastoma multiforme. Nat. Genet. 2018, 50, 1553–1564. [Google Scholar] [CrossRef] [PubMed]
  27. Bleeker, F.E.; Lamba, S.; Zanon, C.; van Tilborg, A.A.; Leenstra, S.; Troost, D.; Hulsebos, T.; Vandertop, W.P.; Bardelli, A. Absence of AKT1 mutations in glioblastoma. PLoS ONE 2009, 4, e5638. [Google Scholar] [CrossRef]
  28. Turner, K.M.; Sun, Y.; Ji, P.; Granberg, K.J.; Bernard, B.; Hu, L.; Cogdell, D.E.; Zhou, X.; Yli-Harja, O.; Nykter, M.; et al. Genomically amplified Akt3 activates DNA repair pathway and promotes glioma progression. Proc. Natl. Acad. Sci. USA 2015, 112, 3421–3426. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Mooney, K.L.; Choy, W.; Sidhu, S.; Pelargos, P.; Bui, T.T.; Voth, B.; Barnette, N.; Yang, I. The role of CD44 in glioblastoma multiforme. J. Clin. Neurosci. 2016, 34, 1–5. [Google Scholar] [CrossRef]
  30. Crespo, I.; Vital, A.L.; Gonzalez-Tablas, M.; Patino Mdel, C.; Otero, A.; Lopes, M.C.; de Oliveira, C.; Domingues, P.; Orfao, A.; Tabernero, M.D. Molecular and genomic alterations in glioblastoma multiforme. Am. J. Pathol. 2015, 185, 1820–1833. [Google Scholar] [CrossRef]
  31. Ohgaki, H.; Kleihues, P. The definition of primary and secondary glioblastoma. Clin. Cancer Res. 2013, 19, 764–772. [Google Scholar] [CrossRef] [PubMed]
  32. Tchoghandjian, A.; Baeza, N.; Colin, C.; Cayre, M.; Metellus, P.; Beclin, C.; Ouafik, L.; Figarella-Branger, D. A2B5 cells from human glioblastoma have cancer stem cell properties. Brain Pathol. 2010, 20, 211–221. [Google Scholar] [CrossRef]
  33. Sugiarto, S.; Persson, A.I.; Munoz, E.G.; Waldhuber, M.; Lamagna, C.; Andor, N.; Hanecker, P.; Ayers-Ringler, J.; Phillips, J.; Siu, J.; et al. Asymmetry-defective oligodendrocyte progenitors are glioma precursors. Cancer Cell 2011, 20, 328–340. [Google Scholar] [CrossRef] [PubMed]
  34. Plaks, V.; Kong, N.; Werb, Z. The cancer stem cell niche: How essential is the niche in regulating stemness of tumor cells? Cell Stem Cell 2015, 16, 225–238. [Google Scholar] [CrossRef]
  35. Jin, X.; Kim, L.J.Y.; Wu, Q.; Wallace, L.C.; Prager, B.C.; Sanvoranart, T.; Gimple, R.C.; Wang, X.; Mack, S.C.; Miller, T.E.; et al. Targeting glioma stem cells through combined BMI1 and EZH2 inhibition. Nat. Med. 2017, 23, 1352–1361. [Google Scholar] [CrossRef]
  36. Shingu, T.; Ho, A.L.; Yuan, L.; Zhou, X.; Dai, C.; Zheng, S.; Wang, Q.; Zhong, Y.; Chang, Q.; Horner, J.W.; et al. Qki deficiency maintains stemness of glioma stem cells in suboptimal environment by downregulating endolysosomal degradation. Nat. Genet. 2017, 49, 75–86. [Google Scholar] [CrossRef] [PubMed]
  37. Bengtsson, H.; Jonsson, G.; Vallon-Christersson, J. Calibration and assessment of channel-specific biases in microarray data with extended dynamical range. BMC Bioinform. 2004, 5, e177. [Google Scholar] [CrossRef]
  38. Bhat, K.P.; Salazar, K.L.; Balasubramaniyan, V.; Wani, K.; Heathcock, L.; Hollingsworth, F.; James, J.D.; Gumin, J.; Diefes, K.L.; Kim, S.H.; et al. The transcriptional coactivator TAZ regulates mesenchymal differentiation in malignant glioma. Genes. Dev. 2011, 25, 2594–2609. [Google Scholar] [CrossRef] [Green Version]
  39. Pearce, M.; Cullinan, A.; Hogg, G.; Hosseini, D.; Ehrich, M. Mutation profiling in tumor samples using the Sequenom OncoCarta™ Panel. Nat. Methods 2009, 6, 7–8. [Google Scholar] [CrossRef]
  40. Liu, J.L.; Sheng, X.; Hortobagyi, Z.K.; Mao, Z.; Gallick, G.E.; Yung, W.K. Nuclear PTEN-mediated growth suppression is independent of Akt down-regulation. Mol. Cell Biol. 2005, 25, 6211–6224. [Google Scholar] [CrossRef]
Cancers 11 01564 g001 550
Figure 1. (a,b): Supervised clustering (The Cancer Genome Atlas (TCGA) 1461 probe sets) and unsupervised clustering (MTA 1600 probe sets) classified 25 glioma-initiating cell (GIC) cell lines into mesenchymal, neural, proneural, and classical subtypes. Defined GIC subtypes could be distinguished by their distinct patterns of gene expression, as related to TCGA subclass. (c) There were 275 probe sets (~20% of total) common to these two methods. (d) PCAof gene expression of the four GIC subtypes indicated that the proneural and classical GICs were most closely related and that neural and mesenchymal GICs deviated in a direction opposite to them, which supports the hypothesis that the “proneural” group has two GIC subtypes.
Figure 1. (a,b): Supervised clustering (The Cancer Genome Atlas (TCGA) 1461 probe sets) and unsupervised clustering (MTA 1600 probe sets) classified 25 glioma-initiating cell (GIC) cell lines into mesenchymal, neural, proneural, and classical subtypes. Defined GIC subtypes could be distinguished by their distinct patterns of gene expression, as related to TCGA subclass. (c) There were 275 probe sets (~20% of total) common to these two methods. (d) PCAof gene expression of the four GIC subtypes indicated that the proneural and classical GICs were most closely related and that neural and mesenchymal GICs deviated in a direction opposite to them, which supports the hypothesis that the “proneural” group has two GIC subtypes.
Cancers 11 01564 g001
Cancers 11 01564 g002 550
Figure 2. (a) GIC subtypes differentially expressed lineage markers for neural stem/progenitor cells (neuroepithelial cells, radial glia, glial progenitors, neuronal progenitors, oligodendrocyte progenitor cells, and SVZ astrocytes). (b) GIC subtypes exhibited distinct lineage profiles that recapitulated neural development. Earlier lineage markers (Nestin, Sox2, Olig2) are highly expressed in proneural and classical GICs, while differentiation markers (GFAP: glial fibrillary acidic protein, TuJ1: Neuron-specific beta-III Tubulin) are abundant in neural GICs. YKL-40: Chitinase-3-like protein 1 is enriched in mesenchymal GICs. (c): Lineage characteristics of GIC subtypes are retained in their intracranial xenografts. One of representative results of each subtypes are illustrated in Figure 2C. Scale bars: 100 μm.
Figure 2. (a) GIC subtypes differentially expressed lineage markers for neural stem/progenitor cells (neuroepithelial cells, radial glia, glial progenitors, neuronal progenitors, oligodendrocyte progenitor cells, and SVZ astrocytes). (b) GIC subtypes exhibited distinct lineage profiles that recapitulated neural development. Earlier lineage markers (Nestin, Sox2, Olig2) are highly expressed in proneural and classical GICs, while differentiation markers (GFAP: glial fibrillary acidic protein, TuJ1: Neuron-specific beta-III Tubulin) are abundant in neural GICs. YKL-40: Chitinase-3-like protein 1 is enriched in mesenchymal GICs. (c): Lineage characteristics of GIC subtypes are retained in their intracranial xenografts. One of representative results of each subtypes are illustrated in Figure 2C. Scale bars: 100 μm.
Cancers 11 01564 g002
Cancers 11 01564 g003 550
Figure 3. GIC subtypes cultured under different conditions (1% (fetal bovine serum) FBS + 1 µM RA) exhibited varied potential for differentiation into neural lineages. Most proneural and classical GICs displayed trilineage differentiation potential, while neural GICs rarely differentiated into the oligodendrocytic lineage. Mesenchymal GICs were less likely to differentiate into astrocytic, neuronal, and oligodendrocytic lineages. One of representative results of each subtypes are illustrated in Figure 3. Scale bars: 100 μm.
Figure 3. GIC subtypes cultured under different conditions (1% (fetal bovine serum) FBS + 1 µM RA) exhibited varied potential for differentiation into neural lineages. Most proneural and classical GICs displayed trilineage differentiation potential, while neural GICs rarely differentiated into the oligodendrocytic lineage. Mesenchymal GICs were less likely to differentiate into astrocytic, neuronal, and oligodendrocytic lineages. One of representative results of each subtypes are illustrated in Figure 3. Scale bars: 100 μm.
Cancers 11 01564 g003
Cancers 11 01564 g004 550
Figure 4. (a,b): GIC subtypes exhibit distinct biological behaviors in relation to self-renewal capacity. Representative results of each subtypes are illustrated in Figure 4a,b. Scale bars: 250 μm.
Figure 4. (a,b): GIC subtypes exhibit distinct biological behaviors in relation to self-renewal capacity. Representative results of each subtypes are illustrated in Figure 4a,b. Scale bars: 250 μm.
Cancers 11 01564 g004
Cancers 11 01564 g005 550
Figure 5. (a) The in vivo biological behavior of GIC subtypes was studied by orthotopic injection of cells into mouse brain. The neural xenograft subtype was the least angiogenic, as indicated by vWF staining. Ki67 staining showed that classical GICs were most active in proliferation in vivo. Scale bars: 100 μm. (b) Neural GICs produced significantly lower amounts of VEGF and MMP-9.
Figure 5. (a) The in vivo biological behavior of GIC subtypes was studied by orthotopic injection of cells into mouse brain. The neural xenograft subtype was the least angiogenic, as indicated by vWF staining. Ki67 staining showed that classical GICs were most active in proliferation in vivo. Scale bars: 100 μm. (b) Neural GICs produced significantly lower amounts of VEGF and MMP-9.
Cancers 11 01564 g005
Cancers 11 01564 g006 550
Figure 6. GIC subtypes exhibited divergent patterns of signaling pathway activation. Multiple pathways, such as TGF-β, Notch, VEGF, and Wnt, were identified in GIC subgroups.
Figure 6. GIC subtypes exhibited divergent patterns of signaling pathway activation. Multiple pathways, such as TGF-β, Notch, VEGF, and Wnt, were identified in GIC subgroups.
Cancers 11 01564 g006
Table
Table 1. Mutation analysis in GICs.
Table 1. Mutation analysis in GICs.
SubtypeCell LineTP53PTENPIK3CAPIK3R1METIDH1N-RasAKT3
NeuralGICT14S241Y R132H E17K
GICT12G105DDel exons 2–9
GICT46 Y88 *
GICT20P98L
GICT24R306 stopcodon
GICT35G105DDel exons 2–9
ProneuralGICT11 Q546P
GICT22 Q546P
GICT14 R88Q
GICT25R273C
GICT18 Exon 6Del AGAA
GICT28 Y68H M326I
ClassicalGICT15M237IDel exons 1
GICT12 Del exons 1–2 M326I
GICT32 M326IR988C
GICT23 R233 * N564K
GICT40C238Y N375S
GICT19 Deletion all exons
MesenchymalGICT17H179RM199del Q61K
GICT38 Del exons 3–9
GICT29C238Y N375S
GICT44P98L
GICT22 R233 *
GICT16 G132D
GICT21 Del exons 2–9
TP53: tumor protein p53; PTEN: phosphatase and tensin homolog; PIK3CA: phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha PIK3R1: phosphoinositide-3-kinase regulatory subunit 1; MET: MET proto-oncogene, receptor tyrosine kinase; IDH1: isocitrate dehydrogenase 1; N-Ras: NRAS proto-oncogene; AKT3: AKT serine/threonine kinase 3; *:Nonsense mutation.
Table
Table 2. The histological phenotypes of tumors formed differed in relation to GIC subtype.
Table 2. The histological phenotypes of tumors formed differed in relation to GIC subtype.
Histological PhenotypeNeuralProneuralClassicalMesenchymal
In vivo invasionLess invasiveinvasiveinvasiveHighly invasive
In vivo angiogenesisLess angiogenicLess angiogenicAngiogenicHighly angiogenic
necrosisFocal, palisading cellsFocal, palisading cellsExtensive, extravasating blood cellsMild
Proliferation (Ki67 staining)LowLowhighLow
Table
Table 3. Primary antibodies list.
Table 3. Primary antibodies list.
Antigen (Clone/Code)SourceAntibody
GLASTAbcamRabbit
NestinCell SignalingMouse
Olig2AbcamRabbit
Sox2Sigma-AldrichRabbit
PDGFRaSanta CruzRabbit
S100bAbcamRabbit
GFAPCell SignalingRabbit
TuJ1Sigma-AldrichRabbit
PSA-NCAMSigma-AldrichMouse
O4Sigma-AldrichMouse
NG2AbcamRabbit
YKL-40Cell SignalingRabbit
A2B5AbcamMouse
CD133AbcamRabbit
Ki67Cell SignalingRabbit
vWFAbcamRabbit

Share and Cite

MDPI and ACS Style

Saito, N.; Hirai, N.; Aoki, K.; Sato, S.; Suzuki, R.; Hiramoto, Y.; Fujita, S.; Nakayama, H.; Hayashi, M.; Sakurai, T.; et al. Genetic and Lineage Classification of Glioma-Initiating Cells Identifies a Clinically Relevant Glioblastoma Model. Cancers 2019, 11, 1564. https://doi.org/10.3390/cancers11101564

AMA Style

Saito N, Hirai N, Aoki K, Sato S, Suzuki R, Hiramoto Y, Fujita S, Nakayama H, Hayashi M, Sakurai T, et al. Genetic and Lineage Classification of Glioma-Initiating Cells Identifies a Clinically Relevant Glioblastoma Model. Cancers. 2019; 11(10):1564. https://doi.org/10.3390/cancers11101564

Chicago/Turabian Style

Saito, Norihiko, Nozomi Hirai, Kazuya Aoki, Sho Sato, Ryo Suzuki, Yu Hiramoto, Satoshi Fujita, Haruo Nakayama, Morito Hayashi, Takatoshi Sakurai, and et al. 2019. "Genetic and Lineage Classification of Glioma-Initiating Cells Identifies a Clinically Relevant Glioblastoma Model" Cancers 11, no. 10: 1564. https://doi.org/10.3390/cancers11101564

APA Style

Saito, N., Hirai, N., Aoki, K., Sato, S., Suzuki, R., Hiramoto, Y., Fujita, S., Nakayama, H., Hayashi, M., Sakurai, T., & Iwabuchi, S. (2019). Genetic and Lineage Classification of Glioma-Initiating Cells Identifies a Clinically Relevant Glioblastoma Model. Cancers, 11(10), 1564. https://doi.org/10.3390/cancers11101564

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop