Skip to main content

Integrated analysis of microRNA-target interactions with clinical outcomes for cancers



Clinical statement alone is not enough to predict the progression of disease. Instead, the gene expression profiles have been widely used to forecast clinical outcomes. Many genes related to survival have been identified, and recently miRNA expression signatures predicting patient survival have been also investigated for several cancers. However, miRNAs and their target genes associated with clinical outcomes have remained largely unexplored.


Here, we demonstrate a survival analysis based on the regulatory relationships of miRNAs and their target genes. The patient survivals for the two major cancers, ovarian cancer and glioblastoma multiforme (GBM), are investigated through the integrated analysis of miRNA-mRNA interaction pairs.


We found that there is a larger survival difference between two patient groups with an inversely correlated expression profile of miRNA and mRNA. It supports the idea that signatures of miRNAs and their targets related to cancer progression can be detected via this approach.


This integrated analysis can help to discover coordinated expression signatures of miRNAs and their target mRNAs that can be employed for therapeutics in human cancers.


As patterns of gene expression correlate with disease phenotype and patient outcome, mRNA expression profiling has been used to classify disease risks as well as prediction of outcome[14]. In addition, survival analysis with gene expression profiles is beneficial for identifying new prognostic targets of diverse diseases.

Many miRNAs have been also found to be correlated with clinical outcome in specific cancer types [510]. miRNAs are a class of small and endogenous RNA molecules that regulate their target mRNAs through translational repression or mRNA degradation [11]. In tumors, many miRNAs can be aberrantly expressed, leading to potentially abnormal regulation of their target mRNAs. Although over 1,000 human miRNAs may be encoded in human genome, the potential therapeutic markers provided by miRNAs for a diverse spectrum of diseases are still unexplored.

Recently, several investigations have put emphasis on the integrated analysis of miRNAs and mRNAs in clinical outcomes [1217]. In general, there are different approaches for the joint analysis of miRNA and mRNA data. For example, several miRNAs associated with survival rate can be extracted by survival analysis and then their relationship of inverse correlation can be identified based on analyzing miRNA and mRNA expression profiles. Most approaches do not test the clinical outcome by considering miRNA expression and mRNAs expression simultaneously. At minimum, this analysis requires the size of the cohort to be large enough for statistically significant measurement outcomes and the paired samples with both mRNA and miRNA expressions should be given. The Cancer Genome Atlas (TCGA) [18] provides different types of genomic datasets and we systematically integrated multi-omics data for cancer clinical outcome prediction [19].

Here, we present the survival analysis considering the regulatory relationships of miRNAs and their target genes. We tested clinical outcomes with the patient survival information, miRNA expression profiles and gene expression profiles for ovarian cancer and glioblastoma multiforme (GBM) through the integrated analysis of miRNA-mRNA interaction pairs. We found the miRNA-mRNA pairs with an inversely correlated expression profile that have significant survival differences between two patient groups. The results presented here suggest that this analysis can help to discover expression signatures of miRNAs and their target mRNAs that can be employed for therapeutics in human cancers.


miRNA target dataset

We obtained the largest collection of human miRNA-mRNA target relationships from the TarBase 6.0 [20]. All targets in this collection were manually curated and experimentally validated. We extracted a total of 12,879 interactions consisting of 229 miRNAs and 6,699 target genes.

Gene and miRNA expression profile

The miRNA and gene expression data sets in ovarian cancer and GBM are acquired from the Cancer Genome Atlas (TCGA) [18]. Ovarian cancer and GBM data sets consist of mRNA expression profiles and miRNA expression profiles for 496 and 425 patients, respectively. We transformed the expression values into a Z-score for each miRNA or each gene. If the Z-score of miRNA (or mRNA) is greater (less) than zero, its expression level is defined to be high (low). We excluded miRNA-mRNA pairs from the expression datasets that were not listed in the TarBase in order to avoid unnecessary calculations. Expression matrices for ovarian cancer composed of 137 miRNAs × 496 patients and 5,707 target genes × 496 patients. Among all possible 776,016 miRNAs and mRNA pairs, the number of validated interactions presented in TarBase 6.0 is 10,574 composed from 137 miRNAs and 5,707 mRNAs. Expression matrices for GBM consist of 144 miRNAs × 425 patients and 6,700 target genes × 425 patients. Among all possible 964,800 miRNA-mRNA pairs, the number of validated interactions is 9,073 from 144 miRNAs and 6,700 mRNAs.

