The impact of realistic axonal shape on axon diameter estimation using diffusion MRI

TL;DR

This paper develops a theory linking axonal shape variations, such as beading and undulations, to diffusion MRI signals, revealing how these features confound accurate axon diameter estimation and proposing methods to mitigate these effects.

Contribution

The study introduces a theoretical framework connecting axonal shape irregularities to diffusion MRI metrics, validated with simulations and microscopy data, highlighting biases in diameter estimation.

Findings

Inner diameter overestimated by about twofold in narrow pulse limit due to shape variations.

Kurtosis deviates from ideal cylinder predictions because of caliber variations.

Undulations significantly affect diameter estimates in the wide pulse limit, but can be reduced by directional averaging.

Abstract

To study axonal microstructure with diffusion MRI, axons are typically modeled as straight impermeable cylinders, whereby the transverse diffusion MRI signal can be made sensitive to the cylinder's inner diameter. However, the shape of a real axon varies along the axon direction, which couples the longitudinal and transverse diffusion of the overall axon direction. Here we develop a theory of the intra-axonal diffusion MRI signal based on coarse-graining of the axonal shape by 3d diffusion. We demonstrate how the estimate of the inner diameter is confounded by the diameter variations (beading), and by the local variations in direction (undulations) along the axon. We analytically relate diffusion MRI metrics, such as time-dependent radial diffusivity D(t) and kurtosis K(t), to the axonal shape, and validate our theory using Monte Carlo simulations in synthetic undulating axons with randomly positioned beads, and in realistic axons reconstructed from electron microscopy images of mouse brain white matter. We show that (i) In the narrow pulse limit, the inner diameter from D(t) is overestimated by about twofold due to a combination of axon caliber variations and undulations (each contributing a comparable effect size); (ii) The narrow-pulse kurtosis K$_{t\to\infty}$ deviates from that in an ideal cylinder due to caliber variations; we also numerically calculate the fourth-order cumulant for an ideal cylinder in the wide pulse limit, which is relevant for inner diameter overestimation; (iii) In the wide pulse limit, the axon diameter overestimation is mainly due to undulations at low diffusion weightings b; and (iv) The effect of undulations can be considerably reduced by directional averaging of high-b signals, with the apparent inner diameter given by a combination of the axon caliber (dominated by the thickest axons), caliber variations, and the residual contribution of undulations.