使用自动化计算机视觉模型研究腰椎椎旁肌健康与年龄、BMI、性别、体力活动和背痛之间的关联:一项UK Biobank研究

Investigating the associations between lumbar paraspinal muscle health and age, BMI, sex, physical activity, and back pain using an automated computer-vision model: a UK Biobank study

作者信息Evert Onno Wesselink, Annelies Pool-Goudzwaard, Benjamin De Leener, Christine Sze Wan Law, Meredith Blair Fenyo, Gabriella Marie Ello, Michel Willem Coppieters, James Matthew Elliott, Sean Mackey, Kenneth Arnold Weber 2nd
PMID38417587
期刊Spine J
发布时间2024-07
DOI10.1016/j.spinee.2024.02.013

实验完整度

观察性横断面研究,包括MRI图像采集、CNN分割模型训练与验证、统计分析及相关性检验。缺乏功能验证和机制实验。

主要模型

UK Biobank参与者(N=9,564) Dixon脂肪-水MRI腰椎图像

重点核对

是否纳入UK Biobank参与者并排除不符合标准的个体(如无MRI、缺BMI、体力活动、癌症或严重疾病等) MRI扫描参数:1.5T Siemens Aera,VIBE序列,TE1=2.39 ms, TE2=4.77 ms, TR=6.69 ms,翻转角10° CNN训练与测试数据集划分(训练n=100,测试n=30)及匹配标准(性别、年龄、BMI、BP状态) 统计模型:偏相关、ANCOVA,控制变量(年龄、BMI、性别、体力活动)及α=0.05、Bonferroni校正

摘要

背景背景:腰椎椎旁肌健康在背痛(BP)中的作用并不明确。该领域的挑战包括缺乏工具和大型、异质性数据集来探究肌肉健康与BP之间的关联。计算机视觉模型在这一领域具有变革性,能够实现肌肉健康的自动化量化和大型数据集的处理。目的:使用自动化计算机视觉模型,在一个大型、异质性数据集中研究腰椎椎旁肌健康与年龄、性别、BMI、体力活动和BP之间的关联。设计:横断面研究。患者样本:纳入UK Biobank中接受腹部Dixon水脂MRI的参与者(N=9,564)(41.8%女性,平均[SD]年龄:63.5 [7.6]岁,BMI:26.4 [4.1] kg/m2),其中6,953人报告无疼痛,930人急性BP,1,681人慢性BP。结局指标:使用计算机视觉模型自动推导左侧和右侧腰椎多裂肌(LM)、竖脊肌(ES)和腰大肌(PM)从L1至L5椎体水平的肌内脂肪(IMF)和平均横截面积(aCSA)。方法:对每块肌肉进行双尾偏Pearson相关分析,以评估肌肉测量值(IMF和aCSA)与年龄(控制BMI、性别和体力活动)、BMI(控制年龄、性别和体力活动)和体力活动(控制年龄、性别和BMI)之间的关系。使用单因素ANCOVA检验每块肌肉IMF和aCSA的性别差异,同时控制年龄、BMI和体力活动。同样,使用单因素ANCOVA检验无疼痛、急性BP和慢性BP组间差异,同时控制年龄、BMI、性别和体力活动(α=0.05)。结果:女性在所有肌肉中具有更高的IMF(LM平均差[MD]=11.1%,ES MD=10.2%,PM MD=0.3%,p<.001)和更低的aCSA(LM MD=47.6 mm2,ES MD=350.0 mm2,PM MD=321.5 mm2,p<.001)。年龄较高与所有肌肉较高的IMF和较低的aCSA相关(r≥0.232,p<.001),但LM的aCSA除外(r≤0.013,p≥.267)。BMI较高与所有肌肉较高的IMF和aCSA相关(r≥0.174,p<.001)。体力活动较高与所有肌肉较低的IMF和较高的aCSA相关(r≥0.036,p≤.002),但LM的aCSA除外(r≤0.010,p≥.405)。慢性BP患者的IMF高于无疼痛者(IMF MD≤1.6%,p<.001),aCSA低于无疼痛者(aCSA MD≤27.4 mm2,p<.001),且IMF高于急性BP患者(IMF MD≤1.1%,p≤.044)。BP患者与无疼痛者之间的差异并非局限于下腰椎水平,而是广泛分布于整个腰椎。结论:椎旁肌健康与年龄、BMI、性别和体力活动相关,但LM aCSA与年龄和体力活动之间的关联除外。BP患者(慢性>急性)比无疼痛者具有更高的IMF和更低的aCSA。差异并非局部化,而是广泛分布于整个腰椎。在科研或临床环境中解读椎旁肌健康测量值时,应考虑与年龄、BMI、性别和体力活动的关联。关键词:脂肪组织,人工智能,背部肌肉,脂肪浸润,下背痛,磁共振成像

