# Transcriptomic Analysis Reveals an NRF2-Mediated Redox and Metabolic Reprogramming in Sorafenib-Resistant Hepatocellular Carcinoma Cells

**Authors:** Angelo Michilli, Cristian Bassi, Farzaneh Moshiri, Bruno De Siena, Rosaria Marinaro, Elisa Callegari, Massimo Negrini, Silvia Sabbioni

PMC · DOI: 10.3390/biotech15010018 · BioTech · 2026-02-11

## TL;DR

This study finds that resistance to sorafenib in liver cancer is linked to a change in cell metabolism and redox balance, driven by a protein called NRF2.

## Contribution

The study identifies a novel NRF2-driven metabolic and redox reprogramming mechanism in sorafenib-resistant hepatocellular carcinoma.

## Key findings

- NRF2-regulated genes like Gpx4, Txn1, and Hmox1 are upregulated in resistant cells, enhancing antioxidant defense.
- Pharmacological inhibition of NRF2 restores sorafenib sensitivity in resistant cells.
- Resistant cells show metabolic and redox reprogramming that suppresses ferroptosis and promotes survival.

## Abstract

Despite the advent of immune checkpoint inhibitor-based regimens, sorafenib remains an important therapeutic option for patients with advanced hepatocellular carcinoma (HCC) who are ineligible for immunotherapy. However, its clinical efficacy is limited by the emergence of drug resistance, whose underlying molecular mechanisms remain incompletely understood. To investigate these mechanisms, we established a murine model of acquired sorafenib resistance and performed comparative RNA sequencing of sorafenib-sensitive versus -resistant Hep55.1C hepatoma cells. Transcriptomic profiling revealed a distinct resistance-associated signature comprising 1264 significantly deregulated genes (adjusted p < 0.03, fold change > 3.0). Pathway analysis and Gene Set Enrichment Analyses (GSEA) indicated a coordinated downregulation of metabolic and intercellular signaling pathways, accompanied by marked upregulation of redox-regulatory, mitochondrial and cellular stress-response programs. Genes transcriptionally regulated by nuclear factor erythroid 2-related factor 2 (NRF2) including Gpx4, Txn1, Txnrd1, Hmox1, Fth1, Taldo1, Phgdh, and MafG, involved in antioxidant defense, ferroptosis suppression and metabolic rewiring, were all upregulated in resistant cells. Pharmacological inhibition of NRF2 activity using brusatol restored sensitivity to sorafenib, functionally implicating NRF2-dependent pathways in the maintenance of the resistant phenotype. These findings demonstrate that acquired sorafenib resistance in HCC is associated with a stable NRF2-driven transcriptional and metabolic reprogramming that enhances antioxidant capacity, suppresses ferroptosis and promotes tumor cell survival. Targeting NRF2-regulated redox metabolism may therefore represent a promising strategy to overcome therapeutic resistance in HCC.

## Linked entities

- **Genes:** GABPA (GA binding protein transcription factor subunit alpha) [NCBI Gene 2551], GPX4 (glutathione peroxidase 4) [NCBI Gene 2879], TXN (thioredoxin) [NCBI Gene 7295], TXNRD1 (thioredoxin reductase 1) [NCBI Gene 7296], HMOX1 (heme oxygenase 1) [NCBI Gene 3162], FTH1 (ferritin heavy chain 1) [NCBI Gene 2495], TALDO1 (transaldolase 1) [NCBI Gene 6888], PHGDH (phosphoglycerate dehydrogenase) [NCBI Gene 26227], MAFG (MAF bZIP transcription factor G) [NCBI Gene 4097]
- **Chemicals:** sorafenib (PubChem CID 216239)
- **Diseases:** hepatocellular carcinoma (MONDO:0007256), HCC (MONDO:0007256)
- **Species:** Mus musculus (taxon 10090)

## Full-text entities

