- 1Department of Radiology, Juntendo University School of Medicine, Tokyo, Japan
- 2Department of Radiology, Toho University, Tokyo, Japan
- 3Department of Neurology, Juntendo University School of Medicine, Tokyo, Japan
- 4Siemens Healthcare K.K., Tokyo, Japan
- 5Siemens Healthcare GmbH, Erlangen, Germany
Microstructure imaging by means of multidimensional diffusion encoding is increasingly applied in clinical research, with expectations that it yields a parameter that better correlates with clinical disability than current methods based on single diffusion encoding. Under the assumption that diffusion within a voxel can be well described by a collection of diffusion tensors, several parameters of this diffusion tensor distribution can be derived, including mean size, variance of sizes, orientational dispersion, and microscopic anisotropy. The information provided by multidimensional diffusion encoding also enables us to decompose the sources of the conventional fractional anisotropy and mean kurtosis. In this study, we explored the utility of the diffusion tensor distribution approach for characterizing white-matter degeneration in aging and in Parkinson disease by using double diffusion encoding. Data from 23 healthy older subjects and 27 patients with Parkinson disease were analyzed. Advanced age was associated with greater mean size and size variances, as well as smaller microscopic anisotropy. By analyzing the parameters underlying diffusion kurtosis, we found that the reductions of kurtosis in aging and Parkinson disease reported in the literature are likely driven by the reduction in microscopic anisotropy. Furthermore, microscopic anisotropy correlated with the severity of motor impairment in the patients with Parkinson disease. The present results support the use of multidimensional diffusion encoding in clinical studies and are encouraging for its future clinical implementation.
Introduction
Parkinson disease (PD) is a neurodegenerative disorder characterized by motor symptoms (akinesia, resting tremor, and rigidity) and a wide range of cognitive, neuropsychiatric, and autonomic dysfunctions (Poewe et al., 2017). Advanced age is a major risk factor for the development of PD and is also associated with faster motor decline (Levy, 2007; Collier et al., 2017). Pathologically, PD is characterized by widespread aggregation of α-synuclein-immunoreactive inclusions in the form of Lewy pathology within both the neuronal cytoplasm (Lewy bodies) and axons (Lewy neurites) (Braak et al., 2003; Kanazawa et al., 2012; Poewe et al., 2017). Neuropathological studies have indicated that Lewy pathologies evolve along major fiber pathways, beginning in the brain stem and eventually advancing to the neocortical regions (Braak et al., 2003). Accumulating evidence has suggested that axonal degeneration is an early event in the process of neurodegeneration that is common to PD and other age-related neurological diseases (Kurowska et al., 2016; Salvadores et al., 2017). Non-invasive characterization of the neurodegeneration underlying the pathogenesis and progression of PD is of high clinical demand, because it will aid in the development of novel therapeutic strategies and in monitoring the effects of treatment.
Diffusion MRI (dMRI) is uniquely sensitive to tissue features on the micrometer scale and therefore has been widely used to study neurodegeneration in aging (Madden et al., 2012; Coutu et al., 2014; Billiet et al., 2015; Benitez et al., 2018; Guerreri et al., 2019) and diseases (Goveas et al., 2015; Andica et al., 2019b). In PD, the majority of studies have applied diffusion tensor imaging (DTI) (Basser, 1995) and have consistently reported smaller fractional anisotropy (FA) and greater mean diffusivity (MD) in the white matter of patients than of controls (for a recent review and meta-analysis see Atkinson-Clement et al., 2017). Several groups (Wang et al., 2011; Kamagata et al., 2013, 2014; Surova et al., 2016, 2018) further explored the utility of dMRI acquisition with higher b-values than used in DTI, which is typically analyzed by means of diffusion kurtosis imaging (DKI) (Jensen et al., 2005). The works by Kamagata et al. (2013, 2014) showed that kurtosis in the white matter is reduced in patients and that DKI is more sensitive to white-matter degeneration than is DTI. Of note, several studies demonstrated that white-matter degeneration as detected by dMRI precedes macroscopic gray matter atrophy (Agosta et al., 2013; Duncan et al., 2016; Rektor et al., 2018), suggesting the potential of dMRI parameters as early biomarkers of PD.
Regarding acquisition, both DTI and DKI use single diffusion encoding (SDE), which uses one pair of diffusion-sensitizing gradients. Multidimensional diffusion encoding (Mitra, 1995; Topgaard, 2017), that can be realized by using either double diffusion encoding (DDE) (Cory et al., 1990; Callaghan and Xia, 1991), triple diffusion encoding (Mori and Van Zijl, 1995), or continuous gradient waveforms (Caprihan et al., 1996; Callaghan, 1997; Eriksson et al., 2013), has recently gained attention in clinical studies because these methods provide more specific information about the tissue microstructure. For example, cumulant expansion of the DDE signal up to the fourth-order term of the gradient amplitude (Jespersen, 2012) shows DDE provides unique information that is not contained in the standard diffusion and kurtosis tensors available with SDE. In addition to DTI/DKI parameters, at least two new properties of biological interest can be obtained: the microscopic anisotropy (Cory et al., 1990; Cheng and Cory, 1999; Callaghan and Komlosh, 2002; Özarslan and Basser, 2008; Shemesh et al., 2009; Lawrenz et al., 2010; Jespersen et al., 2013; Lasič et al., 2014) and the variance of isotropic diffusivities among individual microenvironments (compartments) (Szczepankiewicz et al., 2016; Westin et al., 2016). Several works have demonstrated that the microscopic anisotropy, a measure of diffusion anisotropy that is not confounded by orientational dispersion of axons or fibers, can be estimated reliably in nervous tissues (Komlosh et al., 2007, 2008; Shemesh and Cohen, 2011; Shemesh et al., 2012; Jespersen et al., 2013; Lawrenz and Finsterbusch, 2013; Lawrenz et al., 2015). Microscopic anisotropy has been shown to be useful for characterizing white-matter degeneration in aging (Lawrenz et al., 2016) and in multiple sclerosis (Yang et al., 2018; Andersen et al., 2020). The variance of isotropic diffusivities has also been suggested to be promising for deducing the microstructural underpinnings of the diffusion changes in diseases like schizophrenia (Westin et al., 2016). Moreover, multidimensional diffusion encoding enables us to decompose the sources of DTI/DKI parameters (Lasič et al., 2014; Szczepankiewicz et al., 2016; Westin et al., 2016; Henriques et al., 2020). Such analyses are expected to improve our understanding of neurodegeneration and the sensitivity of imaging to pathology. For example, if the pathology affects two underlying sources of a particular DTI/DKI parameter in opposite directions, as has been suggested in Douaud et al. (2011); Lawrenz et al. (2016), Lampinen et al. (2019), the sensitivity and interpretability of that parameter would be limited, whereas separation of each source may provide useful information.
In this study, we analyzed DDE data by using the covariance tensor framework (Westin et al., 2016) that works under the diffusion tensor distribution (DTD) model (Jian et al., 2007). Given the DTD model, continuous waveforms would offer more efficient diffusion weighting and faster acquisition than DDE (Nilsson et al., 2020). We adopted DDE here because this study was planned before the definition of diffusion time for the continuous waveform was made by Lundell et al. (2019). As DTD is a model for the long diffusion time limit where time-dependence is negligible (Novikov et al., 2019), we needed to compare the diffusion time of our measurement with the previous studies (Clark et al., 2001; Portnoy et al., 2013) to rationalize our assumption. The potential limitation arising from neglecting the effects of diffusion time is detailed in “Discussion.”
The aim of this study was to identify the sources of DTI/DKI changes previously reported in aging and PD in terms of the DTD model parameters. We first analyzed correlation with age in healthy older subjects and then examined group differences between the healthy subjects and patients with PD. Finally, we investigated the correlations with the severity of motor impairment in the patients.
Materials and Methods
Theory
Model and Assumptions
Linking the MRI signal to specific tissue properties involves modeling, simplification of the complex reality relying on a few assumptions. By adopting the DTD model, we assume that the voxel consists of multiple, non-exchanging Gaussian compartments. More precisely, we model the diffusion within a voxel by a distribution of diffusion tensors (free, anisotropic diffusion) (Jian et al., 2007; Topgaard, 2017). In such case, the measurement can be fully characterized by the b-tensor (i.e., b-tensor encoding). In this regard, SDE is linear tensor encoding (LTE), whereas DDE provides planar tensor encoding (PTE) if we apply two diffusion encodings in non-colinear directions (Westin et al., 2016; Topgaard, 2017). Although the DTD model is certainly a crude assumption and not fully validated yet in human brain, describing the tissue by a combination of several diffusion tensors is a common starting point in most dMRI models currently used in clinical studies (Jelescu and Budde, 2017).
Under the DTD model, the signal can be written as:
where S(B) is the signal, S0 is the signal without diffusion weighting, B is the b-tensor, and Dc is the diffusion tensor of each compartment. The bracket < > denotes the ensemble average over the voxel. The colon denotes a generalized scalar product between the two tensors, B:D = ∑i∑jBijDij.
Diffusion Parameters
Multidimensional diffusion encoding provides us means to extract summarizing statistics of the distribution of diffusion tensors. Eq. 1 can be expanded (Westin et al., 2016; Topgaard, 2017) to:
where ⊗ denotes a tensor outer product. Here, D is the well-known, voxel-averaged diffusion tensor, and C is a fourth-order tensor called the covariance tensor. From D and C, several scalar parameters that summarize microstructural features, including DTI/DKI parameters, can be computed (Westin et al., 2016; Topgaard, 2017; Henriques et al., 2020). Below, we briefly describe the meanings of the parameters used in this study and their relations to each other. First, the DTI parameters MD and FA are defined (Basser, 1995; Westin et al., 2016; Henriques et al., 2019) as:
Vλ(D) denote the variance of eigenvalues of D, defined by using the eigenvalues λ1, λ2, and λ3 as . MD and FA can be understood as the size and shape of the voxel-averaged diffusion tensor. The definition of μFA has a similar form to that of FA (Westin et al., 2016; Henriques et al., 2019):
Unlike FA, μFA is not influenced by orientational dispersion and purely reflects the anisotropy (shape) of microstructural environments. In addition, a measure of orientational dispersion (the orientational order parameter, OP) is defined (Lasič et al., 2014; Westin et al., 2016) as:
OP equals 1 for perfectly aligned orientations and equals 0 for fully isotropic dispersion. FA is influenced by both μFA and OP (Lasič et al., 2014):
Furthermore, under the DTD assumption, mean kurtosis (MK), as defined in Hansen et al. (2013), can be decomposed into two kurtosis sources (Westin et al., 2016; Henriques et al., 2020):
where
are the isotropic and anisotropic kurtosis sources. Here, V(Dc) is the variance of isotropic diffusivities among compartments (variance of sizes). The last term in Eq. 8 is a factor related to orientational dispersion and can be expressed as (Westin et al., 2016; Henriques et al., 2020)
Using Eq. 11, Eq. 8 can be re-written as:
In this study, we report MD, FA, MK, Kiso, Kaniso, μFA, OP, and V(Dc). Although the information represented by some parameters overlaps (e.g., both μFA and Kaniso are measures of microscopic anisotropy), having multiple forms of expression helps us to understand how the DTD parameters affect FA and MK (Eqs 7 and 12).
Participants
This study was carried out in accordance with the Declaration of Helsinki for experiments involving humans. The Institutional Review Board approved this study, and all subjects gave written informed consent prior to participation. Twenty-three healthy older subjects (63.1 ± 7.2 years old) and 27 patients with PD (66.1 ± 6.6 years old) were enrolled. Patients with PD were diagnosed by neurologists on the basis of clinical diagnostic criteria of the Movement Disorder Society (Postuma et al., 2015). Motor function and disease stage were evaluated with the Unified Parkinson’s Disease Rating Scale motor part (UPDRS-III) (Goetz et al., 2008) and the Hoehn and Yahr staging scale (Hoehn and Yahr, 1967). White-matter T2 hyperintensities (WMH) were rated according to the Fazekas scale (Fazekas et al., 1987). The clinical and demographic characteristics of the participants are summarized in Table 1.
Because DDE is a relatively new technique for clinical research, we also examined the stability of the dMRI parameters in terms of scan–rescan repeatability in a separate group of 4 young healthy subjects (3 male and 1 female, 25–37 years old), who underwent two scans on different days within a week. In these subjects, we also acquired conventional SDE for comparison.
Image Acquisition
The subjects were scanned by using a clinical 3T scanner (MAGNETOM Prisma, Siemens Healthcare, Erlangen, Germany) equipped with a 64-channel head coil using a prototype sequence. DDE data were acquired using a monopolar spin-echo type acquisition. We used parallel (LTE) and perpendicular (PTE) pairs of diffusion sensitizing gradient blocks, with 30 uniformly distributed directions. Note that, under the DTD assumption, the direction of PTE can be characterized by the normal vector. For both LTE and PTE, we used two shells, which had b-values of 1000 and 2000 s/mm2. One volume without diffusion weighting was also obtained. Thus, a total of 121 volumes were acquired. Other settings were image resolution = 3 × 3 × 3 mm3, TE = 93 ms, TR = 6300 ms, Δ = 18.3 ms, δ = 16.8 ms, mixing time = 24.9 ms, in-plane GRAPPA with acceleration factor 2, through-plane GRAPPA (SMS) with multiband factor 2, partial Fourier 6/8, 44 axial slices, and scan time 14 min. In addition, anatomical T1- and T2-weighted images were acquired and inspected by a neuroradiologist to rate Fazekas grade and to check for any other co-existing pathologies. SDE acquisition in the young subjects used the same b-values and 30 directions, resulting in a total of 61 volumes. The acquisition settings were identical to those of DDE, except Δ = 43.3 ms and δ = 36.4 ms. The scan time for SDE was 7 min.
Image Processing
Images were pre-processed by using FSL 6.0.1 (Jenkinson et al., 2012) and MRTrix3 (Tournier et al., 2019). Raw images were denoised (Veraart et al., 2016) and corrected for Gibbs artifact (Kellner et al., 2016), eddy currents and motion (Andersson and Sotiropoulos, 2016), and B1 inhomogeneity (Tustison et al., 2010). For eddy-current correction of DDE, the second diffusion direction was used as the input to eddy in FSL (Yang et al., 2018). The dMRI parameter maps were computed using the multidimensional diffusion MRI toolbox (Nilsson et al., 2018). In a small number of voxels where the fitting resulted in negative values of V(Dc), the fitting was repeated by using data smoothed with an isotropic 3D Gaussian kernel with sigma = 0.7 × voxel size.
Stability of Measurements
Using the scan-rescan data from the 4 young subjects, we computed the within-subject coefficient of variation (CVws), defined as . Here, μ is the grand mean and σws is the within-subject standard deviation. To compute CVws across the white-matter voxels, each subject’s FA map was non-linearly registered into the Johns Hopkins University (JHU) template in FSL. The other dMRI parameter maps were also transferred into the standard space by using the same deformation. A white-matter mask was generated by thresholding the FA template at FA > 0.2. To mitigate partial volume effects with the gray matter and cerebrospinal fluid (CSF), the mask was further eroded by one voxel. We also report CVws for the mean value of each white-matter region of interest (ROI) that was used for the analyses of aging and PD (Section “Statistical Analysis”).
Statistical Analysis
We examined the effects of the subjects’ characteristics (age, diagnosis, and UPDRS-III score) on the dMRI parameters by means of whole-brain voxel-wise analyses using tract-based spatial statistics (TBSS) (Smith et al., 2006) and atlas-based ROI analyses. First, all subjects’ FA images were aligned into a common space by means of non-linear registration, followed by creation of a mean FA skeleton. The threshold for creating the FA skeleton was set at FA > 0.2. Then, the aligned FA map of each subject was projected onto the FA skeleton. The other parameter maps were projected onto the mean FA skeleton by using the same transformation. Subsequently, the TBSS-processed skeletons were subjected to ROI analyses. Twelve white-matter ROIs were defined as the intersection of the white-matter skeleton and the JHU and Harvard–Oxford atlases in FSL as in Guerreri et al. (2019) (Figure 1), and the mean value within each ROI was extracted. Specifically, the deep white-matter ROIs (the anterior corona radiata [ACR], posterior corona radiata [PCR], superior corona radiata [SCR], anterior and posterior limb of the internal capsule [ALIC, PLIC], and genu and splenium of the corpus callosum [GCC, SCC]) were defined by using the JHU atlas, and the subcortical white-matter ROIs (frontal, sensory-motor, parietal, occipital, and temporal cortices) were defined by using the Harvard–Oxford atlas.
Figure 1. Atlas-based ROIs. Twelve ROIs were defined as the intersections between the TBSS-processed white-matter skeleton and the JHU and Harvard–Oxford atlases in FSL (ACR/SCR/PCR, anterior/superior/posterior corona radiata; ALIC/PLIC, anterior/posterior limb of the internal capsule; GCC/SCC, genu/splenium of the corpus callosum).
We opted not to exclude voxels with WMH from the analyses for the following two reasons. First, evidence suggests associations between WMH and the severity of motor and non-motor symptoms in PD (McDonald et al., 2016; Veselý et al., 2016). The pathophysiology of WMH is complex, and whether WMH are present in the patients because of incidental cerebrovascular disease or are a consequence of PD remains unclear (McDonald et al., 2016). Second, although TBSS allows the exclusion of voxels by using subject-specific lesion masks, this process causes differences in the degrees of freedom across voxels and may complicate the interpretation of statistical results (Benitez et al., 2018).
For both the TBSS and ROI-based analyses, we fit linear models where the dependent variable was each dMRI parameter. First, we examined the effect of age in the healthy subjects, using sex as a nuisance covariate. Then, we examined the differences between the healthy subjects and the patients with PD, using age and sex as nuisance covariates. Finally, we examined the correlation with UPDRS-III score in the patients, using age and sex as nuisance covariates. Because motor impairment might be correlated with age, we checked for multi-collinearity using the variance inflating factor (VIF) (Hair et al., 2010). Although we used a stringent criterion of VIF < 4.0 (Hair et al., 2010), VIF was below this threshold for all predictors in the linear models. Statistical inference of the linear models was conducted by using permutation analysis (Winkler et al., 2014). For computational efficiency, we used 2000 permutation and tail approximations (Winkler et al., 2016). For TBSS, threshold-free cluster enhancement (Smith and Nichols, 2009) was applied. The significance threshold was set at P < 0.05, corrected for family-wise error (FWE) for multiple comparisons across voxels or ROIs. For the ROI-based analyses, we also report results without FWE correction (Puncorrected < 0.05) because these data may be informative for future studies with larger samples.
Results
Scan-Rescan Repeatability
MD, FA, and MK showed excellent repeatability, with CVws below 6–8% for most of the white-matter voxels (Figure 2). The values of CVws were similar between DDE and SDE. μFA and OP also showed CVws below 8% for most of the voxels. The parameters of size variance [Kiso and V(Dc)] showed greater CVws, falling in the range of 10–25% for most of the voxels. As for ROI-based measurements, CVws was below 2% in all ROIs for MD, FA, μFA, and OP (Table 2). MK and Kaniso exhibited slightly greater values of CVws ranging from 0.5 to 4.8%. Although Kiso and V(Dc) showed CVws below 6% for most of the ROIs, poorer repeatability was found in the internal capsule and corpus callosum with CVws around 10%, reflecting the challenges of dMRI in these regions due to partial volume effects with the CSF and the gray matter, Gibbs artifact, and CSF pulsation.
Figure 2. Scan–rescan repeatability of the dMRI parameters derived from DDE. Top: the diffusion parameter maps from a representative subject. Middle: scatterplots comparing the first and the second sessions. Values from the white-matter voxels of all 4 subjects are shown. Warmer color represents higher density of that point. Bottom: within-subject coefficient of variation (CVws) in the white-matter voxels. The blue and orange histograms represent CVws of the dMRI parameters derived from DDE and SDE, respectively.
Correlation With Age in the Healthy Subjects
To display the results of ROI-based correlations, we adopted a figure format used in Billiet et al. (2015) and Guerreri et al. (2019) (Figure 3). For MD and V(Dc), positive correlation with age was seen in extensive regions of the white matter in both ROI-based analyses and TBSS (Figures 3–5). FA, MK, Kaniso, and μFA showed negative correlation with age within these regions. OP demonstrated a unique behavior in that both positive and negative correlations were observed, depending on anatomical locations. In particular, although the correlation was negative in most of the white-matter regions, positive correlation was observed in the ALIC and PLIC, external capsule, and SCR. Kiso showed positive correlation with age, although in fewer regions compared with the other diffusion parameters.
Figure 3. Correlation with age in the healthy control subjects. Pearson’s correlation coefficients are shown where a correlation between a diffusion parameter and age was found (Puncorrected < 0.05). Red denotes positive correlation; blue indicates negative correlation. Asterisks indicate significant correlations after FWE correction. ACR/SCR/PCR, anterior/superior/posterior corona radiata; ALIC/PLIC, anterior/posterior limb of the internal capsule; GCC/SCC, genu/splenium of the corpus callosum.
Figure 4. Correlation with age in two representative ROIs. Red lines represent linear fitting where a significant correlation was found (Puncorrected < 0.05). The shaded area represents the 95% confidence interval. ALIC, anterior limb of the internal capsule; SCC, splenium of the corpus callosum.
Figure 5. TBSS results for correlation with age in the healthy control subjects. Voxels with significant correlation after the FWE correction (P < 0.05) are shown. Hot colors denote positive correlation, whereas cold colors denote negative correlation. To aid visualization, the results were thickened by using the fill function in FSL.
Difference Between the Healthy Control and the Patients With PD
The ROI-based analyses revealed that the patients with PD had greater values of MD, Kiso, and V(Dc) in the subcortical white-matter regions than the controls (Figure 6), although only the differences in Kiso and V(Dc) of the parietal ROI remained significant after the FWE correction. The group differences were below the significance threshold for Kaniso and μFA for all ROIs. Slightly greater values of MK were observed in the patients in the PLIC and SCC. Smaller values of FA and OP were observed in the ALIC and PLIC, respectively. In TBSS, only OP revealed statistically significant results: the patients exhibited smaller values in the left internal capsule, the left cerebral peduncle, the left thalamus, and the right external capsule (Figure 7).
Figure 6. Differences between the healthy control subjects (light blue) and the patients with PD (dark blue) in the ROI-based analyses. The red asterisks denote statistically significant differences (*Puncorrected < 0.05, **P < 0.05). ACR/SCR/PCR, anterior/superior/posterior corona radiata; ALIC/PLIC, anterior/posterior limb of the internal capsule; GCC/SCC, genu/splenium of the corpus callosum.
Figure 7. TBSS results for group differences between the healthy control subjects and the patients with PD. Voxels with significant differences after the FWE correction (P < 0.05) are shown. Voxels where the patients exhibited smaller values are represented in cold colors. In no voxels did patients exhibit greater values. To aid visualization, the results were thickened by using the fill function in FSL.
Correlation With Motor Impairment in PD
The ROI-based analyses revealed positive correlations with UPDRS-III score for MD, Kiso, and V(Dc), and negative correlations for FA, MK, Kaniso, and μFA in GCC and SCC (Figures 8, 9). MK, Kaniso, and μFA exhibited negative correlations also in the subcortical white-matter ROIs. TBSS showed negative correlations for MK, Kaniso, and μFA in extensive areas in the frontal and parietal lobes and the corpus callosum, and positive correlation for V(Dc) in SCC (Figure 10). For MD, FA, and OP, the correlations were below the significance threshold in TBSS.
Figure 8. Correlation with UPDRS-III score in the patients with PD. The Pearson correlation coefficients shown indicate significant correlation between a diffusion parameter and UPDRS-III score (Puncorrected < 0.05). Red denotes positive correlations; blue indicates negative correlations. Asterisks indicate significant correlations after FWE correction. ACR/SCR/PCR, anterior/superior/posterior corona radiata; ALIC/PLIC, anterior/posterior limb of the internal capsule; GCC/SCC, genu/splenium of the corpus callosum.
Figure 9. Correlation with UPDRS-III score in two representative ROIs. Red lines represent linear fitting where significant correlation was found (Puncorrected < 0.05). The shaded area represents the 95% confidence interval. SCC, splenium of the corpus callosum.
Figure 10. TBSS results for correlation with UPDRS-III score in the patients with PD. Voxels with significant correlation after the FWE correction (P < 0.05) are shown. Hot colors denote positive correlation, whereas cold colors denote negative correlation. To aid visualization, the results were thickened by using the fill function in FSL.
Discussion
This study investigated the utility of DTD parameters derived from DDE to characterize white-matter degeneration in aging and PD. Advanced age was associated with greater MD and smaller FA, MK, Kaniso, and μFA, in agreement with previous studies that used SDE (Madden et al., 2012; Coutu et al., 2014; Billiet et al., 2015; Benitez et al., 2018; Guerreri et al., 2019) and DDE (Lawrenz et al., 2016). OP decreased with age for most of the voxels, which, together with the decrease of μFA, led to the age-related decrease of FA. At the same time, a positive correlation between OP and age was observed in some white matter regions including the ALIC, PLIC, and SCR. Our results of age-related changes in μFA and OP replicated the findings by the earlier study (Lawrenz et al., 2016). Age-related OP increase in some regions can be understood as selective degeneration of the secondary crossing fibers, thus increasing the proportion of aligned fibers. The work by Lawrenz et al. (2016) employed a different approach than the DTD model, i.e., microscopic anisotropy estimation proposed in Lawrenz et al. (2010) which is based on the signal equation of restricted diffusion provided by Mitra (1995). Though it is beyond the scope of this study to examine which model should be preferred, it seems that the observed trends in age-related changes of μFA and OP are robust against this difference in model assumptions. We also found an age-related increase in the measures of size variance [Kiso and V(Dc)]. The effects of aging and PD were in the same directions for most of the diffusion parameters and anatomical locations, such as increases of MD, Kiso, and V(Dc) and decreases of FA, Kaniso, and μFA. This is biologically plausible, because neurodegeneration in aging and PD have many features in common, including a declined ability to homeostatically regulate proteostasis, neuroinflammation, mitochondrial dysfunction, oxidative stress, degeneration of the myelin sheath, accumulation of cellular debris, alteration of axonal transport, axonal swelling/beading, and loss of axons (Tagliaferro et al., 2015; Collier et al., 2017; Salvadores et al., 2017; Calabrese et al., 2018; Datar et al., 2019). The present results also suggest that the measures of microscopic anisotropy (Kaniso and μFA) might be useful to track white-matter degeneration related to the motor impairment in PD.
Through the analysis of underlying kurtosis sources via the DTD model, we found that the reductions of MK reported in aging (Coutu et al., 2014; Benitez et al., 2018; Guerreri et al., 2019) and PD (Kamagata et al., 2013, 2014) are likely driven by the reduction of microscopic anisotropy. The decomposition of kurtosis sources (Eq. 12) suggests that either an increase of Kiso or a decrease of OP would lead to an increase of MK. We indeed observed slightly greater MK in the PLIC and SCC in the patients than in the controls. Although this appears to contradict previous reports of smaller MK in patients (Kamagata et al., 2013, 2014), it may instead indicate that the increase of Kiso and decrease of Kaniso have different time trajectories during disease progression. If the increase of Kiso is related to early neuroinflammation and the decrease of Kaniso represents the subsequent degeneration, as speculated by Andica et al. (2019a), MK may possibly have a non-monotonical trajectory, with an initial increase followed by decrease, assuming constant OP (Eq. 12). Future studies to investigate the trajectories of the DTD parameters in relation to disease progression are warranted. We also demonstrated that microscopic anisotropy can be estimated with excellent repeatability comparable to that of SDE-derived DTI/DKI parameters in both the present study and the literature (Shahim et al., 2017; Ades-Aron et al., 2018). Taken together, the present results are encouraging that multidimensional diffusion encoding could be implemented clinically. Although we used a relatively long acquisition for this exploratory study, μFA can be estimated reliably with a reduced DDE acquisition (Yang et al., 2018; Kerkelä et al., 2020) to make the scan time clinically feasible.
Our results support several hypothetical interpretations of previous SDE studies that relied on more restrictive assumptions. (Chad et al., 2018) and (Andica et al., 2019a) applied regularized fitting of the bi-tensor free-water DTI model (Pasternak et al., 2009) to study aging and PD, respectively. Their results showed that aging and PD are associated with an increase of the free-water fraction in the white matter. A larger free-water fraction would lead to greater size variance and smaller microscopic anisotropy (Szczepankiewicz et al., 2015; Westin et al., 2016), as observed in our study. Certainly, an increase in free water is not the only possible mechanism, and further investigation is required to elucidate specific pathological features. For example, patchy axonal loss and demyelination may lead to similar results via heterogenous increase of the radial diffusivities of compartmental diffusion tensors. Changes in axon morphology such as swelling and beading, which would reduce axial diffusivities of compartmental diffusion tensors (Budde and Frank, 2010; Palombo et al., 2018) and hence μFA, also occur heterogeneously within a voxel and may lead to an increase of size variance. (Ikenouchi et al., 2020) estimated μFA by using SDE under the assumption of uniform size and shape of compartmental diffusion tensors (Kaden et al., 2016) and reported a reduction of μFA in patients with PD and correlation with UPDRS-III score, in line with our observations. Although the methods of Pasternak et al. (2009) and Kaden et al. (2016) were designed to be practical for use with clinical SDE data, the dependence on regularization or strongly restrictive assumptions may bias the output parameters (Molina-Romero et al., 2018; Henriques et al., 2019) and therefore confirmation with extended acquisition, as we used in the current study, is important.
Although the differences between the controls and the patients were less prominent than the effects of age, smaller OP in the patients was identified in the left cerebral peduncle and left internal capsule (Figure 7). This overlaps the anatomical location of the medial forebrain bundle (Coenen et al., 2018), a fiber bundle that is involved early in PD (Tagliaferro et al., 2015). Our results are in line with those of previous studies that observed reduced FA within this region (Planetta et al., 2013; Zhang et al., 2015), as well as with reports of left-predominant involvement (Prakash et al., 2012; Scherfler et al., 2012).
Interestingly, we observed Kaniso and μFA were not very sensitive to the group differences between the normal subjects and the patients with PD (Figure 6), despite their correlation with the motor impairment. In contrast, MD, Kiso, and V(Dc) showed some sensitivity to the group differences in several ROIs but their correlation with motor impairment was limited to the callosal ROIs (Figure 8). Although the mechanisms behind these observations remain unknown, our speculation is that these two groups of diffusion parameters weigh different facets of neurodegeneration. It has been hypothesized that age-related changes form pathological foundation on which PD-related neurodegeneration build (Collier et al., 2017), and that there is a biological interaction between aging and PD (Levy, 2007). We speculate that microscopic anisotropy has relatively greater weight on the accumulation of age-related neuronal damage that is perhaps accelerated in PD and exacerbates the clinical expression of symptoms, while reduced microscopic anisotropy is less specific to the diagnosis of PD. On the other hand, MD, Kiso, and V(Dc) may be influenced more by the factors that differentiate PD from normal aging but are not proportional to disease progression.
Several possible confounding effects with regard to acquisition and parameter estimation need to be mentioned. First, by adopting the DTD model, which does not feature the effects of timing parameters (diffusion gradient duration, separation between the gradients, and the mixing time), we assumed the effects of timing parameters were negligible. Though this has some support by earlier studies observing the effects of diffusion time were small in brain tissue for the range above 10 ms (Clark et al., 2001; Portnoy et al., 2013), more recent studies suggest observable time-dependence in the human white matter for even longer times (Fieremans et al., 2016; Lee et al., 2020). A recent observation of intra-compartmental kurtosis in mouse brains (Henriques et al., 2020) also indicate the contribution of restricted compartments might be non-negligible. Taken together, the DTD model employed in this study may not fully capture the microstructural complexities and a more adequate model would need to feature the effects that explain the reported time-dependence, like restricted compartments (Yolcu et al., 2016; Özarslan et al., 2017), structural disorder (Novikov et al., 2014), or inter-compartmental exchange (Nilsson et al., 2013). Incorporating the timing parameters as additional measurement dimensions of DDE is expected to identify a suitable extension of the model in future. Though dependence on the DTD model is an essential limitation to our study, in Kerkelä et al. (2020), μFA estimated under this assumption approximated the gold-standard model-free method (Ianuş et al., 2018) very well in human white matter. Therefore, we believe that any bias introduced by the model assumption was not large, at least for μFA and Kaniso. Second, the expansion Eq. 2 is valid only in the vicinity of b = 0 (Westin et al., 2016; Topgaard, 2017), as in the case of DTI/DKI (Kiselev, 2013). Although the covariance tensor framework is computationally efficient, a possible bias in case of finite signal-to-noise ratio has been pointed out (Reymbaut et al., 2020). In addition, our scan–rescan analyses revealed problematic instability of the V(Dc) and Kiso estimation, which needs to be addressed before clinical application of these metrics. Third, the acquisition used in this study was not controlled for the concomitant fields (transverse magnetic field components accompanying the applied gradient to satisfy Maxwell’s equations) (Baron et al., 2012), which cause bias in the signal and hence the DTD parameters (Szczepankiewicz et al., 2019). This issue would be mitigated by the use of multiple refocusing pulses to make the diffusion encoding symmetric (Callaghan and Komlosh, 2002).
Our study has several limitations. First, PD is a heterogenous disease whose clinical manifestations vary widely among individual patients (Fereshtehnejad et al., 2017). Robust determination of PD subtypes based on cross-sectional data remains a challenge (Simuni et al., 2016) and needs to be addressed in future work. Also, though we merged the left and right ROIs for this exploratory study, neurodegeneration in PD has some asymmetry which correlates with the laterality of symptoms (Riederer et al., 2018). Further investigation with a large sample will be necessary to elucidate the relationships between white-matter degeneration and detailed clinical features, including both motor and non-motor symptoms. Second, we did not prove clinical benefit of DDE as compared to the standard SDE. Though μFA showed better correlation with age and motor impairment than FA, similar correlation could be observed with MK. Given the longer scan time of DDE, DTI/DKI based on standard SDE is still a powerful option for clinical studies. However, the strength of this study is that we showed the sensitivity of MK to neurodegeneration is likely attributed to its link to microscopic anisotropy. Such knowledge can be informative for future researches in selecting dMRI acquisition and parameter to study particular diseases. Third, we applied simplistic linear models for statistical analyses, but the effects of age might differ between the patients and the controls (Zhang et al., 2016; Pozorski et al., 2018). Also, the effect of PD might differ between males and females (Gillies et al., 2014). Although adding age × disease and sex × disease interactions in the linear model did not reveal any significant interactions for the present data (not shown), several imaging studies have suggested that these interactions are non-negligible (Dean et al., 2016; Yadav et al., 2016). The absence of these interaction effects in the present study may be attributed to the relatively small sample size which was not optimally matched regarding male/female ratio. Fourth, the majority of patients had relatively long disease duration and were already under medication, and therefore we could not quantify the effects of medication (Atkinson-Clement et al., 2017) and/or the non-linear (U-shaped) trajectories presumably related to neuronal compensation (Sanjari Moghaddam et al., 2020). Lastly, the scan–rescan repeatability should ideally be measured in a cohort similar to the patients because motion artifacts might be severer in these subjects than in young and healthy subjects. Also, the within-subject variability needs to be evaluated relative to the between-subject variability in the patients.
Conclusion
In this study, we explored the utility of DDE-derived parameters for characterizing white-matter degeneration in aging and PD. Advanced age was associated with greater mean size and size variance of the compartmental diffusion tensors and with smaller microscopic anisotropy. We found the reductions of MK in aging and PD reported in the literature are likely driven by the reduction of microscopic anisotropy. Furthermore, microscopic anisotropy correlated with the severity of motor impairment in PD. We further showed that microscopic anisotropy can be estimated with excellent repeatability by using modern clinical scanners. In conclusion, multidimensional diffusion encoding can provide more comprehensive and clinically relevant information about the white-matter degeneration in aging and PD than conventional SDE-based methods.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.
Ethics Statement
The studies involving human participants were reviewed and approved by the Institutional Review Board of Juntendo University School Hospital. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
KKami and KKama contributed to conception and design of the study. KKama, KO, TH, TO, and HT-A organized the data acquisition and quality control. KKami, SM, and CA performed the image processing. KKami performed the statistical analysis and wrote the first draft of the manuscript. KM, TF, MH, NH, and SA revised the subsections of the manuscript. All authors contributed to the final manuscript revision, read and approved the submitted version.
Funding
This work was supported by the Brain/MINDS Beyond program from AMED (Grant Nos. JP19dm0307024 and JP19dm0307101) and JSPS KAKENHI (Grant Nos. 18K07729 and 19K17244).
Conflict of Interest
KM is an employee of Siemens Healthcare K.K. TF is an employee of Siemens Healthcare GmbH.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
Ades-Aron, B., Veraart, J., Kochunov, P., McGuire, S., Sherman, P., Kellner, E., et al. (2018). Evaluation of the accuracy and precision of the diffusion parameter EStImation with Gibbs and NoisE removal pipeline. Neuroimage 183, 532–543. doi: 10.1016/j.neuroimage.2018.07.066
Agosta, F., Canu, E., Stojković, T., Pievani, M., Tomić, A., Sarro, L., et al. (2013). The topography of brain damage at different stages of Parkinson’s disease. Hum. Brain Mapp. 34, 2798–2807. doi: 10.1002/hbm.22101
Andersson, J. L. R., and Sotiropoulos, S. N. (2016). An integrated approach to correction for off-resonance effects and subject movement in diffusion MR imaging. Neuroimage 125, 1063–1078. doi: 10.1016/j.neuroimage.2015.10.019
Andersen, K. W., Lasić, S., Lundell, H., Nilsson, M., Topgaard, D., Sellebjerg, F., et al. (2020). Disentangling white-matter damage from physiological fiber orientation dispersion in multiple sclerosis. Brain Commun. 2:fcaa077. doi: 10.1093/braincomms/fcaa077
Andica, C., Kamagata, K., Hatano, T., Saito, A., Uchida, W., Ogawa, T., et al. (2019a). Free-water imaging in white and gray matter in Parkinson’s disease. Cells 8:839. doi: 10.3390/cells8080839
Andica, C., Kamagata, K., Hatano, T., Saito, Y., Ogaki, K., Hattori, N., et al. (2019b). MR biomarkers of degenerative brain disorders derived from diffusion imaging. J. Magn. Reson. Imaging. doi: 10.1002/jmri.27019 [Epub ahead of print].
Atkinson-Clement, C., Pinto, S., Eusebio, A., and Coulon, O. (2017). Diffusion tensor imaging in Parkinson’s disease: review and meta-analysis. Neuroimage Clin. 16, 98–110. doi: 10.1016/j.nicl.2017.07.011
Baron, C. A., Lebel, R. M., Wilman, A. H., and Beaulieu, C. (2012). The effect of concomitant gradient fields on diffusion tensor imaging. Magn. Reson. Med. 68, 1190–1201. doi: 10.1002/mrm.24120
Basser, P. J. (1995). Inferring microstructural features and the physiological state of tissues from diffusion-weighted images. NMR Biomed. 8, 333–344. doi: 10.1002/nbm.1940080707
Benitez, A., Jensen, J. H., Falangola, M. F., Nietert, P. J., and Helpern, J. A. (2018). Modeling white matter tract integrity in aging with diffusional kurtosis imaging. Neurobiol. Aging 70, 265–275. doi: 10.1016/j.neurobiolaging.2018.07.006
Billiet, T., Vandenbulcke, M., Mädler, B., Peeters, R., Dhollander, T., Zhang, H., et al. (2015). Age-related microstructural differences quantified using myelin water imaging and advanced diffusion MRI. Neurobiol. Aging 36, 2107–2121. doi: 10.1016/j.neurobiolaging.2015.02.029
Braak, H., Del Tredici, K., Rüb, U., de Vos, R. A., Jansen Steur, E. N., and Braak, E. (2003). Staging of brain pathology related to sporadic Parkinson’s disease. Neurobiol. Aging 24, 197–211. doi: 10.1016/S0197-4580(02)00065-9
Budde, M. D., and Frank, J. A. (2010). Neurite beading is sufficient to decrease the apparent diffusion coefficient after ischemic stroke. Proc. Natl. Acad. Sci. U.S.A. 107, 14472–14477. doi: 10.1073/pnas.1004841107
Calabrese, V., Santoro, A., Monti, D., Crupi, R., Di Paola, R., Latteri, S., et al. (2018). Aging and Parkinson’s disease: inflammaging, neuroinflammation and biological remodeling as key factors in pathogenesis. Free Radic. Biol. Med. 115, 80–91. doi: 10.1016/j.freeradbiomed.2017.10.379
Callaghan, P. T. (1997). A simple matrix formalism for spin echo analysis of restricted diffusion under generalized gradient waveforms. J. Magn. Reson. 129, 74–84. doi: 10.1006/jmre.1997.1233
Callaghan, P. T., and Komlosh, M. E. (2002). Locally anisotropic motion in a macroscopically isotropic system: displacement correlations measured using double pulsed gradient spin-echo NMR. Magn. Reson. Chem. 40, S15–S19. doi: 10.1002/mrc.1122
Callaghan, P. T., and Xia, Y. (1991). Velocity and diffusion imaging in dynamic NMR microscopy. J. Magn. Reson. 91, 326–352. doi: 10.1016/0022-2364(91)90196-Z
Caprihan, A., Wang, L. Z., and Fukushima, E. (1996). A multiple-narrow-pulse approximation for restricted diffusion in a time-varying field gradient. J. Magn. Reson. Ser. A 118, 94–102. doi: 10.1006/jmra.1996.0013
Chad, J. A., Pasternak, O., Salat, D. H., and Chen, J. J. (2018). Re-examining age-related differences in white matter microstructure with free-water corrected diffusion tensor imaging. Neurobiol. Aging 71, 161–170. doi: 10.1016/j.neurobiolaging.2018.07.018
Cheng, Y., and Cory, D. G. (1999). Multiple scattering by NMR. J. Am. Chem. Soc. 121, 7935–7936. doi: 10.1021/ja9843324
Clark, C. A., Hedehus, M., and Moseley, M. E. (2001). Diffusion time dependence of the apparent diffusion tensor in healthy human brain and white matter disease. Magn. Reson. Med. 45, 1126–1129. doi: 10.1002/mrm.1149
Coenen, V. A., Schumacher, L. V., Kaller, C., Schlaepfer, T. E., Reinacher, P. C., Egger, K., et al. (2018). The anatomy of the human medial forebrain bundle: Ventral tegmental area connections to reward-associated subcortical and frontal lobe regions. Neuroimage Clin. 18, 770–783. doi: 10.1016/j.nicl.2018.03.019
Collier, T. J., Kanaan, N. M., and Kordower, J. H. (2017). Aging and Parkinson’s disease: different sides of the same coin? Mov. Disord. 32, 983–990. doi: 10.1002/mds.27037
Cory, D., Garroway, A., and Miller, J. (1990). Applications of spin transport as a probe of local geometry. Polymer Prepr. 31, 149–150.
Coutu, J.-P., Chen, J. J., Rosas, H. D., and Salat, D. H. (2014). Non-gaussian water diffusion in aging white matter. Neurobiol. Aging 35, 1412–1421. doi: 10.1016/j.neurobiolaging.2013.12.001
Datar, A., Ameeramja, J., Bhat, A., Srivastava, R., Mishra, A., Bernal, R., et al. (2019). The roles of microtubules and membrane tension in axonal beading, retraction, and atrophy. Biophys. J. 117, 880–891. doi: 10.1016/j.bpj.2019.07.046
Dean, D. C., Sojkova, J., Hurley, S., Kecskemeti, S., Okonkwo, O., Bendlin, B. B., et al. (2016). Alterations of Myelin Content in Parkinson’s disease: a cross-sectional neuroimaging study. PLoS One 11:e0163774. doi: 10.1371/journal.pone.0163774
Douaud, G., Jbabdi, S., Behrens, T. E. J., Menke, R. A., Gass, A., Monsch, A. U., et al. (2011). DTI measures in crossing-fibre areas: increased diffusion anisotropy reveals early white matter alteration in MCI and mild Alzheimer’s disease. Neuroimage 55, 880–890. doi: 10.1016/j.neuroimage.2010.12.008
Duncan, G. W., Firbank, M. J., Yarnall, A. J., Khoo, T. K., Brooks, D. J., Barker, R. A., et al. (2016). Gray and white matter imaging: a biomarker for cognitive impairment in early Parkinson’s disease? Mov. Disord. 31, 103–110. doi: 10.1002/mds.26312
Eriksson, S., Lasic, S., and Topgaard, D. (2013). Isotropic diffusion weighting in PGSE NMR by magic-angle spinning of the q-vector. J. Magn. Reson. 226, 13–18. doi: 10.1016/j.jmr.2012.10.015
Fazekas, F., Chawluk, J. B., Alavi, A., Hurtig, H. I., and Zimmerman, R. A. (1987). MR signal abnormalities at 1.5 T in Alzheimer’s dementia and normal aging. Am. J. Roentgenol. 149, 351–356. doi: 10.2214/ajr.149.2.351
Fereshtehnejad, S.-M., Zeighami, Y., Dagher, A., and Postuma, R. B. (2017). Clinical criteria for subtyping Parkinson’s disease: biomarkers and longitudinal progression. Brain 140, 1959–1976. doi: 10.1093/brain/awx118
Fieremans, E., Burcaw, L. M., Lee, H. H., Lemberskiy, G., Veraart, J., and Novikov, D. S. (2016). In vivo observation and biophysical interpretation of time-dependent diffusion in human white matter. Neuroimage 129, 414–427. doi: 10.1016/j.neuroimage.2016.01.018
Gillies, G. E., Pienaar, I. S., Vohra, S., and Qamhawi, Z. (2014). Sex differences in Parkinson’s disease. Front. Neuroendocrinol. 35, 370–384. doi: 10.1016/j.yfrne.2014.02.002
Goetz, C. G., Tilley, B. C., Shaftman, S. R., Stebbins, G. T., Fahn, S., Martinez-Martin, P., et al. (2008). Movement Disorder Society-sponsored revision of the Unified Parkinson’s Disease Rating Scale (MDS-UPDRS): scale presentation and clinimetric testing results. Mov. Disord. 23, 2129–2170. doi: 10.1002/mds.22340
Goveas, J., O’Dwyer, L., Mascalchi, M., Cosottini, M., Diciotti, S., De Santis, S., et al. (2015). Diffusion-MRI in neurodegenerative disorders. Magn. Reson. Imaging 33, 853–876. doi: 10.1016/j.mri.2015.04.006
Guerreri, M., Palombo, M., Caporale, A., Fasano, F., Macaluso, E., Bozzali, M., et al. (2019). Age-related microstructural and physiological changes in normal brain measured by MRI γ-metrics derived from anomalous diffusion signal representation. Neuroimage 188, 654–667. doi: 10.1016/j.neuroimage.2018.12.044
Hair, J., Black, W. C., Babin, B. J., and Anderson, R. E. (2010). Multivariate Data Analysis, 7th Edn. Upper Saddle River, NJ: Pearson Education.
Hansen, B., Lund, T. E., Sangill, R., and Jespersen, S. N. (2013). Experimentally and computationally fast method for estimation of a mean kurtosis. Magn. Reson. Med. 69, 1754–1760. doi: 10.1002/mrm.24743
Henriques, R. N., Jespersen, S. N., and Shemesh, N. (2019). Microscopic anisotropy misestimation in spherical-mean single diffusion encoding MRI. Magn. Reson. Med. 81, 3245–3261. doi: 10.1002/mrm.27606
Henriques, R. N., Jespersen, S. N., and Shemesh, N. (2020). Correlation tensor magnetic resonance imaging. Neuroimage 211:116605. doi: 10.1016/j.neuroimage.2020.116605
Hoehn, M. M., and Yahr, M. D. (1967). Parkinsonism: onset, progression, and mortality. Neurology 17, 427–442. doi: 10.1212/WNL.17.5.427
Ianuş, A., Jespersen, S. N., Serradas Duarte, T., Alexander, D. C., Drobnjak, I., and Shemesh, N. (2018). Accurate estimation of microscopic diffusion anisotropy and its time dependence in the mouse brain. Neuroimage 183, 934–949. doi: 10.1016/j.neuroimage.2018.08.034
Ikenouchi, Y., Kamagata, K., Andica, C., Hatano, T., Ogawa, T., Takeshige-Amano, H., et al. (2020). Evaluation of white matter microstructure in patients with Parkinson’s disease using microscopic fractional anisotropy. Neuroradiology 62, 197–203. doi: 10.1007/s00234-019-02301-1
Jelescu, I. O., and Budde, M. D. (2017). Design and validation of diffusion MRI models of white matter. Front. Phys. 28:61. doi: 10.3389/fphy.2017.00061
Jenkinson, M., Beckmann, C. F., Behrens, T. E. J., Woolrich, M. W., and Smith, S. M. (2012). FSL. Neuroimage 62, 782–790. doi: 10.1016/j.neuroimage.2011.09.015
Jensen, J. H., Helpern, J. A., Ramani, A., Lu, H., and Kaczynski, K. (2005). Diffusional kurtosis imaging: the quantification of non-Gaussian water diffusion by means of magnetic resonance imaging. Magn. Reson. Med. 53, 1432–1440. doi: 10.1002/mrm.20508
Jespersen, S. N. (2012). Equivalence of double and single wave vector diffusion contrast at low diffusion weighting. NMR Biomed. 25, 813–818. doi: 10.1002/nbm.1808
Jespersen, S. N., Lundell, H., Sønderby, C. K., and Dyrby, T. B. (2013). Orientationally invariant metrics of apparent compartment eccentricity from double pulsed field gradient diffusion experiments. NMR Biomed. 26, 1647–1662. doi: 10.1002/nbm.2999
Jian, B., Vemuri, B. C., Özarslan, E., Carney, P. R., and Mareci, T. H. (2007). A novel tensor distribution model for the diffusion-weighted MR signal. Neuroimage 37, 164–176. doi: 10.1016/j.neuroimage.2007.03.074
Kaden, E., Kruggel, F., and Alexander, D. C. (2016). Quantitative mapping of the per-axon diffusion coefficients in brain white matter. Magn. Reson. Med. 75, 1752–1763. doi: 10.1002/mrm.25734
Kamagata, K., Tomiyama, H., Hatano, T., Motoi, Y., Abe, O., Shimoji, K., et al. (2014). A preliminary diffusional kurtosis imaging study of Parkinson disease: comparison with conventional diffusion tensor imaging. Neuroradiology 56, 251–258. doi: 10.1007/s00234-014-1327-1
Kamagata, K., Tomiyama, H., Motoi, Y., Kano, M., Abe, O., Ito, K., et al. (2013). Diffusional kurtosis imaging of cingulate fibers in Parkinson disease: comparison with conventional diffusion tensor imaging. Magn. Reson. Imaging 31, 1501–1506. doi: 10.1016/j.mri.2013.06.009
Kanazawa, T., Adachi, E., Orimo, S., Nakamura, A., Mizusawa, H., and Uchihara, T. (2012). Pale neurites, premature α-synuclein aggregates with centripetal extension from axon collaterals. Brain Pathol. 22, 67–78. doi: 10.1111/j.1750-3639.2011.00509.x
Kellner, E., Dhital, B., Kiselev, V. G., and Reisert, M. (2016). Gibbs-ringing artifact removal based on local subvoxel-shifts. Magn. Reson. Med. 76, 1574–1581. doi: 10.1002/mrm.26054
Kerkelä, L., Henriques, R. N., Hall, M. G., Clark, C. A., Shemesh, N., Hall, M. G., et al. (2020). Validation and noise robustness assessment of microscopic anisotropy estimation with clinically feasible double diffusion encoding MRI. Magn. Reson. Med. 83, 1698–1710. doi: 10.1002/mrm.28048
Kiselev, V. G. (2013). “The cumulant expansion: an overarching mathematical framework for understanding diffusion NMR,” in Diffusion MRI, ed. D. K. Jones (Oxford: Oxford University Press). doi: 10.1093/med/9780195369779.003.0010
Komlosh, M. E., Horkay, F., Freidlin, R. Z., Nevo, U., Assaf, Y., and Basser, P. J. (2007). Detection of microscopic anisotropy in gray matter and in a novel tissue phantom using double Pulsed Gradient Spin Echo MR. J. Magn. Reson. 189, 38–45. doi: 10.1016/j.jmr.2007.07.003
Komlosh, M. E., Lizak, M. J., Horkay, F., Freidlin, R. Z., and Basser, P. J. (2008). Observation of microscopic diffusion anisotropy in the spinal cord using double-pulsed gradient spin echo MRI. Magn. Reson. Med. 59, 803–809. doi: 10.1002/mrm.21528
Kurowska, Z., Kordower, J. H., Stoessl, A. J., Burke, R. E., Brundin, P., Yue, Z., et al. (2016). Is axonal degeneration a key early event in Parkinson’s disease? J. Parkinsons Dis. 6, 703–707. doi: 10.3233/JPD-160881
Lampinen, B., Szczepankiewicz, F., Novén, M., van Westen, D., Hansson, O., Englund, E., et al. (2019). Searching for the neurite density with diffusion MRI: challenges for biophysical modeling. Hum. Brain Mapp. 40, 2529–2545. doi: 10.1002/hbm.24542
Lasič, S., Szczepankiewicz, F., Eriksson, S., Nilsson, M., and Topgaard, D. (2014). Microanisotropy imaging: quantification of microscopic diffusion anisotropy and orientational order parameter by diffusion MRI with magic-angle spinning of the q-vector. Front. Phys. 2:11. doi: 10.3389/fphy.2014.00011
Lawrenz, M., Brassen, S., and Finsterbusch, J. (2015). Microscopic diffusion anisotropy in the human brain: reproducibility, normal values, and comparison with the fractional anisotropy. Neuroimage 109, 283–297. doi: 10.1016/j.neuroimage.2015.01.025
Lawrenz, M., Brassen, S., and Finsterbusch, J. (2016). Microscopic diffusion anisotropy in the human brain: age-related changes. Neuroimage 141, 313–325. doi: 10.1016/j.neuroimage.2016.07.031
Lawrenz, M., and Finsterbusch, J. (2013). Double-wave-vector diffusion-weighted imaging reveals microscopic diffusion anisotropy in the living human brain. Magn. Reson. Med. 69, 1072–1082. doi: 10.1002/mrm.24347
Lawrenz, M., Koch, M. A., and Finsterbusch, J. (2010). A tensor model and measures of microscopic anisotropy for double-wave-vector diffusion-weighting experiments with long mixing times. J. Magn. Reson. 202, 43–56. doi: 10.1016/j.jmr.2009.09.015
Lee, H.-H., Papaioannou, A., Kim, S.-L., Novikov, D. S., and Fieremans, E. (2020). A time-dependent diffusion MRI signature of axon caliber variations and beading. Commun. Biol. 3:354. doi: 10.1038/s42003-020-1050-x
Levy, G. (2007). The relationship of Parkinson disease with aging. Arch. Neurol. 64, 1242–1246. doi: 10.1001/archneur.64.9.1242
Lundell, H., Nilsson, M., Dyrby, T. B., Parker, G. J. M., Cristinacce, P. L. H., Zhou, F.-L., et al. (2019). Multidimensional diffusion MRI with spectrally modulated gradients reveals unprecedented microstructural detail. Sci. Rep. 9:9026.
Madden, D. J., Bennett, I. J., Burzynska, A., Potter, G. G., Chen, N. K., et al. (2012). Diffusion tensor imaging of cerebral white matter integrity in cognitive aging. Biochim. Biophys. Acta 1822, 386–400. doi: 10.1016/j.bbadis.2011.08.003
McDonald, C., Newton, J. L., and Burn, D. J. (2016). Orthostatic hypotension and cognitive impairment in Parkinson’s disease: causation or association? Mov. Disord. 31, 937–946. doi: 10.1002/mds.26632
Mitra, P. P. (1995). Multiple wave-vector extensions of the NMR pulsed-field-gradient spin-echo diffusion measurement. Phys. Rev. B Condens. Matter 51, 15074–15078. doi: 10.1103/PhysRevB.51.15074
Molina-Romero, M., Menze, B. H., Pedro, A. G., Czisch, M., Sämann, P. G., Jones, D. K., et al. (2018). A diffusion model-free framework with echo time dependence for free-water elimination and brain tissue microstructure characterization. Magn. Reson. Med. 80, 2155–2172. doi: 10.1002/mrm.27181
Mori, S., and Van Zijl, P. C. M. (1995). Diffusion weighting by the trace of the diffusion tensor within a single scan. Magn. Reson. Med. 33, 41–52. doi: 10.1002/mrm.1910330107
Nilsson, M., Lätt, J., Van Westen, D., Brockstedt, S., Lasič, S., Ståhlberg, F., et al. (2013). Noninvasive mapping of water diffusional exchange in the human brain using filter-exchange imaging. Magn. Reson. Med. 69, 1573–1581. doi: 10.1002/mrm.24395
Nilsson, M., Szczepankiewicz, F., Lampinen, B., Ahlgren, A., de Almeida Martins, J. P., Lasic, S., et al. (2018). “An open-source framework for analysis of multidimensional diffusion MRI data implemented in MATLAB,” in Proceedings of the 26th Annual Meeting of ISMRM, Paris, France.
Nilsson, M., Szczepankiewicz, F., Brabec, J., Taylor, M., Westin, C.-F. C., Golby, A., et al. (2020). Tensor-valued diffusion MRI in under 3 minutes: an initial survey of microscopic anisotropy and tissue heterogeneity in intracranial tumors. Magn. Reson. Med. 83, 608–620. doi: 10.1002/mrm.27959
Novikov, D. S., Fieremans, E., Jespersen, S. N., Kiselev, V. G., Mri, D., and Bird, T. A. (2019). Quantifying brain microstructure with diffusion MRI: theory and parameter estimation. NMR Biomed. 32:e3998. doi: 10.1002/nbm.3998
Novikov, D. S., Jensen, J. H., Helpern, J. A., and Fieremans, E. (2014). Revealing mesoscopic structural universality with diffusion. Proc. Natl. Acad. Sci. U.S.A. 111, 5088–5093. doi: 10.1073/pnas.1316944111
Özarslan, E., and Basser, P. J. (2008). Microscopic anisotropy revealed by NMR double pulsed field gradient experiments with arbitrary timing parameters. J. Chem. Phys. 128, 154511. doi: 10.1063/1.2905765
Özarslan, E., Yolcu, C., Herberthson, M., Westin, C.-F., and Knutsson, H. (2017). Effective potential for magnetic resonance measurements of restricted diffusion. Front. Phys. 5:68. doi: 10.3389/fphy.2017.00068
Palombo, M., Ligneul, C., Hernandez-Garzon, E., and Valette, J. (2018). Can we detect the effect of spines and leaflets on the diffusion of brain intracellular metabolites? Neuroimage 182, 283–293. doi: 10.1016/j.neuroimage.2017.05.003
Pasternak, O., Sochen, N., Gur, Y., Intrator, N., and Assaf, Y. (2009). Free water elimination and mapping from diffusion MRI. Magn. Reson. Med. 62, 717–730. doi: 10.1002/mrm.22055
Planetta, P. J., Schulze, E. T., Geary, E. K., Corcos, D. M., Goldman, J. G., Little, D. M., et al. (2013). Thalamic projection fiber integrity in de novo Parkinson Disease. Am. J. Neuroradiol. 34, 74–79. doi: 10.3174/ajnr.A3178
Poewe, W., Seppi, K., Tanner, C. M., Halliday, G. M., Brundin, P., Volkmann, J., et al. (2017). Parkinson disease. Nat. Rev. Dis. Primers 3:17013. doi: 10.1038/nrdp.2017.13
Portnoy, S., Flint, J. J., Blackband, S. J., and Stanisz, G. J. (2013). Oscillating and pulsed gradient diffusion magnetic resonance microscopy over an extended b-value range: implications for the characterization of tissue microstructure. Magn. Reson. Med. 69, 1131–1145. doi: 10.1002/mrm.24325
Postuma, R. B., Berg, D., Stern, M., Poewe, W., Olanow, C. W., Oertel, W., et al. (2015). MDS clinical diagnostic criteria for Parkinson’s disease. Mov. Disord. 30, 1591–1601. doi: 10.1002/mds.26424
Pozorski, V., Oh, J. M., Adluru, N., Merluzzi, A. P., Theisen, F., Okonkwo, O., et al. (2018). Longitudinal white matter microstructural change in Parkinson’s disease. Hum. Brain Mapp. 39, 4150–4161. doi: 10.1002/hbm.24239
Prakash, B. D., Sitoh, Y.-Y., Tan, L. C. S., and Au, W. L. (2012). Asymmetrical diffusion tensor imaging indices of the rostral substantia nigra in Parkinson’s disease. Parkinsonism Relat. Disord. 18, 1029–1033. doi: 10.1016/j.parkreldis.2012.05.021
Rektor, I., Svátková, A., Vojtíšek, L., Zikmundová, I., Vaníček, J., Király, A., et al. (2018). White matter alterations in Parkinson’s disease with normal cognition precede grey matter atrophy. PLoS One 13:e0187939. doi: 10.1371/journal.pone.0187939
Reymbaut, A., Mezzani, P., Almeida Martins, J. P., and Topgaard, D. (2020). Accuracy and precision of statistical descriptors obtained from multidimensional diffusion signal inversion algorithms. NMR Biomed. 17:e4267. doi: 10.1002/nbm.4267
Riederer, P., Jellinger, K. A., Kolber, P., Hipp, G., Sian-Hülsmann, J., and Krüger, R. (2018). Lateralisation in Parkinson disease. Cell Tissue Res. 373, 297–312. doi: 10.1007/s00441-018-2832-z
Salvadores, N., Sanhueza, M., Manque, P., and Court, F. A. (2017). Axonal degeneration during aging and its functional role in neurodegenerative disorders. Front. Neurosci. 11:451. doi: 10.3389/fnins.2017.00451
Sanjari Moghaddam, H., Dolatshahi, M., Mohebi, F., and Aarabi, M. H. (2020). Structural white matter alterations as compensatory mechanisms in Parkinson’s disease: a systematic review of diffusion tensor imaging studies. J. Neurosci. Res. 98, 1398–1416. doi: 10.1002/jnr.24617
Scherfler, C., Seppi, K., Mair, K. J., Donnemiller, E., Virgolini, I., Wenning, G. K., et al. (2012). Left hemispheric predominance of nigrostriatal dysfunction in Parkinson’s disease. Brain 135, 3348–3354. doi: 10.1093/brain/aws253
Shahim, P., Holleran, L., Kim, J. H., and Brody, D. L. (2017). Test-retest reliability of high spatial resolution diffusion tensor and diffusion kurtosis imaging. Sci. Rep. 7:11141. doi: 10.1038/s41598-017-11747-3
Shemesh, N., Barazany, D., Sadan, O., Bar, L., Zur, Y., Barhum, Y., et al. (2012). Mapping apparent eccentricity and residual ensemble anisotropy in the gray matter using angular double-pulsed-field-gradient MRI. Magn. Reson. Med. 68, 794–806. doi: 10.1002/mrm.23300
Shemesh, N., and Cohen, Y. (2011). Microscopic and compartment shape anisotropies in gray and white matter revealed by angular bipolar double-PFG MR. Magn. Reson. Med. 65, 1216–1227. doi: 10.1002/mrm.22738
Shemesh, N., Özarslan, E., Basser, P. J., and Cohen, Y. (2009). Measuring small compartmental dimensions with low-q angular double-PGSE NMR: the effect of experimental parameters on signal decay. J. Magn. Reson. 198, 15–23. doi: 10.1016/j.jmr.2009.01.004
Simuni, T., Caspell-Garcia, C., Coffey, C., Lasch, S., Tanner, C., and Marek, K. (2016). How stable are Parkinson’s disease subtypes in de novo patients: analysis of the PPMI cohort? Parkinsonism Relat. Disord. 28, 62–67. doi: 10.1016/j.parkreldis.2016.04.027
Smith, S. M., Jenkinson, M., Johansen-Berg, H., Rueckert, D., Nichols, T. E., Mackay, C. E., et al. (2006). Tract-based spatial statistics: voxelwise analysis of multi-subject diffusion data. Neuroimage 31, 1487–1505. doi: 10.1016/j.neuroimage.2006.02.024
Smith, S. M., and Nichols, T. E. (2009). Threshold-free cluster enhancement: addressing problems of smoothing, threshold dependence and localisation in cluster inference. Neuroimage 44, 83–98. doi: 10.1016/j.neuroimage.2008.03.061
Surova, Y., Lampinen, B., Nilsson, M., Lätt, J., Hall, S., Widner, H., et al. (2016). Alterations of diffusion kurtosis and neurite density measures in deep grey matter and white matter in Parkinson’s disease. PLoS One 11:e0157755. doi: 10.1371/journal.pone.0157755
Surova, Y., Nilsson, M., Lampinen, B., Lätt, J., Hall, S., Widner, H., et al. (2018). Alteration of putaminal fractional anisotropy in Parkinson’s disease: a longitudinal diffusion kurtosis imaging study. Neuroradiology 60, 247–254. doi: 10.1007/s00234-017-1971-3
Szczepankiewicz, F., Lasi, S., Westen, D., Van, Sundgren, P. C., Englund, E., et al. (2015). Quantification of microscopic diffusion anisotropy disentangles effects of orientation dispersion from microstructure: applications in healthy volunteers and in brain tumors. Neuroimage 104, 241–252. doi: 10.1016/j.neuroimage.2014.09.057
Szczepankiewicz, F., van Westen, D., Englund, E., Westin, C. F., Ståhlberg, F., Lätt, J., et al. (2016). The link between diffusion MRI and tumor heterogeneity: Mapping cell eccentricity and density by diffusional variance decomposition (DIVIDE). Neuroimage 142, 522–532. doi: 10.1016/j.neuroimage.2016.07.038
Szczepankiewicz, F., Westin, C., and Nilsson, M. (2019). Maxwell-compensated design of asymmetric gradient waveforms for tensor-valued diffusion encoding. Magn. Reson. Med. 82, 1424–1437. doi: 10.1002/mrm.27828
Tagliaferro, P., Kareva, T., Oo, T. F., Yarygina, O., Kholodilov, N., and Burke, R. E. (2015). An early axonopathy in a hLRRK2(R1441G) transgenic model of Parkinson disease. Neurobiol. Dis. 82, 359–371. doi: 10.1016/j.nbd.2015.07.009
Topgaard, D. (2017). Multidimensional diffusion MRI. J. Magn. Reson. 275, 98–113. doi: 10.1016/j.jmr.2016.12.007
Tournier, J.-D., Smith, R., Raffelt, D., Tabbara, R., Dhollander, T., Pietsch, M., et al. (2019). MRtrix3: a fast, flexible and open software framework for medical image processing and visualisation. Neuroimage 202:116137. doi: 10.1016/j.neuroimage.2019.116137
Tustison, N. J., Avants, B. B., Cook, P. A., Zheng, Y, Egan, A., Yushkevich, P. A., et al. (2010). N4ITK: improved N3 bias correction. IEEE Trans. Med. Imaging 29, 1310–1320. doi: 10.1109/TMI.2010.2046908
Veraart, J., Fieremans, E., and Novikov, D. S. (2016). Diffusion MRI noise mapping using random matrix theory. Magn. Reson. Med. 76, 1582–1593. doi: 10.1002/mrm.26059
Veselý, B., Antonini, A., and Rektor, I. (2016). The contribution of white matter lesions to Parkinson’s disease motor and gait symptoms: a critical review of the literature. J. Neural Transm. 123, 241–250. doi: 10.1007/s00702-015-1470-9
Wang, J.-J., Lin, W.-Y., Lu, C., Weng, Y.-H., Ng, S.-H., Wang, C.-H., et al. (2011). Parkinson disease: diagnostic utility of diffusion kurtosis imaging. Radiology 261, 210–217. doi: 10.1148/radiol.11102277
Westin, C. F., Knutsson, H., Pasternak, O., Szczepankiewicz, F., Özarslan, E., van Westen, D., et al. (2016). Q-space trajectory imaging for multidimensional diffusion MRI of the human brain. Neuroimage 135, 345–362. doi: 10.1016/j.neuroimage.2016.02.039
Winkler, A. M., Ridgway, G. R., Douaud, G., Nichols, T. E., and Smith, S. M. (2016). Faster permutation inference in brain imaging. Neuroimage 141, 502–516. doi: 10.1016/j.neuroimage.2016.05.068
Winkler, A. M., Ridgway, G. R., Webster, M. A., Smith, S. M., and Nichols, T. E. (2014). Permutation inference for the general linear model. Neuroimage 92, 381–397. doi: 10.1016/j.neuroimage.2014.01.060
Yadav, S. K., Kathiresan, N., Mohan, S., Vasileiou, G., Singh, A., Kaura, D., et al. (2016). Gender-based analysis of cortical thickness and structural connectivity in Parkinson’s disease. J. Neurol. 263, 2308–2318. doi: 10.1007/s00415-016-8265-2
Yang, G., Tian, Q., Leuze, C., Wintermark, M., and McNab, J. A. (2018). Double diffusion encoding MRI for the clinic. Magn. Reson. Med. 80, 507–520. doi: 10.1002/mrm.27043
Yolcu, C., Memiç, M., Şimşek, K., Westin, C.-F., and Özarslan, E. (2016). NMR signal for particles diffusing under potentials: from path integrals and numerical methods to a model of diffusion anisotropy. Phys. Rev. E 93:052602. doi: 10.1103/PhysRevE.93.052602
Zhang, Y., Wu, I.-W., Buckley, S., Coffey, C. S., Foster, E., Mendick, S., et al. (2015). Diffusion tensor imaging of the nigrostriatal fibers in Parkinson’s disease. Mov. Disord. 30, 1229–1236. doi: 10.1002/mds.26251
Keywords: aging, Parkinson disease, diffusion MRI, double diffusion encoding, microstructure
Citation: Kamiya K, Kamagata K, Ogaki K, Hatano T, Ogawa T, Takeshige-Amano H, Murata S, Andica C, Murata K, Feiweier T, Hori M, Hattori N and Aoki S (2020) Brain White-Matter Degeneration Due to Aging and Parkinson Disease as Revealed by Double Diffusion Encoding. Front. Neurosci. 14:584510. doi: 10.3389/fnins.2020.584510
Received: 17 July 2020; Accepted: 22 September 2020;
Published: 15 October 2020.
Edited by:
Yu-Chien Wu, Indiana University Bloomington, United StatesReviewed by:
Evren Özarslan, Linköping University, SwedenMatthew ManHin Cheung, CUHK Medical Centre, Hong Kong
Copyright © 2020 Kamiya, Kamagata, Ogaki, Hatano, Ogawa, Takeshige-Amano, Murata, Andica, Murata, Feiweier, Hori, Hattori and Aoki. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Kouhei Kamiya, kkamiya-tky@umin.ac.jp