Comprehensive analysis of autophagy-related prognostic genes in breast cancer

A Prognostic Risk Signature of Two Autophagy-Related Genes for Predicting Triple-Negative Breast Cancer Outcomes

Background

Triple-negative breast cancer (TNBC) is recognized as the most aggressive molecular subtype of breast cancer. Recent studies have highlighted the complex role of autophagy in the pathogenesis of TNBC.

Methods

In this study, we evaluated 18,330 genes, including 1111 autophagy-related genes, (ARGs), across 579 TNBC samples from online databases. Differentially expressed ARGs in TNBC were identified using high-throughput RNA-seq data from the Cancer Genome Atlas (TCGA). Prognostic factors were examined through Cox regression and multivariate Cox analyses, with predictive efficacy assessed using receiver operating characteristic (ROC) curves. A nomogram integrating the risk signature with clinicopathological factors, such as TNM stage, was developed. Immunohistochemical analysis of clinical samples was also conducted.

Results

EIF4EBP1 and NPAS3 were significantly correlated with prognostic outcomes in patients with TNBC. Multivariate Cox regression analysis demonstrated that the expression levels of these two genes were accurate predictors of disease progression in TNBC samples from TCGA and the GSE31519 dataset. The efficacy of this predictive model was validated using ROC curve analysis and calibration plots, confirming its ability to accurately estimate the 1-, 2-, and 3-year survival rates for individuals with TNBC. Additionally, EIF4EBP1 and NPAS3 expression influenced drug sensitivity in TNBC cell lines, with notably lower NPAS3 expression in TNBC tissues, particularly in Stage III cases. This study is the first to report NPAS3 expression in patients with TNBC.

Conclusion

The autophagy-related genes EIF4EBP1 and NPAS3 may serve as independent prognostic factors for individuals with TNBC.

Proteomic analysis of breast cancer based on immune subtypes

BACKGROUND

Immunotherapy is applied to breast cancer to resolve the limitations of survival gain in existing treatment modalities. With immunotherapy, a tumor can be classified into immune-inflamed, excluded and desert based on the distribution of immune cells. We assessed the clinicopathological features, each subtype's prognostic value and differentially expressed proteins between immune subtypes.

METHODS

Immune subtyping and proteomic analysis were performed on 56 breast cancer cases with neoadjuvant chemotherapy. The immune subtyping was based on the level of tumor-infiltrating lymphocytes (TILs) and Klintrup criteria. If the level of TILs was ≥ 10%, it was classified as immune-inflamed type without consideration of the Klintrup criteria. In cases of 1-9% TIL, Klintrup criteria 1-3 were classified as the immune-excluded subtype and Klintrup criteria not available (NA) was classified as NA. Cases of 1% TILs and Klintrup 0 were classified as the immune-desert subtype. Mass spectrometry was used to identify differentially expressed proteins in formalin-fixed paraffin-embedded biopsy tissues.

RESULTS

Of the 56 cases, 31 (55%) were immune-inflamed, 21 (38%) were immune-excluded, 2 (4%) were immune-desert and 2 (4%) were NA. Welch's t-test revealed two differentially expressed proteins between immune-inflamed and immune-excluded/desert subtypes. Coronin-1A was upregulated in immune-inflamed tumors (adjusted p = 0.008) and α-1-antitrypsin was upregulated in immune-excluded/desert tumors (adjusted p = 0.008). Titin was upregulated in pathologic complete response (pCR) than non-pCR among immune-inflamed tumors (adjusted p = 0.036).

CONCLUSIONS

Coronin-1A and α-1-antitrypsin were upregulated in immune-inflamed and immune-excluded/desert subtypes, respectively. Titin's elevated expression in pCR within the immune-inflamed subtype may indicate a favorable prognosis. Further studies involving large representative cohorts are necessary to validate these findings.

In Situ Insonation of Alkaline Buffer Containing Liposomes Leads to a Net Improvement of the Therapeutic Outcome in a Triple Negative Breast Cancer Murine Model

Breast cancer is characterized by an acidic micro-environment. Acidic extracellular pH gives cancer cells an evolutionary advantage, hence, neutralization of the extracellular pH has been considered as a potential therapeutic strategy. To address the issue of systemic pH alteration, an approach based on the targeted delivery of the buffering solution to the tumor region is investigated. The method relies on the use of low frequency ultrasound and sono-sensitive liposomes loaded with buffers at alkaline pH (LipHUS). After the i.v. injection of LipHUS, the application of ultrasound (US) at the sites of the pathology induces a local increase of pH that results highly effective in i) inhibiting primary tumor growth, ii) reducing tumor recurrence after surgery, and iii) suppressing metastases' formation. The experiments are carried out on a triple negative breast cancer mouse model. The results obtained demonstrate that localized and triggered release of bicarbonate or PBS buffer from sonosensitive liposomes represents an efficient therapeutic tool for treating triple-negative breast cancer. This approach holds promise for potential clinical translation.

