Skip to main content
Download PDF

Two precision medicine predictive tools for six malignant solid tumors: from gene-based research to clinical application

Journal of Translational Medicine volume 17, Article number: 405 (2019) Cite this article

Abstract

Background

The current study aimed to construct competing endogenous RNA (ceRNA) regulation network and develop two precision medicine predictive tools for colorectal cancer (CRC).

Methods

Differentially expressed (DE) analyses were performed between CRC tissues and normal tissues. A ceRNA regulation network was constructed based on DElncRNAs, DEmiRNAs, and DEmRNAs.

Results

Fifteen mRNAs (ENDOU, MFN2, FASLG, SHOC2, VEGFA, ZFPM2, HOXC6, KLK10, DDIT4, LPGAT1, BEX4, DENND5B, PHF20L1, HSP90B1, and PSPC1) were identified as prognostic biomarkers for CRC by multivariate Cox regression. Then a Fifteen-mRNA signature was developed to predict overall survival for CRC patients. Concordance indexes were 0.817, 0.838, and 0.825 for 1-, 2- and 3-year overall survival. Patients with high risk scores have worse OS compared with patients with low risk scores.

Conclusion

The current study provided deeper understanding of prognosis-related ceRNA regulatory network for CRC. Two precision medicine predictive tools named Smart Cancer Survival Predictive System and Gene Survival Analysis Screen System were constructed for CRC. These two precision medicine predictive tools can provide valuable precious individual mortality risk prediction before surgery and improve the individualized treatment decision-making.

Introduction

Colorectal cancer (CRC) is one of the most prevalent malignant tumors worldwide. There were around 1.1 million newly diagnosed CRC patients and 0.55 million CRC patients died in 2018 [1]. The 5-year overall survival (OS) of CRC patients was less than 12% [2, 3]. Due to the poor overall survival, early detection of CRC patients with high mortality risk has important significance for improving the individualized treatment decision-making. Several prognostic models have been developed for CRC patients [4, 5]. However, the computational formulas of these prognostic models were too complex for clinical application by patients without calculation tools and medical knowledge. Additional, these prognostic models provided overall survival prediction for different groups, but not individual mortality risk prediction.

From the perspective of precise medicine, a good prognostic model should be able to provide individual mortality risk prediction for specific patient at the individual level. Considering the clinical need of precious individual mortality risk prediction for patients with different tumors, our research team has developed several precision medicine predictive tools for gastric cancer [6] and hepatocellular carcinoma [7]. For individual mortality risk prediction, our precision medicine predictive tools have the following advantages: full-time individual risk prediction, visual illustration, numerical presentation, customizable subgroups, and on-line computing.

So far the molecular biological regulatory mechanism of development and prognosis of tumor remains unclear. Several researches have explored underlying molecular biological regulatory mechanism for different tumors [8,9,10,11]. Salmena and collaborators presented an interesting molecular regulatory mechanism named competing endogenous RNAs (ceRNAs) [12]. The lncRNAs can indirectly regulate the expression of mRNAs through binding shared miRNA response elements [13]. Several researches have explored potential molecular regulatory mechanism for CRC [14,15,16]. Therefore the current research aimed to depict prognosis related ceRNA regulatory network and develop two precision medicine predictive tools for CRC patients.

Materials and methods

Study cohort

RNA sequencing and miRNA sequencing data were obtained from TCGA database. RNA sequencing data contained 14,449 lncRNAs and 20,337 mRNAs whereas miRNA sequencing data contained 1881 miRNAs. Four hundred and twenty-eight CRC patients were included after removing patients without complete survival information. GSE17538 dataset were downloaded from GEO database. GSE17538 dataset involved 231 CRC patients and 23,328 RNAs (GPL570 platform). The original read count values in TCGA dataset were normalized by log2 transformation. The gene background file (Gencode.v29 supplied by The European Bioinformatics Institute of The European Molecular Biology Laboratory (EMBL-EBI) database) was used for gene symbol annotation.

Differentially expressed analyses and regulatory network

