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 kidneys
1,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 organs
3,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 mm
2, matrix=340x340, slice thickness 3 mm,
7 slices, b-values 0 and 500 s/mm
2 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 mm
2, matrix=640x640, slice thickness 5 mm in identical orientation.
Geometric distortion
correction of DTI was carried out according to the method described in
5
using CMTK
6. 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
(λ
1,λ
2 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 on
7. 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.