Dataset Information

The effect of gradient nonlinearities on fiber orientation estimates from spherical deconvolution of diffusion magnetic resonance imaging data.

ABSTRACT: Gradient nonlinearities in magnetic resonance imaging (MRI) cause spatially varying mismatches between the imposed and the effective gradients and can cause significant biases in rotationally invariant diffusion MRI measures derived from, for example, diffusion tensor imaging. The estimation of the orientational organization of fibrous tissue, which is nowadays frequently performed with spherical deconvolution techniques ideally using higher diffusion weightings, can likewise be biased by gradient nonlinearities. We explore the sensitivity of two established spherical deconvolution approaches to gradient nonlinearities, namely constrained spherical deconvolution (CSD) and damped Richardson-Lucy (dRL). Additionally, we propose an extension of dRL to take into account gradient imperfections, without the need of data interpolation. Simulations show that using the effective b-matrix can improve dRL fiber orientation estimation and reduces angular deviations, while CSD can be more robust to gradient nonlinearity depending on the implementation. Angular errors depend on a complex interplay of many factors, including the direction and magnitude of gradient deviations, underlying microstructure, SNR, anisotropy of the effective response function, and diffusion weighting. Notably, angular deviations can also be observed at lower b-values in contrast to the perhaps common assumption that only high b-value data are affected. In in vivo Human Connectome Project data and acquisitions from an ultrastrong gradient (300?mT/m) scanner, angular differences are observed between applying and not applying the effective gradients in dRL estimation. As even small angular differences can lead to error propagation during tractography and as such impact connectivity analyses, incorporating gradient deviations into the estimation of fiber orientations should make such analyses more reliable.

SUBMITTER: Guo F

PROVIDER: S-EPMC7776002 | biostudies-literature | 2020 Oct

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

The effect of gradient nonlinearities on fiber orientation estimates from spherical deconvolution of diffusion magnetic resonance imaging data.

Guo Fenghua F de Luca Alberto A Parker Greg G Jones Derek K DK Viergever Max A MA Leemans Alexander A Tax Chantal M W CMW

Human brain mapping 20201009 2

Gradient nonlinearities in magnetic resonance imaging (MRI) cause spatially varying mismatches between the imposed and the effective gradients and can cause significant biases in rotationally invariant diffusion MRI measures derived from, for example, diffusion tensor imaging. The estimation of the orientational organization of fibrous tissue, which is nowadays frequently performed with spherical deconvolution techniques ideally using higher diffusion weightings, can likewise be biased by gradie ...[more]

PMID: 33035372

Similar Datasets

Project description:Diffusion MRI (dMRI) represents one of the few methods for mapping brain fiber orientations non-invasively. Unfortunately, dMRI fiber mapping is an indirect method that relies on inference from measured diffusion patterns. Comparing dMRI results with other modalities is a way to improve the interpretation of dMRI data and help advance dMRI technologies. Here, we present methods for comparing dMRI fiber orientation estimates with optical imaging of fluorescently labeled neurofilaments and vasculature in 3D human and primate brain tissue cuboids cleared using CLARITY. The recent advancements in tissue clearing provide a new opportunity to histologically map fibers projecting in 3D, which represents a captivating complement to dMRI measurements. In this work, we demonstrate the capability to directly compare dMRI and CLARITY in the same human brain tissue and assess multiple approaches for extracting fiber orientation estimates from CLARITY data. We estimate the three-dimensional neuronal fiber and vasculature orientations from neurofilament and vasculature stained CLARITY images by calculating the tertiary eigenvector of structure tensors. We then extend CLARITY orientation estimates to an orientation distribution function (ODF) formalism by summing multiple sub-voxel structure tensor orientation estimates. In a sample containing part of the human thalamus, there is a mean angular difference of 19o±15o between the primary eigenvectors of the dMRI tensors and the tertiary eigenvectors from the CLARITY neurofilament stain. We also demonstrate evidence that vascular compartments do not affect the dMRI orientation estimates by showing an apparent lack of correspondence (mean angular difference = 49o±23o) between the orientation of the dMRI tensors and the structure tensors in the vasculature stained CLARITY images. In a macaque brain dataset, we examine how the CLARITY feature extraction depends on the chosen feature extraction parameters. By varying the volume of tissue over which the structure tensor estimates are derived, we show that orientation estimates are noisier with more spurious ODF peaks for sub-voxels below 30 µm3 and that, for our data, the optimal gray matter sub-voxel size is between 62.5 µm3 and 125 µm3. The example experiments presented here represent an important advancement towards robust multi-modal MRI-CLARITY comparisons.

Project description:Diffusion magnetic resonance imaging (dMRI) is widely used to probe tissue microstructure, and is currently the only non-invasive way to measure the brain's fiber architecture. While a large number of approaches to recover the intra-voxel fiber structure have been utilized in the scientific community, a direct, 3D, quantitative validation of these methods against relevant histological fiber geometries is lacking. In this study, we investigate how well different high angular resolution diffusion imaging (HARDI) models and reconstruction methods predict the ground-truth histologically defined fiber orientation distribution (FOD), as well as investigate their behavior over a range of physical and experimental conditions. The dMRI methods tested include constrained spherical deconvolution (CSD), Q-ball imaging (QBI), diffusion orientation transform (DOT), persistent angular structure (PAS), and neurite orientation dispersion and density imaging (NODDI) methods. Evaluation criteria focus on overall agreement in FOD shape, correct assessment of the number of fiber populations, and angular accuracy in orientation. In addition, we make comparisons of the histological orientation dispersion with the fiber spread determined from the dMRI methods. As a general result, no HARDI method outperformed others in all quality criteria, with many showing tradeoffs in reconstruction accuracy. All reconstruction techniques describe the overall continuous angular structure of the histological FOD quite well, with good to moderate correlation (median angular correlation coefficient > 0.70) in both single- and multiple-fiber voxels. However, no method is consistently successful at extracting discrete measures of the number and orientations of FOD peaks. The major inaccuracies of all techniques tend to be in extracting local maxima of the FOD, resulting in either false positive or false negative peaks. Median angular errors are ?10° for the primary fiber direction and ?20° for the secondary fiber, if present. For most methods, these results did not vary strongly over a wide range of acquisition parameters (number of diffusion weighting directions and b value). Regardless of acquisition parameters, all methods show improved successes at resolving multiple fiber compartments in a voxel when fiber populations cross at near-orthogonal angles, with no method adequately capturing low to moderate angle (<60°) crossing fibers. Finally, most methods are limited in their ability to capture orientation dispersion, resulting in low to moderate, yet statistically significant, correlation with histologically-derived dispersion with both HARDI and NODDI methodologies. Together, these results provide quantitative measures of the reliability and limitations of dMRI reconstruction methods and can be used to identify relative advantages of competing approaches as well as potential strategies for improving accuracy.

Dataset Information

The effect of gradient nonlinearities on fiber orientation estimates from spherical deconvolution of diffusion magnetic resonance imaging data.

Publications

The effect of gradient nonlinearities on fiber orientation estimates from spherical deconvolution of diffusion magnetic resonance imaging data.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets