Non-parametric intravoxel incoherent motion analysis in patients with intracranial lesions: Test-retest reliability and correlation with arterial spin labeling

Intravoxel Incoherent Motion Quantification Dependent on Measurement SNR and Tissue Perfusion: A Simulation Study

Background

The intravoxel incoherent motion (IVIM) model extracts both functional and structural information of a tissue using motion-sensitizing gradients.

Objective

The Objective of the present work is to investigate the impact of signal to noise ratio (SNR) and physiologic conditions on the validity of IVIM parameters.

Material and Methods

This study is a simulation study, modeling IVIM at a voxel, and also done 10,000 times for every single simulation. Complex noises with various standard deviations were added to signal in-silico to investigate SNR effects on output validity. Besides, some blood perfusion situations for different tissues were considered based on their physiological range to explore the impacts of blood fraction at each voxel on the validity of the IVIM outputs. Coefficient variation (CV) and bias of the estimations were computed to assess the validity of the IVIM parameters.

Results

This study has shown that the validity of IVIM output parameters highly depends on measurement SNR and physiologic characteristics of the studied organ.

Conclusion

IVIM imaging could be useful if imaging parameters are correctly selected for each specific organ, considering hardware limitations.

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.

Signal to noise and b-value analysis for optimal intra-voxel incoherent motion imaging in the brain

Intravoxel incoherent motion (IVIM) is a method that can provide quantitative information about perfusion in the human body, in vivo, and without contrast agent. Unfortunately, the IVIM perfusion parameter maps are known to be relatively noisy in the brain, in particular for the pseudo-diffusion coefficient, which might hinder its potential broader use in clinical applications. Therefore, we studied the conditions to produce optimal IVIM perfusion images in the brain. IVIM imaging was performed on a 3-Tesla clinical system in four healthy volunteers, with 16 b values 0, 10, 20, 40, 80, 110, 140, 170, 200, 300, 400, 500, 600, 700, 800, 900 s/mm2, repeated 20 times. We analyzed the noise characteristics of the trace images as a function of b-value, and the homogeneity of the IVIM parameter maps across number of averages and sub-sets of the acquired b values. We found two peaks of noise of the trace images as function of b value, one due to thermal noise at high b-value, and one due to physiological noise at low b-value. The selection of b value distribution was found to have higher impact on the homogeneity of the IVIM parameter maps than the number of averages. Based on evaluations, we suggest an optimal b value acquisition scheme for a 12 min scan as 0 (7), 20 (4), 140 (19), 300 (9), 500 (19), 700 (1), 800 (4), 900 (1) s/mm2.

Brain perfusion imaging in neonates

•

MRI is the modality of choice to image and quantify cerebral perfusion.

•

Imaging of neonatal brain perfusion is possible using MRI and ultrasound.

•

Novel ultrafast ultrasound imaging allows for excellent spatiotemporal resolution.

•

Understanding cerebral hemodynamic changes of neonatal adaptation is key.

Abnormal variations of the neonatal brain perfusion can result in long-term neurodevelopmental consequences and cerebral perfusion imaging can play an important role in diagnostic and therapeutic decision-making. To identify at-risk situations, perfusion imaging of the neonatal brain must accurately evaluate both regional and global perfusion. To date, neonatal cerebral perfusion assessment remains challenging. The available modalities such as magnetic resonance imaging (MRI), ultrasound imaging, computed tomography (CT), near-infrared spectroscopy or nuclear imaging have multiple compromises and limitations. Several promising methods are being developed to achieve better diagnostic accuracy and higher robustness, in particular using advanced MRI and ultrasound techniques.

The objective of this state-of-the-art review is to analyze the methodology and challenges of neonatal brain perfusion imaging, to describe the currently available modalities, and to outline future perspectives.

Impact of b-Value Sampling Scheme on Brain IVIM Parameter Estimation in Healthy Subjects

Purpose:

Intravoxel incoherent motion (IVIM) analysis has attracted the interest of the clinical community due to its close relationship with microperfusion. Nevertheless, there is no clear reference protocol for its implementation; one of the questions being which b-value distribution to use. This study aimed to stress the importance of the sampling scheme and to show that an optimized b-value distribution decreases the variance associated with IVIM parameters in the brain with respect to a regular distribution in healthy volunteers.

Methods:

