Two observables of one wall: how surface relaxivity can bias the diffusion intra-axonal fraction and the myelin water fraction

技术摘要：表面弛豫率偏置扩散与弛豫测量中的微结构比例

问题陈述
表面弛豫率 ($\rho$) 和随时间变化的扩散是源于同一物理过程的两个可观测变量：水分子与有髓轴突壁的碰撞。虽然扩散磁共振成像 (dMRI) 通常由于 $b=0$ 归一化而被认为对“弛豫率不敏感”，但本文认为，这一假设仅在单室模型中成立。实际上，轴内和轴外水分子位于具有不同表面积/体积比 ($S/V$) 的不同壁面旁。因此，它们经历不同的横向弛豫率 ($R_2 = 1/T_2 + \rho S/V$)。这种差异化的 $T_2$ 在用于归一化的回波时间 ($TE$) 处对信号组分进行了重新加权，从而在微结构估计中引入了系统性偏差。本文研究了这一未观察到的几何项如何偏置两个关键指标：扩散轴内信号比例 ($f_{intra}$) 和髓鞘水比例 (MWF)。

方法论
本研究结合了解析推导与基于第一原理的蒙特卡洛 (MC) 模拟，以验证物理机制并量化偏差。

1. 解析推导：

   • 内部速率： 推导了针对 Gamma 分布轴径群的 Brownstein–Tarr 轴内水表面速率。该速率是 $\rho(4/d)$ 的面积加权平均值。
   • 外部速率： 推导了针对轴外水的 Novikov–Burcaw 表面速率，该速率在运动窄化极限下饱和于由平均外部表面密度 ($S_{ext}/V_{ext}$) 决定的平台值。
   • 闭式比率： 证明了内外 $S/V$ 的比率独立于直径分布，仅取决于纤维体积分数 ($V_F$) 和 $g$-ratio，即 $(1-V_F)/(g V_F)$。
   • 偏差公式： 推导了关于 $f_{intra}$ 偏差 (Eq. 9) 以及基于细轴 $T_2$ 值跌落至髓鞘窗口以下的 Apparent MWF 偏移的闭式表达式。

2. 蒙特卡洛验证：

   • 开发了一个壁计数蒙特卡洛模拟器（使用实现的过冲边界局部时间估计量），通过显式堆积几何结构传播随机游走者，而不使用解析弛豫率模型。
   • 验证了 MC 表面衰减与闭式 Brownstein–Tarr 及 Novikov–Burcaw 速率的吻合度在 $\sim0.2$ 个百分点以内。
   • 在 MC 模拟中确认了细轴表现出跌落至 25 ms 髓鞘窗口以下的 $T_2$ 值，从而验证了 MWF 偏差的机制。

3. 前向建模：

   • 扩散： 模拟了启用表面弛豫率的多壳层 dMRI 信号，并使用标准的球面平均技术 (SMT) 进行拟合，以量化 $f_{intra}$ 偏差。
   • 弛豫测量： 模拟了具有直径分布的轴内水池的多回波 CPMG 衰减，并使用固定的髓鞘窗口通过标准非负最小二乘法 (NNLS) 估计 MWF。
   • 纵向轨迹： 模拟了一个在患者内的脱髓鞘轨迹（保持轴管容积，仅髓鞘变薄），以测试该偏差是否在纵向变化中抵消。

关键结果

• 轴内信号比例 ($f_{intra}$) 偏差：

  • 偏差由内外壁之间的差异化 $S/V$ 驱动。
  • 符号规律： 偏差符号由纤维体积分数 ($V_F$) 决定。在交叉点 $V_F^* = 1/(1+g) \approx 0.59$ 以下，轴内水弛豫更快（低估）。在 $V_F^*$ 以上，轴外水弛豫更快，导致轴内信号被过度加权。
  • 量级： 在生理高密度白质中（$V_F \approx 0.65$），在使用引用的弛豫率值 ($\rho = 1.16$ $\mu$m/s) 时，临床回波时间 ($TE=80$ ms) 下的轴内比例会被高估约 12%。该偏差随 $\rho$ 线性缩放，并随 $TE$ 单调增长。
  • 性质： 这是一个一阶的、空间结构化的系统误差，而非零均值噪声，这意味着它不会在不同体素间抵消，而是会移动区域和队列的均值。

