Understanding the heterogeneity in liver hepatocellular carcinoma with a special focus on malignant cell through single-cell analysis

Abstract

Abstract Introduction Hepatocellular carcinoma (HCC) is the most common form of liver cancer globally and remains a major cause of cancer-related deaths. HCC exhibits significant intra-tumoral and interpatient heterogeneity, impacting treatment efficacy and patient prognosis. Methods We acquired transcriptome data from the TCGA and ICGC databases, as well as liver cancer chip data from the GEO database, and processed the data for subsequent analysis. We also obtained single cell data from the GEO database and performed data analysis using the Seurat package. To further investigate epithelial cell subgroups and their copy number variations, we used the Seurat workflow for subgroup classification and the InferCNV software for CNV analysis, utilizing endothelial cells as a reference. Pseudo-time analysis and transcription factor analysis of epithelial cells were performed using the monocle2 and SCENIC software, respectively. To assess intercellular communication, we employed the CellChat package to identify potential ligand-receptor interactions. We also analyzed gene expression differences and conducted enrichment analysis using the limma and clusterProfiler packages. Additionally, we established tumor-related risk characteristics using Cox analysis and Lasso regression, and predicted immunotherapy response using various datasets. Results The samples were classified into 23 clusters, with malignant epithelial cells being the majority. Trajectory analysis revealed the differentiation states of the malignant epithelial cells, with cluster 1 being in the terminal state. Functional analysis revealed higher aggressiveness and epithelial-mesenchymal transition (EMT) scores in cluster 1, indicating a higher propensity for metastasis. RBP4+ tumor cells were highly enriched with hypoxia process and intensive cell-to-cell communication. A prognostic model was established, and immune infiltration analysis showed increased infiltration in the high-risk group. TP53 demonstrated significant differences in mutation rate between the two risk groups. Validation analysis confirmed the up-regulation of model genes, including AKR1B10, ARL6IP4, ATP6V0B, and BSG in tumor tissues. Conclusion A prognostic model was established based on HCC malignant cell associated gene signature, displaying decent prognosis guiding effectiveness in the multiple cohorts. The study provided comprehensive insights into the heterogeneity and potential therapeutic targets of LIHC.