The original RNA data were processed by “edgeR” package, with a defined P of 0.05 and a ratio of 1.5 times between tumor and non-tumor tissues [17]. The original miRNA sequencing data were processed by “limma” package [18]. First, the interaction associations between lncRNAs and miRNAs were identified in miRcode database [19]. Second, miRTarBase [20], miRDB [21], and TargetScan [22] were searched for miRNA-targeted mRNAs. The ceRNA network was visualized by Cytoscape v3.6.1 [23].

Statistical analyses and artificial intelligence algorithms

Random survival forest, Multi-task logistic regression, and Cox survival regression algorithms were carried out according to the algorithms suggested in the original articles [24,25,26,27,28,29]. Statistical analyses were carried out through SPSS Statistics 19.0 (SPSS Inc.,USA). Other analyses were carried out by R version 3.5.2 with corresponding packages. P value < 0.05 was defined as statistically significant.

Results

Baseline characteristics

TCGA dataset contained 428 CRC patients and GSE17538 dataset contained 231 CRC patients. The clinical information of included patients was shown in Table 1. The mortality of GSE17538 dataset was 40.3% (93/231), which was significantly higher than 22.9% (98/428) of TCGA dataset. There were significant differences in terms of survival time and pathological stage, whereas there were no significant differences in terms of age and gender.

Table 1 Clinical features of colorectal cancer patients
Full size table

Differentially expressed analyses

Differential expression analysis could identify genes with significant differences in expression levels between normal samples and tumor samples. Compared with normal tissues, 3005 lncRNAs (2224 up-regulated and 781 down-regulated), 332 miRNAs (246 up-regulated and 86 down-regulated), and 6713 mRNAs (4087 up-regulated and 2626 down-regulated) were identified in CRC tissues. The volcano plots for differentially expressed RNAs were shown in Additional file 1: Figure S1.

Screening of prognostic mRNAs

Univariate Cox regression was used to explore potential prognostic biomarkers for CRC. Out of previous differentially expressed mRNAs, there were 2504 mRNAs identified as prognostic biomarkers for CRC. Out of 2504 potential prognostic biomarkers, there were 1371 risk factors and 1133 protective factors.

Development of ceRNA network

The miRNA- targeted mRNAs that could be searched in three above databases were defined as the miRNA-targeted mRNAs. Third, these miRNA-targeted mRNAs were intersected with previous prognostic mRNAs for development of ceRNA network. Finally, the ceRNA network, consisting of 14 lncRNAs, 29 miRNAs, and 79 mRNAs, were constructed for CRC. The ceRNA network was visualized in Fig. 1 by Cytoscape v3.6.1. This ceRNA network depicted potential regulatory relations among lncRNAs, miRNAs, and mRNAs, and was helpful to understand the potential mechanisms of tumor prognosis.

Fig. 1
figure 1

Competitive endogenous RNA network chart

Full size image

Functional enrichment analyses

Functional enrichment analyses were performed based on previous prognostic mRNAs in ceRNA network and identified 26 enriched Gene Ontology (GO) terms. The top 15 enriched GO terms were shown in Fig. 2a. The prognostic mRNAs were mainly enriched in transcription factor activity, RNA polymerase II proximal promoter sequence-specific DNA binding, regulation of MAP kinase activity, regulation of protein serine/threonine kinase activity, inactivation of MAPK activity, MAP kinase tyrosine/serine, MAP kinase phosphatase activity, protein tyrosine/serine, E-box binding, protein tyrosine phosphatase activity, and negative regulation of MAP kinase activity. The top KEGG pathways (Fig. 2b) were mainly enriched in MAPK signaling pathways, Proteoglycans in cancer, Breast cancer, IL-17 signaling pathway, TNF signaling pathway, Osteaclast differentiation, MicorRNAs in cancer, Human immunodeficiency virus 1 infection, Fluid shear stress and atherosclesosis, Ras signaling pathway, Hepatitis B, Regulation of action cytoskeleton, PI3K-Akt signaling pathway. These biological functions were helpful to understand the roles of these genes in tumor prognosis. The further circular visualization chart was presented in Fig. 3.

