Identification of cuproptosis-related subtypes, construction of a prognosis model, and tumor microenvironment landscape in gastric cancer.

Cuproptosis is a novel identified regulated cell death (RCD), which is correlated with the development, treatment response and prognosis of cancer. However, the potential role of cuproptosis-related genes (CRGs) in the tumor microenvironment (TME) of gastric cancer (GC) remains unknown. Transcriptome profiling, somatic mutation, somatic copy number alteration and clinical data of GC samples were downloaded from the Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) database to describe the alterations of CRGs from genetic and transcriptional fields. Differential, survival and univariate cox regression analyses of CRGs were carried out to investigate the role of CRGs in GC. Cuproptosis molecular subtypes were identified by using consensus unsupervised clustering analysis based on the expression profiles of CRGs, and further analyzed by GO and KEGG gene set variation analyses (GSVA). Genes in distinct molecular subtypes were also analyzed by GO and KEGG gene enrichment analyses (GSEA). Differentially expressed genes (DEGs) were screened out from distinct molecular subtypes and further analyzed by GO enrichment analysis and univariate cox regression analysis. Consensus clustering analysis of prognostic DEGs was performed to identify genomic subtypes. Next, patients were randomly categorized into the training and testing group at a ratio of 1:1. CRG Risk scoring system was constructed through logistic least absolute shrinkage and selection operator (LASSO) cox regression analysis, univariate and multivariate cox analyses in the training group and validated in the testing and combined groups. Real-time quantitative polymerase chain reaction (RT-qPCR) was used to evaluate the expression of key Risk scoring genes. Sensitivity and specificity of Risk scoring system were examined by using receiver operating characteristic (ROC) curves. pRRophetic package in R was used to investigate the therapeutic effects of drugs in high- and low- risk score group. Finally, the nomogram scoring system was developed to predict patients' survival through incorporating the clinicopathological features and CRG Risk score. Most CRGs were up-regulated in tumor tissues and showed a relatively high mutation frequency. Survival and univariate cox regression analysis revealed that LIAS and FDX1 were significantly associated with GC patients' survival. After consensus unsupervised clustering analysis, GC patients were classified into two cuproptosis molecular subtypes, which were significantly associated with clinical features (gender, age, grade and TNM stage), prognosis, metabolic related pathways and immune cell infiltration in TME of GC. GO enrichment analyses of 84 DEGs, obtained from distinct molecular subtypes, revealed that DEGs primarily enriched in the regulation of metabolism and intracellular/extracellular structure in GC. Univariate cox regression analysis of 84 DEGs further screened out 32 prognostic DEGs. According to the expression profiles of 32 prognostic DEGs, patients were re-classified into two gene subtypes, which were significantly associated with patients' age, grade, T and N stage, and survival of patients. Nest, the Risk score system was constructed with moderate sensitivity and specificity. A high CRG Risk score, characterized by decreased microsatellite instability-high (MSI-H), tumor mutation burden (TMB) and cancer stem cell (CSC) index, and high stromal and immune score in TME, indicated poor survival. Four of five key Risk scoring genes expression were dysregulated in tumor compared with normal samples. Moreover, CRG Risk score was greatly related with sensitivity of multiple drugs. Finally, we established a highly accurate nomogram for promoting the clinical applicability of the CRG Risk scoring system. Our comprehensive analysis of CRGs in GC demonstrated their potential roles in TME, clinicopathological features, and prognosis. These findings may improve our understanding of CRGs in GC and provide new perceptions for doctors to predict prognosis and develop more effective and personalized therapy strategies.

Wang J ，Qin D ，Tao Z ，Wang B ，Xie Y ，Wang Y ，Li B ，Cao J ，Qiao X ，Zhong S ，Hu X ... -  《Frontiers in Immunology》

Cuproptosis-related risk score predicts prognosis and characterizes the tumor microenvironment in colon adenocarcinoma.

