Multi-echo acquisition

Reduced lateralization of multiple functional brain networks in autistic males

BACKGROUND

Autism spectrum disorder has been linked to a variety of organizational and developmental deviations in the brain. One such organizational difference involves hemispheric lateralization, which may be localized to language-relevant regions of the brain or distributed more broadly.

METHODS

In the present study, we estimated brain hemispheric lateralization in autism based on each participant's unique functional neuroanatomy rather than relying on group-averaged data. Additionally, we explored potential relationships between the lateralization of the language network and behavioral phenotypes including verbal ability, language delay, and autism symptom severity. We hypothesized that differences in hemispheric asymmetries in autism would be limited to the language network, with the alternative hypothesis of pervasive differences in lateralization. We tested this and other hypotheses by employing a cross-sectional dataset of 118 individuals (48 autistic, 70 neurotypical). Using resting-state fMRI, we generated individual network parcellations and estimated network asymmetries using a surface area-based approach. A series of multiple regressions were then used to compare network asymmetries for eight significantly lateralized networks between groups.

RESULTS

We found significant group differences in lateralization for the left-lateralized Language (d = -0.89), right-lateralized Salience/Ventral Attention-A (d = 0.55), and right-lateralized Control-B (d = 0.51) networks, with the direction of these group differences indicating less asymmetry in autistic males. These differences were robust across different datasets from the same participants. Furthermore, we found that language delay stratified language lateralization, with the greatest group differences in language lateralization occurring between autistic males with language delay and neurotypical individuals.

CONCLUSIONS

These findings evidence a complex pattern of functional lateralization differences in autism, extending beyond the Language network to the Salience/Ventral Attention-A and Control-B networks, yet not encompassing all networks, indicating a selective divergence rather than a pervasive one. Moreover, we observed an association between Language network lateralization and language delay in autistic males.

Reduced Lateralization of Multiple Functional Brain Networks in Autistic Males

Background

Autism spectrum disorder has been linked to a variety of organizational and developmental deviations in the brain. One such organizational difference involves hemispheric lateralization, which may be localized to language-relevant regions of the brain or distributed more broadly.

Methods

In the present study, we estimated brain hemispheric lateralization in autism based on each participant's unique functional neuroanatomy rather than relying on group-averaged data. Additionally, we explored potential relationships between the lateralization of the language network and behavioral phenotypes including verbal ability, language delay, and autism symptom severity. We hypothesized that differences in hemispheric asymmetries in autism would be limited to the language network, with the alternative hypothesis of pervasive differences in lateralization. We tested this and other hypotheses by employing a cross-sectional dataset of 118 individuals (48 autistic, 70 neurotypical). Using resting-state fMRI, we generated individual network parcellations and estimated network asymmetries using a surface area-based approach. A series of multiple regressions were then used to compare network asymmetries for eight significantly lateralized networks between groups.

Results

We found significant group differences in lateralization for the left-lateralized Language (d = -0.89), right-lateralized Salience/Ventral Attention-A (d = 0.55), and right-lateralized Control-B (d = 0.51) networks, with the direction of these group differences indicating less asymmetry in autistic individuals. These differences were robust across different datasets from the same participants. Furthermore, we found that language delay stratified language lateralization, with the greatest group differences in language lateralization occurring between autistic individuals with language delay and neurotypical individuals.

Limitations

The generalizability of our findings is restricted due to the male-only sample and greater representation of individuals with high verbal and cognitive performance.

Conclusions

These findings evidence a complex pattern of functional lateralization differences in autism, extending beyond the Language network to the Salience/Ventral Attention-A and Control-B networks, yet not encompassing all networks, indicating a selective divergence rather than a pervasive one. Furthermore, a differential relationship was identified between Language network lateralization and specific symptom profiles (namely, language delay) of autism.

Framewise multi-echo distortion correction for superior functional MRI

Functional MRI (fMRI) data are severely distorted by magnetic field (B0) inhomogeneities which currently must be corrected using separately acquired field map data. However, changes in the head position of a scanning participant across fMRI frames can cause changes in the B0 field, preventing accurate correction of geometric distortions. Additionally, field maps can be corrupted by movement during their acquisition, preventing distortion correction altogether. In this study, we use phase information from multi-echo (ME) fMRI data to dynamically sample distortion due to fluctuating B0 field inhomogeneity across frames by acquiring multiple echoes during a single EPI readout. Our distortion correction approach, MEDIC (Multi-Echo DIstortion Correction), accurately estimates B0 related distortions for each frame of multi-echo fMRI data. Here, we demonstrate that MEDIC's framewise distortion correction produces improved alignment to anatomy and decreases the impact of head motion on resting-state functional connectivity (RSFC) maps, in higher motion data, when compared to the prior gold standard approach (i.e., TOPUP). Enhanced framewise distortion correction with MEDIC, without the requirement for field map collection, furthers the advantage of multi-echo over single-echo fMRI.