Fig. 2
figure 2

Chord chart of prognostic mRNAs: a GO terms; b KEGG pathways

Full size image
Fig. 3
figure 3

Circular chart of prognostic mRNAs

Full size image

Construction of predictive model

Multivariate Cox regression analyses were used to screen independent risk factors for tumor prognosis. Fifteen mRNAs were used to develop a predictive model for CRC (Table 2). The formula for predictive model was as follows: The risk score = (− 4.273 * ENDOU) + (− 1.489 * MFN2) + (− 1.243 * FASLG) + (− 0.904 * SHOC2) + (− 0.834 * VEGFA) + (− 0.690 * ZFPM2) + (0.249 * HOXC6) + (0.446 * KLK10) + (0.672 * DDIT4) + (0.705 * LPGAT1) + (0.711 * BEX4) + (1.038 * DENND5B) + (1.065 * PHF20L1) + (1.093 * HSP90B1) + (1.146 * PSPC1). The prognostic nomogram chart was shown in Fig. 4.

Table 2 Model information of prognostic mRNAs
Full size table
Fig. 4
figure 4

Prognostic nomogram for overall survival in colorectal cancer

Full size image

Survival curve analyses

Survival curve analyses were performed to explore the survival influences of included mRNAs (Additional file 1: Figure S2). Kaplan–Meier curves indicated that patients with high expression had significantly poor overall survival than that with low expression for these 15 mRNAs (P < 0.05). Comparison of Kaplan–Meier curves supported that these 15 genes were associated with overall survival in CRC patients.

Predictive performance in model dataset

CRC patients were divided into high risk subgroup and low risk subgroup according to median risk score. Figure 5a demonstrated that there was significant difference between two subgroups for OS (P < 0.001). Concordance indexes were 0.817, 0.838, and 0.825 respectively for 1-year, 2-year, and 3-year OS (Fig. 5b). Calibration curves for OS were presented in Fig. 6, indicating a good agreement between predicted mortality and actual mortality for 1-year, 2-year, and 3-year OS.

Fig. 5
figure 5

Clinical performance in model cohort: a survival curve plot; b time-dependent receiver operating characteristic curve plot

Full size image
Fig. 6
figure 6

Calibration curves for model cohort. a Calibration curve for 1-year survival; b calibration curve for 2-year survival; c calibration curve for 3-year survival

Full size image

Predictive performance in validation dataset

Kaplan–Meier plot (Fig. 7a) demonstrated that OS in high risk subgroup was significantly worse than that in low risk subgroup (P < 0.05). Concordance indexes were 0.773, 0.824, and 0.801 for 1-year, 2-year, and 3-year OS (Fig. 7b). Calibration curves for OS were depicted in Fig. 8, demonstrating that the predicted mortality was in good agreement with the actual mortality.

Fig. 7
figure 7

Clinical performance in validation cohort: a survival curve plot; b time-dependent receiver operating characteristic curve plot

Full size image
Fig. 8
figure 8

Calibration curves for validation cohort. a Calibration curve for 1-year survival; b calibration curve for 2-year survival; c calibration curve for 3-year survival

Full size image

Independence assessment of predictive model

In model cohort, prognostic signature and pathological stage were independent influence factors for OS (Table 3). In validation cohort, prognostic signature, pathological stage, and age were identified as independent influence factors for OS. Decision curves and clinical impact curve for OS were presented in Additional file 1: Figure S3, suggesting that the clinical efficacy of the current prognostic model was superior to pathological stage and age.

Table 3 Independence assessment of prognostic model
Full size table

Smart Cancer Survival Predictive System

We developed a precision medicine predictive tool named Smart Cancer Survival Predictive System, providing a novel convenient on-line calculator for prediction of OS. Smart Cancer Survival Predictive System (Fig. 9) is available at: https://zhangzhiqiao6.shinyapps.io/Smart_cancer_predictive_system_11_CRC_B1003/. Smart Cancer Survival Predictive System provided three individual mortality risk curves predicted by Random survival forest algorithm, Multitask logistic regression algorithm, and Cox survival regression algorithm according to the calculation formula in ogininal articles.

