Reproducibility of liver ADC measurements using first moment optimized diffusion imaging

PURPOSE

Cardiac-induced liver motion can bias liver ADC measurements and compromise reproducibility. The purpose of this work was to enable motion-robust DWI on multiple MR scanners and assess reproducibility of the resulting liver ADC measurements.

METHODS

First moment-optimized diffusion imaging (MODI) was implemented on three MR scanners with various gradient performances and field strengths. MODI-DWI and conventional Stejskal-Tanner monopolar (MONO) DWI were acquired in eight (N = 8) healthy volunteers on each scanner, and DWI repetitions were combined using three different averaging methods. For each combination of scanner, acquisition, and averaging method, ADC measurements from each liver segment were collected. Systematic differences in ADC values between scanners and methods were assessed with linear mixed effects modeling, and reproducibility was quantified via reproducibility coefficients.

RESULTS

MODI reduced left-right liver lobe ADC bias from 0.43 × 10-3 mm2/s (MONO) to 0.19 × 10-3 mm2/s (MODI) when simple (unweighted) repetition averaging was used. The bias was reduced from 0.23 × 10-3 mm2/s to 0.06 × 10-3 mm2/s using weighted averaging, and 0.14 × 10-3 mm2/s to 0.01 × 10-3 mm2/s using squared weighted averaging. There was no significant difference in ADC measurements between field strengths or scanner gradient performance. MODI improved reproducibility coefficients compared to MONO: 0.84 × 10-3 mm2/s vs. 0.63 × 10-3 mm2/s (MODI vs. MONO) for simple averaging, 0.66 × 10-3 mm2/s vs. 0.50 × 10-3 mm2/s for weighted averaging, and 0.61 × 10-3 mm2/s vs. 0.47 × 10-3 mm2/s for squared weighted averaging.

CONCLUSION

The feasibility of motion-robust liver DWI using MODI was demonstrated on multiple MR scanners. MODI improved interlobar agreement and reproducibility of ADC measurements in a healthy cohort.

The effects of maternal flow on placental diffusion-weighted MRI and intravoxel incoherent motion parameters

PURPOSE

To investigate and explain observed features of the placental DWI signal in healthy and compromised pregnancies using a mathematical model of maternal blood flow.

METHODS

Thirteen healthy and nine compromised third trimester pregnancies underwent pulse gradient spin echo DWI MRI, with the results compared to MRI data simulated from a 2D mathematical model of maternal blood flow through the placenta. Both sets of data were fitted to an intravoxel incoherent motion (IVIM) model, and a rebound model (defined within text), which described voxels that did not decay monotonically. Both the in vivo and simulated placentas were split into regions of interest (ROIs) to analyze how the signal varies and how IVIM and rebounding parameters change across the placental width.

RESULTS

There was good agreement between the in vivo MRI data, and the data simulated from the mathematical model. Both sets of data included voxels showing a rebounding signal and voxels showing fast signal decay focused near the maternal side of the placenta. In vivo we found higherf IVIM $$ {f}_{IVIM} $$ in the uterine wall and near the maternal side of the placenta, with the slow diffusion coefficientD $$ D $$ reduced in all ROIs in compromised pregnancy.

CONCLUSION

A simulation based entirely on maternal blood explains key features observed in placental DWI, indicating the importance of maternal blood flow in interpreting placental MRI data, and providing potential new metrics for understanding changes in compromised placentas.

Reproducibility of intravoxel incoherent motion quantification in the liver across field strengths and gradient hardware

PURPOSE

To evaluate reproducibility and interlobar agreement of intravoxel incoherent motion (IVIM) quantification in the liver across field strengths and MR scanners with different gradient hardware.

METHODS

Cramer-Rao lower bound optimization was performed to determine optimized monopolar and motion-robust 2D (b-value and first-order motion moment [M1]) IVIM-DWI acquisitions. Eleven healthy volunteers underwent diffusion MRI of the liver, where each optimized acquisition was obtained five times across three MRI scanners. For each data set, IVIM estimates (diffusion coefficient (D), pseudo-diffusion coefficients (d 1 * $$ {d}_1^{\ast } $$ andd 2 * $$ {d}_2^{\ast } $$ ), blood velocity SDs (Vb1 and Vb2), and perfusion fractions [f1 and f2]) were obtained in the right and left liver lobes using two signal models (pseudo-diffusion and M1-dependent physical) with and without T2 correction (fc1 and fc2) and three fitting techniques (tri-exponential region of interest-based full and segmented fitting and blood velocity SD distribution fitting). Reproducibility and interlobar agreement were compared across methods using within-subject and pairwise coefficients of variation (CVw and CVp), paired sample t-tests, and Bland-Altman analysis.

RESULTS

Using a combination of motion-robust 2D (b-M1) data acquisition, M1-dependent physical signal modeling with T2 correction, and blood velocity SD distribution fitting, multiscanner reproducibility with median CVw = 5.09%, 11.3%, 9.20%, 14.2%, and 12.6% for D, Vb1, Vb2, fc1, and fc2, respectively, and interlobar agreement with CVp = 8.14%, 11.9%, 8.50%, 49.9%, and 42.0%, respectively, was achieved.

CONCLUSION

Recently proposed advanced IVIM acquisition, signal modeling, and fitting techniques may facilitate reproducible IVIM quantification in the liver, as needed for establishment of IVIM-based quantitative biomarkers for detection, staging, and treatment monitoring of diseases.

Signal drift in diffusion MRI of the brain: effects on intravoxel incoherent motion parameter estimates

OBJECTIVES

Signal drift has been put forward as one of the fundamental confounding factors in diffusion MRI (dMRI) of the brain. This study characterizes signal drift in dMRI of the brain, evaluates correction methods, and exemplifies its impact on parameter estimation for three intravoxel incoherent motion (IVIM) protocols.

MATERIALS AND METHODS

dMRI of the brain was acquired in ten healthy subjects using protocols designed to enable retrospective characterization and correction of signal drift. All scans were acquired twice for repeatability analysis. Three temporal polynomial correction methods were evaluated: (1) global, (2) voxelwise, and (3) spatiotemporal. Effects of acquisition order were simulated using estimated drift fields.

RESULTS

Signal drift was around 2% per 5 min in the brain as a whole, but reached above 5% per 5 min in the frontal regions. Only correction methods taking spatially varying signal drift into account could achieve effective corrections. Altered acquisition order introduced both systematic changes and differences in repeatability in the presence of signal drift.

DISCUSSION

Signal drift in dMRI of the brain was found to be spatially varying, calling for correction methods taking this into account. Without proper corrections, choice of protocol can affect dMRI parameter estimates and their repeatability.

Pulsation artifact reduction using compressed sensitivity encoding in Gd-EOB-DTPA contrast-enhanced magnetic resonance imaging

In Gd-EOB-DTPA-enhanced MRI, cardiac pulsation artifacts in the left lobe often hinder diagnosis, the image quality need to improve. This study aimed to reduce cardiac pulsation artifacts in Gd-EOB-DTPA-enhanced three-dimensional (3D) T1-weighted turbo-field echo (3D-T1TFE) using compressed sensitivity encoding (CS).For phantom evaluation, the cardiac phantom was manually operated using a metronome-synchronized apparatus, comprising a bag-valve mask, a breathing circuit, and a Jackson-Rees system. Transverse images of a liver phantom were acquired using enhanced T1 high-resolution isotropic volumetric excitation with CS (CS-eTHRIVE) and sensitivity encoding (S-eTHRIVE). For evaluation, images obtained during cardiac phantom operation were subtracted from those obtained when the phantom was stationary. Standard deviation (SD) of the difference images was used as the evaluation metric, and assessments were conducted based on changes in heart rate and TFE factor. For clinical image evaluation, artifacts in hepatobiliary phase images acquired 15 min after Gd-EOB-DTPA injection in the order of S-eTHRIVE and CS-eTHRIVE were visually evaluated at four levels. In heart-rate evaluation (40, 60, and 80 beats/min), CS-eTHRIVE revealed significantly lower SD values compared to S-eTHRIVE across all heart rates (P < 0.01), with no significant differences between heart rates. For TFE factor evaluation, CS-eTHRIVE with a factor of 35 exhibited the lowest SD, which was significantly different from all other groups (P < 0.01). In clinical image evaluation, CS-eTHRIVE achieved higher visual scores (mean ± SD: 3.72 ± 0.46) compared with S-eTHRIVE (2.72 ± 0.98, P < 0.01).CS reduced pulsation artifacts in Gd-EOB-DTPA-enhanced 3D-T1TFE.