Scene Perception and Visuospatial Memory Converge at the Anterior Edge of Visually Responsive Cortex

To fluidly engage with the world, our brains must simultaneously represent both the scene in front of us and our memory of the immediate surrounding environment (i.e., local visuospatial context). How does the brain's functional architecture enable sensory and mnemonic representations to closely interface while also avoiding sensory-mnemonic interference? Here, we asked this question using first-person, head-mounted virtual reality and fMRI. Using virtual reality, human participants of both sexes learned a set of immersive, real-world visuospatial environments in which we systematically manipulated the extent of visuospatial context associated with a scene image in memory across three learning conditions, spanning from a single FOV to a city street. We used individualized, within-subject fMRI to determine which brain areas support memory of the visuospatial context associated with a scene during recall (Experiment 1) and recognition (Experiment 2). Across the whole brain, activity in three patches of cortex was modulated by the amount of known visuospatial context, each located immediately anterior to one of the three scene perception areas of high-level visual cortex. Individual subject analyses revealed that these anterior patches corresponded to three functionally defined place memory areas, which selectively respond when visually recalling personally familiar places. In addition to showing activity levels that were modulated by the amount of visuospatial context, multivariate analyses showed that these anterior areas represented the identity of the specific environment being recalled. Together, these results suggest a convergence zone for scene perception and memory of the local visuospatial context at the anterior edge of high-level visual cortex.SIGNIFICANCE STATEMENT As we move through the world, the visual scene around us is integrated with our memory of the wider visuospatial context. Here, we sought to understand how the functional architecture of the brain enables coexisting representations of the current visual scene and memory of the surrounding environment. Using a combination of immersive virtual reality and fMRI, we show that memory of visuospatial context outside the current FOV is represented in a distinct set of brain areas immediately anterior and adjacent to the perceptually oriented scene-selective areas of high-level visual cortex. This functional architecture would allow efficient interaction between immediately adjacent mnemonic and perceptual areas while also minimizing interference between mnemonic and perceptual representations.

A single case neuroimaging study of tickertape synesthesia

Reading acquisition is enabled by deep changes in the brain's visual system and language areas, and in the links subtending their collaboration. Disruption of those plastic processes commonly results in developmental dyslexia. However, atypical development of reading mechanisms may occasionally result in ticker-tape synesthesia (TTS), a condition described by Francis Galton in 1883 wherein individuals "see mentally in print every word that is uttered (…) as from a long imaginary strip of paper". While reading is the bottom-up translation of letters into speech, TTS may be viewed as its opposite, the top-down translation of speech into internally visualized letters. In a series of functional MRI experiments, we studied MK, a man with TTS. We showed that a set of left-hemispheric areas were more active in MK than in controls during the perception of normal than reversed speech, including frontoparietal areas involved in speech processing, and the Visual Word Form Area, an occipitotemporal region subtending orthography. Those areas were identical to those involved in reading, supporting the construal of TTS as upended reading. Using dynamic causal modeling, we further showed that, parallel to reading, TTS induced by spoken words and pseudowords relied on top-down flow of information along distinct lexical and phonological routes, involving the middle temporal and supramarginal gyri, respectively. Future studies of TTS should shed new light on the neurodevelopmental mechanisms of reading acquisition, their variability and their disorders.

Evaluating the efficacy of multi-echo ICA denoising on model-based fMRI