Fig. 9
figure 9

Home page of Smart Cancer Survival Predictive System. a Survival curves display page; b quantitative display page

Full size image

Gene Survival Analysis Screen System

To further explore survival curves of previous prognostic genes in different gender and pathological stage subgroups, we developed a new online program named Gene Survival Analysis Screen System (Fig. 10). Gene Survival Analysis Screen System is available at: https://zhangzhiqiao5.shinyapps.io/Gene_Survival_Subgroup_Analysis_B1003/. Gene Survival Analysis Screen System provided seven tumor datasets for exploration research and allowed users to select tumor type, gender, and stage, which were the important influence factors to the prognosis.

Fig. 10
figure 10

Home page of Gene Survival Analysis Screen System. a Survival curves display page; b Cox survival analysis display page

Full size image

Clinical application in other tumors

The current study download five different tumor datasets from TCGA database as external validation datasets to explore the clinical application of the current prognostic model. Figure 11 displayed the diagnostic accuracy of the current prognostic model in five tumors, including hepatocellular carcinoma (n = 348), breast cancer (n = 1030), gastric cancer (n = 265), lung cancer (n = 494), and ovarian cancer (n = 370).

Fig. 11
figure 11

Time-dependent receiver operating characteristic curve plots. a Hepatocellular carcinoma; b breast cancer; c gastric cancer; d lung cancer; e ovarian cancer

Full size image

Predictive performance in five malignant solid tumors

To further explore the predictive performance of the current prognostic model, we externally validated the current prognostic model in a super merge dataset (n = 2507), which including five malignant solid tumors (hepatocellular carcinoma, breast cancer, gastric cancer, lung cancer, and ovarian cancer). Kaplan–Meier plot (Fig. 12a) demonstrated that OS in high risk subgroup was significantly worse than that in low risk subgroup (P < 0.05). Concordance indexes were 0.663, 0.639, and 0.655 for 1-year, 2-year, and 3-year OS (Fig. 12b). Calibration curves for OS were depicted in Fig. 13, demonstrating that the predicted mortality was in good agreement with the actual mortality.

Fig. 12
figure 12

Clinical performance in five malignant solid tumor datasets: a survival curves for high risk group and low risk group. b Time-dependent receiver operating characteristic curves

Full size image
Fig. 13
figure 13

Calibration curves in five malignant solid tumor datasets: a calibration curve for 1-year survival; b calibration curve for 2-year survival; c calibration curve for 3-year survival

Full size image

Discussion

The current research constructed a prognosis-related ceRNA network for CRC. A prognostic nomogram was developed for OS based on prognostic mRNAs in ceRNA network. The current prognostic model could discriminate high risk patients from low risk patients for different time points. Meanwhile, two novel precision medicine predictive tools were developed to provide convenient on-line calculation for prediction of OS in CRC patients. Smart Cancer Survival Predictive System provided three individual mortality risk curves predicted by Random survival forest algorithm, Multitask logistic regression algorithm, and Cox survival regression. Gene Survival Analysis Screen System provided survival curve comparison and multivariable survival analysis results. Smart Cancer Survival Predictive System and Gene Survival Analysis Screen System could provide precious individualized mortality risk prediction for CRC patients.

The current study screened differentially expressed RNAs between CRC tissues and normal tissues and then constructed a prognosis-related ceRNA regulatory network for CRC. Based on mRNAs in ceRNA regulatory network, the current study further identified independent prognostic biomarkers for overall survival in CRC patients. Through ceRNA regulatory network, the regulatory relationship among lncRNA-miRNA-mRNA was depicted, providing post-transcriptional biological regulatory pathway information for CRC. GO term and KEGG pathway analyses were helpful to further understand the biological function and molecular regulatory pathway of prognostic mRNAs in ceRNA network. The fifteen prognostic mRNAs were identified as independent prognostic biomarkers for overall survival in CRC patients by multivariable Cox regression. The ceRNA regulatory network and functional enrichment analyses provided potential biological prognostic indicators and therapeutic targets for future researches.

