Evaluation of fiber tracking from subsampled q-space data in diffusion spectrum imaging

Research Progress in Diffusion Spectrum Imaging

Studies have demonstrated that many regions in the human brain include multidirectional fiber tracts, in which the diffusion of water molecules within image voxels does not follow a Gaussian distribution. Therefore, the conventional diffusion tensor imaging (DTI) that hypothesizes a single fiber orientation within a voxel is intrinsically incapable of revealing the complex microstructures of brain tissues. Diffusion spectrum imaging (DSI) employs a pulse sequence with different b-values along multiple gradient directions to sample the diffusion information of water molecules in the entire q-space and then quantitatively estimates the diffusion profile using a probability density function with a high angular resolution. Studies have suggested that DSI can reliably observe the multidirectional fibers within each voxel and allow fiber tracking along different directions, which can improve fiber reconstruction reflecting the true but complicated brain structures that were not observed in the previous DTI studies. Moreover, with increasing angular resolution, DSI is able to reveal new neuroimaging biomarkers used for disease diagnosis and the prediction of disorder progression. However, so far, this method has not been used widely in clinical studies, due to its overly long scanning time and difficult post-processing. Within this context, the current paper aims to conduct a comprehensive review of DSI research, including the fundamental principles, methodology, and application progress of DSI tractography. By summarizing the DSI studies in recent years, we propose potential solutions towards the existing problem in the methodology and applications of DSI technology as follows: (1) using compressed sensing to undersample data and to reconstruct the diffusion signal may be an efficient and promising method for reducing scanning time; (2) the probability density function includes more information than the orientation distribution function, and it should be extended in application studies; and (3) large-sample study is encouraged to confirm the reliability and reproducibility of findings in clinical diseases. These findings may help deepen the understanding of the DSI method and promote its development in clinical applications.

Generalized diffusion spectrum magnetic resonance imaging (GDSI) for model-free reconstruction of the ensemble average propagator

Diffusion spectrum MRI (DSI) provides model-free estimation of the diffusion ensemble average propagator (EAP) and orientation distribution function (ODF) but requires the diffusion data to be acquired on a Cartesian q-space grid. Multi-shell diffusion acquisitions are more flexible and more commonly acquired but have, thus far, only been compatible with model-based analysis methods. Here, we propose a generalized DSI (GDSI) framework to recover the EAP from multi-shell diffusion MRI data. The proposed GDSI approach corrects for q-space sampling density non-uniformity using a fast geometrical approach. The EAP is directly calculated in a preferable coordinate system by multiplying the sampling density corrected q-space signals by a discrete Fourier transform matrix, without any need for gridding. The EAP is demonstrated as a way to map diffusion patterns in brain regions such as the thalamus, cortex and brainstem where the tissue microstructure is not as well characterized as in white matter. Scalar metrics such as the zero displacement probability and displacement distances at different fractions of the zero displacement probability were computed from the recovered EAP to characterize the diffusion pattern within each voxel. The probability averaged across directions at a specific displacement distance provides a diffusion property based image contrast that clearly differentiates tissue types. The displacement distance at the first zero crossing of the EAP averaged across directions orthogonal to the primary fiber orientation in the corpus callosum is found to be larger in the body (5.65 ± 0.09 μm) than in the genu (5.55 ± 0.15 μm) and splenium (5.4 ± 0.15 μm) of the corpus callosum, which corresponds well to prior histological studies. The EAP also provides model-free representations of angular structure such as the diffusion ODF, which allows estimation and comparison of fiber orientations from both the model-free and model-based methods on the same multi-shell data. For the model-free methods, detection of crossing fibers is found to be strongly dependent on the maximum b-value and less sensitive compared to the model-based methods. In conclusion, our study provides a generalized DSI approach that allows flexible reconstruction of the diffusion EAP and ODF from multi-shell diffusion data and data acquired with other sampling patterns.

White matter changes in patients with mild traumatic brain injury: MRI perspective

