Skip to main content

Integration of transcriptome and proteome profiles in glioblastoma: looking for the missing link



Glioblastoma (GB) is the most common and aggressive tumor of the brain. Genotype-based approaches and independent analyses of the transcriptome or the proteome have led to progress in understanding the underlying biology of GB. Joint transcriptome and proteome profiling may reveal new biological insights, and identify pathogenic mechanisms or therapeutic targets for GB therapy. We present a comparison of transcriptome and proteome data from five GB biopsies (TZ) vs their corresponding peritumoral brain zone (PBZ). Omic analyses were performed using RNA microarray chips and the isotope-coded protein label method (ICPL).


As described in other cancers, we found a poor correlation between transcriptome and proteome data in GB. We observed only two commonly deregulated mRNAs/proteins (neurofilament light polypeptide and synapsin 1) and 12 altered biological processes; they are related to cell communication, synaptic transmission and nervous system processes. This poor correlation may be a consequence of the techniques used to produce the omic profiles, the intrinsic properties of mRNA and proteins and/or of cancer- or GB-specific phenomena. Of interest, the analysis of the transcription factor binding sites present upstream from the open reading frames of all altered proteins identified by ICPL method shows a common binding site for the topoisomerase I and p53-binding protein TOPORS. Its expression was observed in 7/11 TZ samples and not in PBZ. Some findings suggest that TOPORS may function as a tumor suppressor; its implication in gliomagenesis should be examined in future studies.


In this study, we showed a low correlation between transcriptome and proteome data for GB samples as described in other cancer tissues. We observed that NEFL, SYN1 and 12 biological processes were deregulated in both the transcriptome and proteome data. It will be important to analyze more specifically these processes and these two proteins to allow the identification of new theranostic markers or potential therapeutic targets for GB.


Glioblastoma (GB) is the most common and aggressive primary tumor in the adult brain. Despite years of research and numerous clinical trials, survival remains poor [1]. Progresses have been made in understanding the underlying biology of GB thank to work involving genotype-based approaches and proteome analyses. The cancer genome atlas (TCGA) analysis identified genetic events that appear to be important in human GBs, including (i) deregulation of growth factor signaling via amplification and mutational activation of receptor tyrosine kinase (RTK) genes; (ii) activation of the phosphatidyl inositol 3-kinase pathway; and (iii) inactivation of the p53 and retinoblastoma tumor suppressor pathways [2]. Genome-wide profiling studies have highlighted the existence of molecular subtypes of GB with distinct biological features and clinical correlates [3,4,5]. For example, Verhaak et al. [3] described four subtypes of GB: proneural, neural, classical and mesenchymal characterized by abnormalities in PDGFRA, IDH1, EGFR and NF1. However, the definition of a Verhaak subtype for a whole tumor has been questioned because GBs are very heterogeneous tumors and recent studies showed that different samples from the same tumor can be of different Verhaak subtypes [6, 7]. Proteome profiling of human GB samples also revealed a protein cluster (Huntingtin, HNF4α, c-Myc and 14-3-3ζ) that is differentially expressed in GB and might also serve as a diagnostic marker [8, 9]. The methylation status of MGMT has also been identified through omic analyses: this feature predicts sensitivity to temozolomide, an alkylating agent that is the current standard treatment for GB patients [10]. Another identified biomarker is the isocitrate dehydrogenase 1 (IDH1) mutation, that has been identified and has diagnostic applications as it helps in distinguishing primary from secondary GB [2].

Until recently, the behavior of GB has been studied through independent analyses of the transcriptome or of the proteome [11,12,13,14,15,16]. Joint transcriptome and proteome profiling may reveal new biological insight, and identify pathogenic mechanisms or therapeutic targets for GB therapy.

We report the analysis of GB biopsies from five patients involving RNA microarray and isotope-coded protein label (ICPL) technologies, part of the Grand Ouest Glioma Project, a translational project aiming to study the intratumoral heterogeneity in GB [11, 12, 15,16,17,18,19,20]. The transcriptome and the proteome of the GB tumor zone (TZ) were defined by comparison with the corresponding peritumoral brain zone (PBZ). The integrated transcriptome and proteome analysis was based on the four different approaches described by Haider and Pal [21]: (1) intersection of transcriptome and proteome data, (2) identification of the common biological processes altered in the two datasets (3) identification of the common functional pathways altered in the two datasets, (4) topological network methods, with the analysis of the transcription factor binding sites (TFBSs) present upstream from the open reading frames of the altered proteins identified by ICPL method.


Patient recruitment

The entire project was approved by the local institutional review board (CPP Ouest II) and the Direction Générale de la Santé (DGS). All patients included in this study were diagnosed for de novo GB (WHO 2007 classification) by a central committee of neuropathologists and gave their written informed consent prior to their enrolment. Five patients (Table 1) with both proteome and transcriptome analysis of their TZ and PBZ tissues were selected from the databank of the “Grand Ouest Glioma Project”. More detailed information on tissue samples characteristics can be found in our previous publications [11, 15, 16].