Cuproptosis is a novel copper-dependent regulatory cell death (RCD), which is closely related to the occurrence and development of multiple cancers. However, the potential role of cuproptosis-related genes (CRGs) in the tumor microenvironment (TME) of colon adenocarcinoma (COAD) remains unclear. Transcriptome, somatic mutation, somatic copy number alteration and the corresponding clinicopathological data of COAD were downloaded from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus database (GEO). Difference, survival and correlation analyses were conducted to evaluate the characteristics of CRGs in COAD patients. Consensus unsupervised clustering analysis of CRGs expression profile was used to classify patients into different cuproptosis molecular and gene subtypes. TME characteristics of different molecular subtypes were investigated by using Gene set variation analysis (GSVA) and single sample gene set enrichment analysis (ssGSEA). Next, CRG Risk scoring system was constructed by applying logistic least absolute shrinkage and selection operator (LASSO) cox regression analysis and multivariate cox analysis. Real-time quantitative polymerase chain reaction (RT-qPCR) and immunohistochemistry (IHC) were used to exam the expression of key Risk scoring genes. Our study indicated that CRGs had relatively common genetic and transcriptional variations in COAD tissues. We identified three cuproptosis molecular subtypes and three gene subtypes based on CRGs expression profile and prognostic differentially expressed genes (DEGs) expression profile, and found that changes in multilayer CRGs were closely related to the clinical characteristics, overall survival (OS), different signaling pathways, and immune cell infiltration of TME. CRG Risk scoring system was constructed according to the expression of 7 key cuproptosis-related risk genes (GLS, NOX1, HOXC6, TNNT1, GLS, HOXC6 and PLA2G12B). RT-qPCR and IHC indicated that the expression of GLS, NOX1, HOXC6, TNNT1 and PLA2G12B were up-regulated in tumor tissues, compared with those in normal tissues, and all of GLS, HOXC6, NOX1 and PLA2G12B were closely related with patient survival. In addition, high CRG risk scores were significantly associated with high microsatellite instability (MSI-H), tumor mutation burden (TMB), cancer stem cell (CSC) indices, stromal and immune scores in TME, drug susceptibility, as well as patient survival. Finally, a highly accurate nomogram was constructed to promote the clinical application of the CRG Risk scoring system. Our comprehensive analysis showed that CRGs were greatly associated with TME, clinicopathological characteristics, and prognosis of patient with COAD. These findings may promote our understanding of CRGs in COAD, providing new insights for physicians to predict prognosis and develop more precise and individualized therapy strategies.

Wang J ，Tao Z ，Wang B ，Xie Y ，Wang Y ，Li B ，Cao J ，Qiao X ，Qin D ，Zhong S ，Hu X ... -  《Frontiers in Oncology》

Identification of a novel cuproptosis-related gene signature and integrative analyses in patients with lower-grade gliomas.

Cuproptosis is a newly discovered unique non-apoptotic programmed cell death distinguished from known death mechanisms like ferroptosis, pyroptosis, and necroptosis. However, the prognostic value of cuproptosis and the correlation between cuproptosis and the tumor microenvironment (TME) in lower-grade gliomas (LGGs) remain unknown. In this study, we systematically investigated the genetic and transcriptional variation, prognostic value, and expression patterns of cuproptosis-related genes (CRGs). The CRG score was applied to quantify the cuproptosis subtypes. We then evaluated their values in the TME, prognostic prediction, and therapeutic responses in LGG. Lastly, we collected five paired LGG and matched normal adjacent tissue samples from Sun Yat-sen University Cancer Center (SYSUCC) to verify the expression of signature genes by quantitative real-time PCR (qRT-PCR) and Western blotting (WB). Two distinct cuproptosis-related clusters were identified using consensus unsupervised clustering analysis. The correlation between multilayer CRG alterations with clinical characteristics, prognosis, and TME cell infiltration were observed. Then, a well-performed cuproptosis-related risk model (CRG score) was developed to predict LGG patients' prognosis, which was evaluated and validated in two external cohorts. We classified patients into high- and low-risk groups according to the CRG score and found that patients in the low-risk group showed significantly higher survival possibilities than those in the high-risk group (P<0.001). A high CRG score implies higher TME scores, more significant TME cell infiltration, and increased mutation burden. Meanwhile, the CRG score was significantly correlated with the cancer stem cell index, chemoradiotherapy sensitivity-related genes and immune checkpoint genes, and chemotherapeutic sensitivity, indicating the association with CRGs and treatment responses. Univariate and multivariate Cox regression analyses revealed that the CRG score was an independent prognostic predictor for LGG patients. Subsequently, a highly accurate predictive model was established for facilitating the clinical application of the CRG score, showing good predictive ability and calibration. Additionally, crucial CRGs were further validated by qRT-PCR and WB. Collectively, we demonstrated a comprehensive overview of CRG profiles in LGG and established a novel risk model for LGG patients' therapy status and prognosis. Our findings highlight the potential clinical implications of CRGs, suggesting that cuproptosis may be the potential therapeutic target for patients with LGG.

Bao JH ，Lu WC ，Duan H ，Ye YQ ，Li JB ，Liao WT ，Li YC ，Sun YP ... -  《Frontiers in Immunology》

Cuproptosis-Related Signature Predicts the Prognosis, Tumor Microenvironment, and Drug Sensitivity of Hepatocellular Carcinoma.