fMRI is an indispensable tool for neuroscience investigation, but this technique is limited by multiple sources of physiological and measurement noise. These noise sources are particularly problematic for analysis techniques that require high signal-to-noise ratio for stable model fitting, such as voxel-wise modeling. Multi-echo data acquisition in combination with echo-time dependent ICA denoising (ME-ICA) represents one promising strategy to mitigate physiological and hardware-related noise sources as well as motion-related artifacts. However, most studies employing ME-ICA to date are resting-state fMRI studies, and therefore we have a limited understanding of the impact of ME-ICA on complex task or model-based fMRI paradigms. Here, we addressed this knowledge gap by comparing data quality and model fitting performance of data acquired during a visual population receptive field (pRF) mapping (N = 13 participants) experiment after applying one of three preprocessing procedures: ME-ICA, optimally combined multi-echo data without ICA-denoising, and typical single echo processing. As expected, multi-echo fMRI improved temporal signal-to-noise compared to single echo fMRI, with ME-ICA amplifying the improvement compared to optimal combination alone. However, unexpectedly, this boost in temporal signal-to-noise did not directly translate to improved model fitting performance: compared to single echo acquisition, model fitting was only improved after ICA-denoising. Specifically, compared to single echo acquisition, ME-ICA resulted in improved variance explained by our pRF model throughout the visual system, including anterior regions of the temporal and parietal lobes where SNR is typically low, while optimal combination without ICA did not. ME-ICA also improved reliability of parameter estimates compared to single echo and optimally combined multi-echo data without ICA-denoising. Collectively, these results suggest that ME-ICA is effective for denoising task-based fMRI data for modeling analyzes and maintains the integrity of the original data. Therefore, ME-ICA may be beneficial for complex fMRI experiments, including voxel-wise modeling and naturalistic paradigms.

Advances in resting state fMRI acquisitions for functional connectomics

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.

Improving the Breath-Holding CVR Measurement Using the Multiband Multi-Echo EPI Sequence

Blood oxygen level-dependent (BOLD) functional MRI (fMRI) is commonly used to measure cerebrovascular reactivity (CVR), which can convey insightful information about neurovascular health. Breath-holding (BH) has been shown to be a practical vasodilatory stimulus for measuring CVR in clinical settings. The conventional BOLD fMRI approach has some limitations, however, such as susceptibility-induced signal dropout at air tissue interfaces and low BOLD sensitivity especially in areas of low T 2 * . These drawbacks can potentially be mitigated with multi-echo sequences, which acquire several images at different echo times in one shot. When combined with multiband techniques, high temporal resolution images can be acquired. This study compared an advanced multiband multi-echo (MBME) echo planar imaging (EPI) sequence with an existing multiband single-echo (MB) sequence to evaluate the repeatability and sensitivity of BH activation and CVR mapping. Images were acquired from 28 healthy volunteers, of which 18 returned for repeat imaging. Both MBME and MB data were pre-processed using both standard and advanced denoising techniques. The MBME data was further processed by combining echoes using a T 2 * -weighted approach and denoising using multi-echo independent component analysis. BH activation was calculated using a general linear model and the respiration response function. CVR was computed as the percent change related to the activation. To account for differences in CVR related to TE, relative CVR (rCVR) was computed and normalized to the mean gray matter CVR. Test-retest metrics were assessed with the Dice coefficient, rCVR difference, within subject coefficient of variation, and the intraclass correlation coefficient. Our findings demonstrate that rCVR for MBME scans were significantly higher than for MB scans across most of the gray matter. In areas of high susceptibility-induced signal dropout, however, MBME rCVR was significantly less than MB rCVR due to artifactually high rCVR for MB scans in these regions. MBME rCVR showed improved test-retest metrics compared with MB. Overall, the MBME sequence displayed superior BOLD sensitivity, improved specificity in areas of signal dropout on MBME scans, enhanced reliability, and reduced variability across subjects compared with MB acquisitions. Our results suggest that the MBME EPI sequence is a promising tool for imaging CVR.

Rapid Precision Functional Mapping of Individuals Using Multi-Echo fMRI

Resting-state functional magnetic resonance imaging (fMRI) is widely used in cognitive and clinical neuroscience, but long-duration scans are currently needed to reliably characterize individual differences in functional connectivity (FC) and brain network topology. In this report, we demonstrate that multi-echo fMRI can improve the reliability of FC-based measurements. In four densely sampled individual humans, just 10 min of multi-echo data yielded better test-retest reliability than 30 min of single-echo data in independent datasets. This effect is pronounced in clinically important brain regions, including the subgenual cingulate, basal ganglia, and cerebellum, and is linked to three biophysical signal mechanisms (thermal noise, regional variability in the rate of T₂* decay, and S₀-dependent artifacts) with spatially distinct influences. Together, these findings establish the potential utility of multi-echo fMRI for rapid precision mapping using experimentally and clinically tractable scan times and will facilitate longitudinal neuroimaging of clinical populations.

Lynch et al. demonstrate that the test-retest reliability of resting-state connectivity measurements can be improved using multi-echo fMRI. This effect is pronounced in clinically important brain regions and could help facilitate precision mapping of functional brain networks in healthy people and patient populations.

Distraction decreases rIFG-putamen connectivity during goal-directed effort for food rewards