Survival analysis of miRNA and target mRNA

We performed the survival analysis for patient groups using the Kaplan-Meier estimator. Each group was constructed based on the expression profiles of a miRNA-mRNA pair. There are four possible combinations of expression profiles (Figure 1(A)). One is a group with low expression of a miRNA and high expression of a mRNA. We defined this group as LH. On the other hand, a second group is HL that is high expression of miRNA and low expression of mRNA. The other groups, HH and LL, denote both high and both low expression of miRNA and mRNA, respectively.

Figure 1

Scheme of survival analysis with miRNA and mRNA. Scheme of survival analysis with miRNA and mRNA A. The grouping method according to miRNA and mRNA expression profiles for survival analysis. The survival analysis based on the combination of expression levels of miRNA and mRNA. B. The survival curve between LH group and HL group.

We carried out the survival analysis of each miRNA-mRNA pair for six possible pairwise sets (HL:LH, HH:LL, HL:LH, HH:HL, LL:LH and LL:HL). Among these sets, our focus was on the survival test between LH and HL (Figure 1(B)) as they are inversely related to the level of expression. The group HL indicates that high expression of miRNA causes low expression of mRNA through the mechanism that a miRNA represses or cleaves its target mRNA, while the group LH states that the low expression of miRNA results in the high expression of mRNA.


We investigated if the comparison between the two groups, HL:LH, showed any distinguishing characteristics compared to the other relationships between pairs in survival analysis. Thus, we performed the survival analysis of each miRNA-mRNA for six possible pairwise sets. Here, HH, HL, LH or LL indicates a patient group classified according to low or high expression level of a miRNA and mRNA. For example, HL means a patient group with high expression of a miRNA and low expression of mRNA. We performed the log-rank test of survival scores for all possible miRNA-mRNA pairs and then extracted the validated pairs whose interactions were shown in TarBase6.0. Table 1 shows the number of significant pairs with p-value of log-rank test <0.01 for both total pairs and validated pairs in each pairwise set. As the table indicates, there was not a substantial difference between the results. For example, in ovarian cancer, although the number of significant validated pairs in HL:LH is 261, which is higher than other pairs, the background number, 22,575 is higher as well. As a result, the number of pairs associated with survival in HL:LH might be less than we expected.

Table 1 Comparison of the number of significant pairs.

Therefore, we extracted miRNA-mRNA interaction pairs that showed significant difference between the HL group and LH group in survival rate using more stringent criteria. In Figure 2(A), we plotted the distribution of p-values corresponding to survival tests for 10,574 validated pairs in the ovarian cancer dataset. Each test was conducted from two groups containing different patients based on the expression values of miRNA and mRNA (Figure 2(B)). We obtained the threshold (corresponding to p-value < 3.66e-3) by calculating the false discovery rate (0.01) from a null model with a total of 776,016 pairs including most non-target relations. A total of 163 significant pairs with validated interactions were selected and they contained 22 miRNAs and 156 mRNAs. Expression patterns of all pairs are shown in Figure 2(C). The list of all significant interacting pairs for ovarian cancer is described in Additional file 1 Supplemental Table 1 and the list for GBM is also shown in Additional file 1 Supplemental Table 2. The survival difference between HL group and LH group suggests that the selected miRNA and its target mRNA might affect the progression of cancer in a coordinated fashion.

Figure 2

Selecting significant miRNA-mRNA pairs from HL:LH (high miRNA & low mRNA expression group vs. low miRNA & high mRNA expression group) type. Selecting significant miRNA-mRNA pairs from HL:LH (high miRNA & low mRNA expression group vs. low miRNA & high mRNA expression group) type. A. Distribution of significances by Log-rank test with HL:LH type in ovarian cancer. The dashed line indicates the threshold that the false discovery ratio is 0.01. B. The scatter plot of expression values for miR-98 and PLAGL1 pair among miRNA and mRNA pairs. Each point indicates a patient with HL:LH expression pattern. C. The heatmap of average expression values of miRNA and mRNAs for all significant pairs. Each cell indicates average expression value of a miRNA or mRNA in a group (i.e. HL or LH)

Table 2 Top 20 ranked miRNA-mRNA interaction pairs (OV).

