Identification and Verification of Key Genes in Colorectal Cancer Liver Metastases Through Analysis of Single-Cell Sequencing Data and TCGA Data

Abstract Background Colorectal cancer (CRC) is highly prevalent worldwide, with more patients experiencing colorectal cancer liver metastases (CRLM). This study aimed to identify key genes in CRLM through single-cell sequencing data reanalysis and experimental validation. Methods The study analyzed single-cell RNA-sequencing (scRNA-seq) data from the Gene Expression Omnibus (GEO) database. Gene ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) were used for gene functional enrichment analysis. The Cancer Genome Atlas (TCGA) data enabled bulk-RNA expression and survival prognosis analysis. Real-time polymerase chain reaction (qPCR) detected mRNA expression, whereas Western blot determined protein levels. Cell function experiments assessed SPARC’s impact on CRC cell behavior. Results Cluster analysis showed 23 classes among 17 CRLM samples, representing six cell types. A GO and KEGG analysis identified interleukin-1 beta (IL1B), CD2 molecule (CD2), and C-X-C motif chemokine ligand 8 (CXCL8) as significant prognostic factors in CRC. Secreted protein acidic and cysteine rich (SPARC) was one of the top differentially expressed genes (DEGs) in tissue stem cells, confirmed in primary and metastatic lesions. Metastatic lesions showed higher expression of SPARC and CRC stem cell marker leucine-rich repeat containing G protein-coupled receptor 5 (LGR5), which was significantly correlated positively with LGR5 expression. Knockdown of SPARC reduced CRC cell sphere- and colony-formation, invasion, and migration abilities. Overexpression of SPARC significantly increased the malignancy of CRC cells. Conclusions Several key genes were identified in the process of CRLM. In CRLM samples and those corresponding to CRC stem cells, SPARC was significantly upregulated. In the therapy of CRLM, SPARC might be a potential target.