灰质生物物理模型的比较系统分析基于公共数据集

《MAGNETIC RESONANCE IN MEDICINE》：Comparative Systematic Analysis of Gray Matter Biophysical Models on a Public Dataset

　　

目的：旨在比较四种灰质扩散模型（NEXI、SMEX、SANDI和SANDIX）的表现。方法：研究人员使用单个公开可用的人类体内数据集，即连接组扩散微结构数据集（CDMD），该数据集采集了两个扩散时间。利用开源工具箱Gray Matter Swiss Knife在26名健康受试者中估计了皮层微结构指标，并评估了拟合优度、解剖模式以及与先前研究的一致性。结果：CDMD数据产生的灰质参数估计值与先前研究在所有四个模型中报告的值一致。NEXI和SMEX产生了相似的皮层解剖模式，在扩散时间之间具有一致的区域分布。拟合优度在不同模型间变化，NEXI显示最佳拟合，其次是SMEX，最后是SANDIX。SANDI参数估计显示对扩散时间选择和拟合算法的强烈依赖性。结论：这项回顾性跨模型分析证明了仅从两个扩散时间估计交换模型的可行性，并强调了生物特异性、模型复杂性和拟合鲁棒性之间的权衡，这些是选择未来临床和研究应用模型时的关键考虑因素。

扩散磁共振成像（dMRI）能够非侵入性地探测脑组织微观结构，为理解灰质（GM）微结构特征提供了独特窗口。近年来，针对GM的生物物理模型迅速发展，包括神经突起交换成像（NEXI）、交换标准模型（SMEX）、体细胞和神经突起密度成像（SANDI）以及体细胞和神经突起密度成像结合交换（SANDIX），它们分别考虑了体细胞贡献、隔室间交换等不同微结构特征。然而，这些模型通常需要多扩散时间采集，且缺乏在公共数据集上的系统比较，这限制了其临床转化潜力。因此，本研究旨在利用公开的CDMD数据集（MGH Connectome Diffusion Microstructure Dataset），比较四种GM模型在仅两个扩散时间下的表现，评估其可行性、可靠性及解剖一致性，为未来模型选择提供依据。该研究发表在《MAGNETIC RESONANCE IN MEDICINE》。

研究人员使用了CDMD数据集，该数据集来自26名健康受试者（36.8±14.6岁，17名女性），在3T Connectom扫描仪上采集，梯度强度300 mT/m，包含两个扩散时间（Δ=19 ms和49 ms），b值范围0.05–6 ms/mm2（Δ=19 ms）和0.2–17.8 ms/mm2（Δ=49 ms），体素分辨率2 mm各向同性，总扫描时间55分钟。采用开源Gray Matter Swiss Knife工具箱，通过非线性最小二乘（NLS）结合L-BFGS-B算法进行参数估计，初始值通过网格搜索确定。对26名受试者的幅度数据应用Rician均值校正（RM），对14名受试者的实值数据则采用复杂去噪。SANDI仅使用短扩散时间（Δ=19 ms）进行拟合，以符合其无交换假设。

**4.1 所有GM生物物理模型的定量比较**
**4.1.1 参数估计总结**
通过森林图比较了皮层带参数估计值。结果显示，所有模型的参数估计在实值与幅度数据去噪策略之间以及与之前研究报告的值均较为一致。但交换时间（τ）在实值数据中显著高于幅度数据（NEXI: 40–70 ms vs. 25–30 ms; SMEX: 40–50 ms vs. 30–40 ms）。神经突起内扩散率（D_n）是所有模型中最不稳健的参数，常触及上界。SANDI的体细胞信号分数（f_s）在NLS拟合下显著低于先前使用随机森林（RF）算法的研究，而神经突起信号分数（f_n）则更高，表明拟合算法对参数估计有显著影响。
**4.1.2 拟合优度**
通过校正Akaike信息准则（AICc）比较模型拟合质量。NEXI获得最佳拟合（AICc = -1.29），其次为SMEX（-1.24），SANDIX最差（-1.14）。SANDI仅用单扩散时间拟合，其AICc（-1.31）因数据点较少而不可直接比较。当用SANDI参数预测另一扩散时间信号时，预测信号与实测信号趋势相反，表明模型无法捕捉多扩散时间依赖。
**4.1.3 群体平均参数图**
参数图显示，SANDIX参数图噪声最大，f_n图在所有模型中噪声最显著。SANDI的f_s图缺乏预期的灰白质对比，而NEXI和SMEX的交换时间（τ）在脑室周围白质区域高于灰质，符合髓鞘轴突渗透性低的特征。
**4.2 NEXI与SMEX的聚焦比较**
**4.2.1 模拟**
通过生成SMEX合成信号（与实验信号相同的脉冲序列参数），并用NEXI和SMEX进行拟合。结果显示，NEXI在参数估计准确性上优于SMEX，尤其在τ和D_n方面，SMEX产生了更大的偏差。这可能归因于SMEX的数值积分求解器精度不足。
**4.2.2 表面图**
NEXI和SMEX的皮层表面图显示几乎相同的解剖趋势，但尺度略有差异。τ在感觉运动皮层（前中央、旁中央区域）以及枕叶和颞叶中较高；f_n在枕叶较高，前颞叶较低；D_n分布与f_n相似；细胞外扩散率（D_e）在枕叶较低，岛叶较高。这些模式与已知的细胞构筑和髓鞘构筑特征一致。

讨论部分指出，噪声建模至关重要，Rician均值校正显著影响τ估计，即使在实值数据中也不能忽略。NEXI在本次实验中表现优于SMEX，主要归因于数值实现差异而非模型本身局限，因为本次采集的脉冲持续时间（δ=8 ms）较短，更接近窄脉冲近似。SANDI由于缺乏交换隔室，在多扩散时间数据中表现不佳；SANDIX因参数过多导致拟合稳定性下降。研究结论部分翻译如下：研究人员回顾性评估并比较了四种先进生物物理模型（NEXI、SMEX、SANDI和SANDIX）在公共扩散MRI数据集CDMD上的表现，使用Gray Matter Swiss Knife包进行拟合。这些模型在拟合优度上有所不同，NEXI最佳，SMEX次之，SANDIX最差。本研究中获得的参数估计与先前研究报告一致。对于交换模型如NEXI和SMEX，研究人员证明使用仅两个扩散时间的采集协议可以获得准确且解剖上有意义的参数估计。关键的是，研究人员强调了显式建模噪声的重要性，即使处理实值数据时也是如此，因为它显著影响交换时间参数估计（τ）。在幅度数据情况下，这意味着在NEXI模型中纳入Rician均值。此外，研究人员显示SANDI参数估计，特别是f_n，对拟合算法的选择敏感，并且由于缺乏交换隔室，该模型无法捕捉多扩散时间依赖性。未来的工作应专注于进一步优化采集参数以减少扫描时间，并使灰质模型能够临床转化到患者群体，同时确保皮层微结构估计的可解释性和可靠性。