Differences in slow drift among echoes in multiband multiecho EPI data compromise TE-dependent analysis
Kun Lu1, Javier Gonzalez-Castillo 2, Matthew Middione3, Brice Fernandez4, Valur Olafsson5, Prantik Kundu6, Ajit Shankaranarayanan7, and Thomas Liu1

1Center for Functional Magnetic Resonance Imaging, University of California, San Diego, La Jolla, CA, United States, 2Laboratory of Brain and Cognition, Section on Functional Imaging Methods, National Institutes of Health, Bethesda, MD, United States, 3Applied Sciences Laboratory West, GE Healthcare, Menlo Park, CA, United States, 4Applications and Workflow, GE Healthcare, Munich, Germany, 5Neuroscience Imaging Center, University of Pittsburgh, Pittsburgh, PA, United States, 6Brain Imaging Center, Icahn Institute of Medicine at Mt. Sinai, New York, NY, United States, 7Applications and Workflow, GE Healthcare, Menlo Park, CA, United States

Synopsis

We recently observed differences in slow signal drift between multiple echoes in multiband multiecho EPI data. Such differences in drift compromise the TE dependence model and could impact the TE dependence analysis (e.g Multiecho Independent Component Analysis ME-ICA). We first designed a metric MEQA to quantitatively measure the observed drift differences, then investigated the effects of the drift differences on ME-ICA analysis.

Introduction

Multiband Multiecho (MBME) EPI has the advantage of higher temporal resolution and improved sensitivity to BOLD signal compared to conventional single band single echo EPI. A few studies have examined the benefits of MBME in resting state and task fMRI [1,2]. It has also been demonstrated that TE-dependent analysis methods (e.g Multiecho Independent Component Analysis ME-ICA [3]) can be applied to MBME data to effectively remove non-BOLD signals [1]. However, we recently observed differences in the slow signal drift between multiple echoes in MBME data, which could compromise the TE-dependence model and impact the ME-ICA analysis. Figure 1C. shows an example of such drift differences in an MBME data acquired on an agar phantom. For better illustration, we show data from an acquisition with a significant amount of aliasing artifact in the reconstructed images (Figure 1A). We have noted that areas of aliasing artifacts often exhibit large drift differences. Our goal is to design a metric to quantify the differences in drift among multiple echoes, and to probe the effects of the drift differences on TE dependent analysis such as ME-ICA.

Method

The previous MBME studies have used in-plane acceleration for faster data acquisition and reduced off resonance effects. Gomez et al [4] also indicated that in-plane acceleration combined with lower multiband factor (2 or 3) increases BOLD temporal SNR. In this study, we acquired both phantom and in-vivo MBME data with and without in-plane acceleration (Auto-calibrating Reconstruction for Cartesian Imaging: ARC) so we can compare the effects of in-plane acceleration on the signal drift. We used a GE 3T Discovery MR750 (GE Healthcare, Waukesha, WI) and a prototype MBME pulse sequence with blipped CAIPI [5]. The parameters were: MB=3, 250 time points, 3 echoes, 24cm FOV, and 3.75x3.75x3.8mm voxel size. For data without ARC: TR = 1.2s, TE = 11, 32.2, and 53.4ms; and with ARC=2: TR = 1s, TE = 11, 25, and 39ms. An offline Matlab recon was used for image reconstruction. To quantify the drift differences, we first performed a fit of the MBME data to obtain T2* and S0 time courses (Figure 2). The residual time course from this fit was then modeled voxel-wise as a 10th order Legendre polynomial to remove high frequency noise components that did not contribute to slow drifts. From this fitted residual, we calculated the standard deviation and used it as a quantitative measure of the drift difference (dubbed as MEQA, Figure 2.) For the in-vivo data, we also performed ME-ICA analysis.

Results

Figure 1B shows an example of the MEQA map. High MEQA values correspond well to voxels with remaining aliasing artifacts and large echo drift differences. Figure 3 shows calculated MEQA maps from the MBME phantom and in-vivo data. The data with in-plane acceleration (ARC=2) has significantly more MEQA voxels than the data without ARC, indicating more prevalent deviation from the TE dependence model in the ARC=2 data. This is true for both the phantom and in-vivo data.

The ME-ICA analysis shows that the data acquired without ARC has higher number of both total ICA components and BOLD-like components than the data with ARC (Table 1.)

Discussion

Our results show that the MEQA metric is a sensitive measure of slow drift differences among multiple echoes in MBME data. Such differences are typically observed in voxels with imaging artifacts due to imperfect image reconstruction, susceptibility, and motion (not shown). Therefore MEQA could potentially be used as a quality assurance metric for MBME data. We also show that data with higher MEQA voxel counts, and thus more prevalent echo drift differences, has lower ME-ICA performance. Although very preliminary, this result warrants further investigations of the source of the drift differences and how they affect ME-ICA analysis.

Acknowledgements

This work was supported in part by a research grant from GE Healthcare.

References

1. Olafsson, V., Kundu, P., Wong, E.C., Bandettini, P.A., Liu, T.T., Enhanced identification of BOLD-like components with multi-echo simultaneous multi-slice (MESMS) fMRI and multi-echo ICA. Neuroimage, 2015; 112: 43–51.

2. Boyacioglu, R., Schulz, J., Koopmans, P., Barth, M., Norris, D.G., Improving sensitivity and specificity for RS fMRI using multiband multi-echo EPI at 7 T. Neuroimage, 2015; 119: 352–361.

3. Kundu, P., Inati, S.J., Evans, J.W., Luh, W.-M., Bandettini, P. A, Differentiating BOLD and non-BOLD signals in fMRI time series using multi-echo EPI. Neuroimage. 2011; 60: 1759–1770.

4. Gomez, ED. , Schulz, J., Boyacioglu R., Norris, DG., Poser, BA., Multiband Multiecho 2D-EPI: Maximizing BOLD CNR for fMRI at 3T. Proc. Intl. Soc. Mag. Reson. Med 2015; 23: 2045.

5. Setsompop, K., Gagoski, B.A., Polimeni, J.R., Witzel, T., Wedeen, V.J., Wald, L.L., Blipped-controlled aliasing in parallel imaging for simultaneous multislice echo planar imaging with reduced g-factor penalty. Magn. Reson. Med. 2012; 67:1210–1224.

Figures

Figure 1. A: Raw MBME images of an agar phantom showing remaining aliasing artifacts. B: Calculated MEQA maps overlaid on top of the MBME image. C: Example time courses from selected voxels (in the green box) showing large differences in slow signal drift over time between the three echoes (black: echo 1; red: echo 2; blue: echo 3).

Figure 2. MEQA calculation flow diagram

Figure 3. Calculated MEQA maps from the phantom and in-vivo MBME data. Data with ARC =2 (right column) has significantly more voxels with large echo drift differences then data without ARC (left column).

Table 1. Results from the ME-ICA analysis



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