Mitofusin 2 (MFN2) promotes cell proliferation and invasiveness in gastric cancer [30]. Lin et al. reported that miRNA-30d can regulate cancer angiogenesis and cancer proliferation through MYPT1/c-JUN/VEGFA pathway [31]. Liu et al. reported that ZFPM2 Antisense RNA 1 (ZFPM2-AS1) promotes tumor proliferation and inhibits apoptosis through regulating miR-137 in renal cell cancer [32]. Li et al. reported that Homeobox C6 (HOXC6) promotes invasion via EMT pathway in hepatocellular carcinoma [33]. Petraki et al. reported that Kallikrein Related Peptidase 10 (KLK10) expression is associated with overall survival in CRC [34]. Xu et al. reported that MiR-199b-5p can promote tumor proliferation through regulating KLK10 in cervical cancer [35]. Liu et al. reported that Heat Shock Protein 90 Beta Family Member 1 (HSP90B1) is significantly related with worse overall survival in lung cancer [36]. Cawthorn et al. reported that up-regulated HSP90B1 is related with distant metastasis [37]. Kessler et al. reported that Paraspeckle Component 1 (PSPC1) is significantly related with poor prognosis in hepatocellular carcinoma [38]. The results in current study further supported the credibility of these findings above in previous researches.

High quality samples and standardized RNA extraction technology are helpful to improve the quality and reliability of research data. According to the Standard Operating Procedure suggested by TCGA database, for gene sequencing sample, conventional methods of fixation and heating of biopsies may lead to inactivation of antigens and genetic material in tissues, therefore frozen sections was recommended in TCGA database. Samples should be snap-frozen and stored in cryovials in liquid nitrogen vapor and should not allow tissues thaw until sectioned. The mirVana isolation technology uses organic extraction and then fixes RNA on glass fiber filter to purify total RNA, so it can prepare high-purity and high-quality small RNA molecules within 2 h. RNA analytes should not undergo multiple freeze–thaw cycles to maintain biological activity. RNA is easily degraded by ribonuclease, so special methods should be taken to prevent RNA degradation. All reagents must be made of ribonuclease free materials and all products and disposable materials used must be free of ribonuclease. In order to prevent contamination by ribonuclease carried by the skin, gloves must be worn before handling biological samples. RNA analytes should be placed on wet ice and RNA quantification should be performed before freezing.

Advantages of the current study: Firstly, our study team developed a precision medicine tool named Smart Cancer Survival Predictive System, providing full-time individual mortality risk prediction with visual illustration and numerical presentation. Secondly, a novel precision medicine tool named Gene Survival Analysis Screen System was developed to explore the associations between prognostic genes (including clinical parameters) and overall survival. Users are free to choose the appropriate subgroup, gender, and stage. Meanwhile, users can upload their own dataset to explore and validate the research result.

Shortcomings of the current research: first, detection platforms were different between model dataset and validation dataset. The influence of different detection platforms on gene expression read counts should be taken into account while evaluating clinical application of the current prognostic model. Second, several important prognostic factors were not included in the current study, including surgical, radiotherapy and chemotherapy regimens. Third, the prognostic model was developed and validated by using datasets downloaded from public databases without research team’s study cohort. External applicability of prognostic model needs to be validated by study cohorts from different population.

Conclusions

In conclusion, the current study provided deeper understanding of prognosis-related ceRNA regulatory network for CRC. Two precision medicine predictive tools named Smart Cancer Survival Predictive System and Gene Survival Analysis Screen System were constructed for CRC. These two precision medicine predictive tools can provide individual prognostic information before surgery and improve the individualized treatment decision-making.

Availability of data and materials

The study data is available at: https://zhangzhiqiao5.shinyapps.io/Gene_Survival_Subgroup_Analysis_B1003/.

Abbreviations

CRC:

colorectal cancer

TCGA:

