Modified triexponential analysis of intravoxel incoherent motion for brain perfusion and diffusion

Presentation of microstructural diffusion components by color schemes in abdominal organs

PURPOSE

Development of a color scheme representation to facilitate the interpretation of tri-exponential DWI data from abdominal organs, where multi-exponential behavior is more pronounced.

METHODS

Multi-exponential analysis of DWI data provides information about the microstructure of the tissue under study. The tri-exponential signal analysis generates numerous parameter images that are difficult to analyze individually. Summarized color images can simplify at-a-glance analysis. A color scheme was developed in which the slow, intermediate, and fast diffusion components were each assigned to a different red, green, and blue color channel. To improve the appearance of the image, histogram equalization, gamma correction, and white balance were used, and the processing parameters were adjusted. Examples of the resulting color maps of the diffusion fractions of healthy and pathological kidney and prostate are shown.

RESULTS

The color maps obtained by the presented method show the merged information of the slow, intermediate, and fast diffusion components in a single view. A differentiation of the different fractions becomes clearly visible. Fast diffusion regimes, such as in the renal hilus, can be clearly distinguished from slow fractions, such as in dense tumor tissue.

CONCLUSION

Combining the diffusion information from tri-exponential DWI analysis into a single color image allows for simplified interpretation of the diffusion fractions. In the future, such color images may provide additional information about the microstructural nature of the tissue under study.

Evaluation of deep learning reconstruction on diffusion-weighted imaging quality and apparent diffusion coefficient using an ice-water phantom

This study assessed the influence of deep learning reconstruction (DLR) on the quality of diffusion-weighted images (DWI) and apparent diffusion coefficient (ADC) using an ice-water phantom. An ice-water phantom with known diffusion properties (true ADC = 1.1 × 10-3 mm2/s at 0 °C) was imaged at various b-values (0, 1000, 2000, and 4000 s/mm2) using a 3 T magnetic resonance imaging scanner with slice thicknesses of 1.5 and 3.0 mm. All DWIs were reconstructed with or without DLR. ADC maps were generated using combinations of b-values 0 and 1000, 0 and 2000, and 0 and 4000 s/mm2. Based on the quantitative imaging biomarker alliance profile, the signal-to-noise ratio (SNRs) in DWIs was calculated, and the accuracy, precision, and within-subject parameter variance (wCV) of the ADCs were evaluated. DLR improved the SNR in DWIs with b-values ranging from 0 to 2000s/mm2; however, its effectiveness was diminished at 4000 s/mm2. There was no noticeable difference in the ADCs of images generated with or without implementing DLR. For a slice thickness of 1.5 mm and combined b-values of 0 and 4000 s/mm2, the ADC values were 0.97 × 10-3and 0.98 × 10-3mm2/s with and without DLR, respectively, both being lower than the true ADC value. Furthermore, DLR enhanced the precision and wCV of the ADC measurements. DLR can enhance the SNR, repeatability, and precision of ADC measurements; however, it does not improve their accuracies.

Tri- and bi-exponential diffusion analyses of the kidney: effect of respiratory-controlled acquisition on diffusion parameters

This study examined whether respiratory-controlled acquisition influences diffusion parameters obtained with intravoxel incoherent motion (IVIM) analysis using tri-exponential and bi-exponential models. Ten healthy volunteers were examined on a 3.0 T MRI system to obtain coronal diffusion-weighted images of both kidneys. The participants were scanned twice using respiratory-triggering (RT) and free-breathing (FB) acquisition to assess the repeatability of the measurements. We determined mean signal intensities in the renal cortex at each b value. Then, perfusion-related diffusion coefficient (Dp), fast-free diffusion coefficient (Df), slow-restricted diffusion coefficient (Ds), and their corresponding fractions (Fp, Ff, and Fs, respectively) were calculated using tri-exponential function. Moreover, perfusion-related diffusion coefficient (D*), the fraction (F), and perfusion-independent diffusion coefficient (D) were calculated using bi-exponential function. Normalized root-mean-square errors for the tri- and bi-exponential analyses (nRMSEtri and nRMSEbi, respectively) were determined to assess the deviation of the fitted to measured data, i.e., the fitting accuracy. Additionally, repeatability coefficients (RCs) were calculated from Bland-Altman plots to evaluate the repeatability of each diffusion parameter. These values were compared between the RT and FB groups. Dp and D* in the RT group were significantly lower than those in the FB group (P < 0.05). In addition, the RT group showed significantly lower nRMSEtri and nRMSEbi values than those in the FB group (P < 0.05). Moreover, Dp, Ds, Fs, and D* at RT showed lower RC values than those at FB. Respiratory-controlled acquisition affects perfusion-related diffusion parameters of the kidney obtained using tri-exponential and bi-exponential analyses.

A three-component multi-b-value diffusion-weighted imaging might be a useful biomarker for detecting microstructural features in gliomas with differences in malignancy and IDH-1 mutation status

OBJECTIVES

