Abstract
Background:
Depression plays a significant role in mediating breast cancer recurrence and metastasis. However, a precise risk model is lacking to evaluate the potential impact of depression on breast cancer prognosis. In this study, we established a depression-related gene (DRG) signature that can predict overall survival (OS) and elucidate its correlation with pathological parameters and sensitivity to therapy in breast cancer.
Methods:
The model training and validation assays were based on the analyses of 1,096 patients from The Cancer Genome Atlas (TCGA) database and 2,969 patients from GSE96058. A risk signature was established through univariate and multivariate Cox regression analyses.
Results:
Ten DRGs were determined to construct the risk signature. Multivariate analysis revealed that the signature was an independent prognostic factor for OS. Receiver operating characteristic (ROC) curves indicated good performance of the model in predicting 1-, 3-, and 5-year OS, particularly for patients with triple-negative breast cancer (TNBC). In the high-risk group, the proportion of immunosuppressive cells, including M0 macrophages, M2 macrophages, and neutrophils, was higher than that in the low-risk group. Furthermore, low-risk patients responded better to chemotherapy and endocrine therapy. Finally, a nomogram integrating risk score, age, tumor-node-metastasis (TNM) stage, and molecular subtypes were established, and it showed good agreement between the predicted and observed OS.
Conclusion:
The 10-gene risk model not only highlights the significance of depression in breast cancer prognosis but also provides a novel gene-testing tool to better prevent the potential adverse impact of depression on breast cancer prognosis.
Introduction
As a formidable health problem for women worldwide, breast cancer ranks first in incidence and mortality among women’s malignancies. According to the World Health Statistics 2020, breast cancer incidence accounted for 11.7% of new cases in 2020 and 6.9% of mortalities worldwide (1). Despite advances in early detection and drug development, significant barriers remain in the prediction of metastasis or recurrence risk (2). Given the long latent period and relatively young age of breast cancer, the development of prognosis-risk prediction models is of great value to improve treatment strategies and overall survival (OS).
Prognostic risk factors of breast cancer have attracted increasing interest for more than a decade. In 2007, the St. Gallen Expert Consensus considered that the risk of breast cancer recurrence was mainly related to age, tumor size, histological grade, peritumor intravascular cancer emboli, and expression of the estrogen receptor (ER)/progesterone receptor (PR) and human epidermal growth factor receptor 2 (HER2) (3). In St. Gallen 2011, breast cancer was divided into four types, namely, luminal A, luminal B, HER2+, and triple-negative breast cancer (TNBC) (4), of which TNBC was considered to have a poor prognosis. In fact, only 77% of TNBC patients survive 5 years after diagnosis, while the survival rate of other subtypes reached 93% (5). The 5-year survival rate in the metastatic TNBC is less than 30% and almost all patients will ultimately die of their disease (6). However, breast cancer is a highly heterogenous disease (7), and some TNBC patients have a low recurrence risk and relatively long survival period. It is believed that more genes are involved in determining breast cancer prognosis in addition to ER, PR and HER2. For example, the luminal A type was further classified into low-, intermediate-, and high-risk recurrence subgroups by a 21-gene recurrence score system (8). In addition, risk prediction tools, such as MammaPrint, Breast Cancer Index, PAM50 and EPclin, have also been developed to address the needs for precise diagnosis and treatment (9–12). Nevertheless, these tools are greatly limited by specific platforms and pathological subtype. Thus, it is significant and urgent to establish novel gene panels to assist with risk prediction and anti-cancer drug development.
Numerous studies have emphasized the close correlation between depression and breast cancer (13–16). The meta-analysis showed that 32.3% of breast cancer patients suffer from depression (17). Nearly half of breast cancer patients at early stage would have depression, anxiety, or both after diagnosis (18). Depression is related to a 24% increased risk of breast cancer recurrence, 30% increased risk of all-cause mortality, as well as 29% increased risk of cancer-specific mortality (19). Animal studies have demonstrated that psychological stress could accelerate tumor growth and promote lung metastasis (20–22). Thus, it is crucial to explore the role of depression-related gene (DRGs) for the prognosis of breast cancer.
Herein, a 10-DRG risk signature was determined and it presented a positive prediction with OS in The Cancer Genome Atlas (TCGA) and GSE96058 cohorts, particularly for TNBC patients. A high-risk score was related to age, tumor size, tumor stage, metastasis, and immunosuppressive cells. Moreover, low-risk patients might be highly sensitive to chemotherapy and endocrine therapy. Finally, a nomogram was successfully established by integrating the risk score and other clinical parameters. Taken together, our study not only provides a novel risk prediction model for evaluating breast cancer prognosis based on DRGs analysis but also highlights the clinical significance of depression evaluation in breast cancer treatment.
Materials and Methods
Data Source and Collection
Breast cancer related clinical information and the gene expression profiles were collected from TCGA database, which was chosen as testing dataset containing 1,096 patient samples and 113 non-tumor tissues. For further verification, the dataset GSE96058 (n=2969) downloaded from the GEO database was selected for validation.
Selection of DEGs in Breast Cancer
A total of 8,479 DRGs were acquired from the Genecard database, which supplies a comprehensive, up-to-date list of human genes. We obtained the expression values of 8,479 DRGs from TCGA. Based on TCGA dataset, DEGs were identified via Wilcoxon test after averaging replicate probes. Those genes with |log2-fold change (FC)| > 2 and adjusted P value (FDR) < 0.05 were considered DEGs.
Subsequently, PPIs of the DEGs were constructed based on the STRING (The Search Tool for the Retrieval of Interacting Genes) online tool. Then Cytoscape software 3.6.1 was applied to visualized the network. Additionally, the key modules in network were analyzed by MCODE in Cytoscape software. The options were set as degree cutoff = 2, K-Core = 4, and Node Score Cutoff = 0.3.
Identification of Prognostic Genes and Construction of a Risk Model
The univariate Cox regression analysis was performed to explore the correlation between genes and OS in the testing dataset. DEGs with p <0.05 were considered as candidate genes. Then LASSO regression was applied to shrink scope of gene screening. Finally, a multivariate Cox analysis was constructed to screen highly prognosis-associated genes and generate the prognostic risk model.
The risk score was calculated as follows:
(where “Exp” represents gene expression level and “β” represents the regression coefficient from the multivariate Cox analysis).
According to the risk score, the optimal cut-off value obtained by Receiver operating characteristic (ROC) analysis was performed for patient risk stratification.
Validation of the Constructed Model
To validate this model, Kaplan-Meier survival analysis along with log-rank test was employed to analyze survival difference between patients with different scores. To confirm whether the signature can be applied as an independent factor predictive of survival, univariate and multivariate Cox regression analyses were carried out.
GSEA
To explore the potential biological pathways, GSEA was carried out during the enrichment of the MSigDB Collection (c2.cp.kegg.v7.1.symbols.gmt). Pathways with |NES| > 1 and p < 0.05 were considered significantly enriched.
Exploration of Tumor-Infiltrating Immune Cells
As the most frequently cited tool for analyzing immune cells infiltration, CIBERSORT was used to quantify the immune cell proportion based on their RNA sequencing data (23). In the study, the tool was used to estimate the proportion of immune cell fractions in each breast patients in TCGA cohort.
Investigation of the Significance of Risk Signature for Drug Treatment
Genomics of Drug Sensitivity in Cancer website (24) was employed to discuss differences between patients with low and high score in response to chemotherapy, endocrine therapy, and targeted therapy. The IC50 of common therapeutic drugs in the TCGA breast cancer cohort was calculated.
Identification of Potential Small Molecule Drugs Targeting High Risk Patients
CMap is an online pharmacogenomic database based on gene expression data of cultured human cells treated with bioactive substance. Those up- and down-regulated target genes were uploaded to CMap. The connectivity score ranges from -1 to 1, which was used to reveal the closeness between expression profiles. The drugs with negative scores were generally considered as potential therapeutic molecules.
Construction and Evaluation of the Nomogram
The rms R package was performed to build a nomogram containing the model and clinicopathological features. The C-index and calibration plot were applied to assess predicted probabilities.
Statistical Analysis
The R software (version 3.6.3) was used to analyze the statistical data. Kaplan-Meier survival curves at 1-, 3-, and 5-years were illustrated with the survminer package. The prognostic effect of the risk model and clinicopathological features were assessed through univariate Cox analysis and multivariate Cox analysis. A two-sided p < 0.05 was considered statistically significant.
Results
Identification of DRGs and Establishment of a 10-Gene Risk Model for Predicting Breast Cancer Prognosis
To construct a risk model based on DRGs, TCGA-BC cohort and GSE96058 cohort were selected as testing and validation datasets, respectively. The flow chart of bioinformatic analysis is illustrated in Figure 1. This study included 1,027 cases from TCGA cohort in which primary breast cancer had been followed up for more than 1 month. Clinical characteristics are summarized in Table S1. In total, 358 upregulated and 400 downregulated DRGs were finally identified (Figures 2A, B). To narrow the gene set and determine the hub genes, a protein–protein interaction (PPI) network based on the differentially expressed genes (DEGs) was built. Subsequently, module analyses of the network were conducted using Molecular Complex Detection (MCODE) to obtain hub genes (Figure 2C). Five significant modules were screened out and genes identified as hub genes were further analyzed.
Figure 1
Figure 2
In univariate Cox analysis, 34 DRGs were found significantly related to OS (Figure 2D). Through lasso regression analyses, 15 genes were further identified as prognosis-related genes (Figure 2E). Multivariate Cox analysis confirmed the following 10 DRGs for constructing the risk model: MT3, SORBS1, IGFALS, AMH, IL12B, TP53AIP1, PXDNL, MC5R, FOXD1 and LHX1 (Table S2). The formula was developed as below:
Prognostic Effect of the Model on Predicting OS of Breast Cancer
To explore the predictive effect of the model, risk scores in TCGA-Breast cancer and GSE96058 cohorts were calculated according to the formula. Subsequently, patients were divided into either low- or high-risk group based on the cut-off point of −1.281, as determined by MedCalc version 15.2.2 (MedCalc Software bvba, Ostend, Belgium). Survival analysis demonstrated a shorter survival time in high-risk patients (p < 0.001), which was verified by the GSE96058 dataset (p < 0.001) (Figure 3A). As illustrated in the Figures 3B, C, the higher the score, the higher proportion of patients dying in TCGA and GSE96058 datasets. Heatmap analysis showed that MT3, SORBS1, IGFALS, IL12B and TP53AIP1 in the low-risk group were overexpressed, whereas AMH, PXDNL, MC5R, FOXD1 and LHX1 in the high-risk group were overexpressed. The outcomes in GSE96058 dataset were in line with those in TCGA (Figure 3D). As displayed in Figure 3E, the maximum area under the ROC curve (AUC) values reached 0.789, indicating the good sensitivity and specificity of the model. Consistently, in GSE96058 cohort the model achieved the maximum AUC of 0.725. These results revealed that the model was effective in predicting the survival probability.
Figure 3
Considering tumor heterogeneity, survival and ROC analyses of different molecular subtypes of breast cancer were carried out, respectively. The risk model had good prognostic performance, particularly for TNBC. As shown in Figure S1, the survival analyses proved the prognostic value for TNBC in both TCGA and GSE96058 cohorts (p < 0.001). ROC analyses further revealed that in the TCGA cohort, the AUC of TNBC was 0.709, 0.732, and 0.730 at 1, 3 and 5 years, respectively. The results were validated in the GSE96058 dataset in which AUC at 1, 3, and 5 years reached 0.747, 0.657, and 0.609 for TNBC, respectively (Figure S2).
Correlation Analysis Between the Model and Pathological Parameters
The relationship between the risk model and pathological parameters was further investigated. As displayed in Figures 4A, B, the risk model in the TCGA testing cohort was closely associated with age, molecular subtypes, clinical stage, T stage, and M stage, whereas N stage was not. Similarly, significant differences were also observed in GSE96058 cohorts, including age, molecular subtypes, tumor size, and positive nodes (Figure S3A, B). The risk score was significant for prognosis prediction when evaluated by univariate Cox analysis (HR = 2.926, 95% CI: 2.306–3.711) (Figure 4C). Furthermore, multivariable Cox analysis revealed that the model remained statistically significant (HR = 1.075, 95% CI:1.058–1.092), indicating that the model is an independent prognostic factor (Figure 4D). Likewise, in GSE96058 cohort, the univariate Cox regression analysis demonstrated that risk score (HR=2.329, 95% CI:1.879–2.887) is significantly related to survival (Figure S3C). The multivariate Cox regression analysis indicated that risk score (HR=1.648, 95% CI:1.324–2.049) is still an independent prognostic factor for breast cancer (Figure S3D).
Figure 4
Molecular Mechanism Analysis of the Model in Predicting Prognosis
We continued to explore the potential biological mechanisms of the risk model. GSEA revealed that group at high-risk was enriched in terpenoid backbone biosynthesis, cell cycle, pentose, and glucuronate interconversions, DNA replication, steroid biosynthesis, fructose, and mannose metabolism as well as mismatch repair. By contrast, the group at low-risk was enriched in pathways such as cell adhesion molecules, intestinal immune network for IgA production, cytokine-cytokine receptor interaction, JAK-STAT signaling pathway, chemokine signaling pathway, T cell receptor signaling pathway, and leukocyte transendothelial migration (Figure 5A).
Figure 5
Growing evidence shows that depression could inhibit the immune effector cells and facilitate immune escape, thus protecting cancer cells from immunological killing (14, 25). Moreover, the preliminary GSEA analyses in this study illustrated that the signature is closely associated with immune pathways. Herein, the abundance of 22 immune cells in low-risk and high-risk cases was analyzed via the CIBERSORT method (Figure 5B). The results revealed that high-risk patients had higher abundances of immunosuppressive cells including M0 macrophages, M2 macrophages, and neutrophils, but significantly lower abundances of M1 macrophages, naive B cells, CD8 T cells, resting memory CD4 T cells, regulatory T cells, resting dendritic cells, and resting mast cells (Figure 5C). In summary, the immunosuppressive microenvironment might be responsible for worse prognosis in high-risk cases.
Then the expression of the immune checkpoint regulators was analyzed in patients with different risk score. As shown in Figure 5D, the expressions of PD-1, PD-L1, and CTLA-4 increased significantly in patients with low-scores (p < 0.05), whereas LAG3 showed no statistical difference. These results implied that low-risk patients might be sensitive to immune checkpoint blockade therapies.
Additionally, it is well established that cytokines play a crucial role in the tumor immune microenvironment. Therefore, the results in TCGA cohort demonstrated that GZMA, GZMB, and IFNG in the low-risk group were highly expressed, while TGFB1 was elevated in high-risk breast cancer patients (p < 0.001) (Figure 5E).
Clinical Effects of the Model in Predicting Treatment Sensitivity
We further identified associations between the model and treatment effect of chemotherapeutics, endocrine drugs, and targeted agents. Low-risk patients had lower half-maximal inhibitory concentration (IC50) of chemotherapeutics such as cisplatin, doxorubicin, cytoxan, methotrexate, and vinorelbine except for paclitaxel, indicating the predictive potential of the model for chemosensitivity (Figure 6A). Regarding endocrine therapy, low-risk patients had lower IC50 of endocrine therapies, such as tamoxifen, fulvestrant, and letrozole (Figure 6B). Low-risk patients had lower IC50 of targeted therapies, such as trastuzumab, olaparib, and temsirolimus, whereas high-risk patients had lower IC50 for targeted therapies such as palbociclib and sorafenib. The two groups showed no statistical differences for lapatinib (Figure 6C). Moreover, the DEGs were subjected to Connectivity Map (CMap) analysis and a total of six drugs, including pronetalol, puromycin, chlorphenamine, megestrol, semustine, and acacetin, were indicated for high-risk breast patients (Figure 6D).
Figure 6
Nomogram Construction Based on the 10-Gene Risk Model
To create a quantitative method to predict OS, the nomogram for predicting 1-year, 3-year and 5-year OS was developed. Four risk factors, including age, stage, molecular subtypes, and risk score, were included in the nomogram (Figure 7A). The C-index reached 0.782 (95%CI: 0.733–0.831). As shown by the calibration curves, the actual and predicted OS matched well (Figures 7B–D).
Figure 7
Discussion
In this study, a novel prognostic model consisting of 10-DRGs was established to enhance predictive capability of OS in breast cancer. The findings demonstrated that the model can effectively distinguish high-risk from low-risk patients and the model was effective in predicting OS in both the testing and validation sets, particularly for TNBC patients. Additionally, when pathological factors were included in multivariate Cox analysis, the model remained as an independent prognostic factor. These results implied that the risk model is a robust prognostic tool. Moreover, it presented that a high score was related to poor survival in TNBC patients. Given that TNBC is the most challenging subtype of breast cancer, predicting prognosis of TNBC by this risk model is of great significance.
Among the 10 genes, five (PXDNL, MC5R, AMH, FOXD1, and LHX1) were high-risk genes and five (IGFALS, SORBS1, IL12B, MT3, and TP53AIP1) were protective factors. PXDNL is reported as a susceptibility gene in patients with depression (26). MC5R, a melanocortin receptor, could mediate the axis responsiveness to integrated signals from diurnal rhythms and cortisol negative feedback. It was reported that MC5R antagonist could treat the depressive or generalized anxiety disorder (27). MT3 may be related to the antidepressant-like activity of eugenol, resulting in increased expression of the brain-derived neurotrophic factor (28). Furthermore, compared with those of the wild-type mice, Koumura et al. (29) found that MT3 KO mice had shorter social interactions duration as well as diminished prepulse inhibition for the acoustic startle response, indicating abnormal psychological behavior in schizophrenia, anxiety, autism, and phencyclidine-induced psychosis. In general, the 10 genes are related to depression but are not solely related to depression. The mental diseases including depression, anxiety, and panic are usually interrelated. Patients with depression, anxiety, or panic experienced declines in quality of life and showed poor treatment adherence as compared with patients without mental disorders (30). Our previous meta-analysis indicated that depression and anxiety would increase the risk of breast cancer recurrence and mortality (19). Likewise, the 10 genes included in this study are not only depression-related but also related to breast cancer progression. Gomulkiewicz et al. (31) demonstrated that patients with ductal breast cancer had lower expression of MT3 than that in non-malignant breast tissue and the level of MT3 in breast cancer patients with lymph node metastasis decreased compared to patients without metastasis. Moreover, the expression of MT3 in breast cancer cell lines was lower than that in the normal human breast epithelial cell lines. These findings indicated that MT3 might be correlated with the malignant transformation of breast epithelial cells and tumor progression. In addition, PXDNL and FOXD1 were reportedly involved in breast cancer pathogenesis and were significant in predicting prognosis. FOXD1 was found to be significantly related to the prognosis of basal-like breast cancer (32). Zhao et al. (33) demonstrated that FOXD1 could increase cell proliferation and enhance chemoresistance of cancer cells by targeting p27. It was demonstrated that decreased levels of SORBS1 have significant correlation with worse survival, distant metastasis, and more malignant phenotype. SORBS1 could inhibit tumorigenesis and metastasis through preventing JNK activation and attenuate cisplatin chemotherapy through p53 accumulation in breast cancer (34). TP53AIP1 inhibits proliferation and growth of breast cancer cells through cell-cycle arrest, apoptosis induction, elevation of the expression of cleaved-caspase-3, cleaved-caspase-9, Bax, p53, and the inhibition of Bcl-2, Ki67, and PI3K/Akt pathways (35). Although reports on the prognostic effect of the remaining genes in breast cancer are limited, it is worth further evaluating their potential as biomarkers of breast cancer progression.
Furthermore, the model was associated with tumor-infiltrating immune cells, as the abundance of M2 macrophages and neutrophils in high-risk patients significantly increased. The abundance of M1 macrophages, naive B cells, CD8 T cells, resting memory CD4 T cells, resting dendritic cells, and resting mast cells was lower in the high-risk group and was related to better OS. It is well established that cytokines play a crucial role in the tumor immune microenvironment. As high-risk patients in our study had higher abundances of immunosuppressive cells, but significantly lower proportions of M1 macrophages, naive B cells, CD8 T cells, resting memory CD4 T cells, resting dendritic cells, resting mast cells, and immune inhibitory cytokines such as TGFB1 and biomarkers of T cell activation such as GZMA, GZMB, and IFNG are selected to explore the differences between high-risk and low-risk groups. Notably, the results in TCGA cohort demonstrated that GZMA, GZMB, and IFNG in the low-risk cases were highly expressed, while TGFB1 was elevated in the high-risk group. Future research should be conducted on the mechanisms of this model in influencing immune cells in tumor microenvironment. Accumulating evidence also demonstrated the potential of immune checkpoint inhibitors (ICIs) in breast cancer treatment (36); however, only a few patients obtained therapeutic benefit from ICIs. It was found that low-risk patients had increased expression of PD-1, PD-L1, as well as CTLA-4, suggesting that low-risk patients are promising candidates for ICIs. In addition, the model provided an opportunity to choose the optimal therapy for breast cancer patients. Our study revealed that low-risk patients seem to be more sensitive to chemotherapy and endocrine therapy, which is in line with better OS of the low-risk patients. These results indicate the potential effect of the model in predicting therapy sensitivity, which is also beneficial for selecting optimal therapeutic strategies and improving breast cancer patients’ prognosis.
Furthermore, our study provides potential therapeutic options for high-risk patients. Among them, pronethalol, a β-adrenoceptor antagonist, is used to treat coronary heart disease and arrhythmias. Clinical studies have demonstrated that the application of β-blocker could improve the prognosis of patients with lung cancer, ovarian cancer, and melanoma (37–39). Although the use of pronethalol in clinical practice was restricted because it proved to be a carcinogen in mice (40), other β-blocker agents, such as propranolol, may be a valid approach for patients at high risk. It was reported that propranolol inhibits metastasis, reduces recurrence, and cancer-specific mortality in several clinical studies (41, 42). In addition, the safety of propranolol in breast cancer patients has been proven in a phase II pilot study (43). Chlorpheniramine, an alkylamine antihistamine, is used to prevent allergic symptoms. Although not generally recognized as an antidepressant or anxiolytic drug, chlorphenamine seems to have these properties as well. Shoko et al. (44) found that chlorpheniramine exerted antidepressant-like effects by activating dopamine D1 and α1-adrenoceptors. This result indicates that antidepressants may be another choice for breast cancer patients. It was proved that application of antidepressants for breast cancer may reduce mortality risk (45). However, antidepressants, such as selective serotonin re-uptake inhibitor, can reduce tamoxifen’s effectiveness and increase death risk (46). Thus, further studies are warranted for validation. Acacetin is a natural flavonoid, which can inhibit the secretion of carcinogenic estrogen metabolites, as well as tumor cell proliferation, invasion, and migration (47, 48). It is worth exploring its potential application in regulating DRGs during breast cancer treatment.
Conclusion
This study provides a novel 10-DRGs risk model to predict breast cancer survival. The model is particularly significant for TNBC and is promising for predicting drug sensitivity and is thus helpful for designing optimal therapeutic strategies to improve clinical prognosis. Our findings further highlight the significance of monitoring and treating psychological stress in preventing recurrence and metastasis of breast cancer. Further experimental validation and prospective clinical trials are worth conducting.
Funding
This work was supported by the National Natural Science Foundation of China [81873306, 81973526, 82004373, 82104869]; the State Key Laboratory of Dampness Syndrome of Chinese Medicine [SZ2021ZZ19]; Science and Technology Planning Project of Guangdong Province [2021A0505030059, 2017B030314166]; Department of Education of Guangdong Province [2018KZDXM022, A1-2606-19-111-009, 2019KQNCX019]; The 2020 Guangdong Provincial Science and Technology Innovation Strategy Special Fund (Guangdong-Hong Kong-Macau Joint Lab) [2020B1212030006]; Guangdong traditional Chinese medicine bureau project [20201132, 20211114]; Guangzhou science and technology project [202102010316, 202102010241, 201904010407]; The Specific Research Fund for TCM Science and Technology of Guangdong provincial Hospital of Chinese Medicine [YN2018MJ07, YN2018QJ08], and the Foundation for Young Scholars of Guangzhou University of Chinese Medicine [QNYC20190101].
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Author contributions
ZW and XW conceived and designed the study. XW and NW conducted the bioinformatic analysis and wrote the manuscript. LZ, WY, and KS participated in data interpretation and revised the manuscript. SW, YZ, BY, JZ, and BP contributed to data collection and discussion.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2022.879563/full#supplementary-material
References
1
SungHFerlayJSiegelRLLaversanneMSoerjomataramIJemalAet al. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin (2021) 71:209–49. doi: 10.3322/caac.21660
2
KimbungSLomanNHedenfalkI. Clinical and Molecular Complexity of Breast Cancer Metastases. Semin Cancer Biol (2015) 35:85–95. doi: 10.1016/j.semcancer.2015.08.009
3
GoldhirschAWoodWCGelberRDCoatesASThürlimannBSennHJ. Progress and Promise: Highlights of the International Expert Consensus on the Primary Therapy of Early Breast Cancer 2007. Ann Oncol (2007) 18:1133–44. doi: 10.1093/annonc/mdm271
4
GoldhirschAWoodWCCoatesASGelberRDRlimannBTSennH-J. Strategies for Subtypes—Dealing With the Diversity of Breast Cancer: Highlights of the St Gallen International Expert Consensus on the Primary Therapy of Early Breast Cancer. Ann Oncol (2011) 22:1736–47. doi: 10.1093/annonc/mdr304
5
BauerKRBrownMCressRDPariseCACaggianoV. Descriptive Analysis of Estrogen Receptor (ER)-Negative, Progesterone Receptor (PR)-Negative, and HER2-Negative Invasive Breast Cancer, the So-Called Triple-Negative Phenotype: A Population-Based Study From the California Cancer Registry. Cancer (2010) 109:1721–8. doi: 10.1002/cncr.22618
6
LiedtkeCMazouniCHessKRAndréFTordaiAMejiaJAet al. Response to Neoadjuvant Therapy and Long-Term Survival in Patients With Triple-Negative Breast Cancer. J Clin Oncol (2008) 26:1275–81. doi: 10.1200/jco.2007.14.4147
7
ZardavasDIrrthumASwantonCPiccartM. Clinical Management of Breast Cancer Heterogeneity. Nat Rev Clin Oncol (2015) 12:381–94. doi: 10.1038/nrclinonc.2015.73
8
SparanoJAGrayRJRavdinPMMakowerDFPritchardKIAlbainKSet al. Clinical and Genomic Risk to Guide the Use of Adjuvant Therapy for Breast Cancer. N Engl J Med (2019) 380:2395–405. doi: 10.1056/NEJMoa1904819
9
CardosoFvan't VeerLJBogaertsJSlaetsLVialeGDelalogeSet al. 70-Gene Signature as an Aid to Treatment Decisions in Early-Stage Breast Cancer. N Engl J Med (2016) 375:717–29. doi: 10.1056/NEJMoa1602253
10
FilipitsMRudasMJakeszRDubskyPFitzalFSingerCFet al. A New Molecular Predictor of Distant Recurrence in ER-Positive, HER2-Negative Breast Cancer Adds Independent Information to Conventional Clinical Risk Factors. Clin Cancer Res (2011) 17:6012–20. doi: 10.1158/1078-0432.Ccr-11-0926
11
SestakIBuusRCuzickJDubskyPKronenwettRDenkertCet al. Comparison of the Performance of 6 Prognostic Signatures for Estrogen Receptor-Positive Breast Cancer: A Secondary Analysis of a Randomized Clinical Trial. JAMA Oncol (2018) 4:545–53. doi: 10.1001/jamaoncol.2017.5524
12
PuMMesserKDaviesSRVickeryTLPittmanEParkerBAet al. Research-Based PAM50 Signature and Long-Term Breast Cancer Survival. Breast Cancer Res Treat (2020) 179:197–206. doi: 10.1007/s10549-019-05446-y
13
ReichMLesurAPerdrizet-ChevallierC. Depression, Quality of Life and Breast Cancer: A Review of the Literature. Breast Cancer Res Treat (2008) 110:9–17. doi: 10.1007/s10549-007-9706-5
14
BorgiMCollacchiBOrtonaECirulliF. Stress and Coping in Women With Breast Cancer:Unravelling the Mechanisms to Improve Resilience. Neurosci Biobehav Rev (2020) 119:406–21. doi: 10.1016/j.neubiorev.2020.10.011
15
CarreiraHWilliamsRMüllerMHarewoodRStanwaySBhaskaranK. Associations Between Breast Cancer Survivorship and Adverse Mental Health Outcomes: A Systematic Review. J Natl Cancer Inst (2018) 110:1311–27. doi: 10.1093/jnci/djy177
16
JohnsSAStutzPVTalibTLCoheeAABeck-CoonKABrownLFet al. Acceptance and Commitment Therapy for Breast Cancer Survivors With Fear of Cancer Recurrence: A 3-Arm Pilot Randomized Controlled Trial. Cancer (2020) 126:211–8. doi: 10.1002/cncr.32518
17
PilevarzadehMAmirshahiMAfsargharehbaghRRafiemaneshHHashemiSMBalouchiA. Global Prevalence of Depression Among Breast Cancer Patients: A Systematic Review and Meta-Analysis. Breast Cancer Res Treat (2019) 176:519–33. doi: 10.1007/s10549-019-05271-3
18
BurgessCCorneliusVLoveSGrahamJRamirezAJB. Depression and Anxiety in Women With Early Breast Cancer: Five Year Observational Cohort Study. BMJ (2005) 330:702. doi: 10.1136/bmj.38343.670868.D3
19
WangXWangNZhongLWangSZhengYYangBet al. Prognostic Value of Depression and Anxiety on Breast Cancer Recurrence and Mortality: A Systematic Review and Meta-Analysis of 282,203 Patients. Mol Psychiatry (2020) 25:3186–97. doi: 10.1038/s41380-020-00865-6
20
ChenHLiuDGuoLChengXGuoNShiM. Chronic Psychological Stress Promotes Lung Metastatic Colonization of Circulating Breast Cancer Cells by Decorating a Pre-Metastatic Niche Through Activating β-Adrenergic Signaling. J Pathol (2018) 244:49–60. doi: 10.1002/path.4988
21
FlahertyRLIntabliHFalcinelliMBuccaGHeskethAPatelBAet al. Stress Hormone-Mediated Acceleration of Breast Cancer Metastasis is Halted by Inhibition of Nitric Oxide Synthase. Cancer Lett (2019) 459:59–71. doi: 10.1016/j.canlet.2019.05.027
22
CuiBLuoYTianPPengFLuJYangYet al. Stress-Induced Epinephrine Enhances Lactate Dehydrogenase A and Promotes Breast Cancer Stem-Like Cells. J Clin Invest (2019) 129:1030–46. doi: 10.1172/jci121685
23
NewmanAMLiuCLGreenMRGentlesAJFengWXuYet al. Robust Enumeration of Cell Subsets From Tissue Expression Profiles. Nat Methods (2015) 12:453–7. doi: 10.1038/nmeth.3337
24
YangWSoaresJGreningerPEdelmanEJLightfootHForbesSet al. Genomics of Drug Sensitivity in Cancer (GDSC): A Resource for Therapeutic Biomarker Discovery in Cancer Cells. Nucleic Acids Res (2013) 41:D955–961. doi: 10.1093/nar/gks1111
25
BlumeJDouglasSDEvansDL. Immune Suppression and Immune Activation in Depression. Brain Behav Immun (2011) 25:221–9. doi: 10.1016/j.bbi.2010.10.008
26
ChenXLongFCaiBChenXQinLChenG. A Novel Relationship for Schizophrenia, Bipolar, and Major Depressive DisorderPart 8: A Hint From Chromosome 8 High Density Association Screen. Mol Neurobiol (2016) 16:5868–82. doi: 10.1007/s12035-016-0102-1
27
HrubyVJCaiMMorganC. Methods for the Treatment of Depression and Anxiety by a Melanocortin 5 Receptor Antagonist, PG-20n. United States Patent US 10653743 (2020).
28
IrieYItokazuNAnjikiNIshigeAWatanabeKKeungWM. Eugenol Exhibits Antidepressant-Like Activity in Mice and Induces Expression of Metallothionein-III in the Hippocampus. Brain Res (2004) 1011:243–6. doi: 10.1016/j.brainres.2004.03.040
29
KoumuraAKakefudaKHondaAItoYHaraHJNL. Metallothionein-3 Deficient Mice Exhibit Abnormalities of Psychological Behaviors. Neurosci Lett (2009) 467:11–4. doi: 10.1016/j.neulet.2009.09.051
30
LawrenceWRKuliszewskiMGHoslerASLeinungMCZhangXZhangWet al. Association Between Preexisting Mental Illnesses and Mortality Among Medicaid-Insured Women Diagnosed With Breast Cancer. Soc Sci Med (2021) 270:113643. doi: 10.1016/j.socscimed.2020.113643
31
GomulkiewiczAJablonskaKPulaBGrzegrzolkaJBorskaSPodhorska-OkolowMet al. Expression of Metallothionein3 in Ductal Breast Cancer. Int J Oncol (2016) 49:2487–97. doi: 10.3892/ijo.2016.3759
32
LiuQSongXLiuZYuZ. Investigation of Candidate Genes and Pathways in Basal/TNBC Patients by Integrated Analysis. Technol Cancer Res Treat (2021) 20:15330338211019506. doi: 10.1177/15330338211019506
33
ZhaoYFZhaoJYYueHHuKSShenHGuoZGet al. FOXD1 Promotes Breast Cancer Proliferation and Chemotherapeutic Drug Resistance by Targeting P27. Biochem Biophys Res Commun (2015) 456:232–7. doi: 10.1016/j.bbrc.2014.11.064
34
SongLChangRDaiCWuYGuoJQiMet al. SORBS1 Suppresses Tumor Metastasis and Improves the Sensitivity of Cancer to Chemotherapy Drug. Oncotarget (2017) 8:9108–22. doi: 10.18632/oncotarget.12851
35
LiangYWangSLiuJ. Overexpression of Tumor Protein P53-Regulated Apoptosis-Inducing Protein 1 Regulates Proliferation and Apoptosis of Breast Cancer Cells Through the PI3K/Akt Pathway. J Breast Cancer (2019) 22:172–84. doi: 10.4048/jbc.2019.22.e21
36
GaynorNCrownJCollinsDM. Immune Checkpoint Inhibitors: Key Trials and an Emerging Role in Breast Cancer. Semin Cancer Biol (2020) 79:44–57. doi: 10.1016/j.semcancer.2020.06.016
37
De GiorgiVGrazziniMGandiniSBenemeiSLottiTMarchionniNet al. Treatment With β-Blockers and Reduced Disease Progression in Patients With Thick Melanoma. Arch Intern Med (2011) 171:779–81. doi: 10.1001/archinternmed.2011.131
38
DiazESKarlanBYLiAJ. Impact of Beta Blockers on Epithelial Ovarian Cancer Survival. Gynecol Oncol (2012) 127:375–8. doi: 10.1016/j.ygyno.2012.07.102
39
WangHMLiaoZXKomakiRWelshJWO'ReillyMSChangJYet al. Improved Survival Outcomes With the Incidental Use of Beta-Blockers Among Patients With Non-Small-Cell Lung Cancer Treated With Definitive Radiation Therapy. Ann Oncol (2013) 24:1312–9. doi: 10.1093/annonc/mds616
40
ShanksRGWoodTMDornhorstACClarkML. Some Pharmacological Properties of a New Adrenergic Beta-Receptor Antagonist. Nature (1966) 212:88–90. doi: 10.1038/212088a0
41
BarronTIConnollyRMSharpLBennettKVisvanathanK. Beta Blockers and Breast Cancer Mortality: A Population- Based Study. J Clin Oncol (2011) 29:2635–44. doi: 10.1200/jco.2010.33.5422
42
Melhem-BertrandtAChavez-MacgregorMLeiXBrownENLeeRTMeric-BernstamFet al. Beta-Blocker Use is Associated With Improved Relapse-Free Survival in Patients With Triple-Negative Breast Cancer. J Clin Oncol (2011) 29:2645–52. doi: 10.1200/jco.2010.33.4441
43
ShaashuaLShabat-SimonMHaldarRMatznerPZmoraOShabtaiMet al. Perioperative COX-2 and β-Adrenergic Blockade Improves Metastatic Biomarkers in Breast Cancer Patients in a Phase-II Randomized Trial. Clin Cancer Res (2017) 23:4651–61. doi: 10.1158/1078-0432.Ccr-17-0152
44
ShokoHMiyataSOnoderaKKameiJ. Involvement of Dopamine D1 Receptors and Alpha1-Adrenoceptors in the Antidepressant-Like Effect of Chlorpheniramine in the Mouse Tail Suspension Test. Eur J Pharmacol (2007) 562:72–6. doi: 10.1016/j.ejphar.2007.01.063
45
ShimEJLeeJWChoJJungHKKimNHLeeJEet al. Association of Depression and Anxiety Disorder With the Risk of Mortality in Breast Cancer: A National Health Insurance Service Study in Korea. Breast Cancer Res Treat (2020) 179:491–8. doi: 10.1007/s10549-019-05479-3
46
KellyCMJuurlinkDNGomesTDuong-HuaMPritchardKIAustinPCet al. Selective Serotonin Reuptake Inhibitors and Breast Cancer Mortality in Women Receiving Tamoxifen: A Population Based Cohort Study. Bmj (2010) 340:c693. doi: 10.1136/bmj.c693
47
RenHMaJSiLRenBChenXWangDet al. Low Dose of Acacetin Promotes Breast Cancer MCF-7 Cells Proliferation Through the Activation of ERK/ PI3K /AKT and Cyclin Signaling Pathway. Recent Pat Anticancer Drug Discov (2018) 13:368–77. doi: 10.2174/1574892813666180420154012
48
JonesAAGehlerS. Acacetin and Pinostrobin Inhibit Malignant Breast Epithelial Cell Adhesion and Focal Adhesion Formation to Attenuate Cell Migration. Integr Cancer Ther (2020) 19:1534735420918945. doi: 10.1177/1534735420918945
Summary
Keywords
breast cancer, depression, predictive model, overall survival, nomogram
Citation
Wang X, Wang N, Zhong LLD, Su K, Wang S, Zheng Y, Yang B, Zhang J, Pan B, Yang W and Wang Z (2022) Development and Validation of a Risk Prediction Model for Breast Cancer Prognosis Based on Depression-Related Genes. Front. Oncol. 12:879563. doi: 10.3389/fonc.2022.879563
Received
19 February 2022
Accepted
01 April 2022
Published
10 May 2022
Volume
12 - 2022
Edited by
Shaoquan Zheng, Sun Yat-sen University Cancer Center (SYSUCC), China
Reviewed by
Sen Li, Beijing University of Chinese Medicine, China; Weihuan Huang, Jinan University, China
Updates
Copyright
© 2022 Wang, Wang, Zhong, Su, Wang, Zheng, Yang, Zhang, Pan, Yang and Wang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhiyu Wang, wangzhiyu976@126.com
†These authors have contributed equally to this work
This article was submitted to Molecular and Cellular Oncology, a section of the journal Frontiers in Oncology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.