Reasearch Awards nomination

Email updates

Keep up to date with the latest news and content from JTM and BioMed Central.

Open Access Research

Interactions of miR-323/miR-326/miR-329 and miR-130a/miR-155/miR-210 as prognostic indicators for clinical outcome of glioblastoma patients

Shuwei Qiu12, Sheng Lin3, Dan Hu24, Yimin Feng2, Yang Tan2 and Ying Peng1*

Author Affiliations

1 Department of Neurology, The Sun Yat-sen Memorial Hospital, Sun Yat-sen University, Guangzhou, 510120, China

2 Department of Internal Medicine, College of Medicine, East Tennessee State University, Johnson City, 37604, TN, USA

3 Laboratory of Medical Genetics, Key Laboratory of Birth Defects Surveillance and Intervention of Guangdong Province, Shenzhen Research Institute of Population and Family Planning, Shenzhen, 518040, China

4 Department of Neurology, Renmin Hospital of Wuhan University, Wuhan, 430060, China

For all author emails, please log on.

Journal of Translational Medicine 2013, 11:10  doi:10.1186/1479-5876-11-10


The electronic version of this article is the complete one and can be found online at: http://www.translational-medicine.com/content/11/1/10


Received:3 November 2012
Accepted:7 January 2013
Published:9 January 2013

© 2013 Qiu et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Glioblastoma multiforme (GBM) is the most common and aggressive brain tumor with poor clinical outcome. Identification and development of new markers could be beneficial for the diagnosis and prognosis of GBM patients. Deregulation of microRNAs (miRNAs or miRs) is involved in GBM. Therefore, we attempted to identify and develop specific miRNAs as prognostic and predictive markers for GBM patient survival.

Methods

Expression profiles of miRNAs and genes and the corresponding clinical information of 480 GBM samples from The Cancer Genome Atlas (TCGA) dataset were downloaded and interested miRNAs were identified. Patients’ overall survival (OS) and progression-free survival (PFS) associated with interested miRNAs and miRNA-interactions were performed by Kaplan-Meier survival analysis. The impacts of miRNA expressions and miRNA-interactions on survival were evaluated by Cox proportional hazard regression model. Biological processes and network of putative and validated targets of miRNAs were analyzed by bioinformatics.

Results

In this study, 6 interested miRNAs were identified. Survival analysis showed that high levels of miR-326/miR-130a and low levels of miR-323/miR-329/miR-155/miR-210 were significantly associated with long OS of GBM patients, and also showed that high miR-326/miR-130a and low miR-155/miR-210 were related with extended PFS. Moreover, miRNA-323 and miRNA-329 were found to be increased in patients with no-recurrence or long time to progression (TTP). More notably, our analysis revealed miRNA-interactions were more specific and accurate to discriminate and predict OS and PFS. This interaction stratified OS and PFS related with different miRNA levels more detailed, and could obtain longer span of mean survival in comparison to that of one single miRNA. Moreover, miR-326, miR-130a, miR-155, miR-210 and 4 miRNA-interactions were confirmed for the first time as independent predictors for survival by Cox regression model together with clinicopathological factors: Age, Gender and Recurrence. Plus, the availability and rationality of the miRNA-interaction as predictors for survival were further supported by analysis of network, biological processes, KEGG pathway and correlation analysis with gene markers.

Conclusions

Our results demonstrates that miR-326, miR-130a, miR-155, miR-210 and the 4 miRNA-interactions could serve as prognostic and predictive markers for survival of GBM patients, suggesting a potential application in improvement of prognostic tools and treatments.

Keywords:
Glioblastoma multiforme; microRNA; Prognostic marker; Overall survival; Progression-free survival; Interaction

Background

Glioblastoma multiforme (GBM) is the most common and aggressive primary adult brain tumor. Despite advances in treatment modalities, the prognosis of GBM patients is very poor [1]. Therefore, it is urgent to develop new diagnostic and prognostic tools and treatments which may be beneficial for improving the clinical management of GBM. Currently, tumor stratifications relying on molecular profiles are increasingly prevalent and important. Furthermore, molecular and genetic profiling studies have identified several prognostic and predictive markers for GBM [2,3].

MicroRNAs (miRNAs or miRs), which are endogenous non-coding small RNAs, post-transcriptionally regulate gene expression through inhibition of translation or degradation of target mRNAs [4]. MiRNAs are aberrantly expressed in a variety of tumor types and exert important regulations on tumor biology via acting as oncogenes or tumor suppressors [5]. Recently, several studies indicate that expressions of miRNAs are associated with patients’ survival and are able to function as prognostic and predictive indicators [6,7]. Moreover, it has been confirmed that miRNA expression profiles are more accurate to classify tumors than mRNA profiles [8]. However, in the selection of miRNA markers for GBM prognosis, the applying of following aspects, such as small dataset, explanatory variables, single miRNA analysis, pre-selection of miRNAs and use of approaches, finally lead to a variety set of different miRNA markers.

