Breast Cancer Prognosis Prediction and Immune Pathway Molecular Analysis Based on Mitochondria-Related Genes

Luo, Weixu; Han, Yuanshan; Li, Xin; Liu, Zhuo; Meng, Pan; Wang, Yuhong

doi:https://doi.org/10.1155/2022/2249909

Genetics Research

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusions Data Availability Disclosure Conflicts of Interest Authors’ Contributions Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2022 | Article ID 2249909 | https://doi.org/10.1155/2022/2249909

Breast Cancer Prognosis Prediction and Immune Pathway Molecular Analysis Based on Mitochondria-Related Genes

Weixu Luo,^1,2Yuanshan Han,³Xin Li,²Zhuo Liu,⁴Pan Meng,¹and Yuhong Wang¹

Academic Editor: Giuseppe Damante

Received09 Feb 2022

Accepted18 May 2022

Published31 May 2022

Abstract

Background. Mitochondria play an important role in breast cancer (BRCA). We aimed to build a prognostic model based on mitochondria-related genes. Method. Univariate Cox regression analysis, random forest, and the LASSO method were performed in sequence on pretreated TCGA BRCA datasets to screen out genes from a Gene Set Enrichment Analysis, Gene Ontology: biological process gene set to build a prognosis risk score model. Survival analyses and ROC curves were performed to verify the model by using the GSE103091 dataset. The BRCA datasets were equally divided into high- and low-risk score groups. Comparisons between clinical features and immune infiltration related to different risk scores and gene mutation analysis and drug sensitivity prediction were performed for different groups. Result. Four genes, MRPL36, FEZ1, BMF, and AFG1L, were screened to construct our risk score model in which the higher the risk score, the poorer the prognosis. Univariate and multivariate analyses showed that the risk score was significantly associated with age, M stage, and N stage. The gene mutation probability in the high-risk score group was significantly higher than that in the low-risk score group. Patients with higher risk scores were more likely to die. Drug sensitivity prediction in different groups indicated that PF-562271 and AS601245 might be new inhibitors of BRCA. Conclusion. We developed a new workable risk score model based on mitochondria-related genes for BRCA prognosis and identified new targets and drugs for BRCA research.

1. Introduction

Breast cancer (BRCA) is one of the most common cancers worldwide; it is also a complex disease with different types and molecular characteristics [1]. According to the expression status of hormone receptors (estrogen and progesterone receptors) and HER2, BRCAs are mainly classified into four subtypes: luminal A, luminal B, HER2-enriched, and triple-negative [2]. Although there are many different subtypes of BRCAs, it has been reported that some of these subtypes might convert to other under specific conditions [3]. Because BRCA is a complex heterogeneous disease, the pathogenesis and clinical manifestations of different patients may differ. Diagnostic testing of patients with cancer has not yet been fully integrated into the clinical practice [4]. In cases where the subtype of BRCA cannot be accurately determined, identifying prognostic risk factors may be an effective method for diagnosis and treatment [5]. Therefore, in this study, we aimed to construct a new prognostic model for BRCA.

Mitochondria are among the most important organelles. Except for highly specialized mammalian mature red blood cells, most eukaryotic cells have their own mitochondria. Oxidative phosphorylation during aerobic respiration occurs in the mitochondria; this process is the main source of cellular energy [6]. Mitochondria play an important role in BRCA. For example, metabolic patterns in cancer change according to the different needs of various solid tumors [7]. Gathering mitochondria, by altering their subcellular localization, results in an increase production of ROS, which are toxic to cells, weakening the invasive ability of BRCA cells [8]. In addition, as a cofactor of the electron transport chain in mitochondrial oxidative phosphorylation, an increase in heme synthesis could inhibit glycolysis and oxidative metabolism, reducing the proliferation of BRCA cells [9]. Besides, many mitochondria-related genes play important roles in cancer. For example, BMF is a proapoptotic gene. In BRCA cell lines MCF-7 and MDA-MB-231, BMF could bind Bcl-2, instead of Bax, leading to mitochondrial outer membrane permeabilization and, finally, to apoptosis [10]. Based on the importance of mitochondria and related genes, we wondered whether we could construct a prognostic model for BRCA using mitochondria-related genes.

The mutation risk of moderate-penetrance genes in patients with BRCA is two to four times that in the general population [11]. Identifying high-risk germline mutations and implementing strategies to reduce risk could increase the survival rate of patients with BRCA [12]. Hence, we analyzed the mutation patterns of different risk score groups. It is difficult to screen useful drugs directly on patients with BRCA because of the high cost, great difficulty, and long periods of treatment needed. Fortunately, the differences between the expression patterns of tumor and normal cells can be used to predict the sensitivity of different cells to a certain drug. Additionally, it has been reported that pRRophetic is a workable tool for predicting drug effects in BRCA; therefore, we used these tools in our study [13]. There are various cell populations in the breast ductal epithelial layer and the normal breast tissue microenvironment, including immune and stromal cells. As their interactions play a major role in early BRCA [14], we performed immune infiltration analyses of patients with different risk scores.

In this study, we used TCGA BRCA datasets to screen four mitochondria-related genes, namely, MRPL36, FEZ1, BMF, and AFG1L, and constructed and verified a risk score model to predict the prognosis of BRCA. We also compared the features of the different risk score groups, such as mutation, drug sensitivity prediction, and immune infiltration.

2. Materials and Methods

2.1. Accessing and Preprocessing of BRCA Datasets

