- Open Access
Pan-cancer patterns of cuproptosis markers reveal biologically and clinically relevant cancer subtypes
Biomarker Research volume 11, Article number: 13 (2023)
Cuproptosis is a newly discovered type of cell death triggered by copper accumulation. Here we exhibited the genetic profiles of 10 cuproptosis-associated genes (CuAGs) across 21 cancer types. Only 8.0% (627/7839) of tumors possessed at least 1 mutation on CuAGs, while the copy number amplifications or deletions on the alleles of CuAGs were ubiquitous. Generally, the expression of CuAGs showed heterogeneity across cancer types and the expression of CuAGs showed different correlations with MSI, TMB, immune and stromal features in different cancer types. Therefore, CuAGs were ubiquitously and heterogeneously dysregulated in pan-cancer. With a Non-negative Matrix Factorization method, we divided patients of each cancer type into cuproptosis-based subtypes, which showed a close but heterogeneous correlation with different biological and clinical features. Accordingly, we summarized all cancer types into four categories. The cancers in which cuproptosis subtypes correlated with MSI and TMB were annotated as Genomic disturbed. Those correlated with stromal scores were categorized as Stromal remolded. The others only associated with immune infiltration were labeled as Immune inhibited. A minor fraction of cancers not correlated with any biological indicators were marked as Cuproptosis inert. Together, we provided a pan-cancer overview of cuproptosis markers which revealed biologically and clinically relevant cancer subtypes in different cancers.
To the editor
Cuproptosis, a unique type of cell death triggered by copper accumulation, was newly characterized . The inner mechanism of cuproptosis involved mitochondrial protein lipoylation, tricarboxylic acid cycle, and reactive oxygen species generation, which were strongly associated with tumor microenvironment (TME) homeostasis [1,2,3]. However, whether the cuproptosis-associated genes (CuAGs) were associated with cancer progression and the TME features were still unknown. In this study, we collected CuAGs which encoded the components of the lipoic acid pathway (FDX1, LIAS, LIPT1, DLD) and the pyruvate dehydrogenase complex (DLAT, PDHA1, PDHB, MTF1, GLS, CDKN2A), including 3 negative regulators (i.e., MTF1, GLS, CDKN2A) and other 7 positive regulators. Then we employed gene expression data of those 10 CuAGs from the TCGA database to reveal the cuproptosis-associated distinct regulatory patterns in different cancers.
Here we found that the overall somatic mutational frequency of the 10 CuAGs was at a moderate level (8.0%, 627/7839), and CDKN2A exhibited the highest somatic mutation rates (Fig. 1A, Additional file 1: Figure S1A). Compared with the mutation status, the CNV alterations are much more prevalent, and the CNV status of each CuAG was comparably consistent among different cancer types (Fig. 1B). The expression of CuAGs exhibited high heterogeneity among different cancers (Fig. 1C), among which CDKN2A was overexpressed in most of the cancers, even with a deletion in CNV (Fig. 1D, Additional file 1: Figure S1B-K). Generally, all of the CuAGs showed expressional dysregulation and prognostic relevance (Fig. 1D-E), among which DLAT was the most widely applicable prognostic CuAG that correlated with the prognosis of 7 different cancers (Additional file 2: Table S1). KIRC was the cancer type most affected by CuAGs in prognosis, in which 7/10 CuAGs were prognostic (Additional file 2: Table S 1). Moreover, CDKN2A, DLAT, GLS, LIAS, and PDHB were also identified as potential chemotherapy response predictors (Additional file 1: Figure S2).
Then the correlations between the CuAGs and biological features of different cancers were evaluated. Most CuAGs exhibited close and consistent associations with the immune score and the stromal score (Fig. 1F-J, Additional file 2: Table S2 and S3). The 7 positive cuproptosis regulators exhibited negative correlations with both immune and stromal scores in the majority of the cancer types, while the 3 negative cuproptosis regulators shared positive associations with TME scores (Additional file 1: Figure S3 and S4). Comparably, the regulation of CuAGs on TMB and MSI was heterogeneous (Additional file 1: Figure S5 and S6, Additional file 2: Table S4 and S5), suggesting that cuproptosis might affect genomic stability only in some cancer types.
To evaluate the overall correlation between cuproptosis and different biological features in various cancers, we employed the non-negative matrix (NMF) factorization to classify patients based on all CuAGs . For each cancer, patients were assorted into 2-4 subgroups (Fig. 1K, Additional file 1: Figure S7). The 2D t-SNE plot verified the mathematical validity of the unsupervised clustering results by NMF (Fig. 1L, Additional file 1: Figure S8). In most cancer types, the NMF subgroups exhibited a strong correlation with immune (15/21) and stromal scores (11/21)  (Fig. 1M-N, Additional file 1: Figure S9A-B). Correspondently, correlations between genomic stability indicators [6, 7] and NMF subgroups were only identified in a small portion of cancers (9/21 for TMB score, Fig. 1O, Additional file 1: Figure S9C; 6/21 for MSI scores, Fig. 1P, Additional file 1: Figure S9D). Additionally, in most cancer types, subgroups showed close association with cuproptosis-related biological processes, such as hypoxia and reactive oxygen species (ROS) pathways [8, 9] (Additional file 1: Figure S10). Meanwhile, in COAD, HNSC, KIRC, LIHC, LUAD, and UCEC, different NMF subgroups also exhibited distinct overall survivals (Fig. 1Q-V, Additional file 1: Figure S11).
Since CuAGs showed a heterogeneous correlation with different biological and clinical features, we further summarized all cancer types into four categories: Genomic disturbed, Stromal remolded, Immune inhibited, and Cuproptosis inert (Fig. 2A, Additional file 2: Table S6). The Genomic disturbed category refers to the cancers most dramatically affected by CuAGs, and their NMF subgroups were generally associated with both genomic instability indicators and tumor microenvironment (e.g., STAD, Fig. 2B-E). The Stromal remolded category refers to the cancers in which the NMF subgroups correlated with the TME-related features but not the genomic-stability ones. Moreover, the Stromal remolded category was also characterized by closely associated with prognosis (e.g., COAD, Fig. 2F-K, Additional file 1: Figure S12). The Immune inhibited category refers to the cancers in which the NMF subtype exhibited correlations mainly with immune scores. We believe that those cancers were only slightly affected by CuAGs, in which cuproptosis only affected the immune microenvironment (e.g., OV, Fig. 2L-O, Additional file 1: Figure S13). The Cuproptosis inert category, which was annotated for the cancers barely affected by cuproptosis, was not associated with any biological and clinical features (e.g., SARC, Additional file 1: Figure S14).
In conclusion, we demonstrated the widespread genetic dysregulation and distinct biological relevance of cuproptosis regulators in pan-cancer. Moreover, we provided a pan-cancer overview of cuproptosis regulating patterns which revealed biologically and clinically relevant cancer subtypes.
Cancers with abbreviations
ACC, adrenocortical carcinoma; BLCA, bladder urothelial carcinoma; BRCA, breast cancer; CESC, cervical squamous cell carcinoma, and endocervical adenocarcinoma; COAD, colon adenocarcinoma; ESCC, esophageal carcinoma; GBM, glioblastoma multiforme; HNSC, head, and neck squamous carcinoma; KIRC, kidney renal clear cell carcinoma; KIRP, kidney renal papillary cell carcinoma; LAML, acute myeloid leukemia; LIHC, liver hepatocellular carcinoma; LUAD, lung adenocarcinoma; LUSC, lung squamous cell carcinoma; OV, ovarian serous cystadenocarcinoma; PAAD, pancreatic adenocarcinoma; PRAD, prostate adenocarcinoma; READ, rectum adenocarcinoma; SARC, sarcoma; STAD, stomach adenocarcinoma; UCEC, uterine corpus endometrial carcinoma.
Availability of data and materials
All sequencing data were public available datasets (TCGA https://portal.gdc.cancer.gov/cart). All other data supporting the conclusions of this article are presented within the article and its supplementary files.
Copy number variations
Cuproptosis associated gene
Gene set enrichment analysis
Non-negative matrix factorization
Reactive oxygen species
Single-sample gene-set enrichment analysis
The Cancer Genome Atlas
Tumor mutation burden
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. https://doi.org/10.1126/science.abf0529.
Saikolappan S, Kumar B, Shishodia G, Koul S, Koul HK. Reactive oxygen species and cancer: A complex interaction. Cancer Lett. 2019;452:132–43. https://doi.org/10.1016/j.canlet.2019.03.020.
Mukha D, Fokra M, Feldman A, Sarvin B, Sarvin N, Nevo-Dinur K, Besser E, Hallo E, Aizenshtein E, Schug ZT, Shlomi T. Glycine decarboxylase maintains mitochondrial protein lipoylation to support tumor growth. Cell Metab. 2022;34:775–82. https://doi.org/10.1016/j.cmet.2022.04.006.
Yu N, Wu MJ, Liu JX, Zheng CH, Xu Y. Correntropy-Based Hypergraph Regularized NMF for Clustering and Feature Selection on Multi-Cancer Integrated Data. IEEE Trans Cybern. 2021;51:3952–63. https://doi.org/10.1109/TCYB.2020.3000799.
Wu J, Li L, Zhang H, Zhao Y, Zhang H, Wu S, Xu B. A risk model developed based on tumor microenvironment predicts overall survival and associates with tumor immunity of patients with lung adenocarcinoma. Oncogene. 2021;40:4413–24. https://doi.org/10.1038/s41388-021-01853-y.
Jardim DL, Goodman A, de Melo GD, Kurzrock R. The Challenges of Tumor Mutational Burden as an Immunotherapy Biomarker. Cancer Cell. 2021;39:154–73. https://doi.org/10.1016/j.ccell.2020.10.001.
Vilar E, Gruber SB. Microsatellite instability in colorectal cancer-the stable evidence. Nat Rev Clin Oncol. 2010;7:153–62. https://doi.org/10.1038/nrclinonc.2009.237.
Burr SP, Costa AS, Grice GL, Timms RT, Lobb IT, Freisinger P, Dodd RB, Dougan G, Lehner PJ, Frezza C, Nathan JA. Mitochondrial Protein Lipoylation and the 2-Oxoglutarate Dehydrogenase Complex Controls HIF1alpha Stability in Aerobic Conditions. Cell Metab. 2016;24:740–52. https://doi.org/10.1016/j.cmet.2016.09.015.
Feng D, Witkowski A, Smith S. Down-regulation of mitochondrial acyl carrier protein in mammalian cells compromises protein lipoylation and respiratory complex I and results in cell death. J Biol Chem. 2009;284:11436–45. https://doi.org/10.1074/jbc.M806991200.
We acknowledge the owners of TCGA and GSCA for providing the valuable platforms for making all those meaningful data available. All the contributors to those public datasets are deeply appreciated.
This work was supported by grants from the Beijing Nova Program of Science and Technology (Z191100001119128); National Natural Science Foundation of China (82073390); Beijing Municipal Science and Technology Project (Z191100006619081); The Digestive Medical Coordinated Development Center of Beijing Municipal Administration of Hospitals (XXZ02). The study sponsors had no role in the design and preparation of this manuscript.
Ethics approval and inform consent
This study was conducted based on the open-source data from the TCGA database. The patients involved in this database have already obtained ethical approval. All users were permitted to download relevant data for research and publish relevant articles. Therefore, there are no ethical issues and other conflicts of interest regarding this article.
Consent for publication
The authors have no competing interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1:
Fig. S1. Somatic mutation and expression of cuproptosis markers in cancers. Fig. S2. GDSC data revealing the sensitivity to different drugs of cuproptosis markers. Fig. S3. The correlation plot exhibiting the correlations between cuproptosis markers and immune scores in 21 TCGA cancers. Fig. S4. The correlation plot exhibiting the correlations between cuproptosis markers and stromal scores in 21 TCGA cancers. Fig. S5. The correlation plot exhibiting the correlations between cuproptosis markers and TMB in 21 TCGA cancers. Fig. S6. The correlation plot exhibiting the correlations between cuproptosis markers and MSI in 21 TCGA cancers. Fig. S7. NMF clustering of different cancers based on cuproptosis markers. Fig. S8. Visualization of NMF subtypes via two-dimensional t-SNE. Fig. S9. Differences in biological features between different NMF subgroups in pan-cancer. Fig. S10. Boxplots showing the hallmark scores of Hypoxia and Reactive oxygen species among different subgroups. Fig. S11. Prognosis analysis comparing different NMF subgroups in pan-cancer. Fig. S12. Analysis on different biological features between different NMF subgroups in COAD patients. Fig. S13. Analysis on different biological features between different NMF subgroups in OV patients. Fig. S14. Analysis on different biological features between different NMF subgroups in SARC patients.
Additional file 2:
Table S1. Prognostic value of CuAGs in pan-cancer. Table S2. Pearson correlation between CuAGs expressions and immune scores in pan-cancer. Table S3. Pearson correlation between CuAGs expressions and stromal scores in pan-cancer. Table S4. Pearson correlation between CuAGs expressions and TMB in pan-cancer. Table S5. Pearson correlation between CuAGs expressions and MSI in pan-cancer. Table S6. An ultimate cancer classification system based on distinct regulatory patterns of CuAGs. Table S7. The gene lists used for marking DNA repairing-related pathways. Table S8. The gene lists used for marking TME-related cells. Table S9. The gene lists used for marking TME-related pathways.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Bu, F., Li, X., Zhao, Y. et al. Pan-cancer patterns of cuproptosis markers reveal biologically and clinically relevant cancer subtypes. Biomark Res 11, 13 (2023). https://doi.org/10.1186/s40364-022-00446-5
- Tumor microenvironment