Comprehensive evaluation of cell death-related genes as novel diagnostic biomarkers for breast cancer

Background

Breast cancer (BRCA) ranks first among cancers in terms of incidence and mortality rates in women, primarily owing to metastasis, chemo-resistance, and heterogeneity. To predict long-term prognosis and design novel therapies for BRCA, more sensitive markers need to be explored.

Methods

Data from 1089 BRCA patients were downloaded from TCGA database. Pearson's correlation analysis and univariate and multivariate Cox regression analyses were performed to assess the role of cell death-related genes (CDGs) in predicting BRCA prognosis. Kaplan-Meier survival curves were generated to compare the overall survival in the two subgroups. A nomogram was constructed using risk scores based on the five CDGs and other clinicopathological features. CCK-8, EdU incorporation, and colony formation assays were performed to verify the inhibitory effect of NFKBIA on BRCA cell proliferation. Transwell assay, flow cytometry, and immunohistochemistry analyses were performed to ascertain the biological function of NFKBIA.

Results

Five differentially expressed CDGs were detected among 156 CDGs. The risk score for each patient was then calculated based on the expression levels of the five CDGs. Distinct differences in immune infiltration, expression of immune-oncological targets, mutation status, and half-maximal inhibitory concentration values of some targeted drugs were observed between the high- and low-risk groups. Finally, in vitro cell experiments verified that NFKBIA overexpression suppresses the proliferation and migration of BRCA cells.

Conclusions

Our study revealed that some CDGs, especially NFKBIA, could serve as sensitive markers for predicting the prognosis of patients with BRCA and designing more personalized clinical therapies.

NOL12 as an Oncogenic Biomarker Promotes Hepatocellular Carcinoma Growth and Metastasis

Hepatocellular carcinoma (HCC) is a common malignancy with a poor prognosis worldwide. However, the pathogenesis of HCC remains poorly understood. In this study, we found that NOL12 was significantly overexpressed in independent HCC datasets from TCGA database. We confirmed that the expression level of NOL12 was upregulated in human HCC tissues and cell lines by RT-qPCR. High expression of NOL12 is associated with worse reduced overall survival (OS), high pathological grade, node metastasis, and advanced clinical stage in patients with HCC. Moreover, knockdown of NOL12 dramatically inhibits the proliferation and metastasis of HCC cells in vitro and in vivo. CIBERSORTx analysis revealed that twelve types of tumor-infiltrating immune cells (TICs) are correlated with NOL12 expression. The risk signature based on 8 NOL12-related genes is an independent prognostic factor for patients with HCC. The OS rate of patients in the low-risk score group was better than that in the high-risk score group. In addition, the total tumor mutation burden (TMB) in the high-risk score group increased significantly, and the risk scores could be used as an alternative indicator of immune checkpoint inhibitor (ICI) response. In conclusion, our findings indicated that NOL12 might be involved in the progression of HCC and can be used as a potential therapeutic target. Moreover, the NOL12-related risk signature may have predictive relevance with regard to ICI therapy.

In search of autophagy biomarkers in breast cancer: Receptor status and drug agnostic transcriptional changes during autophagy flux in cell lines

