Differential analysis of glioblastoma multiforme proteome by a 2D-DIGE approach

Background Genomics, transcriptomics and proteomics of glioblastoma multiforme (GBM) have recently emerged as possible tools to discover therapeutic targets and biomarkers for new therapies including immunotherapy. It is well known that macroscopically complete surgical excision, radiotherapy and chemotherapy have therapeutic limitations to improve survival in these patients. In this study, we used a differential proteomic-based technique (2D-Difference Gel Electrophoresis) coupled with matrix-assisted laser desorption/ionization-time of flight (MALDI-TOF) mass spectrometry to identify proteins that may serve as brain tumor antigens in new therapeutic assays. Five samples of patients presenting a GBM and five samples of microscopically normal brain tissues derived from brain epileptic surgery specimen were labeled and run in 2D-PAGE (Two-Dimensional Polyacrylamide Gel Electrophoresis) with an internal pool sample on each gel. Five gels were matched and compared with DIA (Difference In-gel Analysis) software. Differential spots were picked, in-gel digested and peptide mass fingerprints were obtained. Results From 51 protein-spots significantly up-regulated in GBM samples, mass spectrometry (MS) identified twenty-two proteins. The differential expression of a selected protein set was first validated by western-blotting, then tested on large cohorts of GBM specimens and non-tumor tissues, using immunohistochemistry and real-time RT-PCR. Conclusions Our results confirmed the importance of previously described proteins in glioma pathology and their potential usefulness as biological markers but also revealed some new interesting targets for future therapies.


Introduction
Malignant gliomas are the most common human primary brain tumors, glioblastoma multiforme (GBM) being the most aggressive and lethal form. Median survival of patients affected by GBM is around 15 months and nearly 75% of them will die within 18 months following diagnosis [1] despite treatment combining complete excision then chemo-and radiotherapy, pointing to the need for non-conventional therapies. Research has thus focused on specific strategies targeting intracellular signaling pathways [2], specific surface molecules, antiangiogenesis therapy [3], immunotherapy [4] and combinatorial approaches. Many reports identifying genes [5] and proteins which can be used to distinguish between GBM and either other gliomas or non tumorous brain tissue have increased our knowledge of this pathology, but a combination of genomic, proteomic and transcriptomic data will be needed to identify new therapeutic targets or biomarkers.
Proteomic approaches were most widely based on methods using differential expression on 2D-PAGE gels [6] or, more recently, two dimensional chromatography [7] followed by mass spectrometry protein identification. The 2D-DIGE (two dimensional difference gel electrophoresis) technology, using a mixed-sample internal standard, is now recognized as an accurate method to determine and quantify human proteins, reducing inter-gel variability and simplifying gel analysis. In the present study, GBM (T) and non tumorous (Nt) brain samples were labeled and fractionated using 2D-DIGE. Mass fingerprints of differentially expressed proteins were acquired with a MALDI-TOF/TOF spectrometer. Western-blot experiments, or immunohistochemical analysis for the largest series, were applied to confirm the up-regulation of selected proteins. The relationship with mRNA expression of coding genes for some over expressed proteins was analyzed by quantitative RT-PCR.