The Cancer Genome Atlas

GEO:

The Gene Expression Omnibus

ROC:

receiver operating characteristic

OS:

overall survival

lncRNA:

long non-coding RNA

miRNA:

microRNAs

HR:

hazard ratio

CI:

confidence interval

AJCC:

The American Joint Committee on Cancer

SD:

standard deviation

DCA:

decision curve analysis

ceRNA:

competitive endogenous RNA

References

  1. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.

    Article  PubMed  Google Scholar 

  2. Kalyan A, Kircher S, Shah H, Mulcahy M, Benson A. Updates on immunotherapy for colorectal cancer. J Gastrointest Oncol. 2018;9(1):160–9.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Qi P, Xu MD, Ni SJ, Shen XH, Wei P, Huang D, Tan C, Sheng WQ, Zhou XY, Du X. Down-regulation of ncRAN, a long non-coding RNA, contributes to colorectal cancer cell migration and invasion and predicts poor overall survival for colorectal cancer patients. Mol Carcinog. 2015;54(9):742–50.

    Article  CAS  PubMed  Google Scholar 

  4. Wang X, Zhou J, Xu M, Yan Y, Huang L, Kuang Y, Liu Y, Li P, Zheng W, Liu H, et al. A 15-lncRNA signature predicts survival and functions as a ceRNA in patients with colorectal cancer. Cancer Manag Res. 2018;10:5799–806.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Zhang H, Wang Z, Wu J, Ma R, Feng J. Long noncoding RNAs predict the survival of patients with colorectal cancer as revealed by constructing an endogenous RNA network using bioinformation analysis. Cancer Med. 2019;8(3):863–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Cheng C, Wang Q, Zhu M, Liu K, Zhang Z. Integrated analysis reveals potential long non-coding RNA biomarkers and their potential biological functions for disease free survival in gastric cancer patients. Cancer Cell Int. 2019;19:123.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Zhang Z, Li J, He T, Ouyang Y, Huang Y, Liu Q, Wang P, Ding J. The ceRNA regulatory network reveals potential prognostic biomarkers for overall survival in hepatocellular carcinoma. Cancer Sci. 2019;110:2905–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Zeng J, Cai X, Hao X, Huang F, He Z, Sun H, Lu Y, Lei J, Zeng W, Liu Y, et al. LncRNA FUNDC2P4 down-regulation promotes epithelial–mesenchymal transition by reducing E-cadherin expression in residual hepatocellular carcinoma after insufficient radiofrequency ablation. Int J Hyperthermia. 2018;34(6):802–11.

    Article  CAS  PubMed  Google Scholar 

  9. Zhong X, Long Z, Wu S, Xiao M, Hu W. LncRNA-SNHG7 regulates proliferation, apoptosis and invasion of bladder cancer cells assurance guidelines. J BUON. 2018;23(3):776–81.

    PubMed  Google Scholar 

  10. Shi X, Zhao Y, He R, Zhou M, Pan S, Yu S, Xie Y, Li X, Wang M, Guo X, et al. Three-lncRNA signature is a potential prognostic biomarker for pancreatic adenocarcinoma. Oncotarget. 2018;9(36):24248–59.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Huang Y, Xiang B, Liu Y, Wang Y, Kan H. LncRNA CDKN2B-AS1 promotes tumor growth and metastasis of human hepatocellular carcinoma by targeting let-7c-5p/NAP1L1 axis. Cancer Lett. 2018;437:56–66.

    Article  CAS  PubMed  Google Scholar 

  12. Salmena L, Poliseno L, Tay Y, Kats L, Pandolfi PP. A ceRNA hypothesis: the Rosetta Stone of a hidden RNA language? Cell. 2011;146(3):353–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  13. Thomson DW, Dinger ME. Endogenous microRNA sponges: evidence and controversy. Nat Rev Genet. 2016;17(5):272–83.

    Article  CAS  PubMed  Google Scholar 

  14. Yuan W, Li X, Liu L, Wei C, Sun D, Peng S, Jiang L. Comprehensive analysis of lncRNA-associated ceRNA network in colorectal cancer. Biochem Biophys Res Commun. 2019;508(2):374–9.

    Article  CAS  PubMed  Google Scholar 

  15. Liang Y, Zhang C, Ma MH, Dai DQ. Identification and prediction of novel non-coding and coding RNA-associated competing endogenous RNA networks in colorectal cancer. World J Gastroenterol. 2018;24(46):5259–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Pan H, Pan J, Song S, Ji L, Lv H, Yang Z. Identification and development of long non-coding RNA-associated regulatory network in colorectal cancer. J Cell Mol Med. 2019;23:5200–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.

    Article  CAS  PubMed  Google Scholar 

  18. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Jeggari A, Marks DS, Larsson E. miRcode: a map of putative microRNA target sites in the long non-coding transcriptome. Bioinformatics. 2012;28(15):2062–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Chou CH, Shrestha S, Yang CD, Chang NW, Lin YL, Liao KW, Huang WC, Sun TH, Tu SJ, Lee WH, et al. miRTarBase update 2018: a resource for experimentally validated microRNA-target interactions. Nucleic Acids Res. 2018;46(D1):D296–302.

    Article  CAS  PubMed  Google Scholar 

  21. Wong N, Wang X. miRDB: an online resource for microRNA target prediction and functional annotations. Nucleic Acids Res. 2015;43(Database issue):D146–52.

    Article  CAS  PubMed  Google Scholar 

  22. Agarwal V, Bell GW, Nam JW, Bartel DP. Predicting effective microRNA target sites in mammalian mRNAs. Elife. 2015. https://doi.org/10.7554/eLife.05005.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Fisher LD, Lin DY. Time-dependent covariates in the Cox proportional-hazards regression model. Ann Rev Public Health. 1999;20(undefined):145–57.

    Article  CAS  Google Scholar 

  25. Katzman JL, Shaham U, Cloninger A, Bates J, Jiang T, Kluger Y. DeepSurv: personalized treatment recommender system using a Cox proportional hazards deep neural network. BMC Med Res Methodol. 2018;18(1):24.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Xu H, Gu X, Tadesse MG, Balasubramanian R. A modified random survival forests algorithm for high dimensional predictors and self-reported outcomes. J Comput Graph Stat. 2018;27(4):763–72.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Nasejje JB, Mwambi H. Application of random survival forests in understanding the determinants of under-five child mortality in Uganda in the presence of covariates that satisfy the proportional and non-proportional hazards assumption. BMC Res Notes. 2017;10(1):459.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Alaeddini A, Hong SH. A multi-way multi-task learning approach for multinomial logistic regression. An application in joint prediction of appointment miss-opportunities across multiple clinics. Methods Inf Med. 2017;56(4):294–307.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Bisaso KR, Karungi SA, Kiragga A, Mukonzo JK, Castelnuovo B. A comparative study of logistic regression based machine learning techniques for prediction of early virological suppression in antiretroviral initiating HIV patients. BMC Med Inform Decis Mak. 2018;18(1):77.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Fang CL, Sun DP, Chen HK, Lin CC, Hung ST, Uen YH, Lin KY. Overexpression of mitochondrial GTPase MFN2 represents a negative prognostic marker in human gastric cancer and its inhibition exerts anti-cancer effects. J Cancer. 2017;8(7):1153–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Lin ZY, Chen G, Zhang YQ, He HC, Liang YX, Ye JH, Liang YK, Mo RJ, Lu JM, Zhuo YJ, et al. MicroRNA-30d promotes angiogenesis and tumor growth via MYPT1/c-JUN/VEGFA pathway and predicts aggres sive outcome in prostate cancer. Mol Cancer. 2017;16(1):48.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Liu JG, Wang HB, Wan G, Yang MZ, Jiang XJ, Yang JY. Long noncoding RNA ZFPM2-AS1 promotes the tumorigenesis of renal cell cancer via targeting miR-137. Eur Rev Med Pharmacol Sci. 2019;23(13):5675–81.

    PubMed  Google Scholar 

  33. Li PD, Chen P, Peng X, Ma C, Zhang WJ, Dai XF. HOXC6 predicts invasion and poor survival in hepatocellular carcinoma by driving epithelial–mesenchymal transition. Aging. 2018;10(1):115–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Petraki C, Youssef YM, Dubinski W, Lichner Z, Scorilas A, Pasic MD, Komborozos V, Khalil B, Streutker C, Diamandis EP, et al. Evaluation and prognostic significance of human tissue kallikrein-related peptidase 10 (KLK10) in col orectal cancer. Tumour Biol. 2012;33(4):1209–14.

    Article  CAS  PubMed  Google Scholar 

  35. Xu LJ, Duan Y, Wang P, Yin HQ. MiR-199b-5p promotes tumor growth and metastasis in cervical cancer by down-regulating KLK10. Biochem Biophys Res Commun. 2018;503(2):556–63.

    Article  CAS  PubMed  Google Scholar 

  36. Liu K, Kang M, Li J, Qin W, Wang R. Prognostic value of the mRNA expression of members of the HSP90 family in non-small cell lung cancer. Exp Ther Med. 2019;17(4):2657–65.

    CAS  PubMed  PubMed Central  Google Scholar 

  37. Cawthorn TR, Moreno JC, Dharsee M, Tran-Thanh D, Ackloo S, Zhu PH, Sardana G, Chen J, Kupchak P, Jacks LM, et al. Proteomic analyses reveal high expression of decorin and endoplasmin (HSP90B1) are associated with br east cancer metastasis and decreased survival. PLoS ONE. 2012;7(2):e30992.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Kessler SM, Hosseini K, Hussein UK, Kim KM, List M, Schulthei CS, Schulz MH, Laggai S, Jang KY, Kiemer AK. Hepatocellular carcinoma and nuclear paraspeckles: induction in chemoresistance and prediction for poor survival. Cell Physiol Biochem. 2019;52(4):787–801.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We would like to express sincere thanks to Dr. Gary S. Collins (University of Oxford), Dr. Manali Rupji (Emory University), Mrs. Qingmei Liu for inspirations, suggestions and assistances in development of our precision medicine tools.

