Jed Wingrove1, Stephen J Wastling1, Prantik Kundu2, Gareth Barker1, Donal Hill1, Owen O'Daly1, and Fernando Zelaya1
1Department of Neuroimaging, King's College London, London, United Kingdom, 2Brain Imaging Center, Mount Sinai Hospital, New York, NY, United States
Synopsis
Resting state fMRI data is highly susceptible to low frequency noise fluctuations from motion and pulsatile physiological movement leading to inaccuracies in observed connectivity. In this study we directly compare conventional single echo EPI against a multi echo EPI (ME-EPI) acquisition and analysis method which denoises and cleans acquired fMRI time series. Task evoked and rs-fMRI data was acquired in 8 volunteers and showed improved spatial localisation and discrete clustering in ME-EPI analysis in comparison to single echo data. Furthermore, this work demonstrates the benefits of using ME-EPI for both functional activation and resting state connectivity investigations. Introduction
Resting state fMRI (rs-fMRI) enables the investigation of functional
connectivity, identifying integrated brain networks at rest (RSNs), as well as revealing
specific functional connections to a pre-defined region (‘seed-based’ analysis)
1.
However, rs-fMRI is highly affected by artefactual low frequency components, originating
from head motion, scanner drift and physiological pulsatile motion; leading to
activation error and overestimation. Multi-Echo-EPI (ME-EPI) acquisition,
combined with independent component (ICA) and echo time (TE) dependent
2, is an emerging methodology to de-noise and remove such artefacts
3.
In this investigation we assess the advantages of ME-EPI over conventional
single echo protocols, by acquiring task-evoked and resting functional
information in the same scanning session in 8 healthy volunteers, with both
protocols.
Methods
The study was approved by the KCL human research ethics
committee. All scans were acquired in a 3T
GE-MR750 scanner in 8 healthy volunteers (2 female). rs-fMRI and motor-task fMRI data were
collected in all scans, with single echo and ME-EPI sequences. The order was counterbalanced across subjects,
but both rs-fMRI scans were collected first to avoid carry over effects of
activation.
For all fMRI acquisitions (denoted SETASK, METASK,
SERS and MERS), identical parameters were used: FA = 80°,
TR = 2500ms, TE = 30ms (single echo) and TE = 12, 28, 44, 60ms, voxel-size:
3.3x3.3x4mm. A single, 1mm isotropic T1-weighted MPRAGE images was also
acquired.
Motor Task
Participants were required to squeeze and release two soft
balls in each hand upon a visual cue; once every 2.5s for 40s periods,
alternating with 40s of rest (5 repeats), total acquisition time = 6:40 mins,
160 images.
Resting State
Participants were scanned in the supine position whilst
observing a cross-hair on a dark screen for a total of 8.00mins, 192 images.
Image Processing
Motion correction and DARTEL-based spatial normalisation
were performed with SPM-12, and a ME-ICA de-noising toolkit2. Motion
parameter regressors were used in the first level analysis. A random-effects
analysis was employed to derive group activation maps for the motor task and a paired
t-test was conducted across sequences.
Resting state data required further pre-processing; motion
regression, regression of WM and CSF signal, data de-trending and band-pass filtering
(>0.1Hz) prior to seed based connectivity analysis. After TE dependent analysis
(which removes non-exponential decay of the ME signal) and ICA, a de-noised,
optimally combined time series using all four echoes was obtained3. Two 6mm3 spherical seeds, identified from peak significance in the SETASK
group analysis, were used to produce functional connectivity maps for both SERS
and MERS.
Results
Both
ME-EPI and single echo data from the motor task revealed significant areas (pvoxel<0.001,
pcluster < 0.05) of
activation in bilateral primary motor cortices, supplementary motor area (SMA),
thalamus and cerebellum. However (as clearly shown in Fig. 1), ME-EPI derived
activation was better localised and constrained in comparison to the single
echo data.
One-sample t-test of the seed based rs-fMRI, using a seed in
the SMA (obtained from the SETASK activation map) revealed significant (pvoxel
< 0.001, pcluster < 0.05) connectivity with the bilateral
primary motor areas, visual cortex and superior cerebellum, using both
methodologies. But significant connectivity with the insula cortex and the
brainstem were only identified with the ME-EPI analysis (Seen from Fig. 2).
Discussion
This is (to our knowledge), one of the first within-subject,
within-session evaluations of conventional single echo fMRI acquisition against
a combined ME-EPI sequence. We employed
ME-EPI with TE-dependent and ICA analysis, to remove non-BOLD signal components
in both task and resting state fMRI. We demonstrate that in both instances ME-EPI
produces more spatially defined regions of activation vs single echo, and
better represents the underlying functional anatomy from which the signal
originates. We recommend the use of ME-EPI over single echo protocols for both task evoked and resting state functional studies.
Acknowledgements
No acknowledgement found.References
1. S,M. Smith, D. Vidaurre, C,F. Beckmann, M,F. Glasser, M. Jenkinson, K,L. Miller, T,E. Nichols, E,C. Robinson, G. Salimi-Khorshidi, M,W. Woolrich, D,M. Barch, K. Ugurbil and D,C. Van Essen. Functional Connectomics from Resting-State fMRI. Cell. 2013;17(12):666-682.
2. P. Kundu, S,J Inati, J,W. Evans, W. Luh, P,A. Bandettini. Differentiating BOLD and non-BOLD Signals in fMRI Time Series Using Multi-Echo EPI. Neuroimage. 2011;60(3): 1759-1770.
3. P. Kundu, N,D. Brenowitz, V. Voon, et al. Integrated Strategy for Improving Functional Connectivity Mapping using Multiecho fMRI. PNAS. 2013;110(40):16187-16192.