Magnetic resonance imaging-based approaches for detecting the efficacy of combining therapy following VEGFR-2 and PD-1 blockade in a colon cancer model.

BACKGROUND: Angiogenesis inhibitors have been identified to improve the efficacy of immunotherapy in recent studies. However, the delayed therapeutic effect of immunotherapy poses challenges in treatment planning. Therefore, this study aims to explore the potential of non-invasive imaging techniques, specifically intravoxel-incoherent-motion diffusion-weighted imaging (IVIM-DWI) and blood oxygenation level-dependent magnetic resonance imaging (BOLD-MRI), in detecting the anti-tumor response to the combination therapy involving immune checkpoint blockade therapy and anti-angiogenesis therapy in a tumor-bearing animal model. METHODS: The C57BL/6 mice were implanted with murine MC-38 cells to establish colon cancer xenograft model, and randomly divided into the control group, anti-PD-1 therapy group, and combination therapy group (VEGFR-2 inhibitor combined with anti-PD-1 antibody treatment). All mice were imaged before and, on the 3rd, 6th, 9th, and 12th day after administration, and pathological examinations were conducted at the same time points. RESULTS: The combination therapy group effectively suppressed tumor growth, exhibiting a significantly higher tumor inhibition rate of 69.96% compared to the anti-PD-1 group (56.71%). The f value and D* value of IVIM-DWI exhibit advantages in reflecting tumor angiogenesis. The D* value showed the highest correlation with CD31 (r = 0.702, P = 0.001), and the f value demonstrated the closest correlation with vessel maturity (r = 0.693, P = 0.001). While the BOLD-MRI parameter, R2* value, shows the highest correlation with Hif-1α(r = 0.778, P < 0.001), indicating the capability of BOLD-MRI to evaluate tumor hypoxia. In addition, the D value of IVIM-DWI is closely related to tumor cell proliferation, apoptosis, and infiltration of lymphocytes. The D value was highly correlated with Ki-67 (r = - 0.792, P < 0.001), TUNEL (r = 0.910, P < 0.001) and CD8a (r = 0.918, P < 0.001). CONCLUSIONS: The combination of VEGFR-2 inhibitors with PD-1 immunotherapy shows a synergistic anti-tumor effect on the mouse colon cancer model. IVIM-DWI and BOLD-MRI are expected to be used as non-invasive approaches to provide imaging-based evidence for tumor response detection and efficacy evaluation.

Comparison of methods for intravoxel incoherent motion parameter estimation in the brain from flow-compensated and non-flow-compensated diffusion-encoded data.

PURPOSE: Joint analysis of flow-compensated (FC) and non-flow-compensated (NC) diffusion MRI (dMRI) data has been suggested for increased robustness of intravoxel incoherent motion (IVIM) parameter estimation. For this purpose, a set of methods commonly used or previously found useful for IVIM analysis of dMRI data obtained with conventional diffusion encoding were evaluated in healthy human brain. METHODS: Five methods for joint IVIM analysis of FC and NC dMRI data were compared: (1) direct non-linear least squares fitting, (2) a segmented fitting algorithm with estimation of the diffusion coefficient from higher b-values of NC data, (3) a Bayesian algorithm with uniform prior distributions, (4) a Bayesian algorithm with spatial prior distributions, and (5) a deep learning-based algorithm. Methods were evaluated on brain dMRI data from healthy subjects and simulated data at multiple noise levels. Bipolar diffusion encoding gradients were used with b-values 0-200 s/mm2 and corresponding flow weighting factors 0-2.35 s/mm for NC data and by design 0 for FC data. Data were acquired twice for repeatability analysis. RESULTS: Measurement repeatability as well as estimation bias and variability were at similar levels or better with the Bayesian algorithm with spatial prior distributions and the deep learning-based algorithm for IVIM parameters D $$ D $$ and f $$ f $$ , and for the Bayesian algorithm only for v d $$ {v}_d $$ , relative to the other methods. CONCLUSION: A Bayesian algorithm with spatial prior distributions is preferable for joint IVIM analysis of FC and NC dMRI data in the healthy human brain, but deep learning-based algorithms appear promising.

Velocity-compensated intravoxel incoherent motion of the human calf muscle.

