Characterizing diffusion heterogeneity changes after acute ischemic stroke
Ona Wu1, Arne Lauer2, Gregoire Boulouis2, Lisa Cloonan2, Mark Etherton2, Abigail S. Cohen2, Pedro T. Cougo-Pinto2, Katherine Mott1, William A. Copen3, and Natalia S. Rost2

1Athinoula A Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, United States, 2Department of Neurology, Massachusetts General Hospital, Boston, MA, United States, 3Department of Radiology, Massachusetts General Hospital, Boston, MA, United States

Synopsis

Diffusion kurtosis imaging (DKI) has been suggested to be a more sensitive marker for microstructural injury than diffusion tensor imaging (DTI). To investigate this hypothesis, we analyzed DKI data from acute ischemic stroke patients enrolled in a prospective serial MRI study (N=18). Axial diffusivity and axial kurtosis values within the ischemic core were significantly correlated with time-to-MRI. Regional differences in both diffusivity and kurtosis were observed as a function of tissue outcome suggesting DKI may provide complementary information to that obtained from DTI.

Purpose

Diffusional kurtosis is a measure of the non-Gaussianity of the diffusion process of water molecules in tissue.1 DKI may be a more sensitive marker for pathophysiologic changes in cellular microstructure after acute ischemic stroke (AIS) than diffusion tensor imaging (DTI) metrics such as fractional anisotropy (FA).1, 2 Furthermore, in a subset of AIS patients who undergo successful revascularization treatment, lesions on acute trace diffusion-weighed MRI (DWI) have been noted to “reverse”.3, 4 Coupled with changes in diffusivity, DKI may provide complementary insight into the extent of neuronal injury.2 We therefore investigated DKI metrics as a function of tissue outcome in AIS patients.

Methods

DKI from 18 AIS patients enrolled in a prospective serial MRI study were analyzed. The first MRI was acquired for clinical purposes within 12 h from when the patient was last-known-well on a 1.5T GE scanner, and the second MRI prior to discharge on a 3T Siemens system. MR perfusion-weighted imaging (MRP) was acquired at both sessions with gradient echo echo-planar imaging (TR/TE=1500/40 ms at 1.5T and TR/TE=1500/35 ms at 3T, 80 timepoints) during the first pass of an intravenous bolus injection of a gadolinium-based contrast agent. Perfusion maps (CBF, CBV, MTT and Tmax) were calculated using deconvolution5 with an automatically selected arterial input function.6 DKI was acquired at the second time point using 30 directions with b-value=1000 s/mm2, 2000 s/mm2 and 10 b-value=0 s/mm2 images (3x3x3 mm3) using simultaneous multislice image acquisition.7 Mean kurtosis (MK), axial kurtosis (AK), radial kurtosis (RK), mean diffusivity (MD), axial diffusivity (AD), radial diffusivity (RD) and fractional anisotropy (FA) maps were calculated using the Diffusion Kurtosis Estimator.8 For delineating region of interests (ROI), acute DWI, DKI, MRP and the follow-up infarct (FU) from the second MRI’s FLAIR images were co-registered to one another (MNI autoreg9), and to the MNI 152 1mm atlas. Abnormal perfusion was defined as tissue with Tmax values greater than 6 seconds. DWI and DKI were compared in the following regions: Core (abnormal acute DWI and FU), Growth (normal acute DWI, abnormal FU), and Salvaged (normal acute DWI, abnormal acute perfusion, normal FU). Contralateral normal tissue (CNL) was defined as tissue with no apparent abnormalities on acute and follow-up imaging. DWI and DKI metrics in the Core were correlated against time-to-MRI (nonparametric Spearman’s correlation analysis). To minimize errors due to poor co-registration, only patients >1 cm3 Growth or Saved ROIs were included for analysis. Each ROI was further segmented into gray matter (GM), white matter (WM) and cerebral spinal fluid. Mean acute DTI (MD, AD, RD, FA) and DKI values (MK, AK, RK) in the four ROIs: Core, Growth, Saved and CNL as a whole and in GM and WM areas were evaluated using one-way ANOVA with repeated measures followed by a post-hoc Tukey HSD test.

Results