The main purpose of this study is to identify specific miRNA markers that are closely associated with tumor progression and survivals for GBM patients by analyzing significantly altered miRNAs in a large dataset. Another goal is to investigate the availability and rationality of interactions of interested miRNAs as prognostic and predictive indictors for clinical outcome of GBM patients. In this study, we found that miR-326, miR-130a, miR-155, miR-210 and 4 miRNA-interactions could function as prognostic and predictive markers for survival of GBM patient.

Materials and methods

TCGA miRNA dataset and patient information

Expression data of miRNAs and genes and the corresponding clinical data for glioblastoma patients were downloaded from The Cancer Genome Atlas data portal (July 2012) [9]. 480 GBM patients with full annotation of Age, Gender, Survival time, Vital status, Time to progression/recurrence and miRNA values were identified in this study. There were 186 females and 294 male patients with ages 56.7 ± 15.9 and 58.1 ±13.6 years, respectively. Among the whole set, 318 patients suffered from tumor recurrence while 162 patients were kept away from progression/recurrence. Besides, miRNA expression data from 10 normal brain tissues (NBT) were also collected. The collection of the original material and data of TCGA was conducted in compliance with all applicable laws, regulations and policies for the protection of human subjects, and necessary IRB approvals were obtained [9]. Totally, Expression levels of 534 human microRNAs were detected using the Agilent 8 × 15K Human microRNA platform. The data was quantile-normalized, collapsed within miRNAs, and log2 transformed.

Analysis of expression levels of miRNA data in GBM samples

To reflect the real expression profiles of miRNAs in GBM and normal brain tissues, and to yield the detection error, we analyzed three separate batches, batch 5, batch 16 and batch 20, from TCGA dataset with the biggest numbers of patients, incorporating 63, 47 and 46 patient samples respectively. The level 3 data were directly used to evaluate the relative expression levels of miRNAs in each sample by the Z-score method [10]. The mean Z-score of each miRNA from GBM and NBT samples was calculated and sorted in each batch according to Z-score values in GBM samples, and the top 20 most down-regulated and up-regulated miRNAs were illustrated in the heatmaps pattern. Only those microRNAs overlapping for two or three times among three batches were selected.

Statistical computing methods

MiRNA expression profiles related to glioblastoma survival were identified using the Kaplan-Meier survival analysis and statistical significances of overall survival (OS) and Progression-Free survival (PFS) were determined using the Log-Rank test. Survival analysis was performed on SPSS (version 17.0; SPSS Inc.) and the survival curve was generated by GraphPad Prism (version 5.04; GraphPad Software, Inc.). For survival analysis, patients with survival time lesser than 30 days were excluded, since these patients might have died for reasons other than the disease itself. A total of 458 patients fitting this criterion were included for survival analysis. For stratification analysis of survival, expression levels of down-regulated miRNAs were sorted by ascending order, while the up-regulated miRNAs were sorted by descending order. Then, quartiles of 25%, 50% and 75% of the sorted miRNA values were set as cutoffs for low/high expressions of each miRNA. The survival time was expressed as mean ± SE. To determine whether expression levels of miRNAs were associated with Time to Progression (TTP), discrepancies of miRNA levels between groups were tested by student’s t-test, in terms of TTP of 9, 12 and 15 months. Likewise, according to whether there was recurrence or not, differences of miRNA levels between recurrent and non-recurrent groups were also tested using student’s t-test.

Next, the Cox proportional hazard regression model was used to determine the influences of miRNA expressions as well as clinicopathological factors (age, gender and recurrence) on patient survival. To adjust this potential effect that may be confounded by age, gender and recurrence, a multivariate Cox proportional hazard regression analysis using all these clinicopathological factors was performed. For analysis of interactions of two miRNAs, different quartile stratifications of the expression levels of miRNAs were set as cut-offs for high and/or low levels. The meaning of interaction in this text was defined as combined effect or coaction of high and/or low levels of two different miRNAs from one patient, e.g. high/low level of miR-A and high/low level of miR-B.

Functional Gene Ontology (GO) biological processes terms of the putative targets of candidate miRNAs were performed and the GO enrichment scores of target genes were presented in the form of heatmaps [11,12]. Putative targets of miRNAs were predicted by TargetScan Human 6.2 software [13]. Then the Functional Gene Ontology (GO) biological processes terms of the putative targets of candidate miRNAs, as well as the statistical analysis, were performed by DAVID Bioinformatics Resources 6.7 [14]. The network of validated targets of interested miRNAs was created by BisoGenet software, showing the interactions of validated targets [15]. The validated targets of interested miRNAs were obtained from Diana TarBase v6.0 [16] and miRecords [17] and extensive Pubmed publications search. The differences were considered statistically significant at p < 0.05.

Results

Screening of the most altered miRNAs expressed in GBM samples

