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

Sintesi Tecnica: La rilassività superficiale distorce le frazioni microstrutturali di diffusione e rilassometria

Definizione del Problema
La rilassività superficiale ($\rho$) e la diffusione tempo-dipendente sono due osservabili derivanti dallo stesso processo fisico: la collisione delle molecole d'acqua con le pareti degli assoni mielinizzati. Sebbene la dMRI sia spesso assunta come "insensibile alla rilassività" a causa della normalizzazione $b=0$, questo articolo sostiene che tale assunzione sia valida solo per un singolo compartimento. Nella realtà, l'acqua intra-assonale ed extra-assonale risiede contro pareti diverse con diversi rapporti superficie/volume ($S/V$). Di conseguenza, esse sperimentano tassi di rilassamento trasversale differenti ($R_2 = 1/T_2 + \rho S/V$). Questo $T_2$ differenziale ri-pesa i segnali dei compartimenti al tempo di eco ($TE$) utilizzato per la normalizzazione, introducendo un bias sistematico nelle stime microstrutturali. Il documento investiga come questo termine geometrico non osservato distorca due metriche chiave: la frazione di segnale intra-assonale di diffusione ($f_{intra}$) e la frazione di acqua mielinica (MWF).

Metodologia
Lo studio combina derivazioni analitiche con simulazioni Monte Carlo (MC) basate sui primi principi per validare la fisica e quantificare il bias.

1. Derivazioni Analitiche:

   • Tasso Interiore: Derivato il tasso superficiale di Brownstein–Tarr per l'acqua intra-assonale su una popolazione di calibro assonale distribuita secondo una funzione Gamma. Il tasso è la media ponderata sull'area di $\rho(4/d)$.
   • Tasso Esterno: Derivato il tasso superficiale di Novikov–Burcaw per l'acqua extra-assonale, che satura verso un plateau determinato dalla densità superficiale esterna media ($S_{ext}/V_{ext}$) nel limite del restringimento motorio (motional-narrowing).
   • Rapporto in Forma Chiusa: Dimostrato che il rapporto tra $S/V$ interiore ed esterno è indipendente dalla distribuzione del calibro e dipende esclusivamente dalla frazione di volume delle fibre ($V_F$) e dal rapporto $g$: $(1-V_F)/(g V_F)$.
   • Formule del Bias: Derivate espressioni in forma chiusa per il bias in $f_{intra}$ (Eq. 9) e lo spostamento della MWF apparente in base all'incrocio dei valori di $T_2$ degli assoni sottili al di sotto della finestra della mielina.

2. Validazione Monte Carlo:

   • Sviluppato un simulatore MC di conteggio delle pareti (utilizzando un estimatore del tempo locale del bordo a realizzazione-overshoot) per propagare i random walker attraverso geometrie esplicite impaccate senza modelli di rilassività analitici.
   • Validato che l'attenuazione superficiale MC corrisponde ai tassi Brownstein–Tarr e Novikov–Burcaw in forma chiusa con un errore di circa $\sim 0,2$ punti percentuali.
   • Confermato che gli assoni sottili nella simulazione MC esibiscono valori di $T_2$ che scendono al di sotto della finestra della mielina di 25 ms, validando il meccanismo per il bias della MWF.

