Elsevier

NeuroImage

Volume 154, 1 July 2017, Pages 59-80
NeuroImage

Multi-echo fMRI: A review of applications in fMRI denoising and analysis of BOLD signals

https://doi.org/10.1016/j.neuroimage.2017.03.033Get rights and content

Highlights

  • An emerging fMRI technology, multi-echo (ME)-fMRI, makes fMRI more interpretable by enabling the detection of BOLD or artifact origins of acquired signals.

  • This review covers recent multi-echo fMRI acquisition methods, and the analysis steps for this data that make fMRI at once more principled, straightforward, and powerful.

  • After an overview of history and theory of ME-fMRI, T2* modeling and applications is discussed.

  • A focus is on the use of ME-fMRI data to extract and classify components from spatial ICA, called multi-echo ICA (ME-ICA).

  • New imaging methods such as multi-band multi-echo fMRI and imaging at 7 T are demonstrated throughout the review, and a practical analysis pipeline is described.

Abstract

In recent years the field of fMRI research has enjoyed expanded technical abilities related to resolution, as well as use across many fields of brain research. At the same time, the field has also dealt with uncertainty related to many known and unknown effects of artifact in fMRI data. In this review we discuss an emerging fMRI technology, called multi-echo (ME)-fMRI, which focuses on improving the fidelity and interpretability of fMRI. Where the essential problem of standard single-echo fMRI is the indeterminacy of sources of signals, whether BOLD or artifact, this is not the case for ME-fMRI. By acquiring multiple echo images per slice, the ME approach allows T2* decay to be modeled at every voxel at every time point. Since BOLD signals arise by changes in T2* over time, an fMRI experiment sampling the T2* signal decay can be analyzed to distinguish BOLD from artifact signal constituents. While the ME approach has a long history of use in theoretical and validation studies, modern MRI systems enable whole-brain multi-echo fMRI at high resolution. This review covers recent multi-echo fMRI acquisition methods, and the analysis steps for this data to make fMRI at once more principled, straightforward, and powerful. After a brief overview of history and theory, T2* modeling and applications will be discussed. These applications include T2* mapping and combining echoes from ME data to increase BOLD contrast and mitigate dropout artifacts. Next, the modeling of fMRI signal changes to detect signal origins in BOLD-related T2* versus artifact-related S0 changes will be reviewed. A focus is on the use of ME-fMRI data to extract and classify components from spatial ICA, called multi-echo ICA (ME-ICA). After describing how ME-fMRI and ME-ICA lead to a general model for analysis of fMRI signals, applications in animal and human imaging will be discussed. Applications include removing motion artifacts in resting state data at subject and group level. New imaging methods such as multi-band multi-echo fMRI and imaging at 7 T are demonstrated throughout the review, and a practical analysis pipeline is described. The review culminates with evidence from recent studies of major boosts in statistical power from using multi-echo fMRI for detecting activation and connectivity in healthy individuals and patients with neuropsychiatric disease. In conclusion, the review shows evidence that the multi-echo approach expands the range of experiments that is practicable using fMRI. These findings suggest a compelling future role of the multi-echo approach in subject-level and clinical fMRI.

Introduction

Blood oxygenation level dependent (BOLD) functional MRI (fMRI) is widely used to study brain activity based on hemodynamic signals (Bandettini et al., 1992, Bullmore et al., 1996, Buxton et al., 1998, Friston et al., 1995, Ogawa et al., 1990). However, recent studies show that fMRI data can be severely affected by artifacts (Power et al., 2012). These artifacts relate to subject head motion, cardiac and respiratory effects, and hardware (Glover et al., 2000, Jo et al., 2010). Studies on the effects of fMRI artifacts have brought into question many of the compelling findings on brain function based on fMRI, for example, relating to human brain development (Fair et al., 2009). fMRI artifacts also reduce the statistical power of fMRI studies and lead to spurious findings, which has been associated with a crisis of confidence in fMRI research (Button et al., 2013). Thus, while in recent years the field of fMRI research has both enjoyed advanced technology and expanded use, it also deals with a deep discomfort related to many known and unknown effects of artifact. In this review, we discuss an emerging fMRI approach, called multi-echo (ME)-fMRI, which focuses on improving the fidelity of fMRI signals through a physically-driven determination of the origins of fMRI signals as arising from either BOLD contrast or artifact. We also discuss how ME methods can be combined with emerging multi-band acceleration methods to create fMRI strategies with both high resolution and fidelity. Altogether, it is shown that ME-fMRI involves a few relatively small changes from standard fMRI acquisition technique, but enables analyses that support major improvements in fMRI data quality.

