- Open Access
A pyroptosis-related gene signature provides an alternative for predicting the prognosis of patients with hepatocellular carcinoma
BMC Medical Genomics volume 16, Article number: 2 (2023)
Hepatocellular Carcinoma (HCC) is a common malignant neoplasm with limited treatment options and poor outcomes. Thus, there is an urgent need to find sensitive biomarkers for HCC.
Gene expression and clinicopathological information were obtained from public databases, based on which a pyroptosis-related gene signature was constructed by the least absolute shrinkage and selection operator Cox regression. The applicability of the signature was evaluated via Kaplan–Meier curve and time-dependent ROC curve. TIMER, QUANTISEQ, MCPCOUNTER, EPIC, CIBERSORT, ssGSEA, and ESTIMATE were employed to assess the immune status. Comparisons between groups were analyzed with Wilcoxon test. Pearson and Spearman correlation analyses were adopted for linear correlation analysis. Genetic knockdown was conducted using siRNA transfection and the mRNA expression levels of interest genes were measured using quantitative reverse transcription PCR. Finally, protein levels in 10 paired tumor tissues and adjacent non-tumor tissues from HCC patients were measured using immunohistochemistry.
A pyroptosis-related gene signature was established successfully to calculate independent prognostic risk scores. It was found that survival outcomes varied significantly between different risk groups. In addition, an attenuated antitumor immune response was found in the high-risk group. Meanwhile, multiple immune checkpoints were up-regulated in high-risk score patients. Cell cycle-related genes, angiogenesis-related genes and tumor drug resistance genes were also markedly elevated. Knockdown of prognostic genes in the signature significantly inhibited the expression of immune checkpoint genes and angiogenesis-related genes. Besides, each prognostic gene was expressed at a higher level in HCC tissues than in adjacent normal tissues.
We successfully established a novel pyroptosis-related gene signature which could help predict the overall survival and assess the immune status of HCC patients.
Hepatocellular carcinoma (HCC) is one of the most frequently diagnosed malignant tumors, accounting for more than 90% of all liver cancer cases , and ranking fifth of the incidence and third in mortality of all malignancies worldwide . Risk factors for HCC include chronic hepatitis B and C virus infections, alcohol abuse, nonalcoholic fatty liver disease, and exposure to dietary toxins such as aflatoxins . Patients are often diagnosed with HCC at advanced stages with poor prognosis , and unclear molecular mechanisms lead to the poor understanding of HCC prognosis. The overall survival (OS) of HCC varies across the world, with a 5-year survival rate of 18% in the United States, and 12% in China . Therefore, it is primarily important to clarify the molecular mechanisms underlying the poor prognosis of HCC and explore new prognostic biomarkers for HCC.
Pyroptosis, an inflammatory form of cell death triggered by certain inflammasomes, is primarily mediated by caspase 1/4/5 (11)-cleaved GSDMD (gasdermin D). Afterwards, GSDMD-N oligomerizes form membrane pores, leading to the release of activated cytokines and pyroptotic cell death [5,6,7]. Some studies reported that pyroptosis played a double-edged role in cancer . While still an ambiguous process in cancer, pyroptosis could not only play a potent and persuasive role to conquer apoptosis resistance but also a crucial role in immunity . Conversely, pyroptosis can induce damage to healthy tissues and establish a microenvironment suitable for tumor growth and metastatic progression . Obviously, pyroptosis involves complex and contentious processes, so that its role and detailed mechanism in oncogenesis deserves extensive exploration. Although the exact link between pyroptosis and HCC is presently unclear, increasing research has begun to focus on HCC and pyroptosis [11, 12].
In this study, we attempted to explore the association between pyroptosis and HCC by establishing a predictive signature based on pyroptosis-related genes strongly related to poor outcomes of HCC, hoping that our findings could provide a new signature that could predict the clinical survival outcome and help design individualized treatment for HCC patients.
The data were obtained from public databases. HCC samples collected before 2011 were excluded because of the excessive length of storage. Finally, data of 209 HCC patients were obtained from The Cancer Genome Atlas hepatocellular carcinoma (TCGA-LIHC) portal (https://portal.gdc.cancer. gov/repository). Additional 231 tumor samples and 199 adjacent normal tissue samples were collected from the International Cancer Genome Consortium hepatocellular carcinoma (ICGC-LIRI-JP) portal (https://dcc.icgc.org/projects/LIRI-JP). Also, an independent cohort containing 10 paired HCC samples was recruited. Clinical data for these patients are presented in Table 1. The present study follows access policies and publication guidelines. Then, 69 pyroptosis-related genes were obtained from the previous literature [8, 13,14,15,16,17] (Additional file 1: Table S1).
Identification of prognostic genes
The differentially expressed genes (DEGs) between HCC and adjacent non-tumorous tissues were screened with an FDR < 0.05 and log2 (fold-change) > 0 from the whole genes for candidate genes by the "limma" R package. Univariate Cox analysis was implemented to identify significant prognostic genes from the DEGs (P < 0.05). Correlation networks were performed through “igraph” and “reshape2” packages. Protein–protein interaction (PPI) networks were constructed with the STRING database.
Gene signature establishment and validation
When constructing the prognostic gene model, the collected data were processed based on the least absolute shrinkage and selection operator (LASSO) penalized Cox regression analysis, and the penalty parameters were obtained by cross validation with the R package "glmnet" tool. When the dataset was large, K-Fold Cross Validation analysis was used, so that the optimal model and parameters can be selected for function evaluation. When the dataset was small (N < 50), the leave-one-out Cross Validation analysis was used because it uses more training samples in each iteration. The LASSO analysis was mainly used to select variables and determine the model that met the interpretation requirements according to the regression coefficient. The standardized expression matrix of candidate prognostic genes was the independent variable in the regression equation, and the dependent variables were OS and patient status. Based on the analysis of the expression value of each gene and the corresponding regression coefficient, the risk score of the patients was determined: score = esum (expression of each gene × Correspondence coefficient). All patients were divided into high- and low-risk groups based on the median score. The "Rtsne" R software was used for principal component analysis to determine the distribution of each group. The survival status of patients was analyzed based on the Kaplan–Meier test. The ROC curve was drawn by processing the relevant data, and then the AUC value was calculated by the R package “survivalROC” so as to determine the sensitivity of variables.
Identification of independent prognostic factors for OS in HCC
To evaluate the predictive power of the risk score, each prognostic gene in the signature (CASP3, IRAK1, MAPK1, MAPK3 and YWHAB) and the clinical risk factors, including age (< 60 vs. ≥ 60), gender (male vs. female), tumor grade (G1/2 vs. G3/4), tumor stage (I/II vs. III/IV), univariate and multivariate Cox regression analyses were performed to identify independent prognostic factors.
Immune status in distinct risk groups
Five algorithms (TIMER, QUANTISEQ, MCPCOUNTER, EPIC and CIBERSORT) were applied to evaluate the correlation between immune cells and risk scores by “immunedeconv” R package . The differences between these five algorithms are shown in Additional file 2: Table S2. Single-sample gene set enrichment analysis (ssGSEA), ESTIMATE and CIBERSORT were performed to assess the immune status in high- and low-risk groups by the "GSVA", “estimate” and “limma” R package based on all of the expressed genes. The expression levels of immune checkpoint genes in different risk groups were analyzed by Wilcoxon test.
Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses
To detect the potential biological functions and pathways in different risk groups, GO and KEGG pathway analyses were carried out using the “clusterProfiler” and “enrichplot” packages based on all of the expressed genes [19,20,21].
Cell cycle and angiogenesis-related genes expression in different risk groups
The expression of cell cycle and angiogenesis-related genes in different risk groups was compared using the Wilcoxon test.
Tumor drug resistance genes and chemotherapy sensitivity analysis
The NCI-60 database containing 60 different cancer cell lines from 9 different types of tumors was accessed through the CellMiner interface (https://discover.nci.nih.gov/cellminer/). Altogether 263 drugs on clinical trials or approved by the US Food and Drug Administration, were used to evaluate the correlation between prognostic gene expression and drug sensitivity, and that between the prognostic model and tumor drug resistance genes by Pearson and Spearman correlation analyses respectively.
Cell lines, cell culture and cell transfection
HCC cell line Huh7 was purchased from the Cell Bank of the Chinese Academy of Sciences (Shanghai, China). To reduce the risk of microbial contamination, cells were cultured in a suitable medium containing 10% fetal bovine serum (FBS; Gibco) and supplements of penicillin (100 U/mL) and streptomycin (100 g/mL) in a humidified environment with 5% CO2/95% air at 37 °C. Gene silencing was achieved using siRNA. siRNA was transfected into cells using Lipofectamine 3000 (Life Technologies), according to the manufacturer’s instructions. Transfected cells were analyzed 48 h after transfection with siRNA. siRNA sequences are listed in Additional file 3: Table S3.
Quantitative reverse transcription PCR (qRT-PCR)
Total RNA was extracted from the cells using TRIzol reagent (Thermo Scientific, Cat# 15596018). Then, the RNA was quantified using Nanodrop. The RNA was reverse-transcribed into DNA using a Prime Script RT Reagent Kit (Takara, Cat#RR036A) following the manufacturer’s instructions. Isolated DNA was subjected to qRT-PCR analysis using the CellAmp™ Direct TB Green® RT-qPCR Kit (Takara, Cat# 3735A). Primer sequences are listed in Additional file 4: Table S4.
Ten pairs of HCC and adjacent non-tumorous tissues collected by the First Affiliated Hospital of Wenzhou Medical University (Wenzhou, China) were selected as samples for IHC detection, so as to determine the protein expression of prognostic genes in these samples and evaluate the prognostic value of these genes. The research content was approved by the ethics review committee of the university, and all the subjects were informed and agreed. The samples were collected and then fixed in 10% formalin, placed at room temperature for half an hour, subsequently embedded in paraffin, processed into 4-μm slices, deparaffinized, rehydrated, and boiled for 10 min. Then immerse the slice in 3% hydrogen peroxide solution for 10 min, which could achieve the purpose of inactivation and avoid non-specific binding. Incubate the slice with 1% FBS in PBS for half an hour. After incubation, the sections were stained by the secondary antibody combined with the primary antibody and HRP. The antibodies are detailed in Additional file 5: Table S5. Then the sections were stained continuously with 3,3′-diaminobenzidine and hematoxylin. After dehydration, the samples were sealed, observed and photographed. Representative pictures of prognostic genes were displayed in our study. Finally, a quantitative analysis of IHC staining was conducted by ImageJ software.
Group comparison was made by the Wilcoxon test. Chi-square test was used to check the difference of results of categorical variables, and Kaplan Meier test was used for the difference of survival rates between the two groups. Cox regression analysis was conducted on the collected data to determine the factors closely related to OS, and the corresponding regression analysis model was established based on the results. The relationship between risk scores and drug resistance genes was determined by the Spearman method. The relationship between prognostic gene expression and drug sensitivity was judged by the Pearson coefficient. R4.0.2 and SPSS23.0 software were used for data statistical analysis, and the correlation curves were drawn. P < 0.05 indicates that the difference of results is statistically significant.
Screening of pyroptosis-related prognostic genes
The flowchart for the screening of candidate genes was exhibited in Fig. 1. Among the 69 pyroptosis-related genes, 45 of them were expressed differentially between tumor and adjacent non-tumorous tissues, and 19 of them were linked to OS (Fig. 2A). Ultimately, 17 genes were survival-related and discovered to have significant differential expression, based on which a pyroptosis-related gene signature was constructed. Figure 2B–C shows the prognostic genes and their hazard ratio. PPI networks and gene correlation networks revealed the interactions among these prognostic genes (Fig. 2D–E).
Construction of the prognostic signature
Five genes with maximum prognostic value (CASP3, IRAK1, MAPK1, MAPK3 and YWHAB) were found using LASSO Cox regression analysis (Fig. 2F–G). A pyroptosis-related prognostic model was constructed using the genes mentioned above in the TCGA cohort. The risk score = 0.97 × 10–2 × expression level of CASP3 + 0.65 × 10–2 × expression level of IRAK1 + 0.53 × 10–2 × expression level of MAPK1 + 0.83 × 10–2 × expression level of MAPK3 + 0.18 × 10–2 × expression level of YWHAB. According to the median cut-off value, the patients were grouped into either a high-risk group or a low-risk group (Fig. 3A). It was clearly demonstrated that patients in distinct risk groups were separated in discrete directions (Fig. 3B). In addition, patient survival in high-risk score group was significantly poorer than that in low-risk score group (Fig. 3E). Besides, time-dependent ROC analysis showed excellent performance of the prognostic signature (TCGA: 1-year AUC = 0.716, 2-year AUC = 0.665, 3-year OS = 0.707) (Fig. 3F).
Validation of the prognostic signature
We further examined the prognostic signature in the ICGC cohort for additional independent validation. Patients were separated into different risk groups using the same formula and median value from the TCGA cohort (Fig. 3C). The results in the ICGC cohort were virtually consistent with those in the TCGA cohort. Patients in different risk groups were separated in two directions and OS of patients in the high-risk group was significantly worse than that in the low-risk group. (Fig. 3D, G). The AUC of the pyroptosis-related signature for 1-, 2-, and 3-year OS were 0.755, 0.721 and 0.707, respectively (Fig. 3H).
Clinical characteristics of patients in different risk groups
The proportion of tumor stage I-II patients in the high-risk score group was significantly higher than that in the low-risk score group (P < 0.05), and vice versa for the proportion of tumor stage III-IV (Fig. 4D). Likewise, a similar trend was seen in the ICGC cohort (Fig. 4G), though the difference was not statistically significant (P = 0.055). In addition, no significant difference in age, gender and tumor grade was observed between the two groups (Fig. 4A–C, E–F). All these findings suggest that the tumor stage tended to be higher in patients with high-risk scores.
Independent prognostic analysis of the pyroptosis-related signature
The independent prognostic value of the risk score, each prognostic gene in the signature and clinical features was evaluated by univariate and multivariate Cox analyses. Tumor stages, YWHAB and risk scores in the TCGA cohort and gender, tumor stages, CASP3, IRAK1, MAPK1, MAPK3, YWHAB and risk scores in the ICGC cohort significantly correlated with the patient survival prognosis in univariate analysis were subjected to multivariate analysis, and the results showed that the risk score and the tumor stage were statistically significant (P < 0.05) (Fig. 5A, C). Clearly, the risk score was a prominent risk predictor after adjustment for known clinical and pathologic factors, but no single gene showed an independent prognostic value. Furthermore, the combination of the risk score with the tumor stage could provide a more accurate prediction of 1-, 2-, 3- year OS in HCC (TCGA: 1-year AUC = 0.746, 2-year AUC = 0.745, 3-year OS = 0.771; ICGC: 1-year AUC = 0.855, 2-year AUC = 0.747, 3-year OS = 0.756) (Fig. 5B, D).
Association between the signature and immune status
Next, five algorithms were performed to evaluate the correlation between immune cells and risk scores. It was found that the risk score was positively correlated with the infiltration of macrophages, myeloid dendritic cells, T cell regulatory (Treg), macrophage/monocytes and cancer-associated fibroblasts (Additional file 6: Fig. S1A, C). The expression of immune checkpoint genes including PDCD1 (PD-1), PDCD1LG2 (PD-L2), CTLA4, CD80, CD86, HAVCR2, LGALS9, CD276 and VTCN1 were also up‐regulated in the high‐risk group (Additional file 6: Fig. S1B, D). To further explore the immune status in different risk groups, we performed ssGSEA algorithm. In both the TCGA and the ICGC cohorts, the immune cell subpopulations including aDCs, iDCs, macrophages, Tfh as well as Treg showed high infiltration in the high-risk group (P < 0.05) (Additional file 6: Fig. S2A, B). Moreover, the scores of immune-related functions (CCR, Check point, MHC class I and parainflammation) in the high-risk group were significantly higher than those in the low-risk group (Additional file 6: Fig. S2C, D). Besides, CIBERSORT and ESTIMATE were also performed to distinguish immune status in high- and low-risk groups (Additional file 6: Fig. S2E–H).
Function enrichment analysis
GO terms, such as cell cycle checkpoint, cell cycle G2/M phase transition, cell cycle G1/S phase transition, cell cycle DNA replication, signal transduction involved in cell cycle checkpoint and cellular response to hypoxia, were significantly enriched by GO enrichment analysis (Additional file 6: Fig. S3A, C). In addition, KEGG pathway terms including cell cycle, hepatitis B, PD-L1 expression and PD-1 checkpoint pathway in cancer, T cell receptor signaling pathway, VEGF signaling pathway, HIF-1 signaling pathway, EGFR tyrosine kinase inhibitor resistance and hepatocellular carcinoma were significantly enriched by KEGG analysis (Additional file 6: Fig. S3B, D). Intriguingly, the KEGG HCC signaling pathway was significantly enriched, whether in the TCGA database or the ICGC database.
Expression levels of genes related to cell cycle and tumor angiogenesis
Wilcoxon tests showed that the expressions of cell cycle (CCNA2, CCNB1, CCNB2, CCND2, CCND3, CDC20, CDC23, CDC25A, CDC25B, CDC25C, CDK1, CDK2, CDK4, CDK7, CHEK1, CHEK2, E2F1, E2F3, E2F4 and GSK3B) and angiogenesis-related genes (HIF1A, FDGFRA, FDGFRB, FDGFA, FDGFB, NRP1, NRP2, ANGPT2, VEGFB, FGFR1, FGFR2, FGFR3, RCBO1, RCBC3, SLIT1 and SLIT2) were increased in the high-risk group than those in the low-risk group (Additional file 6: Fig. S4).
Tumor drug resistance analysis
So far, tumor drug resistance remains a main problem in tumor therapy, which urged us to further explore the correlation between the risk score and tumor drug resistance. It was found in our study that the expression levels of MRP1, MRP4 and MRP5 in the high-risk group were higher than those in the low-risk group. Furthermore, the expression levels of MRP1, MRP4 and MRP5 were positively correlated with the risk score (Additional file 6: Fig. S5A–B). Besides, Pearson correlation analysis revealed an inverse correlation (P < 0.05) between the expression levels of the prognostic genes and chemosensitivity (Additional file 6: Fig. S5C).
Knockdown of prognostic genes in HCC cell lines
To further substantiate the functions of the prognostic genes in the signature, siRNA knock-down experiments were performed. siRNA was transfected into HCC cells to inhibit the expression of CASP3, IRAK1, MAPK1, MAPK3, and YWHAB, respectively. Subsequently, we examined the changes in mRNA expression of immune checkpoint genes and angiogenesis-related genes. qRT-PCR results showed that the majority of immune checkpoint genes (PDCD1, HAVCR2, LGALS9 and VTCN1) and angiogenesis-related genes (PDGFRA, PDGFB, VEGFB and FGFR1) had significantly lower levels of mRNA expression in the GMPS knockdown group than in the control group (Additional file 6: Fig. S6).
The protein expression of the target genes was determined by IHC. The results showed that the expression of CASP3, IRAK1, MAPK1, MAPK3 and YWHAB in HCC tissues was significantly higher than that in adjacent non-tumorous tissues, which is consistent with RNA-sequencing data from the public databases (Additional file 6: Fig. S7).
Despite considerable advances toward the understanding of the molecular mechanism of HCC [22,23,24], it remains a major public health problem, especially in China [25, 26]. The existing prognostic staging system still has many limitations in accurate prognosis prediction and individualized precision therapy [27, 28]. Therefore, continued efforts are needed to find better prognostic signatures to guide individualized treatment so that patients can profit more from precision therapy.
In this study, we built a pyroptosis-related gene signature, knowing that patients with similar clinical characteristics such as the tumor stage may have different outcomes due to the heterogeneity of different epigenetic and genetic backgrounds in tumor subtypes . Our study identified that the risk score was an independent prognostic factor of OS and it was superior to the tumor stage in predicting the OS of HCC patients. Combined with the tumor stage, the prognostic signature performed better in predicting OS. Therefore, the combined use of the five-gene signature and tumor stage may be more conducive to predicting the prognosis of HCC.
CASP3 used to be assumed as an executioner of apoptosis. However, the latest viewpoint proposes that caspase-3 can cause GSDME-mediated pyroptosis [30, 31]. Also, it can promote cancer cell growth, cellular migration, invasiveness, and tumor angiogenesis [32,33,34,35]. IRAK1 is a critical mediator of toll-like receptor and interleukin-1 (IL-1) signaling pathways, playing a crucial role in innate immunity and inflammation. Disruption of these pathways is associated with numerous diseases, including malignancies including HCC [36, 37]. MAPK1 (also known as ERK2) and MAPK3 (also known as ERK1) belong to the MAP kinase family. Both of them are associated with the development and progression of multiple tumors including HCC [38,39,40]. YWHAB encodes a number of 14–3-3 family proteins, of which 14–3-3β regulates multiple signaling pathways in normal and cancer cells . Up-regulation of the 14–3-3β enhances HCC cell migration and proliferation [42, 43]. Collectively, all these prognostic genes in the model have been reported to be involved in cancer initiation and progression and most of them show a close link with HCC, suggesting that these genes can potentially be used as prognostic biomarkers for HCC. Nevertheless, whether these genes affect the prognosis of HCC remains to be elucidated.
Subsequently, GO and KEGG pathway analyses were used to identify the potential role of DEGs. These genes were primarily manifested in cancer-associated pathways, such as cell cycle , focal adhesion [45, 46], ras protein signal transduction , canonical wnt signaling pathway [48, 49], cellular response to hypoxia , Notch signaling pathway [51, 52], VEGF signaling pathway , AMPK signaling pathway , HIF-1 signaling pathway [55, 56], NF-kappa B signaling pathway , mTOR signaling pathway [58, 59], MAPK signaling pathway  and HCC. Of these, HCC is in line with our study subject. Based on these results, we speculated that these DEGs may contribute to the poor prognosis of HCC via activating the above biological pathways. Simultaneously, immune-related pathways, including T cell receptor signaling pathway, B cell receptor signaling pathway, regulation of T cell activation, regulation of T cell differentiation and antigen processing and presentation were significantly enriched, implying that immune dysfunction may also account for poor outcomes in HCC patients.
Additionally, we also noted that cell cycle-related pathways (cell cycle checkpoint, cell cycle G1/S phase transition, cell cycle DNA replication and cell cycle G2/M phase transition) and tumor angiogenesis-related pathways (HIF-1 signaling pathway, VEGF signaling pathway, mTOR signaling pathway and MAPK signaling pathway) were markedly enriched [53, 55, 61,62,63,64]. Aberrant cell cycle is known as a common feature of tumorigenesis . Also, tumor angiogenesis is a pivotal step in tumor growth, invasion and migration . So, we performed a further study on the cell cycle and tumor angiogenesis and found that cell cycle and angiogenesis genes were aberrantly up-regulated in the high-risk group. Furthermore, we discovered that the knockout of prognostic genes in the signature significantly reduced the expression of angiogenesis genes in HCC cells. Accordingly, high-risk scores are likely to be associated with cancer cell proliferation and angiogenesis through regulating the aforementioned genes and pathways. Meanwhile, inhibiting the expression of pyroptosis-related prognostic genes can suppress tumor angiogenesis, thereby improving the prognosis of high-risk populations.
The tumor microenvironment (TME) plays an important role in tumor progression and metastasis . It was found in this study that the infiltration of macrophages and Treg in the high-risk group was significantly higher than that in the low-risk group and the risk scores were positively correlated with the infiltration of macrophages, Treg and cancer associated-fibroblasts, which are known as important components of TME [68, 69]. Macrophages as well as Treg are known to suppress anti-tumor immunity and facilitate tumor progression [69,70,71,72]. Cancer-associated fibroblasts have also been reported to contribute to cancer progression [73, 74]. Intriguingly, immune checkpoint molecules of PD-1, PD-L2, CTLA4, CD80, CD86, HAVCR2, LGALS9, CD276 and VTCN1 were up-regulated in the high-risk group. Moreover, our study demonstrated a significant decrease in the expression of PDCD1, HAVCR2, LGALS9 and VTCN1 after silencing the prognostic genes in the signature. It is known that immune checkpoint inhibitory molecules can facilitate immune escape of cancer cells . Hence, together with the immune dysfunction discussed above, we conclude that induction of immunosuppressive microenvironment seems to be associated with an unfavorable prognosis of high-risk patients, and inhibiting the expression of pyroptosis-related prognostic genes may improve tumor immunosuppression.
Unexpectedly, we discovered that aDCs, iDCs and MHC class I were enriched in the high-risk group, and all of them are related to antigen presentation. Notably, antigen processing and presentation was enriched in the high-risk group suggesting the prognostic genes in our study may change the TME and immune status through affecting antigen presentation, but specific mechanisms warrant further investigation.
Additionally, PD1/PDL1 checkpoint attracted our special attention. To the best of our knowledge, tumor cells can mediate tumor immune escape by utilizing the PD1 / PDL1 checkpoint [76, 77]. Interestingly, recent studies suggested that pyroptosis-induced inflammation could activate anti-tumor immune responses and sensitize cancer cells to anti-PD-1 therapy [78, 79]. More importantly, KEGG analysis revealed that PD-L1 expression and PD-1 checkpoint pathway in cancer was enriched in the high-risk group. Thus, the synergistic effect of pyroptosis induction and PD-1 inhibitors might produce potent anti-tumor effects in patients with high-risk scores, though additional studies are required to verify our supposition.
HCC is extremely resistant to traditional chemotherapeutics, with only certain drugs yielding effective clinical responses . Knowing that chemoresistance is a difficult problem in HCC treatment, we examined the impact of pyroptosis-related genes on chemosensitivity. Of note, the expression levels of target genes were inversely related to chemosensitivity and tumor drug resistance genes (MRP1, MRP4 and MRP5) were overexpressed in the high-risk group. Based on the aforementioned results, it is reasonable to believe that chemotherapy resistance in high-risk patients may be caused by up-regulated expression of MRP1, MRP4 and MRP5. In addition, target genes were significantly enriched in HIF-1 signaling pathway and cellular response to hypoxia in our study. Hypoxia is known to closely contribute to chemoresistance in cancers and HIF1A is an important target for hypoxia-driven drug resistance [81, 82]. Therefore, we speculate that hypoxia is one of the possible causes of chemoresistance in the high-risk group. Chemotherapy drugs can activate caspase 3, thereby specifically cleaving GSDME and ultimately inducing pyroptosis [31, 83]. Recently, a combined therapy consisting of DNA demethylation (promote GSDME expression) and chemotherapy (trigger caspase-3-involved pyroptosis of tumor cells) attracts our sight. This combination strategy could stimulate immune responses through pyroptosis-induced cytokine release and suppress tumor growth, metastasis, and recurrence . It is clear that caspase-3 overexpression played a significant role in this process in the high-risk group, indicating that this combined therapy strategy can better improve the poor outcome in the high-risk group.
In this study, we successfully generated a strong prognostic signature which we believe can help further refine the prognostic predictive power of HCC. This pyroptosis-related signature may provide new insights into the immunity of HCC and suggest a possible direction for individualized treatment of HCC in future.
Availability of data and materials
The datasets analysed during the current study are available in the TCGA repository (https://portal.gdc.cancer.gov/repository) and ICGC repository (https://dcc.icgc.org/projects/LIRI-JP). Additional data can be requested from the corresponding author.
The Cancer Genome Atlas
International Cancer Genome Consortium
The least absolute shrinkage and selection operator
Differentially expressed genes
Receiver operating characteristic
Area under the curve
Kyoto Encyclopedia of Genes and Genomes
Gene set enrichment analysis
Quantitative reverse transcription PCR
T cell regulatory
Anwanwan D, Singh SK, Singh S, Saikam V, Singh R. Challenges in liver cancer and possible treatment approaches. Biochim Biophys Acta Rev Cancer. 2020;1873(1): 188314.
Villanueva A. Hepatocellular carcinoma. N Engl J Med. 2019;380(15):1450–62.
Yang JD, Hainaut P, Gores GJ, Amadou A, Plymoth A, Roberts LR. A global view of hepatocellular carcinoma: trends, risk, prevention and management. Nat Rev Gastroenterol Hepatol. 2019;16(10):589–604.
Craig AJ, von Felden J, Garcia-Lezana T, Sarcognato S, Villanueva A. Tumour evolution in hepatocellular carcinoma. Nat Rev Gastroenterol Hepatol. 2020;17(3):139–52.
He WT, Wan H, Hu L, Chen P, Wang X, Huang Z, et al. Gasdermin D is an executor of pyroptosis and required for interleukin-1β secretion. Cell Res. 2015;25(12):1285–98.
Kayagaki N, Stowe IB, Lee BL, O’Rourke K, Anderson K, Warming S, et al. Caspase-11 cleaves gasdermin D for non-canonical inflammasome signalling. Nature. 2015;526(7575):666–71.
Shi J, Zhao Y, Wang K, Shi X, Wang Y, Huang H, et al. Cleavage of GSDMD by inflammatory caspases determines pyroptotic cell death. Nature. 2015;526(7575):660–5.
Xia X, Wang X, Cheng Z, Qin W, Lei L, Jiang J, et al. The role of pyroptosis in cancer: pro-cancer or pro-"host"? Cell Death Dis. 2019;10(9):650.
Ruan J, Wang S, Wang J. Mechanism and regulation of pyroptosis-mediated in cancer cell death. Chem Biol Interact. 2020;323: 109052.
Loveless R, Bloomquist R, Teng Y. Pyroptosis at the forefront of anticancer immunity. J Exp Clin Cancer Res. 2021;40(1):264.
Hage C, Hoves S, Strauss L, Bissinger S, Prinz Y, Poschinger T, et al. Sorafenib induces pyroptosis in macrophages and triggers natural killer cell-mediated cytotoxicity against hepatocellular carcinoma. Hepatology. 2019;70(4):1280–97.
Zhang X, Zhang P, An L, Sun N, Peng L, Tang W, et al. Miltirone induces cell death in hepatocellular carcinoma cell through GSDME-dependent pyroptosis. Acta Pharm Sin B. 2020;10(8):1397–413.
Fang Y, Tian S, Pan Y, Li W, Wang Q, Tang Y, et al. Pyroptosis: a new frontier in cancer. Biomed Pharmacother. 2020;121: 109595.
Shi J, Gao W, Shao F. Pyroptosis: gasdermin-mediated programmed necrotic cell death. Trends Biochem Sci. 2017;42(4):245–54.
Xu YJ, Zheng L, Hu YW, Wang Q. Pyroptosis and its relationship to atherosclerosis. Clin Chim Acta. 2018;476:28–37.
Xue Y, Enosi Tuipulotu D, Tan WH, Kay C, Man SM. Emerging activators and regulators of inflammasomes and pyroptosis. Trends Immunol. 2019;40(11):1035–52.
Zhaolin Z, Guohua L, Shiyuan W, Zuo W. Role of pyroptosis in cardiovascular disease. Cell Prolif. 2019;52(2): e12563.
Sturm G, Finotello F, Petitprez F, Zhang JD, Baumbach J, Fridman WH, et al. Comprehensive evaluation of transcriptome-based cell-type quantification methods for immuno-oncology. Bioinformatics. 2019;35(14):i436–45.
Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28(11):1947–51.
Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49(D1):D545–51.
Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Bruix J, da Fonseca LG, Reig M. Insights into the success and failure of systemic therapy for hepatocellular carcinoma. Nat Rev Gastroenterol Hepatol. 2019;16(10):617–30.
Park JW, Chen M, Colombo M, Roberts LR, Schwartz M, Chen PJ, et al. Global patterns of hepatocellular carcinoma management from diagnosis to death: the BRIDGE study. Liver Int. 2015;35(9):2155–66.
Zheng R, Qu C, Zhang S, Zeng H, Sun K, Gu X, et al. Liver cancer incidence and mortality in China: temporal trends and projections to 2030. Chin J Cancer Res. 2018;30(6):571–9.
Wang Y, Liao J, Qi W, Xie L, Li Y. Predictive value of conventional ultrasound and contrast-enhanced ultrasound in early recurrence of hepatocellular carcinoma after surgical resection. Ultrasound Med Biol. 2016;42(5):1042–8.
Yuan C, Yuan M, Chen M, Ouyang J, Tan W, Dai F, et al. Prognostic implication of a novel metabolism-related gene signature in hepatocellular carcinoma. Front Oncol. 2021;11: 666199.
Hoshida Y, Moeini A, Alsinet C, Kojima K, Villanueva A. Gene signatures in the management of hepatocellular carcinoma. Semin Oncol. 2012;39(4):473–85.
Huitzil-Melendez FD, Capanu M, O’Reilly EM, Duffy A, Gansukh B, Saltz LL, et al. Advanced hepatocellular carcinoma: which staging systems best predict prognosis? J Clin Oncol. 2010;28(17):2889–95.
Ouyang G, Yi B, Pan G, Chen X. A robust twelve-gene signature for prognosis prediction of hepatocellular carcinoma. Cancer Cell Int. 2020;20:207.
McIlwain DR, Berger T, Mak TW. Caspase functions in cell death and disease. Cold Spring Harb Perspect Biol. 2015;7(4):a008656.
Wang Y, Gao W, Shi X, Ding J, Liu W, He H, et al. Chemotherapy drugs induce pyroptosis through caspase-3 cleavage of a gasdermin. Nature. 2017;547(7661):99–103.
Feng X, Yu Y, He S, Cheng J, Gong Y, Zhang Z, et al. Dying glioma cells establish a proangiogenic microenvironment through a caspase 3 dependent mechanism. Cancer Lett. 2017;385:12–20.
Lauber K, Bohn E, Krober SM, Xiao YJ, Blumenthal SG, Lindemann RK, et al. Apoptotic cells induce migration of phagocytes via caspase-3-mediated release of a lipid attraction signal. Cell. 2003;113(6):717–30.
Mukai M, Kusama T, Hamanaka Y, Koga T, Endo H, Tatsuta M, et al. Cross talk between apoptosis and invasion signaling in cancer cells through caspase-3 activation. Cancer Res. 2005;65(20):9121–5.
Zhao X, Wang D, Zhao Z, Xiao Y, Sengupta S, Xiao Y, et al. Caspase-3-dependent activation of calcium-independent phospholipase A2 enhances cell migration in non-apoptotic ovarian cancer cells. J Biol Chem. 2006;281(39):29357–68.
Cheng BY, Lau EY, Leung HW, Leung CO, Ho NP, Gurung S, et al. IRAK1 augments cancer stemness and drug resistance via the AP-1/AKR1B10 signaling cascade in hepatocellular carcinoma. Cancer Res. 2018;78(9):2332–42.
Singer JW, Fleischman A, Al-Fayoumi S, Mascarenhas JO, Yu Q, Agarwal A. Inhibition of interleukin-1 receptor-associated kinase 1 (IRAK1) as a therapeutic strategy. Oncotarget. 2018;9(70):33416–39.
Lin J, Cao S, Wang Y, Hu Y, Liu H, Li J, et al. Long non-coding RNA UBE2CP3 enhances HCC cell secretion of VEGFA and promotes angiogenesis by activating ERK1/2/HIF-1alpha/VEGFA signalling in hepatocellular carcinoma. J Exp Clin Cancer Res. 2018;37(1):113.
Ng KY, Chan LH, Chai S, Tong M, Guan XY, Lee NP, et al. TP53INP1 downregulation activates a p73-dependent DUSP10/ERK signaling pathway to promote metastasis of hepatocellular carcinoma. Cancer Res. 2017;77(17):4602–12.
Yan Z, Ohuchida K, Fei S, Zheng B, Guan W, Feng H, et al. Inhibition of ERK1/2 in cancer-associated pancreatic stellate cells suppresses cancer-stromal interaction and metastasis. J Exp Clin Cancer Res. 2019;38(1):221.
Xu C, Du Z, Ren S, Liang X, Li H. MiR-129-5p sensitization of lung cancer cells to etoposide-induced apoptosis by reducing YWHAB. J Cancer. 2020;11(4):858–66.
Liu TA, Jan YJ, Ko BS, Chen SC, Liang SM, Hung YL, et al. Increased expression of 14-3-3β promotes tumor progression and predicts extrahepatic metastasis and worse survival in hepatocellular carcinoma. Am J Pathol. 2011;179(6):2698–708.
Tang Y, Lv P, Sun Z, Han L, Zhou W. 14-3-3β Promotes migration and invasion of human hepatocellular carcinoma cells by modulating expression of MMP2 and MMP9 through PI3K/Akt/NF-kappaB pathway. PLoS ONE. 2016;11(1): e0146070.
Liu S, Yang TB, Nan YL, Li AH, Pan DX, Xu Y, et al. Genetic variants of cell cycle pathway genes predict disease-free survival of hepatocellular carcinoma. Cancer Med. 2017;6(7):1512–22.
Bijian K, Lougheed C, Su J, Xu B, Yu H, Wu JH, et al. Targeting focal adhesion turnover in invasive breast cancer cells by the purine derivative reversine. Br J Cancer. 2013;109(11):2810–8.
Nagano M, Hoshino D, Koshikawa N, Akizawa T, Seiki M. Turnover of focal adhesions and cancer cell migration. Int J Cell Biol. 2012;2012: 310616.
Simanshu DK, Nissley DV, McCormick F. RAS proteins and their regulators in human disease. Cell. 2017;170(1):17–33.
Anastas JN, Moon RT. WNT signalling pathways as therapeutic targets in cancer. Nat Rev Cancer. 2013;13(1):11–26.
Clevers H, Nusse R. Wnt/beta-catenin signaling and disease. Cell. 2012;149(6):1192–205.
Mennerich D, Kubaichuk K, Kietzmann T. DUBs, hypoxia, and cancer. Trends Cancer. 2019;5(10):632–53.
Bolos V, Grego-Bessa J, de la Pompa JL. Notch signaling in development and cancer. Endocr Rev. 2007;28(3):339–63.
Yuan X, Wu H, Xu H, Xiong H, Chu Q, Yu S, et al. Notch signaling: an emerging therapeutic target for cancer treatment. Cancer Lett. 2015;369(1):20–7.
Sakurai T, Kudo M. Signaling pathways governing tumor angiogenesis. Oncology. 2011;81(Suppl 1):24–9.
Chhipa RR, Fan Q, Anderson J, Muraleedharan R, Huang Y, Ciraolo G, et al. AMP kinase promotes glioblastoma bioenergetics and tumour growth. Nat Cell Biol. 2018;20(7):823–35.
Hong M, Shi H, Wang N, Tan HY, Wang Q, Feng Y. Dual Effects of Chinese herbal medicines on angiogenesis in cancer and ischemic stroke treatments: role of HIF-1 network. Front Pharmacol. 2019;10:696.
Vadde R, Vemula S, Jinka R, Merchant N, Bramhachari PV, Nagaraju GP. Role of hypoxia-inducible factors (HIF) in the maintenance of stemness and malignancy of colorectal cancer. Crit Rev Oncol Hematol. 2017;113:22–7.
Taniguchi K, Karin M. NF-kappaB, inflammation, immunity and cancer: coming of age. Nat Rev Immunol. 2018;18(5):309–24.
Hay N. The Akt-mTOR tango and its relevance to cancer. Cancer Cell. 2005;8(3):179–83.
Polivka J Jr, Janku F. Molecular targets for cancer therapy in the PI3K/AKT/mTOR pathway. Pharmacol Ther. 2014;142(2):164–75.
Dhillon AS, Hagan S, Rath O, Kolch W. MAP kinase signalling pathways in cancer. Oncogene. 2007;26(22):3279–90.
Li W, Tan D, Zhang Z, Liang JJ, Brown RE. Activation of Akt-mTOR-p70S6K pathway in angiogenesis in hepatocellular carcinoma. Oncol Rep. 2008;20(4):713–9.
Mehta VB, Besner GE. HB-EGF promotes angiogenesis in endothelial cells via PI3-kinase and MAPK signaling pathways. Growth Factors. 2007;25(4):253–63.
Munoz-Chapuli R, Quesada AR, Angel MM. Angiogenesis and signal transduction in endothelial cells. Cell Mol Life Sci. 2004;61(17):2224–43.
Saxton RA, Sabatini DM. mTOR signaling in growth, metabolism, and disease. Cell. 2017;168(6):960–76.
Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell. 2011;144(5):646–74.
Albini A, Tosetti F, Li VW, Noonan DM, Li WW. Cancer prevention by targeting angiogenesis. Nat Rev Clin Oncol. 2012;9(9):498–509.
Quail DF, Joyce JA. Microenvironmental regulation of tumor progression and metastasis. Nat Med. 2013;19(11):1423–37.
Curry JM, Sprandio J, Cognetti D, Luginbuhl A, Bar-ad V, Pribitkin E, et al. Tumor microenvironment in head and neck squamous cell carcinoma. Semin Oncol. 2014;41(2):217–34.
Gajewski TF, Schreiber H, Fu YX. Innate and adaptive immune cells in the tumor microenvironment. Nat Immunol. 2013;14(10):1014–22.
Degroote H, Van Dierendonck A, Geerts A, Van Vlierberghe H, Devisscher L. Preclinical and clinical therapeutic strategies affecting tumor-associated macrophages in hepatocellular carcinoma. J Immunol Res. 2018;2018:7819520.
Mantovani A, Marchesi F, Malesci A, Laghi L, Allavena P. Tumour-associated macrophages as treatment targets in oncology. Nat Rev Clin Oncol. 2017;14(7):399–416.
Smyth MJ, Teng MW, Swann J, Kyparissoudis K, Godfrey DI, Hayakawa Y. CD4+CD25+ T regulatory cells suppress NK cell-mediated immunotherapy of cancer. J Immunol. 2006;176(3):1582–7.
Kato T, Noma K, Ohara T, Kashima H, Katsura Y, Sato H, et al. Cancer-associated fibroblasts affect intratumoral CD8(+) and FoxP3(+) T cells via IL6 in the tumor microenvironment. Clin Cancer Res. 2018;24(19):4820–33.
Madar S, Goldstein I, Rotter V. ‘Cancer associated fibroblasts’–more than meets the eye. Trends Mol Med. 2013;19(8):447–53.
Pardoll DM. The blockade of immune checkpoints in cancer immunotherapy. Nat Rev Cancer. 2012;12(4):252–64.
Li CW, Lim SO, Xia W, Lee HH, Chan LC, Kuo CW, et al. Glycosylation and stabilization of programmed death ligand-1 suppresses T-cell activity. Nat Commun. 2016;7:12632.
Tang H, Wang Y, Chlewicki LK, Zhang Y, Guo J, Liang W, et al. Facilitating T cell infiltration in tumor microenvironment overcomes resistance to PD-L1 blockade. Cancer Cell. 2016;30(3):500.
Xiao Y, Zhang T, Ma X, Yang QC, Yang LL, Yang SC, et al. Microenvironment-responsive prodrug-induced pyroptosis boosts cancer immunotherapy. Adv Sci (Weinh). 2021;8(24): e2101840.
Zhang S, Zhang Y, Feng Y, Wu J, Hu Y, Lin L, et al. Biomineralized two-enzyme nanoparticles regulate tumor glycometabolism inducing tumor cell pyroptosis and robust antitumor immunotherapy. Adv Mater. 2022;34:e2206851.
Lohitesh K, Chowdhury R, Mukherjee S. Resistance a major hindrance to chemotherapy in hepatocellular carcinoma: an insight. Cancer Cell Int. 2018;18:44.
Lu J, Tan M, Cai Q. The Warburg effect in tumor progression: mitochondrial oxidative metabolism as an anti-metastasis mechanism. Cancer Lett. 2015;356(2 Pt A):156–64.
Maiso P, Huynh D, Moschetta M, Sacco A, Aljawai Y, Mishima Y, et al. Metabolic signature identifies novel targets for drug resistance in multiple myeloma. Cancer Res. 2015;75(10):2071–82.
Rogers C, Fernandes-Alnemri T, Mayes L, Alnemri D, Cingolani G, Alnemri ES. Cleavage of DFNA5 by caspase-3 during apoptosis mediates progression to secondary necrotic/pyroptotic cell death. Nat Commun. 2017;8:14128.
Fan JX, Deng RH, Wang H, Liu XH, Wang XN, Qin R, et al. Epigenetics-based tumor cells pyroptosis for enhancing the immunological effect of chemotherapeutic nanocarriers. Nano Lett. 2019;19(11):8049–58.
The project was supported by the National Natural Science Foundation of China (No. 81970527). The funders had no role in the study design, data collection, data analysis, interpretation and writing of the report.
Ethics approval and consent to participate
All procedures performed in this study involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards. This study was approved by the Review of Ethics Committee in Clinical Research of the First Affiliated Hospital of Wenzhou Medical University. Written informed consents were obtained from all patients for the use of the biospecimens for research purposes.
Consent for publication
The authors have declared that no competing interest exists.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1.
The description of pyroptosis related genes.
Additional file 2.
Five kinds of immune infiltration algorithm principles.
Additional file 3.
The siRNA sequences used in this study.
Additional file 4.
The primer sequences used in this study.
Additional file 5.
Information of the antibody used in IHC.
Additional file 6.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Lin, D., Chen, Z., Zeng, Y. et al. A pyroptosis-related gene signature provides an alternative for predicting the prognosis of patients with hepatocellular carcinoma. BMC Med Genomics 16, 2 (2023). https://doi.org/10.1186/s12920-023-01431-z
- Hepatocellular carcinoma (HCC)
- Gene signature
- Overall survival
- Immune status
- Drug sensitivity