Next Article in Journal
Understanding the Biological Activities of Vitamin D in Type 1 Neurofibromatosis: New Insights into Disease Pathogenesis and Therapeutic Design
Next Article in Special Issue
Enhanced Malignant Phenotypes of Glioblastoma Cells Surviving NPe6-Mediated Photodynamic Therapy are Regulated via ERK1/2 Activation
Previous Article in Journal
CTEN Induces Tumour Cell Invasion and Survival and Is Prognostic in Radiotherapy-Treated Head and Neck Cancer
Previous Article in Special Issue
Pembrolizumab Activity in Recurrent High-Grade Gliomas with Partial or Complete Loss of Mismatch Repair Protein Expression: A Monocentric, Observational and Prospective Pilot Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Immunohistochemically Characterized Intratumoral Heterogeneity Is a Prognostic Marker in Human Glioblastoma

1
Department of Neuropathology, School of Medicine, Institute of Pathology, Technical University Munich, Trogerstraße 18, 81675 München, Germany
2
Department of Diagnostic and Interventional Neuroradiology, School of Medicine, Technical University Munich, Ismaninger Str. 22, 81675 München, Germany
3
Department of Neurology, School of Medicine, Technical University Munich, Ismaninger Str. 22, 81675 München, Germany
4
Institute of Pathology, School of Medicine, Technical University Munich, Trogerstraße 18, 81675 München, Germany
5
Department of Neurosurgery, School of Medicine, Technical University Munich, Ismaninger Str. 22, 81675 München, Germany
6
Department of RadiationOncology, School of Medicine, Technical University Munich, Ismaninger Str. 22, 81675 München, Germany
7
TranslaTUM (Zentralinstitut für translationale Krebsforschung der Technischen Universität München), Einsteinstraße 25, 81675 München, Germany
*
Author to whom correspondence should be addressed.
These authors contributed to this paper equally.
These authors contributed to this paper equally.
Cancers 2020, 12(10), 2964; https://doi.org/10.3390/cancers12102964
Submission received: 7 September 2020 / Revised: 2 October 2020 / Accepted: 9 October 2020 / Published: 13 October 2020
(This article belongs to the Special Issue Recurrent Glioblastoma)

Abstract

:

Simple Summary

Intratumoral heterogeneity is believed to contribute to the immense therapy resistance and recurrence rate of glioblastoma. The aim of this retrospective study was to analyze the heterogeneity of 36 human glioblastoma samples on a morphological level by immunohistochemistry. We confirmed that this method is valid for heterogeneity detection. 115 Areas of Interest were labelled. By cluster analysis, we defined two subtypes (“classical” and “mesenchymal”). The results of epigenomic analyses corroborated the findings. Interestingly, patients with tumors that consisted of both subtypes (“subtype-heterogeneous”) showed a shorter overall survival compared to patients with tumor that were dominated by one subtype (“subtype-dominant”). Furthermore, the analysis of 21 corresponding pairs of primary and recurrent glioblastoma demonstrated that, additionally to an intratumoral heterogeneity, there is also a chronological heterogeneity with dominance of the mesenchymal subtype in recurrent tumors. Our study confirms the prognostic impact of intratumoral heterogeneity in glioblastoma and makes this hallmark assessable by routine diagnostics.

Abstract

Tumor heterogeneity is considered to be a hallmark of glioblastoma (GBM). Only more recently, it has become apparent that GBM is not only heterogeneous between patients (intertumoral heterogeneity) but more importantly, also within individual patients (intratumoral heterogeneity). In this study, we focused on assessing intratumoral heterogeneity. For this purpose, the heterogeneity of 38 treatment-naïve GBM was characterized by immunohistochemistry. Perceptible areas were rated for ALDH1A3, EGFR, GFAP, Iba1, Olig2, p53, and Mib1. By clustering methods, two distinct groups similar to subtypes described in literature were detected. The classical subtype featured a strong EGFR and Olig2 positivity, whereas the mesenchymal subtype displayed a strong ALDH1A3 expression and a high fraction of Iba1-positive microglia. 18 tumors exhibited both subtypes and were classified as “subtype-heterogeneous”, whereas the areas of the other tumors were all assigned to the same cluster and named “subtype-dominant”. Results of epigenomic analyses corroborated these findings. Strikingly, the subtype-heterogeneous tumors showed a clearly shorter overall survival compared to subtype-dominant tumors. Furthermore, 21 corresponding pairs of primary and recurrent GBM were compared, showing a dominance of the mesenchymal subtype in the recurrent tumors. Our study confirms the prognostic impact of intratumoral heterogeneity in GBM, and more importantly, makes this hallmark assessable by routine diagnostics.

1. Introduction

Tumor heterogeneity is now considered as a hallmark of glioblastoma (GBM). It is believed to strongly contribute to therapy resistance and accordingly, to the poor prognosis of this tumor entity [1]. To understand the various pheno- and genotypical characteristics and their biology, numerous approaches to subclassify glioblastomas exist. Originally, they were developed to describe intertumoral heterogeneity. Subtypes, called “classical”, “mesenchymal”, “proneural”, and “neural” have been consistently discussed since Verhaak et al. stated that the subtypes show different clinical courses and biology [2]. However, it has become apparent that, even more importantly, heterogeneity within a single tumor exists: Several studies have demonstrated that one tumor can harbor multiple subclones, which are assigned to different subtypes by their molecular characteristics [3,4]. Similar to evolutional processes, diversity leads to advantage. There are different hypotheses regarding the development of heterogeneity in tumors. Based on Charles Darwin’s theory of evolution [5], the hypothesis of clonal evolution sees heterogeneity as a result of natural selection [6]. Genetic instability of tumor cells results in accumulated mutations leading to genetic diversity and heterogeneous morphology. By selective stress, e.g., caused by chemo- or radiotherapy, only adequately adapted cell clones survive [7]. In contrast, the stem cell model acts on the assumption of a hierarchical organization of tumor cells [8]. By self-renewal of stem cell like neoplastic cells, genetically and phenotypically diverse daughter cells develop, from which different intratumoral subtypes arise [9,10]. Besides these two main hypotheses, it is assumed that tumor heterogeneity is a consequence of a multifactorial process, including epigenetic alterations [11], intercellular communication, and interaction with the surrounding microenvironment [10,12]. In addition to regional heterogeneity, a chronological heterogeneity can also be observed when comparing pairs of primary and recurrent GBM. The mesenchymal subtype seems to be most the therapy resistant, since its occurrence increases in recurrent tumors [13]. Most studies on heterogeneity of GBM are based on large-scale genomic characterization. This is a powerful tool for discovery and in-depth tumor analysis, but it has a limited availability. There have been previous attempts to capture heterogeneity in GBM by immunohistochemistry, but these studies focused mainly on intertumoral heterogeneity, though [14,15]. Just recently, we published our first study with a morphological approach by using immunohistochemistry. We defined different tumor regions including region of hypoxia and stem cell region [16]. In this current study, we, again, chose the broadly applicable technique immunohistochemistry, but focused on applying the established subtypes on human tumor tissue with the aim to prove that immunohistochemistry is a valid method for detecting these diverse subtypes in an individual tumor. Furthermore, we hypothesized that the detection of different subtypes within one tumor has impact on its biological and clinical behavior. For this study, markers were chosen that have already been proposed for the recognition of different subtypes. Alterations of the epidermal growth factor receptor (EGFR) are very common in GBM [17,18]. The status of EGFR amplification correlates with the tumor’s potential to migrate [19]. The upregulation of EGFR is a characteristic of the classical subtype, according to Verhaak et al. [2]. The glial fibrillary acidic protein (GFAP), an astrocytic intermediate filament, is associated with migration and motility of astrocytes [18]. The mouse models showed that GFAP-positive tumors act out a more aggressive growth [20]. Oligodendrocyte lineage factor 2 (Olig2) is a transcription factor regulating proliferation of stem cells in the central nervous system (CNS) [21]. In tumors, Olig2 abrogates the proliferation inhibition of tumor suppressor p21 [22]. Experimental Olig2 deletion led to a shift from a proneural to the mesenchymal GBM subtype with the abrogation of EGFR [23]. 25% of primary GBM exhibit alterations in the function of p53 [24]. Mutations of this transcription factor were considered to be characteristic for the proneural subtype of GBM by Verhaak et al. [2]. It has to be noted, though, that the proneural subtype was also defined by mutations of isocitrate dehydrogenase (IDH), which leads to the challenge of whether this assignment is still contemporary. Nevertheless, p53 is of high interest and it was included in this study. In addition, expression of the enzyme aldehyde dehydrogenase 1A3 (ALDH1A3) was analyzed. By its catalytic activity, which leads to oxidation from all-trans retinal to retinoic acid [25], it influences cell proliferation, differentiation, and apoptosis [26]. Furthermore, enzymes of the ALDH family counteract oxidative stress and therefore, protect from cell damage by aldehyde oxidation [27]. As tumor marker, ALDH1A3 is associated with poor outcome in a diversity of malignant tumors, amongst others also in high-grade gliomas [28]. It was shown that the enzyme is associated with the mesenchymal subtype in GBM [26]. Because not only tumor cells constitute tumor tissue, we chose to include ionized calcium binding adaptor molecule 1 (Iba1) as marker of microglial cells, which, on average, reach a fraction of 30 to 40% of all cells in the tumor area [29]. Besides immunosuppressive effects [29], microglial cells promote cell proliferation and migration by the secretion of growth factors [28]. Furthermore, a high amount of microglial cells is associated with the mesenchymal subtype of GBM [13]. Lastly, proliferation marker, molecular immunology borstel 1 (Mib1), was included in this study in order to examine whether proliferation activity was associated with certain GBM subtypes.
By staining the mentioned markers immunohistochemically, this study demonstrates the intratumoral heterogeneity in human glioblastoma samples on a regional level and by also comparing pairs of primary and recurrent GBM on a chronological level.