Advanced Diffusion-Weighted MRI for Cancer Microstructure Assessment in Body Imaging, and Its Relationship With Histology

Diffusion-weighted magnetic resonance imaging (DW-MRI) aims to disentangle multiple biological signal sources in each imaging voxel, enabling the computation of innovative maps of tissue microstructure. DW-MRI model development has been dominated by brain applications. More recently, advanced methods with high fidelity to histology are gaining momentum in other contexts, for example, in oncological applications of body imaging, where new biomarkers are urgently needed. The objective of this article is to review the state-of-the-art of DW-MRI in body imaging (ie, not including the nervous system) in oncology, and to analyze its value as compared to reference colocalized histology measurements, given that demonstrating the histological validity of any new DW-MRI method is essential. In this article, we review the current landscape of DW-MRI techniques that extend standard apparent diffusion coefficient (ADC), describing their acquisition protocols, signal models, fitting settings, microstructural parameters, and relationship with histology. Preclinical, clinical, and in/ex vivo studies were included. The most used techniques were intravoxel incoherent motion (IVIM; 36.3% of used techniques), diffusion kurtosis imaging (DKI; 16.7%), vascular, extracellular, and restricted diffusion for cytometry in tumors (VERDICT; 13.3%), and imaging microstructural parameters using limited spectrally edited diffusion (IMPULSED; 11.7%). Another notable category of techniques relates to innovative b-tensor diffusion encoding or joint diffusion-relaxometry. The reviewed approaches provide histologically meaningful indices of cancer microstructure (eg, vascularization/cellularity) which, while not necessarily accurate numerically, may still provide useful sensitivity to microscopic pathological processes. Future work of the community should focus on improving the inter-/intra-scanner robustness, and on assessing histological validity in broader contexts. LEVEL OF EVIDENCE: NA TECHNICAL EFFICACY: Stage 2.

Probing Renal Microstructure and Function with Advanced Diffusion MRI: Concepts, Applications, Challenges, and Future Directions

Diffusion measurements in the kidney are affected not only by renal microstructure but also by physiological processes (i.e., glomerular filtration, water reabsorption, and urine formation). Because of the superposition of passive tissue diffusion, blood perfusion, and tubular pre-urine flow, the limitations of the monoexponential apparent diffusion coefficient (ADC) model in assessing pathophysiological changes in renal tissue are becoming apparent and motivate the development of more advanced diffusion-weighted imaging (DWI) variants. These approaches take advantage of the fact that the length scale probed in DWI measurements can be adjusted by experimental parameters, including diffusion-weighting, diffusion gradient directions and diffusion time. This forms the basis by which advanced DWI models can be used to capture not only passive diffusion effects, but also microcirculation, compartmentalization, tissue anisotropy. In this review, we provide a comprehensive overview of the recent advancements in the field of renal DWI. Following a short introduction on renal structure and physiology, we present the key methodological approaches for the acquisition and analysis of renal DWI data, including intravoxel incoherent motion (IVIM), diffusion tensor imaging (DTI), non-Gaussian diffusion, and hybrid IVIM-DTI. We then briefly summarize the applications of these methods in chronic kidney disease and renal allograft dysfunction. Finally, we discuss the challenges and potential avenues for further development of renal DWI. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY: Stage 2.

Calculation of intravoxel incoherent motion parameter maps using a kernelized total difference-based method

Quantitative analysis of diffusion-weighted magnetic resonance imaging (DW-MRI) has been explored for many clinical applications since its development. In particular, the intravoxel incoherent motion (IVIM) model for DW-MRI has been commonly utilized in various organs. However, because of the presence of excessive noise, the IVIM parameter maps obtained from pixel-wise fitting are often unreliable. In this study, we propose a kernelized total difference-based curve-fitting method to estimate the IVIM parameters. Simulated DW-MRI data at five signal-to-noise ratios (i.e., 10, 20, 30, 50, and 100) and real abdominal DW-MRI data acquired on a 1.5-T MRI scanner with nine b-values (i.e., 0, 10, 25, 50, 100, 200, 300, 400, and 500 s/mm2) and six diffusion-encoding gradient directions were used to evaluate the performance of the proposed method. The results were compared with those obtained by three existing methods: trust-region reflective (TRR) algorithm, Bayesian probability (BP), and deep neural network (DNN). Our simulation results showed that the proposed method outperformed the other three comparing methods in terms of root-mean-square error. Moreover, the proposed method could preserve small details in the estimated IVIM parameter maps. The experimental results showed that, compared with the TRR method, the proposed method as well as the BP (and DNN) method could reduce the overestimation of the pseudodiffusion coefficient and improve the quality of IVIM parameter maps. For all studied abdominal organs except the pancreas, both the proposed method and the BP method could provide IVIM parameter estimates close to the reference values; the former had higher precision. The kernelized total difference-based curve-fitting method has the potential to improve the reliability of IVIM parametric imaging.

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

BACKGROUND

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.

PURPOSE

To minimize variability and maximize repeatably in abdominal DWI in terms of IVIM parameter estimates.

STUDY TYPE

Prospective test-retest and image quality comparison.

SUBJECTS

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.

FIELD STRENGTH/SEQUENCE

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.

ASSESSMENT

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.

STATISTICAL TESTS

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.

RESULTS

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).

DATA CONCLUSION

Simultaneous multislice acquisitions had significantly less variability and higher ICCs of Ds, higher eSNR, less distortion, and reduced scan time compared to EPI.

EVIDENCE LEVEL

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

OBJECTIVES

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.

MATERIALS AND METHODS

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.

RESULTS

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).

CONCLUSION

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.

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.

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.

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 parametersD $$ D $$ andf $$ f $$ , and for the Bayesian algorithm only forv 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.

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.

Feature-guided deep learning reduces signal loss and increases lesion CNR in diffusion-weighted imaging of the liver

PURPOSE

This research aims to develop a feature-guided deep learning approach and compare it with an optimized conventional post-processing algorithm in order to enhance the image quality of diffusion-weighted liver images and, in particular, to reduce the pulsation-induced signal loss occurring predominantly in the left liver lobe.

METHODS

Data from 40 patients with liver lesions were used. For the conventional approach, the best-suited out of five examined algorithms was chosen. For the deep learning approach, a U-Net was trained. Instead of learning "gold-standard" target images, the network was trained to optimize four image features (lesion CNR, vessel darkness, data consistency, and pulsation artifact reduction), which could be assessed quantitatively using manually drawn ROIs. A quality score was calculated from these four features. As an additional quality assessment, three radiologists rated different features of the resulting images.

RESULTS

The conventional approach could substantially increase the lesion CNR and reduce the pulsation-induced signal loss. However, the vessel darkness was reduced. The deep learning approach increased the lesion CNR and reduced the signal loss to a slightly lower extent, but it could additionally increase the vessel darkness. According to the image quality score, the quality of the deep-learning images was higher than that of the images obtained using the conventional approach. The radiologist ratings were mostly consistent with the quantitative scores, but the overall quality ratings differed among the readers.

CONCLUSION

Unlike the conventional algorithm, the deep-learning algorithm increased the vessel darkness. Therefore, it may be a viable alternative to conventional algorithms.

In vivo disentanglement of diffusion frequency-dependence, tensor shape, and relaxation using multidimensional MRI