Distracted eating can lead to increased food intake, but it is unclear how. We aimed to assess how distraction affects motivated, goal-directed responses for food reward after satiation. Thirty-eight healthy normal-weight participants (28F; 10M) performed a visual detection task varying in attentional load (high vs. low distraction) during fMRI. Simultaneously, they exerted effort for sweet and savory food rewards by repeated button presses. Two fMRI runs were separated by sensory-specific satiation (outcome devaluation) of one of the (sweet or savory) reward outcomes, to assess outcome-sensitive, goal-directed, responses (valued vs. devalued reward, post vs. pre satiation). We could not verify our primary hypothesis that more distraction leads to less activation in ventromedial prefrontal cortex (vmPFC) during goal-directed effort. Behaviorally, distraction also did not affect effort for food reward following satiation across subjects. For our secondary hypothesis, we assessed whether distraction affected other fronto-striatal regions during goal-directed effort. We did not obtain such effects at our whole-brain corrected threshold, but at an exploratory uncorrected threshold (p < 0.001), distraction decreased goal-directed responses (devalued vs. valued) in the right inferior frontal gyrus (rIFG). We continued with this rIFG region for the next secondary hypothesis; specifically, that distraction would reduce functional connectivity with the fronto-striatal regions found in the previous analyses. Indeed, distraction decreased functional connectivity between the rIFG and left putamen for valued versus devalued food rewards (pFWE(cluster) < 0.05). In an exploratory brain-behavior analysis, we showed that distraction-sensitive rIFG-responses correlated negatively (r = - 0.40; p = 0.014) with the effect of distraction on effort. Specifically, decreased distraction-related rIFG-responses were associated with increased effort for food reward after satiation. We discuss the absence of distraction effects on goal-directed responses in vmPFC and in behavior across participants. Moreover, based on our significant functional connectivity and brain-behavior results, we suggest that distraction might attenuate the ability to inhibit responses for food reward after satiation by affecting the rIFG and its connection to the putamen.

Causes and Consequences of Diagnostic Heterogeneity in Depression: Paths to Discovering Novel Biological Depression Subtypes

Depression is a highly heterogeneous syndrome that bears only modest correlations with its biological substrates, motivating a renewed interest in rethinking our approach to diagnosing depression for research purposes and new efforts to discover subtypes of depression anchored in biology. Here, we review the major causes of diagnostic heterogeneity in depression, with consideration of both clinical symptoms and behaviors (symptomatology and trajectory of depressive episodes) and biology (genetics and sexually dimorphic factors). Next, we discuss the promise of using data-driven strategies to discover novel subtypes of depression based on functional neuroimaging measures, including dimensional, categorical, and hybrid approaches to parsing diagnostic heterogeneity and understanding its biological basis. The merits of using resting-state functional magnetic resonance imaging functional connectivity techniques for subtyping are considered along with a set of technical challenges and potential solutions. We conclude by identifying promising future directions for defining neurobiologically informed depression subtypes and leveraging them in the future for predicting treatment outcomes and informing clinical decision making.

More than BOLD: Dual-spin populations create functional contrast

PURPOSE

Functional MRI contrast has generally been associated with changes in transverse relaxivity caused by blood oxygen concentration, the so-called blood oxygen level dependent contrast. However, this interpretation of fMRI contrast has been called into question by several recent experiments at high spatial resolution. Experiments were conducted to examine contrast dependencies that cannot be explained only by differences in relaxivity in a single-spin population.

METHODS

Measurements of functional signal and contrast were obtained in human early visual cortex during a high-contrast visual stimulation over a large range of TEs and for several flip angles. Small voxels (1.5 mm) were used to restrict the measurements to cortical gray matter in early visual areas identified using retinotopic mapping procedures.

RESULTS

Measurements were consistent with models that include 2 spin populations. The dominant population has a relatively short transverse lifetime that is strongly modulated by activation. However, functional contrast is also affected by volume changes between this short-lived population and the longer-lived population.

CONCLUSION

Some of the previously observed "nonclassical" behaviors of functional contrast can be explained by these interacting dual-spin populations.

A deconvolution algorithm for multi-echo functional MRI: Multi-echo Sparse Paradigm Free Mapping

