Methods for cleaning the BOLD fMRI signal

Available online 9 December 2016 http://www.sciencedirect.com/science/article/pii/S1053811916307418?via%3Dihubhttp://www.sciencedirect.com/science/article/pii/S1053811916307418?via%3DihubBlood 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.This work was supported by the Spanish Ministry of Economy and Competitiveness [Grant PSI 2013–42343 Neuroimagen Multimodal], the Severo Ochoa Programme for Centres/Units of Excellence in R & D [SEV-2015-490], and the research and writing of the paper were supported by the NIMH and NINDS Intramural Research Programs (ZICMH002888) of the NIH/HHS

Towards a better understanding of the impact of heart rate on the BOLD signal: a new method for physiological noise correction and its applications

Functional magnetic resonance imaging (fMRI) based on blood oxygenation level-dependent (BOLD) contrast allows non-invasive examination of brain activity and is widely used in the neuroimaging field. The BOLD contrast mechanism reflects hemodynamic changes resulting from a complex interplay of blood flow, blood volume, and oxygen consumption. Heart rate (HR) variations are the most intriguing and less understood physiological processes affecting the BOLD signal, as they are the result of a wide variety of interacting factors. The use of the response function that best models HR-induced signal changes, called cardiac response function (CRF), is an effective method to reduce HR noise in fMRI. However, current models of physiological noise correction based on CRF, i.e. canonical and individual, either do not take into account variations in HR between subjects, and are thus inadequate for cohorts with varying HR, or require time-consuming quality control of individual physiological recordings and derived CRFs. By analyzing a large cohort of healthy individuals, the results presented in this thesis show that different HRs influence the BOLD signal and their corresponding spectra differently. A further finding is that HR plays an essential role in determining the shape of the CRF. Slower HRs produce a smoothed CRF with a single well-defined maximum, while faster HRs cause a second maximum. Taking advantage of this dependence of the CRF on HR, a novel method is proposed to model HR-induced fluctuations in the BOLD signal more accurately than current approaches of physiological noise correction. This method, called HR-based CRF, consists of two CRFs: one for HRs below 68 bpm and one for HRs above this value. HR-based CRFs can be directly applied to the fMRI data without the time-consuming task of deriving a CRF for each subject while accounting for inter-subject variability in HR response

Identifying Respiration-Related Aliasing Artifacts in the Rodent Resting-State fMRI

Resting-state functional magnetic resonance imaging (rs-fMRI) combined with optogenetics and electrophysiological/calcium recordings in animal models is becoming a popular platform to investigate brain dynamics under specific neurological states. Physiological noise originating from the cardiac and respiration signal is the dominant interference in human rs-fMRI and extensive efforts have been made to reduce these artifacts from the human data. In animal fMRI studies, physiological noise sources including the respiratory and cardiorespiratory artifacts to the rs-fMRI signal fluctuation have typically been less investigated. In this article, we demonstrate evidence of aliasing effects into the low-frequency rs-fMRI signal fluctuation mainly due to respiration-induced B0 offsets in anesthetized rats. This aliased signal was examined by systematically altering the fMRI sampling rate, i.e., the time of repetition (TR), in free-breathing conditions and by adjusting the rate of ventilation. Anesthetized rats under ventilation showed a significantly narrower frequency bandwidth of the aliasing effect than free-breathing animals. It was found that the aliasing effect could be further reduced in ventilated animals with a muscle relaxant. This work elucidates the respiration-related aliasing effects on the rs-fMRI signal fluctuation from anesthetized rats, indicating non-negligible physiological noise needed to be taken care of in both awake and anesthetized animal rs-fMRI studies

Sampling Rate Effects on Resting State fMRI Metrics

Low image sampling rates used in resting state functional magnetic resonance imaging (rs-fMRI) may cause aliasing of the cardiorespiratory pulsations over the very low frequency (VLF) BOLD signal fluctuations which reflects to functional connectivity (FC). In this study, we examine the effect of sampling rate on currently used rs-fMRI FC metrics. Ultra-fast fMRI magnetic resonance encephalography (MREG) data, sampled with TR 0.1 s, was downsampled to different subsampled repetition times (sTR, range 0.3–3 s) for comparisons. Echo planar k-space sampling (TR 2.15 s) and interleaved slice collection schemes were also compared against the 3D single shot trajectory at 2.2 s sTR. The quantified connectivity metrics included stationary spatial, time, and frequency domains, as well as dynamic analyses. Time domain methods included analyses of seed-based functional connectivity, regional homogeneity (ReHo), coefficient of variation, and spatial domain group level probabilistic independent component analysis (ICA). In frequency domain analyses, we examined fractional and amplitude of low frequency fluctuations. Aliasing effects were spatially and spectrally analyzed by comparing VLF (0.01–0.1 Hz), respiratory (0.12–0.35 Hz) and cardiac power (0.9–1.3 Hz) FFT maps at different sTRs. Quasi-periodic pattern (QPP) of VLF events were analyzed for effects on dynamic FC methods. The results in conventional time and spatial domain analyses remained virtually unchanged by the different sampling rates. In frequency domain, the aliasing occurred mainly in higher sTR (1–2 s) where cardiac power aliases over respiratory power. The VLF power maps suffered minimally from increasing sTRs. Interleaved data reconstruction induced lower ReHo compared to 3D sampling (p < 0.001). Gradient recalled echo-planar imaging (EPI BOLD) data produced both better and worse metrics. In QPP analyses, the repeatability of the VLF pulse detection becomes linearly reduced with increasing sTR. In conclusion, the conventional resting state metrics (e.g., FC, ICA) were not markedly affected by different TRs (0.1–3 s). However, cardiorespiratory signals showed strongest aliasing in central brain regions in sTR 1–2 s. Pulsatile QPP and other dynamic analyses benefit linearly from short TR scanning