Diffusion MRI with free gradient waveforms, combined with simultaneous relaxation encoding, referred to as multidimensional MRI (MD-MRI), offers microstructural specificity in complex biological tissue. This approach delivers intravoxel information about the microstructure, local chemical composition, and importantly, how these properties are coupled within heterogeneous tissue containing multiple microenvironments. Recent theoretical advances incorporated diffusion time dependency and integrated MD-MRI with concepts from oscillating gradients. This framework probes the diffusion frequency,ω $$ \omega $$ , in addition to the diffusion tensor,D $$ \mathbf{D} $$ , and relaxation,R 1 $$ {R}_1 $$ ,R 2 $$ {R}_2 $$ , correlations. AD ω - R 1 - R 2 $$ \mathbf{D}\left(\omega \right)-{R}_1-{R}_2 $$ clinical imaging protocol was then introduced, with limited brain coverage and 3 mm3 voxel size, which hinder brain segmentation and future cohort studies. In this study, we introduce an efficient, sparse in vivo MD-MRI acquisition protocol providing whole brain coverage at 2 mm3 voxel size. We demonstrate its feasibility and robustness using a well-defined phantom and repeated scans of five healthy individuals. Additionally, we test different denoising strategies to address the sparse nature of this protocol, and show that efficient MD-MRI encoding design demands a nuanced denoising approach. The MD-MRI framework provides rich information that allows resolving the diffusion frequency dependence into intravoxel components based on theirD ω - R 1 - R 2 $$ \mathbf{D}\left(\omega \right)-{R}_1-{R}_2 $$ distribution, enabling the creation of microstructure-specific maps in the human brain. Our results encourage the broader adoption and use of this new imaging approach for characterizing healthy and pathological tissues.

Image denoising and model-independent parameterization for IVIM MRI

Objective. To improve intravoxel incoherent motion imaging (IVIM) magnetic resonance Imaging quality using a new image denoising technique and model-independent parameterization of the signal versusb-value curve.Approach. IVIM images were acquired for 13 head-and-neck patients prior to radiotherapy. Post-radiotherapy scans were also acquired for five of these patients. Images were denoised prior to parameter fitting using neural blind deconvolution, a method of solving the ill-posed mathematical problem of blind deconvolution using neural networks. The signal decay curve was then quantified in terms of several area under the curve (AUC) parameters. Improvements in image quality were assessed using blind image quality metrics, total variation (TV), and the correlations between parameter changes in parotid glands with radiotherapy dose levels. The validity of blur kernel predictions was assessed by the testing the method's ability to recover artificial 'pseudokernels'. AUC parameters were compared with monoexponential, biexponential, and triexponential model parameters in terms of their correlations with dose, contrast-to-noise (CNR) around parotid glands, and relative importance via principal component analysis.Main results. Image denoising improved blind image quality metrics, smoothed the signal versusb-value curve, and strengthened correlations between IVIM parameters and dose levels. Image TV was reduced and parameter CNRs generally increased following denoising.AUCparameters were more correlated with dose and had higher relative importance than exponential model parameters.Significance. IVIM parameters have high variability in the literature and perfusion-related parameters are difficult to interpret. Describing the signal versusb-value curve with model-independent parameters like theAUCand preprocessing images with denoising techniques could potentially benefit IVIM image parameterization in terms of reproducibility and functional utility.

Facilitating diffusion tensor imaging of the brain during continuous gross head motion with first and second order motion compensating diffusion gradients

PURPOSE

To demonstrate the feasibility of motion compensating diffusion gradient schemes in the acquisition of quality diffusion tensor images (DTI) of the brain during continuous gross head motion.

METHODS

Five healthy subjects were scanned using a clinical 3 T MRI with and without continuous head motion. For one volunteer, DTI data was acquired using standard (M0) diffusion-weighted (DW) gradients, and first (M1) and second (M2) order gradient schemes that were previously developed for use in cardiac DTI. In four additional volunteers, DTI data was acquired with M0 and M2 gradients. DTI parameters were calculated and compared with established retrospective motion corrections.

RESULTS

In the absence of motion, DTI parameters calculated from M0, M1, and M2 data were consistent. In the presence of motion, up to 44% of DW images acquired with M0 gradients were corrupted by signal dropout, compared to 0% of the M2 images. In voxelwise comparisons, DTI parameters calculated using motion-M0 data were elevated compared to reference data. Retrospective corrections for extreme motion applied to motion-M0 data did not improve consistency with reference data in cases where motion corrupted >15% of DW images. In contrast, DTI parameters calculated with motion-M2 data were consistent with reference data.

CONCLUSION

This proof-of-principle study demonstrates that motion compensating diffusion gradients can mitigate artifacts because of continuous motion in DTI of the brain and offers promise for improved DTI accessibility. Further study will be necessary to determine the robustness of the approach in patient populations with high susceptibility to head motion.

In vivo disentanglement of diffusion frequency-dependence, tensor shape, and relaxation using multidimensional MRI

Diffusion MRI with free gradient waveforms, combined with simultaneous relaxation encoding, referred to as multidimensional MRI (MD-MRI), offers microstructural specificity in complex biological tissue. This approach delivers intravoxel information about the microstructure, local chemical composition, and importantly, how these properties are coupled within heterogeneous tissue containing multiple microenvironments. Recent theoretical advances incorporated diffusion time dependency and integrated MD-MRI with concepts from oscillating gradients. This framework probes the diffusion frequency, ω , in addition to the diffusion tensor, D , and relaxation, R 1 , R 2 , correlations. A D ( ω ) - R 1 - R 2 clinical imaging protocol was then introduced, with limited brain coverage and 3 mm3 voxel size, which hinder brain segmentation and future cohort studies. In this study, we introduce an efficient, sparse in vivo MD-MRI acquisition protocol providing whole brain coverage at 2 mm3 voxel size. We demonstrate its feasibility and robustness using a well-defined phantom and repeated scans of five healthy individuals. Additionally, we test different denoising strategies to address the sparse nature of this protocol, and show that efficient MD-MRI encoding design demands a nuanced denoising approach. The MD-MRI framework provides rich information that allows resolving the diffusion frequency dependence into intravoxel components based on their D ( ω ) - R 1 - R 2 distribution, enabling the creation of microstructure-specific maps in the human brain. Our results encourage the broader adoption and use of this new imaging approach for characterizing healthy and pathological tissues.

Impact of velocity- and acceleration-compensated encodings on signal dropout and black-blood state in diffusion-weighted magnetic resonance liver imaging at clinical TEs

PURPOSE

The study aims to develop easy-to-implement concomitant field-compensated gradient waveforms with varying velocity-weighting (M1) and acceleration-weighting (M2) levels and to evaluate their efficacy in correcting signal dropouts and preserving the black-blood state in liver diffusion-weighted imaging. Additionally, we seek to determine an optimal degree of compensation that minimizes signal dropouts while maintaining blood signal suppression.

METHODS

Numerically optimized gradient waveforms were adapted using a novel method that allows for the simultaneous tuning of M1- and M2-weighting by changing only one timing variable. Seven healthy volunteers underwent diffusion-weighted magnetic resonance imaging (DWI) with five diffusion encoding schemes (monopolar, velocity-compensated (M1 = 0), acceleration-compensated (M1 = M2 = 0), 84%-M1-M2-compensated, 67%-M1-M2-compensated) at b-values of 50 and 800 s/mm2 at a constant echo time of 70 ms. Signal dropout correction and apparent diffusion coefficients (ADCs) were quantified using regions of interest in the left and right liver lobe. The blood appearance was evaluated using two five-point Likert scales.

RESULTS

Signal dropout was more pronounced in the left lobe (19%-42% less signal than in the right lobe with monopolar scheme) and best corrected by acceleration-compensation (8%-10% less signal than in the right lobe). The black-blood state was best with monopolar encodings and decreased significantly (p < 0.001) with velocity- and/or acceleration-compensation. The partially M1-M2-compensated encoding schemes could restore the black-blood state again. Strongest ADC bias occurred for monopolar encodings (difference between left/right lobe of 0.41 μm2/ms for monopolar vs. < 0.12 μm2/ms for the other encodings).

CONCLUSION

All of the diffusion encodings used in this study demonstrated suitability for routine DWI application. The results indicate that a perfect value for the level of M1-M2-compensation does not exist. However, among the examined encodings, the 84%-M1-M2-compensated encodings provided a suitable tradeoff.

Realistic numerical simulations of diffusion tensor cardiovascular magnetic resonance: The effects of perfusion and membrane permeability

PURPOSE

