Coexpression Network Analysis of Genes Related to the Characteristics of Tumor Stemness in Triple-Negative Breast Cancer

Abstract

Cancer stem cells (CSCs) are subsets of cells with the ability of self-renewal and differentiation in neoplasm, which are considered to be related to tumor heterogeneity. It has been reported that CSCs act on tumorigenesis and tumor biology of triple-negative breast cancer (TNBC). However, the key genes that cause TNBC showing stem cell characteristics are still unclear. We combined the RNA sequencing (RNA-seq) data from The Cancer Genome Atlas (TCGA) database and mRNA expression-based stemness index (mRNAsi) to further analyze mRNAsi with regard to molecular subtypes, tumor depth, and pathological staging characteristics of breast cancer (BC). Secondly, we extract the differential gene expression of tumor vs. normal group and TNBC vs. other subtypes of BC group, respectively, and intersect them to achieve precise results. We used a weighted gene coexpression network analysis (WGCNA) to screen significant gene modules and the functions of selected genes including BIRC5, CDC25A, KIF18B, KIF2C, ORC1, RAD54L, and TPX2 were carried out through gene ontology (GO) functional annotation. The Oncomine, bc-GenExMiner v4.4, GeneMANIA, Kaplan-Meier Plotter (KM-plotter), and GEPIA were used to verify the expression level and functions of key genes. In this study, we found that TNBC had the highest stem cell characteristics in BC compared with other subtypes. The lower the mRNAsi score, the better the overall survival and treatment outcome. Seven key genes of TNBC were screened and functional annotation indicated that there were strong correlations between them, relating to nuclear division, organelle fission, mitotic nuclear division, and other events that determine cell fate. Among these genes, we found four genes that were highly associated with adverse survival events. Seven key genes identified in this study were found to be closely related to the maintenance of TNBC stemness, and the overexpression of four showed earlier recurrence. The overall survival (OS) curves of all key genes between differential expression level crossed at around nine-year follow-up, which was consistent with the trend of the OS curve related to mRNAsi. These findings may provide new ideas for screening therapeutic targets in order to depress TNBC stemness.