This work introduces a novel algorithm for deconvolution of the BOLD signal in multi-echo fMRI data: Multi-echo Sparse Paradigm Free Mapping (ME-SPFM). Assuming a linear dependence of the BOLD percent signal change on the echo time (TE) and using sparsity-promoting regularized least squares estimation, ME-SPFM yields voxelwise time-varying estimates of the changes in the apparent transverse relaxation (ΔR₂^⁎) without prior knowledge of the timings of individual BOLD events. Our results in multi-echo fMRI data collected during a multi-task event-related paradigm at 3 Tesla demonstrate that the maps of R₂^⁎ changes obtained with ME-SPFM at the times of the stimulus trials show high spatial and temporal concordance with the activation maps and BOLD signals obtained with standard model-based analysis. This method yields estimates of ΔR₂^⁎ having physiologically plausible values. Owing to its ability to blindly detect events, ME-SPFM also enables us to map ΔR₂^⁎ associated with spontaneous, transient BOLD responses occurring between trials. This framework is a step towards deciphering the dynamic nature of brain activity in naturalistic paradigms, resting-state or experimental paradigms with unknown timing of the BOLD events.

Cerebral Autoregulation Evidenced by Synchronized Low Frequency Oscillations in Blood Pressure and Resting-State fMRI

Resting-state functional magnetic resonance imaging (rs-fMRI) is a widely used technique for mapping the brain’s functional architecture, so delineating the main sources of variance comprising the signal is crucial. Low frequency oscillations (LFO) that are not of neural origin, but which are driven by mechanisms related to cerebral autoregulation (CA), are present in the blood-oxygenation-level-dependent (BOLD) signal within the rs-fMRI frequency band. In this study we use a MR compatible device (Caretaker, Biopac) to obtain a non-invasive estimate of beat-to-beat mean arterial pressure (MAP) fluctuations concurrently with rs-fMRI at 3T. Healthy adult subjects (n = 9; 5 male) completed two 20-min rs-fMRI scans. MAP fluctuations were decomposed into different frequency scales using a discrete wavelet transform, and oscillations at approximately 0.1 Hz show a high degree of spatially structured correlations with matched frequency fMRI fluctuations. On average across subjects, MAP fluctuations at this scale of the wavelet decomposition explain ∼2.2% of matched frequency fMRI signal variance. Additionally, a simultaneous multi-slice multi-echo acquisition was used to collect 10-min rs-fMRI at three echo times at 7T in a separate group of healthy adults (n = 5; 5 male). Multiple echo times were used to estimate the R₂^∗ decay at every time point, and MAP was shown to strongly correlate with this signal, which suggests a purely BOLD (i.e., blood flow related) origin. This study demonstrates that there is a significant component of the BOLD signal that has a systemic physiological origin, and highlights the fact that not all localized BOLD signal changes necessarily reflect blood flow supporting local neural activity. Instead, these data show that a proportion of BOLD signal fluctuations in rs-fMRI are due to localized control of blood flow that is independent of local neural activity, most likely reflecting more general systemic autoregulatory processes. Thus, fMRI is a promising tool for studying flow changes associated with cerebral autoregulation with high spatial resolution.

Improving the Assessment of Breath-Holding Induced Cerebral Vascular Reactivity Using a Multiband Multi-echo ASL/BOLD Sequence

Breath holding (BH) is a viable vasodilatory stimulus for calculating functional MRI-derived cerebral vascular reactivity (CVR). The BH technique suffers from reduced repeatability compared with gas inhalation techniques; however, extra equipment is needed to perform gas inhalation techniques, and this equipment is not available at all institutions. This study aimed to determine the sensitivity and repeatability of BH activation and CVR using a multiband multi-echo simultaneous arterial spin labelling/blood oxygenation level dependent (ASL/BOLD) sequence. Whole-brain images were acquired in 14 volunteers. Ten subjects returned for repeat imaging. Each subject performed four cycles of 16 s BH on expiration interleaved with paced breathing. Following standard preprocessing, the echoes were combined using a T2*-weighted approach. BOLD and ASL BH activation was computed, and CVR was then determined as the percent signal change related to the activation. The "M" parameter from the Davis Model was also computed by incorporating the ASL signal. Our results showed higher BH activation strength, volume, and repeatability for the combined multi-echo (MEC) data compared with the single-echo data. MEC CVR also had higher repeatability, sensitivity, specificity, and reliability compared with the single-echo BOLD data. These data support the usefulness of an MBME ASL/BOLD acquisition for BH CVR and M measurements.

Multiband multi-echo simultaneous ASL/BOLD for task-induced functional MRI

