A fast schema for parameter estimation in diffusion kurtosis imaging

Optimization of quasi-diffusion magnetic resonance imaging for quantitative accuracy and time-efficient acquisition

PURPOSE

Quasi-diffusion MRI (QDI) is a novel quantitative technique based on the continuous time random walk model of diffusion dynamics. QDI provides estimates of the diffusion coefficient, <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mspace/> <mml:msub><mml:mi>D</mml:mi> <mml:mrow><mml:mn>1</mml:mn> <mml:mo>,</mml:mo> <mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {D}_{1,2} $$</mml:annotation></mml:semantics> </mml:math> in mm²  s^-1 and a fractional exponent, <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>α</mml:mi></mml:mrow> <mml:annotation>$$ \upalpha $$</mml:annotation></mml:semantics> </mml:math> , defining the non-Gaussianity of the diffusion signal decay. Here, the b-value selection for rapid clinical acquisition of QDI tensor imaging (QDTI) data is optimized.

METHODS

Clinically appropriate QDTI acquisitions were optimized in healthy volunteers with respect to a multi-b-value reference (MbR) dataset comprising 29 diffusion-sensitized images arrayed between <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>b</mml:mi> <mml:mo>=</mml:mo> <mml:mn>0</mml:mn></mml:mrow> <mml:annotation>$$ b=0 $$</mml:annotation></mml:semantics> </mml:math> and 5000 s mm^-2 . The effects of varying maximum b-value ( <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>b</mml:mi> <mml:mi>max</mml:mi></mml:msub> </mml:mrow> <mml:annotation>$$ {b}_{\mathrm{max}} $$</mml:annotation></mml:semantics> </mml:math> ), number of b-value shells, and the effects of Rician noise were investigated.

RESULTS

QDTI measures showed <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>b</mml:mi> <mml:mi>max</mml:mi></mml:msub> </mml:mrow> <mml:annotation>$$ {b}_{\mathrm{max}} $$</mml:annotation></mml:semantics> </mml:math> dependence, most significantly for <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>α</mml:mi></mml:mrow> <mml:annotation>$$ \upalpha $$</mml:annotation></mml:semantics> </mml:math> in white matter, which monotonically decreased with higher <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>b</mml:mi> <mml:mi>max</mml:mi></mml:msub> </mml:mrow> <mml:annotation>$$ {b}_{\mathrm{max}} $$</mml:annotation></mml:semantics> </mml:math> leading to improved tissue contrast. Optimized 2 b-value shell acquisitions showed small systematic differences in QDTI measures relative to MbR values, with overestimation of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mspace/> <mml:mspace/> <mml:msub><mml:mi>D</mml:mi> <mml:mrow><mml:mn>1</mml:mn> <mml:mo>,</mml:mo> <mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ \kern0.50em {D}_{1,2} $$</mml:annotation></mml:semantics> </mml:math> and underestimation of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>α</mml:mi></mml:mrow> <mml:annotation>$$ \upalpha $$</mml:annotation></mml:semantics> </mml:math> in white matter, and overestimation of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>D</mml:mi> <mml:mrow><mml:mn>1</mml:mn> <mml:mo>,</mml:mo> <mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {D}_{1,2} $$</mml:annotation></mml:semantics> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>α</mml:mi></mml:mrow> <mml:annotation>$$ \upalpha $$</mml:annotation></mml:semantics> </mml:math> anisotropies in gray and white matter. Additional shells improved the accuracy, precision, and reliability of QDTI estimates with 3 and 4 shells at <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>b</mml:mi> <mml:mi>max</mml:mi></mml:msub> <mml:mo>=</mml:mo> <mml:mn>5000</mml:mn></mml:mrow> <mml:annotation>$$ {b}_{\mathrm{max}}=5000 $$</mml:annotation></mml:semantics> </mml:math>  s mm^-2 , and 4 b-value shells at <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>b</mml:mi> <mml:mi>max</mml:mi></mml:msub> <mml:mo>=</mml:mo> <mml:mn>3960</mml:mn></mml:mrow> <mml:annotation>$$ {b}_{\mathrm{max}}=3960 $$</mml:annotation></mml:semantics> </mml:math>  s mm^-2 , providing minimal bias in <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>D</mml:mi> <mml:mrow><mml:mn>1</mml:mn> <mml:mo>,</mml:mo> <mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {D}_{1,2} $$</mml:annotation></mml:semantics> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>α</mml:mi></mml:mrow> <mml:annotation>$$ \upalpha $$</mml:annotation></mml:semantics> </mml:math> compared to the MbR.

CONCLUSION

A highly detailed optimization of non-Gaussian dMRI for in vivo brain imaging was performed. QDI provided robust parameterization of non-Gaussian diffusion signal decay in clinically feasible imaging times with high reliability, accuracy, and precision of QDTI measures.