2. Results

2.1. Areas of Interest (AoI)

For the primary glioblastomas, a total of 115 AoI were defined, ranging from 1 to 8 AoI per tumor. A total of 80 AoI were determined in the recurrent tumors, also ranging from 1 to 8 per tumor. The sample size of primary tumors varied between 12.5 mm2 and 658.0 mm2 (mean: 376.0 mm2); in the recurrent tumors, the range was from 17.8 mm2 to 658.1 mm2 (mean: 372.8 mm2).

2.2. Correlation Analysis

Spearman-rho correlation analysis showed a significant positive correlation between ALDH1A3 expression and GFAP (r = 0.228; p = 0.014) and Iba1 positivity (r = 0.468; p = 0.000). GFAP showed strong negative correlations with EGFR (r = −0.291; p = 0,002), Olig2 (r = −0.345; p = 0.000), and proliferation marker Mib1 (r = −0.413; p = 0.000). Furthermore, Mib1 showed a positive correlation with Olig2 (r = 0.415; p = 0.000). Table 1 shows all of the values. For the actual values of immunohistochemistry, also see Table S1.

2.3. Tumor Cells Express ALDH1A3

Correlation analysis showed a strong correlation between Iba1 positivity and ALDH1A3 expression. We performed an immunofluorescence co-staining for Iba1 and ALDH1A3 to exclude that only Iba1-positive microglial cells express ALDH1A3. It can be nicely shown that many ALDH1A3-positive cells do not express Iba1 (Figure 1) and, hence, are not microglial, but assumingly tumor cells.

2.4. Cluster Analysis Defines Two Immunohistochemical Subtypes

Two clustering methods—hierarchical clustering and partitioning around medoids clustering (PAM)—were used. For both methods, best division was observed for a number of two clusters. Figure 2 shows the hierarchical clustering with the mean scores of all AoI, which were assigned to the two clusters.
The results were corroborated with PAM analysis (Figure 3). This method also detects a medoid for each cluster, i.e., the sample with the most characteristic immunohistochemical profile of each cluster. These profiles were very similar to those of the mean scores of the two clusters that were achieved by hierarchical clustering. We chose the PAM results for further analysis as (a) calculating the medoids gives an intuitive insight into the “most typical” member of each group and (b) through distance calculation, assigning new samples to one of the two clusters is very straightforward.
The two clustering methods assigned 100 of the 115 (87.0%) defined to the concordant cluster, whereas, for 15 (13.0%) areas, the attribution deviated. Figure 2 marks the divergent AoI.
Through this cluster analysis, two immunohistochemical clusters or subtypes can be described:

2.4.1. Cluster A = Classical/Proliferating Subtype

Forty-six AoI of 23 tumors were assigned to cluster A by PAM clustering. The areas of this cluster show a high proliferation index (mean Mib1 value 38.8%, range 20–90% compared to mean Mib1 value of Cluster B of 15.7%, range 5–90%), which justifies the additional name proliferating subtype in our sample collective. The cluster is further characterized by high expression of EGFR, p53, and Olig2. The positive correlation of these three markers collaborate the profile observation (see also Table 1). EGFR is a well-known marker of the classical subtype of the TCGA classification by Verhaak et al. [2]. In cluster A, the mean immune reactive score (IRS) adds up to 7.0 as compared to an average of 5.0, taken all AoI together. The strong EGFR positivity of the classical/proliferating subtype suggests a close relation to the classical subtype described by Verhaak et al. [2]. Furthermore, Olig2 shows a higher IRS in this cluster (IRS 5.9 compared to an average of 3.8). For the expression profile of the classical/proliferating subtype, see also Figure 4.

2.4.2. Cluster B = Mesenchymal/Microglial-Dominant Subtype

PAM clustering assigned 69 AoI originating from 33 tumors to cluster B, which is defined by the high expression of biomarkers ALDH1A3 (mean IRS of 4.5 as compared to an average of 4.0) and GFAP (mean IRS of 6.9 compared to an average of 5.6). The high amount of Iba1-positive microglial cells (mean IRS of 6.3 as compared to an average of 6.1) led to the additional name microglial-dominant subtype. While in cluster A, proliferation processes seem to prevail, cluster B shows microenvironmental changes. Because ALDH1A3 is a mesenchymal marker, this cluster bears the name mesenchymal subtype, as known from Verhaak’s TCGA classification [2]. For the expression profile of the microglial-dominant/mesenchymal subtype, also see Figure 4.

2.5. Cluster Analysis Shows Intratumoral Heterogeneity

