Characterization of spinal cord white matter by suppressing signal from hindered space. A Monte Carlo simulation and an ex vivo ultrahigh-b diffusion-weighted imaging study

Ultrahigh-b diffusion-weighted imaging for quantitative evaluation of myelination in shiverer mouse spinal cord

PURPOSE

To perform a quantitative evaluation of myelination on WT and myelin-deficient (shiverer) mouse spinal cords using ultrahigh-b diffusion-weighted imaging (UHb-DWI).

METHODS

UHb-DWI of ex vivo on spinal cord specimens of two shiverer (C3HeB/FeJ-shiverer, homozygous genotype for MbP_shi ) and six WT (Black Six, C3HeB/FeJ) mice were acquired using 3D multishot diffusion-weighted stimulated-echo EPI, a homemade RF coil, and a small-bore 7T MRI system. Imaging was performed in transaxial plane with 75 × 75 μm² in-plane resolution, 1-mm-slice thickness, and radial DWI using b_max = 42,890 s/mm² . Histological evaluation was performed on upper thoracic sections using optical and transmission electron microscopy. Numerical Monte Carlo simulations (MCSs) of water diffusion were performed to facilitate interpretation of UHb-DWI signal-b curves.

RESULTS

The white matter ultrahigh-b radial DWI (UHb-rDWI) signal-b curves of WT mouse cords behaved biexponentially with high-b diffusion coefficient D_H < 0.020 × 10^-3 mm² /s. However, as expected with less myelination, the signal-b of shiverer mouse cords behaved monoexponentially with significantly greater D_H = 0.162 × 10^-3 , 0.142 × 10^-3 , and 0.164 × 10^-3 mm² /s at anterodorsal, posterodorsal, and lateral columns, respectively. The axial DWI signals of all mouse cords behaved monoexponentially with D = (0.718-1.124) × 10^-3 mm² /s. MCS suggests that these elevated D_H are mainly induced by increased water exchange at the myelin sheath. Microscopic results were consistent with the UHb-rDWI findings.

CONCLUSION

UHb-DWI provides quantitative differences in myelination of spinal cords from myelin-deficit shiverer and WT mice. UHb-DWI may become a powerful tool to evaluate myelination in demyelinating disease models that may translate to human diseases, including multiple sclerosis.

High-b diffusivity of MS lesions in cervical spinal cord using ultrahigh-b DWI (UHb-DWI)

PURPOSE

The purpose of this study was to investigate UHb-rDWI signal in white matter tracts of the cervical spinal cord (CSC) and compare quantitative values between healthy control WM with both MS NAWM and MS WM lesions.

METHODS

UHb-rDWI experiments were performed on (a) 7 MS patients with recently active or chronic lesions in CSC and on (b) 7 healthy control of similar age range and gender distribution to MS subjects. All MRI data were acquired using clinical 3T MRI system. Axial high-b diffusion images were acquired using 2D single-shot DW stimulated EPI with reduced FOV and a CSC-dedicated 8 channel array coil. High-b diffusion coefficient DH was estimated by fitting the signal-b curve to a double or single-exponential function.

RESULTS

The high-b diffusivity DH values were measured as (0.767 ± 0.297) × 10-3 mm2/s in the posterior column lesions, averaged over 6 MS patients, and 0.587 × 10-3 mm2/s in the corticospinal tract for another patient. The averaged DH values of the 7 healthy volunteers from the posterior and lateral column were (0.0312 ± 0.0306) × 10-3 and (0.0505 ± 0.0205) × 10-3 mm2/s, respectively. UHb-rDWI signal-b curves of the MS patients revealed to noticeably behave differently to that of the healthy controls. The patient signal-b curves decayed with greater high-b decay constants to reach lower signal intensities relative to signal-b curves of the healthy controls.

CONCLUSION

