Identification of serum prognostic marker miRNAs and construction of microRNA-mRNA networks of esophageal cancer.

Abstract

Esophageal cancer is a common tumor of the digestive system with poor prognosis. This study was to gain a better understanding of the mechanisms involved in esophageal cancer and to identify new prognostic markers. We downloaded the esophageal cancer miRNA expression profile microarray data (GSE113740, GSE112264, GSE122497, GSE113486, and GSE106817) from the GEO database, extracted the esophageal cancer miRNA sequencing data from The Cancer Genome Atlas (TCGA) database, and then used a bioinformatics approach to select common differentially expressed miRNAs (DEMs). Differentially expressed genes (DEGs) were selected by predicting DEM target genes using the miRWalk database and intersecting with differential genes obtained from TCGA database for esophageal cancer. The STRING database was used to obtain protein-protein interaction (PPI) relationships to construct the DEM-DEG network. Furthermore, we selected core genes and core miRNAs associated with esophageal cancer prognosis by performing survival and univariate/multivariate COX analysis on DEMs and DEGs in the network and performed GSEA analysis on core genes alone, and finally the expression of the markers was verified by qPCR in esophageal cancer cell lines Eca109, SKGT-4 and normal esophageal epithelial cells HEEC. Nine DEMs were obtained, of which three were upregulated and six were downregulated, and 326 DEGs were obtained, of which 105 were upregulated and 221 were downregulated. Survival univariate/multivariate COX analysis revealed that five genes, ZBTB16, AQP4, ADCYAP1R1, PDGFD, and VIPR2, and two microRNAs, miR-99a-5p, and miR-508-5p, were related to esophageal cancer prognosis. GSEA analysis showed that the following genes may be involved in esophageal cancer prognosis: ZBTB16 may through the MTOR signaling pathway, AQP4 through the GNRH signaling pathway, ADCYAP1R1 through the PPAR signaling pathway, VIPR2 through the P53 signaling pathway and PDGFD through the PENTOSE-PHOSPHATE signaling pathway.