Hepatocellular carcinoma (HCC) is a fatal primary malignancy characterized by high invasion and migration. We aimed to explore the underlying metastasis-related mechanism supporting the development of HCC.
The dataset of single cell RNA-seq (GSE149614) were collected for cell clustering by using the Seurat R package, the FindAllMarkers function was used to find the highly expression and defined the cell cluster. The WebGestaltR package was used for the GO and KEGG function analysis of shared genes, the Gene Set Enrichment Analysis (GSVA) was performed by clusterProfiler R package, the hTFtarget database was used to identify the crucial transcription factors (TFs), the Genomics of Drug Sensitivity in Cancer (GDSC) database was used for the drug sensitivity analysis. Finally, the overexpression and trans-well assay was used for gene function analysis.
We obtained 9 cell clusters from the scRNA-seq data, including the nature killer (NK)/T cells, Myeloid cells, Hepatocytes, Epithelial cells, Endothelial cells, Plasma B cells, Smooth muscle cells, B cells, Liver bud hepatic cells. Further cell ecological analysis indicated that the Hepatocytes and Endothelial cell cluster were closely related to the cancer metastasis. Subsequently, the NDUFA4L2-Hepatocyte, GTSE1-Hepatocyte, ENTPD1-Endothelial and NDUFA4L2-Endothelial were defined as metastasis-supporting cell clusters, in which the NDUFA4L2-Hepatocyte cells was closely related to angiogenesis, while the NDUFA4L2-Endothelial was related with the inflammatory response and complement response. The overexpression and trans-well assay displayed that NDUFA4L2 exhibited clearly metastasis-promoting role in HCC progression.