一种基于棕榈酰化的新型分子标签揭示COX6A1在代谢功能障碍相关脂肪性肝病中的关键调控作用

A novel palmitoylation-based molecular signature reveals COX6A1 as a key regulator in metabolic dysfunction-associated steatotic liver disease.

作者信息Tianqi Yu, Zhihao Fang, Yue Cheng, Yongxu Zhou, Yanchao Ji, Chang Liu
PMID41185013
发布时间2025-11-03
DOI10.1186/s12967-025-07253-0

实验完整度

整合多组学分析、机器学习建模、单细胞分析、动物模型(MCD/HFD)和体外细胞实验(HepG2)等多层级验证,包含功能验证和机制探索。

主要模型

HepG2细胞(棕榈酸诱导的脂肪变性模型) MCD饮食小鼠模型 HFD饮食小鼠模型 MASLD患者肝脏组织(单细胞和批量转录组)

重点核对

HepG2细胞中COX6A1 siRNA敲低效率(qRT-PCR和Western blot验证) 棕榈酸处理浓度(200 µM)和处理时间(48小时) MCD饮食喂养8周和HFD饮食喂养20周的小鼠模型 线粒体功能检测(JC-1、MitoSOX)和凋亡检测(Annexin V/PI) 外部验证队列(GSE126848、GSE130970、GSE135251)

摘要

背景:代谢功能障碍相关脂肪性肝病(MASLD)患病率很高,但棕榈酰化、线粒体功能和免疫重塑之间的分子联系仍不清楚。方法:我们整合了四个GEO批量RNA-seq队列(GSE126848、GSE130970、GSE135251、GSE213621)和一个单细胞RNA-seq队列(GSE136103)。表达矩阵被归一化为TPM,进行log转换,并使用ComBat进行批次校正。从GeneCards组装了一个精选的棕榈酰化相关基因(PRG)集,并使用共识聚类(ConsensusClusterPlus)定义MASLD亚型。进行了差异表达(limma;|log2FC|>0.5,P<0.05)、功能富集(GO/KEGG的clusterProfiler)和GSEA。WGCNA鉴定了PRG相关模块;使用机器学习流程(12种算法;113种模型组合)优先考虑枢纽基因,在GSE213621中进行10倍交叉验证,并在GSE126848、GSE130970和GSE135251中进行外部验证。使用xCell、EPIC、MCP-counter、QUANTISEQ、CIBERSORT和TIMER推断免疫浸润。在体外,暴露于棕榈酸的HepG2细胞用COX6A1 siRNA转染;评估了脂质积累(油红O/尼罗红)、线粒体膜电位(JC-1)、线粒体ROS(MitoSOX)、凋亡(Annexin V/PI流式细胞术)以及COX6A1/NDUFA4表达(qRT-PCR/Western blot)。性能通过ROC AUC和决策曲线分析进行评估。结果:确定了两个MASLD亚型,Cluster 1和Cluster 2,其中Cluster 1富集代谢途径,Cluster 2显示免疫激活和线粒体代谢途径。确定了枢纽基因,如COX6A1、COX7A2和NDUFA4,具有诊断潜力(AUC>0.75)。单细胞分析揭示了差异免疫细胞浸润,Cluster 2显示免疫活性增加。体外COX6A1敲低减轻了棕榈酸诱导的线粒体功能障碍和凋亡,表明其作为治疗靶点的潜力。结论:我们的研究发现COX6A1与MASLD状态和分子亚型密切相关,并且在独立队列中显示出强大的诊断价值。COX6A1与异常线粒体氧化还原代谢和免疫调节(特别是巨噬细胞信号传导)相关,其调节直接影响肝细胞模型中的脂毒性损伤(mtROS、ΔΨm、凋亡)。总之,这些发现提名COX6A1作为MASLD的机制基础的生物标志物和潜在治疗靶点。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
棕榈酰化相关基因如何影响MASLD的分子亚型和免疫微环境,以及COX6A1在MASLD发病机制中的作用是什么?
核心机制
COX6A1通过调节线粒体呼吸链复合物IV的活性,影响线粒体ROS产生和p53依赖性凋亡信号,从而参与MASLD的代谢-免疫串扰。
主要证据
体外敲低COX6A1减轻棕榈酸诱导的脂质积累、线粒体膜电位去极化、ROS增加和凋亡;WGCNA和机器学习识别COX6A1为枢纽基因,诊断AUC>0.75,并在单细胞和动物模型中验证表达上调。
研究意义
COX6A1作为MASLD的潜在诊断生物标志物和治疗靶点,为针对棕榈酰化和线粒体功能的治疗策略提供了依据。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

数据整合与预处理

整合多个GEO队列,确保数据可比性,为后续分析奠定基础。

从GEO获取4个批量RNA-seq和1个单细胞RNA-seq数据集,进行TPM归一化、log转换和ComBat批次校正。

2

棕榈酰化相关基因的差异表达和功能富集分析

探究PRG在MASLD中的表达变化及其涉及的生物学功能。

使用limma识别MASLD与对照之间的DEGs,进行GO/KEGG富集和GSEA。

