Identification of candidate biomarkers and analysis of prognostic values in ovarian cancer by integrated bioinformatics analysis

Ovarian cancer is the first leading cause of mortality in gynecological malignancies. To identify key genes and microRNAs in ovarian cancer, mRNA microarray dataset GSE36668, GSE18520, GSE14407 and microRNA dataset GSE47841 were downloaded from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) and microRNAs (DEMs) were obtained using GEO2R. Functional and pathway enrichment analysis were performed for DEGs using DAVID database. Protein-protein interaction (PPI) network was established by STRING and visualized by Cytoscape. Following, overall survival (OS) analysis of hub genes was performed by the Kaplan-Meier plotter online tool. Module analysis of the PPI network was performed using MCODE. Moreover, miRecords was applied to predict the targets of the DEMs. A total of 345 DEGs were obtained, which were mainly enriched in the terms related to cell cycle, mitosis, and ovulation cycle process. A PPI network was constructed, consisting of 141 nodes and 296 edges. Sixteen genes had high degrees in the network. High expression of four genes of the 16 genes was associated with worse OS of patients with ovarian cancer, including CCNB1, CENPF, KIF11, and ZWINT. A significant module was detected from the PPI network. The enriched functions and pathways included cell cycle, nuclear division, and oocyte meiosis. Additionally, a total of 36 DEMs were identified. The expression of KIF11 was negatively correlated with that of has-miR-424 and has-miR-381, and it was also the potential target of two microRNAs. In conclusion, these results identified key genes, which could provide potential targets for ovarian cancer diagnosis and treatment.