- **Genes:** Gsr (glutathione reductase) [NCBI Gene 14782] {aka D8Ertd238e, Gr-1, Gr1}, Keap1 (kelch-like ECH-associated protein 1) [NCBI Gene 50868] {aka INRF2, mKIAA0132}, Mafg (Maf bZIP transcription factor G) [NCBI Gene 17134] {aka C630022N07Rik}, Hmox1 (heme oxygenase 1) [NCBI Gene 15368] {aka D8Wsu38e, HO-1, HO1, Hemox, Hmox, Hsp32}, Nr3c1 (nuclear receptor subfamily 3, group C, member 1) [NCBI Gene 14815] {aka GR, Grl-1, Grl1}, Txn1 (thioredoxin 1) [NCBI Gene 22166] {aka ADF, Trx1, Txn}, Nfe2l2 (nuclear factor, erythroid derived 2, like 2) [NCBI Gene 18024] {aka Nrf2}, Actb (actin, beta) [NCBI Gene 11461] {aka Actx, E430023M04Rik, beta-actin}, Slc7a11 (solute carrier family 7 (cationic amino acid transporter, y+ system), member 11) [NCBI Gene 26570] {aka 9930009M05Rik, sut, xCT}, PDGFRB (platelet derived growth factor receptor beta) [NCBI Gene 5159] {aka CD140B, IBGC4, IMF1, JTK12, KOGS, OPDKD}, KIT (KIT proto-oncogene, receptor tyrosine kinase) [NCBI Gene 3815] {aka C-Kit, CD117, MASTC, PBT, SCFR}, Psat1 (phosphoserine aminotransferase 1) [NCBI Gene 107272] {aka D8Ertd814e, EPIP, PSA, Psat}, NFE2L2 (NFE2 like bZIP transcription factor 2) [NCBI Gene 4780] {aka IMDDHH, NRF2, Nrf-2}, Gpx4 (glutathione peroxidase 4) [NCBI Gene 625249] {aka GPx-4, GSHPx-4, PHGPx, mtPHGPx, snGPx}, Fth1 (ferritin heavy polypeptide 1) [NCBI Gene 14319] {aka FHC, Fth, HFt, MFH}, Phgdh (3-phosphoglycerate dehydrogenase) [NCBI Gene 236539] {aka 3-PGDH, 3PGDH, 4930479N23, A10, PGAD, PGD}, Txnrd1 (thioredoxin reductase 1) [NCBI Gene 50493] {aka TR, TR1, TrxR1}, Shmt2 (serine hydroxymethyltransferase 2 (mitochondrial)) [NCBI Gene 108037] {aka 2700043D08Rik, SHMT}, Slc3a2 (solute carrier family 3 (activators of dibasic and neutral amino acid transport), member 2) [NCBI Gene 17254] {aka 4F2, 4F2HC, Cd98, Ly-10, Ly-m10, Ly10}, Taldo1 (transaldolase 1) [NCBI Gene 21351], Gapdh (glyceraldehyde-3-phosphate dehydrogenase) [NCBI Gene 14433] {aka Gapd}, Psph (phosphoserine phosphatase) [NCBI Gene 100678] {aka PSP, PSPase}
- **Diseases:** Tumors (MESH:D009369), hyperprogressive disease (MESH:D004194), injury to (MESH:D014947), Mitochondrial dysfunction (MESH:D028361), infection (MESH:D007239), cytotoxic (MESH:D064420), viral infection (MESH:D014777), HCC (MESH:D006528)
- **Chemicals:** CO (MESH:D002248), paclitaxel (MESH:D017239), Glycosphingolipid (MESH:D006028), pentose phosphate (MESH:D010428), streptomycin (MESH:D013307), carbon (MESH:D002244), lipid peroxides (MESH:D008054), methionine (MESH:D008715), SDS (MESH:D012967), Bay 43-9006 (MESH:D000077157), bevacizumab (MESH:D000068258), glycine (MESH:D005998), 3-phosphohydroxypyruvate (MESH:C012488), iron (MESH:D007501), 3-phosphoglycerate (MESH:C005156), atezolizumab (MESH:C000594389), thiols (MESH:D013438), heme (MESH:D006418), Amino acid (MESH:D000596), NADP+ (MESH:D009249), serine (MESH:D012694), GSSG (MESH:D019803), ribose-5-phosphate (MESH:C031626), cisplatin (MESH:D002945), 3PG (-), alpha-ketoglutarate (MESH:D007656), fructose-6-phosphate (MESH:C027618), PUFA (MESH:D005231), doxorubicin (MESH:D004317), penicillin (MESH:D010406), cystine (MESH:D003553), Sphingolipid (MESH:D013107), Brusatol (MESH:C020237), PVDF (MESH:C024865), DMSO (MESH:D004121), ROS (MESH:D017382), 3-phosphoserine (MESH:D010768), L-glutamine (MESH:D005973), CO2 (MESH:D002245), GSH (MESH:D005978), cysteine (MESH:D003545), Lipid (MESH:D008055)
- **Species:** Mus musculus (house mouse, species) [taxon 10090], Homo sapiens (human, species) [taxon 9606]
- **Cell lines:** Hep53.4 — Mus musculus (Mouse), Hepatocellular carcinoma of the mouse, Cancer cell line (CVCL_5765), C57BL/6 — Mus musculus (Mouse), Transformed cell line (CVCL_C0MU), Hep55-RES — Mus musculus (Mouse), Hepatocellular carcinoma of the mouse, Cancer cell line (CVCL_5766), C57BL/6J — Mus musculus (Mouse), Transformed cell line (CVCL_C0MW), H55-RES — Homo sapiens (Human), Medulloblastoma, Cancer cell line (CVCL_VU84), H55 — Homo sapiens (Human), Diabetic retinopathy, Induced pluripotent stem cell (CVCL_B3NU)

## Full text

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

## Figures

9 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12921911/full.md

## References

58 references — full list in the complete paper: https://tomesphere.com/paper/PMC12921911/full.md

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