Novel cell subtypes of SPP1 + S100P+, MS4A1-SPP1 + S100P+ were key subpopulations in intrahepatic cholangiocarcinoma

In this study, we integrated single-cell RNA sequencing (scRNA-seq) data to investigate cell heterogeneity and utilized MSigDB and CIBERSORTx to explore the pathways of major cell types and the relationships between different cell subtypes. Subsequently, we explored the correlation of cell subtypes with survival and used Gene Set Enrichment Analysis (GSEA) analyses to assess the pathways associated with the infiltration of specific cell subtypes. Finally, multiplex immunohistochemistry in tissue microarray cohort were performed to validate differences in protein level and their correlation with survival. iCCA presented a unique immune ecosystem, with increased proportions of Epi (epithelial)-SPP1–2, Epi-S100P-1, Epi-DN (double negative for SPP1 and S100P expression)-1, Epi-DN-2, Epi-DP (double positive for SPP1 and S100P expression)-1, Plasma B-3, Plasma B-2, B-HSPA1A-1, B-HSPA1A-2 cells, and decreased proportions of B-MS4A1. High level of Epi-DN-2, Epi-SPP1–1, Epi-SPP1–2, B-MS4A1, and low level of Epi-DB-1, Epi-S100P-1, and Epi-S100P-2 was significantly associated with longer overall survival (OS), and high level of B-MS4A1_Low_Epi-DN-2_Low was associated with the shortest OS. Moreover, the results of MsigDB and GSEA suggest that bile acid metabolism is a crucial process in iCCA. Finally, we found that S100P+, SPP1+, SPP1 + S100P+, and MS4A1-SPP1 + S100P+ were highly expressed, whereas MS4A1 was lowly expressed in iCCA, and patients with high level of S100P+, SPP1 + S100P+, and MS4A1-SPP1 + S100P+ exhibited shorter survival. We identified the cell heterogeneity of iCCA, found that iCCA is a unique immune ecosystem with many cell subtypes, and showed that the novel cell subtypes of SPP1 + S100P+ and MS4A1-SPP1 + S100P+ were key subpopulations in iCCA.