White matter microstructural impairments in amyotrophic lateral sclerosis: A mean apparent propagator MRI study

Background

White matter (WM) impairment is a hallmark of amyotrophic lateral sclerosis (ALS). This study evaluated the capacity of mean apparent propagator magnetic resonance imaging (MAP-MRI) for detecting ALS-related WM alterations.

Methods

Diffusion images were obtained from 52 ALS patients and 51 controls. MAP-derived indices [return-to-origin/-axis/-plane probability (RTOP/RTAP/RTPP) and non-Gaussianity (NG)/perpendicular/parallel NG (NG_⊥/NG_||)] were computed. Measures from diffusion tensor/kurtosis imaging (DTI/DKI) and neurite orientation dispersion and density imaging (NODDI) were also obtained. Voxel-wise analysis (VBA) was performed to determine differences in these parameters. Relationship between MAP parameters and disease severity (assessed by the revised ALS Functional Rating Scale (ALSFRS-R)) was evaluated by Pearson’s correlation analysis in a voxel-wise way. ALS patients were further divided into two subgroups: 29 with limb-only involvement and 23 with both bulbar and limb involvement. Subgroup analysis was then conducted to investigate diffusion parameter differences related to bulbar impairment.

Results

The VBA (with threshold of P < 0.05 after family-wise error correction (FWE)) showed that ALS patients had significantly decreased RTOP/RTAP/RTPP and NG/ NG_⊥/NG_|| in a set of WM areas, including the bilateral precentral gyrus, corona radiata, posterior limb of internal capsule, midbrain, middle corpus callosum, anterior corpus callosum, parahippocampal gyrus, and medulla. MAP-MRI had the capacity to capture WM damage in ALS, which was higher than DTI and similar to DKI/NODDI. RTOP/RTAP/NG/NG_⊥/NG_|| parameters, especially in the bilateral posterior limb of internal capsule and middle corpus callosum, were significantly correlated with ALSFRS-R (with threshold of FWE-corrected P < 0.05). The VBA (with FWE-corrected P < 0.05) revealed the significant RTAP reduction in subgroup with both bulbar and limb involvement, compared with those with limb-only involvement.

Conclusions

Microstructural impairments in corticospinal tract and corpus callosum represent the consistent characteristic of ALS. MAP-MRI could provide alternative measures depicting ALS-related WM alterations, complementary to the common diffusion imaging methods.

MRI phantom for tissue simulation with respect to diffusion coefficient and kurtosis - Validation with injection of liposomal theranostics

This study presents gelatine-based and agar-based phantoms with an addition of glycerol, safflower oil, silicone oil and cellulose microcrystalline with a potential to cover the entire range of tissue diffusion coefficients and kurtosis values. Forty types of phantoms were prepared and examined for NMR relaxation times T₁ and T₂ and diffusional metrics D, K and ADC. Wide ranges of values of D (0.0003-0.0031 mm²s^-1), K (0.00-7.24) and ADC (0.0002-0.0031 mm²s^-1) were observed. Two of the phantoms closely mimic muscle and cortical gray matter with respect to water diffusion parameters. Although many of the presented phantoms display both D and K values within the range of human tissues, they match different tissues with respect to D and K. The imaging results for the gray matter simulating phantom injected with the liposomal solution, bear a resemblance to the particle size effect described in the literature. The phantoms presented in this work are simple in preparation and affordable tissue-simulating materials to be used primarily in development of diffusion kurtosis-based MRI methods and possibly in a preliminary assessment of MRI contrast agents. Further adjustments of the chemical compositions could potentially lead to development of new types of phantoms mimicking diffusional properties of more kinds of soft tissues.

Evaluation of non-local means based denoising filters for diffusion kurtosis imaging using a new phantom

Image denoising has a profound impact on the precision of estimated parameters in diffusion kurtosis imaging (DKI). This work first proposes an approach to constructing a DKI phantom that can be used to evaluate the performance of denoising algorithms in regard to their abilities of improving the reliability of DKI parameter estimation. The phantom was constructed from a real DKI dataset of a human brain, and the pipeline used to construct the phantom consists of diffusion-weighted (DW) image filtering, diffusion and kurtosis tensor regularization, and DW image reconstruction. The phantom preserves the image structure while minimizing image noise, and thus can be used as ground truth in the evaluation. Second, we used the phantom to evaluate three representative algorithms of non-local means (NLM). Results showed that one scheme of vector-based NLM, which uses DWI data with redundant information acquired at different b-values, produced the most reliable estimation of DKI parameters in terms of Mean Square Error (MSE), Bias and standard deviation (Std). The result of the comparison based on the phantom was consistent with those based on real datasets.