Robust, fast and accurate mapping of diffusional mean kurtosis

  1. Megan E Farquhar
  2. Qianqian Yang  Is a corresponding author
  3. Viktor Vegh
  1. School of Mathematical Sciences, Faculty of Science, Queensland University of Technology, Australia
  2. Centre for Data Science, Queensland University of Technology, Australia
  3. Centre for Biomedical Technologies, Queensland University of Technology, Australia
  4. Centre for Advanced Imaging, The University of Queensland, Australia
  5. ARC Training Centre for Innovation in Biomedical Imaging Technology, Australia
10 figures, 5 tables and 1 additional file

Figures

Simulation results on the effect of the number of diffusion times involved in the fitting of the sub-diffusion model (Equation 3) parameters (Dβ,β) and computing K following (Equation 9) at various SNR levels.

The true values for (Dβ, β, K) are set to (3 × 10–4, 0.75, 0.8125) to represent white matter (blue) and (5 × 10–4, 0.85, 0.4733) to represent grey matter (orange). Rows 1–3: Distributions of fitted parameter values using different number of diffusion times. 2′ represents an additional simulation using two diffusion times but set to be the same, so it has the same number of data points in the fitting as for using two different diffusion times. Row 4: Coefficient of variation (CV) of the parameter values fitted using different number of diffusion times.

Surface plots of R2 values achieved with fitting simulated data with two diffusion times, Δ1 and Δ2, to the sub-diffusion model (Equation 3) at various SNR levels.

R2 values were computed by comparing the estimated mean kurtosis with the true kurtosis. R2 contours at the 0.85, 0.90, 0.95, and 0.99 levels have been provided for visualisation purposes.

Scatter plots of simulated K values versus fitted K values for simulated data with different number of diffusion times at various SNR levels.

The simulated data is created using the sub-diffusion model with random normal noise (Equation 10). Blue dots represent kurtosis based on fitting the sub-diffusion model (Equation 3). Orange dots represent kurtosis based on fitting the traditional diffusional kurtosis imaging (DKI) model (Equation 6). Black line is a reference line for R2=1.00, indicating fitted kurtosis values are 100% matching the simulated ones.

Time-dependence in diffusional kurtosis imaging (DKI) metrics using simulated data at different diffusion times (Δ¯).

The true values for (Dβ, β, K) used in the simulations are set to (3 × 10–4, 0.75, 0.8125) to represent white matter (blue dotted lines) and (5 × 10–4, 0.85, 0.4733) to represent grey matter (orange dotted lines). (A) and (B) use data with added random Gaussian noise (SNR = 20) to estimate the parameters DDKI and KDKI. (C) and (D) use noiseless data to obtain estimates for large Δ¯ values. Shaded regions in (A) and (B) represent the 95% confidence intervals of the estimates.

R2 values for the b-value sampling optimisation based on DW-MRI data with SNR = 5, 10, and 20.

The specifically investigated b-value combinations using two, three, and four non-zero b-values have been ordered by the size of the R2 value. The colour bar depicts the proportion of Δ=19 ms and Δ=49 ms b-values needed to produce the corresponding R2 value. Note, the different b-value combinations were assigned a unique identifier, and these appear along the abscissa for each of the three non-zero b-value cases. The b-value combinations achieving the highest R2 values are displayed in the inset pictures and the b-values are provided in Table 1.

Spatially resolved maps of mean kurtosis shown for two example slices and two different subjects, Subject 3 rescan slice 71 (Panel A) and Subject 5 slice 74 (Panel B) from the Connectome 1.0 DW-MRI data.

Individual maps were generated using the sub-diffusion model framework (K), as well as using the traditional approach (KDKI). The diffusion times, Δ, used to generate each plot are provided for each case. We consider the mean kurtosis maps using two diffusion times (Δ=19,49 ms) as the benchmarks.

Representative error map and sample parameter fits for Subject 5 slice 74.

The DW-MRI data with two diffusion times was fitted to the sub-diffusion model in both q-space (A–D) and b-space (E–H), following (Equation 3) and (Equation 4), respectively. The first and second columns are voxels in white matter (30,20,74) and (45,56,74), respectively. The third and fourth columns are voxels in grey matter (58,35,74) and (34,78,74), respectively.

Distributions of the estimated parameter pair (D,K) in different regions of the brain of all subjects, coloured by the probability density.

Yellow indicates high probability density, light blue indicates low probability density. Panel (A): distributions of (D,K), generated using the sub-diffusion model (Equation 3) with both Δ=19,49ms. Panel (B): distributions of (DDKI,KDKI), generated using the standard diffusional kurtosis imaging (DKI) model (Equation 6) with Δ=19ms. Kurtosis is dimensionless and diffusivity is in units of × 10–3 mm2/s.

