Diffusion Tensor Imaging (DTI) of the kidneys incorporating advanced geometric distortion correction using reversed phase encoding images.
Jose Teruel1,2, Jeremy C. Lim3, Eric E. Sigmund4, Elissa Botterill5, Jas-mine Seah6, Shawna Farquharson7, Elif E. Ekinci6,8, and Ruth P. Lim5,9

1Department of Circulation and Medical Imaging, Norwegian University of Science and Technology, Trondheim, Norway, 2St. Olavs University Hospital, Trondheim, Norway, 3Department of Radiology, The Royal Melbourne Hospital, Melbourne, Australia, 4Department of Radiology, NYU Langone Medical Center, New York, NY, United States, 5Department of Radiology, Austin Health, Melbourne, Australia, 6Department of Endocrinology, Austin Health, Melbourne, Australia, 7Florey Neuroscience Institute, Melbourne, Australia, 8Department of Endocrinology, The University of Melbourne, Melbourne, Australia, 9Departments of Radiology and Surgery, The University of Melbourne, Melbourne, Australia

Synopsis

Diffusion tensor imaging is emerging as a promising technique for structural and functional evaluation of the kidneys. However, diffusion sequences employing echo planar imaging readout are prone to geometric distortions due to static field inhomogeneities arising from different magnetic susceptibilities from adjacent tissues and bowel gas. In this study, we evaluated the efficacy of distortion correction using a reversed phase encoding approach for diffusion tensor imaging of healthy controls and patients with Type 1 diabetes.

Purpose

Diffusion tensor imaging (DTI) is emerging as a promising technique for structural and functional evaluation of the kidneys1,2. Diffusion sequences employing echo planar imaging (EPI) readout are prone to geometric distortions due to static field inhomogeneities arising from different susceptibilities. Geometric distortion correction has its main application in brain MRI, but recent studies have shown its applicability to body organs3,4. There is little experience with this technique for renal imaging, which is subject to distortions due to surrounding organs and adjacent bowel gas. The purpose of this work is to evaluate efficacy of distortion correction using EPI-based diffusion images acquired with opposite phase encoding (PE) polarities, in healthy controls and patients with Type 1 diabetes (T1DM).

Methods

