Abstract
Background:
Polyamine metabolism is closely associated with tumorigenesis, progression, and the tumor microenvironment (TME). This study aimed to determine whether polyamine metabolism-related genes (PMRGs) could predict prognosis and immunotherapy efficacy in Breast Cancer (BC).
Methods:
We conducted a comprehensive multi-omics analysis of PMRG expression profiles in BC. Consensus cluster analysis was used to identify PMRG expression subtypes in the METABRIC cohort. Univariate and multivariate Cox regression analyses were performed to identify independent prognostic genes, which were subsequently used to construct a predictive model for BC, along with a novel nomogram based on PMRGs. The model was validated using an independent cohort (GSE86166). Independent prognostic genes were further verified in BC tissues using quantitative real-time PCR (qRT-PCR), Semi-quantitative Western blot, and immunohistochemistry. Additionally, we analyzed the immune microenvironment and enriched pathways across different subtypes using multiple algorithms. Finally, the “oncoPredict” R package was used to assess potential drug sensitivities in high-risk and low-risk groups.
Results:
Seventeen polyamine metabolism genes were identified. PMRGs were abundantly expressed in tumor cells, with 12 survival-related genes being selected. In the METABRIC cohort, two PMRG expression subtypes were identified, with cancer- and immune-related pathways being more active in cluster B, which was associated with a worse prognosis. Six genes were used to construct a prognostic model through univariate and multivariate Cox regression analyses. The predictive performance of the polyamine metabolism model was validated by ROC curve analysis (training cohort: METABRIC, AUC3years=0.684; validation cohort: GSE86166, AUC3years=0.682). A nomogram combining risk scores and clinicopathological features was constructed. Decision Curve Analysis (DCA) demonstrated that the model could guide clinical treatment strategies. Four high-risk independent prognostic factors (OAZ1, SRM, SMOX, and SMS) were validated as being upregulated in breast cancer tissues. The model successfully stratified BC patients into high-risk and low-risk groups, with the high-risk group exhibiting poorer clinical outcomes. Functional analysis revealed significant differences in immune status and drug sensitivity between high-risk and low-risk groups.
Conclusions:
This study elucidated the biological characteristics of PMRG expression subtypes in BC, identifying a polyamine-related prognostic signature and four novel biomarkers to accurately predict prognosis and immunotherapy response in BC patients.
1 Introduction
Polyamines (PAs) are small polycationic alkyl amines, including putrescine, spermidine, and spermine. Polyamines fulfill important cellular functions not only in eukaryotes, but in virtually all organisms, including prokaryotes as well (). Polyamines exist in mammalian cells at millimolar concentrations. These molecules contain multiple amino groups, primarily produced by the decarboxylation of specific amino acids, and are essential for normal cell growth and development in eukaryotic organisms. PAs are involved in various cellular activities through interactions with negatively charged DNA, RNA, or proteins (, ), and their depletion leads to cell stagnation. They participate in stabilizing cell structures, binding nucleic acids, and biosynthesizing proteins (, ). Mutations in polyamine metabolism (PM) enzymes or transporters are linked to Snyder-Robinson syndrome and other diseases (–).
Breast cancer is the second most common cancer worldwide, affecting approximately 42,000 women annually and representing the leading cause of cancer-related deaths among females (). Due to the high heterogeneity of BC treatments are increasingly targeted toward subtypes, stages, and grades (). Conventional therapies for BC include surgery, chemotherapy, and radiotherapy (). Currently, immunotherapy has gained significant attention, encompassing immune checkpoint blockade, vaccines, immune-oncolytic virus drugs, and adoptive cell therapy (). The advancement of targeted and immunotherapy approaches has expanded treatment options, particularly for advanced cases. However, many patients fail to respond effectively, highlighting the need to identify novel biomarkers that can accurately predict immunotherapy response.
Many evidence strengthens the hypothesis that a rise in intracellular PA concentrations, mainly through an up-regulation of PA biosynthetic enzymes, is associated with increased cell proliferation and is usually linked to tumorigenesis (–). Dysregulated PM has been observed in various cancers, including breast, colorectal, prostate, skin, renal, and lung cancers (–). Elevated PA levels have also been detected in body fluids of cancer patients (, ). PM functions downstream of major oncogenic pathways (), and inhibitors of key biosynthetic enzymes in PM pathways can significantly impact tumor progression. Thus, PM is closely associated with tumorigenesis and represents a key target for cancer therapy. While targeted PMRG therapies offer a promising approach for BC, it remains unclear how PM affects the tumor microenvironment (TME) and its role in immunotherapy.
In this study, we identified 17 PMRGs from the literature and combined multi-omics data from the METABRIC, TCGA, and GEO datasets, including transcriptomics, single-cell sequencing, and copy number variation analyses. We classified two PMRG expression subtypes in the METABRIC cohort, uncovering significant differences in immune cell infiltration and pathways, suggesting that PM strongly influences TME characteristics. Using univariate and multivariate Cox analyses, we identified six prognostic genes to develop a risk score model for predicting patient survival, which was validated in the independent GSE86166 cohort. Furthermore, we confirmed the expression of key prognostic genes (OAZ1, SRM, SMOX, SMS) in clinical tissues through qRT-PCR, Western blot, and immunohistochemistry assays. The PMRG model revealed distinct immune features and immunotherapy responses across BC subtypes, enabling the prediction of clinical responses to chemotherapy. This is the first study to comprehensively investigate the role of PMRGs in BC patients from the perspectives of prognosis and immunotherapy response.
2 Materials and methods
2.1 Data Collection
The RNA sequencing expression profile data for BC samples were collected from the TCGA database (https://portal.gdc.cancer.gov/), including 113 normal and 1113 cancer cases. The training dataset was obtained from the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) database, comprising 1,992 cases (https://www.cbioportal.org/). The validation dataset (GSE86166) included 305 cases from the Gene Expression Omnibus (GEO) database (http://www.ncbi.nlm.nih.gov/geo/). The single-cell dataset GSE161529 was sourced from the Tumor Immune Single-Cell Hub (TISCH; http://tisch.comp-genomics.org) ().
2.2 Expression of genes related to polyamine metabolism and protein-protein interaction network
We collected 17 genes related to polyamine metabolism from the literature () (Supplementary Table S1) and compared their expression levels in tumor tissues and adjacent normal tissues using the TCGA database. A protein-protein interaction (PPI) network was constructed using the STRING platform (https://string-db.org/) to analyze the interconnections among these genes ().
2.3 Copy number variation analysis
Copy number variation (CNV) is common in cancers and often serves as a driving event. Chromosomal regions are frequently lost or gained in cancer patients, and CNVs are significantly associated with BC risk (). Breast cancer CNV data were downloaded from XENA (https://xena.ucsc.edu/) to investigate alterations in polyamine metabolism genes and their chromosomal locations.
2.4 Consensus clustering
Consensus cluster analysis was performed using the “ConsensusClusterPlus” R package. The K-means method was employed to identify distinct patterns related to polyamine metabolism gene expression. Principal component analysis (PCA) was conducted to analyze clustering of the two subtypes, and the “ggplot2” R package was used to validate clustering reliability with the uniform manifold approximation and projection (UMAP) method.
2.5 Functional enrichment analysis
We downloaded “c2.cp.kegg.symbols.gmt” and “c5.go.symbols.gmt” data from the MSigDB database (https://www.gsea-msigdb.org/gsea/msigdb) and performed gene set variation analysis (GSVA) using the “GSVA” R package ().
2.6 Development and validation of prognostic features based on polyamine metabolism-related genes
Univariate Cox regression analysis was conducted on the training METABRIC cohort to identify survival-associated genes. A multivariate Cox regression model was subsequently applied to determine key genes and calculate their coefficients for both the training cohort (METABRIC) and the validation cohort (GSE86166). Six polyamine metabolism-related genes closely associated with overall survival (OS) were identified. The risk score was calculated as follows:
The model’s predictive performance was assessed using Kaplan-Meier (KM) survival curves and receiver operating characteristic (ROC) curve analysis over time.
2.7 Construction and evaluation of prediction Nomogram
The “survival” R package was used to compare the prognosis of high-risk and low-risk groups in the METABRIC and GSE86166 BC cohorts. A prognostic nomogram based on the six independent prognostic genes was developed using the “rms” R package. This nomogram predicted 1-, 3-, and 5-year survival rates for BC patients. DCA was performed to evaluate the clinical net benefit ().
2.8 Relationship between risk score and immune cell infiltration
CIBERSORT and ssGSEA R scripts were used to quantify risk scores and immune cell infiltration (). CIBERSORT was applied to estimate the proportions of immune cell types in high-risk and low-risk groups, with the sum of immune cell type scores equaling 1 for each sample. Spearman correlation analysis was performed to investigate the relationship between risk scores and immune-infiltrating cells.
2.9 Immune characteristic analysis
Tumor mutation burden (TMB) for each patient in the METABRIC cohort was evaluated using the “ESTIMATE” software package (). Differences in TMB between high-risk and low-risk subgroups, including Stromal Score, Immune Score, and ESTIMATE Score, were analyzed.
2.10 Chemotherapy response prediction
The “oncoPredict” package was used to predict the sensitivity of therapeutic agents for the different subgroups ().
2.11 Patient samples
Human BC specimens and nearby nontumorigenic specimens were collected from 40 diagnosed BC patients without preoperative treatment during surgery from June 2018 to October 2022 at the First Affiliated Hospital of China Medical University. Among them, 20 samples were embedded in paraffin for immunohistochemical analysis; The remaining samples were all rapidly frozen in liquid nitrogen and then stored at -80°C for future use. This study was approved by the Ethics Committee of the First Affiliated Hospital of China Medical University.
2.12 qRT-PCR
Fifteen paired cancer and adjacent tissue samples from BC patients were collected at our institution and stored at -80°C. Total RNA was extracted using Trizol (Sigma) following the manufacturer’s protocol, and RNA quality was assessed using a Nanodrop (Thermo). cDNA was synthesized via reverse transcription using Hiscript QRT supermix for qPCR (Vazyme). QRT-PCR was performed using SYBR Green Mastermix (Vazyme), with β-actin serving as the reference gene. Expression levels were quantified using the 2-ΔΔCt method. Primer sequences for target genes are listed in Supplementary Table S2.
2.13 Western blot
Five paired cancer and adjacent tissue samples from BRCA patients were collected and stored at -80°C. Protein concentration was measured using Beyotime reagents (China). Twenty micrograms of protein were mixed with 5× SDS loading buffer, separated via 10% sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE), and transferred to a PVDF membrane (Millipore). The membrane was blocked with 3% BSA buffer in PBS at room temperature for 1 hour, followed by overnight incubation at 4°C with primary antibodies diluted in 3% BSA buffer. The next day, membranes were washed with TBST and incubated with specific secondary antibodies at room temperature for 2 hours. Protein expression 33was detected using an ECL chemiluminescence kit (Beyotime).
Antibodies used in this study:
Anti-OAZ1 (Ornithine Decarboxylase Antizyme 1; 1:1000; Rabbit; cat. no. ab223481; Abcam).
Anti-SMOX (Spermine Oxidase; 1:1000; Rabbit; cat. no. ab213631; Abcam).
Anti-SMS (Spermine Synthase; 1:1000; Rabbit; cat. no. ab241496; Abcam).
Anti-SRM (Spermidine Synthase; 1:1000; Rabbit; cat. no. ab156879; Abcam).
Anti-β-actin (Beta Actin; 1:3000; Rabbit; cat. no. AF7018; Affinity).
Goat Anti-Rabbit IgG (H+L) HRP (1:3000; cat. no. S0001; Affinity).
2.14 Immunohistochemistry
Paraffin-embedded tumor tissue sections were deparaffinized, rehydrated, and incubated overnight at 4°C with primary antibodies against OAZ1, SMOX, SRM, and SMS. The sections were then incubated with the corresponding anti-rabbit/mouse secondary antibodies (Zhongshan, China) at 37°C for 2 hours. Sections were treated with ABC-peroxidase and diaminobenzidine (DAB) (Zhongshan, China), counterstained with hematoxylin, and visualized using light microscopy.
2.15 Statistical analysis
The data were analyzed by R software (version 4.3.1, https://www.r-project.org/). R packages (ESTIMATE, ggplot2, GSVA, limma, survminer, and survival) were applied for data analysis and graph plotting. The median value of risk scores was treated as the cutoff value for the two subgroups. The qRT-PCR results were analyzed by GraphPad Prism (version 10.1.2). Student’s T-test was employed to compare the statistical differences between the two groups. The Kaplan– Meier method was performed for prognosis among groups. Multivariate COX analysis was used to screen prognostic related genes. The Pearson test was used for correlation analysis. A P value < 0.05 was considered statistically significant (*, p < 0.05; **, p < 0.01; ***, p < 0.001), and the false detection rate (FDR) q < 0.05 was used for multiple testing correction.
3 Result
3.1 Multiomics analysis of polyamine metabolism-related genes in Breast caner
The study was conducted according to the summary flowchart (Figure 1). A total of 17 polyamine metabolism-related genes (PMRGs) were identified from the literature. In the TCGA-BRCA cohort, 16 of these genes exhibited statistically significant differences in expression. Most polyamine metabolism genes (ATP13A2, AZIN1, AZIN2, OAZ1, OAZ2, OAZ3, AOC1, PAOX, SAT1, SMOX, SRM, SMS, and AGMAT) were highly expressed in BC In contrast, AMD1, ODC1, SAT2, and ARG1 showed lower expression in cancer tissues (Figure 2A).
Figure 1
Figure 2
We further analyzed the expression levels of PMRGs at the single-cell level using the GSE161529 dataset. After quality control, 32,168 cells were annotated, including B cells, CD4Tconv, CD8T, endothelial, epithelial, fibroblasts, monocytes/macrophages (Mono/Macro), NK cells, pericytes, and plasma cells (Figure 2B). PMRGs were expressed explicitly in BC cells. In PMRGs, OAZ1 and SAT1 are highly expressed across all cell types, especially in malignant cells and immune cells (Figures 2C, D). SMS and AMD1 are mainly expressed in malignant cells, epithelial cells and immune cells (CD8T, Mono/Macro). SRM is mainly expressed on epithelial cells, endothelial cells and fibroblasts, and is also expressed to a certain extent in immune cells. However, ARG1 and AOC1 are almost undetectable in the TME. These findings indicate that PMRGs play an important roles in BC and are closely linked to immunity (Supplementary Figures S1, S2).
To explore potential interactions among PMRGs, we performed a correlation analysis and constructed a protein-protein interaction (PPI) network (Figure 2E). ATP13A2 is a late endolysosomal transporter and is involved in polyamine transport (), but it did not directly interact with other proteins in our PPI network. ATP13A2 islinked to genes in the related network, indicating that intermediate proteins maybe involved in tumorigenesis. This result is consistent with that of Tang in oral squamous cell carcinoma (). The regulatory network provided a comprehensive view of the interconnections among the 17 PMRGs and their prognostic significance in BC. The network diagram (Figure 2F) illustrated the relationships among 12 survival-related genes, comprising 8 risk factors and 4 favorable factors. AZIN2 exhibited the highest number of negative correlations (blue lines), while ODC1 and SRM displayed the most positive correlations (pink lines), indicating their central regulatory roles within the PMRG network.
Copy number variation (CNV) analysis revealed frequent alterations of PMRGs in BC patients. Amplifications were most significant for OAZ3 and AZIN1, located on chromosomes 1 and 8, respectively. In contrast, ATP13A2 and SRM exhibited the most extensive copy number deletions, both located on chromosome 1 (Figures 2G, H).
3.2 12 polyamine metabolism-related genes were used for consistent clustering of breast cancer molecular subsets
To better understand the role of PMRGs in BC we performed consensus clustering on 12 prognostic PMRGs using the “ConsensusClusterPlus” R package. As shown in Figure 3A, when k=2, the cohort was effectively divided into two distinct subtypes. Principal component analysis (PCA) revealed significant differences between the two subtypes (Figure 3B).
Figure 3
Overall survival analysis revealed significant prognostic differences between the two subtypes (P<0.001) (Figure 3C). The PMRGs-A subtype demonstrated a significant survival advantage, while PMRGs-B exhibited poorer prognosis. Expression patterns of PMRGs in the two subgroups were visualized using boxplots (Figure 3D), showing that most PMRGs, such as SMS, SMOX, and SAT1, were highly expressed in PMRGs-B. Only ATP13A2, PAOX, and SAT2 displayed low expression, suggesting their potential as therapeutic targets.
3.3 Immune cell infiltration and pathways of two PMRG subtype clusters
A heatmap of PMRG expression and corresponding clinicopathological features for the two subtypes is shown in Figure 4A. To investigate immune cell infiltration differences between the subtypes, we used ssGSEA to visualize and compare 23 immune-infiltrating cell subtypes (Figure 4B). Significant differences were observed in immune cell infiltration between the two groups. In the high-risk PMRGs-B subgroup, there was a marked increase in immune cell infiltration, including activated B cells, activated CD4+ T cells, activated CD8+ T cells, activated dendritic cells, regulatory T cells, follicular helper cells, and type 1/17/2 helper T cells.
Figure 4
To explore the functional differences between the clusters, GSVA was applied to analyze KEGG pathway enrichment (Figure 4C). The PMRGs-B cluster, which exhibited poor prognosis, was enriched in pathways associated with the cell cycle, P53 signaling, and the pentose phosphate pathway— all of which are closely related to tumor development. GSEA enrichment analysis further confirmed that the high-risk PMRGs-B subtype was primarily enriched in autoimmune thyroid disease pathway (Figure 4D). There were significant differences in risk scores between the two PMRGs subtypes (Figure 4E), and alluvial maps (Figure 4F) showed changes in polyamine metabolism-related subtypes, risk scores, and life status.
These findings indicate that the two polyamine modification patterns display distinct immune infiltration characteristics and functional enrichment profiles, with PMRGs-B linked to poor prognosis and immune dysregulation.
3.4 Construction and validation of a well-performing prognostic signal related to polyamine metabolism
Univariate Cox regression analysis in the METABRIC cohort (Figure 5A) identified 12 PMRGs significantly correlated with survival (P<0.05). Of these, AGMAT, AMD1, ATP13A2, OAZ1, ODC1, SAT1, SMOX, SRM, and SMS were risk factors (HR > 1, P<0.05), whereas AZIN2, ATP13A2, PAOX, and SAT2 were favorable factors (HR < 1, P<0.05) for patient prognosis. To assess the clinical relevance of polyamine metabolism-related genes (PMRGs), we developed and validated a prognostic model using the independent cohort GSE86166. A multivariate Cox regression analysis was conducted to evaluate the effectiveness of the risk score model in predicting BC patient prognosis. Based on the median risk value, patients were divided into high- and low-risk groups. The final model consisted of six PMRGs, of which OAZ1, SMOX, SRM, and SMS were significantly correlated with overall survival (OS) and identified as risk factors (HR > 1, P<0.05) (Figure 5B).
Figure 5
Kaplan-Meier survival curves demonstrated that high-risk groups had significantly poorer outcomes in both the METABRIC training cohort and the GSE86166 validation cohort (Figures 5C, D). Time-dependent ROC curve analysis for OS at 1, 3, and 5 years showed strong predictive performance for the model (Figures 5E, F).
To account for clinicopathological factors, the PMRG-based risk score was integrated with clinical data to construct a nomogram (Figure 6A). The calibration plot confirmed the nomogram’s predictive accuracy (Figure 6B). DCA further demonstrated that the nomogram provides substantial benefits for predicting short- and long-term survival in BC patients (Figure 6C). These results highlight that the PMRG-based nomogram is a reliable and effective tool for predicting BC patient prognosis.
Figure 6
3.5 Identification of independent prognostic factors
Using univariate and multivariate Cox regression analyses, OAZ1, SMOX, SRM, and SMS were identified as independent prognostic factors. We further investigated the expression levels of these four genes and their diagnostic and prognostic relevance in BC tissues.
The expression levels of OAZ1, SMOX, SRM, and SMS were significantly elevated in cancer tissues compared to normal tissues (Supplementary Figures S3A–D), as confirmed by qRT-PCR (Figures 7A–D), Western blot (Figures 7E, F), and immunohistochemistry (Figure 8). Survival analysis revealed that patients with high expression of these four genes had a significantly poorer prognosis (Figures 7G–J).
Figure 7
Figure 8
These findings establish OAZ1, SMOX, SRM, and SMS as potential biomarkers and independent prognostic factors in BC providing critical insights for future clinical applications.
3.6 Gene set enrichment analysis and immune activity of different risk scores
The immune microenvironment plays a crucial role in tumor progression and response to immunotherapy. We used ssGSEA to analyze the expression of 22 immune-infiltrating cell types in association with risk scores. BC samples were ranked from low to high-risk scores, showing the proportion of different immune cells across the risk spectrum (Figure 9A). Macrophage M0 composition accounted for the largest proportion of immune cells in the high-risk group (Figure 9B), suggesting that macrophages M0 may contribute significantly to the poor prognosis of BC patients.
Figure 9
Further investigation of the regulatory relationships between immune cells (Figure 9C) revealed a significant negative correlation (cor = -0.42) between macrophages M0 and CD8+ T cells. The nine-gene signatures used to construct the PMRG-based model showed distinct expression patterns in high- and low-risk populations (Figure 9D) and were significantly associated with various immune cell infiltrations (Figure 9E).
We next analyzed the tumor microenvironment (TME) of BC patients in the high- and low-risk groups. By estimating expression-based scores, we obtained the stromal score, immune score, and ESTIMATE score for both groups (Figure 9F). The tumor mutation burden (TMB) was significantly higher in the high-risk group compared to the low-risk group, indicating that higher risk scores correlate with higher TMB.
Chemotherapy remains a cornerstone of cancer treatment. Using the “oncoPredict” R package, we explored potential differences in chemotherapy sensitivity between the high- and low-risk groups. High-risk patients exhibited greater sensitivity to common chemotherapy drugs, such as Cisplatin and Cyclophosphamide (Figures 10A, B). Additionally, small molecule drugs targeting PARP inhibitors (e.g., Talazoparib) and PI3K inhibitors (e.g., Alpelisib) were found to be more effective in high-risk patients (Figures 10C, D).
Figure 10
These findings suggest that BC patients in the high-risk group not only exhibit distinct immune microenvironment characteristics but also show promising sensitivity to chemotherapy and targeted therapies, offering potential clinical benefits for treatment optimization.
4 Discussion
Breast cancer is a heterogeneous and highly aggressive disease, ranking among the leading causes of global morbidity and mortality () and representing a significant disease burden for women worldwide (). Polyamines (PAs) play a critical role in the proliferation of breast cancer cells (). Although intracellular PA concentrations are tightly regulated, dysregulation frequently occurs in breast cancer cells, making polyamine metabolism (PM) a potential target for BC intervention. PM is closely associated with the tumor microenvironment (TME) and is involved in anti-tumor immunity (, ), with intricate crosstalk between PM and the TME. However, the role of PM in the BC TME remains underexplored (). Therefore, this study investigated the clinical significance of polyamine metabolism-related genes (PMRGs) in BC.
In this study, we analyzed PMRG expressions in BC at transcriptome and single-cell levels. Most PMRGs were overexpressed in BC samples, and a majority (12/17) were associated with survival. Given the heterogeneity of BC, patients were stratified based on the expression of 17 PMRGs using consensus clustering. Two distinct PMRG expression subtypes were identified, showing significant differences in prognosis and biological pathways. Patients with high-risk scores have the worse OS related to poor survival prognosis and high TMB than those with low risk scores and inhibitory TME status. These reveal that the risk score can offer an innovative approach to assessing the TME status and prognosis of breast cancer.
Our analysis revealed that the PMRG expression cluster B subtype had a poorer prognosis and was significantly enriched in tumor- and immune-related pathways. These included the cell cycle, natural killer cell-mediated cytotoxicity, T cell receptor signaling pathway, chemokine signaling pathway, P53 signaling pathway, and the pentose phosphate pathway (PPP). The P53 signaling pathway plays a pivotal role in regulating breast cancer progression. For example, activation of P53 by Salt-inducible kinase 1 (SIK1) promotes oxidative phosphorylation, thereby inhibiting aerobic glycolysis and suppressing cell proliferation in breast cancer (). Conversely, inhibition of the P53 signaling pathway contributes to breast cancer progression by increasing cell proliferation, migration, and invasion ().The pentose phosphate pathway (PPP) also plays an essential role in breast cancer metabolism, contributing to oxidative stress regulation, nucleotide synthesis, and maintenance of the REDOX state (, ). It has been implicated in various cancer cell processes, including proliferation, apoptosis, drug resistance, invasiveness, metastasis, and senescence (–47). By maintaining a high proliferative state, the PPP supports cancer cell viability (48, 49).
We developed a risk score model established by six genes (OAZ1, SMOX, SRM, SMS, ATP13A2, and PAOX) from the METABRIC and GEO database via multivariate COX correlation analysis, predicting the different prognosis of the high- and low- risk patients. Survival analysis indicates that this risk model can effectively distinguish high-risk and low-risk patients. To improve the accuracy of model prediction, this study combined prognostic risk scores with clinical characteristics to construct the nomogram for OS prediction. Through analysis of ROC curves and DCA, we found that Nomogram had a better predictive ability to predict BC prognosis than other indicators. We observed a significant positive correlation between risk scores and chemotherapy, as well as ER_IHC. This finding enhances the prediction accuracy and clinical significance of the risk model.
The occurrence and metastasis of tumors are closely related to the tumor microenvironment (TME) (50). In breast cancer TME, tumor-associated macrophages (TAMs) play a crucial role, constituting more than 50% of the tumor volume (51). TAMs are typically activated within tumors, where they exert tumor-promoting (52, 53) and immunosuppressive effects (54), ultimately leading to poor prognosis and chemotherapy resistance (55, 56).
Given the established relationship between polyamines and immunity, we analyzed immune cell infiltration in high- and low-risk groups using ssGSEA, CIBERSORT, and ESTIMATE algorithms. Immune cell infiltration was elevated in the high-risk group, particularly in macrophages M0, M1, T cells CD4 memory activated, and T cells follicular helper (Tfh) cells. Monocytes differentiate into three macrophage subtypes: non-activated M0 macrophages, pro-inflammatory M1 macrophages, and anti-inflammatory M2 macrophages. M1 macrophages release potent pro-inflammatory cytokines, such as TNF-α, IL-1, IL-6, IL-12, and iNOS, contributing to chronic inflammation and fibrous capsule formation. In contrast, M2 macrophages promote anti-inflammatory responses, tissue repair, and growth through the release of factors like IL-4, IL-10, and TGF-β (57–59).
Numerous previous studies have shown that M2 macrophages promote proliferation of breast cancer cells (54), renal tubular cells (60), colon cancer cells (61). Tfh cells provide essential help to B cells for effective antibody-mediated immune responses. In various solid organ tumor types of non-lymphocytic origin, their presence frequently coincides with a better prognosis. Existing studies confirm that Tfh cells is a key to the success of the immune checkpoint blockade (ICB) determinants and predictors (62).
Additionally, the high-risk group demonstrated a relatively favorable immunotherapy response based on the TIDE score. These results suggest that PMRG signatures may serve as predictors of immunotherapy responses. Given the dependency of tumor cells on polyamines and the critical physiological roles of polyamines in various immune cell types, targeting polyamine metabolic pathways may enhance immunotherapy efficacy (, , , 63).
This study first developed and validated the important biological function of PMRG signatures in determining the prognosis of BC patients. More importantly, candidate genes OAZ1, SMOX, SRM, and SMS were identified as independent prognostic factors. We further validated their expression levels in BC cells and explored their prognostic value using survival analysis.
OAZ1 is a key member of the ornithine decarboxylase enzyme family involved in polyamine metabolism (64). OAZ1 exhibits tumor inhibitory activity by affecting cell proliferation, apoptosis, and differentiation in oral cancer cell lines, leukemia, and non-small cell lung cancer (65–67). However, its role in breast cancer remains poorly understood.
SMOX (Spermine Oxidase) plays a major role in the catabolism of mammalian polyamines by oxidizing spermine to spermidine, producing reactive oxygen species (ROS) in the process (68). ROS induced by oxidative stress can lead to epithelial cell apoptosis but also cause DNA damage, thereby increasing the risk of tumorigenesis (69, 70). High SMOX expression has been implicated in gastric cancer (71), hepatocellular carcinoma (72), and colorectal cancer (73). SMOX is also associated with drug responses and cellular reactions to stress stimuli.
SRM (Spermidine Synthase) converts putrescine to spermidine in the polyamine biosynthesis pathway. SRM is overexpressed in prostate cancer and clear cell renal cell carcinoma (ccRCC), serving as a reliable biomarker and therapeutic target (74, 75). Overexpression of SRM may increase bladder cancer resistance to pirarubicin, while SRM knockdown improves chemotherapy efficacy (64). SRM also regulates the immune microenvironment, as its knockdown inhibits fibroblast proliferation (76). Recent studies have shown that targeting SRM can enhance the sensitivity of FGFR-mutant bladder cancer cells to erdafitinib treatment (77).
SMS (Spermine Synthase) is the only enzyme responsible for spermine (SPM) synthesis in mammalian cells (). As the final step in the polyamine biosynthesis pathway, SMS catalyzes the transfer of aminopropyl from decarboxylated S-adenosylmethionine (dcSAM) to spermidine (SPD) to produce spermine (78). Mutations in SMS cause Snyder-Robinson Syndrome (SRS), a rare X-linked recessive disorder characterized by intellectual disability, developmental delays, skeletal abnormalities, and seizures (78–80). Research on SMS in cancers is limited; however, existing studies indicate that SMS overexpression promotes colon and pancreatic cancer progression (81, 82). Targeting SMS with inhibitors reduces polyamine levels, thereby suppressing tumor cell proliferation (83).
In summary, this study, for the first time, developed and validated the critical biological function of polyamine metabolic gene markers in predicting BC patient prognosis. The nomogram based on this model offers a valuable tool for clinicians to develop personalized treatment plans for BC patients in clinical practice. Importantly, candidate genes OAZ1, SMOX, SRM, and SMS were identified as independent prognostic factors. Future research into the molecular mechanisms underlying these markers, along with prospective randomized clinical trials, will have significant clinical implications and provide a roadmap for precision medicine in breast cancer.
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.
Ethics statement
The studies involving humans were approved by the Ethics Committee of the First Affiliated Hospital of China Medical University (Approval Number: [2023] 643). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
RW: Conceptualization, Writing – original draft, Writing – review & editing. SC: Conceptualization, Data curation, Writing – original draft. QG: Formal Analysis, Investigation, Writing – original draft. YC: Methodology, Writing – original draft. XH: Methodology, Writing – original draft. FS: Resources, Writing – review & editing. CL: Software, Writing – review & editing. GZ: Writing – review & editing. BC: Writing – original draft, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. The research was supported by department of Science and Technology of Liaoning Province (2023JH2/101300048).
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.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2025.1613458/full#supplementary-material
Supplementary Figure 1-2Expression of SAT2, AZIN1, SMS, OAZ2, SRM, AMD1, ODC1, ATP13A2, SMOX, AOC1, OAZ3, AGMAT, ARG1, PAOX, AZIN2 in single cell data set GSE161529.
Supplementary Figure 3The expression of independent factors was upregulated in breast cancer tissues in TCGA dataset. (A)OAZ1, (B)SMOX, (C)SRM, and (D)SMS.
References
1
KrysenkoSWohllebenW. Polyamine and ethanolamine metabolism in bacteria as an important component of nitrogen assimilation for survival and pathogenicity. Med Sci (Basel). (2022) 10:40–72. doi: 10.3390/medsci10030040
2
PeggAE. Functions of polyamines in mammals. J Biol Chem. (2016) 291:14904–12. doi: 10.1074/jbc.R116.731661
3
MadeoFEisenbergTPietrocolaFKroemerG. Spermidine in health and disease. Science. (2018) 359:eaan2788. doi: 10.1126/science.aan2788
4
CervelliMPietropaoliSSignoreFAmendolaRMariottiniP. Polyamines metabolism and breast cancer: state of the art and perspectives. Breast Cancer Res Treat. (2014) 148:233–48. doi: 10.1007/s10549-014-3156-7
5
GernerEWBruckheimerECohenA. Cancer pharmacoprevention: Targeting polyamine metabolism to manage risk factors for colon cancer. J Biol Chem. (2018) 293:18770–8. doi: 10.1074/jbc.TM118.003343
6
LiCBrazillJMLiuSBelloCZhuYMorimotoMet al. Spermine synthase deficiency causes lysosomal dysfunction and oxidative stress in models of Snyder-Robinson syndrome. Nat Commun. (2017) 8:1257. doi: 10.1038/s41467-017-01289-7
7
BuppCPSchultzCRUhlKLRajasekaranSBachmannAS. Novel de novo pathogenic variant in the ODC1 gene in a girl with developmental delay, alopecia, and dysmorphic features. Am J Med Genet A. (2018) 176:2548–53. doi: 10.1002/ajmg.a.40523
8
van VeenSMartinSVan den HauteCBenoyVLyonsJVanhoutteRet al. ATP13A2 deficiency disrupts lysosomal polyamine export. Nature. (2020) 578:419–24. doi: 10.1038/s41586-020-1968-7
9
RodanLHAnyane-YeboaKChongKKlein Wassink-RuiterJSWilsonASmithLet al. Gain-of-function variants in the ODC1 gene cause a syndromic neurodevelopmental disorder associated with macrocephaly, alopecia, dysmorphic features, and neuroimaging abnormalities. Am J Med Genet A. (2018) 176:2554–60. doi: 10.1002/ajmg.a.60677
10
MomenimovahedZSalehiniyaH. Epidemiological characteristics of and risk factors for breast cancer in the world. Breast Cancer (Dove Med Press). (2019) 11:151–64. doi: 10.2147/bctt.S176070
11
HarbeckNGnantM. Breast cancer. Lancet. (2017) 389:1134–50. doi: 10.1016/s0140-6736(16)31891-8
12
MarusykAPolyakK. Tumor heterogeneity: causes and consequences. Biochim Biophys Acta. (2010) 1805:105–17. doi: 10.1016/j.bbcan.2009.11.002
13
AdamsSGatti-MaysMEKalinskyKKordeLASharonEAmiri-KordestaniLet al. Current landscape of immunotherapy in breast cancer: A review. JAMA Oncol. (2019) 5:1205–14. doi: 10.1001/jamaoncol.2018.7147
14
NowotarskiSLWosterPMCaseroRAJr. Polyamines and cancer: implications for chemotherapy and chemoprevention. Expert Rev Mol Med. (2013) 15:e3. doi: 10.1017/erm.2013.3
15
PeggAE. Mammalian polyamine metabolism and function. IUBMB Life. (2009) 61:880–94. doi: 10.1002/iub.230
16
AmendolaRCervelliMFratiniEPolticelliFSallustioDEMariottiniP. Spermine metabolism and anticancer therapy. Curr Cancer Drug Targets. (2009) 9:118–30. doi: 10.2174/156800909787580935
17
ManniAGroveRKunselmanSAldazCM. Involvement of the polyamine pathway in breast cancer progression. Cancer Lett. (1995) 92:49–57. doi: 10.1016/0304-3835(95)03763-m
18
CervelliMBellaviaGFratiniEAmendolaRPolticelliFBarbaMet al. Spermine oxidase (SMO) activity in breast tumor tissues and biochemical analysis of the anticancer spermine analogues BENSpm and CPENSpm. BMC Cancer. (2010) 10:555. doi: 10.1186/1471-2407-10-555
19
WallaceHM. The physiological role of the polyamines. Eur J Clin Invest. (2000) 30:1–3. doi: 10.1046/j.1365-2362.2000.00585.x
20
UppJRJr.SaydjariRTownsendCMJr.SinghPBarrancoSCThompsonJC. Polyamine levels and gastrin receptors in colon cancers. Ann Surg. (1988) 207:662–9. doi: 10.1097/00000658-198806000-00004
21
GuptaSAhmadNMarengoSRMacLennanGTGreenbergNMMukhtarH. Chemoprevention of prostate carcinogenesis by alpha-difluoromethylornithine in TRAMP mice. Cancer Res. (2000) 60:5125–33.
22
GilmourSK. Polyamines and nonmelanoma skin cancer. Toxicol Appl Pharmacol. (2007) 224:249–56. doi: 10.1016/j.taap.2006.11.023
23
ChenMNieZHuangDGaoYCaoHZhengLet al. Development of a polyamine gene expression score for predicting prognosis and treatment response in clear cell renal cell carcinoma. Front Immunol. (2022) 13:1048204. doi: 10.3389/fimmu.2022.1048204
24
WangNChaiMZhuLLiuJYuCHuangX. Development and validation of polyamines metabolism-associated gene signatures to predict prognosis and immunotherapy response in lung adenocarcinoma. Front Immunol. (2023) 14:1070953. doi: 10.3389/fimmu.2023.1070953
25
LiuRLinXLiZLiQBiK. Quantitative metabolomics for investigating the value of polyamines in the early diagnosis and therapy of colorectal cancer. Oncotarget. (2018) 9:4583–92. doi: 10.18632/oncotarget.22885
26
NakajimaTKatsumataKKuwabaraHSoyaREnomotoMIshizakiTet al. Urinary polyamine biomarker panels with machine-learning differentiated colorectal cancers, benign disease, and healthy controls. Int J Mol Sci. (2018) 19:756–70. doi: 10.3390/ijms19030756
27
CaseroRAJr.Murray StewartTPeggAE. Polyamine metabolism and cancer: treatments, challenges and opportunities. Nat Rev Cancer. (2018) 18:681–95. doi: 10.1038/s41568-018-0050-3
28
HanYWangYDongXSunDLiuZYueJet al. TISCH2: expanded datasets and new tools for single-cell transcriptome analyses of the tumor microenvironment. Nucleic Acids Res. (2023) 51:D1425–d31. doi: 10.1093/nar/gkac959
29
TangJWuXChengBLuY. Identification of a polyamine-related signature and six novel prognostic biomarkers in oral squamous cell carcinoma. Front Mol Biosci. (2023) 10:1073770. doi: 10.3389/fmolb.2023.1073770
30
JiangPGuSPanDFuJSahuAHuXet al. Signatures of T cell dysfunction and exclusion predict cancer immunotherapy response. Nat Med. (2018) 24:1550–8. doi: 10.1038/s41591-018-0136-1
31
KumaranMKrishnanPCassCEHubauxRLamWYasuiYet al. Breast cancer associated germline structural variants harboring small noncoding RNAs impact post-transcriptional gene regulation. Sci Rep. (2018) 8:7529. doi: 10.1038/s41598-018-25801-1
32
HänzelmannSCasteloRGuinneyJ. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinf. (2013) 14:7. doi: 10.1186/1471-2105-14-7
33
VickersAJCroninAMElkinEBGonenM. Extensions to decision curve analysis, a novel method for evaluating diagnostic tests, prediction models and molecular markers. BMC Med Inform Decis Mak. (2008) 8:53. doi: 10.1186/1472-6947-8-53
34
NewmanAMLiuCLGreenMRGentlesAJFengWXuYet al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. (2015) 12:453–7. doi: 10.1038/nmeth.3337
35
R Core T. A Language and Environment for Statistical Computing. R Foundation for Statistical Computing. (2023). Available online at: https://www.R-project.org/.
36
MaeserDGruenerRFHuangRS. oncoPredict: an R package for predicting in vivo or cancer patient drug response and biomarkers from cell line screening data. Brief Bioinform. (2021) 22:1–7. doi: 10.1093/bib/bbab260
37
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
38
HolbertCECullenMTCaseroRAJr.StewartTM. Polyamines in cancer: integrating organismal metabolism and antitumour immunity. Nat Rev Cancer. (2022) 22:467–80. doi: 10.1038/s41568-022-00473-2
39
NakkinaSPGittoSBPandeyVParikhJGGeertsDMaurerHCet al. Differential expression of polyamine pathways in human pancreatic tumor progression and effects of polyamine blockade on tumor microenvironment. Cancers (Basel). (2021) 13:6391–406. doi: 10.3390/cancers13246391
40
PonnusamyLManoharanR. Distinctive role of SIK1 and SIK3 isoforms in aerobic glycolysis and cell growth of breast cancer through the regulation of p53 and mTOR signaling pathways. Biochim Biophys Acta Mol Cell Res. (2021) 1868:118975. doi: 10.1016/j.bbamcr.2021.118975
41
GuanXShiAZouYSunMZhanYDongYet al. EZH2-Mediated microRNA-375 Upregulation Promotes Progression of Breast Cancer via the Inhibition of FOXO1 and the p53 Signaling Pathway. Front Genet. (2021) 12:633756. doi: 10.3389/fgene.2021.633756
42
TrachoothamDAlexandreJHuangP. Targeting cancer cells by ROS-mediated mechanisms: a radical therapeutic approach? Nat Rev Drug Discov. (2009) 8:579–91. doi: 10.1038/nrd2803
43
GaoWXuYChenTDuZLiuXHuZet al. Targeting oxidative pentose phosphate pathway prevents recurrence in mutant Kras colorectal carcinomas. PloS Biol. (2019) 17:e3000425. doi: 10.1371/journal.pbio.3000425
44
VizánPAlcarraz-VizánGDíaz-MoralliSSolovjevaONFrederiksWMCascanteM. Modulation of pentose phosphate pathway during cell cycle progression in human colon adenocarcinoma cell line HT29. Int J Cancer. (2009) 124:2789–96. doi: 10.1002/ijc.24262
45
SukhatmeVPChanB. Glycolytic cancer cells lacking 6-phosphogluconate dehydrogenase metabolize glucose to induce senescence. FEBS Lett. (2012) 586:2389–95. doi: 10.1016/j.febslet.2012.05.052
46
GeTYangJZhouSWangYLiYTongX. The role of the pentose phosphate pathway in diabetes and cancer. Front Endocrinol (Lausanne). (2020) 11:365. doi: 10.3389/fendo.2020.00365
47
RigantiCGazzanoEPolimeniMAldieriEGhigoD. The pentose phosphate pathway: an antioxidant defense and a crossroad in tumor cell fate. Free Radic Biol Med. (2012) 53:421–36. doi: 10.1016/j.freeradbiomed.2012.05.006
48
JinLZhouY. Crucial role of the pentose phosphate pathway in Malignant tumors. Oncol Lett. (2019) 17:4213–21. doi: 10.3892/ol.2019.10112
49
PatraKCHayN. The pentose phosphate pathway and cancer. Trends Biochem Sci. (2014) 39:347–54. doi: 10.1016/j.tibs.2014.06.005
50
HirataESahaiE. Tumor microenvironment and differential responses to therapy. Cold Spring Harb Perspect Med. (2017) 7:a026781. doi: 10.1101/cshperspect.a026781
51
VitaleIManicGCoussensLMKroemerGGalluzziL. Macrophages and metabolism in the tumor microenvironment. Cell Metab. (2019) 30:36–50. doi: 10.1016/j.cmet.2019.06.001
52
LuHClauserKRTamWLFröseJYeXEatonENet al. A breast cancer stem cell niche supported by juxtacrine signalling from monocytes and macrophages. Nat Cell Biol. (2014) 16:1105–17. doi: 10.1038/ncb3041
53
KomoharaYFujiwaraYOhnishiKTakeyaM. Tumor-associated macrophages: Potential therapeutic targets for anti-cancer therapy. Adv Drug Delivery Rev. (2016) 99:180–5. doi: 10.1016/j.addr.2015.11.009
54
WilliamsCBYehESSoloffAC. Tumor-associated macrophages: unwitting accomplices in breast cancer Malignancy. NPJ Breast Cancer. (2016) 2:15025–. doi: 10.1038/npjbcancer.2015.25
55
DeNardoDGBarretoJBAndreuPVasquezLTawfikDKolhatkarNet al. CD4(+) T cells regulate pulmonary metastasis of mammary carcinomas by enhancing protumor properties of macrophages. Cancer Cell. (2009) 16:91–102. doi: 10.1016/j.ccr.2009.06.018
56
XuanQJWangJXNandingAWangZPLiuHLianXet al. Tumor-associated macrophages are correlated with tamoxifen resistance in the postmenopausal breast cancer patients. Pathol Oncol Res. (2014) 20:619–24. doi: 10.1007/s12253-013-9740-z
57
OrihuelaRMcPhersonCAHarryGJ. Microglial M1/M2 polarization and metabolic states. Br J Pharmacol. (2016) 173:649–65. doi: 10.1111/bph.13139
58
AtriCGuerfaliFZLaouiniD. Role of human macrophage polarization in inflammation during infectious diseases. Int J Mol Sci. (2018) 19:1801–16. doi: 10.3390/ijms19061801
59
VinogradovSWarrenGWeiX. Macrophages associated with tumors as potential targets and therapeutic intermediates. Nanomed (Lond). (2014) 9:695–707. doi: 10.2217/nnm.14.13
60
HuangLChenWTanZHuangYGuXLiuLet al. Mrc1(+) macrophage-derived IGF1 mitigates crystal nephropathy by promoting renal tubule cell proliferation via the AKT/Rb signaling pathway. Theranostics. (2024) 14:1764–80. doi: 10.7150/thno.89174
61
CuiZLiWWangYZhaoMLiuKYangYet al. M2 macrophage-derived exosomal ferritin heavy chain promotes colon cancer cell proliferation. Biol Trace Elem Res. (2023) 201:3717–28. doi: 10.1007/s12011-022-03488-w
62
Gutiérrez-MeloNBaumjohannD. T follicular helper cells in cancer. Trends Cancer. (2023) 9:309–25. doi: 10.1016/j.trecan.2022.12.007
63
SamalKZhaoPKendzickyAYcoLPMcClungHGernerEet al. AMXT-1501, a novel polyamine transport inhibitor, synergizes with DFMO in inhibiting neuroblastoma cell proliferation by targeting both ornithine decarboxylase and polyamine transport. Int J Cancer. (2013) 133:1323–33. doi: 10.1002/ijc.28139
64
ZhuYPiaoCZhangZJiangYKongC. The potential role of c-MYC and polyamine metabolism in multiple drug resistance in bladder cancer investigated by metabonomics. Genomics. (2022) 114:125–37. doi: 10.1016/j.ygeno.2021.11.028
65
WangXJiangL. Effects of ornithine decarboxylase antizyme 1 on the proliferation and differentiation of human oral cancer cells. Int J Mol Med. (2014) 34:1606–12. doi: 10.3892/ijmm.2014.1961
66
WuBWangXMaWZhengWJiangL. Assay of OAZ1 mRNA levels in chronic myeloid leukemia combined with application of leukemia PCR array identified relevant gene changes affected by antizyme. Acta Haematol. (2014) 131:141–7. doi: 10.1159/000353406
67
Murray-StewartTRWosterPMCaseroRAJr. Targeting polyamine metabolism for cancer therapy and prevention. Biochem J. (2016) 473:2937–53. doi: 10.1042/bcj20160383
68
WangYDevereuxWWosterPMStewartTMHackerACaseroRAJr. Cloning and characterization of a human polyamine oxidase that is inducible by polyamine analogue exposure. Cancer Res. (2001) 61:5370–3.
69
ChaturvediRAsimMRomero-GalloJBarryDPHogeSde SabletTet al. Spermine oxidase mediates the gastric cancer risk associated with Helicobacter pylori CagA. Gastroenterology. (2011) 141:1696–708.e1-2. doi: 10.1053/j.gastro.2011.07.045
70
XuHChaturvediRChengYBussiereFIAsimMYaoMDet al. Spermine oxidation induced by Helicobacter pylori results in apoptosis and DNA damage: implications for gastric carcinogenesis. Cancer Res. (2004) 64:8521–5. doi: 10.1158/0008-5472.Can-04-3511
71
ChaturvediRde SabletTAsimMPiazueloMBBarryDPVerriereTGet al. Increased Helicobacter pylori-associated gastric cancer risk in the Andean region of Colombia is mediated by spermine oxidase. Oncogene. (2015) 34:3429–40. doi: 10.1038/onc.2014.273
72
HuTSunDZhangJXueRJanssenHLATangWet al. Spermine oxidase is upregulated and promotes tumor growth in hepatocellular carcinoma. Hepatol Res. (2018) 48:967–77. doi: 10.1111/hepr.13206
73
KimSKimDRohSHongIKimHAhnTSet al. Expression of spermine oxidase is associated with colorectal carcinogenesis and prognosis of patients. Biomedicines. (2022) 10:626–37. doi: 10.3390/biomedicines10030626
74
UrabeFKosakaNSawaYItoKKimuraTEgawaSet al. The miR-1908/SRM regulatory axis contributes to extracellular vesicle secretion in prostate cancer. Cancer Sci. (2020) 111:3258–67. doi: 10.1111/cas.14535
75
KwonOKHaYSNaAYChunSYKwonTGLeeJNet al. Identification of novel prognosis and prediction markers in advanced prostate cancer tissues based on quantitative proteomics. Cancer Genomics Proteomics. (2020) 17:195–208. doi: 10.21873/cgp.20180
76
ForshellTPRimpiSNilssonJA. Chemoprevention of B-cell lymphomas by inhibition of the Myc target spermidine synthase. Cancer Prev Res (Phila). (2010) 3:140–7. doi: 10.1158/1940-6207.Capr-09-0166
77
YuYGaoXZhaoHSunJWangMXiongXet al. A genome-wide synthetic lethal screen identifies spermidine synthase as a target to enhance erdafitinib efficacy in FGFR-mutant bladder cancer. Cancer Res. (2025) 85:2288–301. doi: 10.1158/0008-5472.Can-24-3217
78
PeggAEMichaelAJ. Spermine synthase. Cell Mol Life Sci. (2010) 67:113–21. doi: 10.1007/s00018-009-0165-5
79
CasonALIkeguchiYSkinnerCWoodTCHoldenKRLubsHAet al. X-linked spermine synthase gene (SMS) defect: the first polyamine deficiency syndrome. Eur J Hum Genet. (2003) 11:937–44. doi: 10.1038/sj.ejhg.5201072
80
StewartTMFoleyJRHolbertCEKhomutovMRastkariNTaoXet al. Difluoromethylornithine rebalances aberrant polyamine ratios in Snyder-Robinson syndrome. EMBO Mol Med. (2023) 15:e17833. doi: 10.15252/emmm.202317833
81
GuoYYeQDengPCaoYHeDZhouZet al. Spermine synthase and MYC cooperate to maintain colorectal cancer cell survival by repressing Bim expression. Nat Commun. (2020) 11:3243. doi: 10.1038/s41467-020-17067-x
82
GuoZZhangXLinCHuangYZhongYGuoHet al. METTL3-IGF2BP3-axis mediates the proliferation and migration of pancreatic cancer by regulating spermine synthase m6A modification. Front Oncol. (2022) 12:962204. doi: 10.3389/fonc.2022.962204
83
LiLLiJRaoJNLiMBassBLWangJY. Inhibition of polyamine synthesis induces p53 gene expression but not apoptosis. Am J Physiol. (1999) 276:C946–54. doi: 10.1152/ajpcell.1999.276.4.C946
Summary
Keywords
polyamine metabolism-related genes, breast cancer, multi-omics, METABRIC, prognostic, tumor microenvironment
Citation
Wang R, Cai S, Gao Q, Chen Y, Han X, Shang F, Liang C, Zhu G and Chen B (2025) Polyamine metabolism related gene index prediction of prognosis and immunotherapy response in breast cancer. Front. Oncol. 15:1613458. doi: 10.3389/fonc.2025.1613458
Received
17 April 2025
Accepted
16 July 2025
Published
31 July 2025
Volume
15 - 2025
Edited by
Ali Ata Moazzami, Swedish University of Agricultural Sciences, Sweden
Reviewed by
Sergii Krysenko, Valent BioSciences LLC, United States
Cassandra E Holbert, Loyola University Maryland, United States
Updates
Copyright
© 2025 Wang, Cai, Gao, Chen, Han, Shang, Liang, Zhu and Chen.
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: Guolian Zhu, 490052343@qq.com; Bo Chen, bochen@cmu.edu.cn
†These authors have contributed equally to this work
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.