Ten volunteers were included in this study; images were acquired on a 1.5T MR scanner. Two distributions of 16 b-values were used: one considered ‘regular’ due to its close association with that used in other studies, and the other considered ‘optimized’ according to previous studies. IVIM parameters were adjusted according to the bi-exponential model, using two-step method. Analysis was undertaken in ROI defined using in the Automated Anatomical Labeling atlas, and parameters distributions were compared in a total of 832 ROI.

Results:

Maps with fewer speckles were obtained with the ‘optimized’ distribution. Coefficients of variation did not change significantly for the estimation of the diffusion coefficient D but decreased by approximately 39% for the pseudo-diffusion coefficient estimation and by 21% for the perfusion fraction. Distributions of adjusted parameters were found significantly different in 50% of the cases for the perfusion fraction, in 80% of the cases for the pseudo-diffusion coefficient and 17% of the cases for the diffusion coefficient. Observations across brain areas show that the range of average values for IVIM parameters is smaller in the ‘optimized’ case.

Conclusion:

Using an optimized distribution, data are sampled in a way that the IVIM signal decay is better described and less variance is obtained in the fitted parameters. The increased precision gained could help to detect small variations in IVIM parameters.

Impact of hypertension on cerebral microvascular structure in CPAP-treated obstructive sleep apnoea patients: a diffusion magnetic resonance imaging study

PURPOSE

Obstructive sleep apnoea (OSA) is a highly prevalent sleep-related breathing disorder associated with hypertension, impaired peripheral vascular function and an increased risk of stroke. Evidence suggests that abnormalities of the cerebral microcirculation, such as capillary rarefication, may be present in these patients. We evaluated whether the presence of hypertension may affect the cerebral capillary architecture and function assessed by Intravoxel Incoherent Motion (IVIM) magnetic resonance imaging (MRI) in patients with continuous positive airway pressure (CPAP)-treated OSA.

METHODS

Forty-one patients (88% male, mean age 57 ± 10 years) with moderate-to-severe OSA were selected and divided into two groups (normotensive vs. hypertensive). All hypertensive OSA patients were adherent with their antihypertensive medication. Cerebral microvascular structure was assessed in grey (GM) and white matter (WM) using an echo-planar diffusion imaging sequence with 14 different b values. A step-wise IVIM analysis algorithm was applied to compute true diffusion (D), perfusion fraction (f) and pseudo-diffusion (D*) values. Group comparisons were performed with the Wilcoxon-Mann-Whitney-Test. Regression analysis was adjusted for age.

RESULTS

Diffusion- and perfusion-related indexes in middle-aged OSA normotensive patients were quantified in both tissue types (D [10^-3 mm²/s]: GM = 0.83 ± 0.03; WM = 0.72 ± 0.03; f (%) GM = 0.09 ± 0.01; WM = 0.06 ± 0.01; D* [10^-3 mm²/s]: GM = 7.72 ± 0.89; WM = 7.38 ± 0.98). In the examined tissue types, hypertension did not result in changes on the estimated MRI IVIM index values.

CONCLUSION

Based on IVIM analysis, cerebral microvascular structure and function showed no difference between hypertensive and normotensive patients with moderate-to-severe OSA treated with CPAP. Treatment adherence with antihypertensive drug regime and, in turn, controlled hypertension seems not to affect microvascular structure and perfusion of the brain.

TRIAL REGISTRATION

ClinicalTrials.gov Identifier: NCT02493673.

Automated pixel-wise brain tissue segmentation of diffusion-weighted images via machine learning

The diffusion-weighted (DW) MR signal sampled over a wide range of b-values potentially allows for tissue differentiation in terms of cellularity, microstructure, perfusion, and T₂ relaxivity. This study aimed to implement a machine learning algorithm for automatic brain tissue segmentation from DW-MRI datasets, and to determine the optimal sub-set of features for accurate segmentation. DWI was performed at 3 T in eight healthy volunteers using 15 b-values and 20 diffusion-encoding directions. The pixel-wise signal attenuation, as well as the trace and fractional anisotropy (FA) of the diffusion tensor, were used as features to train a support vector machine classifier for gray matter, white matter, and cerebrospinal fluid classes. The datasets of two volunteers were used for validation. For each subject, tissue classification was also performed on 3D T₁ -weighted data sets with a probabilistic framework. Confusion matrices were generated for quantitative assessment of image classification accuracy in comparison with the reference method. DWI-based tissue segmentation resulted in an accuracy of 82.1% on the validation dataset and of 82.2% on the training dataset, excluding relevant model over-fitting. A mean Dice coefficient (DSC) of 0.79 ± 0.08 was found. About 50% of the classification performance was attributable to five features (i.e. signal measured at b-values of 5/10/500/1200 s/mm² and the FA). This reduced set of features led to almost identical performances for the validation (82.2%) and the training (81.4%) datasets (DSC = 0.79 ± 0.08). Machine learning techniques applied to DWI data allow for accurate brain tissue segmentation based on both morphological and functional information.