The purpose of the study was to explore the performance of a three-component diffusion model in evaluating the degree of malignancy and isocitrate dehydrogenase 1 (IDH-1) gene type of gliomas.

METHODS

Overall, 60 patients with gliomas were enrolled. The intermediate and perfusion-related diffusion coefficients (D_int and D_p) and fractions of strictly limited, intermediate, and perfusion-related diffusion (F_very-slow, F_int, and F_p) were obtained with a three-component diffusion model. Parameters were also obtained from a diffusion kurtosis model and mono- and biexponential models. All parameters were compared between different tumor grades and IDH-1 gene types. Diagnostic performance and logistic regression analyses were performed.

RESULTS

High-grade gliomas (HGGs) had significantly higher F_int, F_very-slow, and D_p values but significantly lower F_p and D_int values than low-grade gliomas (LGGs), and F_int and F_p differed significantly among grade II, III, and IV gliomas (p < 0.05 for all). F_int achieved the highest AUC of 0.872 in differentiating between LGGs and HGGs. Logistic regression analysis revealed that in each model, F_int, diffusion coefficient (D), apparent diffusion coefficient (ADC), mean diffusivity (MD), and mean kurtosis (MK) were associated with glioma grading. After multiple regression analysis, F_int remained the only differentiator. Additionally, F_int and F_p showed significant differences between IDH-1 mutated and IDH-1 wild-type gliomas (p = 0.007 and 0.01, respectively).

CONCLUSIONS

The three-component DWI model served as a useful biomarker for detecting microstructural features in gliomas with different grades and IDH-1 mutation statuses.

KEY POINTS

• The three-component model enables the estimation of an intermediate diffusion component. • The three-component model performed better than the other models in glioma grading and genotyping. • F_int was useful in evaluating the grade and genotype of gliomas.

Precision of region of interest-based tri-exponential intravoxel incoherent motion quantification and the role of the Intervoxel spatial distribution of flow velocities

PURPOSE

The purpose of this work was to obtain precise tri-exponential intravoxel incoherent motion (IVIM) quantification in the liver using 2D (b-value and first-order motion moment [M1 ]) IVIM-DWI acquisitions and region of interest (ROI)-based fitting techniques.

METHODS

Diffusion MRI of the liver was performed in 10 healthy volunteers using three IVIM-DWI acquisitions: conventional monopolar, optimized monopolar, and optimized 2D (b-M1 ). For each acquisition, bi-exponential and tri-exponential full, segmented, and over-segmented ROI-based fitting and a newly proposed blood velocity SDdistribution (BVD) fitting technique were performed to obtain IVIM estimates in the right and left liver lobes. Fitting quality was evaluated using corrected Akaike information criterion. Precision metrics (test-retest repeatability, inter-reader reproducibility, and inter-lobar agreement) were evaluated using Bland-Altman analysis, repeatability/reproducibility coefficients (RPCs), and paired sample t-tests. Precision was compared across acquisitions and fitting methods.

RESULTS

High repeatability and reproducibility was observed in the estimations of the diffusion coefficient (Dtri  = [1.03 ± 0.11] × 10-3  mm2 /s; RPCs ≤ 1.34 × 10-4  mm2 /s), perfusion fractions (F1  = 3.19 ± 1.89% and F2  = 16.4 ± 2.07%; RPCs ≤ 2.51%), and blood velocity SDs (Vb,1  = 1.44 ± 0.14 mm/s and Vb,2  = 3.62 ± 0.13 mm/s; RPCs ≤ 0.41 mm/s) in the right liver lobe using the 2D (b-M1 ) acquisition in conjunction with BVD fitting. Using these methods, significantly larger (p < 0.01) estimates of Dtri and F1 were observed in the left lobe in comparison to the right lobe, while estimates of Vb,1 and Vb,2 demonstrated high interlobar agreement (RPCs ≤ 0.45 mm/s).

CONCLUSIONS

The 2D (b-M1 ) IVIM-DWI data acquisition in conjunction with BVD fitting enables highly precise tri-exponential IVIM quantification in the right liver lobe.

[Feasibility of Cerebrovascular Reserve Assessment Using Stretched Exponential Model in Major Cerebral Artery Steno-occlusive Disease: Comparison with SPECT]

PURPOSE

To clarify whether diffusion-weighted imaging using stretched exponential model can assess cerebrovascular reserve (CVR) in patients with major cerebral artery steno-occlusive disease, we compared stretched exponential parameters and single-photon emission computed tomography (SPECT).

METHODS

Twenty-nine patients with unilateral major cerebral artery steno-occlusive disease (25 men and 4 women; age, 69±11 years) were analyzed in this study. The patients were divided into three groups: normal CVR (CVR≥30%), moderate CVR (10%≤CVR<30%), and severe CVR (CVR<10%). The distributed diffusion coefficient (DDC) and heterogeneity index (α) from the stretched exponential model, apparent diffusion coefficient (ADC) from the monoexponential model, and CVR and resting cerebral blood flow (CBF) from SPECT were measured in the bilateral middle cerebral artery territories, and ipsilateral-to-contralateral ratios (rDDC, rα, rADC, and rCBF) were obtained.

