肝细胞气球样变的转录组轨迹揭示MASLD中一个脂毒性不可逆点并启动致癌性调控重连

《Biochemistry and Biophysics Reports》：Transcriptomic trajectory of hepatocyte ballooning reveals a lipotoxic point of no return and initiates oncogenic regulatory rewiring in MASLD

　　

摘要 背景： 肝细胞气球样变是进行性脂肪性肝炎的一个定义性组织学标志，但其分子基础仍知之甚少。 目的： 在MASLD至肝细胞癌谱系中，鉴定肝细胞气球样变背后的分子特征和机制轨迹。 方法： 将来自43名按组织学气球样变评分分层的患者的肝脏微阵列数据，与机器学习特征选择和公开可用的单核RNA测序数据（87,244个肝细胞核）相结合，推导出一个基因特征，并在一个独立的批量转录组队列（GSE130970，n = 75）中进行测试；拟时间、转录因子调控子和细胞–细胞通讯分析表征了其潜在生物学机制，并在棕榈酸处理的Huh7细胞中进行体外验证，以及在510名TCGA-LIHC患者中进行预后检验。 结果： 一个八基因特征（SPP1、MFSD2A、PRAP1、ANXA2、KRT8、MACROH2A2（H2AFY2）、JUN 和 FABP4）表征了气球样变肝细胞，并以AUC为0.768（95% CI 0.654–0.883）区分中重度与无气球样变。其中六个基因构成一个肝细胞气球样变–肝细胞损伤（HBHi）评分，该评分在肝硬化（F4）肝脏的细胞核中最高。拟时间分析解析了一个四阶段分子崩溃过程，在拟时间≈25处达到顶点，以染色质稳定因子MACROH2A2的缺失、KRT8和SPP1的共诱导，以及HNF4A/PPARA驱动的程序被应激因子JUN、ATF3和EGR1所取代为标志。气球样变肝细胞通过FN1–CD44和FN1–(ITGAV + ITGB1)轴传递促纤维化信号至星状细胞和巨噬细胞；KRT8、JUN和SPP1的高表达均与肝细胞癌中较差的生存率相关。 结论： 肝细胞气球样变是一个转录活跃的、逐步进行的过程，具有一个脂毒性不可逆点，其预后相关性延伸至肝细胞癌。

论文解读：肝细胞气球样变的转录组轨迹揭示MASLD中的脂毒性不可逆点与致癌性调控重连

一、研究背景与科学问题

代谢功能障碍相关脂肪性肝病（MASLD）是全球最常见的慢性肝病，一项涵盖92项人群研究、超过930万个体的荟萃分析估计其全球患病率为30.05%，且从1990–2006年的25.26%上升至2016–2019年的38.00%。部分单纯性脂肪变性患者会进展为代谢功能障碍相关脂肪性肝炎（MASH），后者可能进一步发展为肝硬化和肝细胞癌（HCC）。肝细胞气球样变（HB）是MASH的一个核心特征，指肝细胞的一种退化形式，表现为细胞增大、肿胀、变圆以及特征性的网状细胞质。气球样变细胞通常比正常肝细胞更大更圆，在苏木精–伊红染色下呈现苍白的网状细胞质，即“蜘蛛网”样外观。组织病理学研究已显示脂肪变性肝细胞和窦周纤维化区域与气球样变肝细胞共定位，提示HB参与MASH进展。肝细胞气球样变的存在与更严重的肝病以及更高的肝脏、心血管和慢性肾脏并发症发生率相关，且活检显示气球样变、Mallory–Denk小体或纤维化的患者肝脏相关死亡频率更高。

