3103

To Evaluate Effect of Interpolation Methods for B0 inhomogeneity correction on computation of CEST MRI Contrast in Human Brain at 3T and 7T
Anup Singh1,2, Ayan Debnath1,3, Rakesh Kumar Gupta4, and Ravinder Reddy3
1Centre for Biomedical Engineering, Indian Institute of Technology Delhi, New Delhi, India, 2Biomedical Engineering, AIIMS, New Delhi, India, 3Radiology, University of Pennsylvania, Philadelphia, PA, United States, 4Radiology, Fortis Memorial Research Institute, Gurugram, India

Synopsis

CEST MRI provides high resolution mapping of molecules such as Glutamate, Creatine, labile proteins/peptide, etc. CEST MRI contrast computation using asymmetry analysis require interpolation of data at different frequency offsets for B0 inhomogeneity corrections. In this study, different interpolation methods such as polynomial, cubic, spline and smoothingspline were compared for B0 inhomogeneity correction of various CEST contrasts(GluCEST-w, APT-w, CrCEST-w) at different field strength (3T, 7T). The 2nd and 3rd degree polynomial interpolations provided better B0 inhomogeneity correction for in vivo data from human brain. Polynomial interpolations for APT-w also improved differentiation of high-grade-glioma and low-grade-glioma tumors at 3T.

Introduction

Chemical-Exchange-Saturation-Transfer(CEST)1 MRI is a non-invasive molecular imaging technique that has shown potential applications in pre-clinical2 and clinical studies3. Applications of Amide-proton-transfer-weighted(APT-w)4-6 CEST, Glutamate-weighted CEST (GluCEST-w)7-8, Creatine-weighted CEST (CrCEST-w)9,10, etc. have been reported in different studies4-10. In ideal scenario, only two CEST-weighted (w) images at positive and negative offset frequency are required to compute CEST contrast. But in practical scenario, due to B0 inhomogeneity CEST-w image contains signal intensities from shifted offset frequencies rather than at desired frequency offset. Generally, CEST-w data is acquired at multiple offset-frequencies around frequency of interest for the purpose of B0 inhomogeneity corrections. For B0 inhomogeneity correction, pixel-wise CEST-w intensities obtained within this certain range are interpolated and value at offset-frequency of interest is replaced by the fitted/interpolated data value at adjusted offset-frequency corresponding to frequency shift in the B0 map. Some of the reported studies11,12 on CEST have used cubic interpolation for B0 inhomogeneity correction. Davis, et al13 reported spline interpolation, while Bagga P, et al14 reported polynomial function or linear fitting. Different studies11-14 used various interpolation functions to correct for B0 inhomogeneity. The objective of the current study is to evaluate the effect of interpolation methods to obtain accurate CEST asymmetry contrast using CEST-w data from simulations and in vivo human brain at 3T and 7T.

Methods

Simulations: The modified Bloch-McConnell equation based simulations were used to generate CEST data at 7T with saturation parameters root-mean-square B1(B1rms)=2.9µT, duration=800ms at frequency offsets -5ppm to +5ppm with step-size of 0.2ppm. B0 inhomogeneity of 0.1ppm was added to simulated dataset and GluCEST contrast was computed using asymmetry analysis as reported earlier8. Similarly, simulations were carried out to generate APT-w at B1rms=2µT, duration=800ms as well as CrCEST-w data at B1rms=1.4µT, duration=2000ms. Choice of these saturation parameters is based upon reported studies for corresponding CEST contrast. CEST contrast was computed using asymmetry analysis with normalization using control image M0 as well as Mneg.
Human brain GluCEST MRI at 7T: Six asymptomatic human volunteers(M:F=5:1;35±14.5years) were scanned at whole body 7T MRI scanner(Siemens). The study protocol included WASSR15, CEST and B1 map data acquired using pulse sequence reported previously7. CEST data were acquired at B1rms=2.9µT and duration=800ms at offset-frequencies ranging from ±1.8ppm to ±4.2ppm in steps of 0.1ppm. CEST data were also acquired at multiple B1rms of 0.7µT, 1.4µT and 2.2µT for B1 inhomogeneity correction.
Human brain APT-w MRI at 3T: A total of 24 treatment naive patients (14 low-grade-glioma (LGG); 10 high-grade-glioma (HGG)) were scanned at whole body 3T MRI scanner (Philips). A single slice APT-w images were acquired (positioned through center of lesion) with 64 frequency points (±0 ppm to ±14 ppm with step-size of 0.5 ppm, and control image at 100 ppm) using pulse sequence reported previously6. APT-w contrast corresponding to region-of-interest (ROI)-1 (entire tumor) and ROI-2 (contrast-enhancing tumor) were compared between LGG and HGG using Independent Student’s T-test with two tails. The difference was considered statistically significant if p<0.05.
CEST MRI data for each case was interpolated using various methods such as linear, polynomial (degree=1,2,3), spline, cubic spline and smoothing spline for B0 inhomogeneity correction. Results were compared using goodness-of-fit (R2).

