ObjectiveConventional magnetic resonance imaging (MRI) of the multiple sclerosis spinal cord is limited by low specificity regarding the underlying pathological processes, and new MRI metrics assessing microscopic damage are required. We aim to show for the first time that neurite orientation dispersion (i.e., variability in axon/dendrite orientations) is a new biomarker that uncovers previously undetected layers of complexity of multiple sclerosis spinal cord pathology. Also, we validate against histology a clinically viable MRI technique for dispersion measurement (neurite orientation dispersion and density imaging, NODDI), to demonstrate the strong potential of the new marker.MethodsWe related quantitative metrics from histology and MRI in four post mortem spinal cord specimens (two controls; two progressive multiple sclerosis cases). The samples were scanned at high field, obtaining maps of neurite density and orientation dispersion from NODDI and routine diffusion tensor imaging (DTI) indices. Histological procedures provided markers of astrocyte, microglia, myelin and neurofilament density, as well as neurite dispersion.ResultsWe report from both NODDI and histology a trend toward lower neurite dispersion in demyelinated lesions, indicative of reduced neurite architecture complexity. Also, we provide unequivocal evidence that NODDI‐derived dispersion matches its histological counterpart (P < 0.001), while DTI metrics are less specific and influenced by several biophysical substrates.InterpretationNeurite orientation dispersion detects a previously undescribed and potentially relevant layer of microstructural complexity of multiple sclerosis spinal cord pathology. Clinically feasible techniques such as NODDI may play a key role in clinical trial and practice settings, as they provide histologically meaningful dispersion indices.

This work introduces a compartment-based model for apparent cell body (namely soma) and neurite density imaging (SANDI) using non-invasive diffusion-weighted MRI (DW-MRI). The existing conjecture in brain microstructure imaging trough DW-MRI presents water diffusion in white (WM) and grey (GM) matter as restricted diffusion in neurites, modelled by infinite cylinders of null radius embedded in the hindered extra-neurite water. The extra-neurite pool in WM corresponds to water in the extra-axonal space, but in GM it combines water in the extra-cellular space with water in soma. While several studies showed that this microstructure model successfully describe DW-MRI data in WM and GM at b3 ms/m 2 , it has been also shown to fail in GM at high b values (b>>3 ms/m 2 ). Here we hypothesize that the unmodelled soma compartment may be responsible for this failure and propose SANDI as a new model of brain microstructure where soma is explicitly included. We assess the effects of size and density of soma on the direction-averaged DW-MRI signal at high b values and the regime of validity of the model using numerical simulations and comparison with experimental data from mouse (bmax = 40 ms/m 2 ) and human (bmax = 10 ms/m 2 ) brain. We show that SANDI defines new contrasts representing new complementary information on the brain cyto-and myelo-architecture. Indeed, we show for the first-time maps from 25 healthy human subjects of MR soma and neurite signal fractions, that remarkably mirror contrasts of histological images of brain cyto-and myelo-architecture. Although still under validation, SANDI might provide new insight into tissue architecture by introducing a new set of biomarkers of potential great value for biomedical applications and pure neuroscience.

PurposeTo identify optimal pulsed gradient spin‐echo (PGSE) and oscillating gradient spin‐echo (OGSE) sequence settings for maximizing sensitivity to axon diameter in idealized and practical conditions.MethodsSimulations on a simple two‐compartment white matter model (with nonpermeable cylinders) are used to investigate a wide space of clinically plausible PGSE and OGSE sequence parameters with trapezoidal diffusion gradient waveforms. Signal sensitivity is measured as a derivative of the signal with respect to axon diameter. Models of parallel and dispersed fibers are investigated separately to represent idealized and practical conditions.ResultsSimulations show that, for the simple case of gradients perfectly perpendicular to straight parallel fibers, PGSE always gives maximum sensitivity. However, in real‐world scenarios where fibers have unknown and dispersed orientation, low‐frequency OGSE provides higher sensitivity. Maximum sensitivity results show that on current clinical scanners (G max = 60 mT/m, signal to noise ratio (SNR) = 20) axon diameters below 6 µm are indistinguishable from zero. Scanners with stronger gradient systems such as the Massachusetts General Hospital (MGH) Connectom scanner (G max = 300 mT/m) can extend this sensitivity limit down to 2–3 µm, probing a much greater proportion of the underlying axon diameter distribution.ConclusionLow‐frequency OGSE provides additional sensitivity to PGSE in practical situations. OGSE is particularly advantageous for systems with high performance gradients. Magn Reson Med 75:688–700, 2016. © 2015 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Microscopic diffusion anisotropy (μA) has been recently gaining increasing attention for its ability to decouple the average compartment anisotropy from orientation dispersion. Advanced diffusion MRI sequences, such as double diffusion encoding (DDE) and double oscillating diffusion encoding (DODE) have been used for mapping μA, usually using measurements from a single b shell. However, the accuracy of μA estimation vis-à-vis different b-values was not assessed. Moreover, the time-dependence of this metric, which could offer additional insights into tissue microstructure, has not been studied so far. Here, we investigate both these concepts using theory, simulation, and experiments performed at 16.4T in the mouse brain, ex-vivo. In the first part, simulations and experimental results show that the conventional estimation of microscopic anisotropy from the difference of D(O)DE sequences with parallel and orthogonal gradient directions yields values that highly depend on the choice of b-value. To mitigate this undesirable bias, we propose a multi-shell approach that harnesses a polynomial fit of the signal difference up to third order terms in b-value. In simulations, this approach yields more accurate μA metrics, which are similar to the ground-truth values. The second part of this work uses the proposed multi-shell method to estimate the time/frequency dependence of μA. The data shows either an increase or no change in μA with frequency depending on the region of interest, both in white and gray matter. When comparing the experimental results with simulations, it emerges that simple geometric models such as infinite cylinders with either negligible or finite radii cannot replicate the measured trend, and more complex models, which, for example, incorporate structure along the fibre direction are required. Thus, measuring the time dependence of microscopic anisotropy can provide valuable information for characterizing tissue microstructure.