UHb-DWI of the CSC reveals a marked difference in signal-b-curves and DH values in MS lesions compared to NAWM and healthy control WM. Based on physical principles, we interpret these altered observations of quantitative diffusion values to be indicative of demyelination. Further studies in animal models will be required to fully interpret UHb-DWI quantitative diffusion values during demyelination and remyelination.

Monte Carlo Simulations of Water Exchange Through Myelin Wraps: Implications for Diffusion MRI

Diffusion magnetic resonance imaging (dMRI) yields parameters sensitive to brain tissue microstructure. A structurally important aspect of this microstructure is the myelin wrapping around the axons. This paper investigated the forward problem concerning whether water exchange via the spiraling structure of the myelin can meaningfully contribute to the signal in dMRI. Monte Carlo simulations were performed in a system with intra-axonal, myelin, and extra-axonal compartments. Diffusion in the myelin was simulated as a spiral wrapping the axon, with a custom number of wraps. Exchange (or intra-axonal residence) times were analyzed for various number of wraps and axon diameters. Pulsed gradient sequences were employed to simulate the dMRI signal, which was analyzed using different methods. Diffusional kurtosis imaging analysis yielded the radial diffusivity (RD) and radial kurtosis (RK), while the two-compartment Kärger model yielded estimates the intra-axonal volume fraction ( ν _ic ) and exchange time ( τ ). Results showed that τ was on the sub-second level for geometries with axon diameters below 1.0 μ m and less than eight wraps. Otherwise, exchange was negligible compared to typical experimental durations, with τ of seconds or longer. In situations where exchange influenced the signal, estimates of RK and ν _ic increased with the number of wraps, while RD decreased. τ estimates from simulated signals were in agreement with predicted ones. In conclusion, exchange through spiraling myelin permits sub-second τ for small diameters and low number of wraps. Such conditions may arise in the developing brain or in neurodegenerative disease, and thus the results could aid the interpretation of dMRI studies.

Ultra-high-b radial diffusion-weighted imaging (UHb-rDWI) of human cervical spinal cord

BACKGROUND

Injury in the cervical spinal cord (CSC) can lead to varying degrees of neurologic deficit and persistent disability. Diffusion tensor imaging (DTI) is a promising method to evaluate white matter integrity and pathology. However, the conventional DTI results are limited with respect to the specific details of neuropathology and microstructural architecture. In this study we used ultrahigh-b radial-DWI (UHb-rDWI) with b-values ranging from 0 to ∼7500 s/mm² and calculated decay constant (D_H ) at the high b-values, which gives much deeper insight about the microscopic environment of CSC white matter.

PURPOSE

To evaluate a novel diffusion MRI, UHb-rDWI technique for imaging of the CSC.

STUDY TYPE

Longitudinal.

SUBJECTS

Four healthy controls, each scanned twice.

FIELD STRENGTH/SEQUENCE

3T/2D single shot diffusion-weighted stimulated echo planar imaging with reduced field of view.

ASSESSMENT

The signal from each pixel of b₀ (b = 0) and b-value (b ≠ 0) images were fitted to a biexponential function and normalized. The signal-b curve is obtained by dividing the latter curve by the former. D_H was obtained from the curve at b >4000 s/mm² . A Monte-Carlo Simulation (MCS) was performed to investigate how D_H changes upon the increased water-exchange at the CSC.

RESULTS

The signal-b curves plotted at multiple levels of healthy CSC are almost identical on two successive scans and show a biexponential decay behavior: fast exponential decay at lower b-values and much slower decay at UHb-values. The mean values of D_H were measured as (0.0607 ± 0.02531) ×10^-3 and (0.0357 ± 0.02072) ×10^-3 s/mm² at the lateral funiculus and posterior column, respectively. MCS of diffusion MRI shows that the D_H is elevated by increased water exchange between the intra- and extraaxonal spaces.

DATA CONCLUSION

UHb-rDWI signal-b plots of the normal CSC were highly reproducible on successive scans and their biexponential decay behavior can be used to characterize normal spinal white matter.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2019;49:204-211.