Abstract
Colorectal cancer (CRC) is a common malignant tumor with a high mortality rate worldwide. Advanced CRC often leads to liver metastasis, which has a poor prognosis, highlighting the need to investigate the underlying mechanisms. Omics, encompassing genomics, epigenomics, transcriptomics, proteomics, metabolomics, and microbiomics, enables comprehensive molecular analysis of cells and tissues. Tumor-omics research has advanced rapidly, with growing attention on CRC-related omics. However, systematic reviews on omics research specific to colorectal cancer liver metastasis (CRLM) are limited. This review summarizes the current status and progress of multi-omics research on CRLM and discusses the application of multi-omics technologies in basic research and the significant clinical implications.
keywords
- Colorectal cancer liver metastasis
- genomics
- transcriptomics
- epigenomics
- proteomics
- metabolomics
- microbiomics
Introduction
Global Cancer Observatory (GLOBOCAN) data from 2022 report 1.9 million colorectal cancer (CRC) cases and 904,000 deaths globally, ranking 3rd in incidence and 2nd in mortality worldwide1. Between 25% and 30% of patients with CRC develop liver metastases with a poor prognosis and a 5-year survival rate < 10%2,3. Studies on isolated cellular pathways and tissue metabolism offer limited insight into the pathogenesis of colorectal cancer liver metastasis (CRLM), diagnostic accuracy, recurrence, metastasis, therapeutic efficacy, and patient survival. Therefore, omics research significantly contributes to advancing our understanding of liver metastasis in patients with CRC.
Omics analysis involves statistical processing, and comparative and correlation analyses of bulk data across various molecular levels, such as genomics, epigenomics, transcriptomics, proteomics, metabolomics, and microbiomics. Omics analysis establishes relationships among molecular data and reveals the interconnections4. For example, multi-omics research on CRLM integrates omics data to identify key genes, metabolites, proteins, and microorganisms involved in the process for further experimental analysis and potential applications.
Although there are numerous reports on omics research in CRC, systematic reviews on omics in liver metastasis are lacking. This study provides a comprehensive review of genomics, epigenomics, transcriptomics, proteomics, metabolomics, microbiomics, and multi-omics CRLM studies and explores the current state of omics research in this context (Figure 1).
Omics framework of colorectal cancer liver metastasis. Genomic, epigenetic, transcriptomic, proteomic, metabolomic, and microbiomic analyses were performed on tissue from primary colorectal cancer lesions (T), liver metastases (LM), lymph node metastases (NM), blood circulation (BC), and fecal microbiota (FM) samples.
Genomic insights into CRLM
Genomics explores the structure, dynamics, and functions of the genome. Ki et al.5 utilized complementary (c)DNA microarrays to identify 46 genes potentially involved in CRLM progression. Similarly, Yamasaki et al.6 showed that the potential for liver metastasis in CRC is encoded within the primary tumor and detected by gene expression profiling. Studies in China identified frequent mutations in patients with CRC with liver metastasis, including tumor protein (TP) 53 (82%), adenomatous polyposis coli (APC) (76%), Kirsten rat sarcoma viral oncogene homolog (KRAS) (42%), SMAD family member (SMAD) 4 (14%), filaggrin (FLG) (13%), and f-box and wd repeat domain-containing (FBXW) 7 (11%)7. These findings link high-frequency mutations to poor patient prognosis, thereby highlighting the clinical implications of genomic research for CRLM treatment and prognosis.
Genomic differences between primary CRC lesions and liver metastasis
Although significant research has focused on the genomics of primary CRC tumors, studies involving liver metastasis are limited. Evidence suggests that differential gene expression between primary and metastatic sites is a key genomic feature of liver metastasis. Lin et al.8 performed gene ontology (GO) analysis and showed that genes associated with tissue remodeling and immune responses are upregulated in metastatic tumors, whereas genes related to proliferation and oxidative phosphorylation are downregulated. The study also identified 10 genes with higher expression in CRLM, 5 of which encoded “matrix proteins” involved in extracellular matrix (ECM) interactions, indicating a potential association with liver metastasis9. Notably, genomic profiles in patients with CRC vary significantly by metastasis site, with greater genomic divergence between primary tumors and liver metastases than between primary tumors and lung metastases10. Comparative genomic hybridization (CGH) and other techniques demonstrated a greater number of mutations in liver metastases than in primary tumors, which likely reflects increased genetic instability11. Additionally, patients with sporadic CRLM exhibited specific transcriptional dysregulation in liver metastases, whereas primary tumor expression remained normal. This finding indicates increased genomic instability in metastatic cells or adaptive changes to the liver microenvironment12,13.
Application of genomics in the tumor microenvironment (TME) of CRLM
Traditional sequencing methods treat tumor samples as a whole, whereas single-cell technologies differentiate cellular components, particularly immune cells. This shift has redirected cancer genomics from focusing solely on cancer cells to exploring the TME. Genomic analyses have revealed that immune cell changes within the TME may contribute to tumor progression and liver metastasis. Jackstadt et al.14 found that intestinal epithelial neurogenic locus notch homolog protein (NOTCH) 1 signaling drives metastasis through transforming growth factor (TGF)-dependent neutrophil recruitment, highlighting the role of NOTCH1 in CRC progression. Sathe et al.15 identified distinct intercellular communication programs in which macrophages and fibroblasts regulate the pre-metastatic niche in the liver, thereby forming a “soil” for tumor seeding. Specifically, TME-associated macrophages expressing secreted phosphoprotein (SPP) 1 promote primary CRC progression. Furthermore, Rong et al.16 used single-cell RNA sequencing to identify matrix gla protein (MGP) as a driver of CD8+ T cell exhaustion via NF-κB activation for CRLM. MGP knockout combined with anti-PD1 therapy synergistically resist liver metastasis. These findings highlight the role of genomics in understanding CRLM mechanisms and treatment outcomes.
Application of genomics in imaging diagnosis of CRLM
Studies linking genomics and imaging have proven valuable in diagnosing CRLM. Using next-generation sequencing (NGS) technology, a study analyzed weak target tumor enhancement (TTE) and normal tissues using CRC magnetic resonance imaging (MRI) and revealed correlations between TTE and genetic mutations in CRLM17. Wang et al.18 identified six genes (NRAS, KRAS, BRAF, PIK3CA, MAPK1, and STAT1) linked to specific imaging features, such as shape-elongation and gradient, thus providing a foundation for non-invasive diagnosis of liver metastasis.
Application of genomics in the treatment of CRLM
Surgery remains a viable option for patients with CRLM who meet resectability criteria. However, tumor molecular behavior, rather than surgical technique, ultimately determines patient prognosis19. Chan et al.20 reported no statistically significant differences in patient outcomes across surgical approaches. Chan et al.20 emphasized the importance of integrated management strategies based on cancer genomics to assess prognosis. Genomic analysis can guide patient selection, surgical timing, and technique for metastatic lesions21. Sun et al.22 demonstrated that P-cadherin regulates CRLM. This finding suggested P-cadherin inhibition as a potential therapeutic strategy. Liu et al.23 identified ankyrin repeat domain 42 (ANKRD42) as a key regulatory gene promoting CRLM through genome-wide sgRNA library screening. Knockdown of ANKRD42 significantly attenuates tumor cell migration and invasion, highlighting its potential as a therapeutic target. Additionally, combining oxaliplatin with mTOR inhibitors shows synergistic effects, highlighting the potential of integrating chemotherapy with targeted therapy24. While genomic analysis offers significant insights for CRC treatment, further research is warranted to advance clinical application.
Application of genomics in prognosis of CRLM
Genomic research can reveal the clonality and timing of metastatic CRC25. Recent studies have focused on identifying biomarkers in bodily fluids, such as blood and urine, including serum proteins (peptides), plasma miRNAs, circulating tumor cells, and nucleic acids. When combined with technologies, such as whole-genome sequencing and circulating tumor cell DNA analysis, genomic analysis can help predict CRLM occurrence. Notably, undetectable circulating tumor DNA following liver metastasis resection suggests a better prognosis for recurrence-free survival26.
The core of tumor genomics is mapping cancer-driving factors and providing personalized therapeutic strategies27. Genomics is crucial in understanding CRLM origins, molecular differences between primary and metastatic lesions, genetic profiling, treatment decisions, and survival prognosis.
Epigenomics of CRLM
Although every human body cell contains the same genome, distinct cellular phenotypes arise from different gene expression patterns, driven by cell type-specific transcription factors and epigenome structures28. Epigenomics primarily focuses on DNA and histone modifications that regulate chromatin structure, gene activity, and gene expression without altering the DNA sequence4. Studies have suggested that epigenetic regulation may contribute to tumor metastasis29, and DNA methylation and histone modifications are reversible. Investigating epigenetics in CRLM could reveal new therapeutic targets for CRLM diagnosis and treatment (Figure 2).
Epigenetic regulation in colorectal cancer (CRC) and CRC liver metastasis (CRLM). (A) Histone modifications as potential diagnostic biomarkers for CRC. Trimethylation of H3K9 and H3K27 enhances the expression of immune checkpoint genes, such as CTLA-4, TIGIT, PD-1, and TIM-3, which may serve as epigenetic markers for CRC. Methylation signatures on circulating nucleosomal histones, including H3K27me3 and H4K20me3, are detectable in the bloodstream and have potential diagnostic value. Differences in HIST2H3AK19Ac and H2BLK121Ac acetylation between primary CRC tumors and liver metastases may serve as biomarkers of CRC progression. (B) Schematic representation of histone modifications. Histone methylation, acetylation, and phosphorylation are dynamically regulated by histone methyltransferases (HMTs), histone demethylases (HDMs), histone acetyltransferases (HATs), histone deacetylases (HDACs), protein kinases (PKs), and protein phosphatases (PPs), resulting in chromatin remodeling toward either an open (active) or closed (repressive) state. (C) Role of DNA methylation in CRLM. Aberrant DNA methylation alters gene expression by promoting either transcriptional activation or repression, contributing to the progression of liver metastasis in CRC. (D) Schematic of epigenetic regulation via chromatin remodeling. Epigenetic mechanisms, including histone modifications, DNA methylation, and non-coding RNAs (ncRNAs), mediate chromatin structure reorganization, thereby influencing gene accessibility and transcriptional activity. B3GNT7, beta-1,3-N-acetylglucosaminyltransferase 7; CAB39, calcium-binding protein 39; C2, complement component 2; CEP112, centrosomal protein 112; CFLAR, CASP8 and FADD-like apoptosis regulator; CTLA-4, cytotoxic T-lymphocyte-associated protein 4; CTSC, cathepsin C; EVI2B, ecotropic viral integration site 2B; GUF1, GUF1 homolog, GTPase; KLF4, Kruppel-like factor 4; LZTS1, leucine zipper tumor suppressor 1; MAPT, microtubule-associated protein tau; MUC2, mucin 2; PAX8, paired box 8; PD-1, programmed cell death protein 1; PTH1R, parathyroid hormone 1 receptor; RASSF1A, ras association domain family member 1 isoform A; SHANK2, SH3 and multiple ankyrin repeat domains 2; SLC10A1, solute carrier family 10 member 1; THBS1, thrombospondin 1; TIGIT, T-cell immunoreceptor with Ig and ITIM domains; TIM3, T-cell immunoglobulin and mucin-domain containing-3; TNNI2, troponin I2, fast skeletal type; TRAPPC3, trafficking protein particle complex subunit 3; UPK3A, uroplakin 3A.
Relationship between DNA methylation and the occurrence, development, treatment, and prognosis of CRLM
DNA methylation occurs at cytosine residues preceding guanine in humans, forming CpG dinucleotides (C-phosphate-G bonds)30. While most CpG dinucleotides are methylated, CpG islands, found in the promoter regions of 40%–60% of tumor suppressor genes, remain unmethylated in healthy cells. These islands, spanning 200–2,000 base pairs with a CG content ≥ 50%, contribute to gene expression regulation30–32. Epigenetic mutations silence tumor suppressor genes or activate oncogenes in human cancers.
The methylation of various tumor suppressor genes is associated with CRLM. Zhang et al.33 noted that long non-coding RNA (lnc-LALC) promotes CRLM by epigenetically silencing leucine zipper tumor suppressor gene 1 (LZTS1). Upregulation of genes, such as SLC10A1, MAPT, SHANK2, PTH1R, and C2, and downregulation of genes, including CAB39, CFLAR, CTSC, THBS1, and TRAPPC3, are linked to DNA methylation and may act as biomarkers for CRLM34. Differential DNA methylation in primary CRC can predict risk of liver metastasis, with several regions, including TNNI2, PAX8, GUF1, KLF4, EVI2B, CEP112, and lncRNA AC011298, identified as prognostic biomarkers35,36.
Numerous studies have observed hypermethylation of genes in metastatic liver lesions of CRC. Ebert et al.37 found frequent hypermethylation of TPEF, which encodes a transmembrane protein with epidermal growth factor and follistatin domains, in primary CRC lesions and liver metastases. This finding suggested that epigenetic alterations are common in the early and late stages of carcinogenesis. Therefore, TPEF methylation may act as a diagnostic target. Several studies have identified different gene methylation mechanisms associated with CRLM. ALX4 is hypermethylated in primary colorectal adenocarcinomas and liver metastases38, and MUC2 promoter methylation is linked to liver and lymph node metastases39. Lu et al.40 showed that B3GNT7 methylation enhances CRC cell metastasis to the liver. Tommasi et al.41 found that RASSF1A methylation contributes to liver metastasis progression. Ju et al.42 suggested that DNA methylation-induced silencing of UPK3A enhances liver metastasis CRLM. Rodger et al.43 observed differential CpG methylation patterns, with hypomethylation in normal colon and primary CRC, but hypermethylation in liver metastases, which correlated with methylation levels and chromatin activity44.
DNA methylation is a valuable tool for the prognostic assessment of CRLM. For example, multiplex DNA methylation profiling of circulating free DNA (cfDNA) enables patient stratification and treatment monitoring in advanced CRC with liver metastasis45. Furthermore, methylation profiling has been used to identify effective methylation biomarkers, such as the significant increase in trichohyaline (TCHH) gene methylation, which correlates strongly with liver metastasis tumor volume. This finding suggested the potential of TCHH as a biomarker for liver metastasis46.
Epigenetic regulation is reversible and restoring normal epigenomic profiles could serve as a potential therapeutic target for cancer. CXCL14, a tumor suppressor gene in CRC, is activated but silenced during progression, primarily because of promoter hypermethylation. Therefore, reversing CXCL14 hypermethylation may offer an epigenetic therapeutic strategy for CRC treatment47. However, research on epigenetic therapies to intervene in CRLM is limited and further exploration is needed.
Relationship between histone modifications and the occurrence, development, treatment, and prognosis of CRLM
DNA wraps around nucleosome particles to form the organizational foundation for eukaryotic genomes. Core histones (H2A, H2B, H3, and H4) undergo post-translational modifications that help define specific chromatin regions and states48. These modifications include methylation, acetylation, ubiquitination, and phosphorylation with methylation and acetylation being most common49. CRLM has been reported to be closely linked to histone modifications50.
Histone methylation and demethylation in CRLM
Histone methylation involves adding methyl groups to lysine (K) residues on H3 and H4, and is a key post-translational modification51. Regulation of histone methylation is controlled by two enzyme groups [histone methyltransferases (HMTs), which add methylation marks, and histone demethylases (HDMs), which remove methylation marks]52,53. Each lysine residue can undergo three distinct methylation states (monomethyl, dimethyl, or trimethyl)54.
Histone methylation modifications are reportedly involved in CRC initiation and progression55–57. However, research on the role of histone methylation modifications in liver metastasis is limited. Tamagawa et al.50 found that liver metastasis in CRC is associated with histone modifications with H3K27me2 and H3K4me2 methylation levels potentially serving as independent prognostic factors58. Specifically, lower H3K4me2 levels correlate with poorer survival outcomes50. Further research is required to determine the role of histone methylation in CRLM.
Additionally, histone methylation modifications can serve as diagnostic markers for CRC59 and predictors of chemotherapy sensitivity60. Nair et al.61 compared the epigenetic modification mechanisms of immune checkpoints in CRC tissues and normal intestinal mucosa. Nair et al.61 found high expression of PD-1, CTLA-4, TIM-3, TIGIT, PD-L1, and Galectin-9 in CRC tissues. Nair et al.61 also observed a significant reduction in H3K9me3 binding to the PD-1 and TIGIT promoters and H3K27me3 binding to the CTLA-4 promoter in CRC tissues, as well as lower binding of both repressive histones to the TIM-3 promoter. These findings suggested that H3K9me3 and H3K27me3 are involved in upregulating CTLA-4, TIGIT, and PD-1, making these modifications potential diagnostic biomarkers for CRC.
Along with the gene methylation detected in tissues, histone methylation of circulating nucleosomes (H3K27me3 and H4K20me3) can serve as a biomarker for CRC diagnosis59. Kang et al.60 found that 5-fluorouracil (5-FU)-resistant colon cancer cells exhibit oxidative stress, which induces histone methylation-related proteins and DNA demethylases, thereby contributing to 5-FU resistance. However, the role of histone methylation modifications in predicting liver metastasis and survival prognosis remains unexplored.
Histone acetylation and deacetylation in CRLM
Histone acetylation is regulated by a balance between histone acetyltransferases (HATs) and histone deacetylases (HDACs)62. Shen et al.63 compared acetylation patterns in metastatic liver tumors and primary CRC tissues, and identified 31 downregulated acetylation sites on 22 proteins and 40 upregulated acetylation sites on 32 proteins. The most significantly altered acetylated histones were HIST2H3AK19Ac and H2BLK121Ac, which provides new insights into CRC metastatis. HDAC6-dependent deacetylation of A-kinase-anchoring protein 12 (AKAP12) leads to its ubiquitin-mediated degradation and promotes colon cancer metastasis64. These studies suggest that histone acetylation and deacetylation are involved in CRC metastasis. However, liver metastasis-specific research has been limited.
CRLM transcriptomics
Although each distinct cell type shares an identical genome in multicellular organisms, not all genes are transcriptionally active. Different cell types exhibit unique gene expression profiles. Transcriptomics classifies all types of transcripts, including mRNAs, ncRNAs, and miRNAs. The transcriptional structure of genes is determined by factors, such as the transcriptional start sites, 5′ and 3′ ends, splicing patterns, and other post-transcriptional modifications. Additionally, transcriptomics quantifies changes in the level of each transcript expression during development and under various conditions65.
Transcriptomic study on the occurrence and development of CRLM
Researchers have elucidated the mechanisms driving CRLM based on transcriptomic analysis, which has provided insights into potential diagnostic biomarkers and therapeutic targets. Sha et al.66 demonstrated that single-cell transcriptomic analysis can identify the potential cellular origin and driving factors of CRLM. Additionally, whole exome and RNA sequencing of CRC tumors and liver metastases revealed that copy number variations primarily follow several evolutionary patterns, such as clonal–clonal evolution 20q amplification (amp), 17p deletion (del), 18q del, 8p del, subclonal–clonal evolution (8q amp, 13q amp, and 8p del), and metastasis-specific evolution (8q amp).
These findings suggested a novel evolutionary process for CRLM67 and offer deeper insights into the mechanisms underlying liver metastasis through transcriptomics. The liver-specific transcription factors, FOXA2 and HNF1A, can reportedly bind to acquired enhancers and activate the transcription of liver-specific genes, thereby driving CRLM68. Cai et al.69 reported that polyclonal dissemination and transcriptional repressor GATA binding 1 (TRPS1) gene mutations may be key drivers of CRLM through single-cell exome sequencing. Hur et al.70 concluded that upregulated miR-885 and miR-122 and downregulated miR-10b in CRLM tissues suggests the potential as metastasis-specific biomarkers71. Furthermore, high expression of miR-20a is closely associated with the liver metastatic state of CRC72. Xu et al.73 found that circRNA_0001178 and circRNA_0000826 are significantly upregulated in CRLM tissues, suggesting that these circRNAs may serve as potential biomarkers for predicting CRLM.
Hematogenous dissemination is the primary route for CRLM, followed by two histopathologic growth patterns [sprouting angiogenesis and vessel co-option (VCO)]. Transcriptomic research has revealed distinct mechanisms for both metastatic patterns. Most liver metastases promote hematogenous spread through angiogenesis. Fleischer et al.74 found specific metabolic changes and WNT signaling activation in VCO-related metastatic cancer cells, providing potential therapeutic targets for VCO-based CRLM. Tumor pericytes (TPCs) play a pivotal role in this invasion. Li et al.75 used single-cell RNA sequencing (scRNA-seq) to identify 13 TPC subpopulations and a novel TCF21high “matrix–pericytes” phenotype, which promotes extracellular matrix stiffness, collagen remodeling, and metastasis. Rad et al.76 found that runt-related transcription factor-1 (RUNX1) overexpression in replacement lesions drives cancer cell motility in VCO-type CRLM.
Transcriptional differences between primary CRC lesions and liver metastases
Various researchers have compared the transcriptomic profiles of primary CRC lesions and liver metastases to determine transcriptional heterogeneity and cellular diversity. Liver metastases exhibit lower expression of miR-99b-5p, miR-377, miR-143, miR-10b, miR-28-5p, and miR-200c compared to primary tumor tissues, while miR-196b-5p, miR-122, miR-122*, and miR-885-5p are upregulated71,77. Four lncRNAs (GAS5, H19, MEG3, and Yiya) also exhibit significant differences78. Wang et al.79 identified increased CD8_CXCL13 and CD4_CXCL13 subpopulations in liver metastases, which are linked to better prognosis. Wang et al.79 also observed that distinct MCAM(+) fibroblasts in liver metastases may activate CD8 + CXCL13 cells through Notch signaling. Other studies showed that intestinal epithelial Notch1 signaling drives metastasis formation in genomic research14. This finding suggested a dual role of the Notch pathway in CRC progression and patient survival by regulating immune cells. These studies highlight gene expression differences between primary and metastatic CRC lesions.
Application of transcriptomics in CRLM treatment
Researchers have conducted mRNA expression profiling studies to analyze the correlation between gene expression and therapeutic strategies80. Simultaneously, differences in drug sensitivity can be detected at the transcriptomic level, potentially leading to the development of therapeutic agents based on variations in gene expression81. Transcriptomic sequencing of endothelial cells from CRLM and adjacent normal liver tissues revealed molecular differences, one of which was upregulation of the checkpoint kinase, WEE1. Subsequent experiments confirmed that WEE1 inhibitors could be promising therapeutic options for the treatment of CRLM82. Using single-cell transcriptomic analysis, Zhang et al.83 found that secreted protein acidic and rich in cysteine (SPARC) was significantly upregulated in CRLM. SPARC knockdown reduced CRC cell spheroid and colony formation, invasion, and migration. Sun et al.84 analyzed transcriptomic data from primary adjacent tissues, colorectal tumors, and metastases in patients, which revealed that SLITRK4 inhibition suppressed tumor growth and liver metastasis. Pharmacologic inhibition using lipid-polymer hybrid nanoparticles (NPs) to deliver siRNA effectively inhibited CRLM. A transcriptomic analysis is valuable for developing new therapeutic approaches for patients with CRC, with and without liver metastasis. However, current research on the application of transcriptomics in the treatment of liver metastasis is limited.
Application of transcriptomics in the prognostic analysis of CRLM
Transcriptomics and other methods have identified factors that are negatively correlated with CRLM prognosis, as follows:
high exosomal CXCL10 expression is linked to poor prognosis and shorter progression-free survival85;
overexpression of lncRNA-SNHG15 is associated with lower survival rates in patients with CRLM86;
high miR-221 or miR-222 expression in cancer stroma correlates with liver metastasis, distant metastasis, and short overall survival87;
miR-200c and miR-196b-5p expression are associated with short overall survival in patients with CRLM77;
miR-125, miR-127, miR-145, miR-192, miR-194, miR-199-5p, miR-215, and miR-429 expression is associated with poor survival in liver metastases88;
increased miR-210 and miR-133b expression correlates with reduced survival89; and
high CIB1 expression is associated with poor progression-free survival and OS90.
Additionally, transcriptomic studies have identified the following factors that correlate with CRLM prognosis:
Proteomic studies in CRLM
The concept of a proteome was first proposed by Wilkins et al.92 in 1995. A proteome refers to the entire set of proteins expressed in the genome of a cell, tissue, or organism with the respective forms and functional states93. This chapter analyzes the occurrence and progression of liver metastasis in CRC from a proteomic perspective, highlighting the differences between primary colorectal tumors and metastatic liver lesions. We also discuss the implications of proteomic findings in guiding the diagnosis, treatment, and prognosis of CRC.
Proteomic studies on CRLM occurrence and progression
The development and progression of CRC are closely associated with the ECM, which is composed of various proteins, including collagen, proteoglycans, and glycoproteins. Proteomic analyses have identified upregulated collagen types (COL10A1, COL12A1, COL14A1, and COL15A1) in CRLM, which are also present in primary tumors, suggesting ECM alterations involved in tumor proliferation and liver metastasis94. In a study involving 3 recurrent liver metastases from 1 patient, Voß et al.95 identified 56 upregulated ECM-related proteins in the 3rd metastatic lesion, 26 of which were linked to poor prognosis in CRC. This finding indicated varying malignancy levels in liver metastases.
Epithelial-mesenchymal transition (EMT) is a prerequisite for metastasis. Thrombospondin-1 (THBS1) was shown to promote CRLM by facilitating EMT based on a proteomic analysis of primary CRC specimens96. Unstable epithelial/mesenchymal (E/M)-type tumor cells exhibit the strongest potential for liver metastasis. The proportion of E/M-type circulating tumor cells is associated with distant metastasis. Furthermore, EMT-related protein clusters (ERBB2, COL6A1, and CAVIN1) are likely implicated in EMT-mediated metastasis97.
Other studies utilizing proteomic approaches have identified key proteins associated with CRC and liver metastasis, including heat shock protein family D member 1 (HSPD1), eukaryotic translation elongation factor 1 gamma, heterogeneous nuclear ribonucleoprotein A2/B1, fibrinogen beta chain (FGB), Talin 1, adaptor related protein complex 2 subunit alpha 2, serrate RNA effector molecule homolog, apolipoprotein C3, phosphoglucomutase 598, filamin A-interacting protein 1-like (FILIP1L), and plasminogen (PLG)99, as well as overexpression of apolipoprotein E100. However, the relationship between these proteins and CRLM requires further validation.
Overall, studies employing proteomic approaches have unveiled the mechanisms underlying CRLM, which are important for advancing our understanding of the initiation and progression of liver metastasis.
Proteomic differences between CRLM and primary tumors
Retrospective studies have revealed molecular differences between primary and metastatic CRC lesions, including the CAM-EMT, RAS-RAF-MEK, PI3K-AKT-PTEN, TGFBRII/SMAD4, apoptotic, and transcriptional error regulation pathways. These findings suggest that overall protein marker expression profiles may exhibit significant differences101. Farhner et al.102 compared the proteomes of seven specimens from primary and liver metastatic CRC lesions and revealed several upregulated proteins in liver metastasis, particularly those involved in glucose metabolism, including pyruvate carboxylase, fructose-bisphosphate aldolase B, and fructose-1,6-bisphosphatase 1. Moreover, Kim et al.103 performed a proteomic analysis of primary and liver metastatic CRC lesions. Several proteins were significantly upregulated in liver metastases, including SERPINA1, APOA1, ITLN1, DES, DBI, SDHA, and CA1 1, as well as arginase and glutathione S-transferase A3104. Fahrner et al.102 discovered that liver metastases exhibit a loss of several structural proteins associated with muscle contraction and cell junction assembly, such as vinculin, synaptopodin, and filamin C, compared to primary lesions. This study highlighted the proteomic differences between primary and metastatic liver CRC, providing insights into the clinical challenges of treating liver metastases.
Application of proteomic analysis in diagnosing CRLM
Some studies have identified CRLM predictive factors using proteomic analyses. For example, immunoproteomic profiling of tumor tissues revealed that OLFM4, CD11b, and ITGA2 are overexpressed in both primary colon tumors and liver metastases, suggesting that OLFM4, CD11b, and ITGA2 are potential CRC biomarkers with autoantigenic properties105. Additionally, proteomic analysis of urine samples identified the following two promising peptides derived from type I collagen: AGPP(-OH)GEAGKP(-OH)GEQGVP(-OH)GDLGAP(-OH)GP; and KGNSGEP(-OH)GAPGSKGDTGAKGEP(-OH)GPVG. This study provided a technical reference for the non-invasive detection of CRLM106. Similarly, other researchers have found that the combination of collagen-derived urinary AGP peptides and serum CEA levels may significantly improve the detection of patients at high risk of CRLM107. Furthermore, based on a study of collagen biomarkers in urine, the following two naturally occurring peptides (NOPs) were identified: GPP [collagen α-1(I)]; and GND [collagen α-1(III)]. The development of a molecular model combining these NOPs with CEA has enhanced the sensitivity and specificity of diagnosing high-risk patients for CRLM108.
Application of proteomic analysis in treating CRLM
Proteomic analysis may aid in identifying effective individualized cancer therapeutic targets. Receptor tyrosine kinases (RTKs), which are cell surface receptors involved in the regulation of key biological pathways, are critical targets for cancer treatment. One study quantified RTKs using proteomics to assess the protein abundance of 21 RTKs and demonstrated significant differences between metastatic and normal liver tissues. RET kinase was the most abundant in non-tumor tissues (approximately 35%), while PGFRB was the most abundant RTK in tumors (approximately 47%)109, highlighting the potential of RTK as a therapeutic target for tyrosine kinase inhibition. Turtoi et al.110 found that the proteomic characteristics of different regions of CRC liver metastases vary. Among the antigens that can uniformly target the entire liver lesion, LTBP2, and IGFBI are particularly suitable for liver metastasis-targeted therapy. Wertenbroek et al.111 conducted thermal ablation treatment (either cryoablation or radiofrequency ablation) in patients with CRC with liver metastases, followed by proteomic analysis. The results indicated that radiofrequency ablation is superior to cryoablation, with longer survival and a lower incidence of complications. Therefore, proteomic-based phenotypic analysis provides valuable support for the localized treatment of CRLM.
Application of proteomic analysis in the prognosis of CRLM
Proteomic analysis has been shown to predict the prognosis of patients with liver metastases. Upregulation of gankyrin (a small molecular protein) is positively correlated with disease progression and liver metastasis in patients with CRC112. Additionally, positive expression of the polymeric immunoglobulin receptor (pIgR) is significantly associated with poor prognosis in patients with CRLM113. In patients with CRC and liver metastasis, an elevated level of extracellular vesicle protein is significantly correlated with shortened OS both before and after surgery114. Zhu et al.115 developed a decision-tree model using eight proteomic features (m/z 3315, 6637, 1207, 1466, 4167, 4210, 2660, and 4186) and demonstrated the potential value in predicting liver metastasis in patients undergoing radical CRC surgery. Kirana et al.116 showed that the proteins, HLAB, ADAMTS 2, LTBP 3, JAG 2, and NME 2, are significantly associated with tumor progression, vascular invasion, and distant liver metastasis based on differential protein expression analysis. Spatial mapping of the tumor proteome in CRLM revealed that three specific peptides (histone H4, hemoglobin subunit alpha, and inosine-5′-monophosphate dehydrogenase 2) and two unknown m/z values (1305.840 and 1661.060) are significantly elevated in patients with shorter survival. These molecules may serve as potential prognostic biomarkers117.
Overall, proteomic analysis to study CRLM can assist in the diagnosis, treatment, and improvement of prognosis. Proteomics can be used to identify new drug targets, enabling the design of novel therapies targeting the different pathways involved in CRLM.
Applications of metabolomics in the progression, treatment, and prognosis of CRLM
Metabolomics involves qualitative and quantitative analyses of all metabolites within an organism to determine the relationships with physiologic or pathologic changes4. Metabolomics can identify potential biomarkers for CRLM, monitor metastatic progression, and assess treatment responses.
Serum metabolomic profiling can serve as a tool for distinguishing different stages of CRC. Notably, serum metabolomic features differ between patients with locally advanced CRC and patients with liver-only metastases. Furthermore, significant differences in serum metabolomic profiles were observed between patients with liver-only metastases and patients with extrahepatic metastases, indicating that serum metabolomics can be used for cancer staging and guiding treatment strategies118. The literature is replete with studies showing that metabolomics contributes to the non-invasive diagnosis and staging of CRC119,120. However, research on gut metabolomics with liver metastasis is limited and requires further investigation.
Numerous genes contribute to CRC progression through metabolite alterations and targeting these pathways can prolong survival in CRLM. Research has indicated that the EFNB2/EPHB4 axis increases cholesterol levels in CRLM with hypercholesterolemia linked to tumor growth. A high-cholesterol diet significantly shortens survival in mice, whereas blocking the EFNB2/EPHB4 axis extends the OS121,122. Metabolomic analysis suggests that chondroitin sulfate synthase 1 (CHSY1) gene activation promotes CD8+ T cell exhaustion via succinate metabolism, leading to liver metastasis. Artemisinin significantly inhibits CHSY1 and combined with anti-PD1 therapy, enhances treatment efficacy123. Aerobic glycolysis, which is a key metabolic pathway in tumor metabolism, alters pyruvate entry into mitochondria, which is targeted in cancer research124. Mitochondrial uncouplers, such as niclosamide ethanolamine (NEN) and oxyclozanide, show effective anticancer activity in mouse models of metastatic liver tumors124.
Metabolomics has also been explored in the treatment of liver metastasis in CRC. Metabolomic analysis of circulating tumor cells (CTCs) derived from CRLM revealed that the metabolite profiles and associated gene expression in CTC-derived metastatic tissue, particularly the specific upregulation of mRNAs involved in metabolic pathways, such as the folate one-carbon pool, folate biosynthesis, and histidine metabolism, differed significantly from those in parental cells. Additionally, tumors and metastatic sites derived from CTCs exhibit significantly higher invasiveness and chemoresistance than tumors and metastatic sites from cell lines. Therefore, identifying different metabolites could provide valuable insight into cancer progression and metastasis, thereby offering potential targets for treating metastatic CRC125. Mass spectrometry-based metabolomic analysis was used to characterize the metabolic profile of tumor-associated macrophages (TAMs), which showed that riboflavin is positively correlated with the morphology of larger TAMs (L-TAMs). Riboflavin regulates the expression and activity of lysine-specific demethylase 1A (LSD1) in L-TAMs and in vitro M2-polarized macrophages. The riboflavin-LSD1 axis influences TAM morphology in CRLM, making LSD1 a potential target for reprogramming TAM subtypes and offering new directions for anti-tumor therapy126.
Metabolomics is valuable for assessing the prognosis of liver metastases in CRC. Jonas et al.127 identified preoperative plasma metabolites, such as lysophosphatidylcholines (lysoPCs) and phosphatidylcholines (PCs), as potential predictors for postoperative disease recurrence in patients with CRLM. Costantini et al. 128 demonstrated that combined metabolic and lipidomic biomarkers help predict the prognosis of patients with metastatic CRC undergoing liver metastasis resection after induction therapy. Wu et al.129 categorized the metabolic phenotypes of CAFs within the TME and showed that scoring this classification could be used to evaluate chemotherapy sensitivity and survival prognosis in patients with liver metastases.
Applications of microbiomics in the diagnosis, treatment, progression, and prognosis of CRLM
Microbiomics focuses on the structure and function of microbial communities, interactions within these communities, and relationships between microbes and their environment or host. Microbiomics also explores regulation of microbial community growth, metabolism, and cooperative processes4. The intestinal microbiota significantly contributes to waste degradation, nutrient transport, and gut homeostasis maintenance. Dysbiosis and ecological imbalance in the colon can lead to physiologic disorders130, including the development of CRC131,132.
Microbiomics has been extensively explored in the context of pathogenesis, diagnosis, treatment, and prognosis. However, the role of CRLM is the subject of active research. Na et al. 133 reported that an increase in Proteus mirabilis and a decrease in Bacteroides vulgatus significantly influence CRLM, potentially due to a reduction in Kupffer cells in the liver. Studies have shown that cytolethal distending toxin (CDT) produced by Campylobacter jejuni promotes CRC metastasis to the liver and lungs via the JAK2-STAT3-MMP9 signaling pathway134. However, the specific mechanisms linking the gut microbiome to CRLM, as well as potential applications in diagnosing and predicting liver metastasis, require further investigation.
Studies have shown that microbiomics can be applied in treating liver metastasis. Antibiotic therapy is a potential strategy for CRC. Studies have suggested that changes in the oral microbiota may affect treatment outcomes and prognosis of primary CRC and CRLM radiotherapy. Fusobacterium nucleatum (Fn) in the oral cavity may migrate and colonize CRC lesions, contributing to radiotherapy resistance. However, combining antibiotics [e.g., metronidazole (MTZ)] can prevent this resistance135. Fn in CRC liver metastatic tissues correlates with a lower density of CD8(+) T cells and more myeloid-derived suppressor cells. These findings offer critical insights for the development of microbiota- and immune-targeted therapeutic strategies aimed at preventing and treating colorectal cancer metastasis to the liver.136. Wang et al.137 reported that MTZ inhibits tumor growth and reduces CRLM malignancy by modulating gut Fn and altering microbiota composition. Additionally, augmenting the microbiome has potential therapeutic value. Huang et al. 138 found fecal microbiota transplantation enhances the efficacy of PD-1 blockade in murine models of colorectal cancer. Xu et al.139 developed a nanodrug-engineered bacterium (LR-S-CD/CpG@LNP) in which oral administration of CD/CpG@LNPs induced the production of cytotoxic reactive oxygen species (ROS) at the CRC tumor site through photothermal and photodynamic effects, triggering immunogenic cell death (ICD) in CRC cells. This finding offered a new strategy for oral combination therapy in CRLM. Therefore, gut microbiota transplantation may serve as an adjunctive therapy for advanced CRC, particularly for patients with liver metastasis, by potentially enhancing treatment efficacy through modulation of the TME (Figure 3).
Microbiomics of colorectal cancer liver metastasis (CRLM). Microbiome-related techniques have identified the presence of Proteus mirabilis, Bacteroides vulgatus, and Fusobacterium nucleatum in the primary colorectal cancer (CRC) lesions. Fusobacterium nucleatum is also detected within CRLM sites, where its presence is associated with reduced infiltration of CD8+ T cells and increased accumulation of myeloid-derived suppressor cells (MDSCs). Therapeutic approaches, such as oral administration of metronidazole (MNZ) and a nanoengineered probiotic delivery system (LR-S-CD/CpG@LNP), have demonstrated efficacy in treating both primary CRC and CRLM. Fecal microbiota transplantation (FMT) has shown potential to synergize with PD-1 blockade therapy in CRLM.
Integrated application of multi-omics in the mechanisms, diagnosis, and treatment of CRLM
The integration of multi-omic approaches, including genomics, transcriptomics, epigenomics, microbiomics, metabolomics, and proteomics, offers promise for enhancing biomarker sensitivity and addressing population heterogeneity (Figure 4). Chen et al.140 used genomics, proteomics, and phosphoproteomics to investigate CRC cohorts and found that metastatic tissues are genetically like primary tumors, while significant differences emerged at the proteomic level. Proteomic and phosphoproteomic analyses of primary tumors effectively distinguished metastatic cases with kinase network analysis revealing heterogeneity between primary tumors and liver metastases. Researchers analyzing multi-omics data (genome, epigenome, and transcriptome) from cfDNA and cfRNA samples observed that multi-omics approaches outperformed single-omic data in identifying cancer-associated genes. Additionally, plasma cfRNA was strongly correlated with tumor RNA for specific cancer- and immune-related features, suggesting that cfRNAs is a non-invasive approach for the clinical monitoring cancer progression and patient status141. Therefore, multi-omics research is crucial for advancing the diagnosis and treatment of CRC.
Multi-omics analysis of colorectal cancer liver metastasis (CRLM). Integrated genomic, transcriptomic, and proteomic analyses have suggested a potential role of breast cancer gene 1 (BRCA1) in the early development of CRLM, indicating its promise as a tumor microenvironment (TME)-associated biomarker. Genomic and transcriptomic profiling has identified SMAD family member 4 (SMAD4) R361H/C mutations, which may confer resistance to bevacizumab and 5-fluorouracil (5-FU) via activation of the STAT3 signaling pathway. Four-jointed box 1 (FJX1) is considered to be a potential therapeutic target for CRLM through multi-omics analyses. Combined genomic and proteomic studies have also revealed that SMAD4, calponin-2, and glutathione peroxidase 3 may serve as predictive markers for early recurrence of CRLM.
Numerous studies have used multi-omics approaches to demonstrate the differences and mechanisms underlying primary and metastatic lesions in CRLM. Ham-Karim et al.101 conducted an integrated multi-omics (genomics, transcriptomics, and proteomics) analysis and found that the mutation status of 54 of 60 loci was identical between primary and metastatic tumors, whereas 40 of 58 gene loci exhibited heterogeneity. Furthermore, miRNA expression profiles differed between primary and metastatic tumors. Integrated genomics, transcriptomics, and proteomics analysis identified BRCA 1 as 1 of 20 transcripts simultaneously upregulated in all 3 types of TME liver cells during CRLM. The most likely sequence of cellular activation during metastasis involves endothelial (Ito) followed by Kupffer cells. Immunohistochemical analysis of human liver metastases revealed BRCA1 protein co-localization with Ito, Kupffer, and endothelial cells in 81.8% of early or synchronous metastases. This finding suggested that BRCA1 co-expression in these cells may contribute to the early stages of CRLM, thus highlighting the potential of BRCA1 as a TME biomarker142.
Multi-omics analysis holds promise for elucidating the mechanisms underlying treatment resistance in CRC. Shi et al.143 revealed changes in genes and the microenvironment during chemotherapy (bevacizumab combined with chemotherapy) for CRLM based on multi-omics analysis. Shi et al.143 found that 92% of cases exhibited minimal change at the genomic level but significant differences at the transcriptomic level. Specifically, the SMAD4 R361H/C mutation mediates resistance to bevacizumab and 5-FU via the STAT3 signaling pathway. Combining the STAT3 inhibitor, GB201, with the initial treatment restored the therapeutic sensitivity of SMAD4 R361H/C mutated cancer cells. Thus, multi-omics analysis reveals resistance mechanisms during treatment and enables effective intervention.
Studies on genomics and transcriptomics as prognostic biomarkers for CRLM have shown that CTAG1A, CSTL1, FJX1, IER5L, and KLHL35 expression is higher in metastatic lesions than primary tumors. These genes may contribute to cancer cell metastasis. Specifically, FJX1 mRNA expression is elevated in CRC tissues at various T stages and in metastatic lymph nodes than in normal mucosa or lymphatic tissues. Gene expression analysis suggested that FJX1 regulates chromatin-modifying enzymes, the Notch signaling pathway, cellular senescence, and other signaling pathways. Therefore, FJX1 may be a crucial target for CRLM and a potential novel biomarker for its prediction and treatment144.
Multi-omics analysis aids in prognosis monitoring for patients with CRLM. Wong et al.145 used genomics and proteomics to accurately identify patients at risk of early relapse after curative treatment. Wong et al.145 found the somatic mutations most associated with early relapse were TP53 (88%), APC (71%), KRAS (38%), and SMAD4 (21%). Additionally, proteins associated with early relapse included calponin-2, versican core protein, glutathione peroxidase 3, fibulin-5, and amyloid-β precursor protein. Further analysis indicated that SMAD4, calponin-2, and glutathione peroxidase 3 may have predictive potential for early relapse, which aids in optimizing prognosis evaluation for CRLM.
Overall, multi-omics data related to tumor cells and tissues provide a more comprehensive understanding of the pathogenesis, diagnostic and therapeutic strategies, drug resistance mechanisms, and prognostic evaluation of the corresponding diseases compared to single-omics data. However, research applying multi-omicss approaches to CRLM is limited.
Summary and future directions
Advances in genomics, transcriptomics, epigenomics, proteomics, metabolomics, and microbiomics have paved a solid foundation for precise oncologic research. CRLM is a key factor for poor prognosis in CRC. Applying omics approaches to CRLM can guide the diagnosis, treatment, and prognosis of advanced CRC. This strategy offers more effective management and improved survival outcomes for patients with liver metastasis (Table 1).
Multi-omics summary of colorectal cancer liver metastases
Proteomic-related research on CRLM has focused extensively on pathogenesis, diagnosis, metastatic vs. primary tumors, treatment, and prognosis. Future studies should explore clinical translation to provide more therapeutic options for patients with liver metastasis.
Current omics research on CRLM requires further basic studies to explore underlying mechanisms. Genomic studies mostly focus on differences between primary tumors and metastatic sites but miss mechanisms directly linked to liver metastasis onset and therapeutic strategies. Epigenomic research, especially on DNA methylation, lacks studies on treatments for CRLM. Research on histone methylation and demethylation is limited with few studies investigating specific mechanisms. Although transcriptomic studies on CRLM are abundant, progress in improving the diagnosis of liver metastasis is limited. Metabolomics research in CRC mainly targets lipid and gut microbiota metabolites, with few studies on liver metastasis. Microbiological research has focused on CRC initiation and progression, treatment, and prognosis. However, research into the role of the microbiota in CRLM has been minimal. These gaps in knowledge present significant opportunities for further exploration.
Conflict of interest statement
No potential conflicts of interest are disclosed.
Author contributions
Conceived and designed the analysis: Kexue Zhou, Yan Li.
Collected the data: Kexue Zhou.
Contributed data or analysis tools: Kexue Zhou, Chengxiang Yang.
Performed the analysis: Kexue Zhou, Yan Li.
Wrote the paper: Kexue Zhou.
- Received February 11, 2025.
- Accepted May 14, 2025.
- Copyright: © 2025, The Authors
This work is licensed under the Creative Commons Attribution-NonCommercial 4.0 International License.
References
- 1.↵
- 2.↵
- 3.↵
- 4.↵
- 5.↵
- 6.↵
- 7.↵
- 8.↵
- 9.↵
- 10.↵
- 11.↵
- 12.↵
- 13.↵
- 14.↵
- 15.↵
- 16.↵
- 17.↵
- 18.↵
- 19.↵
- 20.↵
- 21.↵
- 22.↵
- 23.↵
- 24.↵
- 25.↵
- 26.↵
- 27.↵
- 28.↵
- 29.↵
- 30.↵
- 31.
- 32.↵
- 33.↵
- 34.↵
- 35.↵
- 36.↵
- 37.↵
- 38.↵
- 39.↵
- 40.↵
- 41.↵
- 42.↵
- 43.↵
- 44.↵
- 45.↵
- 46.↵
- 47.↵
- 48.↵
- 49.↵
- 50.↵
- 51.↵
- 52.↵
- 53.↵
- 54.↵
- 55.↵
- 56.
- 57.↵
- 58.↵
- 59.↵
- 60.↵
- 61.↵
- 62.↵
- 63.↵
- 64.↵
- 65.↵
- 66.↵
- 67.↵
- 68.↵
- 69.↵
- 70.↵
- 71.↵
- 72.↵
- 73.↵
- 74.↵
- 75.↵
- 76.↵
- 77.↵
- 78.↵
- 79.↵
- 80.↵
- 81.↵
- 82.↵
- 83.↵
- 84.↵
- 85.↵
- 86.↵
- 87.↵
- 88.↵
- 89.↵
- 90.↵
- 91.↵
- 92.↵
- 93.↵
- 94.↵
- 95.↵
- 96.↵
- 97.↵
- 98.↵
- 99.↵
- 100.↵
- 101.↵
- 102.↵
- 103.↵
- 104.↵
- 105.↵
- 106.↵
- 107.↵
- 108.↵
- 109.↵
- 110.↵
- 111.↵
- 112.↵
- 113.↵
- 114.↵
- 115.↵
- 116.↵
- 117.↵
- 118.↵
- 119.↵
- 120.↵
- 121.↵
- 122.↵
- 123.↵
- 124.↵
- 125.↵
- 126.↵
- 127.↵
- 128.↵
- 129.↵
- 130.↵
- 131.↵
- 132.↵
- 133.↵
- 134.↵
- 135.↵
- 136.↵
- 137.↵
- 138.↵
- 139.↵
- 140.↵
- 141.↵
- 142.↵
- 143.↵
- 144.↵
- 145.↵