This review focuses on white matter (WM) changes in mild traumatic brain injury (mTBI) as assessed by multimodal MRI. All the peer reviewed publications on WM changes in mTBI from January 2011 through September 2016 are included in this review. This review is organized as follows: introduction to mTBI, the basics of multimodal MRI techniques that are potentially useful for probing the WM integrity, summary and critical evaluation of the published literature on the application of multimodal MRI techniques to assess the changes of WM in mTBI, and correlation of MRI measures with behavioral deficits. The MRI–pathology correlation studies based on preclinical models of mTBI are also reviewed. Finally, the author's perspective of future research directions is described.

Q-space truncation and sampling in diffusion spectrum imaging

PURPOSE

To characterize the q-space truncation and sampling on the spin-displacement probability density function (PDF) in diffusion spectrum imaging (DSI).

METHODS

DSI data were acquired using the MGH-USC connectome scanner (Gmax  = 300 mT/m) with bmax  = 30,000 s/mm2 , 17 × 17 × 17, 15 × 15 × 15 and 11 × 11 × 11 grids in ex vivo human brains and bmax  = 10,000 s/mm2 , 11 × 11 × 11 grid in vivo. An additional in vivo scan using bmax =7,000 s/mm2 , 11 × 11 × 11 grid was performed with a derated gradient strength of 40 mT/m. PDFs and orientation distribution functions (ODFs) were reconstructed with different q-space filtering and PDF integration lengths, and from down-sampled data by factors of two and three.

RESULTS

Both ex vivo and in vivo data showed Gibbs ringing in PDFs, which becomes the main source of artifact in the subsequently reconstructed ODFs. For down-sampled data, PDFs interfere with the first replicas or their ringing, leading to obscured orientations in ODFs.

CONCLUSION

The minimum required q-space sampling density corresponds to a field-of-view approximately equal to twice the mean displacement distance (MDD) of the tissue. The 11 × 11 × 11 grid is suitable for both ex vivo and in vivo DSI experiments. To minimize the effects of Gibbs ringing, ODFs should be reconstructed from unfiltered q-space data with the integration length over the PDF constrained to around the MDD. Magn Reson Med 76:1750-1763, 2016. © 2016 International Society for Magnetic Resonance in Medicine.

Mapping the Orientation of White Matter Fiber Bundles: A Comparative Study of Diffusion Tensor Imaging, Diffusional Kurtosis Imaging, and Diffusion Spectrum Imaging

BACKGROUND AND PURPOSE

White matter fiber tractography relies on fiber bundle orientation estimates from diffusion MR imaging. However, clinically feasible techniques such as DTI and diffusional kurtosis imaging use assumptions, which may introduce error into in vivo orientation estimates. In this study, fiber bundle orientations from DTI and diffusional kurtosis imaging are compared with diffusion spectrum imaging as a criterion standard to assess the performance of each technique.

MATERIALS AND METHODS

For each subject, full DTI, diffusional kurtosis imaging, and diffusion spectrum imaging datasets were acquired during 2 independent sessions, and fiber bundle orientations were estimated by using the specific theoretic assumptions of each technique. Angular variability and angular error measures were assessed by comparing the orientation estimates. Tractography generated with each of the 3 reconstructions was also examined and contrasted.

RESULTS

Orientation estimates from all 3 techniques had comparable angular reproducibility, but diffusional kurtosis imaging decreased angular error throughout the white matter compared with DTI. Diffusion spectrum imaging and diffusional kurtosis imaging enabled the detection of crossing-fiber bundles, which had pronounced effects on tractography relative to DTI. Diffusion spectrum imaging had the highest sensitivity for detecting crossing fibers; however, the diffusion spectrum imaging and diffusional kurtosis imaging tracts were qualitatively similar.

CONCLUSIONS

