Identification of key genes and long non-coding RNA associated ceRNA networks in hepatocellular carcinoma

Abstract

Background Hepatocellular carcinoma (HCC) is one of the leading causes of cancer-related deaths worldwide. Although multiple efforts have been made to understand the development of HCC, morbidity, and mortality rates remain high. In this study, we aimed to discover the mRNAs and long non-coding RNAs (lncRNAs) that contribute to the progression of HCC. We constructed a lncRNA-related competitive endogenous RNA (ceRNA) network to elucidate the molecular regulatory mechanism underlying HCC. Methods A microarray dataset (GSE54238) containing information about both mRNAs and lncRNAs was downloaded from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) and lncRNAs (DElncRNAs) in tumor tissues and non-cancerous tissues were identified using the limma package of the R software. The miRNAs that are targeted by DElncRNAs were predicted using miRcode, while the target mRNAs of miRNAs were retrieved from miRDB, miRTarBas, and TargetScan. Functional annotation and pathway enrichment of DEGs were performed using the EnrichNet website. We constructed a protein–protein interaction (PPI) network of DEGs using STRING, and identified the hub genes using Cytoscape. Survival analysis of the hub genes and DElncRNAs was performed using the gene expression profiling interactive analysis database. The expression of molecules with prognostic values was validated on the UALCAN database. The hepatic expression of hub genes was examined using the Human Protein Atlas. The hub genes and DElncRNAs with prognostic values as well as the predictive miRNAs were selected to construct the ceRNA networks. Results We found that 10 hub genes (KPNA2, MCM7, CKS2, KIF23, HMGB2, ZWINT, E2F1, MCM4, H2AFX, and EZH2) and four lncRNAs (FAM182B, SNHG6, SNHG1, and SNHG3) with prognostic values were overexpressed in the hepatic tumor samples. We also constructed a network containing 10 lncRNA–miRNA–mRNA pathways, which might be responsible for regulating the biological mechanisms underlying HCC. Conclusion We found that the 10 significantly overexpressed hub genes and four lncRNAs were negatively correlated with the prognosis of HCC. Further, we suggest that lncRNA SNHG1 and the SNHG3-related ceRNAs can be potential research targets for exploring the molecular mechanisms of HCC.