• 髓鞘水比例 (MWF) 偏差：

  • 机制： 表面弛豫率缩短了轴内水的表观 $T_2$。对于最细的轴（内径 $d < d^* \approx 0.17$ $\mu$m，基于引用 $\rho$），其 $T_2$ 会跌落至标准的 25 ms 髓鞘窗口以下。这些水被误分类为髓鞘水。
  • 方向： 细白质（高 $S/V$）读出的髓鞘含量比粗白质更丰富，即使两者的真实髓鞘体积相同。
  • 量级： 在引用的 $\rho$ 下，偏差较小（MWF 中约为 0.33 个百分点，或推断的髓鞘含量约为 0.82 pp）。然而，该偏差相对于 $\rho$ 是超线性的；如果 $\rho$ 高出一个数量级，偏差将变得显著（约 3.4 pp）。
  • 可检测性： 该偏差被掩盖在单体素噪声之下，但在对具有不同直径分布的感兴趣区域 (ROI) 进行平均时，会表现为系统性偏移。

• 纵向脱髓鞘：

  • 在保持轴管容积的脱髓辖轨迹（原发性脱髓鞘）中，MWF 偏差由固定的轴内直径决定。因此，该偏差几乎是恒定的，在计算两个时间点之间的变化 ($\Delta$MWF) 时会相互抵消。
  • 相反，$f_{intra}$ 偏差取决于内外 $S/V$ 的差异。随着髓鞘变薄，外部壁退缩，改变了外部 $S/V$ 而内部保持不变。因此，$f_{intra}$ 偏差在脱髓鞘过程中会发生显著漂移（例如，从 -5% 变为 -21%），在纵向扩散研究中产生巨大的、依赖于 $TE$ 的伪影。

意义与主张
本文声称识别了一种此前未被报道的、源于多回波测量中体相 $T_2$ 与几何表面弛豫率不可分割性的扩散与弛豫微结构成像偏差来源。

1. 统一物理学： 它确立了表面弛豫率和随时间变化的扩散是同一壁碰撞过程的双重读数。弛豫测量中测得的“$T_2$”受到了依赖于几何形状的速率的污染，而扩散 MRI 试图通过归一化消除此影响，但由于各组分间的 $S/V$ 差异，无法完全消除。
2. 系统误差： $f_{intra}$ 偏差被描述为一种一阶的、确定性的系统误差，随填充密度和回波时间缩放，如果回波时间不同，可能会混淆跨中心或跨队列的比较。
3. MWF 解释： MWF 偏差被界定为一个结构性混淆，即“髓鞘水”实际上测量的是“高 $S/V$ 捕获水”。这意味着细白质区域（如胼径前部）可能比粗白质区域（如中体）表现出更高的髓鞘含量，这独立于实际的髓鞘体积。
4. 纵向对比： 一个关键发现是这些偏差在纵向行为上的分歧。在原发性脱髓鞘中，MWF 的变化对该偏差具有鲁棒性（因为偏移量会抵消）；而基于扩散的 $f_{intra}$ 变化则受到移动外部壁几何形状的严重污染。
5. 局限性： 作者谦虚地指出，偏差的绝对量级取决于尚不明确的 $\rho$ 值（其可能存在数量级的差异），且理想化几何结构（非接触圆柱体）可能会低估在现实接触纤维包中的效应。他们并未提出新的采集方案来解决此问题，而是指出这种简并性是结构性的，需要目前临床难以实现的正交约束（如高梯度振荡梯度）。