To study the sensitivity of diffusion tensor cardiovascular magnetic resonance (DT-CMR) to microvascular perfusion and changes in cell permeability.

METHODS

Monte Carlo (MC) random walk simulations in the myocardium have been performed to simulate self-diffusion of water molecules in histology-based media with varying extracellular volume fraction (ECV) and permeable membranes. The effect of microvascular perfusion on simulations of the DT-CMR signal has been incorporated by adding the contribution of particles traveling through an anisotropic capillary network to the diffusion signal. The simulations have been performed considering three pulse sequences with clinical gradient strengths: monopolar stimulated echo acquisition mode (STEAM), monopolar pulsed-gradient spin echo (PGSE), and second-order motion-compensated spin echo (MCSE).

RESULTS

Reducing ECV intensifies the diffusion restriction and incorporating membrane permeability reduces the anisotropy of the diffusion tensor. Widening the intercapillary velocity distribution results in increased measured diffusion along the cardiomyocytes long axis when the capillary networks are anisotropic. Perfusion amplifies the mean diffusivity for STEAM while the opposite is observed for short diffusion encoding time sequences (PGSE and MCSE).

CONCLUSION

The effect of perfusion on the measured diffusion tensor is reduced using an increased reference b-value. Our results pave the way for characterization of the response of DT-CMR to microstructural changes underlying cardiac pathology and highlight the higher sensitivity of STEAM to permeability and microvascular circulation due to its longer diffusion encoding time.

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.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

An unsupervised convolutional neural network method for estimation of intravoxel incoherent motion parameters

Objective. Intravoxel incoherent motion (IVIM) imaging obtained by fitting a biexponential model to multiple b-value diffusion-weighted magnetic resonance imaging (DW-MRI) has been shown to be a promising tool for different clinical applications. Recently, several deep neural network (DNN) methods were proposed to generate IVIM imaging. Approach. In this study, we proposed an unsupervised convolutional neural network (CNN) method for estimation of IVIM parameters. We used both simulated and real abdominal DW-MRI data to evaluate the performance of the proposed CNN-based method, and compared the results with those obtained from a non-linear least-squares fit (TRR, trust-region reflective algorithm) and a feed-forward backward-propagation DNN-based method. Main results. The simulation results showed that both the DNN- and CNN-based methods had lower coefficients of variation than the TRR method, but the CNN-based method provided more accurate parameter estimates. The results obtained from real DW-MRI data showed that the TRR method produced many biased IVIM parameter estimates that hit the upper and lower parameter bounds. In contrast, both the DNN- and CNN-based methods yielded less biased IVIM parameter estimates. Overall, the perfusion fraction and diffusion coefficient obtained from the DNN- and CNN-based methods were close to literature values. However, compared with the CNN-based method, both the TRR and DNN-based methods tended to yield increased pseudodiffusion coefficients (55%–180%). Significance. Our preliminary results suggest that it is feasible to estimate IVIM parameters using CNN.

Reduction of the cardiac pulsation artifact and improvement of lesion conspicuity in flow‐compensated diffusion images in the liver—A quantitative evaluation of postprocessing algorithms

PURPOSE

To enhance image quality of flow-compensated diffusion-weighted liver MRI data by increasing the lesion conspicuity and reducing the cardiac pulsation artifact using postprocessing algorithms.

METHODS

Diffusion-weighted image data of 40 patients with liver lesions had been acquired at 1.5 T. These data were postprocessed with 5 different algorithms (weighted averaging, p-mean, percentile, outlier exclusion, and exception set). Four image properties of the postprocessed data were evaluated for optimizing the algorithm parameters. These properties were the lesion to tissue contrast-to-noise ratio (CNR), the reduction of the cardiac pulsation artifact, the data consistency, and the vessel darkness. They were combined into a total quality score ( Q total , $$ {Q}_{\mathrm{total}}, $$ set to 1 for the trace-weighted reference image), which was used to rate the image quality objectively.

RESULTS

The weighted averaging algorithm performed best according to the total quality score ( Q total = 1.111 ± 0.067 $$ {Q}_{\mathrm{total}}=1.111\pm 0.067 $$ ). The further ranking was outlier exclusion algorithm ( Q total = 1.086 ± 0.061 $$ {Q}_{\mathrm{total}}=1.086\pm 0.061 $$ ), p-mean algorithm ( Q total = 1.045 ± 0.049 $$ {Q}_{\mathrm{total}}=1.045\pm 0.049 $$ ), percentile algorithm ( Q total = 1.012 ± 0.049 $$ {Q}_{\mathrm{total}}=1.012\pm 0.049 $$ ), and exception set algorithm ( Q total = 0.957 ± 0.027 $$ {Q}_{\mathrm{total}}=0.957\pm 0.027 $$ ). All optimized algorithms except for the exception set algorithm corrected the pulsation artifact and increased the lesion CNR. Changes in Q total $$ {Q}_{\mathrm{total}} $$ were significant for all optimized algorithms except for the percentile algorithm. Liver ADC was significantly reduced (except for the exception set algorithm), particularly in the left lobe.

CONCLUSION

Postprocessing algorithms should be used for flow-compensated liver DWI. The proposed weighted averaging algorithm seems to be suited best to increase the image quality of artifact-corrupted flow-compensated diffusion-weighted liver data.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Flow-compensated diffusion encoding in MRI for improved liver metastasis detection

Magnetic resonance (MR) diffusion-weighted imaging (DWI) is often used to detect focal liver lesions (FLLs), though DWI image quality can be limited in the left liver lobe owing to the pulsatile motion of the nearby heart. Flow-compensated (FloCo) diffusion encoding has been shown to reduce this pulsation artifact. The purpose of this prospective study was to intra-individually compare DWI of the liver acquired with conventional monopolar and FloCo diffusion encoding for assessing metastatic FLLs in non-cirrhotic patients. Forty patients with known or suspected multiple metastatic FLLs were included and measured at 1.5 T field strength with a conventional (monopolar) and a FloCo diffusion encoding EPI sequence (single refocused; b-values, 50 and 800 s/mm²). Two board-certified radiologists analyzed the DWI images independently. They issued Likert-scale ratings (1 = worst, 5 = best) for pulsation artifact severity and counted the difference of lesions visible at b = 800 s/mm² separately for small and large FLLs (i.e., < 1 cm or > 1 cm) and separately for left and right liver lobe. Differences between the two diffusion encodings were assessed with the Wilcoxon signed-rank test. Both readers found a reduction in pulsation artifact in the liver with FloCo encoding (p < 0.001 for both liver lobes). More small lesions were detected with FloCo diffusion encoding in both liver lobes (left lobe: six and seven additional lesions by readers 1 and 2, respectively; right lobe: five and seven additional lesions for readers 1 and 2, respectively). Both readers found one additional large lesion in the left liver lobe. Thus, flow-compensated diffusion encoding appears more effective than monopolar diffusion encoding for the detection of liver metastases.

Intravoxel Incoherent Motion Magnetic Resonance Imaging in Skeletal Muscle: Review and Future Directions

Throughout the body, muscle structure and function can be interrogated using a variety of noninvasive magnetic resonance imaging (MRI) methods. Recently, intravoxel incoherent motion (IVIM) MRI has gained momentum as a method to evaluate components of blood flow and tissue diffusion simultaneously. Much of the prior research has focused on highly vascularized organs, including the brain, kidney, and liver. Unique aspects of skeletal muscle, including the relatively low perfusion at rest and large dynamic range of perfusion between resting and maximal hyperemic states, may influence the acquisition, postprocessing, and interpretation of IVIM data. Here, we introduce several of those unique features of skeletal muscle; review existing studies of IVIM in skeletal muscle at rest, in response to exercise, and in disease states; and consider possible confounds that should be addressed for muscle-specific evaluations. Most studies used segmented nonlinear least squares fitting with a b-value threshold of 200 sec/mm2 to obtain IVIM parameters of perfusion fraction (f), pseudo-diffusion coefficient (D*), and diffusion coefficient (D). In healthy individuals, across all muscles, the average ± standard deviation of D was 1.46 ± 0.30 × 10-3  mm2 /sec, D* was 29.7 ± 38.1 × 10-3  mm2 /sec, and f was 11.1 ± 6.7%. Comparisons of reported IVIM parameters in muscles of the back, thigh, and leg of healthy individuals showed no significant difference between anatomic locations. Throughout the body, exercise elicited a positive change of all IVIM parameters. Future directions including advanced postprocessing models and potential sequence modifications are discussed. LEVEL OF EVIDENCE: 2 TECHNICAL EFFICACY: Stage 2.

