Quasi-diffusion magnetic resonance imaging (QDI): A fast, high b-value diffusion imaging technique

Superficial white matter microstructural imaging method based on time-space fractional-order diffusion

Objective. Microstructure imaging based on diffusion magnetic resonance signal is an advanced imaging technique that enablesin vivomapping of the brain's microstructure. Superficial white matter (SWM) plays an important role in brain development, maturation, and aging, while fewer microstructure imaging methods address the SWM due to its complexity. Therefore, this study aims to develop a diffusion propagation model to investigate the microstructural characteristics of the SWM region.Approach. In this paper, we hypothesize that the effect of cell membrane permeability and the water exchange between soma and dendrites cannot be neglected for typical clinical diffusion times (20 ms Collapse

Key Words

• anomalous diffusion
• diffusion magnetic resonance imaging
• fractional diffusion equation
• microstructure imaging
• superficial white matter

Collapse

MESH Headings Collapse

Grants Collapse

Affiliation(s)

Jianglin He School of Health Science and Engineering, University of Shanghai for Science and Technology, Shanghai 200093, People's Republic of China
Yuanjun Wang School of Health Science and Engineering, University of Shanghai for Science and Technology, Shanghai 200093, People's Republic of China

Collapse

Identifying potential imaging markers for diffusion property changes in a mouse model of amyotrophic lateral sclerosis: Application of the continuous time random walk model to ultrahigh b-value diffusion-weighted MR images of spinal cord tissue

Diffusion MRI (dMRI) explores tissue microstructures by analyzing diffusion-weighted signal decay measured at different b-values. While relatively low b-values are used for most dMRI models, high b-value diffusion-weighted imaging (DWI) techniques have gained interest given that the non-Gaussian water diffusion behavior observed at high b-values can yield potentially valuable information. In this study, we investigated anomalous diffusion behaviors associated with degeneration of spinal cord tissue using a continuous time random walk (CTRW) model for DWI data acquired across an extensive range of ultrahigh b-values. The diffusion data were acquired in situ from the lumbar level of spinal cords of wild-type and age-matched transgenic SOD1G93A mice, a well-established animal model of amyotrophic lateral sclerosis (ALS) featuring progressive degeneration of axonal tracts in this tissue. Based on the diffusion decay behaviors at low and ultrahigh b-values, we applied the CTRW model using various combinations of b-values and compared diffusion metrics calculated from the CTRW model between the experimental groups. We found that diffusion-weighted signal decay curves measured with ultrahigh b-values (up to 858,022 s/mm2 in this study) were well represented by the CTRW model. The anomalous diffusion coefficient obtained from lumbar spinal cords was significantly higher in SOD1G93A mice compared with control mice (14.7 × 10-5  ± 5.54 × 10-5  vs. 7.87 × 10-5  ± 2.48 × 10-5  mm2 /s, p = 0.01). We believe this is the first study to illustrate the efficacy of the CTRW model for analyzing anomalous diffusion regimes at ultrahigh b-values. The CTRW modeling of ultrahigh b-value dMRI can potentially present a novel approach for noninvasively evaluating alterations in spinal cord tissue associated with ALS pathology.

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.

Generalisation of continuous time random walk to anomalous diffusion MRI models with an age-related evaluation of human corpus callosum

Diffusion MRI measures of the human brain provide key insight into microstructural variations across individuals and into the impact of central nervous system diseases and disorders. One approach to extract information from diffusion signals has been to use biologically relevant analytical models to link millimetre scale diffusion MRI measures with microscale influences. The other approach has been to represent diffusion as an anomalous transport process and infer microstructural information from the different anomalous diffusion equation parameters. In this study, we investigated how parameters of various anomalous diffusion models vary with age in the human brain white matter, particularly focusing on the corpus callosum. We first unified several established anomalous diffusion models (the super-diffusion, sub-diffusion, quasi-diffusion and fractional Bloch-Torrey models) under the continuous time random walk modelling framework. This unification allows a consistent parameter fitting strategy to be applied from which meaningful model parameter comparisons can be made. We then provided a novel way to derive the diffusional kurtosis imaging (DKI) model, which is shown to be a degree two approximation of the sub-diffusion model. This link between the DKI and sub-diffusion models led to a new robust technique for generating maps of kurtosis and diffusivity using the sub-diffusion parameters β_SUB and D_SUB. Superior tissue contrast is achieved in kurtosis maps based on the sub-diffusion model. 7T diffusion weighted MRI data for 65 healthy participants in the age range 19-78 years was used in this study. Results revealed that anomalous diffusion model parameters α and β have shown consistent positive correlation with age in the corpus callosum, indicating α and β are sensitive to tissue microstructural changes in ageing.

