Maximum b-value dependence of Diffusion kurtosis imaging sensitivity in detecting white matter microstructure
Miao Sha1, Yuanyuan Chen1, Xin Zhao1, Man Sun2, Weiwei Wang1, Hongyan Ni2, and Dong Ming1

1Tianjin University, Tianjin, China, People's Republic of, 2Tianjin First Center Hospital, Tianjin, China, People's Republic of

Synopsis

Diffusion kurtosis imaging is a powerful technique to measure the non-gaussion diffusion as well as the complicated microstructure. In this paper, we conducted a comparison between different acquisitions with different maximum b-value on normal volunteers. We found that the outcome of diffusion kurtosis imaging was influenced by the maximum b-value in the acquisition. And this influence was highly associated with the microstructure, including both radial profile and angular profile in the structure reconstruction, which indicated the mechanism of non-gaussion under high b-value.

Purpose

Diffusion kurtosis imaging (DKI) was suggested to describe the exact non-Gaussian diffusion in tissue and also provided accurate and abundant information of the microstructure1-3. This paper is to investigate the critical b-value in the acquisition and reveal the influence on the kurtosis representation of white matter microstructure with diffusion kurtosis imaging (DKI).

Methods

Nine normal healthy volunteers (4 male, 4 female) aging from 23 to 26 participated and were scanned in a Siemens 3.0 T MRI Scanner System using a 32-channel head coil. A six b-values DWI acquisition including 0, 500, 1000, 1500, 2000, 2500 sec/mm2 with 30 diffusion encoded orientations was conducted on each individal. Other acquisition parameters: TR=6800ms, TE=101ms, voxel size=2.0×2.0×2.0mm3, FOV=256×256 mm2, metrics=128×128, number of sagittal slices=52 with no gap. And a high resolution 3D T1-weighted MPRAGE was acquired with the following imaging parameters: TR = 1900 ms, TE = 2.52 ms, matrix = 240×256, FOV = 240×256 mm2, and slice thickness = 1 mm. After image quality assessment, one participate was excluded due to poor quality of imaging data. After eddy current and head motion correction were finished in FSL sofrware, Diffusional Kurtosis Estimator (DKE) was ultilized to calculated kurtosis indexes radial kurtosis(RK), axial kurtosis(AK) and mean kurtosis(MK) from the three combinations of b-value (BL: 0,1000 sec/mm2,2000 sec/mm2; BH:0,1000sec/mm2,2500sec/mm2 and TA: 0, 500, 1000, 1500, 2000, 2500 sec/mm2) DWIs data. Tract-based spatial statistics analysis for whole brain was firstly performed on RK and AK under significance of FDR corrected p<0.05 in FSL software. Then the 3D geometrical shape analysis of kurtosis for corpus callosum were performed, in which the corpus callosum was divided into four parts. For each part, the averaged 3D spatial kurtosis distribution was reconstructed and graphied for detailed microstructure.

Results

Visually from Fig.1, the post cingulate white matter and splenium of corpus callosum revealed differences in two kurtosis indexes, and more significant differences were observed in AK map than MK maps, especially in the corpus callosum, posterior cingulate white matter and bilateral posterior thalamic radiation (especially left) . Predominantly, two-paired t-test of the differences between BL and TA contributes more than that BH and TA. Note the MK statistic-analysis, in this case, almost none differences have been factored out. As diffusion tensor had little difference between these protocols, Fig.2 is the 3D surface of the kurtosis tensors’ cross section in CC over all orientations. All the blue lines directed the principal diffusion orientation which represented the fiber direction. Some detailed difference was obvious that although the orientation and radius profile all showed some characteristic, the Fig.2-EFGH was more similar to the Fig.2-IJKL, while the Fig.2-ABCD indicated higher variability compared to the other two lines in Fig.2. In summary, the radial kurtosis and axial kurtosis from different maximum b-value acquisitions lead a significant difference to kurtosis estimation in white matter mainly on corpus callosum and corona radiate. The specific microstructure and the geometric shape of white matter revealed specific sensitivity and variability with the varied b-values.

