Diffusion–relaxation correlation multidimensional MRI (MD-MRI) replaces voxel-averaged diffusion tensor quantities and R1 and R2 relaxation rates with their multidimensional distributions, enabling the selective extraction and mapping of specific diffusion–relaxation spectral ranges that correspond to different cellular features. This approach has the potential of achieving high sensitivity and specificity in detecting subtle changes that would otherwise be averaged out. Here, the whole brain characterization of MD-MRI distributions and derived parameters is presented and the intrascanner test–retest reliability, repeatability, and reproducibility are evaluated to promote the further development of these quantities as neuroimaging biomarkers. We compared white matter tracts and cortical and subcortical gray matter regions, revealing notable variations in their diffusion–relaxation profiles, indicative of unique microscopic morphological characteristics. We found that the reliability and repeatability of MD-MRI-derived diffusion and relaxation mean parameters were comparable with values expected in conventional diffusion tensor imaging and relaxometry studies. Importantly, the estimated signal fractions of intravoxel spectral components in the MD-MRI distribution, corresponding to white matter, gray matter, and cerebrospinal fluid, were found to be reproducible. This underscores the viability of employing a spectral analysis approach to MD-MRI data. Our results show that a clinically feasible MD-MRI protocol can reliably deliver information of the rich structural and chemical variety that exists within each imaging voxel, creating potential for new MRI biomarkers with enhanced sensitivity and specificity.

Diffusion–relaxation multidimensional MRI (MD-MRI) acquisition integrates meso- and microstructural probes (Callaghan, 2011) with chemical composition sensitivity (Tofts, 2003). While relaxation values can be transformed into estimates of, for example, myelin fraction (Mackay et al., 1994; Vasilescu et al., 1978) and diffusion tensors into estimates of, for example, general orientation of nerve fiber bundles (Basser et al., 1994) within individual millimeter-scale image voxels, achieving a more precise interpretation is hindered by the complex reality that each voxel comprises numerous cells with varying sizes, shapes, and orientations (Livet et al., 2007). This heterogeneity creates considerable uncertainty when establishing connections between the observable metrics and the specific microscopic details we are often interested to discern.

While the classical diffusion tensor imaging (DTI) model (Basser et al., 1994) has proven sensitive to probing brain microstructure in a range of scenarios and applications (Soares et al., 2013), it is encumbered by inherent limitations that impede its applicability and the precision of its interpretations (Novikov et al., 2018). To enhance specificity, various multicomponent biophysical models have been developed and successfully applied to study brain white matter (WM) (Assaf & Basser, 2005; Benjamini et al., 2016; Stanisz et al., 1997; Veraart et al., 2020; Zhang et al., 2012). While these models offer advantages in specificity and interpretability, they rely on strict assumptions and incorporate fixed parameters to reduce fitting ambiguities (Novikov et al., 2018), which may not be valid in disease states or even in heterogeneous gray matter (GM) regions (Jespersen et al., 2019; Lampinen et al., 2019; Novikov et al., 2018; Veraart et al., 2019).

An alternative approach is to avoid a biophysical model and apply a data-driven estimation of distributions of quantitative metrics. This strategy applies a mathematical representation that balances physical accuracy, simplicity, and usefulness in addressing scientific questions without encouraging overinterpretation. First introduced in 1986 to obtain T2 distributions from muscle tissue (Kroeker & Mark Henkelman, 1986), this approach allows to intuitively expand the dimensionality of the data by simultaneously encoding both diffusion and relaxation. Over the years, it has been used to resolve individual diffusion and relaxation properties of distinct water populations, and their correlations (Hürlimann & Venkataramanan, 2002; Silva et al., 2002). The obtained information has been shown invaluable to the study of tissue microstructure and brain connectivity (de Almeida Martins et al., 2020, 2021; Kundu et al., 2023; Peled et al., 1999; Reymbaut et al., 2021; Stanisz & Henkelman, 1998; Yon et al., 2020) as well as pathology (Benjamini et al., 2020, 2022; Kim et al., 2017; Slator et al., 2021; Wei et al., 2022).

An important aspect of translational motion of particles in biological tissue is its time dependency (Callaghan, 2011; Price, 2009). This ensemble translational motion of particles can be analyzed in two ways: in the time domain using metrics like mean-square displacements (Einstein, 1905), apparent diffusivities (Woessner, 1963), and velocity autocorrelation functions (Uhlenbeck & Ornstein, 1930), or in the frequency domain using tensor-valued diffusion spectra D(ω), derived via Fourier transformation (Callaghan & Stepišnik, 1995; Stepišnik, 1981). Deciding between time- and frequency-domain analysis for a specific application is challenging, but recent studies suggest a preference for the frequency-domain approach in in vivo human research (Arbabi et al., 2020; Hamilton et al., 2024; Hennel et al., 2021; Maekawa et al., 2020; Tan et al., 2020; Tétreault et al., 2020; Xu et al., 2020, 2021).

The ω-dependence of D(ω) is often approximated as Lorentzian, reflecting an exponential velocity autocorrelation function (Stepišnik, 1998; Uhlenbeck & Ornstein, 1930), but more precise expressions exist for simple pore shapes such as parallel planes, cylinders, spheres (Stepišnik, 1993), and for random permeable barriers (Novikov et al., 2014). The latter model, which captures structural disorder without assuming specific geometry, has been increasingly used to analyze in vivo human brain data in recent studies (Arbabi et al., 2020; Fieremans et al., 2016; Lee et al., 2020). Further advancing the field, the employment of diffusion acquisition strategies that incorporate tensor-valued diffusion encoding with free gradient waveforms (Westin et al., 2014) contains information about higher order diffusion propagator covariances not present in linear diffusion encoding. The ω-dependence of D(ω) components can be accessed by using modulated gradients g(t) (Stepišnik, 1981), which produce tensor-valued encoding spectra b(ω) that concentrate spectral power on specific tensor elements and ω values (Lasič et al., 2022; Lundell & Lasič, 2020). The diffusion tensor distribution (DTD) provides a framework for separating and decoupling distinct water pools within complex and heterogeneous systems (Basser & Pajevic, 2003). Parametric DTD models resolve the mathematical instability in the estimation process; however, they rely on preassumed models for the underlining distribution (Jian et al., 2007; Magdoom et al., 2021). An alternative approach is the nonparametric DTD, which can be reliably estimated from tensor-valued diffusion encoding data, and allows for resolving WM, GM, and CSF from isotropic-anisotropic correlations (de Almeida Martins & Topgaard, 2016; Topgaard, 2019). Recently, a nonparametric tensor-valued Lorentzian diffusion spectra framework has been proposed as a flexible representation to capture MRI signal responses from diverse model systems and biological tissues (Narvaez et al., 2024), distinguishing between restricted and anisotropic diffusion effects. It was shown to reproduce the synthetic data generated by the random permeable barrier model (Novikov et al., 2014).

Demonstrations of nonparametric diffusion–relaxation correlations of MD-MRI included D- R1- R2 in vivo in humans (Martin et al., 2021), and D(ω)- R1- R2 ex vivo in rat (Narvaez et al., 2022). The latter was recently translated to an efficient and sparse in vivo frequency-dependent MD-MRI acquisition protocol that provides whole brain coverage at 2-mm isotropic resolution (Johnson et al., 2024). This proof-of-concept study tested the feasibility of estimating voxel-wise distributions of frequency-dependent diffusion tensors, D(ω), and longitudinal and transverse relaxation rates, R1 and R2, from this sparsely acquired MD-MRI data. The rich information contained within the high-dimensional D(ω)- R1- R2 distributions can be distilled into first (means) and second (variances and covariances) moment statistics of, for example, frequency-dependent isotropic diffusivity and diffusion anisotropy, R1, and R2, in each voxel. Further, having voxel-wise distributions provides access to distinct water populations within a voxel, and, therefore, enables improved specificity toward WM, GM, and cerebrospinal fluid (CSF) (de Almeida Martins et al., 2020). In turn, this unique intravoxel information allows to map microstructure-specific diffusion and relaxation properties without relying on biophysical assumptions or models.

In recent studies, the specificity of estimated intravoxel components using MD-MRI was evaluated in both physical (Benjamini & Basser, 2016; de Almeida Martins & Topgaard, 2016, 2018) and synthetic (Reymbaut et al., 2020) phantoms, and via comparisons with microscopy of fixed tissue samples from the rat and ferret brain (Benjamini & Basser, 2017; Narvaez et al., 2022). The applicability of this framework for in vivo human MRI was confirmed by comparing MD-MRI-derived values in the human brain with those reported in the literature (de Almeida Martins et al., 2021; Martin et al., 2021; Naranjo et al., 2021). In this context, our aims are to present a characterization of D(ω)- R1- R2 estimates throughout the human brain and to evaluate their intrascanner test–retest reliability and repeatability. These efforts contribute to the ongoing development of the MD-MRI framework as a neuroimaging modality.

2.1 Participants

Ten healthy participants (mean age 48 years, standard deviation 14.4 years; 4 women) were each scanned twice, a few weeks apart (i.e., a total of 20 scans). Participants were systematically drawn from ongoing healthy cohorts of the National Institute on Drug Abuse (NIDA). Experimental procedures were performed in compliance with our local Institutional Review Board, and participants provided written informed consent. Prior to each scan, NIDA clinical and nursing units conducted COVID-19 testing, urine drug tests, a physical examination, and a questionnaire on pre-existing conditions and daily habits. Exclusion criteria included major medical illness or current medication use, a history of neurological or psychiatric disorders or substance abuse.

2.2 Data acquisition

The MRI data acquisition sequences and procedures used were the same as in a previous study (Johnson et al., 2024). Data were acquired using a 3T scanner (MAGNETOM Prisma, Siemens Healthcare AG, Erlangen, Germany) with a 32-channel head coil. MD-MRI data were acquired with 2-mm isotropic voxel size using a single-shot spin-echo EPI sequence modified for tensor-valued diffusion encoding with free gradient waveforms (Martin et al., 2021; Wetscherek et al., 2015). The acquisition parameters were set as follows: FOV = 228 × 228 × 110 mm3, voxel size = 2 × 2 × 2 mm3, acquisition bandwidth = 1512 Hz/pixel, in-plane acceleration factor 2 using GRAPPA reconstruction with 24 reference lines, effective echo spacing of 0.8 ms, phase-partial Fourier factor of 0.75, and axial slice orientation.

The acquisition protocol was designed based on previously described heuristic principles (Martin et al., 2021) used to achieve an efficient and sparse MD-MRI acquisition (Johnson et al., 2024). Briefly, in addition to a b=0 ms/µm2 volume, numerically optimized gradient waveforms (Sjölund et al., 2015) that produce linear, planar, and spherical b-tensors were employed with b-tensor magnitude, b, ranging between 0.1 and 3 ms/µm2, and b-tensor normalized anisotropy, bΔ, of – 0.5 (planar), 0 (spherical), and 1 (linear). The spectral content of the employed diffusion gradient waveforms was in the range of 6.6–21 Hz centroid frequencies, ωcent/2π. We note that each b-tensor contains a range of frequencies and that the centroid frequency is used only as its representative value. All data inversion and metric calculations use frequency as such and do not rely on reducing the spectral content of b(ω) to a single number. The datasets were acquired with a single phase encoding direction (anterior to posterior, AP), and an additional b=0 ms/µm2 volume with reversed phase encoding direction (PA). Sensitivity to R1 and R2 was achieved by acquiring data with different combinations of repetition times, TR = (0.62, 1.75, 3.5, 5, 7, 7.6) s and echo times, TE = (40, 63, 83, 150) ms. The number of concatenations and preparation scans was increased to allow values of TR below 5 s. A total of 139 individual measurements were recorded over 40 min. More details of the acquisition protocol are available in Johnson et al. (2024), along with a detailed summary of the acquisition protocol shown in Supplementary Figure 1. In addition to MD-MRI, fat-suppressed T1-weighted MPRAGE images with TE=3.42 ms, TR=1900 ms, and 1 mm isotropic voxel size were also acquired as structural images. The structural data were utilized in image registrations and region of interest (ROI) segmentations.

2.3 MD-MRI data preprocessing

We followed recent recommendations for denoising and preprocessing of sparse MD-MRI data (Johnson et al., 2024). In short, all data volumes were first concatenated and underwent denoising with the MPPCA technique (Veraart et al., 2016). Then, Gibbs ringing correction (Kellner et al., 2016) for partial k-space acquisitions (Lee et al., 2021) was performed, followed by motion and eddy currents distortion corrections using TORTOISE’s DIFFPREP (Rohde et al., 2004) with a physically based parsimonious quadratic transformation model and a normalized mutual information metric. For susceptibility distortion correction, a T1-weighted image was initially converted into a T2-weighted image with b=0 s/mm2 contrast (Schilling et al., 2019), which was fed into the DR-BUDDI software (Irfanoglu et al., 2015) for AP-PA distortion correction. The final preprocessed images from the MD-MRI sequence were output with a single interpolation in the space of the anatomical T1-weighted MPRAGE image at native MD-MRI in-plane voxel size (2 mm isotropic).

2.4 Estimation of MD-MRI parameters

Solving the joint D(ω)- R1- R2 probability density distribution from the MD-MRI measurements is an ill-conditioned problem. To solve that inversion problem, the estimation of the MD-MRI parameters used in this study followed the procedures described in Narvaez et al. (2022). The multidimensional diffusion MRI toolbox for MATLAB (Nilsson et al., 2018) was used to perform a Monte Carlo inversion (de Almeida Martins & Topgaard, 2018; Prange & Song, 2009) for the preprocessed MD-MRI data. In the inversion, the MD-MRI signal S(b(ω),TE,TR) is represented as a weighted sum of discrete components c of the joint probability density distribution D(ω)- R1- R2:

