Identification of Novel circRNA-Based ceRNA Network Involved in the Pathogenesis of Gastric Cancer

Abstract

Objective. Evidence increasingly shows that circular RNAs (circRNAs) are closely associated with tumorigenesis and cancer progression. However, the roles of circRNAs and the underlying mechanism behind these circRNAs in gastric cancer (GC) remain to be elucidated. This study is aimed at conferring a better understanding of GC pathogenesis with a specific focus on circRNA-based ceRNA action. Methods. circRNA expression profiles were downloaded from two Gene Expression Omnibus (GEO) microarray datasets, GSE152309 and GSE121445. Expression profiles of miRNAs and mRNAs were collected from The Cancer Genome Atlas (TCGA) database. The ceRNA network was constructed based on circRNA-miRNA pairs and miRNA-mRNA pairs. Functional and pathway enrichment analyses were performed to evaluate functional pathways of differentially expressed mRNAs (DEmRNAs). The PPI network was constructed by mapping DEmRNAs into the STRING database to identify hub genes, and then the DEcircRNA-DEmiRNA-hub gene subnetwork was constructed. The expression levels of candidate differentially expressed circRNAs (DEcircRNAs) in cancerous and matched noncancerous gastric tissues surgically resected from 52 GC patients were determined by the RT-qPCR analysis. Results. Differential expression analysis with Venn diagram analysis showed 11 overlapped DEcircRNAs (7 upregulated circRNAs and 4 downregulated ones) between cancerous tissues and noncancerous gastric tissues. The DEcircRNA-DEmiRNA-DEmRNA network was constructed, consisting of 2 DEcircRNAs, 7 DEmiRNAs, and 104 DEmRNAs. GO and KEGG pathway analyses revealed that 104 DEmRNAs might be associated with GC development and progression. The PPI network was constructed, yielding 16 hub genes, TFDP1, KRAS, LMNB1, MET, MYBL2, CDC25A, E2F5, HMGA1, HMGA2, CBFB, CBX3, CDC7, IGF2BP3, KIF11, PDGFB, and SMC1A, which were all upregulated in GC tissues compared with adjacent tumor-free gastric tissues. Based on the above hub genes in GC, the DEcircRNA-DEmiRNA-hub gene subnetwork was reconstructed based on hsa_circ_0000384 and hsa_circ_0000043, including 22 pairs of the upcircRNA-downmiRNA-upmRNA network. The expression levels of hsa_circ_0000384 and hsa_circ_0000043 were remarkably higher in GC tissues than those in matched adjacent tumor-free gastric tissues ( $p<0.001$ ), concurring with the bioinformatics results. Conclusion. Our study offers a better understanding of circRNA-related ceRNA regulatory mechanism in the pathogenesis of GC, highlighting two ceRNA networks based on hsa_circ_0000384 and hsa_circ_0000043.