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

Technische Zusammenfassung: Oberflächenrelaxivität verzerrt fraktionelle Werte der Diffusions- und Relaxometrie-Mikrostruktur

Problemstellung
Oberflächenrelaxivität ($\rho$) und zeitabhängige Diffusion sind zwei beobachtbare Größen, die aus demselben physikalischen Prozess resultieren: dem Kollidieren von Wassermolekülen mit den Wänden myelinisiertem Axone. Während die Diffusions-MRT (dMRI) oft als „relaxivitätsblind“ angenommen wird, da sie durch eine $b=0$-Normalisierung genormt ist, argumentiert dieses Paper, dass diese Annahme nur für eine einzelne Kompartimentierung gilt. In der Realität befinden sich intra-axonale und extra-axonale Wasser gegen unterschiedliche Wände mit unterschiedlichen Oberflächen-zu-Volumen-Verhältnissen ($S/V$). Folglich erfahren sie unterschiedliche transversale Relaxationsraten ($R_2 = 1/T_2 + \rho S/V$). Diese differenzielle $T_2$-Verkürzung gewichtet die Signal-Kompartimente zum Zeitpunkt des Echos ($TE$), der für die Normalisierung verwendet wird, um und führt so zu einer systematischen Verzerrung der Mikrostruktur-Schätzungen ein. Das Paper untersucht, wie dieser unbeobachtete geometrische Term zwei Schlüsselmetriken verzerrt: den intra-axialen Diffusionssignalanteil ($f_{intra}$) und den Myelinwasseranteil (MWF).

Methodik
Die Studie kombiniert analytische Herleitungen mit Monte-Carlo-Simulationen (MC) auf Basis von Grundprinzipien, um die Physik zu validieren und die Verzerrung zu quantifizieren.

1. Analytische Herleitungen:

   • Interne Rate: Ableitung der Brownstein–Tarr Oberflächenrate für intra-axonales Wasser über eine Gamma-verteilte Axonkaliber-Population. Die Rate ist der flächengewichtete Durchschnitt von $\rho(4/d)$.
   • Externe Rate: Ableitung der Novikov–Burcaw Oberflächenrate für extra-axonales Wasser, die im Grenzfall der Bewegungsverengung (motional narrowing) zu einem Plateau sättigt, das durch die mittlere externe Oberflächendichte ($S_{ext}/V_{ext}$) bestimmt wird.
   • Geschlossene Form des Verhältnisses: Nachweis, dass das Verhältnis von internem zu externem $S/V$ unabhängig von der Kaliberverteilung ist und ausschließlich vom Faser-Volumenanteil ($V_F$) und dem $g$-Ratio abhängt: $(1-V_F)/(g V_F)$.
   • Bias-Formeln: Ableitung geschlossener Ausdrücke für die Verzerrung von $f_{intra}$ (Gl. 9) und die Verschiebung des scheinbaren MWF durch das Unterschreiten des Myelin-Fensters durch die $T_2$-Werte dünner Axone.

2. Monte-Carlo-Validierung:

   • Entwicklung eines Wandenzählenden MC-Simulators (unter Verwendung eines realisierten Overshoot-Boundary-Local-Time-Schätzers), um Random Walker durch explizite gepackte Geometrien ohne analytische Relaxivitätsmodelle zu propagieren.
   • Validierung, dass die MC-Oberflächenattenuierung den geschlossenen Brownstein–Tarr- und Novikov–Burcaw-Raten innerhalb von $\sim 0,2$ Prozentpunkten entspricht.
   • Bestätigung, dass dünne Axone in der MC-Simulation $T_2$-Werte aufweisen, die unter das 25 ms Myelin-Fenster fallen, was den Mechanismus für den MWF-Bias validiert.

3. Forward Modeling:

   • Diffusion: Simulation von Multi-Shell dMRI-Signalen mit aktivierter Oberflächenrelaxivität und Anpassung mittels Standard-Spherical-Mean-Techniken (SMT), um die $f_{intra}$-Verzerrung zu quantifizieren.
   • Relaxometrie: Simulation von Multi-Echo-CPMG-Zerfällen mit einer Kaliber-verteilten intra-axialen Population und Schätzung des MWF mittels Standard-NNLS mit einem fixierten Myelin-Fenster.
   • Longitudinaler Verlauf: Modellierung eines longitudinalen Demyelinisierungs-Trajektoriums innerhalb eines Subjekts (Erhalt des Axonlumens, Ausdünnung des Myelins), um zu testen, ob sich der Bias bei longitudinalen Änderungen aufhebt.

Wichtigste Ergebnisse

• Verzerrung des intra-axialen Signalanteils ($f_{intra}$):

  • Die Verzerrung wird durch das differenzielle $S/V$ zwischen den internen und externen Wänden angetrieben.
  • Vorzeichen-Gesetz: Das Vorzeichen der Verzerrung wird durch den Faser-Volumenanteil ($V_F$) bestimmt. Unterhalb eines Crossover-Punktes $V_F^* = 1/(1+g) \approx 0,59$ relaxiert das intra-axonale Wasser schneller (Unterschätzung). Oberhalb von $V_F^*$ relaxiert das extra-axonale Wasser schneller, was dazu führt, dass das intra-axonale Signal überbewertet wird.
  • Größenordnung: In physiologisch dichtem weißem Substanzgewebe ($V_F \approx 0,65$) wird der intra-axonale Anteil bei einer klinischen Echozeit ($TE=80$ ms) und einem zitierten Relaxivitätswert ($\rho = 1,16$ $\mu$m/s) um etwa 12 % überschätzt. Diese Verzerrung skaliert linear mit $\rho$ und wächst monoton mit $TE$.
  • Natur: Es handelt sich um einen primären, räumlich strukturierten systematischen Fehler, kein Nullmittelwert-Rauschen; er gleicht sich also nicht über Voxel hinweg aus, sondern verschiebt regionale und Kohorten-Mittelwerte.