Typical simultaneous blood oxygenation-level dependent (BOLD) and arterial spin labeling (ASL) sequences acquire two echoes, one perfusion-sensitive and one BOLD-sensitive. However, for ASL, spatial resolution and brain coverage are limited due to the T1 decay of the labeled blood. This study applies a sequence combining a multiband acquisition with four echoes for simultaneous BOLD and pseudo-continuous ASL (pCASL) echo planar imaging (MBME ASL/BOLD) for block-design task-fMRI. A multiband acceleration of four was employed to increase brain coverage and reduce slice-timing effects on the ASL signal. Multi-echo independent component analysis (MEICA) was implemented to automatically denoise the BOLD signal by regressing non-BOLD components. This technique led to increased temporal signal-to-noise ratio (tSNR) and BOLD sensitivity. The MEICA technique was also modified to denoise the ASL signal by regressing artifact and BOLD signals from the first echo time-series. The MBME ASL/BOLD sequence was applied to a finger-tapping task functional MRI (fMRI) experiment. Signal characteristics and activation were evaluated using single echo BOLD, combined ME BOLD, combined ME BOLD after MEICA denoising, perfusion-weighted (PW), and perfusion-weighted after MEICA denoising time-series. The PW data was extracted using both surround subtraction and high-pass filtering followed by demodulation. In addition, the CBF/BOLD response ratio and CBF/BOLD coupling were analyzed. Results showed that the MEICA denoising procedure significantly improved the BOLD signal, leading to increased BOLD sensitivity, tSNR, and activation statistics compared to conventional single echo BOLD data. At the same time, the denoised PW data showed increased tSNR and activation statistics compared to the non-denoised PW data. CBF/BOLD coupling was also increased using the denoised ASL and BOLD data. Our preliminary data suggest that the MBME ASL/BOLD sequence can be employed to collect whole-brain task-fMRI with improved data quality for both BOLD and PW time series, thus improving the results of block-design task fMRI.

Methods for cleaning the BOLD fMRI signal

Blood oxygen-level-dependent functional magnetic resonance imaging (BOLD fMRI) has rapidly become a popular technique for the investigation of brain function in healthy individuals, patients as well as in animal studies. However, the BOLD signal arises from a complex mixture of neuronal, metabolic and vascular processes, being therefore an indirect measure of neuronal activity, which is further severely corrupted by multiple non-neuronal fluctuations of instrumental, physiological or subject-specific origin. This review aims to provide a comprehensive summary of existing methods for cleaning the BOLD fMRI signal. The description is given from a methodological point of view, focusing on the operation of the different techniques in addition to pointing out the advantages and limitations in their application. Since motion-related and physiological noise fluctuations are two of the main noise components of the signal, techniques targeting their removal are primarily addressed, including both data-driven approaches and using external recordings. Data-driven approaches, which are less specific in the assumed model and can simultaneously reduce multiple noise fluctuations, are mainly based on data decomposition techniques such as principal and independent component analysis. Importantly, the usefulness of strategies that benefit from the information available in the phase component of the signal, or in multiple signal echoes is also highlighted. The use of global signal regression for denoising is also addressed. Finally, practical recommendations regarding the optimization of the preprocessing pipeline for the purpose of denoising and future venues of research are indicated. Through the review, we summarize the importance of signal denoising as an essential step in the analysis pipeline of task-based and resting state fMRI studies.

Multi-echo EPI of human fear conditioning reveals improved BOLD detection in ventromedial prefrontal cortex

Standard T₂^* weighted functional magnetic resonance imaging (fMRI) performed with echo-planar imaging (EPI) suffers from signal loss in the ventromedial prefrontal cortex (vmPFC) due to macroscopic field inhomogeneity. However, this region is of special interest to affective neuroscience and psychiatry. The Multi-echo EPI (MEPI) approach has several advantages over EPI but its performance against EPI in the vmPFC has not yet been examined in a study with sufficient statistical power using a task specifically eliciting activity in this region. We used a fear conditioning task with MEPI to compare the performance of MEPI and EPI in vmPFC and control regions in 32 healthy young subjects. We analyzed activity associated with short (12ms), standard (29ms) and long (46ms) echo times, and a voxel-wise combination of these three echo times. Behavioral data revealed successful differentiation of the conditioned versus safety stimulus; activity in the vmPFC was shown by the contrast "safety stimulus > conditioned stimulus" as in previous research and proved significantly stronger with the combined MEPI than standard single-echo EPI. Then, we aimed to demonstrate that the additional cluster extent (ventral extension) detected in the vmPFC with MEPI reflects activation in a relevant cluster (i.e., not just non-neuronal noise). To do this, we used resting state data from the same subjects to show that the time-course of this region was both connected to bilateral amygdala and the default mode network. Overall, we demonstrate that MEPI (by means of the weighted sum combination approach) outperforms standard EPI in vmPFC; MEPI performs always at least as good as the best echo time for a given brain region but provides all necessary echo times for an optimal BOLD sensitivity for the whole brain. This is relevant for affective neuroscience and psychiatry given the critical role of the vmPFC in emotion regulation.

