FA阈值选择对骺-干骺端复合体DTI纤维束成像与年化身高速度相关性的敏感性研究

《European Radiology Experimental》：Sensitivity of physeal-metaphyseal tractography correlations with height velocity for fractional anisotropy threshold selection

　　

目的：骺-干骺端复合体（PMC）的扩散张量成像（DTI，一种基于水分子扩散方向性表征组织微结构的磁共振成像方法）纤维束成像可提供骨骼成熟信息。各向异性分数（FA，量化水扩散方向性的指标）最小和最大阈值分别为0.15和1.0常用于重建纤维束；这些值源于高度有序组织，可能排除驱动生长的较无序肥大区和过渡区。降低FA阈值是否更好捕捉生长敏感特征并与年化身高速度相关尚不清楚。材料与方法：45名健康青少年（20名男孩；12–17岁）接受股骨远端PMC的3 T DTI。采用九个先前未测试的FA阈值组合（FA_min：0.05、0.10、0.15与FA_max：0.40、0.50、0.60配对）和常规阈值（0.15–1.0）进行确定性纤维束成像。计算提取指标（纤维束FA、表观扩散系数[ADC，反映水扩散平均速率的指标]、体积）与年化身高速度之间的偏Pearson相关，并控制Tanner分期（Tanner stage，青春期发育分期）。结果：纤维束ADC在所有阈值下与身高速度呈强正相关（r=0.55–0.63）。纤维束FA呈一致负相关（r=-0.38至-0.45）。纤维束体积相关具有阈值依赖性：在较低FA_min阈值下为正（FA_min=0.05时，r=0.33–0.37），但在较高FA_min阈值下减弱或反转（0.15–1.0时，r=-0.26）。结论：纤维束ADC和FA是生长板活动的稳健指标，对FA阈值选择不敏感。纤维束体积与身高速度的关联具有阈值依赖性，较低FA_min阈值可捕捉生长敏感的PMC组织，较低FA_max可排除无活性PMC组织。这些发现强调标准化纤维束成像参数在生长板DTI研究中的重要性。关键点：问题：如何防止常规DTI FA阈值低估儿童PMC中生长敏感的纤维束体积？发现：较低FA阈值准确捕捉纤维束体积与身高速度的相关性，而纤维束ADC和FA在所有阈值下均可靠相关。相关性声明：优化骺-干骺端复合体DTI中的FA阈值可改善微结构成像与身高速度的相关性，确保儿童骨骼生长研究数据的可重复性，增强青少年无创骨骼生长监测，并可能补充传统骨龄X线片。

研究背景与目的：扩散张量成像（DTI，一种测量水分子扩散方向性的磁共振成像方法）可通过各向异性分数（FA，量化水扩散方向性的指标）和表观扩散系数（ADC，反映水扩散平均速率的指标）表征组织微结构，并用于纤维束成像（tractography，基于DTI重建微结构通路的技术）。骺-干骺端复合体（PMC，包含骺板及邻近干骺端的区域）的纤维束成像可反映与生长相关的微结构活动，其纤维束数量、长度和体积等指标依赖FA阈值。文献中PMC纤维束成像常用FA_min=0.15、FA_max=1.0，该阈值旨在兼顾高度有序PMC组织与噪声/各向同性组织排除，但可能排除驱动生长的较无序肥大区和过渡区。阈值选择会影响纤维束重建，存在敏感性与特异性权衡；然而，PMC纤维束成像对生物学生长速度的敏感性是否受FA阈值影响尚未评估。标准化参数是多中心可重复性和临床实施的前提。本研究旨在评估不同FA阈值对PMC DTI及纤维束成像指标与年化身高速度（height velocity，一年内身高变化）相关性的影响。

关键方法概括：本研究为两所三级儿科中心前瞻性队列的方法学子分析，受试者于2020年7月至2022年7月从美国东北部两城市普通儿科诊所招募。45名健康青少年（20名男孩，12–17岁）在基线和12个月随访接受左膝3 T DTI（20个非共线方向，b=0和600 s/mm2，体素12 mm3）。研究人员手动分割股骨远端骺板，用Diffusion Toolkit拟合扩散张量，Trackvis以FACT算法（Fiber Assignment by Continuous Tracking，连续追踪纤维分配算法）进行确定性纤维束成像，测试10种FA阈值组合（FA_min 0.05、0.10、0.15与FA_max 0.40、0.50、0.60，以及常规0.15–1.0），提取纤维束FA、ADC和体积，并与年化身高速度行控制Tanner分期的偏Pearson相关。