PURPOSE: To determine whether intravoxel incoherent motion (IVIM) describes the blood perfusion in muscles better, assuming pseudo diffusion (Bihan Model 1) or ballistic motion (Bihan Model 2). METHODS: IVIM parameters were measured in 18 healthy subjects with three different diffusion gradient time profiles (bipolar with two diffusion times and one with velocity compensation) and 17 b-values (0-600 s/mm2) at rest and after muscle activation. The diffusion coefficient, perfusion fraction, and pseudo-diffusion coefficient were estimated with a segmented fit in the gastrocnemius medialis (GM) and tibialis anterior (TA) muscles. RESULTS: Velocity-compensated gradients resulted in a decreased perfusion fraction (6.9% ± 1.4% vs. 4.4% ± 1.3% in the GM after activation) and pseudo-diffusion coefficient (0.069 ± 0.046 mm2/s vs. 0.014 ± 0.006 in the GM after activation) compared to the bipolar gradients with the longer diffusion encoding time. Increased diffusion coefficients, perfusion fractions, and pseudo-diffusion coefficients were observed in the GM after activation for all gradient profiles. However, the increase was significantly smaller for the velocity-compensated gradients. A diffusion time dependence was found for the pseudo-diffusion coefficient in the activated muscle. CONCLUSION: Velocity-compensated diffusion gradients significantly suppress the IVIM effect in the calf muscle, indicating that the ballistic limit is mostly reached, which is supported by the time dependence of the pseudo-diffusion coefficient.

Intravoxel incoherent motion and diffusion kurtosis imaging and their machine-learning-based texture analysis for detection and assessment of prostate cancer severity at 3 T.

OBJECTIVES: To evaluate the role of combined intravoxel incoherent motion and diffusion kurtosis imaging (IVIM-DKI) and their machine-learning-based texture analysis for the detection and assessment of severity in prostate cancer (PCa). MATERIALS AND METHODS: Eighty-eight patients underwent MRI on a 3 T scanner after giving informed consent. IVIM-DKI data were acquired using 13 b values (0-2000 s/mm2) and analyzed using the IVIM-DKI model with the total variation (TV) method. PCa patients were categorized into two groups: clinically insignificant prostate cancer (CISPCa) (Gleason grade ≤ 6) and clinically significant prostate cancer (CSPCa) (Gleason grade ≥ 7). One-way analysis-of-variance, t test, and receiver operating characteristic analysis was performed to measure the discriminative ability to detect PCa using IVIM-DKI parameters. A chi-square test was used to select important texture features of apparent diffusion coefficient (ADC) and IVIM-DKI parameters. These selected texture features were used in an artificial neural network for PCa detection. RESULTS: ADC and diffusion coefficient (D) were significantly lower (p < 0.001), and kurtosis (k) was significantly higher (p < 0.001), in PCa as compared with benign prostatic hyperplasia (BPH) and normal peripheral zone (PZ). ADC, D, and k showed high areas under the curves (AUCs) of 0.92, 0.89, and 0.88, respectively, in PCa detection. ADC and D were significantly lower (p < 0.05) as compared with CISPCa versus CSPCa. D for detecting CSPCa was high, with an AUC of 0.63. A negative correlation of ADC and D with GS (ADC, ρ = -0.33; D, ρ = -0.35, p < 0.05) and a positive correlation of k with GS (ρ = 0.22, p < 0.05) were observed. Combined IVIM-DKI texture showed high AUC of 0.83 for classification of PCa, BPH, and normal PZ. CONCLUSION: D, f, and k computed using the IVIM-DKI model with the TV method were able to differentiate PCa from BPH and normal PZ. Texture features of combined IVIM-DKI parameters showed high accuracy and AUC in PCa detection.

Feasibility of Multiparametric Perfusion Assessment in Diabetic Foot Ulcer Using Intravoxel Incoherent Motion and Blood Oxygenation-Level Dependent MRI.

BACKGROUND: Patients with type-2 diabetes (T2DM) are at increased risk of developing diabetic foot ulcers (DFU) and experiencing impaired wound healing related to underlying microvascular disease. PURPOSE: To evaluate the sensitivity of intra-voxel incoherent motion (IVIM) and blood oxygen level dependent (BOLD) MRI to microvascular changes in patients with DFUs. STUDY TYPE: Case-control. POPULATION: 20 volunteers who were age and body mass index matched, including T2DM patients with DFUs (N = 10, mean age = 57.5 years), T2DM patients with controlled glycemia and without DFUs (DC, N = 5, mean age = 57.4 years) and healthy controls (HC, N = 5, mean age = 52.8 years). FIELD STRENGTH/SEQUENCE: 3T/multi-b-value IVIM and dynamic BOLD. ASSESSMENT: Resting IVIM parameters were obtained using a multi-b-value diffusion-weighted imaging sequence and two IVIM models were fit to obtain diffusion coefficient (D), pseudo-diffusion coefficient (D*), perfusion fraction (f) and microvascular volume fraction (MVF) parameters. Microvascular reactivity was evaluated by inducing an ischemic state in the foot with a blood pressure cuff during dynamic BOLD imaging. Perfusion indices were assessed in two regions of the foot: the medial plantar (MP) and lateral plantar (LP) regions. STATISTICAL TESTS: Effect sizes of group mean differences were assessed using Hedge's g adjusted for small sample sizes. RESULTS: DFU participants exhibited elevated D*, f, and MVF values in both regions (g ≥ 1.10) and increased D (g = 1.07) in the MP region compared to DC participants. DC participants showed reduced f and MVF compared to HC participants in the MP region (g ≥ 1.06). Finally, the DFU group showed reduced tolerance for ischemia in the LP region (g = -1.51) and blunted reperfusion response in both regions (g < -2.32) compared to the DC group during the cuff-occlusion challenge. DATA CONCLUSION: The combined use of IVIM and BOLD MRI shows promise in differentiating perfusion abnormalities in the feet of diabetic patients and suggests hyperperfusion in DFU patients. LEVEL OF EVIDENCE: 1 TECHNICAL EFFICACY: Stage 1.

