Biexponential Intravoxel Incoherent Motion Research Articles

Modeling cerebrospinal fluid dynamics across the entire intracranial space through integration of four-dimensional flow and intravoxel incoherent motion magnetic resonance imaging

BackgroundBidirectional reciprocal motion of cerebrospinal fluid (CSF) was quantified using four-dimensional (4D) flow magnetic resonance imaging (MRI) and intravoxel incoherent motion (IVIM) MRI. To estimate various CSF motions in the entire intracranial region, we attempted to integrate the flow parameters calculated using the two MRI sequences. To elucidate how CSF dynamics deteriorate in Hakim’s disease, an age-dependent chronic hydrocephalus, flow parameters were estimated from the two MRI sequences to assess CSF motion in the entire intracranial region.MethodsThis study included 127 healthy volunteers aged ≥ 20 years and 44 patients with Hakim’s disease. On 4D flow MRI for measuring CSF motion, velocity encoding was set at 5 cm/s. For the IVIM MRI analysis, the diffusion-weighted sequence was set at six b-values (i.e., 0, 50, 100, 250, 500, and 1000 s/mm2), and the biexponential IVIM fitting method was adapted. The relationships between the fraction of incoherent perfusion (f) on IVIM MRI and 4D flow MRI parameters including velocity amplitude (VA), absolute maximum velocity, stroke volume, net flow volume, and reverse flow rate were comprehensively evaluated in seven locations in the ventricles and subarachnoid spaces. Furthermore, we developed a new parameter for fluid oscillation, the Fluid Oscillation Index (FOI), by integrating these two measurements. In addition, we investigated the relationship between the measurements and indices specific to Hakim’s disease and the FOIs in the entire intracranial space.ResultsThe VA on 4D flow MRI was significantly associated with the mean f-values on IVIM MRI. Therefore, we estimated VA that could not be directly measured on 4D flow MRI from the mean f-values on IVIM MRI in the intracranial CSF space, using the following formula; e0.2(f−85) + 0.25. To quantify fluid oscillation using one integrated parameter with weighting, FOI was calculated as VA × 10 + f × 0.02. In addition, the FOIs at the left foramen of Luschka had the strongest correlations with the Evans index (Pearson’s correlation coefficient: 0.78). The other indices related with Hakim’s disease were significantly associated with the FOIs at the cerebral aqueduct and bilateral foramina of Luschka. FOI at the cerebral aqueduct was also elevated in healthy controls aged ≥ 60 years.ConclusionsWe estimated pulsatile CSF movements in the entire intracranial CSF space in healthy individuals and patients with Hakim’s disease using FOI integrating VA from 4D flow MRI and f-values from IVIM MRI. FOI is useful for quantifying the CSF oscillation.

Quantitative characterization of breast lesions and normal fibroglandular tissue using compartmentalized diffusion-weighted model: comparison of intravoxel incoherent motion and restriction spectrum imaging

BackgroundTo compare the compartmentalized diffusion-weighted models, intravoxel incoherent motion (IVIM) and restriction spectrum imaging (RSI), in characterizing breast lesions and normal fibroglandular tissue.MethodsThis prospective study enrolled 152 patients with 157 histopathologically verified breast lesions (41 benign and 116 malignant). All patients underwent a full-protocol preoperative breast MRI, including a multi-b-value DWI sequence. The diffusion parameters derived from the mono-exponential model (ADC), IVIM model (Dt, Dp, f), and RSI model (C1, C2, C3, C1C2, F1, F2, F3, F1F2) were quantitatively measured and then compared among malignant lesions, benign lesions and normal fibroglandular tissues using Kruskal-Wallis test. The Mann-Whitney U-test was used for the pairwise comparisons. Diagnostic models were built by logistic regression analysis. The ROC analysis was performed using five-fold cross-validation and the mean AUC values were calculated and compared to evaluate the discriminative ability of each parameter or model.ResultsAlmost all quantitative diffusion parameters showed significant differences in distinguishing malignant breast lesions from both benign lesions (other than C2) and normal fibroglandular tissue (all parameters) (all P < 0.0167). In terms of the comparisons of benign lesions and normal fibroglandular tissues, the parameters derived from IVIM (Dp, f) and RSI (C1, C2, C1C2, F1, F2, F3) showed significant differences (all P < 0.005). When using individual parameters, RSI-derived parameters-F1, C1C2, and C2 values yielded the highest AUCs for the comparisons of malignant vs. benign, malignant vs. normal tissue and benign vs. normal tissue (AUCs = 0.871, 0.982, and 0.863, respectively). Furthermore, the combined diagnostic model (IVIM + RSI) exhibited the highest diagnostic efficacy for the pairwise discriminations (AUCs = 0.893, 0.991, and 0.928, respectively).ConclusionsQuantitative parameters derived from the three-compartment RSI model have great promise as imaging indicators for the differential diagnosis of breast lesions compared with the bi-exponential IVIM model. Additionally, the combined model of IVIM and RSI achieves superior diagnostic performance in characterizing breast lesions.

