Introduction

Patients diagnosed with Glioblastoma (GBM, WHO-IV) have extremely poor median survival times, despite modern microsurgery, chemoradiotherapy, reoperation and experimental therapies [14]. To improve GBM treatment and patient median survival times, fresh insight into the molecular pathogenesis of GBM is essential. Proteomics can define molecular pathways and cellular functions altered in GBM [5]. Genomic studies, although important, are limited by the fact that normal, upregulated or mutated genes may not be transcribed for a number of epigenetic reasons [6]. Multiple discrepancies between mRNA and proteomic expression profiles in differential analyses of gliomas highlight the importance of studying protein expression [7].

Age is a powerful individual prognostic indicator [810]. Long term survivors of GBM are invariably younger patients [810], and in one randomised clinical trial, median survival for GBM cohorts aged <45 years was 48 weeks compared to 19 weeks for those >65 years; and at 18 months 23 % of the younger cohort was alive compared to 3 % of the older cohort [8]. Numerous randomised controlled trials and hospital series have excluded differences in access to health care as the cause for this differential outcome [1, 912].

The biological basis of this powerful age-related effect is not understood. The histological features of GBM, cellular proliferative indices, epidermal growth factor receptor amplification and p53 expression are very similar irrespective of age [9]. Although younger patients are more likely to have secondary GBM than the elderly [3, 9, 13], there is no difference in outcomes between primary and secondary GBM once the diagnosis is made and patients have been aged matched with controls [13]. Although genetic differences associated with short survival (6q loss, 10q loss, 19q gain, sodium ion channel mutations) and longer survival (TP53 mutations and the combination of LOH1p and LOH19q, MGMT status, mutation of IDH1) have been identified, these differences have not, with the exception of IDH1 mutations where the mutation occurred in younger patients, been analysed with respect to patient age [1420].

Analyses of GBM samples from older patients has begun to provide a coherent view of the proteomic response in GBM but interpretation is complicated by differences in experimental design and proteomic technology [5]. In this study we provide the first systematic proteomic analysis in young GBM (versus age-matched peritumoural-control brain) to gain insight into the basis of the importance of age on prognosis. For the purpose of comparison, we performed a parallel, contemporaneous study (using the same experimental design and technology) in old GBM.

Materials and methods

Clinical material

Glioblastoma and peritumoural-control brain samples were obtained from young (<45 years) and old patients (>60 years) undergoing resective brain tumour surgery (Ethical approval: LREC/2004/4/16). The sampling procedure and clinical details of experimental samples are described in Supplementary Tables 1 and 2 and Supplementary methods.

The experimental group sizes used for the primary proteomic analysis were young GBM (n = 7) and young peritumoural-control (n = 12) (based on a priori power calculations to detect significant changes of ≥35 % with power ≥0.8). Tissue was collected for two comparison groups: old GBM (n = 13) and old peritumoural-control (n = 10). The median co-efficients of variation were similar in each experimental group: young GBM 33.55 %, young peritumoural-control 27.38 %, old GBM 33.99 % and old peritumoural-control 26.39 % (Supplementary Fig. 1).

Total protein extracts were separated by isoelectric point and molecular mass using 2DGE (see Supplementary-methods). 2D-gel images were captured using a FluorChem Image Analyser and aligned in a single study using at least four manual alignment vectors followed by automatic placement of further alignment vectors by the software (~200–400 vectors per gel). The mean protein levels of each protein were analysed using Student’s t test (p ≤ 0.003, equivalent to p ≤ 0.01 with Bonferroni correction factor 3 for each comparison). Significant data are presented in Table 1 (uncorrected for multiple comparisons) and all data are presented in Supplementary Table 3. Protein spots differentially expressed in young GBM versus young peritumoural-control were manually excised and proteins identified using LC–MS [21]. LC–MS runs of each sample were combined using Maxquant, assuming a false positive rate of 0.01 [22]. An identical approach was applied in parallel comparing old GBM versus old peritumoural-control, and young GBM versus old GBM.

Table 1 Proteins altered in young GBM