Results

Figure-1 shows comparison of different interpolation methods for GluCEST-w signal intensity(SI) over the frequency range of +1.8ppm to +4.2ppm corresponding to ROI in WM. Accuracy of the linear interpolation(R2=0.98) was less compared to 2nd degree(R2=0.99) and 3rd degree(R2=0.99) polynomial interpolation. The cubic and spline interpolant fits to all the points with R2 =1, while smoothingspline(R2 =0.99) was similar to polynomial interpolant(degree=2, 3). Similar results were observed for APT-w and CrCEST-w SI. Figure 2 presents GluCESTM0-w and GluCESTNeg-w maps for different interpolations. Polynomial provides less noisy maps preserving the distribution of glutamate in WM and GM. Cubic and spline resulted in noisy and poor quality of maps. Smoothingspline provided less noisy maps compared to cubic and spline interpolation. The cubic and spline interpolation resulted in large variations in the values, as can be seen close to mid brain line on APT-w maps in Fig.3 b5 and b6. For the current study data, Polynomial function provided better quality of map compared with other types of interpolations. There was a significant difference(p<0.05) between APT-w contrast of LGG and HGG(Fig.4) corresponding to ROI-2 considering polynomial interpolations (p<0.01 (Fig.4b1); p=0.02 (Fig.4b2); p=0.02 (Fig.4b3)) for correcting B0 inhomogeneity. Table-1 show ROC analysis results and p values of APT-w for differentiation between LGG and HGG.

Discussion

Polynomial interpolation fits better to CEST MRI data over a limited range at 7T and 3T. With the further increase in degree of polynomial, interpolation gets better, but it becomes more prone to noise fluctuations. The cubic and spline interpolations are highly sensitive to noise. Polynomial(degree=2, 3) provided more accurate CEST contrast that is less sensitive to noise.

Conclusion

The current study showed that for in vivo CEST data(GluCEST-w, APT-w, CrCEST-w), polynomial(2nd & 3rd degree) interpolations provided better B0 inhomogeneity corrected CEST maps at 7T and polynomial interpolation, particularly linear showed improved differentiation of HGG and LGG tumors using APT-w MRI at 3T.

Acknowledgements

This study was supported by Indian Institute of Technology Delhi, Fortis Memorial Research Institute Gurugram and University of Pennsylvania. This study was partially supported by funding support from MATRICS, SERB-DST Grant number: MTR_2017_001021. The authors thank Dr. Jinyuan Zhou and Dr. Peter C.M. van Zijl for APT-w pulse sequence. The Authors thanks Dr. Hari Hariharan for CEST pulse sequence for 7T data acquisition.

References

1. Ward KM, Aletras AH, Balaban RS. A New Class of Contrast Agents for MRI Based on Proton Chemical Exchange Dependent Saturation Transfer (CEST). J Magn Reson. 2000;143(1):79-87.

2. Sun PZ, Wang E, Cheung JS, et al Simulation and optimization of pulsed radio frequency irradiation scheme for chemical exchange saturation transfer (CEST) MRI—demonstration of pH‐weighted pulsed‐amide proton CEST MRI in an animal model of acute cerebral ischemia. J Magn Reson. 2018; doi:10.1002/mrm.22894

3. Jones KM, Pollard AC, Pagel MD. Clinical Applications of Chemical Exchange Saturation Transfer (CEST) MRI. J Magn Reson. 2018; 47(1): 11-27.

4. Zhou J, Payen JF, Wilson DA, et al. Using the amide proton signals of intracellular proteins and peptides to detect pH effects in MRI. Nat Med. 2003;9:1085-1090.