研究结果：研究人群基线性别及青春期早晚期分布均匀。纤维束ADC：控制Tanner分期的偏相关分析显示，10种FA阈值下纤维束ADC均与年化身高速度呈强正相关（r=0.55–0.63，名义p<0.001），最大关联见于FA 0.05–0.60（r=0.63，95%CI：0.42–0.78）；FA_min由0.05升至0.15仅轻度衰减（Δr≈0.05），FA_max由0.40至1.00影响较小。纤维束FA：所有阈值下纤维束FA与身高速度呈一致负相关（r=-0.38至-0.45），最大关联见于FA 0.15–0.60（r=-0.45，95%CI：-0.66至-0.18，名义p=0.002）。纤维束体积：控制Tanner分期的偏相关分析显示，纤维束体积与身高速度的关联高度依赖阈值；FA_min=0.05时呈名义正相关，FA 0.05–0.40为r=0.37（95%CI：0.08–0.60，名义p=0.013），FA 0.05–0.50为r=0.34（95%CI：0.05–0.58，名义p=0.021），FA 0.05–0.60为r=0.33（95%CI：0.04–0.57，名义p=0.026）；FA_min=0.10时相关减弱并失去名义显著性（r=0.10–0.20）；FA_min=0.15时相关系数转为负值但未达统计显著，常规FA 0.15–1.00为r=-0.26（95%CI：-0.52–0.03）。代表性高生长速度受试者（8.3 cm/年，早青春期）在FA 0.05–0.40下纤维束体积为25.5 mm3，而常规FA 0.15–1.0下仅6.1 mm3。

讨论部分总结：讨论指出，PMC扩散指标与身高速度相关，且ADC和FA对FA阈值选择稳健；纤维束FA随身高速度增加而降低，ADC随身高速度增加而升高，与兔模型结果互补，提示快速生长时PMC微结构限制较少，ADC可能反映水合增加、细胞外基质周转或肥大软骨细胞区增大。Raya et al显示轴向扩散与胶原纤维组织一致；兔生长板MRI-DTI与组织学和微计算机断层扫描（μCT）比较提示高ADC反映过渡区扩散。低阈值下纤维束体积与身高速度显著相关，而高阈值下无显著相关，提示宽松的FA_min和FA_max可能捕捉更弥漫的生长敏感纤维束并排除较少反映活跃纵向生长的组织；建议未来评估PMC纤维束体积时以FA_min=0.05、FA_max=0.40作为优化起点。低FA_min可能纳入各向同性噪声，但FA_min=0.05时体积与身高速度的相关性提示额外体积为有意义的生物信号而非随机噪声；臂丛研究显示FA_min由0.05升至0.10可使检出纤维束减少20%，提示常规阈值可能低估外周或过渡结构。限制FA_max同样重要，因为接近骺闭合时高FA矿化组织增加，保留最高各向异性体素会偏向成熟组织而非生长活跃区。未来可探索双阈值分割低FA（0.05–0.15）和高FA（0.15–0.60）区域以区分肥大区与增殖区；ADC可作为稳健的骺成熟和生物年龄标志。局限性包括样本量小、两平台（3 T GE Signa Premier和3 T Siemens Magnetom Prisma）关键扩散参数协调但重复/回波时间独立优化且无旅行受试者校准、年龄范围窄（12–17岁）、p值为描述性且未校正多重比较、样本量不足以估计总体效应量；优势为纵向设计及每个受试者时间点采集协议一致。

研究结论部分翻译：PMC DTI指标显著受纤维束成像参数选择影响。尽管ADC和FA在各阈值下保持稳健，纤维束体积与身高速度的关联被常规FA阈值减弱，并在该队列中于宽松FA_min=0.05和FA_max=0.40时保留最佳。常规阈值0.15和1.0对生长相关低各向异性PMC成分敏感性较低，可能排除对骨骼发育关键的较低各向异性过渡区并纳入非骺板组织。这些描述性方法学发现支持标准化报告纤维束成像参数，作为开发儿童肌肉骨骼影像可重复DTI生物标志物的前提；FA_min=0.05和FA_max=0.40被提出作为优化起点而非确定临床标准，尚待在专门有把握估计这些相关性的更大队列中进行推断性确认。