Xiao-Wei Fu1,2,3,4
Chun-Qing Song2,3,4*- 1Fudan University, Shanghai, China
- 2Key Laboratory of Growth Regulation and Translational Research of Zhejiang Province, School of Life Sciences, Westlake University, Hangzhou, China
- 3Westlake Laboratory of Life Sciences and Biomedicine, Hangzhou, China
- 4Laboratory of Gene Therapeutic Biology, Institute of Basic Medical Sciences, Westlake Institute for Advanced Study, Hangzhou, China
Background: Hepatocellular carcinoma (HCC) is characterized by a poor prognosis and accounts for the fourth common cause of cancer-related deaths. Recently, pyroptosis has been revealed to be involved in the progression of multiple cancers. However, the role of pyroptosis in the HCC prognosis remains elusive.
Methods: The clinical information and RNA-seq data of the HCC patients were collected from the TCGA-LIHC datasets, and the differential pyroptosis-related genes (PRG) were firstly explored. The univariate Cox regression and consensus clustering were applied to recognize the HCC subtypes. The prognostic PRGs were then submitted to the LASSO regression analysis to build a prognostic model in the TCGA training cohort. We further evaluated the predictive model in the TCGA test cohort and ICGC validation cohort (LIRI-JP). The accuracy of prediction was validated using the Kaplan—Meier (K-M) and receiver operating characteristic (ROC) analyses. The single-sample gene set enrichment analysis (ssGSEA) was used to determine the differential immune cell infiltrations and related pathways. Finally, the expression of the prognostic genes was validated by qRT-PCR in vivo and in vitro.
Results: We identified a total of 26 differential PRGs, among which three PRGs comprising GSDME, GPX4, and SCAF11 were subsequently chosen for constructing a prognostic model. This model significantly distinguished the HCC patients with different survival years in the TCGA training, test, and ICGC validation cohorts. The risk score of this model was confirmed as an independent prognostic factor. A nomogram was generated indicating the survival years for each HCC patient. The ssGSEA demonstrated several tumor-infiltrating immune cells to be remarkably associated with the risk scores. The qRT-PCR results also showed the apparent dysregulation of PRGs in HCC. Finally, the drug sensitivity was analyzed, indicating that Lenvatinib might impact the progression of HCC via targeting GSDME, which was also validated in human Huh7 cells.
Conclusion: The PRG signature comprised of GSDME, GPX4, and SCAF11 can serve as an independent prognostic factor for HCC patients, which would provide further evidence for more clinical and functional studies.
Introduction
Statistics indicate that by 2025, there will be more than 1 million new cases of liver cancer annually, posing a significant challenge for the global medical field (Llovet et al., 2021). The HCC ranks as the fourth leading cause of cancer-related mortality accounting for 90% of liver cancer cases (Llovet et al., 2016). Recent studies indicated the 5-years survival rate of HCC patients to be lessened by 20% globally and as low as 12% in Asian countries (Craig et al., 2020). Epidemiological studies reported the hepatitis B virus, hepatitis C virus, alcoholism, and aflatoxin as common inducing factors of HCC (Yang et al., 2019).
HCC develops rapidly as well as stealthily such that the patients are diagnosed only when the disease has progressed to the middle and late stage. Alpha-fetoprotein (AFP) is currently a gold standard for the diagnosis and prognosis of HCC. However, its sensitivity and specificity are not very high owing to the interference in the expression by many non-HCC related factors (Cai et al., 2019). In addition, the accuracy of prognosis is also affected by the heterogeneity of HCC (Liu et al., 2016). Therefore, it is urgent to develop a novel prognostic model to improve the accuracy of the prognostic judgment for HCC patients.
Cell death is one of the critical aspects of anti-tumor drug research. It may involve various patterns, including pyroptosis, apoptosis, necrosis, necroptosis, and ferroptosis (Shojaie et al., 2020). Pyroptosis is a novel lytic and pro-inflammatory programmed cell death. It is mediated by the cysteine aspartate specific protein kinases (caspases) 1, 4, 5, and 11 (Galluzzi et al., 2018), and the final step is dependent on the activity of the gasdermin (GSDM) family proteins to form pores in the cell membrane (Kovacs and Miao, 2017; Humphries et al., 2020). Pyroptosis is morphologically characterized by cell swelling, plasma membrane permeability, and the gradual release of the inflammatory factors (Rathinam and Fitzgerald, 2016), while necrosis makes the cytoplasmic membrane rupture with a blast. Pyroptosis is also different from apoptosis because apoptosis does not induce cell membrane breakdown and inflammatory response (Bertheloot et al., 2021). In addition, pyroptosis differs from necroptosis, another inflammatory programmed cell death executed by the mixed lineage kinase domain-like protein (MLKL), in that pyroptosis maintains mitochondria integrity but necroptosis does not (Bertheloot et al., 2021). Furthermore, pyroptosis differs from ferroptosis because ferroptosis is featured by iron-dependent oxidative perturbations, increased membrane density, and small mitochondria (Sun et al., 2021). As studies continue, pyroptosis has been found to plays a vital role in tumor formation and development. Due to its dual functions of resisting infection and inducing pathological inflammation, pyroptosis has dual roles in promoting tumors and changing tumor immune microenvironments (Kong et al., 2015; Rathinam and Fitzgerald, 2016; Xia et al., 2019). Specifically, the role of pyroptosis in the HCC development and prognosis remains elusive.
This study systematically analyzed the differential expressions of PRGs between the HCC and normal samples; explored the clinical prognostic value of these genes through Cox expression analysis; established an independent prognostic model based on PRGs; investigated the relationship between the pyroptosis and tumor immune microenvironments; validated the mRNA expressions of PRGs in vivo and in vitro; and evaluated the drug sensitivity of these prognostic factors. Therefore, this study provides potential targets for the prognosis and treatment of HCC patients.
Materials and Methods
Data Acquisition
The RNA sequence data and related clinical information of 374 liver cancer patients (TCGA-LIHC) were acquired from the TCGA website (https://portal.gdc.cancer.gov/repository). The gene expression data were normalized by scale method using the “limma” package (Yuan C. et al., 2021). After excluding the missing clinical information of patients, 370 HCC patients were randomly separated into the training and the test groups by the “caret” package. Besides, transcriptomics data with clinical features of 231 HCC patients (LIRI-JP) were downloaded from the ICGA database (https://dcc.icgc.org/projects/LIRI-JP). These HCC patients were HBV or HCV carriers from Japan. The gene read count values of these patients were also normalized, and both the TCGA and ICGC data were public.
Analysis of Differential PRGs
A total of thirty-three PRGs were retrieved from the previous literature and are listed in Supplementary Table S1 (Man and Kanneganti, 2015; Tang et al., 2020; Ye J. et al., 2021). The differentially expressed genes (DEGs) were analyzed by the “limma” package in the R software and visualized by the heatmap and volcano plot with adjusted p-value < 0.05 in the TCGA cohort. In addition, a protein-protein interaction network for the PRGs was generated using the software GENEMANIA (http://genemania.org/).
Establishment and Validation of the Prognostic Model Based on the PRGs
The Univariate Cox analysis was employed to screen the PRGs with the prognostic value, and the cut-off p-value was set at 0.05, and 17 survival-related genes were used for further study. The prognostic model was established to minimize overfitting using the LASSO-penalized Cox regression analysis via the “glmnet” R package (Du et al., 2021). Eventually, the three genes and their coefficients were retained, and the minimum criteria determined the penalty parameter (λ). The risk score was obtained using the formula: risk score = (βA × Gene A expression) + (βB × Gene B expression) ⋯ + (βN × Gene N expression), in which β represents regression coefficient (Yang et al., 2021). The HCC patients were separated into the high- and low-risk groups based on the median risk score. Then, the principal component analysis (PCA) for the two risk groups was constructed using the “limma” and “scatterplot3d” R packages in terms of gene expressions in the prognostic model (Yuan M. et al., 2021). The survival analysis between the two risk groups was carried out using the “survminer” R package, and the ROC curve analysis was performed via the “survival” and “timeROC” R packages (Fang et al., 2021). Besides, the univariate and multivariate Cox regression was utilized to determine the independent prognostic value of the 3-gene signature.
To validate the efficiency of our model, the patients in the TCGA internal test cohort or ICGC external validation cohort (LIRI-JP) were applied. The mRNA levels of PRGs were normalized according to the “scale” function, and the risk score was calculated using the same formula applied in the TCGA training cohort. The TCGA test or ICGC cohort patients were separated into the low- and high-risk groups using the median risk score obtained from the TCGA training cohort.
Development of a Predictive Nomogram
Based on the risk score and different clinical features (gender, age, histologic grade, and pathological stage), a nomogram model was established to predict the survival years for the HCC patients using the “rms” and “survival” packages (Dai et al., 2021).
Gene Set Enrichment Analyses
The GSEA analysis was carried out between the two risk groups using the GSEA software 4.0.1 to identify the differential KEGG pathways. The normalized enrichment scores and nominal p-values were determined for analyzing the enrichment levels and statistical significance. In addition, the infiltrating scores of the 16 immune cells and activities of 13 immune-related pathways were analyzed using ssGSEA provided in the “GSVA” R package (Zhao et al., 2021).
Drug Sensitivity
The NCI-60 dataset covering nine cancer categories was available on the CellMiner homepage (https://discover.nci.nih.gov/cellminer) (Shankavaram et al., 2007; Shankavaram et al., 2009). In addition, Pearson correlation analysis was performed to indicate the relevance between the independent prognostic PRGs and drug sensitivity. Drugs used in this sensitivity analysis are those approved by the FDA or those in clinical tests.
Construction of the Hepatocellular Carcinoma Mouse Model
The FVB mice (8 weeks) were purchased from Westlake University (Zhejiang, China). Plasmids for injection were harvested using the Endotoxin Free Maxiprep kit (Qiagen) and delivered to the FVB mice as a mixture using the hydrodynamic tail vein injection for HCC formation as described before (Song et al., 2017). Dosages of these plasmids were: px330-U6-sgP53 (mouse) 20 μg, pT3-EF1α-c-Myc (human) 5 μg, pCMV-sleeping beauty transposase 2 μg. The daily abdominal palpation of mice was observed, and the mice were then sacrificed when they suffered from high burdens of liver tumors. The mice were raised according to the protocols approved by the Institutional Animal Care and Use Committee of Westlake University. The researcher followed the standard biosecurity and institutional safety procedures in this study.
Cell Culture and Drug Treatment
The human normal liver cell line HL-7702 and HCC cell lines (SK-Hep1 and Huh7) were cultured in a DMEM medium containing 10% FBS in a 5% CO2 incubator at 37°C. In addition, huh7 cells were treated by Lenvatinib (HY-10981, MCE) with different concentrations (0, 0.5, 1, 5, 10 μM) for 48 h based on the previous study (Jin et al., 2021).
RNA extraction and quantitative real-time reverse transcriptase-polymerase chain reaction (qRT-PCR).
The total RNA in the liver tissues or liver cell lines was extracted using the TRIZOL (Invitrogen) reagent. Then, the cDNA was obtained by reverse transcription using the cDNA Synthesis Mix (Novoprotein, E047-01B) and analyzed using quantitative PCR (Novoprotein, E096-01B). Finally, the mRNA expression levels were normalized with the ACTB gene. All the primers used in this study are listed in Supplementary Table S2.
Western Blotting
The total proteins in Huh7 cells treated with Lenvatinib were extracted by RIPA buffer containing 1% PMSF (Beyotime) and quantified by the BCA method (Beyotime). Twenty microgram proteins were used for western blotting against GSDME (1:2000, ab215191, Abcam) and Hsp90 (1:2000, 610419, BD Bioscience).
Statistical Analysis
The one-way ANOVA was adopted to compare the DEGs between the HCC and normal liver tissues. The Mann—Whitney test using the BH method adjusted p-value was adopted to measure the ssGSEA scores. The entire statistical analyses were achieved using the R software 4.0.1. The statistical significance was considered as a P -value less than 0.05 if not otherwise specified.
Results
Differential Expression of PRGs Between the Normal and Tumor Tissues
We performed the differential expression analysis of 33 PRGs involved between the normal liver and HCC samples in the TCGA database (Figure 1A). The 33 PRGs were chosen based on their roles in pyroptosis according to previous studies (Karki and Kanneganti, 2019; Ye Y. et al., 2021; Ju et al., 2021). The results were presented in the heatmaps and volcano plots (Figures 1B,C). We identified 26 differential genes (p < 0.05), among which the three genes were significantly downregulated in the HCC samples, including IL6, IL-1β, and NLRP3. Conversely, another 23 genes were significantly upregulated, including GSDMC, PLCG1, PYCARD, GSDME, NLRP1, GSDMB, GSDMD, CASP8, NOD1, CASP8, NOD1, CASP3, GSDMA, PJVK, TIRAP, NOD2, AIM2, GPX4, NLRP7, CASP9, CASP6, PRKACA, SCAF11, CASP4, and NLRP6. The interactions within these PRGs were identified using GeneMANIA and constituted three interconnected subnetworks: the NLR superfamily regulatory network, the caspase family regulatory network, and the GSDM family regulatory network (Figure 1D). Finally, the correlation for these PRGs was analyzed and indicated in Figure 1E.
FIGURE 1. Screening of the PRGs in the HCC patients from the TCGA database. (A) Flow chart of the study. (B) Heatmap of the differential gene expression between the normal and the tumor tissues. (C) Volcano plot of the differentially expressed PRGs. (D) PPI network of the interactions among PRGs. (E) The correlation network (Redline: positive correlation; Blueline: negative correlation. The depth of the colors reflects the strength of the relevance).
Tumor Classification Based on the Prognostic Pyroptosis Regulators
We first analyzed the prognostic values of PRGs in the TCGA cohort with the Univariate Cox regression analysis. The high expression of CASP1, CASP3, CASP5, CASP6, CASP8, GPX4, GSDMA, GSDME, NLRC4, NLRP3, NLRP7, NOD1, NOD2, PLCG1, and SCAF11 correlated with the poor survival of the HCC patients, as indicated by the Hazard ratio (HR) > 1 (Figure 2A). Then, we performed the consensus clustering analysis to investigate the relationship between these prognostic genes and HCC subtypes. According to the CDF value, we classified the 370 HCC patients into two clusters (k = 2, Figures 2B–D), and we found that the patients from cluster 1 tended to survive longer than the patients from cluster 2 (Figure 2E), implying a significant prognostic value of these PRGs. Furthermore, the two clusters did not differ in the clinical parameters such as stage, grade, gender, age, and TMN (Figure 2F).
FIGURE 2. Identification of the molecular subtypes of the HCC patients using the PRGs associated with prognosis. (A) The univariate Cox Regression Analysis based forest plot in PRGs. (B) The HCC patients were stratified into 2 clusters based on the consensus clustering matrix (k = 2). (C, D) Consensus clustering model with cumulative distribution function (CDF) by k from 2 to 9. (E) The Kaplan—Meier curves (KM) of the overall survival in the two HCC clusters. (F) Heatmap with the correlation between the two groups and their clinicopathologic characters.
Establishment of an Independent Prognostic Risk Model in the TCGA Training Cohort
We split the samples in the TCGA-LIHC dataset into two equal cohorts: training cohort and test cohort at random. We found no significant difference between the training and test cohorts from TCGA-LIHC in the leading clinical indicators (Supplementary Table S3). We performed lasso regression analysis using 17 prognostic genes in the training cohort to build the prognostic model. According to the minimum criteria, a risk model consisting of GSDME, GPX4, and SCAF11 was built (Figures 3A,B). The risk score was obtained by the formula: risk score = (0.0182*GSDME exp.) + (0.0005*GPX4 exp.) + (0.0188*SCAF11 exp.)
FIGURE 3. Construction of the pyroptosis genes-based independent prognostic model in the TCGA training cohort. (A, B) Construction of the LASSO regression model based on the 17 predictive genes in the TCGA training cohort. (C–E) Distribution of the risk scores, survival status, and expression of the three pyroptosis-related risk genes. (F) The KM analysis of the overall survival in the high-risk and low-risk groups. (G) The ROC analysis to evaluate the predictive efficiency. (H) The univariate Cox analysis to assess the independence of three pyroptosis-related risk genes. (I) Multivariate Cox analysis to assess the independence of three pyroptosis-related risk genes. (J–L) PCA analysis of the high-risk and low-risk groups based on all genes, pyroptosis genes, and three risk genes.
We separated the patients in the training cohort into the high- and low-risk groups based on the median risk score (Figure 3C). The high-risk group had more deaths and shorter survival years (Figure 3D). The heatmap analysis indicated that high-risk patients had increased expression levels of three risk genes (Figure 3E). In addition, the Kaplan-Meier curve indicated that the high-risk patients had worse overall survival (OS) than low-risk patients (p < 0.001). The subsequent ROC analysis demonstrated that this 3-gene risk model could robustly evaluate and predict the survival of the HCC patients (AUC = 0.763, Figure 3G). Furthermore, the univariate and multivariate Cox regression analyses determined whether the risk score derived from the prognostic risk model could act as an independent prognostic indicator. In the univariate Cox regression analysis, the risk score (p < 0.001, HR = 3.996, 95% CI: 2.253–7.088, Figure 3H) was a potential hazard factor. The multivariate analysis also indicated that the risk score could serve as an independent prognostic factor (p < 0.001, HR = 3.282, 95% CI: 1.787–6.026, Figure 3I). Finally, the PCA plot indicated that the high- and low-risk groups could be well-separated with the three risk genes, but not all the genes or all the pyroptosis genes (Figures 3J–L).
Internal and External Validation of the Risk Signature
The efficiency of the risk model was validated in a TCGA internal test cohort which included 184 patients and an ICGC external validation cohort. In the TCGA internal test cohort, we divided patients into high or low-risk groups according to the median risk score calculated by the formula in the training cohort. The risk score is an independent prognostic factor in the test cohort as indicated by univariate and multivariate Cox regression analysis (p = 0.001, HR = 3.225, 95% CI = 1.566–6.643 for univariate, Figure 4A; p = 0.01, HR = 2.843, 95% CI = 1.288–6.277 for multivariate, Figure 4B). The high-risk group tended to have more deaths and higher expression of risk genes (Figures 4C–E). The KM curve indicated that overall survival was lower in the high-risk group (p = 0.015, Figure 4F).
FIGURE 4. Internal validation of the risk model in the TCGA test cohort. (A, B) Forest plots of Univariate Cox analysis and Multivariate Cox analysis. (C–E) Distribution of risk scores, survival status, and expression of the three pyroptosis-related risk genes. (F) The KM analysis of the overall survival in the high-risk and low-risk groups. (G) The ROC analysis to estimate the predictive efficiency. (E) Univariate analysis. (H) PCA analysis of the high-risk and low-risk groups based on the three risk genes.
Moreover, the ROC curve suggested that the model exhibited an excellent predictive capability (AUC = 0.677, Figure 4G). PCA analysis indicated that the expression levels of the 3-risk gene could separate high-risk patients from low-risk in the TCGA test cohort (Figure 4H). We split the ICGC external validation cohort into the high- and low-risk groups based on the risk score (Figure 5A). The low-risk group had fewer deaths and lower expression of the risk genes (Figures 5B,C). The KM curve and ROC analysis suggested that the overall survival of the low-risk group was higher (p = 0.02, Figure 5D), and the model was reliable (AUC = 0.638, Figure 5E). Finally, the PCA plot indicated that the risk genes were well able to separate the two risk groups (Figure 5F).
FIGURE 5. External validation of the risk model in the ICGC cohort. (A–C) Distribution of the risk scores, survival status, and expression of the three pyroptosis-related risk genes in the ICGC cohort. (D) The KM analysis of the overall survival in the high-risk and low-risk groups. (E) The ROC analysis to estimate the predictive efficiency. (F) The PCA analysis of the high-risk and low-risk groups based on the three risk genes.
Stratification Analysis of the Independent Prognostic Signature
We separated the patients in the TCGA-LIHC dataset into several subgroups according to the different clinical parameters. First, we investigated whether the high- and low-risk patients determined differences in survival. The clinical stratifications for the study included age (>65 vs ≤ 65), gender (female vs male), tumor grade (G3/4 vs G1/2), and AJCC stage (I/II vs III/IV). The KM curve showed that the high-risk patients had a poorer survival probability than the low-risk patients under the condition of age >65, female, male, G1–G2, G3–G4, or Stage I–II, except for age ≤65 or stage III-IV (Figures 6A–H).
FIGURE 6. Identification of the HCC patients suitable for the risk model in the TCGA cohort. KM analysis of the overall survival in the high-risk and low-risk groups based on age (A, B), gender (C, D), grade (E, F), and stage (G, H).
Establishment of a Prognostic Nomogram for Hepatocellular Carcinoma
We developed a novel prognostic nomogram to offer a reliable and quantifiable method for predicting the survival of the HCC patients based on the risk scores and clinical features, such as age, gender, grade, and stage (Figure 7). The nomogram could effectively predict the probability of the 1, 3, and 5-years overall survival in the HCC patients.
FIGURE 7. Nomogram containing the risk score to predict the overall survival in the HCC patients.
Gene Set Enrichment Analysis and Immune Activity Between the Subgroups
Given clinical information, the high- and low-risk groups were significantly correlated with the immune scores (Figure 8A). Specifically, the patients with lower immune scores favored high-risk scores compared to those with higher immune scores. In addition, the patients with higher stages picked higher risk scores (Figures 8B–E). Interestingly, the high-risk group had a high expression level of PD-L1 (Figure 8F). Furthermore, the GSEA analysis significantly enriched the immune-related signaling pathways, such as the B cell and T cell receptor signaling pathways (Figure 8G).
FIGURE 8. GSEA and immune correlation analysis of the immune cells and related immune pathways in the TCGA cohort. (A) Heatmap with correlations between clinical features and risk groups (*p < 0.05). For T, M, N, stage, and grade, the higher the numbers, the more advanced the tumors. For cluster, cluster 1 and cluster 2 were divided based on consensus clustering model. The risk and ImmuneScore were divided into high and low based on the median of all patients. (B–E) Risk scores of the HCC patients are classified by cluster, stage, Tstage, and immune scores. (F) PD-L1 expression in the low-risk and high-risk groups. (G) GSEA of the significantly enriched KEGG pathways in the TCGA cohort classified by independent prognostic genes. (H, I) The ssGSEA scores of 16 kinds of immune cells and 13 immune pathways (*p < 0.05, **p < 0.01, ***p < 0.001).
Further, the ssGSEA analyzed the differences in 16 kinds of immune cell infiltrations and 13 types of immune signal pathways. The high-risk group had decreased infiltrations of the B cells, mast cells, NK cells, and TIL and increased infiltrations of the DCS and macrophage compared to the low-risk group (p < 0.05, Figure 8H). Besides, the high-risk group exhibited suppression in the immune pathways, including the cytolytic activity and type II IFN response (p < 0.05, Figure 8I).
Validation of the Differential Expression of the Independent Prognostic Genes
We analyzed the differential mRNA levels of the three independent prognostic genes in the normal liver tissues (from TCGA and GTEx) and HCC tissues (from TCGA) by GEPIA. The HCC tissues favored an increased mRNA level of GSDME, GPX4, and SCAF11 (Figures 9A–C). In addition, we examined the protein levels of these three genes using the HPA database. The results showed that the HCC tissues had higher protein levels of GSDME, GPX4, and SCAF11 (Figures 9D–I).
FIGURE 9. Expression of the independent prognostic genes. (A–C) The gene mRNA expressions of GPX4, GSDME, SCAF11 in the normal and tumor groups (*p < 0.05). (D–I) Immunohistochemistry of the GPX4, GSDME, SCAF11 in the normal and tumor groups from the HPA database.
We further validated the expression levels of GSDME, GPX4, and SCAF11 in a mouse HCC model, which was constructed by knocking out the p53 and overexpressing the myc in the mouse liver (Figure 10A). The HCC model formed multiple tumors as shown in Figure 10B. As expected, the mouse HCC livers had remarkably increased GSDME, GPX4, and SCAF11 mRNA levels compared to the normal livers (Figures 10C–E).
FIGURE 10. Validation of the independent prognostic genes in vivo and in vitro. (A) HCC Model design. FVB mice were injected with normal saline or sgP53/c-Myc/SB plasmids, respectively. (B) Representative liver tissues of normal (left) and HCC mice (middle) and tumor numbers (right). (C–E) The mRNA expressions of GSDME, GPX4, and SCAF11 in the normal and HCC livers. (F–H) The mRNA expressions of GSDME, GPX4, and SCAF11 in the human normal liver cell HL-7701 and liver cancer cell SK-Hep1. *p < 0.05, **p < 0.01, ***p < 0.001.
In addition, we confirmed the mRNA levels of these independent prognostic genes in human culture cells. As shown in Figures 10F–H, the HCC cell line SK-Hep1 expressed significantly higher mRNA levels of GSDME, GPX4, and SCAF11 compared to the normal liver cell line HL-7702.
Drug Sensitivity Analysis
By analyzing the CellMiner database, the potential drugs were found to be correlated to these independent prognostic genes (Supplementary Table S4). Among the top 16 gene-drug correlations, 15 correlations pointed to the GSDME; the other correlation was SCAF11 (Figure 11). Surprisingly, the HCC drug lenvatinib positively correlated with the expression of GSDME (Cor = 0.453, p < 0.001, Figure 11K). To confirm this correlation, we treated HCC cell line Huh7 with Lenvatinib. As was expected, Lenvatinib treatment upregulated both the mRNA and protein levels of GSDME. Besides, Lenvatinib was able to induce the active form of GSDME (GSDME N-terminal).
FIGURE 11. Sensitivity correlation analysis between the independent prognostic genes and drugs based on the CellMiner Database. (A–P) Correlation analysis. (Q) Relative mRNA expression of GSDME following treatment with Lenvatinib for 48 h in Huh7. (R) Analysis of GSDME and GSDME-N terminal in Huh7 treated with Lenvatinib for 48 h by western blotting. *p < 0.05, **p < 0.01, ***p < 0.001.
Discussion
In recent years, the increasing incidence of HCC has necessitated early diagnosis and prognosis as being incredibly significant for the survival of the patients. At present, the diagnosis and prognosis of HCC patients are mainly based on the pathological evaluation, AJCC TNM, and BCLC stage (Bruix et al., 2016). However, these diagnostic and prognostic methods are not sensitive enough. Therefore, finding efficient diagnostic and prognostic markers is necessary for helping the HCC patients in improving their clinical outcomes. Moreover, the mechanisms underlying the pathogenesis of HCC have not been fully elucidated (Ye J. et al., 2021). To devise precision medicine, improving the prognosis and survival rate of HCC patients, and discovering the molecular mechanisms underlying the development and progression of HCC, accurate and reliable prognostic models should be developed based on the novel biomarkers (Yan et al., 2021). Although the role of pyroptosis in tumors has been gradually explored (Wang et al., 2017; Blasco and Gomis, 2020), it remains unclear mainly in HCC.
This study systematically compared the altered expressions of PRGs in HCC patients and revealed that three genes were significantly downregulated while 23 genes were upregulated. The consensus clustering analysis identified 17 prognostic genes, which helped the HCC patients be divided into two subtypes. Interestingly, the two subtypes were indistinguishable by existing clinical criteria but had significant differential survival rates, proving instrumental for clinical typing. Subsequently, to evaluate the prognosis of the HCC patients, we established an independent prognostic model by performing the LASSO and Cox analysis of the prognosis-related genes, which was well validated in the internal test and external validation cohorts. The model yielded significant survival differences among the patients with different clinical characteristics, except for those with age ≤65 and in stage III-IV. Based on this risk model, a nomogram was drawn to predict the HCC patients’ overall survival. Besides, the high-risk group had a lower immune score and higher expression level of PD-L1 than the low-risk group, indicating the potential differences in the immune function between the two risk groups. Later, the GSEA analysis of the two risk groups revealed that the T cell receptor and B cell receptor signaling pathways were enriched. Furthermore, the ssGSEA analysis suggested that the high-risk group is characterized by the lower levels of immune cell infiltration (B cells, Mast cells, NK cells, and TIL cells) and immune pathways (Cytolytic activity and Type II IFN response). Considering the independent prognostic model was constructed using GSDME, GPX4, and SCAF11, we then validated their differential expressions in the normal and HCC tissues. Using the RNA-seq data from the TCGA and GTEx, and protein expressions from HPA, all the three genes showed higher mRNA and protein levels in the HCC tissues than those in the normal tissues. Furthermore, the expression of these three risk genes was validated in the HCC mouse model and human HCC cell lines. Finally, we analyzed the correlation between FDA-approved drugs and these three targets using CellMiner database. The results revealed that among the 16 most significant correlations, seven drugs showed significantly positive correlations with GSDME, eight drugs conversely correlated with GSDME, and one drug positively correlated with SCAF11. We finally validated the correlation of Lenvatinib and GSDME in mRNA and protein levels.
Since the term “pyroptosis” was raised in 2001 by D’ Souza et al., many studies have focused on this novel pro-inflammatory programmed cell death (Fink and Cookson, 2005). In recent years, pyroptosis has gained increasing prominence in tumor research (Wang et al., 2017; Yu et al., 2019; Wu et al., 2021). However, pyroptosis demonstrates both pro-tumor and anti-tumor functions (Xia et al., 2019). The induction of the cancer cells towards pyroptosis can suppress tumor development, proving to be a promising target for drug discovery. Also, the inflammatory molecules released by the cancer cells undergoing pyroptosis can gradually transform the surrounding normal tissues into cancer cells by changing microenvironments. However, little research has focused on the function and prognosis of multiple PRGs in HCC development (Chu et al., 2016). An independent prognosis model was constructed with three PRGs, including GSDME, GPX4, and SCAF11. GSDME, formerly called DFNA5, was first confirmed in an extended Dutch family with autosomal dominant nonsyndromic hereditary hearing loss (Van Laer et al., 2004). Then, in 2017, it was identified as a new executor of pyroptosis (Wang et al., 2017). GSDME is cleaved by caspase-3 producing the N-terminal fragment (GSDME-NT) that converts the death pathway from apoptosis to pyroptosis. Recently, new studies have recognized GSDME as a conduit for the release of IL-1β to the surrounding microenvironment independent of its capability of triggering cell death (Zhou and Abbott, 2021). Ever since the mechanism of GSDME mediated pyroptosis was revealed, an increasing number of studies have been focused on its role in cancer. In most cancer types, such as breast cancer, colorectal cancer, gastric cancer, and bladder cancer, GSDME has a higher expression level in the normal tissues and is often considered as the tumor suppressor gene by inducing pyroptosis in the cancer cells or by acting on the T lymphocytes through Granzyme B (Xia et al., 2019; De Schutter et al., 2021). Interestingly, ten percent of tumor cells have a higher GSDME expression level (De Schutter et al., 2021). Surprisingly, in this study, the high-risk groups showed increased GSDME levels than the low-risk group. This result demonstrated that the conventional tumor suppressor, GSDME, may act as an oncogene in the HCC microenvironment. It is reasonable to consider GSDME in promoting tumor development. The GSDME -induced pyroptosis can release many IL-1β and other inflammatory factors into the surrounding normal liver cells, inducing the normal cells into the tumor cells in inflammatory conditions (Zhou and Abbott, 2021). Coincidentally, the GSDME-mediated pyroptosis has been verified as the probable key point to respond to the toxic side effects of chemotherapeutic agents (Shen et al., 2021). Besides, some studies have demonstrated that the development of drugs targeting GSDME could be a promising treatment for HCC (Hu et al., 2019; Zhang et al., 2020; Liang et al., 2021; Shangguan et al., 2021). Hence, the function of GSDME in the HCC development and drug treatment requires further investigation. GPX4, a classical selenoprotein belonging to the glutathione peroxidase families, can reduce the membrane peroxidized phospholipids by transferring GSH (Liang et al., 2007), and it is considered to own unique lipid peroxidation inhibitory properties (Dabkowski et al., 2008). Subsequent studies have identified GPX4 as a classical negative regulator of ferroptosis and have roles in various cancers such as clear-cell carcinomas (CCCs), breast cancer, colon cancer (Yang et al., 2014; Seibt et al., 2019; Yang et al., 2020; Lee et al., 2021). Accordingly, recent studies indicated that GPX4 could also promote HCC development via inhibiting ferroptosis (Kim et al., 2020; Alves et al., 2021; Asperti et al., 2021; Chang et al., 2021). For non-cancer liver diseases, we aim to protect live cells from cell death, lipid peroxidation, and ROS release; induction of GPX4 expression was reported to be helpful to inhibit ferroptosis (Mao et al., 2020). However, for HCC, we aim at killing these cancer cells by inducing cell death, so it is beneficial to inhibit GPX4 and then activate ferroptosis (Jin et al., 2020). The latest study by Kang et al. indicated that GPX4 could suppress macrophagic pyroptosis in mice (Guerriero et al., 2015; Zhu et al., 2019). GPX4 knockout activated the lipid peroxidation-dependent caspase-11, which triggered the GSDMD cleavage to induce pyroptosis during polymicrobial sepsis. Here, GPX4 was found to be overexpressed in the HCC tissues, and the high-risk group showed an increased level of GPX4. Given the role of GPX4 in pyroptosis and ferroptosis, developing inhibitors targeting GPX4 might promote the GSDMD-mediated pyroptosis as well as ferroptosis, thus suppressing the survival of the HCC cells and drug resistance. SCAF11 has been previously reported to be involved in pyroptosis (Ye Y. et al., 2021). However, its role in cancer has not been explicitly studied so far. This study showed that the high expression of SCAF11 is related to the poor prognosis in the HCC patients, revealing that the inhibition of SCAF11 should be considered as a target to treat HCC.
Pyroptosis is always accompanied by inflammation and tumor immunity (Tan et al., 2021). However, chronic inflammation exerts essential functions in tumor initiation, progression, and invasion via suppressing the anti-tumor immune responses mediated by the immune cells such as the Natural Killer cells (NK) and M1 macrophages (Raposo et al., 2015). This study found decreased levels of immune cells such as the B cells, NK cells, TIL cells, and Mast cells, and immune pathways such as type II IFN response and cytolytic activity in the high-risk group, indicating that poor prognosis may result due to the decreased levels of anti-tumor immunity. Therefore, promoting anti-tumor immune responses are of great importance for effective clinical therapies.
Analysis of the NCI-60 cell line set in the CellMiner database indicated that the increased levels of the prognosis-related genes are positively correlated to drug resistance, such as Tamoxifen, Vinorelbine, AMG-900, and Litronesib. Notably, an increased level of GSDME is also positively correlated with the sensitivity to Lenvatinib, the first-line drug approved by FDA in 2018 for treating HCC (Yi et al., 2021). In 2019, a study showed that sorafenib, a classic FDA-approved drug for treating HCC, can induce macrophage pyroptosis and promote the NK cell responses against HCC (Hage et al., 2019). This research, combined with the correlation of Lenvatinib and GSDME, made us hypothesize that the mechanisms of Lenvatinib in HCC treatment might also be involved in inducing pyroptosis. We further confirmed this hypothesis by treating Huh7 cells with Lenvatinib, and the results showed that Lenvatinib upregulated the levels of both total GSDME and active GSDME N-terminal. Analyzing the correlation suggested that the pyroptosis-related prognostic genes are promising for anti-tumor drug development.
Conclusion
Our study revealed that developing HCC is inextricably linked to pyroptosis. Furthermore, the functional analysis, immune microenvironment, and drug correlation analysis established a basis for investigating the role of pyroptosis in the HCC development, determining the prognosis for the HCC patients, and providing clinical treatment. Based on the PRGs, an independent prognostic model was constructed for HCC, predicting the OS of the patients in the TCGA test cohort and ICGC external validation cohort. Moreover, our results confirmed the mRNA expression of three independent prognostic PRGs in vivo and in vitro. Our work will further assist in understanding the role of pyroptosis in HCC prognosis and drug sensitivity, thereby providing support for precision medicine.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Ethics Statement
The animal study was reviewed and approved by Institutional Animal Care and Use Committee of Westlake University.
Author Contributions
C-QS conceived and designed the study. X-WF performed bioinformatics analysis and experiments. C-QS and X-WF wrote the manuscript. All authors reviewed and approved the final manuscript.
Funding
This work is supported by grants from the National Natural Science Foundation of China (No. 92068103), Westlake Laboratory of Life Sciences and Biomedicine HRHI program (NO. W101266022101), and Westlake Education Foundation.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2021.748039/full#supplementary-material
References
Alves, A. d. F., Moura, A. C. d., Andreolla, H. F., Veiga, A. B. G. d., Fiegenbaum, M., Giovenardi, M., et al. (2021). Gene Expression Evaluation of Antioxidant Enzymes in Patients with Hepatocellular Carcinoma: RT-qPCR and Bioinformatic Analyses. Genet. Mol. Biol. 44, e20190373. doi:10.1590/1678-4685-GMB-2019-0373
Asperti, M., Bellini, S., Grillo, E., Gryzik, M., Cantamessa, L., Ronca, R., et al. (2021). H-ferritin Suppression and Pronounced Mitochondrial Respiration Make Hepatocellular Carcinoma Cells Sensitive to RSL3-Induced Ferroptosis. Free Radic. Biol. Med. 169, 294–303. doi:10.1016/j.freeradbiomed.2021.04.024
Bertheloot, D., Latz, E., and Franklin, B. S. (2021). Necroptosis, Pyroptosis and Apoptosis: an Intricate Game of Cell Death. Cell Mol Immunol. 18, 1106–1121. doi:10.1038/s41423-020-00630-3
Blasco, M. T., and Gomis, R. R. (2020). PD-L1 Controls Cancer Pyroptosis. Nat. Cel Biol. 22, 1157–1159. doi:10.1038/s41556-020-00582-w
Bruix, J., Reig, M., and Sherman, M. (2016). Evidence-Based Diagnosis, Staging, and Treatment of Patients with Hepatocellular Carcinoma. Gastroenterology 150, 835–853. doi:10.1053/j.gastro.2015.12.041
Cai, J., Chen, L., Zhang, Z., Zhang, X., Lu, X., Liu, W., et al. (2019). Genome-wide Mapping of 5-hydroxymethylcytosines in Circulating Cell-free DNA as a Non-invasive Approach for Early Detection of Hepatocellular Carcinoma. Gut 68, 2195–2205. doi:10.1136/gutjnl-2019-318882
Chang, W.-T., Bow, Y.-D., Fu, P.-J., Li, C.-Y., Wu, C.-Y., Chang, Y.-H., et al. (2021). A Marine Terpenoid, Heteronemin, Induces Both the Apoptosis and Ferroptosis of Hepatocellular Carcinoma Cells and Involves the ROS and MAPK Pathways. Oxid. Med. Cell Longev. 2021, 1–12. doi:10.1155/2021/7689045
Chu, Q., Jiang, Y., Zhang, W., Xu, C., Du, W., Tuguzbaeva, G., et al. (2016). Pyroptosis Is Involved in the Pathogenesis of Human Hepatocellular Carcinoma. Oncotarget 7, 84658–84665. doi:10.18632/oncotarget.12384
Craig, A. J., von Felden, J., Garcia-Lezana, T., Sarcognato, S., and Villanueva, A. (2020). Tumour Evolution in Hepatocellular Carcinoma. Nat. Rev. Gastroenterol. Hepatol. 17, 139–152. doi:10.1038/s41575-019-0229-4
Dabkowski, E. R., Williamson, C. L., and Hollander, J. M. (2008). Mitochondria-specific Transgenic Overexpression of Phospholipid Hydroperoxide Glutathione Peroxidase (GPx4) Attenuates Ischemia/Reperfusion-Associated Cardiac Dysfunction. Free Radic. Biol. Med. 45, 855–865. doi:10.1016/j.freeradbiomed.2008.06.021
Dai, Y., Qiang, W., Lin, K., Gui, Y., Lan, X., and Wang, D. (2021). An Immune-Related Gene Signature for Predicting Survival and Immunotherapy Efficacy in Hepatocellular Carcinoma. Cancer Immunol. Immunother. 70, 967–979. doi:10.1007/s00262-020-02743-0
De Schutter, E., Croes, L., Ibrahim, J., Pauwels, P., Op de Beeck, K., Vandenabeele, P., et al. (2021). GSDME and its Role in Cancer: From behind the Scenes to the Front of the Stage. Int. J. Cancer 148, 2872–2883. doi:10.1002/ijc.33390
Du, H., Xie, S., Guo, W., Che, J., Zhu, L., Hang, J., et al. (2021). Development and Validation of an Autophagy-Related Prognostic Signature in Esophageal Cancer. Ann. Transl Med. 9, 317. doi:10.21037/atm-20-4541
Fang, Y., Huang, S., Han, L., Wang, S., and Xiong, B. (2021). Comprehensive Analysis of Peritoneal Metastasis Sequencing Data to Identify LINC00924 as a Prognostic Biomarker in Gastric Cancer. Cancer Manag. Res. 13, 5599–5611. doi:10.2147/CMAR.S318704
Fink, S. L., and Cookson, B. T. (2005). Apoptosis, Pyroptosis, and Necrosis: Mechanistic Description of Dead and Dying Eukaryotic Cells. Infect. Immun. 73, 1907–1916. doi:10.1128/IAI.73.4.1907-1916.2005
Galluzzi, L., Vitale, I., Aaronson, S. A., Abrams, J. M., Adam, D., Agostinis, P., et al. (2018). Molecular Mechanisms of Cell Death: Recommendations of the Nomenclature Committee on Cell Death 2018. Cell Death Differ. 25, 486–541. doi:10.1038/s41418-017-0012-4
Guerriero, E., Capone, F., Accardo, M., Sorice, A., Costantini, M., Colonna, G., et al. (2015). GPX4 and GPX7 Over-expression in Human Hepatocellular Carcinoma Tissues. Eur. J. Histochem. 59, 2540. doi:10.4081/ejh.2015.2540
Hage, C., Hoves, S., Strauss, L., Bissinger, S., Prinz, Y., Pöschinger, T., et al. (2019). Sorafenib Induces Pyroptosis in Macrophages and Triggers Natural Killer Cell-Mediated Cytotoxicity against Hepatocellular Carcinoma. Hepatology 70, 1280–1297. doi:10.1002/hep.30666
Hu, J., Dong, Y., Ding, L., Dong, Y., Wu, Z., Wang, W., et al. (2019). Local Delivery of Arsenic Trioxide Nanoparticles for Hepatocellular Carcinoma Treatment. Sig Transduct Target. Ther. 4, 28. doi:10.1038/s41392-019-0062-9
Humphries, F., Shmuel-Galia, L., Ketelut-Carneiro, N., Li, S., Wang, B., Nemmara, V. V., et al. (2020). Succination Inactivates Gasdermin D and Blocks Pyroptosis. Science 369, 1633–1637. doi:10.1126/science.abb9818
Jin, M., Shi, C., Li, T., Wu, Y., Hu, C., and Huang, G. (2020). Solasonine Promotes Ferroptosis of Hepatoma Carcinoma Cells via Glutathione Peroxidase 4-induced Destruction of the Glutathione Redox System. Biomed. Pharmacother. 129, 110282. doi:10.1016/j.biopha.2020.110282
Jin, H., Shi, Y., Lv, Y., Yuan, S., Ramirez, C. F. A., Lieftink, C., et al. (2021). EGFR Activation Limits the Response of Liver Cancer to Lenvatinib. Nature 595, 730–734. doi:10.1038/s41586-021-03741-7
Ju, X., Yang, Z., Zhang, H., and Wang, Q. (2021). Role of Pyroptosis in Cancer Cells and Clinical Applications. Biochimie 185, 78–86. doi:10.1016/j.biochi.2021.03.007
Karki, R., and Kanneganti, T.-D. (2019). Diverging Inflammasome Signals in Tumorigenesis and Potential Targeting. Nat. Rev. Cancer 19, 197–214. doi:10.1038/s41568-019-0123-y
Kim, D. H., Kim, W. D., Kim, S. K., Moon, D. H., and Lee, S. J. (2020). TGF-β1-mediated Repression of SLC7A11 Drives Vulnerability to GPX4 Inhibition in Hepatocellular Carcinoma Cells. Cell Death Dis. 11, 406. doi:10.1038/s41419-020-2618-6
Kong, H., Wang, Y., Zeng, X., Wang, Z., Wang, H., and Xie, W. (2015). Differential Expression of Inflammasomes in Lung Cancer Cell Lines and Tissues. Tumor Biol. 36, 7501–7513. doi:10.1007/s13277-015-3473-4
Kovacs, S. B., and Miao, E. A. (2017). Gasdermins: Effectors of Pyroptosis. Trends Cel Biol. 27, 673–684. doi:10.1016/j.tcb.2017.05.005
Lee, N., Carlisle, A. E., Peppers, A., Park, S. J., Doshi, M. B., Spears, M. E., et al. (2021). xCT-Driven Expression of GPX4 Determines Sensitivity of Breast Cancer Cells to Ferroptosis Inducers. Antioxidants 10, 317. doi:10.3390/antiox10020317
Liang, H., Remmen, H. V., Frohlich, V., Lechleiter, J., Richardson, A., and Ran, Q. (2007). Gpx4 Protects Mitochondrial ATP Generation against Oxidative Damage. Biochem. Biophys. Res. Commun. 356, 893–898. doi:10.1016/j.bbrc.2007.03.045
Liang, W.-F., Gong, Y.-X., Li, H.-F., Sun, F.-L., Li, W.-L., Chen, D.-Q., et al. (2021). Curcumin Activates ROS Signaling to Promote Pyroptosis in Hepatocellular Carcinoma HepG2 Cells. In Vivo 35, 249–257. doi:10.21873/invivo.12253
Liu, P.-H., Hsu, C.-Y., Hsia, C.-Y., Lee, Y.-H., Su, C.-W., Huang, Y.-H., et al. (2016). Prognosis of Hepatocellular Carcinoma: Assessment of Eleven Staging Systems. J. Hepatol. 64, 601–608. doi:10.1016/j.jhep.2015.10.029
Llovet, J. M., Zucman-Rossi, J., Pikarsky, E., Sangro, B., Schwartz, M., Sherman, M., et al. (2016). Hepatocellular Carcinoma. Nat. Rev. Dis. Primers 2, 16018. doi:10.1038/nrdp.2016.18
Llovet, J. M., Kelley, R. K., Villanueva, A., Singal, A. G., Pikarsky, E., Roayaie, S., et al. (2021). Hepatocellular Carcinoma. Nat. Rev. Dis. Primers 7, 6. doi:10.1038/s41572-020-00240-3
Man, S. M., and Kanneganti, T.-D. (2015). Regulation of Inflammasome Activation. Immunol. Rev. 265, 6–21. doi:10.1111/imr.12296
Mao, L., Zhao, T., Song, Y., Lin, L., Fan, X., Cui, B., et al. (2020). The Emerging Role of Ferroptosis in Non-cancer Liver Diseases: Hype or Increasing hope? Cel Death Dis. 11, 518. doi:10.1038/s41419-020-2732-5
Raposo, T. P., Beirão, B. C. B., Pang, L. Y., Queiroga, F. L., and Argyle, D. J. (2015). Inflammation and Cancer: till Death Tears Them Apart. Vet. J. 205, 161–174. doi:10.1016/j.tvjl.2015.04.015
Rathinam, V. A. K., and Fitzgerald, K. A. (2016). Inflammasome Complexes: Emerging Mechanisms and Effector Functions. Cell 165, 792–800. doi:10.1016/j.cell.2016.03.046
Seibt, T. M., Proneth, B., and Conrad, M. (2019). Role of GPX4 in Ferroptosis and its Pharmacological Implication. Free Radic. Biol. Med. 133, 144–152. doi:10.1016/j.freeradbiomed.2018.09.014
Shangguan, F., Zhou, H., Ma, N., Wu, S., Huang, H., Jin, G., et al. (2021). A Novel Mechanism of Cannabidiol in Suppressing Hepatocellular Carcinoma by Inducing GSDME Dependent Pyroptosis. Front. Cel Dev. Biol. 9, 697832. doi:10.3389/fcell.2021.697832
Shankavaram, U. T., Reinhold, W. C., Nishizuka, S., Major, S., Morita, D., Chary, K. K., et al. (2007). Transcript and Protein Expression Profiles of the NCI-60 Cancer Cell Panel: an Integromic Microarray Study. Mol. Cancer Ther. 6, 820–832. doi:10.1158/1535-7163.MCT-06-0650
Shankavaram, U. T., Varma, S., Kane, D., Sunshine, M., Chary, K. K., Reinhold, W. C., et al. (2009). CellMiner: a Relational Database and Query Tool for the NCI-60 Cancer Cell Lines. BMC Genomics 10, 277. doi:10.1186/1471-2164-10-277
Shen, X., Wang, H., Weng, C., Jiang, H., and Chen, J. (2021). Caspase 3/GSDME-dependent Pyroptosis Contributes to Chemotherapy Drug-Induced Nephrotoxicity. Cel Death Dis. 12, 186. doi:10.1038/s41419-021-03458-5
Shojaie, L., Iorga, A., and Dara, L. (2020). Cell Death in Liver Diseases: A Review. Int. J. Mol. Sci. 21, 9682. doi:10.3390/ijms21249682
Song, C.-Q., Li, Y., Mou, H., Moore, J., Park, A., Pomyen, Y., et al. (2017). Genome-Wide CRISPR Screen Identifies Regulators of Mitogen-Activated Protein Kinase as Suppressors of Liver Tumors in Mice. Gastroenterology, 152:1161, 1173. doi:10.1053/j.gastro.2016.12.002
Sun, S., Guo, W., Lv, F., Zhang, G., Wang, J., Li, R., et al. (2021). Comprehensive Analysis of Ferroptosis Regulators in Lung Adenocarcinomas Identifies Prognostic and Immunotherapy-Related Biomarkers. Front. Mol. Biosci. 8, 587436. doi:10.3389/fmolb.2021.587436
Tan, Y., Sun, R., Liu, L., Yang, D., Xiang, Q., Li, L., et al. (2021). Tumor Suppressor DRD2 Facilitates M1 Macrophages and Restricts NF-κB Signaling to Trigger Pyroptosis in Breast Cancer. Theranostics 11, 5214–5231. doi:10.7150/thno.58322
Tang, R., Xu, J., Zhang, B., Liu, J., Liang, C., Hua, J., et al. (2020). Ferroptosis, Necroptosis, and Pyroptosis in Anticancer Immunity. J. Hematol. Oncol. 13, 110. doi:10.1186/s13045-020-00946-7
Van Laer, L., Vrijens, K., Thys, S., Van Tendeloo, V. F., Smith, R. J., Van Bockstaele, D. R., et al. (2004). DFNA5: Hearing Impairment Exon Instead of Hearing Impairment Gene? J. Med. Genet. 41, 401–406. doi:10.1136/jmg.2003.015073
Wang, Y., Gao, W., Shi, X., Ding, J., Liu, W., He, H., et al. (2017). Chemotherapy Drugs Induce Pyroptosis through Caspase-3 Cleavage of a Gasdermin. Nature 547, 99–103. doi:10.1038/nature22393
Wu, D., Wang, S., Yu, G., and Chen, X. (2021). Cell Death Mediated by the Pyroptosis Pathway with the Aid of Nanotechnology: Prospects for Cancer Therapy. Angew. Chem. Int. Ed. 60, 8018–8034. doi:10.1002/anie.202010281
Xia, X., Wang, X., Cheng, Z., Qin, W., Lei, L., Jiang, J., et al. (2019). The Role of Pyroptosis in Cancer: Pro-Cancer or Pro-"Host"? Cel Death Dis. 10, 650. doi:10.1038/s41419-019-1883-8
Yan, J., Cao, J., and Chen, Z. (2021). Mining Prognostic Markers of Asian Hepatocellular Carcinoma Patients Based on the Apoptosis-Related Genes. BMC Cancer 21, 175. doi:10.1186/s12885-021-07886-6
Yang, W. S., SriRamaratnam, R., Welsch, M. E., Shimada, K., Skouta, R., Viswanathan, V. S., et al. (2014). Regulation of Ferroptotic Cancer Cell Death by GPX4. Cell 156, 317–331. doi:10.1016/j.cell.2013.12.010
Yang, J. D., Hainaut, P., Gores, G. J., Amadou, A., Plymoth, A., and Roberts, L. R. (2019). A Global View of Hepatocellular Carcinoma: Trends, Risk, Prevention and Management. Nat. Rev. Gastroenterol. Hepatol. 16, 589–604. doi:10.1038/s41575-019-0186-y
Yang, Y., Yang, L., Jiang, S., Yang, T., Lan, J., Lei, Y., et al. (2020). HMGB1 Mediates Lipopolysaccharide-Induced Inflammation via Interacting with GPX4 in colon Cancer Cells. Cancer Cel Int. 20, 205. doi:10.1186/s12935-020-01289-6
Yang, Y., Wu, G., Li, Q., Zheng, Y., Liu, M., Zhou, L., et al. (2021). Angiogenesis-Related Immune Signatures Correlate with Prognosis, Tumor Microenvironment, and Therapeutic Sensitivity in Hepatocellular Carcinoma. Front. Mol. Biosci. 8, 690206. doi:10.3389/fmolb.2021.690206
Ye, J., She, X., Liu, Z., He, Z., Gao, X., Lu, L., et al. (2021a). Eukaryotic Initiation Factor 4A-3: A Review of its Physiological Role and Involvement in Oncogenesis. Front. Oncol. 11, 712045. doi:10.3389/fonc.2021.712045
Ye, Y., Dai, Q., and Qi, H. (2021b). A Novel Defined Pyroptosis-Related Gene Signature for Predicting the Prognosis of Ovarian Cancer. Cell Death Discov. 7, 71. doi:10.1038/s41420-021-00451-x
Yi, C., Chen, L., Lin, Z., Liu, L., Shao, W., Zhang, R., et al. (2021). Lenvatinib Targets FGF Receptor 4 to Enhance Antitumor Immune Response of Anti-programmed Cell Death‐1 in HCC. Hepatology 74, 2544–2560. doi:10.1002/hep.31921
Yu, J., Li, S., Qi, J., Chen, Z., Wu, Y., Guo, J., et al. (2019). Cleavage of GSDME by Caspase-3 Determines Lobaplatin-Induced Pyroptosis in colon Cancer Cells. Cel Death Dis 10, 193. doi:10.1038/s41419-019-1441-4
Yuan, C., Yuan, M., Chen, M., Ouyang, J., Tan, W., Dai, F., et al. (2021a). Prognostic Implication of a Novel Metabolism-Related Gene Signature in Hepatocellular Carcinoma. Front. Oncol. 11, 666199. doi:10.3389/fonc.2021.666199
Yuan, M., Wang, Y., Sun, Q., Liu, S., Xian, S., Dai, F., et al. (2021b). Identification of a Nine Immune-Related lncRNA Signature as a Novel Diagnostic Biomarker for Hepatocellular Carcinoma. Biomed. Res. Int. 2021, 1–10. doi:10.1155/2021/9798231
Zhang, X., Zhang, P., An, L., Sun, N., Peng, L., Tang, W., et al. (2020). Miltirone Induces Cell Death in Hepatocellular Carcinoma Cell through GSDME-dependent Pyroptosis. Acta Pharm. Sin. B 10, 1397–1413. doi:10.1016/j.apsb.2020.06.015
Zhao, W.-j., Ou, G.-y., and Lin, W.-w. (2021). Integrative Analysis of Neuregulin Family Members-Related Tumor Microenvironment for Predicting the Prognosis in Gliomas. Front. Immunol. 12, 682415. doi:10.3389/fimmu.2021.682415
Zhou, B., and Abbott, D. W. (2021). Gasdermin E Permits Interleukin-1 Beta Release in Distinct Sublytic and Pyroptotic Phases. Cel Rep. 35, 108998. doi:10.1016/j.celrep.2021.108998
Keywords: pyroptosis, prognosis, biomarker, immune infiltration, drug sensitivity, hepatocellular carcinoma
Citation: Fu X-W and Song C-Q (2021) Identification and Validation of Pyroptosis-Related Gene Signature to Predict Prognosis and Reveal Immune Infiltration in Hepatocellular Carcinoma. Front. Cell Dev. Biol. 9:748039. doi: 10.3389/fcell.2021.748039
Received: 27 July 2021; Accepted: 22 October 2021;
Published: 08 November 2021.
Edited by:
Fan Feng, The 302nd Hospital of PLA, ChinaReviewed by:
QIsi Lu, Southern Medical University, ChinaTao Huang, Sun Yat-sen University Cancer Center (SYSUCC), China
Copyright © 2021 Fu and Song. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Chun-Qing Song, songchunqing@westlake.edu.cn