Spatially resolved maps of mean kurtosis shown for two example slices and two different subjects, Subject 3 rescan slice 71 (Panel A) and Subject 5 slice 74 (Panel B), based on SNR reduction of the Connectome 1.0 DW-MRI data.

Individual maps were generated using the sub-diffusion model framework (K), considering optimal and sub-optimal four non-zero b-value sampling schemes. Here, two b-values with Δ=19 ms and two b-values with Δ=49 ms were selected for each case. The optimal b-values were chosen as the best for each SNR shown in Table 1. The sub-optimal b-values were chosen to have an R2 = 0.3, 0.45, 0.5 to be about half of the maximum R2, for SNR = 6 (b=800,1500,200,2300 s/mm2), SNR=10 (b=1500,3450,6750,13500 s/mm2), and SNR = 20 (b=3450,4750,2300,4250 s/mm2), respectively. The benchmark kurtosis map is provided in Figure 6.

Interclass correlation coefficient (ICC) results for mean kurtosis are depicted for the 12 brain regions analysed.

The mean (μ) and standard deviation (σ) computed based on all the Connectome 1.0 DW-MRI data (subscript A), and the reduced data achieving an SNR = 10 with optimal four non-zero b-value sampling (subscript O), are provided for each brain region. Histograms were generated using all data. Mean kurtosis based on the optimised protocol was computed using the sub-diffusion framework using DW-MRI data with the four non-zero b-values suggested in Table 1 and diffusion encoding directions down sampled to achieve an SNR = 10.

Tables

Table 1
A selection of the best b-value sampling regimes to achieve the highest R2 value in the three cases considered.

The various categories correspond with two, three and four non-zero b-value sampling schemes, with Δ1 and Δ2 denoting the diffusion time setting used to generate the b-values. Note, entries are b-values in unit of s/mm2, and Δ1=19 ms and Δ2=49 ms were used to match the Connectome 1.0 DW-MRI data collection protocol. The entries listed at the bottom row are suggested optimal non-zero b-values for clinical practice.

Two b-valuesThree b-valuesFour b-values
Δ1Δ2R2Δ1Δ2R2Δ1Δ2R2
35067500.01350, 80023000.46350, 4750950, 42500.61
80023000.11350950, 42500.47350, 6000950, 67500.62
SNR = 535042500.183502300, 42500.47350, 4750950, 67500.62
3509500.2650, 80023000.48350, 2400950, 67500.63
35023000.30350950, 67500.49350, 2400950, 42500.63
35042500.63350, 475023000.83350, 2400950, 98500.91
80042500.6450, 80023000.84350, 3450950, 67500.91
SNR = 103509500.67350, 80023000.84350, 6000950, 98500.91
35023000.72350, 150023000.84350, 4750950, 98500.91
80023000.73350, 240023000.85350, 3450950, 98500.91
150023000.77350, 600023000.89350, 47502300, 13,5000.96
3509500.7850, 150042500.90350, 47502300, 67500.96
SNR = 2035023000.80350, 150023000.90350, 60002300, 98500.96
80042500.82350, 475023000.90350, 47502300, 98500.96
80023000.8550, 240042500.90350, 4750950, 67500.96
80023000.85350, 150023000.90350, 1500950, 42500.92
Table 2
Benchmark kurtosis values estimated using the Connectome 1.0 DW-MRI data for different regions of the human brain.

Results are provided for the traditional mean kurtosis (KDKI) at two distinct diffusion times, and values (K) obtained based on fitting the sub-diffusion model across both diffusion times. Results are for grey matter (GM) and white matter (WM) brain regions, in categories of sub-cortical (sc) and cortical (c), and CC stands for corpus callosum. The pooled means and standard deviations across participants have been tabulated, along with the coefficient of variation in parentheses.

KDKIK
Δ=19 msΔ=49 msΔ=19,49 ms
scGM0.57 ± 0.23(40%)0.50 ± 0.24(48%)0.60 ± 0.21(35%)
Thalamus0.65 ± 0.19(30%)0.58 ± 0.22(39%)0.70 ± 0.17(25%)
Caudate0.41 ± 0.24(58%)0.37 ± 0.23(64%)0.39 ± 0.14(35%)
Putamen0.54 ± 0.21(40%)0.45 ± 0.22(49%)0.49 ± 0.13(27%)
Pallidum0.68 ± 0.25(37%)0.64 ± 0.26(41%)0.93 ± 0.18(19%)
cGM0.53 ± 0.24(46%)0.46 ± 0.23(51%)0.40 ± 0.16(39%)
Fusiform0.55 ± 0.22(40%)0.44 ± 0.22(49%)0.40 ± 0.15(37%)
Lingual0.59 ± 0.21(35%)0.53 ± 0.22(42%)0.47 ± 0.16(34%)
WM0.78 ± 0.20(25%)0.76 ± 0.19(26%)0.87 ± 0.22(25%)
Cerebral WM0.77 ± 0.19(25%)0.75 ± 0.19(26%)0.85 ± 0.22(26%)
Cerebellum WM0.99 ± 0.18(18%)0.95 ± 0.19(20%)1.07 ± 0.22(21%)
CC0.65 ± 0.22(35%)0.65 ± 0.22(34%)0.95 ± 0.25(26%)
Table 3
Benchmark Dβ values estimated using the Connectome 1.0 DW-MRI data for different regions of the human brain.