With cluster analysis, we could confirm the existence of intratumoral heterogeneity, as of the 38 primary GBM, 18 (47.4%) contained AoI, which were assigned to different clusters by PAM. These tumors were named subtype-heterogeneous (ST-het). Of 15 (39.5%) GBM, all intratumoral defined AoI were assigned to the microglial-dominant/mesenchymal subtype. Five (13.2%) tumors consisted of only AoI of the classical/proliferating subtype. Tumors, whose AoI were all assigned to one cluster, were summarized to a subtype-dominant (ST-dom) group.
The mean area of ST-het tumors was not significantly larger than the mean area of ST-dom tumors (433.84 mm2 compared to 323.85 mm2, p = 0.096).

2.6. Patients with Subtype-Heterogeneous Tumors May Have a Poorer Survival Than Patients with Subtype-Dominant Tumors

The clinical relevance of the two immunohistochemical subtypes and their coexistence in one tumor is shown by the results regarding correlation with clinical data. For survival analysis, only glioblastomas with non-methylated promotor of the O-6-methylguanin-DNA methyltransferase (MGMT) were taken into account to allow for a comparable precondition. The tumor samples tested for their MGMT promotor status are marked in Figure 2. Including only tumors with non-methylated MGMT promotor status led to an analysis of 20 tumors, whereof 10 were assigned to ST-het and 10 to ST-dom. Figure 5a shows the result of the Kaplan–Meier method. Even if it did not reach level of significance (p = 0.166), a difference of the two groups can be seen. Patients with ST-het tumors show a shorter overall survival (OS) with a mean OS of 18.6 months (95 % CI 13.9 and 23.2 months) when compared to patients with ST-dom tumors with a mean OS of 25.3 months (95 % CI 16.7 and 33.8).

2.7. Epigenetic Profiles Confirm Existence of Intratumoral Heterogeneity

Immunohistochemistry demonstrated intratumoral heterogeneity on protein level. Subsequently, analyses of methylation profiles of tumor samples were conducted in order to confirm this observation. For that, DNA from each two different AoI of three tumors was extracted and examined for their methylation profile using an epigenome-wide EPIC array. For one tumor, the AoI were assigned to the different clusters defined by PAM clustering (sample Het02), whereas, for the other two, both AoI were assigned to the same cluster (Het01 to mesenchymal/microglial-dominant subtype; Het03 to classical/proliferating subtype). The included samples are also marked in Figure 2. There, it can be noted that, by hierarchical clustering, the chosen AoI of Het03 were assigned to different clusters. Despite that, this particular tumor had been chosen for epigenetic profiling due to the appropriate size of AoI for sufficient DNA extraction. Principal component analysis (PCA) of the six samples revealed an expected proximity of the both samples of tumors Het01 and Het03. Interestingly, the methylation profiling confirms the allocation of the samples of tumor Het02 to the two different clusters, as one AoI shows a distinct nearer proximity to the samples of Het01 as to the other AoI of the same tumor sample (Figure 5b).

2.8. Dominance of Mesenchymal/Microglial-Dominant Subtype in Case of Recurrence

Intratumoral heterogeneity is a local phenomenon, but it has also a temporal component. This can be demonstrated by comparing primary tumor and its relapse. Twenty-one corresponding pairs of primary and relapse were analyzed and compared regarding their cluster assignment. For the recurrent tumors, a clear dominance of the mesenchymal/microglial-dominant subtype was observed.
Seventy-four (92.5%) of the 80 AoI were assigned to this subtype by PAM clustering, whereas only six (7.5%) were grouped into the classical/proliferating subtype. When only taking in account the recurrent tumors with corresponding primary, 69 of 72 AoI (95.8%) were clustered into the mesenchymal/microglial-dominant subtype, whereas only three (4.2%) were assigned to the classical/proliferating subtype. This led to an assignment of 19 of 21 (90%) recurrent tumors to ST-dom, whereof all were mesenchymal/microglial-dominant. The remaining two (10%) tumors belonged to the ST-het group. Figure 6 demonstrates the expression levels of ALDH1A3 (a), EGFR (b), Iba1 (c), and Mib1 (d) in AoI of primary tumors when compared to recurrent tumors. Mean IRS for ALDH1A3 and Iba1, both markers for the mesenchymal/microglial-dominant subtype were much higher in the recurrent tumors compared to the primary tumors (mean IRS for ALDH1A3: 7.19 vs. 3.48; mean IRS for Iba1: 7.46 vs. 6.12). In contrast, EGFR and Mib1, both markers for the classical/proliferating subtype, decreased between primary and recurrent tumors (mean IRS for EGFR: 4.82 vs. 3.08; mean score for Mib1: 5.77 vs. 3.75). T-test analyses showed a significant increase of ALDH1A3 and Iba1 in the progression (p = 0.000 and p = 0.001) and significant decrease of EGFR and amount of proliferating cells (p = 0.005 and p = 0.000).

3. Discussion

This study confirms the existence of intratumoral heterogeneity in GBM and its influence on the clinical outcome. So far, mostly single cell analyses identified intratumoral heterogeneity on genomic level [4], but, as this study shows, morphological examination, including immunohistochemical protein labeling, also serves as valid heterogeneity detection system. Therefore, our results make assessing this central oncogenetic property of GBM broadly clinically available. Most previous immunohistochemical studies focused on intertumoral heterogeneity. Conroy et al. [14] for instance, defined three subtypes with high expression of EGFR defining the classical subtype. Popova et al. [15] also defined a mesenchymal and a classical subtype. With immunohistochemical analyses of EGFR, PDGFRA, and p53, Le Mercier et al. [30] succeeded with a division in a proneural, characterized by PDGFRA and p53 positivity and a classical subtype with high EGFR expression, as well as “others”, which do not express any of the markers. All of these results are in accordance with our findings, which go beyond that by showing that these subtypes can not only be found in different GBM, but also coexist in one tumor. A third subtype, the proneural subtype, which was proposed by several groups, was mostly defined by IDH mutation [14,15,31] and it should not be compared to IDH wildtype GBM according to the current state of knowledge. In early heterogeneity studies, a neural subtype was named [2]. Newer studies propose the influence of contamination with surrounding normal brain tissue, since this subtype was mostly found at the tumor border [13]. In a recent study, we began to apply the concept of intratumoral heterogeneity by immunohistochemistry and defined different tumor regions, e.g., region of hypoxia, proliferative region and stem cell region [16]. Although, the current study focused on applying the established subtypes on human tumor tissue. By cluster analysis, two subgroups were defined. The first subtype is comparable to the classical subtype by Verhaak et al. [2]. Our results show that it is immunohistochemically characterized by a strong positivity for EGFR and Olig2, as well as a high proliferation activity. High EGFR expression leads to angiogenesis and invasion [9,32]. A heterogenous protein expression has already been described for GBM [33], and it may influence the tumor cells’ potential to migrate [19]. Stem cell marker Olig2, which is also associated with the classical subtype in this study, is integrated in regulating stem cell proliferation [22] and driving tumor growth [34]. Contrary to our results, the transcription factor was shown to be associated with the proneural subtype [30,35], but, as mentioned before, it is questionable if this subtype exists. This might explain why an association between EGFR and Olig2 and their belonging to the same subtype have not been described before. Similar applies for p53 that was also mentioned as characteristic for a proneural subtype before [2,14], and that belongs to the classical subtype in our sample collective. The other subtype, which was defined in this study, is comparable to the mesenchymal subtype by Verhaak et al. [2]. Here, a strong positivity for ALDH1A3 and GFAP was observed. Furthermore, there was an above-average amount of microglial cells in this subtype. An association of ALDH1A3 to this subtype has been described before [26,36]. The enzyme’s influence on cell adhesion and tumor invasion and its capability to reduce oxidative stress could contribute to the poorer outcome of this subtype [26,37]. GFAP is also associated with a more aggressive tumor growth and it was mentioned as a characteristic of the mesenchymal subtype before [35]. The high amount of microglial cells in this subtype has a strong influence on the tumor’s microenvironment. Because their number exceed the amount of tumor cells in some tumors or tumor regions, their effect on therapy resistance must be considered. The observed strong negative correlation of EGFR and GFAP in our sample group is notable. Both of the markers have been mentioned as characteristics for miscellaneous subtypes, EGFR for the classical subtype [2] and GFAP for the mesenchymal subtype [35]. We hypothesize that the strong negative correlation could be due to different “strategical focuses” of the tumor cells with strong expression of EGFR in areas of cell proliferation and the strong expression of GFAP in areas of tumor invasion, as proposed before [38].
After defining a mesenchymal and classical subtype by cluster analysis, we could show that in some cases both subtypes can be found in the same tumor. Epigenetic examination confirmed these observations in one specimen. Additionally, survival analysis suggested that tumors that consist of both subtypes, defined as subtype-heterogeneous tumors, have a poorer outcome when compared to subtype-dominant tumors. This underlines the clinical importance of our study. The regional occurrence also has an impact on biopsy planning for therapeutic decision making, as a biopsy only reflects a small part of the tumor and it may not display the whole tumor’s morphology.
Furthermore, this study dealt with tumor development between primary and recurrent GBM. A chronological heterogeneity was also observed. The mesenchymal subtype dominated in the relapses, which suggests that cells of this subtype may have a higher therapy resistance and, hence, are responsible for relapse occurrence. Previous immunohistochemical studies used the term “mesenchymal transition” for this phenomenon [39,40]. However, as shown in our study, the mesenchymal subtype is also found in areas of the primary GBM, which leads to the preference of using the term “mesenchymal dominance” instead.
Even if our results give a clear picture of intratumoral and temporal heterogeneity in GBM, the results demand confirmatory studies with a larger of samples size, especially to substantiate the clinical impact. Broader examination of the epigenetic heterogeneity on a regional level is needed. The advantage of the use of immunohistochemistry is its broad availability and capacity to reflect the tumor’s morphology. This might clear the way to implementation of regional heterogeneity analyses into standard diagnostics. For this purpose, a standardized evaluation, e.g., by means of a defined heterogeneity index, should be developed.