TCGA BRCA datasets were downloaded from the UCSC Xena platform (https://xenabrowser.net/), including RNA-seq data and clinical features of the patients, such as gender, age, stage, grade, survival status, and survival time. The external verification dataset GSE103091 () was downloaded from https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE103091. All probes in the microarray datasets were renamed and re-annotated. For probes with the same gene name, the mean value of the expression levels was defined as the gene’s expression value. Data were normalized before performing any other analyses.

2.2. The Establishment of a Risk Score Model Based on Mitochondria-Related Genes

The mitochondria-related genes were collected from the Gene Set Enrichment Analysis, Gene Ontology: biological process (GSEA GO BP) gene set. Taking the survival time as one of the factors and the expression of the mitochondria-related genes as the other, univariate Cox analyses were carried out to calculate both the value and hazard ratio. The , value indicated whether the gene was significantly related to the patients’ survival time, whereas the hazard ratio indicated their relationships. A preliminary screening of the genes whose value was >5 was performed, and further screening was done using random survival forest. After the error rate stabilized, we determined the importance of the genes in the upcoming prognostic model. After 1000 cross-validations using the LASSO method, the minimum lambda value was selected to determine how many genes would be used to establish the prognostic model. Finally, a group of genes was chosen based on their LASSO regression coefficients. The risk score consisted of the sum of the gene expression values multiplied by their LASSO regression coefficients.

2.3. Relationship between Clinical Features and Patient Risk Scores, and Mutation Analysis

All patient risk scores were calculated using the formula described above. Taking the risk score as one of the factors, univariate and multivariate Cox analyses were carried out with different clinical features, such as age or stage, as the other factor, to investigate the relationships between the risk score and each of these features. Patients in the BRCA datasets were equally divided into two groups according to their risk scores: low- and high-risk score groups. In addition, the mutation information from the patients was separately displayed by group to compare the mutation rates and types between them.

2.4. Immune Infiltration Analysis of Patients with Different Risk Scores

We performed four algorithms [15], ESTIMATE, McCounter, single-sample gene set enrichment analysis (ssGSEA), and TIMER, to estimate the relationship between the risk score and cell components or immune response. Heat maps were used to show differences in immune reactions under different algorithms.

2.5. Gene Set Enrichment Analysis (GSEA) and Drug Sensitivity Prediction

The gene sets used in this study were downloaded from the MSigDB (http://www.gsea-msigdb.org/gsea/msigdb). clusterProfiler (an R package) was used to perform the GSEA. pRRophetic (an R package) was used to predict the IC50, representing drug sensitivity [16], in the low- and high-risk score groups.

2.6. Statistical Analysis

All the statistical analyses in this study were performed using R 3.6.1 (https://www.r-project.org/). The Shapiro–Wilk test was used to check the normal distribution of variables. Unpaired Student’s t-test was used to check the differences in variables that conformed to the normal distribution, whereas the Wilcoxon test was used to check the differences in variables that did not conform to the normal distribution. For multiple groups, one-way ANOVA was used to compare the mean values. The Kaplan–Meier method was used to generate survival curves using the R package “survminer.” All heat maps were generated using “pheatmap.”

3. Results

3.1. Establishment of a Risk Score Model to Predict the Prognosis of BRCA

In the GSEA GO BP gene set, there were 345 genes related to mitochondria, 63 of which were screened out by univariate Cox analysis. The hazard ratios of these 63 genes are shown in Figure 1(a). Using the random forest algorithm (Figure 1(b)) and LASSO regression analysis (Figure 1(c)), we determined that the four genes with the highest variable relative importance were MRPL36, FEZ1, BMF, and AFG1L. Therefore, the LASSO coefficient of these four genes was used as the weight of the risk score model. The risk score of the prognostic model was built using the following formula:

(a)

(b)

(c)

Figure 1

Screening of mitochondria-related genes and establishment of a risk score model. (a) Differential (increased or reduced) mitochondria-related gene expression in the hazard group compared with that in the normal group. (b) Error rates of randomly generated trees (upper panel). Variable relative importance of the four selected mitochondria-related genes. (c) LASSO analysis: partial likelihood deviance values were plotted against log (λ) (upper panel). The relative abundance of the selected genes varies with the risk score.

3.2. Prognosis Prediction of Different Risk Scores

To verify whether the risk score model worked, we evaluated the prognosis of different risk scores in the BRCA and GSE103091 datasets. We found that the increase in the risk score was proportional to the increase in patient mortality (Figures 2(a), 2(d). In addition, the survival curves shown in Figures 2(b) and 2(e) show that the higher the risk score, the worse the prognosis of the patients (). Moreover, the ROC curves in Figures 2(c) and 2(f) show that the AUC values for one and three years were higher than 0.75, indicating that the risk score model was workable.

(a)

(b)

(c)

(d)

(e)

(f)

3.3. Relationships between Risk Scores and Clinical Characteristics, and Mutations in Different Risk Score Groups

The results from univariate and multivariate analyses of risk scores with clinical characteristics of the subjects included in TCGA BRCA datasets are shown in Figure 3(a). The univariate analysis showed that the risk score was significantly associated with age, stage, T stage, M stage, and N stage (). In addition, the multivariate analysis showed that the risk score was significantly related to age, M stage, and N stage (). All hazard ratios were greater than 1.

(a)

(b)

We performed mutation analysis of the low- and high-risk score groups, and the results were sorted by mutation patterns (Figure 3(b)). Comparing the upper and lower panels of Figure 3(b), it was found that patients in the high-risk score group had more mutations than those in the low-risk score group. Comparing the gene name on the left and the mutation rate on the right of Figure 3(b), the top-ranked gene mutations in the high- and low-risk score groups were roughly the same, in which the mutation probability of TP53, PIK3CA, TTN, CDH1, and GATA3 exceeded 10%. Remarkably, the mutation probability of HMCN1 increased from 5% in the low-risk score group to 7% in the high-risk score group. Comparing only the information on the right side of the panels in Figure 3(b), the pattern of mutations changed in the different groups. In the high-risk score group, there were more nonstop mutations in PIK3CA, fewer nonsense mutations in TP53, and more splice site and frameshift deletions in TP53 than in the low-risk score group. In addition, compared to the low-risk score group, fewer missense mutations and multihit mutations in CDH1, fewer frameshift insertions in GTAT3, fewer frameshift insertions and deletions in HMCN1, and more nonsense mutations in HMCN1 occurred in the high-risk score group. In particular, the in-frame insertion mutations in MUC4 were more abundantly found in the high-risk score group than in the low-risk score group. Thus, the mutation patterns differed between the high- and low-risk score groups.

3.4. Relationships between Risk Scores and Clinical Features, and Drug Sensitivity Prediction

As shown in Figure 4(a), as the risk scores increased, the death rate of the patients also increased significantly. In addition, an increase in the risk score corresponded to significant increases in the rate of stage IV and stage M1 BRCAs and to a significant decreased in stage T1 BRCA. The risk score of the patients was not related to age or N stage.

(a)

(b)

To identify potentially effective drugs against BRCA, we used pRRophetic to predict the different clinical drug responses in the high- or low-risk score groups (Figure 4(b)). Compared with the low-risk score group, the high-risk score group seemed more sensitive to cisplatin, crizotinib (PF-02341066), CHIR99021, (-)-parthenolide, AS601245, and PF-562271, indicating that these drugs could be used to treat BRCA.

3.5. Relationship between Risk Scores and Immune Infiltration

To analyze the differences in immune infiltration among patients with different risk scores, we used four methods: ESTIMATE, McCounter, single-sample gene set enrichment analysis (ssGSEA), and the TIMER algorithm (Figure 5). The epidemic immunity in BRCA was analyzed according to different gene expression patterns. According to the ESTIMATE analysis, tumor purity was negatively correlated with the risk score. In other words, the degree of immune cell infiltration increased at higher risk score values. Stromal and immune ESTIMATE scores significantly increased with increasing risk scores. According to the McCounter analysis, with an increase in risk scores, T cells, CD8 T cells, cytotoxic lymphocytes, monocytic lineage, myeloid dendritic cells, neutrophils, endothelial cells, and fibroblasts significantly increase as well. According to the ssGSEA analysis, with an increase in risk score, activated CD4 T cells, activated dendritic cells, CD 56 bright natural killer cells, CD56 dim natural killer cells, central memory CD4 T cells, central memory CD8 T cells, effector memory CD4 T cells, effector memory CD8 T cells, gamma delta T cells, immature dendritic cells, macrophages, mast cells, myeloid-derived suppressor cells, memory B cells, monocytes, natural killer cells, natural killer T cells, neutrophils, plasmacytoid dendritic cells, regulatory T cells, T follicular helper cells, type 1 T helper cells, type 17 T helper cells, and type 2 T helper cells were significantly increased. According to the TIMER analysis, with an increase in the risk score, CD4 T cells, CD8 T cells, neutrophils, macrophages, and dendritic cells were significantly more abundant. These results were consistent and showed that the risk scores were related to immune infiltration.

3.6. GSEA

We then chose several GO processes and KEGG pathways and performed GSEA. The results showed that the NES values of DNA-dependent DNA replication, the HIPPO signaling pathway, regulation of mRNA processing, the TNF-mediated signaling pathway, cell cycle G2 M phase transition, signal transduction by p53 class mediator (Figure 6(a)), the mTOR signaling pathway, the TGFβ signaling pathway, and the VEGF signaling pathway (Figure 6(b)) were greater than 0 (), indicating that these processes or pathways were activated in the high-risk score group.

(a)

(b)

4. Discussion

In recent years, the combination of the random survival forest model and the LASSO regression method has been widely used to establish various disease prognostic models from different perspectives. Not only the gene expression patterns, imaging data, clinical characteristics [17], radiomics results [18], or other variables can also be used to establish prognostic models. Mitochondria play an important role in different cancers; mitochondria-related genes had also been used to build prognostic models of other diseases, for example, acute myeloid leukemia [19]. But this thinking perspective seems did not attract researchers’ attention in previous breast cancer studies. We filled this research gap by constructing a prognosis risk score model based on four mitochondria-related genes. Besides its role in BRCA, BMF is also an important prognostic marker for patients with other cancers, such as colon cancer [20] and hepatocellular carcinoma [21]. MRPL36 is a mitochondrial ribosomal protein that has not been studied in BRCA, but MRPL36 has been correlated with poor progression-free survival in ovarian cancer [22]. The translation product of AFG1L, LACE1, interacts with p53 and mediates its mitochondrial translocation and apoptosis [23]. FEZ1 is an anti-invasive factor [24], and it is related to mitochondrial anterograde movement [25]. There is no research on FEZ1 in BRCA, but it is known that changes in the mitochondrial subcellular position in cells could affect the invasion ability of BRCA [8]. These four genes may be new useful targets for BRCA treatment or prognosis.

It is known that, after accumulating a certain number of mutations, somatic cells will develop into cancer cells that will be out of control and proliferate abnormally without limitation. In this study, we analyzed and compared the mutation patterns of the top 30 genes in the high- and low-risk score groups. TP53 is a well-known tumor suppressor gene, and its translation product, p53, is associated with apoptosis in cancer [26]. The loss of p53 leads to systemic neutrophil inflammation, which promotes BRCA metastasis [27]. PIK3CA is an important molecule in the PI3K signaling pathway, and its translation product, p110α, is a catalytic subunit of the PI3Kα complex [2]. A PIK3CA-activating mutation often occurs in BRCA and causes disorders in the PI3K/Akt pathway. Therefore, PI3K inhibitors have been used to target PIK3CA for the treatment of different subtypes of advanced BRCA [28]. As the translation product of CDH1, E-cadherin is a structural molecule for cell adhesion and an important factor in maintaining epithelial characteristics. The regulation of cell migration is related to changes in E-cadherin subcellular localization [29, 30]. In some genetic contexts, AKT inhibitors may cause the accumulation of mutations in CDH1, thereby accelerating the metastasis and progression of BRCA [29]. Additionally, GATA3 is a transcription factor with a zinc-finger structure that regulates T-cell receptors in the immune system [31]. GATA3 is one of the most frequently mutated genes in BRCA; it can affect chromatin localization of FOXA1 and ER-α [32]. In BRCA, splice, frameshift, truncation, and extension are the most common mutations found in GATA3 [31]. The main function of GATA3 in the breast is to alter the fate of luminal cells [33]. The translation product of MUC4 is Mucin-4, a highly glycosylated protein, located on the surface of the cell membrane, that normally lubricates and protects vascular and epithelial surfaces [34]. It also promotes the combination of blood and tumor cells [35]. Therefore, Mucin-4 plays an important role in tumor progression and metastasis. In addition, an increased expression of Mucin-4 masks the HER2 epitope on the cell surface, resulting in trastuzumab resistance of cancer cells in patients with HER2-overexpressed BRCA [36]. Little research has focused on HMCN1 in breast cancer, but it is known that this gene is related to tumor heterogeneity and poor prognosis of BRCA [37]. The different mutations analyzed in our risk score model were consistent with those reported in other studies.

The assessment of tumor-infiltrating lymphocytes has been proven to be a reliable surrogate indicator and a powerful independent prognostic biomarker for immune antitumor activity in patients with BRCA [38]. Noteworthy, lung metastasis is a major cause of death in patients with BRCA. Neutrophil extracellular traps in the tumor microenvironment can promote this metastasis [39]. In BRCA, activated M2 macrophages can produce cytokines and promote tumor cell proliferation, metastasis, survival, and tumor angiogenesis [40]. Some cancer-associated fibroblast subpopulations can alter the immune microenvironment and lead to BRCA cell metastasis [41]. In some NRF2-positive breast cancers, immune cells, such as CD8+ T, CD4+ T, dendritic, and stromal cells (such as adipocytes, fibroblasts, and keratinocytes), are highly infiltrated [42]. In this study, we used several algorithms to evaluate the degree of immune infiltration. We found that, as the risk score increased, the infiltration of neutrophils and macrophages also increased. The results of the different algorithms in our study reached a consistent conclusion and coincided with those of previous reports.

Cisplatin is one of the most commonly used chemotherapeutic agents. It can block BRCA progression by blocking early epithelial-mesenchymal transition [43], preventing cancer cell division, and inducing apoptosis [44] or other mechanisms. Crizotinib (PF-02341066) is an FDA-approved small-molecule protein kinase inhibitor [45]. It can be used to treat BRCA by targeting the inhibition of MET and many other kinases [46]. CHIR99021 is an inhibitor of GSK3. Studies have shown that CHIR99021 inhibits the growth of the human BRCA cell line MDA-MB-231 [47]. (-)-Parthenolide is a sesquiterpene lactone isolated from Chrysanthemum morifolium. It can target FAK1 and affect the proliferation, survival, and movement of BRCA cells [48]. PF-562271 is a FAK inhibitor. Although it has not been reported in BRCA, it can inhibit the viability and migration ability of embryonal rhabdomyosarcoma cells [49]. Activation of JNK signaling in BRCA promotes the survival of cancer cells, and AS601245 is an inhibitor of JNK [50]. The screening done for all these drugs, according to the different risk score groups, is consistent with the results from other studies and suggests that PF-562271 and AS601245 may be used as new drugs for the treatment of BRCA.

Different biological processes and signaling pathways are involved in BRCA. For example, cancer stem-like cells in BRCA under DNA replication stress induce drug resistance and recurrence [51]. mRNA processing is usually disordered in BRCA [52]. In addition, the HIPPO signaling pathway plays an important role in BRCA metastasis through crosstalk with other signaling pathways [53]. Moreover, TNF-α is an important component of the tumor microenvironment of patients with BRCA; it is a proinflammatory cytokine secreted by macrophages or tumor cells. It can induce macrophages to differentiate into the M1 type, which is mainly used to kill tumors and promote chronic inflammation. TNF-α in peripheral blood has different implications for prognosis, depending on different expression patterns, together with other cytokines. The effect of TNF-α on BRCA cell proliferation or apoptosis depends on whether other pathways such as NF-κB are active [54]. mTOR, downstream of the PI3K/AKT pathway, consists of two complexes with similar structures but different functions: mTORC1 and mTORC2. mTORC1 enhances inner cellular anabolism to promote cell growth, while mTORC2 participates in AKT phosphorylation. Inhibitors of mTORC1 or mTORC2, for example, everolimus, have been proven effective in treating the ER+/HER2+ subtype of BRCA [2]. Additionally, TGFβ is a cytokine that promotes tissue repair and inhibits adaptive immunity, and plays an important role in epithelial-mesenchymal transition and tumor immune evasion [55]. The blockade of TGFβ signaling in CD4+ T cells could restore the responses of helper T cells and inhibit cancer progression, which is mainly mediated by the induction of vascular system reorganization [56]. VEGF not only limits antitumor immunity but also promotes pathological angiogenesis in internal tumor tissues [57]. The above-mentioned studies are consistent with our predictions based on patients with different risk scores, according to our prognostic model.

Our prognostic model found that, in the high-risk score group, mutations, such as PI3KCA and TP53, often occur together with the activation of signaling pathways, such as the Hippo pathway or the PI3K/AKT/mTOR pathway. The wild-type p53 protein cooperates with the Hippo pathway to inhibit tumor cell growth by promoting senescence or apoptosis. This function could be reversed by the mutation status of the TP53 gene (which causes the p53 protein to gain GOF activity) [58]. The hotspot mutation H1047R in PI3KCA, frequently detected in BRCA, activates the PI3K/AKT/mTOR pathway and promotes BRCA tumor formation [59]. The described relationships between mutations and activated pathways are consistent with our results. Moreover, our risk score prognostic model may provide new insights into the specific mechanisms of other mutations in BRCA.

According to our model, patients with high-risk scores are sensitive to drugs such as cisplatin. Similarly, immune cells such as macrophages in patients with high-risk scores are affected. Some studies have shown that regulation of the immune microenvironment can enhance the immunological benefits of chemotherapy drugs such as cisplatin [60] and that cisplatin and other drugs can also regulate the proportion of immune cells in BRCA [61]. Interestingly, macrophage blockade enhances cisplatin response by releasing intratumoral type I interferons [62]. The described relationships between drugs and immune cells are consistent with our results. In addition, our risk score prognostic model may provide possible ideas for mechanistic studies of other chemotherapy drugs for BRCA.

Our study has some limitations. The same gene may have different prognostic significances in different BRCA subtypes. For example, mutations in TP53 may contribute to the mortality rate of luminal B or HER2-enriched subtype BRCAs but make no effort to treat luminal A subtype BRCA [63]. In the future, we intend to validate our model using different subtypes of BRCA datasets and explore the different expressions between patients with different risk scores in different datasets. This may provide new ideas and targets for BRCA treatment and prognosis of different subtypes.

5. Conclusions

In this study, we constructed a risk score model for BRCA prognosis based on four mitochondria-related genes. This model was evaluated from different perspectives, such as mutation, invasion, and sensitivity to drugs prediction, providing new insights for BRCA prognosis.

Data Availability

TCGA Breast Cancer (BRCA) datasets (https://xenabrowser.net/) and GSE103091 were used to support this study.

Disclosure

Weixu Luo and Yuanshan Han are co-first authors.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Authors’ Contributions

Weixu Luo and Yuanshan Han contributed equally to this manuscript.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (82174357 and 81874464), the Youth Science Fund Project of National Natural Science Foundation of China (82104839, 82104846), and Hunan University Innovation Platform Project (19K066).

References

S. Loibl, P. Poortmans, M. Morrow, C. Denkert, and G. Curigliano, “Breast cancer,” The Lancet, vol. 397, no. 10286, pp. 1750–1769, 2021.
View at: Publisher Site | Google Scholar
D. Miricescu, A. Totan, I. I. Stanescu-Spinu, S. C. Badoiu, C. Stefani, and M. Greabu, “PI3K/AKT/mTOR signaling pathway in breast cancer: from molecular landscape to clinical aspects,” International Journal of Molecular Sciences, vol. 22, no. 1, p. 173, 2020.
View at: Publisher Site | Google Scholar
S. K. Yeo and J. L. Guan, “Breast cancer: multiple subtypes within a tumor?” Trends in Cancer, vol. 3, no. 11, pp. 753–760, 2017.
View at: Publisher Site | Google Scholar
P. S. Hegde and D. S. Chen, “Top 10 challenges in cancer immunotherapy,” Immunity, vol. 52, no. 1, pp. 17–35, 2020.
View at: Publisher Site | Google Scholar
F. Yu, F. Quan, J. Xu et al., “Breast cancer prognosis signature: linking risk stratification to disease subtypes,” Briefings in Bioinformatics, vol. 20, no. 6, pp. 2130–2140, 2019.
View at: Publisher Site | Google Scholar
P. Ghosh, C. Vidal, S. Dey, and L. Zhang, “Mitochondria targeting as an effective strategy for cancer therapy,” International Journal of Molecular Sciences, vol. 21, no. 9, p. 3363, 2020.
View at: Publisher Site | Google Scholar
A. Avagliano, M. R. Ruocco, F. Aliotta et al., “Mitochondrial flexibility of breast cancers: a growth advantage and a therapeutic opportunity,” Cells, vol. 8, no. 5, p. 401, 2019.
View at: Publisher Site | Google Scholar
Y. Onodera, J.-M. Nam, M. Horikawa, H. Shirato, and H. Sabe, “Arf6-driven cell invasion is intrinsically linked to TRAK1-mediated mitochondrial anterograde trafficking to avoid oxidative catastrophe,” Nature Communications, vol. 9, no. 1, p. 2682, 2018.
View at: Publisher Site | Google Scholar
P. Kaur, S. Nagar, M. Bhagwat et al., “Activated heme synthesis regulates glycolysis and oxidative metabolism in breast and ovarian cancer cells,” PLoS One, vol. 16, no. 11, Article ID e0260400, 2021.
View at: Publisher Site | Google Scholar
X. Guo, C. Xiang, Z. Zhang, F. Zhang, T. Xi, and L. Zheng, “Displacement of Bax by BMF mediates STARD13 3′UTR-induced breast cancer cells apoptosis in an miRNA-depedent manner,” Molecular Pharmaceutics, vol. 15, no. 1, pp. 63–71, 2018.
View at: Publisher Site | Google Scholar
B. Chen, G. Zhang, X. Li et al., “Comparison of BRCA versus non-BRCA germline mutations and associated somatic mutation profiles in patients with unselected breast cancer,” Aging, vol. 12, no. 4, pp. 3140–3155, 2020.
View at: Publisher Site | Google Scholar
O. M. Valencia, S. E. Samuel, R. K. Viscusi, T. S. Riall, L. A. Neumayer, and H. Aziz, “The role of genetic testing in patients with breast cancer,” JAMA Surgery, vol. 152, no. 6, pp. 589–594, 2017.
View at: Publisher Site | Google Scholar
Y. Peng, H. Yu, Y. Jin et al., “Construction and validation of an immune infiltration-related gene signature for the prediction of prognosis and therapeutic response in breast cancer,” Frontiers in Immunology, vol. 12, Article ID 666137, 2021.
View at: Publisher Site | Google Scholar
S. L. Goff and D. N. Danforth, “The role of immune cells in breast tissue and immunotherapy for the treatment of breast cancer,” Clinical Breast Cancer, vol. 21, no. 1, pp. e63–e73, 2021.
View at: Publisher Site | Google Scholar
Y. Tang, C. Li, Y.-J. Zhang, and Z.-H. Wu, “Ferroptosis-related long non-coding RNA signature predicts the prognosis of head and neck squamous cell carcinoma,” International Journal of Biological Sciences, vol. 17, no. 3, pp. 702–711, 2021.
View at: Publisher Site | Google Scholar
P. Geeleher, N. Cox, and R. S. Huang, “pRRophetic: an R package for prediction of clinical chemotherapeutic response from tumor gene expression levels,” PLoS One, vol. 9, no. 9, Article ID e107468, 2014.
View at: Publisher Site | Google Scholar
C.-F. Ni, S.-M. Zhao, L.-W. Qiu et al., “Prognostic nomogram for hepatocellular carcinoma patients after transarterial chemoembolization based on des-γ-carboxy prothrombin reactivity and modified Response Evaluation Criteria in Solid Tumors,” Journal of Cancer Research and Therapeutics, vol. 17, no. 3, pp. 707–714, 2021.
View at: Publisher Site | Google Scholar
J. Wang, L. Bravo, J. Zhang et al., “Radiomics analysis derived from LGE-MRI predict sudden cardiac death in participants with hypertrophic cardiomyopathy,” Frontiers in Cardiovascular Medicine, vol. 8, Article ID 766287, 2021.
View at: Publisher Site | Google Scholar
N. Jiang, X. Zhang, Q. Chen et al., “Identification of a mitochondria-related gene signature to predict the prognosis in AML,” Frontiers in Oncology, vol. 12, Article ID 823831, 2022.
View at: Publisher Site | Google Scholar
E. A. Badr, M. F. Assar, A. M. A. Eltorgoman, A. Z. Labeeb, G. A. Breaka, and E. A. Elkhouly, “A correlation between BCL-2 modifying factor, p53 and livin gene expressions in cancer colon patients,” Biochemistry and Biophysics Reports, vol. 22, Article ID 100747, 2020.
View at: Publisher Site | Google Scholar
R. Liu, G. Wang, C. Zhang, and D. Bai, “A prognostic model for hepatocellular carcinoma based on apoptosis-related genes,” World Journal of Surgical Oncology, vol. 19, no. 1, p. 70, 2021.
View at: Publisher Site | Google Scholar
H. Xu, R. Zou, F. Li et al., “MRPL15 is a novel prognostic biomarker and therapeutic target for epithelial ovarian cancer,” Cancer Medicine, vol. 10, no. 11, pp. 3655–3673, 2021.
View at: Publisher Site | Google Scholar
J. Cesnekova, J. Spacilova, H. Hansikova, J. Houstek, J. Zeman, and L. Stiburek, “LACE1 interacts with p53 and mediates its mitochondrial translocation and apoptosis,” Oncotarget, vol. 7, no. 30, pp. 47687–47698, 2016.
View at: Publisher Site | Google Scholar
A. P. Sudheesh and R. S. Laishram, “Nuclear phosphatidylinositol-phosphate type I kinase α-coupled star-PAP polyadenylation regulates cell invasion,” Molecular and Cellular Biology, vol. 38, no. 5, 2018.
View at: Publisher Site | Google Scholar
M. B. Teixeira, M. R. Alborghetti, and J. Kobarg, “Fasciculation and elongation zeta proteins 1 and 2: from structural flexibility to functional diversity,” World Journal of Biological Chemistry, vol. 10, no. 2, pp. 28–43, 2019.
View at: Publisher Site | Google Scholar
J. Zhang, C. Li, Q. Xue et al., “An efficient carbon‐based drug delivery system for cancer therapy through the nucleus targeting and mitochondria mediated apoptotic pathway,” Small Methods, vol. 5, no. 12, Article ID e2100539, 2021.
View at: Publisher Site | Google Scholar
M. D. Wellenstein, S. B. Coffelt, D. E. M. Duits et al., “Loss of p53 triggers WNT-dependent systemic inflammation to drive breast cancer metastasis,” Nature, vol. 572, no. 7770, pp. 538–542, 2019.
View at: Publisher Site | Google Scholar
B. Verret, J. Cortes, T. Bachelot, F. Andre, and M. Arnedos, “Efficacy of PI3K inhibitors in advanced breast cancer,” Annals of Oncology, vol. 30, pp. x12–x20, 2019.
View at: Publisher Site | Google Scholar
S. P. Gao, A. J. Kiliti, K. Zhang et al., “AKT1 E17K inhibits cancer cell migration by abrogating β-catenin signaling,” Molecular Cancer Research, vol. 19, no. 4, pp. 573–584, 2021.
View at: Publisher Site | Google Scholar
Y. Yang, B. Zhao, L. Lv, Y. Yang, S. Li, and H. Wu, “FBXL10 promotes EMT and metastasis of breast cancer cells via regulating the acetylation and transcriptional activity of SNAI1,” Cell Death Discovery, vol. 7, no. 1, p. 328, 2021.
View at: Publisher Site | Google Scholar
E. M. Martin, K. A. Orlando, K. Yokobori, and P. A. Wade, “The estrogen receptor/GATA3/FOXA1 transcriptional network: lessons learned from breast cancer,” Current Opinion in Structural Biology, vol. 71, pp. 65–70, 2021.
View at: Publisher Site | Google Scholar
M. Takaku, S. A. Grimm, B. De Kumar, B. D. Bennett, and P. A. Wade, “Cancer-specific mutation of GATA3 disrupts the transcriptional regulatory network governed by Estrogen Receptor alpha, FOXA1 and GATA3,” Nucleic Acids Research, vol. 48, no. 9, pp. 4756–4768, 2020.
View at: Publisher Site | Google Scholar
F. Bai, L.-H. Zhang, X. Liu et al., “GATA3 functions downstream of BRCA1 to suppress EMT in breast cancer,” Theranostics, vol. 11, no. 17, pp. 8218–8233, 2021.
View at: Publisher Site | Google Scholar
C. A. Dreyer, K. VanderVorst, S. Free, A. Rowson-Hodel, and K. L. Carraway, “The role of membrane mucin MUC4 in breast cancer metastasis,” Endocrine-Related Cancer, vol. 29, no. 1, p. R17, 2022.
View at: Publisher Site | Google Scholar
A. R. Rowson-Hodel, J. H. Wald, J. Hatakeyama et al., “Membrane Mucin Muc4 promotes blood cell association with tumor cells and mediates efficient metastasis in a mouse model of breast cancer,” Oncogene, vol. 37, no. 2, pp. 197–207, 2018.
View at: Publisher Site | Google Scholar
M. F. Mercogliano, M. De Martino, L. Venturutti et al., “TNFα-Induced mucin 4 expression elicits trastuzumab resistance in HER2-positive breast cancer,” Clinical Cancer Research, vol. 23, no. 3, pp. 636–648, 2017.
View at: Publisher Site | Google Scholar
C. Kikutake, M. Yoshihara, T. Sato, D. Saito, and M. Suyama, “Intratumor heterogeneity of HMCN1 mutant alleles associated with poor prognosis in patients with breast cancer,” Oncotarget, vol. 9, no. 70, Article ID 33337, 2018.
View at: Publisher Site | Google Scholar
G. Pruneri, A. Vingiani, and C. Denkert, “Tumor infiltrating lymphocytes in early breast cancer,” The Breast, vol. 37, pp. 207–214, 2018.
View at: Publisher Site | Google Scholar
Y. Xiao, M. Cong, J. Li et al., “Cathepsin C promotes breast cancer lung metastasis by modulating neutrophil infiltration and neutrophil extracellular trap formation,” Cancer Cell, vol. 39, no. 3, pp. 423–437, 2021.
View at: Publisher Site | Google Scholar
T. Z. Tien, J. N. L. W. Lee, J. C. T. Lim et al., “Delineating the breast cancer immune microenvironment in the era of multiplex immunohistochemistry/immunofluorescence,” Histopathology, vol. 79, no. 2, pp. 139–159, 2021.
View at: Publisher Site | Google Scholar
F. Pelon, B. Bourachot, Y. Kieffer et al., “Cancer-associated fibroblast heterogeneity in axillary lymph nodes drives metastases in breast cancer through complementary mechanisms,” Nature Communications, vol. 11, no. 1, p. 404, 2020.
View at: Publisher Site | Google Scholar
M. Oshi, F. A. Angarita, Y. Tokumaru et al., “High expression of NRF2 is associated with increased tumor-infiltrating lymphocytes and cancer immunity in ER-positive/HER2-negative breast cancer,” Cancers, vol. 12, no. 12, p. 3856, 2020.
View at: Publisher Site | Google Scholar
H. Wang, S. Guo, S.-J. Kim et al., “Cisplatin prevents breast cancer metastasis through blocking early EMT and retards cancer growth together with paclitaxel,” Theranostics, vol. 11, no. 5, pp. 2442–2459, 2021.
View at: Publisher Site | Google Scholar
H. Fallahi and R. Godini, “System-level responses to cisplatin in pro-apoptotic stages of breast cancer MCF-7 cell line,” Computational Biology and Chemistry, vol. 83, Article ID 107155, 2019.
View at: Publisher Site | Google Scholar
R. Roskoski, “Properties of FDA-approved small molecule protein kinase inhibitors: a 2022 update,” Pharmacological Research, vol. 175, Article ID 106037, 2022.
View at: Publisher Site | Google Scholar
N. M. Ayoub, D. R. Ibrahim, A. E. Alkhalifa, and B. A. Al-Husein, “Crizotinib induced antitumor activity and synergized with chemotherapy and hormonal drugs in breast cancer cells via downregulating MET and estrogen receptor levels,” Investigational New Drugs, vol. 39, no. 1, pp. 77–88, 2021.
View at: Publisher Site | Google Scholar
C. Tu, R. Xu, M. Koleti, and J. Zoldan, “Glycogen synthase kinase-3 inhibition sensitizes human induced pluripotent stem cells to thiol-containing antioxidants induced apoptosis,” Stem Cell Research, vol. 23, pp. 182–187, 2017.
View at: Publisher Site | Google Scholar
C. A. Berdan, R. Ho, H. S. Lehtola et al., “Parthenolide covalently targets and inhibits focal adhesion kinase in breast cancer cells,” Cell Chemical Biology, vol. 26, no. 7, pp. 1027–1035, 2019.
View at: Publisher Site | Google Scholar
A. Al-Ghabkari, D. O. Qasrawi, M. Alshehri, and A. Narendran, “Focal adhesion kinase (FAK) phosphorylation is a key regulator of embryonal rhabdomyosarcoma (ERMS) cell viability and migration,” Journal of Cancer Research and Clinical Oncology, vol. 145, no. 6, pp. 1461–1469, 2019.
View at: Publisher Site | Google Scholar
Q. Wu, W. Wu, V. Jacevic, T. C. C. Franca, X. Wang, and K. Kuca, “Selective inhibitors for JNK signalling: a potential targeted therapy in cancer,” Journal of Enzyme Inhibition and Medicinal Chemistry, vol. 35, no. 1, pp. 574–583, 2020.
View at: Publisher Site | Google Scholar
T. Murayama, Y. Takeuchi, K. Yamawaki et al., “MCM10 compensates for Myc‐induced DNA replication stress in breast cancer stem‐like cells,” Cancer Science, vol. 112, no. 3, pp. 1209–1224, 2021.
View at: Publisher Site | Google Scholar
N. Malik, H. Yan, N. Moshkovich et al., “The transcription factor CBFB suppresses breast cancer through orchestrating translation and transcription,” Nature Communications, vol. 10, no. 1, p. 2071, 2019.
View at: Publisher Site | Google Scholar
C. Wei, Y. Wang, and X. Li, “The role of Hippo signal pathway in breast cancer metastasis,” OncoTargets and Therapy, vol. 11, pp. 2185–2193, 2018.
View at: Publisher Site | Google Scholar
D. Cruceriu, O. Baldasici, O. Balacescu, and I. Berindan-Neagoe, “The dual role of tumor necrosis factor-alpha (TNF-α) in breast cancer: molecular insights and therapeutic approaches,” Cellular Oncology, vol. 43, no. 1, pp. 1–18, 2020.
View at: Publisher Site | Google Scholar
A. Bagati, S. Kumar, P. Jiang et al., “Integrin αvβ6-TGFβ-SOX4 pathway drives immune evasion in triple-negative breast cancer,” Cancer Cell, vol. 39, no. 1, pp. 54–67, 2021.
View at: Publisher Site | Google Scholar
S. Li, M. Liu, M. H. Do et al., “Cancer immunotherapy via targeted TGF-β signalling blockade in TH cells,” Nature, vol. 587, no. 7832, pp. 121–125, 2020.
View at: Publisher Site | Google Scholar
M. Schmittnaegel, N. Rigamonti, E. Kadioglu et al., “Dual angiopoietin-2 and VEGFA inhibition elicits antitumor immunity that is enhanced by PD-1 checkpoint blockade,” Science Translational Medicine, vol. 9, no. 385, Article ID eaak9670, 2017.
View at: Publisher Site | Google Scholar
M. Ferraiuolo, L. Verduci, G. Blandino, and S. Strano, “Mutant p53 protein and the Hippo transducers YAP and TAZ: a critical oncogenic node in human cancers,” International Journal of Molecular Sciences, vol. 18, no. 5, p. 961, 2017.
View at: Publisher Site | Google Scholar
X.-L. Liu, J.-L. Liu, Y.-C. Xu et al., “Membrane metallo-endopeptidase mediates cellular senescence induced by oncogenic PIK3CA H1047R accompanied with pro-tumorigenic secretome,” International Journal of Cancer, vol. 145, no. 3, pp. 817–829, 2019.
View at: Publisher Site | Google Scholar
A. Hanoteau, J. M. Newton, R. Krupar et al., “Tumor microenvironment modulation enhances immunologic benefit of chemoradiotherapy,” Journal for Immuno Therapy of Cancer, vol. 7, no. 1, p. 10, 2019.
View at: Publisher Site | Google Scholar
J. Á. Balog, L. Hackler, A. K. Kovács et al., “Single cell mass cytometry revealed the immunomodulatory effect of cisplatin via downregulation of splenic CD44+, IL-17A+ MDSCs and promotion of circulating IFN-γ+ myeloid cells in the 4T1 metastatic breast cancer model,” International Journal of Molecular Sciences, vol. 21, no. 1, p. 170, 2019.
View at: Publisher Site | Google Scholar
C. Salvagno, M. Ciampricotti, S. Tuit et al., “Therapeutic targeting of macrophages enhances chemotherapy efficacy by unleashing type I interferon response,” Nature Cell Biology, vol. 21, no. 4, pp. 511–521, 2019.
View at: Publisher Site | Google Scholar
L. Silwal-Pandit, A. Langerød, and A. L. Børresen-Dale, “TP53 mutations in breast and ovarian cancer,” Cold Spring Harb Perspect Med, vol. 7, no. 1, 2017.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Weixu Luo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

540

Downloads

731

Citations