Discussion

Two diffusion weighted b-values acquisition was popular in DKI application since that one low b-value and one high b-value was demonstrated to be better fitted the diffusion decay4. Also 2000 sec/mm2 and 2500 sec/mm2 were mostly used in DKI research as maximum b-value5,6. We compared BH and BL to explore the influence of this variation of maximum b-value on the kurtosis estimation. As an control experiment, TA including all b-values was compared with the considered variation of b-values' number and also this six b-value acquisition was popularly used2,3,6. Based on this consideration, the results above might confirmed that the critical b-value was the maximum. The significant regions were located in the corpus callosum and corona radiata which reflected that this maximum b-value influence was associated with the microstruture. For the detailed graphic in Fig 2, radial profile in the fiber direction and angular profile in the myelin direction were showed difference caused by the different maximum b-value. This was agree with the previous research7 that the maximum b-value could not impact the fiber orientaion, but the fiber radius profile and the fiber longitudinal profile which was found in this paper. High b-value was highly sensitive to the non-Gaussion as well as the microstructure8,9, and the relation between them can obtain a further insight into the non-Gaussion diffusion mechanism and the intrinsic principles of DKI4-6 .

Conclusion

What can be concluded that maximum b-value is a key factor that decides the kurtosis representation of white matter, and the proper value should vary from specific microstructure.

Acknowledgements

This paper was supported by National Natural Science Foundation of China (No. 81571762, 81222021, 31500865), National Key Technology R&D Program of the Ministry of Science and Technology of China (No. 2012BAI34B02),the Tianjin Research Program of Application Foundation and Advanced Technology (13JCQNJC14400) and the Tianjin Bureau of Public Health Foundation (No. 09KY10, 11KG108).

References

[1]Jensen J H, Helpern J A, Ramani A, et al. Diffusional kurtosis imaging: The quantification of non-gaussian water diffusion by means of magnetic resonance imaging[J]. Magnetic Resonance in Medicine, 2005, 53(6): 1432-1440.

[2]Jensen J H, Helpern J A. MRI quantification of non-Gaussian water diffusion by kurtosis analysis[J]. NMR in Biomedicine, 2010, 23(7): 698-710.

[3]Hui E S, Cheung M M, Qi L, et al. Towards better MR characterization of neural tissues using directional diffusion kurtosis analysis[J]. Neuroimage, 2008, 42(1): 122-134.

[4]Zhu J, Zhuo C, Qin W, et al. Performances of diffusion kurtosis imaging and diffusion tensor imaging in detecting white matter abnormality in schizophrenia[J]. NeuroImage: Clinical, 2015, 7: 170-176.

[5]Merisaari H, Jambor I. Optimization of b-value distribution for four mathematical models of prostate cancer diffusion-weighted imaging using b values up to 2000 s/mm2: Simulation and repeatability study[J]. Magnetic Resonance in Medicine, 2015, 73(5): 1954-1969.

[6]Hansen B, Lund T E, Sangill R, et al. Experimentally and computationally fast method for estimation of a mean kurtosis [J]. Magnetic Resonance in Medicine, 2013, 69(6): 1754-1760.

[7]Neto H R, Correia M M, Nunes R G, et al. Exploring the 3D geometry of the diffusion kurtosis tensor--impact on the development of robust tractography procedures and novel biomarkers.[J]. Neuroimage, 2015, 111:85-99.

[8]Glenn G R, Helpern J A, Tabesh A, et al. Optimization of white matter fiber tractography with diffusional kurtosis imaging[J]. NMR in Biomedicine, 2015, 28(10): 1245-1256.

Figures

Figure 1. Regions of significant MK, AK difference between the three protocols BL, BH and TA. The columns a and b shows F-test results of the three, the column c and d shows T-test results between BL and TA and column e and f shows T-test results between BH and TA (FDR corrected p<0.05).

Figure 2. Morphological contrast of kurtosis distribution along the corpus callosum between BL (plane A-D), BH (plane E-H) and TA (plane I-L).



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