4. Materials and Methods

4.1. Material

Tumor samples from 38 treatment-naïve IDH wildtype GBM patients (median age at diagnosis: 60 years; 29 male) plus material of 21 corresponding recurrent tumors were included in this retrospective study. According to Bayerisches Krankenhausgesetz, Artikel 27 it is allowed to use patient data for research given that they are anonymous. In the approval of the ethics committee (Bayerisches Krankenhausgesetz) for our study, it is stated that we do not need the patients’ consent for this retrospective study. Additionally, the material of two recurrent glioblastomas alone (median age at diagnosis: 60 years; two male) were analyzed. For patient data, also see Table 2.
The tissue samples, all formalin-fixed and paraffin-embedded (FFPE), derived from the Department of Neuropathology of the Institute of Pathology, Technical University Munich with a period of surgical resection at the Clinic and Polyclinic of Neurosurgery at Klinikum rechts der Isar from 2011 to 2017. Histopathological diagnosis was performed by neuropathologists and it was re-evaluated for this study according to WHO classification of tumors of the central nervous system, 2016 [41]. The study was performed according to the standards of the Helsinki Declaration of 1975 (as revised in 1983) and approved by the local ethics committee (reference number 164/19 S). Tumor samples of 38 patients were analyzed. Of 26 patients, the MGMT promotor status was known. The overall survival was only calculated for patients with non-methylated MGMT promotor.

4.2. Immunohistochemistry

For immunohistochemistry, 2 µm thick slides were cut with a standard microtome and dried at 76 °C for 30 min. EGFR, GFAP, Iba1, Olig2, p53, and Mib1 immunostaining was performed using a fully-automated staining system (Ventana BenchMark ULTRA; Ventana Medical Systems; Tucson, AZ, USA). In brief, the slides were exposed to heat-induced epitope uncovering in pH 8.4 buffer at 95 °C for 32 min. The tissue was incubated with H2O2 as an inhibitor of endogenous peroxidase in order to prevent unspecific bindings of the primary antibody. Afterwards, the slides were charged with anti-EGFR (monoclonal, mouse, dilution 1:50; Clone E30; DakoCytomation Denmark A/S, Glostrup, Denmark), anti-GFAP (monoclonal, mouse, dilution 1:100; Clone 6F2; DakoCytomation Denmark A/S, Denmark), anti-Iba1 (polyclonal, rabbit, dilution 1:500; Wako Pure Chemical Industries, Japan), anti-Olig2 (monoclonal, mouse, dilution 1:100; Clone 211F1.1; Cell Marque, Rocklin, CALIF, USA), anti-p53 (monoclonal, mouse, dilution 1:200; Clone DO-7; DakoCytomation Denmark A/S, Denmark), or anti-Mib1 (monoclonal, mouse, dilution 1:500; Clone MIB-1; VWR, Radnor, PA, USA) antibodies. For antibody detection, 3,3′-diaminobenzidine- (DAB-) based OptiView DAB IHC Detection Kit (Ventana Medical Systems) was used subsequently.
ALDH1A3 immunohistochemistry was manually performed, starting with epitope uncovering in pH 6.0 citrate buffer at 95 °C for 30 min., followed by H2O2 incubation. Anti-ALDH1A3 antibody (polyclonal, rabbit, dilution 1:600; Thermo Fisher Scientific, Waltham, MA, USA) was incubated overnight at 4 °C. Biotinylated secondary anti-rabbit IgG antibody (Vector Laboratories, Burlingame, CA, USA) in a dilution of 1:400 and subsequently, ABC-reagents (Vector Laboratories, Burlingame, CA, USA) were incubated for 30 min each., followed by DAB-reagent (Agilent, Santa Clara, CA USA).
For all immunostainings, counterstaining with haematoxylin was conducted and positive controls used as quality assurance.

4.3. Staining Evaluation

After slide digitalization with Aperio AT2 scanner (Leica biosystems, Wetzlar, Germany), evaluation of immunohistochemistry was performed by neuropathologists using Aperio ImageScope (version 12.3.0.5056, Leica biosystems, Germany). For every immunohistochemical marker, the areas of interest (AoI) were defined for each slide. For that, perceptible either particularly high or low protein expression was labeled for every marker. Afterwards, marked sectors were matched among all stainings of one tumor sample and overlapping zones were defined as AoI. Subsequently, scores for each marker in every AoI were given. Regarding Iba1, EGFR, GFAP, and Olig2, a variant of the established semiquantitative IRS by Remmele and Stegner was used, which consists of a product of scores of staining intensity (0 = no staining, 1 = weak, 2 = moderate, 3 = strong positivity) and percentage of positive cells (0 = 0%, 1 = 1–4%, 2 = 5–50%, 3 = 51–75%, 4 = 76–100%) [35,42]. For ALDH1A3, the score was slightly modified, since the maximum of stained cells only rarely exceeded 50 % (0 = 0%, 1 = 1–10%, 2 = 11–20%, 3 = 21–50%, 4 = 51–100%). For Mib1, only the percentage of positive cells was counted and translated into a score from 0–12 (0 = 0%, 3 = 1–19%, 6 = 20–29%, 9 = 30–39%, 12 = ≥ 40%). For p53, only strong positive cells were counted and scored while taking into account ranges that were suggested by Takami et al. [43] (0 = 0%, 1 = 1–10%, 2 = 11–89%, 3 = 90–100%).

