Complete partial volume solution for ASL brain perfusion data applied to relapsing-remitting multiple sclerosis patients
Ruth Oliver1,2, Linda Ly1,2, Chenyu Wang1,2, Heidi Beadnall2, Ilaria Boscolo Galazzo3,4, Michael Chappell5,6, Xavier Golay7, Enrico De Vita7, David Thomas7, and Michael Barnett1,2

1Sydney Neuroimaging Analysis Centre, Sydney, Australia, 2University of Sydney, Sydney, Australia, 3Institute of Nuclear Medicine, University College London, London, United Kingdom, 4Department of Neuroradiology, University Hospital Verona, Verona, Italy, 5Institute of Biomedical Engineering, University of Oxford, Oxford, United Kingdom, 6FMRIB Centre, University of Oxford, Oxford, United Kingdom, 7Institute of Neurology, University College London, London, United Kingdom

Synopsis

ASL is a low resolution imaging modality that suffers from the partial volume effect, leading to an underestimation of GM perfusion. This effect has two principle causes; blurring from the point spread function in the slice direction, and inadequate resolution due to the need for large voxels to achieve sufficient SNR. Both may act as confounders for measurement of GM CBF abnormalities. Decreased GM perfusion could reflect neuronal loss or metabolic dysfunction; PV correction allows a decoupling of structure and function. We present the first application of a complete PV correction solution for ASL to a cohort of MS patients.

Purpose

To apply partial volume correction to ASL data from patients with MS to improve the accuracy of grey matter (GM) CBF estimates and facilitate the measurement of perfusion changes as a biomarker of treatment response.

Introduction

ASL is a low resolution imaging modality known to suffer from the partial volume (PV) effect, leading to an underestimation of GM perfusion. The PV effect has two principle causes: blurring from the point spread function (PSF) in the slice direction, due to the long echo train employed, and inadequate sampling resolution due to the need for large voxels to achieve sufficient SNR1. Both these effects may act as confounders for the measurement of GM CBF abnormalities. A decrease in GM perfusion could reflect neuronal loss or metabolic dysfunction; PV correction techniques allow decoupling of structure and function. It is hypothesized that reduced GM perfusion after PV correction may serve as a biomarker of reversible neuronal dysfunction prior to substantive tissue loss2. In this work, we present the first application of a complete PV correction solution for ASL to a cohort of MS patients.

Methods

Six RRMS patients (5F, 1M aged 24-56Y, mean=43.5Y, EDSS scores ranged from 1.0 to 3.5; mean=2.3) were imaged on a GE MR750 3T scanner with 8 channel head coil. Acquisitions included: T2 FLAIR (TE/TR=162/8000ms, TI=2181ms, matrix=512x512x480, resolution=0.47x0.47x0.6mm3) for lesion characterization and T1-weighted anatomical IR-SPGR for tissue concentration (TE/TR=2.8/7.1ms, TI=450ms, FA=12°, matrix 512x512x248, resolution=0.47x0.47x0.7mm3). The ASL sequence was pCASL labeling with 8-arm 3D stack-of-spirals read-out (inflow time=1525ms, bolus length=1525ms, matrix 128x128x32, resolution=1.88x1.88x5mm3). MS lesions were manually segmented on the FLAIR image using Jim3 software to create lesion masks. FSL's4 lesion infilling was used to in-paint the lesions on the anatomical image, prior to segmentation into GM, WM, CSF components using FSL FAST. The segmentations were downsampled to perfusion space using FSL’s flirt and applywarp functions to ensure the PV maps within each ASL voxel represented the average PV estimate across the high-resolution region5. The 1D PSF in the slice direction was estimated using the Extended Phase Graph algorithm for the difference and proton density images6. These images were both deblurred using the calculated PSF with a Richardson-Lucy deconvolution algorithm, in order to restore signal to the voxel of origin. The images were PV corrected using a 3 x 3 x 3 kernel in a linear regression algorithm before combining to produce separate maps of GM and WM CBF7. Regions of interest (ROI) were created from a Freesurfer parcellation8 which was transformed to ASL space in the same manner as the PV estimates and the mean GM or WM CBF was calculated for each of these.

Results

The 1D PSF in the slice direction was estimated to be 1.6 and 1.3 voxels at full width half maximum for the ASL difference and proton density images respectively. ASL and proton density images for patient 002NN before and after deblurring are shown in fig. 1,2. Middle slice CBF and PV-corrected GM and WM CBF are shown in fig. 3. Fig. 4 shows the mean GM and WM CBF with and without deblurring for all 6 subjects, which shows a broad trend of increased GM and reduced WM CBF. There was a consistent increase in GM CBF with the inclusion of deblurring compared to without, ranging from 5-12%, mean 8%. In 4 out of 6 subjects, the WM CBF decreased by -2 to -11%, mean -5%. For subject 004SS, there was no change, and a 1% increase for 006JF. Fig. 5 shows the mean regional CBF and GM CBF post-correction, as well as the WM CBF.

Discussion

We observed significant through-plane blurring in 3D ASL images, even in highly-segmented acquisitions. The combination of deblurring and PV-correction produces an increase in mean GM CBF and a decrease in WM CBF over and above that produced by either correction alone. GM CBF values are modified to a greater extent than WM CBF, due to the cortex being significantly thinner than the width of the PSFs and GM voxels therefore suffering more signal attenuation. WM CBF values for these MS patients were higher than usually observed in healthy cohorts. Elevated WM CBF preceding lesion development has previously been observed in MS patients, which suggests that CBF may be sensitive to inflammation and serve as a early stage biomarker of lesion formation9,10.

Conclusion

More accurate estimations of GM and WM CBF can be obtained by applying PV-correction solutions that correct for both PSF blurring and tissue heterogeneity. This will be particularly important for SPMS patients, who exhibit tissue atrophy over the disease course.

Acknowledgements

Thanks to Fernando Zelaya for his assistance with the GE ASL sequence.

We acknowledge partial funding from Biogen for this study.

References

1. Asllani, I. MRM 2008 2. Debernard, L. JNNP 2014 3. www.xinapse.com 4. http://fsl.fmrib.ox.ac.uk/fsl 5. Chappell, MA. MRM 2011 6. Hennig, J. J Mag Res 1988 7 Oliver, RA. PhD thesis 2015. 8. http://freesurfer.net 9. Paling, D. J Cere Blood Flow & Metab. 2013 1-9 10. Wuerfel, J. Brain 2004 127:111-119

Figures

Fig. 1 Original, acquired difference images (LHS) and after deblurring with 1D PSF in slice direction of FWHM 1.6 voxels.

Fig. 2 Original, acquired proton density images (LHS) and after deblurring with 1D PSF in slice direction of FWHM 1.3 voxels.

Fig. 3 LHS: Uncorrected CBF; Middle: GM CBF; RHS: WM CBF for middle slice from patient 002NN. Complete PV correction with deblurring component is applied to difference and proton density images, GM and WM CBF masked at PVF < 0.1.

Fig. 4 Comparison of the effect of including deblurring of the PSF to PV correction technique. GM CBF is consistently elevated, with WM CBF reducing in 4 out of 6 cases. No change for one subject, and 1% increase for other.

Fig. 5 Regional mean CBF and GM and WM CBF after partial volume correction for all six RRMS patients.



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