Pyroptosis-related genes features on prediction of the prognosis in liver cancer: An integrated analysis of bulk and single-cell RNA sequencing

Abstract

ObjectiveThis study explores the impact of pyroptosis-related genes (PRG) on the prognosis of liver cancer (LC). Methods421 samples (371 tumor samples and 50 normal samples) from the Cancer Genome Atlas (TCGA) were included in this study. GSE14520 dataset (data of RNA expression and relevant clinicopathological features), GSE125449 dataset (single-cell data in LC) and HCCDB18 dataset (validation on the reliability of the model) were downloaded as appropriate. Download the PRG and its corresponding pathway information from the gene set enrichment analysis (GSEA) website. The consensus clustering was performed by ConsensusClusterPlus package. Differentially expressed genes (DEGs) were identified using limma package, and prognostic features were constructed using un/multivariate and Lasso Cox regression. Pathway enrichment analysis was conducted by ssGSEA method. Receiver Operating Characteristic and the survival analysis were conducted by timeROC and Survminer packages. The Seurat package was used for single-cell RNA sequencing (scRNA-seq) analysis. For cellular validation, following the quantification on the key genes via reverse-transcription quantitative PCR, the Transwell and scratch assays were applied to evaluate the in-vitro invasion and migration of LC cells Huh-7. Results12 prognosis-related genes were identified to be related to the progression of LC. Three subtypes including C1, C2 and C3 were categorized using the 12 prognosis-related genes and PRGs significantly related to the prognosis of LC patients. The worst and best prognosis was seen in C3 subtype and C2 subtype, respectively. Hallmark pathway enrichment analysis has shown the concurrent immunoactivation and immune escape in C3 subtype. A RiskScore model was constructed using 8 key genes (KPNA2, UCK2, FTCD, CBX2, RAB32, HMMR, S100A9 and ANXA10) from the DEGs of three subtypes. The RiskScore system as an independent prognostic factor dividing the patients into high and low risk groups, and patients of the high-risk group had poor prognosis in both test set and validation set. A nomogram model combining the risk score had the extreme higher benefit. Further, 6 subclusters were identified from scRNA-seq analysis, where the highest PYROPTOSIS score was seen in Monocytic-Macrophages. The quantification on the key genes has suggested the high expressions of KPNA2, UCK2, CBX2, RAB32, HMMR and S100A9 and the low expressions of FTCD and ANXA10 in LC cells Huh-7. Particularly, UCK2 knockdown evidently diminished the number of invaded and migrated LC cells in vitro. ConclusionThe risk model associated with pyproptosis is crucial for the tumor immunity of LC and may serve as a prognostic indicator for patients suffering from LC. Our findings will offer new perspectives for immunotherapies targeting LC.