RESULTS

The rDDC values in severe CVR were significantly higher than those in normal CVR (P=0.003). The rDDC values were significantly negatively correlated with ipsilateral CVR (rho=-0.31, P=0.009). The rDDC values were not significantly correlated with rCBF (P=0.34).

CONCLUSION

We have shown that elevated rDDC values are associated with impaired CVR. Our results suggest that diffusion-weighted imaging using stretched exponential model has a potential to evaluate hemodynamic impairment.

b value and first-order motion moment optimized data acquisition for repeatable quantitative intravoxel incoherent motion DWI

PURPOSE

To design a b value and first-order motion moment (M1 ) optimized data acquisition for repeatable intravoxel incoherent motion (IVIM) quantification in the liver.

METHODS

Cramer-Rao lower bound optimization was performed to determine optimal monopolar and optimal 2D samplings of the b-M1 space based on noise performance. Monte Carlo simulations were used to evaluate the bias and variability in estimates obtained using the proposed optimal samplings and conventional monopolar sampling. Diffusion MRI of the liver was performed in 10 volunteers using 3 IVIM acquisitions: conventional monopolar, optimized monopolar, and b-M1 -optimized gradient waveforms (designed based on the optimal 2D sampling). IVIM parameter maps of diffusion coefficient, perfusion fraction, and blood velocity SD were obtained using nonlinear least squares fitting. Noise performance (SDs), stability (outlier percentage), and test-retest or scan-rescan repeatability (intraclass correlation coefficients) were evaluated and compared across acquisitions.

RESULTS

Cramer-Rao lower bound and Monte Carlo simulations demonstrated improved noise performance of the optimal 2D sampling in comparison to monopolar samplings. Evaluating the designed b-M1 -optimized waveforms in healthy volunteers, significant decreases (p < 0.05) in the SDs and outlier percentages were observed for measurements of diffusion coefficient, perfusion fraction, and blood velocity SD in comparison to measurements obtained using monopolar samplings. Good-to-excellent repeatability (intraclass correlation coefficients ≥ 0.77) was observed for all 3 parameters in both the right and left liver lobes using the b-M1 -optimized waveforms.

CONCLUSIONS

2D b-M1 -optimized data acquisition enables repeatable IVIM quantification with improved noise performance. 2D acquisitions may advance the establishment of IVIM quantitative biomarkers for liver diseases.

Triexponential Diffusion Analysis of Diffusion-weighted Imaging for Breast Ductal Carcinoma in Situ and Invasive Ductal Carcinoma

Purpose

To obtain detailed information in breast ductal carcinoma in situ (DCIS) and invasive ductal carcinoma (IDC) using triexponential diffusion analysis.

Methods

Diffusion-weighted images (DWI) of the breast were obtained using single-shot diffusion echo-planar imaging with 15 b-values. Mean signal intensities at each b-value were measured in the DCIS and IDC lesions and fitted with the triexponential function based on a two-step approach: slow-restricted diffusion coefficient (D_s) was initially determined using a monoexponential function with b-values > 800 s/mm². The diffusion coefficient of free water at 37°C was assigned to the fast-free diffusion coefficient (D_f). Finally, the perfusion-related diffusion coefficient (D_p) was derived using all the b-values. Furthermore, biexponential analysis was performed to obtain the perfusion-related diffusion coefficient (D^*) and the perfusion-independent diffusion coefficient (D). Monoexponential analysis was performed to obtain the apparent diffusion coefficient (ADC). The sensitivity and specificity of the aforementioned diffusion coefficients for distinguishing between DCIS and IDC were evaluated using the pathological results.

Results

The D_s, D, and ADC of DCIS were significantly higher than those of IDC (P < 0.01 for all). There was no significant correlation between D_p and D_s, but there was a weak correlation between D^* and D. The combination of D_p and D_s showed higher sensitivity and specificity (85.9% and 71.4%, respectively), compared to the combination of D^* and D (81.5% and 33.3%, respectively).

Conclusion

Triexponential analysis can provide detailed diffusion information for breast tumors that can be used to differentiate between DCIS and IDC.

Noncontrast Pediatric Brain Perfusion: Arterial Spin Labeling and Intravoxel Incoherent Motion

Noncontrast magnetic resonance imaging techniques for measuring brain perfusion include arterial spin labeling (ASL) and intravoxel incoherent motion (IVIM). These techniques provide noninvasive and repeatable assessment of cerebral blood flow or cerebral blood volume without the need for intravenous contrast. This article discusses the technical aspects of ASL and IVIM with a focus on normal physiologic variations, technical parameters, and artifacts. Multiple pediatric clinical applications are presented, including tumors, stroke, vasculopathy, vascular malformations, epilepsy, migraine, trauma, and inflammation.

