Spectral diffusion analysis of kidney intravoxel incoherent motion MRI in healthy volunteers and patients with renal pathologies

Investigation of diffusion time dependence of apparent diffusion coefficient and intravoxel incoherent motion parameters in the human kidney

PURPOSE

To characterize the diffusion time (Δeff) dependence of apparent diffusion coefficient (ADC) and intravoxel incoherent motion-related parameters in the human kidney at 3 T.

METHODS

Sixteen healthy volunteers underwent an MRI examination at 3 T including diffusion-weighted imaging at different Δeff ranging from 24.1 to 104.1 ms. The extended mono-exponential ADC and intravoxel incoherent motion models were fitted to the data for each Δeff and the medullary and cortical ADC, (pseudo-)diffusion coefficients (D* and D) and flow-related signal fraction (f) were calculated.

RESULTS

When all the data were used for fitting, a significant trend toward higher ADC with increasing Δeff was observed between 24.1 and 104.1 ms (median and interquartile range: 2.38 [2.19, 2.47] to 2.84 [2.36, 2.90] × 10-3 mm2/s for cortex, and 2.28 [2.18, 2.37] to 2.82 [2.58, 3.11] × 10-3 mm2/s for medulla). In contrast, no significant differences in ADC were found when only the data acquired at b-values higher than 200 s/mm2 were used for fitting. When the intravoxel incoherent motion model was applied, cortical and medullary f increased significantly (cortex: 0.21 [0.15 0.27] to 0.37 [0.32, 0.49] × 10-3 mm2/s; medulla: 0.15 [0.13 0.29] to 0.41 [0.36 0.51] × 10-3 mm2/s). No significant changes in cortical and medullary D and D* were observed as diffusion time increased.

CONCLUSION

Renal perfusion and tubular flow substantially contribute to the observed increase in ADC over a wide range of Δeff between 24 and 104 ms.

Presentation of microstructural diffusion components by color schemes in abdominal organs

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

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.

MRI of kidney size matters

OBJECTIVE

To highlight progress and opportunities of measuring kidney size with MRI, and to inspire research into resolving the remaining methodological gaps and unanswered questions relating to kidney size assessment.

MATERIALS AND METHODS

This work is not a comprehensive review of the literature but highlights valuable recent developments of MRI of kidney size.

RESULTS

The links between renal (patho)physiology and kidney size are outlined. Common methodological approaches for MRI of kidney size are reviewed. Techniques tailored for renal segmentation and quantification of kidney size are discussed. Frontier applications of kidney size monitoring in preclinical models and human studies are reviewed. Future directions of MRI of kidney size are explored.

CONCLUSION

MRI of kidney size matters. It will facilitate a growing range of (pre)clinical applications, and provide a springboard for new insights into renal (patho)physiology. As kidney size can be easily obtained from already established renal MRI protocols without the need for additional scans, this measurement should always accompany diagnostic MRI exams. Reconciling global kidney size changes with alterations in the size of specific renal layers is an important topic for further research. Acute kidney size measurements alone cannot distinguish between changes induced by alterations in the blood or the tubular volume fractions-this distinction requires further research into cartography of the renal blood and the tubular volumes.

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.

Noninvasive Assessment of Diabetic Kidney Disease With MRI: Hype or Hope?

Owing to the increasing prevalence of diabetic mellitus, diabetic kidney disease (DKD) is presently the leading cause of chronic kidney disease and end-stage renal disease worldwide. Early identification and disease interception is of paramount clinical importance for DKD management. However, current diagnostic, disease monitoring and prognostic tools are not satisfactory, due to their low sensitivity, low specificity, or invasiveness. Magnetic resonance imaging (MRI) is noninvasive and offers a host of contrast mechanisms that are sensitive to pathophysiological changes and risk factors associated with DKD. MRI tissue characterization involves structural and functional information including renal morphology (kidney volume (TKV) and parenchyma thickness using T1- or T2-weighted MRI), renal microstructure (diffusion weighted imaging, DWI), renal tissue oxygenation (blood oxygenation level dependent MRI, BOLD), renal hemodynamics (arterial spin labeling and phase contrast MRI), fibrosis (DWI) and abdominal or perirenal fat fraction (Dixon MRI). Recent (pre)clinical studies demonstrated the feasibility and potential value of DKD evaluation with MRI. Recognizing this opportunity, this review outlines key concepts and current trends in renal MRI technology for furthering our understanding of the mechanisms underlying DKD and for supplementing clinical decision-making in DKD. Progress in preclinical MRI of DKD is surveyed, and challenges for clinical translation of renal MRI are discussed. Future directions of DKD assessment and renal tissue characterization with (multi)parametric MRI are explored. Opportunities for discovery and clinical break-through are discussed including biological validation of the MRI findings, large-scale population studies, standardization of DKD protocols, the synergistic connection with data science to advance comprehensive texture analysis, and the development of smart and automatic data analysis and data visualization tools to further the concepts of virtual biopsy and personalized DKD precision medicine. We hope that this review will convey this vision and inspire the reader to become pioneers in noninvasive assessment and management of DKD with MRI. LEVEL OF EVIDENCE: 1 TECHNICAL EFFICACY: Stage 2.

