酮体代谢与钠-葡萄糖协同转运蛋白2抑制剂（SGLT2）在缺血性病因心力衰竭中的疗效关系：一种整合代谢组学与全球疾病负担人群数据、跨越不同遗传背景的精准医学研究方法

《Journal of the American Heart Association》：Ketone Body Metabolism and Sodium Glucose Cotransporter 2 Inhibitor Response in Heart Failure With Ischemic Cause: A Precision Medicine Approach Using Metabolomics Integration With Global Burden of Disease Population Data Across Diverse Genetic Backgrounds

　　

背景：缺血性病因心力衰竭（heart failure with ischemic cause, HFIC）与显著心血管死亡率相关。钠-葡萄糖协同转运蛋白2抑制剂（sodium glucose cotransporter 2 inhibitor, SGLT2抑制剂）虽显示心血管获益，但个体间疗效差异的机制尚不清楚。研究人员调查了基线酮体代谢与SGLT2抑制剂在跨不同遗传背景的HFIC患者中疗效的关系。方法：研究人员分析了来自23个国家（2020–2024年）3847例HFIC患者的代谢组学数据。酮体代谢物（β-羟基丁酸〔β-hydroxybutyrate, β-HB〕、乙酰乙酸〔acetoacetate〕、丙酮〔acetone〕）通过液相色谱-质谱定量。SGLT2抑制剂疗效通过复合终点评估，包括心血管死亡、心力衰竭住院和肾功能下降。分析包括多变量Cox回归、机器学习（随机森林、极端梯度提升〔Extreme Gradient Boosting, XGB〕）和孟德尔随机化（Mendelian randomization, MR），并整合全球疾病负担（Global Burden of Disease, GBD）2021年数据跨越5个遗传祖源组。结果：基线β-HB与SGLT2抑制剂结局呈负相关（r=?0.67，P<0.001）。最低酮体三分位显示更优结局（风险比〔hazard ratio, HR〕0.58〔95% CI 0.51–0.66〕，P<0.001），最高三分位风险升高（HR 1.58〔95% CI 1.39–1.79〕，P<0.001）。东亚人群基线酮体水平较欧洲祖源高34.24%（2.47±0.83对1.84±0.61 mmol/L，P<0.001），且治疗获益减弱。机器学习模型基于基线代谢组学谱预测个体结局的受试者工作特征曲线下面积（area under the receiver operating characteristic curve, AUC-ROC）达0.8245（95% CI 0.8012–0.8478）。结论：基线酮体代谢与HFIC中SGLT2抑制剂疗效密切相关，且存在显著的祖源间差异。代谢组学分析可为精准医学指导治疗决策提供依据，尚需前瞻性验证。

该研究发表于《Journal of the American Heart Association》。研究背景方面，缺血性病因心力衰竭（HFIC）是全球重大公共卫生问题，据全球疾病负担（GBD）2021研究约影响5650万人，占心力衰竭亚型约33.8%，预后差且5年死亡率近49%。病理生理涉及心肌能量代谢失调、慢性炎症和心室重构。钠-葡萄糖协同转运蛋白2抑制剂（SGLT2抑制剂）虽在DAPA-HF、EMPEROR-Reduced等试验中显示心血管死亡和心衰住院降低25%至30%，并被纳入指南导向药物治疗，但仍存在约20%至35%患者疗效不佳，个体反应异质性机制不明。近年证据提示心肌能量代谢特别是酮体利用在心衰病理和药物反应中关键作用；SGLT2抑制剂可通过促进内源性酮生成等机制发挥心脏保护，但基线酮体代谢与疗效关系尚未明确。此外 cardio-kidney-metabolic综合征概念强调代谢功能障碍跨器官交互，酮体代谢异常可能是统一机制。精准医学中代谢组学可提供功能表型读取，且药物反应存在遗传祖源差异（多数试验以欧洲祖源为主），因此研究人员通过开展跨国家、跨祖源的大型代谢组学整合GBD人群数据研究，假设基线酮体代谢谱可预测SGLT2抑制剂在HFIC中的疗效并存在祖源差异，目标包括：刻画不同人群HFIC酮体代谢模式；明确基线代谢谱与SGLT2抑制剂疗效关系；整合个体代谢组学与GBD人群负担数据；开发基于代谢表型的精准治疗选择算法。