To decrease the detection error within and cross microarrays, mean Z-score values of each miRNA in bath5, 16 and 20, which incorporate the biggest numbers of GBM samples in TCGA dataset, were calculated and sorted, and then the top 20 most altered miRNAs in GBM samples were illustrated in the form of heatmaps. After overlap analysis, expression levels of 19 miRNAs were found to be decreased in GBM sam-ples, including miR-124a, miR-128a, miR-128b, miR-129, miR-132, miR-137, miR-139, miR-203, miR-218, miR-323, miR-326, miR-329, miR-330, miR-383, miR-433, miR-485-5p, miR-491, miR-628, and miR-769-5p, while 19 miRNAs increased in GBM samples, involving let-7i, miR-106a, miR-106b, miR-130a, miR-130b, miR-155, miR-15a, miR-15b, miR-16, miR-195, miR-21, miR-210, miR-23a, miR-25, miR-27a, miR-605, miR-92, miR-92b, and miR-93 (Figure 1). Expression levels of all these candidate miRNAs were validated by other independent miRNAs detection studies with RNA sequencing, real-time PCR or microarray [18-22].

thumbnailFigure 1. Screening of miRNAs differentially expressed in GBM samples. Shown are top 20 most down-regulated and up-regulated miRNAs in GBM samples within batch 5, 16 and 20 from TCGA dataset. MiRNAs overlapped no less than two times cross 3 batches are selected as candidates for further studies. Red denotes low levels of miRNAs in GBM samples, while green denotes low levels of miRNAs in normal brain tissues.

Correlations between miRNA expression and survival of GBM patients

Next, using Kaplan-Meier survival analysis, correlations between miRNA expression profiles and patient survival were performed. Expression levels of the 38 aberrantly expressed miRNAs were stratified by quartiles 25%, 50% and 75%. Totally, 16 miRNAs were shown to correlate with OS of GBM patients at different quartile stratifications (Additional file 1: Table S1). Only miR-323, miR-326, miR-329, miR-130a, miR-155 and miR-210 were selected as interested candidates for further analysis; this was because, according to our knowledge, no publications on patients survival associated with these miRNAs were reported. Our analysis showed that low expression levels of miR-323, miR-329, miR-155 and miR-210 significantly correlated with long OS survival (p = 0.0043, 0.0182, 0.0191 and 0.0077, respectively), while high expression levels of miR-326 and miR-130a were associated with long OS survival (p = 0.0377 and 0.0099; Figure 2A). We also showed that high levels of miR-326 and miR-130a were associated with extended survival without tumor progression (p = 0.036 and 0.0098; Figure 2B), while low expression levels of miR-155 and miR-210 were related with long PFS survival (p = 0.0055 and 0.0212; Figure 2B). These results suggested that miR-326 and miR-130a may serve as tumor suppressors, while miR-323, miR-329, miR-155 and miR-210 serve as oncogenes.

Additional files 1. Table S1. Parameter estimates of microRNAs associated with survival time at different quartile stratifications.

Format: DOC Size: 108KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

thumbnailFigure 2. Kaplan-Meier survival curves of the interested miRNAs in GBM patients. A. Analysis of the interested miRNAs on patients’ Overall survival (OS). B. Analysis of the interested miRNAs on Progression-free survival (PFS). Survival time displayed here is mean ± SE.

The association among miRNAs expression and recurrence and time to progression in GBM patients

Both tumor recurrence and time to progression (TTP) are closely related to clinical outcome. Therefore, miRNA expression levels were evaluated to identify the differently expressed miRNAs according to tumor recurrence as well as long/short time to progression. Initially, we found that, among the 38 aberrant miRNAs, levels of miR-323 and miR-329 from the non-recurrent group (n = 162) were significantly increased in comparison to the recurrent group (n = 318; p = 0.0314 and 0.0364, respectively), whereas miR-132, miR-433, miR-628 and miR-769-5p were marginally increased in the non-recurrent group with miR-92b marginally decreased (Table 1). Then we analyzed expression discrepancies of miRNAs between long and short TTP groups. The results showed that significant up-regulations of miR-128b and miR-323 were observed in the group of TTP longer than 15 months (p = 0.0420 and 0.0451), with marginal increase of miR-329 level (p = 0.0981; Table 2). Likewise, among groups of TTP more than 9, 12 or 15months, obvious elevations of miR-106a, miR-106b and miR-92 were exhibited (p < 0.05), while levels of miR-21 was significantly reduced in the group of TTP more than 9 and 12months (p < 0.01; Table 2). However, no significant associations were observed between expression levels of miR-326/miR-130/miR-155/miR-210 and recurrence/TTP time.

Table 1. Differential expression levels of microRNAs associated with progression/recurrence or not in GBM patients

Table 2. Differential expression levels of microRNAs associated with time to progression in patients with GBM

The association between interactions of miRNAs and survival of GBM patients