Assessment of fitting methods and variability of IVIM parameters in muscles of the lumbar spine at rest.

Intravoxel incoherent motion (IVIM) MRI provides insight into tissue diffusion and perfusion. Here, estimates of perfusion fraction ( ), pseudo-diffusion coefficient ( ), and diffusion coefficient ( ) obtained via different fitting methods are compared to ascertain (1) the optimal analysis strategy for muscles of the lumbar spine and (2) repeatability of IVIM parameters in skeletal muscle at rest. Diffusion-weighted images were acquired in the lumbar spine at rest in 15 healthy participants. Data were fit to the bi-exponential IVIM model to estimate and using three variably segmented approaches based on non-linear least squares fitting, and a Bayesian fitting method. Assuming that perfusion and diffusion are temporally stable in skeletal muscle at rest, and spatially uniform within a spinal segment, the optimal analysis strategy was determined as the approach with the lowest temporal or spatial variation and smallest residual between measured and fit data. Inter-session repeatability of IVIM parameters was evaluated in a subset of 11 people. Finally, simulated IVIM signal at varying signal to noise ratio were evaluated to understand precision and bias. Experimental results showed that IVIM parameter values differed depending on the fitting method. A three-step non-linear least squares fitting approach, where , and were estimated sequentially, generally yielded the lowest spatial and temporal variation. Solving all parameters simultaneously yielded the lowest residual between measured and fit data, however there was substantial spatial and temporal variability. Results obtained by Bayesian fitting had high spatial and temporal variability in addition to a large residual between measured and fit data. Simulations showed that all fitting methods can fit the IVIM data at signal to noise ratios >35, and that was the most challenging to accurately obtain. Overall, this study motivates use of a three-step non-linear least squares fitting strategy to quantify IVIM parameters in skeletal muscle.

Reproducibility and Repeatability of Intravoxel Incoherent Motion MRI Acquisition Methods in Liver.

Intravoxel incoherent motion (IVIM) diffusion weighted MRI (DWI) has potential for evaluating hepatic fibrosis but image acquisition technique influence on diffusion parameter estimation bears investigation. To minimize variability and maximize repeatably in abdominal DWI in terms of IVIM parameter estimates. Prospective test-retest and image quality comparison. Healthy volunteers (3F/7M, 29.9 ± 12.9 years) and Family Study subjects (18F/12M, 51.7 ± 16.7 years), without and with liver steatosis. Abdominal single-shot echo-planar imaging (EPI) and simultaneous multi-slice (SMS) DWI sequences with respiratory triggering (RT), breath-holding (BH), and navigator echo (NE) at 3 Tesla. SMS-BH, EPI-NE, and SMS-RT data from twice-scanned healthy volunteers were analyzed using 6 × b-values (0-800 s⋅mm-2 ) and lower (LO) and higher (HI) b-value ranges. Family Study subjects were scanned using SMS and standard EPI sequences. The biexponential IVIM model was used to estimate fast-diffusion coefficient (Df ), fraction of fast diffusion (f), and slow-diffusion coefficient (Ds ). Scan time, estimated signal-to-noise ratio (eSNR), eSNR per acquisition, and distortion ratio were compared. Coefficients of variation (CoV) and Bland Altman analyses were performed for test-retest repeatability. Interclass correlation coefficient (ICC) assessed interobserver agreement with P < 0.05 deemed significant. Within-subject CoVs among volunteers (N = 10) for f and Ds were lowest in EPI-NE-LO (11.6%) and SMS-RT-HI (11.1%). Inter-observer ICCs for f and Ds were highest for EPI-NE-LO (0.63) and SMS-RT-LO (0.76). Df could not be estimated for most subjects. Estimated eSNR (EPI = 21.9, SMS = 4.7) and eSNR time (EPI = 6.7, SMS = 16.6) were greater for SMS, with less distortion in the liver region (DR-PE: EPI = 23.6, SMS = 13.1). Simultaneous multislice acquisitions had significantly less variability and higher ICCs of Ds , higher eSNR, less distortion, and reduced scan time compared to EPI. 2 TECHNICAL EFFICACY: Stage 1.

Effect of simultaneous multislice imaging, slice properties, and repetition time on the measured magnetic resonance biexponential intravoxel incoherent motion in the liver.

