Analysis and validation of abnormal signaling pathways and immune cell infiltration characteristics in digestive system cancers based on peroxisome-related genes.

Abstract

BACKGROUND: Although emerging evidence suggests a role for peroxisomes in tumorigenesis, their functions in digestive cancers remain unclear. This study aims to investigate the association between peroxisomes and digestive tract tumors. METHODS: To systematically investigate peroxisomal functions in digestive cancers, we first constructed and validated tumor-specific prognostic signatures based on peroxisome-related genes (PRGs) through univariate Cox, least absolute shrinkage and selection operator (LASSO), and multivariate Cox regression analyses. We then characterized the tumor immune microenvironment (TIME) with CIBERSORT, X-CELL, and EPIC algorithms, and identified tumor-specific and common signalings via Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and gene set enrichment analysis (GSEA). Focusing on hepatocellular carcinoma (HCC), we experimentally validated peroxisome-related therapeutic responses by profiling signature genes in radioresistant cells and an orthotopic transarterial chemoembolization (TACE) rat model. PEX13 knockdown further assessed peroxisomal role in radiosensitivity and targeted therapy response. Clinical relevance of PEX13 was evaluated in HCC cohort. Single-cell RNA sequencing dataset and lipidomics further revealed peroxisomal mechanisms in HCC progression. Finally, peroxisomal function in colorectal cancer (CRC) was validated in vitro. RESULTS: Novel peroxisome-related prognostic signatures demonstrated strong predictive power in HCC, colon adenocarcinoma, rectal adenocarcinoma, pancreatic adenocarcinoma, gastric adenocarcinoma, esophageal adenocarcinoma, esophageal squamous cell carcinoma, and cholangiocarcinoma. High-risk patients displayed an immunosuppressive microenvironment, characterized by increased infiltration of regulatory T cells, M2 macrophages, Th2 cells, or cancer-associated fibroblasts, or Th1 cells' reduction. Peroxisomes engaged in several distinct yet convergent pathways, most notably "positive regulation of response to stimuli". HCC prognostic genes were dynamically regulated in response to therapeutic stimuli, including radiotherapy, targeted therapy, and TACE. Clinically, the expression of PEX13 was markedly upregulated in tumor tissues from therapy-resistant HCC patients. Mechanistically, peroxisomal dysfunction induced by silencing PEX13 in HCC or UBE2D2 in CRC may overcome therapeutic resistance (radiotherapy/ lenvatinib resistance in HCC, radioresistance in CRC) through reprogramming lipid metabolism. CONCLUSIONS: Peroxisomes act as pivotal regulators of digestive cancer progression by modulating signaling pathways, the TIME, therapeutic resistance, and lipid metabolism. Targeting peroxisomal function, particularly in high-risk subgroups of HCC and CRC, warrants further exploration as a promising therapeutic strategy.