- 1Department of Breast and Thyroid Surgery, The Second Affiliated Hospital of Chongqing Medical University, Chongqing, China
- 2Department of Cardiology, The Second Affiliated Hospital, Chongqing Medical University, Chongqing, China
Background: Compelling evidence has demonstrated the pivotal role of autophagy in the prognosis of breast cancer. Breast cancer (BC) patients with early relapse consistently exhibited worse survival.
Methods: The autophagy-related genes were derived from the Human Autophagy Database (HADb) and high-sequencing data were obtained from The Cancer Genome Atlas (TCGA). Discrepantly expressed autophagy genes (DEAGs) between early relapse and long-term survival groups were performed using the Linear Models for Microarray data (LIMMA) method. Lasso Cox regression analysis was conducted for the selection of the 4-gene autophagy-related gene signature. GSE42568 and GSE21653 databases were enrolled in this study for the external validation of the signature. Then patients were divided into high and low-risk groups based on the specific score formula. GSEA was used to discover the related signaling pathway. The Kaplan-Meier curves and the receiver operating characteristic (ROC) curves were used to evaluate the discrimination and accuracy of the 4-gene signature.
Results: A signature composed of four autophagy-related mRNA including APOL1, HSPA8, SIRT1, and TP73, was identified as significantly associated with the early relapse in BC patients. Time-dependent receiver-operating characteristic at 1 year suggested remarkable accuracy of the signature [area under the curve (AUC = 0.748)]. The risk score model based on the autophagy-related signature showed favorable predicting value in 1-, 2-, and 3-year relapse-free survival (RFS) in training and two validating cohorts. The GSEA displayed gene sets were remarkably enriched in carcinogenic activation pathways and autophagy-related pathways. The nomogram involving three variables (progesterone receptor status, T stage, and 4-gene signature) exhibited relatively good discrimination with a C-index of 0.766.
Conclusions: Our study establishes an autophagy-related 4-gene signature that can effectively stratify the high-risk and low-risk BC patients for early relapse. Combined with the clinicopathological variables, the signature could significantly help oncologists tailor more efficient treatment strategies for BC patients.
Introduction
Breast cancer (BC) is currently the most frequent malignancy and one of the leading causes of cancer death in the United States (estimated 279,100 new cases and 42,690 death) (1) and China mainland (estimated 304,000 new cases and 70,000 death) (2). Although the long-term survival of patients with BC has been significantly increased in the past years with the application of targeted therapy (3), endocrine therapy (4), and even immunotherapy (5, 6), early relapse (2 years after initial treatment) with metastasis could reverse this favorable outcome (7, 8). Regardless of the prognosis, all women with BC are at risk for early recurrence. According to a recent review report, nearly 50% of early recurrences occur within 5 years of surgery, and they peak at 2 years after surgery in women treated with adjuvant tamoxifen (9). Besides, early relapse in BC patients is frequently associated with poor clinicopathological features [such as young age (10), late TNM stage, poor differentiation grade, and worse histopathological type (11, 12)] and resistance to adjuvant chemotherapy or endocrine therapy (13–16). Those cases who developed early relapse consistently tended to have poorer long-term survival rates. Notably, a recent study has demonstrated that BC patients experienced altered hormone receptor and HER2 status throughout tumor progression, which significantly influences survival (17). Thus, for the great heterogeneity of BC, the prognosis varies significantly in BC patients with the same stage and comparable clinicopathological features. For this reason, hall markers and other biological indicators could help to predict the recurrence of BC (18).
Autophagy is a routine physiological process associated with aging and human disease via guiding the degradation of damaged, denatured, or senescent proteins and organelles in lysosomes (19, 20). Accordingly, compelling evidence has demonstrated that autophagy plays a pivotal role in tumor growth, metastasis, and recurrence of BC, which could maintain the homeostasis and the survival of BC cells by removing dysfunctional or unnecessary substances (21–24). On the other hand, accumulating evidence showed that autophagy-related genes were significantly involved in the regulations of the autophagy process. For instance, recent two basic research demonstrated that MTA1 (metastasis-associated 1) (13) and long noncoding RNA H19 (14) were the regulators of autophagy in resistance to the endocrine therapy (tamoxifen). However, Marsh et al. (25) discovered autophagy could inhibit the metastasis of BC cells by accumulating the autophagy cargo receptor (ACR), neighbor to BRCA1. A number of coding RNA (mRNA) and non-coding RNAs (microRNA, lncRNA, and circRNA) signatures have been identified for predicting the proliferation and prognosis of BC patients (26–30). Nonetheless, most of these signatures focused on overall survival, there is still a lack of work on investigating the impact of mRNA on the relapse-free survival of BC, and none of the previous studies have concentrated on early relapse. Therefore, identifying autophagy-related mRNA signature could not only easily help oncologists classify the BC patients with a high risk of early relapse but also make more efficient therapeutic modalities at an earlier stage of a patient’s treatment (31).
In the present study, we conducted an autophagy-related 4-gene signature to predict the early relapse of BC patients and construct a nomogram for predicting the 1-, 2-, and 3-year RFS probability during clinical practice.
Materials and Methods
Data Source Collection
The messenger RNA (mRNA)-seq expression and clinicopathological characteristics of 1,025 BC patients were obtained from the TCGA program website (https://cancergenome.nih.gov/). Meanwhile, the autophagy-related genes were derived from the human autophagy database (HADb, http://www.autophagy.lu/). After excluding BC patients with incomplete clinicopathological medical records and patients initially diagnosed with metastasis, there were 785 BC patients included for further analysis. The data from the TCGA database were assigned as a training cohort (Figure 1). Moreover, two Gene Expression Omnibus (GEO) cohorts including the GSE42568 and GSE21653 datasets (detailed clinical information was summarized in Table S1) were obtained from the GEO database (https://www.ncbi.nlm.nih.gov/geo/) and assigned as the validation cohorts. All GEO datasets were produced by the Affymetrix HG-U133 plus 2.0. Raw microarray was normalized using Robust Multichip Average (32). When multiple probes were mapped to the same Entrez Gene ID, we used the mean value to represent its average expression level.
Ethics Approval
The protocol for this study was approved by Chongqing Medical University. Ethical approval was waived by the local Ethics Committee of the Chongqing Medical University in view of the retrospective nature of the study and all the procedures being performed were part of the routine care.
Identification of Autophagy-Related mRNA Signature for the Early Relapse of BC
Recurrence of BC patients was frequently occurred within 5 years after the initial treatment, while the first 2 years were the peak of recurrence (9). However, the definition of early relapse in BC patients is still ambiguous in recently published literature (33–36). Thus, early relapse in the present study was defined as the locoregional recurrence or distant metastasis within a short-term of 2 years follow-up after the initial primary resection. Samples in the training set were selected and divided into early relapse group and long-term survival group (no relapse at least 5 years follow-up). The calculations of differentially expressed genes (DEGs) between early relapse and long-term survival BC patients were conducted using the linear models for microarray data (LIMMA) method. The threshold for identification of DEGs was set as P value < 0.1. Besides, the LASSO Cox regression model (37) was used to select the most significantly relapse-associated mRNA of all the DEGs. A risk score model containing both coefficients and mRNA expression levels was established to generate the risk score for all BC patients in the training cohort. Based on the risk score, patients were divided into high-risk and low-risk groups with the median risk score as the cut-off point.
Gene Set Enrichment Analysis
The Gene Set Enrichment Analysis (GSEA, http://www.broadinstitute.org/gsea/index.jsp) was applied to evaluate differences between the low-risk and high-risk groups. Namely, Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were conducted to differentially expressed genes between these two groups. Normal P values <0.05 were regarded as statistically significantly enriched.
Validation Analysis
To further confirm the classification reliability and prognosis value of the 4-gene signature analyzed by TCGA, similar analyses were performed on GSE42568 and GSE216533 datasets to validate the prognostic significance of this autophagy-related signature.
Statistical Analysis
Survival differences between the low-risk and high-risk groups in each set were assessed by the Kaplan–Meier estimate and compared via the log-rank test. Baseline characteristics between low-risk and high-risk groups in each set were compared using the Pearson-chi square test (minimal expected value > 5). Multivariate Cox regression analysis and data stratification analysis were exploited to evaluate the independent prognostic significance of risk score and clinicopathological factors in predicting the RFS of BC patients. Time-dependent receiver-operating characteristic (ROC) analysis was used to investigate the prognostic and predictive accuracy of the signature. To access the probability of RFS survival in BC patients, a nomogram was subsequently developed based on the risk score and clinical features by using the “rms” R package. And the predictive feasibility of the nomogram was weighed by the Harrell concordance indexes (C-index) and calibration curves. All statistical analyses were performed with the use of R (version 4.0.3, www.r-project.org). All statistical tests were two-tailed, and P values < 0.05 were considered statistically significant.
Results
Identification of an Autophagy-Related Gene Signature for the Early Relapse of BC
Generally, we take the intersection of mRNA from the TCGA database with 222 autophagy genes in HADb. 46 differentially expressed autophagy-related genes were identified between the early relapse group and long-term survival group by using the “limma” package in the R software. These genes were subsequently included for LASSO analysis (Figure 2). Based on the LASSO analysis, four genes including the APOL1, HSPA8, SIRT1, and TP73, were regarded as the independent prognostic factor in early relapse BC patients. A gene-based prognostic model was further established to evaluate the RFS risk for each patient. The results are as follows: Risk score= (-0.209* status of APOL1) +(0.387* status of HSPA8) +(-1.073* status of SIRT1) +(-0.233* status of TP73). Thus, BC patients were divided into high-risk and low-risk groups with the median risk score as the cut-off point (Table 1).
Figure 2 (A) The heat map demonstrates forty-six differentially expressed mRNA in breast cancer among early relapse and long-term survival group both in training cohort; (B) LASSO coefficient profiles of the 4 early relapse-associated mRNA. A vertical line is drawn at the value chosen by 10-fold cross-validation; (C) X-tile analysis of the 4 selected GRGs.
Table 1 The demographic characteristics of breast cancer patients among high-risk and low-risk groups according the autophagy-related 4-gene signature.
The Prognostic Value of 4-Gene Signature in Training and Validating Cohorts
Among the high-risk group and low-risk group, the distribution of risk score and relapse status of BC patients were displayed. The results showed that the higher the risk score, the higher the morbidity rate was observed in the training group and two validating groups (Figure 3). Similarly, the Kaplan-Meier survival analysis demonstrated that the RFS of the high-risk group was significantly inferior to those of the low-risk group in the training group and two validating groups (log-rank p=0.002 in the training set, log-rank p=0.013 in the validation set I, log-rank p=0.003 in the validation set II, respectively). Moreover, the time-dependent ROC analyses at 1-, 2-, and 3-year were also conducted to evaluate the accuracy of the 4-gene classifier. In the training cohort, the AUC of 1-, 2-, and 3-year was 0.748, 0.696, and 0.651, respectively. Accordingly, a relatively promising AUC value was also observed in the validating sets (validating set I: the AUC of 1-, 2-, and 3-year was 0.611, 0.649, and 0.655, respectively; validating set II: the AUC of 1-, 2-, and 3-year was 0.524, 0.640, and 0.654, respectively).
Figure 3 (A) Distribution of risk score, time-dependent ROC curves at 1, 2, and 3 years and Kaplan–Meier survival analysis between patients at low and high risks of relapse in training cohort; (B) first external validation cohort; (C) second external validation cohort.
Establishment of a Predictive Nomogram
To access the independence and accuracy of the 4-gene signature in predicting the RFS of BC patients. The univariate and multivariate Cox analyses integrated with the clinicopathological characteristics were performed (Table 2). At univariate analysis, PR status, primary tumor size, regional lymph node status, and the 4-gene based signature were significantly associated with the early relapse of BC patients. During the multivariate analysis, larger tumor size (T2: HR=1.82, 95%CI: 0.51- 6.38, T3: HR=3.06, 95%CI: 0.50- 18.61, T4: HR=19.99, 95%CI: 3.88- 102.76, p=0.001) and high-risk BC patients derived from the 4-gene classifier (HR=5.73, 95%CI: 1.63- 20.16, p=0.006) were identified as the independent risk factors in promoting the early relapse of BC. PR status reached marginal significance (HR=2.34, 95%CI: 0.96- 5.73, p=0.063). A novel nomogram (Figure 4) was subsequently established with the three variables involvement (PR status, tumor size, and 4-gene signature). Optimally, the model contained a satisfying C-index of 0.766 (95%CI: 0.604-0.927). Moreover, three calibration curves for evaluating the accuracy of the predictive ability in short-term RFS were also performed via 1000 bootstrap repetitions (Figure 4). The curves (apparent, ideal, and bias-corrected lines) suggested a promising agreement in the training model.
Table 2 The univariate and multivariate Cox regression analysis in the early relapse of breast cancer in the training group.
Figure 4 (A) The nomogram for predicting the 1-, 2-, and 3- year relapse-free survival in breast cancer patients, based on the autophagy-related 4-gene signature selection and clinical factors. (B) The 1-, 2- and 3-year calibration curves were derived from the nomogram, respectively.
Gene Set Enrichment Analysis
The KEGG pathway analysis was conducted to discover the associated biological signaling pathway of 4 autophagy-related mRNA sets. Notably, differentially expressed genes between high-risk and low-risk groups were determined. Namely, the GSEA results indicated that the genes enriched in the high-risk group were related to the regulation of homologous recombination, N-glycan biosynthesis, oxidative phosphorylation, protein export, and RNA polymerase (Figure 5). On the contrary, in the low-risk group, the autophagy-related gene sets were involved in pathways related to dilated cardiomyopathy, cardiomyopathy HCM, phosphatidylinositol signaling system, proximal tubule bicarbonate reclamation, and vascular smooth muscle contraction (Figure 5).
Figure 5 Gene Set Enrichment Analysis (GSEA). (A) GSEA shows a significant enrichment of cancer-related pathways in the high-risk group based on the training cohort. (B) GSEA shows a significant enrichment of cancer-related pathways in the low-risk group based on the training cohort.
Discussion
To date, BC has become the leading malignancy among women worldwide (1, 2, 38) with a promising relatively higher 5- year survival rate, compared with other invasive cancers. Nevertheless, survivors can experience early recurrence with resistance to the initial treatments paralleled by highly invasive metastasis (9, 35). TNM stage and immunohistochemical indicators like ER, PR, Her-2, and Ki-67 index were frequently used to access the prognosis of BC patients. Chen et al. (12) determined that the late-stage (p< 0.001), poor differentiated grade (p = 0.002), PR-negative status (p = 0.014), and HER2-negative status (p = 0.033) were significant associated with the early relapse of BC. However, Huang et al. (11) determined that the cancer TNM stage was significantly associated with the early-relapse in BC patients, while clinical variables including age, tumor location, ER status, PR status, or HER2 status were not. In addition, a different survival pattern has been observed in BC patients with a relatively similar condition during clinical practice. These results indicated that genetic biomarkers also played a pivotal role in regulating tumor cell cycle progression and metastasis.
Regarding the gene signatures, previous works highlighted that the imbalance of cell proliferation and apoptosis, as well as autophagy regulation disorder, might also be attributed to the occurrence and development of BC. Namely, autophagy is a pivotal process in control of cell fate and significantly correlates with apoptosis via inactivating the mammalian target of rapamycin (mTOR) signaling pathway or directly activates the initiation step of autophagy by phosphorylating unc-51-like autophagy activating kinase 1 (ULK1). In terms of cancer initiation, autophagy is considered tumor-suppressive due to its cytoprotective role (23, 25, 39, 40). Notably, Marsh et al. (25) discovered autophagy could inhibit the metastasis of BC cells by accumulating the autophagy cargo receptor (ACR), neighbor to BRCA1. Moreover, recent two basic research demonstrated that MTA1 (metastasis-associated 1) (13) and long noncoding RNA H19 (14) were the regulators of autophagy in resistance to the endocrine therapy (tamoxifen). On the contrary, several studies (22, 41) demonstrated the autophagy was positively associated with the tumor growth, metastasis, and recurrence of BC, which could maintain the homeostasis and the survival of BC cells by removing dysfunctional or unnecessary substances. Therefore, autophagy is a powerful but double-edged sword, which had an essential impact on the prognosis of BC (41).
In the present study, 4 autophagy-related mRNA including APOL1, HSPA8, SIRT1, and TP73 were pivotal genes in the RFS of BC. Of these genes, APOL1 (apolipoprotein-L1) has been observed significantly associated with kidney disease, especially in terms of HIV-related chronic renal disease (42, 43). In regulating the proliferation and metastasis of cancer cells, recent studies speculated the phenotype of APOLs was involved in several cancers’ metastasis via the strong reduction of cellular adherence and increased in cell motility, together with an important reduction of the capacity for apoptosis (44–47). Besides, members of the heat-shock protein 70 (HSPA) family gained plenty of attention as a potential target for tumor therapy, which could promote cancer cell growth by different mechanisms (48–50). For instance, Rohde et al. (50) demonstrated the suppression role of HSPA in “HeLa” cells, namely, depletion of HSPA and HSPA2 arrested cancer cells in G2/M and G1, respectively. Regarding the SIRT1 (Sirtuin-1) mRNA, it significantly participated in gene regulation, genome stability maintenance, apoptosis, autophagy, and tumorigenesis (51). As recent studies reported, a downregulation of SIRT1 has already been described in gastric cancer (52) and breast cancer (53, 54). Zhang et al. demonstrated the activation of SIRT1 could suppress gastric cancer cells proliferation and metastasis via STAT3/MMP-13 signaling pathway (52). Meanwhile. Latifkar et al. (54) reported that inhibition of SIRT1 would impair the lysosomal function, resulting in the enhanced secretion of pro-tumorigenic exosomes which might reconstruct the extracellular matrix and enhance the invasive properties of cultured BC cells. Additionally, accumulating evidence has proved the dysfunction of TP73 (tumor protein p73) was associated with the proliferation and prognosis of different cancers (55). Notably, in vitro study, Sharif et al. (56) demonstrated that high expression levels of TP73 suppressed the proliferation of BC via enhanced autophagy and cell death. Alternatively, knockdown of TP73 decreased NAMPT (nicotinamide phosphoribosyl transferase) inhibition-induced autophagy and cell death.
Regarding the clinicopathological characteristics of BC patients, only primary tumor size was significantly associated with the early relapse of BC patients after stepwise multivariate analysis (adjust p=0.001). Interestingly, negative progesterone receptor status trended towards significantly increasing the risk of early relapse of BC patients (adjust p=0.063). Previously, compare to the prognostic value of estrogen receptor and HER2 status, progesterone receptor status was not so important. However, recent studies have demonstrated that progesterone receptor-negative tumors have generally been shown to have a poorer prognosis than progesterone receptor-positive tumors (57). Notably, evidence from one large population-based study, negative progesterone receptor status was associated with higher differentiation grade and subsequent recurrence score (58). Meanwhile, Zhang et al. demonstrated that the expression levels of progesterone receptor (cutoff point: 55%) played a pivotal role in predicting the relapse of hormone receptor-positive BC patients (59). However, the underlying mechanism and potential signaling pathway are still not clear but worth further investigation.
To our knowledge, we first discovered this autophagy-related 4-gene signature involved in the early relapse of BC. Based on the 4-gene signature, a risk score model was successfully established. And it was externally validated by two cohorts of GSE42568 and GSE21653, suggesting the favorable reproducibility of this signature in BC. However, the underlying molecular mechanism and signaling pathways of this signature are still inadequately clarified in BC. Nonetheless, the GSEA showed that the genes enriched in the high-risk group were related to the regulation of homologous recombination, N-glycan biosynthesis, oxidative phosphorylation, protein export, and RNA polymerase cancer-related signaling conduction. Alternatively, among the low-risk population, the autophagy-related gene sets were involved in pathways related to dilated cardiomyopathy, cardiomyopathy HCM, phosphatidylinositol signaling system, proximal tubule bicarbonate reclamation, and vascular smooth muscle contraction. Thus, further investigation of the underlying mechanisms may be meaningful. Additionally, constructing a convenient while reliable autophagy-related mRNA signature for identifying the risk biomarkers in promoting early relapse of BC would make up for the deficiency of clinicopathological classification, and further assist oncologists in formulating more efficient treatment modalities at an earlier stage of patients’ management. For this reason, we constructed a nomogram combined with two clinicopathological prognostic factors to predict the 1-, 2-, and 3-year RFS of BC patients in an effective quantitative approach. An optimal C-index of 0.766 was achieved which indicated the feasibility of identifying the high-risk BC patients with early relapse during clinical practice.
Indeed, there are some limitations in the current study needed to be mentioned and addressed in future works. First, this is a retrospective-designed study, and all BC samples were identified from the public database which inevitably weakened the findings we determined. Second, further basic research in our department and other medical centers is merited to external validate our conclusions and elucidate the functional roles of autophagy-related mRNA signature involved in the early relapse of BC. Moreover, with a significant improvement of overall survival in BC patients, longer follow-up (like 10 years) time could better help oncologists predict the clinical outcome in these patients. Last, the risk score model and nomogram can only be applied to predict early relapse in BC patients, and its prognostic role in the different molecular subtypes of BC warrants further evaluation.
Conclusion
In summary, our works demonstrate that 4 autophagy-related mRNA (APOL1, HSPA8, SIRT1, and TP73) are significantly associated with the early relapse of breast cancer during the postoperative follow-up. Based on the autophagy-related mRNA signature risk score classifier, good discrimination in identifying the BC patients with a high risk of early relapse is achieved. Moreover, we successfully establish and validate a utility nomogram derived from the risk scores combining tumor size and PR status for clinically predicting the 1-, 2-,3-year RFS probability in BC patients after initial surgical intervention. Future prospective clinical trials could verify the clinical significance of our autophagy-related mRNA signature in stratifying early relapse in BC patients postoperatively. The mechanisms and underlying signaling of the identified genes on the early relapse of BC are also needed to be further explored.
Author’s Note
The software application generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
(I) Conception and design: YM and YF. (II) Administrative support: SF and GY. (III) Provision of study materials or patients: YF, XW, and YM. (IV) Collection and assembly of data: YM and HL. (V) Data analysis and interpretation: YM and XW. (VI) Manuscript writing: All authors. All authors contributed to the article and approved the submitted version.
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.
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.
Acknowledgments
We thank the GEO and TCGA database for providing their platforms and contributors for their valuable datasets.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2022.824362/full#supplementary-material
References
1. Siegel RL, Miller KD, Jemal A. Cancer Statistics, 2020. CA Cancer J Clin (2020) 70(1):7–30. doi: 10.3322/caac.21590
2. Zhang S, Sun K, Zheng R, Zeng H, Wang S, Chen R, et al. Cancer Incidence and Mortality in China, 2015. J Natl Cancer Center (2020) 1:2–11. doi: 10.1016/j.jncc.2020.12.001
3. von Minckwitz G, Huang CS, Mano MS, Loibl S, Mamounas EP, Untch M, et al. Trastuzumab Emtansine for Residual Invasive HER2-Positive Breast Cancer. N Engl J Med (2019) 380(7):617–28. doi: 10.1056/NEJMoa1814017
4. Harbeck N, Gnant M. Breast Cancer. Lancet (2017) 389(10074):1134–50. doi: 10.1016/s0140-6736(16)31891-8
5. DeSantis C, Siegel R, Bandi P, Jemal A. Breast Cancer Statistics, 2011. CA Cancer J Clin (2011) 61(6):409–18. doi: 10.3322/caac.20134
6. Goff SL, Danforth DN. The Role of Immune Cells in Breast Tissue and Immunotherapy for the Treatment of Breast Cancer. Clin Breast Cancer (2021) 21(1):e63–73. doi: 10.1016/j.clbc.2020.06.011
7. Kennecke H, Yerushalmi R, Woods R, Cheang MC, Voduc D, Speers CH, et al. Metastatic Behavior of Breast Cancer Subtypes. J Clin Oncol (2010) 28(20):3271–7. doi: 10.1200/jco.2009.25.9820
8. Weigelt B, Peterse JL, van ‘t Veer LJ. Breast Cancer Metastasis: Markers and Models. Nat Rev Cancer (2005) 5(8):591–602. doi: 10.1038/nrc1670
9. Markopoulos CJ. Minimizing Early Relapse and Maximizing Treatment Outcomes in Hormone-Sensitive Postmenopausal Breast Cancer: Efficacy Review of AI Trials. Cancer Metastasis Rev (2010) 29(4):581–94. doi: 10.1007/s10555-010-9248-x
10. Yau C, Wang Y, Zhang Y, Foekens JA, Benz CC. Young Age, Increased Tumor Proliferation and FOXM1 Expression Predict Early Metastatic Relapse Only for Endocrine-Dependent Breast Cancers. Breast Cancer Res Treat (2011) 126(3):803–10. doi: 10.1007/s10549-011-1345-1
11. Huang MY, Wang YH, Chen FM, Lee SC, Fang WY, Cheng TL, et al. Multiple Genetic Polymorphisms of GSTP1 313ag, MDR1 3435CC, and MTHFR 677cc Highly Correlated With Early Relapse of Breast Cancer Patients in Taiwan. Ann Surg Oncol (2008) 15(3):872–80. doi: 10.1245/s10434-007-9719-7
12. Chen L, Romond E, Chokshi S, Saeed H, Hodskins J, Stevens M, et al. A Prognostic Model of Early Breast Cancer Relapse After Standard Adjuvant Therapy and Comparison With Metastatic Disease on Initial Presentation. Breast Cancer Res Treat (2012) 136(2):565–72. doi: 10.1007/s10549-012-2265-4
13. Lee MH, Koh D, Na H, Ka NL, Kim S, Kim HJ, et al. MTA1 Is a Novel Regulator of Autophagy That Induces Tamoxifen Resistance in Breast Cancer Cells. Autophagy (2018) 14(5):812–24. doi: 10.1080/15548627.2017.1388476
14. Wang J, Xie S, Yang J, Xiong H, Jia Y, Zhou Y, et al. The Long Noncoding RNA H19 Promotes Tamoxifen Resistance in Breast Cancer via Autophagy. J Hematol Oncol (2019) 12(1):81. doi: 10.1186/s13045-019-0747-0
15. Hoppe R, Achinger-Kawecka J, Winter S, Fritz P, Lo WY, Schroth W, et al. Increased Expression of MiR-126 and MiR-10a Predict Prolonged Relapse-Free Time of Primary Oestrogen Receptor-Positive Breast Cancer Following Tamoxifen Treatment. Eur J Cancer (2013) 49(17):3598–608. doi: 10.1016/j.ejca.2013.07.145
16. Filipits M, Rudas M, Heinzl H, Jakesz R, Kubista E, Lax S, et al. Low P27 Expression Predicts Early Relapse and Death in Postmenopausal Hormone Receptor-Positive Breast Cancer Patients Receiving Adjuvant Tamoxifen Therapy. Clin Cancer Res (2009) 15(18):5888–94. doi: 10.1158/1078-0432.Ccr-09-0728
17. Lindström LS, Karlsson E, Wilking UM, Johansson U, Hartman J, Lidbrink EK, et al. Clinically Used Breast Cancer Markers Such as Estrogen Receptor, Progesterone Receptor, and Human Epidermal Growth Factor Receptor 2 Are Unstable Throughout Tumor Progression. J Clin Oncol (2012) 30(21):2601–8. doi: 10.1200/jco.2011.37.2482
18. Hanahan D, Weinberg RA. Hallmarks of Cancer: The Next Generation. Cell (2011) 144(5):646–74. doi: 10.1016/j.cell.2011.02.013
19. Kroemer G. Autophagy: A Druggable Process That Is Deregulated in Aging and Human Disease. J Clin Invest (2015) 125(1):1–4. doi: 10.1172/jci78652
20. Yun CW, Lee SH. The Roles of Autophagy in Cancer. Int J Mol Sci (2018) 19(11):3466. doi: 10.3390/ijms19113466
21. Wang B, Mao JH, Wang BY, Wang LX, Wen HY, Xu LJ, et al. Exosomal MiR-1910-3p Promotes Proliferation, Metastasis, and Autophagy of Breast Cancer Cells by Targeting MTMR3 and Activating the NF-κb Signaling Pathway. Cancer Lett (2020) 489:87–99. doi: 10.1016/j.canlet.2020.05.038
22. Vera-Ramirez L, Vodnala SK, Nini R, Hunter KW, Green JE. Autophagy Promotes the Survival of Dormant Breast Cancer Cells and Metastatic Tumour Recurrence. Nat Commun (2018) 9(1):1944. doi: 10.1038/s41467-018-04070-6
23. Marsh T, Debnath J. Autophagy Suppresses Breast Cancer Metastasis by Degrading Nbr1. Autophagy (2020) 16(6):1164–5. doi: 10.1080/15548627.2020.1753001
24. Du WW, Yang W, Li X, Awan FM, Yang Z, Fang L, et al. A Circular RNA Circ-DNMT1 Enhances Breast Cancer Progression by Activating Autophagy. Oncogene (2018) 37(44):5829–42. doi: 10.1038/s41388-018-0369-y
25. Marsh T, Kenific CM, Suresh D, Gonzalez H, Shamir ER, Mei W, et al. Autophagic Degradation of NBR1 Restricts Metastatic Outgrowth During Mammary Tumor Progression. Dev Cell (2020) 52(5):591–604.e6. doi: 10.1016/j.devcel.2020.01.025
26. Nassar FJ, Nasr R, Talhouk R. MicroRNAs as Biomarkers for Early Breast Cancer Diagnosis, Prognosis and Therapy Prediction. Pharmacol Ther (2017) 172:34–49. doi: 10.1016/j.pharmthera.2016.11.012
27. Hamam R, Hamam D, Alsaleh KA, Kassem M, Zaher W, Alfayez M, et al. Circulating MicroRNAs in Breast Cancer: Novel Diagnostic and Prognostic Biomarkers. Cell Death Dis (2017) 8(9):e3045. doi: 10.1038/cddis.2017.440
28. Sun X, Hu Y, Wu J, Shi L, Zhu L, Xi PW, et al. RBMS2 Inhibits the Proliferation by Stabilizing P21 MRNA in Breast Cancer. J Exp Clin Cancer Res (2018) 37(1):298. doi: 10.1186/s13046-018-0968-z
29. Zhao W, Geng D, Li S, Chen Z, Sun M. LncRNA HOTAIR Influences Cell Growth, Migration, Invasion, and Apoptosis via the MiR-20a-5p/HMGA2 Axis in Breast Cancer. Cancer Med (2018) 7(3):842–55. doi: 10.1002/cam4.1353
30. Tokumaru Y, Asaoka M, Oshi M, Katsuta E, Yan L, Narayanan S, et al. High Expression of MicroRNA-143 Is Associated With Favorable Tumor Immune Microenvironment and Better Survival in Estrogen Receptor Positive Breast Cancer. Int J Mol Sci (2020) 21(9):842–55. doi: 10.3390/ijms21093213
31. Ma L, Liang Z, Zhou H, Qu L. Applications of RNA Indexes for Precision Oncology in Breast Cancer. Genomics Proteomics Bioinf (2018) 16(2):108–19. doi: 10.1016/j.gpb.2018.03.002
32. Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, et al. Exploration, Normalization, and Summaries of High Density Oligonucleotide Array Probe Level Data. Biostatistics (2003) 4(2):249–64. doi: 10.1093/biostatistics/4.2.249
33. Speers C, Chang SL, Pesch A, Ritter C, Olsen E, Chandler B, et al. A Signature That may be Predictive of Early Versus Late Recurrence After Radiation Treatment for Breast Cancer That may Inform the Biology of Early, Aggressive Recurrences. Int J Radiat Oncol Biol Phys (2020) 108(3):686–96. doi: 10.1016/j.ijrobp.2020.05.015
34. Ogiya A, Yamazaki K, Horii R, Shien T, Horimoto Y, Masuda N, et al. Post-Relapse Survival in Patients With the Early and Late Distant Recurrence in Estrogen Receptor-Positive HER2-Negative Breast Cancer. Breast Cancer (2017) 24(3):473–82. doi: 10.1007/s12282-016-0730-3
35. Geurts SM, de Vegt F, Siesling S, Flobbe K, Aben KK. Van Der Heiden-Van Der Loo MPattern of Follow-Up Care and Early Relapse Detection in Breast Cancer Patients. Breast Cancer Res Treat (2012) 136(3):859–68. doi: 10.1007/s10549-012-2297-9
36. Germain DR, Graham K, Glubrecht DD, Hugh JC, Mackey JR, Godbout R. DEAD Box 1: A Novel and Independent Prognostic Marker for Early Recurrence in Breast Cancer. Breast Cancer Res Treat (2011) 127(1):53–63. doi: 10.1007/s10549-010-0943-7
37. Tibshirani R. The Lasso Method for Variable Selection in the Cox Model. Stat Med (1997) 16(4):385–95. doi: 10.1002/(sici)1097-0258(19970228)16:4<385::aid-sim380>3.0.co;2-3
38. Jung KW, Won YJ, Hong S, Kong HJ, Lee ES. Prediction of Cancer Incidence and Mortality in Korea, 2020. Cancer Res Treat (2020) 52(2):351–8. doi: 10.4143/crt.2020.203
39. Tavera-Mendoza LE, Westerling T, Libby E, Marusyk A, Cato L, Cassani R, et al. Vitamin D Receptor Regulates Autophagy in the Normal Mammary Gland and in Luminal Breast Cancer Cells. Proc Natl Acad Sci USA (2017) 114(11):E2186–94. doi: 10.1073/pnas.1615015114
40. Niklaus NJ, Tokarchuk I, Zbinden M, Schläfli AM, Maycotte P, Tschan MP. The Multifaceted Functions of Autophagy in Breast Cancer Development and Treatment. Cells (2021) 10(6):1447. doi: 10.3390/cells10061447
41. Tyutyunyk-Massey L, Gewirtz DA. Roles of Autophagy in Breast Cancer Treatment: Target, Bystander or Benefactor. Semin Cancer Biol (2020) 66:155–62. doi: 10.1016/j.semcancer.2019.11.008
42. Goyal R, Singhal PC. APOL1 Risk Variants and the Development of HIV-Associated Nephropathy. FEBS J (2020) 288(19):5586–97. doi: 10.1111/febs.15677
43. Friedman DJ, Pollak MR. APOL1 and Kidney Disease: From Genetics to Biology. Annu Rev Physiol (2020) 82:323–42. doi: 10.1146/annurev-physiol-021119-034345
44. Ren L, Yi J, Li W, Zheng X, Liu J, Wang J, et al. Apolipoproteins and Cancer. Cancer Med (2019) 8(16):7032–43. doi: 10.1002/cam4.2587
45. Hu CA, Klopfer EI, Ray PE. Human Apolipoprotein L1 (Apol1) in Cancer and Chronic Kidney Disease. FEBS Lett (2012) 586(7):947–55. doi: 10.1016/j.febslet.2012.03.002
46. Vanhollebeke B, Pays E. The Function of Apolipoproteins L. Cell Mol Life Sci (2006) 63(17):1937–44. doi: 10.1007/s00018-006-6091-x
47. Pays E. The Function of Apolipoproteins L (Apols): Relevance for Kidney Disease, Neurotransmission Disorders, Cancer and Viral Infection. FEBS J (2021) 288(2):360–81. doi: 10.1111/febs.15444
48. Shevtsov M, Balogi Z, Khachatryan W, Gao H, Vígh L, Multhoff G. Membrane-Associated Heat Shock Proteins in Oncology: From Basic Research to New Theranostic Targets. Cells (2020) 9(5):1263. doi: 10.3390/cells9051263
49. Hao Y, Kacal M, Ouchida AT, Zhang B, Norberg E, Vakifahmetoglu-Norberg H. Targetome Analysis of Chaperone-Mediated Autophagy in Cancer Cells. Autophagy (2019) 15(9):1558–71. doi: 10.1080/15548627.2019.1586255
50. Rohde M, Daugaard M, Jensen MH, Helin K, Nylandsted J, Jäättelä M. Members of the Heat-Shock Protein 70 Family Promote Cancer Cell Growth by Distinct Mechanisms. Genes Dev (2005) 19(5):570–82. doi: 10.1101/gad.305405
51. Alves-Fernandes DK, Jasiulionis MG. The Role of SIRT1 on DNA Damage Response and Epigenetic Alterations in Cancer. Int J Mol Sci (2019) 20(13):3153. doi: 10.3390/ijms20133153
52. Zhang S, Yang Y, Huang S, Deng C, Zhou S, Yang J, et al. SIRT1 Inhibits Gastric Cancer Proliferation and Metastasis via STAT3/MMP-13 Signaling. J Cell Physiol (2019) 234(9):15395–406. doi: 10.1002/jcp.28186
53. Moore RL, Faller DV. SIRT1 Represses Estrogen-Signaling, Ligand-Independent Erα-Mediated Transcription, and Cell Proliferation in Estrogen-Responsive Breast Cells. J Endocrinol (2013) 216(3):273–85. doi: 10.1530/joe-12-0102
54. Latifkar A, Ling L, Hingorani A, Johansen E, Clement A, Zhang X, et al. Loss of Sirtuin 1 Alters the Secretome of Breast Cancer Cells by Impairing Lysosomal Integrity. Dev Cell (2019) 49(3):393–408.e7. doi: 10.1016/j.devcel.2019.03.011
55. Rodríguez N, Peláez A, Barderas R, Domínguez G. Clinical Implications of the Deregulated TP73 Isoforms Expression in Cancer. Clin Transl Oncol (2018) 20(7):827–36. doi: 10.1007/s12094-017-1802-3
56. Sharif T, Ahn DG, Liu RZ, Pringle E, Martell E, Dai C, et al. The NAD(+) Salvage Pathway Modulates Cancer Cell Viability via P73. Cell Death Differ (2016) 23(4):669–80. doi: 10.1038/cdd.2015.134
57. Cork DM, Lennard TW, Tyson-Capper AJ. Alternative Splicing and the Progesterone Receptor in Breast Cancer. Breast Cancer Res (2008) 10(3):207. doi: 10.1186/bcr2097
58. Wu SG, Zhang WW, Wang J, Lian CL, Sun JY, Chen YX, et al. Progesterone Receptor Status and Tumor Grade Predict the 21-Gene Recurrence Score of Invasive Lobular Breast Cancer. Biomark Med (2019) 13(12):1005–12. doi: 10.2217/bmm-2019-0209
Keywords: breast cancer, early relapse, autophagy, signature, nomogram
Citation: Min Y, Feng Y, Luo H, Hu D, Wei X, He D, Yin G and Fan S (2022) Identifying and Validating of an Autophagy-Related Gene Signature for the Prediction of Early Relapse in Breast Cancer. Front. Endocrinol. 13:824362. doi: 10.3389/fendo.2022.824362
Received: 29 November 2021; Accepted: 26 January 2022;
Published: 16 February 2022.
Edited by:
Claire Perks, University of Bristol, United KingdomReviewed by:
Shengchun Liu, First Affiliated Hospital of Chongqing Medical University, ChinaYu Jiang, Sichuan University, China
Copyright © 2022 Min, Feng, Luo, Hu, Wei, He, Yin and Fan. 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: Shenghao Fan, ZnNoMTAxM0Bob3NwaXRhbC5jcW11LmVkdS5jbg==; orcid.org/0000-0002-0929-2199
†These authors have contributed equally to this work and share the first authorship