Cuproptosis engages in c-Myc-mediated breast cancer stemness
Journal of Translational Medicine volume 21, Article number: 409 (2023)
Intra-tumoral heterogeneity (ITH) is a distinguished hallmark of cancer, and cancer stem cells (CSCs) contribute to this malignant characteristic. Therefore, it is of great significance to investigate and even target the regulatory factors driving intra-tumoral stemness. c-Myc is a vital oncogene frequently overexpressed or amplified in various cancer types, including breast cancer. Our previous study indicated its potential association with breast cancer stem cell (BCSC) biomarkers.
In this research, we performed immunohistochemical (IHC) staining on sixty breast cancer surgical specimens for c-Myc, CD44, CD24, CD133 and ALDH1A1. Then, we analyzed transcriptomic atlas of 1533 patients with breast cancer from public database.
IHC staining indicated the positive correlation between c-Myc and BCSC phenotype. Then, we used bioinformatic analysis to interrogate transcriptomics data of 1533 breast cancer specimens and identified an intriguing link among c-Myc, cancer stemness and copper-induced cell death (also known as “cuproptosis”). We screened out cuproptosis-related characteristics that predicts poor clinical outcomes and found that the pro-tumoral cuproptosis-based features were putatively enriched in MYC-targets and showed a significantly positive correlation with cancer stemness.
In addition to previous reports on its oncogenic roles, c-Myc showed significant correlation to stemness phenotype and copper-induced cell toxicity in breast cancer tissues. Moreover, transcriptomics data demonstrated that pro-tumoral cuproptosis biomarkers had putative positive association with cancer stemness. This research combined clinical samples with large-scale bioinformatic analysis, covered description and deduction, bridged classic oncogenic mechanisms to innovative opportunities, and inspired the development of copper-based nanomaterials in targeting highly heterogeneous tumors.
Breast cancer is the most common malignant tumor in women, leading to over 684,000 deaths according to latest report by CA Cancer J Clin [1, 2]. Though long-term prognosis of breast cancer is much better than other highly malignant tumors, such as melanoma, challenges including recurrence, distant metastasis and drug resistance still limit therapeutic effects and survival. Cancer stem cells (CSCs) are a self-renewal population critical to tumorigenesis, intra-tumoral heterogeneity, therapeutic resistance and recurrence . Plasticity, quiescence, differentiative potential, and even its low mutational state means that CSCs have considerable research value [4,5,6].
High expression of CD44, CD133, ALDH and low expression of CD24 are well-acknowledged biomarkers in maintaining the stemness of CSC in solid tumors, including breast cancer and lung cancer [7,8,9]. For breast cancer cases, the CD44 high/CD24 low phenotype exacerbates aggressiveness, accelerates antitumor drug resistance and facilitates tumor progression [10, 11]. Moreover, even CSCs present varying phenotypes and spatial characteristics. Our previous research investigated clinical specimens of triple-negative breast cancer (TNBC), and termed this observation intra-tumoral stemness heterogeneity (ITSH) and concluded that ITSH is a promising biomarker of poor clinical outcomes . CSCs differentiate into heterogeneous cancer cell clusters by asymmetrical division, decrease energy consumption in nutrition-deprived tumor microenvironment (TME), or maintain in a low-mutant state to survive immune attack . However, molecular progenitors and driving mechanisms underlying cancer stemness are still elusive.
c-Myc elevation or hyperactivation exists in various cancers . c-Myc exerts pro-tumoral effects via regulating cell cycle and activating downstream proliferative signaling pathways [14, 15]. In vitro experiments showed that c-Myc overexpression transforms normal breast epithelial cells into breast cancer cells and induces epithelial-mesenchymal transition (EMT) process . For metastasis-initiating cells, cancer cell colonization is accompanied by increased proliferation and elevated c-Myc expression . Clinically, c-Myc amplification is associated with greater benefit from trastuzumab for patients with HER2-positive breast cancer . Our previous bioinformatic analysis proposed that c-Myc dysfunction is a potential regulator of cancer stemness and participates in rewiring intracellular molecular interactions and metabolisms [19, 20]. Therefore, in this research, we aimed to investigate the mechanisms between c-Myc, breast cancer stemness and clinical outcomes of patients.
Cuproptosis is a unique form of cell death induced by copper overloading that targets the tricarboxylic acid (TCA) cycle . Some studies reported its role in cancer progression [22,23,24], but it remains to be seen whether this new concept offers druggable targets for use in clinical practice. In this study, we used transcriptomic profiling from large-scale databases of breast cancer patients, conducted rounds of screening and identified that cuproptosis was associated with poor clinical outcomes of breast cancer and cuproptosis engaged in c-Myc mediated breast cancer stemness. Thus, we proposed that cuproptosis-based risk scoring evaluation is conducive to predicting clinical outcomes of breast cancer. This research bridged the classic oncogene c-Myc and the new concept cuproptosis, proposed for the first time that cuproptosis may engaged in c-Myc-mediated breast cancer stemness, and identified a new cuproptosis-based scoring system that could be used for clinical applications.
Materials and methods
Slices of breast cancer tissue and corresponding lymph nodes with cancer cell infiltration were collected from sixty patients after breast cancer surgery. Histological grading was performed based on the criteria of the World Health Organization (2019 World Health Organization classification of tumours of the breast) . Pathologic staging was evaluated by pathologists according to the current International Union against Cancer Tumor Lymph Node Metastasis.
IHC staining of breast tumor tissue samples was performed using antibodies against c-Myc, CD44, CD24, CD133 and ALDH1A1. Pathologic staging was determined by the current International Union against Cancer Tumor Lymph Node metastasis classification. Scanning files or images on paraffin-embedded specimens were collected by digital tissue scanner or imaging system, and the tissue measurement area was performed by pathologists. The number of weak, medium and strong positive cells in the measurement area was analyzed and calculated respectively (negative without staining, 0 points; weak positive light yellow, count 1 point; medium positive brown yellow, count 2 points; strong positive tan count 3 points), total cell count, positive cumulative optical density IOD value, positive pixel area, tissue area mm. The following results were calculated to reflect the degree of positivity .
Data sources and normalization
Clinical information and gene expression profiling data (fragments per kilobase million, FPKM) of patients with Breast invasive carcinoma (BRCA) were downloaded from Gene Expression Omnibus (GEO) database and The Cancer Genome Atlas (TCGA) database. Three GEO BRCA cohorts (GSE42568, GSE48390, and GSE88770) and TCGA cohorts were obtained for subsequent analyses (Table 1) [27,28,29,30]. The FPKM values of TCGA-Breast invasive carcinoma (BRCA) were transformed into transcripts per kilobase million (TPM), as previously described, and were believed to be identical to those from microarrays. Four datasets were combined, and batch effect removal were performed by applying the “combat” algorithm of the SVA package, and the filtered expression matrix was normalized using Seurat’s NormalizeData function [31,32,33,34]. After data normalization, highly variable genes were identified and used for the following principal component analysis (PCA). A total of 1533 patients with BRCA from four data sets were enrolled in this study, including 104 patients from GSE42568, 81 patients from GSE48390, 117 patients from GSE88770, and 1091 patients from TCGA database [27,28,29,30]. The clinical variables included age, gender, grade, tumor stage, TNM Stage, HER2 mutation, ER mutation, follow-up time, and survival status.
Nonnegative matrix factorization (NMF) clustering analysis of cuproptosis biomarkers
Ten cuproptosis biomarkers (FDX1, LIAS, LIPT1, DLD, DLAT, PDHA1, PDHB, MTF1, GLS, CDKN2A) were retrieved from previous publications . R package “NMF” was employed to divide all samples into two CuRGcluster 1 and CuRGcluster 2 according to expression of cuproptosis biomarkers . This clustering was performed based on the following criteria. First, the cumulative distribution function (CDF) curve increased gradually and smoothly. Second, no groups had a small sample size. Lastly, after clustering, the intra-group correlation increased, while the inter-group correlation decreased.
Relationship between molecular subtypes and the clinical features and prognosis of breast cancer
To examine the clinical value of the two subtypes identified by NMF clustering, we compared the relationships between molecular subtypes, clinicopathological characteristics, and prognosis. Clinicopathological features included age, gender, grade, tumor stage, TNM stage, HER2 mutation, ER status, follow-up time, and survival status. Furthermore, the differences in overall survival (OS) among different subtypes were assessed using Kaplan–Meier curves generated by R packages “survival” and “survminer” .
Correlations of cuproptosis-based clustering with TME in breast cancer
The ESTIMATE algorithm was employed to evaluate the immune and stromal scores of each patient . In addition, the proportion of 22 human immune cell subsets of every breast cancer sample were calculated by the CIBERSORT algorithm . Furthermore, immune cell infiltration in the breast cancer TME were also determined using a single-sample gene set enrichment analysis (ssGSEA) algorithm .
Differentiated expressed genes (DEGs) identification and functional annotation
DEGs between CuRGcluster C1 and C2 were identified using R package “limma” with a fold-change of 0.05 and an adjusted p-value of < 0.05. In order to further explore the potential functions of cuproptosis pattern-related DEGs and identify the related gene functions and enriched pathways, functional enrichment analyses were executed on the DEGs using the R package “clusterprofiler”.
Construction of the cuproptosis-related prognostic CuRG_score
To investigate whether cuproptosis was associated with the outcomes of breast cancer patients, we aimed to establish a CuRG_score. First, the top three highly expressed cuproptosis-related genes in geneclusterC were exposed to the univariate Cox regression analysis and two significant genes were subsequently exposed to a LASSO cox regression analysis to compute the exact coefficient values of each identified association. Then, all BRCA patients were randomly categorized into training (n = 686) and testing (n = 685) sets at a ratio of 1:1, then the former was used to construct the cuproptosis-related prognostic CuRG_score. Eventually, two candidate genes were established to get a prognostic CuRG_score in the training set. The CuRG_score was calculated as: CuRG_score = -0.469160009286577*(LIPT1) + 0.322953604620369*(DLAT). Based on the median risk score, a total of 686 patients in the training set were divided into low-risk (CuRG_score < median value) and high-risk (CuRG_score > median value) groups and then subjected to Kaplan–Meier survival analysis. Afterwards, principal component analysis (PCA) was performed using R package “ggplot2”. Similarly, the testing and all sets were divided into low- and high-risk groups, each of which was subjected to Kaplan–Meier survival analysis and the generation of receiver operating characteristic (ROC) curves.
Clinical correlation and stratification analyses of the CuRG_score. Chi-square tests were used to explore the relationships between the CuRG_score and the clinical characteristics (age, gender, TNM stage, HER2 status, and ER status). To assess whether risk scores were independent of other available clinicopathological features, we subjected the training and testing sets to univariate and multivariate analyses. In addition, we performed a stratified analysis to determine whether the CuRG_score retained its predictive ability in different subgroups according to age, gender, T stage, N stage, M stage, tumor stage, HER2 status, and ER status.
Evaluation of immune status, microsatellite instability (MSI), and stemness index between the high- and low-risk groups
To evaluate the proportions of tumor infiltrating immune cells in the TME, CIBERSORT was employed to quantify the abundance of 22 infiltrating immune cells in heterogeneous samples in the low- and high-risk groups. We explored the associations between the fractions of 22 types of infiltrating immune cells and eight genes in the CuRG_score. Boxplots were used to examine the differential expression levels of immune checkpoints between the low- and high-score groups. Furthermore, we analyzed the relationships between the two risk groups and MSI and CSC.
Mutation and drug susceptibility analysis
To determine the somatic mutations of breast cancer patients between high- and low-risk groups, the mutation annotation format from the TCGA database was generated using the “maftools” R package. We also calculated the tumor mutation burden (TMB) score for each patient with BRCA in the two groups. To explore differences in the therapeutic effects of chemotherapeutic drugs in patients in the two groups, we calculated the semi-inhibitory concentration (IC50) values of chemotherapeutic drugs commonly used to treat BRCA using the “pRRophetic” package.
Establishment and validation of a nomogram scoring system
The clinical characteristics and risk score were used to develop a predictive nomogram using the “rms” package based on the outcome of the independent prognosis analysis. In the nomogram scoring system, each variable was matched with a score, and the total score was obtained by adding the scores across all variables of each sample. Time-dependent ROC curves for 1-, 3-, and 5-year survivals were used to assess the nomogram. Calibration plots of the nomogram were used to depict the predictive value between the predicted 1-, 3-, and 5-year survival events and the virtually observed outcomes.
All statistical analyses were performed using R version 4.1.0. Statistical significance was set at p < 0.05.
c-Myc was a potential regulator of breast cancer intra-tumoral heterogeneity
We first investigated the expression of the oncoprotein c-Myc and four breast cancer stem cell (BCSC) biomarkers: CD44, CD24, CD133 and ALDH1A1 . We performed immunohistochemical (IHC) staining on surgical specimens of the primary cancer tissue (PC) and an axillary lymph node with cancer evasion (LN +). Paired surgical specimens were derived from 60 patients diagnosed with various breast cancer subtypes (Table S1). For 60 breast cancer cases, paraffin-embedded tissue from primary cancers and lymph nodes was available for a total of 120 samples. Protein expression of each sample was assessed based on mean expression of three slices. Then, we evaluated the expression of c-Myc, CD44, CD24, CD133, and ALDH1A1 of sixty primary cancers and paired lymph nodes with cancer evasion (Fig. 1A). In both PC and LN + , c-Myc expression positively correlated with CD44, CD133, and ALDH1A1, and negatively correlated with CD24 (Fig. 1B and Table 2). This finding verified the positive correlation between c-Myc and BCSC biomarkers [41,42,43,44]. Since these tissues were collected from paired PC and LN + of the same patients, we also performed Spearman’s correlation analysis on stemness biomarkers in respective tumor sites and identified that a more prominent positive correlation between c-Myc and BCSC phenotype in PC compared with LN + (Table 2), suggesting a potentially more prominent regulation of c-Myc to stemness in the primary tissue. This finding was consistent with our previous report that PC-derived BCSC differentiated into cancer cell subclusters with better metastatic potential . Therefore, based on IHC staining of 120 breast cancer surgical specimens, we identified the positive correlation between c-Myc and BCSC phenotypes, and supposed that the oncoprotein c-Myc was a driving regulator of cancer stemness in breast cancer.
Cuproptosis was dysregulated in breast cancer tissues
Copper-induced cell death (cuproptosis) is a newly reported form of cell death. Excessive copper binds to constituents of the TCA cycle in mitochondria, leading to lipoylated protein aggregation, proteotoxic stress and ultimately cell death . Seven genes (FDX1, LIAS, LIPT1, DLD, DLAT, PDHA1, and PDHB) were verified to rescue cells from copper toxicity, and three genes (MTF1, GLS, and CDKN2A) exacerbated cell death induced by copper overloading. Since dysregulation of intracellular homeostasis and cell cycle contribute to multiple pathologies including cancer , we interrogated the expression of cuproptosis biomarkers in breast cancer specimens (Fig. 2A). We conducted transcriptome analysis on a large-scale sample derived from four open databases (TCGA, GSE42568, GSE48390, and GSE88770). Nine out of ten cuproptosis biomarkers were differentially expressed between breast cancer tissues and normal control. PDHB and CDKN2A expression were significantly higher in breast cancer tissues compared with control. PDHB encodes pyruvate dehydrogenase, which links glycolysis to the TCA cycle and can be mediated by the oncogenic AMPK signaling pathway . Moreover, murine models demonstrated that ectopic activation of the pyruvate dehydrogenase complex by exogenous expression of PDHB increased metastatic potential and survival of cancer cells . CDKN2A is a tumor suppressor gene and its target, p16, inactivates the cyclinD-CDK4/6 complex and leads to cell cycle arrest . CDKN2A alteration is a contributor of tumorigenesis, however, its overexpression was also reported in several types of cancers due to its association with aberrant apoptosis, senescence, angiogenesis and cancer migration [49, 50]. Inverse correlation between Rb and p16 was reported in several cancers, including breast cancer. Additionally, p16 overexpression and heterozygous Rb loss are predictors to CDK4/6 resistance in hormone receptor-positive breast cancer . Survival analysis demonstrated that among the cuproptosis biomarkers, FDX1, PDHA1, DLAT, and DLD predicts poor prognosis (Fig. 2B–E).
Identifying predictive biomarkers of breast cancer based on cuproptosis
To further identify the potential correlation between cuproptosis and breast cancer prognosis, we conducted unsupervised clustering on all breast cancer cases and grouped all patients into two CuRGClusters (C1 and C2, Fig. 2F). Bioinformatic analysis of the survival curve for the two CuRGClusters identified that patients from C2 had worse prognosis than C1 (Fig. 2G). We then investigated the expression of cuproptosis-related genes in the two CuRGClusters. Expression of FDX1, DLD, DLAT, PDHA1, GLS and CDKN2A were more enriched in the cluster with worse prognosis (Fig. 3A). C1 showed significantly higher expression of LIAS, LIPT1, and PDHB. C2 showed significantly higher expression of FDX1, DLD, DLAT, PDHA1, GLS, and CDKN2A. For C2, we identified increased enrichment of pathways including ROS pathway, glycolysis and mTORC1 signaling. Because cuproptosis-mediated cellular damage mainly targets mitochondrial TCA, and c-Myc was reported to maintain breast cancer stemness by increasing oxidative phosphorylation (OXPHOS) and producing reactive oxygen species (ROS) , we interrogated the enrichment of c-Myc-associated pathways in these breast cancer specimens. Notably, we found that tissue samples of C2 cluster showed putative enrichment of c-Myc-targeted pathways (Fig. 3B). This finding prompted us to speculate that cuproptosis was the potential downstream of c-Myc-mediated tumor progression. Additional file 1.
To investigate whether cuproptosis exerts effects on the TME, we evaluated the immune cell infiltration of breast cancer patients. In the breast cancer subgroup with worse prognosis, inhibitory immune response was identified by a higher proportion of inhibitory T cells and immunosuppressive components represented by myeloid-derived suppressor cells (MDSC) and regulatory T (Treg) cells (Fig. 3C). A larger proportion of γδT also suggested immune tolerance in the breast cancer samples from C2. γδT is a CD4-negative and CD8-negative T-cell subgroup having both anti- and pro-tumoral functions during cancer progression . On the one hand, γδT present antigens and exert direct cytotoxicity, on the other, certain γδT subgroups express exhausted ligands, facilitate pro-tumoral chronic inflammation and attenuate adaptive immunity . Based on transcriptomics and the enriched functional pathways, both CuRGClusters (C1 and C2) are active in mitochondrial biological functions such as acetyl-CoA metabolism, purine nucleoside bisphosphate biosynthesis, and reducing equivalents production (Fig. 3D).
Cuproptosis was predictive to clinical outcomes of breast cancer
To further identify the pro-tumoral effects of cuproptosis-related mechanisms and compute the exact coefficient values of each identified association, we performed univariate Cox regression analysis on the differentially expressed cuproptosis biomarkers of C1 and C2. The univariate Cox regression demonstrated that LIPT1, PDHA1, and DLAT as predictor to poor survival and further classified breast cancer samples into three geneClusters (geneClusterA, geneClusterB and geneClusterC, Fig. 4A, B). Then, we performed LASSO regression analysis on the three genes to establish a CuRG_score that evaluate the degree of cuproptosis in breast cancer samples (Fig. 4C). Eventually, we got the CuRG_score calculated as: CuRG_score = − 0.469160009286577* (LIPT1) + 0.322953604620369* (DLAT). Breast cancer samples were divided into the high- and low-risk group based on CuRG_score (Fig. 4D). We checked the correspondence between CuRGclusters classified by differentially expressed cuproptosis biomarkers (C1 and C2), geneClusters classified by univariate Cox regression (A, B, C), and the risk groups classified by CuRG_scores (Fig. 4E). It turned out that previous clusters with comparatively poor outcomes (C2 compared with C1, geneClusterC compared with geneClusterA and geneClusterB) showed putatively higher scores of cuproptosis. Furthermore, survival curve further verified that higher CuRG_score is associated with poorer clinical prognosis (Fig. 4F). The high-risk group had higher expression of DLAT, DLD, FDX1, GLS, MTF1, and PDHA1, and lower expression of LIAS and LIPT1 compared with the low-risk group (Fig. 4G). Based on the above analysis of large-scale transcriptomics, we identified that cuproptosis was predictive to poorer clinical outcomes of patients with breast cancer.
Cuproptosis engaged in c-Myc-mediated cancer stemness
c-Myc is an oncogene regulating cancer cell proliferation through multiple mechanisms and modulating tissue-specific chemotherapy sensitivity via mitochondrial apoptosis . Moreover, c-Myc contributes to maintaining BCSC by increasing OXPHOS and producing ROS . The above bioinformatic analysis identified that cuproptosis predicted poor clinical outcomes of patients with breast cancer, and suggested that this pro-tumoral effect was significantly associated with c-Myc-mediated downstream targets, we inferred that cuproptosis engaged in c-Myc-mediated tumor progression. To further investigate whether the relationship between c-Myc and cuproptosis was also associated to cancer stemness, an embodiment of c-Myc-mediated carcinogenic property, we performed another correlation analysis on stemness and cuproptosis. We found a significant positive correlation between CuRG_score and breast cancer stemness phenotype (Fig. 5A). These findings collectively indicated that cuproptosis  was engaged in the c-Myc-mediated breast cancer stemness and malignancy.
It is acknowledged that cancer stem cell is a regulator of intra-tumoral heterogeneity and a marker of chemo- and radio-resistance. To test whether this CuRG_score predicted therapeutic effects, we investigated the sensitivity to several anti-tumor drugs that were frequently used. The high-risk group had better responses to AKT inhibitors and vinorelbine than the low-risk group (Fig. 5B). AKT inhibitors in combination with first-line chemotherapy  and vinorelbine-based combinatorial therapy [57, 58] both showed therapeutic effects in breast cancer cases, and we demonstrated that high CuRG_score was significantly associated with higher IC50 of these drugs. TMB has promising predictive ability to identify patients who derive clinically significant improvement in survival from immunotherapy [59, 60]. Bioinformatic analysis demonstrated that high CuRG_score was correlated with lower TMB (Fig. 5C). These findings suggested that cuproptosis engaged in c-Myc-mediated breast cancer stemness, therefore rendering breast cancer cells in a dormient and less mutative state that is insensitive to chemotherapy and immunotherapy.
Since the first report of copper-induced cell death in March 2022, there has been a great deal of interest in this novel concept. So far, various bioinformatic analysis have linked this new form of cell death to multiple diseases, including cancers and non-malignant diseases such as rheumatoid arthritis [61, 62]. Moreover, revelation of the exact mechanisms also boosted copper-based nanomaterials and enriched optional carrier for cancer treatment . Though various studies have established prognostic models or simulated the tumor immune microenvironment based on public database, few focused on the molecular mechanisms and the therapeutic potential of this newly identified cell death. The Science report verified ten genes related to cuproptosis among which seven genes (LIAS, LIPT1, DLD, DLAT, PDHA1, PDHB, FDX1) rescued cells from copper overloading and mitochondrial failure, while the other three (MTF1, GLS, CDKN2A) offers synergy and promotes cuproptosis . Another form of metal-dependent cell death that arouses enormous interest is ferroptosis. After ferroptosis was reported, researchers expanded their study into the molecular pathways and utilized cell death resulted from ferroptosis-mediated lipid peroxide accumulation to target cancerous proliferation. The key enzymes and regulatory mechanisms involved in cuproptosis have been identified. However, its role in pathological conditions as well as the potential in clinical practice remains unclear. Inspired by ferroptosis-mediated cancer biology, researchers devoted to the potential value of cuproptosis in cancer. By combining copper toxicity with RNA methylation, a report on hepatocellular carcinoma demonstrates the RNA regulators of cuproptosis biomarkers and shows correlation with tumor mutation which partly predicts sensitivity to immunotherapy . Moreover, studies in combination with lncRNA network indicate that cuproptosis-related lncRNAs can predict the recruitment and infiltration of immune components in TME . These interesting findings motivate us to take cuproptosis as an opportunity to extend our research in oncology.
Heterogeneity is one of the challenges of cancer treatment. Intratumoral heterogeneity, demonstrated in spatiotemporal evolution and instability, greatly contributed to failure of monotherapy, therefore leading to tumor progression and metastasis. CSCs are a unique cancer cell population comparatively immune-privileged than other subclusters that drive the initiation, metastasis and therapeutic resistance of most neoplasms [65, 66]. Investigation into heterogeneity, as well as the contributing factors underlying this intratumoral heterogeneity, is highly valuable to research and clinical practice. c-Myc has been recognized as a “grand orchestrator” of carcinogenesis that mediates tumor progression. In addition to being an important transcriptional factor of cancer cells, c-Myc dysregulates the TME and remodels antitumor immunity .
Our previous bioinformatic analysis suggested that c-Myc was a driving factor of breast cancer stemness. Therefore, we aimed to further specify the specific mechanisms by which c-Myc mediates breast cancer stemness and intratumoral heterogeneity. Start with IHC staining of breast cancer surgical specimens, we verified that c-Myc is putatively related to BCSC phenotype. Our analysis based on large-scale transcriptomics profiles supported that cuproptosis-related poor prognosis is a downstream target of c-Myc-mediated mechanisms. Further statistical analysis and screening established a CuRG_score as an indicator to poor prognosis and showed putative correlation with breast cancer stemness. The above findings proposed a novel mechanism of cuproptosis in breast cancer, that is cuproptosis engaged in breast cancer stemness via c-Myc-related mechanisms. We proposed that cuproptosis was a potential target of the classic oncogene c-Myc, the intermediator of c-Myc-mediated breast cancer stemness and a predictor to poor prognosis of breast cancer.
In conclusion, based on a sample size of 1,533 breast cancer cases, we combined large-scale transcriptomics analysis and pathological examination and revealed that cuproptosis engaged in c-Myc-mediated breast cancer stemness. This finding suggests the potential of innovative nanomaterials that pack c-Myc-targeted inhibitors with copper polymers.
Data used for bioinformatic analysis in this manuscript was acquired from the GEO database (index number: GSE42568, GSE48390, and GSE88770).
Siegel RL, Miller KD, Fuchs HE, Jemal A. Cancer statistics, 2021. CA Cancer J Clin. 2021;71:7–33.
Siegel RL, Miller KD, Fuchs HE, Jemal A. Cancer statistics, 2022. CA Cancer J Clin. 2022;72:7–33.
Batlle E, Clevers H. Cancer stem cells revisited. Nat Med. 2017;23:1124–34.
Xu K, Wang R, Xie H, Hu L, Wang C, Xu J, Zhu C, Liu Y, Gao F, Li X, et al. Single-cell RNA sequencing reveals cell heterogeneity and transcriptome profile of breast cancer lymph node metastasis. Oncogenesis. 2021;10:66.
Beck B, Blanpain C. Unravelling cancer stem cell potential. Nat Rev Cancer. 2013;13:727–38.
Najafi M, Mortezaee K, Ahadi R. Cancer stem cell (a) symmetry and plasticity: tumorigenesis and therapy relevance. Life Sci. 2019;231: 116520.
Zhang H, Brown RL, Wei Y, Zhao P, Liu S, Liu X, Deng Y, Hu X, Zhang J, Gao XD, et al. CD44 splice isoform switching determines breast cancer stem cell state. Genes Dev. 2019;33:166–79.
Zhang R, Tu J, Liu S. Novel molecular regulators of breast cancer stem cell plasticity and heterogeneity. Semin Cancer Biol. 2022;82:11–25.
Wei W, Hu H, Tan H, Chow LW, Yip AY, Loo WT. Relationship of CD44+CD24-/low breast cancer stem cells and axillary lymph node metastasis. J Transl Med. 2012;10(Suppl 1):S6.
Jafari SM, Joshaghani HR, Panjehpour M, Aghaei M, Zargar Balajam N. Apoptosis and cell cycle regulatory effects of adenosine by modulation of GLI-1 and ERK1/2 pathways in CD44 (+) and CD24 (−) breast cancer stem cells. Cell Prolif. 2017. https://doi.org/10.1111/cpr.12345.
Ryoo IG, Choi BH, Ku SK, Kwak MK. High CD44 expression mediates p62-associated NFE2L2/NRF2 activation in breast cancer stem cell-like cells: Implications for cancer stem cell resistance. Redox Biol. 2018;17:246–58.
Yang F, Cao L, Sun Z, Jin J, Fang H, Zhang W, Guan X. Evaluation of breast cancer stem cells and intratumor stemness heterogeneity in triple-negative breast cancer as prognostic factors. Int J Biol Sci. 2016;12:1568–77.
Hsu TY, Simon LM, Neill NJ, Marcotte R, Sayad A, Bland CS, Echeverria GV, Sun T, Kurley SJ, Tyagi S, et al. The spliceosome is a therapeutic vulnerability in MYC-driven cancer. Nature. 2015;525:384–8.
Xu Z, Xu M, Liu P, Zhang S, Shang R, Qiao Y, Che L, Ribback S, Cigliano A, Evert K, et al. The mTORC2-Akt1 cascade is crucial for c-Myc to promote hepatocarcinogenesis in mice and humans. Hepatology. 2019;70:1600–13.
Huang X, Sun Q, Fu H, Zhou X, Guan X, Wang J. Both c-Myc and Ki-67 expression are predictive markers in patients with extranodal NK/T-cell lymphoma, nasal type: a retrospective study in China. Pathol Res Pract. 2014;210:351–6.
Cho KB, Cho MK, Lee WY, Kang KW. Overexpression of c-myc induces epithelial mesenchymal transition in mammary epithelial cells. Cancer Lett. 2010;293:230–9.
Lawson DA, Bhakta NR, Kessenbrock K, Prummel KD, Yu Y, Takai K, Zhou A, Eyob H, Balakrishnan S, Wang CY, et al. Single-cell analysis reveals a stem-cell program in human metastatic breast cancer cells. Nature. 2015;526:131–5.
Perez EA, Jenkins RB, Dueck AC, Wiktor AE, Bedroske PP, Anderson SK, Ketterling RP, Sukov WR, Kanehira K, Chen B, et al. C-MYC alterations and association with patient outcome in early-stage HER2-positive breast cancer from the north central cancer treatment group N9831 adjuvant trastuzumab trial. J Clin Oncol. 2011;29:651–9.
Liu Y, Zhu C, Tang L, Chen Q, Guan N, Xu K, Guan X. MYC dysfunction modulates stemness and tumorigenesis in breast cancer. Int J Biol Sci. 2021;17:178–87.
Gao FY, Li XT, Xu K, Wang RT, Guan XX. c-MYC mediates the crosstalk between breast cancer cells and tumor microenvironment. Cell Commun Signal. 2023;21:28.
Tsvetkov P, Coy S, Petrova B, Dreishpoon M, Verma A, Abdusamad M, Rossen J, Joesch-Cohen L, Humeidi R, Spangler RD, et al. Copper induces cell death by targeting lipoylated TCA cycle proteins. Science. 2022;375:1254–61.
Hu Q, Wang R, Ma H, Zhang Z, Xue Q. Cuproptosis predicts the risk and clinical outcomes of lung adenocarcinoma. Front Oncol. 2022;12: 922332.
Song S, Zhang M, Xie P, Wang S, Wang Y. Comprehensive analysis of cuproptosis-related genes and tumor microenvironment infiltration characterization in breast cancer. Front Immunol. 2022;13: 978909.
Huang T, Liu Y, Li J, Shi B, Shan Z, Shi Z, Yang Z. Insights into prognosis and immune infiltration of cuproptosis-related genes in breast cancer. Front Immunol. 2022;13:1054305.
Tan PH, Ellis I, Allison K, Brogi E, Fox SB, Lakhani S, Lazar AJ, Morris EA, Sahin A, Salgado R, et al. The 2019 World Health Organization classification of tumours of the breast. Histopathology. 2020;77:181–5.
Xie P, Zhang M, He S, Lu K, Chen Y, Xing G, Lu Y, Liu P, Li Y, Wang S, et al. The covalent modifier Nedd8 is critical for the activation of Smurf1 ubiquitin ligase in tumorigenesis. Nat Commun. 2014;5:3733.
Clarke C, Madden SF, Doolan P, Aherne ST, Joyce H, O’Driscoll L, Gallagher WM, Hennessy BT, Moriarty M, Crown J, et al. Correlating transcriptional networks to breast cancer survival: a large-scale coexpression analysis. Carcinogenesis. 2013;34:2300–8.
Huang CC, Tu SH, Lien HH, Jeng JY, Huang CS, Huang CJ, Lai LC, Chuang EY. Concurrent gene signatures for han chinese breast cancers. PLoS ONE. 2013;8: e76421.
Chen, YJ., Huang, CS., Phan, NN., Lu, TP., Liu, CY., Huang, CJ., Chiu, JH., Tseng, LM. & Huang, CC. Molecular subtyping of breast cancer intrinsic taxonomy with oligonucleotide microarray and NanoString nCounter, Biosci Rep, https://doi.org/10.1042/BSR20211428 (2021).
Metzger-Filho O, Michiels S, Bertucci F, Catteau A, Salgado R, Galant C, Fumagalli D, Singhal SK, Desmedt C, Ignatiadis M, et al. Genomic grade adds prognostic value in invasive lobular carcinoma. Ann Oncol. 2013;24:377–84.
Hao Y, Hao S, Andersen-Nissen E, Mauck WM 3rd, Zheng S, Butler A, Lee MJ, Wilk AJ, Darby C, Zager M, et al. Integrated analysis of multimodal single-cell data. Cell. 2021;184:3573-3587.e3529.
Stuart T, Butler A, Hoffman P, Hafemeister C, Papalexi E, Mauck WM 3rd, Hao Y, Stoeckius M, Smibert P, Satija R. Comprehensive integration of single-cell data. Cell. 2019;177:1888-1902.e1821.
Butler A, Hoffman P, Smibert P, Papalexi E, Satija R. Integrating single-cell transcriptomic data across different conditions, technologies, and species. Nat Biotechnol. 2018;36:411–20.
Satija R, Farrell JA, Gennert D, Schier AF, Regev A. Spatial reconstruction of single-cell gene expression data. Nat Biotechnol. 2015;33:495–502.
Gaujoux R, Seoighe C. A flexible R package for nonnegative matrix factorization. BMC Bioinformatics. 2010;11:367.
Wickham H: ggplot2: Elegant Graphics for Data Analysis. 2016. https://doi.org/10.1007/978-3-319-24277-4.
Yoshihara K, Shahmoradgoli M, Martínez E, Vegesna R, Kim H, Torres-Garcia W, Treviño V, Shen H, Laird PW, Levine DA, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. 2013;4:2612.
Newman AM, Steen CB, Liu CL, Gentles AJ, Chaudhuri AA, Scherer F, Khodadoust MS, Esfahani MS, Luca BA, Steiner D, et al. Determining cell type abundance and expression from bulk tissues with digital cytometry. Nat Biotechnol. 2019;37:773–82.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102:15545–50.
Vora P, Venugopal C, Salim SK, Tatari N, Bakhshinyan D, Singh M, Seyfrid M, Upreti D, Rentas S, Wong N, et al. The rational development of CD133-targeting immunotherapies for glioblastoma. Cell Stem Cell. 2020;26:832-844.e836.
Liu M, Liu Y, Deng L, Wang D, He X, Zhou L, Wicha MS, Bai F, Liu S. Transcriptional profiles of different states of cancer stem cells in triple-negative breast cancer. Mol Cancer. 2018;17:65.
Ginestier C, Hur MH, Charafe-Jauffret E, Monville F, Dutcher J, Brown M, Jacquemier J, Viens P, Kleer CG, Liu S, et al. ALDH1 is a marker of normal and malignant human mammary stem cells and a predictor of poor clinical outcome. Cell Stem Cell. 2007;1:555–67.
Wang C, Xu K, Wang R, Han X, Tang J, Guan X. Heterogeneity of BCSCs contributes to the metastatic organotropism of breast cancer. J Exp Clin Cancer Res. 2021;40:370.
Yang F, Xu J, Tang L, Guan X. Breast cancer stem cell: the roles and therapeutic implications. Cell Mol Life Sci. 2017;74:951–66.
Thuy LT, Lee S, Dongquoc V, Choi JS. Nanoemulsion composed of α-tocopherol succinate and dequalinium shows mitochondria-targeting and anticancer effects. Antioxidants. 2023. https://doi.org/10.3390/antiox12020437.
Li N, Li H, Wang Y, Cao L, Zhan X. Quantitative proteomics revealed energy metabolism pathway alterations in human epithelial ovarian carcinoma and their regulation by the antiparasite drug ivermectin: data interpretation in the context of 3P medicine. Epma J. 2020;11:661–94.
Cai Z, Li CF, Han F, Liu C, Zhang A, Hsu CC, Peng D, Zhang X, Jin G, Rezaeian AH, et al. Phosphorylation of PDHA by AMPK drives TCA cycle to promote cancer metastasis. Mol Cell. 2020;80:263-278.e267.
Liggett WH Jr, Sidransky D. Role of the p16 tumor suppressor gene in cancer. J Clin Oncol. 1998;16:1197–206.
Lee CT, Capodieci P, Osman I, Fazzari M, Ferrara J, Scher HI, Cordon-Cardo C. Overexpression of the cyclin-dependent kinase inhibitor p16 is associated with tumor recurrence in human prostate cancer. Clin Cancer Res. 1999;5:977–83.
Romagosa C, Simonetti S, López-Vicente L, Mazo A, Lleonart ME, Castellvi J. Ramon y Cajal S: p16(Ink4a) overexpression in cancer: a tumor suppressor gene associated with senescence and high-grade tumors. Oncogene. 2011;30:2087–97.
Palafox M, Monserrat L, Bellet M, Villacampa G, Gonzalez-Perez A, Oliveira M, Brasó-Maristany F, Ibrahimi N, Kannan S, Mina L, et al. High p16 expression and heterozygous RB1 loss are biomarkers for CDK4/6 inhibitor resistance in ER (+) breast cancer. Nat Commun. 2022;13:5258.
Lee KM, Giltnane JM, Balko JM, Schwarz LJ, Guerrero-Zotano AL, Hutchinson KE, Nixon MJ, Estrada MV, Sánchez V, Sanders ME, et al. MYC and MCL1 cooperatively promote chemotherapy-resistant breast cancer stem cells via regulation of mitochondrial oxidative phosphorylation. Cell Metab. 2017;26:633-647.e637.
Wu D, Wu P, Qiu F, Wei Q, Huang J. Human γδT-cell subsets and their involvement in tumor immunity. Cell Mol Immunol. 2017;14:245–53.
Daley D, Zambirinis CP, Seifert L, Akkad N, Mohan N, Werba G, Barilla R, Torres-Hernandez A, Hundeyin M, Mani VRK, et al. γδ T cells support pancreatic oncogenesis by restraining αβ T cell activation. Cell. 2016;166:1485-1499.e1415.
Sarosiek KA, Fraser C, Muthalagu N, Bhola PD, Chang W, McBrayer SK, Cantlon A, Fisch S, Golomb-Mello G, Ryan JA, et al. Developmental regulation of mitochondrial apoptosis by c-Myc governs age- and tissue-specific sensitivity to cancer therapeutics. Cancer Cell. 2017;31:142–56.
Pascual J, Turner NC. Targeting the PI3-kinase pathway in triple-negative breast cancer. Ann Oncol. 2019;30:1051–60.
Harbeck N, Huang CS, Hurvitz S, Yeh DC, Shao Z, Im SA, Jung KH, Shen K, Ro J, Jassem J, et al. Afatinib plus vinorelbine versus trastuzumab plus vinorelbine in patients with HER2-overexpressing metastatic breast cancer who had progressed on one previous trastuzumab treatment (LUX-Breast 1): an open-label, randomised, phase 3 trial. Lancet Oncol. 2016;17:357–66.
Martín M, Ruiz A, Muñoz M, Balil A, García-Mata J, Calvo L, Carrasco E, Mahillo E, Casado A, García-Saenz JA, et al. Gemcitabine plus vinorelbine versus vinorelbine monotherapy in patients with metastatic breast cancer previously treated with anthracyclines and taxanes: final results of the phase III Spanish breast cancer research group (GEICAM) trial. Lancet Oncol. 2007;8:219–25.
Chowell D, Yoo SK, Valero C, Pastore A, Krishna C, Lee M, Hoen D, Shi H, Kelly DW, Patel N, et al. Improved prediction of immune checkpoint blockade efficacy across multiple cancer types. Nat Biotechnol. 2022;40:499–506.
Gandara DR, Paul SM, Kowanetz M, Schleifman E, Zou W, Li Y, Rittmeyer A, Fehrenbacher L, Otto G, Malboeuf C, et al. Blood-based tumor mutational burden as a predictor of clinical benefit in non-small-cell lung cancer patients treated with atezolizumab. Nat Med. 2018;24:1441–8.
Zhao J, Guo S, Schrodi SJ, He D. Cuproptosis and cuproptosis-related genes in rheumatoid arthritis: Implication, prospects, and perspectives. Front Immunol. 2022;13: 930278.
Xu Y, Liu SY, Zeng L, Ma H, Zhang Y, Yang H, Liu Y, Fang S, Zhao J, Xu Y, et al. An enzyme-engineered nonporous copper(I) coordination polymer nanoplatform for cuproptosis-based synergistic cancer therapy. Adv Mater. 2022;34: e2204733.
Li D, Shi Z, Liu X, Jin S, Chen P, Zhang Y, Chen G, Fan X, Yang J, Lin H. Identification and development of a novel risk model based on cuproptosis-associated RNA methylation regulators for predicting prognosis and characterizing immune status in hepatocellular carcinoma. Hepatol Int. 2023;17:112–30.
Li D, Jin S, Chen P, Zhang Y, Li Y, Zhong C, Fan X, Lin H. Comprehensive analysis of cuproptosis-related lncRNAs for prognostic significance and immune microenvironment characterization in hepatocellular carcinoma. Front Immunol. 2022;13:991604.
Black JRM, McGranahan N. Genetic and non-genetic clonal diversity in cancer evolution. Nat Rev Cancer. 2021;21:379–92.
Vitale I, Shema E, Loi S, Galluzzi L. Intratumoral heterogeneity in cancer progression and response to immunotherapy. Nat Med. 2021;27:212–24.
Dhanasekaran R, Deutzmann A, Mahauad-Fernandez WD, Hansen AS, Gouw AM, Felsher DW. The MYC oncogene—the grand orchestrator of cancer growth and immune evasion. Nat Rev Clin Oncol. 2022;19:23–36.
This work was supported by a Key International Cooperation of the National Natural Science Foundation of China (81920108029) and the Key Foundation for Social Development Project of the Jiangsu Province of China (BE2021741). Authors would like to thank the experiment center of Nanjing Medical University for instrumental assistance.
Competing of interests
The authors declare no competing of interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Wang, R., Xu, K., Chen, Q. et al. Cuproptosis engages in c-Myc-mediated breast cancer stemness. J Transl Med 21, 409 (2023). https://doi.org/10.1186/s12967-023-04204-5
- Intra-tumoral heterogeneity
- Cancer stem cell
- Copper-induced cell death
- Tumor microenvironment