Nonparametric D-R1-R2 distribution MRI of the living human brain

Diffusion-relaxation correlation NMR can simultaneously characterize both the microstructure and the local chemical composition of complex samples that contain multiple populations of water. Recent developments on tensor-valued diffusion encoding and Monte Carlo inversion algorithms have made it possible to transfer diffusion-relaxation correlation NMR from small-bore scanners to clinical MRI systems. Initial studies on clinical MRI systems employed 5D D-R₁ and D-R₂ correlation to characterize healthy brain in vivo. However, these methods are subject to an inherent bias that originates from not including R₂ or R₁ in the analysis, respectively. This drawback can be remedied by extending the concept to 6D D-R₁-R₂ correlation. In this work, we present a sparse acquisition protocol that records all data necessary for in vivo 6D D-R₁-R₂ correlation MRI across 633 individual measurements within 25 min-a time frame comparable to previous lower-dimensional acquisition protocols. The data were processed with a Monte Carlo inversion algorithm to obtain nonparametric 6D D-R₁-R₂ distributions. We validated the reproducibility of the method in repeated measurements of healthy volunteers. For a post-therapy glioblastoma case featuring cysts, edema, and partially necrotic remains of tumor, we present representative single-voxel 6D distributions, parameter maps, and artificial contrasts over a wide range of diffusion-, R₁-, and R₂-weightings based on the rich information contained in the D-R₁-R₂ distributions.

Echo time dependence of biexponential and triexponential intravoxel incoherent motion parameters in the liver

PURPOSE

Intravoxel incoherent motion (IVIM) studies are performed with different acquisition protocols. Comparing them requires knowledge of echo time (TE) dependencies. The TE-dependence of the biexponential perfusion fraction f is well-documented, unlike that of its triexponential counterparts f₁ and f₂ and the biexponential and triexponential pseudodiffusion coefficients D^* , D 1 ∗ , and D 2 ∗ . The purpose was to investigate the TE-dependence of these parameters and to check whether the triexponential pseudodiffusion compartments are associated with arterial and venous blood.

METHODS

Fifteen healthy volunteers (19-58 y; mean: 24.7 y) underwent diffusion-weighted imaging of the abdomen with 24 b-values (0.2-800 s/mm² ) at TEs of 45, 60, 75, and 90 ms. Regions of interest (ROIs) were manually drawn in the liver. One set of bi- and triexponential IVIM parameters per volunteer and TE was determined. The TE-dependence was assessed with the Kruskal-Wallis test.

RESULTS

TE-dependence was observed for f (P < .001), f₁ (P = .001), and f₂ (P < .001). Their median values at the four measured TEs were: f: 0.198/0.240/0.274/0.359, f₁ : 0.113/0.139/0.146/0.205, f₂ : 0.115/0.155/0.182/0.194. D, D^* , D 1 ∗ , and D 2 ∗ showed no significant TE-dependence. Their values were: diffusion coefficient D (10^-4 mm² /s): 9.45/9.63/9.75/9.41, biexponential D^* (10^-2 mm² /s): 5.26/5.52/6.13/5.82, triexponential D 1 ∗ (10^-2 mm² /s): 1.73/2.91/2.25/2.51, triexponential D 2 ∗ (mm² /s): 0.478/1.385/0.616/0.846.

CONCLUSION

f₁ and f₂ show similar TE-dependence as f, ie, increase with rising TE; an effect that must be accounted for when comparing different studies. The diffusion and pseudodiffusion coefficients might be compared without TE correction. Because of the similar TE-dependence of f₁ and f₂ , the triexponential pseudodiffusion compartments are most probably not associated to venous and arterial blood.

Characterisation of microvessel blood velocity and segment length in the brain using multi-diffusion-time diffusion-weighted MRI

Multi-diffusion-time diffusion-weighted MRI can probe tissue microstructure, but the method has not been widely applied to the microvasculature. At long diffusion-times, blood flow in capillaries is in the diffusive regime, and signal attenuation is dependent on blood velocity (v) and capillary segment length (l). It is described by the pseudo-diffusion coefficient (D*=vl/6) of intravoxel incoherent motion (IVIM). At shorter diffusion-times, blood flow is in the ballistic regime, and signal attenuation depends on v, and not l. In theory, l could be estimated using D* and v. In this study, we compare the accuracy and repeatability of three approaches to estimating v, and therefore l: the IVIM ballistic model, the velocity autocorrelation model, and the ballistic approximation to the velocity autocorrelation model. Twenty-nine rat datasets from two strains were acquired at 7 T, with b-values between 0 and 1000 smm^-2 and diffusion times between 11.6 and 50 ms. Five rats were scanned twice to assess scan-rescan repeatability. Measurements of l were validated using corrosion casting and micro-CT imaging. The ballistic approximation of the velocity autocorrelation model had lowest bias relative to corrosion cast estimates of l, and had highest repeatability.

Short-Term Repeatability of in Vivo Cardiac Intravoxel Incoherent Motion Tensor Imaging in Healthy Human Volunteers

BACKGROUND

Intravoxel incoherent motion (IVIM) tensor imaging is a promising technique for diagnosis and monitoring of cardiovascular diseases. Knowledge about measurement repeatability, however, remains limited.

PURPOSE

To evaluate short-term repeatability of IVIM tensor imaging in normal in vivo human hearts.

STUDY TYPE

Prospective.

POPULATION

Ten healthy subjects without history of heart diseases.

FIELD STRENGTH/SEQUENCE

Balanced steady-state free-precession cine sequence and single-shot spin-echo echo planar IVIM tensor imaging sequence (9 b-values, 0-400 seconds/mm² and six diffusion-encoding directions) at 3.0 T.

ASSESSMENT

Subjects were scanned twice with an interval of 15 minutes, leaving the scanner between studies. The signal-to-noise ratio (SNR) was evaluated in anterior, lateral, septal, and inferior segments of the left ventricle wall. Fractional anisotropy (FA), mean diffusivity (MD), mean fraction (MF), and helix angle (HA) in the four segments were independently measured by five radiologists.

STATISTICAL TESTS

IVIM tensor indexes were compared between observers using a one-way analysis of variance or between scans using a paired t-test (normal data) or a Wilcoxon rank-sum test (non-normal data). Interobserver agreement and test-retest repeatability were assessed using the intraclass correlation coefficient (ICC), within-subject coefficient of variation (WCV), and Bland-Altman limits of agreements.

RESULTS

SNR of inferior segment was significantly lower than the other three segments, and inferior segment was therefore excluded from repeatability analysis. Interobserver repeatability was excellent for all IVIM tensor indexes (ICC: 0.886-0.972; WCV: 0.62%-4.22%). Test-retest repeatability was excellent for MD of the self-diffusion tensor (D) and MF of the perfusion fraction tensor (f_p ) (ICC: 0.803-0.888; WCV: 1.42%-9.51%) and moderate for FA and MD of the pseudo-diffusion tensor (D^* ) (ICC: 0.487-0.532; WCV: 6.98%-10.89%). FA of D and f_p and HA of D presented good test-retest repeatability (ICC: 0.732-0.788; WCV: 3.28%-8.71%).

DATA CONCLUSION

The D and f_p indexes exhibited satisfactory repeatability, but further efforts were needed to improve repeatability of D^* indexes.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY: Stage 1.

Simulation of intravoxel incoherent perfusion signal using a realistic capillary network of a mouse brain

PURPOSE

To simulate the intravoxel incoherent perfusion magnetic resonance magnitude signal from the motion of blood particles in three realistic vascular network graphs from a mouse brain.

METHODS