实验结论

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

研究问题
腰椎椎旁肌的健康(肌内脂肪和肌肉大小)是否与年龄、BMI、性别、体力活动和背痛存在关联?
核心机制
文中提出了可能的机制解释,但未验证:年龄相关的纤维类型转变(II型纤维萎缩)和性别差异可能解释年龄和性别与脂肪浸润的关联;体力活动可能通过运动预防年龄相关的脂肪增加;慢性BP相比急性BP的更显著差异可能与长期疼痛和失用引起的肌肉形态变化有关。然而,这些机制属于推测,未通过实验验证。
主要证据
来自UK Biobank的9,564名参与者的Dixon MRI数据,使用自动化计算机视觉模型量化肌内脂肪(IMF)和平均横截面积(aCSA),通过偏相关和ANCOVA分析发现:女性、年龄增加、BMI增加与更高的IMF相关,体力活动与更低IMF相关;慢性BP患者比无疼痛者有更高的IMF和更低的aCSA;腰椎多裂肌的IMF差异不局限于下腰椎,而是广泛分布。
研究意义
研究提供了与年龄、BMI、性别和体力活动相关的椎旁肌健康关联的参考数据,为开发参考值和临床评估提供了基础。

研究路径

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

1

参与者筛选与分组

从UK Biobank中识别符合纳入标准的无疼痛、急性背痛和慢性背痛参与者,确保各组特征具有可比性。

从UK Biobank初始影像评估(应用号67450)中筛选参与者,排除退组、缺BMI或体力活动、癌症、药物滥用、重大疾病、无Dixon MRI、视野不足、水脂交换伪影、图像质量差及分割质量差的参与者。根据疼痛报告(数据字段6159和3571)将参与者分为无疼痛(6,953人)、急性背痛(BP<3个月,930人)和慢性背痛(BP≥3个月,1,681人)。

2

MRI图像采集与肌肉分割

使用Dixon脂肪-水MRI获取腰椎图像,并采用计算机视觉模型自动分割腰椎椎旁肌肉,以量化肌内脂肪和肌肉大小。

使用1.5T Siemens Aera MR扫描仪进行Dixon脂肪-水MRI。由训练有素的评估者(EW)手动分割130名参与者的左右腰椎多裂肌、竖脊肌和腰大肌(L1-L5),作为训练和测试数据。训练了一个基于MONAI框架的改进2D U-Net CNN用于自动分割。

3

CNN模型评估

验证CNN自动分割的准确性、可靠性和速度,确保自动量化肌肉参数的有效性。

评估CNN的分割准确性(Sørensen-Dice指数)、IMF和aCSA的准确性(平均绝对误差MAE和组内相关系数ICC),以及分割速度。

4

统计分析

检验椎旁肌健康(IMF和aCSA)与年龄、BMI、性别、体力活动及背痛状态之间的关联,并评估差异的空间分布。

在无疼痛组中,使用双尾偏相关分析检验IMF和aCSA与年龄、BMI、体力活动的关联(控制其他变量);使用单因素ANCOVA检验性别差异和组间差异(无疼痛、急性BP、慢性BP);在每个轴位切片上使用ANCOVA评估组间差异的分布。

5

空间分布分析

确定椎旁肌IMF和CSA沿腰椎头尾侧的分布,并比较不同背痛组之间的差异是否局限于特定节段。

在每个轴位切片计算IMF和CSA,使用样条插值(interp1d, SciPy)将肌肉长度归一化为0%(最下)至100%(最上),然后在每个归一化切片上使用ANCOVA比较组间差异。

研究方法

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

产品清单

实验环节名称品牌货号
1.5T 西门子 Aera 磁共振扫描仪SiemensSyngo MR D13
MONAI框架----

关键环节

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

环节核对要点
参与者筛选
参与者纳入和排除标准,特别是疼痛报告(Data-Field 6159和3571)、排除条件(癌症、药物滥用等)以及样本量(9,564)
阅读提示:Methods 部分:UK Biobank 和 Inclusion and exclusion criteria
MRI采集与分割
MRI扫描参数(提供商、型号、序列、TR/TE、翻转角、矩阵、分辨率、层厚)以及CNN训练和测试数据划分和参数(训练/测试样本量、匹配变量、训练迭代次数)
阅读提示:Methods 部分:Dixon fat-water MRI 和 CNN training and testing
统计分析与关联评估
统计检验类型(偏相关、ANCOVA)、控制变量、显著性水平(α=0.05)、校正方法(Bonferroni)以及软件版本(R-studio 4.2.2)
阅读提示:Methods 部分:Statistical analyses
空间分布分析
样条插值参数(归一化方法)和每个切片ANCOVA的模型
阅读提示:Methods 部分:Statistical analyses 和 Figure legends (Fig. 5-6)