Abstract:Background: Metabolic reprogramming plays a crucial role in the development of colorectal cancer (CRC), influencing tumor heterogeneity, the tumor microenvironment, and metastasis. While the interaction between metabolism and CRC is critical for developing personalized treatments, gaps remain in understanding how tumor cell metabolism affects prognosis. Our study introduces novel insights by integrating single-cell and bulk transcriptome analyses to explore the metabolic landscape within CRC cells and its mechanisms influencing disease progression. This approach allows us to uncover metabolic heterogeneity and identify specific metabolic genes impacting metastasis, which have not been thoroughly examined in previous studies. Methods: We sourced microarray and single-cell RNA sequencing datasets from the Gene Expression Omnibus (GEO) and bulk sequencing data for CRC from The Cancer Genome Atlas (TCGA). We employed Gene Set Variation Analysis (GSVA) to assess metabolic pathway activity, consensus clustering to identify CRC-specific transcriptome subtypes in bulkseq, and rigorous quality controls, including the exclusion of cells with high mitochondrial gene expression in scRNA seq. Advanced analyses such as AUCcell, infercnvCNV, Non-negative Matrix Factorization (NMF), and CytoTRACE were utilized to dissect the cellular landscape and evaluate pathway activities and tumor cell stemness. The hdWGCNA algorithm helped identify prognosis-related hub genes, integrating these findings using a random forest machine learning model. Results: Kaplan-Meier survival curves identified 21 significant metabolic pathways linked to prognosis, with consensus clustering defining three CRC subtypes (C3, C2, C1) based on metabolic activity, which correlated with distinct clinical outcomes. The metabolic activity of the 13 cell subpopulations, particularly the epithelial cell subpopulation with active metabolic levels, was evaluated using AUCcell in scRNA seq. To further analyze tumor cells using infercnv, NMF disaggregated these cells into 10 cellular subpopulations. Among these, the C2 subpopulation exhibited higher stemness and tended to have a poorer prognosis compared to C6 and C0. Conversely, the C8, C3, and C1 subpopulations demonstrated a higher level of the five metabolic pathways, and the C3 and C8 subpopulations tended to have a more favorable prognosis. hdWGCNA identified 20 modules, from which we selected modules primarily expressed in high metabolic tumor subgroups and highly correlated with clinical information, including blue and cyan. By applying variable downscaling of RF to a total of 50 hub genes, seven gene signatures were obtained. Furthermore, molecules that were validated to be protective in GEO were screened alongside related molecules, resulting in the identification of prognostically relevant molecules such as UQCRFS1 and GRSF1. Additionally, the expression of GRSF1 was examined in colon cancer cell lines using qPCR and phenotypically verified by in vitro experiments. Conclusion: Our findings emphasize that high activity in specific metabolic pathways, including pyruvate metabolism and the tricarboxylic acid cycle, correlates with improved colon cancer outcomes, presenting new avenues for metabolic-based therapies. The identification of hub genes like GRSF1 and UQCRFS1 and their link to favorable metabolic profiles offers novel insights into tumor neovascularization and metastasis, with significant clinical implications for targeting metabolic pathways in CRC therapy.

Single-cell analyses reveal the metabolic heterogeneity and plasticity of the tumor microenvironment during head and neck squamous cell carcinoma progression

Integrative single-cell and bulk transcriptomes analyses reveals heterogeneity of serine-glycine-one-carbon metabolism with distinct prognoses and therapeutic vulnerabilities in HNSCC

Metabolic heterogeneity in clear cell renal cell carcinoma revealed by single-cell RNA sequencing and spatial transcriptomics

Single Cell Analysis Reveals Reciprocal Tumor-Macrophage Intercellular Communications Related with Metabolic Reprogramming in Stem-like Gastric Cancer

Cellular hierarchy framework based on single-cell/multi-patient sample sequencing reveals metabolic biomarker PYGL as a therapeutic target for HNSCC

Single cell atlas reveals multilayered metabolic heterogeneity across tumour types

Single-cell Sequencing of Head and Neck Carcinoma: Transcriptional Landscape and Prognostic Model Based on Malignant Epithelial Cell Features

Metabolic pathway-based subtypes associate glycan biosynthesis and treatment response in head and neck cancer

Single cell deciphering of progression trajectories of the tumor ecosystem in head and neck cancer

Single‐cell RNA sequencing reveals aberrant sphingolipid metabolism in non‐small cell lung cancer impacts tumor‐associated macrophages and stimulates angiogenesis via macrophage inhibitory factor signaling

In Silico analysis and linking of metabolism-related genes with the immune landscape in head and neck squamous cell carcinoma

Identify metabolism-related genes IDO1, ALDH2, NCOA2, SLC7A5, SLC3A2, LDHB, and HPRT1 as potential prognostic markers and correlate with immune infiltrates in head and neck squamous cell carcinoma

Spatial transcriptomics reveals that metabolic characteristics define the tumor immunosuppression microenvironment via iCAF transformation in oral squamous cell carcinoma

Metabolomics, Transcriptome and Single-Cell RNA Sequencing Analysis of the Metabolic Heterogeneity between Oral Cancer Stem Cells and Differentiated Cancer Cells

Unraveling metabolic characteristics and clinical implications in gastric cancer through single-cell resolution analysis

Metabolic Landscape of the Tumor Microenvironment at Single Cell Resolution

Metabolic Heterogeneity of Tumor Cells and its Impact on Colon Cancer Metastasis: Insights from Single-Cell and Bulk Transcriptome Analyses

Targeting Cellular Metabolism in Head and Neck Cancer Precision Medicine Era: A Promising Strategy to Overcome Therapy Resistance

MicroRNAs: Their Role in Metabolism, Tumor Microenvironment, and Therapeutic Implications in Head and Neck Squamous Cell Carcinoma

Single-Cell Transcriptomic Analysis of Primary and Metastatic Tumor Ecosystems in Head and Neck Cancer

Single-cell transcriptomic analyses of tumor microenvironment and molecular reprograming landscape of metastatic laryngeal squamous cell carcinoma