Comprehensive profiling of cell subsets of gastric cancer and liver metastasis based on single cell RNA-sequencing analysis

Background: Liver metastasis (Li) is one of the most common distant metastatic sites for gastric cancer. A deeper understanding of its mechanism of action from a bioinformatics perspective may provide new insights. Therefore, the aim of this study was to use single cell RNA sequencing (scRNA-seq) to evaluate cell subtypes and understand the molecular mechanism of gastric cancer Li. Methods: The scRNA-seq data GSE163558 of gastric cancer and Li were downloaded from Gene Expression Omnibus (GEO). Single cell data were analyzed by various R packages such as Seurat, CellChat, gene set variation analysis (GSVA), monocle, gene set enrichment analysis (GSEA), and survival analysis and the results were plotted by ggplot2, R4.1.1. Protein expression was confirmed by immunohistochemistry in an additional patient cohort. Results: The genes APOD, CXCL5, and JUN are involved in epithelial cell metastasis. The infiltration of cytotoxic CD8 T cells was more frequent in gastric primary tumors (PTs) than in Lis. IL7R high natural killer (NK) cells that had high TXNIP and RIPOR2 expression were located at the site of Li and corresponded to a favorable prognosis. NK cells with high TNFAIP3 expression were located at the PT site and corresponded to a poor prognosis. Furthermore, the gene expression of myeloid cells was different depending on their localization in the PT site or Li. MHC-I signaling pathway was found to be activated in the PT compared to MHC-II at the site of Li, as revealed by cell-cell interaction, and HLA-E-CD94/NKG2A of NK cells was only activated in the PT and not in the Li. Conclusions: The present study revealed the difference between Li and gastric PT by scRNA-seq, demonstrating the impact of partial gene expression on patient prognosis. Our study provides new insights into gastric cancer and Li.