Two-Compartment Perfusion MR IVIM Model to Investigate Normal and Pathological Placental Tissue.

BACKGROUND: Perfusion and diffusion coexist in the placenta and can be altered by pathologies. The two-perfusion model, where f1 and, f2 are the perfusion-fraction of the fastest and slowest perfusion compartment, respectively, and D is the diffusion coefficient, may help differentiate between normal and impaired placentas. PURPOSE: Investigate the potential of the two-perfusion IVIM model in differentiating between normal and abnormal placentas. STUDY-TYPE: Retrospective, case-control. POPULATION: 43 normal pregnancy, 9 fetal-growth-restriction (FGR), 6 small-for-gestational-age (SGA), 4 accreta, 1 increta and 2 percreta placentas. FIELD STRENGTH/SEQUENCE: Diffusion-weighted-echo planar imaging sequence at 1.5 T. ASSESSMENT: Voxel-wise signal-correction and fitting-controls were used to avoid overfitting obtaining that two-perfusion model fitted the observed data better than the IVIM model (Akaike weight: 0.94). The two-perfusion parametric-maps were quantified from ROIs in the fetal and maternal placenta and in the accretion zone of accreta placentas. The diffusion coefficient D was evaluated using a b ≥ 200 sec/mm2 -mono-exponential decay fit. IVIM metrics were quantified to fix f1 + f2 = fIVIM . STATISTICAL-TESTS: ANOVA with Dunn-Sidák's post-hoc correction and Cohen's d test were used to compare parameters between groups. Spearman's coefficient was evaluated to study the correlation between variables. A P-value<0.05 indicated a statistically significant difference. RESULTS: There was a significant difference in f1 between FGR and SGA, and significant differences in f2 and fIVIM between normal and FGR. The percreta + increta group showed the highest f1 values (Cohen's d = -2.66). The f2 between normal and percreta + increta groups showed Cohen's d = 1.12. Conversely, fIVIM had a small effective size (Cohen's d = 0.32). In the accretion zone, a significant correlation was found between f2 and GA (ρ = 0.90) whereas a significant negative correlation was found between fIVIM and D (ρ = -0.37 in fetal and ρ = -0.56 in maternal side) and f2 and D (ρ = -0.38 in fetal and ρ = -0.51 in maternal side) in normal placentas. CONCLUSION: The two-perfusion model provides complementary information to IVIM parameters that may be useful in identifying placenta impairment. LEVEL OF EVIDENCE: 2 TECHNICAL EFFICACY STAGE: 1.

Intravoxel incoherent motion magnetic resonance imaging in the assessment of brain microstructure and perfusion in idiopathic normal-pressure hydrocephalus.

PURPOSE: To determine the relationship between intravoxel incoherent motion (IVIM) MRI parameters and clinical changes post-tap test (TT) in idiopathic normal-pressure hydrocephalus (iNPH) patients. METHODS: Forty-four probable iNPH patients underwent 3 T MRI before and after TT. IVIM parameters were calculated from eight different bilateral regions of interest in basal ganglia, centrum semiovale, and corona radiata. Patients were categorized based on TT response into positive (group 1) and negative (group 2) groups. A Welch two-sample t-test was used to compare differences in D, D*, f, and ADC between the two groups, while a paired t-test was employed to assess the changes within each group before and after TT. These parameters were then correlated with clinical results. RESULTS: In the lenticular and thalamic nuclei, D value was significantly lower in the group 1 compared to group 2 both pre- and post-TT (p = 0.002 and p = 0.007 respectively). Post-TT, the positive response group exhibited a notably reduced D* value (p = 0.012) and significantly higher f values (p = 0.028). In the corona radiata and centrum semiovale, a significant post-TT reduction in D* was observed in the positive response group (p = 0.017). Within groups, the positive response cohort showed a significant post-TT increase in ADC (p < 0.001) and a decrease in D* (p = 0.007). CONCLUSION: IVIM permits the acquisition of important non-invasive information about tissue and vascularization in iNPH patients. Enhanced perfusion in the lenticular and thalamic nuclei may suggest the role of re-established microvascular and glymphatic pathways, potentially elucidating the functional improvement in motor function after TT in iNPH patients.