The Mathematics of Quasi-Diffusion Magnetic Resonance Imaging

Quasi-diffusion imaging (QDI) is a novel quantitative diffusion magnetic resonance imaging (dMRI) technique that enables high quality tissue microstructural imaging in a clinically feasible acquisition time. QDI is derived from a special case of the continuous time random walk (CTRW) model of diffusion dynamics and assumes water diffusion is locally Gaussian within tissue microstructure. By assuming a Gaussian scaling relationship between temporal (α) and spatial (β) fractional exponents, the dMRI signal attenuation is expressed according to a diffusion coefficient, D (in mm2 s−1), and a fractional exponent, α. Here we investigate the mathematical properties of the QDI signal and its interpretation within the quasi-diffusion model. Firstly, the QDI equation is derived and its power law behaviour described. Secondly, we derive a probability distribution of underlying Fickian diffusion coefficients via the inverse Laplace transform. We then describe the functional form of the quasi-diffusion propagator, and apply this to dMRI of the human brain to perform mean apparent propagator imaging. QDI is currently unique in tissue microstructural imaging as it provides a simple form for the inverse Laplace transform and diffusion propagator directly from its representation of the dMRI signal. This study shows the potential of QDI as a promising new model-based dMRI technique with significant scope for further development.

Diffusion in Sephadex Gel Structures: Time Dependency Revealed by Multi-Sequence Acquisition over a Broad Diffusion Time Range

It has been increasingly reported that in biological tissues diffusion-weighted MRI signal attenuation deviates from mono-exponential decay, especially at high b-values. A number of diffusion models have been proposed to characterize this non-Gaussian diffusion behavior. One of these models is the continuous-time random-walk (CTRW) model, which introduces two new parameters: a fractional order time derivative α and a fractional order spatial derivative β. These new parameters have been linked to intravoxel diffusion heterogeneities in time and space, respectively, and are believed to depend on diffusion times. Studies on this time dependency are limited, largely because the diffusion time cannot vary over a board range in a conventional spin-echo echo-planar imaging sequence due to the accompanying T2 decays. In this study, we investigated the time-dependency of the CTRW model in Sephadex gel phantoms across a broad diffusion time range by employing oscillating-gradient spin-echo, pulsed-gradient spin-echo, and pulsed-gradient stimulated echo sequences. We also performed Monte Carlo simulations to help understand our experimental results. It was observed that the diffusion process fell into the Gaussian regime at extremely short diffusion times whereas it exhibited a strong time dependency in the CTRW parameters at longer diffusion times.

Quartile histogram assessment of glioma malignancy using high b-value diffusion MRI with a continuous-time random-walk model

