Mining TCGA database for gene expression in ovarian serous cystadenocarcinoma microenvironment

Abstract

Background Ovarian cancer is one of the leading causes of female deaths worldwide. Ovarian serous cystadenocarcinoma occupies about 90% of it. Effective and accurate biomarkers for diagnosis, outcome prediction and personalized treatment are needed urgently Methods Gene expression profile for OSC patients was obtained from the TCGA database. The ESTIMATE algorithm was used to calculate immune scores and stromal scores of expression data of ovarian serous cystadenocarcinoma samples. Survival results between high and low groups of immune and stromal score were compared and differentially expressed genes (DEGs) were screened out by limma package. The Gene Ontology (GO), the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis and the protein-protein interaction (PPI) network analysis were performed with the g:Profiler database, the Cytoscape and Search Tool for the Retrieval of Interacting Genes (STRING-DB). Survival results between high and low immune and stromal score groups were compared. Kaplan-Meier plots based on TCGA follow up information were generated to evaluate patients’ overall survival. Results Eighty-six upregulated DEGs and one downregulated DEG were identified. Three modules, which included 49 nodes were chosen as important networks. Seven DEGs (VSIG4, TGFBI, DCN, F13A1, ALOX5AP, GPX3, SFRP4) were considered to be correlated with poor overall survival. Conclusion Seven DEGs (VSIG4, TGFBI, DCN, F13A1, ALOX5AP, GPX3, SFRP4) were correlated with poor overall survival in our study. This new set of genes can become strong predictor of survival, individually or combined. Further investigation of these genes is needed to validate the conclusion to provide novel understanding of tumor microenvironment with ovarian serous cystadenocarcinoma prognosis and treatment.