关键方法上，研究人员采用数据库挖掘与整合分析：样本队列来源于4个主要心血管代谢组学数据库（MESA、UK Biobank心血管亚研究、4C Study、ELSA-Brasil）合计3847例确诊HFIC患者（年龄≥18岁、有心肌梗死或血运重建史、左室射血分数LVEF≤49%、基线代谢组学及SGLT2抑制剂治疗≥12个月随访），遗传祖源通过全基因组基因分型主成分分析或自我报告结合混合比例验证分为5组（欧洲、东亚、南亚、非洲、西班牙裔/拉丁裔）。酮体代谢物（β-HB、乙酰乙酸、丙酮）及其他247种代谢物用超高效液相色谱串联质谱等定量，跨队列数据经过严格批次校正（ComBat算法）。SGLT2抑制剂疗效终点为复合终点（心血管死亡、心衰住院、eGFR持续下降≥40%）。统计方法包括Cox比例风险模型（按β-HB三分位分层）、机器学习（随机森林、XGBoost等，数据集按7:3分训练/测试，5折分层交叉验证调参）、孟德尔随机化（两样本MR以酮相关SNVs为工具变量）、多变量调整模型、时间依赖协变量分析（6个月重复酮测量）、通路富集（MetaboAnalyst 5.0）、加权整合GBD 2021人群数据（按国别、年龄、性别权重计算人群归因风险）。后续主要结果均按预设小标题归纳。

结果部分：基线特征与人群人口学：研究队列平均66.34±11.28岁，男性71.36%，LVEF均值34.56±7.82%（81.21%为射血分数降低心衰），糖尿病63.84%，慢性肾病（eGFR<60）40.21%；各祖源组年龄、BMI、血压、eGFR、糖尿病患病率等有显著差异（P<0.001），欧洲祖源年龄最大（68.72±10.45岁），南亚祖源最小（62.18±12.34岁），东亚祖源糖尿病比例最高（71.14%），女性患者年龄更大、高血压更多但β-HB水平与男性无显著差异（P=0.152）。

跨遗传祖源的酮体代谢谱：基线酮体浓度祖源间差异显著（P<0.001），β-HB均值东亚最高（2.47±0.83 mmol/L），其次南亚（2.28±0.79）、非洲（2.01±0.72）、西班牙裔（1.95±0.68）、欧洲最低（1.84±0.61）；乙酰乙酸和总酮体也呈类似趋势；乙酰乙酸/β-HB比值东亚最低（0.167±0.048对欧洲0.175±0.052，P=0.003），提示尽管酮生成更高但利用能力可能受损；跨队列内分层分析确认祖源差异主要为生物学差异而非批次效应。

与GBD 2021全球负担数据的整合：GBD 2021显示2021年全球缺血性心脏病患病人数254.3百万（95% UI 221.4–295.5）、死亡9.0百万（95% UI 8.3–9.5）、DALY 188.4百万（95% UI 177.0–198.1），年龄标准化DALY率地区差异近6倍（中亚最高4864.5/10万，西欧最低843.8/10万）。用GBD权重整合个体代谢组数据估算全球约58.73%（95% UI 52.14%–65.32%）HFIC病例存在酮代谢功能障碍（定义为β-HB>2.41 mmol/L且乙酰乙酸/β-HB<0.15），相当于约3247万例；人群归因风险测算酮代谢异常贡献23.45%（95% UI 18.67%–28.92%）HFIC心血管死亡和31.78%（95% UI 26.34%–37.23%）心衰住院。

基线酮体与临床结局的关联：中位随访28.47个月，复合终点发生率32.42%（心血管死亡10.99%、心衰住院19.39%、eGFR降≥40%为2.03%）。按基线β-HB三分位（T1最低0.28–1.67，T2中间1.68–2.41，T3最高2.42–4.89），Kaplan–Meier分析显示事件率T1（21.34%）、T2（30.67%）、T3（43.89%）对数秩P<0.001；完全调整Cox模型以T2为参照，T1的HR=0.58（95% CI 0.51–0.66，P<0.001），T3的HR=1.58（95% CI 1.39–1.79，P<0.001）；各单项终点也呈一致趋势（心血管死亡T1 HR=0.61、T3 HR=1.42；心衰住院T1 HR=0.64、T3 HR=1.41；eGFR降幅≥40%趋势类似但统计显著性较弱）。逆概率加权Kaplan–Meier和调整曲线支持稳健性。

遗传祖源特异的治疗相关结局：SGLT2抑制剂疗效与酮三分位存在显著祖源交互（交互P<0.001）。欧洲祖源T1的HR=0.63（95% CI 0.56–0.71，P<0.001）、T2 HR=0.84（0.75–0.94，P=0.003）、T3 HR=0.88（0.75–1.03，P=0.112），36个月绝对风险降低（ARR）T1为18.67%、T3仅4.23%；东亚祖源T1 HR=0.78（0.66–0.92，P=0.003，ARR 11.34%），T2 HR=0.93（0.81–1.07，P=0.289），T3 HR=1.02（0.84–1.24，P=0.834，ARR ?0.89%提示潜在危害）；南亚、非洲呈中间模式，西班牙裔近似欧洲但样本小精度低。