This study aims to investigate the previously reported dependency of intravoxel incoherent motion (IVIM) parameters on simultaneous multislice (SMS) acquisition and repetition time (TR). This includes the influence of slice thickness, slice gaps, and slice order on measured IVIM parameters. Diffusion-weighted imaging (DWI) of the liver was performed on 10 healthy volunteers (aged 20-30 years) at 3T with a slice thickness of 5 mm, a slice gap of 5 mm, and a linear slice order. Diffusion-weighted images were acquired with 19 b-values (0-800 s/mm2) using both conventional slice excitation with an acceleration factor of one (AF1) and SMS excitation with an acceleration factor of three (AF3). Each of these measurements were carried out with two repetition times (TRs)- 1,300 ms (prefix s) and 4,500 ms (prefix l)-resulting in four different combinations: sAF1, sAF3, lAF1, and lAF3. Five volunteers underwent additional measurements using a 10 mm slice thickness and with AF1. Median signal values in the liver were used to determine the biexponential IVIM parameters. Statistical significances were assessed using the Kruskal-Wallis test, Wilcoxon signed-rank test, and Student's t-test. In-silico investigations were also used to interpret the data. There were no significant differences between the biexponential IVIM parameters acquired from sAF1, sAF3, lAF1, and lAF3. Median values of the perfusion fraction f were as follows: 29.9% (sAF1), 26.9% (sAF3), 28.1% (lAF1), and 27.5% (lAF3). In the 10 mm-thick slices, f decreased from 31.3% (lAF1) to 27.4% (sAF1) (p = 0.141). The slice excitation mode did not appear to have any significant influence on the biexponential IVIM parameters. However, our simulations, as well as values reported from previous publications, show that slice thickness, slice gaps, and slice order are relevant and should thus be reported in IVIM studies.

Multisite MRI Intravoxel Incoherent Motion Repeatability and Reproducibility across 3 T Scanners in a Breast Diffusion Phantom: A BReast Intravoxel Incoherent Motion Multisite (BRIMM) Study.

Monoexponential apparent diffusion coefficient (ADC) and biexponential intravoxel incoherent motion (IVIM) analysis of diffusion-weighted imaging is helpful in the characterization of breast tumors. However, repeatability/reproducibility studies across scanners and across sites are scarce. To evaluate the repeatability and reproducibility of ADC and IVIM parameters (tissue diffusivity (Dt), perfusion fraction (Fp) and pseudo-diffusion (Dp)) within and across sites employing MRI scanners from different vendors utilizing 16-channel breast array coils in a breast diffusion phantom. Phantom repeatability. A breast phantom containing tubes of different polyvinylpyrrolidone (PVP) concentrations, water, fat, and sponge flow chambers, together with an MR-compatible liquid crystal (LC) thermometer. Bipolar gradient twice-refocused spin echo sequence and monopolar gradient single spin echo sequence at 3 T. Studies were performed twice in each of two scanners, located at different sites, on each of 2 days, resulting in four studies per scanner. ADCs of the PVP and water were normalized to the vendor-provided calibrated values at the temperature indicated by the LC thermometer for repeatability/reproducibility comparisons. ADC and IVIM repeatability and reproducibility within and across sites were estimated via the within-system coefficient of variation (wCV). Pearson correlation coefficient (r) was also computed between IVIM metrics and flow speed. A P value <0.05 was considered statistically significant. ADC and Dt demonstrated excellent repeatability (<2%; <3%, respectively) and reproducibility (both <5%) at the two sites. Fp and Dp exhibited good repeatability (mean of two sites 3.67% and 5.59%, respectively) and moderate reproducibility (mean of two sites 15.96% and 13.3%, respectively). The mean intersite reproducibility (%) of Fp/Dp/Dt was 50.96/13.68/5.59, respectively. Fp and Dt demonstrated high correlations with flow speed while Dp showed lower correlations. Fp correlations with flow speed were significant at both sites. IVIM reproducibility results were promising and similar to ADC, particularly for Dt. The results were reproducible within both sites, and a progressive trend toward reproducibility across sites except for Fp. 2 TECHNICAL EFFICACY: Stage 1.

Usefulness of intravoxel incoherent motion MRI for visualizing slow cerebrospinal fluid motion