Autophagy drives drug resistance and drug-induced cancer cell cytotoxicity. Targeting the autophagy process could greatly improve chemotherapy outcomes. The discovery of specific inhibitors or activators has been hindered by challenges with reliably measuring autophagy levels in a clinical setting. We investigated drug-induced autophagy in breast cancer cell lines with differing ER/PR/Her2 receptor status by exposing them to known but divergent autophagy inducers each with a unique molecular target, tamoxifen, trastuzumab, bortezomib or rapamycin. Differential gene expression analysis from total RNA extracted during the earliest sign of autophagy flux showed both cell- and drug-specific changes. We analyzed the list of differentially expressed genes to find a common, cell- and drug-agnostic autophagy signature. Twelve mRNAs were significantly modulated by all the drugs and 11 were orthogonally verified with Q-RT-PCR (Klhl24, Hbp1, Crebrf, Ypel2, Fbxo32, Gdf15, Cdc25a, Ddit4, Psat1, Cd22, Ypel3). The drug agnostic mRNA signature was similarly induced by a mitochondrially targeted agent, MitoQ. In-silico analysis on the KM-plotter cancer database showed that the levels of these mRNAs are detectable in human samples and associated with breast cancer prognosis outcomes of Relapse-Free Survival in all patients (RSF), Overall Survival in all patients (OS), and Relapse-Free Survival in ER+ Patients (RSF ER+). High levels of Klhl24, Hbp1, Crebrf, Ypel2, CD22 and Ypel3 were correlated with better outcomes, whereas lower levels of Gdf15, Cdc25a, Ddit4 and Psat1 were associated with better prognosis in breast cancer patients. This gene signature uncovers candidate autophagy biomarkers that could be tested during preclinical and clinical studies to monitor the autophagy process.

Codon Usage Bias in Autophagy-Related Gene 13 in Eukaryotes: Uncovering the Genetic Divergence by the Interplay Between Nucleotides and Codon Usages

Synonymous codon usage bias is a universal characteristic of genomes across various organisms. Autophagy-related gene 13 (atg13) is one essential gene for autophagy initiation, yet the evolutionary trends of the atg13 gene at the usages of nucleotide and synonymous codon remains unexplored. According to phylogenetic analyses for the atg13 gene of 226 eukaryotic organisms at the nucleotide and amino acid levels, it is clear that their nucleotide usages exhibit more genetic information than their amino acid usages. Specifically, the overall nucleotide usage bias quantified by information entropy reflected that the usage biases at the first and second codon positions were stronger than those at the third position of the atg13 genes. Furthermore, the bias level of nucleotide ‘G’ usage is highest, while that of nucleotide ‘C’ usage is lowest in the atg13 genes. On top of that, genetic features represented by synonymous codon usage exhibits a species-specific pattern on the evolution of the atg13 genes to some extent. Interestingly, the codon usages of atg13 genes in the ancestor animals (Latimeria chalumnae, Petromyzon marinus, and Rhinatrema bivittatum) are strongly influenced by mutation pressure from nucleotide composition constraint. However, the distributions of nucleotide composition at different codon positions in the atg13 gene display that natural selection still dominates atg13 codon usages during organisms’ evolution.

Searching for essential genes and drug discovery in breast cancer and periodontitis via text mining and bioinformatics analysis

The primary purpose of the study was (1) to search for the essential genes associated with breast cancer and periodontitis, and (2) to identify candidate drugs targeted to these genes for expanding the potential drug indications. The genes related to both breast cancer and periodontitis were determined by text mining. Gene ontology and Kyoto Encyclopedia of Genes and Genomes analysis were performed on these genes, and protein-protein interaction analysis was carried out to export significant module genes. Drug-gene interaction database was employed for potential drug discovery. We identified 221 genes common to both breast cancer and periodontitis. The top six significant enrichment terms and 15 enriched signal pathways were selected. Among 24 significant genes demonstrated as a gene cluster, we found SERPINA1 and TF were significantly related to poor overall survival between the relatively high and low groups in patients. Using the final two genes, 12 drugs were identified that had potential therapeutic effects. SERPINA1 and TF were screened out as essential genes related to both breast cancer and periodontitis, targeting 12 candidate drugs that may expand drug indications. Drug discovery using text mining and analysis of different databases can promote the identification of existing drugs that have the potential of administration to improve treatment in breast cancer.

CDKN1C as a prognostic biomarker correlated with immune infiltrates and therapeutic responses in breast cancer patients

Breast cancer (BC) prognosis and therapeutic sensitivity could not be predicted efficiently. Previous evidence have shown the vital roles of CDKN1C in BC. Therefore, we aimed to construct a CDKN1C‐based model to accurately predicting overall survival (OS) and treatment responses in BC patients. In this study, 995 BC patients from The Cancer Genome Atlas database were selected. Kaplan‐Meier curve, Gene set enrichment and immune infiltrates analyses were executed. We developed a novel CDKN1C‐based nomogram to predict the OS, verified by the time‐dependent receiver operating characteristic curve, calibration curve and decision curve. Therapeutic response prediction was followed based on the low‐ and high‐nomogram score groups. Our results indicated that low‐CDKN1C expression was associated with shorter OS and lower proportion of naïve B cells, CD8 T cells, activated NK cells. The predictive accuracy of the nomogram for 5‐year OS was superior to the tumour‐node‐metastasis stage (area under the curve: 0.746 vs. 0.634, p < 0.001). The nomogram exhibited excellent predictive performance, calibration ability and clinical utility. Moreover, low‐risk patients were identified with stronger sensitivity to therapeutic agents. This tool can improve BC prognosis and therapeutic responses prediction, thus guiding individualized treatment decisions.

