Slice-acceleration Related Biases in Multiband-EPI Resting State Functional Connectivity
Zahra Faraji-Dana1,2, Ali Golestani3, Yasha Khatamian3, Simon Graham1,2, and J. Jean Chen1,3

1Department of Medical Biophysics, University of Toronto, Toronto, ON, Canada, 2Sunnybrook Research Institute, Sunnybrook Health Science Centre, Toronto, ON, Canada, 3Rotman Research Institute, Baycrest Health Science Centre, Toronto, ON, Canada

Synopsis

Resting-state functional connectivity MRI (rs-fcMRI) is most commonly computed as the temporal dependency amongst blood oxygenation Level dependent (BOLD) signal patterns of different brain regions. It has been shown that rs-fMRI can benefit from faster imaging times offered by simultaneous multi-slice (a.k.a. multiband, referred to as “MB”) slice-acceleration that enables acquiring “groups” of slices at the same time. However, this slice grouping may incur aliasing artifacts, primarily from motion and physiological fluctuations. These spurious time-dependent signals can adversely affect the rs-fcMRI maps in the simultaneously-acquired slices (i.e., in one slice-group). In this work we investigate two hypotheses 1) the simultaneously sampled physiological noises as well as the residual aliasing introduce a slice-group effect in rs-fcMRI maps; 2) this slice-group effect can be mitigated by physiological noise correction.

Introduction

Resting-state functional connectivity MRI (rs-fcMRI) commonly measures brain network activity in the absence of a behavioural task using blood oxygenation level dependent (BOLD) contrast. Typically computed from the temporal dependency of BOLD signal patterns across different brain regions, rs-fcMRI has been used extensively to investigate various neuroscience and clinical questions 1. It has been shown that rs-fcfMRI can benefit from faster imaging times offered by simultaneous multi-slice (a.k.a. multiband, or “MB”) slice-acceleration 2,3 that enables “groups” of slices to be acquired at the same time. However, this slice-grouping may incur aliasing artifacts, primarily from motion 4 and “physiological noise”. The latter spurious fluctuations arise from oscillatory respiration and cardiac signals, as well as respiration volume per time (RVT) 5 and cardiac rate variation (CRV) 6. It is possible that aliasing of physiological noise may introduce unwanted signal components across slices acquired within a given group, with a deleterious effect on functional conductivity estimates. Consequently, two hypotheses are investigated in this work: 1) physiological noise and associated aliasing effects introduce a slice-group effect in rs-fcMRI maps; and 2) this slice-group effect can be mitigated by physiological noise correction by retrospective gating.

Method

All images were acquired using a Siemens TIM Trio 3 T MRI System (Siemens, Erlangen, Germany) employing the 32-channel phased-array head coil receiver. Twelve subjects were each scanned with a functional single-shot MB gradient-echo (GE) echo-planar imaging (EPI) 1 sequence (TR=323ms, TE=30 ms, flip angle=40°, 15 slices, 3.44 mm × 3.44 mm × 4.6 mm, 2230 volumes, acceleration factor=3, phase encoding shift factor=2), and a T1-weighted anatomical sequence (MPRAGE, TR = 2400 ms, TE = 2.43 ms, FOV = 256 mm, TI = 1000 ms, readout bandwidth = 180 Hz/px, voxel size = 1 mm × 1 mm × 1 mm). In the MB-EPI scan, the 15 slices were grouped as 1-6-11, 2-7-12, 3-8-13, 4-9-14, 5-10-15, and acquired in this order. For each subject, a pneumatic belt (BioPac, Goleta, USA) and the vendor supplied pulse oximeter were used to record the respiratory and cardiac signals, respectively. Subjects were instructed to close their eyes and relax. The rs-fcfMRI processing was carried out using the FMRIB Software Library v5 (FSL5) (Analysis Group, FMRIB, Oxford, UK) 7.

The preprocessing pipeline included motion correction, brain extraction, spatial smoothing, and regression of six motion parameters. Five different physiological noise correction schemes were performed on the data of each subject using MATLAB (MathWorks, Natick, MA): a) no correction; b) image-based retrospective correction of cardiac and respiratory waveforms (RETROICOR) 8; c) RETROICOR+RVTcor; d) RETROICOR+CRVcor; and e) RETROICOR+RVTcor+CRVcor, where RVTcor and CRVcor refer to the inclusion of these physiological noise signals in the retrospective correction algorithm. Thus, five post-processed fMRI datasets were considered using a multi-regression model, focusing on rs-fcMRI maps of the default mode network (DMN), as it has extensive brain coverage. A 4 mm-radius spherical seed was generated over the posterior cingulate cortex (PCC) using the coordinates reported by reported by Van Dijk et al 9. For each of the five corrected datasets, the spatial average of the signal from the seed was correlated with all other voxels to generate the connectivity maps. All statistical connectivity maps were then corrected using mixture modelling 10 and the resulting Z-maps were then thresholded at Z>2.3 (p<0.01). Comparing the rs-fcMRI maps to an appropriate connectivity template for the DMN 11, the number of false positive and false negative activations were separately obtained for each slice, each physiological noise correction scheme, and each subject. For each subject and correction scheme, the mean number of false positives and false negatives per slice was then determined. Two repeated-measures analyses of variance (rANOVAs) were subsequently performed using MATLAB to compare the effect of the correction schemes on the mean number of false positives and false negatives. Using the known MB-EPI slice-grouping, the slice-group effect on the number of false positives and false negatives of each slice was independently studied via linear mixed effect (LME) analysis in R (R Foundation for Statistical Computing, Vienna, Austria).