Immunoblot analysis was performed on a subset of the tissue samples used for proteomics. Proteins (10 μg) were separated by SDS-PAGE and transferred to nitrocellulose membrane. Primary antibodies were detected using fluorescently-labelled secondary antibodies and were visualized using an Odyssey Imager.

To assess functional protein–protein interactions between the proteins altered in young GBM (p ≤ 0.01), altered protein identifiers were uploaded to Ingenuity Pathway Analysis (IPA; http://www.ingenuity.com). Networks were algorithmically generated based on direct relationships (physical interactions and/or associations) between eligible proteins. Networks are scored and ranked according to the inclusion of as many proteins inputted as possible. Network scores are putatively a measure of probability [23]. For comparison purposes, network analysis was also performed on proteins differentially regulated in old GBM versus old peritumoural-control. Network analysis is a powerful tool for identifying potential interactions between altered proteins (hypothesis generation) which can be subsequently explored in functional analyses.

Results

The median survival in the young GBM cohort was >39 months (3 of the 7 patients are still alive) and this was significantly greater than the median survival of 9 months in the older reference group (p < 0.02). Performance status (see Supplementary Table 1) was not significantly different between the young (mean age 36 years) and old (mean age 67 years) GBM cohorts.

Young GBM: proteins differentially expressed in young GBM compared to age matched controls

A total of 405 protein spots were matched across every 2Dgel (young GBM and young peritumoural-control gels) and analysed. Logarithmic association of the 405 protein expression levels (mean normalised volumes) highlights multiple protein alterations in young GBM (Fig. 1a, Supplementary Fig. 2). 90 protein spots were altered in young GBM (versus young peritumoural-control; p ≤ 0.01) and the identity of these 90 statistically significant altered spots was established by LCMS.

Fig. 1
figure 1

a Overview of proteomic analysis of young GBM. Over 400 spots were identified by 2D gel electrophoresis. The normalised volume represents the relative amount of protein in the spot. Each point in the graph represents the relative amount of protein in the 400 spots analysed. Table 1 lists the proteins which are significantly altered in young GBM. In contrast this graph emphasises that the levels of the majority (more than 75 %) of proteins are unaltered in young GBM. Because of the dynamic range (300 fold difference from the most abundant to the least abundant protein), data are presented as logs. There is a good correlation between young GBM and young peritumoural control (r2 = 0.85), with 22 % of the spots significantly altered (see Supplementary Table 1). b Overview of proteomic response in young GBM compared to old GBM. Over 400 spots were identified by 2D gel electrophoresis. The normalised volume represents the relative amount of protein in the spot. Each point in the graph represents the relative amount of protein in the 400 spots analysed. There is an excellent correlation between young GBM and old GBM (r2 = 0.95) with only 1 % of the spots significantly altered (5 out of 405; see text for details)

Sixty eight unique proteins were significantly altered in young GBM (Table 1, Supplementary Table 3). 15 of these proteins were identified multiple times in 2–5 spots (ATP6V1B2, OXCT1, ALDOA, GFAP, DCD, DPYSL2, TUBB2A, INA, MBP, ACOT7, VDAC2, UCHL1, PGAM1, PRDX3 and GNB1). Identification of the same protein in several spots is a feature of 2DGE proteomic studies and explains the difference between the number of altered protein spots and number of unique proteins identified. From the 68 altered proteins identified, 29 proteins were up-regulated and 39 proteins were down-regulated. A major fraction of the proteins altered in young GBM (25 %; 16 out of the 68 proteins) are putatively localised to mitochondria (OXCT1, PEBP1, DPYSL2, CKMT1A, ACOT7, CKB, IDH3A, SNAP, VDAC2, PRDX3, HSPD1, CAT, ATP6V1E1, GLUD1, CLIC4 and NDUFS3). 12 of the 68 proteins altered in young GBM have previously been described altered in proteomic studies of glioma (APOA1, GFAP, HSPA5, PDIA3, TUBB2A, GLUD1, GSTP1, PGAM1, UCHL1, HSPB1, HSPD1 and SRI) [5]. Notably, over 50 proteins have been described altered in GBM for the first time.

Ten proteins (DPYSL2, SRI, OXCT1, UCHL1, CAT, SEPT11, IDH3A, PDIA3, ATP6V1B2, PRDX3), altered in young GBM were examined using western blotting. Western blotting of young GBM versus young peritumoural-control tissue, demonstrated that 7 out of the 10 proteins tested were significantly altered (p ≤ 0.01) and that 10 out of the 10 proteins showed the same direction of response as the proteomic analysis (Fig. 2, Supplementary Fig. 3).

Fig. 2
figure 2

Proteomic alterations in young GBM: confirmation with western analysis. Western blotting replicates the alterations in defined proteins in GBM in a subset (determined by tissue availability) from the same subjects as used in the proteomic 2D gel electrophoresis. a 2D gel electrophoresis identified a consistent significant (p = 1.8E−06) reduction in DPYSL2 in young GBM. b Western blot analysis identified a similar consistent reduction in DPYSL2 in young GBM. c 2D gel electrophoresis identified a significant increase (p = 0.0057) in Sorcin (though with inter-subject variability) in young GBM. d Western blot analysis identified a similar increase in Sorcin in young GBM, again with inter-subject variability. There was also good correspondence between 2D gel electrophoresis and western blot analysis in young GBM for all 10 proteins examined with both techniques (see Supplementary Fig. 3)

IPA network analysis was performed on the proteomic dataset referred to as young GBM and included 68 proteins. The young GBM dataset generated multiple functional protein networks (Table 2) including 4 high scoring networks containing 23, 15, 12 and 11 dataset proteins respectively (Table 2).

Table 2 Putative interactions between proteins altered in young GBM and old GBM

The top network generated by IPA (Fig. 3) included multiple structural proteins downregulated in young GBM, for example Strathmin (STMN1) and dihydropyrimidinase-related protein 2 (DPYSL2). The network also contained GFAP, upregulated in young GBM, which has long been considered a fundamental and diagnostic feature of glioma [23]. The network included heat shock proteins (HSPD1, HSPA8, HSPB1), and a group of downregulated proteins involved in ATP homeostasis and energy metabolism (ALDOA, ATP6V1E1, CKB, CKMT1A), consistent with existing evidence but also identifying for the first time specific protein networks that may be involved in the dysregulation of energy metabolism in malignant glioma [24]. Lastly a cluster of upregulated proteins, integral to the immunoproteosome (PSME1, PSME2, 20 s/26 s proteosome, PSMB7), was highlighted in the top network. Network 2 was characterised by a cluster of Septin proteins (GTPase proteins that have been shown to play a role in gliomagenesis [25], and the insertion of a hub protein, TRAF6, a signal transducer in NFkappaB signalling. Network 3 was characterised by the insertion of a hub protein HNF4alpha, a transcription factor recently shown to play a role in other neoplasias [26] and Network 4 was characterised by numerous mitochondrial-localised proteins (CAT, IDH3A, NDUFS3 and other complex 1 proteins, OXCT1 and PRDX3).

Fig. 3
figure 3

Protein-protein interactions in young GBM. Visual representation of putative protein–protein interactions in Network 1 (the highest scoring network; Table 2) generated by IPA in young GBM. Each node (shape) represents a protein and its association with other proteins, is represented by a line. Nodes have different shapes that represent different molecule types, for example, transcription factors, enzymes, kinases and phosphatases (refer to Ingenuity Systems Software for detailed node information). Proteins or ‘nodes’ with a coloured background were regulated in the study (green = downregulated; red = upregulated) whilst other interacting proteins with no background are proteins not detected in this study that have been inserted by IPA to produce a highly connected network. The solid lines represent direct interactions or associations between proteins

Old GBM: proteins differentially expressed in Old GBM compared to age matched control tissue

To allow the extensive protein alterations in young GBM to be compared with those in old GBM, a proteomic evaluation was conducted contemporaneously in old GBM (patients >60 years) using the same technology. A total of 405 protein spots were matched across every 2D gel (old GBM and old peritumoural-control gels) and analysed. Logarithmic association of the 405 protein expression levels (mean normalised volumes) was broadly similar to that seen in young GBM (Supplementary Fig. 2). 70 protein spots were altered in old GBM versus old peritumoural-control (p ≤ 0.01).

55 unique proteins were altered significantly in old GBM (listed in Table 1 and Supplementary Table 3 in full). 8 of these proteins were identified multiple times in 2-4 spots (GNB1, INA, ALDOA, SEPT11, HSPB1, CKB, CLIC and GLUL). From the 55 altered proteins identified, 16 proteins were up-regulated and 39 proteins were down-regulated. 19 of the 55 proteins have been reported to be putatively localised to mitochondria (OXCT1, ATP6V1E1, NAPG, NDUFS3, ACOT7, CKB, CKMT1A, DPYSL2, GLUL, HSPD1, IDH3A, PRDX1, PRDX3, PRDX5, PEBP1, SUCLA2, SOD1, UQCRFSL1 and VDAC2; 13 of these 19 proteins were also altered in young GBM). 9 of the 55 proteins have previously been described altered in proteomic studies of glioma (APOA1, GFAP, HSPA5, PDIA3, GSTP1, UCHL1, HSPB1, HSPD1 and PRDX1; 8 of these 9 proteins were also altered in young GBM). Over 40 proteins have been described altered in GBM for the first time (24 of which were also altered in young GBM).

As in young GBM, protein alterations were confirmed by Western blot analysis: 8 out of 10 proteins tested were significantly altered (p ≤ 0.01) and 10 out of the 10 proteins showed the same direction of response as the proteomic analysis (Supplementary Fig. 4).

The highest scoring networks for old GBM highlighted considerable commonality to the proteomic response in young GBM (Table 2).

Proteomic response in young GBM: a comparison with the proteomic response in old GBM

Five protein spots were differentially altered between young and old GBM (p ≤ 0.01). Only three unique statistically altered proteins were identified (PEBP1, NDUFA10 and PGK1). The proteins in two spots were not identified. This number of altered proteins lies beneath the multiple testing threshold of potential false positive results in the study. There was an excellent correlation (r2 = 0.95) between the level of 405 proteins analysed in GBM from younger patients and their level in GBM from older patients (Fig. 1b). This correlation was similar to that seen in peritumoural-control samples from the two age groups (Supplementary Fig. 2). There were good correlations (r2 = 0.85 and 0.90) between the level of 405 proteins in young GBM relative to young peritumoural-controls (Fig. 1a) and old GBM relative to old peritumoural-control (Supplementary Fig. 2) respectively. 48 unique proteins altered in GBM were common in young and old cohorts. The direction and fold change of all 48 proteins was consistent in both young and old GBM. From the top 25 altered protein spots identified in Young GBM (ranked by p value), 17 were demonstrated altered in old GBM (p ≤ 0.01).

Discussion

The present study provides a powerful example of how proteomics can reliably test a hypothesis (i.e. is the most important prognostic variable in GBM, age, associated with a distinct response) and demonstrates that proteomics can play an important role in understanding GBM pathophysiology. Definition of the proteomic response in samples from patients with a homogeneous and clinically defined age range (18–45 years) addresses one of the design weaknesses in proteomic studies of GBM to date [5].

Young (<45 years of age) and old (>60 years of age) GBM cohorts with a mean age difference of 31 years and a significantly better median survival despite optimal therapy in the younger cohort were recruited. Multiple protein alterations were detected in young and old GBM versus age matched control tissue, and included a mixture of previously well-characterised protein alterations in GBM (for example, GFAP and UCHL1), and the identification of many ‘highly expected’ heat shock proteins (HSPD1, HSPB1, HSPA5, HSPA8) and cytoskeletal proteins (TUBB2A, TUBB2C), which confirm the robustness of our proteomic data.

One cluster of upregulated proteins (Fig. 3) in both young and old GBM comprised PSME1, PSME2, 20 s/26 s proteosome and PSMB7. These interacting proteins are central to the immunoproteosome (i-proteosome). All proteasomes contain a 20 s subunit flanked by either 19 s subunits or 11 s subunits. In the standard proteasome two 19 s subunits enclose a 20 s subunit of 2α rings sandwiching 2β rings with proteolytic subunits (β1, β2, β5). In the i-proteasome these catalytic subunits are substituted by LMP2, MECL, LMP7 and the 20 s is flanked by two 11 s/PA28 subunits. The 11 s contains 3α & 3β alternating subunits regulated by PSME1 and PSME2 respectively [27], two of the proteins upregulated in our young and old GBM analyses. Inhibition of the 20 s/26 s proteasome by drugs, such as carfilzomib leads to a build up of poly-ubiquinated proteins causing cell cycle arrest, apoptosis and inhibition of tumour growth [28].

i-proteasome function is to provide peptides for MHC-class1-antigen presentation. Interferon increases i-proteasome numbers during inflammation and oxidative damage [29]. 26 s proteasomes are ineffective at degrading oxidised proteins, but i-proteasomes can efficiently process these damaged proteins [27]. Increased PSME1 and PSME2 as a result of interferon, would prevent protein build up and apoptosis. Conversely loss of i-proteasome function, through inhibition of 11 s subunit formation or joining of the 11 s subunit to the 20 s, would have the two fold effect of damaged protein aggregation, leading to apoptosis; and the removal of ‘self’ peptides from the cell surface, alerting the immune system to the malignant tumour cells. Elucidating the mechanisms of GBM immune resistance and causes of immunosuppression is currently an area of intense research and therapeutic effort in GBM [30, 31].

The proteomic analyses of young and old GBM also highlighted multiple proteins (PRDX3, UCHL1, PEBP1, DPYSL2, UBE2 N, GSTO) involved in nuclear factor kappaB (NFkB) regulation. NFkB is a transcription factor capable of mediating many cellular responses and modulates oncogenesis, tumour progression and chemotherapy resistance [3235]. In the cytoplasm, NFkB is a small protein complex containing two subunits that bind inhibitory kappa B (IkB). IkB binding prevents NFkB translocation to the nucleus. Activation of NFkB, with subsequent translocation to the nucleus can occur through the canonical (utilising IkB kinase, IKK), the non-canonical or the alternative pathway. In the nucleus NFkB regulates transcription of proteins that down-regulate apoptosis, increase cell invasiveness, increase angiogenesis and increase vascular permeability, thereby promoting tumourgenesis [32, 36]. Proteins that regulate NFkB function were altered in GBM. GSTO1 (upregulated in the young GBM analysis) increases IL1β levels which activates IKK. GSTO1 also increases Akt phosphorylation in cells exposed to the pro-apoptotic drug cisplantin. Phosphorylated Akt inhibits apoptosis via NFkB [37]. UBE2 N (also upregulated in the young GBM analysis) is also vital for the activation of IKK via TRAF6 [38]. TRAF6, a core signal transducer in the NFkB pathway was highlighted as a hub protein in IPA Network 2 of both young and old GBM IPA analyses. PRDX3 (downregulated in young and old GBM analyses) also increases IKK activation [39], and knock-down studies of UCHL1 (also downregulated in young and old GBM analyses), show an increase in NFkB function via IKK activation [40]. PEBP1 (upregulated in young GBM and downregulated in old GBM) antagonises NFkB function by interfering with the TNFα pathway [41], resulting in an increase in NFkB function. NFkB’s role in gliomagenesis is summarised in Supplementary Fig. 5. NFkB inhibitors have shown promise in inducing cell death in GBM [42].

Alignment of protein alterations identified in young and old GBM versus age-matched peritumoural-controls showed considerable commonality in the proteomic response of GBM in different aged patients (and also demonstrated the rigour of our two distinct proteomic analyses of GBM). Our study does not provide a clear explanation as to why young and old patients with GBM have differential prognoses. One of the few proteins putatively altered in expression level between young and old GBM, is Phosphatidyl ethanolamine binding protein 1 (PEBP1; also known as Raf1-kinase inhibitor protein, RKIP). PEBP1 was found significantly upregulated in our young GBM proteomic analysis, significantly downregulated in our old GBM analysis, and significantly downregulated in old GBM compared to young GBM. PEBP1 inhibits the RAF/MEK/ERK pro-oncogenic pathway and also inhibits NFkB (also pro-oncogenic) by antagonising the activity of IKK either directly or via Tumour Necrosis Factor alpha (TNFalpha) [41]. The difference in PEBP1 expression levels between young and old GBM could contribute to their different prognosis.