Fig. 2

Patients with LAR TNBC are subject to epigenetic regulation. A Heatmap showing the top 30 differentially expressed epigenetic-related genes (ERGs) filtered by P  <  0.05 and log₂FC > 1. B–C Enrichment of ERGs in signaling pathways as determined by GO (B) and KEGG analysis (C). D–E Construction of ERG WGCNA network. The demonstration of module feature vector clustering (D). Correlation analysis between each module and TNBC subtype feature (E). Color labels are exclusively employed to differentiate between various gene modules and hold no intrinsic significance. F Functional annotation of hub genes in the LAR subtype processed by GO enrichment analysis. G PPI network showing coexpressed proteins of WCGNA hub genes. H ConsensusCluster determined by ERGs in the FUSCC TNBC cohort. I The Sankey diagram displaying relationships among TNBC subtypes and gene clusters. J Chord plot demonstrating the top 10 enriched signaling pathways. K Copy number alterations (left) and transcriptional expression levels (right) of representative HDAC genes. L Univariate Cox regression analysis of ERGs in the LAR subtype of TNBC. Abbreviations: LAR, luminal androgen receptor; TNBC: Triple-Negative Breast Cancer; WGCNA: Weighted Gene Co-expression Network Analysis; ERGs: Epigenetic-Related Genes; HDAC: Histone Deacetylases. Statistical analysis was performed by the Student’s t test (E, L, K) or Mann–Whitney test (B, C, J). ****Indicates P < 0.0001; ***indicates P < 0.001; **indicates P < 0.01; *indicates P < 0.05; ns indicates no significance