Multiband multi-echo imaging of simultaneous oxygenation and flow timeseries for resting state connectivity

A novel sequence has been introduced that combines multiband imaging with a multi-echo acquisition for simultaneous high spatial resolution pseudo-continuous arterial spin labeling (ASL) and blood-oxygenation-level dependent (BOLD) echo-planar imaging (MBME ASL/BOLD). Resting-state connectivity in healthy adult subjects was assessed using this sequence. Four echoes were acquired with a multiband acceleration of four, in order to increase spatial resolution, shorten repetition time, and reduce slice-timing effects on the ASL signal. In addition, by acquiring four echoes, advanced multi-echo independent component analysis (ME-ICA) denoising could be employed to increase the signal-to-noise ratio (SNR) and BOLD sensitivity. Seed-based and dual-regression approaches were utilized to analyze functional connectivity. Cerebral blood flow (CBF) and BOLD coupling was also evaluated by correlating the perfusion-weighted timeseries with the BOLD timeseries. These metrics were compared between single echo (E2), multi-echo combined (MEC), multi-echo combined and denoised (MECDN), and perfusion-weighted (PW) timeseries. Temporal SNR increased for the MECDN data compared to the MEC and E2 data. Connectivity also increased, in terms of correlation strength and network size, for the MECDN compared to the MEC and E2 datasets. CBF and BOLD coupling was increased in major resting-state networks, and that correlation was strongest for the MECDN datasets. These results indicate our novel MBME ASL/BOLD sequence, which collects simultaneous high-resolution ASL/BOLD data, could be a powerful tool for detecting functional connectivity and dynamic neurovascular coupling during the resting state. The collection of more than two echoes facilitates the use of ME-ICA denoising to greatly improve the quality of resting state functional connectivity MRI.

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

Multi-echo fMRI, particularly the multi-echo independent component analysis (ME-ICA) algorithm, has previously proven useful for increasing the sensitivity and reducing false positives for functional MRI (fMRI) based resting state connectivity studies. Less is known about its efficacy for task-based fMRI, especially at the single subject level. This work, which focuses exclusively on individual subject results, compares ME-ICA to single-echo fMRI and a voxel-wise T2(⁎) weighted combination of multi-echo data for task-based fMRI under the following scenarios: cardiac-gated block designs, constant repetition time (TR) block designs, and constant TR rapid event-related designs. Performance is evaluated primarily in terms of sensitivity (i.e., activation extent, activation magnitude, percent detected trials and effect size estimates) using five different tasks expected to evoke neuronal activity in a distributed set of regions. The ME-ICA algorithm significantly outperformed all other evaluated processing alternatives in all scenarios. Largest improvements were observed for the cardiac-gated dataset, where ME-ICA was able to reliably detect and remove non-neural T1 signal fluctuations caused by non-constant repetition times. Although ME-ICA also outperformed the other options in terms of percent detection of individual trials for rapid event-related experiments, only 46% of all events were detected after ME-ICA; suggesting additional improvements in sensitivity are required to reliably detect individual short event occurrences. We conclude the manuscript with a detailed evaluation of ME-ICA outcomes and a discussion of how the ME-ICA algorithm could be further improved. Overall, our results suggest that ME-ICA constitutes a versatile, powerful approach for advanced denoising of task-based fMRI, not just resting-state data.

Negative BOLD in default-mode structures measured with EEG-MREG is larger in temporal than extra-temporal epileptic spikes

Introduction: EEG-fMRI detects BOLD changes associated with epileptic interictal discharges (IED) and can identify epileptogenic networks in epilepsy patients. Besides positive BOLD changes, negative BOLD changes have sometimes been observed in the default-mode network, particularly using group analysis. A new fast fMRI sequence called MREG (Magnetic Resonance Encephalography) shows increased sensitivity to detect IED-related BOLD changes compared to the conventional EPI sequence, including frequent occurrence of negative BOLD responses in the DMN. The present study quantifies the concordance between the DMN and negative BOLD related to IEDs of temporal and extra-temporal origin.

Methods: Focal epilepsy patients underwent simultaneous EEG-MREG. Areas of overlap were calculated between DMN regions, defined as precuneus, posterior cingulate, bilateral inferior parietal and mesial prefrontal cortices according to a standardized atlas, and significant negative BOLD changes revealed by an event-related analysis based on the timings of IED seen on EEG. Correlation between IED number/lobe of origin and the overlap were calculated.