In three networks generated from the cortex of a mouse scanned by two-photon laser microscopy, blood flow in each vessel was simulated using Poiseuille's law. The trajectories, flow speeds and phases acquired by a fixed number of simulated blood particles during a Stejskal-Tanner bipolar pulse gradient scheme were computed. The resulting magnitude signal was obtained by integrating all phases and the pseudo-diffusion coefficient D* was estimated by fitting an exponential signal decay. To better understand the anatomical source of the intravoxel incoherent motion (IVIM) perfusion signal, the above was repeated restricting the simulation to various types of vessel.

RESULTS

The characteristics of the three microvascular networks were respectively vessel lengths (mean ± std. dev.) 67.2 ± 53.6 μm, 59.8 ± 46.2 μm and 64.5 ± 50.9 μm, diameters 6.0 ± 3.5 μm, 5.7 ± 3.6 μm and 6.1 ± 3.7 μm and simulated blood velocity 0.9 ± 1.7 μm/ms, 1.4 ± 2.5 μm/ms and 0.7 ± 2.1 μm/ms. Exponential fitting of the simulated signal decay as a function of b-value resulted in the following D*-values [10-3 mm2 /s]: 31.7, 40.4 and 33.4. The signal decay for low b-values was the largest in the larger vessels, but the smaller vessels and the capillaries accounted for more of the total volume of the networks.

CONCLUSION

This simulation improves the theoretical understanding of the IVIM perfusion estimation method by directly linking the MR IVIM perfusion signal to an ultra-high resolution measurement of the microvascular network and a realistic blood flow simulation.

Understanding diffusion-weighted MRI analysis: Repeatability and performance of diffusion models in a benign breast lesion cohort

Diffusion-weighted MRI (DWI) is an important tool for oncology research, with great clinical potential for the classification and monitoring of breast lesions. The utility of parameters derived from DWI, however, is influenced by specific analysis choices. The purpose of this study was to critically evaluate repeatability and curve-fitting performance of common DWI signal representations, for a prospective cohort of patients with benign breast lesions. Twenty informed, consented patients with confirmed benign breast lesions underwent repeated DWI (3 T) using: sagittal single-shot spin-echo echo planar imaging, bipolar encoding, TR/TE: 11,600/86 ms, FOV: 180 x 180 mm, matrix: 90 x 90, slices: 60 x 2.5 mm, iPAT: GRAPPA 2, fat suppression, and 13 b-values: 0-700 s/mm² . A phase-reversed scan (b = 0 s/mm² ) was acquired for distortion correction. Voxel-wise repeat-measures coefficients of variation (CoVs) were derived for monoexponential (apparent diffusion coefficient [ADC]), biexponential (intravoxel incoherent motion: f, D, D*) and stretched exponential (α, DDC) across the parameter histograms for lesion regions of interest (ROIs). Goodness-of-fit for each representation was assessed by Bayesian information criterion. The volume of interest (VOI) definition was repeatable (CoV 13.9%). Within lesions, and across both visits and the cohort, there was no dominant best-fit model, with all representations giving the best fit for a fraction of the voxels. Diffusivity measures from the signal representations (ADC, D, DDC) all showed good repeatability (CoV < 10%), whereas parameters associated with pseudodiffusion (f, D*) performed poorly (CoV > 50%). The stretching exponent α was repeatable (CoV < 12%). This pattern of repeatability was consistent over the central part of the parameter percentiles. Assumptions often made in diffusion studies about analysis choices will influence the detectability of changes, potentially obscuring useful information. No single signal representation prevails within or across lesions, or across repeated visits; parameter robustness is therefore a critical consideration. Our results suggest that stretched exponential representation is more repeatable than biexponential, with pseudodiffusion parameters unlikely to provide clinically useful biomarkers.

Bayesian inference using hierarchical and spatial priors for intravoxel incoherent motion MR imaging in the brain: Analysis of cancer and acute stroke

The intravoxel incoherent motion (IVIM) model allows to map diffusion (D) and perfusion-related parameters (F and D*). Parameter estimation is, however, error-prone due to the non-linearity of the signal model, the limited signal-to-noise ratio (SNR) and the small volume fraction of perfusion in the in-vivo brain. In the present work, the performance of Bayesian inference was examined in the presence of brain pathologies characterized by hypo- and hyperperfusion. In particular, a hierarchical and a spatial prior were combined. Performance was compared relative to conventional segmented least squares regression, hierarchical prior only (non-segmented and segmented data likelihoods) and a deep learning approach. Realistic numerical brain IVIM simulations were conducted to assess errors relative to ground truth. In-vivo, data of 11 central nervous system cancer patients and 9 patients with acute stroke were acquired. The proposed method yielded reduced error in simulations for both the cancer and acute stroke scenarios compared to other methods across the whole investigated SNR range. The contrast-to-noise ratio of the proposed method was better or on par compared to the other techniques in-vivo. The proposed Bayesian approach hence improves IVIM parameter estimation in brain cancer and acute stroke.

Characterization and correction of cardiovascular motion artifacts in diffusion-weighted imaging of the pancreas

PURPOSE

To assess the effects of cardiovascular-induced motion on conventional DWI of the pancreas and to evaluate motion-robust DWI methods in a motion phantom and healthy volunteers.

METHODS

3T DWI was acquired using standard monopolar and motion-compensated gradient waveforms, including in an anatomically accurate pancreas phantom with controllable compressive motion and healthy volunteers (n = 8, 10). In volunteers, highly controlled single-slice DWI using breath-holding and cardiac gating and whole-pancreas respiratory-triggered DWI were acquired. For each acquisition, the ADC variability across volunteers, as well as ADC differences across parts of the pancreas were evaluated.

RESULTS

In motion phantom scans, conventional DWI led to biased ADC, whereas motion-compensated waveforms produced consistent ADC. In the breath-held, cardiac-triggered study, conventional DWI led to heterogeneous DW signals and highly variable ADC across the pancreas, whereas motion-compensated DWI avoided these artifacts. In the respiratory-triggered study, conventional DWI produced heterogeneous ADC across the pancreas (head: 1756 ± 173 × 10^-6 mm² /s; body: 1530 ± 338 × 10^-6 mm² /s; tail: 1388 ± 267 × 10^-6 mm² /s), with ADCs in the head significantly higher than in the tail (P < .05). Motion-compensated ADC had lower variability across volunteers (head: 1277 ± 102 × 10^-6 mm² /s; body: 1204 ± 169 × 10^-6 mm² /s; tail: 1235 ± 178 × 10^-6 mm² /s), with no significant difference (P ≥ .19) across the pancreas.

CONCLUSION

Cardiovascular motion introduces artifacts and ADC bias in pancreas DWI, which are addressed by motion-robust DWI.

Effect of compression garments on muscle perfusion in delayed-onset muscle soreness: A quantitative analysis using intravoxel incoherent motion MR perfusion imaging

The aim of this prospective cohort study was to evaluate the effect of compression garments under resting conditions and after the induction of delayed-onset muscle soreness (DOMS) by MR perfusion imaging using intravoxel incoherent motion (IVIM). Magnetic resonance imaging of both lower legs of 16 volunteers was performed before and after standardized eccentric exercises that induced DOMS. A compression garment (21-22 mmHg) was worn during and for 6 h after exercise on one randomly selected leg. IVIM MR imaging, represented as total muscle perfusion D*f, perfusion fraction f and tissue diffusivity D, were compared between baseline and directly, 30 min, 6 h and 48 h after exhausting exercise with and without compression. Creatine kinase levels and T2-weighted images were acquired at baseline and after 48 h. DOMS was induced in the medial head of the gastrocnemius muscle (MGM) in all volunteers. Compression garments did not show any significant effect on IVIM perfusion parameters at any time point in the MGM or the tibialis anterior muscle (p > 0.05). Microvascular perfusion in the MGM increased significantly in both the compressed and noncompressed leg between baseline measurements and those taken directly after and 30 min after the exercise: the relative median f increased by 31.5% and 24.7% in the compressed and noncompressed leg, respectively, directly after the exercise compared with the baseline value. No significant change in tissue perfusion occurred 48 h after the induction of DOMS compared with baseline. It was concluded that compression garments (21-22 mmHg) do not alter microvascular muscle perfusion at rest, nor do they have any significant effect during the regeneration phase of DOMS. In DOMS, only a short-term effect of increased muscle perfusion (30 min after exercise) was observed, with normalization occurring during regeneration after 6-48 h. The normalization of perfusion independently of compression after 6 h may have implications for diagnostic and therapeutic strategies and for the better understanding of pathophysiological pathways in DOMS.

