# Integrative spatial and single-cell transcriptomics elucidate programmed cell death-driven tumor microenvironment dynamics in hepatocellular carcinoma

**Authors:** Kai Lei, Yutong Zhao, Shumin Li, Jiawei Liu, Wenhao Chen, Caihong Zhou, Yi Zhang, Jinmei Tan, Jian Wu, Qi Zhou, Jiehui Tan

PMC · DOI: 10.3389/fimmu.2025.1589563 · Frontiers in Immunology · 2025-07-16

## TL;DR

This study uses gene and cell data to predict liver cancer outcomes and understand tumor environments linked to cell death.

## Contribution

A novel PCD score model for HCC prognosis and tumor microenvironment insights using integrative transcriptomic approaches.

## Key findings

- Seventeen PCD-related genes identified as significant prognostic indicators for HCC.
- High PCD scores correlate with poor survival and distinct tumor cell and immune microenvironment features.
- UBE2E1 is identified as a key oncogenic gene in HCC linked to the PCD model.

## Abstract

Programmed cell death (PCD) mechanisms play crucial roles in cancer progression and treatment response. This study aims to develop a PCD scores prediction model to evaluate the prognosis of hepatocellular carcinoma (HCC) and elucidate the tumor microenvironment differences.

We analyzed transcriptomic data from 363 HCC patients in the TCGA database and 221 patients in the GEO database to develop a PCD prediction model. Single-cell RNA sequencing (scRNA-seq) and spatial transcriptomics sequencing (ST-seq) data from HCC patients were analyzed to investigate the tumor microenvironment and functional disparities. The oncogenic role of the key gene UBE2E1 in the model was explored in HCC through various in vitro experiments.

Seventeen PCD-related genes were identified as significant prognostic indicators, forming the basis of our PCD prediction model. High-PCD scores correlated with poorer overall survival (OS) and exhibited significant predictive capabilities. scRNA-seq analysis revealed distinct tumor cell characteristics and immune microenvironment differences between high- and low-PCD groups. High-PCD tumors showed increased cell proliferation and malignancy-associated gene expression. T cells in high-PCD patients were more likely to be exhausted, with elevated expression of exhaustion markers. ST-seq data also confirmed these results. Among the genes associated with the PCD prognostic model, UBE2E1 was identified as a key oncogenic marker in HCC.

The PCD prediction model effectively predicts prognosis in HCC patients and reveals critical insights into the tumor microenvironment and immune cell exhaustion. This study underscores the potential of PCD-related biomarkers in guiding personalized treatment strategies for HCC.

## Linked entities

- **Genes:** UBE2E1 (ubiquitin conjugating enzyme E2 E1) [NCBI Gene 7324]
- **Diseases:** hepatocellular carcinoma (MONDO:0007256), HCC (MONDO:0007256)

## Full-text entities