Quantification of Regional Cerebral Blood Flow Using Diffusion Imaging With Phase Contrast

BACKGROUND

The perfusion-related diffusion coefficient obtained from triexponential diffusion analysis is closely correlated with regional cerebral blood flow (rCBF), as assessed by arterial spin labeling (ASL) methods. However, this provides only a semiquantitative measure of rCBF, thereby making absolute rCBF quantification challenging.

PURPOSE

To obtain rCBF in a noninvasive manner using a novel diffusion imaging method with phase contrast (DPC), in which the total CBF from phase-contrast (PC) MRI was utilized to convert perfusion-related diffusion coefficients to rCBF values.

STUDY TYPE

Prospective.

SUBJECTS

Eleven healthy volunteers (nine men and two women; mean age, 23.9 years) participated in this study.

FIELD STRENGTH/SEQUENCE

A 3.0 T, single-shot diffusion echo-planar imaging with multiple b-values (0-3000 s/mm² ), PC-MRI, pulsed continuous ASL, and 3D T₁ -weighted fast field echo.

ASSESSMENT

rCBF and its correlations in the gray matter (GM) and white matter (WM) were compared between DPC and ASL methods. rCBF in the GM and WM and the GM/WM ratio were compared with the literature values obtained using [¹⁵ O]-water positron emission tomography (¹⁵ O-H₂ O PET).

STATISTICAL TESTS

Spearman's correlation coefficient and Wilcoxon signed-rank test were used. Significance was set at P < 0.05.

RESULTS

A significant positive correlation between DPC and ASL in terms of rCBF was observed in GM (R = 0.9), whereas the correlation between the two methods was poor in WM (R = 0.09). The rCBF in GM and WM and the GM/WM ratio obtained using DPC were consistent with the literature values assessed using ¹⁵ O-H₂ O PET. The rCBF value obtained using DPC was significantly higher in the GM and WM than that using ASL.

DATA CONCLUSION

DPC enabled noninvasive quantification of rCBF.

EVIDENCE LEVEL

2 TECHNICAL EFFICACY: Stage 1.

Pattern recognition analysis of directional intravoxel incoherent motion MRI in ischemic rodent brains

This study aimed to demonstrate a reliable automatic segmentation method for independently separating reduced diffusion and decreased perfusion areas in ischemic stroke brains using constrained nonnegative matrix factorization (cNMF) pattern recognition in directional intravoxel incoherent motion MRI (IVIM-MRI). First, the feasibility of cNMF-based segmentation of IVIM signals was investigated in both simulations and in vivo experiments. The cNMF analysis was independently performed for S₀ -normalized and scaled (by the difference between the maximum and minimum) IVIM signals, respectively. Segmentations of reduced diffusion (from S₀ -normalized IVIM signals) and decreased perfusion (from scaled IVIM signals) areas were performed using the corresponding cNMF pattern weight maps. Second, Monte Carlo simulations were performed for directional IVIM signals to investigate the relationship between the degree of vessel alignment and the direction of the diffusion gradient. Third, directional IVIM-MRI experiments (x, y and z diffusion-gradient directions, 20 b values at 7 T) were performed for normal (n = 4), sacrificed (n = 1, no flow) and ischemic stroke models (n = 4, locally reduced flow). The results showed that automatic segmentation of the hypoperfused lesion using cNMF analysis was more accurate than segmentation using the conventional double-exponential fitting. Consistent with the simulation, the double-exponential pattern of the IVIM signals was particularly strong in white matter and ventricle regions when the directional flows were aligned with the applied diffusion-gradient directions. cNMF analysis of directional IVIM signals allowed robust automated segmentation of white matter, ventricle, vascular and hypoperfused regions in the ischemic brain. In conclusion, directional IVIM signals were simultaneously sensitive to diffusion and aligned flow and were particularly useful for automatically segmenting ischemic lesions via cNMF-based pattern recognition.

Acute renal impairment characterization using diffusion magnetic resonance imaging: Validation by histology

Diffusion magnetic resonance imaging has been demonstrated to be a simple, noninvasive and accurate method for the detection of renal microstructure and microcirculation, which are closely linked to renal function. Moreover, serum endothelin-1 (ET-1) was also reported as a good indicator of early renal injury. The aim of this study was to evaluate the feasibility and capability of diffusion MRI and ET-1 to detect acute kidney injury by an operation simulating high-pressure renal pelvic perfusion, which is commonly used during ureteroscopic lithotripsy. Histological findings were used as a reference. Fourteen New Zealand rabbits in an experimental group and 14 in a control group were used in this study. Diffusion tensor imaging and intravoxel incoherent motion diffusion-weighted imaging were acquired by a 3.0 T MRI scanner. Significant corticomedullary differences were found in the values of the apparent diffusion coefficient (ADC), pure tissue diffusion, volume fraction of pseudo-diffusion (fp) and fractional anisotropy (FA) (P < 0.05 for all) in both preoperation and postoperation experimental groups. Compared with the control group, the values of cortical fp_mean , medullary ADC_mean and FA_mean decreased significantly (P < 0.05) after the operation in the experimental group. Also, the change rate of medullary ADC_mean in the experimental group was more pronounced than that in the control group (P = 0.018). No significant change was found in serum ET-1 concentration after surgery in either the experimental (P = 0.80) or control (P = 0.17) groups. In the experimental group, histological changes were observed in the medulla, while no visible change was found in the cortex. This study demonstrated the feasibility of diffusion MRI to detect the changes of renal microstructure and microcirculation in acute kidney injury, with the potential to evaluate renal function. Moreover, the sensitivity of diffusion MRI to acute kidney injury appears to be superior to that of serum ET-1.

