Establishment of a 7-gene prognostic signature based on oxidative stress genes for predicting chemotherapy resistance in pancreatic cancer

Oxidative stress is involved in regulating various biological processes in human cancers. However, the effect of oxidative stress on pancreatic adenocarcinoma (PAAD) remained unclear. Pancreatic cancer expression profiles from TCGA were downloaded. Consensus ClusterPlus helped classify molecular subtypes based on PAAD prognosis-associated oxidative stress genes. Limma package filtered differentially expressed genes (DEGs) between subtypes. A multi-gene risk model was developed using Lease absolute shrinkage and selection operator (Lasso)-Cox analysis. A nomogram was built based on risk score and distinct clinical features. Consistent clustering identified 3 stable molecular subtypes (C1, C2, C3) based on oxidative stress-associated genes. Particularly, C3 had the optimal prognosis with the greatest mutation frequency, activate cell cycle pathway in an immunosuppressed status. Lasso and univariate cox regression analysis selected 7 oxidative stress phenotype-associated key genes, based on which we constructed a robust prognostic risk model independent of clinicopathological features with stable predictive performance in independent datasets. High-risk group was found to be more sensitive to small molecule chemotherapeutic drugs including Gemcitabine, Cisplatin, Erlotinib and Dasatinib. The 6 of 7 genes expressions were significantly associated with methylation. Survival prediction and prognostic model was further improved through a decision tree model by combining clinicopathological features with RiskScore. The risk model containing seven oxidative stress-related genes may have a greater potential to assist clinical treatment decision-making and prognosis determination.