Identification of prognostic markers for hepatocellular carcinoma based on miRNA expression profiles

Abstract Aims The aim of the study was to identify key miRNAs related to hepatocellular carcinoma (HCC) and then to explore their potential function and clinical significance. Materials and methods The miRNA expression profiles of 387 HCC and 62 normal liver tissues were obtained from the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) databases. GEO2R tool and edgeR package in R/Bioconductor were used to screen out HCC-related miRNAs. VennDiagram package was used to identify key miRNAs related to HCC. The miRWalk tool and multiple R packages, such as pROC and survival, were used to explore potential function and clinical significance of these key miRNAs. Key findings A total of 17 and 300 HCC-related human miRNAs were identified in GEO dataset and TCGA, respectively. Thereinto seven miRNAs including hsa-miR-199a-3p, hsa-miR-199b-3p, hsa-miR-139-5p, hsa-miR-139-3p, hsa-miR-424-3p, hsa-miR-1269b and hsa-miR-1269a were key miRNAs related to HCC. Functional enrichment analysis showed that these key miRNAs were involved in multiple biological processes, such as telomere maintenance via telomerase, protein sumoylation, histone mRNA metabolic process and angiotensin maturation. Cox regression analysis indicated that hsa-miR-139-5p expression was associated with the prognosis of HCC patients. ROC curve analysis suggested that survival prediction model developed based on tumor stage and hsa-miR-139-5p exhibited good performance in predicting 3-year overall survival of HCC patients. Significance The present study identified several HCC-related miRNAs, which might serve as new diagnostic markers and therapeutic targets for HCC. In addition, hsa-miR-139-5p might act as a promising prognostic indicator for HCC patients.