Construction of circRNA-Based ceRNA Network to Reveal the Role of circRNAs in the Progression and Prognosis of Hepatocellular Carcinoma

Front Genet. 2021 Feb 26:12:626764. doi: 10.3389/fgene.2021.626764. eCollection 2021.

Abstract

Background: Circular RNAs (circRNAs) are now under hot discussion as novel promising biomarkers for patients with hepatocellular carcinoma (HCC). The purpose of our study is to identify several competing endogenous RNA (ceRNA) networks related to the prognosis and progression of HCC and to further investigate the mechanism of their influence on tumor progression.

Methods: First, we obtained gene expression data related to liver cancer from The Cancer Genome Atlas (TCGA) database (http://www.portal.gdc.cancer.gov/), including microRNA (miRNA) sequence, RNA sequence, and clinical information. A co-expression network was constructed through the Weighted Correlation Network Analysis (WGCNA) software package in R software. The differentially expressed messenger RNAs (DEmRNAs) in the key module were analyzed with the Database for Annotation Visualization and Integrated Discovery (DAVID) (https://david.ncifcrf.gov/summary.jsp) to perform functional enrichment analysis including Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO). The data of miRNA expression and clinical information downloaded from TCGA were utilized for survival analysis to detach the prognostic value of the DEmiRNAs of the key module.

Results: The 201 differentially expressed miRNAs (DEmiRNAs) and 3,783 DEmRNAs were preliminarily identified through differential expression analysis. The co-expression networks of DEmiRNAs and DEmRNAs were constructed with WGCNA. Further analysis confirmed four miRNAs in the most significant module (blue module) were associated with the overall survival (OS) of patients with liver cancer, including hsa-miR-92b-3p, hsa-miR-122-3p, hsa-miR-139-5p, and hsa-miR-7850-5p. DAVID was used for functional enrichment analysis of 286 co-expressed mRNAs. The GO analysis results showed that the top enriched GO terms were oxidation-reduction process, extracellular exosome, and iron ion binding. In KEGG pathway analysis, the top three enriched terms included metabolic pathways, fatty acid degradation, and valine, leucine, and isoleucine degradation. In addition, we intersected the miRNA-mRNA interaction prediction results with the differentially expressed and prognostic mRNAs. We found that hsa-miR-92b-3p can be related to CPEB3 and ACADL. By overlapping the data of predicted circRNAs by circBank and differentially expressed circRNAs of GSE94508, we screened has_circ_0077210 as the upstream regulatory molecule of hsa-miR-92b-3p. Hsa_circ_0077210/hsa-miR-92b-3p/cytoplasmic polyadenylation element binding protein-3 (CPEB3) and acyl-Coenzyme A dehydrogenase, long chain (ACADL) were validated in HCC tissue.

Conclusion: Our research provides a mechanistic elucidation of the unknown ceRNA regulatory network in HCC. Hsa_circ_0077210 might serve a momentous therapeutic role to restrain the occurrence and development of HCC.

Keywords: bio-markers; ceRNA; circRNAs; hepatocellular carcinoma; prognosis.