Measurement and scan reproducibility of parameters of intravoxel incoherent motion in renal tumor and normal renal parenchyma: a preliminary research at 3.0 T MR

PURPOSE

To prospectively estimate measurement and scan reproducibility of parameters of intravoxel incoherent motion (IVIM) in renal tumors, normal renal cortex, and medulla.

METHODS

Twenty-four consecutive patients (twelve males and twelve females; median age 56.7 years, range 32-71 years) with 25 renal tumors (20 renal cell carcinomas, one urothelium carcinoma, three angiomyolipomas, and one oncocytoma) were examined twice using IVIM₁ and IVIM₂ with 9 and 16 b values, respectively, at 3.0 T. All the patients were re-scanned in 24-48 h. Regions of interest (ROIs) were placed in solid part of tumor, normal cortex, and medulla to derive IVIM parameters D (true diffusion coefficient), D* (pseudodiffusion coefficient), and f (perfusion fraction of pseudodiffusion). Differences in parameters between two IVIM sets and intra-observer, inter-observer, and scan-rescan differences were assessed using paired t tests. Intra-observer, inter-observer, and scan-rescan reproducibility were assessed by measuring coefficient of variation and Bland-Altman limits of agreements.

RESULTS

Intra-observer reproducibility of renal tumors, normal renal cortex, and medulla was excellent for apparent diffusion coefficient (ADC; CV: 3.45%-5.34%, BA-LA: -14% to 18%) and D (CV: 3.65% to 6.04%, BA-LA: -18% to 19%), good for f (CV: 11.96%-16.08%, BA-LA: -76.4% to 92.1% except f of medulla with CV of 32.59% and BA-LA of -76.4% to 92.1% in IVIM₁), and poor for D* (CV: 25.0% to 75.4%, BA-LA: -111% to 150%). The same order was in inter-observer reproducibility analysis. Scan-rescan reproducibility was the worst of the three parameters. Renal medulla showed worse reproducibility than renal tumors and the normal cortex. The metrics of IVIM₂ had better reproducibility than IVIM₁.

CONCLUSION

Excellent reproducibility evaluation for ADC and D, good for f, and poor for D* derived from IVIM was performed in renal tumors, normal renal cortex, and medulla. D* has limited reliability and scan-rescan reproducibility should be improved.

Intravoxel incoherent motion diffusion-weighted magnetic resonance imaging in characterization of axillary lymph nodes: Preliminary animal experience

PURPOSE

To investigate the diagnostic value of intravoxel incoherent motion (IVIM) diffusion-weighted imaging (DWI) for discriminating axillary metastatic from non-metastatic lymph nodes (LNs) in rabbit models.

MATERIALS AND METHODS

The institutional animal care and use committee approved this study. Forty New Zealand white rabbits were randomly divided into two groups. The axillary LN models were created by inoculating VX2 cell suspension and complete Freund's adjuvant in the mammary glands of 20 female rabbits of each group, respectively. Conventional MRI and IVIM DWI were performed after animal models successfully established. Images of axillary LNs were analyzed with regard to long-axis diameter (L), short-axis diameter (S), apparent diffusion coefficient (ADC) and IVIM parameters (D, D*, f). Receiver operating characteristic analyses were conducted to determine the diagnostic performance of aforementioned criteria.

RESULTS

A total of 42 metastatic and 30 non-metastatic LNs were successfully isolated. ADC and D of metastatic LNs were significantly lower than those of non-metastatic ones (all P < 0.001), whereas D* was statistically higher (P = 0.033). L, S, and f showed no significant difference between the two groups (P = 0.089, 0.058, 0.054, respectively). Optimal cutoff values, area under the curve, sensitivity, and specificity for differentiation were as follows: ADC = 1.101 × 10^-3 mm²/s, 0.886, 78.6%, 90.0%; D = 0.938 × 10^-3 mm²/s, 0.927, 83.3%, 93.3%; and D* = 12.635 × 10^-3 mm²/s, 0.657, 52.4%, 80.0%.

CONCLUSION

