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

技術要約：表面緩和率が拡散およびリラクソメトリの微細構造分率に及ぼすバイアス

問題提起
表面緩和率（$\rho$）と時間依存的な拡散は、同一の物理的プロセス、すなわち水分子が髄鞘化された軸索の壁に衝突するという現象から生じる2つの観測量である。拡散MRI（dMRI）は、$b=0$による正規化により「緩和率に対して盲目（relaxivity-blind）」であるとしばしば仮定されるが、本論文は、この仮定は単一のコンパートメントに対してのみ成立すると論じている。実際には、軸索内（intra-axonal）の水と軸索外（extra-axonal）の水は、異なる表面積対体積比（$S/V$）を持つ異なる壁に接している。その結果、両者は異なる横緩和率（$R_2 = 1/T_2 + \rho S/V$）を経験する。この差のある$T_2$は、正規化に使用されるエコータイム（$TE$）における信号コンパートメントの重み付けを変化させ、微細構造推定に系統的なバイアスを導入する。本論文では、この観測不可能な幾何学的項が、2つの主要な指標である拡散軸索内信号分率（$f_{intra}$）およびミエリン水分率（MWF）にどのようにバイアスを与えるかを調査している。

手法
本研究は、物理現象の検証とバイアスの定量化のために、解析的な導出と第一原理に基づくモンテカルロ（MC）シミュレーションを組み合わせている。

1. 解析的導出:

   • 内部レート: ガンマ分布に従う軸索径集団に対する、ブラウンシュタイン・タール（Brownstein–Tarr）の表面レートを導出した。このレートは、$\rho(4/d)$の面積加重平均である。
   • 外部レート: 軸索外の水に対するノビコフ・バーカウ（Novikov–Burcaw）の表面レートを導出した。これは、運動型狭窄限界における平均外部表面密度（$S_{ext}/V_{ext}$）によって決定されるプラトーへと飽和する。
   • 閉形式の比: 内部と外部の$S/V$の比は、径分布に依存せず、線維体積分率（$V_F$）と$g$比のみに依存することを示した：$(1-V_F)/(g V_F)$。
   • バイアス公式: $f_{intra}$におけるバイアス（式9）および、細い軸索の$T_2$値がミエリン窓を下回ることによる見かけのMWFのシフトに関する閉形式の式を導出した。

2. モンテカルロによる検証:

   • 解析的な緩和率モデルを用いずに、明示的な充填幾何学構造内をランダムウォーカーを伝播させるための、壁カウント型MCシミュレータ（実現されたオーバーシュート境界局所時間推定器を使用）を開発した。
   • MCによる表面減衰が、解析的なブラウンシュタイン・タールおよびノビコフ・バーカウのレートと$\sim0.2$パーセントポイント以内で一致することを確認し、MCシミュレータを検証した。
   • 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$が1桁大きければ、バイアスは顕著になる（約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$の値（これは1桁程度変動する可能性がある）に依存すること、および理想的な幾何学（接触しない円柱）が、現実の接触している線維パックにおける影響を過小評価している可能性があることを控えめに述べている。彼らは、この縮退を解決するための新しい撮像プロトコルを提案しておらず、それは現在の臨床技術を超えた、高勾配振動勾配などの直交する制約を必要とするものであるとしている。