Results: 15 patients were analyzed, some showing IED over more than one location resulting in 30 different IED types. The average overlap between negative BOLD and DMN was significantly larger in temporal (23.7 ± 19.6 cm³) than extra-temporal IEDs (7.4 ± 5.1 cm³, p = 0.008). There was no significant correlation between the number of IEDs and the overlap between DMN structures and negative BOLD areas.

Discussion: MREG results in an increased sensitivity to detect negative BOLD responses related to focal IED in single patients, with responses often occurring in DMN regions. In patients with high overlap with the DMN, this suggests that epileptic IEDs may be associated with a brief decrease in attention and cognitive ability. Interestingly this observation was not dependent on the frequency of IED but more common in IED of temporal origin.

High-speed real-time resting-state FMRI using multi-slab echo-volumar imaging

We recently demonstrated that ultra-high-speed real-time fMRI using multi-slab echo-volumar imaging (MEVI) significantly increases sensitivity for mapping task-related activation and resting-state networks (RSNs) compared to echo-planar imaging (Posse et al., 2012). In the present study we characterize the sensitivity of MEVI for mapping RSN connectivity dynamics, comparing independent component analysis (ICA) and a novel seed-based connectivity analysis (SBCA) that combines sliding-window correlation analysis with meta-statistics. This SBCA approach is shown to minimize the effects of confounds, such as movement, and CSF and white matter signal changes, and enables real-time monitoring of RSN dynamics at time scales of tens of seconds. We demonstrate highly sensitive mapping of eloquent cortex in the vicinity of brain tumors and arterio-venous malformations, and detection of abnormal resting-state connectivity in epilepsy. In patients with motor impairment, resting-state fMRI provided focal localization of sensorimotor cortex compared with more diffuse activation in task-based fMRI. The fast acquisition speed of MEVI enabled segregation of cardiac-related signal pulsation using ICA, which revealed distinct regional differences in pulsation amplitude and waveform, elevated signal pulsation in patients with arterio-venous malformations and a trend toward reduced pulsatility in gray matter of patients compared with healthy controls. Mapping cardiac pulsation in cortical gray matter may carry important functional information that distinguishes healthy from diseased tissue vasculature. This novel fMRI methodology is particularly promising for mapping eloquent cortex in patients with neurological disease, having variable degree of cooperation in task-based fMRI. In conclusion, ultra-high-real-time speed fMRI enhances the sensitivity of mapping the dynamics of resting-state connectivity and cerebro-vascular pulsatility for clinical and neuroscience research applications.

The neural processing of masked speech

Spoken language is rarely heard in silence, and a great deal of interest in psychoacoustics has focused on the ways that the perception of speech is affected by properties of masking noise. In this review we first briefly outline the neuroanatomy of speech perception. We then summarise the neurobiological aspects of the perception of masked speech, and investigate this as a function of masker type, masker level and task. This article is part of a Special Issue entitled "Annual Reviews 2013".

Removing motion and physiological artifacts from intrinsic BOLD fluctuations using short echo data

Differing noise variance across study populations has been shown to cause artifactual group differences in functional connectivity measures. In this study, we investigate the use of short echo time functional MRI data to correct for these noise sources in blood oxygenation level dependent (BOLD)-weighted time series. A dual‐echo sequence was used to simultaneously acquire data at both a short (TE = 3.3 ms) and a BOLD-weighted (TE = 35 ms) echo time. This approach is effectively “free,” using dead-time in the pulse sequence to collect an additional echo without affecting overall scan time or temporal resolution. The proposed correction method uses voxelwise regression of the short TE data from the BOLD-weighted data to remove noise variance. In addition to a typical resting state scan, non-compliant behavior associated with patient groups was simulated via increased head motion or physiological fluctuations in 10 subjects. Short TE data showed significant correlations with the traditional motion-related and physiological noise regressors used in current connectivity analyses. Following traditional preprocessing, the extent of significant additional variance explained by the short TE data regressors was significantly correlated with the average head motion across the scan in the resting data (r² = 0.93, p < 0.0001). The reduction in data variance following the inclusion of short TE regressors was also correlated with scan head motion (r² = 0.48, p = 0.027). Task-related data were used to demonstrate the effects of the short TE correction on BOLD activation time series with known temporal structure; the size and strength of the activation were significantly decreased, but it is not clear whether this reflects BOLD contamination in the short TE data or correlated changes in blood volume. Finally, functional connectivity maps of the default mode network were constructed using a seed correlation approach. The effects of short TE correction and low-pass filtering on the resulting correlations maps were compared. Results suggest that short TE correction more accurately differentiates artifactual correlations from the correlations of interest in conditions of amplified noise.