Four Metabolic Genes Predict Survival in Aggressive Liver Cancer

Intrahepatic cholangiocarcinoma, the second most common primary malignant tumor of the liver, is one of oncology’s most formidable adversaries. Its incidence has surged by 140 percent over the past decade, yet the five-year overall survival rate remains stubbornly below 10 percent. Surgery and liver transplantation are still the only definitive treatments, and for patients whose tumors cannot be removed, gemcitabine-based chemotherapy remains the mainstay despite limited efficacy. A major reason for this grim picture is the tumor’s remarkable heterogeneity, which drives rapid progression and resistance to therapy. Now, a team of researchers reporting in the open-access journal Heliyon has taken a step toward taming this disease with a computational and experimental strategy that reads the tumor’s metabolic fingerprint, distilling it into a four-gene signature that outperforms conventional clinical markers in forecasting patient survival.

The study, led by Feng Li and Wenfeng Zhang along with colleagues in China, rests on a deceptively simple biological premise: cancer is not just a disease of uncontrolled growth, but of rewired metabolism. Tumor cells undergo metabolic reprogramming, altering how they consume glucose, synthesize lipids, and process amino acids to fuel their expansion. The genes that govern these processes, collectively termed metabolism-related genes, sit upstream of many of the changes clinicians currently measure. The authors argue that shifts in these regulatory genes precede changes in established tumor markers such as carbohydrate antigen 19-9, carcinoembryonic antigen, and lactic acid, making them potentially more sensitive and reliable indicators of disease course.

To test this idea, the team mined the Kyoto Encyclopedia of Genes and Genomes database for a comprehensive catalog of 3,216 metabolism-related genes. They then compared the expression of these genes between normal liver tissue and intrahepatic cholangiocarcinoma specimens using the TCGA-CHOL collection, which provided mRNA profiles and clinical data for 32 tumor samples and 8 normal controls. Rigorous statistical filtering, employing the limma framework with a false discovery rate threshold below 0.01 and absolute log2 fold changes exceeding 2, revealed a striking landscape of metabolic disruption: 1,002 differentially expressed metabolic genes, of which 717 were elevated and 285 were suppressed in tumors. Hierarchical clustering showed that these expression patterns alone were sufficient to cleanly separate cancerous from healthy tissue.

Functional enrichment analysis added biological texture to the numbers. Gene Ontology and KEGG pathway analyses indicated that the dysregulated genes clustered around the metabolism and synthesis of energy-bearing substances, catalytic activity, the chemical carcinogenesis pathway, and signaling through peroxisome proliferator-activated receptors, or PPARs. That last finding is particularly intriguing given prior evidence that PPARγ can suppress cholangiocarcinoma growth through p53-dependent pathways and that a microRNA targeting PPARγ contributes to gemcitabine resistance. The enrichment results suggest that metabolic gene dysregulation in this cancer is not random noise but converges on pathways with known roles in tumor initiation and treatment failure.

The next challenge was to compress more than a thousand candidate genes into a clinically usable predictor. The researchers first applied univariate Cox regression with a stringent significance threshold of p less than 0.005, which flagged six genes associated with elevated mortality risk: SLC16A3, B4GALNT1, ACSL4, MTHFD1L, PYGB, and AGPAT4. They then turned to LASSO Cox regression, a machine-learning technique that shrinks the coefficients of less informative variables toward zero to guard against overfitting in small datasets. With three-fold cross-validation and an optimal penalty parameter of 0.1196, the model settled on four genes. The resulting risk score is a weighted sum of their expression levels, with MTHFD1L carrying the largest coefficient at 0.059857, followed by ACSL4, PYGB, and SLC16A3.

The performance of this compact signature was remarkable. When patients in the TCGA cohort were split into high-risk and low-risk groups at the median score, Kaplan-Meier analysis showed dramatically worse overall survival in the high-risk category. More striking still, receiver operating characteristic analysis yielded an area under the curve of 0.994 for the risk score, dwarfing the predictive power of age (0.819), sex (0.623), histological grade (0.609), pathologic stage (0.545), and nodal or metastatic staging, which fell below 0.5. In both univariate and multivariate Cox regression, the risk score was the only factor significantly associated with overall survival, marking it as an independent prognostic variable rather than a proxy for tumor burden.

Crucially, the model survived contact with external data. In an independent validation cohort of 30 patients from the GEO database (dataset GSE107943), high-risk patients again died at markedly higher rates, and the risk score achieved an area under the curve of 0.789, the best among all tested variables. In this cohort, only vascular invasion and the risk score emerged as significant predictors in multivariate analysis. Subgroup analyses reinforced the biological plausibility of the signature: patients with better-prognosis molecular subclass A tumors showed significantly lower expression of MTHFD1L, SLC16A3, and PYGB, and correspondingly lower risk scores, than those in the poorer-prognosis subclass B.

The team did not stop at computational modeling. Quantitative reverse transcription polymerase chain reaction confirmed that all four genes were expressed at significantly higher levels in two cholangiocarcinoma cell lines, CCLP1 and HuCCT-1, than in immortalized normal biliary epithelial cells. Immunohistochemical data from the Human Protein Atlas showed that the proteins SLC16A3, ACSL4, and PYGB were also elevated in tumor tissue, while MTHFD1L protein was paradoxically lower in tumors despite elevated mRNA, a discrepancy the authors acknowledge and that may reflect post-transcriptional regulation. Each gene carries mechanistic intrigue: SLC16A3, a lactate transporter, is co-activated by mutant KRAS and CK2 to drive tumor progression; ACSL4 participates in ferroptosis, an iron-dependent form of cell death; MTHFD1L is a folate-cycle enzyme implicated in hepatocellular carcinoma; and PYGB, brain-type glycogen phosphorylase, fuels aerobic glycolysis. Notably, PYGB expression showed the strongest and most consistent association with survival across both cohorts, and prior work has shown that inhibiting PYGB slows cholangiocarcinoma progression by reprogramming glycolysis, making it a compelling therapeutic target.

The authors are candid about limitations. The training and validation cohorts were small, retrospective, and drawn from public repositories, and bulk transcriptomic data are vulnerable to technical and biological confounding. Prospective validation in larger, randomized populations or patient-derived xenograft models will be essential before the signature can inform clinical decisions. Even so, the study delivers what the field has lacked: a parsimonious, independently validated metabolic gene model that surpasses traditional staging in predictive accuracy. If larger trials confirm its performance, the four-gene score could help clinicians identify high-risk patients earlier, stratify them for intensified therapy, and, through targets like PYGB and SLC16A3, open new pharmacological avenues in a cancer that urgently needs them.

Subject of Research: A metabolism-related gene signature for predicting prognosis in intrahepatic cholangiocarcinoma

Article Title: A new predictive model for intrahepatic cholangiocarcinoma based on metabolism-related genes

Article References: Li, F., Su, D., Deng, X., Wang, M., Tan, J., Chen, B., He, W., Miao, C., & Zhang, W. (2026). A new predictive model for intrahepatic cholangiocarcinoma based on metabolism-related genes. Heliyon, 12(15), Article e45494. https://doi.org/10.1016/j.heliyon.2026.e45494

Image Credits: AI Generated

DOI: 10.1016/j.heliyon.2026.e45494

Keywords: intrahepatic cholangiocarcinoma, metabolic reprogramming, prognostic model, LASSO Cox regression, PYGB, SLC16A3, ACSL4, MTHFD1L, biomarker, TCGA, gene expression, liver cancer