A Novel Transcription Factor-Based Prognostic Signature in Endometrial Cancer: Establishment and Validation

Background

Endometrial cancer (EC) is a common malignancy of the female reproductive system worldwide. Increasing evidence has suggested that many transcription factors are aberrantly expressed in various cancers. This study aimed to develop a transcription factor-based prognostic signature for EC.

Methods

Gene expression data and clinical data of EC patients were downloaded from The Cancer Genome Atlas (TCGA) database. Univariate Cox regression and Multivariate Cox regression analysis was used to construct a prognostic signature. Then, the efficacy of the prognostic signature was validated in a training cohort, testing cohort and then the entire cohort. Correlations between clinical features and the model were also analyzed, and a nomogram based on the multivariate Cox analysis was developed. Furthermore, we verified the effect of a key transcription factor, E2F1, on biological functions of EC in vitro.

Results

We developed a nine-transcription factor (MSX1, HOXB9, E2F1, DLX4, BNC2, DLX2, PDX1, POU3F2, and FOXP3) prognostic signature. Compared with those in the low-risk group, patients in the high-risk group had worse clinical outcomes. The area under the curve (AUC) of this prognostic signature for 5-year survival was 0.806 in the training cohort, 0.710 in the testing cohort and 0.761 in the entire cohort. Gene set enrichment analysis (GSEA) revealed a correlation between the prognostic signature and various cancer signaling pathways, and a hub transcription factor regulatory network was constructed. The prognostic signature was confirmed to have independent predictive value. Finally, a nomogram based on the prognostic signature and clinical independent prognostic factors was also established and performed well according to the calibration curves. Further, knockdown of E2F1 inhibited invasion and metastasis of EC cells.

Conclusion

Our study developed and validated a transcription factor-based prognostic signature that accurately predicts prognosis of EC patients. Moreover, E2F1 may represent a potential target for the treatment of EC.

Identification and development of an independent immune-related genes prognostic model for breast cancer

Background

Breast cancer is one of the main malignant tumors that threaten the lives of women, which has received more and more clinical attention worldwide. There are increasing evidences showing that the immune micro-environment of breast cancer (BC) seriously affects the clinical outcome. This study aims to explore the role of tumor immune genes in the prognosis of BC patients and construct an immune-related genes prognostic index.

Methods

The list of 2498 immune genes was obtained from ImmPort database. In addition, gene expression data and clinical characteristics data of BC patients were also obtained from the TCGA database. The prognostic correlation of the differential genes was analyzed through Survival package. Cox regression analysis was performed to analyze the prognostic effect of immune genes. According to the regression coefficients of prognostic immune genes in regression analysis, an immune risk scores model was established. Gene set enrichment analysis (GSEA) was performed to probe the biological correlation of immune gene scores. P < 0.05 was considered to be statistically significant.

Results

In total, 556 immune genes were differentially expressed between normal tissues and BC tissues (p < 0. 05). According to the univariate cox regression analysis, a total of 66 immune genes were statistically significant for survival risk, of which 30 were associated with overall survival (P < 0.05). Finally, a 15 immune genes risk scores model was established. All patients were divided into high- and low-groups. KM survival analysis revealed that high immune risk scores represented worse survival (p < 0.001). ROC curve indicated that the immune genes risk scores model had a good reliability in predicting prognosis (5-year OS, AUC = 0.752). The established risk model showed splendid AUC value in the validation dataset (3-year over survival (OS) AUC = 0.685, 5-year OS AUC = 0.717, P = 0.00048). Moreover, the immune risk signature was proved to be an independent prognostic factor for BC patients. Finally, it was found that 15 immune genes and risk scores had significant clinical correlations, and were involved in a variety of carcinogenic pathways.

Conclusion

In conclusion, our study provides a new perspective for the expression of immune genes in BC. The constructed model has potential value for the prognostic prediction of BC patients and may provide some references for the clinical precision immunotherapy of patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-08041-x.

Role of Hypoxia-Mediated Autophagy in Tumor Cell Death and Survival