Among top 20 ranked miRNA-mRNA interaction pairs, the roles of four targets including PLAGL1, MTUS1, MEF and IFNGR1 have been reported in ovarian cancer (Table 2). miR-148a was down-regulated in ovarian cancer cell lines and might be involved in the carcinogenesis of ovarian cancer [21, 22]. Previous research reported that overexpression of miR-148b in ovarian cancer tissues was not associated with any of the pathological features of patients with ovarian cancer. It suggested that miR-148b might be involved in the early stage of ovarian carcinogenesis [23]. Although other miRNAs including miR-196a, miR-374b, miR-124, and miR-98 have not been reported for associations with ovarian cancer, they are related with the oncogenic phenotype or their expression of other cancers [2427].

Several miRNAs are dominated in significant miRNA-mRNA pairs. miR-98 and miR-148b-3p have more than 446 and 187 target genes, respectively. miR-98 or miR-148b-3p itself shows a significant survival difference between high and low expression (pvalue < 0.0077 and pvalue < 0.0014), while miR-124-3p with 979 targets genes shows a borderline significance (pvalue < 0.0051).

LOT1(PLAGL1/ZAC1) is known to possess anti-proliferative effects and is frequently silenced in ovarian cancer and breast cancer [28]. Previous studies suggest that a shortage of the PLAGL1 protein might impair its role in regulating the cell cycle and interfere with apoptosis. Consequently, cells may grow and divide too quickly in an uncontrolled manner. Mitochondrial tumor suppressor 1 (MTUS1) is a newly identified candidate tumor suppressor gene [29]. Previous studies have shown that MTUS1 expression levels are down-regulated in cancers of the colon, ovary, pancreas, head and neck, and breast cancer [30]. MEF (myeloid ELF1-like factor, also known as ELF4) is expressed in a significant proportion of ovarian carcinomas [31]. The oncogenic activity of MEF was shown by the ability of MEF to transform NIH3T3 cells, and induce the formation of tumors in nude mice. The expression level of IFNGR1 in a typical ovarian cancer population varies, with 22% of them displaying a complete loss of the IFNγ receptor [32]. Low levels of receptor expression seem to have a negative effect on survival and are unrelated to other pathologic variables. Therefore, low expression of IFNGR1 could be regarded as an independent prognostic marker in ovarian cancer. Although we have found only several previously reported functional roles related to ovarian cancer, they hint at the possibility that other target genes might be associated with ovarian cancer development and progression.

From the analysis of GBM dataset, we could find functional evidences of four targets including CDKN1A, MTAP, KIT and ATM among top 20 ranked miRNA-mRNA interaction pairs have been reported in GBM cancer (Table 3). The roles of most top-ranked miRNAs in GBM have been reported in previous literatures. miR-34a is a transcriptional target of p53 and it also suppresses brain tumor growth by targeting c-Met and Notch[33]. The significant correlation between miR-106a expression levels and overall survival was observed in a large set of FFPE GBM specimen and it could be used as an independent prognostic biomarker in GBM patients[34]. Previous study has reported that miR-221/222 cluster is involved in regulation of cell cycle progression and cell proliferation by targeting p27 and p57 [35]. Real-time PCR experiment revealed that the expression level of miR-26b was inversely correlated with the grade of glioma[36]. Thus it suggests that miR-26b may act as a tumor suppressor in GBM. miR-17 expression was higher in glioma tissue and it is significantly related to poor overall survival[37].

Table 3 Top 20 ranked miRNA-mRNA interaction pairs (GBM).

Previous quantitative real-time PCR revealed overexpression of CDKN1A in primary GBM[38]. This result suggests that CDKN1A expression is regarded as a putative marker to distinguish primary GBM from secondary GBM. The gene for methylthioadenosine phosphorylase (MTAP) is located closely to the gene CDKN2A. MTAP-deficiency in many tumors that have been most resistant to treatment occurs commonly. Especially 70% of glioblastoma lack MTAP. Amplification of KIT in 17 (4.4%) glioblastomas was reveled from screening of 390 glioblastomas[40]. A borderline positive association (p=0.0579) between KIT amplification and TP53 mutation was also observed. ATM expression is correlated with radioresistance in primary GBM cells in culture. Genes encoding components of the DNA-damage response (DDR) pathway are frequently altered in human GBM patients and the ATM/Chk2/p53 cascade suppresses GBM formation[41].