where Nc is the number of fitted components, fc is the weight of the estimated component c, : denotes a generalized dot product, and D(ω) is an axially symmetric diffusion tensor that depends on the angular frequency of the diffusion gradient waveform. The diffusion tensors are approximated as axisymmetric Lorentzians parameterized by the zero-frequency axial and radial diffusivities [D, D], azimuthal and polar angles [θ, ϕ], high frequency isotropic diffusivity, D0, axial and radial transition frequencies, [Γ, Γ], along with longitudinal and transversal relaxation rates [R1, R2]. In this study, these parameters were sampled in the ranges 0.05<D,,0<5 µm2/ms, 0<θ<π, 0<ϕ<2π, 0.2<R1<2 s-1, 1<R2<30 s-1, and 0.01<Γ,<10000 s-1. The Monte Carlo inversion used non-negative least squares with quasi-genetic filtering, bootstrapping, and constraints for the MD-MRI parameter estimates (Benjamini, 2020; de Almeida Martins & Topgaard, 2018). As a result, for each imaging voxel, 300 solutions (bootstraps) of the D(ω)- R1- R2 probability density distribution with up to 10 weighted components (for each bootstrap) were estimated. Voxel-wise D(ω)- R1- R2 distributions in the primary analysis space [D, D, θ, ϕ, D0, Γ, Γ, R1, R2] were evaluated at selected values of ω/2π (Narvaez et al., 2024) within the narrow 6.6–21 Hz experimental window, giving a set of ω-dependent distributions in the [D(ω), D(ω), θ, ϕ, R1, R2] space. It should be noted that the highest frequencies reached for b-values of 0.5, 1.5, and 3 ms/µm2 were 21, 15, and 11 Hz, respectively. We further chose to express the axial and radial diffusivities as isotropic diffusivity, Diso, and squared normalized anisotropy DΔ2 (Conturo et al., 1996).

Two-dimensional projections of the full 6D D(ω)- R1- R2 distributions were constructed to facilitate visualization and to allow investigation of certain ROIs. To obtain the distributions in each voxel, the bootstrap solutions first have to be consolidated into a single representative distribution. This was done using the following procedure: (1) grouping D(ω)- R1- R2 components across different bootstraps solutions using k-means clustering over [Diso(ω), DΔ2(ω), θ, ϕ, R1, R2] as features (after normalization according to their respective maximal values) with the L1 distance metric. (2) Once grouped, each cluster is transformed into a single D(ω)- R1- R2 component by taking the weighted mean of each parameter (e.g., Diso), and the median of the weights. This process then allowed to project and map the median weights of the discrete components from the consolidated distributions onto 64 × 64 meshes in the Diso- DΔ2, Diso- R1, and Diso- R2 planes. These voxel-wise 2D projections were then summed over entire ROIs within each scan, normalized, and finally averaged across all 20 scans to produce characteristic distributions.

The estimated MD-MRI parameters assessed for reliability and repeatability in this study included the expectation E[x] (mean) values of isotropic diffusivity Diso (E[Diso]), squared diffusion anisotropy DΔ2 (E[DΔ2]), R1 relaxation rate (E[R1]), and R2 relaxation rate (E[R2]). The effects of restricted diffusion were quantified by a finite difference approximation of the rate of change of the diffusivity metrics with frequency ω/2π within the investigated window (6.6 to 21 Hz), expressed as Δω/2π[x] of Diso and DΔ2 statistics (Δω/2π E[Diso], Δω/2π E[DΔ2], respectively) (Johnson et al., 2024; Narvaez et al., 2022). The two-dimensional diffusivity – diffusion anisotropy space can be used to delineate subvoxel components corresponding to different microstructural environments (Pierpaoli et al., 1996; Topgaard, 2019). Subvoxel components were expressed as bin signal fractions, fbin1, fbin2, and fbin3 (Diso<2.5 µm2/ms and DΔ2>0.25; Diso<2.5 µm2/ms and DΔ2<0.25; and Diso>2.5 µm2/ms and 0<DΔ2<1, respectively), which roughly represent WM, GM, and CSF signal fractions, respectively. Intravoxel information from diffusion–relaxation correlations can be directly imaged by resolving the relaxation properties according to the above-defined bins (de Almeida Martins et al., 2020). We, therefore, included in this study the bin-resolved means of the R1 and R2 relaxation, for example, E[R1]bin1, E[R2]bin3. The expectation values E[x] provide estimates of voxel averages of isotropic diffusivity Diso, squared diffusion anisotropy DΔ2, and relaxation rates R1 and R2, and are the closest analogues of mean diffusivity and fractional anisotropy in diffusion tensor imaging and monoexponential R1 and R2 relaxation rates in relaxation mapping. The frequency dependence parameters (Δω/2π[x]) describe the diffusion spectrum, its rate of change, and how the diffusion parameters depend on the frequency of the diffusion encoding gradients. Specifically, Δω/2πE[Diso] is more commonly known as ΔfADC (Aggarwal et al., 2012) or the diffusion dispersion rate (Xu et al., 2016). The bin fractions fbin1, fbin2, fbin3 are computed from the discrete Diso and DΔ2 components and, thus, depend on their estimation accuracy. All of the above parameters are computed by taking the median over bootstrap replicas of the estimates.

2.5 Assessment of reliability and repeatability

2.5.1 Intraclass correlation coefficient (ICC)

We used the intraclass correlation coefficient (ICC) as a measure of reliability. We used a single-measurement, absolute-agreement, two-way mixed-effects model (ICC(A, 1) in the convention of McGraw and Wong) for the computation (McGraw & Wong, 1996). MATLAB was used for the computation (MATLAB Release R2023a, The MathWorks, Inc., Natick, MA, USA). ICC estimates the proportion of variance due to differences between the participants compared with the total variance of the data. In other words, it is an estimate of the fraction of the total variance that is not due to measurement errors (or errors introduced by the processing pipeline, such as image registration inaccuracy). It characterizes whether differences between the participants can be detected from the noisy measurements. ICC can be improved by either decreasing the measurement error or by increasing the biological variability in the sample, such as by selecting participants with larger differences in age. The ICC ranges from 0 to 1, and while negative values are possible, they are treated as zeros. We use the classification used by Koo and Li (2016) to interpret the ICC values: poor reliability (ICC <0.5), moderate reliability (0.5 ≤ICC <0.75), good reliability (0.75 ≤ICC <0.90), and excellent reliability (ICC ≥0.9).

2.5.2 Within-subject coefficient of variation (CVws)

We used the within-subject coefficient of variation (CVws) as a measure of repeatability (Luque Laguna et al., 2020). It is the ratio of within-subject standard deviation σws to the mean over subjects μ. We estimate CVws as

where MSws is the within-subject mean squares and μ¯ is the mean over participants and repetitions. Note that we take the absolute value of the mean because some of the MD-MRI parameters we use can take negative values. A within-subject coefficient of variation CVws=0 indicates perfect repeatability.

2.6 Voxel-wise analysis

Discounting super-resolution and biophysical modeling-based approaches, voxel-wise analysis enables the most spatially specific analyses of MR images. To facilitate such an analysis, we registered (transformed) the MR images of the different subjects and time points to a common space (a study-specific template). Because the images are aligned, differences in MR signals between subjects (or time points) can be evaluated at each imaging voxel. However, this approach of analysis is sensitive to image registration inaccuracies and imaging noise. Accuracy of the image registration depends on the resolution, contrast, and signal-to-noise ratio of the images, as well as on the degree of anatomical differences between the subjects. We applied spatial filtering after image registration to alleviate the confounds of image registration inaccuracies and measurement noise on the voxel-wise analyses. Although commonly used in voxel-wise test–retest neuroimaging studies (Caceres et al., 2009; Shou et al., 2013; Somandepalli et al., 2015), spatial filtering of images significantly reduces the effective resolution of the resulting data.

2.6.1 Image registration