BackgroundIn the cerebrospinal fluid (CSF) dynamics, the pulsations of cerebral arteries and brain is considered the main driving force for the reciprocating bidirectional CSF movements. However, measuring these complex CSF movements on conventional flow-related MRI methods is difficult. We tried to visualize and quantify the CSF motion by using intravoxel incoherent motion (IVIM) MRI with low multi-b diffusion-weighted imaging.MethodsDiffusion-weighted sequence with six b values (0, 50, 100, 250, 500, and 1000 s/mm2) was performed on 132 healthy volunteers aged ≥ 20 years and 36 patients with idiopathic normal pressure hydrocephalus (iNPH). The healthy volunteers were divided into three age groups (< 40, 40 to < 60, and ≥ 60 years). In the IVIM analysis, the bi-exponential IVIM fitting method using the Levenberg–Marquardt algorithm was adapted. The average, maximum, and minimum values of ADC, D, D*, and fraction of incoherent perfusion (f) calculated by IVIM were quantitatively measured in 45 regions of interests in the whole ventricles and subarachnoid spaces.ResultsCompared with healthy controls aged ≥ 60 years, the iNPH group had significantly lower mean f values in all the parts of the lateral and 3rd ventricles, whereas significantly higher mean f value in the bilateral foramina of Luschka. In the bilateral Sylvian fossa, which contain the middle cerebral bifurcation, the mean f values increased gradually with increasing age, whereas those were significantly lower in the iNPH group. In the 45 regions of interests, the f values in the bilateral foramina of Luschka were the most positively correlated with the ventricular size and indices specific to iNPH, whereas that in the anterior part of the 3rd ventricle was the most negatively correlated with the ventricular size and indices specific to iNPH. Other parameters of ADC, D, and D* were not significantly different between the two groups in any locations.ConclusionsThe f value on IVIM MRI is useful for evaluating small pulsatile complex motion of CSF throughout the intracranial CSF spaces. Patients with iNPH had significantly lower mean f values in the whole lateral ventricles and 3rd ventricles and significantly higher mean f value in the bilateral foramina of Luschka, compared with healthy controls aged ≥ 60 years.

Influence of saturation effects on biexponential liver intravoxel incoherent motion.

Studies on intravoxel incoherent motion (IVIM) imaging in the liver have been carried out with different acquisition protocols. The number of acquired slices and the distances between slices can influence IVIM measurements due to saturation effects, but these effects have often been disregarded. This study investigated differences in biexponential IVIM parameters between two slice settings. Fifteen healthy volunteers (21-30 years) were examined at a field strength of 3T. Diffusion-weighted images of the abdomen were acquired with 16 b values (0-800 s/mm2 ), with four slices for the few slices setting and 24-27 slices for the many slices setting. Regions of interest were manually drawn in the liver. The data were fitted with a monoexponential signal curve and a biexponential IVIM curve, and biexponential IVIM parameters were determined. The dependence on the slice setting was assessed with Student's t test for paired samples (normally distributed IVIM parameters) and the Wilcoxon signed-rank test (non-normally distributed parameters). None of the parameters were significantly different between the settings. For few slices and many slices, respectively, the mean values (SDs) for were ( ) and ( ); for they were 29.7% (6.2%) and 27.7% (3.6%); and for they were ( ) and ( ). Biexponential IVIM parameters in the liver are comparable among IVIM studies that use different slice settings, with mostly negligible saturation effects. However, this may not hold for studies that use much shorter TR.

Noninvasively evaluating the grade and IDH mutation status of gliomas by using mono-exponential, bi-exponential diffusion-weighted imaging and three-dimensional pseudo-continuous arterial spin labeling

ObjectivesTo noninvasively assess the diagnostic performance of diffusion-weighted imaging (DWI), bi-exponential intravoxel incoherent motion imaging (IVIM) and three-dimensional pseudo-continuous arterial spin labeling (3D pCASL) in differentiating lower-grade gliomas (LGGs) from high-grade gliomas (HGGs), and predicting the isocitrate dehydrogenase (IDH) mutation status. Materials and methodsNinety-five patients with pathologically confirmed grade 2–4 gliomas with preoperative DWI, IVIM and 3D pCASL were enrolled in this study. The Student’s t test and Mann-Whitney U test were used to evaluate differences in parameters of DWI, IVIM and 3D pCASL between LGG and HGG as well as between mutant and wild-type IDH in grade 2 and 3 diffusion astrocytoma; receiver operator characteristic (ROC) analysis was used to assess the diagnostic performance. ResultsThe value of ADCmean, ADCmin, Dmean and Dmin in HGGs were lower than in LGGs, while the value of CBFmean and CBFmax in HGGs were higher than in LGGs. In ROC analysis, the AUC values of Dmean, Dmin and CBFmax were 0.827, 0.878 and 0.839, respectively. The combination of CBFmax and Dmin displayed the highest diagnostic performance to distinguish LGGs from HGGs, with AUC 0.906, sensitivity 82.4 %, and specificity 86.4 %. In grades 2 and 3 diffusion astrocytoma patients, ADCmin, Dmean, Dmin, CBFmean and CBFmax showed significant differences between IDHmut and IDHwt group (p < 0.05, 0.001, 0.001, 0.01 and 0.001, respectively) and the AUC values were 0. 709, 0.849, 0.919, 0.755 and 0.873, respectively. Similarly, the combination of CBFmax and Dmin demonstrated the highest AUC value (0.938) in prediction IDH mutation status, with sensitivity 92.9 %, and specificity 95.5 %. ConclusionThe combination of IVIM and 3D pCASL can be used in prediction histologic grade and IDH mutation status of glioma noninvasively.