3

共识聚类识别MASLD亚型

基于PRG表达谱将MASLD患者分为不同的分子亚型。

使用ConsensusClusterPlus对GSE213621中的PRG表达进行共识聚类。

4

WGCNA和机器学习筛选枢纽基因

识别与PRG显著相关的基因模块并筛选关键枢纽基因,评估诊断价值。

基于ssGSEA评分进行WGCNA分析,结合DEGs,利用12种机器学习算法构建诊断模型并用外部数据集验证。

5

免疫浸润和单细胞分析

比较MASLD亚型的免疫微环境差异,并解析关键基因的细胞类型特异性表达。

使用六种算法评估免疫浸润,对单细胞数据GSE136103进行聚类、细胞注释、轨迹和细胞通讯分析。

6

动物模型验证

在体内验证COX6A1和NDUFA4在MASLD中的表达变化。

建立MCD和HFD小鼠模型,通过组织学染色、血清生化、RT-qPCR、Western blot和免疫组化检测。

7

体外细胞功能实验

验证COX6A1在肝细胞脂质代谢、线粒体功能和凋亡中的作用。

在HepG2细胞中用siRNA敲低COX6A1,棕榈酸处理诱导脂肪变性,检测脂质积累、线粒体膜电位、ROS、凋亡和蛋白表达。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
DMEM培养基Gibco--
胎牛血清Invitrogen--
青霉素-链霉素Gibco--
棕榈酸KunChuang Technology Development Co.--
Lipofectamine 2000转染试剂Invitrogen--
TRIzol试剂Invitrogen--
cDNA合成试剂盒TOYOBO--
SYBR Green预混液Sevenbio--
RIPA裂解液Solarbio--
抗COX6A1抗体AbmartPS03703
抗NDUFA4抗体AbmartPA1787
抗p53抗体HUABIOET1601-13
抗BCL2抗体ZEN-BIOSCIENCE381702
抗BAX抗体ZEN-BIOSCIENCER22708
抗caspase-3抗体ZEN-BIOSCIENCER23315
抗caspase-9抗体ZEN-BIOSCIENCER22844
抗cleaved caspase-3抗体CST9661
抗cleaved caspase-9抗体CST9505
抗GAPDH抗体AbmartP60037
抗β-actin抗体ABclonalAC026
HRP标记二抗ABclonalAS014
JC-1染料BeyotimeC2003S
MitoSOX Red探针BeyotimeS0061S
油红OSolarbio--
尼罗红Solarbio--
Annexin V-APC/PI凋亡检测试剂盒Procell--
BriCyte E6流式细胞仪Mindray--
醋酸铀Sigma--
柠檬酸铅Sigma--
JEOL-1230透射电子显微镜JEOL--
高脂饮食Xiao Shu You Tai--
蛋氨酸胆碱缺乏饮食Xiao Shu You Tai--
总胆固醇检测试剂盒Beijing Solarbio Science & Technology Co., Ltd.--
甘油三酯检测试剂盒Beijing Solarbio Science & Technology Co., Ltd.--
AST检测试剂盒Nanjing Jiancheng Bioengineering Institute--
ALT检测试剂盒Nanjing Jiancheng Bioengineering Institute--
OCT包埋剂Sakura Finetek--
Masson三色染色Solarbio--
兔抗COX6A1抗体AbmartPS03703
兔抗NDUFA4抗体AbmartPA1787
AlphaTSA多重免疫组化试剂盒AlphaxbioAXT37100031
鼠抗ONECUT1抗体AbmartMG520196
DAPI----
Pannoramic MIDI II扫描仪----
SlideViewer软件----

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
数据获取与预处理
GEO数据集ID(GSE126848, GSE130970, GSE135251, GSE213621, GSE136103);归一化方法(TPM, log2);批次校正方法(ComBat)
阅读提示:Methods部分 'Dataset source'
共识聚类
PRG基因集(GeneCards,相关性评分>12,50个基因);聚类参数(k=2,50次重采样,80%样本)
阅读提示:Methods部分 'Consensus clustering analysis'
机器学习模型
训练集GSE213621;10倍交叉验证;113种组合;最佳模型(Stepwise GLM+LDA);外部验证集
阅读提示:Methods部分 'Machine learning-based identification'
细胞实验
HepG2细胞来源(Procell);PA浓度200 µM,处理48小时;siRNA序列(Supplementary Table S3);转染试剂(Lipofectamine 2000)
阅读提示:Methods部分 'Cell culture and small interfering RNA transfection'
动物实验
小鼠品系(C57BL/6,雄性,6周龄);饮食类型(NCD, HFD, MCD);喂养周期(MCD 8周,HFD 20周);分组样本量(NCD n=14, HFD n=7, MCD n=7)
阅读提示:Methods部分 'Animals and treatment'
线粒体功能检测
JC-1浓度2 µM,孵育20分钟;MitoSOX Red浓度1 µM,孵育30分钟;细胞密度(8×10^5或1×10^6每孔)
阅读提示:Methods部分 'Mitochondrial membrane potential' 和 'Measurement of mitochondrial ROS'