A Novel Ferroptosis-Related Gene Signature Predicts Recurrence in Patients With Pancreatic Ductal Adenocarcinoma

Abstract

Background: Recurrence after surgery is largely responsible for the extremely poor outcomes for patients with pancreatic ductal adenocarcinoma (PDAC). Ferroptosis is implicated in chemotherapy sensitivity and tumor recurrence, we aimed to find out survival-associated ferroptosis-related genes and use them to build a practical risk model with the purpose to predict PDAC recurrence.

Methods: Univariate Cox regression analysis was conducted to obtain prognostic ferroptosis-related genes in The Cancer Genome Atlas (TCGA, N = 140) cohort. Multivariate Cox regression analysis was employed to construct a reliable and credible gene signature. The prognostic performance was verified in a MTAB-6134 (N = 286) validation cohort and a PACA-CA (N = 181) validation cohort. The stability of the signature was tested in TCGA and MTAB-6134 cohorts by ROC analyses. Pathway enrichment analysis was adopted to preliminary illuminate the biological relevance of the gene signature.

Results: Univariate and multivariate Cox regression analyses identified a 5-gene signature that contained CAV1, DDIT4, SLC40A1, SRXN1 and TFAP2C. The signature could efficaciously stratify PDAC patients with different recurrence-free survival (RFS), both in the training and validation cohorts. Results of subgroup receiver operating characteristic curve (ROC) analyses confirmed the stability and the independence of this signature. Our signature outperformed clinical indicators and previous reported models in predicting RFS. Moreover, the signature was found to be closely associated with several cancer-related and drug response pathways.

Conclusion: This study developed a precise and concise prognostic model with the clinical implication in predicting PDAC recurrence. These findings may facilitate individual management of postoperative recurrence in patients with PDAC.

Introduction

Pancreatic ductal adenocarcinoma (PDAC) is one of the most aggressive gastrointestinal tumors with a 5-years survival rate not exceeding 9% (Siegel et al., 2019). Surgery combined with adjuvant chemotherapy is the standard treatment for resectable PDAC (Garrido-Laguna and Hidalgo, 2015), and it significantly elevates the 5-years survival rate to 20–25% (Hackert et al., 2016). Unfortunately, majority of PDAC patients miss the chance for this regimen owing to atypical and unspecific symptoms at an early, resectable stage (Heinemann and Boeck, 2008; Sharma et al., 2011; Rombouts et al., 2015). Furthermore, even after radical resection, most patients will develop recurrence (Vincent et al., 2011), and patients with tumor recurrence have significantly decreased overall survival rates compared with patients without tumor recurrence (Breidert et al., 2012). Under this grim circumstance, novel development of recurrence risk prediction models is urgently needed for clinicians to manage patients by more tailored therapeutic strategies and postoperative surveillance.

Ferroptosis is an iron-dependent form of non-apoptotic cell death, which is featured with lipid peroxidation (Dixon et al., 2012a). Triggering ferroptosis has shown potential therapeutic value in oncology (Shan et al., 2020), especially for eradicating aggressive malignancies that are refractory to conventional treatments (Roh et al., 2017; Shin et al., 2018; Gao et al., 2019). In PDAC, activation of ferroptosis by knockdown ARF6 gene could enhance the sensitivity to gemcitabine (Ye et al., 2020), the cornerstone drug of chemotherapy for PDAC (Heinemann, 2001). Anticancer effects of ferroptosis inducers on PDAC have also been reported (Eling et al., 2015). In addition, Tang et al. found that diminished ferroptosis sensitivity is associated with decreased immune activity and robust gemcitabine resistance in PDAC (Tang et al., 2020). These previous studies revealed that ferroptosis represents an under-explored source of prognostic factors that could be used to predict tumor recurrence.

In this study, we analyzed the correlation between ferroptosis-related genes with recurrence-free survival (RFS) of patients, and developed a 5-gene risk prediction model with satisfactory predictive performance in both training and validation datasets. The proposed signature outperformed clinical parameters and previous reported models in predicting tumor recurrence. In addition, this signature is positively correlated with cell division and DNA replication. The 5-gene signature is promising to help to make individual treatments and has potential to be implemented for clinical practice.

Materials and Methods

Data Collection