Non-Invasive Evaluation of Cerebral Microvasculature Using Pre-Clinical MRI: Principles, Advantages and Limitations

Alterations to the cerebral microcirculation have been recognized to play a crucial role in the development of neurodegenerative disorders. However, the exact role of the microvascular alterations in the pathophysiological mechanisms often remains poorly understood. The early detection of changes in microcirculation and cerebral blood flow (CBF) can be used to get a better understanding of underlying disease mechanisms. This could be an important step towards the development of new treatment approaches. Animal models allow for the study of the disease mechanism at several stages of development, before the onset of clinical symptoms, and the verification with invasive imaging techniques. Specifically, pre-clinical magnetic resonance imaging (MRI) is an important tool for the development and validation of MRI sequences under clinically relevant conditions. This article reviews MRI strategies providing indirect non-invasive measurements of microvascular changes in the rodent brain that can be used for early detection and characterization of neurodegenerative disorders. The perfusion MRI techniques: Dynamic Contrast Enhanced (DCE), Dynamic Susceptibility Contrast Enhanced (DSC) and Arterial Spin Labeling (ASL), will be discussed, followed by less established imaging strategies used to analyze the cerebral microcirculation: Intravoxel Incoherent Motion (IVIM), Vascular Space Occupancy (VASO), Steady-State Susceptibility Contrast (SSC), Vessel size imaging, SAGE-based DSC, Phase Contrast Flow (PC) Quantitative Susceptibility Mapping (QSM) and quantitative Blood-Oxygenation-Level-Dependent (qBOLD). We will emphasize the advantages and limitations of each strategy, in particular on applications for high-field MRI in the rodent's brain.

The value of intravoxel incoherent motion imaging in predicting the survival of patients with astrocytoma

BACKGROUND

The evaluation of the prognosis of gliomas may have great value in individualized treatment.

PURPOSE

To evaluate the value of intravoxel incoherent motion (IVIM) in predicting the survival of patients with astrocytoma and comparing it to apparent diffusion coefficients (ADC).

MATERIAL AND METHODS

Sixty patients with pathologically confirmed cerebral astrocytomas underwent IVIM scans before any treatment was performed. Patients were divided into death group and survival group according to a two-year follow-up. ADC and quantitative parameters of IVIM including D, D*, and f were measured. Independent sample t test was used to compare the two groups of parameters. The accuracy of each parameter for two-year survival rate was analyzed by receiver operating characteristic (ROC) curve and Kaplan-Meier survival curves. The correlation between quantitative parameters and survival days was analyzed by Pearson correlation analysis.

RESULTS

The ADC, D*, and f values were statistically significant different between the death and the survival groups (P < 0.05). The AUC of the ADC, D*, and f were 0.811, 0.858, and 0.892, respectively. The ADC cut-off value of 0.668 × 10^-3 mm²/s corresponded to 82.6% sensitivity and 73% specificity. The D* cut-off value of 3.913 × 10^-3 mm²/s corresponded to 78.4% sensitivity and 87% specificity. The f cut-off value of 0.487 corresponded to 83.8% sensitivity and 87% specificity. Significant log rank test was performed for each parameter to predict overall survival (P < 0.05). There was a correlation between ADC (r = 0.625, P = 0.023), D* (r = -0.655, P = 0.012), f (r = -0.725, P = 0.000) and survival days.

CONCLUSION

The D* and f values demonstrated great potential in predicting the two-year survival rate for patients with astrocytoma.

Repeatability of IVIM biomarkers from diffusion-weighted MRI in head and neck: Bayesian probability versus neural network

Purpose

The intravoxel incoherent motion (IVIM) model for DWI might provide useful biomarkers for disease management in head and neck cancer. This study compared the repeatability of three IVIM fitting methods to the conventional nonlinear least‐squares regression: Bayesian probability estimation, a recently introduced neural network approach, IVIM‐NET, and a version of the neural network modified to increase consistency, IVIM‐NET_mod.

Methods

Ten healthy volunteers underwent two imaging sessions of the neck, two weeks apart, with two DWI acquisitions per session. Model parameters (ADC, diffusion coefficient Dt, perfusion fraction fp, and pseudo‐diffusion coefficient Dp) from each fit method were determined in the tonsils and in the pterygoid muscles. Within‐subject coefficients of variation (wCV) were calculated to assess repeatability. Training of the neural network was repeated 100 times with random initialization to investigate consistency, quantified by the coefficient of variance.

Results

The Bayesian and neural network approaches outperformed nonlinear regression in terms of wCV. Intersession wCV of Dt in the tonsils was 23.4% for nonlinear regression, 9.7% for Bayesian estimation, 9.4% for IVIM‐NET, and 11.2% for IVIM‐NET_mod. However, results from repeated training of the neural network on the same data set showed differences in parameter estimates: The coefficient of variances over the 100 repetitions for IVIM‐NET were 15% for both Dt and fp, and 94% for Dp; for IVIM‐NET_mod, these values improved to 5%, 9%, and 62%, respectively.

Conclusion

Repeatabilities from the Bayesian and neural network approaches are superior to that of nonlinear regression for estimating IVIM parameters in the head and neck.

Assessment of tissue perfusion of pancreatic cancer as potential imaging biomarker by means of Intravoxel incoherent motion MRI and CT perfusion: correlation with histological microvessel density as ground truth

Background/objectives

The aim of this study was to compare intravoxel incoherent motion (IVIM) diffusion weighted (DW) MRI and CT perfusion to assess tumor perfusion of pancreatic ductal adenocarcinoma (PDAC).

Methods

In this prospective study, DW-MRI and CT perfusion were conducted in nineteen patients with PDAC on the day before surgery. IVIM analysis of DW-MRI was performed and the parameters perfusion fraction f, pseudodiffusion coefficient D*, and diffusion coefficient D were extracted for tumors, upstream, and downstream parenchyma. With a deconvolution-based analysis, the CT perfusion parameters blood flow (BF) and blood volume (BV) were estimated for tumors, upstream, and downstream parenchyma. In ten patients, intratumoral microvessel density (MVD_tumor) and microvessel area (MVA_tumor) were analyzed microscopically in resection specimens. Correlation coefficients between IVIM parameters, CT perfusion parameters, and histological microvessel parameters in tumors were calculated. Receiver operating characteristic (ROC) analysis was performed for differentiation of tumors and upstream parenchyma.

Results

f_tumor significantly positively correlated with BF_tumor (r = 0.668, p = 0.002) and BV_tumor (r = 0.672, p = 0.002). There were significant positive correlations between f_tumor and MVD_tumor/ MVA_tumor (r ≥ 0.770, p ≤ 0.009) as well as between BF_tumor and MVD_tumor/ MVA_tumor (r ≥ 0.697, p ≤ 0.025). Correlation coefficients between f_tumor and MVD_tumor/ MVA_tumor were not significantly different from correlation coefficients between BF_tumor and MVD_tumor/ MVA_tumor (p ≥ 0.400). Moreover, f, BF, BV, and permeability values (PEM) showed excellent performance in distinguishing tumors from upstream parenchyma (area under the ROC curve ≥0.874).

Conclusions

The study shows that IVIM derived f_tumor and CT perfusion derived BF_tumor similarly reflect vascularity of PDAC and seem to be comparably applicable for the evaluation of tumor perfusion for tumor characterization and as potential quantitative imaging biomarker.

Trial registration

DRKS, DRKS00022227, Registered 26 June 2020, retrospectively registered. https://www.drks.de/drks_web/navigate.do?navigationId=trial. HTML&TRIAL_ID=DRKS00022227.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40644-021-00382-x.

Renal Diffusion-Weighted Imaging (DWI) for Apparent Diffusion Coefficient (ADC), Intravoxel Incoherent Motion (IVIM), and Diffusion Tensor Imaging (DTI): Basic Concepts