As shown above, Low levels of miR-323 and miR-329 were associated with long survival, whereas high levels of both miRNA were present in patients without recurrence or long TTP time. This inconsistence highlighted the complication and importance of these two miRNAs in GBM tumor progression. Moreover, interactions of miRNAs may be more specific and feasible to discriminate and predict the potential survival time of GBM patients. Therefore, interactions of miR-323/miR-326/miR-329 and miR-130a/miR-155/miR-210 on survival were initially performed using Kaplan-Meier analysis with the aim at obtaining more specific predictors. The OS curves showed that interactions of two miRNAs were more potent to discriminate survival time of the same patients (Figure 3A). For instance, the OS of patients with low miR-323 and high miR-130a were longer than that of patients with high miR-323 and low miR-130a, with the survival gap of 16.5 months (p = 0.0007; Figure 3 A1). Also, the interaction of high miR-326 and high miR-130a was associated with long OS while that of low miR-326 and low miR-130a was not, with the survival span of 21.5 months (p = 0.0003; Figure 3 A4). Furthermore, PFS analysis showed that several interactions of two miRNAs were more effective to distinguish survival times of patients (Figure 3B). For example, the interaction of high miR-326 and high miR-130a was associated with much longer PFS survival than that of low miR-326 and low miR-130a, with the survival gap of 15.8 months (Figure 3 B2). Collectively, these data indicated that interactions of different levels of miRNAs significantly correlated with OS and/or PFS, and it was potent and effective to discriminate and predict the survival time for GBM patients by this method.

thumbnailFigure 3. Kaplan-Meier survival analysis of the miRNA-interactions. A. Analysis of the miRNA-interaction on patients’ OS. B. Analysis of the miRNA-interaction on patients’ PFS.

Functional analyses of the interested miRNAs in GBM

The availability of interaction of miRNAs for prognosis survival was confirmed as above. Subsequently, the rationality was further analyzed. First, GO biological processes of putative targets of miR-323 and miR-130a were enriched in apoptosis, cell proliferation, cell cycle, and cell adhesion/migration, indicating their complicated biological functions during tumor progression, while the enrichment of biological processes in cell adhesion of targets of miR-329 suggested that it mainly took part in tumor migration. Biological processes of target genes of miR-326 and miR-155 were dominantly associated with apoptosis (Figure 4A). Likewise, analysis of KEGG pathway showed that target genes of miR-323, miR-326, miR-130a and miR-155 were all enriched in Pathways in cancer, while miR-329 was mainly related with Adherens junction (Additional file 2: Table S2). However, there was no enrichment of biological process and KEGG pathway of target genes of miR-210 due to less predictive targets according to TargetScan (Figure 4A and Additional file 2: Table S2). Second, according to Diana TarBase and miRecords, a total of 45 validated target genes of the four miRNAs were observed, whereas the network was built with only 18 real target genes of miR-130a, miR-155 and miR-210, which exhibited direct protein-protein interactions (Figure 4B). This simple net was composed of molecules associated with apoptosis, cell proliferation, inflammation or carcinogenesis, which illustrated that targets of the interested miRNAs were not independent and they interacted with each other. All in all, the complication and cross-talk of the six miRNAs suggested the rationality of interactions of miRNAs to predict the clinical outcome.

Additional files 2. Table S2. KEGG pathway analysis of putative targets of candidate miRNAs.

Format: DOCX Size: 20KB Download fileOpen Data

thumbnailFigure 4. Functional analysis of the interested miRNAs in GBM. A. Functional GO biological processes of putative target genes of the 6 miRNAs enriched in cell proliferation and growth, cell cycle, apoptosis and cell adhesion and migration are illustrated. Red denotes highly significant correlation with miRNAs; black denotes lowly significant correlation with miRNAs; grey denotes no correlation. All the GO terms were selected at the criteria of p value < 0.01. B. Network of validated target genes of miRNAs. Only the direct interactions between these target genes are exhibited. pp: protein-protein interaction; violet: miR-155 targets; wine red: miR-210 targets; cyan: miR-130a targets. Targets of miR-323, miR-326 and miR-329 are not involved, because of non-interaction among their few validated target genes. C. Correlation analysis of the interested miRNAs with confirmed GBM gene markers.

Genomic analysis of human GBM showed that EGFR, PTEN and IDH1 were among the most altered genes [23], which are used as commonly monitored markers [24]. Then we tried to determine whether the 6 miRNAs were associated with these confirmed GBM markers, and the analysis showed that miR-326 and miR-329 were negatively correlated with EGFR expression levels, while miR-155 was positively related with PTEN expression levels (Figure 4C). Meanwhile, miR-329 was also inversely associated with IDH1 expression (Figure 4C). These correlation analysis results suggested the potential application of the interested miRNAs in predicting outcome of GBM patients.

Interactions of miRNAs as prognostic and predictive indicators for survival of patients with GBM