Figure 3 displays the network of miRNA-target mRNA interactions that are associated with the difference of clinical outcome of OV cancer patients. It shows that many target genes in the network, indicated by highlighted nodes, are associated with neoplasm or cancer-related conditions. Some are interaction pairs where HL groups show good outcomes while others are when LH groups have good outcome. This network implies that the increased expression of miR-148b is significantly associated with survival rate for ovarian cancer, which agrees with the previous observation [23]. In contrast, the low level of expression of miR-124 is shown to have strong association with survival rate. The complete network of interactions can be found in Additional file 1 Supplementary Figure 1.

Figure 3

Network visualization of miRNA-mRNA interaction pairs in terms of clinical outcome in ovarian cancer patients. Network visualization of miRNA-mRNA interaction pairs in terms of clinical outcome in ovarian cancer patients. Highlighted nodes indicate genes that are associated with Neoplasm/Cancer-related condition.


We have introduced the integrated analysis of survival test with datasets of miRNA and mRNA expression profiles, as well as clinical information in ovarian cancer and GBM. We have seen that the combined expression patterns between miRNAs and mRNAs can distinguish between risk groups related to co-regulation of both. In addition, we have presented supporting evidence for functional roles of miRNA and their targets in specific cancer from the literature. Our approach can be utilized to detect clinical and therapeutic miRNAs and their targets related to outcome of several cancers.


  1. 1.

    Vasselli JR, Shih JH, Iyengar SR, Maranchie J, Riss J, Worrell R, Torres-Cabala C, Tabios R, Mariotti A, Stearman R, et al: Predicting survival in patients with metastatic kidney cancer by gene-expression profiling in the primary tumor. Proc Natl Acad Sci USA. 2003, 100 (12): 6958-6963. 10.1073/pnas.1131754100.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  2. 2.

    van de Vijver MJ, He YD, van't Veer LJ, Dai H, Hart AA, Voskuil DW, Schreiber GJ, Peterse JL, Roberts C, Marton MJ, et al: A gene-expression signature as a predictor of survival in breast cancer. The New England journal of medicine. 2002, 347 (25): 1999-2009. 10.1056/NEJMoa021967.

    Article  CAS  PubMed  Google Scholar 

  3. 3.

    van 't Veer LJ, Dai H, van de Vijver MJ, He YD, Hart AA, Mao M, Peterse HL, van der Kooy K, Marton MJ, Witteveen AT, et al: Gene expression profiling predicts clinical outcome of breast cancer. Nature. 2002, 415 (6871): 530-536. 10.1038/415530a.

    Article  PubMed  Google Scholar 

  4. 4.

    West M, Blanchette C, Dressman H, Huang E, Ishida S, Spang R, Zuzan H, Olson JA, Marks JR, Nevins JR: Predicting the clinical status of human breast cancer by using gene expression profiles. Proc Natl Acad Sci USA. 2001, 98 (20): 11462-11467. 10.1073/pnas.201162998.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  5. 5.

    Lawrie CH, Chi J, Taylor S, Tramonti D, Ballabio E, Palazzo S, Saunders NJ, Pezzella F, Boultwood J, Wainscoat JS, et al: Expression of microRNAs in diffuse large B cell lymphoma is associated with immunophenotype, survival and transformation from follicular lymphoma. Journal of cellular and molecular medicine. 2009, 13 (7): 1248-1260. 10.1111/j.1582-4934.2008.00628.x.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  6. 6.

    Lawrie CH, Gal S, Dunlop HM, Pushkaran B, Liggins AP, Pulford K, Banham AH, Pezzella F, Boultwood J, Wainscoat JS, et al: Detection of elevated levels of tumour-associated microRNAs in serum of patients with diffuse large B-cell lymphoma. British journal of haematology. 2008, 141 (5): 672-675. 10.1111/j.1365-2141.2008.07077.x.

    Article  PubMed  Google Scholar 

  7. 7.

    Lawrie CH, Soneji S, Marafioti T, Cooper CD, Palazzo S, Paterson JC, Cattan H, Enver T, Mager R, Boultwood J, et al: MicroRNA expression distinguishes between germinal center B cell-like and activated B cell-like subtypes of diffuse large B cell lymphoma. International journal of cancer Journal international du cancer. 2007, 121 (5): 1156-1161. 10.1002/ijc.22800.

    Article  CAS  PubMed  Google Scholar 

  8. 8.

    Iorio MV, Ferracin M, Liu CG, Veronese A, Spizzo R, Sabbioni S, Magri E, Pedriali M, Fabbri M, Campiglio M, et al: MicroRNA gene expression deregulation in human breast cancer. Cancer research. 2005, 65 (16): 7065-7070. 10.1158/0008-5472.CAN-05-1783.

    Article  CAS  PubMed  Google Scholar 

  9. 9.

    Calin GA, Ferracin M, Cimmino A, Di Leva G, Shimizu M, Wojcik SE, Iorio MV, Visone R, Sever NI, Fabbri M, et al: A MicroRNA signature associated with prognosis and progression in chronic lymphocytic leukemia. The New England journal of medicine. 2005, 353 (17): 1793-1801. 10.1056/NEJMoa050995.

    Article  CAS  PubMed  Google Scholar 

  10. 10.

    Takamizawa J, Konishi H, Yanagisawa K, Tomida S, Osada H, Endoh H, Harano T, Yatabe Y, Nagino M, Nimura Y, et al: Reduced expression of the let-7 microRNAs in human lung cancers in association with shortened postoperative survival. Cancer research. 2004, 64 (11): 3753-3756. 10.1158/0008-5472.CAN-04-0637.

    Article  CAS  PubMed  Google Scholar 

  11. 11.

    Bartel DP, MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116 (2): 281-297. 10.1016/S0092-8674(04)00045-5.

    Article  CAS  PubMed  Google Scholar 

  12. 12.

    Volinia S, Croce CM: Prognostic microRNA/mRNA signature from the integrated analysis of patients with invasive breast cancer. Proc Natl Acad Sci USA. 2013, 110 (18): 7413-7417. 10.1073/pnas.1304977110.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  13. 13.

    Buffa FM, Camps C, Winchester L, Snell CE, Gee HE, Sheldon H, Taylor M, Harris AL, Ragoussis J: microRNA-associated progression pathways and potential therapeutic targets identified by integrated mRNA and microRNA expression profiling in breast cancer. Cancer research. 2011, 71 (17): 5635-5645. 10.1158/0008-5472.CAN-11-0489.

    Article  CAS  PubMed  Google Scholar 

  14. 14.

    Gade S, Porzelius C, Falth M, Brase JC, Wuttig D, Kuner R, Binder H, Sultmann H, Beissbarth T: Graph based fusion of miRNA and mRNA expression data improves clinical outcome prediction in prostate cancer. BMC Bioinformatics. 2011, 12: 488-10.1186/1471-2105-12-488.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  15. 15.

    Creighton CJ, Hernandez-Herrera A, Jacobsen A, Levine DA, Mankoo P, Schultz N, Du Y, Zhang Y, Larsson E, Sheridan R, et al: Integrated analyses of microRNAs demonstrate their widespread influence on gene expression in high-grade serous ovarian carcinoma. PLoS One. 2012, 7 (3): e34546-10.1371/journal.pone.0034546.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  16. 16.

    Enerly E, Steinfeld I, Kleivi K, Leivonen SK, Aure MR, Russnes HG, Ronneberg JA, Johnsen H, Navon R, Rodland E, et al: miRNA-mRNA integrated analysis reveals roles for miRNAs in primary breast tumors. PLoS One. 2011, 6 (2): e16915-10.1371/journal.pone.0016915.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  17. 17.

    Bair E, Tibshirani R: Semi-supervised methods to predict patient survival from gene expression data. PLoS biology. 2004, 2 (4): E108-10.1371/journal.pbio.0020108.

    PubMed Central  Article  PubMed  Google Scholar 

  18. 18.

    TCGA Network: Integrated genomic analyses of ovarian carcinoma. Nature. 2011, 474 (7353): 609-615. 10.1038/nature10166.

    Article  Google Scholar 

  19. 19.

    Kim D, Shin H, Song YS, Kim JH: Synergistic effect of different levels of genomic data for cancer clinical outcome prediction. J Biomed Inform. 2012, 45 (6): 1191-1198. 10.1016/j.jbi.2012.07.008.

    Article  CAS  PubMed  Google Scholar 

  20. 20.

    Vergoulis T, Vlachos IS, Alexiou P, Georgakilas G, Maragkakis M, Reczko M, Gerangelos S, Koziris N, Dalamagas T, Hatzigeorgiou AG: TarBase 6.0: capturing the exponential growth of miRNA targets with experimental support. Nucleic acids research. 2012, 40 (Database): D222-229.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  21. 21.

    Zhou X, Zhao F, Wang ZN, Song YX, Chang H, Chiang Y, Xu HM: Altered expression of miR-152 and miR-148a in ovarian cancer is related to cell proliferation. Oncology reports. 2012, 27 (2): 447-454.

    CAS  PubMed  Google Scholar 

  22. 22.

    Delfino KR, Rodriguez-Zas SL: Transcription factor-microRNA-target gene networks associated with ovarian cancer survival and recurrence. PLoS One. 2013, 8 (3): e58608-10.1371/journal.pone.0058608.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  23. 23.

    Chang H, Zhou X, Wang ZN, Song YX, Zhao F, Gao P, Chiang Y, Xu HM: Increased expression of miR-148b in ovarian carcinoma and its clinical significance. Molecular medicine reports. 2012, 5 (5): 1277-1280.

    CAS  PubMed  Google Scholar 

  24. 24.

    Furuta M, Kozaki KI, Tanaka S, Arii S, Imoto I, Inazawa J: miR-124 and miR-203 are epigenetically silenced tumor-suppressive microRNAs in hepatocellular carcinoma. Carcinogenesis. 2010, 31 (5): 766-776. 10.1093/carcin/bgp250.

    Article  CAS  PubMed  Google Scholar 

  25. 25.

    Dalmay T, Edwards DR: MicroRNAs and the hallmarks of cancer. Oncogene. 2006, 25 (46): 6170-6175. 10.1038/sj.onc.1209911.

    Article  CAS  PubMed  Google Scholar 

  26. 26.

    Lui WO, Pourmand N, Patterson BK, Fire A: Patterns of known and novel small RNAs in human cervical cancer. Cancer research. 2007, 67 (13): 6031-6043. 10.1158/0008-5472.CAN-06-0561.

    Article  CAS  PubMed  Google Scholar 

  27. 27.

    Hoffman AE, Zheng T, Yi C, Leaderer D, Weidhaas J, Slack F, Zhang Y, Paranjape T, Zhu Y: microRNA miR-196a-2 and breast cancer: a genetic and epigenetic association study and functional analysis. Cancer research. 2009, 69 (14): 5970-5977. 10.1158/0008-5472.CAN-09-0236.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  28. 28.

    Abdollahi A, Pisarcik D, Roberts D, Weinstein J, Cairns P, Hamilton TC: LOT1 (PLAGL1/ZAC1), the candidate tumor suppressor gene at chromosome 6q24-25, is epigenetically regulated in cancer. The Journal of biological chemistry. 2003, 278 (8): 6041-6049. 10.1074/jbc.M210361200.

    Article  CAS  PubMed  Google Scholar 

  29. 29.

    Xiao J, Chen JX, Zhu YP, Zhou LY, Shu QA, Chen LW: Reduced expression of MTUS1 mRNA is correlated with poor prognosis in bladder cancer. Oncology letters. 2012, 4 (1): 113-118.

    PubMed Central  CAS  PubMed  Google Scholar 

  30. 30.

    Califano D, Pignata S, Pisano C, Greggi S, Laurelli G, Losito NS, Ottaiano A, Gallipoli A, Pasquinelli R, De Simone V, et al: FEZ1/LZTS1 protein expression in ovarian cancer. Journal of cellular physiology. 2010, 222 (2): 382-386. 10.1002/jcp.21962.

    Article  CAS  PubMed  Google Scholar 

  31. 31.

    Yao JJ, Liu Y, Lacorazza HD, Soslow RA, Scandura JM, Nimer SD, Hedvat CV: Tumor promoting properties of the ETS protein MEF in ovarian cancer. Oncogene. 2007, 26 (27): 4032-4037. 10.1038/sj.onc.1210170.

    Article  CAS  PubMed  Google Scholar 

  32. 32.

    Duncan TJ, Rolland P, Deen S, Scott IV, Liu DT, Spendlove I, Durrant LG: Loss of IFN gamma receptor is an independent prognostic factor in ovarian cancer. Clinical cancer research : an official journal of the American Association for Cancer Research. 2007, 13 (14): 4139-4145. 10.1158/1078-0432.CCR-06-2833.

    Article  CAS  Google Scholar 

  33. 33.

    Li Y, Guessous F, Zhang Y, Dipierro C, Kefas B, Johnson E, Marcinkiewicz L, Jiang J, Yang Y, Schmittgen TD, et al: MicroRNA-34a inhibits glioblastoma growth by targeting multiple oncogenes. Cancer Res. 2009, 69 (19): 7569-7576. 10.1158/0008-5472.CAN-09-0529.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  34. 34.

    Zhao S, Yang G, Mu Y, Han D, Shi C, Chen X, Deng Y, Zhang D, Wang L, Liu Y, et al: MiR-106a is an independent prognostic marker in patients with glioblastoma. Neuro Oncol. 2013, 15 (6): 707-717. 10.1093/neuonc/not001.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  35. 35.

    Zhang C, Kang C, You Y, Pu P, Yang W, Zhao P, Wang G, Zhang A, Jia Z, Han L, et al: Co-suppression of miR-221/222 cluster suppresses human glioma cell growth by targeting p27kip1 in vitro and in vivo. Int J Oncol. 2009, 34 (6): 1653-1660.

    CAS  PubMed  Google Scholar 

  36. 36.

    Wu N, Zhao X, Liu M, Liu H, Yao W, Zhang Y, Cao S, Lin X: Role of microRNA-26b in glioma development and its mediated regulation on EphA2. PLoS One. 2011, 6 (1): e16264-10.1371/journal.pone.0016264.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  37. 37.

    Lu S, Wang S, Geng S, Ma S, Liang Z, Jiao B: Increased expression of microRNA-17 predicts poor prognosis in human glioma. J Biomed Biotechnol. 2012, 2012: 970761-

    PubMed Central  PubMed  Google Scholar 

  38. 38.

    Ernst A, Hofmann S, Ahmadi R, Becker N, Korshunov A, Engel F, Hartmann C, Felsberg J, Sabel M, Peterziel H, et al: Genomic and expression profiling of glioblastoma stem cell-like spheroid cultures identifies novel tumor-relevant genes associated with survival. Clin Cancer Res. 2009, 15 (21): 6541-6550. 10.1158/1078-0432.CCR-09-0695.

    Article  CAS  PubMed  Google Scholar 

  39. 39.

    Lubin M, Lubin A: Selective killing of tumors deficient in methylthioadenosine phosphorylase: a novel strategy. PLoS One. 2009, 4 (5): e5735-10.1371/journal.pone.0005735.

    PubMed Central  Article  PubMed  Google Scholar 

  40. 40.

    Nobusawa S, Stawski R, Kim YH, Nakazato Y, Ohgaki H: Amplification of the PDGFRA, KIT and KDR genes in glioblastoma: a population-based study. Neuropathology. 2011, 31 (6): 583-588. 10.1111/j.1440-1789.2011.01204.x.

    Article  PubMed  Google Scholar 

  41. 41.

    Tribius S, Pidel A, Casper D: ATM protein expression correlates with radioresistance in primary glioblastoma cells in culture. Int J Radiat Oncol Biol Phys. 2001, 50 (2): 511-523. 10.1016/S0360-3016(01)01489-4.

    Article  CAS  PubMed  Google Scholar 

