白质扩散MRI信号高阶球谐功率的b值依赖符合轴突Stick模型预测：独立TDE测量Da的严格验证

《MAGNETIC RESONANCE IN MEDICINE》：The b-Value Dependence for Higher-Degree Spherical Harmonic Powers of the Diffusion MRI Signal Follows Stick Model Predictions in White Matter

　　

目的：比较白质扩散MRI（dMRI）信号球谐（SH）功率对b值的依赖与轴突内水扩散Stick模型的理论预测。
方法：在3名健康成人中于3 T采集dMRI数据，包含10个b值壳层（b=1000–10 000 s/mm²），每壳层64个方向。对每个b值，将信号展开为球谐（SH），包含所有偶数阶l至l=6，并计算各阶谐波功率。在选定白质感兴趣区（ROI）中，将这些功率与Stick模型理论预测比较；该预测仅依赖两个可调参数：轴突内扩散率（D_a）和每阶振幅因子。采用三重扩散编码（TDE）MRI独立估计D_a，使振幅因子成为唯一拟合参数。
结果：对于b ? 4000 s/mm²（预期轴突外水信号被强烈抑制的阈值），在总大脑白质以及胼胝体膝部、体部和压部中，实验与理论之间均观察到密切一致。这些发现扩展了既往仅考虑l=0谐波功率或将D_a作为拟合参数的研究。
结论：在白质中，对于b ? 4000 s/mm²，扩散MRI信号偶数阶球谐功率（至l=6）的b值依赖可被轴突Stick模型相对较好地描述。

扩散MRI（dMRI）通过测量信号随扩散编码强度（b值）和方向的变化，为脑白质微结构组织提供非侵入性观察窗口。白质dMRI模型常将轴突视为一维不可渗透的“Stick”，并据此预测轴突内水方向平均信号在高b值下随b值增加而下降；该行为已在b值约4000–6000 s/mm²以上的脑白质中得到实验支持。球谐（SH）展开中，方向平均信号对应l=0球谐功率，而Stick模型还对l≥2偶数阶球谐功率给出由轴突内扩散率（D_a）和阶数特异振幅因子控制的b值依赖。既往研究多只检验l=0，或将D_a作为拟合参数、采用文献代表值，高阶球谐预测缺乏严格实验验证。高阶谐波恰是决定纤维取向密度函数（fODF）角变化的关键，对约束球谐反卷积（CSD）、神经突方向离散度与密度成像（NODDI）、纤维球成像（FBI）、球平均技术（SMT）和标准模型成像（SMI）等基于Stick模型的方法具有重要意义；同时，高阶谐波功率对测量噪声和有限方向采样更敏感，可能产生系统偏倚。因此，研究人员需要在独立测量D_a的条件下，逐阶检验l=0、2、4、6的Stick模型预测，并评估噪声与采样稳健性。

研究人员开展了一项技术验证研究：采用三名健康成人（29岁女性、33岁男性、49岁女性）的多壳高b值dMRI数据，计算总白质及胼胝体膝部、体部、压部的偶数阶SH功率，并与由三重扩散编码（TDE）独立估计的D_a所参数化的Stick模型预测比较；同时通过AICc比较固定TDE D_a、共享自由D_a和阶数特异D_a模型，并用蒙特卡洛模拟评估有限方向采样和Rician幅度噪声影响。研究结论是：在高b值区，l=0、2、4的实测b值依赖与Stick模型预测高度一致，l=6一致性较弱但总体支持；固定或共享D_a模型获得几乎全部模型支持，阶数特异D_a支持很少。其意义在于为Stick模型描述白质dMRI信号高阶角结构、进而支持基于SH的fODF估计和纤维追踪提供了更严格验证，同时提示l=6结果需谨慎解释。

关键技术方法不超过250字：研究纳入三名健康成人，数据在南卡罗来纳医科大学机构审查委员会批准方案下采集；主dMRI在3 T扫描仪上采集10个b值壳层（b=1000–10 000 s/mm²）、每壳层64方向，TDE数据在b=4000和4600 s/mm²采集。图像经MP-PCA去噪、Gibbs振铃去除、涡流与运动校正及Rician噪声偏倚校正；在总白质和胼胝体三个亚区ROI中，将各b值归一化方向信号展开至l=6的复数SH，计算偶数阶谐波功率。以TDE逐体素估计D_a并固定于模型，仅拟合每阶振幅；用几何均方根误差（gRMSE）和Lin一致性相关系数（CCC）评价实验与理论一致性；用AICc比较固定、共享和阶数特异D_a模型；并以实际64方向采集几何和经验校准信噪比（SNR）开展蒙特卡洛模拟。