Finally, a univariate Cox proportional hazard regression model was carried out to determine the influence of the 6 miRNAs as well as clinicopathological factors (gender, age and recurrence) on patient survival. This univariate analysis indicated that age, recurrence, expression levels of miR-323, miR-326, miR-329, miR-130a, miR-155 and miR-210 were significantly associated with survival (Table 3). Furthermore, the interactions of miR-323 and miR-130a, miR-326 and miR-155, miR-326 and miR-210, and miR-329 and miR-130 were more sensitively related with survival (Table 3). Then the multivariate Cox regression model was performed to adjust the potentially confounded effects by age, gender and recurrence. The result showed that miR-326, miR-130a, miR-155, miR-210 and the four miRNA-interactions were found to be still significantly associated with survivals, whereas miR-323 and miR-329 had marginal impacts on survival (Table 3). These results indicated that miR-326, miR-130a, miR-155, miR-210 and the 4 miRNA-interactions could serve as prognostic and predictive indicators for GBM patients, which were independent of clinical variables.

Table 3. Cox regression analysis of GBM patients in relation to clinicopathological factors and miRNA expression

Discussion

In this study, we identified 38 differentially expressed miRNAs from the most significantly altered miRNAs using data from TCGA dataset. Kaplan-Meier survival and Cox multivariate proportional hazard model confirmed that the expression of miR-326, miR-130a, miR-155 and miR-210 were correlated with OS and PFS of GBM patients and were verified for the first time as independent predictors for GBM patient survival. More importantly, interactions between miR-323/miR-326/miR-329 and miR-130a/miR-155/miR-210 were also significantly related with clinical outcome and were more sensitive to discriminate and predict survival time of patients. Moreover, interactions of miR-323 and miR-130a, miR-326 and miR-155, miR-326 and miR-210 and of miR-329 and miR-130 were also confirmed as independent prognostic indicators to clinical outcome of GBM patients. In addition, the availability and rationality of these interactions as independent prognostic and predictive indicators were supported by integrated analysis of network, biological processes and correlation analysis with confirmed GBM gene markers. Our results suggest a potential application of miRNA profiles and their interactions in development and improvement of prognostic tools and treatments.

Presently, except that presurgical prognosis relies largely on age and Karnofsky Performance Status (KPS), no convincing prognostic and predictive factors have been prevalent in clinical management of GBM patients, although several prognostic and predictive markers or models have been proposed or developed, such as MGMT promoter methylation [25], BRAF fusions and IDH1 mutations [26], subclassification based on gene expression [2], Immunohistochemical analysis [27] and Volume-Age-KPS (VAK) prognostic model related with MR-imaging [28]. Notably, due to development of molecular and gene profiles, molecular stratification for patients’ outcome are increasingly emphasized, which leads to the extensive investigation and exploration of molecular markers.

MicroRNAs, as a family of small non-coding RNAs which are negatively involved in gene regulations, have been recognized as important intervention targets and predictive tools for several diseases because of the stability and convenience of miRNA detection [29-31]. Actually, several study groups have identified a pool of miRNA signatures for clinical outcome prediction. Through screening expression profiles of 200 miRNAs from 84 astrocytoma samples, miR-106a, miR-181b and miR-21 were identified as diagnostic and prognostic markers in defining the signature of astrocytomas and predicting the post-surgical outcome [6]. In another study including 38 GBM samples, miR-21, miR-181c, miR-195, and miR-196b were associated with survival of GBM patients [32]. Using TCGA dataset with 253 individuals, 23 and 19 miRNAs were defined to be associated with OS and PFS, respectively [7]. Also, in another publication with 222 GBM samples, a risk score, formulated on the basis of expression signatures of 10 miRNAs, was associated with GBM patient survival, which was suggested to predict GBM patient survival [33]. On one hand, all these studies indicated that miRNAs were thoroughly involved in GBM biology and several miRNAs could act as predictive and classified indicators for GBM clinical outcome. On the other hand, one concern has been aroused that all the identified miRNAs were almost totally different among these publications, which may be due to different uses of approaches or pre-selections of target miRNAs and so on. In this study, through calculating, sorting and overlapping mean Z-score values in GBM samples from three separate batches, we obtained the most altered miRNAs, which ensured that these candidate miRNAs were more specific and accurate to distinguish expression differences between GBM and normal brain tissues. Herein, the candidate miRNAs in this article were more convincing and feasible for further potential application in clinical practice.

This study did not follow the conventional training and validation test analysis. However, selection bias was yielded and validation of our findings was supported though corroborations as follows. First of all, all miRNAs were selected from the top most altered and overlapped miRNAs which were sorted according to mean Z-scores originated from 3 independent batches. Then, all expression levels of interested miRNAs have been validated on miRNAMap and other independent miRNA detections, which could be considered as external validations. Furthermore, biological function of the interested miRNAs and their target genes were analyzed, which may directly reflect the roles of miRNAs in tumor progression.