4.4. Immunofluorescence Double Staining

Because the AoI mostly showed similar values for Iba1 and ALDH1A3, double staining with immunofluorescence was conducted in order to evaluate whether tumor cells of Iba1-positive microglial cells are ALDH1A3-positive. The pretreatment was conducted analogous to the above described ALDH1A3 immunohistochemistry. Afterwards, a mixture of the primary antibodies (anti-Iba1: monoclonal, mouse, dilution 1:100; Wako Pure Chemical Industries, Japan; anti-ALDH1A3: polyclonal, rabbit, dilution 1:800; Abcam, UK) was incubated over night at 4 °C followed by the secondary antibody (for Iba1: anti-mouse IgG Alexa Fluor 568, donkey, dilution 1:2000; Thermo Fischer Scientific, USA; for ALDH1A3: anti-mouse IgG Alexa Fluor 488, donkey, dilution 1:2000; Thermo Fischer Scientific, USA) for 45 min. Lastly, counterstaining with DAPI (Roche Diagnostics GmbH, Germany) was conducted, followed by covering with Aqua poly Maunt (Polysciences Inc., Warrington, PA, USA). Positive controls served as quality for evaluation, fluorescence microscope Axio Imager.Z2 (Zeiss, Oberkochen, Germany), and software AxioVision (version 4.8, Zeiss, Germany) were used.

4.5. 850k Methylation Array

DNA from marked tumor areas were extracted from FFPE material, followed by measuring the DNA concentration using Quibits dsDNA High sensitivity Assay kit (Invitrogen, Waltham, MA, USA) on a QuBit 4 system. DNA was then applied to Illumina EPIC BeadChip (Illumina, San Diego, CA, USA) for methylation analysis, as previously described [44].
For subsequent analysis of the 850k arrays, the BioConductor package “minfi” [45] was used. After initial quality control, standard processing (SWAN normalization, filtering probes with a detection p > 0.05) was applied. A set of glioma-specific CpG sites from Ceccarelli et al. [46] was selected for further analysis in order to visualize spatial proximity of samples. Principal Component Analysis of the filtered methylation data was performed, and samples were plotted by the two largest eigenvectors.

4.6. Statistical Analysis

Data were analyzed using R (version 3.3.2, R Foundation for Statistical Computing, Vienna, Austria) and IBM SPSS Statistics (version 22.0 IBM, Armonk, NY, USA). Marker correlation was tested by Spearman–Rho correlation analysis. Welch’s t-test and scatterplots served for comparison of primary and recurrent tumors. The connection of tissue size and number of AoI was analyzed by t-test. The Kaplan–Meier method was used for survival probability analysis. Hierarchical and PAM cluster analyses were conducted for the results of immunohistochemistry. Given the non-continuous nature of this data, Gower’s distance was chosen in order to calculate the underlying dissimilarity matrix (using R function “daisy”). Agglomerative hierarchical clustering with Ward’s method for clustering was performed using R function “agnes”. In parallel, PAM was performed for two groups (k = 2) while using R function “pam”. Statistical significance was defined as p < 0.05.

5. Conclusions

Our results indicate that the immunohistochemical detection of regional and temporal heterogeneity in GBM is feasible and it potentially provides information important for prognosis and therapy resistance. Consideration should be given to implementing immunohistochemical evaluation of tumor heterogeneity into standard neuropathological diagnostics.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-6694/12/10/2964/s1, Table S1: Results of immunohistochemistry.

Author Contributions

Conceptualization, F.L.-S., K.M., B.W. and J.S.; Data curation, F.S.-G., J.K., J.G., S.E.C., C.Z. and B.M.; Formal analysis, J.G., C.Z., B.M., B.V. and J.S.; Funding acquisition, B.W.; Investigation, F.L.-S., K.M., B.W. and J.S.; Methodology, F.L.-S., K.M., S.B., G.P., N.P. and W.W.; Project administration, F.L.-S.; Software, K.M., F.K., G.P. and B.W.; Supervision, B.W. and J.S.; Validation, F.L.-S., F.K., F.S.-G., W.W., J.G. and J.S.; Visualization, F.L.-S. and K.M.; Writing— original draft, F.L.-S. and K.M.; Writing— review & editing, S.B., F.S.-G., J.K., G.P., N.P., W.W., J.G., S.E.C., C.Z., B.M., B.W. and J.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the DFG, SFB-824, subproject B12 (grant to B.W.).

Acknowledgments