Table 1 Description of patients’ characteristics

GB and control brain sampling

For each patient, image-guided neuronavigation was used during pre-surgical planning to define TZ and PBZ samplings sites. The TZ sample (volume around 1 cm3) was then collected in the contrast-enhanced area of the tumor by computer-assisted, image-guided brain biopsies before the surgical resection of the tumor (Brainlab®, La Défense, France). Control brain samples were taken from the PBZ to be used as control samples for transcriptomic and proteomic analyses in radiologically normal, non-enhancing brain at least 1 cm from the contrast enhancing tumor. All the PBZ showed minimal genomic alteration (< 1%) and did not show tumor cell infiltration on histopathological analysis except for the PBZ of GB-10 [15].

Samples were transferred to the Department of Pathology of the University hospital of Angers, France, for pathological diagnostic, and to transcriptomic and proteomic platforms in the University hospital of Rennes, France for molecular analyses.

Transcriptome analysis

Transcriptome analyses of the tumor samples were performed as previously described in one of our previous publication on the transcriptomic platform Biosit, Rennes, France [16, 22]. In brief, total RNA was isolated from the GB samples using the NucleoSpin RNAII kit (Macherey–Nagel, Hoerdt, France) and RNA integrity (RNA integrity NC8) was assessed with an Agilent 2100 bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). Extraction of RNA microarray data was performed with an Agilent Whole Human Genome 4 × 44 K Microarray 15 Kit (Agilent technologies), according to the manufacturer’s recommendations.