研究结果如下。4.1 谐波功率对b值的依赖（Harmonic-Power Dependence on b-Value）：通过多壳dMRI计算各偶数阶SH功率，并在总白质及胼胝体亚区与固定TDE D_a后的Stick模型预测比较，发现总白质中l=0、2、4在b≥4000 s/mm²与b≥6000 s/mm²均一致较好；l=0的gRMSE为3.45%和2.97%，CCC均为0.987；l=2的gRMSE为3.88%和3.37%，CCC约0.982–0.983；l=4在更严格阈值下gRMSE由3.44%降至1.76%，CCC由0.910升至0.967；l=6一致性较弱且CCC在b≥6000 s/mm²降至0.596。胼胝体体部和压部对l=0、2、4显示强一致，膝部l=4在更严格阈值下改善明显；各区域l=6均较弱且变异较大。结论是l=0、2、4支持预测的高b值依赖，l=6支持较弱。4.2 拟合模型对D_a的依赖（Dependence of the Fitted Model on D_a）：通过逐体素AICc比较三种D_a参数化，发现在b≥5000 s/mm²时固定TDE D_a模型在52.6%有效体素中占优，共享自由D_a模型占46.4%，阶数特异D_a仅1.1%；平均Akaike权重分别为0.465、0.520和0.015，约98.5%模型支持归于固定或共享D_a。共享D_a中位数约2.67 μm²/ms，较TDE估计高约16.9%；b≥6000 s/mm²时固定TDE模型支持增至62.4%。ROI层面共享D_a模型在12个受试者×ROI拟合中占11个。阶数特异D_a的可识别性差，部分参数触及拟合边界。结论是研究支持跨SH阶使用共同D_a，不支持阶数特异D_a，但自由共享D_a与TDE值并不完全数值等同。4.3 噪声与角度采样偏倚的蒙特卡洛评价（Monte Carlo Evaluation of Noise and Angular-Sampling Bias）：通过蒙特卡洛模拟实际64方向采集几何，发现无噪声时有限方向采样误差随SH阶数增加，b≥5000 s/mm²下l=0、2、4、6的中位绝对采样误差约1.2%、2.0%、3.4%、5.9%，最坏方向相关误差达9.7%、16.9%、18.6%、34.3%。Rician幅度噪声在SNR=40时明显改变绝对谐波功率，l=0中位偏倚约+8.4%至+16.4%，l=2约?8.4%至?12.1%，l=4约?5.9%至?10.3%，l=6约?2.2%至+0.5%；但重新拟合振幅后，高b值曲线形状的gRMSE约0.8%–3.2%。幅度噪声还使自由共享D_a上偏约4.7%–9.6%。结论是噪声主要影响谐波功率幅度，而有限角度采样对高阶尤其l=6影响更大。

讨论部分总结：TDE所得D_a（2.26–2.34 μm²/ms）与既往报告范围一致；固定D_a后，l<6的谐波功率在高b值区符合预测，b≥6000 s/mm²时更紧密，l=4改善最明显，而l=6一致性和稳定性较弱。模型比较支持共同D_a，阶数特异D_a缺乏支持且估计受限；自由共享D_a高于TDE值，可能受幅度噪声影响，但研究未要求二者数值完全等同。蒙特卡洛表明噪声诱导的偏倚多为可被振幅吸收的乘性变化，不能解释l=0、2、4的主要一致性；l=6的弱一致性可能与有限角度采样、高阶角结构截断以及其谐波峰值接近采集上限有关。由于谐波功率旋转不变并仅依赖信号SH系数，观察到的渐近b值行为可约束任何基于SH的角结构表示，包括CSD等fODF方法；高阶谐波对fODF多峰分辨和纤维追踪方向估计具有实际意义。研究为技术验证而非群体推断，三名受试者提供独立生物学重复，体素分析提供大量受试者内空间评估但体素非独立生物样本；样本量限制个体间变异推断。Stick模型为理想化描述，结果验证的是观察到的Stick样信号成分的b值依赖，而非唯一确定其细胞来源；残余轴突外信号、有限轴突半径、水交换及其他受限结构可能影响谐波功率，需未来更密角度采样和更敏感采集进一步界定Stick近似的定量适用区间。

研究结论（翻译）：研究人员表明，白质dMRI信号至l=6的偶数阶SH功率在高b值区与Stick模型预测的阶数特异b值依赖一致。通过TDE MRI独立估计并在主要谐波功率分析中固定D_a，仅每阶拟合单一振幅，构成了对模型的严格验证。对于b ? 4000 s/mm²，l=0、2、4显示实验与理论强一致，并在更高b值进一步改善；l=6一致性降低，提示其对噪声和拟合不稳定性更敏感。额外模型比较几乎不支持阶数特异D_a，固定TDE D_a和单一共享D_a表述占 nearly all model support；蒙特卡洛模拟进一步显示幅度噪声可偏倚绝对谐波功率，而拟合的高b值曲线形状相对稳健；有限角度采样效应随谐波阶数增加，在l=6最大。因此，l=6结果支持结论但定量稳健性低于l=0、2、4。总体而言，这些发现进一步支持将Stick模型用于描述大b值下扩散信号的角结构，并为基于SH的白质表示（包括用于纤维追踪和微结构建模的fODF）提供依据。