The purpose of this study is to investigate the feasibility of using a continuous-time random-walk (CTRW) diffusion model, together with a quartile histogram analysis, for assessing glioma malignancy by probing tissue heterogeneity as well as cellularity. In this prospective study, 91 patients (40 females, 51 males) with histopathologically proven gliomas underwent MRI at 3 T. The cohort included 42 grade II (GrII), 19 grade III (GrIII) and 29 grade IV (GrIV) gliomas. Echo-planar diffusion-weighted imaging was conducted using 17 b-values (0-4000 s/mm² ). Three CTRW model parameters, including an anomalous diffusion coefficient D_m , and two parameters related to temporal and spatial diffusion heterogeneity α and β, respectively, were obtained. The mean parameter values within the tumor regions of interest (ROIs) were computed by utilizing the first quartile of the histograms as well as the full ROI for comparison. A Bonferroni-Holm-corrected Mann-Whitney U-test was used for the group comparisons. Individual and combinations of the CTRW parameters were evaluated for the characterization of gliomas with a receiver operating characteristic analysis. All first-quartile mean CTRW parameters yielded significant differences (p-values < 0.05) between pair-wise comparisons of GrII (D_m : 1.14 ± 0.37 μm² /ms; α: 0.904 ± 0.03, β: 0.913 ± 0.06), GrIII (D_m : 0.88 ± 0.21 μm² /ms; α: 0.888 ± 0.01, β: 0.857 ± 0.06) and GrIV gliomas (D_m : 0.73 ± 0.22 μm² /ms; α: 0.878 ± 0.01; β: 0.791 ± 0.07). The highest sensitivity, specificity, accuracy and area-under-the-curve of using the combinations of the first-quartile parameters were 84.2%, 78.5%, 75.4% and 0.76 for GrII and GrIII classification; 86.2%, 89.4%, 75% and 0.76 for GrIII and GrIV classification; and 86.2%, 85.7%, 84.5% and 0.90 for GrII and GrIV classification, respectively. Quartile-based analysis produced higher accuracy and area-under-the-curve than the full ROI-based analysis in all classifications. The CTRW diffusion model, together with a quartile-based histogram analysis, offers a new way for probing tumor structural heterogeneity at a subvoxel level, and has potential for in vivo assessment of glioma malignancy to complement histopathology.

In vivo assessment of Lauren classification for gastric adenocarcinoma using diffusion MRI with a fractional order calculus model

OBJECTIVES

To evaluate the performance of a fractional order calculus (FROC) diffusion model for imaging-based assessment of Lauren classification in gastric adenocarcinoma.

METHODS

In this study, 43 patients (15 females, 28 males) with gastric adenocarcinoma underwent MRI at 1.5 T. According to pathology-based Lauren classification, 10 patients had diffuse-type, 20 had intestinal-type, and 13 had mixed-type lesions. The diffuse and mixed types were combined as diffuse-and-mixed type to be differentiated from the intestinal type using diffusion MRI. Diffusion-weighted images were acquired by using eleven b-values (0-2000 s/mm²). Three FROC model parameters comprising diffusion coefficient D, intravoxel diffusion heterogeneity β, and a microstructural quantity μ, together with a conventional apparent diffusion coefficient (ADC), were estimated. The mean parameter values in the tumour were computed by using a percentile histogram analysis. Individual or linear combinations of the mean parameters in the tumour were used to differentiate the diffuse-and-mixed type from the intestinal type using descriptive statistics and receiver operating characteristic (ROC) analyses.

RESULTS

Significant differences were observed between diffuse-and-mixed-type and intestinal-type lesions in D (0.99 ± 0.20 μm²/ms vs. 1.11 ± 0.23 μm²/ms; p = 0.036), β (0.37 ± 0.08 vs. 0.43 ± 0.11; p = 0.043), μ (7.92 ± 2.79 μm vs. 9.87 ± 1.52 μm; p = 0.038), and ADC (0.81 ± 0.34 μm²/ms vs. 0.96 ± 0.19 μm²/ms; p = 0.033). Among the individual parameters, μ produced the largest area under the ROC curve (0.739). The combinations of (D, β, μ) and (β and μ) produced the best overall performance with a sensitivity of 0.739, specificity of 0.750, accuracy of 0.744, and area under the curve of 0.793 (95% confidence interval: 0.657-0.929).

CONCLUSION

Diffusion MRI with the FROC model holds promise for non-invasive assessment of Lauren classification for gastric adenocarcinoma.

KEY POINTS

• High b-value diffusion MRI with a FROC model that is sensitive to tissue microstructures can differentiate the diffuse-and-mixed type from intestinal type of gastric adenocarcinoma. • The combination of FROC parameters produced the best result for distinguishing the diffuse-and-mixed type from the intestinal type with an area under the receiver operating characteristic curve of 0.793. • The FROC model parameters, individually or conjointly, hold promise for repeated, non-invasive evaluations of gastric adenocarcinoma at various time points throughout disease progression or regression to complement conventional Lauren classification.