《Scientific Reports》：Metabolomics-based central obesity identifies subtypes musculoskeletal disorders management
编辑推荐：
该研究旨在利用核磁共振（NMR）代谢组学数据识别中枢性肥胖的新型代谢亚型，并评估其与肌少症、椎间盘疾病、骨质疏松和类风湿关节炎四种经典肌肉骨骼疾病发病风险的前瞻性关联。研究纳入英国生物样本库（UK Biobank）247,107名社区居住参与者，中位随访13.
该研究旨在利用核磁共振（NMR）代谢组学数据识别中枢性肥胖的新型代谢亚型，并评估其与肌少症、椎间盘疾病、骨质疏松和类风湿关节炎四种经典肌肉骨骼疾病发病风险的前瞻性关联。研究纳入英国生物样本库（UK Biobank）247,107名社区居住参与者，中位随访13.8年。采用分类提升（Catboost）算法从251种血浆NMR代谢物中筛选对中枢性肥胖判别力最高的代谢特征。通过Shapley值排序，识别出13种关键代谢物。随后应用主成分分析（PCA）和K-means聚类构建肥胖代谢亚型。采用Cox比例风险模型评估这些亚型与新发肌肉骨骼疾病风险之间的关联。研究识别出4种稳定的中枢性肥胖代谢亚型。亚型III以糖蛋白乙酰基（Glycoprotein Acetyls）、缬氨酸（valine）和大高密度脂蛋白（large HDL）中磷脂/脂质水平显著升高为特征，而亚型II表现为白蛋白（albumin）和谷氨酰胺（glutamine）水平显著降低。多变量调整后，与非肥胖人群相比，亚型II和III发生肌少症（风险比，HR=1.43，95%置信区间（CI）1.28–1.61）和椎间盘疾病（HR=1.31，95%CI 1.19–1.45）的风险显著更高。各亚型之间在骨质疏松或类风湿关节炎风险方面未观察到显著差异。基于代谢组学的中枢性肥胖分类识别出具有显著不同肌肉骨骼疾病风险的异质性亚组。高炎症-高分解代谢代谢表型代表肌少症和椎间盘疾病的独立高风险组，为精准预防和管理提供了新的代谢分类框架。
中枢性肥胖已成为全球公共卫生危机，患病率在过去四十年持续上升，目前影响超过三分之一全球成年人群。该代谢异常以 visceral adipose tissue 过度积聚为特征，不仅显著增加心血管疾病和2型糖尿病风险，还与肌肉骨骼系统退行性病变密切相关。骨质疏松、类风湿关节炎、椎间盘疾病和肌少症作为典型肌肉骨骼疾病，共同构成中老年人群残疾的主要原因，并带来显著社会经济负担。然而，其潜在分子机制仍 largely elusive，这一知识缺口使传统基于人体测量学的肥胖分类难以指导肌肉骨骼疾病的精准预防和管理。NMR波谱可同时定量数百种循环代谢物，提供系统代谢稳态全景；既往研究已发现支链氨基酸、糖蛋白乙酰基和脂蛋白亚类分布等代谢特征与中枢性肥胖严重程度相关，并参与骨代谢、软骨基质降解和骨骼肌蛋白稳态调控。例如，缬氨酸等支链氨基酸升高可能通过激活哺乳动物雷帕霉素靶蛋白（mTOR）信号通路并抑制自噬加速肌少症进展；糖蛋白乙酰基作为慢性低度炎症整合标志物，已被确立为类风湿关节炎和椎间盘退变的独立危险因素。但多数现有研究局限于单一代谢物或单个代谢通路关联，未能系统刻画肥胖相关代谢失调的复杂表型异质性。因此，研究人员旨在利用无监督聚类算法构建新型中枢性肥胖代谢组学亚型，并系统评估各亚型与四种典型肌肉骨骼疾病发病风险的关联；研究假设该数据驱动代谢表型策略可揭示被传统肥胖分类掩盖的异质性亚组，从而为肌肉骨骼疾病早期风险识别和个体化干预提供新的代谢分类体系。
研究人员基于UK Biobank 247,107名社区居住参与者、中位随访13.8年的前瞻性队列，利用NMR代谢组学、Catboost结合Shapley加性解释（SHAP）值、主成分分析（PCA）与K-means聚类，以及Cox比例风险模型和Kaplan-Meier曲线，识别出4种中枢性肥胖代谢亚型；亚型II和III与肌少症和椎间盘疾病风险显著升高相关，且独立于年龄、性别、社会经济地位和传统心血管危险因素；骨质疏松和类风湿关节炎在各亚型间无显著差异。该研究提示基于代谢组学的肥胖精细分型可捕捉传统腰围指标未反映的异质性风险信息，为肌少症和椎间盘疾病早期预警及代谢表型导向的精准预防提供新框架。
关键技术方法：研究样本来源于UK Biobank，共纳入247,107名社区居住参与者，中位随访13.8年；采用NMR代谢组学检测约280,000份EDTA血浆样本，获得251项代谢指标，其中170项为绝对定量值、81项为计算浓度比值；以中枢性肥胖状态为输出标签构建Catboost分类模型，按70%训练集和30%验证集划分，五折交叉验证优化超参数，并基于SHAP值排序筛选13种判别力最强代谢物；对13种代谢物行PCA后采用K-means聚类，综合轮廓分析与50次bootstrap重采样稳定性确定k=4，平均调整兰德指数（ARI）为0.9550（SD 0.0217）、平均标准化互信息（NMI）为0.9326（SD 0.0258）；采用Cox比例风险模型和Kaplan-Meier累积风险曲线评估结局；结局由国际疾病分类第十版（ICD-10）诊断编码确定，基线患病者从相应分析中排除；中枢性肥胖按腰围定义，男性>102 cm、女性>88 cm；缺失率低于20%的协变量采用链式方程多重插补。
研究结果如下。
研究人群基线特征：从UK Biobank初始502,401人中，247,107人符合纳入标准并进入分析，其中84,733人被归类为中枢性肥胖。基线描述显示，中枢性肥胖个体更年长、女性比例更高、Townsend剥夺指数更高、非白人比例更高且教育程度更低；其Glycoprotein Acetyls、Large HDL中磷脂/脂质、酪氨酸、缬氨酸和Very Small VLDL中磷脂/脂质水平升高。
中枢性肥胖相关关键特征与新型亚型识别：对251项代谢物数据进行Z-score标准化后采用Catboost算法分析，SHAP分析对最佳训练模型计算SHAP值并排序13个影响最大的特征；Glycoprotein Acetyls为最重要代谢物，其次为Large HDL中磷脂/脂质；SHAP蜂群图进一步确认Glycoprotein Acetyls为主要特征；模型预测性能稳健，训练集曲线下面积（AUC）为0.80，验证集AUC为0.772。对13个所选代谢特征行PCA后进行K-means聚类，最佳聚类数为4；聚类方案稳定性高，ARI为0.955、NMI为0.933；随后将中枢性肥胖参与者分为4个不同亚型进一步分析。
新型中枢性肥胖亚型的基线与代谢特征：亚型I包含21,184人，白蛋白、酪氨酸和谷氨酰胺水平最高；亚型II包含21,265人，以白蛋白和谷氨酰胺水平最低为特征；亚型III包含12,276人，Glycoprotein Acetyls、Large HDL中磷脂/脂质、肌酐、Large HDL中甘油三酯和缬氨酸水平最高；亚型IV包含30,008人，L-HDL中胆固醇酯（CE）/脂质、亚油酸/脂肪酸比值、DHA/总脂肪酸比值及VS-VLDL中游离胆固醇/脂质水平最高。亚型III和IV参与者更年长、白人比例更高，并伴更高收缩压（SBP）水平。
新型肥胖亚型与慢性疾病的关联：采用多变量Cox比例风险模型，调整年龄、性别、Townsend指数、族裔、教育、当前吸烟和饮酒、SBP及舒张压（DBP）后，多数结果仍具统计学显著性。以非肥胖组为参照（HR=1.0），亚型II和III参与者发生肌少症和椎间盘疾病的风险比显著高于亚型I；与非肥胖人群相比，亚型II和III发生肌少症（HR=1.43，95%CI 1.28–1.61）和椎间盘疾病（HR=1.31，95%CI 1.19–1.45）风险显著更高。相反，各肥胖亚型与骨质疏松或类风湿关节炎风险未见显著关联。在平均随访13.8年中，Kaplan-Meier曲线显示四种肌肉骨骼疾病在4个新型肥胖亚型及非肥胖参与者间具有不同累积发病轨迹；椎间盘疾病观察到的累积发病率最高，亚型II和III的累积发病轨迹非常相似，并持续显示肌少症和椎间盘疾病累积发病率更高。
讨论部分总结：研究人员利用13种与中枢性肥胖密切相关的NMR代谢物识别出4种代谢谱显著不同的肥胖亚型，并系统阐明其与新发肌肉骨骼疾病风险的差异关联。主要发现为，在13.8年随访中，亚型II和III较亚型I发生肌少症和椎间盘疾病的风险显著更高，且该关联独立于年龄、性别、社会经济地位和传统心血管危险因素；骨质疏松和类风湿关节炎风险在各亚型间无显著差异。这表明基于代谢组学的肥胖精细分类可有效捕捉传统腰围指标未反映的异质性风险信息，为肌少症和椎间盘疾病早期预警提供新的复合生物标志物策略。讨论指出，肌少症发病涉及脂肪组织浸润骨骼肌所致脂毒性和线粒体功能障碍；椎间盘退变与肥胖状态下椎体机械负荷增加及髓核细胞细胞外基质降解酶上调相关。亚型II和III累积发病率显著升高，提示这些亚组存在更明显的脂解、炎症激活和蛋白分解代谢通路；研究结果支持代谢失调本身在肥胖-肌肉骨骼疾病关联中的核心作用，即使在满足同一中枢性肥胖诊断标准者中，不同代谢表型间疾病风险差异可超过30%。亚型III为高炎症-高分解代谢表型，以Glycoprotein Acetyls、缬氨酸和large HDL磷脂同时升高为特征；亚型II为代谢不足表型，以白蛋白和谷氨酰胺显著降低为特征，提示合成代谢储备受损和抗氧化能力下降；二者代谢标志不同，但均趋同于肌少症和椎间盘疾病风险升高，提示过度炎症和代谢代偿不足是两条不同但同样有害的路径。亚型I以白蛋白和谷氨酰胺最高、风险最低，亚型IV呈中间谱。阴性结果可能与以下因素有关：骨质疏松发病主要受雌激素耗竭、钙磷代谢和年龄相关骨重建调控，而非炎症-代谢通路；类风湿关节炎虽涉及慢性炎症，但发病主要由自身免疫过程驱动，可能较少受NMR panel所捕获代谢特征直接调节；UK Biobank中骨质疏松和类风湿关节炎发生率相对较低，即使样本量大也可能限制检测中等效应量的统计效能；队列平均年龄约57岁，可能不足以充分捕获通常于更老年发生的骨质疏松结局。研究优势包括：超过500,000名社区居住个体的前瞻性队列、近14年随访提供充分统计效能；率先将梯度提升机与SHAP值归因整合用于代谢组学特征选择，克服传统高维降维方法常忽视非线性交互的局限；无监督聚类所得代谢表型经多轮交叉验证和敏感性分析显示高稳定性。局限包括：尽管调整众多协变量，残余混杂不能完全排除；UK Biobank参与者以欧洲血统为主，外推至其他族裔需谨慎；肌肉骨骼疾病结局通过医院住院记录和死亡登记确定，轻症或未住院病例可能存在漏诊。
研究结论翻译：总之，研究人员利用大规模人群队列和NMR代谢组学数据，成功 delineated 4种以中枢性肥胖为中心的新型代谢亚型。其中两种亚型被证实携带肌少症和椎间盘疾病的显著残余风险，该风险不能由传统肥胖指标解释。这些发现拓展了对肥胖相关肌肉骨骼疾病异质性的理解，并提示未来预防策略应从“一刀切”减重范式转向由代谢指纹指导的风险分层。具体而言，被归类为亚型II或III、且肌少症和椎间盘疾病风险升高者，可能从代谢靶向策略中潜在获益，例如抗炎饮食模式或蛋白质-能量补充。后续研究应整合多组学数据和重复测量设计，以阐明从特定代谢亚型到肌肉骨骼结局的动态轨迹，并探索针对高炎症-高分解代谢表型的靶向营养和药物干预策略。
