- Open Access
Combination of ferroptosis and pyroptosis to construct a prognostic classifier and predict immune landscape, chemotherapeutic efficacy and immunosuppressive molecules in hepatocellular carcinoma
BMC Cancer volume 22, Article number: 229 (2022)
The induction of ferroptosis and pyroptosis has been highlighted as a novel approach to decide cancer cell fate. However, few studies have systematically explored the role of combining these two novel cell death modalities in hepatocellular carcinoma (HCC).
Ferroptosis-related genes (FRGs) and pyroptosis-related genes (PRGs) were retrieved and downloaded from FerrDb and GeneCards database, respectively. A prognostic classifier integrating with prognostic differentially expressed FRGs and PRGs was constructed by the least absolute shrinkage and selection operator (LASSO) algorithm in the TCGA-LIHC dataset and verified using the ICGC (LIRI-JP) dataset.
A total of 194 differentially expressed FRGs and PRGs were identified and named as differentially expressed genes (DEGs) and, out of them 79 were found dramatically correlated with prognosis in HCC. Based on 13 key DEGs with prognostic value, a novel expression signature was constructed and used to stratify HCC patients into 2 groups. Kaplan–Meier analysis demonstrated that high-risk patients had a more dismal prognosis. Receiver operating characteristic curve (ROC) and multivariate Cox analysis confirmed its predictive power and independent characteristic. Immune profile analysis demonstrated that high-risk group had prominent upregulation of immunosuppressive cells, including macrophages, Th2_cells and Treg. The correlation analysis between this signature and immunosuppressive molecules, Immunophenoscore (IPS) and chemotherapeutic efficacy demonstrated that low-risk group had a higher IC50 of cisplatin, mitomycin and doxorubicin and negatively related with CTLA4, HAVCR2, LAG3, PDCD1, TIGIT and ICIs treatment represented by CTLA4-/PD-1-, CTLA4 + /PD-1-, CTLA4-/PD-1 + .
In this research, a novel expression signature was identified based on FRGs and PRGs in HCC, and this signature could be used to predict prognosis and select patients potentially benefiting from immunotherapies and chemotherapy.
Hepatocellular carcinoma (HCC) is a prevalent malignancy worldwide which is characterized by increasing incidence and unfavorable prognosis [1, 2]. Although early-stage HCC patients could receive liver resection, transplantation and radiofrequency ablation, many patients still suffer from tumor recurrence . As a novel therapeutic approach, immunotherapies based on immune checkpoint inhibitors (ICIs) have benefited HCC patients in many clinical trials . For unresectable HCC patients, the therapeutic efficacy of atezolizumab plus bevacizumab for overall survival (OS) is superior to that of sorafenib . However, some HCC patients who receive ICIs treatment, such as nivolumab and pembrolizumab fail to show significant improvement in OS [6, 7], which might be due to tumors’ innate resistance to apoptosis . Thus, inducing novel modalities of cell death has become a promising target of antitumor therapeutic strategy. Ferroptosis and pyroptosis are such essential biological processes in HCC [9,10,11].
As an iron-dependent type of regulated cell death, ferroptosis is characterized by accumulation of lipid peroxidation to lethal levels . Currently, genes identified to regulate this novel form of cell death could be classified into 3 categories: drivers of ferroptosis (DOF), suppressors of ferroptosis (SOF) and others, which could either drive or suppress ferroptosis based on the context [13, 14]. Pyroptosis is a lytic form of regulated cell death characterized by release of many proinflammatory mediators. There are 2 major methods by which dead cells could activate pyroptosis: GSDMD-dependent manner regulated by caspase1/4/5/11 and GSDME-dependent manner regulated by caspase 3 [15,16,17,18,19].
Accumulating evidence has identified the induction of ferroptosis and pyroptosis as a novel approach by which CD8 + T cells could inhibit tumor growth. For instance, CD8 + T cells could release IFN-γ to downregulate SLC7A11 expression, resulting in lipid ROS accumulation and tumor cell ferroptosis . The activation of ferroptosis further promotes antitumor immunity. Besides, CD8 + T cells could release GzmA (GSDMB-cleaving enzyme) and GzmB (GSDME-cleaving enzyme) to induce pyroptosis. Induced tumor cell pyroptosis could activate IL-1β, which is derived from macrophages and required for antitumor immunity . The induction of ferroptosis and pyroptosis could enhance anticancer immunity and suppress tumor growth, suggesting a favorable prognosis for HCC patients. However, few studies have systematically discussed the possibility of combining these 2 cell death modalities in HCC.
Thus, our study focuses on the comprehensive analysis of a combined ferroptosis-related genes (FRGs) and pyroptosis related genes (PRGs) for HCC with regard to prognosis, clinicopathological feature, chemotherapeutic efficacy, tumor-infiltrating immune cells and immunosuppressive molecules.
Materials and methods
Acquisition of data, FRGs and PRGs
Gene expression profiling and survival data of 365 HCC patients were obtained from The Cancer Genome Atlas liver hepatocellular carcinoma (TCGA-LIHC) dataset [21, 22]. The scale method provided by R “limma” package was used to normalize gene expression values. Another 231 HCC patients with valid RNA-seq data and survival data from the ICGC (LIRI-JP) dataset were downloaded (Table 1). Gene expression values after read count normalization were used.
Then, 173 FRGs and 120 PRGs were retrieved from the FerrDb and GeneCards website, respectively.
Generation of differentially expressed genes (DEGs) with prognostic value
DEGs between HCC samples and normal ones were identified by R “limma” package in the TCGA dataset and false discovery rate (FDR) < 0.05 was set as the threshold. Then, DEGs were subjected to univariate Cox analysis to screen out FRGs and PRGs with prognostic value and P value < 0.05 was regarded as statistical difference. Venn diagram was plotted in which the interaction between DEGs and prognosis-related genes was displayed. Correlation analysis among these prognosis-related DEGs was conducted and an interaction network was analyzed in the STRING database to identify hub genes .
Construction of a combined ferroptosis and pyroptosis signature and assessment of its clinical utility
The least absolute shrinkage and selection operator (LASSO) algorithm with tenfold cross-validation, which could minimize the risk of overfitting was used to shrink and select variables . Some genes with a regression coefficient of non-zero were identified as the optimal predictors for OS and incorporated into this novel signature, whose risk score was calculated according to the normalized gene expression value and its corresponding regression coefficient. The median risk score was then used as the cutoff value to divide HCC patients into 2 groups. The difference of clinicopathological features between high- and low-risk group was analyzed by Wilcoxon signed-rank test. The relationship between clinicopathological parameters and risk score was investigated by Chi-square test.
Estimation of chemotherapeutic efficacy, ICIs-related molecules and Immunophenoscore (IPS) with this signature
To assess whether this signature was associated with the half inhibitory concentration IC50 of common antitumor drugs and chemotherapeutic efficacy, we applied “pRRophetic” package in R. By constructing the ridge regression model based on Genomics of Drug Sensitivity in Cancer (GDSC) (www.cancerrxgene.org/) cell line expression spectrum and TCGA gene expression profiles, the package could apply pRRophetic algorithm to predict drug IC50. Wilcoxon signed-rank test was implemented to compare the difference of IC50 between different risk groups. To investigate the relationship between this combined ferroptosis and pyroptosis signature and immunosuppressive molecules, we explored the difference of CTLA4, HAVCR2, LAG3, PDCD1 and TIGIT expression between high- and low-risk group using R “limma” package and applied “ggpubr” package to transform the results into a visual violin plot. As a superior biomarker to predict response of anti-PD-1 and CTLA-4 therapies, IPS could calculate the determinants of tumor immunogenicity and depict the cancer antigenomes and intra-tumoral immune profiles. This scoring scheme derived from a panel of immune-related genes, which belong to four classes: suppressor cells, effector cells, immunomodulators or checkpoints, and MHC-related molecules. By averaging the samplewise Z scores of the four classes within the respective category, the sum of the weighted averaged Z score was calculated as the IPS.
Validation of this combined ferroptosis and pyroptosis related signature
The Kaplan–Meier analysis was conducted to analyze the difference of OS between risk groups. R software was used to visualize the distribution of risk score and survival outcome of each HCC patient. R “timeROC” package was used to calculate area under the curve (AUC) of 1-, 2-, 3-year receiver operating characteristic curve (ROC) to evaluate the predictive ability of this novel signature. Principal component analysis (PCA) and t-SNE analysis were conducted to explore whether this signature could differentiate HCC patients between different risk groups. Uni- and multi variate Cox analyses were conducted to confirm whether this signature could serve as an independent predictor for HCC prognosis.
Functional enrichment analysis and immune profile analysis
DEGs between different risk groups were screened out by R “limma” package and we set FDR < 0.05 and |log2 fold change > 1| as the threshold. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis [25,26,27] were then performed to understand the biological function and pathways of these DEGs by using “clusterProfiler” R package. To explore the immune infiltration profiles between different risk groups, we conducted single-sample gene set enrichment analysis (ssGSEA) to calculate the score of 16 immune-cell features and 13 immune-function characteristics .
Identification of prognostic DEGs
A total of 194 DEGs between 374 HCC samples and 50 normal ones were screened out and out of them 79 were identified associated with OS in the univariate Cox analysis (Fig. 1). The protein–protein interaction (PPI) network among these prognostic DEGs were presented in Fig. 2a, in which there were 79 nodes and 230 edges. Genes with the top 15 degree of interaction were identified as hub genes (Fig. 2b). The correlation among 79 prognostic DEGs was displayed in Fig. 2c.
Construction of a combined ferroptosis and pyroptosis signature
Based on the expression profiles of 79 prognostic DEGs mentioned above, we conducted LASSO regression analysis to develop a novel signature, in which a total of 13 genes were identified as the optimal variables. Among them ATG3, FTL3, G6PD, HILPDA, NRAS, PRDX6, SLC1A5, SLC7A11 were FRGs, SQSTM1 participated in both ferroptosis and pyroptosis, and the remaining 4 genes (GLMN, LRPPRC, MKI67, UBE2D2) were PRGs. The risk score for this novel signature was: [ATG3 expression * (0.0599818988151381)] + [FLT3 expression * (-0.321132320389413)] + [G6PD expression * (0.0881814324303116)] + [GLMN expression * (0.130781902193193)] + [HILPDA expression * (0.119282064768739)] + [LRPPRC expression * (0.00792886569542188)] + [MKI67 expression * (0.0165502840606549)] + [NRAS expression * (0.0916391974243284)] + [PRDX6 expression * (0.114398632925529)] + [SLC1A5 expression * (0.0521211560497305)] + [SLC7A11 expression * (0.0616722848553423)] + [SQSTM1 expression * (0.00940765304518399)] + [UBE2D2 expression * (0.04686256927426)]. The median risk score was then used as the cutoff value to stratify HCC samples into 2 groups, in which there were 182 high- and 183 low-risk cases in the TCGA dataset.
Evaluation of the clinical utility of this novel signature
Kaplan–Meier analysis indicated that high-risk group had a more dismal OS (Fig. 3a, P < 0.001). Consistently, as was shown in Fig. 3b, low-risk patients had a lower probability to suffer from earlier death compared with those high-risk counterparts. PCA and t-SNE analysis demonstrated that HCC samples in different risk groups were easily distinguished (Fig. 3c and d). Uni- and multi- variate Cox analyses confirmed that this signature could predict prognosis independent of clinicopathological parameters (Fig. 3e and f). The AUC of 1-, 2-, 3-year ROC for this signature was 0.811, 0.743 and 0.721 (Fig. 3g) and the AUC of this signature was higher than that of clinicopathological indicators (Fig. 3h). Besides, Wilcoxon signed rank test and Chi-square test showed that clinicopathological features (Fig. 4a), including tumor grade (Fig. 4b), clinical stage (Fig. 4c) and T stage (Fig. 4d) were different between risk groups.
Exploration of the relationship between this signature and chemotherapeutic efficacy, immunosuppressive molecules and IPS
By evaluating the role of this signature in predicting the efficacy of common chemotherapeutics, immunosuppressive molecules and IPS, we discovered that low-risk group had a higher IC50 of cisplatin, mitomycin and doxorubicin (Fig. 5a) and was negatively related with CTLA4, HAVCR2, LAG3, PDCD1, TIGIT (Fig. 5b) and ICIs treatment represented by CTLA4-/PD-1-, CTLA4 + /PD-1-, CTLA4-/PD-1 + (Fig. 5c).
Verification of this combined ferroptosis and pyroptosis signature and construction of a nomogram
To assess the predictive ability of this novel signature, ICGC dataset was used as external validation. By calculating the risk score of each HCC sample based on the same formula derived from TCGA dataset, we stratified them into high- or low- risk group according to the cutoff median value (Fig. 6a). Patients with low risk were less susceptible to earlier death and had favorable OS by comparison with those high-risk counterparts (Fig. 6b and c). The AUC of 1-, 2-, 3-year ROC for this signature in the ICGC dataset was 0.750, 0.728 and 0.715 (Fig. 6d). PCA and t-SNE analysis demonstrated that high- and low-risk HCC samples were scattered in two directions (Fig. 6e and f). Uni- and multi-variate Cox analysis confirmed that this signature could predict OS independent of clinicopathological parameters (Fig. 6g and h). To provide a quantitative approach by which clinicians could predict probability of survival in HCC patients, we developed a nomogram by integrating of this novel prognostic signature and clinicopathological characteristics, including age, albumin, total bilirubin, prothrombin time, Child_Pugh classification, histologic grade, TNM stage and liver cirrhosis (Fig. 7a). The predictive accuracy of this nomogram was validated by its discrimination and calibration performance (Fig. 7b and c).
Functional analysis and immune cell infiltration
To explore the biological function and pathways associated with this signature, we performed GO and KEGG analysis based on DEGs between high- and low-risk group. In both TCGA and ICGC dataset, GO terms for biological process were organelle fission, nuclear division and sister chromatid segregation; for cellular component were mainly associated with chromosome, such as spindle, condensed chromosome, centromeric region (Fig. 8a and b). Besides, the KEGG results revealed that these DEGs mainly participated in cell cycle, PI3K-Akt signaling pathway and cellular senescence (Fig. 8c and d).
To further investigate the relationship between this signature and immune profiles, we performed ssGSEA to calculate the score of immune-cell features and immune-function characteristics. The score of macrophages, Th2_cells and Treg was higher in high-risk group, while neutrophils and NK cells earned a higher score in low-risk group (Fig. 9a and b). Besides, MHC_class_I and Type_II_IFN_Response was mainly enriched in high- and low-risk group, respectively (Fig. 9c and d).
As two novel cell death modalities, mounting evidence has demonstrated that ferroptosis and pyroptosis are intimately related with tumor progression [29,30,31]. However, few studies have thoroughly explored the role of combing these 2 novel cell death modalities in HCC. In the present research, we systematically investigated the expression and prognostic significance of 173 FRGs and 120 PRGs in HCC. To our surprise, most of the FRGs (75.2%) and PRGs (59.2%) were found differentially expressed between HCC samples and normal ones, and approximately 40% of them were associated with survival in the univariate Cox analysis. These findings significantly demonstrated the potential role of ferroptosis and pyroptosis in HCC and the possibility of constructing a prognostic expression signature from the perspective of these 2 cell death modalities.
A total of 13 genes were finally incorporated into this novel signature. ATG3, a major autophagy regulator, could promote ferroptosis via degradation of ferritin and accumulation of iron . FLT3 inhibitor could prevent lipid peroxidation to protect cells against oxidative glutamate toxicity, which involves 3 cell death modalities: apoptosis, ferroptosis and necroptosis . G6PD has been considered as a critical pacesetter of pentose phosphate pathway that prevents erastin-induced ferroptosis when it was knocked down . HIF-2a-HILPDA could selectively enrich polyunsaturated fatty acids and acts as a central driver of a ferroptotic cell death . In HL-60/NRASQ61L cells, HMGB1 knockdown could decrease erastin-induced ferroptosis via RAS-JNK/p38 signaling . PRDX6, an essential family member of nonselenium peroxidases, could protect cells against ferroptotic process and serve as a potential target to enhance antitumor activity of ferroptosis-based chemotherapeutics . SLC1A5, which mediates uptake of neutral amino acids, has been reported to cause ferroptotic cell death when it is suppressed by miR-137 . SLC7A11 is a mechanistic determinant which promotes tumoral lipid oxidation and ferroptosis when it is synergistically suppressed by immunotherapy and radiotherapy . SQSTM1 is passively released in the context of GSDMD-mediated pyroptosis and found to mediate various biological processes, especially autophagy and ferroptosis [38, 39]. Downregulation of GLMN levels could activate inflammasome and pyroptotic cell death of macrophages . Tetherin could prevent T cell pyroptosis by interacting with LRPPRC and preventing the formation of LRPPRC-Bcl-2-Beclin 1 ternary complex . MKI67, a well-known marker of proliferation, was found involved in caspase-3-mediated apoptosis and caspase-1-mediated pyroptosis of CD4 T cells in HIV . Knockdown of UBE2D2, a member of ubiquitin-conjugating enzymes, could significantly reduce SQSTM1 recruitment, which is passively released in the context of GSDMD-mediated pyroptosis .
Recent studies have confirmed that ferroptosis and pyroptosis are significantly associated with antitumor activity. Tumor cells undergoing ferroptosis and pyroptosis could recruit tumor-suppressed immune cells and enhance antitumor immunity. To investigate the immune-cell characteristics and immune-function features of this signature, we performed ssGSEA based on DEGs between different risk groups. The results showed that high-risk group of both datasets earned higher scores of macrophages, Th2_cells and Treg, which are thought to promote tumor growth and invasion and tightly associated with dismal prognosis [44,45,46,47]. MHC_class_I, a major antigen presenting process, was also enriched in high-risk group. It is possibly because the pore formation with the pyroptotic cellular plasma membrane could release immune stimulants, which could attract activated dendritic cells and thus, promote anti-tumor T cell activity . Besides, cells undergoing ferroptosis could attract dendritic cells by releasing lipid mediators . The scores of NK_cells and Type_II_IFN_Response were higher in low-risk group by the fact that type II IFN is mainly released by NK cells , which are major components of immune defense against tumorigenesis .
Several shortcomings should be acknowledged in our study. First, this combined ferroptosis and pyroptosis signature was constructed and verified using public dataset. Real-world dataset is required to assess its accuracy and efficacy. Second, the molecular mechanisms between genes identified by this signature and tumor immunity in HCC are not elucidated. Future experimental studies are needed to address this problem.
In conclusion, a combined ferroptosis and pyroptosis signature was constructed for HCC, which was tightly associated with prognosis, clinicopathological features, immune profiles, chemotherapeutic efficacy and immunosuppressive molecules. The molecular mechanisms between this signature and tumor immunity are largely unknown and require further experimental investigation.
Immune checkpoint inhibitors
Drivers of ferroptosis
Suppressors of ferroptosis
Differentially expressed genes
False discovery rate
The least absolute shrinkage and selection operator
- IC50 :
The half inhibitory centration
Area under the curve
Receiver operating characteristic curve
Principal component analysis
Kyoto Encyclopedia of Genes and Genomes
Single-sample gene set enrichment analysis
Hartke J, Johnson M, Ghabril M. The diagnosis and treatment of hepatocellular carcinoma [J]. Semin Diagn Pathol. 2017;34(2):153–9.
Kulik L, El-Serag HB. Epidemiology and Management of Hepatocellular Carcinoma [J]. Gastroenterology. 2019;156(2):477-91.e1.
Tabrizian P, Jibara G, Shrager B, et al. Recurrence of hepatocellular cancer after resection: patterns, treatments, and prognosis [J]. Ann Surg. 2015;261(5):947–55.
Hilmi M, Neuzillet C, Calderaro J, et al. Angiogenesis and immune checkpoint inhibitors as therapies for hepatocellular carcinoma: current knowledge and future research directions [J]. J Immunother Cancer. 2019;7(1):333.
Finn RS, Qin S, Ikeda M, et al. Atezolizumab plus Bevacizumab in Unresectable Hepatocellular Carcinoma [J]. N Engl J Med. 2020;382(20):1894–905.
Finn RS, Ryoo BY, Merle P, et al. Pembrolizumab As Second-Line Therapy in Patients With Advanced Hepatocellular Carcinoma in KEYNOTE-240: A Randomized, Double-Blind, Phase III Trial [J]. J Clin Oncol. 2020;38(3):193–202.
Liu Z, Lin Y, Zhang J, et al. Molecular targeted and immune checkpoint therapy for advanced hepatocellular carcinoma [J]. J Exp Clin Cancer Res. 2019;38(1):447.
Tang R, Xu J, Zhang B, et al. Ferroptosis, necroptosis, and pyroptosis in anticancer immunity [J]. J Hematol Oncol. 2020;13(1):110.
Zhang Z, Zhang Y, Xia S, et al. Gasdermin E suppresses tumour growth by activating anti-tumour immunity [J]. Nature. 2020;579(7799):415–20.
Lang X, Green MD, Wang W, et al. Radiotherapy and Immunotherapy Promote Tumoral Lipid Oxidation and Ferroptosis via Synergistic Repression of SLC7A11 [J]. Cancer Discov. 2019;9(12):1673–85.
Wang W, Green M, Choi JE, et al. CD8(+) T cells regulate tumour ferroptosis during cancer immunotherapy [J]. Nature. 2019;569(7755):270–4.
Dixon SJ, Lemberg KM, Lamprecht MR, et al. Ferroptosis: an iron-dependent form of nonapoptotic cell death [J]. Cell. 2012;149(5):1060–72.
Jiang L, Kon N, Li T, et al. Ferroptosis as a p53-mediated activity during tumour suppression [J]. Nature. 2015;520(7545):57–62.
Tarangelo A, Magtanong L, Bieging-Rolett KT, et al. p53 Suppresses Metabolic Stress-Induced Ferroptosis in Cancer Cells [J]. Cell Rep. 2018;22(3):569–75.
Agard NJ, Maltby D, Wells JA. Inflammatory stimuli regulate caspase substrate profiles [J]. Mol Cell Proteomics. 2010;9(5):880–93.
Kayagaki N, Stowe IB, Lee BL, et al. Caspase-11 cleaves gasdermin D for non-canonical inflammasome signalling [J]. Nature. 2015;526(7575):666–71.
Shi J, Zhao Y, Wang K, et al. Cleavage of GSDMD by inflammatory caspases determines pyroptotic cell death [J]. Nature. 2015;526(7575):660–5.
Crawford ED, Wells JA. Caspase substrates and cellular remodeling [J]. Annu Rev Biochem. 2011;80:1055–87.
Julien O, Wells JA. Caspases and their substrates [J]. Cell Death Differ. 2017;24(8):1380–9.
Stockwell BR, Jiang X. A Physiological Function for Ferroptosis in Tumor Suppression by the Immune System [J]. Cell Metab. 2019;30(1):14–5.
Gu X, Guan J, Xu J, et al. Model based on five tumour immune microenvironment-related genes for predicting hepatocellular carcinoma immunotherapy outcomes [J]. J Transl Med. 2021;19(1):26.
Liu GM, Zeng HD, Zhang CY, et al. Identification of a six-gene signature predicting overall survival for hepatocellular carcinoma [J]. Cancer Cell Int. 2019;19:138.
Szklarczyk D, Franceschini A, Kuhn M, et al. The STRING database in 2011: functional interaction networks of proteins, globally integrated and scored [J]. Nucleic Acids Res. 2011;39(Database issue):D561-8.
Tibshirani R. The lasso method for variable selection in the Cox model [J]. Stat Med. 1997;16(4):385–95.
Kanehisa M. Toward understanding the origin and evolution of cellular organisms [J]. Protein Sci. 2019;28(11):1947–51.
Kanehisa M, Furumichi M, Sato Y, et al. KEGG: integrating viruses and cellular organisms [J]. Nucleic Acids Res. 2021;49(D1):D545–51.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes [J]. Nucleic Acids Res. 2000;28(1):27–30.
Rooney MS, Shukla SA, Wu CJ, et al. Molecular and genetic properties of tumors associated with local immune cytolytic activity [J]. Cell. 2015;160(1–2):48–61.
Cui J, Zhou Z, Yang H, et al. MST1 Suppresses Pancreatic Cancer Progression via ROS-Induced Pyroptosis [J]. Mol Cancer Res. 2019;17(6):1316–25.
Liu HJ, Hu HM, Li GZ, et al. Ferroptosis-Related Gene Signature Predicts Glioma Cell Death and Glioma Patient Progression [J]. Front Cell Dev Biol. 2020;8:538.
Wu Y, Zhang S, Gong X, et al. The epigenetic regulators and metabolic changes in ferroptosis-associated cancer progression [J]. Mol Cancer. 2020;19(1):39.
Zhou B, Liu J, Kang R, et al. Ferroptosis is a type of autophagy-dependent cell death [J]. Semin Cancer Biol. 2020;66:89–100.
Kang Y, Tiziani S, Park G, et al. Cellular protection using Flt3 and PI3Kα inhibitors demonstrates multiple mechanisms of oxidative glutamate toxicity [J]. Nat Commun. 2014;5:3672.
Zou Y, Palte MJ, Deik AA, et al. A GPX4-dependent cancer cell state underlies the clear-cell morphology and confers sensitivity to ferroptosis [J]. Nat Commun. 2019;10(1):1617.
Ye F, Chai W, Xie M, et al. HMGB1 regulates erastin-induced ferroptosis via RAS-JNK/p38 signaling in HL-60/NRAS(Q61L) cells [J]. Am J Cancer Res. 2019;9(4):730–9.
Lu B, Chen XB, Hong YC, et al. Identification of PRDX6 as a regulator of ferroptosis [J]. Acta Pharmacol Sin. 2019;40(10):1334–42.
Luo M, Wu L, Zhang K, et al. miR-137 regulates ferroptosis by targeting glutamine transporter SLC1A5 in melanoma [J]. Cell Death Differ. 2018;25(8):1457–72.
Zou B, Liu J, Klionsky DJ, et al. Extracellular SQSTM1 as an inflammatory mediator [J]. Autophagy. 2020;16(12):2313–5.
Liu J, Yang M, Kang R, et al. Autophagic degradation of the circadian clock regulator promotes ferroptosis [J]. Autophagy. 2019;15(11):2033–5.
Suzuki S, Suzuki T, Mimuro H, et al. Shigella hijacks the glomulin-cIAPs-inflammasome axis to promote inflammation [J]. EMBO Rep. 2018;19(1):89–101.
Zou J, Li W, Misra A, et al. The viral restriction factor tetherin prevents leucine-rich pentatricopeptide repeat-containing protein (LRPPRC) from association with beclin 1 and B-cell CLL/lymphoma 2 (Bcl-2) and enhances autophagy and mitophagy [J]. J Biol Chem. 2015;290(11):7269–79.
Doitsh G, Galloway NL, Geng X, et al. Cell death by pyroptosis drives CD4 T-cell depletion in HIV-1 infection [J]. Nature. 2014;505(7484):509–14.
Geisler S, Vollmer S, Golombek S, et al. The ubiquitin-conjugating enzymes UBE2N, UBE2L3 and UBE2D2/3 are essential for Parkin-dependent mitophagy [J]. J Cell Sci. 2014;127(Pt 15):3280–93.
Zhou SL, Zhou ZJ, Hu ZQ, et al. Tumor-Associated Neutrophils Recruit Macrophages and T-Regulatory Cells to Promote Progression of Hepatocellular Carcinoma and Resistance to Sorafenib [J]. Gastroenterology. 2016;150(7):1646-58.e17.
Zhang Q, He Y, Luo N, et al. Landscape and Dynamics of Single Immune Cells in Hepatocellular Carcinoma [J]. Cell. 2019;179(4):829-45.e20.
Yao F, Yuan Q, Song X, et al. Yupingfeng Granule Improves Th2-Biased Immune State in Microenvironment of Hepatocellular Carcinoma through TSLP-DC-OX40L Pathway [J]. Evid Based Complement Alternat Med. 2020;2020:1263053.
Fu J, Xu D, Liu Z, et al. Increased regulatory T cells correlate with CD8 T-cell impairment and poor survival in hepatocellular carcinoma patients [J]. Gastroenterology. 2007;132(7):2328–39.
Erkes DA, Cai W, Sanchez IM, et al. Mutant BRAF and MEK Inhibitors Regulate the Tumor Immune Microenvironment via Pyroptosis [J]. Cancer Discov. 2020;10(2):254–69.
FriedmannAngeli JP, Krysko DV, Conrad M. Ferroptosis at the crossroads of cancer-acquired drug resistance and immune evasion [J]. Nat Rev Cancer. 2019;19(7):405–14.
di Franco S, Turdo A, Todaro M, et al. Role of Type I and II Interferons in Colorectal Cancer and Melanoma [J]. Front Immunol. 2017;8:878.
Liu P, Chen L, Zhang H. Natural Killer Cells in Liver Disease and Hepatocellular Carcinoma and the NK Cell-Based Immunotherapy [J]. J Immunol Res. 2018;2018:1206737.
All authors thank the family for participating in the study.
No sources of funding were required for this work.
Ethics approval and consent to participate
Investigation has been conducted in accordance with the Declaration of Helsinki and national and international guidelines. All datasets in the present research were downloaded from public databases, including TCGA and ICGC. Thus, the present study was exempted from the approval of ethics review board in Renji hospital, School of Medicine, Shanghai Jiao Tong University. The current study follows the TCGA and ICGC data access policies and publication guidelines.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Xu, L., Zheng, Q. & Liu, W. Combination of ferroptosis and pyroptosis to construct a prognostic classifier and predict immune landscape, chemotherapeutic efficacy and immunosuppressive molecules in hepatocellular carcinoma. BMC Cancer 22, 229 (2022). https://doi.org/10.1186/s12885-022-09301-0
- Hepatocellular carcinoma
- Overall survival
- Immune profiles
- Immunosuppressive molecules