Self-supervised neural network improves tri-exponential intravoxel incoherent motion model fitting compared to least-squares fitting in non-alcoholic fatty liver disease.

Recent literature suggests that tri-exponential models may provide additional information and fit liver intravoxel incoherent motion (IVIM) data more accurately than conventional bi-exponential models. However, voxel-wise fitting of IVIM results in noisy and unreliable parameter maps. For bi-exponential IVIM, neural networks (NN) were able to produce superior parameter maps than conventional least-squares (LSQ) generated images. Hence, to improve parameter map quality of tri-exponential IVIM, we developed an unsupervised physics-informed deep neural network (IVIM3-NET). We assessed its performance in simulations and in patients with non-alcoholic fatty liver disease (NAFLD) and compared outcomes with bi-exponential LSQ and NN fits and tri-exponential LSQ fits. Scanning was performed using a 3.0T free-breathing multi-slice diffusion-weighted single-shot echo-planar imaging sequence with 18 b-values. Images were analysed for visual quality, comparing the bi- and tri-exponential IVIM models for LSQ fits and NN fits using parameter-map signal-to-noise ratios (SNR) and adjusted R 2. IVIM parameters were compared to histological fibrosis, disease activity and steatosis grades. Parameter map quality improved with bi- and tri-exponential NN approaches, with a significant increase in average parameter-map SNR from 3.38 to 5.59 and 2.45 to 4.01 for bi- and tri-exponential LSQ and NN models respectively. In 33 out of 36 patients, the tri-exponential model exhibited higher adjusted R 2 values than the bi-exponential model. Correlating IVIM data to liver histology showed that the bi- and tri-exponential NN outperformed both LSQ models for the majority of IVIM parameters (10 out of 15 significant correlations). Overall, our results support the use of a tri-exponential IVIM model in NAFLD. We show that the IVIM3-NET can be used to improve image quality compared to a tri-exponential LSQ fit and provides promising correlations with histopathology similar to the bi-exponential neural network fit, while generating potentially complementary additional parameters.

Bifurcated Topological Optimization for IVIM.

In this work, we shed light on the issue of estimating Intravoxel Incoherent Motion (IVIM) for diffusion and perfusion estimation by characterizing the objective function using simplicial homology tools. We provide a robust solution via topological optimization of this model so that the estimates are more reliable and accurate. Estimating the tissue microstructure from diffusion MRI is in itself an ill-posed and a non-linear inverse problem. Using variable projection functional (VarPro) to fit the standard bi-exponential IVIM model we perform the optimization using simplicial homology based global optimization to better understand the topology of objective function surface. We theoretically show how the proposed methodology can recover the model parameters more accurately and consistently by casting it in a reduced subspace given by VarPro. Additionally we demonstrate that the IVIM model parameters cannot be accurately reconstructed using conventional numerical optimization methods due to the presence of infinite solutions in subspaces. The proposed method helps uncover multiple global minima by analyzing the local geometry of the model enabling the generation of reliable estimates of model parameters.

IVIM Parameters on MRI Could Predict ISUP Risk Groups of Prostate Cancers on Radical Prostatectomy.

PurposeTo elucidate the usefulness of intravoxel incoherent motion (IVIM)/apparent diffusion coefficient (ADC) parameters in preoperative risk stratification using International Society of Urological Pathology (ISUP) grades.Materials and MethodsForty-five prostate cancer (PCa) patients undergoing radical prostatectomy (RP) after prostate multiparametric magnetic resonance imaging (mpMRI) were included. The ISUP grades were categorized into low-risk (I-II) and high-risk (III-V) groups, and the concordance between the preoperative and postoperative grades was analyzed. The largest region of interest (ROI) of the dominant tumor on each IVIM/ADC image was delineated to obtain its histogram values (i.e., minimum, mean, and kurtosis) of diffusivity (D), pseudodiffusivity (D*), perfusion fraction (PF), and ADC. Multivariable logistic regression analysis of the IVIM/ADC parameters without and with preoperative ISUP grades were performed to identify predictors for the postoperative high-risk group.ResultsThirty-two (71.1%) of 45 patients had concordant preoperative and postoperative ISUP grades. Dmean, D*kurtosis, PFkurtosis, ADCmin, and ADCmean were significantly associated with the postoperative ISUP risk group (all p < 0.05). Dmean and D*kurtosis (model I, both p < 0.05) could predict the postoperative ISUP high-risk group with an area under the curve (AUC) of 0.842 and a 95% confidence interval (CI) of 0.726–0.958. The addition of D*kurtosis to the preoperative ISUP grade (model II) may enhance prediction performance, with an AUC of 0.907 (95% CI 0.822–0.992).ConclusionsThe postoperative ISUP risk group could be predicted by Dmean and D*kurtosis from mpMRI, especially D*kurtosis. Obtaining the biexponential IVIM parameters is important for better risk stratification for PCa.

