Abstract
Due to the high heterogeneity of ovarian cancer (OC), it occupies the main cause of cancer-related death among women. As the most aggressive and frequent subtype of OC, high-grade serous cancer (HGSC) represents around 70 % of all patients. With the booming progress of single-cell RNA sequencing (scRNA-seq), unique and subtle changes among different cell states have been identified including novel risk genes and pathways. Here, our present study aims to identify differentially correlated core genes between normal and tumor status through HGSC scRNA-seq data analysis. R package high-dimension Weighted Gene Co-expression Network Analysis (hdWGCNA) was implemented for building gene interaction networks based on HGSC scRNA-seq data. DiffCorr was integrated for identifying differentially correlated genes between tumor and their adjacent normal counterparts. Software Cytoscape was implemented for constructing and visualizing biological networks. Real-time qPCR (RT-qPCR) was utilized to confirm expression pattern of new genes. We introduced ScHGSC-IGDC (Identifying Genes with Differential Correlations of HGSC based on scRNA-seq analysis), an in silico framework for identifying core genes in the development of HGSC. We detected thirty-four modules in the network. Scores of new genes with opposite correlations with others such as NDUFS5, TMSB4X, SERPINE2 and ITPR2 were identified. Further survival and literature validation emphasized their great values in the HGSC management. Meanwhile, RT-qPCR verified expression pattern of NDUFS5, TMSB4X, SERPINE2 and ITPR2 in human OC cell lines and tissues. Our research offered novel perspectives on the gene modulatory mechanisms from single cell resolution, guiding network based algorithms in cancer etiology field.