To enable the voxel-wise analysis, the MD-MRI parameter maps of the different participants were aligned to a common space using a study-specific template. Advanced Normalization Tools (Avants et al., 2011) (ANTs version 2.4.4, http://stnava.github.io/ANTs/, accessed on 14 June 2023) was used to compute and apply the image registrations. High-resolution (1-mm isotropic) T1-weighted MPRAGE images were used to create the study-specific template. First, repetition 1 and repetition 2 of each participant were registered to a mid-way space between the repetitions with an affine transform and averaged to avoid causing interpolation bias that would differ between each repetition (Reuter et al., 2012). Then, a study-specific template was formed from the test–retest average T1-weighted images of the different participants using a combination of affine and symmetric image normalization (SyN) transforms. The T1-weighted template image (1-mm isotropic resolution) was downsampled to match the resolution of the MD-MRI images (2-mm isotropic resolution). Finally, because the MD-MRI images of each participant had already been aligned to their corresponding T1-weighted image during the preprocessing stage (Section 2.3), the MD-MRI parameter maps of each participant and repetition were transformed to the 2-mm resolution template space by combining and applying the image transforms that were computed to make the T1-weighted template. Only one interpolation step was applied (after the application of all the image transforms in succession) to avoid the accumulation of image interpolation errors.

2.6.2 Spatial filtering

When image misregistration occurs, a brain structure will be misaligned between different participants. Spatial filtering can be applied to blend together the surroundings of each voxel, thus increasing the degree of overlap in the brain structure between the participants at the cost of effective image resolution. As a side effect, spatial filtering also increases the signal-to-noise ratio of each voxel. After the transformation of the parameter maps to the template space, spatial filtering was applied to reduce the influence of inaccuracies caused by image registration. The choice of the optimal filtering function is nontrivial. However, as Caceres et al. (2009) demonstrated, ICC can be used to guide the decision because it quantifies how well differences between the subjects can be detected over the overall level of variability. Moderate amounts of smoothing can decrease the overall variability (image misalignment and imaging noise), thus increasing the ICC. However, smoothing decreases the spatial specificity of the parametric maps. If there is a regionally specific difference between subjects, it may vanish when excessive smoothing is applied. We applied Gaussian spatial filtering at different values of full width at half maximum (FWHM) for the kernel. When applying the filter, signal contributions from voxels outside of the brain mask were excluded. Then, we computed ICC for each voxel and the median of the ICC values of all brain voxels. Supplementary Figure 2 demonstrates varying behaviors of median ICC as a function of the amount of smoothing applied, depending on the MD-MRI parameter map. Some, such as E[Diso], increase monotonically with the amount of smoothing applied. Such behavior indicates global differences in the parameter values between the participants. Other parameters, such as E[DΔ2], display an optimum at a moderate amount of smoothing, indicating spatially localized differences between the participants that we should seek to preserve. A value of FWHM = 2 pixels (4 mm) was chosen for the voxel-wise analyses, as that value appears to retain such localized differences in all of the parameter maps, yet offers a significant boost to ICC values by decreasing the overall variability due to image misalignment and noise. Values of smoothing FWHM larger than 2 pixels decreased the median voxel-wise ICC values for some parameter maps, suggesting that the decrease in spatial specificity was outweighing the benefits associated with spatial smoothing.

2.6.3 Brain-level statistics of voxel-wise metrics

After the application of spatial filtering (FWHM = 2 pixels = 4 mm), we computed reliability (ICC) and repeatability (CVws) measures for each brain voxel. Violin plots were used to visualize the distributions of voxel-wise values of ICC and CVws over the whole brain. For some MD-MRI parameter maps, a small number of voxels had very high values of CVws. The resulting distributions of voxel-wise values were very broad yet very sparse, making it difficult to produce plots that represented the correct shapes of the distributions. For those distributions of CVws that contained values higher than 100%, the MATLAB function “histcounts” (with automatically selected bins) was used to iteratively remove the voxels that were contained within histogram bins with a density of less than 0.1% of the total number of voxels. The resulting violin plots were labeled to indicate that they contained outliers that were removed prior to plotting. This outlier removal process was applied only for the violin plots of CVws.

2.7 ROI-based analysis

In contrast to voxel-wise analysis, ROI-based analysis computes statistical measures within a chosen anatomical region. In our analysis, we compute the mean of each MD-MRI parameter map for each ROI. The resulting parameter estimates are thus spatially averaged based on anatomically delineated brain regions. Compared with voxel-wise analysis, ROI-based analysis is more robust as it is less susceptible to imaging noise and image registration inaccuracies. Voxel-wise analysis has gained popularity as an exploratory method, whereas ROI-based analysis is more common for (anatomically or functionally) hypothesis-driven studies.

2.7.1 Whole brain segmentation

The spatially localized atlas network tiles (SLANT) method was used to perform whole brain segmentation (Huo et al., 2019) and obtain cortical and subcortical GM ROIs. Briefly, SLANT employs multiple independent 3D convolutional networks to segment the brain. Each of the networks is only responsible for a particular spatial region, thus the task of each network is simplified to focus on patches from a similar portion of the brain. Within this end-to-end pipeline, affine registration, N4 bias field correction, and intensity normalization are employed to roughly normalize each brain to the same space before segmentation. After each network performs its duty, the segmentation labels are fused together to form labels for 132 anatomical regions, primarily from subcortical and cortical areas, based on the BrainCOLOR protocol (https://mindboggle.info/braincolor/). The SLANT framework has shown high intra- and interscan protocol reproducibility (Xiong et al., 2019). It was performed using the 1-mm-resolution T1-weighted MPRAGE images to obtain GM segmentations in the native space of each participant and repetition. Finally, GM segmentations were downsampled to 2-mm resolution to match the resolution of the MD-MRI scans.

The Johns Hopkins University DTI-based WM atlas (https://identifiers.org/neurovault.collection:264) was used to generate segmentations of white matter tracts. First, by using the participant-to-cohort average image transforms computed earlier (see Section 2.6.1), we formed a cohort-averaged E[DΔ2] image (2-mm isotropic resolution). Then, we registered that image to the 2-mm resolution fractional anisotropy image of the WM atlas. Finally, we applied the inverses of the computed transforms (participant to cohort average and cohort average to WM atlas) to transform the ROIs of the WM atlas to the native space of each participant and repetition.

The segmented ROIs we focused on in our ROI analysis are shown in study-specific template space in Figure 1. They covered different regions of the brain and included subcortical GM regions (basal ganglia and thalamus), cortical GM regions (middle occipital gyrus and inferior frontal gyrus), commissural tracts (genu, body, and splenium of the corpus callosum (CC)), projection tracts (internal capsule and corona radiata), and association tracts (sagittal stratum and external capsule). For each ROI that had a separate ROI for the left and right hemisphere of the brain, we combined the left and the right side to create a single ROI. To ascertain that each image voxel only belonged to one ROI, we removed from the GM ROIs those voxels that belonged to any of the WM ROIs. Finally, we removed all voxels with fbin3>0.2 (estimated CSF signal fraction > 0.2) from the ROIs. Some of the ROIs, shown for the different repetitions and participants in native space, are shown in Supplementary Figure 3.

Fig. 1.

Regions of interest (ROIs) used in the ROI analysis as shown on the cohort-averaged T1-weighted template. Subcortical (basal ganglia and thalamus) and cortical (middle occipital gyrus and inferior frontal gyrus) gray matter ROIs were produced using the spatially localized atlas network tiles method (SLANT). Commissural (genu, body, and splenium of the corpus callosum (CC)), projection (internal capsule and corona radiata), and association (sagittal stratum and external capsule) white matter tracts were produced using the Johns Hopkins University diffusion tensor imaging-based white matter atlas. The slice coordinates shown are approximate coordinates in the Montreal Neurological Institute (MNI) space. Abbreviations: L-R, left-right; P-A, posterior-anterior; I-S, inferior-superior.

Fig. 1.

Regions of interest (ROIs) used in the ROI analysis as shown on the cohort-averaged T1-weighted template. Subcortical (basal ganglia and thalamus) and cortical (middle occipital gyrus and inferior frontal gyrus) gray matter ROIs were produced using the spatially localized atlas network tiles method (SLANT). Commissural (genu, body, and splenium of the corpus callosum (CC)), projection (internal capsule and corona radiata), and association (sagittal stratum and external capsule) white matter tracts were produced using the Johns Hopkins University diffusion tensor imaging-based white matter atlas. The slice coordinates shown are approximate coordinates in the Montreal Neurological Institute (MNI) space. Abbreviations: L-R, left-right; P-A, posterior-anterior; I-S, inferior-superior.

Close modal

2.7.2 Estimation of ICC and CVws

The estimation of ICC (reliability) and CVws (repeatability) values for each ROI and MD-MRI parameter map was performed in the native space of each participant and repetition to avoid unnecessary image interpolation. We employed an approach typical to ROI analysis, that is, we averaged the parameter value over the voxels within each ROI. Then, ICC and CVws were computed from those ROI-averaged parameter values. Thus, we get estimates for the reliability and repeatability of ROI-averaged MD-MRI parameter maps. Spider (radar) plots were used to visualize differences in MD-MRI parameter values, ICC, and CVws between different ROIs (modified MATLAB function spider_plot, https://github.com/NewGuy012/spider_plot, accessed Aug 28, 2024).

2.7.3 Bland–Altman plots

Bland–Altman plots (Bland & Altman, 1986) were used to assess the presence of a bias between the test and retest measurement for each MD-MRI parameter map. In a Bland–Altman plot, the difference between the test and retest measurement is plotted as a function of the mean of the test and the retest measurement for each data point (in our case, each voxel). To achieve this, the parameter maps and ROIs of each participant and repetition were transformed to the cohort template space. In the cohort template space, a template ROI was formed by the majority vote of the ROIs of the participants and repetitions (at least 50% agreement between the participants and repetitions was needed for a voxel to belong to a certain ROI). Each voxel in the cohort template space that belonged to any of the ROIs included in our ROI analysis was included as a data point in the Bland–Altman plot, for each participant, resulting in 152620 data points (15262 voxels, 10 participants). The MATLAB class “densityScatterChart” (https://github.com/MATLAB-Graphics-and-App-Building/densityScatterChart, accessed July 17, 2023) was used to plot the density scatter plots. The mean of the test–retest difference over all data points (voxels and participants), along with the 95% confidence interval of the observations, was also computed and plotted for each MD-MRI parameter. Finally, we reported the relative mean test–retest difference as a percentage of the grand mean (mean over repetitions, participants, and voxels).

3.1 Profiling of white matter, cortical, and subcortical regions in the human brain

Leveraging the rich information content of the voxel-wise D(ω)- R1- R2 distributions, we first provide robust microstructural profiling across various brain regions, based on the average data from the 20 scans in our study. The ROIs used in this study are shown in Figure 1, and projections from full D(ω)- R1- R2 distributions in six representative ROIs, comprising two cortical and two subcortical regions, along with two WM tracts are shown in Figure 2.

Fig. 2.

Projections of estimated D(ω)- R1- R2 distributions in different regions of interest (ROI). The displayed 2D projections include Diso- DΔ2 (top row), Diso- R1 (middle row), and Diso- R2 (bottom row) distributions. Different grayscale contour lines represent different centroid frequencies for the diffusion encoding gradients (ω/2π=6.621 Hz). Colored areas in the Diso- DΔ2 distribution plots represent bins that approximately correspond to white matter (red bin), gray matter (green bin), and cerebrospinal fluid (blue bin). Each distribution was first clustered across bootstraps, averaged over all voxels within each ROI, then averaged over the 10 participants and 2 repetitions. (A) Middle occipital gyrus (MOG), (B) inferior occipital gyrus (IFG), (C) basal ganglia, (D) thalamus, (E) splenium of the corpus callosum (CC), (F) sagittal stratum.

Fig. 2.

Projections of estimated D(ω)- R1- R2 distributions in different regions of interest (ROI). The displayed 2D projections include Diso- DΔ2 (top row), Diso- R1 (middle row), and Diso- R2 (bottom row) distributions. Different grayscale contour lines represent different centroid frequencies for the diffusion encoding gradients (ω/2π=6.621 Hz). Colored areas in the Diso- DΔ2 distribution plots represent bins that approximately correspond to white matter (red bin), gray matter (green bin), and cerebrospinal fluid (blue bin). Each distribution was first clustered across bootstraps, averaged over all voxels within each ROI, then averaged over the 10 participants and 2 repetitions. (A) Middle occipital gyrus (MOG), (B) inferior occipital gyrus (IFG), (C) basal ganglia, (D) thalamus, (E) splenium of the corpus callosum (CC), (F) sagittal stratum.

Close modal

For each ROI, the distributions are visualized as projections onto the 2D Diso- DΔ2, Diso- R1, and Diso- R2 planes for 5 frequencies ω/2π between 6.6 and 21 Hz. The well-known differences in tissue microstructure are most clearly manifested in the 2D Diso- DΔ2 projection (Topgaard, 2019), which is, therefore, used to define three subvoxel bins, as illustrated by the red, green, and blue overlays in Figure 2. The bin signal fractions, fbin1, fbin2, and fbin3, evaluated at the lowest frequency (ω/2π = 6.6 Hz), are obtained by integrating over these spectral regions. Additional 2D projections that display the 2D DΔ2- Diso, DΔ2- R1, and DΔ2- R2 (i.e., the anisotropy as a function of relaxation) are shown in Supplementary Figure 4 to allow for a more complete investigation of the specific microstructural features of different brain regions.

The distributions shown in Figure 2 reflect the microstructural differences of the respective regions and their subvoxel tissue components, revealing differences in composition between cortical and subcortical GM, and WM. Specifically for the 2D Diso- DΔ2 projections, cortical GM regions such as the inferior frontal gyrus are mainly composed of morphologically isotropic tissue elements with low diffusivity (Fig. 2A-B, green shading). Subcortical regions such as the thalamus present a wider distribution of tissue elements in terms of both their microscopic shape and size, with an additional highly anisotropic component (Fig. 2D, red shading); and WM tracts such as the splenium of the corpus callosum are mainly composed of morphologically anisotropic and low diffusivity tissue elements (Fig. 2E-F).

The R1 and R2 relaxation rates can complement the microstructural picture with compositional information. The 2D Diso- R1 (Fig. 2) and DΔ2- R1 (Supplementary Fig. 4) projections show a gradual increase in the R1 value of the main component from the cortical GM to the subcortical GM to the WM regions. A low R1 peak (~0.4 s-1) was observed in the cortical regions, visible on the 1D R1 projection. In the 2D Diso- R2 and DΔ2- R2 projections, distinct and more heterogeneous trends are observed compared with R1, with a single R2 population in the cortical regions, compared with two R2 populations in the subcortical and WM regions. From the Diso- R2 correlations (Fig. 2), both fast- and slow-relaxing R2 components have qualitatively similar Diso values. On the other hand, from the DΔ2- R2 correlations (Supplementary Fig. 4), regional differences were observed in which a high anisotropy fast-relaxing population was only observed in the CC splenium, while in the basal ganglia, the thalamus, and the sagittal stratum, only low anisotropy fast-relaxing and high anisotropy slow-relaxing components have been observed.

Diffusion frequency dependency is shown as grayscale contours in Figure 2, and can be more readily observed from the 1D projections in each dimension. Although some frequency dependency can be seen for Diso, it is more pronounced when examining DΔ2. In all regions, frequency dependence of the microscopic anisotropy behaves in a similar manner, in which DΔ2 shifts toward lower values with increased frequency.

The voxel-wise multidimensional distributions can be used in different ways to extract information of interest. The most straightforward approach to achieve dimensionality reduction involves computing the per-voxel statistics, such as means E[x], for the different dimensions, assessed at one of the probed diffusion frequencies. This way, the means E[Diso], E[R1], and E[R2] correspond to conventional mean diffusivity (Basser et al., 1994), quantitative R1 and R2 (Weiskopf et al., 2021), respectively. The E[DΔ2] map is comparable with metrics used to quantify microscopic diffusion anisotropy (Lasič et al., 2014; Lawrenz et al., 2010; Shemesh et al., 2016). Expressing the diffusion frequency dependence, rate of change with frequency (within the investigated range of 6.6–21 Hz) for the per-voxel means and variances can be used (Narvaez et al., 2024). Thus, positive values of the Δω/2π E[Diso] parameter indicate diffusion time dependency behavior suggestive of restriction (Aggarwal et al., 2012), while decreased anisotropy with higher frequency results in negative values in the Δω/2π E[DΔ2] map. These MD-MRI diffusion–relaxation parameters, averaged within the ROIs in Figure 2, are presented in Figure 3. The mean values agree with the expected microstructural differences between GM and WM regions, and also correspond with the observed trends of the 2D projections in Figure 2.

Fig. 3.

Mean (solid line) and standard deviation (shaded area covers the range of ±1 standard deviation) of region-of-interest-averaged multidimensional MRI parameters. Abbreviations: BG, basal ganglia; Thal, thalamus; MOG, middle occipital gyrus; IFG, inferior frontal gyrus; CC genu, genu of the corpus callosum; CC body, body of the corpus callosum; CC splenium, splenium of the corpus callosum; IC, internal capsule; CR, corona radiata; SS, sagittal stratum; EC, external capsule.

Fig. 3.

Mean (solid line) and standard deviation (shaded area covers the range of ±1 standard deviation) of region-of-interest-averaged multidimensional MRI parameters. Abbreviations: BG, basal ganglia; Thal, thalamus; MOG, middle occipital gyrus; IFG, inferior frontal gyrus; CC genu, genu of the corpus callosum; CC body, body of the corpus callosum; CC splenium, splenium of the corpus callosum; IC, internal capsule; CR, corona radiata; SS, sagittal stratum; EC, external capsule.

Close modal

An alternative way to extract information from the high dimensional voxel-wise data is to perform integration across a predefined parameter range, akin to a “spectral” ROI. The obtained integral value, ranging from 0 to 1, indicates the signal fraction within a particular multidimensional distribution. This voxel-wise computation enables the creation of an image reflecting the selected spectral component (Benjamini & Basser, 2020; Benjamini, Bouhrara, et al., 2021; de Almeida Martins et al., 2021; Martin et al., 2021). In this study, we used the previously suggested spectral ROIs shown in Figure 2 (de Almeida Martins et al., 2020; Topgaard, 2019), to summarize the WM, GM, and CSF contributions. Voxel-wise signal fractions, averaged within the ROIs in Figure 2, are also presented in Figure 3. We can further leverage the derived subvoxel information from the Diso- DΔ2 plane bins illustrated in Figure 2 (i.e., fbin1, fbin2, and fbin3) and use it to resolve relaxation parameters based on their diffusion characteristics (de Almeida Martins et al., 2020; Reymbaut et al., 2021). Summarized values expressing diffusion–relaxation correlation via the (diffusion-) bin-resolved means of R1 and R2 relaxation properties are captured with the bin-resolved mean relaxation metrics, for example, E[R1]bin1, E[R2]bin2, which can also be found in Figure 3.

To complement the data in Figure 3, numerical values for the means and standard deviations from all ROIs of each MD-MRI parameter, computed over all participants, repetitions, are given in Supplementary Figure 5.

3.2 Test–retest reliability and repeatability

We hypothesized that in vivo human brain D(ω)- R1- R2 MD-MRI parameter estimates acquired using a clinically achievable protocol are repeatable and reliable. Below, we evaluate this hypothesis for two prevalent image analysis approaches: voxel-wise and ROI-based. We employ the intraclass correlation coefficient (ICC) to assess the reliability, the within-subject coefficient of variation (CVws) to quantify the repeatability, and Bland–Altman plots to estimate the bias of the MD-MRI parameter maps.

3.2.1 Voxel-wise analysis

We transformed MD-MRI parameter maps to the cohort template space and estimated ICC (reliability) and CVws (repeatability) values for each brain voxel, after applying spatial filtering to the maps to reduce the effects of image misalignment. Parametric maps of mean values and signal fractions are shown in Figure 4. Figure 4A shows a representative image slice of the MD-MRI parameter maps from a single subject. The maps are shown after transformation to the cohort template space but without the application of spatial filtering to give a representative example of the image quality. Figure 4B shows the cohort averages (10 participants, 2 repetitions) for those MD-MRI parameter maps. A representative slice for the MD-MRI parameter maps from all the participants in the study is given in Supplementary Figure 6 to visualize the intersubject variability. For each MD-MRI parameter, voxel-wise ICC and CVws values are shown in Figure 4C and 4D, respectively. These were computed after applying a Gaussian spatial filter to the maps (FWHM = 2 pixels = 4 mm). Supplementary Figure 7 shows ICC and CVws values computed with no spatial filtering applied to the parameter maps. The subvoxel spectral information can be further utilized to resolve relaxation parameters based on their diffusion characteristics, which are depicted in Figure 5. Each map utilizes two orthogonal scales: the brightness intensity indicates the relative signal fraction, while the color scale denotes the specific relaxation property value. Voxel-wise ICC and CVws values are shown in Figure 5C and 5D, respectively.

Fig. 4.

Voxel-wise MD-MRI parameter means E[x] and signal fraction maps, and their reliability and repeatability measures. Means E[x] and signal fraction maps, shown for (A) a representative single participant, and for a (B) cohort average. Voxel-wise (C) intraclass correlation (ICC; reliability) and (D) within-subject coefficient of variation (CVws; repeatability), computed for each parameter map. Spatial smoothing was applied to each map before the ICC and CVws computations. The ICC and CVws values demonstrate good reliability and repeatability for the first order statistical parameters E[x] of Diso, DΔ2, R1, and R2, and low reliability and repeatability for the dependence on diffusion gradient frequency Δω/2π[x] of the diffusion parameters. The bin signal fractions fbin1 and fbin2 demonstrate high reliability and repeatability. The bin signal fraction fbin3 shows low repeatability in white matter but high reliability and repeatability elsewhere.

Fig. 4.

Voxel-wise MD-MRI parameter means E[x] and signal fraction maps, and their reliability and repeatability measures. Means E[x] and signal fraction maps, shown for (A) a representative single participant, and for a (B) cohort average. Voxel-wise (C) intraclass correlation (ICC; reliability) and (D) within-subject coefficient of variation (CVws; repeatability), computed for each parameter map. Spatial smoothing was applied to each map before the ICC and CVws computations. The ICC and CVws values demonstrate good reliability and repeatability for the first order statistical parameters E[x] of Diso, DΔ2, R1, and R2, and low reliability and repeatability for the dependence on diffusion gradient frequency Δω/2π[x] of the diffusion parameters. The bin signal fractions fbin1 and fbin2 demonstrate high reliability and repeatability. The bin signal fraction fbin3 shows low repeatability in white matter but high reliability and repeatability elsewhere.

Close modal
Fig. 5.

Voxel-wise MD-MRI parameter for bin-resolved relaxation values, and their reliability and repeatability measures. MD-MRI parameter R1 and R2 means, resolved according to their respective diffusion bin fractions, shown for (A) a representative single participant, and for (B) a cohort-average. Voxel-wise (C) intraclass correlation coefficient (ICC; reliability) and (D) within-subject coefficient of variation (CVws; repeatability), computed for each parameter map. Spatial smoothing was applied to each map before the ICC and CVws computations. The ICC and CVws values demonstrate good reliability and repeatability except for fbin3, which shows low repeatability in white matter but high reliability and repeatability elsewhere.

Fig. 5.

Voxel-wise MD-MRI parameter for bin-resolved relaxation values, and their reliability and repeatability measures. MD-MRI parameter R1 and R2 means, resolved according to their respective diffusion bin fractions, shown for (A) a representative single participant, and for (B) a cohort-average. Voxel-wise (C) intraclass correlation coefficient (ICC; reliability) and (D) within-subject coefficient of variation (CVws; repeatability), computed for each parameter map. Spatial smoothing was applied to each map before the ICC and CVws computations. The ICC and CVws values demonstrate good reliability and repeatability except for fbin3, which shows low repeatability in white matter but high reliability and repeatability elsewhere.

Close modal

Figure 6 shows the distributions of voxel-wise ICC and CVws values for each MD-MRI parameter map. With most MD-MRI parameter maps, high ICC values (high reliability) corresponded to low CVws values (high repeatability), as could be expected. E[Diso], E[DΔ2], and E[R2] exhibited both high ICC values and low CVws values. E[R1] showed one of the lowest CVws values but also low ICC values, suggesting low relative measurement error but also small biological variability between the participants. The estimated signal fraction parameter fbin3 had relatively high ICC but high CVws values as well. The high CVws values of fbin3 were concentrated in WM regions. In those regions, fbin3 values are close to 0, and so even minute differences in fbin3 between repetitions lead to high CVws values. We also estimated the distributions of the ICC and CVws values over all brain voxels. Supplementary Figure 8 shows the distributions of ICC and CVws values in the brain with different levels of spatial filtering applied to the parameter maps.

Fig. 6.

Summary of region of interest (ROI)-based and voxel-wise (A) reliability (intraclass correlation coefficient, ICC) and (B) repeatability (within-subject coefficient of variation, CVws) results. For ROI-based results (blue dots), ICC and CVws were computed for each ROI-averaged multidimensional (MD)-MRI parameter map, and a boxplot was plotted for the collection of ROI-wise results. Voxel-wise ICC and CVws values were plotted as violin plots (orange). Crosses (x) indicate that the parameter in question contained values outside of the shown axis limits (CVws>0.20).

Fig. 6.

Summary of region of interest (ROI)-based and voxel-wise (A) reliability (intraclass correlation coefficient, ICC) and (B) repeatability (within-subject coefficient of variation, CVws) results. For ROI-based results (blue dots), ICC and CVws were computed for each ROI-averaged multidimensional (MD)-MRI parameter map, and a boxplot was plotted for the collection of ROI-wise results. Voxel-wise ICC and CVws values were plotted as violin plots (orange). Crosses (x) indicate that the parameter in question contained values outside of the shown axis limits (CVws>0.20).

Close modal

3.2.2 ROI analysis

The ROIs used in the analysis are shown in Figure 1. Within each ROI and for each participant and repetition, the mean MD-MRI parameter map was computed, from which ICC (reliability) and CVws (repeatability) were estimated. As the parameters are averaged over the ROIs, the number of voxels in each ROI affects the quality of the metrics. All ROIs in this study had comparable volume (~1500 voxels), with the exception of the corona radiata (~3500 voxels) and the sagittal stratum (~500 voxels). The number of voxels for each participant and repetition in each ROI is shown in Supplementary Figure 9.

Reliability and repeatability values for different MD-MRI parameter maps and ROIs are plotted in Figure 7, and their numerical values are given in Supplementary Figure 10. Figure 7 allows to visually assess uniformity of the reliability and repeatability across the examined ROIs for individual MD-MRI parameters. Boxplots over all ROIs for the ROI-based ICC and CVws values, shown in Figure 6, indicate, as expected, that ROI-averaged MD-MRI parameters are more reliable and repeatable than their voxel-wise counterparts. The relationship between high reliability (high ICC) and high repeatability (low CVws) in the ROI-based analysis was similar to the results of the voxel-wise analysis, except in the case of fbin3. The reliability (ICC) of fbin3 benefited more from the ROI averaging than the other MD-MRI parameters, achieving the highest median ROI-based reliability (ICC = 0.94) out of all the MD-MRI parameters even though its median repeatability was the 4th worst (CVws = 0.106). As in the voxel-wise analysis, the ROI-averaged values of E[Diso], E[DΔ2], and E[R2] displayed both high reliability and high repeatability.

Fig. 7.

Reliability (intraclass correlation coefficient, ICC; blue line and axis ticks) and repeatability (within-subject coefficient of variation, CVws; red line and axis ticks) of region-of-interest-averaged MD-MRI parameter maps. CVws>0.20values are marked as crosses on the CVws=0.20 line. Generally, ROI-based ICC values were higher and CVws values lower than those of individual voxels (Fig. 6). That was expected due to the increase in signal-to-noise ratio achieved through averaging. Abbreviations: BG, basal ganglia; Thal, thalamus; MOG, middle occipital gyrus; IFG, inferior frontal gyrus; CC genu, genu of the corpus callosum; CC body, body of the corpus callosum; CC splenium, splenium of the corpus callosum; IC, internal capsule; CR, corona radiata; SS, sagittal stratum; EC, external capsule.

Fig. 7.

Reliability (intraclass correlation coefficient, ICC; blue line and axis ticks) and repeatability (within-subject coefficient of variation, CVws; red line and axis ticks) of region-of-interest-averaged MD-MRI parameter maps. CVws>0.20values are marked as crosses on the CVws=0.20 line. Generally, ROI-based ICC values were higher and CVws values lower than those of individual voxels (Fig. 6). That was expected due to the increase in signal-to-noise ratio achieved through averaging. Abbreviations: BG, basal ganglia; Thal, thalamus; MOG, middle occipital gyrus; IFG, inferior frontal gyrus; CC genu, genu of the corpus callosum; CC body, body of the corpus callosum; CC splenium, splenium of the corpus callosum; IC, internal capsule; CR, corona radiata; SS, sagittal stratum; EC, external capsule.

Close modal

We assessed the presence of bias between the test and the retest measurement using Bland–Altman plots comprising the different participants and the voxels included in the analyzed ROIs (Fig. 8). For each of the MD-MRI parameters, the 95% confidence interval for the observations included the zero-bias level. The relative mean differences between the test and the retest measurement for the different parameters were E[Diso]: 0.55%, E[DΔ2]: -0.031%, E[R1]: -0.66%, E[R2]:-0.65%, Δω/2πE[Diso]: 11%, Δω/2πE[DΔ2]: 3.4%, fbin1: -0.60%, fbin2: 1.1%, fbin3: -1.7%, E[R1]bin1: -0.84%, E[R1]bin2: -0.51%, E[R1]bin3: 2.0%, E[R2]bin1: -0.96%, E[R2]bin2: 0.51%, E[R2]bin3: 6.0%. Those biases were, on average, about one order of magnitude smaller (range: 4 to 140 times smaller) than the median voxel-wise within-subject coefficients of variation CVws (Supplementary Fig. 5).

Fig. 8.

Bland–Altman density plots for the different MD-MRI parameter maps. In each plot, all the voxels that were included in the region of interest analysis are included as data points, for all participants. The mean (solid red line) and the 95% confidence interval of the observations (1.96 ×standard deviation; dashed red line) are shown for each plot.

Fig. 8.

Bland–Altman density plots for the different MD-MRI parameter maps. In each plot, all the voxels that were included in the region of interest analysis are included as data points, for all participants. The mean (solid red line) and the 95% confidence interval of the observations (1.96 ×standard deviation; dashed red line) are shown for each plot.

Close modal

Diffusion–relaxation correlation MD-MRI is a method that allows to acquire rich information about tissue microstructure by exploring the correlations between different MRI contrasts. This framework replaces voxel-averaged quantities with multidimensional distributions of those quantities that allows to selectively extract and map specific diffusion–relaxation spectral ranges. As a result, it has the potential of achieving high sensitivity and specificity toward detecting subtle changes that would have been otherwise averaged out. The purpose of this work was to methodologically characterize different brain regions in terms of their multicomponent diffusion–relaxation properties, interpret differences and similarities between them, and to assess the framework’s overall repeatability. We compared WM tracts, cortical GM, and subcortical GM regions, and found pronounced differences in their diffusion–relaxation profile, reflecting distinct diversities in their microscopic morphological content. We then applied a test–retest paradigm to estimate the repeatability and reliability of parameters derived from the 40-min whole brain in vivo MD-MRI framework, and showed them to be comparable with known DTI reproducibility.

Results from comprehensive profiling of representative WM, cortical GM, and subcortical GM regions are shown in Figure 2, Supplementary Figure 4, and Figure 3. Given the extensive history of DTI studies, it is not surprising that the MD-MRI data allowed to easily differentiate between WM and GM regions based on their different diffusion anisotropy distributions. However, we demonstrated here that meaningful differences can be gleaned even between WM tracts: compared with the sagittal stratum, the CC splenium showed a reduced content of low anisotropy – slow diffusivity components (i.e., fbin2, see Figure 3 and Supplementary Fig. 5), and had, in addition, a higher content of a fast-relaxing R2 component. Further, this fast-relaxing R2 population displayed high microscopic anisotropy. The observed differences in the R2 components and in E[R2]bin1 can be attributed to the higher myelin content in the CC splenium compared with the sagittal stratum (Dvorak et al., 2021), which would appear as a fast-relaxing microscopically anisotropic peak. The differences in the diffusive properties may be driven by the layered structure of the sagittal stratum (Di Carlo et al., 2019) that makes it more microstructurally heterogeneous compared with the CC.

Our findings in GM revealed distinct microstructural attributes that differentiate between subcortical and cortical regions. Compared with the cortical regions, the thalamus and basal ganglia contained a higher content of microscopically anisotropic components (i.e., fbin1, see Fig. 3 and Supplementary Fig. 5), and a wider distribution of diffusivities, skewed toward larger values. These differences in microscopic anisotropy and in diffusivity are likely to be driven, respectively, by the presence of fiber tracts and the relatively large subcortical nuclei (Morel, 2007). In addition to the diffusion properties, cortical and subcortical GM regions showed differences in their R2 distributions (see Fig. 2), in which the more microstructurally heterogeneous subcortical regions displayed a fast-relaxing R2 peak that may be modulated by the presence of macromolecules and myelin lipids.

While we were able to apply the MD-MRI framework to perform a microstructural study of the human brain thus demonstrating its potential, there are two factors that may introduce unwanted variability to its outcomes: the acquisition parameter space is sparsely sampled to keep the scan time feasible, and the parameter estimation is an ill-posed inverse problem. Until now, there has been no evaluation of MD-MRI parameter reproducibility in the living human brain. Therefore, we investigated the reliability (ICC) and repeatability (CVws) of estimated MD-MRI parameters using a test–retest paradigm. The CVws quantifies the relative error between repetitions of the same measurement (Luque Laguna et al., 2020). Thus, high repeatability means that the difference between the repetitions is small compared with the measured value. The ICC quantifies the fraction of variance that is not explained by errors, that is, due to differences between the participants (Matheson, 2019). High reliability means that biological differences between participants can be detected reliably. Like CVws, ICC depends on errors but, unlike CVws, also on differences between participants. Thus, our results for reliability are representative of our sample of participants (10 healthy participants with a mean age of 48 years, standard deviation 14.4 years; 4 women). The results could be different for a different age or sex distribution or in the presence of disease or pathology. It should also be noted that the reproducibility results are affected by the method used to estimate the measures from MD-MRI data (i.e., both the signal processing pipeline and the template-based ROI segmentation).

We investigated the reliability and repeatability of both voxel-based and ROI-averaged MD-MRI parameter estimates. The MD-MRI distribution mean parameters all had median voxel-wise CVws values smaller than 5%, exhibiting good repeatability. Voxel-wise E[Diso] and E[DΔ2], demonstrated good reliability (ICC = 0.75) and E[R2] moderate reliability (ICC = 0.72). However, E[R1] had poor reliability (ICC<0.50), which may be due to the relatively narrow range of R1 encoding sensitivity in the acquisition protocol, leading to limited sensitivity. The ROI-averaged distribution means E[Diso], E[DΔ2], and E[R2] had high repeatability (CVws < 2%) and reliability (ICC > 0.85).

The diffusion gradient frequency dependency of mean isotropic diffusivity, Δω/2π E[Diso], exhibited poor voxel-wise and ROI-averaged reliability and repeatability. Δω/2π E[DΔ2] had poor reliability but reasonably good repeatability (voxel-wise CVws = 18%, ROI CVws = 7.5%). We believe that these results have two main reasons. First, diffusion gradient hardware limitations constrain the diffusion frequency range to centroid frequencies of 6.6–21 Hz, with the highest frequency achieved only at a low b-value of 0.5 ms/μm2, and thus restrict the ability to reliably decompose and estimate the expected time scales in the system. This frequency content range, assuming spherical compartment, zero permeability, and intracellular diffusivity 3.5 μm2/ms, is expected to be particularly sensitive to the length scale of approximately 13 µm (Stepišnik, 1993). While this length scale may be characteristic in some brain regions, typical length scales in most of the human brain are significantly smaller. This mismatch between the sensitivity range of the MD-MRI acquisition and the underlying microstructure is a potential driver of the poor reproducibility we found. We believe that the second contributor to these findings is the way in which the frequency-dependence parameters are defined, that is, as a subtraction of two parameters. Thus, they suffer from error propagation as inaccuracies are summated during subtraction. Furthermore, the absolute differences between the highest and lowest frequencies are very small due to the limited frequency range, leading to small numbers around 0, with reduced sensitivity. Based on these results, finding an alternative way to quantify frequency dependence from the D(ω)- R1- R2 distributions is encouraged.

One of the unique features of the MD-MRI framework is the ability to access and visualize intravoxel information within the human brain by selectively integrating regions of the voxel-wise D(ω)- R1- R2 distributions. Consistent with previous studies, we delineated three spectral ROIs in the 2D Diso- DΔ2 distribution space, roughly corresponding to WM, GM, and CSF, denoted as bin1, bin2, and bin3, respectively, as illustrated in Figure 2. Voxel-wise application of partial integration based on these bins yields corresponding signal fraction maps, labeled as fbin1, fbin2, and fbin3. Voxel-wise reliability and repeatability analyses resulted in whole brain median ICC values of 0.67, 0.68, and 0.59, and CVws values of 6%, 15%, and 90%, for fbin1, fbin2, and fbin3, respectively. The ROI-based performance was better, with median ICC values ranging from 0.80 to 0.94, and CVws values in the range of 2–11%. We hypothesize that the relatively poor reproducibility of fbin3 was due to the small amounts of free water in non-CSF regions in the brain, leading to high relative errors in the estimated signal fractions. This effect is particularly exacerbated when reporting median values across the whole brain, with only a small portion of the voxels residing in CSF-containing regions, as can be seen in the almost binary-looking fbin3 map in Figure 4. The overall good reproducibility of the intravoxel fraction maps is encouraging, especially because such subvoxel component maps have demonstrated specificity toward pathology and have been used to characterize subtle processes noninvasively (Benjamini, Iacono, et al., 2021; Benjamini et al., 2022; Slator et al., 2019).

The subvoxel spectral information can be further utilized to resolve relaxation parameters based on their diffusion characteristics, such as the bin-resolved mean relaxation rates E[R1]bin1, E[R1]bin2, E[R1]bin3, E[R2]bin1, E[R2]bin2, and E[R2]bin3, thus directly visualizing diffusion–relaxation correlation. Of those bin-resolved parameters, E[R1]bin2, E[R1]bin3, E[R2]bin1, and E[R2]bin2 had moderate median voxel-wise reliabilities (ICC > 0.5) and good ROI-averaged reliabilities (ICC > 0.75). Their repeatability values were relatively good (voxel-wise CVws 4-12%, ROI CVws 2-5%), except those of E[R1]bin3, which were poorer (voxel-wise CVws = 100%, ROI CVws = 16%). E[R1]bin1 had low voxel-wise reliability (ICC = 0.39), moderate ROI reliability (ICC = 0.54), and good repeatability (voxel-wise CVws = 5.5%, ROI CVws = 3.1%). Interestingly, E[R2]bin3 had poor voxel-wise reliability (ICC = 0.29) and repeatability (CVws = 130%) but good ROI-averaged reliability (ICC = 0.81). Taken together, these results suggest that MD-MRI scalar parameters that combine information from different contrasts (diffusion and relaxation) can be reliable and repeatable. However, parameters related to bin3 (CSF) provide a special case to consider. Such parameters had considerably lower, even very poor, median voxel-wise reliability and repeatability values. Averaging those parameters over an ROI greatly increased their reliability and repeatability values. The voxel-wise estimates could have been greatly influenced by the low amount of CSF in most tissues and by inaccuracies in image registration.

Comparing our reliability and repeatability results with previous studies would provide important context. While the current study is the first to report on reproducibility of MD-MRI, we can draw parallels with previous diffusion MRI variability investigations. Specifically, MD and FA in our study align with E[Diso] and E[DΔ2], respectively. Most of the recent DTI studies investigated ICC and CVws in WM regions. Repeatability and reliability of DTI metrics in the CC reported, for MD and FA, respectively, CVws ranges of 1–2.5% and 2.5–3%, and ICC ranges of 0.54–0.73 and 0.76–0.98 (Coelho et al., 2022; Fan et al., 2021; Grech-Sollars et al., 2015; Luque Laguna et al., 2020). Our investigation in the CC revealed, for E[Diso] and E[DΔ2], both an average CVws of 1%, and average ICC of 0.88 and 0.91, respectively. In addition to E[Diso] and E[DΔ2], the reproducibility of E[R1] and E[R2] can also be compared with relaxometry studies in the CC, in which the CVws was established at under 1% for both T1 and T2 (Hagiwara et al., 2019), comparable with our reported values around 2%. Another rough comparison could be drawn between the free water fraction from the so-called Standard Model (Novikov et al., 2019) and MD-MRI derived fbin3. A recent reproducibility study found CVws values of 9% and 15% in the CC and internal capsule, respectively (Coelho et al., 2022), comparable with CVws average values of 10% and 16% we reported here for fbin3. Notably, the reproducibility of MD-MRI diffusion measures in the subcortical regions we investigated surpassed those of reported DTI parameters. For instance, reliability for MD and FA in the basal ganglia (CVws of 4% and 10%) and in the thalamus (CVws of 10% and 8%) (Grech-Sollars et al., 2015) should be compared with the corresponding values of 2% for E[Diso] and 4% for E[DΔ2] reported in our study.

Our study has several limitations and assumptions. The MD-MRI signal representation used provides nonparametric distributions of diffusion and relaxation components and is based on the Gaussian phase distribution (GPD) approximation (Neuman, 1974; Stepišnik, 1981). Non-Gaussian diffusion effects that violate the GPD assumption (Jespersen et al., 2019) and induce microscopic kurtosis are possible (Novello et al., 2022). In this case, the microscopic kurtosis is not accounted for in our model, and its effect would split between variances of Diso and DΔ2. Considering relaxation, quantitative comparisons of R1 and R2 across studies are hindered by fundamental considerations in pulse sequence and acquisition design. The recorded signal results from intricate interactions among partial excitation, relaxation processes, and proton pool exchange, each possessing distinct MR properties (Manning et al., 2021). These factors make R1 and, to a lesser extent, R2, dependent on specific pulse sequence parameters, slice thickness, and radiofrequency bandwidth. Furthermore, hardware limitations restricted the minimal echo time to 40 ms in this study, which may fully attenuate myelin water (R2~ 100 s-1) (Manning et al., 2021). Future improvements in diffusion gradient frequency range, b-values, and echo time minimization could be achieved by utilizing next-generation diffusion gradients (Dai et al., 2023; Huang et al., 2021; Michael et al., 2022). When assessing the reproducibility results, we should keep in mind that the study included 10 participants, and the results related to reliability are representative of that sample. Although consistent with the median sample size of six reported in technical MRI studies (Hanspach et al., 2021), a larger sample would likely yield more robust MD-MRI reproducibility values. Additionally, our focus on intrascanner reproducibility highlights the need for future investigations into interscanner assessment.

While informative, voxel-wise reproducibility analysis requires image registration to a common space, introducing inaccuracies due to biological variation and image interpolation. Although we applied Gaussian spatial filtering to mitigate these effects, complete elimination cannot be achieved. The level of smoothing we applied for our voxel-wise results is consistent with levels of smoothing typically applied (FWHM: 2–4 pixels; 4–10 mm), indicating that our results are representative of this field of research (Caceres et al., 2009; Shou et al., 2013; Somandepalli et al., 2015). However, smoothing is not always applied in voxel-wise analyses (Coelho et al., 2022; Luque Laguna et al., 2020; Veraart et al., 2021). Similar to the voxel-wise analysis, the ROI-based analysis is also affected by the accuracy of the segmentation results (Supplementary Fig. 3). For ROI-based analysis, we first estimated voxel-wise MD-MRI parameters and then averaged them over the ROI, as this approach reduces tissue heterogeneity in the estimation step. Conversely, averaging signals before parameter estimation would increase the expected number of distinct components, potentially affecting accuracy, and was, therefore, avoided.

We established reliable quantification of relaxation and diffusion properties, including intravoxel information, using a clinically feasible, 40-min in vivo diffusion–relaxation correlation (D(ω)- R1- R2) multidimensional MRI scan in the human brain. This enabled us to methodologically explore brain regions in terms of their diffusion–relaxation multicomponent properties and interpret differences and similarities between them. We were able to observe subtle variations between subcortical and cortical regions, and diversity of WM tracts, demonstrating the microstructural sensitivity of the MD-MRI framework. We further established repeatability and reliability of the pipeline, and showed it to be comparable with known diffusion MRI reproducibility. The possibilities of using the rich spectral information in the MD-MRI data to create biomarkers that show high sensitivity to specific cell types or pathologies should be further investigated. To ensure the reproducibility and adoption of such biomarkers, future studies should also investigate the interscanner reproducibility and the harmonization of MD-MRI datasets.

The data that support the findings of this study are available on request from the corresponding author. The data are not publicly available due to privacy or ethical restrictions. Code to process the MD-MRI data is freely available as implemented in the multidimensional diffusion MRI toolbox (https://github.com/markus-nilsson/md-dmri).

E.M.: Conceptualization, Methodology, Software, Formal Analysis, Writing – Original Draft, Writing – Review & Editing, Visualization; S.B.: Software, Resources, Writing – Review & Editing; B.A.L.: Software, Resources, Writing – Review & Editing; Y.Y.: Resources, Funding Acquisition, Writing – Review & Editing; D.T.: Conceptualization, Methodology, Software, Writing – Review & Editing; D.B.: Conceptualization, Methodology, Software, Investigation, Resources, Writing – Original Draft, Writing – Review & Editing, Visualization, Supervision, Project Administration, Funding Acquisition.

This work was funded by the Intramural Research Programs of the National Institute on Aging (NIA), the National Institute on Drug Abuse (NIDA) of the National Institutes of Health (NIH). D.T. was supported by the Swedish Foundation for Strategic Research (Grant No. ITM17-0267), and the Swedish Research Council (Grant No. 2022-04422_VR).

The authors declare no conflict of interest.

The authors would like to thank Mr. Phil Cholak for facilitating the MRI scans.

Supplementary material for this article is available with the online version here: https://doi.org/10.1162/imag_a_00387.

Aggarwal
,
M.
,
Jones
,
M. V.
,
Calabresi
,
P. A.
,
Mori
,
S.
, &
Zhang
,
J.
(
2012
).
Probing mouse brain microstructure using oscillating gradient diffusion MRI
.
Magnetic Resonance in Medicine
,
67
(
1
),
98
109
. https://doi.org/10.1002/mrm.22981
Arbabi
,
A.
,
Kai
,
J.
,
Khan
,
A. R.
, &
Baron
,
C. A.
(
2020
).
Diffusion dispersion imaging: Mapping oscillating gradient spin-echo frequency dependence in the human brain
.
Magnetic Resonance in Medicine
,
83
(
6
),
2197
2208
. https://doi.org/10.1002/mrm.28083
Assaf
,
Y.
, &
Basser
,
P. J.
(
2005
).
Composite hindered and restricted model of diffusion (CHARMED) MR imaging of the human brain
.
NeuroImage
,
27
(
1
),
48
58
. https://doi.org/10.1016/j.neuroimage.2005.03.042
Avants
,
B. B.
,
Tustison
,
N. J.
,
Song
,
G.
,
Cook
,
P. A.
,
Klein
,
A.
, &
Gee
,
J. C.
(
2011
).
A reproducible evaluation of ANTs similarity metric performance in brain image registration
.
NeuroImage
,
54
(
3
),
2033
2044
. https://doi.org/10.1016/j.neuroimage.2010.09.025
Basser
,
P. J.
,
Mattiello
,
J.
, &
LeBihan
,
D.
(
1994
).
MR diffusion tensor spectroscopy and imaging
.
Biophysical Journal
,
66
(
1
),
259
267
. https://doi.org/10.1016/S0006-3495(94)80775-1
Basser
,
P. J.
, &
Pajevic
,
S.
(
2003
).
A normal distribution for tensor-valued random variables: Applications to diffusion tensor MRI
.
IEEE Transactions on Medical Imaging
,
22
(
7
),
785
794
. https://doi.org/10.1109/TMI.2003.815059
Benjamini
,
D.
(
2020
).
Nonparametric inversion of relaxation and diffusion correlation data
. In
D.
Topgaard
(Ed.),
Advanced diffusion encoding methods in MRI
(pp.
278
316
).
The Royal Society of Chemistry
. https://doi.org/10.1039/9781788019910-00278
Benjamini
,
D.
, &
Basser
,
P. J.
(
2016
).
Use of marginal distributions constrained optimization (MADCO) for accelerated 2D MRI relaxometry and diffusometry
.
Journal of Magnetic Resonance
,
271
,
40
45
. https://doi.org/10.1016/j.jmr.2016.08.004
Benjamini
,
D.
, &
Basser
,
P. J.
(
2017
).
Magnetic resonance microdynamic imaging reveals distinct tissue microenvironments
.
NeuroImage
,
163
,
183
196
. https://doi.org/10.1016/j.neuroimage.2017.09.033
Benjamini
,
D.
, &
Basser
,
P. J.
(
2020
).
Multidimensional correlation MRI
.
NMR in Biomedicine
,
33
(
12
),
e4226
. https://doi.org/10.1002/nbm.4226
Benjamini
,
D.
,
Bouhrara
,
M.
,
Komlosh
,
M. E.
,
Iacono
,
D.
,
Perl
,
D. P.
,
Brody
,
D. L.
, &
Basser
,
P. J.
(
2021
).
Multidimensional MRI for characterization of subtle axonal injury accelerated using an adaptive nonlocal multispectral filter
.
Frontiers in Physics
,
9
,
737374
. https://doi.org/10.3389/fphy.2021.737374
Benjamini
,
D.
,
Hutchinson
,
E. B.
,
Komlosh
,
M. E.
,
Comrie
,
C. J.
,
Schwerin
,
S. C.
,
Zhang
,
G.
,
Pierpaoli
,
C.
, &
Basser
,
P. J.
(
2020
).
Direct and specific assessment of axonal injury and spinal cord microenvironments using diffusion correlation imaging
.
NeuroImage
,
221
,
117195
. https://doi.org/10.1016/j.neuroimage.2020.117195
Benjamini
,
D.
,
Iacono
,
D.
,
Komlosh
,
M. E.
,
Perl
,
D. P.
,
Brody
,
D. L.
, &
Basser
,
P. J.
(
2021
).
Diffuse axonal injury has a characteristic multidimensional MRI signature in the human brain
.
Brain
,
144
(
3
),
800
816
. https://doi.org/10.1093/brain/awaa447
Benjamini
,
D.
,
Komlosh
,
M. E.
,
Holtzclaw
,
L. A.
,
Nevo
,
U.
, &
Basser
,
P. J.
(
2016
).
White matter microstructure from nonparametric axon diameter distribution mapping
.
NeuroImage
,
135
,
333
344
. https://doi.org/10.1016/j.neuroimage.2016.04.052
Benjamini
,
D.
,
Priemer
,
D. S.
,
Perl
,
D. P.
,
Brody
,
D. L.
, &
Basser
,
P. J.
(
2022
).
Mapping astrogliosis in the individual human brain using multidimensional MRI
.
Brain
,
146
(
3
),
1212
1226
. https://doi.org/10.1093/brain/awac298
Bland
,
J. M.
, &
Altman
,
D. G.
(
1986
).
Statistical methods for assessing agreement between two methods of clinical measurement
.
The Lancet
,
327
(
8476
),
307
310
. https://doi.org/10.1016/S0140-6736(86)90837-8
Caceres
,
A.
,
Hall
,
D. L.
,
Zelaya
,
F. O.
,
Williams
,
S. C. R.
, &
Mehta
,
M. A.
(
2009
).
Measuring fMRI reliability with the intra-class correlation coefficient
.
NeuroImage
,
45
(
3
),
758
768
. https://doi.org/10.1016/j.neuroimage.2008.12.035
Callaghan
,
P. T.
(
2011
).
Translational dynamics and magnetic resonance: Principles of pulsed gradient spin echo NMR
.
Oxford University Press
. https://doi.org/10.1093/acprof:oso/9780199556984.001.0001
Callaghan
,
P. T.
, &
Stepišnik
,
J.
(
1995
).
Frequency-domain analysis of spin motion using modulated-gradient NMR
.
Journal of Magnetic Resonance, Series A
,
117
(
1
),
118
122
. https://doi.org/10.1006/jmra.1995.9959
Coelho
,
S.
,
Baete
,
S. H.
,
Lemberskiy
,
G.
,
Ades-Aron
,
B.
,
Barrol
,
G.
,
Veraart
,
J.
,
Novikov
,
D. S.
, &
Fieremans
,
E.
(
2022
).
Reproducibility of the Standard Model of diffusion in white matter on clinical MRI systems
.
NeuroImage
,
257
,
119290
. https://doi.org/10.1016/j.neuroimage.2022.119290
Conturo
,
T. E.
,
McKinstry
,
R. C.
,
Akbudak
,
E.
, &
Robinson
,
B. H.
(
1996
).
Encoding of anisotropic diffusion with tetrahedral gradients: A general mathematical diffusion formalism and experimental results
.
Magnetic Resonance in Medicine
,
35
(
3
),
399
412
. https://doi.org/10.1002/mrm.1910350319
Dai
,
E.
,
Zhu
,
A.
,
Yang
,
G. K.
,
Quah
,
K.
,
Tan
,
E. T.
,
Fiveland
,
E.
,
Foo
,
T. K. F.
, &
McNab
,
J. A.
(
2023
).
Frequency-dependent diffusion kurtosis imaging in the human brain using an oscillating gradient spin echo sequence and a high-performance head-only gradient
.
NeuroImage
,
279
,
120328
. https://doi.org/10.1016/j.neuroimage.2023.120328
de Almeida Martins
,
J. P.
,
Tax
,
C. M. W.
,
Reymbaut
,
A.
,
Szczepankiewicz
,
F.
,
Chamberland
,
M.
,
Jones
,
D. K.
, &
Topgaard
,
D.
(
2021
).
Computing and visualising intra-voxel orientation-specific relaxation–diffusion features in the human brain
.
Human Brain Mapping
,
42
(
2
),
310
328
. https://doi.org/10.1002/hbm.25224
de Almeida Martins
,
J. P.
,
Tax
,
C. M. W.
,
Szczepankiewicz
,
F.
,
Jones
,
D. K.
,
Westin
,
C. F.
, &
Topgaard
,
D.
(
2020
).
Transferring principles of solid-state and Laplace NMR to the field of in vivo brain MRI
.
Magnetic Resonance
,
1
(
1
),
27
43
. https://doi.org/10.5194/mr-1-27-2020
de Almeida Martins
,
J. P.
, &
Topgaard
,
D.
(
2016
).
Two-dimensional correlation of isotropic and directional diffusion using NMR
.
Physical Review Letters
,
116
(
8
),
087601
. https://doi.org/10.1103/PhysRevLett.116.087601
de Almeida Martins
,
J. P.
, &
Topgaard
,
D.
(
2018
).
Multidimensional correlation of nuclear relaxation rates and diffusion tensors for model-free investigations of heterogeneous anisotropic porous materials
.
Scientific Reports
,
8
(
1
),
2488
. https://doi.org/10.1038/s41598-018-19826-9
Di Carlo
,
D. T.
,
Benedetto
,
N.
,
Duffau
,
H.
,
Cagnazzo
,
F.
,
Weiss
,
A.
,
Castagna
,
M.
,
Cosottini
,
M.
, &
Perrini
,
P.
(
2019
).
Microsurgical anatomy of the sagittal stratum
.
Acta Neurochirurgica
,
161
(
11
),
2319
2327
. https://doi.org/10.1007/s00701-019-04019-8
Dvorak
,
A. V.
,
Swift-LaPointe
,
T.
,
Vavasour
,
I. M.
,
Lee
,
L. E.
,
Abel
,
S.
,
Russell-Schulz
,
B.
,
Graf
,
C.
,
Wurl
,
A.
,
Liu
,
H.
,
Laule
,
C.
,
Li
,
D. K. B.
,
Traboulsee
,
A.
,
Tam
,
R.
,
Boyd
,
L. A.
,
MacKay
,
A. L.
, &
Kolind
,
S. H.
(
2021
).
An atlas for human brain myelin content throughout the adult life span
.
Scientific Reports
,
11
(
1
),
269
. https://doi.org/10.1038/s41598-020-79540-3
Einstein
,
A.
(
1905
).
Über die von der molekularkinetischen Theorie der Wärme geforderte Bewegung von in ruhenden Flüssigkeiten suspendierten Teilchen
.
Annalen der Physik
,
322
(
8
),
549
560
. https://doi.org/10.1002/andp.19053220806
Fan
,
Q.
,
Polackal
,
M. N.
,
Tian
,
Q.
,
Ngamsombat
,
C.
,
Nummenmaa
,
A.
,
Witzel
,
T.
,
Klawiter
,
E. C.
, &
Huang
,
S. Y.
(
2021
).
Scan-rescan repeatability of axonal imaging metrics using high-gradient diffusion MRI and statistical implications for study design
.
NeuroImage
,
240
,
118323
. https://doi.org/10.1016/j.neuroimage.2021.118323
Fieremans
,
E.
,
Burcaw
,
L. M.
,
Lee
,
H.-H.
,
Lemberskiy
,
G.
,
Veraart
,
J.
, &
Novikov
,
D. S.
(
2016
).
In vivo observation and biophysical interpretation of time-dependent diffusion in human white matter
.
NeuroImage
,
129
,
414
427
. https://doi.org/10.1016/j.neuroimage.2016.01.018
Grech-Sollars
,
M.
,
Hales
,
P. W.
,
Miyazaki
,
K.
,
Raschke
,
F.
,
Rodriguez
,
D.
,
Wilson
,
M.
,
Gill
,
S. K.
,
Banks
,
T.
,
Saunders
,
D. E.
,
Clayden
,
J. D.
,
Gwilliam
,
M. N.
,
Barrick
,
T. R.
,
Morgan
,
P. S.
,
Davies
,
N. P.
,
Rossiter
,
J.
,
Auer
,
D. P.
,
Grundy
,
R.
,
Leach
,
M. O.
,
Howe
,
F. A.
,
Clark
,
C.
A
. (
2015
).
Multi-centre reproducibility of diffusion MRI parameters for clinical sequences in the brain
.
NMR in Biomedicine
,
28
(
4
),
468
485
. https://doi.org/10.1002/nbm.3269
Hagiwara
,
A.
,
Hori
,
M.
,
Cohen-Adad
,
J.
,
Nakazawa
,
M.
,
Suzuki
,
Y.
,
Kasahara
,
A.
,
Horita
,
M.
,
Haruyama
,
T.
,
Andica
,
C.
,
Maekawa
,
T.
,
Kamagata
,
K.
,
Kumamaru
,
K. K.
,
Abe
,
O.
, &
Aoki
,
S.
(
2019
).
Linearity, bias, intrascanner repeatability, and interscanner reproducibility of quantitative multidynamic multiecho sequence for rapid simultaneous relaxometry at 3 T: A validation study with a standardized phantom and healthy controls
.
Investigative Radiology
,
54
(
1
),
39
47
. https://doi.org/10.1097/rli.0000000000000510
Hamilton
,
J.
,
Xu
,
K.
,
Geremia
,
N.
,
Prado
,
V. F.
,
Prado
,
M. A. M.
,
Brown
,
A.
, &
Baron
,
C. A.
(
2024
).
Robust frequency-dependent diffusional kurtosis computation using an efficient direction scheme, axisymmetric modelling, and spatial regularization
.
Imaging Neuroscience
,
2
,
1
22
. https://doi.org/10.1162/imag_a_00055
Hanspach
,
J.
,
Nagel
,
A. M.
,
Hensel
,
B.
,
Uder
,
M.
,
Koros
,
L.
, &
Laun
,
F. B.
(
2021
).
Sample size estimation: Current practice and considerations for original investigations in MRI technical development studies
.
Magnetic Resonance in Medicine
,
85
(
4
),
2109
2116
. https://doi.org/10.1002/mrm.28550
Hennel
,
F.
,
Michael
,
E. S.
, &
Pruessmann
,
K. P.
(
2021
).
Improved gradient waveforms for oscillating gradient spin-echo (OGSE) diffusion tensor imaging
.
NMR in Biomedicine
,
34
(
2
),
e4434
. https://doi.org/10.1002/nbm.4434
Huang
,
S. Y.
,
Witzel
,
T.
,
Keil
,
B.
,
Scholz
,
A.
,
Davids
,
M.
,
Dietz
,
P.
,
Rummert
,
E.
,
Ramb
,
R.
,
Kirsch
,
J. E.
,
Yendiki
,
A.
,
Fan
,
Q.
,
Tian
,
Q.
,
Ramos-Llordén
,
G.
,
Lee
,
H.-H.
,
Nummenmaa
,
A.
,
Bilgic
,
B.
,
Setsompop
,
K.
,
Wang
,
F.
,
Avram
,
A. V.
,
Rosen
,
B.
R
. (
2021
).
Connectome 2.0: Developing the next-generation ultra-high gradient strength human MRI scanner for bridging studies of the micro-, meso- and macro-connectome
.
NeuroImage
,
243
,
118530
. https://doi.org/10.1016/j.neuroimage.2021.118530
Huo
,
Y.
,
Xu
,
Z.
,
Xiong
,
Y.
,
Aboud
,
K.
,
Parvathaneni
,
P.
,
Bao
,
S.
,
Bermudez
,
C.
,
Resnick
,
S. M.
,
Cutting
,
L. E.
, &
Landman
,
B. A.
(
2019
).
3D whole brain segmentation using spatially localized atlas network tiles
.
NeuroImage
,
194
,
105
119
. https://doi.org/10.1016/j.neuroimage.2019.03.041
Hürlimann
,
M. D.
, &
Venkataramanan
,
L.
(
2002
).
Quantitative measurement of two-dimensional distribution functions of diffusion and relaxation in grossly inhomogeneous fields
.
Journal of Magnetic Resonance
,
157
(
1
),
31
42
. https://doi.org/10.1006/jmre.2002.2567
Irfanoglu
,
M. O.
,
Modi
,
P.
,
Nayak
,
A.
,
Hutchinson
,
E. B.
,
Sarlls
,
J.
, &
Pierpaoli
,
C.
(
2015
).
DR-BUDDI (diffeomorphic registration for blip-up blip-down diffusion imaging) method for correcting echo planar imaging distortions
.
NeuroImage
,
106
,
284
299
. https://doi.org/10.1016/j.neuroimage.2014.11.042
Jespersen
,
S. N.
,
Olesen
,
J. L.
,
Ianuş
,
A.
, &
Shemesh
,
N.
(
2019
).
Effects of nongaussian diffusion on “isotropic diffusion” measurements: An ex-vivo microimaging and simulation study
.
Journal of Magnetic Resonance
,
300
,
84
94
. https://doi.org/10.1016/j.jmr.2019.01.007
Jian
,
B.
,
Vemuri
,
B. C.
,
Özarslan
,
E.
,
Carney
,
P. R.
, &
Mareci
,
T. H.
(
2007
).
A novel tensor distribution model for the diffusion-weighted MR signal
.
NeuroImage
,
37
(
1
),
164
176
. https://doi.org/10.1016/j.neuroimage.2007.03.074
Johnson
,
J. T. E.
,
Irfanoglu
,
M. O.
,
Manninen
,
E.
,
Ross
,
T. J.
,
Yang
,
Y.
,
Laun
,
F. B.
,
Martin
,
J.
,
Topgaard
,
D.
, &
Benjamini
,
D.
(
2024
).
In vivo disentanglement of diffusion frequency-dependence, tensor shape, and relaxation using multidimensional MRI
.
Human Brain Mapping
,
45
(
7
),
e26697
. https://doi.org/10.1002/hbm.26697
Kellner
,
E.
,
Dhital
,
B.
,
Kiselev
,
V. G.
, &
Reisert
,
M.
(
2016
).
Gibbs-ringing artifact removal based on local subvoxel-shifts
.
Magnetic Resonance in Medicine
,
76
(
5
),
1574
1581
. https://doi.org/10.1002/mrm.26054
Kim
,
D.
,
Doyle
,
E. K.
,
Wisnowski
,
J. L.
,
Kim
,
J. H.
, &
Haldar
,
J. P.
(
2017
).
Diffusion-relaxation correlation spectroscopic imaging: A multidimensional approach for probing microstructure
.
Magnetic Resonance in Medicine
,
78
(
6
),
2236
2249
. https://doi.org/10.1002/mrm.26629
Koo
,
T. K.
, &
Li
,
M. Y.
(
2016
).
A guideline of selecting and reporting intraclass correlation coefficients for reliability research
.
Journal of Chiropractic Medicine
,
15
(
2
),
155
163
. https://doi.org/10.1016/j.jcm.2016.02.012
Kroeker
,
R. M.
, &
Mark
Henkelman
, R
. (
1986
).
Analysis of biological NMR relaxation data with continuous distributions of relaxation times
.
Journal of Magnetic Resonance (1969)
,
69
(
2
),
218
235
. https://doi.org/10.1016/0022-2364(86)90074-0
Kundu
,
S.
,
Barsoum
,
S.
,
Ariza
,
J.
,
Nolan
,
A. L.
,
Latimer
,
C. S.
,
Keene
,
C. D.
,
Basser
,
P. J.
, &
Benjamini
,
D.
(
2023
).
Mapping the individual human cortex using multidimensional MRI and unsupervised learning
.
Brain Communications
,
5
(
6
),
fcad258
. https://doi.org/10.1093/braincomms/fcad258
Lampinen
,
B.
,
Szczepankiewicz
,
F.
,
Novén
,
M.
,
van Westen
,
D.
,
Hansson
,
O.
,
Englund
,
E.
,
Mårtensson
,
J.
,
Westin
,
C.-F.
, &
Nilsson
,
M.
(
2019
).
Searching for the neurite density with diffusion MRI: Challenges for biophysical modeling
.
Human Brain Mapping
,
40
(
8
),
2529
2545
. https://doi.org/10.1002/hbm.24542
Lasič
,
S.
,
Szczepankiewicz
,
F.
,
Eriksson
,
S.
,
Nilsson
,
M.
, &
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
.
Frontiers in Physics
,
2
,
11
. https://doi.org/10.3389/fphy.2014.00011
Lasič
,
S.
,
Yuldasheva
,
N.
,
Szczepankiewicz
,
F.
,
Nilsson
,
M.
,
Budde
,
M.
,
Dall’Armellina
,
E.
,
Schneider
,
J. E.
,
Teh
,
I.
, &
Lundell
,
H.
(
2022
).
Stay on the beat with tensor-valued encoding: Time-dependent diffusion and cell size estimation in ex vivo heart
.
Frontiers in Physics
,
10
. https://doi.org/10.3389/fphy.2022.812115
Lawrenz
,
M.
,
Koch
,
M. A.
, &
Finsterbusch
,
J.
(
2010
).
A tensor model and measures of microscopic anisotropy for double-wave-vector diffusion-weighting experiments with long mixing times
.
Journal of Magnetic Resonance
,
202
(
1
),
43
56
. https://doi.org/10.1016/j.jmr.2009.09.015
Lee
,
H.-H.
,
Novikov
,
D. S.
, &
Fieremans
,
E.
(
2021
).
Removal of partial Fourier-induced Gibbs (RPG) ringing artifacts in MRI
.
Magnetic Resonance in Medicine
,
86
(
5
),
2733
2750
. https://doi.org/10.1002/mrm.28830
Lee
,
H.-H.
,
Papaioannou
,
A.
,
Novikov
,
D. S.
, &
Fieremans
,
E.
(
2020
).
In vivo observation and biophysical interpretation of time-dependent diffusion in human cortical gray matter
.
NeuroImage
,
222
,
117054
. https://doi.org/10.1016/j.neuroimage.2020.117054
Livet
,
J.
,
Weissman
,
T. A.
,
Kang
,
H.
,
Draft
,
R. W.
,
Lu
,
J.
,
Bennis
,
R. A.
,
Sanes
,
J. R.
, &
Lichtman
,
J. W.
(
2007
).
Transgenic strategies for combinatorial expression of fluorescent proteins in the nervous system
.
Nature
,
450
(
7166
),
56
62
. https://doi.org/10.1038/nature06293
Lundell
,
H.
, &
Lasič
,
S.
(
2020
).
Diffusion encoding with general gradient waveforms
. In
D.
Topgaard
(Ed.),
Advanced diffusion encoding methods in MRI
(pp.
12
16
).
The Royal Society of Chemistry
. https://doi.org/10.1039/9781788019910-00012
Luque Laguna
,
P. A.
,
Combes
,
A. J. E.
,
Streffer
,
J.
,
Einstein
,
S.
,
Timmers
,
M.
,
Williams
,
S. C. R.
, &
Dell’Acqua
,
F.
(
2020
).
Reproducibility, reliability and variability of FA and MD in the older healthy population: A test-retest multiparametric analysis
.
NeuroImage: Clinical
,
26
,
102168
. https://doi.org/10.1016/j.nicl.2020.102168
Mackay
,
A.
,
Whittall
,
K.
,
Adler
,
J.
,
Li
,
D.
,
Paty
,
D.
, &
Graeb
,
D.
(
1994
).
In vivo visualization of myelin water in brain by magnetic resonance
.
Magnetic Resonance in Medicine
,
31
(
6
),
673
677
. https://doi.org/10.1002/mrm.1910310614
Maekawa
,
T.
,
Hori
,
M.
,
Murata
,
K.
,
Feiweier
,
T.
,
Kamiya
,
K.
,
Andica
,
C.
,
Hagiwara
,
A.
,
Fujita
,
S.
,
Koshino
,
S.
,
Akashi
,
T.
,
Kamagata
,
K.
,
Wada
,
A.
,
Abe
,
O.
, &
Aoki
,
S.
(
2020
).
Differentiation of high-grade and low-grade intra-axial brain tumors by time-dependent diffusion MRI
.
Magnetic Resonance Imaging
,
72
,
34
41
. https://doi.org/10.1016/j.mri.2020.06.018
Magdoom
,
K. N.
,
Pajevic
,
S.
,
Dario
,
G.
, &
Basser
,
P. J.
(
2021
).
A new framework for MR diffusion tensor distribution
.
Scientific Reports
,
11
(
1
),
2766
. https://doi.org/10.1038/s41598-021-81264-x
Manning
,
A. P.
,
MacKay
,
A. L.
, &
Michal
,
C. A.
(
2021
).
Understanding aqueous and non-aqueous proton T1 relaxation in brain
.
Journal of Magnetic Resonance
,
323
,
106909
. https://doi.org/10.1016/j.jmr.2020.106909
Martin
,
J.
,
Reymbaut
,
A.
,
Schmidt
,
M.
,
Doerfler
,
A.
,
Uder
,
M.
,
Laun
,
F. B.
, &
Topgaard
,
D.
(
2021
).
Nonparametric D-R1-R2 distribution MRI of the living human brain
.
NeuroImage
,
245
,
118753
. https://doi.org/10.1016/j.neuroimage.2021.118753
Matheson
,
G. J.
(
2019
).
We need to talk about reliability: Making better use of test-retest studies for study design and interpretation
.
PeerJ
,
7
,
e6918
. https://doi.org/10.7717/peerj.6918
McGraw
,
K. O.
, &
Wong
,
S. P.
(
1996
).
Forming inferences about some intraclass correlation coefficients
.
Psychological Methods
,
1
(
1
),
30
46
. https://doi.org/10.1037/1082-989X.1.1.30
Michael
,
E. S.
,
Hennel
,
F.
, &
Pruessmann
,
K. P.
(
2022
).
Evaluating diffusion dispersion across an extended range of b-values and frequencies: Exploiting gap-filled OGSE shapes, strong gradients, and spiral readouts
.
Magnetic Resonance in Medicine
,
87
(
6
),
2710
2723
. https://doi.org/10.1002/mrm.29161
Morel
,
A.
(
2007
).
Stereotactic atlas of the human thalamus and basal ganglia
(1 ed.).
CRC Press
. https://doi.org/10.3109/9781420016796
Naranjo
,
I. D.
,
Reymbaut
,
A.
,
Brynolfsson
,
P.
,
Lo
Gullo
, R.,
Bryskhe
,
K.
,
Topgaard
,
D.
,
Giri
,
D. D.
,
Reiner
,
J. S.
,
Thakur
,
S. B.
, &
Pinker-Domenig
,
K.
(
2021
).
Multidimensional diffusion magnetic resonance imaging for characterization of tissue microstructure in breast cancer patients: A prospective pilot study
.
Cancers
,
13
(
7
),
1606
. https://doi.org/10.3390/cancers13071606
Narvaez
,
O.
,
Svenningsson
,
L.
,
Yon
,
M.
,
Sierra
,
A.
, &
Topgaard
,
D.
(
2022
).
Massively multidimensional diffusion-relaxation correlation MRI
.
Frontiers in Physics
,
9
,
793966
. https://doi.org/10.3389/fphy.2021.793966
Narvaez
,
O.
,
Yon
,
M.
,
Jiang
,
H.
,
Bernin
,
D.
,
Forssell-Aronsson
,
E.
,
Sierra
,
A.
, &
Topgaard
,
D.
(
2024
).
Nonparametric distributions of tensor-valued Lorentzian diffusion spectra for model-free data inversion in multidimensional diffusion MRI
.
The Journal of Chemical Physics
,
161
(
8
),
084201
. https://doi.org/10.1063/5.0213252
Neuman
,
C. H.
(
1974
).
Spin echo of spins diffusing in a bounded medium
.
The Journal of Chemical Physics
,
60
(
11
),
4508
4511
. https://doi.org/10.1063/1.1680931
Nilsson
,
M.
,
Szczepankiewicz
,
F.
,
Lampinen
,
B.
,
Ahlgren
,
A.
,
de Almeida Martins
,
J. P.
,
Lasic
,
S.
,
Westin
,
C.-F.
, &
Topgaard
,
D.
(
2018
).
An open-source framework for analysis of multidimensional diffusion MRI data implemented in MATLAB
.
Proceedings of the International Society for Magnetic Resonance in Medicine
.
Paris, France
. http://archive.ismrm.org/2018/5355.html
Novello
,
L.
,
Henriques
,
R. N.
,
Ianuş
,
A.
,
Feiweier
,
T.
,
Shemesh
,
N.
, &
Jovicich
,
J.
(
2022
).
In vivo correlation tensor MRI reveals microscopic kurtosis in the human brain on a clinical 3T scanner
.
NeuroImage
,
254
,
119137
. https://doi.org/10.1016/j.neuroimage.2022.119137
Novikov
,
D. S.
,
Fieremans
,
E.
,
Jespersen
,
S. N.
, &
Kiselev
,
V. G.
(
2019
).
Quantifying brain microstructure with diffusion MRI: Theory and parameter estimation
.
NMR in Biomedicine
,
32
(
4
),
e3998
. https://doi.org/10.1002/nbm.3998
Novikov
,
D. S.
,
Jensen
,
J. H.
,
Helpern
,
J. A.
, &
Fieremans
,
E.
(
2014
).
Revealing mesoscopic structural universality with diffusion
.
Proceedings of the National Academy of Sciences
,
111
(
14
),
5088
5093
. https://doi.org/10.1073/pnas.1316944111
Novikov
,
D. S.
,
Kiselev
,
V. G.
, &
Jespersen
,
S. N.
(
2018
).
On modeling
.
Magnetic Resonance in Medicine
,
79
(
6
),
3172
3193
. https://doi.org/10.1002/mrm.27101
Peled
,
S.
,
Cory
,
D. G.
,
Raymond
,
S. A.
,
Kirschner
,
D. A.
, &
Jolesz
,
F. A.
(
1999
).
Water diffusion, T2, and compartmentation in frog sciatic nerve
.
Magnetic Resonance in Medicine
,
42
(
5
),
911
918
. https://doi.org/10.1002/(SICI)1522-2594(199911)42:5<911::AID-MRM11>3.0.CO;2-J
Pierpaoli
,
C.
,
Jezzard
,
P.
,
Basser
,
P. J.
,
Barnett
,
A.
, &
Chiro
,
G. D.
(
1996
).
Diffusion tensor MR imaging of the human brain
.
Radiology
,
201
(
3
),
637
648
. https://doi.org/10.1148/radiology.201.3.8939209
Prange
,
M.
, &
Song
,
Y.-Q.
(
2009
).
Quantifying uncertainty in NMR T2 spectra using Monte Carlo inversion
.
Journal of Magnetic Resonance
,
196
(
1
),
54
60
. https://doi.org/10.1016/j.jmr.2008.10.008
Price
,
W. S.
(
2009
).
NMR studies of translational motion: Principles and applications
.
Cambridge University Press
. https://doi.org/10.1017/CBO9780511770487
Reuter
,
M.
,
Schmansky
,
N. J.
,
Rosas
,
H. D.
, &
Fischl
,
B.
(
2012
).
Within-subject template estimation for unbiased longitudinal image analysis
.
NeuroImage
,
61
(
4
),
1402
1418
. https://doi.org/10.1016/j.neuroimage.2012.02.084
Reymbaut
,
A.
,
Critchley
,
J.
,
Durighel
,
G.
,
Sprenger
,
T.
,
Sughrue
,
M.
,
Bryskhe
,
K.
, &
Topgaard
,
D.
(
2021
).
Toward nonparametric diffusion – Characterization of crossing fibers in the human brain
.
Magnetic Resonance in Medicine
,
85
(
5
),
2815
2827
. https://doi.org/10.1002/mrm.28604
Reymbaut
,
A.
,
Mezzani
,
P.
,
de Almeida Martins
,
J. P.
, &
Topgaard
,
D.
(
2020
).
Accuracy and precision of statistical descriptors obtained from multidimensional diffusion signal inversion algorithms
.
NMR in Biomedicine
,
33
(
12
),
e4267
. https://doi.org/10.1002/nbm.4267
Rohde
,
G. K.
,
Barnett
,
A. S.
,
Basser
,
P. J.
,
Marenco
,
S.
, &
Pierpaoli
,
C.
(
2004
).
Comprehensive approach for correction of motion and distortion in diffusion-weighted MRI
.
Magnetic Resonance in Medicine
,
51
(
1
),
103
114
. https://doi.org/10.1002/mrm.10677
Schilling
,
K. G.
,
Blaber
,
J.
,
Huo
,
Y.
,
Newton
,
A.
,
Hansen
,
C.
,
Nath
,
V.
,
Shafer
,
A. T.
,
Williams
,
O.
,
Resnick
,
S. M.
,
Rogers
,
B.
,
Anderson
,
A. W.
, &
Landman
,
B. A.
(
2019
).
Synthesized b0 for diffusion distortion correction (Synb0-DisCo)
.
Magnetic Resonance Imaging
,
64
,
62
70
. https://doi.org/10.1016/j.mri.2019.05.008
Shemesh
,
N.
,
Jespersen
,
S. N.
,
Alexander
,
D. C.
,
Cohen
,
Y.
,
Drobnjak
,
I.
,
Dyrby
,
T. B.
,
Finsterbusch
,
J.
,
Koch
,
M. A.
,
Kuder
,
T.
,
Laun
,
F.
,
Lawrenz
,
M.
,
Lundell
,
H.
,
Mitra
,
P. P.
,
Nilsson
,
M.
,
Özarslan
,
E.
,
Topgaard
,
D.
, &
Westin
,
C.-F.
(
2016
).
Conventions and nomenclature for double diffusion encoding NMR and MRI
.
Magnetic Resonance in Medicine
,
75
(
1
),
82
87
. https://doi.org/10.1002/mrm.25901
Shou
,
H.
,
Eloyan
,
A.
,
Lee
,
S.
,
Zipunnikov
,
V.
,
Crainiceanu
,
A. N.
,
Nebel
,
M. B.
,
Caffo
,
B.
,
Lindquist
,
M. A.
, &
Crainiceanu
,
C. M.
(
2013
).
Quantifying the reliability of image replication studies: The image intraclass correlation coefficient (I2C2)
.
Cognitive, Affective, & Behavioral Neuroscience
,
13
(
4
),
714
724
. https://doi.org/10.3758/s13415-013-0196-0
Silva
,
M. D.
,
Helmer
,
K. G.
,
Lee
,
J.-H.
,
Han
,
S. S.
,
Springer
,
C. S.
, &
Sotak
,
C. H.
(
2002
).
Deconvolution of compartmental water diffusion coefficients in yeast-cell suspensions using combined T1 and diffusion measurements
.
Journal of Magnetic Resonance
,
156
(
1
),
52
63
. https://doi.org/10.1006/jmre.2002.2527
Sjölund
,
J.
,
Szczepankiewicz
,
F.
,
Nilsson
,
M.
,
Topgaard
,
D.
,
Westin
,
C.-F.
, &
Knutsson
,
H.
(
2015
).
Constrained optimization of gradient waveforms for generalized diffusion encoding
.
Journal of Magnetic Resonance
,
261
,
157
168
. https://doi.org/10.1016/j.jmr.2015.10.012
Slator
,
P. J.
,
Hutter
,
J.
,
Palombo
,
M.
,
Jackson
,
L. H.
,
Ho
,
A.
,
Panagiotaki
,
E.
,
Chappell
,
L. C.
,
Rutherford
,
M. A.
,
Hajnal
,
J. V.
, &
Alexander
,
D. C.
(
2019
).
Combined diffusion-relaxometry MRI to identify dysfunction in the human placenta
.
Magnetic Resonance in Medicine
,
82
(
1
),
95
106
. https://doi.org/10.1002/mrm.27733
Slator
,
P. J.
,
Palombo
,
M.
,
Miller
,
K. L.
,
Westin
,
C.-F.
,
Laun
,
F.
,
Kim
,
D.
,
Haldar
,
J. P.
,
Benjamini
,
D.
,
Lemberskiy
,
G.
,
de Almeida Martins
,
J. P.
, &
Hutter
,
J.
(
2021
).
Combined diffusion-relaxometry microstructure imaging: Current status and future prospects
.
Magnetic Resonance in Medicine
,
86
(
6
),
2987
3011
. https://doi.org/10.1002/mrm.28963
Soares
,
J.
,
Marques
,
P.
,
Alves
,
V.
, &
Sousa
,
N.
(
2013
).
A hitchhiker’s guide to diffusion tensor imaging
.
Frontiers in Neuroscience
,
7
,
31
. https://doi.org/10.3389/fnins.2013.00031
Somandepalli
,
K.
,
Kelly
,
C.
,
Reiss
,
P. T.
,
Zuo
,
X.-N.
,
Craddock
,
R. C.
,
Yan
,
C.-G.
,
Petkova
,
E.
,
Castellanos
,
F. X.
,
Milham
,
M. P.
, &
Di
Martino
, A
. (
2015
).
Short-term test–retest reliability of resting state fMRI metrics in children with and without attention-deficit/hyperactivity disorder
.
Developmental Cognitive Neuroscience
,
15
,
83
93
. https://doi.org/10.1016/j.dcn.2015.08.003
Stanisz
,
G. J.
, &
Henkelman
,
R. M.
(
1998
).
Diffusional anisotropy of T2 components in bovine optic nerve
.
Magnetic Resonance in Medicine
,
40
(
3
),
405
410
. https://doi.org/10.1002/mrm.1910400310
Stanisz
,
G. J.
,
Wright
,
G. A.
,
Henkelman
,
R. M.
, &
Szafer
,
A.
(
1997
).
An analytical model of restricted diffusion in bovine optic nerve
.
Magnetic Resonance in Medicine
,
37
(
1
),
103
111
. https://doi.org/10.1002/mrm.1910370115
Stepišnik
,
J.
(
1981
).
Analysis of NMR self-diffusion measurements by a density matrix calculation
.
Physica B+C
,
104
(
3
),
350
364
. https://doi.org/10.1016/0378-4363(81)90182-0
Stepišnik
,
J.
(
1993
).
Time-dependent self-diffusion by NMR spin-echo
.
Physica B: Condensed Matter
,
183
(
4
),
343
350
. https://doi.org/10.1016/0921-4526(93)90124-O
Stepišnik
,
J.
(
1998
).
Spin echo attenuation of restricted diffusion as a discord of spin phase structure
.
Journal of Magnetic Resonance
,
131
(
2
),
339
346
. https://doi.org/10.1006/jmre.1997.1232
Tan
,
E. T.
,
Shih
,
R. Y.
,
Mitra
,
J.
,
Sprenger
,
T.
,
Hua
,
Y.
,
Bhushan
,
C.
,
Bernstein
,
M. A.
,
McNab
,
J. A.
,
DeMarco
,
J. K.
,
Ho
,
V. B.
, &
Foo
,
T. K. F.
(
2020
).
Oscillating diffusion-encoding with a high gradient-amplitude and high slew-rate head-only gradient for human brain imaging
.
Magnetic Resonance in Medicine
,
84
(
2
),
950
965
. https://doi.org/10.1002/mrm.28180
Tétreault
,
P.
,
Harkins
,
K. D.
,
Baron
,
C. A.
,
Stobbe
,
R.
,
Does
,
M. D.
, &
Beaulieu
,
C.
(
2020
).
Diffusion time dependency along the human corpus callosum and exploration of age and sex differences as assessed by oscillating gradient spin-echo diffusion tensor imaging
.
NeuroImage
,
210
,
116533
. https://doi.org/10.1016/j.neuroimage.2020.116533
Tofts
,
P.
(
2003
).
Quantitative MRI of the brain: Measuring changes caused by disease
.
John Wiley & Sons, Ltd
. https://doi.org/10.1002/0470869526
Topgaard
,
D.
(
2019
).
Diffusion tensor distribution imaging
.
NMR in Biomedicine
,
32
(
5
),
e4066
. https://doi.org/10.1002/nbm.4066
Uhlenbeck
,
G. E.
, &
Ornstein
,
L. S.
(
1930
).
On the theory of the brownian motion
.
Physical Review
,
36
(
5
),
823
841
. https://doi.org/10.1103/PhysRev.36.823
Vasilescu
,
V.
,
Katona
,
E.
,
Simplâceanu
,
V.
, &
Demco
,
D.
(
1978
).
Water compartments in the myelinated nerve. III. Pulsed NMR result
.
Experientia
,
34
(
11
),
1443
1444
. https://doi.org/10.1007/BF01932339
Veraart
,
J.
,
Fieremans
,
E.
, &
Novikov
,
D. S.
(
2019
).
On the scaling behavior of water diffusion in human brain white matter
.
NeuroImage
,
185
,
379
387
. https://doi.org/10.1016/j.neuroimage.2018.09.075
Veraart
,
J.
,
Novikov
,
D. S.
,
Christiaens
,
D.
,
Ades-aron
,
B.
,
Sijbers
,
J.
, &
Fieremans
,
E.
(
2016
).
Denoising of diffusion MRI using random matrix theory
.
NeuroImage
,
142
,
394
406
. https://doi.org/10.1016/j.neuroimage.2016.08.016
Veraart
,
J.
,
Nunes
,
D.
,
Rudrapatna
,
U.
,
Fieremans
,
E.
,
Jones
,
D. K.
,
Novikov
,
D. S.
, &
Shemesh
,
N.
(
2020
).
Noninvasive quantification of axon radii using diffusion MRI
.
eLife
,
9
,
e49855
. https://doi.org/10.7554/eLife.49855
Veraart
,
J.
,
Raven
,
E. P.
,
Edwards
,
L. J.
,
Weiskopf
,
N.
, &
Jones
,
D. K.
(
2021
).
The variability of MR axon radii estimates in the human white matter
.
Human Brain Mapping
,
42
(
7
),
2201
2213
. https://doi.org/10.1002/hbm.25359
Wei
,
X.
,
Zhu
,
L.
,
Zeng
,
Y.
,
Xue
,
K.
,
Dai
,
Y.
,
Xu
,
J.
,
Liu
,
G.
,
Liu
,
F.
,
Xue
,
W.
,
Wu
,
D.
, &
Wu
,
G.
(
2022
).
Detection of prostate cancer using diffusion-relaxation correlation spectrum imaging with support vector machine model – A feasibility study
.
Cancer Imaging
,
22
(
1
),
77
. https://doi.org/10.1186/s40644-022-00516-9
Weiskopf
,
N.
,
Edwards
,
L. J.
,
Helms
,
G.
,
Mohammadi
,
S.
, &
Kirilina
,
E.
(
2021
).
Quantitative magnetic resonance imaging of brain anatomy and in vivo histology
.
Nature Reviews Physics
,
3
(
8
),
570
588
. https://doi.org/10.1038/s42254-021-00326-1
Westin
,
C.-F.
,
Szczepankiewicz
,
F.
,
Pasternak
,
O.
,
Özarslan
,
E.
,
Topgaard
,
D.
,
Knutsson
,
H.
, &
Nilsson
,
M.
(
2014
).
Measurement tensors in diffusion MRI: Generalizing the concept of diffusion encoding
.
Medical Image Computing and Computer-Assisted Intervention – MICCAI
2014
, Cham. https://doi.org/10.1007/978-3-319-10443-0_27
Wetscherek
,
A.
,
Stieltjes
,
B.
, &
Laun
,
F. B.
(
2015
).
Flow-compensated intravoxel incoherent motion diffusion imaging
.
Magnetic Resonance in Medicine
,
74
(
2
),
410
419
. https://doi.org/10.1002/mrm.25410
Woessner
,
D. E.
(
1963
).
N.M.R. spin-echo self-diffusion measurements on fluids undergoing restricted diffusion
.
The Journal of Physical Chemistry
,
67
(
6
),
1365
1367
. https://doi.org/10.1021/j100800a509
Xiong
,
Y.
,
Huo
,
Y.
,
Wang
,
J.
,
Davis
,
L. T.
,
McHugo
,
M.
, &
Landman
,
B.
(
2019
).
Reproducibility evaluation of SLANT whole brain segmentation across clinical magnetic resonance imaging protocols
(Vol. 10949).
SPIE
. https://doi.org/10.1117/12.2512561
Xu
,
J.
,
Jiang
,
X.
,
Devan
,
S. P.
,
Arlinghaus
,
L. R.
,
McKinley
,
E. T.
,
Xie
,
J.
,
Zu
,
Z.
,
Wang
,
Q.
,
Chakravarthy
,
A. B.
,
Wang
,
Y.
, &
Gore
,
J. C.
(
2021
).
MRI-cytometry: Mapping nonparametric cell size distributions using diffusion MRI
.
Magnetic Resonance in Medicine
,
85
(
2
),
748
761
. https://doi.org/10.1002/mrm.28454
Xu
,
J.
,
Jiang
,
X.
,
Li
,
H.
,
Arlinghaus
,
L. R.
,
McKinley
,
E. T.
,
Devan
,
S. P.
,
Hardy
,
B. M.
,
Xie
,
J.
,
Kang
,
H.
,
Chakravarthy
,
A. B.
, &
Gore
,
J. C.
(
2020
).
Magnetic resonance imaging of mean cell size in human breast tumors
.
Magnetic Resonance in Medicine
,
83
(
6
),
2002
2014
. https://doi.org/10.1002/mrm.28056
Xu
,
J.
,
Li
,
H.
,
Li
,
K.
,
Harkins
,
K. D.
,
Jiang
,
X.
,
Xie
,
J.
,
Kang
,
H.
,
Dortch
,
R. D.
,
Anderson
,
A. W.
,
Does
,
M. D.
, &
Gore
,
J. C.
(
2016
).
Fast and simplified mapping of mean axon diameter using temporal diffusion spectroscopy
.
NMR in Biomedicine
,
29
(
4
),
400
410
. https://doi.org/10.1002/nbm.3484
Yon
,
M.
,
de Almeida Martins
,
J. P.
,
Bao
,
Q.
,
Budde
,
M. D.
,
Frydman
,
L.
, &
Topgaard
,
D.
(
2020
).
Diffusion tensor distribution imaging of an in vivo mouse brain at ultrahigh magnetic field by spatiotemporal encoding
.
NMR in Biomedicine
,
33
(
11
),
e4355
. https://doi.org/10.1002/nbm.4355
Zhang
,
H.
,
Schneider
,
T.
,
Wheeler-Kingshott
,
C. A.
, &
Alexander
,
D. C.
(
2012
).
NODDI: Practical in vivo neurite orientation dispersion and density imaging of the human brain
.
NeuroImage
,
61
(
4
),
1000
1016
. https://doi.org/10.1016/j.neuroimage.2012.03.072
This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International (CC BY 4.0) license, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. For a full description of the license, please visit https://creativecommons.org/licenses/by/4.0/legalcode.

Supplementary data