Dynamic Contrast Enhanced MRI and Intravoxel Incoherent Motion to Identify Molecular Subtypes of Breast Cancer with Different Vascular Normalization Gene Expression.

ObjectiveTo assess the expression of vascular normalization genes in different molecular subtypes of breast cancer and to determine whether molecular subtypes with a higher vascular normalization gene expression can be identified using dynamic contrast-enhanced (DCE) magnetic resonance imaging (MRI) and intravoxel incoherent motion (IVIM) diffusion-weighted imaging (DWI).Materials and MethodsThis prospective study evaluated 306 female (mean age ± standard deviation, 50 ± 10 years), recruited between January 2014 and August 2017, who had de novo breast cancer larger than 1 cm in diameter (308 tumors). DCE MRI followed by IVIM DWI studies using 11 different b-values (0 to 1200 s/mm2) were performed on a 1.5T MRI system. The Tofts model and segmented biexponential IVIM analysis were used. For each tumor, the molecular subtype (according to six [I-VI] subtypes and PAM50 subtypes), expression profile of genes for vascular normalization, pericytes, and normal vascular signatures were determined using freshly frozen tissue. Statistical associations between imaging parameters and molecular subtypes were examined using logistic regression or linear regression with a significance level of p = 0.05.ResultsBreast cancer subtypes III and VI and PAM50 subtypes luminal A and normal-like exhibited a higher expression of genes for vascular normalization, pericyte markers, and normal vessel function signature (p < 0.001 for all) compared to other subtypes. Subtypes III and VI and PAM50 subtypes luminal A and normal-like, versus the remaining subtypes, showed significant associations with Ktrans, kep, vp, and IAUGCBN90 on DEC MRI, with relatively smaller values in the former. The subtype grouping was significantly associated with D, with relatively less restricted diffusion in subtypes III and VI and PAM50 subtypes luminal A and normal-like.ConclusionDCE MRI and IVIM parameters may identify molecular subtypes of breast cancers with a different vascular normalization gene expression.

B-Value Optimization in the Estimation of Intravoxel Incoherent Motion Parameters in Patients with Cervical Cancer.

ObjectiveThis study aimed to find the optimal number of b-values for intravoxel incoherent motion (IVIM) imaging analysis, using simulated and in vivo data from cervical cancer patients.Materials and MethodsSimulated data were generated using literature pooled means, which served as reference values for simulations. In vivo data from 100 treatment-naïve cervical cancer patients with IVIM imaging (13 b-values, scan time, 436 seconds) were retrospectively reviewed. A stepwise b-value fitting algorithm calculated optimal thresholds. Feed forward selection determined the optimal subsampled b-value distribution for biexponential IVIM fitting, and simplified IVIM modeling using monoexponential fitting was attempted. IVIM parameters computed using all b-values served as reference values for in vivo data.ResultsIn simulations, parameters were accurately estimated with six b-values, or three b-values for simplified IVIM, respectively. In vivo data showed that the optimal threshold was 40 s/mm2 for patients with squamous cell carcinoma and a subsampled acquisition of six b-values (scan time, 198 seconds) estimated parameters were not significantly different from reference parameters (individual parameter error rates of less than 5%). In patients with adenocarcinoma, the optimal threshold was 100 s/mm2, but an optimal subsample could not be identified. Irrespective of the histological subtype, only three b-values were needed for simplified IVIM, but these parameters did not retain their discriminative ability.ConclusionSubsampling of six b-values halved the IVIM scan time without significant losses in accuracy and discriminative ability. Simplified IVIM is possible with only three b-values, at the risk of losing diagnostic information.

Correction for fast pseudo-diffusive fluid motion contaminations in diffusion tensor imaging

