Skip to main content

Screening prognostic markers for hepatocellular carcinoma based on pyroptosis-related lncRNA pairs

Abstract

Background

Pyroptosis is closely related to cancer prognosis. In this study, we tried to construct an individualized prognostic risk model for hepatocellular carcinoma (HCC) based on within-sample relative expression orderings (REOs) of pyroptosis-related lncRNAs (PRlncRNAs).

Methods

RNA-seq data of 343 HCC samples derived from The Cancer Genome Atlas (TCGA) database were analyzed. PRlncRNAs were detected based on differentially expressed lncRNAs between sample groups clustered by 40 reported pyroptosis-related genes (PRGs). Univariate Cox regression was used to screen out prognosis-related PRlncRNA pairs. Then, based on REOs of prognosis-related PRlncRNA pairs, a risk model for HCC was constructed by combining LASSO and stepwise multivariate Cox regression analysis. Finally, a prognosis-related competing endogenous RNA (ceRNA) network was built based on information about lncRNA–miRNA–mRNA interactions derived from the miRNet and TargetScan databases.

Results

Hierarchical clustering of HCC patients according to the 40 PRGs identified two groups with a significant survival difference (Kaplan–Meier log-rank, p = 0.026). Between the two groups, 104 differentially expressed lncRNAs were identified (|log2(FC)|> 1 and FDR < 5%). Among them, 83 PRlncRNA pairs showed significant associations between their REOs within HCC samples and overall survival (Univariate Cox regression, p < 0.005). An optimal 11-PRlncRNA-pair prognostic risk model was constructed for HCC. The areas under the curves (AUCs) of time-dependent receiver operating characteristic (ROC) curves of the risk model for 1-, 3-, and 5-year survival were 0.737, 0.705, and 0.797 in the validation set, respectively. Gene Set Enrichment Analysis showed that inflammation-related interleukin signaling pathways were upregulated in the predicted high-risk group (p < 0.05). Tumor immune infiltration analysis revealed a higher abundance of regulatory T cells (Tregs) and M2 macrophages and a lower abundance of CD8 + T cells in the high-risk group, indicating that excessive pyroptosis might occur in high-risk patients. Finally, eleven lncRNA–miRNA–mRNA regulatory axes associated with pyroptosis were established.

Conclusion

Our risk model allowed us to determine the robustness of the REO-based PRlncRNA prognostic biomarkers in the stratification of HCC patients at high and low risk. The model is also helpful for understanding the molecular mechanisms between pyroptosis and HCC prognosis. High-risk patients may have excessive pyroptosis and thus be less sensitive to immune therapy.

Peer Review reports

Background

Hepatocellular carcinoma (HCC) is one of the most common human malignancies, with over 800,000 new cases and nearly 700,000 deaths worldwide yearly [1]. HCC is highly heterogeneous and insidious; most patients are diagnosed at advanced stages with poor prognosis [2]. According to statistics, the overall median survival time of patients with advanced HCC is only 9 months, and the 5-year overall survival (OS) is solely 10% [3]. Therefore, exploring the molecular biomarkers associated with HCC prognosis has been a hot issue in HCC research [4].

Pyroptosis is a sort of programmed cell death related to inflammation, which is mediated by the intracellular inflammasome and gasdermins [5, 6]. Appropriate induction of pyroptosis could trigger a moderate inflammatory reaction that might enhance innate immunity and generate an antitumor immune response [7, 8]. In contrast, excessive pyroptosis might excite a hyperinflammatory response that disrupts immune homeostasis and promotes cancer progression [7, 8]. Consequently, some researchers attempt to improve the prognosis of tumor patients by regulating the activation of pyroptosis to produce antitumor immune responses [9]. For instance, migration and invasion of oral squamous cell carcinoma cells could be inhibited via pyroptosis activation by anthocyanidins [10]. High expression of gasdermin E induced by miltirone could be used to provoke pyroptosis in cancer cells [11]. NLRP3 (NLR Family Pyrin Domain Containing 3) inflammasomes could be utilized to mediate pyroptosis to suppress the growth and metastasis of HCC cells [12]. These suggest pyroptosis is closely associated with cancer prognosis.

Long noncoding RNA (lncRNA) is a transcription product of DNA with a length greater than 200 nucleotides, which can regulate gene expression by interacting with proteins, DNA, or other RNAs [13]. It has been reported that lncRNAs are critical regulators of pyroptosis [14]. For example, lncRNA HOTTIP could inhibit pyroptosis and promote ovarian cancer cell proliferation by targeting miR-148a-3p/AKT2 axis [15]. LncRNA MEG3 could inhibit the growth and metastasis of triple-negative breast cancer by activating pyroptosis via NLRP3/caspase-1/GSDMD pathway [16]. These studies demonstrate that aberrant alterations of pyroptosis-related lncRNAs (PRlncRNAs) also impact cancer prognosis.

Prognosis-associated PRlncRNA biomarkers have been identified for multiple cancer types. For HCC, seven and five prognosis-associated PRlncRNAs have been reported [17, 18]. These PRlncRNA prognostic risk models have some efficacy in training and testing datasets. However, their risk thresholds are summarized from the absolute expression levels of PRlncRNAs, which are often data-dependent and unstable, leading to difficulties when applied in clinical settings [19]. In recent years, researchers have found that the within-sample relative expression orderings (REOs) of genes were more robust than the absolute expression levels of genes across samples. Furthermore, the REO-based molecular biomarkers can be easily applied to individual diagnosis, which is more suitable for clinics [20,21,22].

