Abstract
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; ). 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) (; 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 (). 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; ; ; ; ) and diseases (; ). In PD, the majority of studies have applied diffusion tensor imaging (DTI) () 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 ). 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 (; ; 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) (; ), triple diffusion encoding (Mori and Van Zijl, 1995), or continuous gradient waveforms (; ; ), 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 (; ; ; Ö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; ). 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; ). 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 ; 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 (; 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; ). 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 (; Westin et al., 2016; ) 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; ):
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 , can be decomposed into two kurtosis sources (Westin et al., 2016; ):
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; )
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) () and the Hoehn and Yahr staging scale (). White-matter T2 hyperintensities (WMH) were rated according to the Fazekas scale (). The clinical and demographic characteristics of the participants are summarized in Table 1.
TABLE 1
| Healthy controls | Patients with PD | P* | |
| Number | 23 | 27 | – |
| Age | 63.1 ± 7.2 | 66.1 ± 6.6 | 0.14 |
| Sex (male/female) | 8/15 | 16/11 | 0.08 |
| Disease duration (years) | – | 13.4 ± 7.1 | – |
| UPDRS-III | – | 13.6 ± 7.7 | – |
| Hoehn–Yahr stage | – | 1.7 ± 0.8 | – |
| Levodopa equivalent dose (mg) | – | 1040 ± 556 | – |
| Fazekas periventricular white matter | 0.57 ± 0.66 | 0.63 ± 0.69 | 0.74 |
| Fazekas deep white matter | 0.57 ± 0.73 | 0.70 ± 0.67 | 0.49 |
Demographic features of the study participants.
*Group differences between patients and controls were examined by using the χ2 test for sex and Welch’s t-test for age and Fazekas scales. The significance threshold was set at P < 0.05.
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 (), 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 (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
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 ().
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) (). Although we used a stringent criterion of VIF < 4.0 (), 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
TABLE 2
| ROIs | MD | FA | MK | Kaniso | Kiso | μFA | OP | V(Dc) | MDSDE | FASDE | MKSDE |
| Frontal | 0.5 | 1.1 | 1.3 | 2.1 | 0.8 | 0.6 | 0.8 | 1.4 | 1.1 | 3.2 | 0.7 |
| Sensory- motor | 0.5 | 0.8 | 0.5 | 1.8 | 1.9 | 0.5 | 0.9 | 2.9 | 0.4 | 1.4 | 0.6 |
| Parietal | 0.4 | 0.7 | 0.6 | 2.6 | 3.0 | 0.9 | 0.9 | 3.1 | 0.4 | 0.8 | 0.8 |
| Occipital | 0.6 | 1.1 | 1.0 | 3.0 | 3.9 | 1.0 | 0.9 | 4.2 | 0.5 | 0.6 | 0.8 |
| Temporal | 0.6 | 1.0 | 0.7 | 3.3 | 3.3 | 1.5 | 1.0 | 3.6 | 0.6 | 2.5 | 1.5 |
| ALIC | 0.9 | 1.1 | 3.4 | 2.3 | 9.8 | 1.0 | 1.3 | 10.3 | 0.6 | 1.6 | 2.0 |
| PLIC | 0.6 | 0.8 | 1.7 | 1.7 | 6.1 | 0.5 | 0.4 | 6.4 | 1.4 | 1.9 | 1.2 |
| ACR | 1.1 | 1.6 | 1.3 | 3.6 | 3.5 | 1.2 | 0.5 | 4.5 | 0.7 | 1.5 | 1.0 |
| SCR | 0.4 | 1.3 | 1.0 | 2.1 | 3.2 | 0.7 | 0.7 | 3.5 | 0.8 | 1.1 | 1.0 |
| PCR | 0.6 | 0.9 | 1.2 | 2.2 | 1.7 | 0.7 | 0.5 | 2.2 | 1.4 | 0.9 | 0.7 |
| GCC | 0.7 | 1.2 | 4.8 | 4.2 | 10.9 | 0.9 | 0.3 | 7.5 | 3.9 | 2.7 | 2.4 |
| SCC | 0.5 | 0.6 | 1.2 | 2.2 | 4.0 | 0.4 | 0.7 | 3.4 | 2.5 | 1.1 | 2.1 |
Within-subject coefficients of variation of the ROI means (values are percentages).
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.
Correlation With Age in the Healthy Subjects
To display the results of ROI-based correlations, we adopted a figure format used in and (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
FIGURE 4
FIGURE 5
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
FIGURE 7
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
FIGURE 9
FIGURE 10
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; ; ; ; ) 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; ; Salvadores et al., 2017; ; ). 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 (; ; ) 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 , 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; ). 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. () and () 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 (; 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; ) 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 (), 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 (), 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 (; Portnoy et al., 2013), more recent studies suggest observable time-dependence in the human white matter for even longer times (; Lee et al., 2020). A recent observation of intra-compartmental kurtosis in mouse brains () 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) (), 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 ().
Our study has several limitations. First, PD is a heterogenous disease whose clinical manifestations vary widely among individual patients (). 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 (). 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 (; 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 () 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.
Statements
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
1
Ades-AronB.VeraartJ.KochunovP.McGuireS.ShermanP.KellnerE.et al (2018). Evaluation of the accuracy and precision of the diffusion parameter EStImation with Gibbs and NoisE removal pipeline.Neuroimage183532–543. 10.1016/j.neuroimage.2018.07.066
2
AgostaF.CanuE.StojkovićT.PievaniM.TomićA.SarroL.et al (2013). The topography of brain damage at different stages of Parkinson’s disease.Hum. Brain Mapp.342798–2807. 10.1002/hbm.22101
3
AnderssonJ. L. R.SotiropoulosS. N. (2016). An integrated approach to correction for off-resonance effects and subject movement in diffusion MR imaging.Neuroimage1251063–1078. 10.1016/j.neuroimage.2015.10.019
4
AndersenK. W.LasićS.LundellH.NilssonM.TopgaardD.SellebjergF.et al (2020). Disentangling white-matter damage from physiological fiber orientation dispersion in multiple sclerosis.Brain Commun.2:fcaa077. 10.1093/braincomms/fcaa077
5
AndicaC.KamagataK.HatanoT.SaitoA.UchidaW.OgawaT.et al (2019a). Free-water imaging in white and gray matter in Parkinson’s disease.Cells8:839. 10.3390/cells8080839
6
AndicaC.KamagataK.HatanoT.SaitoY.OgakiK.HattoriN.et al (2019b). MR biomarkers of degenerative brain disorders derived from diffusion imaging.J. Magn. Reson. Imaging.10.1002/jmri.27019 [Epub ahead of print].
7
Atkinson-ClementC.PintoS.EusebioA.CoulonO. (2017). Diffusion tensor imaging in Parkinson’s disease: review and meta-analysis.Neuroimage Clin.1698–110. 10.1016/j.nicl.2017.07.011
8
BaronC. A.LebelR. M.WilmanA. H.BeaulieuC. (2012). The effect of concomitant gradient fields on diffusion tensor imaging.Magn. Reson. Med.681190–1201. 10.1002/mrm.24120
9
BasserP. J. (1995). Inferring microstructural features and the physiological state of tissues from diffusion-weighted images.NMR Biomed.8333–344. 10.1002/nbm.1940080707
10
BenitezA.JensenJ. H.FalangolaM. F.NietertP. J.HelpernJ. A. (2018). Modeling white matter tract integrity in aging with diffusional kurtosis imaging.Neurobiol. Aging70265–275. 10.1016/j.neurobiolaging.2018.07.006
11
BillietT.VandenbulckeM.MädlerB.PeetersR.DhollanderT.ZhangH.et al (2015). Age-related microstructural differences quantified using myelin water imaging and advanced diffusion MRI.Neurobiol. Aging362107–2121. 10.1016/j.neurobiolaging.2015.02.029
12
BraakH.Del TrediciK.RübU.de VosR. A.Jansen SteurE. N.BraakE. (2003). Staging of brain pathology related to sporadic Parkinson’s disease.Neurobiol. Aging24197–211. 10.1016/S0197-4580(02)00065-9
13
BuddeM. D.FrankJ. A. (2010). Neurite beading is sufficient to decrease the apparent diffusion coefficient after ischemic stroke.Proc. Natl. Acad. Sci. U.S.A.10714472–14477. 10.1073/pnas.1004841107
14
CalabreseV.SantoroA.MontiD.CrupiR.Di PaolaR.LatteriS.et al (2018). Aging and Parkinson’s disease: inflammaging, neuroinflammation and biological remodeling as key factors in pathogenesis.Free Radic. Biol. Med.11580–91. 10.1016/j.freeradbiomed.2017.10.379
15
CallaghanP. T. (1997). A simple matrix formalism for spin echo analysis of restricted diffusion under generalized gradient waveforms.J. Magn. Reson.12974–84. 10.1006/jmre.1997.1233
16
CallaghanP. T.KomloshM. E. (2002). Locally anisotropic motion in a macroscopically isotropic system: displacement correlations measured using double pulsed gradient spin-echo NMR.Magn. Reson. Chem.40S15–S19. 10.1002/mrc.1122
17
CallaghanP. T.XiaY. (1991). Velocity and diffusion imaging in dynamic NMR microscopy.J. Magn. Reson.91326–352. 10.1016/0022-2364(91)90196-Z
18
CaprihanA.WangL. Z.FukushimaE. (1996). A multiple-narrow-pulse approximation for restricted diffusion in a time-varying field gradient.J. Magn. Reson. Ser. A11894–102. 10.1006/jmra.1996.0013
19
ChadJ. A.PasternakO.SalatD. H.ChenJ. J. (2018). Re-examining age-related differences in white matter microstructure with free-water corrected diffusion tensor imaging.Neurobiol. Aging71161–170. 10.1016/j.neurobiolaging.2018.07.018
20
ChengY.CoryD. G. (1999). Multiple scattering by NMR.J. Am. Chem. Soc.1217935–7936. 10.1021/ja9843324
21
ClarkC. A.HedehusM.MoseleyM. E. (2001). Diffusion time dependence of the apparent diffusion tensor in healthy human brain and white matter disease.Magn. Reson. Med.451126–1129. 10.1002/mrm.1149
22
CoenenV. A.SchumacherL. V.KallerC.SchlaepferT. E.ReinacherP. C.EggerK.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.18770–783. 10.1016/j.nicl.2018.03.019
23
CollierT. J.KanaanN. M.KordowerJ. H. (2017). Aging and Parkinson’s disease: different sides of the same coin?Mov. Disord.32983–990. 10.1002/mds.27037
24
CoryD.GarrowayA.MillerJ. (1990). Applications of spin transport as a probe of local geometry.Polymer Prepr.31149–150.
25
CoutuJ.-P.ChenJ. J.RosasH. D.SalatD. H. (2014). Non-gaussian water diffusion in aging white matter.Neurobiol. Aging351412–1421. 10.1016/j.neurobiolaging.2013.12.001
26
DatarA.AmeeramjaJ.BhatA.SrivastavaR.MishraA.BernalR.et al (2019). The roles of microtubules and membrane tension in axonal beading, retraction, and atrophy.Biophys. J.117880–891. 10.1016/j.bpj.2019.07.046
27
DeanD. C.SojkovaJ.HurleyS.KecskemetiS.OkonkwoO.BendlinB. B.et al (2016). Alterations of Myelin Content in Parkinson’s disease: a cross-sectional neuroimaging study.PLoS One11:e0163774. 10.1371/journal.pone.0163774
28
DouaudG.JbabdiS.BehrensT. E. J.MenkeR. A.GassA.MonschA. 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.Neuroimage55880–890. 10.1016/j.neuroimage.2010.12.008
29
DuncanG. W.FirbankM. J.YarnallA. J.KhooT. K.BrooksD. J.BarkerR. A.et al (2016). Gray and white matter imaging: a biomarker for cognitive impairment in early Parkinson’s disease?Mov. Disord.31103–110. 10.1002/mds.26312
30
ErikssonS.LasicS.TopgaardD. (2013). Isotropic diffusion weighting in PGSE NMR by magic-angle spinning of the q-vector.J. Magn. Reson.22613–18. 10.1016/j.jmr.2012.10.015
31
FazekasF.ChawlukJ. B.AlaviA.HurtigH. I.ZimmermanR. A. (1987). MR signal abnormalities at 1.5 T in Alzheimer’s dementia and normal aging.Am. J. Roentgenol.149351–356. 10.2214/ajr.149.2.351
32
FereshtehnejadS.-M.ZeighamiY.DagherA.PostumaR. B. (2017). Clinical criteria for subtyping Parkinson’s disease: biomarkers and longitudinal progression.Brain1401959–1976. 10.1093/brain/awx118
33
FieremansE.BurcawL. M.LeeH. H.LemberskiyG.VeraartJ.NovikovD. S. (2016). In vivo observation and biophysical interpretation of time-dependent diffusion in human white matter.Neuroimage129414–427. 10.1016/j.neuroimage.2016.01.018
34
GilliesG. E.PienaarI. S.VohraS.QamhawiZ. (2014). Sex differences in Parkinson’s disease.Front. Neuroendocrinol.35370–384. 10.1016/j.yfrne.2014.02.002
35
GoetzC. G.TilleyB. C.ShaftmanS. R.StebbinsG. T.FahnS.Martinez-MartinP.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.232129–2170. 10.1002/mds.22340
36
GoveasJ.O’DwyerL.MascalchiM.CosottiniM.DiciottiS.De SantisS.et al (2015). Diffusion-MRI in neurodegenerative disorders.Magn. Reson. Imaging33853–876. 10.1016/j.mri.2015.04.006
37
GuerreriM.PalomboM.CaporaleA.FasanoF.MacalusoE.BozzaliM.et al (2019). Age-related microstructural and physiological changes in normal brain measured by MRI γ-metrics derived from anomalous diffusion signal representation.Neuroimage188654–667. 10.1016/j.neuroimage.2018.12.044
38
HairJ.BlackW. C.BabinB. J.AndersonR. E. (2010). Multivariate Data Analysis, 7th Edn. Upper Saddle River, NJ: Pearson Education.
39
HansenB.LundT. E.SangillR.JespersenS. N. (2013). Experimentally and computationally fast method for estimation of a mean kurtosis.Magn. Reson. Med.691754–1760. 10.1002/mrm.24743
40
HenriquesR. N.JespersenS. N.ShemeshN. (2019). Microscopic anisotropy misestimation in spherical-mean single diffusion encoding MRI.Magn. Reson. Med.813245–3261. 10.1002/mrm.27606
41
HenriquesR. N.JespersenS. N.ShemeshN. (2020). Correlation tensor magnetic resonance imaging.Neuroimage211:116605. 10.1016/j.neuroimage.2020.116605
42
HoehnM. M.YahrM. D. (1967). Parkinsonism: onset, progression, and mortality.Neurology17427–442. 10.1212/WNL.17.5.427
43
IanuşA.JespersenS. N.Serradas DuarteT.AlexanderD. C.DrobnjakI.ShemeshN. (2018). Accurate estimation of microscopic diffusion anisotropy and its time dependence in the mouse brain.Neuroimage183934–949. 10.1016/j.neuroimage.2018.08.034
44
IkenouchiY.KamagataK.AndicaC.HatanoT.OgawaT.Takeshige-AmanoH.et al (2020). Evaluation of white matter microstructure in patients with Parkinson’s disease using microscopic fractional anisotropy.Neuroradiology62197–203. 10.1007/s00234-019-02301-1
45
JelescuI. O.BuddeM. D. (2017). Design and validation of diffusion MRI models of white matter.Front. Phys.28:61. 10.3389/fphy.2017.00061
46
JenkinsonM.BeckmannC. F.BehrensT. E. J.WoolrichM. W.SmithS. M. (2012). FSL.Neuroimage62782–790. 10.1016/j.neuroimage.2011.09.015
47
JensenJ. H.HelpernJ. A.RamaniA.LuH.KaczynskiK. (2005). Diffusional kurtosis imaging: the quantification of non-Gaussian water diffusion by means of magnetic resonance imaging.Magn. Reson. Med.531432–1440. 10.1002/mrm.20508
48
JespersenS. N. (2012). Equivalence of double and single wave vector diffusion contrast at low diffusion weighting.NMR Biomed.25813–818. 10.1002/nbm.1808
49
JespersenS. N.LundellH.SønderbyC. K.DyrbyT. B. (2013). Orientationally invariant metrics of apparent compartment eccentricity from double pulsed field gradient diffusion experiments.NMR Biomed.261647–1662. 10.1002/nbm.2999
50
JianB.VemuriB. C.ÖzarslanE.CarneyP. R.MareciT. H. (2007). A novel tensor distribution model for the diffusion-weighted MR signal.Neuroimage37164–176. 10.1016/j.neuroimage.2007.03.074
51
KadenE.KruggelF.AlexanderD. C. (2016). Quantitative mapping of the per-axon diffusion coefficients in brain white matter.Magn. Reson. Med.751752–1763. 10.1002/mrm.25734
52
KamagataK.TomiyamaH.HatanoT.MotoiY.AbeO.ShimojiK.et al (2014). A preliminary diffusional kurtosis imaging study of Parkinson disease: comparison with conventional diffusion tensor imaging.Neuroradiology56251–258. 10.1007/s00234-014-1327-1
53
KamagataK.TomiyamaH.MotoiY.KanoM.AbeO.ItoK.et al (2013). Diffusional kurtosis imaging of cingulate fibers in Parkinson disease: comparison with conventional diffusion tensor imaging.Magn. Reson. Imaging311501–1506. 10.1016/j.mri.2013.06.009
54
KanazawaT.AdachiE.OrimoS.NakamuraA.MizusawaH.UchiharaT. (2012). Pale neurites, premature α-synuclein aggregates with centripetal extension from axon collaterals.Brain Pathol.2267–78. 10.1111/j.1750-3639.2011.00509.x
55
KellnerE.DhitalB.KiselevV. G.ReisertM. (2016). Gibbs-ringing artifact removal based on local subvoxel-shifts.Magn. Reson. Med.761574–1581. 10.1002/mrm.26054
56
KerkeläL.HenriquesR. N.HallM. G.ClarkC. A.ShemeshN.HallM. G.et al (2020). Validation and noise robustness assessment of microscopic anisotropy estimation with clinically feasible double diffusion encoding MRI.Magn. Reson. Med.831698–1710. 10.1002/mrm.28048
57
KiselevV. G. (2013). “The cumulant expansion: an overarching mathematical framework for understanding diffusion NMR,” in Diffusion MRI, ed.JonesD. K. (Oxford: Oxford University Press). 10.1093/med/9780195369779.003.0010
58
KomloshM. E.HorkayF.FreidlinR. Z.NevoU.AssafY.BasserP. 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.18938–45. 10.1016/j.jmr.2007.07.003
59
KomloshM. E.LizakM. J.HorkayF.FreidlinR. Z.BasserP. J. (2008). Observation of microscopic diffusion anisotropy in the spinal cord using double-pulsed gradient spin echo MRI.Magn. Reson. Med.59803–809. 10.1002/mrm.21528
60
KurowskaZ.KordowerJ. H.StoesslA. J.BurkeR. E.BrundinP.YueZ.et al (2016). Is axonal degeneration a key early event in Parkinson’s disease?J. Parkinsons Dis.6703–707. 10.3233/JPD-160881
61
LampinenB.SzczepankiewiczF.NovénM.van WestenD.HanssonO.EnglundE.et al (2019). Searching for the neurite density with diffusion MRI: challenges for biophysical modeling.Hum. Brain Mapp.402529–2545. 10.1002/hbm.24542
62
LasičS.SzczepankiewiczF.ErikssonS.NilssonM.TopgaardD. (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. 10.3389/fphy.2014.00011
63
LawrenzM.BrassenS.FinsterbuschJ. (2015). Microscopic diffusion anisotropy in the human brain: reproducibility, normal values, and comparison with the fractional anisotropy.Neuroimage109283–297. 10.1016/j.neuroimage.2015.01.025
64
LawrenzM.BrassenS.FinsterbuschJ. (2016). Microscopic diffusion anisotropy in the human brain: age-related changes.Neuroimage141313–325. 10.1016/j.neuroimage.2016.07.031
65
LawrenzM.FinsterbuschJ. (2013). Double-wave-vector diffusion-weighted imaging reveals microscopic diffusion anisotropy in the living human brain.Magn. Reson. Med.691072–1082. 10.1002/mrm.24347
66
LawrenzM.KochM. A.FinsterbuschJ. (2010). A tensor model and measures of microscopic anisotropy for double-wave-vector diffusion-weighting experiments with long mixing times.J. Magn. Reson.20243–56. 10.1016/j.jmr.2009.09.015
67
LeeH.-H.PapaioannouA.KimS.-L.NovikovD. S.FieremansE. (2020). A time-dependent diffusion MRI signature of axon caliber variations and beading.Commun. Biol.3:354. 10.1038/s42003-020-1050-x
68
LevyG. (2007). The relationship of Parkinson disease with aging.Arch. Neurol.641242–1246. 10.1001/archneur.64.9.1242
69
LundellH.NilssonM.DyrbyT. B.ParkerG. J. M.CristinacceP. L. H.ZhouF.-L.et al (2019). Multidimensional diffusion MRI with spectrally modulated gradients reveals unprecedented microstructural detail.Sci. Rep.9:9026.
70
MaddenD. J.BennettI. J.BurzynskaA.PotterG. G.ChenN. K.et al (2012). Diffusion tensor imaging of cerebral white matter integrity in cognitive aging.Biochim. Biophys. Acta1822386–400. 10.1016/j.bbadis.2011.08.003
71
McDonaldC.NewtonJ. L.BurnD. J. (2016). Orthostatic hypotension and cognitive impairment in Parkinson’s disease: causation or association?Mov. Disord.31937–946. 10.1002/mds.26632
72
MitraP. P. (1995). Multiple wave-vector extensions of the NMR pulsed-field-gradient spin-echo diffusion measurement.Phys. Rev. B Condens. Matter5115074–15078. 10.1103/PhysRevB.51.15074
73
Molina-RomeroM.MenzeB. H.PedroA. G.CzischM.SämannP. G.JonesD. 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.802155–2172. 10.1002/mrm.27181
74
MoriS.Van ZijlP. C. M. (1995). Diffusion weighting by the trace of the diffusion tensor within a single scan.Magn. Reson. Med.3341–52. 10.1002/mrm.1910330107
75
NilssonM.LättJ.Van WestenD.BrockstedtS.LasičS.StåhlbergF.et al (2013). Noninvasive mapping of water diffusional exchange in the human brain using filter-exchange imaging.Magn. Reson. Med.691573–1581. 10.1002/mrm.24395
76
NilssonM.SzczepankiewiczF.LampinenB.AhlgrenA.de Almeida MartinsJ. P.LasicS.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.
77
NilssonM.SzczepankiewiczF.BrabecJ.TaylorM.WestinC.-F. C.GolbyA.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.83608–620. 10.1002/mrm.27959
78
NovikovD. S.FieremansE.JespersenS. N.KiselevV. G.MriD.BirdT. A. (2019). Quantifying brain microstructure with diffusion MRI: theory and parameter estimation.NMR Biomed.32:e3998. 10.1002/nbm.3998
79
NovikovD. S.JensenJ. H.HelpernJ. A.FieremansE. (2014). Revealing mesoscopic structural universality with diffusion.Proc. Natl. Acad. Sci. U.S.A.1115088–5093. 10.1073/pnas.1316944111
80
ÖzarslanE.BasserP. J. (2008). Microscopic anisotropy revealed by NMR double pulsed field gradient experiments with arbitrary timing parameters.J. Chem. Phys.128154511. 10.1063/1.2905765
81
ÖzarslanE.YolcuC.HerberthsonM.WestinC.-F.KnutssonH. (2017). Effective potential for magnetic resonance measurements of restricted diffusion.Front. Phys.5:68. 10.3389/fphy.2017.00068
82
PalomboM.LigneulC.Hernandez-GarzonE.ValetteJ. (2018). Can we detect the effect of spines and leaflets on the diffusion of brain intracellular metabolites?Neuroimage182283–293. 10.1016/j.neuroimage.2017.05.003
83
PasternakO.SochenN.GurY.IntratorN.AssafY. (2009). Free water elimination and mapping from diffusion MRI.Magn. Reson. Med.62717–730. 10.1002/mrm.22055
84
PlanettaP. J.SchulzeE. T.GearyE. K.CorcosD. M.GoldmanJ. G.LittleD. M.et al (2013). Thalamic projection fiber integrity in de novo Parkinson Disease.Am. J. Neuroradiol.3474–79. 10.3174/ajnr.A3178
85
PoeweW.SeppiK.TannerC. M.HallidayG. M.BrundinP.VolkmannJ.et al (2017). Parkinson disease.Nat. Rev. Dis. Primers3:17013. 10.1038/nrdp.2017.13
86
PortnoyS.FlintJ. J.BlackbandS. J.StaniszG. 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.691131–1145. 10.1002/mrm.24325
87
PostumaR. B.BergD.SternM.PoeweW.OlanowC. W.OertelW.et al (2015). MDS clinical diagnostic criteria for Parkinson’s disease.Mov. Disord.301591–1601. 10.1002/mds.26424
88
PozorskiV.OhJ. M.AdluruN.MerluzziA. P.TheisenF.OkonkwoO.et al (2018). Longitudinal white matter microstructural change in Parkinson’s disease.Hum. Brain Mapp.394150–4161. 10.1002/hbm.24239
89
PrakashB. D.SitohY.-Y.TanL. C. S.AuW. L. (2012). Asymmetrical diffusion tensor imaging indices of the rostral substantia nigra in Parkinson’s disease.Parkinsonism Relat. Disord.181029–1033. 10.1016/j.parkreldis.2012.05.021
90
RektorI.SvátkováA.VojtíšekL.ZikmundováI.VaníčekJ.KirályA.et al (2018). White matter alterations in Parkinson’s disease with normal cognition precede grey matter atrophy.PLoS One13:e0187939. 10.1371/journal.pone.0187939
91
ReymbautA.MezzaniP.Almeida MartinsJ. P.TopgaardD. (2020). Accuracy and precision of statistical descriptors obtained from multidimensional diffusion signal inversion algorithms.NMR Biomed.17:e4267. 10.1002/nbm.4267
92
RiedererP.JellingerK. A.KolberP.HippG.Sian-HülsmannJ.KrügerR. (2018). Lateralisation in Parkinson disease.Cell Tissue Res.373297–312. 10.1007/s00441-018-2832-z
93
SalvadoresN.SanhuezaM.ManqueP.CourtF. A. (2017). Axonal degeneration during aging and its functional role in neurodegenerative disorders.Front. Neurosci.11:451. 10.3389/fnins.2017.00451
94
Sanjari MoghaddamH.DolatshahiM.MohebiF.AarabiM. H. (2020). Structural white matter alterations as compensatory mechanisms in Parkinson’s disease: a systematic review of diffusion tensor imaging studies.J. Neurosci. Res.981398–1416. 10.1002/jnr.24617
95
ScherflerC.SeppiK.MairK. J.DonnemillerE.VirgoliniI.WenningG. K.et al (2012). Left hemispheric predominance of nigrostriatal dysfunction in Parkinson’s disease.Brain1353348–3354. 10.1093/brain/aws253
96
ShahimP.HolleranL.KimJ. H.BrodyD. L. (2017). Test-retest reliability of high spatial resolution diffusion tensor and diffusion kurtosis imaging.Sci. Rep.7:11141. 10.1038/s41598-017-11747-3
97
ShemeshN.BarazanyD.SadanO.BarL.ZurY.BarhumY.et al (2012). Mapping apparent eccentricity and residual ensemble anisotropy in the gray matter using angular double-pulsed-field-gradient MRI.Magn. Reson. Med.68794–806. 10.1002/mrm.23300
98
ShemeshN.CohenY. (2011). Microscopic and compartment shape anisotropies in gray and white matter revealed by angular bipolar double-PFG MR.Magn. Reson. Med.651216–1227. 10.1002/mrm.22738
99
ShemeshN.ÖzarslanE.BasserP. J.CohenY. (2009). Measuring small compartmental dimensions with low-q angular double-PGSE NMR: the effect of experimental parameters on signal decay.J. Magn. Reson.19815–23. 10.1016/j.jmr.2009.01.004
100
SimuniT.Caspell-GarciaC.CoffeyC.LaschS.TannerC.MarekK. (2016). How stable are Parkinson’s disease subtypes in de novo patients: analysis of the PPMI cohort?Parkinsonism Relat. Disord.2862–67. 10.1016/j.parkreldis.2016.04.027
101
SmithS. M.JenkinsonM.Johansen-BergH.RueckertD.NicholsT. E.MackayC. E.et al (2006). Tract-based spatial statistics: voxelwise analysis of multi-subject diffusion data.Neuroimage311487–1505. 10.1016/j.neuroimage.2006.02.024
102
SmithS. M.NicholsT. E. (2009). Threshold-free cluster enhancement: addressing problems of smoothing, threshold dependence and localisation in cluster inference.Neuroimage4483–98. 10.1016/j.neuroimage.2008.03.061
103
SurovaY.LampinenB.NilssonM.LättJ.HallS.WidnerH.et al (2016). Alterations of diffusion kurtosis and neurite density measures in deep grey matter and white matter in Parkinson’s disease.PLoS One11:e0157755. 10.1371/journal.pone.0157755
104
SurovaY.NilssonM.LampinenB.LättJ.HallS.WidnerH.et al (2018). Alteration of putaminal fractional anisotropy in Parkinson’s disease: a longitudinal diffusion kurtosis imaging study.Neuroradiology60247–254. 10.1007/s00234-017-1971-3
105
SzczepankiewiczF.LasiS.WestenD.Van, SundgrenP. C.EnglundE.et al (2015). Quantification of microscopic diffusion anisotropy disentangles effects of orientation dispersion from microstructure: applications in healthy volunteers and in brain tumors.Neuroimage104241–252. 10.1016/j.neuroimage.2014.09.057
106
SzczepankiewiczF.van WestenD.EnglundE.WestinC. F.StåhlbergF.LättJ.et al (2016). The link between diffusion MRI and tumor heterogeneity: Mapping cell eccentricity and density by diffusional variance decomposition (DIVIDE).Neuroimage142522–532. 10.1016/j.neuroimage.2016.07.038
107
SzczepankiewiczF.WestinC.NilssonM. (2019). Maxwell-compensated design of asymmetric gradient waveforms for tensor-valued diffusion encoding.Magn. Reson. Med.821424–1437. 10.1002/mrm.27828
108
TagliaferroP.KarevaT.OoT. F.YaryginaO.KholodilovN.BurkeR. E. (2015). An early axonopathy in a hLRRK2(R1441G) transgenic model of Parkinson disease.Neurobiol. Dis.82359–371. 10.1016/j.nbd.2015.07.009
109
TopgaardD. (2017). Multidimensional diffusion MRI.J. Magn. Reson.27598–113. 10.1016/j.jmr.2016.12.007
110
TournierJ.-D.SmithR.RaffeltD.TabbaraR.DhollanderT.PietschM.et al (2019). MRtrix3: a fast, flexible and open software framework for medical image processing and visualisation.Neuroimage202:116137. 10.1016/j.neuroimage.2019.116137
111
TustisonN. J.AvantsB. B.CookP. A.ZhengYEganA.YushkevichP. A.et al (2010). N4ITK: improved N3 bias correction.IEEE Trans. Med. Imaging291310–1320. 10.1109/TMI.2010.2046908
112
VeraartJ.FieremansE.NovikovD. S. (2016). Diffusion MRI noise mapping using random matrix theory.Magn. Reson. Med.761582–1593. 10.1002/mrm.26059
113
VeselýB.AntoniniA.RektorI. (2016). The contribution of white matter lesions to Parkinson’s disease motor and gait symptoms: a critical review of the literature.J. Neural Transm.123241–250. 10.1007/s00702-015-1470-9
114
WangJ.-J.LinW.-Y.LuC.WengY.-H.NgS.-H.WangC.-H.et al (2011). Parkinson disease: diagnostic utility of diffusion kurtosis imaging.Radiology261210–217. 10.1148/radiol.11102277
115
WestinC. F.KnutssonH.PasternakO.SzczepankiewiczF.ÖzarslanE.van WestenD.et al (2016). Q-space trajectory imaging for multidimensional diffusion MRI of the human brain.Neuroimage135345–362. 10.1016/j.neuroimage.2016.02.039
116
WinklerA. M.RidgwayG. R.DouaudG.NicholsT. E.SmithS. M. (2016). Faster permutation inference in brain imaging.Neuroimage141502–516. 10.1016/j.neuroimage.2016.05.068
117
WinklerA. M.RidgwayG. R.WebsterM. A.SmithS. M.NicholsT. E. (2014). Permutation inference for the general linear model.Neuroimage92381–397. 10.1016/j.neuroimage.2014.01.060
118
YadavS. K.KathiresanN.MohanS.VasileiouG.SinghA.KauraD.et al (2016). Gender-based analysis of cortical thickness and structural connectivity in Parkinson’s disease.J. Neurol.2632308–2318. 10.1007/s00415-016-8265-2
119
YangG.TianQ.LeuzeC.WintermarkM.McNabJ. A. (2018). Double diffusion encoding MRI for the clinic.Magn. Reson. Med.80507–520. 10.1002/mrm.27043
120
YolcuC.MemiçM.ŞimşekK.WestinC.-F.ÖzarslanE. (2016). NMR signal for particles diffusing under potentials: from path integrals and numerical methods to a model of diffusion anisotropy.Phys. Rev. E93:052602. 10.1103/PhysRevE.93.052602
121
ZhangY.WuI.-W.BuckleyS.CoffeyC. S.FosterE.MendickS.et al (2015). Diffusion tensor imaging of the nigrostriatal fibers in Parkinson’s disease.Mov. Disord.301229–1236. 10.1002/mds.26251
122
ZhangY.WuI.-W.TosunD.FosterE.SchuffN. (2016). Progression of regional microstructural degeneration in Parkinson’s disease: a multicenter diffusion tensor imaging study.PLoS One11:e0165540. 10.1371/journal.pone.0165540
Summary
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
Volume
14 - 2020
Edited by
Yu-Chien Wu, Indiana University Bloomington, United States
Reviewed by
Evren Özarslan, Linköping University, Sweden; Matthew ManHin Cheung, CUHK Medical Centre, Hong Kong
Updates
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
This article was submitted to Brain Imaging Methods, a section of the journal Frontiers in Neuroscience
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.