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

Resumen Técnico: El sesgo de la relajación superficial en las fracciones de microestructura de difusión y relaxometría

Planteamiento del problema
La relajación superficial ($\rho$) y la difusión dependiente del tiempo son dos observables que surgen de un mismo proceso físico: las colisiones de las moléculas de agua con las paredes de los axones mielinizados. Aunque a menudo se asume que la RM de difusión (dMRI) es "insensible a la relajación" debido a la normalización con $b=0$, este artículo sostiene que dicha suposición solo es válida para un único compartimento. En realidad, el agua intraaxonal y la extraaxonal residen contra paredes distintas con relaciones superficie-volumen ($S/V$) diferentes. En consecuencia, experimentan tasas de relajación transversal distintas ($R_2 = 1/T_2 + \rho S/V$). Esta diferencia en $T_2$ repondera los compartimentos de la señal en el tiempo de eco ($TE$) utilizado para la normalización, introduciendo un sesgo sistemático en las estimaciones microestructurales. El estudio investiga cómo este término geométrico no observado sesga dos métricas clave: la fracción de señal intraaxonal de difusión ($f_{intra}$) y la fracción de agua de mielina (MWF).

Metodología
El estudio combina derivaciones analíticas con simulaciones de Monte Carlo (MC) de principios fundamentales para validar la física y cuantificar el sesgo.

1. Derivaciones analíticas:

   • Tasa interior: Se derivó la tasa de superficie de Brownstein–Tarr para el agua intraaxonal sobre una población de calibres de axones con distribución Gamma. La tasa es el promedio ponderado por el área de $\rho(4/d)$.
   • Tasa exterior: Se derivó la tasa de superficie de Novikov–Burcaw para el agua extraaxonal, que se satura en una meseta determinada por la densidad de superficie exterior media ($S_{ext}/V_{ext}$) en el límite de estrechamiento de movimiento (motional-narrowing).
   • Relación de forma cerrada: Se demostró que la relación entre el $S/V$ interior y el exterior es independiente de la distribución de calibres y depende únicamente de la fracción de volumen de la fibra ($V_F$) y el $g$-ratio: $(1-V_F)/(g V_F)$.
   • Fórmulas de sesgo: Se derivaron expresiones de forma cerrada para el sesgo en $f_{intra}$ (Ec. 9) y el desplazamiento en la MWF aparente basado en el cruce de los valores de $T_2$ de axones finos por debajo de la ventana de la mielina.

2. Validación de Monte Carlo:

   • Se desarrolló un simulador de conteo de paredes de MC (utilizando un estimador de tiempo local de frontera de exceso realizado) para propagar caminantes aleatorios a través de geometrías explícitas empaquetadas sin modelos de relajación analíticos.
   • Se validó que la atenuación de superficie de MC coincide con las tasas de Brownstein–Tarr y Novikov–Burcaw con una precisión de $\sim0.2$ puntos porcentuales.
   • Se confirmó que los axones finos en la simulación de MC exhiben valores de $T_2$ que caen por debajo de la ventana de mielina de 25 ms, validando el mecanismo del sesgo de la MWF.

3. Modelado hacia adelante (Forward Modeling):

   • Difusión: Se simularon señales de dMRI de múltiples capas con relajación superficial habilitada y se ajustaron utilizando técnicas estándar de media esférica (SMT) para cuantificar el sesgo de $f_{intra}$.
   • Relaxometría: Se simularon decaimientos CPMG de múltiples ecos con un grupo intraaxonal de calibre distribuido y se estimó la MWF utilizando mínimos cuadrados no negativos (NNLS) con una ventana de mielina fija.
   • Trayectoria longitudinal: Se modeló una trayectoria de desmielinización intra-sujeto (preservando el lumen del axón, adelgazando la mielina) para probar si el sesgo se cancela en los cambios longitudinales.

Resultados clave

• Sesgo de la fracción de señal intraaxonal ($f_{intra}$):

  • El sesgo es impulsado por la diferencia de $S/V$ entre las paredes interior y exterior.
  • Ley de signo: El signo del sesgo está determinado por la fracción de volumen de la fibra ($V_F$). Por debajo de un punto de cruce $V_F^* = 1/(1+g) \approx 0.59$, el agua intraaxonal se relaja más rápido (subestimación). Por encima de $V_F^*$, el agua extraaxonal se relaja más rápido, causando que la señal intraaxonal se sobreponderre.
  • Magnitud: En la sustancia blanca fisiológicamente densa ($V_F \approx 0.65$), la fracción intraaxonal se sobreestima aproximadamente en un 12% en un tiempo de eco clínico ($TE=80$ ms) utilizando un valor de relajación citado ($\rho = 1.16$ $\mu$m/s). Este sesgo escala linealmente con $\rho$ y crece monótonamente con $TE$.
  • Naturaleza: Este es un error sistemático de primer orden y espacialmente estructurado, no un ruido de media cero, lo que significa que no se promedia entre vóxeles, sino que desplaza las medias regionales y de cohorte.