Results are provided for fitting the data at two distinct diffusion times, and fitting the sub-diffusion model across both diffusion times. Results are for grey matter (GM) and white matter (WM) brain regions, in categories of sub-cortical (sc) and cortical (c), and CC stands for corpus callosum. The pooled means and standard deviations across participants have been tabulated, along with the coefficient of variation in parentheses.

Dβ(×104mm2/sβ)
Δ=19 msΔ=49 msΔ=19,49 ms
scGM3.12 ± 0.96(31%)3.56 ± 1.03(29%)4.13 ± 0.86(21%)
Thalamus2.73 ± 0.83(30%)3.14 ± 0.89(28%)3.89 ± 0.81(21%)
Caudate4.32 ± 0.94(22%)4.80 ± 0.89(19%)5.21 ± 1.03(20%)
Putamen3.36 ± 0.58(17%)3.96 ± 0.54(14%)4.31 ± 0.44(10%)
Pallidum1.84 ± 0.63(34%)1.88 ± 0.67(36%)2.79 ± 0.52(19%)
cGM4.17 ± 0.97(23%)5.19 ± 0.96(19%)5.32 ± 0.90(17%)
Fusiform3.94 ± 0.83(21%)4.93 ± 0.74(15%)5.05 ± 0.69(14%)
Lingual3.78 ± 1.00(27%)4.62 ± 0.99(21%)5.01 ± 1.02(20%)
WM1.84 ± 0.79(43%)2.14 ± 0.93(43%)2.94 ± 0.72(24%)
Cerebral WM1.90 ± 0.77(41%)2.20 ± 0.92(42%)2.98 ± 0.71(24%)
Cerebellum WM0.81 ± 0.44(55%)1.17 ± 0.61(52%)2.03 ± 0.49(24%)
CC2.58 ± 1.84(71%)2.49 ± 1.91(76%)3.78 ± 2.12(56%)
Table 4
Benchmark β values estimated using the Connectome 1.0 DW-MRI data for different regions of the human brain.

Results are provided for fitting the data at two distinct diffusion times, and fitting the sub-diffusion model across both diffusion times. Results are for grey matter (GM) and white matter (WM) brain regions, in categories of sub-cortical (sc) and cortical (c), and CC stands for corpus callosum. The pooled means and standard deviations across participants have been tabulated, along with the coefficient of variation in parentheses.

β
Δ=19 msΔ=49 msΔ=19,49 ms
scGM0.74 ± 0.10(13%)0.69 ± 0.14(19%)0.81 ± 0.06(8%)
Thalamus0.70 ± 0.09(13%)0.63 ± 0.12(20%)0.78 ± 0.05(6%)
Caudate0.83 ± 0.07(8%)0.81 ± 0.07(9%)0.88 ± 0.04(5%)
Putamen0.79 ± 0.07(8%)0.77 ± 0.06(8%)0.85 ± 0.04(5%)
Pallidum0.61 ± 0.11(19%)0.47 ± 0.14(30%)0.72 ± 0.05(7%)
cGM0.81 ± 0.08(10%)0.83 ± 0.07(9%)0.87 ± 0.05(5%)
Fusiform0.81 ± 0.08(9%)0.83 ± 0.07(8%)0.87 ± 0.04(5%)
Lingual0.78 ± 0.09(11%)0.78 ± 0.09(12%)0.85 ± 0.05(6%)
WM0.61 ± 0.12(20%)0.52 ± 0.16(32%)0.74 ± 0.06(9%)
Cerebral WM0.62 ± 0.12(19%)0.53 ± 0.16(30%)0.74 ± 0.06(8%)
Cerebellum WM0.41 ± 0.16(39%)0.32 ± 0.19(60%)0.68 ± 0.06(9%)
CC0.61 ± 0.14(23%)0.43 ± 0.18(43%)0.71 ± 0.07(10%)
Table 5
Kurtosis values (K) under the optimal and sub-optimal b-value sampling regimes for specific brain regions.

