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 7T 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.
•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 7T are demonstrated throughout the review, and a practical analysis pipeline is described.
•Most rs-fMRI is done at 3T with T2*w-EPI of moderate spatiotemporal resolution.•7T rs-fMRI offers higher spatial resolution, for example to study cortical depths.•Multiband EPI and other ...undersampling schemes are used to increase the sampling rate.•Multi-echo EPI allows measurement of connectivity of deep gray matter structures.•Alternative contrasts include VASO-CBV and T2-weighting (for spatial specificity).•Another alternative contrast, CBF, is often used for connectivity measures in clinical populations.
Resting state functional magnetic resonance imaging (rs-fMRI) is based on spontaneous fluctuations in the blood oxygen level dependent (BOLD) signal, which occur simultaneously in different brain regions, without the subject performing an explicit task. The low-frequency oscillations of the rs-fMRI signal demonstrate an intrinsic spatiotemporal organization in the brain (brain networks) that may relate to the underlying neural activity. In this review article, we briefly describe the current acquisition techniques for rs-fMRI data, from the most common approaches for resting state acquisition strategies, to more recent investigations with dedicated hardware and ultra-high fields. Specific sequences that allow very fast acquisitions, or multiple echoes, are discussed next. We then consider how acquisition methods weighted towards specific parts of the BOLD signal, like the Cerebral Blood Flow (CBF) or Volume (CBV), can provide more spatially specific network information. These approaches are being developed alongside the commonly used BOLD-weighted acquisitions. Finally, specific applications of rs-fMRI to challenging regions such as the laminae in the neocortex, and the networks within the large areas of subcortical white matter regions are discussed. We finish the review with recommendations for acquisition strategies for a range of typical applications of resting state fMRI.
A central challenge in the fMRI based study of functional connectivity is distinguishing neuronally related signal fluctuations from the effects of motion, physiology, and other nuisance sources. ...Conventional techniques for removing nuisance effects include modeling of noise time courses based on external measurements followed by temporal filtering. These techniques have limited effectiveness. Previous studies have shown using multi-echo fMRI that neuronally related fluctuations are Blood Oxygen Level Dependent (BOLD) signals that can be characterized in terms of changes in R2* and initial signal intensity (S0) based on the analysis of echo-time (TE) dependence. We hypothesized that if TE-dependence could be used to differentiate BOLD and non-BOLD signals, non-BOLD signal could be removed to denoise data without conventional noise modeling. To test this hypothesis, whole brain multi-echo data were acquired at 3 TEs and decomposed with Independent Components Analysis (ICA) after spatially concatenating data across space and TE. Components were analyzed for the degree to which their signal changes fit models for R2* and S0 change, and summary scores were developed to characterize each component as BOLD-like or not BOLD-like. These scores clearly differentiated BOLD-like “functional network” components from non BOLD-like components related to motion, pulsatility, and other nuisance effects. Using non BOLD-like component time courses as noise regressors dramatically improved seed-based correlation mapping by reducing the effects of high and low frequency non-BOLD fluctuations. A comparison with seed-based correlation mapping using conventional noise regressors demonstrated the superiority of the proposed technique for both individual and group level seed-based connectivity analysis, especially in mapping subcortical–cortical connectivity. The differentiation of BOLD and non-BOLD components based on TE-dependence was highly robust, which allowed for the identification of BOLD-like components and the removal of non BOLD-like components to be implemented as a fully automated procedure.
Display omitted
► ICA components were analyzed for TE-dependence. ► Component-level TE-dependence scores identified BOLD components. ► BOLD networks were identified without spatial templates or frequency thresholds. ► Removing non-BOLD components dramatically denoised data. ► Robust subcortical-cortical connectivity was revealed after denoising.
Abstract Background A recent hypothesis has suggested that core deficits in goal-directed behavior in obsessive-compulsive disorder (OCD) are caused by impaired frontostriatal function. We tested ...this hypothesis in OCD patients and control subjects by relating measures of goal-directed planning and cognitive flexibility to underlying resting-state functional connectivity. Methods Multiecho resting-state acquisition, combined with micromovement correction by blood oxygen level–dependent sensitive independent component analysis, was used to obtain in vivo measures of functional connectivity in 44 OCD patients and 43 healthy comparison subjects. We measured cognitive flexibility (attentional set-shifting) and goal-directed performance (planning of sequential response sequences) by means of well-validated, standardized behavioral cognitive paradigms. Functional connectivity strength of striatal seed regions was related to cognitive flexibility and goal-directed performance. To gain insights into fundamental network alterations, graph theoretical models of brain networks were derived. Results Reduced functional connectivity between the caudate and the ventrolateral prefrontal cortex was selectively associated with reduced cognitive flexibility. In contrast, goal-directed performance was selectively related to reduced functional connectivity between the putamen and the dorsolateral prefrontal cortex in OCD patients, as well as to symptom severity. Whole-brain data-driven graph theoretical analysis disclosed that striatal regions constitute a cohesive module of the community structure of the functional connectome in OCD patients as nodes within the basal ganglia and cerebellum were more strongly connected to one another than in healthy control subjects. Conclusions These data extend major neuropsychological models of OCD by providing a direct link between intrinsically abnormal functional connectivity within dissociable frontostriatal circuits and those cognitive processes underlying OCD symptoms.
The sense of agency (SA) is an established framework that refers to our ability to exert and perceive control over our own actions. Having an intact SA provides the basis for the human perception of ...voluntariness, while impairments in SA are hypothesized to lead to the perception of movements being involuntary that may be seen many neurological or psychiatric disorders. Individuals with functional movement disorders (FMD) experience a lack of control over their movements, yet these movements appear voluntary by physiology. We used fMRI to explore whether alterations in SA in an FMD population could explain why these patients feel their movements are involuntary. We compared the FMD group to a control group that was previously collected using an ecologically valid, virtual-reality movement paradigm that could modulate SA. We found selective dysfunction of the SA neural network, whereby the dorsolateral prefrontal cortex and pre-supplementary motor area on the right did not respond differentially to the loss of movement control. These findings provide some of the strongest evidence to date for a physiological basis underlying these disabling disorders.
Functional connectivity analysis of resting state blood oxygen level–dependent (BOLD) functional MRI is widely used for noninvasively studying brain functional networks. Recent findings have ...indicated, however, that even small (≤1 mm) amounts of head movement during scanning can disproportionately bias connectivity estimates, despite various preprocessing efforts. Further complications for interregional connectivity estimation from time domain signals include the unaccounted reduction in BOLD degrees of freedom related to sensitivity losses from high subject motion. To address these issues, we describe an integrated strategy for data acquisition, denoising, and connectivity estimation. This strategy builds on our previously published technique combining data acquisition with multiecho (ME) echo planar imaging and analysis with spatial independent component analysis (ICA), called ME-ICA, which distinguishes BOLD (neuronal) and non-BOLD (artifactual) components based on linear echo-time dependence of signals—a characteristic property of BOLD $T_{2}^{*}$ signal changes. Here we show for 32 control subjects that this method provides a physically principled and nearly operator-independent way of removing complex artifacts such as motion from resting state data. We then describe a robust estimator of functional connectivity based on interregional correlation of BOLD-independent component coefficients. This estimator, called independent components regression, considerably simplifies statistical inference for functional connectivity because degrees of freedom equals the number of independent coefficients. Compared with traditional connectivity estimation methods, the proposed strategy results in fourfold improvements in signal-to-noise ratio, functional connectivity analysis with improved specificity, and valid statistical inference with nominal control of type 1 error in contrasts of connectivity between groups with different levels of subject motion.
Head motion can be estimated at any point of fMRI image processing. Processing steps involving temporal interpolation (e.g., slice time correction or outlier replacement) often precede motion ...estimation in the literature. From first principles it can be anticipated that temporal interpolation will alter head motion in a scan. Here we demonstrate this effect and its consequences in five large fMRI datasets. Estimated head motion was reduced by 10-50% or more following temporal interpolation, and reductions were often visible to the naked eye. Such reductions make the data seem to be of improved quality. Such reductions also degrade the sensitivity of analyses aimed at detecting motion-related artifact and can cause a dataset with artifact to falsely appear artifact-free. These reduced motion estimates will be particularly problematic for studies needing estimates of motion in time, such as studies of dynamics. Based on these findings, it is sensible to obtain motion estimates prior to any image processing (regardless of subsequent processing steps and the actual timing of motion correction procedures, which need not be changed). We also find that outlier replacement procedures change signals almost entirely during times of motion and therefore have notable similarities to motion-targeting censoring strategies (which withhold or replace signals entirely during times of motion).
Adolescent changes in human brain function are not entirely understood. Here, we used multiecho functional MRI (fMRI) to measure developmental change in functional connectivity (FC) of resting-state ...oscillations between pairs of 330 cortical regions and 16 subcortical regions in 298 healthy adolescents scanned 520 times. Participants were aged 14 to 26 y and were scanned on 1 to 3 occasions at least 6 mo apart. We found 2 distinct modes of age-related change in FC: “conservative” and “disruptive.” Conservative development was characteristic of primary cortex, which was strongly connected at 14 y and became even more connected in the period from 14 to 26 y. Disruptive development was characteristic of association cortex and subcortical regions, where connectivity was remodeled: connections that were weak at 14 y became stronger during adolescence, and connections that were strong at 14 y became weaker. These modes of development were quantified using the maturational index (MI), estimated as Spearman’s correlation between edgewise baseline FC (at 14 y, FC14) and adolescent change in FC (ΔFC14−26), at each region. Disruptive systems (with negative MI) were activated by social cognition and autobiographical memory tasks in prior fMRI data and significantly colocated with prior maps of aerobic glycolysis (AG), AG-related gene expression, postnatal cortical surface expansion, and adolescent shrinkage of cortical thickness. The presence of these 2 modes of development was robust to numerous sensitivity analyses. We conclude that human brain organization is disrupted during adolescence by remodeling of FC between association cortical and subcortical areas.
Artifact removal in resting state fMRI (rfMRI) data remains a serious challenge, with even subtle head motion undermining reliability and reproducibility. Here we compared some of the most popular ...single-echo de-noising methods-regression of Motion parameters, White matter and Cerebrospinal fluid signals (MWC method), FMRIB's ICA-based X-noiseifier (FIX) and ICA-based Automatic Removal Of Motion Artifacts (ICA-AROMA)-with a multi-echo approach (ME-ICA) that exploits the linear dependency of BOLD on the echo time. Data were acquired using a clinical scanner and included 30 young, healthy participants (minimal head motion) and 30 Attention Deficit Hyperactivity Disorder patients (greater head motion). De-noising effectiveness was assessed in terms of data quality after each cleanup procedure, ability to uncouple BOLD signal and motion and preservation of default mode network (DMN) functional connectivity. Most cleaning methods showed a positive impact on data quality. However, based on the investigated metrics, ME-ICA was the most robust. It minimized the impact of motion on FC even for high motion participants and preserved DMN functional connectivity structure. The high-quality results obtained using ME-ICA suggest that using a multi-echo EPI sequence, reliable rfMRI data can be obtained in a clinical setting.