Estimation of diffusion, perfusion and fractional volumes using a multi-compartment relaxation-compensated intravoxel incoherent motion (IVIM) signal model

Compartmental diffusion MRI models that account for intravoxel incoherent motion (IVIM) of blood perfusion allow for estimation of the fractional volume of the microvascular compartment. Conventional IVIM models are known to be biased by not accounting for partial volume effects caused by free water and cerebrospinal fluid (CSF), or for tissue-dependent relaxation effects. In this work, a three-compartment model (tissue, free water and blood) that includes relaxation terms is introduced. To estimate the model parameters, in vivo human data were collected with multiple echo times (TE), inversion times (TI) and b-values, which allowed a direct relaxation estimate alongside estimation of perfusion, diffusion and fractional volume parameters. Compared to conventional two-compartment models (with and without relaxation compensation), the three-compartment model showed less effects of CSF contamination. The proposed model yielded significantly different volume fractions of blood and tissue compared to the non-relaxation-compensated model, as well as to the conventional two-compartment model, suggesting that previously reported parameter ranges, using models that do not account for relaxation, should be reconsidered.

On the Field Strength Dependence of Bi- and Triexponential Intravoxel Incoherent Motion (IVIM) Parameters in the Liver

BACKGROUND

Studies on intravoxel incoherent motion (IVIM) imaging are carried out with different acquisition protocols.

PURPOSE

To investigate the dependence of IVIM parameters on the B₀ field strength when using a bi- or triexponential model.

STUDY TYPE

Prospective.

STUDY POPULATION

20 healthy volunteers (age: 19-28 years).

FIELD STRENGTH/SEQUENCE

Volunteers were examined at two field strengths (1.5 and 3T). Diffusion-weighted images of the abdomen were acquired at 24 b-values ranging from 0.2 to 500 s/mm² .

ASSESSMENT

ROIs were manually drawn in the liver. Data were fitted with a bi- and a triexponential IVIM model. The resulting parameters were compared between both field strengths.

STATISTICAL TESTS

One-way analysis of variance (ANOVA) and Kruskal-Wallis test were used to test the obtained IVIM parameters for a significant field strength dependency.

RESULTS

At b-values below 6 s/mm² , the triexponential model provided better agreement with the data than the biexponential model. The average tissue diffusivity was D = 1.22/1.00 μm² /msec at 1.5/3T. The average pseudodiffusion coefficients for the biexponential model were D^* = 308/260 μm² /msec at 1.5/3T; and for the triexponential model D 1 * = 81.3/65.9 μm² /msec, D 2 * = 2453/2333 μm² /msec at 1.5/3T. The average perfusion fractions for the biexponential model were f = 0.286/0.303 at 1.5/3T; and for the triexponential model f₁ = 0.161/0.174 and f₂ = 0.152/0.159 at 1.5/3T. A significant B₀ dependence was only found for the biexponential pseudodiffusion coefficient (ANOVA/KW P = 0.037/0.0453) and tissue diffusivity (ANOVA/KW: P < 0.001).

DATA CONCLUSION

Our experimental results suggest that triexponential pseudodiffusion coefficients and perfusion fractions obtained at different field strengths could be compared across different studies using different B₀ . However, it is recommended to take the field strength into account when comparing tissue diffusivities or using the biexponential IVIM model. Considering published values for oxygenation-dependent transversal relaxation times of blood, it is unlikely that the two blood compartments of the triexponential model represent venous and arterial blood.

LEVEL OF EVIDENCE

1 Technical Efficacy Stage: 2 J. Magn. Reson. Imaging 2019;50:1883-1892.

[Preparation of a Small Acute-phase Cerebral Infarction Phantom for Diffusion-weighted Imaging]