IVIM DWI is useful to distinguish metastatic from non-metastatic LNs in axilla. D was the most discriminative variable for predicting metastatic LNs.

Cardiac-gated intravoxel incoherent motion diffusion-weighted magnetic resonance imaging for the investigation of intracranial cerebrospinal fluid dynamics in the lateral ventricle: a feasibility study

PURPOSE

Intravoxel incoherent motion (IVIM) in diffusion-weighted magnetic resonance imaging (DW-MRI) attributes the signal attenuation to the molecular diffusion and to a faster pseudo-diffusion. Purpose of the study was to demonstrate the feasibility of IVIM for the investigation of intracranial cerebrospinal fluid (CSF) dynamics.

METHODS

Cardiac-gated DW-MRI images with fifteen b-values (0-1300s/mm²) along three orthogonal directions (mediolateral (ML), anteroposterior (AP), and craniocaudal (CC)) were acquired during maximum systole and diastole in 10 healthy volunteers (6 males, mean age 36 ± 15 years). A pixel-wise bi-exponential fitting with an iterative nonparametric algorithm was carried out to calculate the following parameters: diffusion coefficient (D), fast diffusion coefficient (D*), and fraction of fast diffusion (f). Region of interest measurements were performed in both lateral ventricles. Comparison of IVIM parameters was performed among two cardiac cycle acquisitions and among the diffusion-encoding directions using a paired Student's t test.

RESULTS

f significantly (p < 0.05) depended on the diffusion-encoding direction and on the cardiac cycle (diastole AP 0.30 ± 0.13, ML 0.22 ± 0.12, CC 0.26 ± 0.17; systole AP 0.45 ± 0.17, ML 0.34 ± 0.15, CC 0.40 ± 0.21). Neither a cardiac cycle nor a direction dependency was found among mean D values (which is in line with the expected intraventricular isotropic diffusion) and D* values (p > 0.05 each).

CONCLUSION

The fraction of fast diffusion from IVIM is feasible to detect a direction-dependent and cardiac-dependent pulsatile CSF flow within the lateral ventricles allowing for quantitative monitoring of CSF dynamics. This technique might provide opportunities to further investigate the pathophysiology of various neurological disorders involving altered CSF dynamics.

Assessment of Cervical Cancer with a Parameter-Free Intravoxel Incoherent Motion Imaging Algorithm

Objective

To evaluate the feasibility of a parameter-free intravoxel incoherent motion (IVIM) approach in cervical cancer, to assess the optimal b-value threshold, and to preliminarily examine differences in the derived perfusion and diffusion parameters for different histological cancer types.

Materials and Methods

After Institutional Review Board approval, 19 female patients (mean age, 54 years; age range, 37–78 years) gave consent and were enrolled in this prospective magnetic resonance imaging study. Clinical staging and biopsy results were obtained. Echo-planar diffusion weighted sequences at 13 b-values were acquired at 3 tesla field strength. Single-sliced region-of-interest IVIM analysis with adaptive b-value thresholds was applied to each tumor, yielding the optimal fit and the optimal parameters for pseudodiffusion (D^*), perfusion fraction (F_p) and diffusion coefficient (D). Monoexponential apparent diffusion coefficient (ADC) was calculated for comparison with D.

Results

Biopsy revealed squamous cell carcinoma in 10 patients and adenocarcinoma in 9. The b-value threshold (median [interquartile range]) depended on the histological type and was 35 (22.5–50) s/mm² in squamous cell carcinoma and 150 (100–150) s/mm² in adenocarcinoma (p < 0.05). Comparing squamous cell vs. adenocarcinoma, D^* (45.1 [25.1–60.4] × 10⁻³ mm²/s vs. 12.4 [10.5–21.2] × 10⁻³ mm²/s) and F_p (7.5% [7.0–9.0%] vs. 9.9% [9.0–11.4%]) differed significantly between the subtypes (p < 0.02), whereas D did not (0.89 [0.75–0.94] × 10⁻³ mm²/s vs. 0.90 [0.82–0.97] × 10⁻³ mm²/s, p = 0.27). The residuals did not differ (0.74 [0.60–0.92] vs. 0.94 [0.67–1.01], p = 0.32). The ADC systematically underestimated the magnitude of diffusion restriction compared to D (p < 0.001).

Conclusion

The parameter-free IVIM approach is feasible in cervical cancer. The b-value threshold and perfusion-related parameters depend on the tumor histology type.