Intravoxel incoherent motion (IVIM)-derived perfusion fraction mapping for the visual evaluation of MR-guided high intensity focused ultrasound (MR-HIFU) ablation of uterine fibroids.

BACKGROUND: A method for periprocedural contrast agent-free visualization of uterine fibroid perfusion could potentially shorten magnetic resonance-guided high intensity focused ultrasound (MR-HIFU) treatment times and improve outcomes. Our goal was to test feasibility of perfusion fraction mapping by intravoxel incoherent motion (IVIM) modeling using diffusion-weighted MRI as method for visual evaluation of MR-HIFU treatment progression. METHODS: Conventional and T2-corrected IVIM-derived perfusion fraction maps were retrospectively calculated by applying two fitting methods to diffusion-weighted MRI data (b = 0, 50, 100, 200, 400, 600 and 800 s/mm2 at 1.5 T) from forty-four premenopausal women who underwent MR-HIFU ablation treatment of uterine fibroids. Contrast in perfusion fraction maps between areas with low perfusion fraction and surrounding tissue in the target uterine fibroid immediately following MR-HIFU treatment was evaluated. Additionally, the Dice similarity coefficient (DSC) was calculated between delineated areas with low IVIM-derived perfusion fraction and hypoperfusion based on CE-T1w. RESULTS: Average perfusion fraction ranged between 0.068 and 0.083 in areas with low perfusion fraction based on visual assessment, and between 0.256 and 0.335 in surrounding tissues (all p < 0.001). DSCs ranged from 0.714 to 0.734 between areas with low perfusion fraction and the CE-T1w derived non-perfused areas, with excellent intraobserver reliability of the delineated areas (ICC 0.97). CONCLUSION: The MR-HIFU treatment effect in uterine fibroids can be visualized using IVIM perfusion fraction mapping, in moderate concordance with contrast enhanced MRI. IVIM perfusion fraction mapping has therefore the potential to serve as a contrast agent-free imaging method to visualize the MR-HIFU treatment progression in uterine fibroids.

Spatial profiling of in vivo diffusion-weighted MRI parameters in the healthy human kidney.

OBJECTIVE: Diffusion-weighted MRI is a technique that can infer microstructural and microcirculatory features from biological tissue, with particular application to renal tissue. There is extensive literature on diffusion tensor imaging (DTI) of anisotropy in the renal medulla, intravoxel incoherent motion (IVIM) measurements separating microstructural from microcirculation effects, and combinations of the two. However, interpretation of these features and adaptation of more specific models remains an ongoing challenge. One input to this process is a whole organ distillation of corticomedullary contrast of diffusion metrics, as has been explored for other renal biomarkers. MATERIALS AND METHODS: In this work, we probe the spatial dependence of diffusion MRI metrics with concentrically layered segmentation in 11 healthy kidneys at 3 T. The metrics include those from DTI, IVIM, a combined approach titled "REnal Flow and Microstructure AnisotroPy (REFMAP)", and a multiply encoded model titled "FC-IVIM" providing estimates of fluid velocity and branching length. RESULTS: Fractional anisotropy decreased from the inner kidney to the outer kidney with the strongest layer correlation in both parenchyma (including cortex and medulla) and medulla with Spearman correlation coefficients and p-values (r, p) of (0.42, <0.001) and (0.37, <0.001), respectively. Also, dynamic parameters derived from the three models significantly decreased with a high correlation from the inner to the outer parenchyma or medulla with (r, p) ranges of (0.46-0.55, <0.001). CONCLUSIONS: These spatial trends might find implications for indirect assessments of kidney physiology and microstructure using diffusion MRI.

IDEAL-IQ combined with intravoxel incoherent motion diffusion-weighted imaging for quantitative diagnosis of osteoporosis.