The present study aimed to prepare a small acute-phase cerebral infarction phantom made of gelatin and sucrose to simulate brain parenchymal cells, and a phantom made of collagen peptides and sucrose to simulate cerebral infarction for diffusion-weighted imaging (DWI). During the preparation of gelatin and sucrose mixture (17.0 wt% gelatin, 20.0 wt% sucrose), a cylindrical wooden bar was placed in the center of the phantom and covered with a heat-shrinkable film to ensure space remained after gelling. A mixed solution composed of collagen peptide and sucrose (16.0 wt% collagen peptide, 27.5 wt% sucrose) was then enclosed within the space. The T₂ relaxation time and apparent diffusion coefficient (ADC) of the phantom were set equal to those observed in actual patients with acute-phase cerebral infarction. The mixture was selected based on the signal intensity of both the healthy brain tissue and that subjected to acute cerebral infarction, such that no contrast was observed during T₂-weighted imaging (T₂WI). T₂WI and DWI were performed using a 1.5 T scanner. Although contrast between the mixed gel and mixed solution was obscure on T₂WI, cerebral infarction was clearly visible on DWI. However, the phantom exhibited mono-exponential changes in the ADC value at b values of 0 and 1,000 (s/mm²), and was affected by the proton density and T₁ value depending on the imaging condition.

Hemodynamically self-corrected ΔADC analysis in idiopathic normal pressure hydrocephalus

OBJECTIVE

To clarify the cause of higher water fluctuation of the brain in idiopathic normal pressure hydrocephalus (iNPH), we assessed change in hemodynamic-independent apparent diffusion coefficient during the cardiac cycle (ΔADC) in iNPH.

METHODS

Electrocardiographically synchronized single-shot diffusion echo-planer imaging (b = 0, 500, and 1000 s/mm²) was performed in healthy volunteers, atrophic ventricular dilation group, and iNPH group, respectively. The ΔADC (b = 0 and 1000 s/mm²) and maximum ADC (b = 0 and 500 s/mm²) in the cardiac cycles were measured at the frontal white matter in the brain. Then, self-corrected ΔADC was obtained from the ΔADC divided by the maximum ADC (ADC_peak: perfusion-related diffusion) to correct the blood flow effect.

RESULTS

The ΔADC after correction was significantly higher in the iNPH group than in the other two groups. However, there was no significant difference in ADC_peak values among the groups.

CONCLUSION

Self-corrected ΔADC in iNPH increased because of changes in the biomechanical properties of the brain. Self-corrected ΔADC analysis makes it possible to obtain information on hemodynamically independent water fluctuation as well as perfusion in iNPH.

ADVANCES IN KNOWLEDGE

Analysis self-corrected ΔADC provides simultaneously information on biomechanical properties, perfusion, and water fluctuation in iNPH.

A robust deconvolution method to disentangle multiple water pools in diffusion MRI

The diffusion-weighted magnetic resonance imaging (dMRI) signal measured in vivo arises from multiple diffusion domains, including hindered and restricted water pools, free water and blood pseudo-diffusion. Not accounting for the correct number of components can bias metrics obtained from model fitting because of partial volume effects that are present in, for instance, diffusion tensor imaging (DTI) and diffusion kurtosis imaging (DKI). Approaches that aim to overcome this shortcoming generally make assumptions about the number of considered components, which are not likely to hold for all voxels. The spectral analysis of the dMRI signal has been proposed to relax assumptions on the number of components. However, it currently requires a clinically challenging signal-to-noise ratio (SNR) and accounts only for two diffusion processes defined by hard thresholds. In this work, we developed a method to automatically identify the number of components in the spectral analysis, and enforced its robustness to noise, including outlier rejection and a data-driven regularization term. Furthermore, we showed how this method can be used to take into account partial volume effects in DTI and DKI fitting. The proof of concept and performance of the method were evaluated through numerical simulations and in vivo MRI data acquired at 3 T. With simulations our method reliably decomposed three diffusion components from SNR = 30. Biases in metrics derived from DTI and DKI were considerably reduced when components beyond hindered diffusion were taken into account. With the in vivo data our method determined three macro-compartments, which were consistent with hindered diffusion, free water and pseudo-diffusion. Taking free water and pseudo-diffusion into account in DKI resulted in lower mean diffusivity and higher fractional anisotropy values in both gray and white matter. In conclusion, the proposed method allows one to determine co-existing diffusion compartments without prior assumptions on their number, and to account for undesired signal contaminations within clinically achievable SNR levels.

High-Resolution MR Imaging of Muscular Fat Fraction-Comparison of Three T2-Based Methods and Chemical Shift-Encoded Imaging

Chemical shift-encoded imaging (CSEI) is the most common magnetic resonance imaging fat–water separation method. However, when high spatial resolution fat fraction (FF) images are desired, CSEI might be challenging owing to the increased interecho spacing. Here, 3 T₂-based methods have been assessed as alternative methods for obtaining high-resolution FF images. Images from the calf of 10 healthy volunteers were acquired; FF maps were then estimated using 3 T₂-based methods (2- and 3-parameter nonlinear least squares fit and a Bayesian probability method) and CSEI for reference. In addition, simulations were conducted to characterize the performance of various methods. Here, all T₂-based methods resulted in qualitatively improved high-resolution FF images compared with high-resolution CSEI. The 2-parameter fit showed best quantitative agreement to low-resolution CSEI, even at low FF. The estimated T₂-values of fat and water, and the estimated muscle FF of the calf, agreed well with previously published data. In conclusion, T₂-based methods can provide improved high-resolution FF images of the calf compared with the CSEI method.