Among the 6 interested miRNAs, miR-326 was reported to inhibit GBM cell growth, whereas miR-155 was shown to promote GBM proliferation [34,35], which could be explained by that the target genes of miR-326 and miR-155 were mostly related with apoptosis (Figure 4A). However, according to our knowledge, there is no study reporting the associations between OS and PFS and miR-326/miR-155, while our result for the first time showed that high level of miR-326 and low level of miR-155 were significantly associated with long OS and PFS. Likewise, we first found that low levels of miR-323 and miR-329 correlated with long OS, and high level of miR-130a and low level of miR-210 were linked with extended either OS or PFS. These survival analyses indicated that miR-326 and miR-130a functioned as tumor suppressors while the others as oncogenes. However, it should be noted that expression levels of miR-323/miR-329 were elevated in no-recurrent and longer TTP patients, which were not consistent with oncogenic roles of miR-323/miR-329. Several reasons may be responsible for this inconsistence. Initially, it has been confirmed that on average one miRNA has approximately 100 target sites, regulating a large fraction of protein-coding genes involving in several biological processes, such as cell proliferation, apoptosis, and cell motion etc. [36]. Second, putative targets of miR-323 and miR-329 incorporated a family of molecules associated with cell migration and adhesion as shown in Figure 4A. Furthermore, the recurrence of GBM is related with these migration and adhesion genes [37]. Herein, miR-323/miR-329 may be involved in migration inhibition in non-recurrent patients through elevation of their expression levels. This inconsistence also occurred to miR-130a, which was shown to not only inhibit tumor suppressor RUNX3 in hepatocellular carcinoma [38] but also suppress proto-oncogene MET in lung cancer. This may be due to extensive distribution of predictive targets of miR-130a [39]. To date, there is no functional study related with miR-323, miR-329, miR-130a and miR-210 in GBM.

The complication of biological function of these miRNAs also indicated that it may be more reasonable to study their interactions, because of the multifactorial nature of the disease, and the distinguishing feature of miRNAs that an average miRNA has approximately 100 target sites and regulates a large fraction of protein-coding genes, which form a regulatory network [36,40]. OS and PFS analysis showed that the two-miRNA interaction were more sensitive and accurate to discriminate and predict the survival time in relative to one single miRNA. For instance, the longest gap of mean survival time of OS and PFS occurred on miR-130a, with 8.1 months and 5.4 months, respectively. However, the longest gap of mean survival of OS and PFS was 20.9 months and 15.8 months related with the interaction of miR-326 and miR-130a. Moreover, this interaction effect made the stratification of patients’ survival more detailed and specific. For example, the mean OS of patients with high and low miR-130a was 26.6 months and 18.5 months respectively, whereas the corresponding survival of patients with both low miR-323 and miR-130a was 30.9months, and that with both low miR-323 and high miR-130a was 14.4 months. Therefore, the interaction analysis of miRNAs may provide new views on diagnosis and prognosis of GBM patients.

In summary, we identify miR-326, miR-130a, miR-155 and miR-210 markers related with survival of GBM. More importantly, we determine the availability and rationality of 4 miRNA-interactions as more specific and accurate prognostic and predictive indicators to clinical outcome of GBM patients, implying the application for diagnostic and prognostic tools and treatments.

Competing interests

All authors declared no conflicts of interest.

Authors’ contributions

QSW designed the study, carried out data analysis and drafted the manuscript. LS performed bioinformatics analysis. HD, FYM and TY participated in the collection and analysis of the data. PY conceived of the study and participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.

Acknowledgements

This work was supported by Grants from national natural science foundation of China (NSFC, 81272197, 30973479 and 31070953 to Ying Peng), Oversea collaboration grand of national natural science foundation of China (81228010 to Deling Yin and Ying Peng) and Science & Technology Planning Project of Guangdong Province, China (2009B060700040,2011B031800141 to Ying Peng).