- **Genes:** ERBB2 (erb-b2 receptor tyrosine kinase 2) [NCBI Gene 2064] {aka CD340, HER-2, HER-2/neu, HER2, MLN 19, MLN-19}, CTLA4 (cytotoxic T-lymphocyte associated protein 4) [NCBI Gene 1493] {aka ALPS5, CD, CD152, CELIAC3, CTLA-4, GRD4}, IL7R (interleukin 7 receptor) [NCBI Gene 3575] {aka CD127, CDW127, IL-7R-alpha, IL-7Ralpha, IL7RA, IL7Ralpha}, CD47 (CD47 molecule) [NCBI Gene 961] {aka IAP, MER6, OA3}, CD8A (CD8 subunit alpha) [NCBI Gene 925] {aka CD8, CD8alpha, IMD116, Leu2, p32}, APOE (apolipoprotein E) [NCBI Gene 348] {aka AD2, APO-E, ApoE4, LDLCQ5, LPG}, HSPA1A (heat shock protein family A (Hsp70) member 1A) [NCBI Gene 3303] {aka HEL-S-103, HSP70, HSP70-1, HSP70-1A, HSP70-2, HSP70.1}, LAG3 (lymphocyte activating 3) [NCBI Gene 3902] {aka CD223}, EPCAM (epithelial cell adhesion molecule) [NCBI Gene 4072] {aka Ber-Ep4, BerEp4, DIAR5, EGP-2, EGP314, EGP40}, YWHAQ (tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein theta) [NCBI Gene 10971] {aka 14-3-3, 1C5, HS1}, GZMA (granzyme A) [NCBI Gene 3001] {aka CTLA3, HFSP}, CD3D (CD3 delta subunit of T-cell receptor complex) [NCBI Gene 915] {aka CD3-DELTA, CD3DELTA, IMD19, T3D}, IGFBP3 (insulin like growth factor binding protein 3) [NCBI Gene 3486] {aka BP-53, IBP-3, IBP3, IGFBP-3}, ANXA5 (annexin A5) [NCBI Gene 308] {aka ANX5, CPB-I, ENX2, HEL-S-7, PP4, RPRGL3}, TIGIT (T cell immunoreceptor with Ig and ITIM domains) [NCBI Gene 201633] {aka VSIG9, VSTM3, WUCAM}, CD4 (CD4 molecule) [NCBI Gene 920] {aka CD4mut, IMD79, Leu-3, OKT4D, T4}, KLHDC10 (kelch domain containing 10) [NCBI Gene 23008] {aka PNAS-138, slim}, CACYBP (calcyclin binding protein) [NCBI Gene 27101] {aka GIG5, PNAS-107, S100A6BP, SIP}, SLC4A10 (solute carrier family 4 member 10) [NCBI Gene 57282] {aka NBCn2, NCBE}, CXCL13 (C-X-C motif chemokine ligand 13) [NCBI Gene 10563] {aka ANGIE, ANGIE2, BCA-1, BCA1, BLC, BLR1L}, UBE2E1 (ubiquitin conjugating enzyme E2 E1) [NCBI Gene 7324] {aka UBCH6}, UBB (ubiquitin B) [NCBI Gene 7314] {aka HEL-S-50}, PCK2 (phosphoenolpyruvate carboxykinase 2, mitochondrial) [NCBI Gene 5106] {aka PEPCK, PEPCK-M, PEPCK2, mtPCK2}, STAT6 (signal transducer and activator of transcription 6) [NCBI Gene 6778] {aka D12S1644, HIES6, IL-4-STAT, STAT6B, STAT6C}, ENO1 (enolase 1) [NCBI Gene 2023] {aka ENO1-IT1, ENO1L1, HEL-S-17, MPB1, NNE, PPH}, SPP1 (secreted phosphoprotein 1) [NCBI Gene 6696] {aka BNSP, BSPI, ETA-1, OPN}, GZMB (granzyme B) [NCBI Gene 3002] {aka C11, CCPI, CGL-1, CGL1, CSP-B, CSPB}, NCR [NCBI Gene 4827], NFKB1 (nuclear factor kappa B subunit 1) [NCBI Gene 4790] {aka CVID12, EBP-1, KBF1, NF-kB, NF-kB1, NF-kappa-B1}, HDAC1 (histone deacetylase 1) [NCBI Gene 3065] {aka GON-10, HD1, KDAC1, RPD3, RPD3L1}, NKG7 (natural killer cell granule protein 7) [NCBI Gene 4818] {aka GIG1, GMP-17, p15-TIA-1}, ALB (albumin) [NCBI Gene 213] {aka FDAHT, HSA, PRO0883, PRO0903, PRO1341}, CCR7 (C-C motif chemokine receptor 7) [NCBI Gene 1236] {aka BLR2, CC-CKR-7, CCR-7, CD197, CDw197, CMKBR7}, RPS17 (ribosomal protein S17) [NCBI Gene 6218] {aka DBA4, RPS17L, RPS17L1, RPS17L2, S17, eS17}, PDLIM1 (PDZ and LIM domain 1) [NCBI Gene 9124] {aka CLIM1, CLP-36, CLP36, HEL-S-112, hCLIM1}, S100A11 (S100 calcium binding protein A11) [NCBI Gene 6282] {aka HEL-S-43, MLN70, S100C}, AFP (alpha fetoprotein) [NCBI Gene 174] {aka AFPD, FETA, HPAFP}, CDK4 (cyclin dependent kinase 4) [NCBI Gene 1019] {aka CMM3, MCPH31, PSK-J3}, ANXA1 (annexin A1) [NCBI Gene 301] {aka ANX1, LPC1}, KRT18 (keratin 18) [NCBI Gene 3875] {aka CK-18, CYK18, K18}, FABP1 (fatty acid binding protein 1) [NCBI Gene 2168] {aka FABPL, L-FABP}, CCL5 (C-C motif chemokine ligand 5) [NCBI Gene 6352] {aka D17S136E, RANTES, SCYA5, SIS-delta, SISd, TCP228}, RUNX3 (RUNX family transcription factor 3) [NCBI Gene 864] {aka AML2, CBFA3, PEBP2aC}, USF2 (upstream transcription factor 2, c-fos interacting) [NCBI Gene 7392] {aka FIP, bHLHb12}, HAVCR2 (hepatitis A virus cellular receptor 2) [NCBI Gene 84868] {aka CD366, HAVcr-2, KIM-3, SPTCL, TIM3, TIMD-3}, CD3E (CD3 epsilon subunit of T-cell receptor complex) [NCBI Gene 916] {aka CD3epsilon, IMD18, T3E, TCRE}, ENPEP (glutamyl aminopeptidase) [NCBI Gene 2028] {aka APA, CD249, gp160}, LDHA (lactate dehydrogenase A) [NCBI Gene 3939] {aka GSD11, HEL-S-133P, LDHM, PIG19}, CD79A (CD79a molecule) [NCBI Gene 973] {aka IGA, IGAlpha, MB-1, MB1}, CBS (cystathionine beta-synthase) [NCBI Gene 875] {aka HIP4}
- **Diseases:** OS (MESH:D011475), tumorigenesis (MESH:D063646), cytotoxic (MESH:D064420), MIBC (MESH:D000093284), lung and colorectal cancer (MESH:D015179), LDC (MESH:D016464), HCC (MESH:D006528), ICD (MESH:D003643), necrosis (MESH:D009336), metastasis (MESH:D009362), TNM (MESH:D008207), inflammation (MESH:D007249), Tumor (MESH:D009369)
- **Chemicals:** PI (MESH:D010716), PVDF (MESH:C024865), hematoxylin (MESH:D006416), penicillin (MESH:D010406), puromycin (MESH:D011691), crystal violet (MESH:D005840), DMEM (-), polybrene (MESH:D006583), glycine (MESH:D005998), SDS (MESH:D012967), sorafenib (MESH:D000077157), H&amp;E (MESH:D006371), TBS (MESH:D013725), eosin (MESH:D004801), PBS (MESH:D007854), streptomycin (MESH:D013307), CO2 (MESH:D002245)
- **Species:** Homo sapiens (human, species) [taxon 9606], Mus musculus (house mouse, species) [taxon 10090]
- **Cell lines:** HEK293T — Homo sapiens (Human), Transformed cell line (CVCL_0063), pLKO.1 — Mus musculus (Mouse), Hybridoma (CVCL_C7RB), SNU-449 — Homo sapiens (Human), Adult hepatocellular carcinoma, Cancer cell line (CVCL_0454), Huh-7 — Homo sapiens (Human), Adult hepatocellular carcinoma, Cancer cell line (CVCL_0336)

## Full text

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

## Figures

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

## References

49 references — full list in the complete paper: https://tomesphere.com/paper/PMC12308848/full.md

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