BACKGROUND: Osteoporosis (OP) is a common chronic metabolic bone disease characterized by decreased bone mineral content and microstructural damage, leading to increased fracture risk. Traditional methods for measuring bone density have limitations in accurately distinguishing vertebral bodies and are influenced by vertebral degeneration and surrounding tissues. Therefore, novel methods are needed to quantitatively assess changes in bone density and improve the accurate diagnosis of OP. METHODS: This study aimed to explore the applicative value of the iterative decomposition of water and fat with echo asymmetry and least-squares estimation-iron (IDEAL-IQ) sequence combined with intravoxel incoherent motion diffusion-weighted imaging (IVIM-DWI) for the diagnosis of osteoporosis. Data from 135 patients undergoing dual-energy X-ray absorptiometry (DXA), IDEAL-IQ, and IVIM-DWI were prospectively collected and analyzed. Various parameters obtained from IVIM-DWI and IDEAL-IQ sequences were compared, and their diagnostic efficacy was evaluated. RESULTS: Statistically significant differences were observed among the three groups for FF, R2*, f, D, DDC values, and BMD values. FF and f values exhibited negative correlations with BMD values, with r=-0.313 and - 0.274, respectively, while R2*, D, and DDC values showed positive correlations with BMD values, with r = 0.327, 0.532, and 0.390, respectively. Among these parameters, D demonstrated the highest diagnostic efficacy for osteoporosis (AUC = 0.826), followed by FF (AUC = 0.713). D* exhibited the lowest diagnostic performance for distinguishing the osteoporosis group from the other two groups. Only D showed a significant difference between genders. The AUCs for IDEAL-IQ, IVIM-DWI, and their combination were 0.74, 0.89, and 0.90, respectively. CONCLUSIONS: IDEAL-IQ combined with IVIM-DWI provides valuable information for the diagnosis of osteoporosis and offers evidence for clinical decisions. The superior diagnostic performance of IVIM-DWI, particularly the D value, suggests its potential as a more sensitive and accurate method for diagnosing osteoporosis compared to IDEAL-IQ. These findings underscore the importance of integrating advanced imaging techniques into clinical practice for improved osteoporosis management and highlight the need for further research to explore the full clinical implications of these imaging modalities.

Contralateral renal change in a unilateral ureteral obstruction rat model using intravoxel incoherent motion diffusion-weighted imaging.

OBJECTIVES: Most functional magnetic resonance research has primarily examined alterations in the affected kidney, often neglecting the contralateral kidney. Our study aims to investigate whether imaging parameters accurately depict changes in both the renal cortex and medulla in a unilateral ureteral obstruction rat model, thereby showcasing the utility of intravoxel incoherent motion (IVIM) in evaluating contralateral renal changes. METHODS: Six rats underwent MR scans and were subsequently sacrificed for baseline histological examination. Following the induction of left ureteral obstruction, 48 rats were scanned, and the histopathological examinations were conducted on days 3, 7, 10, 14, 21, 28, 35, and 42. The apparent diffusion coefficient (ADC), pure molecular diffusion (D), pseudodiffusion (D*), and perfusion fraction (f) values were measured using IVIM. RESULTS: On the 10th day of obstruction, both cortical and medullary ADC values differed significantly between the UUO10 group and the sham group (p < 0.01). The cortical D values showed statistically significant differences between UUO3 group and sham group (p < 0.01) but not among UUO groups at other time point. Additionally, the cortical and medullary f values were statistically significant between the UUO21 group and the sham group (p < 0.01). Especially, the cortical f values exhibited significant differences between the UUO21 group and the UUO groups with shorter obstruction time (at time point of 3, 7, 10, 14 day) (p < 0.01). CONCLUSIONS: Significant hemodynamic alterations were observed in the contralateral kidney following renal obstruction. IVIM accurately captures changes in the unobstructed kidney. Particularly, the cortical f value exhibits the highest potential for assessing contralateral renal modifications.

Image downsampling expedited adaptive least-squares (IDEAL) fitting improves intravoxel incoherent motion (IVIM) analysis in the human kidney.

PURPOSE: To improve the reliability of intravoxel incoherent motion (IVIM) model parameter estimation for the DWI in the kidney using a novel image downsampling expedited adaptive least-squares (IDEAL) approach. METHODS: The robustness of IDEAL was investigated using simulated DW-MRI data corrupted with different levels of Rician noise. Subsequently, the performance of the proposed method was tested by fitting bi- and triexponential IVIM model to in vivo renal DWI data acquired on a clinical 3 Tesla MRI scanner and compared to conventional approaches (fixed D* and segmented fitting). RESULTS: The numerical simulations demonstrated that the IDEAL algorithm provides robust estimates of the IVIM parameters in the presence of noise (SNR of 20) as indicated by relatively low absolute percentage bias (maximal sMdPB <20%) and normalized RMSE (maximal RMSE <28%). The analysis of the in vivo data showed that the IDEAL-based IVIM parameter maps were less noisy and more visually appealing than those obtained using the fixed D* and segmented methods. Further, coefficients of variation for nearly all IVIM parameters were significantly reduced in cortex and medulla for IDEAL-based biexponential (coefficients of variation: 4%-50%) and triexponential (coefficients of variation: 7.5%-75%) IVIM modelling compared to the segmented (coefficients of variation: 4%-120%) and fixed D* (coefficients of variation: 17%-174%) methods, reflecting greater accuracy of this method. CONCLUSION: The proposed fitting algorithm yields more robust IVIM parameter estimates and is less susceptible to poor SNR than the conventional fitting approaches. Thus, the IDEAL approach has the potential to improve the reliability of renal DW-MRI analysis for clinical applications.

