Combining Immune-Related Genes For Delineating the Extracellular Matrix and Predicting Hormone Therapy and Neoadjuvant Chemotherapy Benefits In Breast Cancer

Breast cancer (BC) is the most prevalent cancer in women worldwide. A systematic approach to BC treatment, comprising adjuvant and neoadjuvant chemotherapy (NAC), as well as hormone therapy, forms the foundation of the disease’s therapeutic strategy. The extracellular matrix (ECM) is a dynamic network that exerts a robust biological effect on the tumor microenvironment (TME), and it is highly regulated by several immunological components, such as chemokines and cytokines. It has been established that the ECM promotes the development of an immunosuppressive TME. Therefore, while analyzing the ECM of BC, immune-related genes must be considered. In this study, we used bioinformatic approaches to identify the most valuable ECM-related immune genes. We used weighted gene co-expression network analysis to identify the immune-related genes that potentially regulate the ECM and then combined them with the original ECM-related gene set for further analysis. Least absolute shrinkage and selection operator (LASSO) regression and SurvivalRandomForest were used to narrow our ECM-related gene list and establish an ECM index (ECMI) to better delineate the ECM signature. We stratified BC patients into ECMI high and low groups and evaluated their clinical, biological, and genomic characteristics. We found that the ECMI is highly correlated with long-term BC survival. In terms of the biological process, this index is positively associated with the cell cycle, DNA damage repair, and homologous recombination but negatively with processes involved in angiogenesis and epithelial–mesenchymal transition. Furthermore, the tumor mutational burden, copy number variation, and DNA methylation levels were found to be related to the ECMI. In the Metabric cohort, we demonstrated that hormone therapy is more effective in patients with a low ECMI. Additionally, differentially expressed genes from the ECM-related gene list were extracted from patients with a pathologic complete response (pCR) to NAC and with residual disease (RD) to construct a neural network model for predicting the chance of achieving pCR individually. Finally, we performed qRT-PCR to validate our findings and demonstrate the important role of the gene OGN in predicting the pCR rate. In conclusion, delineation of the ECM signature with immune-related genes is anticipated to aid in the prediction of the prognosis of patients with BC and the benefits of hormone therapy and NAC in BC patients.