Multiple Omics Integration Reveals Key Circular RNAs in Hepatocellular Carcinoma

Background HCC is one of the most common malignancies with an increasing incidence worldwide, especially in Asian countries. However, even though targeted cancer therapy drugs such as sorafenib and regorafenib are available, the overall outcome of HCC remains unsatisfactory. Thus, it is urgent to investigate the molecular mechanisms of HCC progression, so as to provide accurate diagnostic criteria and therapeutic targets. Methods RNA-seq data was used to identify and quantify circular RNAs (circRNAs). DESeq2 was used to identify the differentially expressed circRNAs. miRNA binding sites within circRNAs were identified by miRanda. Gene set enrichment analysis (GSEA) was conducted to predict the biological function of circRNAs. Results The differential expression analysis identified 107 upregulated and 95 downregulated circRNAs in HCC tissues. We observed that a differentially expressed circRNA (DE-circRNA), hsa_circ_0141900 was highly negatively correlated with its parental gene RAB1A (PCC < -0.6), which was also closely associated with mTOR signaling pathway. Moreover, we also constructed competing endogenous RNA (ceRNA) network to identify key circRNAs involved in HCC. Notably, hsa_circ_0002130 and hsa_circ_0008774 were highly correlated with the genes involved in gluconeogenesis and HNF3A pathway via the target genes, GOT2 and AR, suggesting that the two circRNAs might regulate these pathways, respectively. Survival analysis revealed that GOT2 was associated with favorable prognosis. Furthermore, high expression of hsa_circ_0002130 was found to inhibit tumor cell growth and promotes GOT2 expression. Conclusion In summary, the circRNAs highlighted by the integrative analysis greatly improved our understanding of the underlying mechanism of circRNAs in HCC.