Construction of a ceRNA network to reveal a vascular invasion associated prognostic model in hepatocellular carcinoma

The aim of this study is to explore the prognostic value of vascular invasion (VI) in hepatocellular carcinoma (HCC) by searching for competing endogenous RNAs (ceRNA) network and constructing a new prognostic model for HCC. The differentially expressed genes (DEGs) between HCC and normal tissues were identified from GEO and TCGA. StarBase and miRanda prediction tools were applied to construct a circRNA-miRNA-mRNA network. The DEGs between HCC with and without VI were also identified. Then, the hub genes were screened to build a prognostic risk score model through the method of least absolute shrinkage and selection operator. The prognostic ability of the model was assessed using the Kaplan-Meier method and Cox regression analysis. In result, there were 221 up-regulated and 47 down-regulated differentially expressed circRNAs (DEcircRNAs) in HCC compared with normal tissue. A circRNA-related ceRNA network was established, containing 11 DEcircRNAs, 12 DEmiRNAs, and 161 DEmRNAs. Meanwhile, another DEG analysis revealed 625 up-regulated and 123 down-regulated DEGs between HCC with and without VI, and then a protein-protein interaction (PPI) network was built based on 122 VI-related DEGs. From the intersection of DEGs within the PPI and ceRNA networks, we obtained seven hub genes to build a novel prognostic risk score model. HCC patients with high-risk scores had shorter survival time and presented more advanced T/N/M stages as well as VI occurrence. In conclusion a novel prognostic model based on seven VI-associated DEGs within a circRNA-related ceRNA network was constructed in this study, with great ability to predict the outcome of HCC patients.