尽管肝细胞气球样变与炎症活动、纤维化进展和不良临床结局相关，是MASLD–MASH连续谱中的一个关键病理转变，但目前两种主要组织学评分系统均仅对气球样变进行半定量分级：NASH临床研究网络NAFLD活动评分（NAS）以0–2分制评分，SAF（脂肪变性–活动–纤维化）系统将气球样变纳入复合活动评分。尽管具有诊断重要性，肝细胞气球样变仍主要依据形态学标准定义，将该表型与恶性进展联系起来的分子机制尚未完全阐明。本研究旨在通过整合批量转录组学、机器学习特征选择和单核RNA测序，重建健康–MASLD–MASH–HCC谱系中肝细胞损伤的分子轨迹，鉴定与肝细胞气球样变相关的转录程序，解析连续的肝细胞疾病轴，并将肝细胞气球样变表征为连接代谢性肝损伤与肝细胞癌的转录调控状态。

二、主要关键技术方法

研究人员首先重新分析了本课题组既往生成的肝脏微阵列数据，经质量过滤后保留43个样本，以软阈值功率5构建加权基因共表达网络分析（WGCNA）共表达网络，并用limma进行差异表达分析（|log?FC| > 0.5，FDR < 0.05）。随后采用LASSO逻辑回归结合10,000次bootstrap迭代进行特征选择与稳定性分析，并在独立批量转录组队列（GSE130970，n = 75）中用岭惩罚逻辑回归模型验证八基因特征，构建列线图并进行决策曲线分析和bootstrap校准。在单核水平，研究人员整合了两个公开可用的单核RNA测序数据集GSE189175和GSE202379（共87,244个肝细胞核），经Seurat v5处理、scDblFinder去双联体、RPCA整合后对肝细胞亚群进行亚聚类，用Monocle 3进行拟时间轨迹推断，用scMetabolism评分KEGG代谢通路活性，用pySCENIC推断转录因子调控子网络，用CellChat分析细胞–细胞通讯。此外，在棕榈酸（PA）处理的Huh7细胞中进行体外验证（RT-qPCR和Western blot），并在510名TCGA-LIHC患者中进行生存分析。

三、研究结果

3.1 从肝脏组织转录组中鉴定肝细胞气球样变基因特征

研究人员分析了43名按组织学气球样变评分（HBS）分类的患者肝脏微阵列数据，其中HBS 0（无气球样变）24例，HBS 1+2（中重度）19例。WGCNA鉴定出MEblack模块（1238个基因），其特征基因与肝细胞气球样变评分的相关性最强（r = 0.57），并与NAFLD活动评分（r = 0.62）、脂肪变性分级（r = 0.58）、小叶炎症（r = 0.49）、纤维化分期（r = 0.46）等均相关。limma差异表达分析鉴定出21个显著差异表达基因，与MEblack模块取交集后获得20个高置信候选基因。LASSO回归经10,000次bootstrap迭代最终选定八个基因：SPP1、MFSD2A、PRAP1、ANXA2、KRT8、MACROH2A2（H2AFY2）、JUN和FABP4。HSPC157（假基因）和C10orf116（ADIRF）虽得分较高但因缺乏肝细胞损伤中的已知功能而被排除。该八基因面板在独立队列GSE130970（n = 75）中经岭回归评估，AUC为0.768（95% CI 0.654–0.883）。模型权重用于构建列线图，将八基因表达谱转换为1+2级气球样变的风险评分，其中贡献最大的三个基因为MACROH2A2、FABP4和ANXA2。决策曲线分析显示在0至0.92阈值范围内相对于治疗全部和不治疗策略具有更优的临床净获益，bootstrap校准（B = 200）的平均绝对误差为0.046。

3.2 拟时间轨迹分析解析四阶段分子崩溃