综合代谢组学谱与通路分析：非应答者相比应答者有187个代谢物显著差异化（FDR q<0.05），前两位为β-HB（非应答/应答倍数1.62，P=3.2×10^?24）和乙酰乙酸（1.52，P=8.7×10^?21）；其余包括支链氨基酸（BCAA：亮氨酸、异亮氨酸、缬氨酸均上调约28%）、中链酰基肉碱（C6–C12上调约35%）、三羧酸（TCA）循环中间物（琥珀酸上调41%、α-酮戊二酸下调23%）、特定神经酰胺、三甲胺N-氧化物（TMAO）、芳香氨基酸等。通路富集显著的有酮体代谢（P=1.2×10^?18，通路影响0.89）、BCAA降解（P=3.4×10^?14，0.76）、TCA循环（P=8.9×10^?12，0.71）、脂肪酸β氧化（P=2.3×10^?10，0.68）、鞘脂代谢（P=5.6×10^?9，0.62）等14条通路。

机器学习预测治疗相关结局：在独立测试集（n=1154）中，极端梯度提升（XGBoost）表现最优，AUC-ROC 0.8245（95% CI 0.8012–0.8478），敏感度78.34%、特异度84.67%、PPV 81.23%、NPV 82.45%、F1 0.7976、Brier 0.134；随机森林次之（AUC 0.8134），神经网络0.7923，支持向量机0.7756，弹性网络0.7589；仅临床人口模型（年龄、性别、祖源、LVEF、NT-proBNP、eGFR、HbA1c、BMI、SBP、糖尿病、SGLT2抑制剂种类）AUC仅0.7198，显著低于代谢组学模型（DeLong P<0.001）。校准曲线良好（Hosmer-Lemeshow χ²=8.45，P=0.391），决策曲线优于全治/全不治策略。XGBoost特征重要性前10为β-HB（100归一化）、乙酰乙酸（87.34）、亮氨酸（72.56）、辛酰肉碱C8（68.92）、神经酰胺C18:0（64.78）、琥珀酸（61.23）、NT-proBNP（58.67）、祖源主成分1（55.89）、LVEF（53.45）、TMAO（51.23）。留一队列交叉验证平均AUC 0.7893支持泛化性。

酮代谢相关基因变异与治疗结局：候选基因分析发现酮代谢基因（HMGCS2、BDH1、OXCT1、ACAT1、SLC16A1等）中23个SNVs与基线β-HB或结局显著相关（P<5×10^?4）；最强为HMGCS2 rs2267623（效应等位G，欧洲EAF 0.18、东亚EAF 0.35），每G等位升β-HB 0.34 mmol/L（P=2.3×10^?8），且弱化治疗获益（HR 1.21〔1.12–1.35〕，P=4.7×10^?5）；其他涉及转运（ABCB1、SLCO1B1）、脂质/糖代谢（PNPLA3、HMGCR、APOA5等）的SNVs也显著。酮代谢多基因风险评分（89个全基因组显著变异）解释基线β-HB方差18.7%（P<0.001），最高五分位比最低五分位β-HB高41.23%、治疗相对获益降28.67%（HR 1.34〔1.18–1.52〕，P<0.001）。祖源主成分1（东欧亚vs欧洲）解释β-HB方差12.34%（P<0.001）；西班牙裔 admixture定位发现染色体1q42、11q13、15q24有祖源特异性酮代谢/结局效应。

孟德尔随机化分析因果关系：两样本MR以34个β-HB相关SNVs为工具变量，逆方差加权估计遗传预测的β-HB每升1 mmol/L关联心衰住院OR=1.47（95% CI 1.21–1.79，P=1.2×10^?4）、心血管死亡OR=1.38（1.14–1.67，P=8.9×10^?4）；MR-Egger截距P=0.234、加权中位数一致、MR-PRESSO无离群、逐一剔除稳定，提示极小水平多效性，支持升高酮水平潜在因果增加不良结局；反向MR（遗传心衰工具→酮）不显著（OR 1.03〔0.91–1.16〕，P=0.634），支持方向为酮升高→不良结局而非反向。

