Figure 4. (A) The correlations of the core genes of each key functional modules to each other and the correlations of the core genes of each key functional modules to cytotoxic resistance genes (HMGB1, docetaxel resistance: SKP2, AXL, KDM5D, MDH2, PIM1, SPHK1, LDHA, SOX2, Hsa-mir-143, Hsa-mir-193a, Hsa-mir-195a, Hsa-mir-204, Hsa-mir-216b, Hsa-mir-323a, Hsa-323b, Hsa-mir-34a, Hsa-mir-375, platinum resistance: Hsa-mir-205, HOTAIR, NEAT1), endocrine therapy resistance genes (AR, FHL2, VAV3, LDHA, AKR1C3, KIF4A, KDM4B) reported the in literature. (B) The core genes or targets of each key functional modules in this work and cytotoxic resistance genes (HMGB1, docetaxel resistance: SKP2, AXL, MDH2, PIM1, SPHK1, LDHA, SOX2), endocrine therapy resistance genes (AR, FHL2, VAV3, LDHA, AKR1C3, KIF4A, KDM4B) reported in the literature, were taken for regulatory network analysis, the genes which regulated well from each other in network were identified as drug metabolism-related core genes for further study. (C, D) Drug metabolism-related core genes differences analysis for RACES as both single gene or total, which were shown in the hot map and box plots. (Notes: mRNA, miRNA, lncRNA, methylation expressed as the mean value of the core genes of each key functional modules in network, miRNA 1 expressed as the mean value of antineoplastic agent response related core miRNAs in network, miRNA 2 expressed as the mean value of platinum resistance related core miRNAs in network).