Influence of saturation effects on biexponential liver intravoxel incoherent motion.

PURPOSE: 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. METHODS: Fifteen healthy volunteers (21-30 years) were examined at a field strength of 3 T. 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). RESULTS: None of the parameters were significantly different between the settings. For few slices and many slices, respectively, the mean values (SDs) for D $$ D $$ were 1.21 µm 2 / ms $$ 1.21{\upmu \mathrm{m}}^2/\mathrm{ms} $$ ( 0.19 µm 2 / ms $$ 0.19\kern0.3em {\upmu \mathrm{m}}^2/\mathrm{ms} $$ ) and 1.20 µm 2 / ms $$ 1.20{\upmu \mathrm{m}}^2/\mathrm{ms} $$ ( 0.11 µm 2 / ms $$ 0.11\kern0.3em {\upmu \mathrm{m}}^2/\mathrm{ms} $$ ); for f $$ f $$ they were 29.7% (6.2%) and 27.7% (3.6%); and for D * $$ {D}^{\ast } $$ they were 8.76 ⋅ 10 - 2 mm 2 / s $$ 8.76\cdot {10}^{-2}{\mathrm{mm}}^2/\mathrm{s} $$ ( 4.54 ⋅ 10 - 2 mm 2 / s $$ 4.54\cdot {10}^{-2}\kern0.3em {\mathrm{mm}}^2/\mathrm{s} $$ ) and 8.71 ⋅ 10 - 2 mm 2 / s $$ 8.71\cdot {10}^{-2}{\mathrm{mm}}^2/\mathrm{s} $$ ( 4.06 ⋅ 10 - 2 mm 2 / s $$ 4.06\cdot {10}^{-2}\kern0.3em {\mathrm{mm}}^2/\mathrm{s} $$ ). CONCLUSION: 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.

An optimized b-value sampling for the quantification of interstitial fluid using diffusion-weighted MRI, a genetic algorithm approach.

PURPOSE: Multi-b-value diffusion-weighted MRI techniques can simultaneously measure the parenchymal diffusivity, microvascular perfusion, and a third, intermediate diffusion component. This component is related to the interstitial fluid in the brain parenchyma. However, simultaneously estimating three diffusion components from multi-b-value data is difficult and has strong dependence on SNR and chosen b-values. As the number of acquired b-values is limited due to scanning time, it is important to know which b-values are most effective to be included. Therefore, this study evaluates an optimized b-value sampling for interstitial fluid estimation. METHOD: The optimized b-value sampling scheme is determined using a genetic algorithm. Subsequently, the performance of this optimized sampling is assessed by comparing it with a linear, logarithmic, and previously proposed sampling scheme, in terms of the RMS error (RMSE) for the intermediate component estimation. The in vivo performance of the optimized sampling is assessed using 7T data with 101 equally spaced b-values ranging from 0 to 1000 s/mm2 . In this case, the RMSE was determined by comparing the fit that includes all b-values. RESULTS: The optimized b-value sampling for estimating the intermediate component was reported to be [0, 30, 90, 210, 280, 350, 580, 620, 660, 680, 720, 760, 980, 990, 1000] s/mm2 . For computer simulations, the optimized sampling had a lower RMSE, compared with the other samplings for varying levels of SNR. For the in vivo data, the voxel-wise RMSE of the optimized sampling was lower compared with other sampling schemes. CONCLUSION: The genetic algorithm-optimized b-value scheme improves the quantification of the diffusion component related to interstitial fluid in terms of a lower RMSE.

Registration and quantification network (RQnet) for IVIM-DKI analysis in MRI.

PURPOSE: A deep learning method is proposed for aligning diffusion weighted images (DWIs) and estimating intravoxel incoherent motion-diffusion kurtosis imaging parameters simultaneously. METHODS: We propose an unsupervised deep learning method that performs 2 tasks: registration and quantification for intravoxel incoherent motion-diffusion kurtosis imaging analysis. A common registration method in diffusion MRI is based on minimizing dissimilarity between various DWIs, which may result in registration errors due to different contrasts in different DWIs. We designed a novel unsupervised deep learning method for both accurate registration and quantification of various diffusion parameters. In order to generate motion-simulated training data and test data, 17 volunteers were scanned without moving their heads, and 4 volunteers moved their heads during the scan in a 3 Tesla MRI. In order to investigate the applicability of the proposed method to other organs, kidney images were also obtained. We compared the registration accuracy of the proposed method, statistical parametric mapping, and a deep learning method with a normalized cross-correlation loss. In the quantification part of the proposed method, a deep learning method that considered the diffusion gradient direction was used. RESULTS: Simulations and experimental results showed that the proposed method accurately performed registration and quantification for intravoxel incoherent motion-diffusion kurtosis imaging analysis. The registration accuracy of the proposed method was high for all b values. Furthermore, quantification performance was analyzed through simulations and in vivo experiments, where the proposed method showed the best performance among the compared methods. CONCLUSION: The proposed method aligns the DWIs and accurately quantifies the intravoxel incoherent motion-diffusion kurtosis imaging parameters.

