Fibroblast Heterogeneity in Hepatocellular Carcinoma and Identification of Prognostic Markers Based on Single-cell Transcriptome Analysis

Background: HCC is a malignant tumor with high morbidity and mortality. Fibroblasts play a key role in the tumor microenvironment (TME). However, the transcriptional regulatory mechanisms of fibroblasts remained unclear in HCC.

Aim: The aim of this study was to explore the complex role of fibroblasts in hepatocellular carcinoma (HCC) and to reveal their transcriptional regulatory mechanisms.

Objective: The goal of this study was to discover potential prognostic markers for HCC by analyzing the genetic variations and differentiation process of fibroblasts.

Methods: Single-cell transcriptome data from the non-tumor liver site and primary tumor site of HCC were acquired from GSE149614, processed, and clustered using the Seurat pipeline. The inferCNV algorithm was applied to infer copy number variations (CNVs) in fibroblasts. Subsequently, the mechanism underlying the interaction between fibroblasts and other cells in the TME of HCC was analyzed using CellChat software. The trajectory of cellular differentiation of fibroblasts from normal state to malignant state was examined using Monocle 2. SCENIC analysis was performed to identify key transcription factors (TFs) in fibroblasts and assess their correlation with HCC prognosis. Finally, qRT-PCR and Transwell assays were carried out to analyze the mRNA expression and cell metastasis.

Results: We identified a total of nine different cell types (B cells, cycling cells, endothelial cells, epithelial cells, fibroblasts, hepatocytes, macrophages, plasma cells, and T cells) based on the single-cell transcriptomic data of HCC. Among them, fibroblasts were highly enriched at the primary tumor site, and their number increased with advanced stages. In addition, significant deletions were detected on chromosome 6p of fibroblasts, and genes in this region were remarkably enriched in pathways associated with antigen processing and presentation. Intercellular communication showed that epithelial cells regulated fibroblasts the most. The differentiation of fibroblasts was mainly accompanied by a transition from normal to malignant state. Importantly, CEBPD and FOSB, the TFs most associated with the putative timing of fibroblasts, were under-expressed in human hepatocytes and showed a significant correlation with HCC prognosis. Overexpressed CEBPD inhibited HCC cell migration and invasion.

Conclusion: In conclusion, our study revealed that fibroblast recruitment and differentiation, as well as copy number loss at chromosome 6p, were associated with a higher degree of malignancy and immune dysfunction in HCC. The current discoveries provided new insights into the clinical treatment and diagnosis of HCC.