MR Perfusion Imaging for Kidney Disease

Renal perfusion reflects overall function of a kidney. As an important indicator of kidney diseases, renal perfusion can be noninvasively measured by multiple methods of MR imaging, such as dynamic contrast-enhanced MR imaging, intravoxel incoherent motion analysis, and arterial spin labeling method. In this article we introduce the principle of the methods, review their recent technical improvements, and then focus on summarizing recent applications of the methods in assessing various renal diseases. By this review, we demonstrate the capability and clinical potential of the imaging methods, with the hope of accelerating their adoption to clinical practice.

Renal MRI: From Nephron to NMR Signal

Renal diseases pose a significant socio-economic burden on healthcare systems. The development of better diagnostics and prognostics is well-recognized as a key strategy to resolve these challenges. Central to these developments are MRI biomarkers, due to their potential for monitoring of early pathophysiological changes, renal disease progression or treatment effects. The surge in renal MRI involves major cross-domain initiatives, large clinical studies, and educational programs. In parallel with these translational efforts, the need for greater (patho)physiological specificity remains, to enable engagement with clinical nephrologists and increase the associated health impact. The ISMRM 2022 Member Initiated Symposium (MIS) on renal MRI spotlighted this issue with the goal of inspiring more solutions from the ISMRM community. This work is a summary of the MIS presentations devoted to: 1) educating imaging scientists and clinicians on renal (patho)physiology and demands from clinical nephrologists, 2) elucidating the connection of MRI parameters with renal physiology, 3) presenting the current state of leading MR surrogates in assessing renal structure and functions as well as their next generation of innovation, and 4) describing the potential of these imaging markers for providing clinically meaningful renal characterization to guide or supplement clinical decision making. We hope to continue momentum of recent years and introduce new entrants to the development process, connecting (patho)physiology with (bio)physics, and conceiving new clinical applications. We envision this process to benefit from cross-disciplinary collaboration and analogous efforts in other body organs, but also to maximally leverage the unique opportunities of renal physiology. LEVEL OF EVIDENCE: 1 TECHNICAL EFFICACY STAGE: 2.

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.

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

PURPOSE

To improve the reliability of intravoxel incoherent motion (IVIM) model parameter estimation for the DWI in the kidney using a novel image downsampling expedited adaptive least-squares (IDEAL) approach.

METHODS

The robustness of IDEAL was investigated using simulated DW-MRI data corrupted with different levels of Rician noise. Subsequently, the performance of the proposed method was tested by fitting bi- and triexponential IVIM model to in vivo renal DWI data acquired on a clinical 3 Tesla MRI scanner and compared to conventional approaches (fixed D* and segmented fitting).

RESULTS

The numerical simulations demonstrated that the IDEAL algorithm provides robust estimates of the IVIM parameters in the presence of noise (SNR of 20) as indicated by relatively low absolute percentage bias (maximal sMdPB <20%) and normalized RMSE (maximal RMSE <28%). The analysis of the in vivo data showed that the IDEAL-based IVIM parameter maps were less noisy and more visually appealing than those obtained using the fixed D* and segmented methods. Further, coefficients of variation for nearly all IVIM parameters were significantly reduced in cortex and medulla for IDEAL-based biexponential (coefficients of variation: 4%-50%) and triexponential (coefficients of variation: 7.5%-75%) IVIM modelling compared to the segmented (coefficients of variation: 4%-120%) and fixed D* (coefficients of variation: 17%-174%) methods, reflecting greater accuracy of this method.

CONCLUSION

The proposed fitting algorithm yields more robust IVIM parameter estimates and is less susceptible to poor SNR than the conventional fitting approaches. Thus, the IDEAL approach has the potential to improve the reliability of renal DW-MRI analysis for clinical applications.

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.

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.

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.