The three PDAC cohorts included in this study for survival analyses were the MTAB-6134 cohort (N = 286), PACA-CA cohort (N = 181) and TCGA cohort (N = 140). TCGA cohort was used as the training set, while MTAB-6134 cohort and PACA-CA cohort were adopted for external validation. For samples in TCGA cohort, gene expression data and corresponding clinical data were downloaded from the TCGA hub at UCSC Xena (https://tcga.xenahubs.net). Microarray data and clinical information of MTAB-6134 cohort was obtained from ArrayExpress (https://www.ebi.ac.uk/arrayexpress/) database. In the cases of PACA-CA cohort, normalized RNA-sequencing data and survival data were downloaded from the International Cancer Genome Consortium (ICGC, https://icgc.org/) database. For TCGA cohort expression profiles, the gene expression profile was measured experimentally using the Illumina HiSeq 2000 RNA Sequencing platform by the University of North Carolina TCGA genome characterization center, and the gene-level transcription was estimated by log2(x+1) transformed RSEM normalized count. For MTAB-6134, Affymetrix Human Genome U219 Array were performed to annotate transcription profiles, and the expression levels can be fetched, too. For PACA-CA, ICGC portal provided us with convenient access to fetch normalized expression profiles. Patients whose clinical information was incomplete or whose histopathological type was not PDAC were strictly removed. Clinical RFS data was available in all of these three cohorts, and patients with a RFS < 1 month were excluded. In addition, we also downloaded gene expression data of GSE15471 cohort and GSE16515 cohort from the Gene Expression Omnibus (GEO, https://www.ncbi.nlm.nih.gov/ geo/) database, aiming to assess the distribution of risk score in PDAC tissues and matched adjacent normal tissues.

Construction of the Ferroptosis-Related Gene Signature for Prognostic Prediction

260 ferroptosis-related gene were retrieved and obtained from the FerrDb database (http://www.zhounan.org/ferrdb) (Zhou and Bao, 2020). In the training TCGA cohort, ferroptosis-related gene that were significantly associated with the RFS of PDAC patients were screened through using the univariate Cox regression analysis. The following multivariate Cox regression analysis was employed to select an optimal prognostic model with the minimum Akaike Information Criterion value. Based on the expression level and corresponding coefficient of each prognostic gene derived from the multivariate Cox regression analysis, the risk score of every patient in all of three cohorts was computed.

Prognostic Performance of the Ferroptosis-Related Gene Signature

According to the best cut-off value determined by the X-tile software (Camp et al., 2004), patients in both the training and validation cohorts were assigned to a low- or a high-risk group. The differences in OS time between the two groups were estimated by the Kaplan–Meier survival curves. To determine whether our signature provided improved survival prediction, ROC analyses were performed on the ferroptosis-related gene signature, clinical indicators and two previously reported gene signatures predicting PDAC recurrence (Shi et al., 2017; Kim et al., 2019a).

Functional Enrichment Analyses

Genes correlated with the risk score were identified using Spearman correlation analysis in MTAB-6134 cohort (p < 0.05) and TCGA cohort (p < 0.05). The top 1,000 positively correlated genes were selected and then submitted to the Metascape database (Zhou et al., 2019) for function annotation and pathway enrichment.

Statistical Analysis

The statistical analysis and graphical work were finished in the R environment (version 3.5.2). Cox regression analyses and K-M survival curves were conducted by the “survival” package. The ROC curves for prognosis were plotted by the “survivalROC” package. Boxplots were generated from the “ggpubr” package. Calibration curves were derived from the “rms” package. A two-sided log-rank p < 0.05 was considered significant.

Results

Construction of the Prognostic Ferroptosis-Related Signature for PDAC

The univariate Cox regression analysis screened 21 of 260 ferroptosis-related genes that were significantly associated with RFS in the TCGA training cohort (p < 0.05). Subsequently, the multivariate Cox regression analysis was implemented to generate a prognostic model. Finally, five ferroptosis-related genes with high correlation with RFS were selected to make up the signature, and the HRs, 95%CIs, and p values of these five genes were illustrated in Figure 1A. Each of the five genes could effectively stratify PDAC patients with different RFS time in the TCGA training cohort (p < 0.05, Figures 1B–F). Among the five genes, CAV1, DDIT4, SRXN1 and TFAP2C were indicative of better survival while SLC40A1 was associated with worse survival. According to the expression of these five genes, we developed a risk-score formula: Risk score = (0.224696×expression value of CAV1)+(0.265016×expression value of DDIT4)-(0.24004×expression value of SLC40A1)+(0.934599×expression value of SRXN1)+(0.144059×expression value of TFAP2C).

FIGURE 1

Distribution of the Risk Score

We first calculated the risk score of samples in GSE15471 and GSE16515 cohorts. These two cohorts contained PDAC and paired adjacent normal samples. As shown in Figures 2A,B, PDAC tissues exhibited a remarkably elevated risk score compared with adjacent normal tissues, suggesting the hazardous role of risk score. Meanwhile, we observed that the risk score increased with the increase of histological grade in MTAB-6134 cohort (Figure 2C). This finding indicated that risk score is associated with tumor malignancy. According to the optimal cut-off value of risk score computed by X-tile software, patients were divided into high- and low-risk groups. The distribution of the risk scores, survival status, and heatmap for five genes in both training and validation samples were illustrated in Figures 2D–F. Low-risk-score patients had a dramatically decreased mortality rate compared with high-risk-score patients in these three cohorts.

FIGURE 2

Prognostic Validation of the Ferroptosis-Related Signature

The Kaplan–Meier curve demonstrated that patients in high-risk group had a significantly shorter RFS than patients in low-risk group in both training and validation cohorts (Figures 3A–C). The calibration curves showed that the RFS predicted by the ferroptosis-related signature were in good accordance with the observed RFS (Figures 3D–F). ROC curves revealed that this signature had moderate accuracy for predicting 1-year RFS, as the area under the curve (AUC) value was 0.787, 0.686 and 0.652 in three cohorts, respectively (Figures 3G–I). In addition, the AUC value of our signature was higher than that of two previous reported gene signatures with recurrence predictive ability. This finding indicated that the proposed ferroptosis-related signature could provide improved outcome prediction.

FIGURE 3

Subgroup Analyses of the Ferroptosis-Related Signature

As we know, several clinical variables are important risk factors in PDAC recurrence. To verify the independence and stability of the ferroptosis-related signature, we performed the subgroup analyses through stratifying patients based on gender, histological grade, T stage and N stage. More than half of the patients in PACA-CA cohort have incomplete clinical information and thus TCGA and MTAB-6134 cohorts were selected for further analyses. As shown in Figure 4, no matter what the group was, our signature still had the capacity to predict prognosis with satisfactory accuracy, as the AUC values were no less than 0.75 in TCGA cohort and 0.65 in MTAB-6134 cohort. More importantly, the AUC value of our model was close to or not lower than that of the previous models in each subgroup.

FIGURE 4

Predictive Comparison of the Ferroptosis-Related Signature With Clinical Indicators

To improve the clinical relevance of this model, we compared the predictive accuracy of the ferroptosis-related signature with that of clinical parameters like grade, T stage and N stage. In TCGA cohort, this signature exhibited more powerful performance in predictive 6-months, 12-months and 18-months recurrence risk compared with traditional parameters (Figure 5A). Similar trends were observed in MTAB-6134 cohort (Figure 5B).

FIGURE 5

Enriched Functions and Pathways Related to the Ferroptosis-Related Signature

With the aim to preliminarily elucidate biological functions and pathways of the ferroptosis-related signature, we identified top 1,000 genes whose expression were correlated with risk score in TCGA cohort and MTAB-6134 cohort, respectively. The results showed that these genes were involved in several recurrence-related pathways like cell cycle, cell division, DNA replication and TP53-related pathways in two cohorts (Figures 6A,B).

FIGURE 6

Discussion

PDAC has an extremely poor outcome even after surgery. Majority of patients develop recurrence after surgical resection (Kim et al., 2019b) and high incidence of recurrence remains the main reason for the unfavorable prognosis of PDAC (Li et al., 2019). In this study, recurrence was observed in 65.7, 80.8 and 61.9% of patients in the TCGA, MTAB-6134 and PACA-CA cohorts, respectively. Accurate recurrence potential prediction of PDAC is of great significance for clinical practitioners to perform specific follow-up strategies and timely intervention of possible recurrence (Daamen et al., 2019). So far, several clinical features have been proved as predictors of recurrence in PDAC, such as resection margin (Ghaneh et al., 2019), tumor size (Ansari et al., 2017), tumor grade (Bilici, 2014) and serum carbohydrate antigen 19‐9 (CA19‐9) (Sugiura et al., 2012). However, few biomarkers and models considering genetic and genomic features of PDAC patients have been identified as postoperative predictors of recurrence in PDAC (Qian et al., 2018).

Ferroptosis, usually regarded as a specific form of regulated cell death, is dependent on the presence of intracellular iron and the accumulation of reactive oxygen species (ROS) (Dixon et al., 2012b). Ferroptosis is mediated by the Fenton reaction, in which Fe2+ reacts with hydrogen peroxide to generate ROS (Xie et al., 2016). ROS are not only a very important secondary signal in cells, but also can also impair the stability of DNA and promote cell death (Dixon and Stockwell, 2014). Variety of molecules have been shown to induce iron degradation and deposition in pancreatic cancer cells, suggesting some new options for the treatment of pancreatic cancer. For example, some researchers have revealed that nuclear receptor coactivator 4(NCOA4)-mediated ferritinophagy, as autophagic process, contributes to ferroptosis via the degradation of ferritin in pancreatic cancer. Degradation of ferritin results in increased intracellular free iron, which triggers ROS generation and consequent ferroptosis in pancreatic cancer (Mancias et al., 2014; Hou et al., 2016). Zhu et al. discovered that heat shock 70 kDa protein 5 (HSPA5) is a negative regulator of ferroptosis in pancreatic cancer cells. HSPA5 is regulated by activating transcription factor 4 (ATF4), and with the activation of ATF4, HSPA5 protein binding to GPX4 increases the stability of GPX4 to protect against ferroptosis. HSPA5-GPX4 pathway make great contribution to regulating ferroptosis in pancreatic cancer (Zhu et al., 2017). In addition to iron metabolism, ROS metabolism also plays a pivotal role in ferroptosis. For instance, ALOX15 can catalyze lipid hydroperoxide generation and ferroptosis in pancreatic cancer (Shintoku et al., 2017), while mitochondrial-targeted nitroxide holds the ability to prevent mitochondrial lipid oxidation on the contrary (Krainz et al., 2016).

In this study, we developed a ferroptosis-related gene signature to predictive recurrence risk in patients with PDAC. This signature exhibited moderate predictive ability with cross-cohort compatibility. In the process of building the signature, we initially performed univariate Cox regression analysis on 260 ferroptosis-related genes and screened 21 survival-related genes. Following multivariate Cox regression analysis identified a 5-gene signature. K-M survival curves confirmed that the signature could effectively capture the RFS differences between low-risk and high-risk groups. ROC analyses demonstrated that the signature had elevated AUC values compared with previously reported models and traditional predictors, indicating a imptoved survival prediction of our model.

Advances in gene chip and high-throughput sequencing have enabled novel development of reliable prognostic models for PDAC patients. To our surprise, there are many models that predict overall survival (Feng et al., 2020; Tan et al., 2020; Wu et al., 2020a; Xu et al., 2021), but few that predict RFS. Meanwhile, few studies have analyzed ferroptosis-related genes expression profiles to construct a risk prediction model. Ferroptosis, a novel form of cell death different from apoptosis, has the potential to promote the efficiency of apoptosis-inducing chemotherapeutic drugs in cell death induction. Thus, ferroptosis inducers are emerging as a promising solution to the problem of tumor drug resistance and combination of ferroptosis inducer and chemotherapy can achieve synergistic reaction and improve the sensitivity of chemotherapy (Stockwell et al., 2017; Shen et al., 2018; Skoupilová et al., 2018). Gemcitabine has become standard of care for postoperative PDAC worldwide and its sensitivity is a key factor affecting tumor recurrence (Liao et al., 2020). Previous studies have proved that the activation of ferroptosis could enhance gemcitabine sensitivity while diminished ferroptosis sensitivity is associated with robust gemcitabine resistance in PDAC (Tang et al., 2020; Ye et al., 2020), indicating that expression of ferroptosis-related gene is a practical tool to predict recurrence.

As for the specific function of these five genes, previous studies have demonstrated that CAV1, DDIT4, SRXN1 and TFAP2C were closely related to PDAC tumorigenesis and progression. CAV1 is overexpressed in PDAC tissues compared with normal tissues and its elevated expression indicates unfavorable prognosis. CAV-1 expression is positively correlated with disease stage, tumor grade and metastatic potential (Huang et al., 2012; Liu et al., 2014). Knockdown of CAV1 results in decreased cell proliferation, invasion, migration and chemotherapy resistance (Chatterjee et al., 2015). DDIT4 is highly expressed in PDAC cells, and DDIT4 inhibition leads to suppressed proliferation, migration, and invasion of PDAC cells (Ma et al., 2020). Expression level of SRXN1 in gemcitabine-resistant SW1990 cells was down-regulated to 6.84-fold change compared with its parental PDAC cell line SW 1990, but the biological role of SRXN1 in PDAC remains unclear (Zhou et al., 2015). Low TFAP2C protein expression is an independent adverse prognostic factor in PDAC (Xiong et al., 2018), while our study proved that high TFAP2C mRNA expression was an unfavorable prognostic indicator. SLC40A1 is a novel iron metabolism associated gene and it regulates iron metabolism to overcome cisplatin resistance in ovarian cancer (Wu et al., 2017; Wu et al., 2020b).

To be honest, there are still many limitations in current research. First, this signature was based on the retrospective data, and need to be verified in more prospective cohorts. Second, all cohorts used in current study have relatively small size (not exceeding 300 patients), partly due to the low opportunities for surgical resection. Our signature needs to be validated in more large-size cohorts in the future. Third, because of the limited number of samples in TCGA and MTAB-6134 cohorts, some subgroup analyses cannot be implemented. In addition, further experiments are needed to clarify the biological implications and relevant mechanisms of selected genes in PDAC progression.

In conclusion, this study is the first to examine ferroptosis-related gene expression profiles in PDAC patients and assess their relationship with tumor recurrence. The proposed signature showed satisfactory performance in predicting PDAC recurrence. The underlying mechanisms of the signature were also been investigated. These results are helpful for improving precision for treatment applications and promoting patient survival and life quality.

References

• 1

  AnsariD.BaudenM.BergströmS.RylanceR.Marko-VargaG.AnderssonR. (2017). Relationship between Tumour Size and Outcome in Pancreatic Ductal Adenocarcinoma. Br. J. Surg.104, 600–607. 10.1002/bjs.10471

  • CrossRef
  • Google Scholar

• 2

  BiliciA. (2014). Prognostic Factors Related with Survival in Patients with Pancreatic Adenocarcinoma. Wjg20, 10802–10812. 10.3748/wjg.v20.i31.10802

  • CrossRef
  • Google Scholar

• 3

  BreidertM.KeckT.MakowiecF.LohrmannC.HarderJ.FischerR. (2012). Early Recurrence of Pancreatic Cancer after Resection and during Adjuvant Chemotherapy. Saudi J. Gastroenterol.18, 118–121. 10.4103/1319-3767.93815

  • CrossRef
  • Google Scholar

• 4

  CampR. L.Dolled-FilhartM.RimmD. L. (2004). X-tile: a New Bio- Informatics tool for Biomarker Assessment and Outcome-Based cut-Point Optimization. Clin. Cancer Res.10, 7252–7259. 10.1158/1078-0432.Ccr-04-0713

  • CrossRef
  • Google Scholar

• 5

  ChatterjeeM.Ben-JosefE.ThomasD. G.MorganM. A.ZalupskiM. M.KhanG.et al (2015). Caveolin-1 Is Associated with Tumor Progression and Confers a Multi-Modality Resistance Phenotype in Pancreatic Cancer. Sci. Rep.5, 10867. 10.1038/srep10867

  • CrossRef
  • Google Scholar

• 6

  DaamenL. A.GrootV. P.IntvenM. P. W.BesselinkM. G.BuschO. R.KoerkampB. G.et al (2019). Postoperative Surveillance of Pancreatic Cancer Patients. Eur. J. Surg. Oncol.45, 1770–1777. 10.1016/j.ejso.2019.05.031

  • CrossRef
  • Google Scholar

• 7

  DixonS. J.StockwellB. R. (2014). The Role of Iron and Reactive Oxygen Species in Cell Death. Nat. Chem. Biol.10 (1), 9–17. 10.1038/nchembio.1416

  • CrossRef
  • Google Scholar

• 8

  DixonS. J.LembergK. M.LamprechtM. R.SkoutaR.ZaitsevE. M.GleasonC. E.et al (2012). Ferroptosis: an Iron-dependent Form of Nonapoptotic Cell Death. Cell149, 1060–1072. 10.1016/j.cell.2012.03.042

  • CrossRef
  • Google Scholar

• 9

  DixonS. J.LembergK. M.LamprechtM. R.SkoutaR.ZaitsevE. M.GleasonC. E.et al (2012). Ferroptosis: an Iron-dependent Form of Nonapoptotic Cell Death. Cell149 (5), 1060–1072. 10.1016/j.cell.2012.03.042

  • CrossRef
  • Google Scholar

• 10

  ElingN.ReuterL.HazinJ.Hamacher-BradyA.BradyN. R. (2015). Identification of Artesunate as a Specific Activator of Ferroptosis in Pancreatic Cancer Cells. Oncoscience2, 517–532. 10.18632/oncoscience.160

  • CrossRef
  • Google Scholar

• 11

  FengZ.ShiM.LiK.MaY.JiangL.ChenH.et al (2020). Development and Validation of a Cancer Stem Cell-Related Signature for Prognostic Prediction in Pancreatic Ductal Adenocarcinoma. J. Transl Med.18, 360. 10.1186/s12967-020-02527-1

  • CrossRef
  • Google Scholar

• 12

  GaoM.DengJ.LiuF.FanA.WangY.WuH.et al (2019). Triggered Ferroptotic Polymer Micelles for Reversing Multidrug Resistance to Chemotherapy. Biomaterials223, 119486. 10.1016/j.biomaterials.2019.119486

  • CrossRef
  • Google Scholar

• 13

  Garrido-LagunaI.HidalgoM. (2015). Pancreatic Cancer: from State-Of-The-Art Treatments to Promising Novel Therapies. Nat. Rev. Clin. Oncol.12, 319–334. 10.1038/nrclinonc.2015.53

  • CrossRef
  • Google Scholar

• 14

  GhanehP.KleeffJ.HalloranC. M.RaratyM.JacksonR.MellingJ.et al (2019). The Impact of Positive Resection Margins on Survival and Recurrence Following Resection and Adjuvant Chemotherapy for Pancreatic Ductal Adenocarcinoma. Ann. Surg.269, 520–529. 10.1097/sla.0000000000002557

  • CrossRef
  • Google Scholar

• 15

  HackertT.UlrichA.BüchlerM. W. (2016). Borderline Resectable Pancreatic Cancer. Cancer Lett.375, 231–237. 10.1016/j.canlet.2016.02.039

  • CrossRef
  • Google Scholar

• 16

  HeinemannV.BoeckS. (2008). Perioperative Management of Pancreatic Cancer. Ann. Oncol.19 (Suppl. 7), vii273–vii278. 10.1093/annonc/mdn450

  • CrossRef
  • Google Scholar

• 17

  HeinemannV. (2001). Gemcitabine: Progress in the Treatment of Pancreatic Cancer. Oncology60, 8–18. 10.1159/000055290

  • CrossRef
  • Google Scholar

• 18

  HouW.XieY.SongX.SunX.LotzeM. T.ZehH. J.3rdet al (2016). Autophagy Promotes Ferroptosis by Degradation of Ferritin. Autophagy12 (8), 1425–1428. 10.1080/15548627.2016.1187366

  • CrossRef
  • Google Scholar

• 19

  HuangC.QiuZ.WangL.PengZ.JiaZ.LogsdonC. D.et al (2012). A Novel FoxM1-Caveolin Signaling Pathway Promotes Pancreatic Cancer Invasion and Metastasis. Cancer Res.72, 655–665. 10.1158/0008-5472.Can-11-3102

  • CrossRef
  • Google Scholar

• 20

  KimJ.JoY. H.JangM.NguyenN. N. Y.YunH. R.KoS. H.et al (2019). PAC-5 Gene Expression Signature for Predicting Prognosis of Patients with Pancreatic Adenocarcinoma. Cancers11, 1749. 10.3390/cancers11111749

  • CrossRef
  • Google Scholar

• 21

  KimY. I.SongK. B.LeeY.-J.ParkK.-M.HwangD. W.LeeJ. H.et al (2019). Management of Isolated Recurrence after Surgery for Pancreatic Adenocarcinoma. Br. J. Surg.106, 898–909. 10.1002/bjs.11144

  • CrossRef
  • Google Scholar

• 22

  KrainzT.GaschlerM. M.LimC.SacherJ. R.StockwellB. R.WipfP. (2016). A Mitochondrial-Targeted Nitroxide Is a Potent Inhibitor of Ferroptosis. ACS Cent. Sci.2 (9), 653–659. 10.1021/acscentsci.6b00199

  • CrossRef
  • Google Scholar

• 23

  LiS.XuH.WuC.WangW.JinW.GaoH.et al (2019). Prognostic Value of γ-glutamyltransferase-to-albumin Ratio in Patients with Pancreatic Ductal Adenocarcinoma Following Radical Surgery. Cancer Med.8, 572–584. 10.1002/cam4.1957

  • CrossRef
  • Google Scholar

• 24

  LiaoW.HuangJ.ZhuG.ZhouJ.WenF.ZhangP.et al (2020). S-1 or Gemcitabine Adjuvant Therapy in Resected Pancreatic Cancer: a Cost-Effectiveness Analysis Based on the JASPAC-01 Trial. Expert Rev. Pharmacoeconomics Outcomes Res.20, 133–138. 10.1080/14737167.2020.1677155

  • CrossRef
  • Google Scholar

• 25

  LiuL.XuH.-X.WangW.-Q.WuC.-T.ChenT.QinY.et al (2014). Cavin-1 Is Essential for the Tumor-Promoting Effect of Caveolin-1 and Enhances its Prognostic Potency in Pancreatic Cancer. Oncogene33, 2728–2736. 10.1038/onc.2013.223

  • CrossRef
  • Google Scholar

• 26

  MaY.LiY.LingS.LiX.KongB.HuM.et al (2020). Loss of Heterozygosity for KrasG12D Promotes REDD1-dependent, Non-canonical Glutamine Metabolism in Pancreatic Ductal Adenocarcinoma. Biochem. Biophys. Res. Commun.526, 880–888. 10.1016/j.bbrc.2020.03.137

  • CrossRef
  • Google Scholar

• 27

  ManciasJ. D.WangX.GygiS. P.HarperJ. W.KimmelmanA. C. (2014). Quantitative Proteomics Identifies NCOA4 as the Cargo Receptor Mediating Ferritinophagy. Nature509 (7498), 105–109. 10.1038/nature13148

  • CrossRef
  • Google Scholar

• 28

  QianZ. R.RubinsonD. A.NowakJ. A.Morales-OyarvideV.DunneR. F.KozakM. M.et al (2018). Association of Alterations in Main Driver Genes with Outcomes of Patients with Resected Pancreatic Ductal Adenocarcinoma. JAMA Oncol.4, e173420. 10.1001/jamaoncol.2017.3420

  • CrossRef
  • Google Scholar

• 29

  RohJ.-L.KimE. H.JangH.ShinD. (2017). Nrf2 Inhibition Reverses the Resistance of Cisplatin-Resistant Head and Neck Cancer Cells to Artesunate-Induced Ferroptosis. Redox Biol.11, 254–262. 10.1016/j.redox.2016.12.010

  • CrossRef
  • Google Scholar

• 30

  RomboutsS. J. E.VogelJ. A.van SantvoortH. C.van LiendenK. P.van HillegersbergR.BuschO. R. C.et al (2015). Systematic Review of Innovative Ablative Therapies for the Treatment of Locally Advanced Pancreatic Cancer. Br. J. Surg.102, 182–193. 10.1002/bjs.9716

  • CrossRef
  • Google Scholar

• 31

  ShanX.LiS.SunB.ChenQ.SunJ.HeZ.et al (2020). Ferroptosis-driven Nanotherapeutics for Cancer Treatment. J. Controlled Release319, 322–332. 10.1016/j.jconrel.2020.01.008

  • CrossRef
  • Google Scholar

• 32

  SharmaC.EltawilK. M.RenfrewP. D.WalshM. J.MolinariM. (2011). Advances in Diagnosis, Treatment and Palliation of Pancreatic Carcinoma: 1990-2010. Wjg17, 867–897. 10.3748/wjg.v17.i7.867

  • CrossRef
  • Google Scholar

• 33

  ShenZ.SongJ.YungB. C.ZhouZ.WuA.ChenX. (2018). Emerging Strategies of Cancer Therapy Based on Ferroptosis. Adv. Mater.30, 1704007. 10.1002/adma.201704007

  • CrossRef
  • Google Scholar

• 34

  ShiG.ZhangJ.LuZ.LiuD.WuY.WuP.et al (2017). A Novel Messenger RNA Signature as a Prognostic Biomarker for Predicting Relapse in Pancreatic Ductal Adenocarcinoma. Oncotarget8, 110849–110860. 10.18632/oncotarget.22861

  • CrossRef
  • Google Scholar

• 35

  ShinD.KimE. H.LeeJ.RohJ.-L. (2018). Nrf2 Inhibition Reverses Resistance to GPX4 Inhibitor-Induced Ferroptosis in Head and Neck Cancer. Free Radic. Biol. Med.129, 454–462. 10.1016/j.freeradbiomed.2018.10.426

  • CrossRef
  • Google Scholar

• 36

  ShintokuR.TakigawaY.YamadaK.KubotaC.YoshimotoY.TakeuchiT.et al (2017). Lipoxygenase‐mediated Generation of Lipid Peroxides Enhances Ferroptosis Induced by Erastin and RSL3. Cancer Sci.108 (11), 2187–2194. 10.1111/cas.13380

  • CrossRef
  • Google Scholar

• 37

  SiegelR. L.MillerK. D.JemalA. (2019). Cancer Statistics, 2019. CA A. Cancer J. Clin.69, 7–34. 10.3322/caac.21551

  • CrossRef
  • Google Scholar

• 38

  SkoupilováH.MichalováE.HrstkaR. (2018). Ferroptosis as a New Type of Cell Death and its Role in Cancer Treatment. Klin Onkol31, 21–26. 10.14735/amko20182S21

  • CrossRef
  • Google Scholar

• 39

  StockwellB. R.Friedmann AngeliJ. P.BayirH.BushA. I.ConradM.DixonS. J.et al (2017). Ferroptosis: A Regulated Cell Death Nexus Linking Metabolism, Redox Biology, and Disease. Cell171, 273–285. 10.1016/j.cell.2017.09.021

  • CrossRef
  • Google Scholar

• 40

  SugiuraT.UesakaK.KanemotoH.MizunoT.SasakiK.FurukawaH.et al (2012). Serum CA19-9 Is a Significant Predictor Among Preoperative Parameters for Early Recurrence after Resection of Pancreatic Adenocarcinoma. J. Gastrointest. Surg.16, 977–985. 10.1007/s11605-012-1859-9

  • CrossRef
  • Google Scholar

• 41

  TanZ.LeiY.XuJ.ShiS.HuaJ.ZhangB.et al (2020). The Value of a Metabolic Reprogramming-Related Gene Signature for Pancreatic Adenocarcinoma Prognosis Prediction. Aging12, 24228–24241. 10.18632/aging.104134

  • CrossRef
  • Google Scholar

• 42

  TangR.HuaJ.XuJ.LiangC.MengQ.LiuJ.et al (2020). The Role of Ferroptosis Regulators in the Prognosis, Immune Activity and Gemcitabine Resistance of Pancreatic Cancer. Ann. Transl Med.8, 1347. 10.21037/atm-20-2554a

  • CrossRef
  • Google Scholar

• 43

  VincentA.HermanJ.SchulickR.HrubanR. H.GogginsM. (2011). Pancreatic Cancer. The Lancet378, 607–620. 10.1016/s0140-6736(10)62307-0

  • CrossRef
  • Google Scholar

• 44

  WuJ.BaoL.ZhangZ.YiX. (2017). Nrf2 Induces Cisplatin Resistance via Suppressing the Iron export Related Gene SLC40A1 in Ovarian Cancer Cells. Oncotarget8, 93502–93515. 10.18632/oncotarget.19548

  • CrossRef
  • Google Scholar

• 45

  WuM.LiX.LiuR.YuanH.LiuW.LiuZ. (2020). Development and Validation of a Metastasis-Related Gene Signature for Predicting the Overall Survival in Patients with Pancreatic Ductal Adenocarcinoma. J. Cancer11, 6299–6318. 10.7150/jca.47629

  • CrossRef
  • Google Scholar

• 46

  WuJ.ZhangL.WuS.YiX.LiuZ. (2020). miR-194-5p Inhibits SLC40A1 Expression to Induce Cisplatin Resistance in Ovarian Cancer. Pathol. - Res. Pract.216, 152979. 10.1016/j.prp.2020.152979

  • CrossRef
  • Google Scholar

• 47

  XieY.HouW.SongX.YuY.HuangJ.SunX.et al (2016). Ferroptosis: Process and Function. Cell Death Differ.23 (3), 369–379. 10.1038/cdd.2015.158

  • CrossRef
  • Google Scholar

• 48

  XiongG.HuangH.FengM.YangG.ZhengS.YouL.et al (2018). MiR-10a-5p Targets TFAP2C to Promote Gemcitabine Resistance in Pancreatic Ductal Adenocarcinoma. J. Exp. Clin. Cancer Res.37, 76. 10.1186/s13046-018-0739-x

  • CrossRef
  • Google Scholar

• 49

  XuD.WangY.LiuX.ZhouK.WuJ.ChenJ.et al (2021). Development and Clinical Validation of a Novel 9-gene Prognostic Model Based on Multi-Omics in Pancreatic Adenocarcinoma. Pharmacol. Res.164, 105370. 10.1016/j.phrs.2020.105370

  • CrossRef
  • Google Scholar

• 50

  YeZ.HuQ.ZhuoQ.ZhuY.FanG.LiuM.et al (2020). Abrogation of ARF6 Promotes RSL3-Induced Ferroptosis and Mitigates Gemcitabine Resistance in Pancreatic Cancer Cells. Am. J. Cancer Res.10, 1182–1193.

  • Google Scholar

• 51

  ZhouN.BaoJ. (2020, FerrDb: a Manually Curated Resource for Regulators and Markers of Ferroptosis and Ferroptosis-Disease Associations, Database (Oxford). 2020. baaa02110.1093/database/baaa021

  • CrossRef
  • Google Scholar

• 52

  ZhouM.YeZ.GuY.TianB.WuB.LiJ. (2015). Genomic Analysis of Drug Resistant Pancreatic Cancer Cell Line by Combining Long Non-coding RNA and mRNA Expression Profling. Int. J. Clin. Exp. Pathol.8, 38–52.

  • Google Scholar

• 53

  ZhouY.ZhouB.PacheL.ChangM.KhodabakhshiA. H.TanaseichukO.et al (2019). Metascape Provides a Biologist-Oriented Resource for the Analysis of Systems-Level Datasets. Nat. Commun.10, 1523. 10.1038/s41467-019-09234-6

  • CrossRef
  • Google Scholar

• 54

  ZhuS.ZhangQ.SunX.ZehH. J.3rdLotzeM. T.KangR.et al (2017). HSPA5 Regulates Ferroptotic Cell Death in Cancer Cells. Cancer Res.77 (8), 2064–2077. 10.1158/0008-5472.Can-16-1979

  • CrossRef
  • Google Scholar