To date, the challenges in controlling fMRI artifacts have been met with generic time series signal processing methods such as regression and frequency-restricting bandpass filters (Satterthwaite et al., 2012). Head motion artifacts, modeled from shifts in brain images over time, are regressed out of fMRI time series. Recordings of physiology during MRI are used to model cardiac and respiratory artifact regressors (Glover et al., 2000). Spontaneous activity is band pass filtered to retain a narrow range of frequencies in order to exclude hardware related signal drifts and high-frequency noise (Carp, 2013). Time series models of noise and usually their temporal derivatives are linearly regressed out of data. Despite the significant reduction in information in data after applying these steps, it is now clear that much artifact remains. Of late, simple deletion of volumes from fMRI datasets has been suggested (Power et al., 2012). This altogether means that artifact signals are not characterized well enough by modeling artifact time courses and regressing them out of data. This reality also suggests that artifacts are themselves various and complex, and may interact in unpredictable ways. At the heart of the issue, however, is that standard fMRI approaches do not have a strong and general ground truth to precisely relate signals to biophysical signal mechanisms versus artifacts. Given information on how signals scale across the echo images of an ME-fMRI experiment, however, valuable insight on the origins of fMRI signals in BOLD contrast or artifact can be gained.

Section snippets

ME-fMRI and fMRI relaxometry

After excitation, standard fMRI uses 2-D echo planar imaging to acquire slice images at a single TE, one slice at a time. At 3 T, this TE is usually 30 ms. ME-fMRI uses a slightly different approach. After a normal excitation pulse, a slice image is acquired at the earliest TE possible. Without exciting again, readout of another image of the same slice is then acquired immediately afterward, at a longer TE, and so forth up to the desired number of images and TEs. This happens for each slice of

TE-dependence and TE-independence of signal changes in ME-fMRI data

Managing the contribution of various artifacts such as from subject head motion to fMRI time series is the central problem of fMRI denoising and analysis. When using single-echo fMRI, the challenge is in distinguishing neurally-related activity from artifacts. However, the scaling of fMRI signals across TEs, based on ME-fMRI data, can be analyzed to infer origins in neurally-related versus artifactual signal processes.

The amplitude of a signal fluctuation from standard single-echo fMRI data is

Classifying statistical components of fMRI datasets as problem and opportunity

TE-dependence and TE-independence models can be applied to classify signal components of fMRI that are found using spatial independent component analysis (ICA) (Kundu et al., 2012). ICA is a powerful method that can factor fMRI image time series into a set of linearly separable component maps and their corresponding time courses, expressed asX(v,t)=A(v,c)w(c,t)where X is a voxel-by-time matrix of time series, A is a voxel-by-component matrix expressing spatial maps, and W is a

A general model for fMRI time series

ME-ICA showed how brain network components could be grouped apart from artifacts of many kinds without spatial or temporal templates of expected functional effects. While ME-ICA was shown first for human resting-state fMRI data from 3 T MRI, it is in principle more broadly applicable. To denoise fMRI, TE-independent component time series could be regressed out of fMRI time series (McKeown et al., 2005). This approach can be used to denoise ME-fMRI of both task and resting-state activity since

Motion artifacts

While motion artifacts have long been known to affect fMRI, only recently has the severity of their impact on studies of connectivity become widely accepted (Power et al., 2012, Satterthwaite et al., 2012). There remains a vital need for an fMRI denoising method that can remove these artifacts in a principled way. Ideally, such a denoising method would lead to improvements in both subject and group levels of study. At the subject level, motion artifacts can either inflate or reduce

ME-fMRI using MB EPI and 7 T MRI

A perceived cost of the ME approach relates to the investment of imaging time to acquire more echoes, instead of using imaging time to acquire at high resolution. On the one hand, high spatial resolution is important for localizing brain function precisely. On the other hand, smaller voxel size leads to lower per-voxel signal to noise ratio, lowering the odds of detecting weaker signals, such as in deep brain areas, despite high spatial resolution. In regards to temporal resolution, acquiring

ME-fMRI of non-normative neuroanatomy

The use of ME-ICA at 7 T could make high-resolution fMRI studies possible for anatomies that are non-standard since ME-ICA needs no comparisons to standard anatomy to identify networks and artifacts. Fig. 17 demonstrates MEMB-fMRI with 2.5 mm isotropic voxel size with 4 TEs using 7 T MRI acquired from controls and epilepsy patients, analyzed with ME-ICA. For a 10-minute resting state dataset of a normal volunteer, representative BOLD networks are shown. Subcortical activity readily sub-divides

Making ME-fMRI practical

It might be surprising that with all the information to be gained from ME-fMRI over standard fMRI, ME-fMRI is still not standard. For example, the spiral and echo planar (SPEP) sequence by Wong et al., available for several years for the GE platform, could flexibly implement ME-fMRI as well as other techniques such as echo relaxation imaging. The main practical problem has been that gradient echo image readout takes too much time without in-plane acceleration so that by even a third TE most

