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
microstructure
1-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/mm
2 with 30 diffusion
encoded orientations was conducted on each individal. Other acquisition parameters: TR=6800ms, TE=101ms, voxel
size=2.0×2.0×2.0mm
3, FOV=256×256 mm
2, 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 mm
2, 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/mm
2,2000 sec/mm
2; BH:0,1000sec/mm
2,2500sec/mm
2
and TA: 0, 500, 1000, 1500, 2000, 2500
sec/mm
2) 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
decay
4. Also 2000 sec/mm
2 and 2500 sec/mm
2 were mostly used in DKI research as maximum b-value
5,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 used
2,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 research
7 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 microstructure
8,9, and the
relation between them can obtain a further insight into the non-Gaussion
diffusion mechanism and the intrinsic principles of DKI
4-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.