A Modified Tri-Exponential Model for Multi-b-value Diffusion-Weighted Imaging: A Method to Detect the Strictly Diffusion-Limited Compartment in Brain

Purpose: To present a new modified tri-exponential model for diffusion-weighted imaging (DWI) to detect the strictly diffusion-limited compartment, and to compare it with the conventional bi- and tri-exponential models.

Methods: Multi-b-value diffusion-weighted imaging (DWI) with 17 b-values up to 8,000 s/mm² were performed on six volunteers. The corrected Akaike information criterions (AICc) and squared predicted errors (SPE) were calculated to compare these three models.

Results: The mean f₀ values were ranging 11.9–18.7% in white matter ROIs and 1.2–2.7% in gray matter ROIs. In all white matter ROIs: the AICcs of the modified tri-exponential model were the lowest (p < 0.05 for five ROIs), indicating the new model has the best fit among these models; the SPEs of the bi-exponential model were the highest (p < 0.05), suggesting the bi-exponential model is unable to predict the signal intensity at ultra-high b-value. The mean ADC_very−slow values were extremely low in white matter (1–7 × 10⁻⁶ mm²/s), but not in gray matter (251–445 × 10⁻⁶ mm²/s), indicating that the conventional tri-exponential model fails to represent a special compartment.

Conclusions: The strictly diffusion-limited compartment may be an important component in white matter. The new model fits better than the other two models, and may provide additional information.

Differentiation Between Luminal-A and Luminal-B Breast Cancer Using Intravoxel Incoherent Motion and Dynamic Contrast-Enhanced Magnetic Resonance Imaging

RATIONALE AND OBJECTIVES

The study aimed to investigate whether intravoxel incoherent motion (IVIM) and dynamic contrast-enhanced (DCE) magnetic resonance imaging (MRI) can differentiate luminal-B from luminal-A breast cancer MATERIALS AND METHODS: Biexponential analyses of IVIM and DCE MRI were performed using a 3.0-T MRI scanner, involving 134 patients with 137 pathologically confirmed luminal-type invasive breast cancers. Luminal-type breast cancer was categorized as luminal-B breast cancer (LBBC, Ki-67 ≧ 14%) or luminal-A breast cancer (LABC, Ki-67 < 14%). Quantitative parameters from IVIM (pure diffusion coefficient [D], perfusion-related diffusion coefficient [D*], and fraction [f]) and DCE MRI (initial percentage of enhancement and signal enhancement ratio [SER]) were calculated. The apparent diffusion coefficient (ADC) was also calculated using monoexponential fitting. We correlated these data with the Ki-67 status.

RESULTS

The D and ADC values of LBBC were significantly lower than those of LABC (P = 0.028, P = 0.037). The SER of LBBC was significantly higher than that of LABC (P = 0.004). A univariate analysis showed that a significantly lower D (<0.847 x 10^-3 mm²/s), lower ADC (<0.960 × 10^-3 mm²/s), and higher SER (>1.071) values were associated with LBBC (all P values <0.01), compared to LABC. In a multivariate analysis, a higher SER (>1.071; odds ratio: 3.0099, 95% confidence interval: 1.4246-6.3593; P = 0.003) value and a lower D (<0.847 × 10^-3 mm²/s; odds ratio: 2.6878, 95% confidence interval: 1.0445-6.9162; P = 0.040) value were significantly associated with LBBC, compared to LABC.

CONCLUSION

The SER derived from DCE MRI and the D derived from IVIM are associated independently with the Ki-67 status in patients with luminal-type breast cancer.

Rapid measurement of intravoxel incoherent motion (IVIM) derived perfusion fraction for clinical magnetic resonance imaging

OBJECTIVE

This study aimed to investigate the reliability of intravoxel incoherent motion (IVIM) model derived parameters D and f and their dependence on b value distributions with a rapid three b value acquisition protocol.

MATERIALS AND METHODS

Diffusion models for brain, kidney, and liver were assessed for bias, error, and reproducibility for the estimated IVIM parameters using b values 0 and 1000, and a b value between 200 and 900, at signal-to-noise ratios (SNR) 40, 55, and 80. Relative errors were used to estimate optimal b value distributions for each tissue scenario. Sixteen volunteers underwent brain DW-MRI, for which bias and coefficient of variation were determined in the grey matter.

RESULTS

Bias had a large influence in the estimation of D and f for the low-perfused brain model, particularly at lower b values, with the same trends being confirmed by in vivo imaging. Significant differences were demonstrated in vivo for estimation of D (P = 0.029) and f (P < 0.001) with [300,1000] and [500,1000] distributions. The effect of bias was considerably lower for the high-perfused models. The optimal b value distributions were estimated to be brain_500,1000, kidney_300,1000, and liver_200,1000.