• Myelinwasseranteil (MWF) Bias:

  • Mechanismus: Oberflächenrelaxivität verkürzt die scheinbare $T_2$ des intra-axialen Wassers. Für die dünnsten Axone (innerer Durchmesser $d < d^* \approx 0,17$ $\mu$m bei zitiertem $\rho$) sinkt das $T_2$ unter das Standard-25-ms-Myelin-Fenster. Dieses Wasser wird fälschlicherweise als Myelinwasser klassifiziert.
  • Richtung: Feines weißes Gewebe (hohes $S/V$) erscheint myelinreicher als grobes weißes Gewebe, selbst wenn das wahre Myelinvolumen identisch ist.
  • Größenordnung: Bei dem zitierten $\rho$ ist die Verzerrung gering (0,33 Prozentpunkte im MWF oder ~0,82 pp im inferierten Myelingehalt). Sie ist jedoch superlinear bezüglich $\rho$: Wenn $\rho$ eine Größenordnung höher ist, wird die Verzerrung signifikant (3,4 pp).
  • Detektierbarkeit: Die Verzerrung liegt unter dem Single-Voxel-Rauschen, tritt aber als systematischer Offset auf, wenn man über Regionen von Interesse (ROIs) mit unterschiedlichen Kaliberverteilungen mittelt.

• Longitudinale Demyelinisierung:

  • In einem Lumen-erhaltenden Demyelinisierungs-Trajektorium (primäre Demyelinisierung) wird der MWF-Bias durch den fixen inneren Axondurchmesser bestimmt. Daher ist der Bias nahezu konstant und hebt sich bei der Berechnung der Änderung ($\Delta$MWF) zwischen den Zeitpunkten auf.
  • Im Gegensatz dazu hängt der $f_{intra}$-Bias vom Unterschied zwischen internem und externem $S/V$ ab. Während das Myelin dünner wird, zieht sich die externe Wand zurück, wodurch sich das externe $S/V$ ändert, während das interne fix bleibt. Infolgedessen driftet der $f_{intra}$-Bias signifikant (z. B. von -5 % auf -21 %) während der Demyelinisierung, was ein großes, TE-abhängiges Artefakt in longitudinalen Diffusionsstudien erzeugt.

Bedeutung und Ansprüche
Das Paper behauptet, eine bisher nicht berichtete Quelle der Verzerrung in der Diffusions- und Relaxometrie-Mikrostruktur-Bildgebung identifiziert zu haben, die daraus resultiert, dass Bulk-$T_2$ und geometrische Oberflächenrelaxivität in Multi-Echo-Messungen nicht trennbar sind.

1. Vereinheitlichte Physik: Es stellt fest, dass Oberflächenrelaxivität und zeitabhängige Diffusion duale Messgrößen desselben Wandkollisionsprozesses sind. Das gemessene „$T_2$“ in der Relaxometrie ist durch eine geometrieabhängige Rate kontaminiert, die die dMRI versucht zu normalisieren, dabei aber aufgrund der unterschiedlichen Kompartiment-$S/V$ nicht vollständig eliminieren kann.
2. Systematische Fehler: Die $f_{intra}$-Verzerrung wird als primärer, deterministischer systematischer Fehler charakterisiert, der mit der Packungsdichte und der Echozeit skaliert und potenziell den Vergleich zwischen Standorten oder Kohorten konfundiert, falls unterschiedliche Echozeiten verwendet werden.
3. MWF-Interpretation: Die MWF-Verzerrung wird als struktureller Confounder gerahmt, bei dem „Myelinwasser“ effektiv „hoch-$S/V$-gefangenes Wasser“ misst. Dies impliziert, dass feine weiße Substanzregionen (z. B. Genu des Corpus Callosum) im Vergleich zu groben Regionen (z. B. Mittelkörper) ein höheres Myelingehalt aufweisen können, einfach aufgrund von Kaliberunterschieden, unabhängig vom tatsächlichen Myelinvolumen.
4. Longitudinaler Kontrast: Eine zentrale Erkenntung ist die Divergenz im Verhalten dieser Biase longitudinal. Während MWF-Änderungen in der primären Demyelinisierung gegenüber diesem Bias robust sind (da der Offset sich aufhebt), sind diffusionsbasierte $f_{intra}$-Änderungen stark durch die sich ändernde Geometrie der externen Wand kontaminiert.
5. Limitierungen: Die Autoren merken bescheiden an, dass die absolute Größe der Verzerrung vom schwer zu bestimmenden Wert von $\rho$ abhängt (der um eine Größenordnung variieren kann) und dass die idealisierte Geometrie (nicht-berührende Zylinder) die Effekte in realistischen, sich berührenden Faserkonstruktionen unterschätzen könnte. Sie schlagen kein neues Akquisitionsprotokoll zur Lösung vor, da die Degenerierung strukturell ist und orthogonale Constraints erfordert (z. B. hochgradig gradienten-oszillierende Gradienten), die derzeit klinisch nicht erreichbar sind.