Does higher sampling rate (multiband + SENSE) improve group statistics - An example from social neuroscience block design at 3T

Multiband (MB) or Simultaneous multi-slice (SMS) acquisition schemes allow the acquisition of MRI signals from more than one spatial coordinate at a time. Commercial availability has brought this technique within the reach of many neuroscientists and psychologists. Most early evaluation of the performance of MB acquisition employed resting state fMRI or the most basic tasks. In this study, we tested whether the advantages of using MB acquisition schemes generalize to group analyses using a cognitive task more representative of typical cognitive neuroscience applications. Twenty-three subjects were scanned on a Philips 3 ​T scanner using five sequences, up to eight-fold acceleration with MB-factors 1 to 4, SENSE factors up to 2 and corresponding TRs of 2.45s down to 0.63s, while they viewed (i) movie blocks showing complex actions with hand object interactions and (ii) control movie blocks without hand object interaction. Data were processed using a widely used analysis pipeline implemented in SPM12 including the unified segmentation and canonical HRF modelling. Using random effects group-level, voxel-wise analysis we found that all sequences were able to detect the basic action observation network known to be recruited by our task. The highest t-values were found for sequences with MB4 acceleration. For the MB1 sequence, a 50% bigger voxel volume was needed to reach comparable t-statistics. The group-level t-values for resting state networks (RSNs) were also highest for MB4 sequences. Here the MB1 sequence with larger voxel size did not perform comparable to the MB4 sequence. Altogether, we can thus recommend the use of MB4 (and SENSE 1.5 or 2) on a Philips scanner when aiming to perform group-level analyses using cognitive block design fMRI tasks and voxel sizes in the range of cortical thickness (e.g. 2.7 ​mm isotropic). While results will not be dramatically changed by the use of multiband, our results suggest that MB will bring a moderate but significant benefit

Depth-Dependent Physiological Modulators of the BOLD Response in the Human Motor Cortex

This dissertation proposes a set of methods for improving spatial localization of cerebral metabolic changes using functional magnetic resonance imaging (fMRI). Blood oxygen level dependent (BOLD) fMRI estabilished itself as the most frequently used technique for mapping brain activity in humans. It is non-invasive and allows to obtain information about brain oxygenation changes in a few minutes. It was discovered in 1990 and, since then, it contributed enormously to the developments in neuroscientific research. Nevertheless, the BOLD contrast suffers from inherent limitations. This comes from the fact that the observed response is the result of a complex interplay between cerebral blood flow (CBF), cerebral blood volume (CBV) and cerebral metabolic rate of oxygen consumption (CMRO2) and has a strong dependency on baseline blood volume and oxygenation. Therefore, the observed response is mislocalized from the site where the metabolic activity takes place and it is subject to high variability across experiments due to normal brain physiology. Since the peak of BOLD changes can be as much as 4 mm apart from the site of metabolic changes, the problem of spatial mislocalization is particularly constraining at submillimeter resolution. Three methods are proposed in this work in order to overcome this limitation and make data more comparable. The first method involves a modification of an estabilished model for calibration of BOLD responses (the dilution model), in order to render it applicable at higher resolutions. The second method proposes a model-free scaling of the BOLD response, based on spatial normalization by a purely vascular response pattern. The third method takes into account the hypothesis that the cortical vasculature could act as a low-pass filter for BOLD fluctuations as the blood is carried downstream, and investigates differences in frequency composition of cortical laminae. All methods are described and tested on a depth-dependent scale in the human motor cortex

Quality data assessment and improvement in pre-processing pipeline to minimize impact of spurious signals in functional magnetic imaging (fMRI)

In the recent years, the field of quality data assessment and signal denoising in functional magnetic resonance imaging (fMRI) is rapidly evolving and the identification and reduction of spurious signal with pre-processing pipeline is one of the most discussed topic. In particular, subject motion or physiological signals, such as respiratory or/and cardiac pulsatility, were showed to introduce false-positive activations in subsequent statistical analyses. Different measures for the evaluation of the impact of motion related artefacts, such as frame-wise displacement and root mean square of movement parameters, and the reduction of these artefacts with different approaches, such as linear regression of nuisance signals and scrubbing or censoring procedure, were introduced. However, we identify two main drawbacks: i) the different measures used for the evaluation of motion artefacts were based on user-dependent thresholds, and ii) each study described and applied their own pre-processing pipeline. Few studies analysed the effect of these different pipelines on subsequent analyses methods in task-based fMRI.The first aim of the study is to obtain a tool for motion fMRI data assessment, based on auto-calibrated procedures, to detect outlier subjects and outliers volumes, targeted on each investigated sample to ensure homogeneity of data for motion. The second aim is to compare the impact of different pre-processing pipelines on task-based fMRI using GLM based on recent advances in resting state fMRI preprocessing pipelines. Different output measures based on signal variability and task strength were used for the assessment