Download references


This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIP) (No. 2010-0028631); Ministry of Health and Welfare (HI13C2164).


Funding for open access charge: National Research Foundation of Korea (NRF).

This article has been published as part of BMC Medical Genomics Volume 7 Supplement 1, 2014: Selected articles from the 3rd Translational Bioinformatics Conference (TBC/ISCB-Asia 2013). The full contents of the supplement are available online at

Author information



Corresponding author

Correspondence to Ju Han Kim.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

JGJ proposed the idea, organized overall procedure and carried out the analysis. DK built the dataset of computational experiments and developed the concept. SYL interpreted the results. HJK improved the manuscript writing and reviewed the manuscript. JHK provided intellectual guidance and mentorship.

Electronic supplementary material


Additional file 1: Survival analysis of six different test set and network of significant miRNA-mRNA interaction pairs. Supplemental Table 1. Survival analysis of six different test sets or validated miRNA-mRNA interactions with ovarian cancer dataset. Supplemental able 2. Survival analysis of six different test sets for validated miRNA-mRNA interactions with GBM dataset. Supplemental Figure 1. Network visualization of all significant iRNA-mRNA interaction pairs in terms of clinical outcome in ovarian cancer patients. (PDF 455 KB)

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit

The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Joung, JG., Kim, D., Lee, S.Y. et al. Integrated analysis of microRNA-target interactions with clinical outcomes for cancers. BMC Med Genomics 7, S10 (2014).

Download citation


  • Survival Analysis
  • miRNA-mRNA interaction
  • miRNA Expression
  • Gene Expression Profile
  • Cancer Genes