研究人员整合了跨越健康–MASLD–MASH–HCC谱系的两个snRNA-seq数据集GSE202379和GSE189175。六个基因（JUN、ANXA2、KRT8、SPP1、MFSD2A和MACROH2A2）因在两个数据集的超过10%肝细胞核中被检测到且均有肝细胞损伤和MASH相关文献支持而被保留用于后续单细胞分析。FABP4因主要表达于肝脏巨噬细胞和脂肪细胞而被排除，PRAP1因核转录本计数极低而被排除。整合后恢复五种主要肝脏细胞类型，肝细胞亚群再聚类得到11个亚群。Monocle 3轨迹推断将肝细胞沿健康至HCC的拟时间轴排序，疾病分期沿该轴按预期时间顺序分布。HBHi评分与拟时间相关（Spearman ρ = 0.40，p < 2.2 × 10?1?），按纤维化分期分组时HBHi评分在肝硬化（F4）肝脏中最高。沿拟时间轴观察到MFSD2A表达进行性下降，提示肝脏空间分区破坏；JUN和ANXA2短暂升高，与AP-1应激信号和膜损伤反应一致；MACROH2A2在拟时间≈25处急剧下降，此时KRT8和SPP1同时被激活。晚期轨迹肝细胞进入以正常肝细胞身份丧失和结构崩溃为特征的状态。

3.3 沿MASLD谱系鉴定脂毒性临界点

研究人员用scMetabolism对每个细胞核的代谢通路活性进行评分。拟时间早期肝细胞保留正常功能，包括宏量营养素和异生物质代谢。随着疾病沿拟时间轴进展，这些通路协调下降。脂肪酸降解能力在轨迹约第30百分位处急剧下降，恰为上升的HBHi损伤负荷超越该能力的点。活性氧和TNF/NF-κB信号是最早上升的损伤通路，而凋亡和未折叠蛋白反应活性上升较晚并在轨迹最后三分之一达到峰值。

3.4 特征基因在体外被棕榈酸驱动的脂毒性急性诱导

研究人员用150 μM棕榈酸（PA）处理Huh7细胞。RT-qPCR显示JUN、KRT8和SPP1转录本呈时间依赖性诱导：JUN和KRT8在24小时已显著上调（分别为1.6倍，p < 0.01和2.3倍，p < 0.01），48小时进一步升高（3.5倍和5.1倍，均p < 0.001）；SPP1在24小时无变化，仅在48小时显著诱导（2.4倍，p < 0.01）。Western blot证实PA处理细胞中骨桥蛋白（SPP1）和KRT8蛋白表达增加。这些结果表明急性脂毒性应激足以在体外激活肝细胞气球样变基因程序。

3.5 应激转录因子取代代谢调控程序

研究人员用pySCENIC绘制snRNA-seq肝细胞亚群的转录因子调控子活性并将其与HBHi评分关联。出现两组对立的转录因子：活性随HBHi升高的一组以EGR1（r = 0.56）、JUN（r = 0.52）和ATF3（r = 0.45）为首，构成应激驱动臂；活性降低的一组包括HNF4A（r = ?0.40）、PPARA（r = ?0.38）和NR1H4（r = ?0.33），构成崩溃的保护臂。沿拟时间轴，CEBPB、FOS、ATF3和JUN逐渐上升，而HNF4A、NR1H3和RXRA下降。值得注意的是，EGR1转录本计数在各疾病阶段稀疏，但其调控子活性在HCC中急剧上升。调控网络分析显示，应激侧ATF3直接靶向JUN（网络中最强调控权重，权重 = 5.2），EGR1也直接调控JUN（权重 = 2.7）；保护侧SOX5维持MACROH2A2（权重 = 3.2），NR1H4维持MFSD2A（权重 = 2.5）。随着HBHi评分升高，两个保护轴逐渐失去活性。

3.6 气球样变肝细胞向周围细胞发送促纤维化信号

研究人员用CellChat建模所有疾病阶段所有细胞类型间的配体–受体通讯。细胞类型比例随疾病进展显著变化，HB_high肝细胞亚群进行性扩增，在晚期NASH和HCC中最显著。从HB_high肝细胞到星状细胞和巨噬细胞的关键配体–受体相互作用包括FN1–CD44和FN1–(ITGAV + ITGB1)，表明活跃的细胞外基质重塑，以及LAMC1–(ITGA1+ITGB1)和LAMC1–(ITGA9+ITGB1)，表明基底膜重编程。