Funding

The present study was funded by Medical Science and Technology Foundation of Guangdong Province (A2016450 and B2018237).

Author information

Author notes

  1. Zhiqiao Zhang, Tingshan He and Liwen Huang were the co-first authors

Authors and Affiliations

  1. Department of Infectious Diseases, Shunde Hospital, Southern Medical University, No. 1 Jiazi Road, Lunjiao, Shunde District, Foshan, 528308, Guangdong Province, China

    Zhiqiao Zhang, Tingshan He, Liwen Huang, Yanling Ouyang, Jing Li, Yiyan Huang, Peng Wang & Jianqiang Ding

Authors
  1. Zhiqiao Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Tingshan He
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Liwen Huang
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Yanling Ouyang
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Jing Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Yiyan Huang
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Peng Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Jianqiang Ding
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

Conceptualization, methodology and resources: ZZ, PW, JD, LJ, LH, TH, YO, and YH; Investigation, data curation, formal analysis, validation, software, project administration, and supervision: ZZ, PW, JD, LJ, LH, TH, YO, and YH; Writing and visualization: ZZ, PW, and JD; Funding acquisition: ZZ. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jianqiang Ding.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

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

Supplementary information

Additional file 1: Figure S1.

Volcano plot: (A) Volcano plot of lncRNAs; (B) Volcano plot of miRNAs; (C) Volcano plot of mRNAs. Figure S2. Survival curves of prognostic genes. Figure S3. Decision curves and clinical impact curves.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), 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 (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, Z., He, T., Huang, L. et al. Two precision medicine predictive tools for six malignant solid tumors: from gene-based research to clinical application. J Transl Med 17, 405 (2019). https://doi.org/10.1186/s12967-019-02151-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12967-019-02151-8

Keywords

Journal of Translational Medicine

ISSN: 1479-5876

Contact us