Skip to main content

A necroptosis-related gene signature to predict prognosis and immune features in hepatocellular carcinoma

Abstract

Background and Aim

Necroptosis plays an important role in hepatocellular carcinoma (HCC) development, recurrence, and immunotherapy tolerance. We aimed to build a new prognostic necroptosis-related gene signature that could be used for survival and immunotherapy prediction in HCC patients.

Methods

We found that necroptosis was associated with HCC progression and survival outcomes and was involved in the immune infiltration of HCC. Multiple bioinformatics methods including WGCNA, LASSO-Cox regression, stepwise Cox regression, and Random Forest and Boruta model analysis, were used to establish a prognostic profile related to necroptosis. The necroptosis-related gene signature was validated in ICGC and GSE14520 datasets.

Results

This five-gene signature showed excellent predictive performance and was an independent risk factor for patients’ overall survival outcome in the three cohorts. Moreover, this signature was an exact predictor using fewer genes than previous gene signatures. Finally, qRT-PCR and immunohistochemical staining investigations were performed in previously collected fresh frozen tumor tissues from HCC patients and their paracancerous normal tissues, and the results were consistent with the bioinformatics results. We found that LGALS3 not only affected the proliferation and migration ability of HepG2 cells but also affected necroptosis and the expression of inflammatory cytokines.

Conclusion

In summary, we established and validated an individualized prognostic profile related to necroptosis to forecast the therapeutic response to immune therapy, which might offer a potential non-apoptotic therapeutic target for HCC patients.

Peer Review reports

Introduction

Necroptosis is a novel form of lytic cell death distinct from apoptosis, and the immunostimulatory molecules released during its activation, including inflammatory cytokines and damage-associated molecular patterns (DAMP), can spark an intense inflammatory response and trigger the death of infected or damaged cells in vivo [1]. There is growing evidence that multiple stimuli including TNFα/TNFR signaling, hypoxia, radiation, and chemotherapeutic agents can trigger the necroptosis pathway [2, 3]. Although the intracellular signaling pathways involved in the regulation of apoptosis and necroptosis vary widely; for example, MLKL (mixed lineage kinase domain-like pseudokinase) and RIPKs (receptor-interacting protein kinases) are thought to be key regulators of the necroptosis process, whereas caspases, such as caspase-8, are key mediators of intracellular apoptosis, crosstalk between apoptosis and necroptosis was found during signal transduction [4]. In addition to its ability to fight pathogenic infections and induce strong inflammatory responses, necroptosis is involved in the regulation of cell death and neuroinflammation during the pathology of several neurodegenerative diseases, including Alzheimer’s disease and Parkinson’s disease [5]. Necroptosis has also been found to be involved in tumor progression as a two-sided joker in various types of cancers: on the one hand, key mediators of the necroptosis pathway promote migration and immunosuppression; on the other hand, necroptosis plays a tumor-suppressive role in cancer when apoptosis is impaired [6,7,8]. The above suggests that further exploration of the mechanisms and functions of the necroptosis pathway is necessary given the critical role of necroptosis in cancer progression and migration as well as in cancer immunosurveillance mechanisms, which will contribute to the development of a new therapeutic approach to eliminate apoptosis-resistant cancer cells for cancer patients.

Previous studies have reported that increased necroptosis is significantly related to nonalcoholic fatty liver disease (NAFLD) in humans [9]. Necroptosis-mediated inflammatory cytokines released in the hepatic microenvironment accelerate liver fibrosis in combination with oxidative stress in a mouse model [10]. The epigenetic hepatic microenvironment is closely related to liver tumorigenesis and hepatocytes containing necroptosis-mediated inflammatory cytokines such as Cxcl13, Ccl6, and Ccl8, which can induce hepatocellular carcinoma (HCC) when surrounded by apoptotic hepatocytes in mice [11]. Furthermore, decreased RIPK3 may enhance the accumulation and polarization of M2 tumor-associated macrophages (TAMs) and participate in the reprogramming of TAM lipid metabolism [12]. We, therefore, speculate that increased hepatic necroptosis plays an important role in the prognosis and tumor immune microenvironment (TME) for HCC, and a deeper insight into it can help clinicians explore new targets for tumor immunotherapy. In this study, we confirmed the speculation of a strong link between necroptosis and poor prognosis and immune cell infiltration in HCC patients, and subsequently established and validated an individualized prognostic profile associated with the necroptosis pathway using multiple bioinformatics methods to forecast the therapeutic response to immune checkpoint blockade (ICB) therapy.

Materials and methods

Acquisition of public datasets and characterization of Z-scores associated with necroptosis