The authors would like to thank Ch. Schustetter and P. Meyer for their excellent technical assistance. This work was supported through the DFG, SFB-824, TPB12 (grant to B.W.)

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lee, J.-K.; Wang, J.; Sa, J.K.; Ladewig, E.; Lee, H.-O.; Lee, I.-H.; Kang, H.J.; Rosenbloom, D.S.; Camara, P.G.; Liu, Z.; et al. Spatiotemporal genomic architecture informs precision oncology in glioblastoma. Nat. Genet. 2017, 49, 594–599. [Google Scholar] [CrossRef] [PubMed]
  2. 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] [PubMed] [Green Version]
  3. Parker, N.R.; Hudson, A.L.; Khong, P.; Parkinson, J.F.; Dwight, T.; Ikin, R.J.; Zhu, Y.; Cheng, Z.J.; Vafaee, F.; Chen, J.; et al. Intratumoral heterogeneity identified at the epigenetic, genetic and transcriptional level in glioblastoma. Sci. Rep. 2016, 6, 22477. [Google Scholar] [CrossRef]
  4. Patel, A.P.; Tirosh, I.; Trombetta, J.J.; Shalek, A.K.; Gillespie, S.M.; Wakimoto, H.; Cahill, D.P.; Nahed, B.V.; Curry, W.T.; Martuza, R.L.; et al. Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma. Science 2014, 344, 1396–1401. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Darwin, C. The Origin of Species: By Means of Natural Selection or the Preservation of Favoured Races in the Struggle for Life; John Murray: London, UK, 1859. [Google Scholar]
  6. Nowell, P.C. The clonal evolution of tumor cell populations. Science 1976, 194, 23–28. [Google Scholar] [CrossRef] [PubMed]
  7. Merlo, L.M.; Pepper, J.W.; Reid, B.J.; Maley, C.C. Cancer as an evolutionary and ecological process. Nat. Rev. Cancer 2006, 6, 924–935. [Google Scholar] [CrossRef]
  8. Shackleton, M.; Quintana, E.; Fearon, E.R.; Morrison, S.J. Heterogeneity in cancer: Cancer stem cells versus clonal evolution. Cell 2009, 138, 822–829. [Google Scholar] [CrossRef] [Green Version]
  9. Bonavia, R.; Inda, M.D.M.; Cavenee, W.K.; Furnari, F.B. Heterogeneity maintenance in glioblastoma: A social network. Cancer Res. 2011, 71, 4055–4060. [Google Scholar] [CrossRef] [Green Version]
  10. Ho, I.A.W.; Shim, W.S.N. Contribution of the microenvironmental niche to glioblastoma heterogeneity. BioMed Res. Int. 2017, 2017, 1–13. [Google Scholar] [CrossRef]
  11. Wenger, A.; Vega, S.F.; Kling, T.; Bontell, T.O.; Jakola, A.S.; Carén, H. Intratumor DNA methylation heterogeneity in glioblastoma: Implications for DNA methylation-based classification. Neuro-Oncology 2019, 21, 616–627. [Google Scholar] [CrossRef]
  12. Dirkse, A.; Golebiewska, A.; Buder, T.; Nazarov, P.V.; Muller, A.; Poovathingal, S.; Brons, N.H.C.; Leite, S.; Sauvageot, N.; Sarkisjan, D.; et al. Stem cell-associated heterogeneity in Glioblastoma results from intrinsic tumor plasticity shaped by the microenvironment. Nat. Commun. 2019, 10, 1787. [Google Scholar] [CrossRef] [PubMed]
  13. Wang, Q.; Hu, B.; Hu, X.; Kim, H.; Squatrito, M.; Scarpace, L.; Decarvalho, A.C.; Lyu, S.; Li, P.; Li, Y.; et al. Tumor evolution of glioma-intrinsic gene expression subtypes associates with immunological changes in the microenvironment. Cancer Cell 2018, 33, 152. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Conroy, S.; Kruyt, F.A.E.; Joseph, J.V.; Balasubramaniyan, V.; Bhat, K.P.; Wagemakers, M.; Enting, R.H.; Walenkamp, A.M.E.; Dunnen, W.F.A.D. Subclassification of newly diagnosed glioblastomas through an immunohistochemical approach. PLoS ONE 2014, 9, e115687. [Google Scholar] [CrossRef] [PubMed]
  15. Popova, S.N.; Bergqvist, M.; Dimberg, A.; Edqvist, P.-H.; Ekman, S.; Hesselager, G.; Pontén, F.; Smits, A.; Sooman, L.; Alafuzoff, I. Subtyping of gliomas of various WHO grades by the application of immunohistochemistry. Histopathology 2013, 64, 365–379. [Google Scholar] [CrossRef] [Green Version]
  16. Bergmann, N.; Delbridge, C.; Gempt, J.; Feuchtinger, A.; Walch, A.; Schirmer, L.; Bunk, W.; Aschenbrenner, T.; Liesche, F.; Schlegel, J. The intratumoral heterogeneity reflects the intertumoral subtypes of glioblastoma multiforme: A regional immunohistochemistry analysis. Front. Oncol. 2020, 10, 494. [Google Scholar] [CrossRef]
  17. Gan, H.K.; Kaye, A.H.; Luwor, R. The EGFRvIII variant in glioblastoma multiforme. J. Clin. Neurosci. 2009, 16, 748–754. [Google Scholar] [CrossRef]
  18. Miconi, G.; Palumbo, P.; Dehcordi, S.R.; La Torre, C.; Lombardi, F.; Evtoski, Z.; Cimini, A.M.; Galzio, R.; Cifone, M.G.; Cinque, B. Immunophenotypic characterization of human glioblastoma stem cells: Correlation with clinical outcome. J. Cell. Biochem. 2015, 116, 864–876. [Google Scholar] [CrossRef]
  19. Parker, J.J.; Canoll, P.; Niswander, L.; Kleinschmidt-DeMasters, B.K.; Foshay, K.; Waziri, A. Intratumoral heterogeneity of endogenous tumor cell invasive behavior in human glioblastoma. Sci. Rep. 2018, 8, 18002. [Google Scholar] [CrossRef] [Green Version]
  20. Irvin, D.M.; McNeill, R.S.; Bash, R.E.; Miller, C.R. Intrinsic astrocyte heterogeneity influences tumor growth in glioma mouse models. Brain Pathol. 2016, 27, 36–50. [Google Scholar] [CrossRef] [Green Version]
  21. Trépant, A.L.; Bouchart, C.; Rorive, S.; Sauvage, S.; Decaestecker, C.; Demetter, P.; Salmon, I. Identification of OLIG2 as the most specific glioblastoma stem cell marker starting from comparative analysis of data from similar DNA chip microarray platforms. Tumor Biol. 2014, 36, 1943–1953. [Google Scholar] [CrossRef] [Green Version]
  22. Tsigelny, I.F.; Kouznetsova, V.L.; Lian, N.; Kesari, S. Molecular mechanisms of OLIG2 transcription factor in brain cancer. Oncotarget 2016, 7, 53074–53101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Kupp, R.; Shtayer, L.; Tien, A.C.; Szeto, E.; Sanai, N.; Rowitch, D.H.; Mehta, S. Lineage-restricted OLIG2-RTK signaling governs the molecular subtype of glioma stem-like cells. Cell Rep. 2016, 16, 2838–2845. [Google Scholar] [CrossRef] [Green Version]
  24. Viotti, J.; Duplan, E.; Caillava, C.; Condat, J.; Goiran, T.; Giordano, C.; Marie, Y.; Idbaih, A.; Delattre, J.-Y.; Honnorat, J.; et al. Glioma tumor grade correlates with parkin depletion in mutant p53-linked tumors and results from loss of function of p53 transcriptional activity. Oncogene 2013, 33, 1764–1775. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Duan, J.J.; Cai, J.; Guo, Y.F.; Bian, X.W.; Yu, S.C. ALDH1A3, a metabolic target for cancer diagnosis and therapy. Int. J. Cancer 2016, 139, 965–975. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Li, G.; Li, Y.; Liu, X.; Wang, Z.; Zhang, C.; Wu, F.; Jiang, H.; Zhang, W.; Bao, Z.; Wang, Y.; et al. ALDH1A3 induces mesenchymal differentiation and serves as a predictor for survival in glioblastoma. Cell Death Dis. 2018, 9, 1190. [Google Scholar] [CrossRef] [PubMed]
  27. Marchitti, S.A.; Brocker, C.; Stagos, D.; Vasiliou, V. Non-P450 aldehyde oxidizing enzymes: The aldehyde dehydrogenase superfamily. Expert Opin. Drug Metab. Toxicol. 2008, 4, 697–720. [Google Scholar] [CrossRef] [Green Version]
  28. Sarkar, S.; Döring, A.; Zemp, F.J.; Silva, C.; Lun, X.; Wang, X.; Kelly, J.; Hader, W.; Hamilton, M.; Mercier, P.; et al. Therapeutic activation of macrophages and microglia to suppress brain tumor-initiating cells. Nat. Neurosci. 2013, 17, 46–55. [Google Scholar] [CrossRef]
  29. Chen, Z.; Hambardzumyan, D. Immune microenvironment in glioblastoma subtypes. Front. Immunol. 2018, 9, 1004. [Google Scholar] [CrossRef] [Green Version]
  30. Le Mercier, M.; Hastir, D.; Lopez, X.M.; De Neve, N.; Maris, C.; Trépant, A.L.; Rorive, S.; Decaestecker, C.; Salmon, I. A simplified approach for the molecular classification of glioblastomas. PLoS ONE 2012, 7, e45475. [Google Scholar] [CrossRef] [Green Version]
  31. Nagy, A.; Garzuly, F.; Padányi, G.; Szűcs, I.; Feldmann, Á; Murnyák, B.; Hortobágyi, T.; Kálmán, B. Molecular subgroups of glioblastoma—An assessment by immunohistochemical markers. Pathol. Oncol. Res. 2017, 25, 21–31. [Google Scholar] [CrossRef]
  32. Keller, S.; Schmidt, M.H.H. EGFR and EGFRvIII promote angiogenesis and cell invasion in glioblastoma: Combination therapies for an effective treatment. Int. J. Mol. Sci. 2017, 18, 1295. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Furnari, F.B.; Cloughesy, T.F.; Cavenee, W.K.; Mischel, P.S. Heterogeneity of epidermal growth factor receptor signalling networks in glioblastoma. Nat. Rev. Cancer 2015, 15, 302–310. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Lu, F.; Chen, Y.; Zhao, C.; Wang, H.; He, D.; Xu, L.; Wang, J.; He, X.; Deng, Y.; Lu, E.E.; et al. Olig2-dependent reciprocal shift in PDGF and EGF receptor signaling regulates tumor phenotype and mitotic growth in malignant glioma. Cancer Cell 2016, 29, 669–683. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Motomura, K.; Natsume, A.; Watanabe, R.; Ito, I.; Kato, Y.; Momota, H.; Nishikawa, R.; Mishima, K.; Nakasu, Y.; Abe, T.; et al. Immunohistochemical analysis-based proteomic subclassification of newly diagnosed glioblastomas. Cancer Sci. 2012, 103, 1871–1879. [Google Scholar] [CrossRef] [PubMed]
  36. Zhang, W.; Liu, Y.; Hu, H.; Huang, H.; Bao, Z.; Yang, P.; Wang, Y.; You, G.; Yan, W.; Jiang, T.; et al. ALDH1A3: A marker of mesenchymal phenotype in gliomas associated with cell invasion. PLoS ONE 2015, 10, e0142856. [Google Scholar] [CrossRef] [PubMed]
  37. Chang, P.M.-H.; Chen, C.-H.; Yeh, C.-C.; Lu, H.-J.; Liu, T.-T.; Chen, M.-H.; Liu, C.-Y.; Wu, A.T.H.; Yang, M.-H.; Tai, S.-K.; et al. Transcriptome analysis and prognosis of ALDH isoforms in human cancer. Sci. Rep. 2018, 8, 2713. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Herpers, M.J.H.M.; Budka, H.; McCormick, D. Production of glial fibrillary acidic protein (GFAP) by neoplastic cells: Adaptation to the microenvironment. Acta Neuropathol. 1984, 64, 333–338. [Google Scholar] [CrossRef]
  39. Schäfer, N.; Gielen, G.H.; Rauschenbach, L.; Kebir, S.; Till, A.; Reinartz, R.; Simon, M.; Niehusmann, P.; Kleinschnitz, C.; Herrlinger, U.; et al. Longitudinal heterogeneity in glioblastoma: Moving targets in recurrent versus primary tumors. J. Transl. Med. 2019, 17, 96. [Google Scholar] [CrossRef] [Green Version]
  40. Wood, M.D.; Reis, G.F.; Reuss, D.E.; Phillips, J.J. Protein analysis of glioblastoma primary and posttreatment pairs suggests a mesenchymal shift at recurrence. J. Neuropathol. Exp. Neurol. 2016, 75, 925–935. [Google Scholar] [CrossRef]
  41. Louis, D.N.; Perry, A.; Reifenberger, G.; Von Deimling, A.; Figarella-Branger, D.; Cavenee, W.K.; Ohgaki, H.; Wiestler, O.D.; Kleihues, P.; Ellison, D.W. The 2016 World Health Organization classification of tumors of the central nervous system: A summary. Acta Neuropathol. 2016, 131, 803–820. [Google Scholar] [CrossRef] [Green Version]
  42. Remmele, W.; Stegner, H.E. Recommendation for uniform definition of an immunoreactive score (IRS) for immunohistochemical estrogen receptor detection (ER-ICA) in breast cancer tissue. Der Pathol. 1987, 8, 138–140. [Google Scholar]
  43. Takami, H.; Yoshida, A.; Fukushima, S.; Arita, H.; Matsushita, Y.; Nakamura, T.; Ohno, M.; Miyakita, Y.; Shibui, S.; Narita, Y.; et al. Revisiting TP53 mutations and immunohistochemistry—A comparative study in 157 diffuse gliomas. Brain Pathol. 2014, 25, 256–265. [Google Scholar] [CrossRef] [PubMed]
  44. Reinhardt, A.; Stichel, D.; Schrimpf, D.; Sahm, F.; Korshunov, A.; Reuss, D.E.; Koelsche, C.; Huang, K.; Wefers, A.K.; Hovestadt, V.; et al. Anaplastic astrocytoma with piloid features, a novel molecular class of IDH wildtype glioma with recurrent MAPK pathway, CDKN2A/B and ATRX alterations. Acta Neuropathol. 2018, 136, 273–291. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Aryee, M.J.; Jaffe, A.E.; Corrada-Bravo, H.; Ladd-Acosta, C.; Feinberg, A.P.; Hansen, K.D.; Irizarry, R.A. Minfi: A flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics 2014, 30, 1363–1369. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Ceccarelli, M.; Barthel, F.P.; Malta, T.M.; Sabedot, T.S.; Salama, S.R.; Murray, B.A.; Morozova, O.; Newton, Y.; Radenbaugh, A.; Pagnotta, S.M.; et al. Molecular profiling reveals biologically discrete subsets and pathways of progression in diffuse glioma. Cell 2016, 164, 550–563. [Google Scholar] [CrossRef] [PubMed]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure 1. Iba1 and ALDH1A3 immunofluorescence double staining. The picture shows that most cells do not exhibit a co-expression of Iba1 and ALDH1A3. Green arrows mark Iba1-positive microglial cells that lack red signal for ALDH1A3. Red arrows point to ALDH1A3-positive cells without green signal for Iba1. Scale bars: 10 µm.
