Transient Anomalous Diffusion MRI in Excised Mouse Spinal Cord: Comparison Among Different Diffusion Metrics and Validation With Histology

Neural tissue is a hierarchical multiscale system with intracellular and extracellular diffusion compartments at different length scales. The normal diffusion of bulk water in tissues is not able to detect the specific features of a complex system, providing nonlocal, diffusion measurement averaged on a 10-20 mm length scale. Being able to probe tissues with sub-micrometric diffusion length and quantify new local parameters, transient anomalous diffusion (tAD) would dramatically increase the diagnostic potential of diffusion MRI (DMRI) in detecting collective and sub-micro architectural changes of human tissues due to pathological damage. In DMRI, the use of tAD parameters quantified using specific DMRI acquisition protocols and their interpretation has often aroused skepticism. Although the derived formulas may accurately fit experimental diffusion-weighted data, the relationships between the postulated dynamical feature and the underlying geometrical structure remains elusive, or at most only suggestive. This work aimed to elucidate and validate the image contrast and information that can be obtained using the tAD model in white matter (WM) through a direct comparison between different diffusion metrics and histology. Towards this goal, we compared tAD metrics extracted from pure subdiffusion (a-imaging) and superpseudodiffusion (g-imaging) in excised mouse spinal cord WM, together with T2 and T2  relaxometry, conventional (normal diffusion-based) diffusion tensor imaging (DTI) and q-space imaging (QSI), with morphologic measures obtained by optical microscopy, to determine which structural and topological characteristics of myelinated axons influenced tAD contrast. Axon diameter (AxDiam), the standard deviation of diameters (SDax:diam), axonal density (AxDens) and effective local density (ELD) were extracted from optical images in several WM tracts. Among all the diffusion parameters obtained at 9.4 T, g-metrics confirmed a strong dependence on magnetic in-homogeneities quantified by R2  = 1/T2  and showed the strongest associations with AxDiam and ELD. On the other hand, a-metrics showed strong associations with SDax:diam and was significantly related to AxDens, suggesting its ability to quantify local heterogeneity degree in neural tissue. These results elucidate the biophysical mechanism underpinning tAD parameters and show the clinical potential of tAD-imaging, considering that both physiologic and pathologic neurodegeneration translate into alterations of WM morphometry and topology.