3. Modellazione Forward:

   • Diffusione: Simulati segnali dMRI multi-shell con rilassività superficiale abilitata e adattati utilizzando tecniche standard di media sferica (SMT) per quantificare il bias di $f_{intra}$.
   • Rilassometria: Simulati decadimenti CPMG multi-eco con un pool intra-assonale a calibro distribuito e stimata la MWF utilizzando l'algoritmo NNLS standard con una finestra di mielina fissa.
   • Traiettoria Longitudinale: Modellata una traiettoria di demielinizzazione intra-soggetto (preservando il lume dell'assone, assottigliando la mielina) per testare se il bias si annulli nei cambiamenti longitudinali.

Risultati Chiave

• Bias della Frazione di Segnale Intra-assonale ($f_{intra}$):

  • Il bias è guidato dal diverso $S/V$ tra le pareti interne ed esterne.
  • Legge del Segno: Il segno del bias è determinato dalla frazione di volume delle fibre ($V_F$). Al di sotto di un punto di incrocio $V_F^* = 1/(1+g) \approx 0,59$, l'acqua intra-assonale rilassa più velocemente (sottostima). Al di sopra di $V_F^*$, l'acqua extra-assonale rilassa più velocemente, causando una sovrastima del segnale intra-assonale.
  • Magnitudo: In una sostanza bianca fisiologicamente densa ($V_F \approx 0,65$), la frazione intra-assonale è sovrastimata di circa il 12% a un tempo di eco clinico ($TE=80$ ms) utilizzando un valore di rilassività citato ($\rho = 1,16$ $\mu$m/s). Questo bias scala linearmente con $\rho$ e cresce monotonicamente con $TE$.
  • Natura: Si tratta di un errore sistematico di primo ordine e spazialmente strutturato, non di rumore a media nulla; pertanto, non si annulla mediando tra i voxel, ma sposta le medie regionali e di coorte.

• Bias della Frazione di Acqua Mielinica (MWF):

  • Meccanismo: La rilassività superficiale accorcia l'apparente $T_2$ dell'acqua intra-assonale. Per gli assoni più sottili (diametro interno $d < d^* \approx 0,17$ $\mu$m al $\rho$ citato), il $T_2$ scende al di sotto della finestra standard della mielina di 25 ms. Quest'acqua viene classificata erroneamente come acqua mielinica.
  • Direzione: La sostanza bianca fine (alto $S/V$) risulta più ricca di mielina rispetto alla sostanza bianca grossolana, anche se il volume reale della mielina è identico.
  • Magnitudo: Al $\rho$ citato, il bias è piccolo (0,33 punti percentuali nella MWF, o ~0,82 pp nel contenuto di mielina inferito). Tuttavia, è super-lineare rispetto a $\rho$; se $\rho$ fosse un ordine di grandezza superiore, il bias diventerebbe significativo (3,4 pp).
  • Rilevabilità: Il bias è sepolto sotto il rumore del singolo voxel, ma emerge come un offset sistematico quando si mediano regioni di interesse (ROI) con diverse distribuzioni di calibro.

• Demielinizzazione Longitudinale:

  • In una traiettoria di demielinizzazione con lume preservato (demielinizzazione primaria), il bias della MWF è determinato dal diametro assonale interno fisso. Pertanto, il bias è quasi costante e si annulla quando si calcola la variazione ($\Delta$MWF) tra i punti temporali.
  • Al contrario, il bias di $f_{intra}$ dipende dalla differenza tra $S/V$ interno ed esterno. Mentre la mielina si assottiglia, la parete esterna recede, alterando il $S/V$ esterno mentre l'interno rimane fisso. Di conseguenza, il bias di $f_{intra}$ deriva significativamente (ad esempio, da -5% a -21%) durante la demielinizzazione, creando un grande artefatto dipendente dal $TE$ negli studi longitudinali sulla diffusione.

Significato e Rivendicazioni
L'articolo rivendica di aver identificato una fonte di bias precedentemente non riportata in entrambe l'imaging della microstruttura di diffusione e rilassometria, che deriva dall'inseparabilità tra $T_2$ bulk e rilassività superficiale geometrica nelle misurazioni multi-eco.

1. Fisica Unificata: Stabilisce che la rilassività superficiale e la diffusione tempo-dipendente sono letture duali dello stesso processo di collisione con la parete. Il "$T_2$" misurato nella rilessometria è contaminato da un tasso dipendente dalla geometria che la dMRI tenta di normalizzare ma non riesce a eliminare completamente a causa delle differenze di $S/V$ tra i compartimenti.
2. Errori Sistematici: Il bias di $f_{intra}$ è caratterizzato come un errore sistematico deterministico di primo ordine che scala con la densità di impacchettamento e il tempo di eco, potenzialmente confondendo i confronti tra siti o tra coorti se i tempi di eco differiscono.
3. Interpretazione della MWF: Il bias della MWF è inquadrato come un confondimento strutturale dove l' "acqua mielinica" misura effettivamente l' "acqua intrappolata ad alto $S/V$". Ciò implica che le regioni di sostanza bianca fine (es. genu del corpo calloso) possono apparire più ricche di mielina rispetto alle regioni di sostanza bianca grossolana (es. corpo centrale) semplicemente a causa delle differenze di calibro, indipendentemente dal volume reale della mielina.
4. Contrasto Longitudinale: Una scoperta chiave è la divergenza nel modo in cui questi bias si comportano longitudinalmente. Mentre i cambiamenti della MWF nella demielinizzazione primaria sono robusti a questo bias (poiché l'offset si annulla), i cambiamenti della $f_{intra}$ basati sulla diffusione sono pesantemente contaminati dalla variazione della geometria della parete esterna.
5. Limitazioni: Gli autori notano con modestia che l'entità assoluta del bias dipende dal valore poco stabilito di $\rho$ (che può variare di un ordine di grandezza) e che la geometria idealizzata (cilindri non a contatto) può sottostimare gli effetti in pacchetti di fibre realistici che si toccano. Non propongono un nuovo protocollo di acquisizione per risolvere il problema, osservando che la degenerazione è strutturale e richiede vincoli ortogonali (es. gradienti oscillanti ad alta intensità) che sono attualmente oltre la portata clinica.