Copper (Cu) metabolism is strongly associated with liver disease. Cuproptosis is a novel format of cell death, and cuproptosis-related genes (CRGs) were identified. However, the role of CRGs in Hepatocellular Carcinoma (HCC) remains unknown. The mRNA transcriptome profiling data, somatic mutation data, and copy number gene level data of The Cancer Genome Atlas-Liver Hepatocellular Carcinoma project (TCGA-LIHC) were downloaded for subsequent analysis. Molecular characterization analysis of CRGs, including differential gene expression analysis, mutation analysis, copy number variation (CNV) analysis, Kaplan-Meier analysis, and immune regulator prioritization analysis, was implemented. The nonnegative matrix factorization (NMF) approach was used to identify the CRG-related molecular subtypes. Principal component analysis was adopted to verify the robustness and reliability of the molecular subtype. The least absolute shrinkage and selection operator regression analysis was performed to construct the prognostic signature based on differentially expressed genes between molecular subtypes. The survival characteristics of the molecular subtype and the signature were analyzed. The Gene Set Variation Analysis was performed for functional annotation. The immune landscape analysis, including immune checkpoint gene analysis, single sample gene set enrichment analysis, tumor immune dysfunction and exclusion (TIDE) analysis, immune infiltration cell, and tumor mutation burden analysis (TMB), was conducted. The ability of the signature to predict conventional anti-HCC agent responses was evaluated. The signature was validated in the LIRI-JP cohort and the IMvigor210 cohort. A total of 13 CRGs are differentially expressed between the tumor and normal samples, while the mutation of CRGs in HCC is infrequent. The expression of CRGs is associated with the CNV level. Fourteen CRGs are associated with the prognosis of HCC. Two clusters were identified and HCC patients were divided into 2 groups with a cutoff risk score value of 1.570. HCC patients in the C1 cluster and high-risk have a worse prognosis. The area under the receiver operating characteristic curve for predicting 1-, 2-, and 3-year overall survival is 0.775, 0.768, and 0.757 in the TCGA-LIHC cohort, and 0.811, 0.741, and 0.775 in the LIRI-JP cohort. Multivariate Cox regression analysis indicates that the signature is an independent prognostic factor. Pathways involved in metabolism and gene stability and immune infiltration cells are significantly enriched. Immune checkpoint genes are highly expressed in the C1 cluster. TMB is positively correlated with the risk score. HCC patients in the high-risk group are more likely to benefit from conventional anti-HCC agents and immune checkpoint inhibitor therapies. The molecular characterization of CRGs in HCC is presented in this study, and a successful prognostic signature for HCC based on the cuproptosis-related molecular subtype was constructed.

Qi X ，Guo J ，Chen G ，Fang C ，Hu L ，Li J ，Zhang C ... -  《-》

Comprehensive analysis of a cuproptosis-related ceRNA network implicates a potential endocrine therapy resistance mechanism in ER-positive breast cancer.

While adjuvant endocrine therapy (ET) may decrease the mortality rate of estrogen receptor-positive (ER+) breast cancer (BC), the likelihood of relapse and metastasis due to ET resistance remains high. Cuproptosis is a recently discovered regulated cell death (RCD), whose role in tumors has yet to be elucidated. Thus, there is a need to study its specific regulatory mechanism in resistance to ET in BC, to identify novel therapeutic targets. The prognostic cuproptosis-related genes (CRGs) in ER+ BC were filtered by undergoing Cox regression and least absolute shrinkage and selection operator (LASSO) regression analyses in TCGA-BRCA, and a CRGs risk signature was constructed using the correlation coefficient. Immune infiltration analysis, immune function analysis, tumor microenvironment (TME) analysis, immune checkpoint analysis, immunotherapy response analysis, drug sensitivity analysis, and pathway activation analysis were carried out among the high- and low-risk groups in turn. The central CRG of cuproptosis in ER+ BC resistance to ET was acquired through the intersection of protein interaction network (PPI) analysis, genes differentially expressed (DEGs) between human BC cells LCC9 and MCF-7 (GSE159968), and CRGs with prognostic significance in TCGA-BRCA ER+ BC. The miRNAs upstream of the core CRGs were predicted based on the intersection of 4 databases, miRDB, RNA22, miRWalk, and RNAlnter. Candidate miRNAs consisted of the intersection of predicted miRNAs and miRNAs differentially expressed in the LCC9 and MCF-7 cell lines (GSE159979). Candidate lncRNAs were the intersection of the differential lncRNAs from the LCC9 and MCF-7 cell lines and the survival-related lncRNAs obtained from a univariate Cox regression analysis. Pearson's correlation analysis was performed between mRNA-miRNA, miRNA-lncRNA, and mRNA-lncRNA expression separately. We constructed A risk signature of 4-CRGs to predict the prognosis of ER+ BC in TCGA-BRCA, a risk score = DLD*0.378 + DBT*0.201 + DLAT*0.380 + ATP7A*0.447 was used as the definition of the formula. There were significant differences between the high- and low-risk groups based on the risk score of 4-CRGs in aspects of immune infiltration, immune function, expression levels of immune checkpoint genes, and signaling pathways. DLD was determined to be the central CRG of cuproptosis in ER+ BC resistance to ET through the intersection of the PPI network analysis, DEGs between LCC9 and MCF-7 and 4-CRGs. Two miRNAs hsa-miR-370-3p and hsa-miR-432-5p were found taking DLD mRNA as a target, and the lncRNA C6orf99 has been hypothesized to be a competitive endogenous RNA that regulates DLD mRNA expression by sponging off hsa-miR-370-3p and hsa-miR-432-5p. This study built a prognostic model based on genes related to cuproptosis in ER+ BC. We considered DLD to be the core gene associated with resistance to ET in ER+ BC via copper metabolism. The search for promising therapeutic targets led to the establishment of a cuproptosis-related ceRNA network C6orf99/hsa-miR-370-3p and hsa-miR-432-5p/DLD.

Zhang D ，Lu W ，Zhuo Z ，Wang Y ，Zhang W ，Zhang M ... -  《BMC Medical Genomics》