Programmed cell death or type I apoptosis has been extensively studied and its contribution to the pathogenesis of disease is well established. However, autophagy functions together with apoptosis to determine the overall fate of the cell. The cross talk between this active self-destruction process and apoptosis is quite complex and contradictory as well, but it is unquestionably decisive for cell survival or cell death. Autophagy can promote tumor suppression but also tumor growth by inducing cancer-cell development and proliferation. In this review, we will discuss how autophagy reprograms tumor cells in the context of tumor hypoxic stress. We will illustrate how autophagy acts as both a suppressor and a driver of tumorigenesis through tuning survival in a context dependent manner. We also shed light on the relationship between autophagy and immune response in this complex regulation. A better understanding of the autophagy mechanisms and pathways will undoubtedly ameliorate the design of therapeutics aimed at targeting autophagy for future cancer immunotherapies.

Sanguisorba officinalis L. Suppresses Triple-Negative Breast Cancer Metastasis by Inhibiting Late-Phase Autophagy via Hif-1α/Caveolin-1 Signaling

Sanguisorba officinalis L. (SA) is a common herb for cancer treatment in the clinic, particularly during the consolidation phase to prevent occurrence or metastasis. Nevertheless, there are limited studies reporting the molecular mechanisms about its anti-metastatic function. It is well demonstrated that autophagy is one of the critical mechanisms accounting for metastasis and anti-cancer pharmacological actions of Chinese herbs. On the threshold, the regulatory effects and molecular mechanisms of SA in suppressing autophagy-related breast cancer metastasis were investigated in this study. In vitro findings demonstrated that SA potently suppressed the proliferation, colony formations well as metastasis process in triple-negative breast cancer. Network and biological analyses predicted that SA mainly targeted caveolin-1 (Cav-1) to induce anti-metastatic effects, and one of the core mechanisms was via regulation of autophagy. Further experiments—including western blotting, transmission electron microscopy, GFP-mRFP-LC3 immunofluorescence, and lysosomal-activity detection—validated SA as a potent late-stage autophagic inhibitor by increasing microtubule-associated light chain 3-II (LC3-II) conversion, decreasing acidic vesicular-organelle formation, and inducing lysosomal dysfunction even under conditions of either starvation or hypoxia. Furthermore, the anti-autophagic and anti-metastatic activity of SA was Cav-1-dependent. Specifically, Cav-1 knockdown significantly facilitated SA-mediated inhibition of autophagy and metastasis. Furthermore, hypoxia inducible factor-1α (Hif-1α) overexpression attenuated the SA-induced inhibitory activities on Cav-1, autophagy, and metastasis, indicating that SA may have inhibited autophagy-related metastasis via Hif-1α/Cav-1 signaling. In both mouse breast cancer xenograft and zebrafish xenotransplantation models, SA inhibited breast cancer growth and inhibited late-phase autophagy in vivo, which was accompanied by suppression of Hif-1α/Cav-1 signaling and the epithelial-mesenchymal transition. Overall, our findings not only indicate that SA acts as a novel late-phase autophagic inhibitor with anti-metastatic activities in triple-negative breast cancer, but also highlight Cav-1 as a regulator in controlling late-phase autophagic activity.

A novel autophagy-related lncRNA prognostic risk model for breast cancer

Long non-coding RNAs (lncRNAs) are well known as crucial regulators to breast cancer development and are implicated in controlling autophagy. LncRNAs are also emerging as valuable prognostic factors for breast cancer patients. It is critical to identify autophagy-related lncRNAs with prognostic value in breast cancer. In this study, we identified autophagy-related lncRNAs in breast cancer by constructing a co-expression network of autophagy-related mRNAs-lncRNAs from The Cancer Genome Atlas (TCGA). We evaluated the prognostic value of these autophagy-related lncRNAs by univariate and multivariate Cox proportional hazards analyses and eventually obtained a prognostic risk model consisting of 11 autophagy-related lncRNAs (U62317.4, LINC01016, LINC02166, C6orf99, LINC00992, BAIAP2-DT, AC245297.3, AC090912.1, Z68871.1, LINC00578 and LINC01871). The risk model was further validated as a novel independent prognostic factor for breast cancer patients based on the calculated risk score by Kaplan-Meier analysis, univariate and multivariate Cox regression analyses and time-dependent receiver operating characteristic (ROC) curve analysis. Moreover, based on the risk model, the low-risk and high-risk groups displayed different autophagy and oncogenic statues by principal component analysis (PCA) and Gene Set Enrichment Analysis (GSEA) functional annotation. Taken together, these findings suggested that the risk model of the 11 autophagy-related lncRNAs has significant prognostic value for breast cancer and might be autophagy-related therapeutic targets in clinical practice.