Transcriptome expression data from TCGA-LIHC, ICGC (LIRI-JP), and GSE14520 datasets were included in this study, and patients’ general information is listed in Table S1. The TCGA-LIHC and ICGC (LIRI-JP) datasets were obtained based on RNA-seq analysis, and we performed log2 (FPKM + 1) transformation to normalize the RNA-seq data. The GSE14520 dataset was obtained based on Affymetrix Human Genome U133A Array, and we used the ‘LIMMA’ package’s normalizeBetweenArrays function to normalize the data. Samples in the three datasets without complete survival data were not included in this study. Genes participating in the necroptosis pathway were identified from the Molecular Signatures Database (https://www.gsea-msigdb.org/gsea/msigdb) and quantified in each HCC sample in the TCGA using a single sample genomic enrichment analysis (ssGSEA) method based on transcriptome analysis data, with Z-score scaling applied to the ssGSEA scores.

Functional annotation

To illustrate the functional annotations implicated with the necroptosis-associated genes, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis and gene ontology (GO) annotation analysis were performed using the “clusterProfiler” package in R software [13,14,15]. The significance level (P-value) and false-positive rate (FDR) of each signal pathway were calculated using Fisher exact test and multiple comparison test [16].

Hub gene selection and gene signature construction

Based on transcriptome expression data and necroptosis Z-scores in TCGA, weighted gene coexpression network analysis (WGCNA) was utilized to build a scale-free coexpression network to determine the gene module most important to necroptosis [17]. Moreover, gene significance (GS) was used to quantify the connections between individual genes and necroptosis Z-scores, while module members reflected the correlations between module characteristic genes and gene expression patterns. Genes identified from the module that was most associated with necroptosis Z-scores were evaluated as candidates using the p value threshold of GS < 0.0001 and the significance level of univariate Cox regression of p < 0.01. The least absolute shrinkage and selection operator (LASSO) Cox regression, stepwise Cox regression, and the Random Forest and Boruta (RFB) model were then used to filter the most reliable prognostic candidates. The multivariate Cox relapse coefficient (β) was utilized to produce a risk score based on the notion of combining the following equation directly with the amount of overlapping gene expression. Risk score=∑iCoefficient (mRNAi)*Expression (mRNAi). We finally separated the patients in the three available datasets into two groups based on the optimum hazard score edge. The predictive control and autonomy of the prognostic signature were assessed using ROC analysis, Kaplan‒Meier survival analysis, and Cox relative risk relapse study.

Tumor-infiltrating immune cells, genetic alterations, and immune checkpoint gene analysis

Multiple approaches were used to determine the abundance ratios of tumor-infiltrating immune cells (TIICs) in the HCC TME, including CIBERSORT [18], TIMER [19], QUANTISEQ [20], and xCELL [21]. The mutation and CNA data of HCC patients were retrieved from TCGA and the changes in genetic variants across various subgroups were examined using the R package “maftools.“ To determine how the prognostic classifier influences immunotherapies, we looked at the relationship between necroptosis Z-scores, risk scores, and three potentially available targeted immune checkpoint genes, PD1, PD-L1, and CTLA4.

Drug susceptibility analysis

With 574 advanced clinical trials and 216 FDA-approved drugs used for follow-up research, the association between anticancer drug sensitivity and molecules in the model was directly explored in the CellMiner database [22]. The cut-off criteria for tumor-sensitive drugs were an adjusted P value of < 0.001 and a Pearson correlation coefficient of > 0.4.

Quantitative real-time polymerase chain reaction (qRT‒PCR) and immunohistochemistry staining assay in clinical samples

Twenty previously collected HCC patients’ fresh frozen tumor tissues and their paracancerous normal tissues were chosen as clinical samples for qRT‒PCR investigation. Table S2 contains the primer sequences. Single-blind and unified criteria methods were used by two experienced pathologists to assess the IHC results. A final score of the sum of the expression score (no positive cells = 0, < 10% = 1, 10–50% = 2, positive staining of > 50% = 3) and the intensity score (negative = 0, weak = 1, moderate = 2, strong = 3) was used to distinguish between low (≤ 4) and high (> 4) expression of LGALS3 in HCC and normal tissues.

Cell culture and siRNA transfection

HepG2 cells were purchased from the American Type Culture Collection and cultured in recommended DMEM (Sangon Biotech, China) with 10% foetal bovine serum (FBS, Sangon Biotech, China) in 100% humidity at 37 °C with 5% CO2. Small interfering RNA (siRNA) (Sense: 5′-GCC ACU GAU UGU GCC UUA UTT-3′; Antisense:5′-AUA AGG CAC AAU CAG UGG CTT-3′ ) and negative control (Sense: 5′- UUC UCC GAA CGU GUC ACG UTT-3′; Antisense: 5′- ACG UGA CAC GUU CGG AGA ATT-3′) were selected to reduce LGALS3 expression in Hep3B cells using lipofectamine™ 3000 transfection reagent (Invitrogen, Carlsbad, USA) according to the manufacturer’s instructions.

Cell proliferation and migration assay

According to the manufacturer’s instructions, the vitality of cells was determined using the cell counting kit-8 (CCK-8, Sangon Biotech, Shanghai, China). In the CCK-8 experiment, cells were seeded at a density of 1 × 105 cells per well in 96-well cell culture clusters and cultivated for 1 day, 2 days, and 3 days. After culture, 10 µL of CCK-8 solution was applied to each well, and the absorbance was measured with a microplate reader at 450 nm within 4 h. Cells were seeded in 6-well plates and cultivated to approximately 80% confluence in serum-free media before being scratched with a sterile pipette tip for the wound-healing experiment. After removing cell debris with PBS washing, cells were grown in DMEM with 10% FBS for the following 48 h. After photographing the wound width of the cell monolayers, the area of the wound width was measured. This test required three independent duplicates to assure accuracy, and the wound closure rate was determined from images as [1 - (wound area/original wound area)]. After covering Transwell filters with Matrigel, cultured cells were resuspended in 200 µL of serum-free DMEM at a density of 1 × 104 cells per mL and plated into Transwell inserts, and the wells were filled with 500 L of DMEM supplemented with 10% FBS for the Transwell assay. After being cultured at 37 °C for 48 h, cells adhered to the bottom of the Transwell filters were stained with 0.1% crystal violet in PBS for 15 min and counted under microscopy at 200 x magnification.

Western blotting

HepG2 cells were harvested in RIPA Lysis Buffer, and were lysed by ultrasound treatment (#DH92-IIN; LAWSON). Supernatants were collected and protein concentrations were determined using Coomassie brilliant blue G250 stain (#C8420, Solarbio). Proteins were separated by SDS-PAGE and transferred to a Nylon membrane. Membranes were blocked in 5% Milk/TBS-T for 2 h at room temperature and incubated overnight at 4 °C with primary antibodies and subsequently with HRP-conjugated secondary antibody [16]. Antibodies against MLKL (#26,539), RIPK1(#3493), RIPK3(#10,188) and β-actin (#93,473) were purchased from CST. Necrostatin-1(Nec-1) (HY-15,760), a necrotic apoptosis inhibitor was purchased from MCE.

Statistical analysis

Quantitative variables were evaluated using an independent-samples t-test. ROC curve analysis and Kaplan‒Meier survival analysis were utilized to examine the prediction performance of survival outcomes using R software (version 4.0.3). A Cox proportional model was used to explore the relationship between a prognostic classifier and survival outcomes, as well as other clinical factors. The results were considered statistically significant when the P value was less than 0.05.

Results

Necroptosis was associated with HCC progression and survival outcomes

A total of 159 necroptosis-associated genes were extracted from prior research and quantified in each sample using ssGSEA analysis in TCGA, as shown in Fig. 1A. HCC patients with a low necroptosis Z-score showed a greater recurrence incidence (Fig. 1B) and more advanced clinical stages during follow-up (Fig. 1C). Furthermore, as demonstrated in Fig. 1D, patients with a greater necroptosis Z-score had considerably better survival outcomes.

Fig. 1
figure 1

Necroptosis was associated with HCC progression and survival outcomes. (A) Overview of the Z-scores associated with necroptosis in TCGA. Differences in the percentage of tumor recurrence (B) and advanced clinical stages (C) between different Z-score groups. (D) Kaplan-Meier survival plot of different Z-scores

Necroptosis was involved in the immune infiltration of HCC

Stromal and immune fractions were estimated by ESTIMATE to assess the abundance of stromal and immune cells within the tumor [23], and the results showed that HCC samples with high necroptosis Z-score defined as having ‘immune hot’ tumors exhibited higher stromal, immune and ESTIMATE scores (Fig. 2A). Then, a heatmap was plotted to assess the relationship between the necroptosis Z-score and the level of immune cell infiltration, and substantial differences were found. We summarized the composition of the 42 immune cells infiltrating in HCC patients from both high and low necroptosis Z-score groups. The results indicated that high necroptosis Z-score group had a distinct pattern of immune infiltration compared to low necroptosis Z-score group. Furthermore, high necroptosis Z-score group had a significantly higher abundance of immune cells—including memory CD8 T cells, memory B cells, CD4 T cells, TAMs, than low necroptosis Z-score group, but had a significantly lower abundance of immune cells—including resting NK cells, mast cells, than low necroptosis Z-score group (Fig. 2B). Finally, we found that the necroptosis Z-score was significantly related to the expression of the three immune checkpoint genes (Figure S1A), and patients with higher necroptosis Z-scores had remarkably lower expression of PD1, PD-L1, and CTLA4 (Figure S1B).

Fig. 2
figure 2

Necroptosis was involved in the immune infiltration of HCC. (A) Differences in stromal, immunes, and estimate scores between different Z-score groups. (B) Heatmap of levels of immune cell infiltration in different Z-score groups

Necroptosis-related gene signature construction in TCGA

Following the removal of outliers (Figure S2A), the top 5000 ranked genes with the greatest variation were chosen for WGCNA. A soft threshold power of 9 (Figure S2B) was used to assure the reasonableness of the scale-free network (scale-free R2 = 0.9). Four nongray modules were created, as illustrated in Fig. 3A. The green module, which had 639 genes, was the most closely associated with necroptosis of the four modules (Fig. 3B&C). The biological importance of these genes in the green module was investigated using enrichment analysis, and the results revealed that they were mostly enriched in immune activation and response (Figure S3). After that, a univariate Cox regression analysis with P values less than 0.01 was used to look for predictive genes, yielding a total of 81 candidate genes identified (Fig. 4A). The LASSO-Cox regression model was used to further filter these candidate genes (Fig. 4B), and 15 genes were ruled out (Fig. 4C). Similarly, the RFB model was used to filter these candidate genes, and the top 30 most important predictors are shown in Fig. 4D. Eighteen candidate genes were filtered by stepwise Cox proportional analysis. Five overlapping candidate genes were discovered (Fig. 4E) and then incorporated into a multivariate Cox regression model to create a predictive five-gene signature. Risk score = HMOX1*0.3318815 + VNN2*0.1493589 + TNFRSF4*0.4052564 - KLRB1*0.6610017 + LGALS3*0.2188195.

Fig. 3
figure 3

WGCNA for hub genes. (A) Four non-gray modules were created in the WGCNA network. (B-C) The green module was the most closely associated with necroptosis of the four modules

Fig. 4
figure 4

Necroptosis-related gene signature construction in TCGA. (A) A total of 81 candidate genes were identified from the green module. (B) Genes screened by the LASSO-Cox regression model. (C) 15 genes were ruled out. (D) The top 30 most important predictors are identified by the RFB model. (E) Five overlapping candidate genes were discovered

The necroptosis-related gene signature was an independent risk factor for patients’ overall survival outcome in TCGA

Using an adequate risk score threshold, patients were split into high- and low-risk categories (Fig. 5A). During follow-up, we discovered that HCC patients with high risk scores had higher mortality (Fig. 5B), and recurrence rates (Fig. 5C) and more advanced clinical stages (Fig. 5D-F). Furthermore, patients with lower risk scores had significantly improved survival results, as shown in Fig. 5G. Additionally, ROC analysis revealed that this necroptosis-related gene signature showed excellent predictive performance, with AUCs of 0.775, 0.745, and 0.777 at one, three, and five years, respectively (Fig. 5H). Finally, after adjusting for other clinical features, multivariate Cox regression analysis revealed that this signature could serve as an independent prognostic factor for HCC patients (HR = 2.278, 95% CI 1.590–3.263, P < 0.001).

Fig. 5
figure 5

Necroptosis-related gene signature was an independent risk factor for patients’ overall survival outcome in TCGA. (A) Distribution of risk scores, OS status, and gene expression profiles. HCC patients with high-risk scores had higher mortality (B), recurrence rate (C), and more advanced clinical stages (D-F). (G) Patients with lower risk scores had significantly improved survival results. (H) ROC analysis revealed that this necroptosis-related gene signature showed excellent predictive performance

The necroptosis-related gene signature was also an independent risk factor for patients’ overall survival outcome in the ICGC and GSE14520 datasets

To validate the signature, the ICGC and GSE14520 datasets were utilized as validation cohorts. In the ICGC cohort, patients were split into high- and low-risk categories by using an adequate risk score threshold (Fig. 6A). During follow-up, we discovered that HCC patients with high-risk scores had higher mortality (Fig. 6B) and more advanced clinical stages (Fig. 6C). As shown in Fig. 6D, patients with lower risk scores had significantly improved survival results. Moreover, ROC analysis revealed that this necroptosis-related gene signature showed good predictive performance (Fig. 6E). Similarly, patients in the GSE14520 dataset were split in to high- and low-risk categories by using an adequate risk score threshold (Fig. 7A). HCC patients with high-risk scores had higher mortality (Fig. 7B) and more advanced clinical stages (Fig. 7C) during follow-up. As seen in Fig. 7D, patients with higher risk scores had significantly poorer survival results. Additionally, ROC analysis showed that this necroptosis-associated genetic signature was an exact predictor of the overall survival rate (Fig. 7E). Most importantly, multivariate Cox regression analysis revealed that this signature could also serve as an independent prognostic factor for HCC patients in ICGC (HR = 2.278, 95% CI 1.590–3.263, P < 0.001) and GSE14520 (HR = 2.278, 95% CI 1.590–3.263, P < 0.001) after adjusting for other clinical features.

Fig. 6
figure 6

Necroptosis-related gene signature was an independent risk factor for patients’ overall survival outcome in the ICGC dataset. (A) Distribution of risk scores, OS status, and gene expression profiles. HCC patients with high-risk scores had higher mortality (B) and more advanced clinical stages (C). (D) Patients with lower risk scores had significantly improved survival results. (E) ROC analysis revealed that this necroptosis-related gene signature showed excellent predictive performance

Fig. 7
figure 7

Necroptosis-related gene signature was an independent risk factor for patients’ overall survival outcome in the GSE14520 dataset. (A) Distribution of risk scores, OS status, and gene expression profiles. HCC patients with high-risk scores had higher mortality (B) and more advanced clinical stages (C). (D) Patients with lower risk scores had significantly improved survival results. (E) ROC analysis revealed that this necroptosis-related gene signature showed excellent predictive performance

Genetic alterations and immune checkpoint gene analysis

As shown in Figure S4, the mutation rates of the top 10 most substantially altered genes were significantly different between the high- and low-risk score subgroups, according to the results of genetic alteration analysis. Furthermore, we found that risk scores were significantly related to the expression of the three immune checkpoint genes (Figure S5A) and patients with higher risk scores had remarkably higher expression of PD1, PD-L1, and CTLA4 (Figure S5B).

The necroptosis-associated gene signature was an exact predictor using fewer genes

Comparison of the predictive power of multiple genetic signatures can help to further explore the prognostic value of these genetic signatures. Therefore, we compared the predictive ability of the following five gene signatures: (1) the necroptosis-associated-gene signature constructed in this study; (2) the six-gene signature constructed by Liu et al. [24]; (3) the ten-gene signature constructed by Liang et al. [25]; (4) the twelve-gene signature constructed by Ouyang et al. [26] and (5) the fourteen-gene signature constructed by Zhang et al. [27]. As shown in Figure S6, our five-gene signature has the best survival prediction ability with fewer genes.

Building a Nomogram Model

The coefficient prediction efficacy of this necroptosis-associated-gene signature was investigated using a nomogram model in the TCGA dataset, and nomogram construction and validation were performed in accordance with the nomogram guidelines [28, 29]. The findings revealed that the nomogram might help us provide a quantitative approach for appropriately predicting the 1-, 3-, and 5-year survival rates (Figure S7A). The calibration curves revealed that the expected and actual likelihood of 1-, 3-, and 5-year survival rates were in good agreement (Figure S7B).

Expression levels of genes in the necroptosis-associated gene signature

Gene Expression Profiling Interactive Analysis (GEPIA) [30] was used first to investigate the expression levels of the five genes, and the results showed that only TNFRSF4 and LGALS3 varied substantially between normal and tumor samples (Figure S8A). Then, using the Human Protein Atlas database (HPA) [31], the protein expression of the five genes was investigated, and differences in protein levels were detected between HMOX1 and LGALS3 (Figure S8B). In addition, the expression levels of the five genes in HCC cell lines (Supplementary Figure S8C) were investigated using the Cancer Cell Line Encyclopedia (CCLE) database [32]. Finally, in clinical samples, qRT‒PCR and IHC assays were performed to corroborate the expression of genes in our model. As shown in Fig. 8A-E, TNFRSF4 and LGALS3 varied substantially between tumor tissues and their paracancerous normal tissues. Moreover, the scores obtained according to the formula in clinical samples could distinguish well between tumor and normal tissue (AUC = 0.882, Fig. 8F). Considering that only TNFRSF4 and LGALS3 differed in clinical samples and that no protein expression of TNFRSF4 was detected in the HPA database, we performed immunohistochemical staining analysis only for LGALS3. As shown in Fig. 8G, LGALS3 was remarkably elevated in tumor tissues.

Fig. 8
figure 8

Expression levels of the five genes analyzed by qRT-PCR and IHC assay. (A-E) Expression levels of the five genes between tumor tissues and their paracancerous normal tissues. (F) Scores in clinical samples could distinguish well between tumor and normal tissue. (G) LGALS3 was remarkably elevated in tumor tissues and increased with the progression of HCC.

Effect of LGALS3 knockdown on the proliferation and migratory capacity of HepG2 cells and the expression of key regulatory genes of necroptosis and inflammatory cytokines

In HepG2 cells, the siRNA significantly reduced LGALS3 expression, as shown in Fig. 9A. CCK-8 analysis revealed that downregulated LGALS3 had a lower proliferation rate (Fig. 9B), and wound-healing experiments revealed that LGALS3 knockdown significantly slowed wound closure (Fig. 9C). LGALS3 knockdown significantly decreased cell invasion in transwell chambers using the Matrigel test (Fig. 9D). Furthermore, knocking down LGALS3 greatly elevated the expression of three important necroptosis regulating genes (MLKL, RIPK1, and RIPK3, Fig. 9E), as well as the levels of inflammatory cytokines such as IL10 and TNF-α (Fig. 9F). Sametime, in order to further determine the protein expression of necroptosis regulating genes, western blot was performed LGALS3 knockdown and co-culture with necrotic apoptosis inhibitor (Nec-1). This was consistent with our expectation that MLKL, RIPK1, and RIPK3 were high expression in HepG2 cell line after LGALS3 knockdown but were significantly decreased after treatment with Nec-1(Figure S10).

Fig. 9
figure 9

Knockdown of LGALS3 inhibited HepG2 cells proliferation and migration and increased the expression of key regulatory genes of Necroptosis and inflammatory cytokines. (A) Decreased LGALS3 in HepG2 cells. (B) CCK-8 assay. (C) Wound-healing assay. (D) Transwell assay. (E) and (F) Increased MLKL, RIPK1, RIPK3, IL6, and TNF-α levels after knockdown of LGALS3.

Drug susceptibility analysis

Based on the results of drug sensitivity studies, 21 drugs were found to be tumor-sensitive, and the top 16 most important tumor-sensitive drugs are depicted in Supplementary Figure S9.

Discussion

Necroptosis is a cellular response to inflammation or infection that is distinct from apoptosis and regulated by the TNF-α receptor system [4]. Cells that die from apoptosis or necrosis have different effects on inflammation within the surrounding microenvironment. The plasma membrane of cells that die due to apoptosis remains intact and has a low inflammatory response around it; however, the plasma membrane of cells dying due to necroptosis ruptures, and the released intracellular contents directly activate and modulate the inflammatory response. There is growing evidence that necroptosis is involved in many states of pathology, including stroke [33], neurodegenerative diseases [5], myocardial infarction [34], liver or kidney injury [35, 36], pancreatitis [37], autoimmune diseases [38] and even cancers [39]. Considering that the expression of the major mediators of necroptosis can be regulated by gene deletion or modification in mouse models with no effect on mouse development [40], necroptosis is expected to be the next target for immunotherapy.

Previous studies have shown that intratumor necroptosis can be induced by hypoxia, malnutrition, chemotherapeutic or biological agents (such as cisplatin and bortezomib), or pro-necrotic stimulating factors secreted by tumor-infiltrating immune cells [41,42,43]. DNA methylation progressively suppresses RIPK3 protein expression during tumor progression to avoid necroptosis, while DNA hypomethylating agents can rescue the necrotic capacity of RIPK3 and promote the efficacy of antitumor immunity [44]. In addition, the release of DAMPs and cytokines in the TME by cancer cells in the necrotic state can facilitate the interaction between dying cancer cells and immune cells, thereby mediating the immune response associated with cancer [45]. Although several studies have confirmed the impact of necroptosis on the prognosis of cancer patients, the conclusions are inconsistent and may be related to tumor heterogeneity, cancer type, and different origin organs of the tumor [46]. Interestingly, methylation-dependent deletion of RIPK3 expression significantly inhibited necroptosis in a variety of common human hepatoma cell lines [44]. Sabira reported that increased necroptosis, which is closely associated with liver ageing, could cause chronic inflammation and promote liver fibrosis [47]. Chen found that heparinase could facilitate HCC metastasis by triggering necroptosis of microvascular endothelial cells [48].In contrast, Silvano reported that increased necroptosis by overexpressing SerpinB3 was beneficial in improving the survival of HCC patients [49]. Given the controversial prognostic role of necroptosis in oncology, there is an urgent need for a specific biomarker or spectrum technology that can consistently and specifically identify necroptosis in complex human tumor tissues.

In this study, after Z-scores of necroptosis were identified by the ssGSEA algorithm, they were shown to be strongly related to patient survival outcomes and involved in the immune infiltration of HCC. Subsequently, we established and validated an individualized prognostic profile associated with the necroptosis pathway using multiple bioinformatics methods to forecast therapeutic responses to ICB therapy. Moreover, when compared with previous signatures, our necroptosis-related gene signature had the best survival prediction ability with fewer genes. Among the five genes, HMOX1 was an unfavorable predictor of HCC prognosis [50], and its expression was also correlated with the efficacy of immunotherapy, such as sorafenib [51, 52]. In addition, HMOX1 was also associated with post-transplant HCC recurrence [53]. VNN2 is a gene related to hydrolase activity that has been identified as a prognostic indicator for HCC patients [54]. The TNFRSF4-NF-κB pathway is considered to be a promising therapeutic target for immune-related gene dysregulation in HCC [55]. High expression of KLRB1 is closely associated with early HCC relapse, and it was also involved in CD8 + T-cell associated immune escape mechanisms [56]. LGALS3 promotes tumor metastasis by activating the β-catenin signaling pathway in HCC, and it is also a reliable biomarker for predicting sorafenib resistance [57]. In this study, we found that LGALS3 not only had an effect on the proliferation and migration ability of HepG2 cells but also affected necroptosis and the expression of inflammatory cytokines.

Our study has certain limitations. First, individual variability in HCC patients may alter the expression of this necroptosis-related gene signature, and therefore future prospective multicenter randomized controlled trials are needed to assess the robustness and clinical usability of the five gene signatures. Furthermore, the biological mechanism between this five-gene signature and the necroptosis pathway in HCC needs to be clarified by biological experimental research in the future.

Conclusions

Overall, we found that necroptosis was associated with HCC progression and survival outcomes and was involved in the immune infiltration of HCC. Subsequently, we established and validated an individualized prognostic profile related to necroptosis using multiple bioinformatics methods to forecast the therapeutic response to immune therapy, and found that LGALS3 affected necroptosis, which might offer a potential nonapoptotic therapeutic target for HCC patients.

Data Availability

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

References

  1. Weinlich R, Oberst A, Beere HM, Green DR. Necroptosis in development, inflammation and disease. Nat Rev Mol Cell Biol. 2017;18(2):127–36.

    Article  CAS  PubMed  Google Scholar 

  2. Lalaoui N, Lindqvist LM, Sandow JJ, Ekert PG. The molecular relationships between apoptosis, autophagy and necroptosis. Semin Cell Dev Biol. 2015;39:63–9.

    Article  CAS  PubMed  Google Scholar 

  3. Huang CY, Kuo WT, Huang YC, Lee TC, Yu LC. Resistance to hypoxia-induced necroptosis is conferred by glycolytic pyruvate scavenging of mitochondrial superoxide in colorectal cancer cells. Cell Death Dis. 2013;4(5):e622.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Khoury MK, Gupta K, Franco SR, Liu B. Necroptosis in the pathophysiology of Disease. Am J Pathol. 2020;190(2):272–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Yuan J, Amin P, Ofengeim D. Necroptosis and RIPK1-mediated neuroinflammation in CNS diseases. Nat Rev Neurosci. 2019;20(1):19–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Strilic B, Yang L, Albarrán-Juárez J, Wachsmuth L, Han K, Müller UC, Pasparakis M, Offermanns S. Tumour-cell-induced endothelial cell necroptosis via death receptor 6 promotes metastasis. Nature. 2016;536(7615):215–8.

    Article  CAS  PubMed  Google Scholar 

  7. Feng X, Song Q, Yu A, Tang H, Peng Z, Wang X. Receptor-interacting protein kinase 3 is a predictor of survival and plays a tumor suppressive role in colorectal cancer. Neoplasma. 2015;62(4):592–601.

    Article  CAS  PubMed  Google Scholar 

  8. Lomphithak T, Akara-Amornthum P, Murakami K, Hashimoto M, Usubuchi H, Iwabuchi E, Unno M, Cai Z, Sasano H, Jitkaew S. Tumor necroptosis is correlated with a favorable immune cell signature and programmed death-ligand 1 expression in cholangiocarcinoma. Sci Rep. 2021;11(1):11743.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Gautheron J, Vucur M, Reisinger F, Cardenas DV, Roderburg C, Koppe C, Kreggenwinkel K, Schneider AT, Bartneck M, Neumann UP, et al. A positive feedback loop between RIP3 and JNK controls non-alcoholic steatohepatitis. EMBO Mol Med. 2014;6(8):1062–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Mohammed S, Nicklas EH, Thadathil N, Selvarani R, Royce GH, Kinter M, Richardson A, Deepa SS. Role of necroptosis in chronic hepatic inflammation and fibrosis in a mouse model of increased oxidative stress. Free Radic Biol Med. 2021;164:315–28.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Seehawer M, Heinzmann F, D’Artista L, Harbig J, Roux PF, Hoenicke L, Dang H, Klotz S, Robinson L, Doré G, et al. Necroptosis microenvironment directs lineage commitment in liver cancer. Nature. 2018;562(7725):69–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Wu L, Zhang X, Zheng L, Zhao H, Yan G, Zhang Q, Zhou Y, Lei J, Zhang J, Wang J, et al. RIPK3 orchestrates fatty acid metabolism in Tumor-Associated Macrophages and Hepatocarcinogenesis. Cancer Immunol Res. 2020;8(5):710–21.

    Article  CAS  PubMed  Google Scholar 

  13. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein science: a publication of the Protein Society. 2019;28(11):1947–51.

    Article  CAS  PubMed  Google Scholar 

  15. Kanehisa M, Furumichi M, Sato Y, Kawashima M, Ishiguro-Watanabe M. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 2023;51(D1):D587–d592.

    Article  CAS  PubMed  Google Scholar 

  16. Zhang K, Zhang L, Mi Y, Tang Y, Ren F, Liu B, Zhang Y, Zheng P. A ceRNA network and a potential regulatory axis in gastric cancer with different degrees of immune cell infiltration. Cancer Sci. 2020;111(11):4041–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Gentles AJ, Newman AM, Liu CL, Bratman SV, Feng W, Kim D, Nair VS, Xu Y, Khuong A, Hoang CD, et al. The prognostic landscape of genes and infiltrating immune cells across human cancers. Nat Med. 2015;21(8):938–45.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, Li B, Liu XS. TIMER: a web server for Comprehensive Analysis of Tumor-Infiltrating Immune cells. Cancer Res. 2017;77(21):e108–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Finotello F, Mayer C, Plattner C, Laschober G, Rieder D, Hackl H, Krogsdam A, Loncova Z, Posch W, Wilflingseder D, et al. Molecular and pharmacological modulators of the tumor immune contexture revealed by deconvolution of RNA-seq data. Genome Med. 2019;11(1):34.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Aran D, Hu Z, Butte AJ. xCell: digitally portraying the tissue cellular heterogeneity landscape. Genome Biol. 2017;18(1):220.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Reinhold WC, Sunshine M, Liu H, Varma S, Kohn KW, Morris J, Doroshow J, Pommier Y. CellMiner: a web-based suite of genomic and pharmacologic tools to explore transcript and drug patterns in the NCI-60 cell line set. Cancer Res. 2012;72(14):3499–511.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Yoshihara K, Shahmoradgoli M, Martínez E, Vegesna R, Kim H, Torres-Garcia W, Treviño V, Shen H, Laird PW, Levine DA, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. 2013;4:2612.

    Article  PubMed  Google Scholar 

  24. Liu GM, Zeng HD, Zhang CY, Xu JW. Identification of a six-gene signature predicting overall survival for hepatocellular carcinoma. Cancer Cell Int. 2019;19:138.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Liang JY, Wang DS, Lin HC, Chen XX, Yang H, Zheng Y, Li YH. A Novel ferroptosis-related gene signature for overall survival prediction in patients with Hepatocellular Carcinoma. Int J Biol Sci. 2020;16(13):2430–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Zhang BH, Yang J, Jiang L, Lyu T, Kong LX, Tan YF, Li B, Zhu YF, Xi AY, Xu X, et al. Development and validation of a 14-gene signature for prognosis prediction in hepatocellular carcinoma. Genomics. 2020;112(4):2763–71.

    Article  CAS  PubMed  Google Scholar 

  28. Balachandran VP, Gonen M, Smith JJ, DeMatteo RP. Nomograms in oncology: more than meets the eye. Lancet Oncol. 2015;16(4):e173–180.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Li Z, Sun Y, Xu J, Yang H, Liu X, Tian Y, Cao S, Zhou Y. Integrin-beta6 serves as a potential prognostic serum biomarker for gastric Cancer. Front Oncol. 2021;11:770997.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45(W1):W98–w102.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Thul PJ, Lindskog C. The human protein atlas: a spatial map of the human proteome. Protein Sci. 2018;27(1):233–44.

    Article  CAS  PubMed  Google Scholar 

  32. Nusinow DP, Szpyt J, Ghandi M, Rose CM, McDonald ER 3rd, Kalocsay M, Jané-Valbuena J, Gelfand E, Schweppe DK, Jedrychowski M, et al. Quantitative proteomics of the Cancer Cell Line Encyclopedia. Cell. 2020;180(2):387–402e316.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Degterev A, Huang Z, Boyce M, Li Y, Jagtap P, Mizushima N, Cuny GD, Mitchison TJ, Moskowitz MA, Yuan J. Chemical inhibitor of nonapoptotic cell death with therapeutic potential for ischemic brain injury. Nat Chem Biol. 2005;1(2):112–9.

    Article  CAS  PubMed  Google Scholar 

  34. Luedde M, Lutz M, Carter N, Sosna J, Jacoby C, Vucur M, Gautheron J, Roderburg C, Borg N, Reisinger F, et al. RIP3, a kinase promoting necroptotic cell death, mediates adverse remodelling after myocardial infarction. Cardiovasc Res. 2014;103(2):206–16.

    Article  CAS  PubMed  Google Scholar 

  35. Ramachandran A, McGill MR, Xie Y, Ni HM, Ding WX, Jaeschke H. Receptor interacting protein kinase 3 is a critical early mediator of acetaminophen-induced hepatocyte necrosis in mice. Hepatology. 2013;58(6):2099–108.

    Article  CAS  PubMed  Google Scholar 

  36. Linkermann A, De Zen F, Weinberg J, Kunzendorf U, Krautwald S. Programmed necrosis in acute kidney injury. Nephrol Dial Transplant. 2012;27(9):3412–9.

    Article  CAS  PubMed  Google Scholar 

  37. Wu J, Huang Z, Ren J, Zhang Z, He P, Li Y, Ma J, Chen W, Zhang Y, Zhou X, et al. Mlkl knockout mice demonstrate the indispensable role of Mlkl in necroptosis. Cell Res. 2013;23(8):994–1006.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Saito N, Honma M, Shibuya T, Iinuma S, Igawa S, Kishibe M, Ishida-Yamamoto A. RIPK1 downregulation in keratinocyte enhances TRAIL signaling in psoriasis. J Dermatol Sci. 2018;91(1):79–86.

    Article  CAS  PubMed  Google Scholar 

  39. Gong Y, Fan Z, Luo G, Yang C, Huang Q, Fan K, Cheng H, Jin K, Ni Q, Yu X, et al. The role of necroptosis in cancer biology and therapy. Mol Cancer. 2019;18(1):100.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Newton K, Dugger DL, Wickliffe KE, Kapoor N, de Almagro MC, Vucic D, Komuves L, Ferrando RE, French DM, Webster J, et al. Activity of protein kinase RIPK3 determines whether cells die by necroptosis or apoptosis. Science. 2014;343(6177):1357–60.

    Article  CAS  PubMed  Google Scholar 

  41. Grivennikov SI, Greten FR, Karin M. Immunity, inflammation, and cancer. Cell. 2010;140(6):883–99.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Sancho D, Joffre OP, Keller AM, Rogers NC, Martínez D, Hernanz-Falcón P, Rosewell I. Reis e Sousa C: identification of a dendritic cell receptor that couples sensing of necrosis to immunity. Nature. 2009;458(7240):899–903.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Jing L, Song F, Liu Z, Li J, Wu B, Fu Z, Jiang J, Chen Z. MLKL-PITPα signaling-mediated necroptosis contributes to cisplatin-triggered cell death in lung cancer A549 cells. Cancer Lett. 2018;414:136–46.

    Article  CAS  PubMed  Google Scholar 

  44. Koo GB, Morgan MJ, Lee DG, Kim WJ, Yoon JH, Koo JS, Kim SI, Kim SJ, Son MK, Hong SS, et al. Methylation-dependent loss of RIP3 expression in cancer represses programmed necrosis in response to chemotherapeutics. Cell Res. 2015;25(6):707–25.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. McCracken MN, Cha AC, Weissman IL. Molecular Pathways: activating T cells after Cancer Cell phagocytosis from Blockade of CD47 “Don’t eat Me” signals. Clin Cancer Res. 2015;21(16):3597–601.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Sprooten J, De Wijngaert P, Vanmeerbeerk I, Martin S, Vangheluwe P, Schlenner S, Krysko DV, Parys JB, Bultynck G, Vandenabeele P et al. Necroptosis in Immuno-Oncology and Cancer Immunotherapy. Cells 2020, 9(8).

  47. Mohammed S, Thadathil N, Selvarani R, Nicklas EH, Wang D, Miller BF, Richardson A, Deepa SS. Necroptosis contributes to chronic inflammation and fibrosis in aging liver. Aging Cell. 2021;20(12):e13512.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Chen X, Cheng B, Dai D, Wu Y, Feng Z, Tong C, Wang X, Zhao J. Heparanase induces necroptosis of microvascular endothelial cells to promote the metastasis of hepatocellular carcinoma. Cell Death Discov. 2021;7(1):33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Fasolato S, Ruvoletto M, Nardo G, Rasola A, Sciacovelli M, Zanus G, Turato C, Quarta S, Terrin L, Fadini GP et al. Low P66shc with high SerpinB3 levels favors necroptosis and better survival in Hepatocellular Carcinoma. Biology (Basel) 2021, 10(5).

  50. Jiang C, He ZL, Hu XH, Ma PY. MiRNA-15a-3p inhibits the metastasis of hepatocellular carcinoma by interacting with HMOX1. Eur Rev Med Pharmacol Sci. 2020;24(24):12694–700.

    CAS  PubMed  Google Scholar 

  51. Wang Y, Xie Y, Ma J, Wang Y, Gong R. Development and validation of a prognostic and immunotherapeutically relevant model in hepatocellular carcinoma. Ann Transl Med. 2020;8(18):1177.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Jiang HY, Ning G, Wang YS, Lv WB. Ahypoxia-related signature enhances the prediction of the prognosis in hepatocellular carcinoma patients and correlates with sorafenib treatment response. Am J Transl Res. 2020;12(12):7762–81.

    CAS  PubMed  PubMed Central  Google Scholar 

  53. Oldani G, Crowe LA, Orci LA, Slits F, Rubbia-Brandt L, de Vito C, Morel P, Mentha G, Berney T, Vallée JP, et al. Pre-retrieval reperfusion decreases cancer recurrence after rat ischemic liver graft transplantation. J Hepatol. 2014;61(2):278–85.

    Article  PubMed  Google Scholar 

  54. Li W, Lu J, Ma Z, Zhao J, Liu J. An Integrated Model based on a six-gene signature predicts overall survival in patients with Hepatocellular Carcinoma. Front Genet. 2019;10:1323.

    Article  CAS  PubMed  Google Scholar 

  55. Li WT, Zou AE, Honda CO, Zheng H, Wang XQ, Kisseleva T, Chang EY, Ongkeko WM. Etiology-specific analysis of Hepatocellular Carcinoma Transcriptome reveals genetic dysregulation in pathways implicated in Immunotherapy Efficacy. Cancers (Basel) 2019, 11(9).

  56. Sun Y, Wu L, Zhong Y, Zhou K, Hou Y, Wang Z, Zhang Z, Xie J, Wang C, Chen D, et al. Single-cell landscape of the ecosystem in early-relapse hepatocellular carcinoma. Cell. 2021;184(2):404–421e416.

    Article  CAS  PubMed  Google Scholar 

  57. Song M, Pan Q, Yang J, He J, Zeng J, Cheng S, Huang Y, Zhou ZQ, Zhu Q, Yang C, et al. Galectin-3 favours tumour metastasis via the activation of β-catenin signalling in hepatocellular carcinoma. Br J Cancer. 2020;123(10):1521–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We thank the patients and their families for their signed informed consent forms in clinical sample experiment.

Funding

This work is supported by PhD research startup foundation of the Third Affiliated Hospital of Zhengzhou University (2021079).

Author information

Authors and Affiliations

Authors

Contributions

KZ designed this work and analyzed the data; JPL performed experiments and analyses; EWY helped for providing tumor samples; KZ and EWY wrote the manuscript.

Corresponding authors

Correspondence to Kai Zhang or Enwu Yuan.

Ethics declarations

Competing interests

The authors declare no competing interests.

Ethics approval and consent to participate

All patients provided an informed consent before undergoing any study-specifc procedures. The study was conducted in accordance with the principles of the Declaration of Helsinki. The study was approved by the Ethics Committee of Zhengzhou University.

Consent for publication

Not Applicable.

Competing interests

The authors declare that the authors have no conflict of interest that might be perceived to influence the results and/or discussion reported in this paper.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Supplementary Material 2

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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, K., Li, J. & Yuan, E. A necroptosis-related gene signature to predict prognosis and immune features in hepatocellular carcinoma. BMC Cancer 23, 660 (2023). https://doi.org/10.1186/s12885-023-11168-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-023-11168-8

Keywords