Characterization of cancer-related fibroblasts (CAF) in hepatocellular carcinoma and construction of CAF-based risk signature based on single-cell RNA-seq and bulk RNA-seq data

• Published in: Frontiers in immunology 2022-09, Vol.13, p.1009789
• Main Authors: Yu, Lianghe, Shen, Ningjia, Shi, Yan, Shi, Xintong, Fu, Xiaohui, Li, Shuang, Zhu, Bin, Yu, Wenlong, Zhang, Yongjie
• Format: Article
• Language: English
• Subjects: Biomarkers, Tumor - metabolism; cancer-associated fibroblasts; Carcinoma, Hepatocellular - pathology; differentially expressed genes; Fibroblasts - metabolism; Gene Expression Regulation, Neoplastic; Humans; Immunology; immunotherapy; Kaplan-Meier Estimate; liver hepatocellular carcinoma; Liver Neoplasms - pathology; nomogram; Reproducibility of Results; RNA-Seq
• ISSN: 1664-3224; 1664-3224
• DOI: 10.3389/fimmu.2022.1009789

Cancer-associated fibroblasts (CAFs) are involved in tumor growth, angiogenesis, metastasis, and resistance to therapy. We sought to explore the CAFs characteristics in hepatocellular carcinoma (HCC) and establish a CAF-based risk signature for predicting the prognosis of HCC patients. The signal-cell RNA sequencing (scRNA-seq) data was obtained from the GEO database. Bulk RNA-seq data and microarray data of HCC were obtained from the TCGA and GEO databases respectively. Seurat R package was applied to process scRNA-seq data and identify CAF clusters according to the CAF markers. Differential expression analysis was performed to screen differentially expressed genes (DEGs) between normal and tumor samples in TCGA dataset. Then Pearson correlation analysis was used to determine the DEGs associated with CAF clusters, followed by the univariate Cox regression analysis to identify CAF-related prognostic genes. Lasso regression was implemented to construct a risk signature based on CAF-related prognostic genes. Finally, a nomogram model based on the risk signature and clinicopathological characteristics was developed. Based on scRNA-seq data, we identified 4 CAF clusters in HCC, 3 of which were associated with prognosis in HCC. A total of 423 genes were identified from 2811 DEGs to be significantly correlated with CAF clusters, and were narrowed down to generate a risk signature with 6 genes. These six genes were primarily connected with 39 pathways, such as angiogenesis, apoptosis, and hypoxia. Meanwhile, the risk signature was significantly associated with stromal and immune scores, as well as some immune cells. Multivariate analysis revealed that risk signature was an independent prognostic factor for HCC, and its value in predicting immunotherapeutic outcomes was confirmed. A novel nomogram integrating the stage and CAF-based risk signature was constructed, which exhibited favorable predictability and reliability in the prognosis prediction of HCC. CAF-based risk signatures can effectively predict the prognosis of HCC, and comprehensive characterization of the CAF signature of HCC may help to interpret the response of HCC to immunotherapy and provide new strategies for cancer treatment.