In this prospective study, we quantified the fast pseudo-diffusion contamination by blood perfusion or cerebrospinal fluid (CSF) intravoxel incoherent movements on the measurement of the diffusion tensor metrics in healthy brain tissue.Diffusion-weighted imaging (TR/TE = 4100 ms/90 ms; b-values: 0, 5, 10, 20, 35, 55, 80, 110, 150, 200, 300, 500, 750, 1000, 1300 s/mm2, 20 diffusion-encoding directions) was performed on a cohort of five healthy volunteers at 3 Tesla. The projections of the diffusion tensor along each diffusion-encoding direction were computed using a two b-value approach (2b), by fitting the signal to a monoexponential curve (mono), and by correcting for fast pseudo-diffusion compartments using the biexponential intravoxel incoherent motion model (IVIM) (bi). Fractional anisotropy (FA) and mean diffusivity (MD) of the diffusion tensor were quantified in regions of interest drawn over white matter areas, gray matter areas, and the ventricles.A significant dependence of the MD from the evaluation method was found in all selected regions. A lower MD was computed when accounting for the fast-diffusion compartments. A larger dependence was found in the nucleus caudatus (bi: median 0.86 10−3 mm2/s, Δ2b: −11.2%, Δmono: −14.4%; p = 0.007), in the anterior horn (bi: median 2.04 10−3 mm2/s, Δ2b: −9.4%, Δmono: −11.5%, p = 0.007) and in the posterior horn of the lateral ventricles (bi: median 2.47 10−3 mm2/s, Δ2b: −5.5%, Δmono: −11.7%; p = 0.007). Also for the FA, the signal modeling affected the computation of the anisotropy metrics. The deviation depended on the evaluated region with significant differences mainly in the nucleus caudatus (bi: median 0.15, Δ2b: +39.3%, Δmono: +14.7%; p = 0.022) and putamen (bi: median 0.19, Δ2b: +3.1%, Δmono: +17.3%; p = 0.015).Fast pseudo-diffusive regimes locally affect diffusion tensor imaging (DTI) metrics in the brain. Here, we propose the use of an IVIM-based method for correction of signal contaminations through CSF or perfusion.

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

Studies on intravoxel incoherent motion (IVIM) imaging are carried out with different acquisition protocols. To investigate the dependence of IVIM parameters on the B0 field strength when using a bi- or triexponential model. Prospective. 20 healthy volunteers (age: 19-28 years). 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/mm2 . 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. One-way analysis of variance (ANOVA) and Kruskal-Wallis test were used to test the obtained IVIM parameters for a significant field strength dependency. At b-values below 6 s/mm2 , the triexponential model provided better agreement with the data than the biexponential model. The average tissue diffusivity was D = 1.22/1.00 μm2 /msec at 1.5/3T. The average pseudodiffusion coefficients for the biexponential model were D* = 308/260 μm2 /msec at 1.5/3T; and for the triexponential model = 81.3/65.9 μm2 /msec, = 2453/2333 μm2 /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 f1 = 0.161/0.174 and f2 = 0.152/0.159 at 1.5/3T. A significant B0 dependence was only found for the biexponential pseudodiffusion coefficient (ANOVA/KW P = 0.037/0.0453) and tissue diffusivity (ANOVA/KW: P < 0.001). Our experimental results suggest that triexponential pseudodiffusion coefficients and perfusion fractions obtained at different field strengths could be compared across different studies using different B0 . 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. 1 Technical Efficacy Stage: 2 J. Magn. Reson. Imaging 2019;50:1883-1892.

Auto-Regressive Discrete Acquisition Points Transformation for Diffusion Weighted MRI Data.

A new method for fitting diffusion-weighted magnetic resonance imaging (DW-MRI) data composed of an unknown number of multi-exponential components is presented and evaluated. The auto-regressive discrete acquisition points transformation (ADAPT) method is an adaption of the auto-regressive moving average system, which allows for the modeling of multi-exponential data and enables the estimation of the number of exponential components without prior assumptions. ADAPT was evaluated on simulated DW-MRI data. The optimum ADAPT fit was then applied to human brain DWI data and the correlation between the ADAPT coefficients and the parameters of the commonly used bi-exponential intravoxel incoherent motion (IVIM) method were investigated. The ADAPT method can correctly identify the number of components and model the exponential data. The ADAPT coefficients were found to have strong correlations with the IVIM parameters. ADAPT(1,1)-β0 correlated with IVIM-D: ρ = 0.708, P < 0.001. ADAPT(1,1)-α1 correlated with IVIM-f: ρ = 0.667, P < 0.001. ADAPT(1,1)-β1 correlated with IVIM-D*: ρ = 0.741, P < 0.001). ADAPT provides a method that can identify the number of exponential components in DWI data without prior assumptions, and determine potential complex diffusion biomarkers. ADAPT has the potential to provide a generalized fitting method for discrete multi-exponential data, and determine meaningful coefficients without prior information.

Differentiation of peripheral nerve sheath tumors in patients with neurofibromatosis type 1 using diffusion-weighted magnetic resonance imaging.