Conclusion

Over the history of the development of fMRI methods, ME-fMRI has been important in the process of validating BOLD-related origins of novel fMRI signal observations, such as resting state time series correlation. The need for validation of novel BOLD phenomena is due to the fact that only a fraction of variance captured in a typical fMRI experiment is related to BOLD contrast. Bandettini et al. and Peltier and Noll, separated by a few years, used highly targeted ME-fMRI applications to conduct

Acknowledgements

PK and PB are supported by the Icahn School of Medicine Capital Campaign, the Translational and Molecular Imaging Institute, Brain Imaging Center and the Department of Radiology at the Icahn School of Medicine. PB is supported by funding from NIH-NCI (R01 CA202911) and NIH-NIMH (R01 MH109544) and Siemens Healthcare. VV is a Wellcome Trust Intermediate Fellow in Clinical Neurosciences (093705/Z/10/Z). We would like to thank Drs. Lara V. Marcuse and Madeline Fields at the Mount Sinai Epilepsy

References (67)

  • J.W. Evans et al.

    Separating slow BOLD from non-BOLD baseline drifts using multi-echo fMRI

    NeuroImage

    (2015)
  • K.J. Friston et al.

    Analysis of fMRI time-series revisited

    Neuroimage

    (1995)
  • J. Gonzalez-Castillo et al.

    Evaluation of multi-echo ICA denoising for task based fMRI studies: block designs, rapid event-related designs, and cardiac-gated fMRI

    NeuroImage

    (2016)
  • A. Ing et al.

    A dual echo approach to motion correction for functional connectivity studies

    NeuroImage

    (2012)
  • H.J. Jo et al.

    Mapping sources of correlation in resting state FMRI, with artifact detection and removal

    NeuroImage

    (2010)
  • E. Kirilina et al.

    The quest for the best: the impact of different EPI sequences on the sensitivity of random effect fMRI group analyses

    NeuroImage

    (2016)
  • P.J. Koopmans et al.

    Multi-echo fMRI of the cortical laminae in humans at 7 T

    NeuroImage

    (2011)
  • P. Kundu et al.

    Differentiating BOLD and non-BOLD signals in fMRI time series using multi-echo EPI

    NeuroImage

    (2012)
  • P. Kundu et al.

    Differentiating BOLD and non-BOLD signals in fMRI time series from anesthetized rats using multi-echo EPI at 11.7 T

    NeuroImage

    (2014)
  • L.S. Morris et al.

    Jumping the gun: mapping neural correlates of waiting impulsivity and relevance across alcohol misuse

    Biol. Psychiatry

    (2016)
  • L.S. Morris et al.

    Fronto-striatal organization: defining functional and microstructural substrates of behavioural flexibility

    Cortex J. Devoted Study Nerv. Syst. Behav.

    (2016)
  • K. Murphy et al.

    The impact of global signal regression on resting state correlations: are anti-correlated networks introduced?

    NeuroImage

    (2009)
  • V. Olafsson et al.

    Enhanced identification of BOLD-like components with multi-echo simultaneous multi-slice (MESMS) fMRI and multi-echo ICA

    NeuroImage

    (2015)
  • S.J. Peltier et al.

    Analysis of fMRI signal and noise component TE dependence

    Neuroimage

    (2000)
  • J.R. Polimeni et al.

    Laminar analysis of 7 T BOLD using an imposed spatial activation pattern in human V1

    NeuroImage

    (2010)
  • J.D. Power et al.

    Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion

    Neuroimage

    (2012)
  • Z.S. Saad et al.

    A new method for improving functional-to-structural MRI alignment using local Pearson correlation

    Neuroimage

    (2009)
  • G. Salimi-Khorshidi et al.

    Automatic denoising of functional MRI data: combining independent component analysis and hierarchical fusion of classifiers

    NeuroImage

    (2014)
  • S. Smith et al.

    Advances in functional and structural MR image analysis and implementation as FSL

    Neuroimage

    (2004)
  • C.G. Thomas et al.

    Noise reduction in BOLD-based fMRI using component analysis

    Neuroimage

    (2002)
  • D.C. Van Essen et al.

    The human connectome project: a data acquisition perspective

    NeuroImage

    (2012)
  • F. Zhao et al.

    Cortical layer-dependent BOLD and CBV responses measured by spin-echo and gradient-echo fMRI: insights into hemodynamic regulation

    Neuroimage

    (2006)
  • P. Bandettini et al.

    Time course EPI of human brain function during task activation

    Magn. Reson. Med.

    (1992)
  • Cited by (0)

    View full text