Comparison of fast and conventional diffusion kurtosis imaging in an anisotropic synthetic phantom
Ganna Blazhenets1,2, Farida Grinberg1,3, Ezequiel Farrher1, Xiang Gao1, Mikheil Kelenjeridze4, Tamo Xechiashvili4, and N. Jon Shah1,3

1Institute of Neuroscience and Medicine - 4, Forschungszentrum Juelich, Juelich, Germany, 2Institute of Nuclear Physics, University of Cologne, Cologne, Germany, 3Department of Neurology, Faculty of Medicine, JARA, RWTH Aachen University, Aachen, Germany, 4Department of Physics, Georgian Technical University, Tbilisi, Georgia

Synopsis

We compare the sensitivity and applicability of two methods for the estimation of mean kurtosis in a multi-sectional, anisotropic diffusion phantom using conventional diffusion kurtosis imaging and a fast protocol for rapid mean kurtosis metric estimation suggested by Hansen et al. (2013). Both methods provide similar image quality and it can be concluded that fast estimation of mean kurtosis is a useful tool that can be used as a fast method for clinical applications. An interesting finding of this work is a stronger dependence of fast computed kurtosis metrics on the orientation of fibres with respect to the static magnetic field than of the conventional method.

Target Audience

Two methods for measuring diffusional kurtosis images were compared regarding their perfomance using an anisotropic synthetic fibre phantom. Our results can be of interest for researchers investigating water diffusion properties in biological tissues.

Purpose

Diffusion magnetic resonance imaging is widely recognized to be a unique, powerful and non-invasive technique for fibrous tissue microstructure investigation, particularly the human brain. Usually, water diffusion in biological tissues is quantified via diffusion tensor imaging (DTI), which uses a Gaussian approximation for the probability distribution of molecular displacements. However, the displacement probability distribution can considerably deviate from the Gaussian form for increasing diffusion weightings ($$$b$$$-values). The propagation of water molecules is impeded by interaction with tissue components, confinements, anisotropy, and is modulated by interfacial interactions with the cell membranes. Due to the heterogeneity and complexity of the tissue microstructure, the various contributions to the average MRI signal in in vivo studies cannot be easily separated from each other. Estimation of the next term in the cumulant expansion of the MRI signal, proportional to the so-called kurtosis, has lately become a popular method to quantify diffusion beyond the Gaussian information offered by DTI. Diffusional mean kurtosis (MK) was suggested to be a sensitive marker for tissue pathology, although, it requires longer acquisition times compared to DTI. A robust and rapid method for estimating MK, with a protocol adapted on commercial scanners, was proposed by Hansen et al.1, and proved to be feasible in clinical use. Several recent studies demonstrated its reproducibility, comparable to the usually obtained MK metrics, in the rat and human brain, as biomarker for glioma subtypes2. The primary purpose of this study was to perform a comparative analysis of the mean kurtosis metrics, obtained with a conventional and a fast protocol for a highly anisotropic medium with a well-known structure, such as a synthetic phantom and to investigate the effect of self-induced susceptibility gradients at different orientations on of the fibres in the static magnetic field.

Materials and methods

The construction of a fibre phantom with a fibre-density gradient region is described in Ref3. Imaging was performed on a whole-body 3T Siemens MAGNETOM Trio Scanner (Siemens, Medical Systems, Erlangen, Germany) using a 12 channel phased-array receive head coil. Diffusion-weighted images were measured with the fibre axis oriented at three angles with the static magnetic field: parallel (0°), perpendicular (90°) and at 45°. The double-refocused diffusion-weighted spin-echo EPI pulse sequence was used for acquiring the signal over 64 field gradient directions. For the fast protocol only 13 diffusion-weighted images (described in Ref.1) were acquired for three $$$b$$$-values (0, 1000 and 2500 s/mm2). The bias in the diffusion-weighted images due to the background noise was corrected using the power-images method4. Conventional DKI parameter maps were computed with ExploreDTI5 for conventional kurtosis imaging and with the help of in-house Matlab scripts (Matlab, The MathWorks) for fast DKI. The number of averages for fast kurtosis was 20 in order to equal the number of diffusion-weighted volumes in conventional DKI experiments.

Results and discussions

Fig.1 shows a comparison of the conventional MK metric and the fast mean kurtosis parameter, $$$\bar{W}$$$, in two structurally different platforms of the phantom. The phantom was positioned in such a way that the fibers were aligned at 45 degrees to the main magnetic field. One can see a similar quality for the images in both methods with minor detectable differences. Fig.1(c) presents a quantitative comparison of the data in both platforms as scatter plots. In voxel by voxel comparison we clearly observe significant correlations between the two maps. However, the spread of points around the $$$MK = \bar{W}$$$ identity line is larger than that reported for brain1, probably due to the larger anisotropic effects of the synthetic fibers. Fig.2 demonstrates the same quantitative and qualitative comparison of both DKI methods for two platforms in the phantom. Here, the phantom was positioned perpendicularly to the direction of a static magnetic field. According to a strong dependence of absorption ability of fibers on the orientation to magnetic field, we see a lower quality of the maps computed by fast kurtosis method in this orientation.

Conclusion

The close correspondence of new measure $$$\bar{W}$$$ of kurtosis with the conventional MK was investigated for an anisotropic synthetic phantom. Fast estimation of diffusion kurtosis metric has the potential to constitute a useful tool for kurtosis imaging and become a beneficial microstructural biomarker. We also observed a strong dependence of the image quality on the orientation of fibers with respect to the magnetic field, that is, the fast method is sensitive to the susceptibility differences arising at interfaces.

Acknowledgements

References

1. B. Hansen, T. Lund, R. Sangill, S. Jespersen, Experimentally and computationally fast method for estimation of a mean kurtosis, Magn. Reson. Med. 00:000-000 (2013); 2. A. Tietze, B. Hansen et al., Diffusional Kurtosis in Patients with Glioma: Initial Results with a Fast Imaging Method in a Clinical Setting, AJNR Am J Neuroradiol. 2015 Aug;36(8):1472-8; 3. E. Farrher, J. Kaffanke, A.A. Celik, T. Stocker, F. Grinberg, N.J. Shah, Novel multisection design of anisotropic diffusion phantoms, Magn. Reson. Imaging, 30 (2012) 518-526; 4. Miller, A.J., Joseph, P.M., The use of power images to perform quantitative analysis on low SNR MR images. Magn. Reson. Imaging 11, 1051–1056. 5.Leemans A, Jeurissen B, Sijbers J, and Jones DK. ExploreDTI: a graphical toolbox for processing, analyzing, and visualizing diffusion MR data.17th Annual Meeting of Intl Soc Mag Reson Med, p. 3537, 2009

Figures

Fig.1 Map of MK with the phantom axes oriented at 45 degree to the magnetic field. Top – region of parallel fibers with a gradient of fiber density along the main fiber direction, bottom – region of parallel fibers of homogeneous fiber density.(a)Selected slice of the mean kurtosis map obtained with a fast protocol for DKI;(b)Selected slice of the mean kurtosis map obtained with conventional DKI;(c)Scatter plots corresponding to the data shown in (a) and (b).The solid line is the identity line.

Fig.2 Map of MK with the phantom axes oriented at 90 degree to the magnetic field. Top – region of parallel fibers with a gradient of fiber density along the main fiber direction, bottom – region of parallel fibers of homogeneous fiber density.(a)Selected slice of the mean kurtosis map obtained with a fast protocol for DKI;(b)Selected slice of the mean kurtosis map obtained with conventional DKI;(c)Scatter plots corresponding to the data shown in (a) and (b).The solid line is the identity line.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
2016