Deep learning intravoxel incoherent motion modeling: Exploring the impact of training features and learning strategies.

PURPOSE: The development of advanced estimators for intravoxel incoherent motion (IVIM) modeling is often motivated by a desire to produce smoother parameter maps than least squares (LSQ). Deep neural networks show promise to this end, yet performance may be conditional on a myriad of choices regarding the learning strategy. In this work, we have explored potential impacts of key training features in unsupervised and supervised learning for IVIM model fitting. METHODS: Two synthetic data sets and one in-vivo data set from glioma patients were used in training of unsupervised and supervised networks for assessing generalizability. Network stability for different learning rates and network sizes was assessed in terms of loss convergence. Accuracy, precision, and bias were assessed by comparing estimations against ground truth after using different training data (synthetic and in vivo). RESULTS: A high learning rate, small network size, and early stopping resulted in sub-optimal solutions and correlations in fitted IVIM parameters. Extending training beyond early stopping resolved these correlations and reduced parameter error. However, extensive training resulted in increased noise sensitivity, where unsupervised estimates displayed variability similar to LSQ. In contrast, supervised estimates demonstrated improved precision but were strongly biased toward the mean of the training distribution, resulting in relatively smooth, yet possibly deceptive parameter maps. Extensive training also reduced the impact of individual hyperparameters. CONCLUSION: Voxel-wise deep learning for IVIM fitting demands sufficiently extensive training to minimize parameter correlation and bias for unsupervised learning, or demands a close correspondence between the training and test sets for supervised learning.

Model-based reconstructions for intravoxel incoherent motion and diffusion tensor imaging parameter map estimations.

Intravoxel incoherent motion (IVIM) imaging and diffusion tensor imaging (DTI) facilitate noninvasive quantification of tissue perfusion and diffusion. Both are promising biomarkers in various diseases and a combined acquisition is therefore desirable. This comes with challenges, including noisy parameter maps and long scan times, especially for the perfusion fraction f and pseudo-diffusion coefficient D*. A model-based reconstruction has the potential to overcome these challenges. As a first step, our goal was to develop a model-based reconstruction framework for IVIM and combined IVIM-DTI parameter estimation. The IVIM and IVIM-DTI models were implemented in the PyQMRI model-based reconstruction framework and validated with simulations and in vivo data. Commonly used voxel-wise nonlinear least-squares fitting was used as the reference. Simulations with the IVIM and IVIM-DTI models were performed with 100 noise realizations to assess accuracy and precision. Diffusion-weighted data were acquired for IVIM reconstruction in the liver (n = 5), as well as for IVIM-DTI in the kidneys (n = 5) and lower-leg muscles (n = 6) of healthy volunteers. The median and interquartile range (IQR) values of the IVIM and IVIM-DTI parameters were compared to assess bias and precision. With model-based reconstruction, the parameter maps exhibited less noise, which was most pronounced in the f and D* maps, both in the simulations and in vivo. The bias values in the simulations were comparable between model-based reconstruction and the reference method. The IQR was lower with model-based reconstruction compared with the reference for all parameters. In conclusion, model-based reconstruction is feasible for IVIM and IVIM-DTI and improves the precision of the parameter estimates, particularly for f and D* maps.

Cardiac Phase and Flow Compensation Effects on REnal Flow and Microstructure AnisotroPy MRI in Healthy Human Kidney.