Results

Fig. 1 illustrates the rs-fcMRI maps for a representative subject, corresponding to the five physiological noise correction schemes (a-e), and the DMN template (f) 11. Using the template as reference, the number of false positives gradually reduces in progression from correction scheme a) to c), with little difference observed between schemes (c-e). Similar results were observed across the cohort, except that the number of false positives continues to decrease across all five correction schemes. The rANOVA results showed a strong trend of physiological correction reducing the number of false positives, that approached statistical significance (p=0.09). No effects were found for false negatives (p=0.81).

Table 1 shows the p-values of the slice-group LME modeling for different physiological noise correction schemes. According to the LME model without physiological noise correction, the number of false positives demonstrates a significant slice-group effect, whereas the number of false negatives does not. Additionally, the slice-group effect (both false positives and false negatives) incrementally diminishes (i.e. p-value increases) progressing from correction scheme a) to scheme e). Specifically, for false positives, RETROICOR was effective at suppressing the slice-group effect so that it was no longer statistically significant.

Discussion and Conclusion

This work shows that in the absence of physiological noise correction, a significant slice-group effect is present in the number of false positives in MB-EPI rs-fcMRI maps. Due to the type of template comparison that was employed, the number of false negatives was large and did not show a significant slice-group effect, although a trend was observed. The work also suggests that slice-group biases can be significantly reduced by correcting for physiological effects in MB-EPI data. In this regard, appropriate physiological correction is potentially even more important for MB-EPI than for regular EPI acquisitions.

Acknowledgements

No acknowledgement found.

References

1. Van den Heuvel, M. P. & Hulshoff Pol, H. E. Exploring the brain network: a review on resting-state fMRI functional connectivity. Eur. Neuropsychopharmacol. 20, 519–34 (2010).

2. Setsompop, K. et al. Blipped-controlled aliasing in parallel imaging for simultaneous multislice echo planar imaging with reduced g-factor penalty. Magn. Reson. Med. 67, 1210–24 (2012).

3. Feinberg, D. a et al. Multiplexed echo planar imaging for sub-second whole brain FMRI and fast diffusion imaging. PLoS One 5, e15710 (2010).

4. Bhat, H. et al. Motion Insensitive ACS Acquisition Method for in-plane and Simultaneous Multi-Slice Accelerated EPI. Proc. Intl. Soc. Mag. Reson. Med 22, 7255 (2014).

5. Birn, R. M., Diamond, J. B., Smith, M. a & Bandettini, P. a. Separating respiratory-variation-related fluctuations from neuronal-activity-related fluctuations in fMRI. Neuroimage 31, 1536–48 (2006).

6. Chang, C., Cunningham, J. P. & Glover, G. H. Influence of heart rate on the BOLD signal: the cardiac response function. Neuroimage 44, 857–69 (2009).

7. Smith, S. M. et al. Advances in functional and structural MR image analysis and implementation as FSL. Neuroimage 23 Suppl 1, S208–19 (2004).

8. Glover, G. H., Li, T. & Ress, D. Image-Based Method for Retrospective Correction of Physiological Motion Effects in fMRI: RETROICOR. Magn. Reson. Med. 167, 162–167 (2000).

9. Van Dijk, K. R. a et al. Intrinsic functional connectivity as a tool for human connectomics: theory, properties, and optimization. J. Neurophysiol. 103, 297–321 (2010).

10. Woolrich, M. W., Behrens, T. E. J., Beckmann, C. F. & Smith, S. M. Mixture models with adaptive spatial regularization for segmentation with an application to FMRI data. IEEE Trans. Med. Imaging 24, 1–11 (2005).

11. Yeo, B. T. T. et al. The organization of the human cerebral cortex estimated by intrinsic functional connectivity. J. Neurophysiol. 106, 1125–65 (2011).

Figures

Fig. 1: rs-fcMRI maps for a representative subject after applying a) no correction b) RETROICOR, c) RETROICOR+RVTcor, d) RETROICOR+CRVcor, e) RETROICOR+ RVTcor+CRVcor, (f) the DMN resting-state connectivity template in the subject coordinates. The arrows point at false positive activations.

Table1: p-values of the slice-group LME analysis after applying different physiological noise correction schemes (a-e).



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