We sought to determine the value of diffusion-weighted (DW) magnetic resonance imaging (MRI) for characterization of benign and malignant peripheral nerve sheath tumors (PNSTs) in patients with neurofibromatosis type 1 (NF1). Twenty-six patients with NF1 and suspicion of malignant transformation of PNSTs were prospectively enrolled and underwent DW MRI at 3T. For a set of benign (n = 55) and malignant (n = 12) PNSTs, functional MRI parameters were derived from both biexponential intravoxel incoherent motion (diffusion coefficient D and perfusion fraction f) and monoexponential data analysis (apparent diffusion coefficients [ADCs]). A panel of morphological MRI features was evaluated using T1- and T2-weighted imaging. Mann-Whitney U-test, Fisher's exact test, and receiver operating characteristic (ROC) analyses were applied to assess the diagnostic accuracy of quantitative and qualitative MRI. Cohen's kappa was used to determine interrater reliability. Malignant PNSTs demonstrated significantly lower diffusivity (P < 0.0001) compared with benign PNSTs. The perfusion fraction f was significantly higher in malignant PNSTs (P < 0.001). In ROC analysis, functional MRI parameters showed high diagnostic accuracy for differentiation of PNSTs (eg, ADCmean, 92% sensitivity with 98% specificity, AUC 0.98; Dmean, 92% sensitivity with 98% specificity, AUC 0.98). By contrast, morphological imaging features had only limited sensitivity (18-94%) and specificity (18-82%) for identification of malignancy. Interrater reliability was higher for monoexponential data analysis. DW imaging shows better diagnostic performance than morphological features and allows accurate differentiation of benign and malignant peripheral nerve sheath tumors in NF1.

Image acquisition for intravoxel incoherent motion imaging of kidneys should be triggered at the instant of maximum blood velocity: evidence obtained with simulations and in vivo experiments.

To demonstrate that diffusion-weighted images should be acquired at the instant of maximum blood velocity in kidneys to extract the perfusion fraction (PF) by the bi-exponential intravoxel incoherent motion model. The PF values were measured in Monte-Carlo simulations corresponding to different blood velocities with a constant known PF. The distribution of the measured PF values (PF-distribution) was characterized quantitatively by 3 markers highlighting the deviation of the measurement from the true PF. Diffusion-weighted images of kidneys were acquired in 10 healthy volunteers at the instant of maximal respectively minimal blood velocity in the renal artery (Vmax versus Vmin acquisition). The PF-distributions measured from the Vmax and Vmin acquisitions were compared mutually and with simulated PF-distributions using the 3 markers. A radiologist evaluated the quality of the PF maps. The PF-distributions measured in the simulations were spread around the true PF value, and spreading was reduced as blood velocity increased. A comparison between simulated and in vivo PF-distributions suggests that a similar phenomenon is plausible in vivo. The quality of the PF maps of the Vmax -acquisition was scored higher by the radiologist than those of the Vmin -acquisition in 95% of cases (19 of 20). The PF maps are of better quality when the Vmax -acquisition is used. We show evidence supporting the hypothesis that the variation of PF along the cardiac cycle is due to oscillations between a poor estimation when the blood velocity is low, and a better estimation when blood velocity is higher.

Preoperative evaluation of pelvine lymph node metastasis in high risk prostate cancer with intravoxel incoherent motion (IVIM) MRI

ObjectivesTo evaluate benign and malignant pelvine lymph nodes in prostate cancer patients with biexponential intravoxel incoherent motion (IVIM) MRI of the prostate prior to radical prostatectomy. MethodsThe ethics committee approved this retrospective study with waiver of informed consent. From February 2012 to November 2013 43 patients with histopathologically proven prostatic cancer were included. All patients were examined applying a standardized MRI protocol including IVIM diffusion weighted imaging with multiple b-values ranging from 0 to 950 s/mm². MR imaging was performed one day prior to radical prostatectomy. Thereafter, extended lymph node resection was performed. For each MRI all visible lymph nodes were registered and calculated as individual regions of interest. These findings were correlated with postoperative pathology. The apparent diffusion coefficient ADC, the diffusion coefficient D and the perfusion fraction f were calculated from IVIM DWI using a biexponential fit. ResultsA total of 120 lymph nodes were detected on MRI. 95 of these were determined as benign and 25 as malignant. The average ADC was significantly lower in malignant compared to benign lymph nodes (0.88 × 10−³ vs 1.67 × 10−³ mm²/s, p < 0.001). Likewise, the average diffusion coefficient D was significantly lower in lymph node metastasis (0.54 × 10−³ vs 1.10 × 10−³ mm²/s, p < .001). The signal rate due to perfusion was significantly higher in malignant compared to benign nodes (33.4% vs. 27.1%, p = 0.02). ConclusionsApplying biexponential IVIM MRI demonstrates significant differences in diffusion parameters ADC and D, as well as in the perfusion fraction f for benign and malignant lymph nodes. Therefore, IVIM might help to further improve the preoperative assessment of lymph nodes in MRI.