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

Technische Samenvatting: Oppervlakte-relaxiviteit vertekent de fracties van diffusie- en relaxometrie-microstructuur

Probleemstelling
Oppervlakte-relaxiviteit ($\rho$) en tijdsafhankelijke diffusie zijn twee observeerbare grootheden die voortkomen uit hetzelfde fysieke proces: de botsingen van watermoleculen met de wanden van gemyeliniseerde axonen. Hoewel diffusie-MRI (dMRI) vaak wordt verondersteld "relaxiviteit-blind" te zijn vanwege de $b=0$ normalisatie, betoogt dit artikel dat deze aanname alleen geldt voor een enkel compartiment. In werkelijkheid bevinden intra-axonaal en extra-axonaal water zich tegen verschillende wanden met verschillende oppervlakte-volumeverhoudingen ($S/V$). Daarom ervaren zij verschillende transversale relaxatiesnelheden ($R_2 = 1/T_2 + \rho S/V$). Deze differentiële $T_2$ herweegt de signaalcompartimenten bij de echo-tijd ($TE$) die voor normalisatie wordt gebruikt, wat een systematische bias introduceert in microstructuur-schattingen. Het artikel onderzoekt hoe deze onwaargenomen geometrische term twee belangrijke metrieken beïnvloedt: de diffusie intra-axonale signaalfractie ($f_{intra}$) en de myeline-waterfractie (MWF).

Methodologie
De studie combineert analytische afleidingen met Monte Carlo (MC) simulaties op basis van eerste principes om de fysica te valideren en de bias te kwantificeren.

1. Analytische Afleidingen:

   • Interne Snelheid: De Brownstein–Tarr oppervlaktesnelheid voor intra-axonaal water over een Gamma-verdeelde axonale diameterpopulatie afgeleid. De snelheid is het oppervlakte-gewogen gemiddelde van $\rho(4/d)$.
   • Externe Snelheid: De Novikov–Burcaw oppervlaktesnelheid voor extra-axonaal water afgeleid, die verzadigt naar een plateau bepaald door de gemiddelde externe oppervlaktedichtheid ($S_{ext}/V_{ext}$) in de motionele vernauwingslimiet.
   • Gesloten Vorm Ratio: Gedemonstreerd dat de ratio van de interne naar de externe $S/V$ onafhankelijk is van de diameterverdeling en uitsluitend afhangt van de vezelvolumefractie ($V_F$) en de $g$-ratio: $(1-V_F)/(g V_F)$.
   • Bias Formules: Afgeleide gesloten vorm expressies voor de bias in $f_{intra}$ (Vergelijking 9) en de verschuiving in de schijnbare MWF gebaseerd op het oversteken van dunne-axon $T_2$-waarden onder het myeline-venster.

2. Monte Carlo Validatie:

   • Een wall-counting MC-simulator ontwikkeld (gebruikmakend van een realized-overshoot boundary local-time estimator) om random walkers te propageren door expliciete gepakte geometrieën zonder analytische relaxiviteitsmodellen.
   • Gevalideerd dat de MC-oppervlakte-attenuatie overeenkomt met de gesloten vorm Brownstein–Tarr en Novikov–Burcaw snelheden tot op $\sim0,2$ procentpunten.
   • Bevestigd dat dunne axonen in de MC-simulatie $T_2$-waarden vertonen die onder de 25 ms myeline-venster kruisen, wat de mechanisme voor MWF-bias valideert.

3. Forward Modeling:

   • Diffusie: Multi-shell dMRI-signalen gesimuleerd met ingeschakelde oppervlakte-relaxiviteit en deze gefit met standaard sferische-gemiddelde technieken (SMT) om de $f_{intra}$-bias te kwantificeren.
   • Relaxometrie: Multi-echo CPMG-vervallen gesimuleerd met een diameter-verdeelde intra-axonale pool en MWF geschat met standaard non-negative least squares (NNLS) met een vast myeline-venster.
   • Longitudinale Trajectorie: Een binnen-subject demyelinisatie-traject gemodelleerd (behoud van axon lumen, dunner wordend myeline) om te testen of de bias annuleert bij longitudinale veranderingen.

Belangrijkste Resultaten

• Intra-axonale Signaalfractie ($f_{intra}$) Bias:

  • De bias wordt gedreven door de differentiële $S/V$ tussen de interne en externe wanden.
  • Tekenwet: Het teken van de bias wordt bepaald door de vezelvolumefractie ($V_F$). Onder een crossover punt $V_F^* = 1/(1+g) \approx 0,59$, relaxeert intra-axonaal water sneller (onderschatting). Boven $V_F^*$ relaxeert extra-axonaal water sneller, wat leidt tot een overschatting van het intra-axonale signaal.
  • Grootte: In fysiologisch dicht wit weefsel ($V_F \approx 0,65$), wordt de intra-axonale fractie met ongeveer 12% overschat bij een klinische echo-tijd ($TE=80$ ms) gebruikmakend van een geciteerde relaxiviteitswaarde ($\rho = 1,16$ $\mu$m/s). Deze bias schaalt lineair met $\rho$ en groeit monotoon met $TE$.
  • Aard: Dit is een eerste-orde, ruimtelijk gestructureerde systematische fout, geen zero-mean ruis, wat betekent dat het niet gemiddeld wordt over voxels maar regionale en cohort-gemiddelden verschuift.