亚组与敏感性分析：酮-疗效负关联在各亚组一致：LVEF≤35% vs >35%（交互P=0.423）、糖尿病有无（P=0.567）、eGFR分层（P=0.289）、NT-proBNP四分位（P=0.378）、年龄<65 vs ≥65（P=0.512）、性别（交互P=0.634，男性T1 vs T2 HR=0.57、女性HR=0.61）；但HbA1c存在交互（P=0.012），HbA1c<7%时关联弱（每mmol/L β-HB HR=1.18〔1.06–1.32〕），HbA1c≥7%时强（HR=1.43〔1.28–1.60〕）。时间依赖分析：有6个月重复酮数据（n=2156）显示持续高酮（>2.41两次）HR=2.34（1.89–2.89，P<0.001）比持续低酮高；SGLT2抑制剂启动后酮下降>20%者结局更好（HR=0.71〔0.62–0.82〕，P<0.001）；酮变化幅度连续预测也显著。排除极端BMI、晚期CKD、近期DKA、随访<12月等敏感性分析均一致；逆概率加权、多重插补与完整病例结果一致；比较排除缺失数据者仅年龄稍大、晚期CKD更多但核心变量无差异且加权分析一致。

讨论部分研究人员总结：本研究整合代谢组学与GBD 2021数据揭示：(1)基线酮体代谢存在显著祖源差异（东亚β-HB较欧洲高约34%）；(2)高基线酮体（最高三分位）关联SGLT2抑制剂疗效减弱（HR 1.58 vs 中间三分位）；(3)全球估测58.73% HFIC病例存在酮代谢功能障碍，贡献近1/4心血管死亡和近1/3心衰住院；(4)代谢组学+XGBoost预测个体疗效AUC-ROC>0.82，支持精准医学应用。衰竭心脏代谢重塑表现为能量低效和底物僵化；高基线酮可能反映心肌缺血/梗死后线粒体氧化受限致酮利用障碍、胰岛素抵抗促脂解/酮前体而外周摄取受阻、神经激素激活促肝酮生成；乙酰乙酸/β-HB降低支持利用受损而非单纯生成增多。祖源差异可由酮代谢酶基因频率不同（如HMGCS2 rs2267623东亚EAF更高）、饮食模式（东亚传统饮食低脂高碳水可能影响酮设定点）、体成分（东亚低BMI但可能潜藏肌肉减少性酮症）等解释；需注意东亚高基线酮部分可能标记亚临床恶液质/肌肉减少，且不同SGLT2抑制剂（本研究含empagliflozin、dapagliflozin、canagliflozin、ertugliflozin）间酮效应异质性未细分。SGLT2抑制剂促酮生成机制明确，但本研究提示并非“酮越高越好”，已达高基线者进一步酮生成可能收益递减甚至有害（天花板效应），需代谢灵活性双向可调。GBD整合显示酮代谢异常为大的可干预靶点但地区负担差异大需因地制宜。机器学习模型AUC 0.82超临床模型，关键特征酮、乙酰乙酸、BCAA、酰基肉碱、神经酰胺等，未来可用于分诊：低酮/有利代谢谱优先SGLT2抑制剂；高酮/高BCAA/神经酰胺等可能需联合生活方式、线粒体靶向、或早考虑器械/高级疗法；系列酮监测（如每3–6个月）可识别无应答调整治疗。遗传学MR支持酮升高因果恶化心衰结局（挑战“酮均为超级燃料”简单观点），短时急性酮可能适应而慢性高酮反映失适应；多基因评分解释有限方差（18.7%）说明环境/生活方式可修饰。局限性：观察性数据挖掘不能确立因果（虽有MR支持）、空腹单次代谢组未捕获餐后/运动应激、某些祖源样本偏小、跨来源批次虽校正但异质性难全消、随访最长4年、男性为主限性别推断、GBD本身有报告偏差、SGLT2抑制剂种类合并分剂未单独分析、缺失数据排除虽检验但潜在选择偏倚不能完全排除。未来需前瞻性随机试验验证代谢组学指导vs标准治疗（可自适应设计）、机制实验（细胞/动物）剖析酮通路与SGLT2抑制剂互作、多组学（蛋白、脂质、转录、表观）整合提升预测、床旁快速代谢检测开发、卫生经济学评估成本效果。

结论部分原文总结为：这项整合个体水平代谢组学与人群水平GBD 2021数据的综合研究确定酮体代谢是与HFIC中SGLT2抑制剂治疗相关结局密切相关的因素，存在显著族裔间差异。高基线酮体关联跨不同遗传背景的治疗相关结局减弱，机器学习算法在结局分类中AUC-ROC达0.8245。这些发现支持实施含代谢组学分析的精准医学方法优化日益扩大的cardio-kidney-metabolic综合征人群治疗选择。未来研究应聚焦前瞻性验证、机制阐明和卫生经济学评价以转化见解为改善患者结局。