# Single-cell RNAseq reveals adverse metabolic transcriptional program in intrahepatic cholangiocarcinoma malignant cells

**Authors:** Christophe Desterke, Raquel Francés, Claudia Monge, Yuanji Fu, Agnès Marchio, Pascal Pineau, Jorge Mata-Garrido

PMC · DOI: 10.1016/j.bbrep.2025.101949 · Biochemistry and Biophysics Reports · 2025-02-15

## TL;DR

This study identifies a harmful metabolic gene pattern in liver cancer cells that predicts worse patient outcomes, independent of other factors.

## Contribution

The study reveals a novel metabolic transcriptional program in ICA malignant cells linked to poor prognosis.

## Key findings

- 16 metabolic genes are overexpressed in ICA patients with poor survival.
- A metabolic gene score predicts worse overall survival independent of fibrosis and tumor stage.
- Seven key metabolic genes are mainly expressed in malignant cells and mark the proliferative ICA subclass.

## Abstract

Intrahepatic cholangiocarcinoma (ICA) is a highly aggressive primary liver cancer, which originates from the epithelial cells of the bile ducts. The transcriptional profile of metabolic enzymes was investigated at both bulk and single-cell levels in tumor samples from distinct ICA cohorts. In a training cohort (TCGA consortium), 16 genes encoding for metabolic enzymes were found overexpressed in cases with poor survival. A computed metabolic gene expression score was significantly associated with worse ICA prognosis at the univariate level (overall survival [OS] log-rank p = 8.2e-4). After adjusting for Ishak fibrosis score and tumor staging, the metabolic expression remained an independent predictor of poor prognosis (multivariate OS log-rank p = 0.01). Seven genes encoding key enzymes (FH, MAT2B, PLOD2, PLOD1, PDE6D, ALDOC, and NT5DC3) were validated as markers of the proliferative subclass of ICA in the GSE32225 dataset, related to poor prognosis. The metabolic score was significantly different between the inflammatory and proliferative subclasses in the validation cohort (p < 2.2e-16). At the single-cell level, in the tumor microenvironment of 10 ICA patients, these seven enzymes were predominantly expressed by malignant cells. The single-cell metabolic score was thus higher in malignant cells. This study identifies a metabolic transcriptional program linked to poor prognosis in ICA, independent of fibrosis and tumor staging.

•Metabolic transcriptional program linked to poor prognosis identified in intrahepatic cholangiocarcinoma (ICA).•16 metabolic enzymes were overexpressed in ICA patients with poor survival in the TCGA cohort.•A metabolic gene score predicts worse OS, independent of fibrosis and tumor stage.•Seven metabolic genes mark the proliferative ICA subclass.•Single-cell analysis shows these enzymes are mainly in malignant tumor cells.•The metabolic gene expression score was significantly higher in proliferative ICA subclass, linking it to poor prognosis.

Metabolic transcriptional program linked to poor prognosis identified in intrahepatic cholangiocarcinoma (ICA).

16 metabolic enzymes were overexpressed in ICA patients with poor survival in the TCGA cohort.

A metabolic gene score predicts worse OS, independent of fibrosis and tumor stage.

Seven metabolic genes mark the proliferative ICA subclass.

Single-cell analysis shows these enzymes are mainly in malignant tumor cells.

The metabolic gene expression score was significantly higher in proliferative ICA subclass, linking it to poor prognosis.

## Linked entities

- **Genes:** FH (fumarate hydratase) [NCBI Gene 2271], MAT2B (methionine adenosyltransferase 2 non-catalytic beta subunit) [NCBI Gene 27430], PLOD2 (procollagen-lysine,2-oxoglutarate 5-dioxygenase 2) [NCBI Gene 5352], PLOD1 (procollagen-lysine,2-oxoglutarate 5-dioxygenase 1) [NCBI Gene 5351], PDE6D (phosphodiesterase 6D) [NCBI Gene 5147], ALDOC (aldolase, fructose-bisphosphate C) [NCBI Gene 230], NT5DC3 (5'-nucleotidase domain containing 3) [NCBI Gene 51559]
- **Diseases:** intrahepatic cholangiocarcinoma (MONDO:0003210), liver cancer (MONDO:0002691)

## Full-text entities

- **Genes:** PLOD1 (procollagen-lysine,2-oxoglutarate 5-dioxygenase 1) [NCBI Gene 5351] {aka EDS6, EDSKCL1, LH, LH1, LLH, PLOD}, PDE6D (phosphodiesterase 6D) [NCBI Gene 5147] {aka JBTS22, PDED}, NT5DC3 (5'-nucleotidase domain containing 3) [NCBI Gene 51559] {aka TU12B1-TY, TU12B1TY}, PLOD2 (procollagen-lysine,2-oxoglutarate 5-dioxygenase 2) [NCBI Gene 5352] {aka BRKS2, LH2, TLH}, MAT2B (methionine adenosyltransferase 2 non-catalytic beta subunit) [NCBI Gene 27430] {aka MAT-II, MATIIbeta, Nbla02999, SDR23E1, TGR}, ALDOC (aldolase, fructose-bisphosphate C) [NCBI Gene 230] {aka ALDC}
- **Diseases:** fibrosis (MESH:D005355), tumor (MESH:D009369), primary liver cancer (MESH:D006528), ICA (MESH:D018281), inflammatory (MESH:D007249)
- **Species:** Homo sapiens (human, species) [taxon 9606]

## Full text

_Full body text omitted from this summary view._ Fetch the complete paper as Markdown: https://tomesphere.com/paper/PMC11872667/full.md

## Figures

7 figures with captions in the complete paper: https://tomesphere.com/paper/PMC11872667/full.md

## References

67 references — full list in the complete paper: https://tomesphere.com/paper/PMC11872667/full.md

---
Source: https://tomesphere.com/paper/PMC11872667