Bioinformatics-Based Identification of a circRNA-miRNA-mRNA Axis in Esophageal Squamous Cell Carcinomas

Esophageal squamous cell carcinoma (ESCC) has a poor prognosis due to the lack of early disease symptoms. Using bioinformatics tools, this study aimed to discover differentially expressed nonprotein-coding RNAs and genes with potential prognostic relevance in ESCC.

Two microRNAs (miRNAs) and one circular RNA (circRNA) microarray datasets were downloaded from the Gene Expression Omnibus (GEO) database. Differential expression of miRNAs (DEMs) and circRNAs (DECs) was, respectively, identified in ESCC tissue and compared to adjacent healthy tissue. Further analysis was performed using the miRNA microarray datasets, where miRTarBase was used to predict which messenger RNAs (mRNAs) was present. This was followed by protein-protein interaction (PPI) network, Kyoto Encyclopedia of Genes and Genomes (KEGG), and Gene Ontology (GO) analyses. Moreover, cytoHubba and UALCAN were used to predict the important nodes and perform patient survival analysis, respectively. The miRNA-associated circRNAs were predicted using the ENCORI website. The interaction between DECs and the predicted circRNAs was also determined. A circRNA-miRNA-mRNA axis was constructed.

Associated with RAP1B and circ_0052867, two miRNAs (miR-133b and miR-139-5p) were identified as being differentially expressed and downregulated across the two datasets. Finally, the circ_0052867/miR-139-5p/RAP1B regulatory axis was established.

This study provides support for the possible mechanisms of disease progression in ESCC.