Cuproptosis-related lncRNA scoring system to predict the clinical outcome and immune landscape in pancreatic adenocarcinoma

Abstract

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.