• Myeline-waterfractie (MWF) Bias:

  • Mechanisme: Oppervlakte-relaxiviteit verkort de schijnbare $T_2$ van intra-axonaal water. Voor de dunste axonen (binnenste diameter $d < d^* \approx 0,17$ $\mu$m bij geciteerde $\rho$) daalt de $T_2$ onder het standaard 25 ms myeline-venster. Dit water wordt verkeerd geclassificeerd als myeline-water.
  • Richting: Fijn wit weefsel (hoge $S/V$) leest als myeline-rijker dan grof wit weefsel, zelfs als de werkelijke myline-volume identiek is.
  • Grootte: Bij de geciteerde $\rho$ is de bias klein (0,33 procentpunt in MWF, of ~0,82 pp in geïnferred myline-gehalte). Echter, de bias is super-lineair ten opzichte van $\rho$; als $\rho$ een orde van grootte hoger is, wordt de bias significant (3,4 pp).
  • Detecteerbaarheid: De bias is begraven onder single-voxel ruis, maar komt naar voren als een systematische offset bij het middelen over regio's van interesse (ROI's) met verschillende diameterverdelingen.

• Longitudinale Demyelinisatie:

  • In een lumen-preserverend demyelinisatie-traject (primaire demyelinisatie) wordt de MWF-bias bepaat door de vaste binnenste axon diameter. Hierdoor is de bias bijna constant en annuleert deze bij het berekenen van de verandering ($\Delta$MWF) tussen tijdstippen.
  • Daarentegen hangt de $f_{intra}$ bias af van het verschil tussen de interne en externe $S/V$. Terwijl het myeline dunner wordt, trekt de externe wand zich terug, waardoor de externe $S/V$ verandert terwijl de interne gelijk blijft. Bijgevolg drift de $f_{intra}$ bias aanzienlijk (bijv. van -5% naar -21%) tijdens demyelinisatie, wat een grote, TE-afhankelijke artefact creëert in longitudinale diffusie-studies.

Betekenis en Claims
Het artikel claimt een voorheen ongerapporteerde bron van bias geïdentificeerd te hebben in zowel diffusie- als relaxometrie-microstructuur imaging, voortvloeiend uit de onscheidbaarheid van bulk $T_2$ en geometrische oppervlakte-relaxiviteit in multi-echo metingen.

1. Verenigde Fysica: Het stelt dat oppervlakte-relaxiviteit en tijdsafhankelijke diffusie duale uitdragen zijn van hetzelfde wand-botsingsproces. De "T2" gemeten in relaxometrie is gecontamineerd door een geometrie-afhankelijke snelheid die diffusie-MRI probeert weg te normaliseren, maar hierin faalt door de verschillen in compartimentele $S/V$.
2. Systematische Fouten: De $f_{intra}$ bias wordt gekarakteriseerd als een eerste-orde, deterministische systematische fout die schaalt met verpakkingsdichtheid en echo-tijd, wat potentieel cohort-vergelijkingen bemoeilijkt als echo-tijden verschillen.
3. MWF Interpretatie: De MWF-bias wordt geframed als een structurele confound waarbij "myeline-water" effectief "hoog-$S/V$ gevangen water" meet. Dit impliceert dat regio's met fijn wit weefsel (bijv. genu van het corpus callosum) meer myline-rijk kunnen lijken dan grove regio's (bijv. mid-body), simpelweg door diameterverschillen, onafhankelijk van het werkelijke myline-volume.
4. Longitudinaal Contrast: Een sleutelfinding is de divergentie in hoe deze biases longitudinaal gedragen. Terwijl MWF-veranderingen in primaire demyelinisatie robuust zijn voor deze bias (omdat de offset annuleert), zijn diffusie-gebaseerde $f_{intra}$ veranderingen zwaar gecontamineerd door de verschuivende externe wand-geometrie.
5. Beperkingen: De auteurs merken bescheiden op dat de absolute grootte van de bias afhangt van de slecht vastgestelde waarde van $\rho$ (die met een orde van grootte kan variëren) en dat de geïdealiseerde geometrie (niet-raakende cilinders) de effecten in realistische, raakende vezelpacks mogelijk onderschat. Ze stellen geen nieuw acquisitieprotocol voor om dit op te lossen, aangezien de degeneratie structureel is en orthogonale beperkingen vereist die momenteel buiten klinisch bereik liggen.