References

  1. Porter KR, McCarthy BJ, Freels S, Kim Y, Davis FG: Prevalence estimates for primary brain tumors in the United States by age, gender, behavior, and histology.

    Neuro Oncol 2010, 12:520-527. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  2. Verhaak RG, Hoadley KA, Purdom E, Wang V, Qi Y, Wilkerson MD, Miller CR, Ding L, Golub T, Mesirov JP, 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  3. Noushmehr H, Weisenberger DJ, Diefes K, Phillips HS, Pujara K, Berman BP, Pan F, Pelloski CE, Sulman EP, Bhat KP, et al.: Identification of a CpG island methylator phenotype that defines a distinct subgroup of glioma.

    Cancer Cell 2010, 17:510-522. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  4. Inui M, Martello G, Piccolo S: MicroRNA control of signal transduction.

    Nat Rev Mol Cell Biol 2010, 11:252-263. PubMed Abstract | Publisher Full Text OpenURL

  5. Bader AG, Brown D, Winkler M: The promise of microRNA replacement therapy.

    Cancer Res 2010, 70:7027-7030. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. Zhi F, Chen X, Wang S, Xia X, Shi Y, Guan W, Shao N, Qu H, Yang C, Zhang Y, et al.: The use of hsa-miR-21, hsa-miR-181b and hsa-miR-106a as prognostic indicators of astrocytoma.

    Eur J Cancer 2010, 46:1640-1649. PubMed Abstract | Publisher Full Text OpenURL

  7. Delfino KR, Serao NV, Southey BR, Rodriguez-Zas SL: Therapy-, gender- and race-specific microRNA markers, target genes and networks related to glioblastoma recurrence and survival.

    Cancer Genomics Proteomics 2011, 8:173-183. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  8. Lu J, Getz G, Miska EA, Alvarez-Saavedra E, Lamb J, Peck D, Sweet-Cordero A, Ebert BL, Mak RH, Ferrando AA, et al.: MicroRNA expression profiles classify human cancers.

    Nature 2005, 435:834-838. PubMed Abstract | Publisher Full Text OpenURL

  9. TCGA: Comprehensive genomic characterization defines human glioblastoma genes and core pathways.

    Nature 2008, 455:1061-1068. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Gabriely G, Yi M, Narayan RS, Niers JM, Wurdinger T, Imitola J, Ligon KL, Kesari S, Esau C, Stephens RM, et al.: Human glioma growth is controlled by microRNA-10b.

    Cancer Res 2011, 71:3563-3572. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  11. Consortium TGO: Gene ontology: tool for the unification of biology.

    Nat Genet 2000, 25:25-29. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Kanehisa M, Goto S: KEGG: kyoto encyclopedia of genes and genomes.

    Nucleic Acids Res 2000, 28:27-30. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Friedman RC, Farh KK, Burge CB, Bartel DP: Most mammalian mRNAs are conserved targets of microRNAs.

    Genome Res 2009, 19:92-105. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. da Huang W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources.

    Nat Protoc 2009, 4:44-57. PubMed Abstract | Publisher Full Text OpenURL

  15. Martin A, Ochagavia ME, Rabasa LC, Miranda J, Fernandez-de-Cossio J, Bringas R: BisoGenet: a new tool for gene network building, visualization and analysis.

    BMC Bioinformatics 2010, 11:91. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  16. Vlachos IS, Kostoulas N, Vergoulis T, Georgakilas G, Reczko M, Maragkakis M, Paraskevopoulou MD, Prionidis K, Dalamagas T, Hatzigeorgiou AG: DIANA miRPath v. 2.0: investigating the combinatorial effect of microRNAs in pathways.

    Nucleic Acids Res 2012, 40:W498-504. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Xiao F, Zuo Z, Cai G, Kang S, Gao X, Li T: miRecords: an integrated resource for microRNA-target interactions.

    Nucleic Acids Res 2009, 37:D105-110. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Landgraf P, Rusu M, Sheridan R, Sewer A, Iovino N, Aravin A, Pfeffer S, Rice A, Kamphorst AO, Landthaler M, et al.: A mammalian microRNA expression atlas based on small RNA library sequencing.

    Cell 2007, 129:1401-1414. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  19. Silber J, Lim DA, Petritsch C, Persson AI, Maunakea AK, Yu M, Vandenberg SR, Ginzinger DG, James CD, Costello JF, et al.: miR-124 and miR-137 inhibit proliferation of glioblastoma multiforme cells and induce differentiation of brain tumor stem cells.

    BMC Med 2008, 6:14. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  20. Huse JT, Brennan C, Hambardzumyan D, Wee B, Pena J, Rouhanifard SH, Sohn-Lee C, le Sage C, Agami R, Tuschl T, Holland EC: The PTEN-regulating microRNA miR-26a is amplified in high-grade glioma and facilitates gliomagenesis in vivo.

    Genes Dev 2009, 23:1327-1337. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  21. Malzkorn B, Wolter M, Liesenberg F, Grzendowski M, Stuhler K, Meyer HE, Reifenberger G: Identification and functional characterization of microRNAs involved in the malignant progression of gliomas.

    Brain Pathol 2009, 20:539-550. PubMed Abstract | Publisher Full Text OpenURL

  22. Zhou X, Ren Y, Moore L, Mei M, You Y, Xu P, Wang B, Wang G, Jia Z, Pu P, et al.: Downregulation of miR-21 inhibits EGFR pathway and suppresses the growth of human glioblastoma cells independent of PTEN status.

    Lab Invest 2010, 90:144-155. PubMed Abstract | Publisher Full Text OpenURL

  23. Parsons DW, Jones S, Zhang X, Lin JC, Leary RJ, Angenendt P, Mankoo P, Carter H, Siu IM, Gallia GL, et al.: An integrated genomic analysis of human glioblastoma multiforme.

    Science 2008, 321:1807-1812. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Camara-Quintana JQ, Nitta RT, Li G: Pathology: commonly monitored glioblastoma markers: EFGR, EGFRvIII, PTEN, and MGMT.

    Neurosurg Clin N Am 2012, 23:237-246. PubMed Abstract | Publisher Full Text OpenURL

  25. Bady P, Sciuscio D, Diserens AC, Bloch J, van den Bent MJ, Marosi C, Dietrich PY, Weller M, Mariani L, Heppner FL, et al.: MGMT methylation analysis of glioblastoma on the Infinium methylation BeadChip identifies two distinct CpG regions associated with gene silencing and outcome, yielding a prediction model for comparisons across datasets, tumor grades, and CIMP-status.

    Acta Neuropathol 2012. OpenURL

  26. von Deimling A, Korshunov A, Hartmann C: The next generation of glioma biomarkers: MGMT methylation, BRAF fusions and IDH1 mutations.

    Brain Pathol 2011, 21:74-87. PubMed Abstract | Publisher Full Text OpenURL

  27. Motomura K, Natsume A, Watanabe R, Ito I, Kato Y, Momota H, Nishikawa R, Mishima K, Nakasu Y, Abe T, et al.: Immunohistochemical analysis-based proteomic subclassification of newly diagnosed glioblastomas.

    Cancer Sci 2012. OpenURL

  28. Zinn PO, Sathyan P, Mahajan B, Bruyere J, Hegi M, Majumder S, Colen RR: A Novel Volume-Age-KPS (VAK) Glioblastoma Classification Identifies a Prognostic Cognate microRNA-Gene Signature.

    PLoS One 2012, 7:e41522. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Kong YW, Ferland-McCollough D, Jackson TJ, Bushell M: microRNAs in cancer management.

    Lancet Oncol 2012, 13:e249-258. PubMed Abstract | Publisher Full Text OpenURL

  30. Osman A: MicroRNAs in health and disease–basic science and clinical applications.

    Clin Lab 2012, 58:393-402. PubMed Abstract OpenURL

  31. Corsini LR, Bronte G, Terrasi M, Amodeo V, Fanale D, Fiorentino E, Cicero G, Bazan V, Russo A: The role of microRNAs in cancer: diagnostic and prognostic biomarkers and targets of therapies.

    Expert Opin Ther Targets 2012, 16(Suppl 2):S103-109. PubMed Abstract | Publisher Full Text OpenURL

  32. Lakomy R, Sana J, Hankeova S, Fadrus P, Kren L, Lzicarova E, Svoboda M, Dolezelova H, Smrcka M, Vyzula R, et al.: MiR-195, miR-196b, miR-181c, miR-21 expression levels and O-6-methylguanine-DNA methyltransferase methylation status are associated with clinical outcome in glioblastoma patients.

    Cancer Sci 2011, 102:2186-2190. PubMed Abstract | Publisher Full Text OpenURL

  33. Srinivasan S, Patric IR, Somasundaram K: A ten-microRNA expression signature predicts survival in glioblastoma.

    PLoS One 2011, 6:e17438. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Kefas B, Comeau L, Floyd DH, Seleverstov O, Godlewski J, Schmittgen T, Jiang J, diPierro CG, Li Y, Chiocca EA, et al.: The neuronal microRNA miR-326 acts in a feedback loop with notch and has therapeutic potential against brain tumors.

    J Neurosci 2009, 29:15161-15168. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  35. Poltronieri P, D'Urso PI, Mezzolla V, D'Urso OF: Potential of anti-cancer therapy based on anti-miR-155 oligonucleotides in glioma and brain tumours.

    Chem Biol Drug Des 2012. OpenURL

  36. Brennecke J, Stark A, Russell RB, Cohen SM: Principles of microRNA-target recognition.

    PLoS Biol 2005, 3:e85. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  37. Monticone M, Daga A, Candiani S, Romeo F, Mirisola V, Viaggi S, Melloni I, Pedemonte S, Zona G, Giaretti W, et al.: Identification of a novel set of genes reflecting different in vivo invasive patterns of human GBM cells.

    BMC Cancer 2012, 12:358. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  38. Xu N, Shen C, Luo Y, Xia L, Xue F, Xia Q, Zhang J: Upregulated miR-130a increases drug resistance by regulating RUNX3 and Wnt signaling in cisplatin-treated HCC cell.

    Biochem Biophys Res Commun 2012, 425:468-472. PubMed Abstract | Publisher Full Text OpenURL

  39. Acunzo M, Visone R, Romano G, Veronese A, Lovat F, Palmieri D, Bottoni A, Garofalo M, Gasparini P, Condorelli G, et al.: miR-130a targets MET and induces TRAIL-sensitivity in NSCLC by downregulating miR-221 and 222.

    Oncogene 2012, 31:634-642. PubMed Abstract | Publisher Full Text OpenURL

  40. Ebert MS, Sharp PA: Roles for microRNAs in conferring robustness to biological processes.

    Cell 2012, 149:515-524. PubMed Abstract | Publisher Full Text OpenURL