Results
Detection and identification of differentially-expressed proteins in GBM versus non-tumorous brain tissue 2D-DIGE was used to identify proteins differentially expressed in GBM versus non tumorous samples. The two-group experimental design ( Figure 1A) using a mixed internal standard enabled us to normalize spot volumes from each sample, to perform inter-gel comparisons and to recognize statistically significant interspot variations. A representative set of overlaid 2D-DIGE images gels is given in Figure 1B.
On each analytical 2D gel, an average of 700 spots corresponding to proteins with a pI between 3 and 10 and a molecular weight from 10 to 250 kDa was detected and 457 spots were matched across the 5 gels. A total of 51 spots showed statistically significant differences (Student's t-test p-value ≤ 0.05) in protein expression between the two populations with an average ratio > 2 and a statistical power of 0.88. These 51 spots were up-regulated in GBM. Only protein spots on silver or Coomassie Blue-stained gels manually matching with the corresponding DIGE images were processed. Thirty protein-spots were identified, corresponding to 22 different proteins (Table 1). This discrepancy is explained by the fact that some proteins, such as hemoglobin, DRP-2 (dihydropyrimidinase related protein 2, or collapsin response mediator protein 2, CRMP-2) or ER-60 (Endoplasmic reticulum resident protein 60 or protein disulfide-isomerase A3), are found in several spots while in other circumstances two proteins can contribute to the same spot. Results of DIGE analysis (average ratio glioblastoma/ normal brain tissue (GBM/NT), T-test), Mascot search (matched peptides, % of coverage and Mascot score) and protein characteristics (accession number, theoretical Mt and pI) are listed for each protein in Table 1.
GBM-overexpressed proteins included several blood proteins. Tumors spots for hemoglobin, apolipoprotein A1, serum albumin, and alpha-2 globulin were five times larger than the same spots from non-tumorous extracts.
Transthyretin was also four times more abundant in GBM. Two members of the dihydropyrimidinase family, DRP-2 and DRP-3 (or CRMP-4), were identified with an average GBM/NT ratio of 6.5 and 4.29 respectively. Three chaperone proteins HSP 27 (heat shock 27 kDa protein), GRP-75 (75 kDa glucose-related protein or heat-shock 70 kDa protein 9 also named mortalin), and ER 60 were also over-expressed in GBM.

Validation of 2D-DIGE results by western blot analysis
To validate the results of the DIGE analysis, five proteins, selected on the basis of interesting biological functions and high-fold changes, were tested in western blotting experiments with specific antibodies available commercially, on tumor protein extracts from GBM and control patients already used in 2D-DIGE. Equal amount of proteins from each sample were loaded. As shown in Figure 2, the expression levels for HSP 27 (29 kDa), ALDH (aldehyde dehydrogenase, 50 kDa), and Mn-SOD (20 kDa) were higher in GBM (n = 4 for Hsp27 and ALDH and n = 6 for MnSOD) in comparison to their expression in non-tumorous samples, thus confirming the DIGE analysis.
Two bands (55 kDa and 65 kDa) were found for DRP-2 and DRP-3. DRP-2 protein was expressed with a strong signal at 55 kDa in GBM (n = 2) and a slight signal at 65 kDa while the contrary was seen for nontumorous samples. DRP-3 protein expression was also different in GBM (n = 2) with a high signal at 65 kDa and a slighter one at 55 kDa and non-tumorous samples where no or faint bands could be found.

In situ validation of proteomic data by immunohistochemistry
Immunohistochemical studies were carried out with the five selected proteins already tested in western blot analysis, namely HSP 27, ALDH, Mn-SOD, DRP-2 and DRP-3, revealed in 25 GBM and three non-tumoral brain samples. The results confirmed their high expression in GBM. Figure 3 shows a representative paraffinembedded section and the respective non-tumorous control for each protein. Immunohistochemical studies gave the distribution of these proteins by cell category of the normal parenchyma: astrocytes, oligodendrocytes, neurons. All neurons (100%) were positive for ALDH, Mn-SOD, DRP-2 and DRP-3 ( Figure 4) except HSP 27 (0%). DRP-family molecules were strongly expressed in normal astrocytes and oligodendrocytes (51% and 52% positive cells respectively for DRP-2 and 68% and 75% for DRP-3) however this expression was higher in GBM cells (87% for DRP-2 and 88% for DRP-3). HSP 27 and ALDH were also more expressed in GBM (68% and 69% respectively) than in astrocytes (17% and 7%) or oligodendrocytes (23% and 2%). used non-tumorous brain samples (or "control") [9]. mRNA expression for GAPDH was identical for control and GBM and none of the nine non-tumorous samples tested against this control was considered as positive in our study (data not shown).
UCH-L1 expression was not increased in the GBM samples. ALDH, GRP-75 and DRP-2 disclosed mild positivity with respectively 4%, 4% and 6% of the GBM samples exhibiting an increase in mRNA expression more than 2-fold greater than the reference. In contrast, HSP 27, MnSOD and DRP-3 were clearly up-regulated with respectively 32%, 44% and 46% of the GBM samples having a 5-fold increase in mRNA expression compared with the pool of non-tumorous brain samples ( Figure 5).

Discussion
Extensive proteomic studies of human GBM have emerged in the last decade allowing comparison between GBM and normal brain proteomes and identification of differences in protein expression and profiles of diverse grades of gliomas (for review see [6,10]).
The aim of the present proteomic work was to highlight potential new glioblastoma tumor antigens thus extending the number of proteins already used and investigated, particularly in immunotherapy assays [9]. The 2D-DIGE technique was used to compare five GBM and five non-tumorous brain samples. Differential protein expression was analyzed by mass spectrometry. 2D-DIGE is a very powerful technique allowing controls and experimental samples to be run on the same gel together with an internal reference, each sample being labeled with a different cyanine. Using the DeCyder software, gel images can be analyzed to identify statistically significant differences in protein expression between  T9  T8  T7  T6   T9  T8  T7  T6   T9  T8  T7  T6  T6  different samples on a limited number of gels [11]. This technique is very useful for applications requiring accurate quantization and direct differential proteomic analysis of normal and abnormal tissues [12]. Twenty-two proteins were statistically highly expressed in GBM comparatively to non-tumorous brain samples. Some of them are involved in key metabolic pathways and suspected to participate in tumorigenic processes.
Blood and cytoskeletal proteins: two major groups upregulated in GBM Blood proteins emerged as a significant group highlighting the high angiogenesis capacity of GBM. GBM are highly vascular tumours and this study analyzed the proteome of the whole GBM sample extract, without removing blood vessels. In a recent review [13], up-regulation of serum albumin and apolipoprotein-A1 were commented as the consequence of the blood brain barrier degradation in malignant gliomas. These results reinforce previous reports describing hemoglobin [8], apolipoprotein A-1 or albumin [8,14] as brain tumor markers up-regulated in GBM.
Transthyretin (TTR), a plasma protein described as thyroxin and retinol transporter, has recently gained interest in neurobiology and seems to represent a key protein in nervous system physiology [15]. Although TTR is synthesized in the choroid plexus, there are several unclear points concerning its other precise localizations and functions. In cerebrospinal fluid, TTR levels are raised and lowered in Parkinson's and Alzheimer's diseases respectively. Furthermore a physiological variant, transthyretin Thr119met chain, is highly expressed in malignant areas of GBM [16]. The biological significance of this GBM-TTR overexpression and either its implication in gliomagenesis or its possible therapeutic role need further investigations.
Another significant class of differentially expressed proteins is represented by cytoskeletal proteins: vimentin, alpha (6,3) and beta (2) tubulin, beta actin and GFAP are three-to six-fold increased in GBM. Beta-tubulins and microtubule components are encoded by a multigene family whose expression patterns are complex. In tumor tissues, most isotypes exhibited an altered expression which might represent markers of sensitivity for drug response [17]. GFAP, a well-known specific astrocytic marker, has been suggested as a serum diagnostic marker for GBM [18]. It has been also disclosed in different forms in high-grade gliomas (grade III and IV) but its overexpression has been linked with the lower grade [19]. In the present work, GFAP was identified in three spots, all of them being up-regulated in GBM.
Mitochondrial and endoplasmic reticulum proteins: enzymatic and chaperone proteins involved in key pathways and dysregulated in GBM In eukaryotic cells, mitochondria and endoplasmic reticulum form two endomembrane networks which interact closely with each other to control metabolic flow, protein transport, intracellular signaling and cell death. In normal glial cells, besides their essential role of adenosine triphosphate (ATP) generation to produce energy, mitochondria are also involved in the regulation of cellular proliferation and apoptosis. The important consequences of mitochondrial dysfunction in glioma in these three areas was highlighted in a recent study [20].
It is not surprising to find some mitochondrial or endoplasmic reticulum proteins differentially expressed in GBM and in non-tumorous samples. In the present work, DIGE analysis and immunohistochemistry study revealed that two mitochondrial enzymes, ALDH and Mn-SOD, are highly expressed in GBM: on average in a series of 50 GBM, 69% tumor cells were ALDH positive and 81% Mn-SOD positive. RT-PCR analysis also revealed a high level of transcripts encoding for MnSOD. Mn-SOD expression has recently been suggested as a potential GBM prognostic marker [21] because the lack of overexpression is found in three years survival group patients compared to short-term survival patients. These results confirmed earlier work which correlated shortest median survival with high Mn-SOD enzyme expression level in GBM patients [22].
The second mitochondrial enzyme, ALDH, has not, to our knowledge, been reported as a brain tumor marker, but ALDH activity, detected commonly with ALDE-FLUOR assay, has been used for identification and isolation of adult stem cells and particularly neural stem cells [23]. Cancer stem cells -or tumor initiating cells (TICs)-, might be also identified in this way. Recently, ALDH expression, characterized by immunohistochemical staining of common epithelial cancers and their corresponding normal tissues, has been correlated with ALDH enzymatic activity. In tumors for which corresponding normal tissues expressed relatively low ALDH levels (breast, lung, ovarian or colonic cancer), ALDH can be proposed as a cancer stem cell marker [24]. Moreover, in breast and pancreatic cancers [25][26][27] the presence of positive cells has been related with poor survival, which was possibly explained by a greatest resistance to chemotherapeutic drugs as demonstrated for ovarian cancer cells [24]. In GBM, there is now much evidence that such tumor stem cells exist and despite several markers proposed for their identification (CD133, SSEA-1 [28], A2B5 [29]), no consensus has been reached. ALDH could therefore be a good candidate for best TIC enrichment in GBM.
GRP-75 (mortalin, MTHSP75), the mitochondrial chaperone Hsp70 kDa protein isoform 9, a major mitochondrial protein, plays a central role in protein import and export. It could also localize in endoplasmic reticulum and contribute to modulation of the stress response. GRP-75 was previously linked with malignant progression of lowgrade astrocytoma [30] and an increased expression of these chaperones contribute to tumorigenesis [31]. Our 2D-DIGE analysis confirmed a higher expression in GBM samples for GRP-75 and also two additional chaperone proteins: the reticulum endoplasmic protein disulfite isomerase (ER-60, ERp57 or GRP-57) and the heat-shock protein HSP 27. In brain, HSP 27 is linked with highgrade gliomas [32] or high malignancy human glioma cell line U87MGΔEGFR [33]. In cancer cells, HSP 27 is currently considered to participate in oncogenesis and in resistance to chemotherapy [34]. Conversely, ER-60 was previously reported as less expressed in GBM relatively to low-grade astrocytomas [32,35], although its key role in the control of newly synthesized glycoproteins during oncogenic transformation has been recently assessed [36].
Additionally, ER-60 is involved in the modulation of STAT-3 signaling and could contribute to a neoplastic diseased state. These findings highlight the importance of mitochondrial and reticulum endoplasmic proteins as targets for new therapies and suggest a particular study of these organelles via specific proteome analysis. In a recent paper, [13] a review of technical limitations in actual glioma proteomics shows the interest of sub-proteomes and particularly mitochondrial proteomes.

Dihydropyrimidinase proteins: importance of isoforms
Two dihydropyrimidinase proteins overexpressed in GBM were selected and retained for validation experiments by western blot, immunochemistry and RT-PCR analysis. DRP-2 (CRMP-2/Ulip2) and DRP-3 (CRMP-4/ Ulip 1) are members of collapsing response mediator proteins (CRMP) involved in regulation of neurite guidance and synapse formation; they are highly expressed during brain development and rarely in adult brains [37] except DRP-2 which can be detected until adulthood in several types of specialized neurons in the hippocampus, cerebellum and dorsal root ganglion. CRMP family members (DRP-1 to -5) have been essentially studied in the context of neurodegenerative diseases [38][39][40] where their expression is often altered. Nonetheless, some reports respectively outlined DRP-2 [8] down-expression and DRP-3 [14] up-regulation in GBM. In agreement with the later, DRP-3 which was identified in 2D-gels in only a single spot (65 kDa and 6.3 for MW and pI respectively) was confirmed up-regulated by western blot (65 kDa band), immunohistochemistry and RQ-PCR analysis. It was not so obvious for DRP-2: DIGEstudy of GBM individualized four overexpressed spots: two spots with an apparent mass of 65 kD and pI around 5.9 and two acidic 55 kDa forms with 5.3 pI. In western blots studies, a 65 kDa-DRP-2 was over expressed in normal brain while in GBM it represented a lighter isoform (55-60 kDa). RQ-PCR results revealed that protein overexpression was not linked to a coding gene up-regulation. In addition, immunochemistry, which could not differentiate DRP-2 isoforms, confirmed that 87% of GBM cells were positive compared to 51% of astrocytes and 52% of oligodendrocytes in normal brain. The DRP family is known to be highly phosphorylated [41,42]. Depending on study design and experimental techniques, the molecular weight of phosphorylated or non-phosphorylated forms of DRP-2 has been reported to vary from 55 kDa to 70 kDa [43,44]. Other post-translational modifications such as myristilation, methylation or esterification can alter the pI of the proteins, and glycosylation and prenylation can alter their molecular weight [8]. In addition, some light forms (55 kDa) could be breakdown products [45]. Our findings revealed that DRP-2 was overexpressed in GBM, and that two isoforms could discriminate GBM from non-tumorous brain tissue. Indeed, additional work is necessary to analyze and determinate the nature of translational modifications responsible for these isoforms and confirm the 55 KDa form as a potential glioma marker.
In our study, we present five proteins that are also regulated in neurodegenerative diseases: UCHL-1, transhtyretin, Apolipoprotein A1, DRP-2 and DRP-3. Altered expression of these proteins in Alzheimer's disease (AD), Parkinson's disease (PD) or Down Syndrome (DS) may reflect their importance in brain physiology or may be the result of stress generated by the disease. UCHL-1 is down-regulated in AD and in PD [46], transthyretin is raised in PD, but lowered in AD. High levels of Apolipoprotein A1 were found to be associated with neurodegeneration in AD. Concerning DRP-family, DRP-2 was increased in paired helical filaments and was related to the loss of neurofibrillary tangles-free neurons in AD ( [40]. In another work [38], DRP-2 was found downregulated at the mRNA level in DS and AD. In contrast, at the proteomic level, an increased ratio of 55 kDa to the 65 kDa forms in AD and DS brains suggested decreased phosphorylation of the DRP-2 protein. Differential expression in spots for DRP-2 and DRP-3 was reported in the study of R. Weitzdoerfer et al [39] in a fetal DS brain analysis: increased levels of one spot assigned to DRP-3 and decreased levels of spots assigned to DRP-2. More studies about these proteins and their isoforms would be necessary to conclude about their use as markers in GBM.
Overall, this work underlines the up-regulation of several GBM-related proteins; this was confirmed by western blot and immunohistochemistry techniques in some instances. These proteins might be attractive clinical biomarkers, linked or not with survival, or a specific GBM immunotherapy target. It could be also interesting to explore specifically the mitochondrial or endoplasmic reticulum GBM proteome which might open the way to novel fields of therapeutic advances.

Materials and methods
Samples GBM samples were obtained from patients admitted in the Neurosurgery Department of Rennes University Hospital and collected in accordance with the French regulations. For the entire study, 66 patients were included, with newly-diagnosed, untreated, primary glioblastoma (7 being classified as giant cell glioblastoma,) and two patients with a secondary glioblastoma. Twenty-three patients received only radiotherapy after surgical resection, when 27 patients received radiotherapy with concomitant temozolomide followed by adjuvant temozolomide (the current standard treatment). Other patients were treated after surgical resection with radiotherapy followed by different chemotherapeutic regimens.
2D-DIGE and western blot analysis were conducted on samples arising from five men (age range 43 -72 years, mean: 60 years), classified as primary glioblastoma. The validation cohort for immunochemistry and RT-PCR analysis was composed of 34 men and 30 women, (age range 35 -75 years mean: 59 years). Some patients were tested with both RT-PCR and immunochemistry. Non-tumorous brain tissues were obtained from normal areas (either grey or white matter) of brain tissues removed from patients undergoing non-tumor epileptic surgery. All samples were either conserved at -80°C for 2D-DIGE, real time RT-PCR and western blotting analysis, or formalin-fixed and paraffin-embedded for immunohistochemistry.

Protein extraction
Cell lysates were prepared from five GBM and five non tumorous brain samples by mechanical disruption in 2.5 volumes of ice-cold lysis buffer (Tris 20 mM, pH7.5, CHAPS 4%, urea 8M (Sigma-Aldrich, St Louis, USA) and antiproteases cocktail (Complete EDTA-free tablets, Roche Diagnostics, Mannheim, Germany)). Samples were sonicated (6 cycles of ten seconds with relapse of 30 seconds in ice-bath) and centrifuged (15 000 g, 30 minutes, 4°C). Supernatants were ultra centrifuged at 108,000 g for 60 minutes at 4°C. Protein concentration, in the resulting supernatants containing cytosolic protein extracts, was determined using the Bradford protein assay. Aliquots (100 μg) were conserved at -80°C.

Protein labeling with cyanin dyes
Cytosolic extracts were labeled with CyDyes DIGE Fluors developed for fluorescence 2-D DIGE technology (GE Healthcare, Bucks, UK) according to the manufacturer's recommended protocol. Briefly, 50 μg of each sample were minimally labeled with 400 pmol of aminereactive cyanine dyes, Cy3 or Cy5, on ice for 30 minutes, in the dark. GBM and non-tumorous samples were labeled with Cy5 or Cy3 in a random manner to avoid dye-specific protein labeling, as described in Figure 1A. An internal pool, labeled with Cy2 fluorescent dye, was generated by combining equal amounts of all cell cytosolic extracts and was included in all the gels run in this study. The labeling reaction was quenched by incubation, for 10 minutes with 1 μL of 10 mM lysine (Sigma-Aldrich, ST Louis, USA) on ice, in a darkroom.
Following the labeling reaction, GBM extracts and their random non-tumorous tissue counterparts were combined together with the internal pool, and Des-treak™ IEF buffer (GE Healthcare) was added to make up the volume to 450 μl prior to IEF (isoelectric focalisation) on five 24 cm gel strips.

Two-dimensional SDS-PAGE
A first focusing isolectric electrophoresis was carried out on IPGphor™ system (GE Healthcare). Pre-cast immobilized pH gradient strips (pH 3-10 NL, 24 cm) were used for this first-dimensional separation for a total focusing time of 60 kV-h. After IEF, the IPG strips were incubated two times under ambient temperature, for 15 minutes in an equilibration solution (0.05M Tris-HCl pH 8.8, 6M Urea, 30% glycerol, 2% SDS and bromophenol blue) containing 65 mM DTT and 250 mM iodoacetamide respectively. Strips were directly applied on top of pre-cast 12% SDS-PAGE gels (GE Healthcare) and run in a vertical Ettan DaltSix system (GE Healthcare) for approximately 5 hours. Five gels were processed simultaneously.

Gel imaging and data analysis
After SDS-PAGE, cyanine-labeled proteins were directly visualized using a Typhoon™ 9400 imager scan (GE Healthcare) in fluorescence mode. Cy2 images were scanned using a 488 nm laser and an emission filter of 520 nm. Cy3 images were scanned using a 532 nm laser and an emission filter of 580 nm. Cy5 images were scanned using a 633 nm laser and an emission filter of 670 nm. Each gel was scanned at 200 μm (pixel size) resolution and was processed using the DeCyder software V5.01 (GE Healthcare) allowing quantification, gel matching and statistical analyses.
The Differential In-gel Analysis module (DIA) was used for pair-wise comparison of the two samples (GBM and non-tumorous sample) on each gel, to exclude artifacts from gel images and differentially quantify the protein spots in the image. The Biological Variation Analysis module (BVA) was used to match the entire set of protein-spot maps from comparable gels simultaneously. Student's test (p < 0.05) was used for statistical analyses. Only spots with at least 2-fold changes in volume after normalization samples were defined as down-or up-regulated. The statistical power of the analysis was calculated similarly to Karp N. et al [47] and Engelen K. et al [48]. The standard deviation of the log 10 (standardized abundance) per condition was calculated for each spots that have been matched across the 5 gels of the analysis. The median of these standard deviations was calculated in each condition to estimate the global variance of the replicates. The statistical power was then calculated for each condition using a tool available online (http://udel.edu/~mcdonald/statttest.html [49]), for a 2-fold change (effect size = 0.301), alpha = 0.05 and 5 replicates per group. After 2D-DIGE imaging and analysis, gels were post-stained with silverstain [50] or Coomassie-blue. Gels were scanned (Image Scanner TM GE Healthcare) and stored in 1% acetic acid at 4°C until spot excision. Matching between silveror Coomassie-blue stained gels and fluorescence maps was performed manually and pick lists were generated using the Image Master™ 2D Elite software (GE Healthcare).
Peptide Mass fingerprints were acquired using a MALDI-ToF/ToF mass spectrometer (Ultraflex™; Bruker Daltonics, GmbH) and processed using the FlexAnalysis™ software (version 2.2; Bruker Daltonics, GmbH) for peak list generation and a first internal calibration with trypsin autodigestion peptides. Peak lists were then transferred to ProteinScape™ software (version 1.3; Bruker Daltonics, GmbH) for another automatic calibration based on a calibration list (related to the sample type and treatment) containing autolysis peaks and contaminants (keratins, polymers and background peaks). After re-calibration, an automatic trypsin and contaminants filtering and removal was performed in order to submit only m/z related to the sample and to obtain higher identification rates (Score-Booster). Only the monoisotopic masses of tryptic peptides were then used to query NCBInr sequence databases (May 2008, 6493741 sequences) using the Mascot search algorithm (Mascot server version 2.1.04; http://www. matrixscience.com). Search conditions were as follows: initial rather open mass window of 70 ppm for an internal calibration, one missed cleavage allowed, modification of cysteines by iodoacetamide and methionine oxidation as variable modifications. Results were scored using the probability-based Mowse score (the protein score is -10 × log (P) were P is the probability that the observed match is a random event. In our conditions, a score greater than 81 indicated a significant identification (p < 0.05).

Immunohistochemical procedure
Immunohistochemistry was performed on formalin-fixed and paraffin-embedded tissues, using 5-μm sections mounted on silanized slides. Antigen retrieval was performed using citrate buffer (pH 6.0) at 80°C Bain-marie (40 min) for DRP-3. The sections were incubated respectively at 20°for 30 and 60 minutes with diluted primary antibodies against HSP 27 and DRP-2 and overnight at 4°C with antibodies against DRP-3, ALDH, Mn-SOD (dilutions of 1:200, 1:25, 1:50, 1:200 and 1:800 for respectively HSP 27, DRP-2, DRP-3, ALDH, and Mn-SOD in antibody diluent of the Dako Cytomation Kit (Trappes, France)). A second incubation (60 minutes) with an anti-goat or anti-rabbit antibody (dilution 1:200, DAKO Cytomation) was followed by a peroxidase-staining procedure using the RTU Vectastain Elite ABC kit (Vector). Sections were counterstained with hematoxylin-eosin-safran. Microscopic analysis was performed using a Leitz-Diaplan microscope (Nuremberg, Germany). Negative controls were obtained in the absence of primary antibody. External positive controls were used for each staining: normal breast for HSP 27; dendate gyrus for DRP-2 and DRP-3; spleen for ALDH and myocard for Mn-SOD. A total of twenty-five GBM and three non-tumorous brain sampling were analyzed. For each antibody, 1000 tumor cells were counted and results were expressed as the percentage of cells with positive cytoplasmic staining in two different and most expressive analyzed areas.

Quantification of gene expression by real-time RT-PCR
Total RNA was extracted from frozen tissues samples: 50 GBM and 9 non-tumorous brain as previously described [9]. In brief, sections of tumor samples, free of necrotic areas, were processed with the Rneasy Plus mini kit (Qiagen, Courtaboeuf, France). cDNA was prepared from 0.5 μg purified RNA (High Capacity cDNA Reverse transcription kit, Applied Biosystems, Courtaboeuf, France). Real-time-PCR was performed with a spectrofluorometric thermal cycler (ABI prism 7900, Applied Biosystems), following the manufacturer's recommendations. Primers and probes for HSP 27 (HSPB1), DRP-2 (DPYSL2), DRP-3 (DPYSL3), ALDH (ALDH2), MnSOD (SOD2),UCHL-1 (UCHL1) and GRP-75 (HSPA9) genes were purchased from Applied Biosystems (Assays-on-demand). Each point of data was run in duplicate. To normalize the data, GAPDH was chosen as an endogenous control. The comparative Ct method was used to determine relative gene copy numbers using the formula 2-ΔΔCt with a pool of non-neoplastic brain samples as reference.