In this study, we tried identifying REO-based prognosis-associated PRlncRNA biomarkers to construct an individualized prognostic risk model for HCC and explore the molecular mechanisms between pyroptosis and HCC prognosis.

Materials and methods

Data collection and preprocessing

The RNA expression data of the 371 HCC and 50 adjacent normal tissue samples analyzed in this study were downloaded from The Cancer Genome Atlas (TCGA) database. The data were obtained using Illumina HiSeq 2000 RNA Sequencing technology. A total of 60,483 RNAs were detected. The data were preprocessed by the following procedures. Firstly, remove the samples with a survival time of less than 30 days or missing survival time. Secondly, normalize expression values to Transcript Per Million (TPM). Thirdly, annotate each examined RNA in the GENCODE database. Finally, exclude mRNAs with count values less than 1 and lncRNAs with count values less than 0.5 in all samples.

For miRNA data downloaded from the TCGA database, sequencing was performed on an Illumina HiSeq miRNASeq platform. After deleting the miRNAs with a count value of 0 in more than 50% of the samples, 578 miRNAs were kept. The miRNA expression values were normalized to Reads Per Million (RPM).

Detection of prognosis-related lncRNA pairs

Clustering HCC samples with the ward linkage algorithm were performed on 40 pyroptosis-related genes (PRGs), which were derived from the pyroptosis pathway in the Molecular Signatures Database (MSigDB) (Additional file 1: Table S1). Then, differentially expressed lncRNAs detected between the two HCC groups were defined as pyroptosis-related lncRNAs. Any two pyroptosis-related lncRNAs can form a lncRNA pair. For a lncRNA pair (lncRNAi|lncRNAj), there are two REO status, lncRNAi < lncRNAj or lncRNAi ≥ lncRNAj, in a sample. The prognosis-related lncRNA pairs were then identified by the following procedures.

  1. (1)

    Combine the lncRNAs of interest two by two to form \({C}_{k}^{2}\) (k is the number of pyroptosis-related lncRNAs) lncRNA pairs.

  2. (2)

    The REO matrix, X, was constructed based on lncRNA pairs for the training set. Xij denoted the REO of the i-th lncRNA pair (lncRNAi1|lncRNAi2) in the j-th sample, taking 1 or 0, with 1 representing lncRNAi1 < lncRNAi2 and 0 representing lncRNAi1 ≥ lncRNAi2.

  3. (3)

    Remove lncRNA pairs for which the percentage of REOs with 1 was less than 20% or greater than 80%. This criterion ensured that the apparent reversal of REOs of lncRNA pairs occurred in a certain amount of HCC samples to facilitate the identification of sample subgroups with different prognoses, as lncRNA pairs with the same score (0 or 1) in more than 80% of samples were considered uninformative [21].

  4. (4)

    For each remaining lncRNA pair, the correlation between its REO values and OS times was evaluated by univariate Cox regression. If the Wald test's p value is less than 0.005, the lncRNA pair was considered as a prognosis-related lncRNA pair.

Construction and evaluation of prognostic lncRNA pair risk model

The LASSO regression and stepwise multivariate Cox regression algorithms were applied to select candidate PRlncRNA pairs as prognostic biomarkers. LASSO regression was adopted to choose the prognosis-related PRlncRNA pairs most predictive of OS. Lambda values corresponding to the smallest partial likelihood deviance were chosen as the optimal parameters after tenfold cross-validation [23, 24]. Then, multivariate Cox regression analysis based on the Akaike information criterion (AIC) method was used to determine the optimal model. The model with the lowest AIC value was considered as the optimal prognostic risk model, with the corresponding PRlncRNA pairs as the eventual predictive risk biomarkers [25].

The time-dependent receiver operating characteristic (ROC) curves were applied to assess the performance, and the Youden index determined the risk threshold. Multivariate Cox proportional hazards regression analysis was employed to evaluate independent prognostic factors associated with OS [26]. Covariates included risk scores for prognostic PRlncRNA pairs, gender, age, tumor stage, and grade.

Enrichment analysis and estimation of immune cell infiltration

Functional enrichment analysis of differentially expressed genes between the high- and low-risk groups was completed by Gene Set Enrichment Analysis (GSEA) based on the Reactome database with annotation information from the MSigDB database (v7.5.1).