• Sesgo de la Fracción de Agua de Mielina (MWF):

  • Mecanismo: La relajación superficial acorta el $T_2$ aparente del agua intraaxonal. Para los axones más finos (diámetro interno $d < d^* \approx 0.17$ $\mu$m con el $\rho$ citado), el $T_2$ cae por debajo de la ventana estándar de mielina de 25 ms. Esta agua es clasificada erróneamente como agua de mielina.
  • Dirección: La sustancia blanca fina (alta $S/V$) se lee como más rica en mielina que la sustancia blanca gruesa, incluso si el volumen de mielina real es idéntico.
  • Magnitud: Con el $\rho$ citado, el sesgo es pequeño (0.33 puntos porcentuales en MWF, o ~0.82 pp en contenido de mielina inferido). Sin embargo, es superlineal respecto a $\rho$; si $\rho$ es un orden de magnitud mayor, el sesgo se vuelve significativo (3.4 pp).
  • Detectabilidad: El sesgo está oculto bajo el ruido de un solo vóxel, pero emerge como un desfase sistemático al promediar regiones de interés (ROI) con diferentes distribuciones de calibre.

• Diferenciación Longitudinal:

  • En una trayectoria de desmielinización de lumen preservado (desmielinización primaria), el sesgo de la MWF está determinado por el diámetro interno fijo del axón. Por lo tanto, el sesgo es casi constante y se cancela al calcular el cambio ($\Delta$MWF) entre puntos temporales.
  • Por el contrario, el sesgo de $f_{intra}$ depende de la diferencia entre el $S/V$ interior y el exterior. A medida que la mielina se adelgaza, la pared exterior se retira, alterando el $S/V$ exterior mientras que el interior permanece fijo. En consecuencia, el sesgo de $f_{intra}$ deriva significativamente (por ejemplo, de -5% a -21%) durante la desmielinización, creando un artefacto de gran magnitud dependiente del $TE$ en estudios de difusión longitudinales.

Significado y Reivindicaciones
El artículo afirma identificar una fuente de sesgo no reportada anteriormente en la imagen de microestructura tanto de difusión como de relaxometría, que surge de la inseparabilidad del $T_2$ de volumen y la relajación superficial geométrica en mediciones de múltiples ecos.

1. Física Unificada: Establece que la relajación superficial y la difusión dependiente del tiempo son lecturas duales del mismo proceso de colisión con la pared. El "$T_2$" medido en la relaxometría está contaminado por una tasa dependiente de la geometría que la RM de difusión intenta normalizar pero no logra eliminar por completo debido a las diferencias de $S/V$ entre compartimentos.
2. Errores Sistemáticos: El sesgo de $f_{intra}$ se caracteriza como un error sistemático de primer orden y determinista que escala con la densidad de empaquetamiento y el tiempo de eco, pudiendo confundir las comparaciones entre centros o cohortes si los tiempos de eco difieren.
3. Interpretación de la MWF: El sesgo de la MWF se plantea como un error estructural donde el "agua de mielina" mide efectivamente "agua atrapada de alto $S/V$". Esto implica que las regiones de sustancia blanca fina (ej. el genu del cuerpo calloso) pueden parecer tener un mayor contenido de mielina que las regiones de sustancia blanca gruesa (ej. el cuerpo medio) simplemente debido a las diferencias de calibre, independientemente del volumen de mielina real.
4. Contraste Longitudinal: Un hallazgo clave es la divergencia en cómo se comportan estos sesgos longitudinalmente. Mientras que los cambios en la MWF son robustos a este sesgo en la desmielinización primaria (ya que el desfase se cancela), los cambios en la $f_{intra}$ basados en difusión están fuertemente contaminados por el cambio en la geometría de la pared exterior.
5. Limitaciones: Los autores señalan modestamente que la magnitud absoluta del sesgo depende del valor de $\rho$ (que puede variar en un orden de magnitud) y que la geometría idealizada (cilindros que no se tocan) puede subestimar los efectos en empaquetamientos de fibras realistas que sí se tocan. No proponen un nuevo protocolo de adquisición para resolver esto, señalando que la degeneración es estructural y requiere restricciones ortogonales (ej. gradientes oscilantes de alto gradiente) que actualmente están fuera del alcance clínico.