Biophysical compartment models for single-shell diffusion MRI in the human brain: a model fitting comparison

Clinically oriented studies commonly acquire diffusion MRI (dMRI) data with a single non-zero b-value (i.e. single-shell) and diffusion weighting of b=1000 s/mm2. To produce microstructural parameter maps, the tensor model is usually used, despite known limitations. Although compartment models have demonstrated improved fits in multi-shell dMRI data, they are rarely used for single-shell parameter maps, where their effectiveness is unclear from the literature. Here, various compartment models combining isotropic balls and symmetric tensors were fitted to single-shell dMRI data to investigate model fitting optimization and extract the most information possible. Full testing was performed in 5 subjects, and 3 subjects with multi-shell data were included for comparison. The results were tested and confirmed in a further 50 subjects. The Markov chain Monte Carlo (MCMC) model fitting technique outperformed non-linear least squares. Using MCMC, the 2-fibre-orientation mono-exponential ball & stick model (BSME 2) provided artifact-free, stable results, in little processing time. The analogous ball & zeppelin model (BZ2) also produced stable, low-noise parameter maps, though it required much greater computing resources (50 000 burn-in steps). In single-shell data, the gamma-distributed diffusivity ball & stick model (BSGD 2) underperformed relative to other models, despite being an often-used software default. It produced artifacts in the diffusivity maps even with extremely long processing times. Neither increased diffusion weighting nor a greater number of gradient orientations improved BSGD 2 fits. In white matter (WM), the tensor produced the best fit as measured by Bayesian information criterion. This result contrasts with studies using multi-shell data. However, in crossing fibre regions the tensor confounded geometric effects with fractional anisotropy (FA): the planar/linear WM FA ratio was 49%, while BZ2 and BSME 2 retained 76% and 83% of restricted fraction, respectively. As a result, the BZ2 and BSME 2 models are strong candidates to optimize information extraction from single-shell dMRI studies.