Raw RNA data were then log2-transformed and normalized (quantile normalization and baseline transformation) using R v3.1.0. ( For the transcriptomic profile identification, we used a non-parametric rank product method to account for hybridization bias, allowing the identification up- or down-regulated genes in pooled GB tumor tissue by comparison to the pooled peritumoral brain samples, using the RankProd R package. RNAs were considered significantly differentially expressed if the false detection rate (FDR) was below 0.05 and the absolute fold-change between pooled data from TZ vs. PBZ was greater than 2.

Proteome analysis

The protocol for proteome analyses has already been described in our previous publications [11, 12]. TZ samples were analyzed using ICPL, that allows a high-throughput identification and quantification of a sample’s protein profile [23, 24]. Intact proteins were labeled with isotopic derivatives of nicotinic acid of different molecular weight, then subject to gel liquid chromatography and tandem mass spectrometry with an Esquire HCT Ultra PTM Discovery mass spectrometer, to identify and quantify proteins.

Peptides were identified by querying the human Swiss-Prot database with the Mascot search engine (v.2.2.07) applying a score above the identity threshold and a FDR < 1%. Differentially expressed proteins were identified in the TZ by comparison to the PBZ samples with a threshold > 1.41 for up-regulated proteins and < 0.71 for down-regulated proteins, which is above the calculated technical variation of the method [11].

Comparison of transcriptome and proteome data

We used four different methods to compare the results from the transcriptome and proteome analyses of the GB tumor tissues:

  1. (1)

    We performed a direct comparison of the intersection of transcripts and proteins found to be deregulated between TZ samples and their corresponding PBZ, to identify the overlap of direct features between transcriptome and proteome data. The comparison of the results from the transcriptome and proteome analyses of the GB tumor tissues was performed on pooled patient data and on paired patient to reduce the uncertainty due to variability between the patients.

  2. (2)

    We conducted an analysis of the biological and functional processes found to be altered in transcriptome and proteome data using DavidGenes ( The probability of alteration of the biological process was calculated using one-sided Fisher exact P-value and the False Discovery Rate (FDR) was calculated using one-sided Fisher exact P-value corrected for multiple comparisons. Biological and functional processes were considered significantly altered in each dataset with P and FDR < 0.05.

  3. (3)

    The functional pathways identified by both transcriptome and proteome data as being altered were identified using KEGG database ( Functional pathways were considered significantly altered in each dataset with P < 0.05 corrected for multiple comparisons using Benjamini–Hochberg method, and a fold-change > 2.

  4. (4)

    We looked for the presence of direct edges between transcripts and proteins with the identification of TFBSs in the regulatory region upstream from the open reading frames of the proteins identified as deregulated in GB. The main objective of this analysis was to identify factors that may bind to (and increase the expression of) the DNA encoding proteins identified in proteome analysis as being over-expressed. The Multi-genome Analysis of Positions and Patterns of Elements of Regulation search engine was used, running the set of proteins found to be deregulated in GB in our study (MAPPER2, A graphical summary of the analysis techniques is available in the Additional file 1: Figure S1.

Western blot analysis

TZ and PBZ samples (n = 11) were lysed in RIPA buffer containing a protease inhibitor cocktail and PMSF (Fisher Scientific, Illkirch, France) at 4 °C for 30 min. The lysates were clarified by centrifugation at 14,000g at 4 °C for 30 min. Protein concentrations were determined using the Pierce™ BCA protein assay kit (Fisher Scientific) with BSA as the standard, and equal samples of proteins (10 μg/lane) from the samples were resolved on a 7.5% SDS–polyacrylamide gel. The proteins were then electrotransferred onto PVDF membranes. After blocking in TBS blocking buffer (Fisher Scientific) at 4 °C overnight, blots were incubated with the respective primary antibodies [anti-actin (Merk Millipore, Guyancourt, France), anti-neurofilament light polypeptide (NEFL), anti-synapsin 1 (SYN1) and anti-topoisomerase I binding, arginine/serine-rich, E3 ubiquitin protein ligase (TOPORS) (CliniSciences, Nanterre, France)] for 2 h at room temperature. Horseradish peroxidase-conjugated secondary antibodies were then used and visualized with an enhanced chemoluminesence (ECL) reagent and the LAS4000 digital imaging system (Fisher Scientific).


Direct comparison of deregulated RNA and proteins

The transcriptome pooled analysis with 41,000 probes of the five TZ vs PBZ samples identified 478 mRNAs differentially expressed between TZ samples and their corresponding PBZ samples. A total of 437 non-redundant genes were identified; 101 genes were over-expressed in the florid TZ, and 300 genes were under-expressed (list of differentially expressed probes in Additional file 2: Table S1). Proteome analysis identified 584 non-redundant proteins, and 259 were quantified: 31 proteins were found to be up-regulated in the TZ in at least 3/5 patients (Full proteome data available in [11]).

The intersection between transcriptome and proteome data consisted of two genes: that for the NEFL and that for SYN1. They did not show the same deregulation in the two omic analyses: transcriptome analysis indicated that they are under-expressed in TZ whereas proteome analysis indicated that they are over-expressed. The Western blot analysis of expression of NEFL and SYN1 in a larger cohort of TZ samples and their corresponding PBZ (n = 11) confirmed the transcriptomic results; in most cases, an under-expression of NEFL and SYN1 proteins was observed in the TZ (9/11 for NEFL and 10/11 for SYN1) (Fig. 1).

Fig. 1

Analysis of NEFL, SYN1 and TOPORS expression in TZ biopsies and their corresponding PBZ (n = 11). a Distribution of densitometry data obtained for NEFL, SYN1 and TOPORS in 11 GB patients. Data are presented as TZ/PBZ ratios. b Example of western blot showing the expression of NEFL, SYN1, TOPORS and actin in TZ and PBZ samples

We also performed paired patient-specific comparisons to reduce the uncertainty due to variability between the patients. We observed similar trends with the pooled data TZ vs. PBZ comparisons with a low correlation rate (about 24%) between transcriptomic and proteomic data (Table 2).

Table 2 Patients-paired comparison of transcriptomic and proteomic data

Comparison of altered biological processes

The transcriptomic data indicated that 149 biological processes were altered in the tumor samples, and the proteome analysis indicated that 23 biological processes were altered (list in Additional file 3: Table S2). Twelve biological processes were found to be altered in both datasets, all repressed in the transcriptome analysis but enriched in the proteome analysis, except for the “regulation of biological quality” process, which was enriched in both datasets (Table 3).

Table 3 Biological processes found to be altered in both transcriptome and proteome datasets

Comparison of altered functional pathways

The analysis of altered functional pathways using the KEGG database found eight upregulated functional pathways in transcriptome analysis and two upregulated pathways in proteome analysis. There was no significant overlap between transcriptome and proteome analyses (Table 4).

Table 4 Functional pathways altered in transcriptome and proteome analyses in KEGG database

Identification of TFBSs

We identified six specific TFBSs present upstream from the open reading frames of the altered proteins identified (Table 5, Additional file 4: Table S3). Of interest, a binding site for the RING finger protein TOPORS was present upstream from the open reading frames of all altered proteins. TOPORS was not deregulated at the mRNA level in the transcriptomic analysis (Additional file 4: Table S3) and it was not identified in the proteome analysis. We performed the expression of TOPORS by Western blot analysis on 11 TZ biopsies and their corresponding PBZ. We didn’t observe TOPORS expression in PBZ while 7/11 TZ samples expressed this protein with a TZ/PBZ ratio of 1.14–2.85 (Fig. 1).

Table 5 Transcription factor binding sites (TFBSs) identified upstream from the open reading frames encoding the deregulated proteins


Although the comparison of transcriptomic and proteomic profiles has been done in several cancers, this study is the first of which we are aware to compare transcriptome and proteome data in a same cohort of GBs. Haider and Pal [21] reviewed the existing major approaches for joint analysis of transcriptome and proteome data. As recommended by this review, we directly compared the deregulated proteins and mRNAs, and then compared the functional processes and regulators identified in these transcriptome and proteome data sets.

There were few common features in the transcriptome and proteome data: they were the deregulation of two mRNAs/proteins (NEFL and SYN1) and 12 biological processes; they are related to cell communication, synaptic transmission and nervous system processes. These findings are consistent with previous reports [11, 16, 38]. No biological processes linked to tumorigenesis, for example cell cycle regulation, cell metabolism or cell motility, were found in both transcriptome and proteome analysis to be altered in GB. Some such processes were found in one dataset to be altered in GB, for example d-glutamine metabolism in the transcriptome and glycolysis/glyconeogenesis in the proteome [39, 40]. We observed the “pathogenic Escherichia coli infection” enrichment in GB through the proteome analysis. Its relation with the GB is presently unknown. Considering the involvement of genes in multiple biological processes, this enrichment could be an artefact. However, recent studies highlight the role of microbiota in gastric and breast cancer developments and may be present in GB [25, 26]. Further studies are needed to reply to this observation.

Although NEFL, SYN1 and twelve biological processes were in common between the two omic analyses, they did not show the same deregulation (except for the “regulation of biological quality” process): transcriptome analysis indicated that they are under-expressed in TZ whereas proteome analysis indicated that they are over-expressed. Mismatches of this type between such pairs of datasets has already been described [21, 25, 26]; there are several possible explanations including translation of mRNAs being up-regulated by RNA binding proteins and/or down-regulation of miRNA targeting these mRNAs. Also, the half-lives of mRNAs are very much shorter than those of proteins and protein stability may be affected by post-translational modifications like phosphorylation, acetylation and glycosylation [21, 27,28,29]. Note that Western blot analysis on more TZ/PBZ samples from other GB patients confirmed the deregulation of NEFL and SYN1 at the protein level but in most cases, an under-expression of these proteins was observed in TZ in line with transcriptomic data. This result highlights the importance to have large cohort of patients to limit the misinterpretation of the overall transcriptomic and proteomic data.

The low correlation between transcriptome and proteome data is not new and not restricted to cancer tissue [28,29,30,31,32,33]. For example, a correlation of only 17% was found between mRNAs and proteins in lung adenocarcinoma [28, 34]; in prostate cancer, the correlation between gene expression and protein levels is also poor to moderate [35]. Song et al. [36] performed proteomic profiling of eight GBs and their paired normal brain tissues and afterwards assessed overlap with RNA gene expression profiling from GEO and TCGA datasets of GBs. They found a correlation of only 2% between the differentially expressed proteins and genes from microarray. The low overlap between transcriptome and proteome data observed in our study can be explained by the general biological phenomena described above but also by analytical bias and GB-specific alterations (Fig. 2). On the analytical side, the choice of the analysis technique is crucial in omic studies and directly affects both the results and the feasibility of comparison between different datasets [37, 38]. RNA microarray techniques are the most widely used methods, allowing fast and accurate identification of mRNAs [39]; however, the number of probes on the microarray chip limits the extent of mRNA detection and probe set identification is a source of error in mRNA identification [40]. The ICPL method used here for proteome analysis allows assessment of protein levels but only a fraction of the proteome corresponding to the abundant proteins is analyzed [24]. Furthermore, only 60% of the identified proteins in one analysis were quantified. Another issue is that correlation coefficients are a crude method of measuring associations. For example, they do not account for interactions. However, because our analysis is based on only five samples it would not be meaningful to apply more advanced statistical methods. More than these analytical biases, GB tumor tissue possesses several specific properties that may influence the comparison between transcriptome and proteome data. GB is by definition an inter- and intra-heterogeneous tumor, which complicates the comparison between the two different analyses. The inter-heterogeneity is well defined at the level of the TZ through the identification of several GB subtypes. Recently, we described that this inter-heterogeneity was also present at the level of the PBZ [15, 18,19,20]. We identified two extratumoral microenvironments that can be encountered in the PBZ of GB patients: an extratumoral microenvironment containing GB-associated stromal cells (GASCs) with procarcinogenic properties and another containing GASCs without such properties [20]. As these cells may have their own specific signatures, this adds a level of complexity to make an inter-individual comparison between transcriptome and proteome data. Furthermore, the use of mirror samples to perform transcriptome and proteome analyses does not guarantee that the two samples are identical due to the intra-heterogeneity of GB [15, 22]. In our previous studies [16, 31], we confirmed this intratumoral heterogeneity at the transcriptomic and proteomic levels by comparing four regions of the GB (necrotic zone, TZ, interface zone between the tumor and the parenchyma and PBZ). The proteomic analysis generated a specific dataset of proteins for which a gradient of over-expression was observed from the periphery to the core of the GB [31]. At the transcriptome level, we observed that the molecular heterogeneity was much more important within tumors than between patients [16, 41]. The molecular definition of this intratumoral heterogeneity of GB is still incomplete, but the data are rapidly growing. For example, Nobusawa et al. [42] observed numerous tumor area-specific genomic imbalances, and our previous study as others reported inherent intratumor molecular subtype heterogeneity in GBs [7, 16, 43]. Other studies showed that RTK amplifications as well as MGMT status are heterogeneously distributed in GB [44,45,46]. More recently, next-generation sequencing techniques were able to highlight this intratumoral heterogeneity, to identify different clonal population of GB cells and understand their role in the recurrence [47, 48].

Fig. 2

Summary of potential bias when comparing transcriptome and proteome data. This figure summarizes various mechanisms that may alter translation, and lead to differences between transcriptome and proteome (Kozak sequence: initiating sequence for translation, located on the mRNA; non-sense read through: misreading of the mRNA in the opposite direction from 3′ to 5′)

Of interest, the search of TFBSs upstream from the open reading frames encoding the deregulated proteins that may explain the modification of their level of expression revealed a common binding site for TOPORS. We observed through Western blot analysis an expression of TOPORS in 7/11 TZ samples while no expression was evidenced in PBZ. The consultation of the Human Protein Atlas is consistent with our data where a higher expression of TOPORS was observed through immunohistochemistry analysis in 5/11 glioma specimens relative to brain tissue [49]. TOPORS is a RING finger protein that was identified originally as a topoisomerase I-binding protein and as p53-binding protein. TOPORS was shown to function as both a ubiquitin and SUMO E3 ligase for p53 [50, 51]. Its overexpression leads to a proteasome-dependent decrease in p53 [50]. Human TOPORS is located on chromosome 9p21, a region found frequently altered in several different malignancies of which GBs [52]. Some findings suggest that TOPORS may function as a tumor suppressor [53, 54]; further studies are needed to clarify the exact clinical significance as well as the exact biological function of TOPORS expression in GBs.


In this study, we showed a low correlation between transcriptome and proteome data for GB samples as described in other cancer tissues. We recognize that the number of studied samples was only five and that only a single sample per patient was used indicating that the results should be considered observational at this time. Future multi-omics studies must be performed on large cohort of patients and on spatially distinct tumor fragments per patient to consider the inter- and intra-heterogeneity of GBs. This may lead to an interpretation accuracy of the overall transcriptomic and proteomic data. We observed that NEFL, SYN1 and 12 biological processes were deregulated in both the transcriptome and proteome data. It will be important to analyze more specifically these processes and these two proteins to allow the identification of new theranostic markers or potential therapeutic targets for GB. Furthermore, a more detailed study of TOPORS which its TFBS was present upstream from the open reading frames of all proteins altered between TZ and PBZ may be promising.



European bioinformatics institute


flase detection rate


glioma-associated stroma cells




isotope-coded protein labeling


O6 methyl guanine methyl transferase


mass spectrometry


neurofilament light polypeptide


peritumoral brain zone


polyvinylidene difluoride


receptor tyrosine kinase


synapsin 1


the cancer genome atlas


transcription factor binding site


TOP1 binding arginine/serine rich protein


tumor zone


  1. 1.

    Stupp R, Mason WP, van den Bent MJ, Weller M, Fisher B, Taphoorn MJB, et al. Radiotherapy plus concomitant and adjuvant temozolomide for glioblastoma. N Engl J Med. 2005;352:987–96.

    CAS  Article  Google Scholar 

  2. 2.

    Cancer Genome Atlas Research Network. Comprehensive genomic characterization defines human glioblastoma genes and core pathways. Nature. 2008;455:1061–8.

    Article  Google Scholar 

  3. 3.

    Verhaak RGW, Hoadley KA, Purdom E, Wang V, Qi Y, Wilkerson MD, et al. Integrated genomic analysis identifies clinically relevant subtypes of glioblastoma characterized by abnormalities in PDGFRA, IDH1, EGFR, and NF1. Cancer Cell. 2010;17:98–110.

    CAS  Article  Google Scholar 

  4. 4.

    Phillips HS, Kharbanda S, Chen R, Forrest WF, Soriano RH, Wu TD, et al. Molecular subclasses of high-grade glioma predict prognosis, delineate a pattern of disease progression, and resemble stages in neurogenesis. Cancer Cell. 2006;9:157–73.

    CAS  Article  Google Scholar 

  5. 5.

    Li A, Walling J, Ahn S, Kotliarov Y, Su Q, Quezado M, et al. Unsupervised analysis of transcriptomic profiles reveals six glioma subtypes. Cancer Res. 2009;69:2091–9.

    CAS  Article  Google Scholar 

  6. 6.

    Eder K, Kalman B. Molecular heterogeneity of glioblastoma and its clinical relevance. Pathol Oncol Res. 2014;20:777–87.

    CAS  Article  Google Scholar 

  7. 7.

    Sottoriva A, Spiteri I, Piccirillo SGM, Touloumis A, Collins VP, Marioni JC, et al. Intratumor heterogeneity in human glioblastoma reflects cancer evolutionary dynamics. Proc Natl Acad Sci USA. 2013;110:4009–14.

    CAS  Article  Google Scholar 

  8. 8.

    Simeone P, Trerotola M, Urbanella A, Lattanzio R, Ciavardelli D, Di Giuseppe F, et al. A unique four-hub protein cluster associates to glioblastoma progression. PLoS ONE. 2014;9:e103030.

    Article  Google Scholar 

  9. 9.

    Malkki H. Neuro-oncology: proteomic profiling could facilitate glioblastoma diagnosis. Nat Rev Neurol. 2014;10:484.

    Article  Google Scholar 

  10. 10.

    Hegi ME, Liu L, Herman JG, Stupp R, Wick W, Weller M, et al. Correlation of O6-methylguanine methyltransferase (MGMT) promoter methylation with clinical outcomes in glioblastoma and clinical strategies to modulate MGMT activity. J Clin Oncol. 2008;26:4189–99.

    CAS  Article  Google Scholar 

  11. 11.

    Com E, Clavreul A, Lagarrigue M, Michalak S, Menei P, Pineau C. Quantitative proteomic Isotope-Coded Protein Label (ICPL) analysis reveals alteration of several functional processes in the glioblastoma. J Proteom. 2012;75:3898–913.

    CAS  Article  Google Scholar 

  12. 12.

    Lemée J-M, Com E, Clavreul A, Avril T, Quillien V, de Tayrac M, et al. Proteomic analysis of glioblastomas: what is the best brain control sample? J Proteom. 2013;85:165–73.

    Article  Google Scholar 

  13. 13.

    Mangiola A, Saulnier N, De Bonis P, Orteschi D, Sica G, Lama G, et al. Gene expression profile of glioblastoma peritumoral tissue: an ex vivo study. PLoS ONE. 2013;8:e57145.

    CAS  Article  Google Scholar 

  14. 14.

    Piccirillo SGM, Dietz S, Madhu B, Griffiths J, Price SJ, Collins VP, et al. Fluorescence-guided surgical sampling of glioblastoma identifies phenotypically distinct tumour-initiating cell populations in the tumour mass and margin. Br J Cancer. 2012;107:462–8.

    CAS  Article  Google Scholar 

  15. 15.

    Lemée J-M, Clavreul A, Aubry M, Com E, de Tayrac M, Eliat P-A, et al. Characterizing the peritumoral brain zone in glioblastoma: a multidisciplinary analysis. J Neurooncol. 2015;122:53–61.

    Article  Google Scholar 

  16. 16.

    Aubry M, de Tayrac M, Etcheverry A, Clavreul A, Saikali S, Menei P, et al. From the core to beyond the margin: a genomic picture of glioblastoma intratumor heterogeneity. Oncotarget. 2015;6:12094–109.

    Article  Google Scholar 

  17. 17.

    Lemée J-M, Clavreul A, Menei P. Intratumoral heterogeneity in glioblastoma: don’t forget the peritumoral brain zone. Neuro-Oncology. 2015;17:1322–32.

    Article  Google Scholar 

  18. 18.

    Clavreul A, Etcheverry A, Chassevent A, Quillien V, Avril T, Jourdan M-L, et al. Isolation of a new cell population in the glioblastoma microenvironment. J Neurooncol. 2012;106:493–504.

    Article  Google Scholar 

  19. 19.

    Clavreul A, Guette C, Faguer R, Tétaud C, Boissard A, Lemaire L, et al. Glioblastoma-associated stromal cells (GASCs) from histologically normal surgical margins have a myofibroblast phenotype and angiogenic properties. J Pathol. 2014;233:74–88.

    CAS  Article  Google Scholar 

  20. 20.

    Clavreul A, Etcheverry A, Tétaud C, Rousseau A, Avril T, Henry C, et al. Identification of two glioblastoma-associated stromal cell subtypes with different carcinogenic properties in histologically normal surgical margins. J Neurooncol. 2014;122:1–10.

    Article  Google Scholar 

  21. 21.

    Haider S, Pal R. Integrated analysis of transcriptomic and proteomic data. Curr Genom. 2013;14:91–110.

    CAS  Article  Google Scholar 

  22. 22.

    de Tayrac M, Etcheverry A, Aubry M, Saïkali S, Hamlat A, Quillien V, et al. Integrative genome-wide analysis reveals a robust genomic glioblastoma signature associated with copy number driving changes in gene expression. Genes Chromosomes Cancer. 2009;48:55–68.

    Article  Google Scholar 

  23. 23.

    Lottspeich F, Kellermann J. ICPL labeling strategies for proteome research. Methods Mol Biol Clifton NJ. 2011;753:55–64.

    CAS  Article  Google Scholar 

  24. 24.

    Schmidt A, Kellermann J, Lottspeich F. A novel strategy for quantitative proteomics using isotope-coded protein labels. Proteomics. 2005;5:4–15.

    CAS  Article  Google Scholar 

  25. 25.

    Xuan C, Shamonki JM, Chung A, Dinome ML, Chung M, Sieling PA, et al. Microbial dysbiosis is associated with human breast cancer. PLoS ONE. 2014;9:e83744.

    Article  Google Scholar 

  26. 26.

    Neto AG, Bradshaw A-D, Pei Z. Microbiome, a new dimension in cancer research. Ann Transl Med. 2015;3:229.

    PubMed  PubMed Central  Google Scholar 

  27. 27.

    Schwanhäusser B, Busse D, Li N, Dittmar G, Schuchhardt J, Wolf J, et al. Global quantification of mammalian gene expression control. Nature. 2011;473:337–42.

    Article  Google Scholar 

  28. 28.

    Maier T, Güell M, Serrano L. Correlation of mRNA and protein in complex biological samples. FEBS Lett. 2009;583:3966–73.

    CAS  Article  Google Scholar 

  29. 29.

    Ghazalpour A, Bennett B, Petyuk VA, Orozco L, Hagopian R, Mungrue IN, et al. Comparative analysis of proteome and transcriptome variation in mouse. PLoS Genet. 2011;7:e1001393.

    CAS  Article  Google Scholar 

  30. 30.

    Gygi SP, Rochon Y, Franza BR, Aebersold R. Correlation between protein and mRNA abundance in yeast. Mol Cell Biol. 1999;19:1720–30.

    CAS  Article  Google Scholar 

  31. 31.

    Com E, Boitier E, Marchandeau J-P, Brandenburg A, Schroeder S, Hoffmann D, et al. Integrated transcriptomic and proteomic evaluation of gentamicin nephrotoxicity in rats. Toxicol Appl Pharmacol. 2012;258:124–33.

    CAS  Article  Google Scholar 

  32. 32.

    Silvera D, Formenti SC, Schneider RJ. Translational control in cancer. Nat Rev Cancer. 2010;10:254–66.

    CAS  Article  Google Scholar 

  33. 33.

    Cuesta R, Gupta M, Schneider RJ. The regulation of protein synthesis in cancer. Prog Mol Biol Transl Sci. 2009;90:255–92.

    CAS  Article  Google Scholar 

  34. 34.

    Chen G, Gharib TG, Huang C-C, Taylor JMG, Misek DE, Kardia SLR, et al. Discordant protein and mRNA expression in lung adenocarcinomas. Mol Cell Proteomics. 2002;1:304–13.

    CAS  Article  Google Scholar 

  35. 35.

    Pascal LE, True LD, Campbell DS, Deutsch EW, Risk M, Coleman IM, et al. Correlation of mRNA and protein levels: cell type-specific gene expression of cluster designation antigens in the prostate. BMC Genom. 2008;9:246.

    Article  Google Scholar 

  36. 36.

    Song Y-C, Lu G-X, Zhang H-W, Zhong X-M, Cong X-L, Xue S-B, et al. Proteogenomic characterization and integrative analysis of glioblastoma multiforme. Oncotarget. 2017;8:97304–12.

    PubMed  PubMed Central  Google Scholar 

  37. 37.

    Louhimo R, Lepikhova T, Monni O, Hautaniemi S. Comparative analysis of algorithms for integration of copy number and expression data. Nat Methods. 2012;9:351–5.

    CAS  Article  Google Scholar 

  38. 38.

    Megger DA, Pott LL, Ahrens M, Padden J, Bracht T, Kuhlmann K, et al. Comparison of label-free and label-based strategies for proteome analysis of hepatoma cell lines. Biochim Biophys Acta. 2014;1844:967–76.

    CAS  Article  Google Scholar 

  39. 39.

    Heller MJ. DNA microarray technology: devices, systems, and applications. Annu Rev Biomed Eng. 2002;4:129–53.

    CAS  Article  Google Scholar 

  40. 40.

    Raghavachari N, Barb J, Yang Y, Liu P, Woodhouse K, Levy D, et al. A systematic comparison and evaluation of high density exon arrays and RNA-seq technology used to unravel the peripheral blood transcriptome of sickle cell disease. BMC Med Genom. 2012;5:28.

    CAS  Article  Google Scholar 

  41. 41.

    Aubry M, de Tayrac M, Etcheverry A, Clavreul A, Saikali S, Menei P, et al. Correction: From the core to beyond the margin: a genomic picture of glioblastoma intratumor heterogeneity. Oncotarget. 2016;7:67685.

    Article  Google Scholar 

  42. 42.

    Nobusawa S, Lachuer J, Wierinckx A, Kim YH, Huang J, Legras C, et al. Intratumoral patterns of genomic imbalance in glioblastomas. Brain Pathol. 2010;20:936–44.

    CAS  PubMed  Google Scholar 

  43. 43.

    Teng J, Carla da Hora C, Kantar RS, Nakano I, Wakimoto H, Batchelor TT, et al. Dissecting inherent intratumor heterogeneity in patient-derived glioblastoma culture models. Neuro-Oncology. 2017;19:820–32.

    CAS  PubMed  PubMed Central  Google Scholar 

  44. 44.

    Little SE, Popov S, Jury A, Bax DA, Doey L, Al-Sarraj S, et al. Receptor tyrosine kinase genes amplified in glioblastoma exhibit a mutual exclusivity in variable proportions reflective of individual tumor heterogeneity. Cancer Res. 2012;72:1614–20.

    CAS  Article  Google Scholar 

  45. 45.

    Szerlip NJ, Pedraza A, Chakravarty D, Azim M, McGuire J, Fang Y, et al. Intratumoral heterogeneity of receptor tyrosine kinases EGFR and PDGFRA amplification in glioblastoma defines subpopulations with distinct growth factor response. Proc Natl Acad Sci USA. 2012;109:3041–6.

    CAS  Article  Google Scholar 

  46. 46.

    Parker NR, Hudson AL, Khong P, Parkinson JF, Dwight T, Ikin RJ, et al. Intratumoral heterogeneity identified at the epigenetic, genetic and transcriptional level in glioblastoma. Sci Rep. 2016;6:22477.

    CAS  Article  Google Scholar 

  47. 47.

    Kim H, Zheng S, Amini SS, Virk SM, Mikkelsen T, Brat DJ, et al. Whole-genome and multisector exome sequencing of primary and post-treatment glioblastoma reveals patterns of tumor evolution. Genome Res. 2015;25:316–27.

    CAS  Article  Google Scholar 

  48. 48.

    Favero F, McGranahan N, Salm M, Birkbak NJ, Sanborn JZ, Benz SC, et al. Glioblastoma adaptation traced through decline of an IDH1 clonal driver and macro-evolution of a double-minute chromosome. Ann Oncol. 2015;26:880–7.

    CAS  Article  Google Scholar 

  49. 49.

    Uhlén M, Fagerberg L, Hallström BM, Lindskog C, Oksvold P, Mardinoglu A, et al. Proteomics. Tissue-based map of the human proteome. Science. 2015;347:1260419.

    Article  Google Scholar 

  50. 50.

    Rajendra R, Malegaonkar D, Pungaliya P, Marshall H, Rasheed Z, Brownell J, et al. Topors functions as an E3 ubiquitin ligase with specific E2 enzymes and ubiquitinates p53. J Biol Chem. 2004;279:36440–4.

    CAS  Article  Google Scholar 

  51. 51.

    Weger S, Hammer E, Heilbronn R. Topors acts as a SUMO-1 E3 ligase for p53 in vitro and in vivo. FEBS Lett. 2005;579:5007–12.

    CAS  Article  Google Scholar 

  52. 52.

    Bredel M, Bredel C, Juric D, Harsh GR, Vogel H, Recht LD, et al. High-resolution genome-wide mapping of genetic alterations in human glial brain tumors. Cancer Res. 2005;65:4088–96.

    CAS  Article  Google Scholar 

  53. 53.

    Saleem A, Dutta J, Malegaonkar D, Rasheed F, Rasheed Z, Rajendra R, et al. The topoisomerase I- and p53-binding protein topors is differentially expressed in normal and malignant human tissues and may function as a tumor suppressor. Oncogene. 2004;23:5293–300.

    CAS  Article  Google Scholar 

  54. 54.

    Oyanagi H, Takenaka K, Ishikawa S, Kawano Y, Adachi Y, Ueda K, et al. Expression of LUN gene that encodes a novel RING finger protein is correlated with development and progression of non-small cell lung cancer. Lung Cancer. 2004;46:21–8.

    Article  Google Scholar 

Download references

Authors’ contributions

MA, EC, MdT and JM analyzed and produced raw data for transcriptomics (MA, MdT, JM) and proteomics (EC). JML, AC and PM treated and analyzed the raw transcriptomics and proteomics data. AC made the western Blot analysis. JML wrote the manuscript. Each authors reviewed the manuscript. JM and PM supervised this work. PM obtained the funding for this work. All authors read and approved the final manuscript.


We gratefully acknowledge the contribution of the Cancéropole Grand Ouest. We also thank neurosurgeons at the University hospitals of Angers, Rennes, Tours and Brest (France) for supplying tissue samples.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

Transcriptomic data have been deposited in the ArrayExpress repository under the E-MTAB-5098, E-MTAB- 5092 and E-MTAB-5109 accession numbers.

Raw MS proteomic data have been converted and submitted to the EBI (European Bioinformatics Institute: and are accessible via the PRIDE repository ( under the project name: “intra-tumoral heterogeneity of glioblastoma multiforme” (Accession Numbers: 19502–19641) with a reviewer access (username: review81913 and password: %tm.jbqN).

Consent for publication

Not applicable.

Ethics approval and consent to participate

The experimental design of this study was approved by the local ethics committee (CPP Ouest II) and the Direction Générale de la Santé (DGS).

All patients gave their written informed consent prior to their enrolment in this study and consent for the use of the anonymized acquired data for research and publication purposes.


This work was funded by the French national health research institute (INSERM) and the Canceropole Grand Ouest, for the cost of omics analyses and the tissue samples conservation.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information



Corresponding author

Correspondence to Jean-Michel Lemée.

Additional files

Additional file 1: Figure S1.

Graphical abstract.

Additional file 2: Table S1.

Differentially expressed probes between TZ and PBZ pooled samples transcriptomic analysis.

Additional file 3: Table S2.

List of altered biological processes in pooled GB transcriptomic and proteomic profiles.

Additional file 4: Table S3.

List of identified TFBS located upstream of deregulated proteins’ coding region.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Lemée, J., Clavreul, A., Aubry, M. et al. Integration of transcriptome and proteome profiles in glioblastoma: looking for the missing link. BMC Molecular Biol 19, 13 (2018).

Download citation


  • Glioblastoma
  • Molecular biology
  • Transcriptomics
  • Proteomics
  • SYN1
  • NEFL