Front Pharmacol. 2026 ;17
1732137
Background: Colorectal cancer (CRC) remains a major cause of cancer-related morbidity and mortality, with high recurrence rates and limited treatment options for metastatic disease. The tumor microenvironment (TME) and metabolic reprogramming are critical drivers of CRC progression, influencing immune responses, therapeutic resistance, and patient outcomes.
Objective: This study explores the interplay between metabolic reprogramming and the TME in CRC using transcriptomic data and bioinformatics approaches to identify metabolically and microenvironmentally defined CRC subtypes and candidate biomarkers.
Methods: Gene expression and clinical data were obtained from TCGA colorectal adenocarcinoma (COAD), rectal adenocarcinoma (READ), and six GEO CRC datasets. Immunohistochemistry (IHC) was performed to validate PDE2A and CKMT2 expression in CRC tissues. Bioinformatic analyses were conducted using R software v4.0.3.
Results: We identified 220 TME- and 40 metabolism-related differentially expressed genes (DEGs) in CRC. Consensus clustering of these TMET genes revealed two distinct subtypes: Cluster 1 (C1), associated with poorer survival, an immune-mesenchymal phenotype, and frequent mutations in TTN and BRAF, and Cluster 2 (C2), characterized by enriched TP53 and APC mutations, classic tumor suppressor pathway activation, and higher genomic instability. Metabolically, C1 was characterized by lipid metabolism and extracellular matrix remodeling, whereas C2 showed enrichment of nucleotide and amino acid metabolism linked to cell cycle progression and DNA repair. Single-cell RNA sequencing confirmed these distinctions, revealing that C1-upregulated genes were predominantly expressed in immune and stromal compartments, whereas C2-upregulated genes were enriched in epithelial and malignant cells. PDE2A, primarily expressed by endothelial cells, was identified as a metabolic biomarker of C1, while CKMT2, expressed in malignant cells, defined C2. These genes serve as key metabolic markers distinguishing CRC subtypes based on molecular heterogeneity and prognosis.
Conclusion: PDE2A and CKMT2 were identified as critical metabolic biomarkers associated with distinct CRC subtypes and TME compositions. These findings highlight the intricate relationship between metabolic reprogramming, the tumor microenvironment, and tumor heterogeneity, providing insights into CRC molecular subtypes and their prognostic significance.
Keywords: RNA sequencing; cancer hallmarks; colorectal cancer; metabolic reprogramming; tumor microenvironment