K was estimated based on fitting the sub-diffusion model to the Connectome 1.0 DW-MRI data with two diffusion times and selected four b-shells. Optimal b-value sampling is considered to have R2 = 0.63, 0.91 and 0.96 for the SNR = 6, 10, and 20 columns, according to Table 1. Sub-optimal b-values are chosen to have R2 = 0.3, 0.45, and 0.5, respectively, as reported in Figure 9. Individual entries are for grey matter (GM) and white matter (WM) brain regions, in categories of sub-cortical (sc) and cortical (c), and CC stands for corpus callosum. A reduction in SNR level was achieved by reducing the number of diffusion encoding directions in each b-shell of the DW-MRI data. The pooled means and standard deviations across participants have been tabulated, along with the coefficient of variation in parentheses. The entries identified in italic under the optimal b-value heading were found to be significantly different from the benchmark mean K reported in Table 2. Sub-optimal result population means were mostly significantly different from the benchmark mean K, and they are not italicised. The average errors (last column) are relative errors compared to the benchmark kurtosis values reported in Table 2.

SNR = 6SNR = 10SNR = 20Average error for SNR = 6/10/20
Optimal b-values
scGM0.47 ± 0.27(57%)0.65 ± 0.22(33%)0.61 ± 0.21(34%)37/23/12%
Thalamus0.57 ± 0.26(46%)0.75 ± 0.20(26%)0.71 ± 0.17(24%)33/20/11%
Caudate0.30 ± 0.20(68%)0.46 ± 0.17(38%)0.40 ± 0.15(38%)43/30/15%
Putamen0.38 ± 0.21(56%)0.58 ± 0.16(28%)0.52 ± 0.14(28%)40/28/13%
Pallidum0.66 ± 0.33(49%)0.90 ± 0.24(26%)0.89 ± 0.20(22%)39/19/12%
cGM0.40 ± 0.22(54%)0.51 ± 0.20(39%)0.45 ± 0.18(40%)30/36/19%
Fusiform0.37 ± 0.21(57%)0.54 ± 0.20(37%)0.45 ± 0.17(38%)35/43/20%
Lingual0.45 ± 0.22(50%)0.59 ± 0.19(33%)0.52 ± 0.18(34%)28/34/18%
WM0.72 ± 0.26(36%)0.91 ± 0.23(26%)0.89 ± 0.22(25%)23/13/8%
Cerebral WM0.71 ± 0.25(35%)0.89 ± 0.23(25%)0.88 ± 0.22(25%)22/13/8%
Sub-optimal b-values
scGM0.31 ± 0.22(71%)0.82 ± 0.29(36)0.58 ± 0.25(43%)55/45/34%
Thalamus0.36 ± 0.23(65%)0.95 ± 0.25(27%)0.72 ± 0.19(26%)54/43/31%
Caudate0.23 ± 0.19(81%)0.57 ± 0.17(31%)0.40 ± 0.20(52%)56/57/40%
Putamen0.25 ± 0.16(67%)0.65 ± 0.17(26%)0.44 ± 0.20(46%)54/42/36%
Pallidum0.44 ± 0.29(66%)1.27 ± 0.32(25%)0.76 ± 0.25(33%)56/42/32%
cGM0.25 ± 0.17(67%)0.44 ± 0.14(33%)0.37 ± 0.17(47%)47/30/36%
Fusiform0.22 ± 0.15(68%)0.47 ± 0.16(35%)0.33 ± 0.17(52%)51/32/38%
Lingual0.29 ± 0.19(66%)0.53 ± 0.16(30%)0.39 ± 0.20(51%)48/29/40%
WM0.42 ± 0.18(44%)1.16 ± 0.36(31%)0.72 ± 0.22(31%)53/38/28%
Cerebral WM0.41 ± 0.18(43%)1.15 ± 0.35(31%)0.74 ± 0.20(28%)53/38/27%
Cerebellum WM0.47 ± 0.21(45%)1.27 ± 0.33(26%)0.35 ± 0.26(74%)57/24/72%
CC0.69 ± 0.42(61%)1.94 ± 0.44(23%)0.92 ± 0.19(20%)47/102/28%
Cerebellum WM0.83 ± 0.29(35%)1.17 ± 0.25(21%)1.15 ± 0.23(20%)27/15/11%
CC0.88 ± 0.36(41%)0.98 ± 0.31(31%)0.86 ± 0.25(29%)26/17/13%

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Megan E Farquhar
  2. Qianqian Yang
  3. Viktor Vegh
(2024)
Robust, fast and accurate mapping of diffusional mean kurtosis
eLife 12:RP90465.
https://doi.org/10.7554/eLife.90465.3