Lecture Notes Juntendo Medical Journal 2017. 63 (4), 285-292
Diffusion Magnetic Resonance Imaging: From Isotropic Diffusion-Weighted Imaging to Diffusion Tensor Imaging and Beyond CHRISTINA ANDICA*, ASAMI SAITO*, SYO MURATA*, AKI HATTORI*, YUTAKA IKENOUCHI*, MASAAKI HORI*, SHIGEKI AOKI* *Department
of Radiology, Juntendo University Graduate School of Medicine, Tokyo, Japan
Diffusion magnetic resonance imaging (dMRI), which was established in the mid 1980ʼs, is an imaging technique that is based on the diffusion of water molecules in tissue. Initially it was an isotropic diffusion-weighted imaging (DWI) technique, and has been widely used mainly for investigating the tissue microstructural changes in neurological abnormalities, especially ischemic stroke. Many studies have been conducted to validate the usefulness of DWI in clinical settings and to improve the technique. Diffusion tensor imaging (DTI) was later introduced to deal with the anisotropic diffusion; the use of an ellipsoid tensor model means that the size and the direction of the water diffusion can be delineated. DTI and its developments are becoming powerful tools, and have become the standard imaging techniques for analyzing white matter fiber structure and connectivity in vivo. Furthermore, to facilitate researchers in applying DTI for multi-subject studies, tract-based spatial statistics (TBSS) was developed. TBSS is a voxelwise statistical analysis that was originally used to evaluate the fractional anisotropy (FA) values in the white matter. Thereafter, diffusion kurtosis imaging (DKI) was developed to overcome the limitations of DTI such as the non-Gaussian characteristics of tissue and the crossing fiber problems. The latest promising technique is neurite orientation dispersion and density imaging (NODDI), which can be used to show the microstructural changes of brain tissue more sensitively and specifically. Key words: diffusion MRI, diffusion tensor imaging
in vivo 4)-6). More advanced diffusion techniques will be also described in this manuscript.
Introduction Diffusion magnetic resonance imaging (dMRI) has been used widely in the clinical setting to analyze the microstructural changes of tissue based on the movement of water molecules 1). Nowadays, the clinical usage of dMRI is mainly to investigate neurological abnormalities, especially for ischemic stroke 1)-3). Diffusion tensor imaging (DTI) and DT tractography, which were developed from dMRI, are now becoming powerful tools and have become the standard imaging techniques for analyzing white matter fiber structure and brain connectivity
The history and basic principles of dMRI Molecular diffusion or Brownian motion was proposed by Einstein in 1905 as a random motion of molecules in a fluid that is induced by thermal energy 1) 6). In the mid 1980ʼ s, Le Bihan proposed that microstructural information of tissue could be delineated by imaging the water diffusion through intravoxel incoherent motion (IVIM) MR Imaging 1)-3). To produce diffusion-weighted images
Corresponding author: Christina Andica Department of Radiology, Juntendo University Graduate School of Medicine 2-1-1 Hongo, Bunkyo-ku, Tokyo, 113-8421, Japan TEL: + 81-3-5802-1230 〔Received
FAX: + 81-3-3816-0958
Jan. 16, 2017〕〔Accepted
E-mail:
[email protected]
Feb. 22, 2017〕
Copyright © 2017 The Juntendo Medical Society. This is an open access article distributed under the terms of Creative Commons Attribution License (CC BY), which permits unrestricted use, distribution, and reproduction in any medium, provided the original source is properly credited. doi: 10.14789/jmj.63.285
285
Andica, et al: Diffusion MRI: isotropic DWI and beyond
90 RF
180 RF
echo
MPG
MPG
G δ A)Stable spins
b=γ2G2δ2(Δ−δ/3)
δ
Δ 180°
Refocus Signal intensity
Phase B)Moving spins
180°
random walk
Failure of refocus
Signal intensity Phase
Figure-1
Sequence chart of diffusion
To produce diffusion weighted images, a pair of motion probing gradients (MPGs) that are equal in magnitude but have an opposite effect in direction are added to the spin echo-echo planar imaging (EPI) sequence. In a voxel of tissue with stable spins (A), the opposite effect of gradient pulses will cancel each other and a strong signal will be produced. However, moving spins (B) will strongly affected by the opposite effect of gradient pulses in the different locations. The two gradients no longer have the same effect due to the different locations, and they cannot cancel each other, and thus little signal intensity is produced 7). (Kore de wakaru Kakusan MRI, 3rd Edition, 2013: 16-55 7))
(DWI), Le Bihan applied the sequence described by Stejskal and Tanner in 1965, which is a gradient spin echo T2-weighted sequence with two extra gradient pulses or motion-probing gradients (MPGs) that are equal in magnitude but have an opposite effect in direction 1) 7)-9). These MPGs were applied along the same directional axis before and after the 180° RF pulse 10). The first of the two gradient pulses acts as a dephasing gradient and the second pulse acts as a rephasing gradient 4) 9)-11). As the result, in a voxel of tissue with stationary water molecules or stable spins (Figure-1A), the opposite effect of the gradient pulses will cancel each other and strong signal intensity will be produced. However, moving water molecules or moving spins (Figure-1B) will be affected by the opposite effect of gradient pulses in the different locations, so the two gradients no longer have the same magnitude and they cannot cancel out each other, and thus little signal intensity is produced 6) 7) 9) 11). The use of dMRI in clinical application was
286
established in 1990 after Michael Moseley discovered that the water diffusion coefficient drops approximately 30-50% in an acute cat brain ischemia model, within several minutes after the occlusion of the middle cerebral artery 1) 2). Shortly after that, Robert Turner applied echo-planar imaging (EPI) (Figure-2A), with the result that diffusion images can be obtained in a short enough time to be feasible for clinical application; the motion artifacts were reduced, and sensitivity to signal changes due to the molecular motion was increased 1) 2) 11) 12). DWI (Figure-2B) provides an image with a unique contrast, and the diffusion coefficient could be quantitatively evaluated with apparent diffusion coefficients (ADC) 9) 13). Moreover, an ADC map (Figure-2C) can be created and ADC values can be displayed by drawing a region of interest (ROI) 9). The signal intensity (SI) of each voxel of tissue in diffusion weighted image is calculated with the following equation 9): SI = SI0×exp (-b×ADC),
Juntendo Medical Journal 63 (4), 2017
A
B
C
D
E
F
Figure-2
Images of diffusion magnetic resonance imaging (dMRI)
A. T2-weighted image (WI)-echo planar imaging (EPI). B. Isotropic diffusion-weighted (DW) image. C. Apparent diffusion coefficient (ADC) map. D. Fractional anisotropy (FA) map. E. FA color map. F. Tractography 7). (Kore de wakaru Kakusan MRI, 3rd Edition, 2013: 16-55 7))
where SI0 is the signal intensity on the T2-weighted (or b = 0 sec/mm 2 or b0 image), b is a diffusion sensitivity factor that consists of γ 2 G 2 δ 2 (Δ-δ/3). Here, G is the strength of the gradients, Δ (large delta) is the time interval between the gradients, δ (small delta) is the duration of the gradient, Δ-δ/3 is the diffusion time, and γ is the gyromagnetic constant 7) 9) 10). Higher ADC values that reflect faster water diffusion appear as lower signal intensity on DW images, while lower ADC values reflect restricted diffusion that is represented as high signal intensity on DW images 14). It is also important that the contrast of DW images is not only due to the differences in ADC but is also influenced by the T2-weighted contrast. Increased signal intensity on DW images that is caused by T2 prolongation is termed the“T2 shine through effect”7) 9) 14). Clinical applications of isotropic diffusion MRI Below are some examples of the utilization of DWI in clinical practice:
1. Acute brain infarction A rapid disruption of cerebral blood flow will interrupt the energy metabolism and ion exchange pumps and cause cytotoxic edema. This pathology will appear as a high signal intensity on DWI with low ADC value (Figure-3A) 9) 15 ). 2. Hemorrhage The imaging characteristics of hemorrhage lesions on dMRI are dependent on the age of the hemorrhage, which contains different blood products. During the hyperacute stage, intracellular oxyhemoglobin will restrict the water movement inside the red blood cells and cause high signal intensity on DWI and lower ADC (Figure-3B) 9) 15). 3. Abscess The high viscosity of pus causes restricted diffusion that is marked with a strong hyperintensity on DWI and reduced ADC 9) 16). The characteristics of an abscess on DWI can be used to differentiate the necrotic portion of a tumor which has a lower signal intensity and higher ADC ratios (Figure-3C) 9).
287
Andica, et al: Diffusion MRI: isotropic DWI and beyond
B
C
D
E
F
ADC
DWI
ADC
DWI
A
Figure-3
Example cases that show the utility of isotropic dMRI
The first row is the diffusion-weighted images (DWI) and the second row is the apparent diffusion coefficient (ADC) maps, all the brain pathologies are indicated with red arrows. A. Acute brain infarction, B. Hemorrhage, C. Abscess, D. Herpes encephalitis, and E. Creutzfeldt-Jakob disease (CJD) show high signal intensity on DWI and low signal intensity on the ADC map, which reflects the restricted diffusion. F. Epidermoid cyst has high signal intensity on DWI but is isointense with the gray matter on the ADC map; this might be caused by the“T2-shine through”effect.
288
Juntendo Medical Journal 63 (4), 2017
4. Herpes encephalitis Herpes encephalitis lesions are shown as hyperintense lesions on DWI and hypointense on ADC map that caused by the cytotoxic edema and high cellularity, that reflect irreversible neuronal damage undergoing necrosis (Figure-3D) 9) 15) 17). 5. Creutzfeldt-Jakob disease (CJD) DWI is proposed as the most sensitive imaging technique in the diagnosis of CJD 18). It will show hyperintense lesions with low ADC value in the cortex and basal ganglia (putamen and caudate nucleus) or thalamus (Figure-3E) 9) 15) 18)-20). Histopathologically, the hyperintensity is correlated with intraneuronal vacuolation, astrocytic gliosis, and the deposition of prion (PrP Sc) 19). 6. Epidermoid tumor DWI is an important imaging technique to differentiate between epidermoid and arachnoid cysts. An epidermoid cyst appears hyperintense on DWI, while an arachnoid cyst will not. Epidermoid cysts have almost similar ADC value as the gray and white matter, suggesting that the high signal intensity in DWI is caused by the“T2 shinethrough”effect (Figure-3F) 9) 13) 15). 7. Multiple sclerosis (MS) On DWI, the signal intensity of MS is variable. MS
A
active lesions that have enhanced contrast usually appear with hyperintensity relative to white matter, while chronic lesions were isointense. Some studies showed increasing ADC values in MS lesions and normal-appearing white matter. The “T2 shine-through”effect is predicted as the cause of this increased signal intensity 15). 8. Other brain lesions that show restricted diffusion Brain tumors with high cellularities, such as glioma, malignant meningioma, and lymphoma might appear as high signal intensity in DWI with low ADC values. However, in many cases these brain tumors might also appear iso- or hypointense 15) 21). Metabolic brain diseases with white matter abnormalities or leukodystrophies, such as metachromatic leukodystrophy, Krabbe disease, Canavan disease, and phenylketonuria also appear as high signal intensity on DW images with low ADC values 22). The development of dMRI 1. Diffusion tensor imaging (DTI) During the development of DTI, Michael Moseley found that water diffusion in a parallel direction to the white matter fibers is faster than in the perpendicular direction 1). This anisotropic phenomenon of
B
Figure-4 The anisotropic property of the diffusion of the brain
C
D
A. The corpus callosum (arrow) is hypointense when the motion-probing gradient (MPG) is applied in the right-to-left direction. B. The peritrigonal white matter tracts, including the optic radiation tract and occipito-thalamic tract (arrow), are hypointense when the MPG is applied in the anterior-to-posterior direction. C. The corticospinal tracts (arrow) are hypointense when the MPG applied in the superiorto-inferior direction. D. The DWI image is the result of the cube root of the multiplication of images A, B, and C. In this image, an acute infarct in the left basal ganglia (arrow) appears as hyperintense 7). (Kore de wakaru Kakusan MRI, 3rd Edition, 2013: 16-55 7))
289
Andica, et al: Diffusion MRI: isotropic DWI and beyond
the diffusion in the brain can be seen by comparing images that are acquired in three different gradient directions 9). Structures that have the same direction as the gradient appear as low signal intensity on the b0 image 9). In the x- (right-to-left) direction, the corpus callosum is hypointense (Figure-4A), in the y- (anterior-to-posterior) direction, the peritrigonal white matter tracts, including the optic radiation tract and occipito-thalamic tract are hypointense (Figure-4B), and in the z- (superior-to-inferior) direction, the corticospinal tracts are hypointense (Figure-4C). The signal intensities on DWI are the cube root of the multiplication of those three images with different gradient directions (Figure-4D) 7) 9). To deal with the anisotropic diffusion, Peter Basser introduced diffusion tensor imaging (DTI) as a technique to determine the true directionality of diffusion 1) 9) 23). Since then, DTI has been widely used for research and clinical applications as it can show the white matter pathologies 4). DTI is an imaging technique that is represented with a 3D ellipsoid shape model (Figure-5) with eigenvalues: λ1 (axial diffusivity) the longest axes, reflects the diffusion that parallel to the white matter fibers; λ2 and λ3, the middle and the shortest axes reflects the diffusion that is perpendicular to white matter fibers. The radial diffusivity is the average value of λ2 and λ3 7) 24). Moreover, eigenvectors; V1, V2, and V3 represent the direction of the diffusion 7) 24). The interpretations of DTI parameters are: 1) fractional anisotropy (FA) that describe the degree of anisotropy of a diffusion process, with values from 0 (isotropic) to 1 (anisotropic) 24); 2) mean diffusivity (MD); 3) radial diffusivity (RD);
V1
and 4) axial diffusivity (AD) 4). FA reflects the axonal density, axonal integrity, and fiber tract complexity 25). Reduced FA was found in many diseases, that is attributed to degradation of myelin sheaths and axonal membranes, or reduced density of axonal fibers 4) 25). MD, which is the mean-squared displacement of water molecules restricted by organelles and membranes, reflects cellular density and extracellular volume and relates to the volume fraction of the interstitial space 25). Moreover, RD will increase in demyelination, edema will increase the MD, and cell proliferation will decrease the MD 4). Below are some developments of DTI: 1) White matter tractography The major eigenvector of DT is assumed to be parallel to the white matter fibers and can be simply visualized using color maps (Figure-2E) 4). The connectivity of white matter on the basis of anisotropic diffusion of water can also be visualized using tractography (Figure-2F) 4). To delineate the white matter tract with tractography, first we have to draw the first region of interest (ROI) in the “seed”point, followed by the second ROI in the “target”point 26). From the seed, the lines were sequentially along the major eigenvector of the diffusion ellipsoid and only the lines that reached the target will be displayed 26). Until now, tractography has been used to reconstruct many white matter pathways, such as the corticospinal tract and corona radiata 4). 2) Tract-based spatial statistics (TBSS) TBSS is a voxelwise statistical analysis, that was originally used to evaluate FA values in the white matter 23). The aims of this method are: 1) to provide an appropriate registration algorithm, and
V3
V2
λ1
Axial diffusivity/D//
λ2
Radial Diffusivity/D⊥ =(λ2+λ3/2)
λ3
Figure-5
Ellipsoid model of diffusion tensor imaging (DTI)
In this model, the size and the direction of the diffusion are represented by the eigenvalues (λ) and eigenvectors (V), respectively. λ1 (the axial diffusivity) with the longest axes, reflects the diffusion that parallel to the white matter fibers, whereas, λ2 and λ3 reflects the diffusion that is perpendicular to the white matter fibers.
290
Juntendo Medical Journal 63 (4), 2017
A
Figure-6
B
C
Tract-based spatial statistics (TBSS) is a voxelwise statistical analysis of the FA values and other diffusion parameters in the white matter
The steps involved in TBSS are: Non-linear registration of native FA image (A) to normalized image (B) and create a mean FA skeleton (C). The mean FA skeleton is representative of the centers of all tracts 7). (Kore de wakaru Kakusan MRI, 3rd Edition, 2013: 16-55 7))
2) to improve the sensitivity, objectivity, and interpretability of multi-subject analysis studies 27). Firstly, FA images were co-registered to each other using nonlinear registration (Figure-6A, B). Next, the group mean FA image was created and thinned to create a mean FA skeleton that represents the centers of all tracts common to the group (Figure-6C). Each subjectʼs FA images (that can now also be used for other DTI parameters) was then projected into the skeleton and the resulting data fed into voxelwise cross-subject statistics 25)-28). As an example, Kamagata, et al. 28) applied TBSS to investigate the white-matter alteration in Parkinson disease with dementia (PDD) and as the results, they found lower FA values in many major tracts in PDD patients. Recent studies also showed that TBSS can be used to analyze the parameters of more advanced DWI techniques such as diffusion kurtosis imaging (DKI) 28) and neurite orientation dispersion and density imaging (NODDI) 29), that will be discussed in the next section. 2. Diffusion kurtosis imaging (DKI) DWI and DTI are in vivo imaging modalities that are based on the Gaussian distribution of water molecules 8) 23); however, biological tissue has a complex microstructure, such that water diffusion is hindered and the distribution of diffusion does not follow a Gaussian distribution or called nonGaussian 23) 30). For that reason, diffusion kurtosis imaging (DKI) was introduced to image the non-Gaussian characteristics of water diffusion and
reflect the heterogeneity of tissue 23). Moreover, DKI is sensitive to the alterations in the area of crossing fibers, which FA was not able to measure 30). The indices of DKI are as follows: 1) mean kurtosis (MK), the average of the diffusion kurtosis along all diffusion directions; 2) axial kurtosis (AK), the kurtosis along the axial direction of the directions of the axons; and 3) radial kurtosis (RK), the kurtosis along the radial direction of the directions of the axons 23). MK reflects the complexity of the tissue, the axonal water fraction, and the tortuosity of the extracellular space 28). Kamiya, et al. 30) showed that MK had a significant correlation with clinical severity and provide more information compared to FA in normal pressure hydrocephalus patients. 3. Neurite orientation dispersion and density imaging (NODDI) NODDI is an imaging method based on dMRI, and was proposed to overcome the limitations of DTI 29) 31) 32). NODDI is able to quantify the direction and the structure of neurites (axons and dendrites) by the orientation-dispersed cylinder model and the Watson distribution 29) 31) 32). Regarding NODDI parameters, the intracellular volume fraction (Vic) represents the neurite density, and the orientation dispersion index (OD) represents the OD of the neurite 32). The indices of NODDI are assumed to describe the changes in FA 31) 32). Kamagata, et al. 31) compared DTI and NODDI parameters in diagnosing PD and found that Vic in the substantia nigra pars compacta is the best index for diagnosing PD. Reduced Vic reflects the neuron loss which is the
291
Andica, et al: Diffusion MRI: isotropic DWI and beyond
histopathological hallmark of PD 31). References 1) Le Bihan D: Diffusion MRI: what water tells us about the brain. EMBO Mol Med, 2014; 6: 569-573. 2) Le Bihan D, Breton E, Lallemand D, Aubin ML, Vignaud J, Laval-Jeantet M: Separation of diffusion and perfusion in intravoxel incoherent motion MR imaging. Radiology, 1988; 168: 497-505. 3) Le Bihan D, Iima M: Diffusion magnetic resonance imaging: what water tells us about biological tissues. PLoS Biol, 2015; 13: e1002203. 4) Alexander AL, Lee JE, Lazar M, Field AS: Diffusion tensor imaging of the brain. Neurotherapeutics, 2007; 4: 316-329. 5) Aoki S, Kamagata K, Hori M, et al: Advances diffusional MR imaging of the brain. Juntendo Medical Journal, 2016; 62: 300-301. 6) Hagmann P, Jonasson L, Maeder P, Thiran JP, Wedeen VJ, Meuli R: Understanding diffusion MR imaging techniques: from scalar diffusion-weighted imaging to diffusion tensor imaging and beyond. Radiographics, 2006; 26 Suppl 1: S205-223. 7) Kore de wakaru Kakusan MRI, 3rd Edition. In: Aoki S, Abe O, Masutani Y, Takahara T, eds. Tokyo: Gakken Medical Shujunsha Co., Ltd. 2013: 16-55. (in Japanese) 8) Le Bihan D, Johansen-Berg H: Diffusion MRI at 25: exploring brain tissue structure and function. Neuroimage, 2012; 61: 324-341. 9) Schaefer PW, Grant PE, Gonzalez RG: Diffusionweighted MR imaging of the brain. Radiology, 2000; 217: 331-345. 10) Nucifora PG, Verma R, Lee SK, Melhem ER: Diffusiontensor MR imaging and tractography: exploring brain microstructure and connectivity. Radiology, 2007; 245: 367-384. 11) Chilla GS, Tan CH, Xu C, Poh CL: Diffusion weighted magnetic resonance imaging and its recent trend-a survey. Quant Imaging Med Surg, 2015; 5: 407-422. 12) Turner R, Le Bihan D, Chesnick AS: Echo-planar imaging of diffusion and perfusion. Magn Reson Med, 1991; 19: 247-253. 13) Chen S, Ikawa F, Kurisu K, Arita K, Takaba J, Kanou Y: Quantitative MR evaluation of intracranial epidermoid tumors by fast fluid-attenuated inversion recovery imaging and echo-planar diffusion-weighted imaging. AJNR Am J Neuroradiol, 2001; 22: 1089-1096. 14) Mukherjee P, Berman JI, Chung SW, Hess CP, Henry RG: Diffusion tensor MR imaging and fiber tractography: theoretic underpinnings. AJNR Am J Neuroradiol, 2008; 29: 632-641. 15) Stadnik TW, Demaerel P, Luypaert RR, et al: Imaging tutorial: differential diagnosis of bright lesions on diffusion-weighted MR images. Radiographics, 2003; 23: e7. 16) Cartes-Zumelzu FW, Stavrou I, Castillo M, Eisenhuber E, Knosp E, Thurnher MM: Diffusion-weighted imaging in the assessment of brain abscesses therapy. AJNR Am J Neuroradiol, 2004; 25: 1310-1317.
292
17) Sener RN: Herpes simplex encephalitis: diffusion MR imaging findings. Comput Med Imaging Graph, 2001; 25: 391-397. 18) Kallenberg K, Schulz-Schaeffer WJ, Jastrow U, et al: Creutzfeldt-Jakob disease: comparative analysis of MR imaging sequences. AJNR Am J Neuroradiol, 2006; 27: 1459-1462. 19) Geschwind MD, Potter CA, Sattavat M, et al: Correlating DWI MRI with pathologic and other features of Jakob-Creutzfeldt disease. Alzheimer Dis Assoc Disord, 2009; 23: 82-87. 20) Ukisu R, Kushihashi T, Tanaka E, et al: Diffusionweighted MR imaging of early-stage Creutzfeldt-Jakob disease: typical and atypical manifestations. Radiographics, 2006; 26 Suppl 1: S191-204. 21) Maier SE, Sun Y, Mulkern RV: Diffusion imaging of brain tumors. NMR Biomed, 2010; 23: 849-864. 22) Patay Z: Diffusion-weighted MR imaging in leukodystrophies. Eur Radiol, 2005; 15: 2284-2303. 23) Steven AJ, Zhuo J, Melhem ER: Diffusion kurtosis imaging: an emerging technique for evaluating the microstructural environment of the brain. AJR Am J Roentgenol, 2014; 202: W26-33. 24) Mori S, Zhang J: Principles of diffusion tensor imaging and its applications to basic neuroscience research. Neuron, 2006; 51: 527-539. 25) Yeh PH, Simpson K, Durazzo TC, Gazdzinski S, Meyerhoff DJ: Tract-Based Spatial Statistics (TBSS) of diffusion tensor imaging data in alcohol dependence: abnormalities of the motivational neurocircuitry. Psychiatry Res, 2009; 173: 22-30. 26) Kunimatsu A, Aoki S, Masutani Y, et al: The optimal trackability threshold of fractional anisotropy for diffusion tensor tractography of the corticospinal tract. Magn Reson Med Sci, 2004; 3: 11-17. 27) Smith SM, Jenkinson M, Johansen-Berg H, et al: Tractbased spatial statistics: voxelwise analysis of multisubject diffusion data. Neuroimage, 2006; 31: 1487-1505. 28) Kamagata K, Motoi Y, Tomiyama H, et al: Relationship between cognitive impairment and white-matter alteration in Parkinsonʼs disease with dementia: tract-based spatial statistics and tract-specific analysis. Eur Radiol, 2013; 23: 1946-1955. 29) Timmers I, Roebroeck A, Bastiani M, Jansma B, RubioGozalbo E, Zhang H: Assessing microstructural substrates of white matter abnormalities: a comparative study using DTI and NODDI. PLoS One, 2016; 11: e0167884. 30) Kamiya K, Kamagata K, Miyajima M, et al: Diffusional kurtosis imaging in idiopathic normal pressure hydrocephalus: correlation with severity of cognitive impairment. Magn Reson Med Sci, 2016; 15: 316-323. 31) Kamagata K, Hatano T, Okuzumi A, et al: Neurite orientation dispersion and density imaging in the substantia nigra in idiopathic Parkinson disease. Eur Radiol, 2016; 26: 2567-2577. 32) Zhang H, Schneider T, Wheeler-Kingshott CA, Alexander DC: NODDI: practical in vivo neurite orientation dispersion and density imaging of the human brain. Neuroimage, 2012; 61: 1000-1016.