Fiber bundle orientation estimates from diffusional kurtosis imaging have less systematic error than those from DTI, which can noticeably affect tractography. Moreover, tractography obtained with diffusional kurtosis imaging is qualitatively comparable with that of diffusion spectrum imaging. Because diffusional kurtosis imaging has a shorter typical scan time than diffusion spectrum imaging, diffusional kurtosis imaging is potentially more suitable for a variety of clinical and research applications.

Strengths and weaknesses of state of the art fiber tractography pipelines--A comprehensive in-vivo and phantom evaluation study using Tractometer

Many different tractography approaches and corresponding isolated evaluation attempts have been presented over the last years, but a comparative and quantitative evaluation of tractography algorithms still remains a challenge, particularly in-vivo. The recently presented evaluation framework Tractometer is the first attempt to approach this challenge in a quantitative, comparative, persistent and open-access way. Tractometer is currently based on the evaluation of several global connectivity and tract-overlap metrics on hardware phantom data. The work presented in this paper focuses on extending Tractometer with a metric that enables the assessment of the local consistency of tractograms with the underlying image data that is not only applicable to phantom dataset but allows the quantitative and purely data-driven evaluation of in-vivo tractography. We furthermore present an extensive reference-based evaluation study of 25,000 tractograms obtained on phantom and in-vivo datasets using the presented local metric as well as all the methods already established in Tractometer. The experiments showed that the presented local metric successfully reflects the behavior of in-vivo tractography under different conditions and that it is consistent with the results of previous studies. Additionally our experiments enabled a multitude of conclusions with implications for fiber tractography in general, including recommendations regarding optimal choice of a local modeling technique, tractography algorithm, and parameterization, confirming and complementing the results of earlier studies.

Simultaneous multi-scale diffusion estimation and tractography guided by entropy spectrum pathways

We have developed a method for the simultaneous estimation of local diffusion and the global fiber tracts based upon the information entropy flow that computes the maximum entropy trajectories between locations and depends upon the global structure of the multi-dimensional and multi-modal diffusion field. Computation of the entropy spectrum pathways requires only solving a simple eigenvector problem for the probability distribution for which efficient numerical routines exist, and a straight forward integration of the probability conservation through ray tracing of the convective modes guided by a global structure of the entropy spectrum coupled with a small scale local diffusion. The intervoxel diffusion is sampled by multi b-shell multi q-angle diffusion weighted imaging data expanded in spherical waves. This novel approach to fiber tracking incorporates global information about multiple fiber crossings in every individual voxel and ranks it in the most scientifically rigorous way. This method has potential significance for a wide range of applications, including studies of brain connectivity.

Function-structure associations of the brain: evidence from multimodal connectivity and covariance studies

Despite significant advances in multimodal imaging techniques and analysis approaches, unimodal studies are still the predominant way to investigate brain changes or group differences, including structural magnetic resonance imaging (sMRI), functional MRI (fMRI), diffusion tensor imaging (DTI) and electroencephalography (EEG). Multimodal brain studies can be used to understand the complex interplay of anatomical, functional and physiological brain alterations or development, and to better comprehend the biological significance of multiple imaging measures. To examine the function-structure associations of the brain in a more comprehensive and integrated manner, we reviewed a number of multimodal studies that combined two or more functional (fMRI and/or EEG) and structural (sMRI and/or DTI) modalities. In this review paper, we specifically focused on multimodal neuroimaging studies on cognition, aging, disease and behavior. We also compared multiple analysis approaches, including univariate and multivariate methods. The possible strengths and limitations of each method are highlighted, which can guide readers when selecting a method based on a given research question. In particular, we believe that multimodal fusion approaches will shed further light on the neuronal mechanisms underlying the major structural and functional pathophysiological features of both the healthy brain (e.g. development) or the diseased brain (e.g. mental illness) and, in the latter case, may provide a more sensitive measure than unimodal imaging for disease classification, e.g. multimodal biomarkers, which potentially can be used to support clinical diagnosis based on neuroimaging techniques.