Molecular characterization of breast cancer: a potential novel immune-related lncRNAs signature

Background

Accumulating evidence has demonstrated that immune-related lncRNAs (IRLs) are commonly aberrantly expressed in breast cancer (BC). Thus, we aimed to establish an IRL-based tool to improve prognosis prediction in BC patients.

Methods

We obtained IRL expression profiles in large BC cohorts (N = 911) from The Cancer Genome Atlas (TCGA) database. Then, in light of the correlation between each IRL and recurrence-free survival (RFS), we screened prognostic IRL signatures to construct a novel RFS nomogram via a Cox regression model. Subsequently, the performance of the IRL-based model was evaluated through discrimination, calibration ability, risk stratification ability and decision curve analysis (DCA).

Results

A total of 52 IRLs were obtained from TCGA. Based on multivariate Cox regression analyses, four IRLs (A1BG-AS1, AC004477.3, AC004585.1 and AC004854.2) and two risk parameters (tumor subtype and TNM stage) were utilized as independent indicators to develop a novel prognostic model. In terms of predictive accuracy, the IRL-based model was distinctly superior to the TNM staging system (AUC: 0.728 VS 0.673, P = 0.010). DCA indicated that our nomogram had favorable clinical practicability. In addition, risk stratification analysis showed that the IRL-based tool efficiently divided BC patients into high- and low-risk groups (P < 0.001).

Conclusions

A novel IRL-based model was constructed to predict the risk of 5-year RFS in BC. Our model can improve the predictive power of the TNM staging system and identify high-risk patients with tumor recurrence to implement more appropriate treatment strategies.

Pathogenic Single Nucleotide Polymorphisms on Autophagy-Related Genes

In recent years, the study of single nucleotide polymorphisms (SNPs) has gained increasing importance in biomedical research, as they can either be at the molecular origin of a determined disorder or directly affect the efficiency of a given treatment. In this regard, sequence variations in genes involved in pro-survival cellular pathways are commonly associated with pathologies, as the alteration of these routes compromises cellular homeostasis. This is the case of autophagy, an evolutionarily conserved pathway that counteracts extracellular and intracellular stressors by mediating the turnover of cytosolic components through lysosomal degradation. Accordingly, autophagy dysregulation has been extensively described in a wide range of human pathologies, including cancer, neurodegeneration, or inflammatory alterations. Thus, it is not surprising that pathogenic gene variants in genes encoding crucial effectors of the autophagosome/lysosome axis are increasingly being identified. In this review, we present a comprehensive list of clinically relevant SNPs in autophagy-related genes, highlighting the scope and relevance of autophagy alterations in human disease.

Comprehensive analysis of autophagy-related prognostic genes in breast cancer

Accumulating evidence revealed that autophagy played vital roles in breast cancer (BC) progression. Thus, the aim of this study was to investigate the prognostic value of autophagy-related genes (ARGs) and develop a ARG-based model to evaluate 5-year overall survival (OS) in BC patients. We acquired ARG expression profiling in a large BC cohort (N = 1007) from The Cancer Genome Atlas (TCGA) database. The correlation between ARGs and OS was confirmed by the LASSO and Cox regression analyses. A predictive model was established based on independent prognostic variables. Thus, time-dependent receiver operating curve (ROC), calibration plot, decision curve and subgroup analysis were conducted to determine the predictive performance of ARG-based model. Four ARGs (ATG4A, IFNG, NRG1 and SERPINA1) were identified using the LASSO and multivariate Cox regression analyses. A ARG-based model was constructed based on the four ARGs and two clinicopathological risk factors (age and TNM stage), dividing patients into high-risk and low-risk groups. The 5-year OS of patients in the low-risk group was higher than that in the high-risk group (P < 0.0001). Time-dependent ROC at 5 years indicated that the four ARG-based tool had better prognostic accuracy than TNM stage in the training cohort (AUC: 0.731 vs 0.640, P < 0.01) and validation cohort (AUC: 0.804 vs 0.671, P < 0.01). The mutation frequencies of the four ARGs (ATG4A, IFNG, NRG1 and SERPINA1) were 0.9%, 2.8%, 8% and 1.3%, respectively. We built and verified a novel four ARG-based nomogram, a credible approach to predict 5-year OS in BC, which can assist oncologists in determining effective therapeutic strategies.