The analysis of gene co-expression network and immune infiltration revealed biomarkers between triple-negative and non-triple negative breast cancer

Triple-negative breast cancer (TNBC) is a heterogeneous disease with a worse prognosis. Despite ongoing efforts, existing therapeutic approaches show limited success in improving early recurrence and survival outcomes for TNBC patients. Therefore, there is an urgent need to discover novel and targeted therapeutic strategies, particularly those focusing on the immune infiltrate in TNBC, to enhance diagnosis and prognosis for affected individuals. The gene co-expression network and gene ontology analyses were used to identify the differential modules and their functions based on the GEO dataset of GSE76275. The Weighted Gene Co-Expression Network Analysis (WGCNA) was used to describe the correlation patterns among genes across multiple samples. Subsequently, we identified key genes in TNBC by assessing genes with an absolute correlation coefficient greater than 0.80 within the eigengene of the enriched module that were significantly associated with breast cancer subtypes. The diagnostic potential of these key genes was evaluated using receiver operating characteristic (ROC) curve analysis with three-fold cross-validation. Furthermore, to gain insights into the prognostic implications of these key genes, we performed relapse-free survival (RFS) analysis using the Kaplan-Meier plotter online tool. CIBERSORT analysis was used to characterize the composition of immune cells within complex tissues based on gene expression data, typically derived from bulk RNA sequencing or microarray datasets. Therefore, we explored the immune microenvironment differences between TNBC and non-TNBC by leveraging the CIBERSORT algorithm. This enabled us to estimate the immune cell compositions in the breast cancer tissue of the two subtypes. Lastly, we identified key transcription factors involved in macrophage infiltration and polarization in breast cancer using transcription factor enrichment analysis integrated with orthogonal omics. The gene co-expression network and gene ontology analyses revealed 19 modules identified using the dataset GSE76275. Of these, modules 5, 11, and 12 showed significant differences between in breast cancer tissue between TNBC and non-TNBC. Notably, module 11 showed significant enrichment in the WNT signaling pathway, while module 12 demonstrated enrichment in lipid/fatty acid metabolism pathways. Subsequently, we identified SHC4/KCNK5 and ABCC11/ABCA12 as key genes in module 11 and module 12, respectively. These key genes proved to be crucial in accu