Patient characteristics were: mean±SD age 66.2±10.4 years, 72.2% (N=13) males, median [IQR] admission NIH Stroke Scale (NIHSS) 6 [2.75-11], time-to-acute MRI 6.2±2.1 h, time-to-F/u MRI 3.0 ± 1.3 days, 33% (N=6) received thrombolysis, acute DWI lesion 3.5 [0.6-20.4] cm3, acute Tmax lesion 10.5 [0-73.1] cm3, FU lesion 13.4 [1.8-56.3] cm3 and 90 day modified Rankin Scale 1 [1-2.25]. Figure 1 shows examples of the maps used for ROI placement, and example DTI and DKI maps. Increased kurtosis is evident, and most conspicuous in AK maps, in which normal GM and WM values are similar. Correlation between time-to-MRI and diffusion metrics were significant only for AD (ρ=-0.54, P=0.022), MK (ρ=0.37, P=0.049), and AK (ρ=0.55, P=0.017). Nine patients exhibited salvaged tissue. Significant differences were found between the ROIs for both GM (Figure 2) and WM (Figure 3). For GM, diffusivity metrics were significantly lower in all ROIs compared to CNL and were significantly with respect to one other. FA was lower in the Core than in Saved tissue for both GM and WM. In WM, Growth ROIs had greater FA values than Core, but smaller values than Saved and CNL. MK and AK values in the Core and Growth ROIs were significantly higher than both CNL and Saved ROI values for GM. For WM, this was true only for AK.

Discussion

The significant inverse correlation between time-to-MRI and axial diffusivity changes in core tissue suggests that ischemia-induced cellular swelling may increase tortuosity of water diffusion paths, imposing direction-dependent restrictions upon diffusion. Coupled with changes in diffusivity, DKI may provide complementary insight into neuronal injury.2 However, future studies are needed in the hyperacute stage to fully understand the role of DTI and DKI in identifying salvageable tissue.

Acknowledgements

We thank Drs. Himanshu Bhat, Kawin Setsompop and Steven Cauley for providing the simultaneous-multislice pulse sequence used for the DKI acquisition.

References

1. Jensen JH, Helpern JA, Ramani A, Lu H, Kaczynski K. Diffusional kurtosis imaging: The quantification of non-gaussian water diffusion by means of magnetic resonance imaging. MRM. 2005;53:1432-1440

2. Hui ES, Fieremans E, Jensen JH, Tabesh A, Feng W, Bonilha L, et al. Stroke assessment with diffusional kurtosis imaging. Stroke. 2012;43:2968-2973

3. Fiehler J, Knudsen K, Kucinski T, Kidwell CS, Alger JR, Thomalla G, et al. Predictors of apparent diffusion coefficient normalization in stroke patients. Stroke. 2004;35:514-519

4. Kidwell CS, Saver JL, Mattiello J, Starkman S, Vinuela F, Duckwiler G, et al. Thrombolytic reversal of acute human cerebral ischemic injury shown by diffusion/perfusion magnetic resonance imaging. Ann Neurol. 2000;47:462-469.

5. Wu O, Ostergaard L, Weisskoff RM, Benner T, Rosen BR, Sorensen AG. Tracer arrival timing-insensitive technique for estimating flow in MR perfusion-weighted imaging using singular value decomposition with a block-circulant deconvolution matrix. Magn Reson Med. 2003;50:164-174

6. Copen WA, Deipolyi AR, Schaefer PW, Schwamm LH, Gonzalez RG, Wu O. Exposing hidden truncation-related errors in acute stroke perfusion imaging. AJNR Am J Neuroradiol. 2015;36:638-645

7. Setsompop K, Gagoski BA, Polimeni JR, Witzel T, Wedeen VJ, Wald LL. Blipped-controlled aliasing in parallel imaging for simultaneous multislice echo planar imaging with reduced g-factor penalty. Magn Reson Med. 2012;67:1210-1224

8. Tabesh A, Jensen JH, Ardekani BA, Helpern JA. Estimation of tensors and tensor-derived measures in diffusional kurtosis imaging. Magn Reson Med. 2011;65:823-836

9. Collins DL, Neelin P, Peters TM, Evans AC. Automatic 3d intersubject registration of MR volumetric data in standardized talairach space. J Comput Assist Tomogr. 1994;18:192-205

Figures

Figure 1: Example of DTI and DKI changes in a patient imaged 9h from LKW and 7 days later with DKI. Despite reperfusion, lesion expansion is noted from acute DWI lesion (red arrow). There are evident increases in kurtosis in the region of the F/u lesion (green arrow). There are also subtle decreases in FA (yellow arrow).

Figure 2: Differences between diffusivity and kurtosis metrics in core, growth and saved regions of interests in gray matter. There were significant differences (P<0.05) across the following regions – MD: All regions; AD: All regions; RD: All regions; FA: Core vs Growth, Core vs Saved; MK: CNL vs Core, CNL vs Growth, Core vs Saved, Growth vs Saved; AK: CNL vs Core, CNL vs Growth, Core vs Saved, Growth vs Saved; RK: CNL vs Core, CNL vs Growth.

Figure 3: Differences between diffusivity and kurtosis metrics in core, growth and saved regions of interests in white matter. There were significant differences (P<0.05) across the following regions – MD: All regions except Core vs Growth; AD: All regions except Core vs Growth; RD: CNL vs Core, CNL vs Growth; FA: All regions except CNL vs Saved; MK: CNL vs Core, CNL vs Growth, CNL vs Saved; AK: All regions except Core vs Growth.



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