Screening of Candidate Inflammatory Markers of Epithelial Cells in Hepatocellular Carcinoma based on Integration Analysis of TCGA/ICGC Databases and Single-cell Sequencing

Background: Hepatocellular Carcinoma (HCC) is closely linked to inflammatory reactions, with chronic liver diseases acting as major risk factors. In the inflammatory microenvironment, repeated damage and repair of liver cells lead to genetic mutations, abnormal proliferation, and tumorigenesis. Objective: This study aimed to investigate the expression profile of specific cell clusters under inflammatory stimulation in HCC and identify potential therapeutic drugs. Methods: Comprehensive analysis of HCC transcriptome data and single-cell sequencing data from TCGA, ICGC, and GEO databases was conducted to explore the specific molecular mechanisms of epithelial cells. Virtual screening of natural compounds in the ZINC database and in vitro cell experiments were performed to identify drugs that regulate the expression of inflammatory factors in epithelial cells. Results: Analysis of the single-cell dataset revealed cell clusters closely associated with HCC, notably Epithelial cells, Hepatocytes, MSC, and iPS cells, with Epithelial cells playing a pivotal role in HCC development. Further investigation of TCGA data unveiled 83 differentially expressed genes (DEGs) related to inflammatory responses in HCC. Intersection analysis of DEGs in epithelial cells and HCC DEGs identified 12 common DEGs, including ADRM1, ATP2B1, FZD5, GPC3, KIF1B, KLF6, LY6E, MET, NAMPT, SERPINE1, SPHK1, and SRI. Prognostic analysis revealed that CCL7, GPR132, ITGB8, PTAFR, SELL, and VIP were influential in the survival prognosis of HCC. A prognostic model based on the expression levels of these genes demonstrated an increased risk of HCC associated with higher differential expression of inflammatory response genes. Additionally, molecular dynamics simulations indicated that compounds NADH and Deferoxamine formed stable docking models with the inflammatory protein VIP, suggesting their potential as candidates for targeted therapy. Conclusion: Inflammatory factors CCL7, GPR132, ITGB8, PTAFR, SELL, and VIP influence the inflammatory cascade response in HCC epithelial cells, and their expression correlates with the survival prognosis of HCC patients. Interfering with VIP expression effectively suppresses proliferation, migration, and invasion of HCC cells, as well as inhibiting the occurrence of inflammatory cascade reactions, thus slowing down the progression of hepatocellular carcinoma. Furthermore, compounds NADH and Deferoxamine have the potential to target and bind to the inflammatory protein VIP, highlighting their relevance in potential HCC treatment.