The estimation of the absolute abundance of tumor-infiltrating (immune cells in HCC samples was achieved by the CIBERSORT algorithm [27].

Construction of prognostic pyroptosis-related competing endogenous RNA (ceRNA) network

The regulatory relationships of lncRNAs, miRNAs, and mRNAs were obtained from the miRNet database and the TargetScan database [28, 29]. The miRNet database was utilized to predict target miRNAs for lncRNAs, and the TargetScan database was used to predict target miRNAs for PRGs.

Statistical analysis

All statistical analysis were completed with R 4.1.0 software. Cluster analysis was finished with the ward linkage algorithm. Differential expression analysis was implemented with the "limma"package. Survival analysis and corresponding plotting were based on the "survival", "glmnet", "MASS" and "survminer" packages. ROC analysis and the determination of risk threshold were completed based on the "survivalROC" package. GSEA was based on the "clusterProfiler" package. Tumor immune infiltration analysis was implemented by the "immunedeconv" package. The Benjamini–Hochberg (BH) method was applied to control the false discovery rate (FDR). Unless otherwise specified, the statistical significance level was set uniformly at 0.05.

Results

Pyroptosis-related lncRNAs

The workflow of this study is illustrated in Fig. 1. After data preprocessing, the TCGA HCC data set included expression measurements of 8477 lncRNAs and 17,596 mRNAs from 343 HCC samples. Firstly, all samples were randomly categorized into training (n = 240) and validation (n = 103) sets. No significant differences have been observed in the clinical features between the training and validation sets (p value < 0.05, Additional file 1: Table S2).

Fig. 1
figure 1

Workflow of this study

Hierarchical clustering was performed on the expression levels of the 40 PRGs in the training set. Samples were clustered into two groups, containing 35 and 205 samples, respectively (Fig. 2A). Kaplan–Meier survival analysis revealed a significant difference in survival between these two groups of patients (log-rank test, p value = 0.026 (Fig. 2B)), which suggested that the expression pattern of PRGs was associated with the prognosis of HCC patients.

Fig. 2
figure 2

Hierarchical clustering and Kaplan–Meier survival analysis. A Hierarchical clustering for the 40 PRGs in the training set. B Kaplan–Meier survival curves for clusters derived from hierarchical clustering

Using the R package "limma", 104 lncRNAs that were differentially expressed between the two groups were identified at |log2(FC)|> 1 and FDR < 5%. These 104 lncRNAs were considered as PRlncRNAs for the following analysis.

Prognosis-related lncRNA pairs and risk model

The 104 PRlncRNAs were paired, and prognosis-related lncRNA pairs were identified based on the within-sample REOs in the training set (see Methods). The REOs of 83 PRlncRNA pairs were observed to be significantly associated with the OS time of HCC by univariate Cox regression analysis. To choose representative prognosis-related lncRNA pairs, we performed the LASSO Cox regression analysis via tenfold cross-validation on these 83 PRlncRNA pairs, and 26 PRlncRNA pairs were selected at the smallest partial likelihood deviance (Fig. 3A, B). Then, the stepwise multivariate Cox regression analysis was performed on 26 PRlncRNA pairs to choose prognosis-related lncRNA pair biomarkers and construct the risk model. Finally, as shown in Fig. 3C, 11 PRlncRNA pairs involving 22 PRlncRNAs were selected at the smallest AIC value. The corresponding risk model is: risk score = 0.5447 × VIM-AS1|AC005083.1 + 0.7797 × LINC01057|RP11-43F13.3 + 1.0279 × TNRC6C-AS1|RP11-395G23.3 − 0.7839 × NRSN2-AS1|LINC01554 + 0.7766 × PCED1B-AS1|AC079466.1 + 1.0191 × LINC00342|CASC9 − 0.7017 × RP1-239B22.5|RP11-344B5.2 + 0.7464 × PSMB8-AS1|ANKRD10-IT1 + 0.6299 × KB-68A7.1|ZFPM2-AS1 − 0.5646 × FOXD2-AS1|AC092580.4 − 1.3199 × LINC00942|RP11-109M17.2 (Fig. 3C). The threshold for high- and low-risk groups was determined by the point with the largest Youden index on the 5-year ROC curve of the training set (Youden index = 0.739, risk score = 0.025), with 105 and 135 patients classified as high- and low-risk samples, respectively. Among the 22 PRlncRNAs, 16 were differentially expressed between high- and low-risk patients. Comparing all HCC samples to adjacent normal samples, 14 of these 16 PRlncRNAs were differentially expressed.

Fig. 3
figure 3

Establishment of the prognostic lncRNA pair risk model. A The plot of partial likelihood deviance versus log(λ). A vertical dotted line marks the smallest partial likelihood deviance; B LASSO coefficient profiles of HCC OS-associated lncRNA pairs; C forest plot depicting associations between lncRNA pairs and risk values determined by multivariate Cox regression analysis. *p < 0.05; **p < 0.01; ***p < 0.001

Validation of prognostic lncRNA pair risk model

According to the risk threshold, the patients in the validation set were divided into a high-risk group (risk score ≥ 0.025) and a low-risk group (risk score < 0.025), respectively. Kaplan–Meier survival analysis showed that OS rates were significantly different between the high- and low-risk groups in both the training and validation sets (log-rank test: p value < 2 × 10–16 and p value = 1.219 × 10–5, Fig. 4A, B). The AUCs of time-dependent ROC curves showed the prediction accuracies of 1-, 3- and 5-year survival were 0.855, 0.891, and 0.902 for the training set, and 0.737, 0.705, and 0.797 for the validation set, respectively (Fig. 4C, D). In addition, the results of multivariate Cox regression indicated that the risk model was an independent prognostic factor for patients with HCC (p value < 0.001, Fig. 4E, F). These results suggest that the risk score model based on lncRNA pairs can be an efficient tool for predicting the prognostic risk of HCC.

Fig. 4
figure 4

Validation of the prognostic lncRNA pair risk model. A, B Kaplan–Meier survival curves for high- and low-risk groups in the training and validation sets; C, D ROC curves based on the predictive efficacy of the lncRNA pair risk model for the training and validation sets; E, F Multivariate Cox regression analysis showed the independence of the risk model for the prediction for the training and validation sets

Analysis of immune-related characteristics of high- and low-risk groups

Among the 40 PRGs, 27 were differentially expressed between high- and low-risk groups in the training set at FDR < 5%. Notably, 25 PRGs were significantly upregulated in the high-risk group. GSEA showed the inflammation-related interleukin-mediated signaling pathways, including the interleukin 4 and interleukin 13 signaling pathways (q-value = 0.021), and signaling by interleukins pathway (q-value = 0.030), were also found to be significantly upregulated in the high-risk group [8]. Furthermore, as shown in Fig. 5, the abundance of regulatory T cells (Tregs) and M2 macrophages was significantly higher in the high-risk group. In comparison, the abundance of CD8 + T cells was significantly lower in the high-risk group. These results suggested that the immune system might have overreacted in the high-risk group due to the upregulation of PRGs expression, disrupting the immune homeostasis and making the prognosis worse.

Fig. 5
figure 5

Immune infiltration analysis between high- and low-risk groups

Establishment of prognostic pyroptosis-related ceRNA network

The target miRNAs of the 22 lncRNAs involved in the 11 PRlncRNA pair biomarkers in the risk score model were predicted by the miRNet database. Finally, 140 target miRNAs of 10 lncRNAs were obtained. The TargetScan database showed that 112 of the 140 miRNAs targeted 39 PRGs. By univariate Cox regression, three of the 112 target miRNAs and eight of the 39 target PRGs were significantly associated with the prognosis of HCC in our data. Moreover, five of the eight PRGs were targets of the three miRNAs. These prognosis-related lncRNAs, miRNAs, and PRGs formed eleven lncRNA–miRNA–mRNA regulatory axes (Fig. 6), involving four lncRNAs, three miRNAs, and five PRGs.

Fig. 6
figure 6

Pyroptosis-associated prognostic ceRNA network

Discussion

In this study, we have constructed a prognostic risk model for HCC by survival analysis of PRlncRNA pairs based on the within-sample REOs of PRlncRNAs. The risk model has good predictive performance at classifying the HCC patients into high- and low-risk groups in the validation set. Given the use of within-sample REOs of PRlncRNA pairs, our prognostic risk model has potential implications for clinical translation and application. Our model is independent of systematic bias and suitable for the individualized clinic and may help to stratify HCC patients at high risk of poor prognosis.

The 240 training samples were clustered by 40 PRGs into two groups of samples with significantly different prognoses, containing 35 cases with relatively poor survival and 205 cases with relatively good survival, respectively. Moreover, these 240 samples were divided by our prognostic prediction model into 105 high-risk cases and 135 low-risk cases. They were significantly overlapped (p = 3.45 × 10–4, hypergeometric test): the 35 poor prognosis samples overlapped with the 105 high-risk samples by 25 and the 205 good prognosis samples overlapped with the 135 low-risk samples by 125. The discrepancy may be due to the different purposes of stratifying samples. The clustering of samples using PRGs was intended to identify differential lncRNAs as PRlncRNAs. The two groups clustered by them may reflect more the similarity of expression patterns of the samples across all PRGs. The high- and low-risk groups were stratified by the constructed prognostic model. They should be more associated with the prognosis of HCC.

We found that more PRGs were upregulated in the high-risk group compared with the low-risk group, indicating more inflammation responses in the high-risk group. Further evidence was provided by the pathway enrichment and immune infiltration analysis results. GSEA showed that interleukin-mediated inflammation-related signaling pathways were upregulated in the high-risk group. CIBERSORT-based analysis showed higher abundance of Tregs and M2 macrophages and a lower abundance of CD8 + T cells in the high-risk group. It has been reported that interleukin 4 and interleukin 13 signaling could induce type 2 inflammatory processes [30]. If the type 2 inflammatory responses were out of control, M2 polarization of macrophages could be promoted, effectively suppressing the cytotoxicity of CD8 + T cells and NK cells [30, 31]. Except for M2 macrophages, Tregs, which function could be enhanced by chronic inflammation, could also efficiently inhibit the function of CD8 + T cells [8, 32]. Relative low abundance of CD8 + T cells has been reported to indirectly induce weaker cytotoxicity, while lower cytotoxicity might induce more insensitivity to immunotherapy [33]. Therefore, we additionally analyzed the cytotoxic-related genes (GZMA, GZMB, GZMK, PRF1) [34] and observed downregulation of these genes in the high-risk group compared to low-risk group (Wilcoxon rank-sum test: p value < 0.05; p value < 0.05; p value < 0.05; p value < 0.001). We then applied the Immune Cell Abundance Identifier (ImmuCellAI) database to predict the immunotherapeutic responses. We found that patients in the high-risk group were likely to have lower scores and be less sensitive to immune checkpoint blockade therapy (Wilcoxon rank-sum test: p value = 0.010) [35]. Therefore, we inferred that excessive pyroptosis might have arisen in high-risk patients, reducing the amount and activity of tumor-infiltrating lymphocytes and worsening tumor prognosis [36].

In clinical practice, due to the simplicity and non-invasiveness, serum markers such as AFP, DCP, and AFP-L3 are often used to diagnose HCC and predict prognosis. Many scoring models have used the three markers for the diagnosis or prognosis of HCC. The GALAD, consisting of age, sex, and the three markers, has been reported to have high predictive accuracy for early HCC in patients with nonalcoholic steatohepatitis (AUC = 0.96) [37] and also accurately classified patients with HCC in stage 0/A of Barcelona Clinic Liver Cancer (AUC = 0.9242) [38]. Studies have found that another scoring model, BALAD-2, consisting of these three markers combined with serum bilirubin and albumin, could stratify the HCC patients into distinct prognostic groups [39]. Furthermore, in 2008, the combined biomarker Japan Integrated Staging was proposed to provide better survival predictions for HCC patients [40]. These studies illustrated the potential of these three markers in the diagnosis and prognosis prediction of HCC. However, in the TCGA data we used, only the serum levels of AFP and DCP were provided. Therefore, we could not directly compare the predictive efficacy of our model with BALAD-2. We compared the AFP and DCP levels between the high- and low-risk groups predicted by our model. Results showed that these two proteins were not significantly differentially expressed between the training set's high- and low-risk groups (Wilcoxon–Mann–Whitney test: p = 0.486 and p = 0.771, respectively,). Elevated serum levels of AFP, AFP-L3, and DCP at baseline had been reported to be associated with a worse prognosis after resection of HCC [41]. We thus compared the corresponding mRNAs in the HCC patients, which showed that the mRNAs of the two markers were significantly up-regulated in the high-risk HCC samples compared to the normal controls (Wilcoxon–Mann–Whitney test: p < 0.05).

To further validate the predictive value of our prognostic PRlncRNA risk model, we compared its performance with three different prognostic models previously reported, which were also constructed based on PRlncRNAs using the same TCGA RNA sequencing data. Zhang et al. constructed a risk-scoring model of 5 PRlncRNAs, with 5-year AUCs of 0.688 for the training and 0.714 for the validation set, respectively [42]. Liu et al. built a prognostic risk-scoring model using 5 PRlncRNAs, with 5-year AUCs of 0.707 for the training set and 0.642 for the validation set, respectively [18]. The 5-year AUCs for the 9-PRlncRNA model built by Zhang et al. were 0.812 for the training set, and 0.722 for the validation set, respectively [43]. The predictive performance of our PRlncRNA prognostic model was higher than the three published models, with 5-year AUCs of 0.902 for the training set and 0.797 for the validation set, respectively.

The prognostic risk model in our study was constructed based on the REOs of PRlncRNAs. Although it is technically simpler to detect serum levels of protein markers such as the commonly used AFP, analysis at the RNA level may provide additional information for understanding cancer mechanistically. Most lncRNAs involved in the lncRNA–miRNA-PRG regulatory axes have been reported to be prognostically relevant in various cancers, including HCC. For example, LINC01554-mediated glucose metabolism reprogramming could suppress the tumorigenicity of HCC through the downregulation of PKM2 expression and inhibition of the Akt/mTOR signaling pathway [44]. NRSN2-AS1 could promote ovarian carcinogenesis through the miR-744-5p/PRKX axis [45]. Upregulation of lncRNA FOXD2-AS1 expression could promote the progression of HCC by causing epigenetic silencing of DKK1 and activating the Wnt/β-catenin signaling pathway [46]. Thus, these eleven lncRNA–miRNA-PRG regulatory axes could be helpful for further understanding the relationship between lncRNAs and PRGs and deserve further investigation.

Conclusion

In summary, we propose an 11-PRlncRNA-pair personalized prognostic risk model for HCC, which let us see the robustness of the REO-based PRlncRNA prognostic biomarkers in the stratification of HCC patients at high and low risk. And the model is also helpful for understanding the molecular mechanisms between pyroptosis and HCC prognosis. High-risk patients may have excessive pyroptosis and thus be less sensitive to immune therapy.

Availability of data and materials

The original contributions presented in the study are included in the article or supplementary material, and further inquiries can be directed to the corresponding authors. The original data for the analysis in this study included RNA-seq data, miRNA-seq data and corresponding clinical data available at TCGA database (https://www.cancer.gov/about-nci/organization/ccg/research/structural-genomics/tcga), annotation information of all pathways available at the MSigDB database (https://www.gsea-msigdb.org/gsea/msigdb), information about miRNA-lncRNA interactions available at the miRNet database (https://www.mirnet.ca/), and information about PRG-miRNA interactions available at the TargetScan database (https://www.targetscan.org/vert_72/).

Abbreviations

AFP:

Alpha fetoprotein

AFP-L3:

Lens culinaris agglutinin-reactive alpha-fetoprotein

AIC:

Akaike information criterion

Akt:

Protein kinase B

AUC:

Area under the curve

BH:

Benjamini–Hochberg

ceRNA:

Competing endogenous RNA

Cox:

Proportional hazards model

DCP:

Des-γ-carboxy prothrombin

DKK1:

Dickkopf WNT signaling pathway inhibitor 1

FDR:

False discovery rate

FOXD2-AS1:

FOXD2 adjacent opposite strand RNA 1

GSDMD:

Gasdermin D

GSEA:

Gene set enrichment analysis

GZMA:

Granzyme A

GZMB:

Granzyme B

GZMK:

Granzyme K

HCC:

Hepatocellular carcinoma

HOTTIP:

HOXA distal transcript antisense RNA

ImmuCellAI:

Immune Cell Abundance Identifier

LASSO:

Least absolute shrinkage and selection operator

LINC01554:

Long intergenic non-protein coding RNA 1554

lncRNA:

Long noncoding RNA

MEG3:

Maternally expressed 3

miR-744-5p:

MicroRNA-744-5p

MSigDB:

Molecular Signatures Database

mTOR:

Mammalian target of rapamycin

NK:

Natural killer

NLRP3:

NLR family pyrin domain containing 3

NRSN2-AS1:

NRSN2 antisense RNA 1

OS:

Overall survival

PKM2:

Pyruvate kinase M2

PRF1:

Perforin 1

PRG:

Pyroptosis-related gene

PRKX:

Protein Kinase CAMP-Dependent X-Linked Catalytic Subunit

PRlncRNA:

Pyroptosis-related lncRNA

REO:

Relative expression ordering

RNA-Seq:

Quantitative transcriptomics data

RPM:

Reads per million

TCGA:

The Cancer Genome Atlas

TPM:

Transcript per million

Treg:

Regulatory T cell

t-ROC:

Time-dependent receiver operating characteristics

Wnt:

Wingless/integrated

References

  1. Sachdeva M, Arora SK. Prognostic role of immune cells in hepatocellular carcinoma. EXCLI J. 2020. https://doi.org/10.17179/excli2020-1455.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Guo Z, Zhong N, Xu X, Zhang Y, Luo X, Zhu H, Zhang X, Wu D, Qiu Y, Tu F. Prediction of hepatocellular carcinoma response to transcatheter arterial chemoembolization: a real-world study based on non-contrast computed tomography radiomics and general image features. J Hepatocell Carcinoma. 2021. https://doi.org/10.2147/JHC.S316117.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Nie Y, Li J, Wu W, Guo D, Lei X, Zhang T, Wang Y, Mao Z, Zhang X, Song W. A novel nine-lncRNA risk signature correlates with immunotherapy in hepatocellular carcinoma. Front Oncol. 2021. https://doi.org/10.3389/fonc.2021.706915.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Yang Q, Xie B, Tang H, Meng W, Jia C, Zhang X, Zhang Y, Zhang J, Li H, Fu B. Minichromosome maintenance 3 promotes hepatocellular carcinoma radioresistance by activating the NF-κB pathway. J Exp Clin Cancer Res. 2019. https://doi.org/10.1186/s13046-019-1241-9.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Ruan S, Han C, Sheng Y, Wang J, Zhou X, Guan Q, Li W, Zhang C, Yang Y. Antcin A alleviates pyroptosis and inflammatory response in Kupffercells of non-alcoholic fatty liver disease by targeting NLRP3. Int Immunopharmacol. 2021. https://doi.org/10.1016/j.intimp.2021.108126.

    Article  PubMed  Google Scholar 

  6. Zhang Z, Zhang Y, Xia S, Kong Q, Li S, Liu X, Junqueira C, Meza-Sosa KF, Mok TMY, Ansara J, Sengupta S, Yao Y, Wu H, Lieberman J. Gasdermin E suppresses tumour growth by activating anti-tumour immunity. Nature. 2020. https://doi.org/10.1038/s41586-020-2071-9.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Khalifeh M, Penson PE, Banach M, Sahebkar A. Statins as anti-pyroptotic agents. Arch Med Sci. 2021. https://doi.org/10.5114/aoms/141155.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Zhao H, Wu L, Yan G, Chen Y, Zhou M, Wu Y, Li Y. Inflammation and tumor progression: signaling pathways and targeted intervention. Signal Transduct Target Ther. 2021. https://doi.org/10.1038/s41392-021-00658-5.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Li Y, Lv J, Shi W, Feng J, Liu M, Gan S, Wu H, Fan W, Shi M. Inflammasome signaling: a novel paradigm of hub platform in innate immunity for cancer immunology and immunotherapy. Front Immunol. 2021. https://doi.org/10.3389/fimmu.2021.710110.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Wang YY, Liu XL, Zhao R. Induction of pyroptosis and its implications in cancer management. Front Oncol. 2019. https://doi.org/10.3389/fonc.2019.00971.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Zhang X, Zhang P, An L, Sun N, Peng L, Tang W, Ma D, Chen J. Miltirone induces cell death in hepatocellular carcinoma cell through GSDME-dependent pyroptosis. Acta Pharm Sin B. 2020. https://doi.org/10.1016/j.apsb.2020.06.015.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Zhang Y, Yang H, Sun M, He T, Liu Y, Yang X, Shi X, Liu X. Alpinumisoflavone suppresses hepatocellular carcinoma cell growth and metastasis via NLRP3 inflammasome-mediated pyroptosis. Pharmacol Rep. 2020. https://doi.org/10.1007/s43440-020-00064-8.

    Article  PubMed  Google Scholar 

  13. Karger A, Nandigama R, Stenzinger A, Grimminger F, Pullamsetti SS, Seeger W, Savai R. Hidden treasures: macrophage long non-coding RNAs in lung cancer progression. Cancers (Basel). 2021. https://doi.org/10.3390/cancers13164127.

    Article  PubMed  Google Scholar 

  14. He D, Zheng J, Hu J, Chen J, Wei X. Long non-coding RNAs and pyroptosis. Clin Chim Acta. 2020. https://doi.org/10.1016/j.cca.2019.11.035.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Tan C, Liu W, Zheng ZH, Wan XG. LncRNA HOTTIP inhibits cell pyroptosis by targeting miR-148a-3p/AKT2 axis in ovarian cancer. Cell Biol Int. 2021. https://doi.org/10.1002/cbin.11588.

    Article  PubMed  Google Scholar 

  16. Yan H, Luo B, Wu X, Guan F, Yu X, Zhao L, Ke X, Wu J, Yuan J. Cisplatin induces pyroptosis via activation of MEG3/NLRP3/caspase-1/GSDMD pathway in triple-negative breast cancer. Int J Biol Sci. 2021. https://doi.org/10.7150/ijbs.60292.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Qu G, Wang D, Xu W, Guo W. Comprehensive analysis of the correlation between pyroptosis-related LncRNAs and tumor microenvironment, prognosis, and immune infiltration in hepatocellular carcinoma. Front Genet. 2022. https://doi.org/10.3389/fgene.2022.867627.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Liu ZK, Wu KF, Zhang RY, Kong LM, Shang RZ, Lv JJ, Li C, Lu M, Yong YL, Zhang C, Zheng NS, Li YH, Chen ZN, Bian H, Wei D. Pyroptosis-related LncRNA signature predicts prognosis and is associated with immune infiltration in hepatocellular carcinoma. Front Oncol. 2022. https://doi.org/10.3389/fonc.2022.794034.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Qi L, Chen L, Li Y, Qin Y, Pan R, Zhao W, Gu Y, Wang H, Wang R, Chen X, Guo Z. Critical limitations of prognostic signatures based on risk scores summarized from gene expression levels: a case study for resected stage I non-small-cell lung cancer. Brief Bioinform. 2016. https://doi.org/10.1093/bib/bbv064.

    Article  PubMed  Google Scholar 

  20. Kim S, Lin CW, Tseng GC. MetaKTSP: a meta-analytic top scoring pair method for robust cross-study validation of omics prediction analysis. Bioinformatics. 2016. https://doi.org/10.1093/bioinformatics/btw115.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Cao K, Liu M, Ma K, Jiang X, Ma J, Zhu J. Prediction of prognosis and immunotherapy response with a robust immune-related lncRNA pair signature in lung adenocarcinoma. Cancer Immunol Immunother. 2022. https://doi.org/10.1007/s00262-021-03069-1.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Guo Y, Jiang W, Ao L, Song K, Chen H, Guan Q, Gao Q, Cheng J, Liu H, Wang X, Guan G, Guo Z. A qualitative signature for predicting pathological response to neoadjuvant chemoradiation in locally advanced rectal cancers. Radiother Oncol. 2018. https://doi.org/10.1016/j.radonc.2018.01.010.

    Article  PubMed  Google Scholar 

  23. Zhou C, Wang S, Shen Z, Shen Y, Li Q, Shen Y, Huang J, Deng H, Ye D, Zhan G, Li J. Construction of an m6A-related lncRNA pair prognostic signature and prediction of the immune landscape in head and neck squamous cell carcinoma. J Clin Lab Anal. 2022. https://doi.org/10.1002/jcla.24113.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Tang R, Wu Z, Rong Z, Xu J, Wang W, Zhang B, Yu X, Shi S. Ferroptosis-related lncRNA pairs to predict the clinical outcome and molecular characteristics of pancreatic ductal adenocarcinoma. Brief Bioinform. 2022. https://doi.org/10.1093/bib/bbab388.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Guo Y, Yang PT, Wang ZW, Xu K, Kou WH, Luo H. Identification of three autophagy-related long non-coding RNAs as a novel head and neck squamous cell carcinoma prognostic signature. Front Oncol. 2021. https://doi.org/10.3389/fonc.2020.603864.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Wu F, Wei H, Liu G, Zhang Y. Bioinformatics profiling of five immune-related lncRNAs for a prognostic model of hepatocellular carcinoma. Front Oncol. 2021. https://doi.org/10.3389/fonc.2021.667904.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, Hoang CD, Diehn M, Alizadeh AA. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015. https://doi.org/10.1038/nmeth.3337.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Chang L, Zhou G, Soufan O, Xia J. miRNet 2.0: network-based visual analytics for miRNA functional analysis and systems biology. Nucleic Acids Res. 2020. https://doi.org/10.1093/nar/gkaa467.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Rawat M, Nighot M, Al-Sadi R, Gupta Y, Viszwapriya D, Yochum G, Koltun W, Ma TY. IL1B increases intestinal tight junction permeability by up-regulation of MIR200C-3p, which degrades occludin mRNA. Gastroenterology. 2020. https://doi.org/10.1053/j.gastro.2020.06.038.

    Article  PubMed  Google Scholar 

  30. McCormick SM, Heller NM. Commentary: IL-4 and IL-13 receptors and signaling. Cytokine. 2015. https://doi.org/10.1016/j.cyto.2015.05.023.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Arvanitakis K, Koletsa T, Mitroulis I, Germanidis G. Tumor-associated macrophages in hepatocellular carcinoma pathogenesis, prognosis and therapy. Cancers (Basel). 2022. https://doi.org/10.3390/cancers14010226.

    Article  PubMed  Google Scholar 

  32. Togashi Y, Shitara K, Nishikawa H. Regulatory T cells in cancer immunosuppression—implications for anticancer therapy. Nat Rev Clin Oncol. 2019. https://doi.org/10.1038/s41571-019-0175-7.

    Article  PubMed  Google Scholar 

  33. Niu X, Chen L, Li Y, Hu Z, He F. Ferroptosis, necroptosis, and pyroptosis in the tumor microenvironment: perspectives for immunotherapy of SCLC. Semin Cancer Biol. 2022. https://doi.org/10.1016/j.semcancer.2022.03.009.

    Article  PubMed  Google Scholar 

  34. Ayers M, Lunceford J, Nebozhyn M, Murphy E, Loboda A, Kaufman DR, Albright A, Cheng JD, Kang SP, Shankaran V, Piha-Paul SA, Yearley J, Seiwert TY, Ribas A, McClanahan TK. IFN-γ-related mRNA profile predicts clinical response to PD-1 blockade. J Clin Investig. 2017. https://doi.org/10.1172/JCI91190.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Miao YR, Zhang Q, Lei Q, Luo M, Xie GY, Wang H, Guo AY. ImmuCellAI: a unique method for comprehensive T-cell subsets abundance prediction and its application in cancer immunotherapy. Adv Sci (Weinh). 2020. https://doi.org/10.1002/advs.201902880.

    Article  PubMed  Google Scholar 

  36. Rizvi S, Wang J, El-Khoueiry AB. Liver cancer immunity. Hepatology. 2021. https://doi.org/10.1002/hep.31416.

    Article  PubMed  Google Scholar 

  37. Best J, Bechmann LP, Sowa JP, Sydor S, Dechêne A, Pflanz K, Bedreli S, Schotten C, Geier A, Berg T, Fischer J, Vogel A, Bantel H, Weinmann A, Schattenberg JM, Huber Y, Wege H, von Felden J, Schulze K, Bettinger D, Thimme R, Sinner F, Schütte K, Weiss KH, Toyoda H, Yasuda S, Kumada T, Berhane S, Wichert M, Heider D, Gerken G, Johnson P, Canbay A. GALAD score detects early hepatocellular carcinoma in an international cohort of patients with nonalcoholic steatohepatitis. Clin Gastroenterol Hepatol. 2020. https://doi.org/10.1016/j.cgh.2019.11.012.

    Article  PubMed  Google Scholar 

  38. Best J, Bilgi H, Heider D, Schotten C, Manka P, Bedreli S, Gorray M, Ertle J, van Grunsven LA, Dechêne A. The GALAD scoring algorithm based on AFP, AFP-L3, and DCP significantly improves detection of BCLC early stage hepatocellular carcinoma. Z Gastroenterol. 2016. https://doi.org/10.1055/s-0042-119529.

    Article  PubMed  Google Scholar 

  39. Berhane S, Toyoda H, Tada T, Kumada T, Kagebayashi C, Satomura S, Schweitzer N, Vogel A, Manns MP, Benckert J, Berg T, Ebker M, Best J, Dechêne A, Gerken G, Schlaak JF, Weinmann A, Wörns MA, Galle P, Yeo W, Mo F, Chan SL, Reeves H, Cox T, Johnson P. Role of the GALAD and BALAD-2 serologic models in diagnosis of hepatocellular carcinoma and prediction of survival in patients. Clin Gastroenterol Hepatol. 2016. https://doi.org/10.1016/j.cgh.2015.12.042.

    Article  PubMed  Google Scholar 

  40. Adhoute X, Penaranda G, Raoul JL, Le Treut P, Bollon E, Hardwigsen J, Castellani P, Perrier H, Bourlière M. Usefulness of staging systems and prognostic scores for hepatocellular carcinoma treatments. World J Hepatol. 2016. https://doi.org/10.4254/wjh.v8.i17.703.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Schütte K, Schulz C, Link A, Malfertheiner P. Current biomarkers for hepatocellular carcinoma: surveillance, diagnosis and prediction of prognosis. World J Hepatol. 2015. https://doi.org/10.4254/wjh.v7.i2.139.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Zhang Z, Shang J, Hu B, Shi H, Cao Y, Li J, Jiao T, Zhang W, Lu S. Prognosis and tumour immune microenvironment of patients with hepatocellular carcinoma by a novel pyroptosis-related lncRNA signature. Front Immunol. 2022. https://doi.org/10.3389/fimmu.2022.836576.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Zhang Z, Xia F, Xu Z, Peng J, Kang F, Li J, Zhang W, Hong Q. Identification and validation of a novel pyroptosis-related lncRNAs signature associated with prognosis and immune regulation of hepatocellular carcinoma. Sci Rep. 2022. https://doi.org/10.1038/s41598-022-13046-y.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Zheng YL, Li L, Jia YX, Zhang BZ, Li JC, Zhu YH, Li MQ, He JZ, Zeng TT, Ban XJ, Yuan YF, Li Y, Guan XY. LINC01554-mediated glucose metabolism reprogramming suppresses tumorigenicity in hepatocellular carcinoma via downregulating PKM2 expression and inhibiting Akt/mTOR signaling pathway. Theranostics. 2019. https://doi.org/10.7150/thno.28992.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Chen Q, Xie J, Yang Y. Long non-coding RNA NRSN2-AS1 facilitates tumorigenesis and progression of ovarian cancer via miR-744-5p/PRKX axis. Biol Reprod. 2022. https://doi.org/10.1093/biolre/ioab212.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Lei T, Zhu X, Zhu K, Jia F, Li S. EGR1-induced upregulation of lncRNA FOXD2-AS1 promotes the progression of hepatocellular carcinoma via epigenetically silencing DKK1 and activating Wnt/β-catenin signaling pathway. Cancer Biol Ther. 2019. https://doi.org/10.1080/15384047.2019.1595276.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors thank the staff members of TCGA, PubMed, GENCODE, MSigDB, miRNet, TargetScan, and ImmuCellAI databases.

Funding

This work was supported in part by the National Natural Science Foundation of China (Grant No. 61961002), and the Thousand Talents Program of Jiangxi for High-level talents in innovation and entrepreneurship (No. Jxsq2020101096).

Author information

Authors and Affiliations

Authors

Contributions

TW and HL designed the research. TW completed all analysis in the research. TW and HL drafted the manuscript. GH participated in the critical revision of the manuscript and provided support in English polishing work. NL, FL, ZC, LM, and TH participated in the manuscript's data collection and literal modification. HL supervised the whole study process. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Guini Hong or Hongdong Li.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Additional information

Publisher's Note

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

Supplementary Information

Additional file 1. Table S1.

The 40 pyroptosis-associated genes identified from the MSigDB database.

Additional file 2. Table S2.

Clinical characteristics of patients in training cohort and validation cohort.

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

Wu, T., Li, N., Luo, F. et al. Screening prognostic markers for hepatocellular carcinoma based on pyroptosis-related lncRNA pairs. BMC Bioinformatics 24, 176 (2023). https://doi.org/10.1186/s12859-023-05299-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12859-023-05299-9

Keywords