Single-nucleus RNA sequencing reveals the pathological process and key transcription factor of sorafenib-induced cardiotoxicity.
6 affiliations
- 1Department of Cardiology, The Second Affiliated Hospital of Harbin Medical University, Harbin, 150086, China.
- 2State Key Laboratory of Frigid Zone Cardiovascular Diseases, Department of Cardiology and Cardiovascular Research Institute, General Hospital of Northern Theater Command, Shenyang, 110016, Liaoning Province, China.
- 3Department of Cardiology, The Second Affiliated Hospital of Harbin Medical University, Harbin, 150086, China. [email protected].
- 4State Key Laboratory of Frigid Zone Cardiovascular Diseases, Department of Cardiology and Cardiovascular Research Institute, General Hospital of Northern Theater Command, Shenyang, 110016, Liaoning Province, China. [email protected].
- 5Department of Cardiology, The Second Affiliated Hospital of Harbin Medical University, Harbin, 150086, China. [email protected].
- 6State Key Laboratory of Frigid Zone Cardiovascular Diseases, Department of Cardiology and Cardiovascular Research Institute, General Hospital of Northern Theater Command, Shenyang, 110016, Liaoning Province, China. [email protected].
Abstract
Background
Sorafenib is a first-line treatment for advanced solid tumors. However, its clinical application is limited by cardiotoxicity. The underlying mechanisms remain unclear. We aimed to identify key mechanisms of sorafenib-induced cardiotoxicity using single-nucleus RNA sequencing (snRNA-seq) and experimental verification.Methods
SnRNA-seq was performed on cardiac tissue from a mouse model of sorafenib-induced cardiotoxicity. Bioinformatics analyses included weighted gene co-expression network analysis (WGCNA), differentially expressed gene (DEG) analysis, intercellular communication inference, and transcriptional regulatory network analysis. We conducted functional verification both in vivo and in vitro.Results
SnRNA-seq demonstrated that cell clustering and component analysis confirmed the reduction in cardiomyocytes (CMs) and endothelial cells (ECs), along with an increase in immune cell populations. DEG enrichment analysis indicated that the upregulated DEGs in CMs, ECs, and fibroblasts (FBs) were consistently enriched in pathways associated with cardiomyopathy, whereas the downregulated DEGs were enriched in myocardial remodeling and repair processes. Gene set scoring revealed activation of fibrosis and inflammation programs in CMs and FBs. Cell communication analysis revealed a widespread disruption of cellular crosstalk in sorafenib-induced cardiotoxicity and identified key ligand-receptor pairs contributing to pro-inflammatory and pro-fibrotic phenotypes. Meanwhile we identified the common transcription factors and main effect genes in the three main effector cells during sorafenib-induced cardiotoxicity. Augur analysis and WGCNA consistently identified CMs as the most significantly affected cell type. Transcriptional regulatory network analysis highlighted melanocyte-inducing transcription factor (MITF) as a core transcription factor involved in sorafenib-induced cardiotoxicity. Functional validation showed that MITF expression was significantly increased in cardiac tissues and CMs after sorafenib treatment and that knockdown of MITF attenuated sorafenib-mediated CM hypertrophy in vitro.Conclusions
By constructing the first snRNA-seq data for sorafenib-induced cardiotoxicity, we identified MITF as a key transcription factor associated with sorafenib-induced cardiotoxicity. Our study provides potential therapeutic targets for mitigating sorafenib-induced cardiotoxicity.
Related
- Anthracycline Cardiotoxicity as Multicellular Network Failure: Cardiac Niche Reprogramming and Pathological Cell–Cell Communication
- Cost-effectiveness analysis of penpulimab plus anlotinib versus sorafenib as first-line therapy for advanced hepatocellular carcinoma
- Efficacy of Sorafenib Combined with Venetoclax and Azacitidine in Newly Diagnosed FLT3-ITD-Mutant Acute Myeloid Leukemia: Remission Efficacy and Prognostic Influence of Post-Remission Consolidation Therapy
- Sorafenib-Induced Hand-Foot Skin Reaction in a Patient With Acute Myeloid Leukemia Following Haploidentical Stem Cell Transplantation: A Case Report.
- COL11A1 promotes Sorafenib resistance in lung adenocarcinoma: from bioinformatic discovery to functional validation.
BibTeX
@article{wang2026singlenucleus,
title = {{Single-nucleus RNA sequencing reveals the pathological process and key transcription factor of sorafenib-induced cardiotoxicity.}},
author = {Yani Wang and Zheming Yang and Jing Wang and Ziqi Liu and Yan Zhang and Xiaolin Zhang and Haixu Song and Chenghui Yan and Dan Liu and Yaling Han},
journal = {BMC medicine},
year = {2026},
month = jul,
doi = {10.1186/s12916-026-05060-y},
url = {https://europepmc.org/article/MED/42897636},
note = {CC BY-NC-ND},
}
RIS
TY - JOUR
TI - Single-nucleus RNA sequencing reveals the pathological process and key transcription factor of sorafenib-induced cardiotoxicity.
AU - Yani Wang
AU - Zheming Yang
AU - Jing Wang
AU - Ziqi Liu
AU - Yan Zhang
AU - Xiaolin Zhang
AU - Haixu Song
AU - Chenghui Yan
AU - Dan Liu
AU - Yaling Han
PY - 2026
DA - 2026/07/01
JO - BMC medicine
DO - 10.1186/s12916-026-05060-y
UR - https://europepmc.org/article/MED/42897636
AB - Background Sorafenib is a first-line treatment for advanced solid tumors. However, its clinical application is limited by cardiotoxicity. The underlying mechanisms remain unclear. We aimed to identify key mechanisms of sorafenib-induced cardiotoxicity using single-nucleus RNA sequencing (snRNA-seq) and experimental verification.Methods SnRNA-seq was performed on cardiac tissue from a mouse model of sorafenib-induced cardiotoxicity. Bioinformatics analyses included weighted gene co-expression network analysis (WGCNA), differentially expressed gene (DEG) analysis, intercellular communication inference, and transcriptional regulatory network analysis. We conducted functional verification both in vivo and in vitro.Results SnRNA-seq demonstrated that cell clustering and component analysis confirmed the reduction in cardiomyocytes (CMs) and endothelial cells (ECs), along with an increase in immune cell populations. DEG enrichment analysis indicated that the upregulated DEGs in CMs, ECs, and fibroblasts (FBs) were consistently enriched in pathways associated with cardiomyopathy, whereas the downregulated DEGs were enriched in myocardial remodeling and repair processes. Gene set scoring revealed activation of fibrosis and inflammation programs in CMs and FBs. Cell communication analysis revealed a widespread disruption of cellular crosstalk in sorafenib-induced cardiotoxicity and identified key ligand-receptor pairs contributing to pro-inflammatory and pro-fibrotic phenotypes. Meanwhile we identified the common transcription factors and main effect genes in the three main effector cells during sorafenib-induced cardiotoxicity. Augur analysis and WGCNA consistently identified CMs as the most significantly affected cell type. Transcriptional regulatory network analysis highlighted melanocyte-inducing transcription factor (MITF) as a core transcription factor involved in sorafenib-induced cardiotoxicity. Functional validation showed that MITF expression was significantly increased in cardiac tissues and CMs after sorafenib treatment and that knockdown of MITF attenuated sorafenib-mediated CM hypertrophy in vitro.Conclusions By constructing the first snRNA-seq data for sorafenib-induced cardiotoxicity, we identified MITF as a key transcription factor associated with sorafenib-induced cardiotoxicity. Our study provides potential therapeutic targets for mitigating sorafenib-induced cardiotoxicity.
SN - 1741-7015
ER -