CONCLUSION

IVIM parameters can be estimated using a rapid DW-MRI protocol, where the optimal b value distribution depends on tissue characteristics and compromise between bias and variability.

Separating blood and water: Perfusion and free water elimination from diffusion MRI in the human brain

The assessment of the free water fraction in the brain provides important information about extracellular processes such as atrophy and neuroinflammation in various clinical conditions as well as in normal development and aging. Free water estimates from diffusion MRI are assumed to account for freely diffusing water molecules in the extracellular space, but may be biased by other pools of molecules in rapid random motion, such as the intravoxel incoherent motion (IVIM) of blood, where water molecules perfuse in the randomly oriented capillary network. The goal of this work was to separate the signal contribution of the perfusing blood from that of free-water and of other brain diffusivities. The influence of the vascular compartment on the estimation of the free water fraction and other diffusivities was investigated by simulating perfusion in diffusion MRI data. The perfusion effect in the simulations was significant, especially for the estimation of the free water fraction, and was maintained as long as low b-value data were included in the analysis. Two approaches to reduce the perfusion effect were explored in this study: (i) increasing the minimal b-value used in the fitting, and (ii) using a three-compartment model that explicitly accounts for water molecules in the capillary blood. Estimation of the model parameters while excluding low b-values reduced the perfusion effect but was highly sensitive to noise. The three-compartment model fit was more stable and additionally, provided an estimation of the volume fraction of the capillary blood compartment. The three-compartment model thus disentangles the effects of free water diffusion and perfusion, which is of major clinical importance since changes in these components in the brain may indicate different pathologies, i.e., those originating from the extracellular space, such as neuroinflammation and atrophy, and those related to the vascular space, such as vasodilation, vasoconstriction and capillary density. Diffusion MRI data acquired from a healthy volunteer, using multiple b-shells, demonstrated an expected non-zero contribution from the blood fraction, and indicated that not accounting for the perfusion effect may explain the overestimation of the free water fraction evinced in previous studies. Finally, the applicability of the method was demonstrated with a dataset acquired using a clinically feasible protocol with shorter acquisition time and fewer b-shells.

Technical Note: Development of a cranial phantom for assessing perfusion, diffusion, and biomechanics

PURPOSE

A novel cranial phantom was developed to simulate the relationships among factors such as blood perfusion, water diffusion, and biomechanics in intracranial tissue.

METHODS

The cranial phantom consisted of a high-density polypropylene filter (mimicking brain parenchyma) with intra- and extrafilter spaces (mimicking cerebral artery and vein, respectively), and a capacitor space (mimicking the cerebrospinal fluid space). Pulsatile and steady flow with different flow rates were applied to the cranial phantom using a programmable pump. On 3.0-T MRI, the measurements of the internal pressure in the phantom, apparent diffusion coefficient (ADC) with monoexponential analysis in the filter, and total simulated cerebral blood flow (tSCBF) into the phantom were synchronized with the pulsatile flow. We obtained their maximum changes during the pulsation period (ΔP, ΔADC, and ΔtSCBF, respectively). Then, the compliance index (CI) was calculated by dividing the volume change (ΔV) by the ΔP in the phantom. Moreover, the same measurements were repeated after the compliance of the phantom was reduced by increasing the water volume in the capacitor space. Under steady flow conditions, we determined the regional SCBF (rSCBF) and perfusion-related and restricted diffusion coefficients (D* and D, respectively) with biexponential analysis in the filter.

RESULTS

The internal pressure, ADC, and tSCBF varied over the pulsation period depending on the input flow. Moreover, the ΔP, ΔADC, ΔtSCBF, and rSCBF increased with the input flow rate. Compared to the high compliance condition, in the low compliance condition, the ΔP and ΔADC were higher by factors of 2.5 and 1.3, respectively, and the CI was smaller by a factor of 2.7, whereas the ΔV was almost unchanged. The D* was strongly affected by the input flow.

CONCLUSION

Our original phantom models the relationships among the blood perfusion, water diffusion, and biomechanics of the intracranial tissue, potentially facilitating the validation of novel MRI techniques and optimization of imaging parameters.

Simple noise reduction for diffusion weighted images

Our purpose in this study was to reduce the noise in order to improve the SNR of Dw images with high b-value by using two correction schemes. This study was performed with use of phantoms made from water and sucrose at different concentrations, which were 10, 30, and 50 weight percent (wt%). In noise reduction for Dw imaging of the phantoms, we compared two correction schemes that are based on the Rician distribution and the Gaussian distribution. The highest error values for each concentration with use of the Rician distribution scheme were 7.3 % for 10 wt%, 2.4 % for 30 wt%, and 0.1 % for 50 wt%. The highest error values for each concentration with use of the Gaussian distribution scheme were 20.3 % for 10 wt%, 11.6 % for 30 wt%, and 3.4 % for 50 wt%. In Dw imaging, the noise reduction makes it possible to apply the correction scheme of Rician distribution.