AUTHOR=Xiang Jiankang , Liu Chuan , He Qingmin , He Pengzhan , Dong Weiguo 

TITLE=Comprehensive analysis of immunogenic cell death associated genes expression, tumor microenvironment, and prognosis in hepatocellular carcinoma

JOURNAL=Frontiers in Pharmacology

VOLUME=14

YEAR=2023

URL=https://www.frontiersin.org/journals/pharmacology/articles/10.3389/fphar.2023.1122011

DOI=10.3389/fphar.2023.1122011

ISSN=1663-9812

ABSTRACT=<p><bold>Background:</bold> Immunogenic cell death (ICD) plays an important role in the development of cancers. This study attempted to explore the role of ICD in the prognosis of hepatocellular carcinoma (HCC).</p><p><bold>Methods:</bold> Gene expression and clinical data were downloaded from The Cancer Genome Alas and Gene Expression Omnibus dataset. The immune/stromal/Estimate scores of the tumor microenvironment (TME) were calculated by ESTIMATE and CIBERSORT algorithms. Kaplan-Meier analysis, functional enrichment analysis, least absolute shrinkage and selection operator (LASSO) analysis, and univariate and multivariate Cox regression analysis were used for prognostic gene screening and prognostic model construction. The correlation of immune cell infiltration and risk scores was analyzed as well. Molecular docking was used to explore the relevance of related genes to anti-cancer drugs.</p><p><bold>Results:</bold> Ten ICD associated differentially expressed genes in HCC were found, and all of them had good predictive ability for HCC. ICD gene high amount of expression group was associated with poor prognosis (<italic>p</italic> = 0.015). The TME, immune cell infiltration and gene expression were different between ICD high and low groups (all <italic>p</italic> &lt; 0.05). Six ICD associated genes (BAX, CASP8, IFNB1, LY96, NT5E and PIK3CA) which could predict the survival status were identified and used to construct the prognostic model for HCC. A risk score was calculated and it could be used as an independent prognostic factor in HCC patients (<italic>p</italic> &lt; 0.001). In addition, the risk score had a positive correlation with macrophage M0 (r = 0.33, <italic>p</italic> = 0.0086). Molecular docking indicated that sorafenib could bind strongly to the target protein, representing that sorafenib may exert anticancer effects through these six ICD associated genes.</p><p><bold>Conclusion:</bold> This study established a prognostic model including six ICD associated genes for HCC, which may deepen our understanding of ICD and guide therapy for HCC patients.</p>