Figure 1. Iba1 and ALDH1A3 immunofluorescence double staining. The picture shows that most cells do not exhibit a co-expression of Iba1 and ALDH1A3. Green arrows mark Iba1-positive microglial cells that lack red signal for ALDH1A3. Red arrows point to ALDH1A3-positive cells without green signal for Iba1. Scale bars: 10 µm.
Cancers 12 02964 g001
Figure 2. Hierarchical clustering divides the tumor areas into two subtypes. The blue frame is used for cluster A, which is comparable to the classical subtype by Verhaak et al. [2], the red frame marks cluster B with characteristics of the mesenchymal subtype. In the lower section of the figure, the characteristic immunohistochemical profiles of both clusters are pictured as bar charts. Red bars are used for markers that exceed the average immune reactive score (IRS). The stripes directly below the hierarchical tree with the same color represent the different AoI of the same tumor. In these stripes, information about O-6-methylguanin-DNA methyltransferase (MGMT) promotor status can be found (M = methylated, NM = non-methylated). The dots under the different colored AoIs samples, which consist of pairs of three colors, mark the AoIs that were analyzed epigenetically (the same colors were used in Figure 5b). The red crosses mark the AoIs, whose cluster assignment differs in partitioning around medoids (PAM).
Figure 2. Hierarchical clustering divides the tumor areas into two subtypes. The blue frame is used for cluster A, which is comparable to the classical subtype by Verhaak et al. [2], the red frame marks cluster B with characteristics of the mesenchymal subtype. In the lower section of the figure, the characteristic immunohistochemical profiles of both clusters are pictured as bar charts. Red bars are used for markers that exceed the average immune reactive score (IRS). The stripes directly below the hierarchical tree with the same color represent the different AoI of the same tumor. In these stripes, information about O-6-methylguanin-DNA methyltransferase (MGMT) promotor status can be found (M = methylated, NM = non-methylated). The dots under the different colored AoIs samples, which consist of pairs of three colors, mark the AoIs that were analyzed epigenetically (the same colors were used in Figure 5b). The red crosses mark the AoIs, whose cluster assignment differs in partitioning around medoids (PAM).
Cancers 12 02964 g002
Figure 3. PAM divides the tumor areas into two subtypes. Similar to Figure 2, the blue frame is used for the classical subtype, the red for the mesenchymal subtype. The bar charts show the immunohistochemical expression profiles of the medoids, the samples with the most characteristic immunohistochemical profile of each cluster. Red bars are used for markers that exceed the average (immune reactive score) IRS.
Figure 3. PAM divides the tumor areas into two subtypes. Similar to Figure 2, the blue frame is used for the classical subtype, the red for the mesenchymal subtype. The bar charts show the immunohistochemical expression profiles of the medoids, the samples with the most characteristic immunohistochemical profile of each cluster. Red bars are used for markers that exceed the average (immune reactive score) IRS.
Cancers 12 02964 g003
Figure 4. Characteristic immunohistochemical profile of both subtypes. The blue frame marks the classical subtype, the red the mesenchymal subtype. For each, the average IRS of all markers are pictured in the bar charts. Basis is the cluster assignment of PAM. Red bars are used for markers that exceed the average IRS. The light grey bars depict the average IRS for all AoIs.
Figure 4. Characteristic immunohistochemical profile of both subtypes. The blue frame marks the classical subtype, the red the mesenchymal subtype. For each, the average IRS of all markers are pictured in the bar charts. Basis is the cluster assignment of PAM. Red bars are used for markers that exceed the average IRS. The light grey bars depict the average IRS for all AoIs.
Cancers 12 02964 g004
Figure 5. Survival analysis and principal component analysis (PCA) of epigenetic examination. (a): By Kaplan–Meier method, a shorter survival of patients with ST-dom tumors compared to patients with ST-het tumors can be observed (p = 0.166). (b): The methylation profile of each two different AoI of three tumors was examined. PCA revealed a proximity of the both samples of tumors Het01 and Het03. Interestingly, the methylation profiling confirms the allocation of the samples of tumor Het02 to the two different clusters by PAM, as one AoI shows a distinct nearer proximity to the samples of Het01 as to the other AoI of the same tumor sample.
Figure 5. Survival analysis and principal component analysis (PCA) of epigenetic examination. (a): By Kaplan–Meier method, a shorter survival of patients with ST-dom tumors compared to patients with ST-het tumors can be observed (p = 0.166). (b): The methylation profile of each two different AoI of three tumors was examined. PCA revealed a proximity of the both samples of tumors Het01 and Het03. Interestingly, the methylation profiling confirms the allocation of the samples of tumor Het02 to the two different clusters by PAM, as one AoI shows a distinct nearer proximity to the samples of Het01 as to the other AoI of the same tumor sample.
Cancers 12 02964 g005
Figure 6. Comparison of marker expression in primary and recurrent GBM. The charts show a strong increase of ALDH1A3 (a) and Iba1 (c) positivity between primary and recurrent tumors. In contrast, EGFR (b) and proliferation (Mib1; d) decrease.
Figure 6. Comparison of marker expression in primary and recurrent GBM. The charts show a strong increase of ALDH1A3 (a) and Iba1 (c) positivity between primary and recurrent tumors. In contrast, EGFR (b) and proliferation (Mib1; d) decrease.
Cancers 12 02964 g006
Table 1. Spearman-Rho correlation analysis.
Table 1. Spearman-Rho correlation analysis.
ALDH1A3GFAPIba1EGFRp53Olig2Mib1
ALDH1A3rs 0.228 *0.468 **−0.0640.111−0.048−0.121
p 0.0140.0000.4950.2400.6120.199
GFAPrs0.228 * 0.085−0.291 **−0.115−0.345 **−0.413 **
p0.014 0.3690.0020.2200.0000.000
Iba1rs0.468 **0.085 −0.0900.176−0.1310.018
p0.0000.369 0.3380.0590.1620.848
EGFRrs−0.064−0.291 **−0.090 −0.1050.326 **0.282 **
p0.4950.0020.338 0.2630.0000.002
p53rs0.111−0.1150.176−0.105 0.0260.298 **
p0.2400.2200.0590.263 0.7830.001
Olig2rs−0.048−0.345 **−0.1310.326 **0.026 0.415 **
p0.6120.0000.1620.0000.783 0.000
Mib1rs−0.121−0.413 **0.0180.282 **0.298 ** 0.415 **
p0.1990.0000.8480.0020.0010.000
One star indicates significant level on a basis of p < 0.05, two stars on a basis of p < 0.01.
Table 2. Patient data.
Table 2. Patient data.
ParameterNumber of PatientsData
Age (in years)n = 38mean59
median60
range27–84
Sexn = 38male29
female9
MGMT promotor statusn = 26methylated6
non-methylated20
Overall survival (in months)n = 20mean (95% CI)22 (17–26)
median (95% CI)19 (13–25)