3.7 气球样变相关基因在肝细胞癌中具有预后价值

研究人员分析了TCGA-LIHC队列510名HCC患者的表达和生存数据，按各基因表达中位数二分后以Kaplan–Meier曲线和Cox比例风险回归评估结局。高KRT8表达与较差的总生存相关（log-rank p = 0.038；HR = 1.5，p = 0.04）。高JUN表达同样预测较差的总生存（log-rank p = 0.02；HR = 1.5，p = 0.021）。高SPP1表达携带所有测试基因中最强的预后信号（log-rank p = 1.6 × 10??；HR = 2.2，p = 2.4 × 10??）。相反，高MFSD2A表达趋向于更好的无病生存（log-rank p = 0.1；HR = 0.75）。所有风险比均为单变量分析。

四、讨论与结论

讨论部分指出，本研究将肝细胞气球样变定义为一种转录活跃的动态细胞状态，而非被动的形态学观察。八基因特征在独立队列中达到AUC 0.768的区分能力，虽非决定性但具有临床参考价值。每个基因均有支持的生物学作用：SPP1（骨桥蛋白）由hedgehog通路激活诱导并通过激活肝星状细胞促进脂肪性肝炎中的纤维化进展；KRT8聚集体形成Mallory–Denk小体，是气球样变肝细胞的结构标志；MACROH2A2是此背景下新增的染色质稳定组蛋白变体，其缺失可能允许应激和致癌基因的去抑制；含c-Jun的AP-1二聚体在饮食性脂肪肝模型中调控脂肪性肝炎和PPARγ信号。拟时间轨迹分析确定了四个连续事件：早期MFSD2A高表达标志门静脉周围身份，其缺失是肝脏空间组织开始瓦解的早期指标；随后JUN和ANXA2短暂激活反映AP-1介导的应激反应和膜损伤信号；在拟时间≈25处三个事件同时发生——MACROH2A2急剧下降、KRT8被诱导、SPP1被激活，这一染色质去抑制、细胞骨架损伤和促纤维化信号的巧合定义了“灾难性连接点”（Catastrophic Nexus），标志从可逆损伤状态向结构和促纤维化崩溃的转变；晚期肝细胞进入正常身份丧失状态。调控网络分析确定了保护臂（HNF4A、PPARA、NR1H4）和应激臂（EGR1、JUN、ATF3）的对立，应激臂内ATF3和EGR1直接驱动JUN形成自我强化环路，而保护臂随HBHi升高逐渐失效，这为晚期MASH趋向不可逆提供了机制解释。细胞通讯分析显示重度气球样变肝细胞是所有疾病阶段的主要信号发送者，HB_high亚群随疾病向HCC进展而扩增，FN1–CD44和FN1–ITGB1轴可能代表可独立于气球样变本身解决的治疔靶点。体外实验表明急性脂毒性应激足以激活核心特征基因。TCGA-LIHC中SPP1、KRT8和JUN高表达均与较差总生存相关，提示气球样变转录状态与致癌进展共享分子特征，但这些关联是观察性的且来自HCC患者而非MASH患者。

研究结论： 本研究提供了肝细胞气球样变作为一种进行性、转录活跃状态的分子模型。八基因特征和HBHi评分为客观的、基于转录组学的分期提供了工具，可补充组织学评分。灾难性连接点——以MACROH2A2的同时缺失、KRT8的诱导和SPP1在拟时间≈25处的激活为定义——及其上游驱动因子ATF3和EGR1可能界定了一个治疗窗口。在不可逆纤维化或恶性转化发生前拦截这一转变可能代表MASH治疗的新策略。未来工作应聚焦于在MASH队列中对HBHi评分进行前瞻性验证、确认所鉴定转录因子因果作用的功能研究，以及靶向FN1–CD44和FN1–ITGB1信号轴的干预措施的临床前测试。