《iScience》：Metabolic exhaustion underlies sarcopenic obesity susceptibility in aging mice as revealed by multi-omics factor analysis
编辑推荐：
肌少性肥胖（sarcopenic obesity, SO）影响高达40%的老年肥胖人群，然而许多人在相当代谢负担下仍保持肌肉健康。这种异质性的分子基础尚不明确。研究人员将12月龄雌性小鼠分别置于高脂饮食（high-fat diet, HFD）或标准饮食（sta
肌少性肥胖（sarcopenic obesity, SO）影响高达40%的老年肥胖人群，然而许多人在相当代谢负担下仍保持肌肉健康。这种异质性的分子基础尚不明确。研究人员将12月龄雌性小鼠分别置于高脂饮食（high-fat diet, HFD）或标准饮食（standard diet, SD）条件下饲养7个月。通过多组学因子分析（multi-omics factor analysis, MOFA）整合肌肉转录组学和代谢组学数据，识别解释表型变异的潜在因子。尽管肥胖发展一致，但HFD喂养的衰老小鼠在肌肉结局方面表现出显著异质性。因子3（Factor 3）独特地区分了饮食组别，并能独立于肥胖严重程度预测肌肉功能障碍。该因子揭示了代谢耗竭状态：脂质过载诱导的脂肪酸氧化（fatty acid oxidation, FAO）基因和氧化标志物上调、长链酰基肉碱（long-chain acylcarnitines, LCACs）积累以及游离左旋肉碱（free L-carnitine）耗竭，表明β-氧化不完全。p53、TNF/NFκB及氧化应激通路的激活指示了全面的代谢失代偿。血浆乙酰肉碱（acetylcarnitine）和泛醇（pantothenol）区分肌少症表型的曲线下面积（area under curve, AUC）达到0.91。个体对SO的易感性反映的是差异性的代谢耗竭而非肥胖严重程度，且可通过循环生物标志物进行检测。
研究背景与科学问题
肌少性肥胖（sarcopenic obesity, SO）——即脂肪过多与肌肉质量和力量下降并存——是威胁健康老龄化的重要风险因素，在60岁以上肥胖人群中患病率可达40%-50%。然而，相当一部分老年肥胖者在相似代谢负荷下却从未发生肌肉退化。无论是肥胖持续时间、严重程度还是传统风险因素均无法解释这一异质性，这构成了关键的认知空白。女性因绝经后雌激素缺乏导致脂肪堆积和肌肉分解加速而尤为脆弱。目前的研究模型难以复现这种复杂性：年轻啮齿动物缺乏塑造个体反应的累积代谢与表观遗传多样性，而遗传性肥胖模型则绕过了营养过剩的渐进适应过程。这些标准化模型未能捕捉SO发展的“二次打击”现象——年龄相关的衰退形成易损基质，代谢应激在此基质上可能产生异质性效应。
研究设计与主要发现
为探究可比代谢应激下肌肉功能障碍的差异性易感性及其分子基础，研究人员采用二次打击实验模型，将代谢挑战叠加于衰老生理之上。12月龄雌性C57BL/6N小鼠（相当于中年人类）在HFD或SD下饲养7个月，跨越其向老年期的转变。通过多组学因子分析（MOFA）整合肌肉转录组学和代谢组学数据，研究人员发现：尽管肥胖发展一致，但HFD喂养的衰老小鼠在肌肉结局方面表现出显著异质性。因子3（Factor 3）独特地区分了饮食组别，并能独立于肥胖严重程度预测肌肉功能障碍。该因子揭示了代谢耗竭状态：脂质过载诱导的FAO基因和氧化标志物上调、LCACs积累以及游离L-肉碱耗竭，表明β-氧化不完全。p53、TNF/NFκB及氧化应激通路的激活指示了全面的代谢失代偿。血浆乙酰肉碱和泛醇区分肌少症表型的AUC达0.91。个体对SO的易感性反映的是差异性的代谢耗竭而非肥胖严重程度，且可通过循环生物标志物进行检测。
关键技术方法
研究人员采用了以下主要技术方法：1）二次打击小鼠模型：将12月龄雌性C57BL/6N小鼠（来自Charles River Laboratories）暴露于HFD（35.5%脂肪）或SD（10%脂肪）7个月，分别在1、3、7个月处死取材；2）多组学因子分析（MOFA）：整合骨骼肌转录组学（RNA-seq，5,000个最可变基因）和代谢组学（UHPLC-MS，251种代谢物）数据，通过25次独立初始化选择最高证据下界模型，识别5个潜在因子；3）蛋白质-蛋白质相互作用网络分析：利用STRING数据库（置信度>0.7）结合马尔可夫聚类算法（MCL，膨胀参数2.5）鉴定功能模块；4）因果中介分析：评估因子3介导HFD对肌肉表型效应的比例；5）血浆-肌肉代谢物相关性分析与受试者工作特征（ROC）曲线评估生物标志物的判别能力。
研究结果
衰老小鼠在HFD下加速但异质性发展肌少症
通过7个月干预，HFD显著改变了衰老小鼠的身体成分和肌肉结构。尽管体重和脂肪量百分比从第1个月开始增加，但肌肉反应更为异质。瘦体重百分比从第3个月起显著低于SD对照组。MRI证实第3个月腿部肌肉体积减少，而组织学显示肌纤维萎缩仅在第7个月才明显。值得注意的是，肌肉功能退化早于可检测的质量变化：跑步距离在整个研究期间进行性受损，握力在HFD起始1个月内下降，悬挂线测试从第2个月起降低。功能缺陷比可检测的质量损失早2个月、比肌纤维萎缩早6个月。HFD小鼠在肌肉参数上表现出显著个体间变异（总体稳健变异系数31.3% vs. SD对照14.7%），而脂肪量百分比变异较低（15.4% vs. 31.0%）且保持稳定。脂肪量百分比与任何肌肉参数均无显著相关，表明肥胖严重程度本身不能预测肌肉功能障碍。表型分类显示HFD小鼠中29%保持非肌少症，40%发展为功能性前肌少症，17%进展为肌少症，9%成为重度肌少症。
潜在因子预测代谢应激下的肌肉功能障碍严重程度
MOFA将多组学数据分解为五个潜在因子，共同解释32.99%的转录组变异和13.03%的代谢组变异。因子3独特地在所有时间点分离HFD与SD组，且值保持稳定。线性回归分析显示只有因子3显著预测肌纤维横截面积（cross-sectional area, CSA）和跑步距离。在HFD小鼠中，因子3独立于时间点强预测CSA和跑步距离，而在SD小鼠中无显著关联。因果中介分析显示，对于跑步距离，因子3呈现经典部分中介作用，间接效应占总HFD效应的58.4%。对于CSA，因子3介导HFD诱导的肌肉萎缩，同时掩盖了HFD潜在的合成代谢效应。因子3值随肌少症表型严重程度逐步增加，从正常小鼠的负值到重度肌少症小鼠的高正值。
因子3揭示代谢重编程下的代谢耗竭
因子3捕获了一种低效的代谢转变，表现为脂肪酸代谢增强但不完全脂质氧化的证据。脂质代谢基因主导正载荷（如Fabp3 0.789、Lpl 0.763、Acadl 0.683），伴随草酰琥珀酸和LCACs积累，指示FAO不完全。同时葡萄糖代谢受抑制，游离L-肉碱、乙酰肉碱和抗氧化剂（GSH/GSSG）耗竭。两个AMPK调节亚基（Prkag2、Prkab2）呈负载荷。LCAC/游离肉碱和LCAC/短链酰基肉碱比值在HFD小鼠中升高2.8倍和2.9倍，而乙酰肉碱/LCAC比值降低2.9倍，且这些比值与因子3相关。基因集富集分析识别188条显著改变通路，包括脂质分解代谢、PPAR信号、氧化磷酸化、氧化应激反应和脂肪生成上调，以及碳水化合物代谢下调。蛋白互作网络揭示两大功能模块：脂质代谢模块（44个基因，主要为正向调控）和肌肉结构与功能模块（19个基因，双向调控）。后者中快缩糖酵解标志物（Myh4、Actn3、Mybpc2等）呈负载荷，而氧化性成分（Myh1、Myh2、慢缩/氧化纤维蛋白）呈正载荷。
血浆代谢组学揭示肌肉因子3特征的循环生物标志物
血浆与骨骼肌间240种共有代谢物的交叉组织相关性分析显示广泛代谢协调，氨基酸、酰基肉碱和三羧酸（tricarboxylic acid, TCA）循环中间体尤为突出。33种血浆代谢物与因子3显著相关，其中三种候选标志物|r|>0.5：泛醇（r 0.56）、乙酰肉碱（r -0.53）和草酰琥珀酸（r 0.51）。泛醇从正常小鼠到重度肌少症小鼠逐步升高，乙酰肉碱呈相反递减模式。逻辑回归分析显示乙酰肉碱对前肌少症和肌少症均有保护效应，而泛醇与肌少症风险增加相关。ROC分析中，乙酰肉碱鉴别前肌少症的AUC为0.79、肌少症为0.83，泛醇鉴别肌少症的AUC为0.92。两者联合模型对前肌少症AUC为0.81，对肌少症AUC达0.91。
讨论与结论
本研究通过二次打击模型成功复现了SO的临床异质性——当代谢应激与衰老相交汇时，个体肌肉命运的分化源于定义明确的分子程序而非随机变异。MOFA的应用代表了方法论上的重要进步，其无监督发现的因子3作为唯一预测因子，证明SO的表型异质性反映了一个协调性代谢耗竭程序的差异性激活。该程序符合以下序列：慢性脂质超载驱动PPAR介导的FAO基因上调，但氧化通量无法匹配底物进入，导致LCACs积累和不完全β-氧化；LCACs积累螯合游离L-肉碱，耗竭线粒体脂质输入所需的肉碱池；由此产生的电子流障碍与抗氧化池耗竭共同造成氧化还原瓶颈，触发下游p53、TNF/NFκB和HMOX1等应激通路。因子3值独立于肥胖严重程度或持续时间解释了从韧性到重度肌少症的连续谱系。差异性的中介模式提供了额外见解：耐力损伤的经典中介作用确认代谢功能障碍直接转化为功能缺陷，而肌肉质量的抑制效应揭示了竞争过程——肥胖的合成代谢驱动力被因子3介导的分解代谢所抵消。因子3载荷中AMPK调节亚基的抑制与PPAR信号的强劲激活提示能量感知与脂质代谢转录调控之间的解耦。在结构层面，因子3捕获了肌纤维特异性基因的双向调控——快缩糖酵解标志物受抑而氧化性成分上调——表明活跃的重塑而非被动退化。这些发现将代谢功能障碍置于最早阶段，先于功能损伤和后续结构恶化，建立了治疗窗口。血浆乙酰肉碱和泛醇作为因子3特征的生物标志物实现了临床转化潜力，两者联合模型对肌少症的判别AUC达0.91。总之，本研究证明个体对SO的易感性反映的是一个代谢耗竭程序的激活，而非随机变异或肥胖严重程度。该系统的框架为解析其他复杂代谢疾病提供了模板，最终实现精准医学策略。