BACKGROUND: Renal diffusion-weighted imaging (DWI) involves microstructure and microcirculation, quantified with diffusion tensor imaging (DTI), intravoxel incoherent motion (IVIM), and hybrid models. A better understanding of their contrast may increase specificity. PURPOSE: To measure modulation of DWI with cardiac phase and flow-compensated (FC) diffusion gradient waveforms. STUDY TYPE: Prospective. POPULATION: Six healthy volunteers (ages: 22-48 years, five females), water phantom. FIELD STRENGTH/SEQUENCE: 3-T, prototype DWI sequence with 2D echo-planar imaging, and bipolar (BP) or FC gradients. 2D Half-Fourier Single-shot Turbo-spin-Echo (HASTE). Multiple-phase 2D spoiled gradient-echo phase contrast (PC) MRI. ASSESSMENT: BP and FC water signal decays were qualitatively compared. Renal arteries and velocities were visualized on PC-MRI. Systolic (peak velocity), diastolic (end stable velocity), and pre-systolic (before peak velocity) phases were identified. Following mutual information-based retrospective self-registration of DWI within each kidney, and Marchenko-Pastur Principal Component Analysis (MPPCA) denoising, combined IVIM-DTI analysis estimated mean diffusivity (MD), fractional anisotropy (FA), and eigenvalues (λi) from tissue diffusivity (Dt ), perfusion fraction (fp ), and pseudodiffusivity (Dp , Dp,axial , Dp,radial ), for each tissue (cortex/medulla, segmented on b0/FA respectively), phase, and waveform (BP, FC). Monte Carlo water diffusion simulations aided data interpretation. STATISTICAL TESTS: Mixed model regression probed differences between tissue types and pulse sequences. Univariate general linear model analysis probed variations among cardiac phases. Spearman correlations were measured between diffusion metrics and renal artery velocities. Statistical significance level was set at P < 0.05. RESULTS: Water BP and FC signal decays showed no differences. Significant pulse sequence dependence occurred for λ1 , λ3 , FA, Dp , fp , Dp,axial , Dp,radial in cortex and medulla, and medullary λ2 . Significant cortex/medulla differences occurred with BP for all metrics except MD (systole [P = 0.224]; diastole [P = 0.556]). Significant phase dependence occurred for Dp , Dp,axial , Dp,radial for BP and medullary λ1 , λ2 , λ3 , MD for FC. FA correlated significantly with velocity. Monte Carlo simulations indicated medullary measurements were consistent with a 34 µm tubule diameter. DATA CONCLUSION: Cardiac gating and flow compensation modulate of measurements of renal diffusion. EVIDENCE LEVEL: 2 TECHNICAL EFFICACY STAGE: 2.

Application of Intravoxel Incoherent Motion in Clinical Liver Imaging: A Literature Review.

Intravoxel incoherent motion (IVIM) modeling is a widely used double-exponential model for describing diffusion-weighted imaging (DWI) signal, with a slow component related to pure molecular diffusion and a fast component associated with microcirculatory perfusion, which compensates for the limitations of traditional DWI. IVIM is a noninvasive technique for obtaining liver pathological information and characterizing liver lesions, and has potential applications in the initial diagnosis and treatment monitoring of liver diseases. Recent studies have demonstrated that IVIM-derived parameters are useful for evaluating liver lesions, including nonalcoholic fatty liver disease (NAFLD), liver fibrosis and liver tumors. However, the results are not stable. Therefore, it is necessary to summarize the current applications of IVIM in liver disease research, identify existing shortcomings, and point out the future development direction. In this review, we searched for studies related to hepatic IVIM-DWI applications over the past two decades in the PubMed database. We first introduce the fundamental principles and influential factors of IVIM, and then discuss its application in NAFLD, liver fibrosis, and focal hepatic lesions. It has been found that IVIM is still unstable in ensuring the robustness and reproducibility of measurements in the assessment of liver fibrosis grade and liver tumors differentiation, due to inconsistent and substantial overlap in the range of IVIM-derived parameters for different fibrotic stages. In the end, the future direction of IVIM-DWI in the assessment of liver diseases is discussed, emphasizing the need for further research on the stability of IVIM-derived parameters, particularly perfusion-related parameters, in order to promote the clinical practice of IVIM-DWI. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY: Stage 3.

Multiparametric MRI of Solid Renal Masses: Principles and Applications of Advanced Quantitative and Functional Methods for Tumor Diagnosis and Characterization.

Solid renal masses (SRMs) are increasingly detected and encompass both benign and malignant masses, with renal cell carcinoma (RCC) being the most common malignant SRM. Most patients with SRMs will undergo management without a priori pathologic confirmation. There is an unmet need to noninvasively diagnose and characterize RCCs, as significant variability in clinical behavior is observed and a wide range of differing management options exist. Cross-sectional imaging modalities, including magnetic resonance imaging (MRI), are increasingly used for SRM characterization. Multiparametric (mp) MRI techniques can provide insight into tumor biology by probing different physiologic/pathophysiologic processes noninvasively. These include sequences that probe tissue microstructure, including intravoxel incoherent motion diffusion-weighted imaging (IVIM-DWI) and T1 relaxometry; oxygen metabolism (blood oxygen level dependent [BOLD-MRI]); as well as vascular flow and perfusion (dynamic contrast-enhanced MRI [DCE-MRI] and arterial spin labeling [ASL]). In this review, we will discuss each mpMRI method in terms of its principles, roles, and discuss the results of human studies for SRM assessment. Future validation of these methods may help to enable a personalized management approach for patients with SRM in the emerging era of precision medicine. EVIDENCE LEVEL: 5. TECHNICAL EFFICACY: 2.