The specialized function of the kidney is reflected in its unique structure, characterized by juxtaposition of disorganized and ordered elements, including renal glomerula, capillaries, and tubules. The key role of the kidney in blood filtration, and changes in filtration rate and blood flow associated with pathological conditions, make it possible to investigate kidney function using the motion of water molecules in renal tissue. Diffusion-weighted imaging (DWI) is a versatile modality that sensitizes observable signal to water motion, and can inform on the complexity of the tissue microstructure. Several DWI acquisition strategies are available, as are different analysis strategies, and models that attempt to capture not only simple diffusion effects, but also perfusion, compartmentalization, and anisotropy. This chapter introduces the basic concepts of DWI alongside common acquisition schemes and models, and gives an overview of specific DWI applications for animal models of renal disease.This chapter is based upon work from the COST Action PARENCHIMA, a community-driven network funded by the European Cooperation in Science and Technology (COST) program of the European Union, which aims to improve the reproducibility and standardization of renal MRI biomarkers. This introduction chapter is complemented by two separate chapters describing the experimental procedure and data analysis.

An optimized b-value distribution for triexponential intravoxel incoherent motion (IVIM) in the liver

PURPOSE

To find an optimized b-value distribution for reproducible triexponential intravoxel incoherent motion (IVIM) exams in the liver.

METHODS

A numeric optimization of b-value distributions was performed using the triexponential IVIM equation and 27 different IVIM parameter sets. Starting with an initially optimized distribution of 6 b-values, the number of b-values was increased stepwise. Each new b-value was chosen from a set of 64 predefined b-values based on the computed summed relative mean error of the fitted triexponential IVIM parameters. This process was repeated for up to 100 b-values. In simulations and in vivo measurements, optimized b-value distributions were compared to 4 representative distributions found in literature.

RESULTS

The first 16 optimized b-values were 0, 0.3, 0.3, 70, 200, 800, 70, 1, 3.5, 5, 70, 1.2, 6, 45, 1.5, and 60 in units of s/mm² . Low b-values were much more frequent than high b-values. The optimized b-value distribution resulted in a higher fit stability compared to distributions used in literature in both, simulation and in vivo measurements. Using more than 6 b-values, ideally 16 or more, increased the fit stability considerably.

CONCLUSION

Using optimized b-values, the fit uncertainty in triexponential IVIM can be largely reduced. Ideally, 16 or more b-values should be acquired.

On the dependence of the cardiac motion artifact on the breathing cycle in liver diffusion-weighted imaging

Purpose

The purpose of this study was to investigate whether the cardiac motion artifact that regularly appears in diffusion-weighted imaging of the left liver lobe might be reduced by acquiring images in inspiration, when the coupling between heart and liver might be minimal.

Materials and methods

43 patients with known or suspected focal liver lesions were examined at 1.5 T with breath hold acquisition, once in inspiration and once in expiration. Data were acquired with a diffusion-weighted echo planar imaging sequence and two b-values (b50 = 50 s/mm² and b800 = 800 s/mm²). The severity of the cardiac motion artifact in the left liver lobe was rated by two experienced radiologists for both b-values with a 5 point Likert scale. Additionally, the normalized signal S(b800)/S(b50) in the left liver lobe was computed. The Wilcoxon signed-rank test was used comparing the scores of the two readers obtained in inspiration and expiration, and to compare the normalized signal in inspiration and expiration.

Results

The normalized signal in inspiration was slightly higher than in expiration (0.349±0.077 vs 0.336±0.058), which would indicate a slight reduction of the cardiac motion artifact, but this difference was not significant (p = 0.24). In the qualitative evaluation, the readers did not observe a significant difference for b50 (reader 1: p = 0.61; reader 2: p = 0.18). For b800, reader 1 observed a significant difference of small effect size favouring expiration (p = 0.03 with a difference of mean Likert scores of 0.27), while reader 2 observed no significant difference (p = 0.62).

Conclusion

Acquiring the data in inspiration does not lead to a markedly reduced cardiac motion artifact in diffusion-weighted imaging of the left liver lobe and is in this regard not to be preferred over acquiring the data in expiration.

Investigation of intravoxel incoherent motion tensor imaging for the characterization of the in vivo human heart

PURPOSE

To investigate intravoxel incoherent motion (IVIM) tensor imaging of the in vivo human heart and elucidate whether the estimation of IVIM tensors is affected by the complexity of pseudo-diffusion components in myocardium.

METHODS

The cardiac IVIM data of 10 healthy subjects were acquired using a diffusion weighted spin-echo echo-planar imaging sequence along 6 gradient directions with 10 b values (0~400 s/mm² ). The IVIM data of left ventricle myocardium were fitted to the IVIM tensor model. The complexity of myocardial pseudo-diffusion components was reduced through exclusion of low b values (0 and 5 s/mm² ) from the IVIM curve-fitting analysis. The fractional anisotropy, mean fraction/mean diffusivity, and Westin measurements of pseudo-diffusion tensors (f_p and D*) and self-diffusion tensor (D), as well as the angle between the main eigenvector of f_p (or D*) and that of D, were computed and compared before and after excluding low b values.

RESULTS

The fractional anisotropy values of f_p and D* without low b value participation were significantly higher (P < .001) than those with low b value participation, but an opposite trend was found for the mean fraction/diffusivity values. Besides, after removing low b values, the angle between the main eigenvector of f_p (or D*) and that of D became small, and both f_p and D* tensors presented significant decrease of spherical components and significant increase of linear components.

CONCLUSION

The presence of multiple pseudo-diffusion components in myocardium indeed influences the estimation of IVIM tensors. The IVIM tensor model needs to be further improved to account for the complexity of myocardial microcirculatory network and blood flow.

Probing the ballistic microcirculation in placenta using flow-compensated and non-compensated intravoxel incoherent motion imaging

PURPOSE

Intravoxel incoherent motion (IVIM) imaging is widely used to evaluate microcirculatory flow, which consists of diffusive and ballistic flow components. We proposed a joint use of flow-compensated (FC) and non-compensated (NC) diffusion gradients to probe the fraction and velocity of ballistic flow in the placenta.

METHODS

Forty pregnant women were included in this study and scanned on a 1.5T clinical scanner. FC and NC diffusion MRI (dMRI) sequences were achieved using a pair of identical or mirrored bipolar gradients. A joint FC-NC model was established to estimate the fraction (f_b ) and velocity (v_b ) of the ballistic flow. Conventional IVIM parameters (f, D, and D*) were obtained from the FC and NC data, separately. The v_b and f·D*, as placental flow velocity measurements, were correlated with the umbilical-artery Doppler ultrasound indices and gestational ages.

RESULTS

The ballistic flow component can be observed from the difference between the FC and NC dMRI signal decay curves. v_b fitted from the FC-NC model showed strong correlations with umbilical-artery impedance indices, the systolic-to-diastolic (SD) ratio and pulsatility index (PI), with correlation coefficients of 0.65 and 0.62. The f·D* estimated from the NC data positively correlated with SD and PI, while the FC-based f·D* values showed weak negative correlations. Significant gestational-age dependence was also found in the flow velocity measurements.

CONCLUSION

Our results demonstrated the feasibility of using FC and NC dMRI to noninvasively measure ballistic flow velocity in the placenta, which may be used as a new marker to evaluate placenta microcirculation.

Quantitative imaging for radiotherapy purposes

Quantitative imaging biomarkers show great potential for use in radiotherapy. Quantitative images based on microscopic tissue properties and tissue function can be used to improve contouring of the radiotherapy targets. Furthermore, quantitative imaging biomarkers might be used to predict treatment response for several treatment regimens and hence be used as a tool for treatment stratification, either to determine which treatment modality is most promising or to determine patient-specific radiation dose. Finally, patient-specific radiation doses can be further tailored to a tissue/voxel specific radiation dose when quantitative imaging is used for dose painting. In this review, published standards, guidelines and recommendations on quantitative imaging assessment using CT, PET and MRI are discussed. Furthermore, critical issues regarding the use of quantitative imaging for radiation oncology purposes and resultant pending research topics are identified.