5. Jones CK, Schlosser MJ, et al. Amide proton transfer imaging of human brain tumors at 3T. Magn Reson Med. 2006;56(3):585-592.

6. Debnath A, Gupta RK, Singh A. Evaluating the Role of Amide-Proton-Transfer Weighted contrast, Optimized for Normalization and Regions of Interest Selection, in Differentiation of Neoplastic and Infective Mass Lesions on 3T MRI. Molecular Imaging and Biology. 2019 May. DOI: 10.1007/s11307-019-01382-x

7. Cai K, Haris M, Singh A, et al. Magnetic resonance imaging of glutamate. Nat Med. 2012;18(2):302-306.

8. Debnath A, Gupta RK, Singh A. To evaluate the role of APT-w contrast, optimized for normalization and ROIs selection, in differentiating Infective and Neo-plastic Mass Lesions. International Society for Magnetic Resonance in Medicine (ISMRM): 27th Annual Meeting & Exhibition-May 2019, Canada, (Proc.Intl.Soc.Mag.Reson. 27 (2019), program number 4347)).

9. Haris M, Singh A, Cai K, et al. A technique for in vivo mapping of myocardial creatine kinase metabolism. Nat Med. 2014;20(2):209-214. doi:10.1038/nm.3436

10. Singh A, Debnath A, Cai K, et al. Evaluating Feasibility of Creatine-weighted CEST MRI in Human Brain at 7T using Z-spectral fitting approach. NMR in Biomedicine. 2019 August. DOI: 10.1002/nbm.4176

11. Neal A, Moffat BA, Stein JM, et al. Glutamate weighted imaging contrast in gliomas with 7 Tesla magnetic resonance imaging. NeuroImage Clin. 2019;22:101694.

12. Pépin J, Francelle L, Carrillo-de Sauvage MA, et al. In vivo imaging of brain glutamate defects in a knock-in mouse model of Huntington’s disease. Neuroimage. 2016;139:53-64.

13. Davis KA, Nanga RPR, Das S, et al. Glutamate imaging (GluCEST) lateralizes epileptic foci in nonlesional temporal lobe epilepsy. Sci Transl Med. 2015;7(309):309ra161.

14. Bagga P, Pickup S, Crescenzi R, et al. In vivo GluCEST MRI: Reproducibility, background contribution and source of glutamate changes in the MPTP model of Parkinson’s disease. Sci Rep. 2018;8(1):1-9. 15. Kim M, Gillen J, Landman BA, Zhou J, et al. Water saturation shift referencing (WASSR) for chemical exchange saturation transfer (CEST) experiments. Magn Reson Med. 2009;61(6):1441-1450.

Figures

Figure-1: Comparison of different interpolation methods for CEST-weighted signal intensity (SI) over the frequency range of +1.8 ppm to +4.2 ppm. Fitting of SI corresponding to an ROI in white matter (a1-a6) with step-size of 0.2 ppm. The columns 1-3 show interpolation methods using polynomial degree 1, 2 and 3 respectively. The columns 4-6 show interpolation methods using cubic, spline and smoothingspline respectively. R2 represents goodness-of-fit.

Figure-2: A representative case of B0 and B1 corrected GluCESTM0 maps (row-1) and GluCESTNeg maps (row-2) of a healthy volunteer for different interpolation methods. The columns 1-3 show interpolation methods using polynomial degree 1, 2 and 3 respectively. The columns 4-6 show interpolation methods using cubic, spline and smoothingspline respectively.

Figure-3: B0 inhomogeneity corrected APT-w maps of representative low-grade glioma (row-1) and high-grade-glioma(row-2) using different interpolation methods. The red arrow represents tumor regions. Since cubic and spline functions are based on piecewise polynomial functions, thus cubic and spline functions are showing large variations near mid brain line. Polynomial interpolation provided better APT-w maps.

Figure-4: Quantitative comparison using box and whisker plots of contrast normalized with APT-w value in normal-appearing-white-matter between low-grade-glioma and high-grade-glioma corresponding to ROI-1 and ROI-2. Polynomial interpolations provide significant difference (p<0.05) between the groups.

Table-1: show ROC analysis results of APT-w for differentiation between LGG and HGG.

Proc. Intl. Soc. Mag. Reson. Med. 28 (2020)
3103