A Cuproptosis-Related LncRNA Risk Model for Predicting Prognosis and Immunotherapeutic Efficacy in Patients with Hepatocellular Carcinoma.
Cuproptosis is a novel programmed cell death pathway that is initiated by direct binding of copper to lipoylated tricarboxylic acid (TCA) cycle proteins. Recent studies have demonstrated that cuproptosis-related genes regulate tumorigenesis. However, the potential role and clinical significance of cuproptosis-related long noncoding RNAs (lncRNAs) in hepatocellular carcinoma (HCC) have not been established. We performed a bioinformatics analyses of RNA-sequencing data of HCC patients extracted from The Cancer Genome Atlas (TCGA) dataset to identify and validate a cuproptosis-related lncRNA prognostic signature. Furthermore, we analyzed the clinical significance of the prognostic signature of cuproptosis-related lncRNA in predicting the immunotherapeutic efficacy and the status of the tumor immune microenvironment. The RNA-sequencing data, genomic mutations, and clinical information were downloaded for 374 HCC samples and 50 normal liver samples from TCGA-Liver Hepatocellular Carcinoma (TCGA-LIHC) dataset. Co-expression analysis of Gene-lncRNA pairs with 49 known cuproptosis-related prognostic genes was used to define cuproptosis-related prognostic lncRNAs. We performed the LASSO algorithm and univariate and multivariate Cox regression analysis, respectively, to gradually identify the prognostic risk models of cuproptosis-related lncRNA based on the TCGA-LIHC dataset. Subsequently, the predictive performance of the model was evaluated using receiver operation characteristic (ROC) curves, Kaplan-Meier survival curves, and prognostic nomogram. The analysis of gene-lncRNA co-expression with 49 known cuproptosis-related genes identified 1359 cuproptosis-related lncRNAs in the TCGA-LIHC data set. A prognostic model was constructed with nine cuproptosis-related prognostic lncRNAs (AC007998.3, AC003086.1, AC009974.2, IQCH-AS1, LINC0256 1, AC105345.1, ZFPM2-AS1, AL353708.1 and WAC-AS1) using LASSO regression and Cox regression analyses. Risk scores were calculated for all HCC patient samples based on the four cuproptosis-related lncRNA prognostic models. All HCC patients were divided into high-risk and low-risk subgroups according to a 1:1 ratio column. The Kaplan-Meier survival curve analysis showed that the overall survival rate (OS) of the high-risk group patients was significantly lower than that of the low-risk group. The principal component analysis (PCA) confirmed that the prognostic lncRNA model accurately distinguished between high- and low-risk HCC patients. Furthermore, regression analysis as well as ROC curves confirmed the prognostic value of the risk score. A nomogram with risk scores and other clinicopathological characteristics was constructed. The nomogram accurately predicted the probability of 1-, 3-, and 5-year OS in HCC patients. Tumor mutation burden (TMB) scores were higher for high-risk patients than for low-risk patients. HCC patients in the low-risk group showed lower TIDE scores and greater sensitivity to antitumor drugs than those in the high-risk group. Tumor immune responses and tumor immune cell infiltration were significantly different between the high-risk and low-risk groups of patients with HCC. Our study identified a 9-cuproptosis-related lncRNA signature that accurately predicted prognosis, immunotherapeutic efficacy, and the status of the tumor immune microenvironment in HCC patients. Therefore, this cuproptosis-related lncRNA risk model is a potential prognostic biometric feature in HCC and shows high clinical value in identifying HCC patients who are potentially responsive to immunotherapy.
Wang S
,Bai H
,Fei S
,Miao B
... -
《-》
Cuproptosis-related lncRNA scoring system to predict the clinical outcome and immune landscape in pancreatic adenocarcinoma.
Cuproptosis is a recently discovered novel programmed cell death pathway that differs from traditional programmed cell death and has an important role in cancer and immune regulation. Long noncoding RNA (lncRNA) is considered new potential prognostic biomarkers in pancreatic adenocarcinoma (PAAD). However, the prognostic role and immune landscape of cuproptosis-related lncRNA in PAAD remain unclear. The transcriptome and clinical data of PAAD were obtained from The Cancer Genome Atlas (TCGA) database. Cuproptosis-related lncRNA was identified using Pearson correlation analysis. The optimal lncRNA was screened by Cox and the Least Absolute Shrinkage and Selection Operator (LASSO) regression mode, and for the construction of risk scoring system. PAAD patients were divided into high- and low-risk groups according to the risk score. Clinicopathological parameter correlation analysis, univariate and multivariate Cox regression, time-dependent receiver operating characteristic (ROC) curves, and nomogram were performed to evaluate the model. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were used to explore differences in biological function between different risk groups. Single-sample gene set enrichment analysis (ssGSEA) and Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data (ESTIMATE) algorithm were used to analyze the differences in tumor immune microenvironment (TIME) in different risk groups of PAAD. Additionally, the Tumor Immune Dysfunction and Exclusion (TIDE) algorithm was used to predict immunotherapy response and identify potential immune beneficiaries. Immune checkpoints and tumor mutation burden (TMB) were also systematically analyzed. Finally, drug sensitivity analysis was used to explore the reactivity of different drugs in high- and low-risk groups to provide a reference for the selection of precise therapeutic drugs. Six cuproptosis-related lncRNAs (AL117335.1, AC044849.1, AL358944.1, ZNF236-DT, Z97832.2, and CASC8) were used to construct risk model. Survival analysis showed that overall survival and progression-free survival in the low-risk group were better than those in the high-risk group, and it is suitable for PAAD patients with different clinical characteristics. Univariate and multifactorial Cox regression analysis showed that risk score was an independent prognostic factor in PAAD patients. ROC analysis showed that the AUC values of the risk score in 1 year, 3 years and 5 years were 0.707,0.762 and 0.880, respectively. Nomogram showed that the total points of PAAD patients at 1 year, 3 years, and 5 years were 0.914,0.648, and 0.543. GO and KEGG analyses indicated that the differential genes in the high- and low-risk groups were associated with tumor proliferation and metastasis and immune regulatory pathway. Immune correlation analysis showed that the amount of pro-inflammatory cells, including CD8+ T cells, was significantly higher in the low-risk group than in the high-risk group, and the expression of immune checkpoint genes, including PD-1 and CTLA-4, was increased in the low-risk group. TIDE analysis suggests that patients in the low-risk group may benefit from immunotherapy. Finally, there was significant variability in multiple chemotherapeutic and targeted drugs across the risk groups, which informs our clinical drug selection. Our cuproptosis-related lncRNA scoring system (CRLss) could predict the clinical outcome and immune landscape of PAAD patients, identify the potential beneficiaries of immunotherapy, and provide a reference for precise therapeutic drug selection.
Huang Y
,Gong P
,Su L
,Zhang M
... -
《Scientific Reports》