Comprehensive analysis of competing endogenous RNA network and 3-mRNA signature predicting survival in papillary renal cell cancer

Long non-coding RNAs (lncRNAs) can act as competing endogenous RNAs (ceRNAs) to exert significant roles in regulating the expression of mRNAs by sequestering and binding miRNAs. To elucidate the functional roles and regulatory mechanism of lncRNAs in papillary renal cell cancer (pRCC), we conducted a comprehensive analysis of ceRNA network and constructed a mRNA signature to predict prognosis of pRCC. We collected mRNAs and lncRNAs expression profiles of 289 pRCC samples and 32 normal renal tissues, and miRNA expression profiles of 292 pRCC samples and 34 normal samples from The Cancer Genome Atlas (TCGA) database. Differential expressions of RNAs were evaluated by the “edgeR” package in R. Functional enrichment analysis of DEmRNA was performed by DAVID 6.8 and KEGG, while PPI network of top 200 DEmRNAs was conducted using the STRING database. The univariate and multivariate Cox regression were conducted to figure out the candidate DEmRNAs with predictive values in prognosis. Receiver operator characteristic (ROC) curve estimation was performed to achieve the area under the curve (AUC) of the ROC curve to judge mRNA-associated prognosic model. A ceRNA network was established relying on the basis of combination of lncRNA-miRNA interactions and miRNA-mRNA interactions. A total of 1928 DEmRNAs, 981 DElncRNAs, and 52 DEmiRNAs were identified at significance level of |log(2)Fold Change |>2 and adjusted P-value < .01. A 3-mRNA signatures consisting of ERG, RRM2, and EGF was constructed to predict survival in pRCC. Moreover, a pRCC-associated ceRNA network was constructed, with 57 lncRNAs, 11 miRNAs, and 28 mRNAs. Our study illustrated the regulatory mechanism of ceRNA network in papillary renal cancer. The identified mRNA signatures could be used to predict survival of pRCC.