Ten T1DM patients (6M/4F, mean age [range] 44[19-62]y, mean eGFR [range] 84 [55-91], and six healthy volunteers (3M/3F, mean age [range] 44[28-63]y) were scanned using a 3T system (Skyra, Siemens) including a respiratory-navigated DTI sequence with the following parameters: TR/TE=1800/81 ms, FOV=400x400 mm2, matrix=340x340, slice thickness 3 mm, 7 slices, b-values 0 and 500 s/mm2 in 12 directions and 4 averages with a feet-to-head [FH] PE and EPI readout in oblique/coronal orientation. Additionally, one b0 acquisition with identical parameters and reversed head-to-feet [HF] PE was acquired. Breath-hold T2-weighted HASTE images (T2WI) were obtained covering the full extent of the kidneys with: TR/TE=1000/95 ms, FOV=420x420 mm2, matrix=640x640, slice thickness 5 mm in identical orientation. Geometric distortion correction of DTI was carried out according to the method described in5 using CMTK6. An experienced abdominal radiologist selected a single representative slice from T2WI and matching slice from DTI datasets for subsequent evaluation. The same radiologist segmented bilateral kidneys in each subject on the selected T2WI slice, to obtain an area as the reference standard for each kidney. After a two week period, the same radiologist and an inexperienced resident radiologist independently segmented the matching uncorrected and corrected DTI-b0 images in random order. Segmentation and parametric maps were produced using non-commercial software developed using Matlab R2014. DTI parametric maps were fitted for both uncorrected and corrected series including the apparent diffusion coefficient (ADC), fractional anisotropy (FA), and the three tensor eigenvalues (λ12 and λ3). Segmented area was compared with Pearson correlation and Wilcoxon signed ranks test. Intraclass correlation (ICC) between readers was assessed.

Results

Geometric distortion correction was successfully performed in all cases (Figure 1). Table 1 presents kidney area segmented on T2WI, uncorrected DTI-b0, and corrected DTI-b0 images. Excellent agreement was found between both readers with an ICC of 0.97 and 0.98 for kidney segmentation in uncorrected DTI-b0 and corrected DTI-b0 respectively. Figure 2 presents a scatter/boxplot for direct comparison of T2WI segmented area with areas obtained by the experienced reader for uncorrected and corrected DTI-b0 images. Pearson correlation between segmented area on T2WI and corrected DTI-b0 was higher (r=0.904; p<0.001), than the correlation between the area segmented in T2WI and uncorrected DTI-b0 (r=0.840; p<0.001). Median [range] of DTI parametric maps for all cases are presented in Table 2 using the experienced reader segmentation. Figure 3 shows ADC and FA parametric maps for one case before and after distortion correction.

Discussion / Conclusion

Our results show there is a significant improvement in the agreement between the area segmented on T2WI (the anatomic reference standard) and the area segmented on corresponding slices in the DTI series after distortion correction. This includes T1DM patients with a spectrum of renal function. This has implications for more reliable segmentation of renal cortex and medulla in DTI, with better correlation to higher resolution anatomic images, of importance when renal pathology is present. In our small cohort, the parameters derived from DTI did not seem to be affected. This observation strengthens its applicability as parameter quantification would not be affected by the processing routine. Nevertheless, it is important to consider that statistical comparison was carried out accounting only for the median value of each parameter in the single-slice kidney segmentation. In this study the PE trajectory for the complete DTI acquisition was chosen as FH, based on7. Therefore, geometric distortion resulted in artefactual elongation of the kidney (compression in the reversed HF b0 image) in the PE direction, as clearly depicted in Figure 1. The choice of PE in this study required that for most cases, the stomach needed to be masked out of the images before correction, as it caused significant artifact affecting the left kidney. Future work to evaluate how DTI parameters might be affected with smaller ROIs to segment renal cortex from medulla using DTI parametric maps with reference to anatomic imaging is planned.

Acknowledgements

This study was supported by a Diabetes Australia Research Trust Grant.

References

1. Notohamiprodjo M, Glaser C, Herrmann KA, et al. Diffusion tensor imaging of the kidney with parallel imaging: initial clinical experience. Invest Radiol 2008;43:677–685.

2. Sigmund EE, Vivier PH, Sui D, et al. Intravoxel incoherent motion and diffusion-tensor imaging in renal tissue under hydration and furosemide ?ow challenges. Radiology 2012;263:758–769.

3. Teruel JR, Fjøsne HE, Østlie A, et al. Inhomogeneous static magnetic field-induced distortion correction applied to diffusion weighted MRI of the breast at 3T.Magn Reson Med 2015;74(4):1138-44.

4. Rakow-Penner RA, White NS, Margolis DJ, et al. Prostate diffusion imaging with distortion correction. Magn Reson Imaging 2015;33:1178-81.

5. Holland D, Kuperman JM, Dale AM. Efficient correction of inhomogeneous static magnetic field-induced distortion in Echo Planar Imaging. NeuroImage 2010;50:175:83.

6. Computational Morphometry Toolkit. https://www.nitrc.org/projects/cmtk/ (Accessed November 10, 2015)

7. Notohamiprodjo M, Dietrich O, Horger W, et al. Diffusion tensor imaging (DTI) of the kidney at 3 tesla-feasibility, protocol evaluation and comparison to 1.5 Tesla. Invest Radiol 2010;45:245-54.

Figures

Figure 1. DTI b0 image (28 year old male, healthy volunteer) before correction with A: feet-head (FH) phase encoding direction; B: head-feet (HF) phase encoding direction. C: DTI b0 corrected image (showing FH corrected image). D: Corresponding slice from T2WI series. The corrected image most closely approximates renal anatomy as shown on T2WI.

Table 1. Median [full range] of the area segmented for each imaging modality divided by reader and laterality. All results expressed in cm2.

Figure 2. Scatter/boxplot showing the area segmented by the experienced reader for the 16 cases (both kidneys) in the three scenarios: T2-weighted HASTE, corresponding slice of uncorrected b0 DTI, and corresponding slice of corrected b0 DTI. Median (solid line). Mean (dashed line). Pearson correlation between T2WI and uncorrected b0 (r=0.840; p<0.001); between T2WI and corrected b0 (r=0.904; p<0.001). Wilcoxon signed ranks test for pairwise related samples comparison.

Table 2. Median [full range] of the DTI parameters obtained for each case. Median value for each kidney analyzed was employed (n=32). No significant differences were found for each of the parameters reported between uncorrected and corrected DTI groups using a non-parametric Mann-Whitney U test.

Figure 3. A: Apparent diffusion coefficient (ADC) and fractional anisotropy (FA) parametric maps of the left kidney for one patient with Type 1 diabetes (50 year old male, eGFR 91) before correction (top), and after correction (bottom) overlaid over DTI b0 images. B: Corresponding T2WI slice.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
0271