Share and Cite

MDPI and ACS Style

Liesche-Starnecker, F.; Mayer, K.; Kofler, F.; Baur, S.; Schmidt-Graf, F.; Kempter, J.; Prokop, G.; Pfarr, N.; Wei, W.; Gempt, J.; et al. Immunohistochemically Characterized Intratumoral Heterogeneity Is a Prognostic Marker in Human Glioblastoma. Cancers 2020, 12, 2964. https://doi.org/10.3390/cancers12102964

AMA Style

Liesche-Starnecker F, Mayer K, Kofler F, Baur S, Schmidt-Graf F, Kempter J, Prokop G, Pfarr N, Wei W, Gempt J, et al. Immunohistochemically Characterized Intratumoral Heterogeneity Is a Prognostic Marker in Human Glioblastoma. Cancers. 2020; 12(10):2964. https://doi.org/10.3390/cancers12102964

Chicago/Turabian Style

Liesche-Starnecker, Friederike, Karoline Mayer, Florian Kofler, Sandra Baur, Friederike Schmidt-Graf, Johanna Kempter, Georg Prokop, Nicole Pfarr, Wu Wei, Jens Gempt, and et al. 2020. "Immunohistochemically Characterized Intratumoral Heterogeneity Is a Prognostic Marker in Human Glioblastoma" Cancers 12, no. 10: 2964. https://doi.org/10.3390/cancers12102964

APA Style

Liesche-Starnecker, F., Mayer, K., Kofler, F., Baur, S., Schmidt-Graf, F., Kempter, J., Prokop, G., Pfarr, N., Wei, W., Gempt, J., Combs, S. E., Zimmer, C., Meyer, B., Wiestler, B., & Schlegel, J. (2020). Immunohistochemically Characterized Intratumoral Heterogeneity Is a Prognostic Marker in Human Glioblastoma. Cancers, 12(10), 2964. https://doi.org/10.3390/cancers12102964

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