Functional magnetic resonance imaging based on SEEP contrast: response function and anatomical specificity

Functional magnetic resonance imaging based on a non-BOLD-contrast mechanism, which we have termed "SEEP" (Signal Enhancement by Extravascular water Protons), has previously been demonstrated. Here the reproducibility of areas of activity identified with both SEEP and BOLD contrast is assessed in duplicate experiments in healthy volunteers, with relatively high resolution (1.6 mm) image data at 1.5 T. These areas of activity are equally well localized to voxels containing primarily gray matter with the two contrast mechanisms. As in previous studies, areas of SEEP activity are observed to be immediately adjacent to areas of BOLD activity, with very little overlap. The response functions were estimated for both SEEP and BOLD contrast, and are observed to be distinct. The peak SEEP response is observed to lag the BOLD response by approximately 1 s and to decay more slowly with no poststimulus undershoot. Average BOLD signal changes (GE-EPI, TE=50 ms) were observed to be 3.4+/-1.2% (mean+/-S.D.), whereas SEEP signal changes (SE-EPI, TE=23 ms) were 1.9+/-0.5%, consistent with previous studies carried out at 0.2 and 3 T. These observations provide further support for the existence of a non-BOLD-contrast mechanism for fMRI, based on changes in extravascular spin density.     

Measurement and characterization of the human spinal cord SEEP response using event-related spinal fMRI

Although event-related fMRI is able to reliably detect brief changes in brain activity and is now widely used throughout systems and cognitive neuroscience, there have been no previous reports of event-related spinal cord fMRI. This is likely attributable to the various technical challenges associated with spinal fMRI (e.g., imaging a suitable length of the cord, reducing image artifacts from the vertebrae and intervertebral discs, and dealing with physiological noise from spinal cord motion). However, with many of these issues now resolved, the largest remaining impediment for event-related spinal fMRI is a deprived understanding of the spinal cord fMRI signal time course. Therefore, in this study, we used a proton density-weighted HASTE sequence, with functional contrast based on signal enhancement by extravascular water protons (SEEP), and a motion-compensating GLM analysis to (i) characterize the SEEP response function in the human cervical spinal cord and (ii) demonstrate the feasibility of event-related spinal fMRI. This was achieved by applying very brief (1 s) epochs of 22°C thermal stimulation to the palm of the hand and measuring the impulse response function. Our results suggest that the spinal cord SEEP response (time to peak ≈8 s; FWHM ≈4 s; and probably lacking pre- and poststimulus undershoots) is slower than previous estimates of SEEP or BOLD responses in the brain, but faster than previously reported spinal cord BOLD responses. Finally, by detecting and mapping consistent signal-intensity changes within and across subjects, and validating these regions with a block-designed experiment, this study represents the first successful demonstration of event-related spinal fMRI.     

Comparison of α-chloralose,medetomidine and isoflurane anesthesia for functional connectivity mapping in the rat

Functional connectivity measures based upon low-frequency blood-oxygenation-level-dependent functional magnetic resonance imaging (BOLD fMRI) signal fluctuations have become a widely used tool for investigating spontaneous brain activity in humans. Still unknown, however, is the precise relationship between neural activity, the hemodynamic response and fluctuations in the MRI signal. Recent work from several groups had shown that correlated low-frequency fluctuations in the BOLD signal can be detected in the anesthetized rat — a first step toward elucidating this relationship. Building on this preliminary work, through this study, we demonstrate that functional connectivity observed in the rat depends strongly on the type of anesthesia used. Power spectra of spontaneous fluctuations and the cross-correlation-based connectivity maps from rats anesthetized with α-chloralose, medetomidine or isoflurane are presented using a high-temporal-resolution imaging sequence that ensures minimal contamination from physiological noise. The results show less localized correlation in rats anesthetized with isoflurane as compared with rats anesthetized with α-chloralose or medetomidine. These experiments highlight the utility of using different types of anesthesia to explore the fundamental physiological relationships of the BOLD signal and suggest that the mechanisms contributing to functional connectivity involve a complicated relationship between changes in neural activity, neurovascular coupling and vascular reactivity.     

Illuminating the BOLD signal: combined fMRI-fNIRS studies

Functional magnetic resonance imaging (fMRI) is currently combined with electrophysiological methods to identify the relationship between neuronal activity and the blood oxygenation level-dependent (BOLD) signal. Several processes like neuronal activity, synaptic activity, vascular dilation, blood volume and oxygenation changes underlie both response modalities, that is, the electrophysiological signal and the vascular response. However, accessing single process relationships is absolutely mandatory when aiming at a deeper understanding of neurovascular coupling and necessitates studies on the individual building blocks of the vascular response. Combined fMRI and functional near-infrared spectroscopy studies have been performed to validate the correlation of the BOLD signal to the hemodynamic changes in the brain. Here we review the current status of the integration of both technologies and judge these studies in the light of recent findings on neurovascular coupling.     

Quantification of venous blood signal contribution to BOLD functional activation in the auditory cortex at 3 T

Most modern techniques for functional magnetic resonance imaging (fMRI) rely on blood-oxygen-level-dependent (BOLD) contrast as the basic principle for detecting neuronal activation. However, the measured BOLD effect depends on a transfer function related to neurophysiological changes accompanying electrical neural activation. The spatial accuracy and extension of the region of interest are determined by vascular effect, which introduces incertitude on real neuronal activation maps. Our efforts have been directed towards the development of a new methodology that is capable of combining morphological, vascular and functional information; obtaining new insight regarding foci of activation; and distinguishing the nature of activation on a pixel-by-pixel basis. Six healthy volunteers were studied in a parametric auditory functional experiment at 3 T; activation maps were overlaid on a high-resolution brain venography obtained through a novel technique. The BOLD signal intensities of vascular and nonvascular activated voxels were analyzed and compared: it was shown that nonvascular active voxels have lower values for signal peak (P<10(-7)) and area (P<10(-8)) with respect to vascular voxels. The analysis showed how venous blood influenced the measured BOLD signals, supplying a technique to filter possible venous artifacts that potentially can lead to misinterpretation of fMRI results. This methodology, although validated in the auditory cortex activation, maintains a general applicability to any cortical fMRI study, as the basic concepts on which it relies on are not limited to this cortical region. The results obtained in this study can represent the basis for new methodologies and tools that are capable of adding further characterization to the BOLD signal properties.     

Changes in MRI signal intensity during hypercapnic challenge under conscious and anesthetized conditions

Most functional magnetic resonance imaging (fMRI) studies in animals are conducted under anesthesia to minimize motion artifacts. However, methods and techniques have been developed recently for imaging fully conscious rats. Functional MRI studies on conscious animals report enhanced BOLD signal changes as compared to the anesthetized condition. In this study, rats were exposed to different concentrations of carbon dioxide (CO(2)) while conscious and anesthetized to test whether cerebrovascular reactivity may be contributing to these enhanced BOLD signal changes. Hypercapnia produced significantly greater increases in MRI signal intensity in fully conscious animals (6.7-13.3% changes) as when anesthetized with 1% isoflurane (3.2-4.9% changes). In addition, the response to hypercapnia was more immediate in the conscious condition (< 30s) with signal risetimes twice as fast as in the anesthetized state (60s). Both cortical and subcortical brain regions showed a robust, dose- dependent increase in MRI signal intensity with hypercapnic challenge while the animals were conscious but little or no change when anesthetized. Baseline variations in MRI signal were higher while animals were conscious but this was off set by greater signal intensity changes leading to a greater contrast-to-noise ratio, 13.1 in conscious animals, as compared to 8.0 in the anesthetized condition. In summary, cerebral vasculature appears to be more sensitive to hypercapnic challenge in the conscious condition resulting in enhanced T2* MRI signal intensity and the potential for better BOLD signal changes during functional imaging.     

Functional phantom for fMRI: a feasibility study

Functional magnetic resonance imaging (fMRI) reveals changes in blood oxygen level-dependent (BOLD) signal after considerable processing. This paper describes the implementation and testing of an fMRI phantom where electric current applied to a thin wire within a proton-rich medium substituted BOLD distortion of the magnetic field; the scanner detects these two distortions as practically identical signal changes. The magnitude of the change depended on the current strength. The phantom has a number of possible applications. Signal changes across sessions, days, instruments and individuals could be monitored. Placing the phantom close to a subject during an fMRI experiment could allow differentiating sensitivity changes in the scanner due to instrumentation from changes in the subject's state and performance during the experiment. The spatial extent of brain activations and effects of various changes in the chain of image formation could be analyzed using current-induced "activations". Furthermore, the phantom could expedite fMRI sequence development by reducing the need to scan human subjects, who introduce uncertainty to the signal. Thus, this fMRI phantom could be useful for both cognitive fMRI studies and scanner calibration.     

Regional homogeneity of resting-state fMRI contributes to both neurovascular and task activation variations

The task induced blood oxygenation level dependent signal changes observed using functional magnetic resonance imaging (fMRI) are critically dependent on the relationship between neuronal activity and hemodynamic response. Therefore, understanding the nature of neurovascular coupling is important when interpreting fMRI signal changes evoked via task. In this study, we used regional homogeneity (ReHo), a measure of local synchronization of the BOLD time series, to investigate whether the similarities of one voxel with the surrounding voxels are a property of neurovascular coupling. FMRI scans were obtained from fourteen subjects during bilateral finger tapping (FTAP), digit–symbol substitution (DSST) and periodic breath holding (BH) paradigm. A resting-state scan was also obtained for each of the subjects for 4 min using identical imaging parameters. Inter-voxel correlation analyses were conducted between the resting-state ReHo, resting-state amplitude of low frequency fluctuations (ALFF), BH responses and task activations within the masks related to task activations. There was a reliable mean voxel-wise spatial correlation between ReHo and other neurovascular variables (BH responses and ALFF). We observed a moderate correlation between ReHo and task activations (FTAP: r = 0.32; DSST: r = 0.22) within the task positive network and a small yet reliable correlation within the default mode network (DSST: r = − 0.08). Subsequently, a linear regression was used to estimate the contribution of ReHo, ALFF and BH responses to the task activated voxels. The unique contribution of ReHo was minimal. The results suggest that regional synchrony of the BOLD activity is a property that can explain the variance of neurovascular coupling and task activations; but its contribution to task activations can be accounted for by other neurovascular factors such as the ALFF.     

Integration of EEG source imaging and fMRI during continuous viewing of natural movies

Electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) are noninvasive neuroimaging tools which can be used to measure brain activity with excellent temporal and spatial resolution, respectively. By combining the neural and hemodynamic recordings from these modalities, we can gain better insight into how and where the brain processes complex stimuli, which may be especially useful in patients with different neural diseases. However, due to their vastly different spatial and temporal resolutions, the integration of EEG and fMRI recordings is not always straightforward. One fundamental obstacle has been that paradigms used for EEG experiments usually rely on event-related paradigms, while fMRI is not limited in this regard. Therefore, here we ask whether one can reliably localize stimulus-driven EEG activity using the continuously varying feature intensities occurring in natural movie stimuli presented over relatively long periods of time. Specifically, we asked whether stimulus-driven aspects in the EEG signal would be co-localized with the corresponding stimulus-driven BOLD signal during free viewing of a movie. Secondly, we wanted to integrate the EEG signal directly with the BOLD signal, by estimating the underlying impulse response function (IRF) that relates the BOLD signal to the underlying current density in the primary visual area (V1). We made sequential fMRI and 64-channel EEG recordings in seven subjects who passively watched 2-min-long segments of a James Bond movie. To analyze EEG data in this natural setting, we developed a method based on independent component analysis (ICA) to reject EEG artifacts due to blinks, subject movement, etc., in a way unbiased by human judgment. We then calculated the EEG source strength of this artifact-free data at each time point of the movie within the entire brain volume using low-resolution electromagnetic tomography (LORETA). This provided for every voxel in the brain (i.e., in 3D space) an estimate of the current density at every time point. We then carried out a correlation between the time series of visual contrast changes in the movie with that of EEG voxels. We found the most significant correlations in visual area V1, just as seen in previous fMRI studies (Bartels A, Zeki, S, Logothetis NK. Natural vision reveals regional specialization to local motion and to contrast-invariant, global flow in the human brain. Cereb Cortex 2008;18(3):705–717), but on the time scale of milliseconds rather than of seconds. To obtain an estimate of how the EEG signal relates to the BOLD signal, we calculated the IRF between the BOLD signal and the estimated current density in area V1. We found that this IRF was very similar to that observed using combined intracortical recordings and fMRI experiments in nonhuman primates. Taken together, these findings open a new approach to noninvasive mapping of the brain. It allows, firstly, the localization of feature-selective brain areas during natural viewing conditions with the temporal resolution of EEG. Secondly, it provides a tool to assess EEG/BOLD transfer functions during processing of more natural stimuli. This is especially useful in combined EEG/fMRI experiments, where one can now potentially study neural-hemodynamic relationships across the whole brain volume in a noninvasive manner.     

An investigation of the relationship between BOLD and perfusion signal changes during epileptic generalised spike wave activity

In pathological conditions interpretation of functional magnetic resonance imaging (fMRI) results can be difficult. This is due to a reliance on the assumed coupling between neuronal activity and changes in cerebral blood flow (CBF) and oxygenation. We wanted to investigate the coupling between blood oxygen level dependant contrast (BOLD) and CBF time courses in epilepsy patients with generalised spike wave activity (GSW) to better understand the underlying mechanisms behind the EEG-fMRI signal changes observed, especially in regions of negative BOLD response (NBR). Four patients with frequent GSW were scanned with simultaneous electroencephalographic (EEG)-fMRI with BOLD and arterial spin labeling (ASL) sequences. We examined the relationship between simultaneous CBF and BOLD measurements by looking at the correlation of the two signals in terms of percentage signal change on a voxel-by-voxel basis. This method is not reliant on coincident activation. BOLD and CBF were positively correlated in patients with epilepsy during background EEG activity and GSW. The subject average value of the Delta CBF/Delta BOLD slope lay between +19 and +36 and also showed spatial variation which could indicate areas with altered vascular response. There was not a significant difference between Delta CBF/Delta BOLD during GSW, suggesting that neurovascular coupling to BOLD signal is generally maintained between states and, in particular, within areas of NBR.     

Effects of mild hypoxic hypoxia on poststimulus undershoot of blood-oxygenation-level-dependent fMRI signal in the human visual cortex

Characteristics of the blood-oxygenation-level-dependent (BOLD) functional magnetic resonance imaging (fMRI) signal poststimulus undershoot in the visual cortex were studied at varying levels of arterial blood oxygen saturation (Ysat). Undershoot with an amplitude of -0.6+/-0.2% appeared after positive BOLD response (+1.7+/-0.5%) under control conditions. Cerebral blood volume (CBV), as determined with vascular-space-occupancy-dependent fMRI, increased by 26-43% during the positive BOLD peak, but the CBV proceeded at baseline level during the BOLD poststimulus undershoot. Mild hypoxic hypoxia (Ysat ranging from 0.82 to 0.89) had no effect on the amplitude or duration of poststimulus undershoot in activated BOLD pixels. Hypoxia did not influence CBV during the BOLD poststimulus undershoot. In contrast, the positive BOLD signal at the level of all activated pixels was smaller in hypoxia than in normoxia. The present results show that the BOLD poststimulus undershoot is not influenced by curtailed oxygen availability and that, during the undershoot, CBV is not different from resting state.     

Neonatal auditory activation detected by functional magnetic resonance imaging

The objective of this study was to detect auditory cortical activation in non-sedated neonates employing functional magnetic resonance imaging (fMRI). Using echo-planar functional brain imaging, subjects were presented with a frequency-modulated pure tone; the BOLD signal response was mapped in 5 mm-thick slices running parallel to the superior temporal gyrus. Twenty healthy neonates (13 term, 7 preterm) at term and 4 adult control subjects. Blood oxygen level-dependent (BOLD) signal in response to auditory stimulus was detected in all 4 adults and in 14 of the 20 neonates. FMRI studies of adult subjects demonstrated increased signal in the superior temporal regions during auditory stimulation. In contrast, signal decreases were detected during auditory stimulation in 9 of 14 newborns with BOLD response. fMRI can be used to detect brain activation with auditory stimulation in human infants.     

Brain oscillatory activity during motor imagery in EEG-fMRI coregistration

The purpose of the present work was to investigate the correlation between topographical changes in brain oscillatory activity and the blood oxygenation level-dependent (BOLD) signal during a motor imagery (MI) task using electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) coregistration.     

Origins of intersubject variability of blood oxygenation level dependent and arterial spin labeling fMRI: implications for quantification of brain activity

Accurate localization of brain activity using blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging (fMRI) has been challenged because of the large BOLD signal within distal veins. Arterial spin labeling (ASL) techniques offer greater sensitivity to the microvasculature but possess low temporal resolution and limited brain coverage. In this study, we show that the physiological origins of BOLD and ASL depend on whether percent change or statistical significance is being considered. For BOLD and ASL fMRI data collected during a simple unilateral hand movement task, we found that in the area of the contralateral motor cortex the centre of gravity (CoG) of the intersubject coefficient of variation (CV) of BOLD fMRI was near the brain surface for percent change in signal, whereas the CoG of the intersubject CV for Z-score was in close proximity of sites of brain activity for both BOLD and ASL. These findings suggest that intersubject variability of BOLD percent change is vascular in origin, whereas the origin of inter-subject variability of Z-score is neuronal for both BOLD and ASL. For longer duration tasks (12 s or greater), however, there was a significant correlation between BOLD and ASL percent change, which was not evident for short duration tasks (6 s). These findings suggest that analyses directly comparing percent change in BOLD signal between pre-defined regions of interest using short duration stimuli, as for example in event-related designs, may be heavily weighted by large-vessel responses rather than neuronal responses.     

The longitudinal changes of BOLD response and cerebral hemodynamics from acute to subacute stroke. A fMRI and TCD study

Background  

By mapping the dynamics of brain reorganization, functional magnetic resonance imaging MRI (fMRI) has allowed for significant progress in understanding cerebral plasticity phenomena after a stroke. However, cerebro-vascular diseases can affect blood oxygen level dependent (BOLD) signal. Cerebral autoregulation is a primary function of cerebral hemodynamics, which allows to maintain a relatively constant blood flow despite changes in arterial blood pressure and perfusion pressure. Cerebral autoregulation is reported to become less effective in the early phases post-stroke.     

Synchronized detection of minute electrical currents with MRI using Lorentz effect imaging

The blood oxygenation level-dependent (BOLD) effect is the most commonly used contrast mechanism in functional magnetic resonance imaging (fMRI), due to its relatively high spatial resolution and sensitivity. However, the ability of BOLD fMRI to accurately localize neuronal activation in space and time is limited by the inherent hemodynamic modulation. There is hence a need to develop alternative MRI methods that can directly image neuroelectric activity, thereby achieving both a high temporal resolution and spatial specificity as compared to conventional BOLD fMRI. In this paper, we extend the Lorentz effect imaging technique, which can detect spatially incoherent yet temporally synchronized minute electrical activity in a strong magnetic field, and demonstrate its feasibility for imaging randomly oriented electrical currents on the order of microamperes with a temporal resolution on the order of milliseconds in gel phantoms. This constitutes a promising step towards its application to direct imaging of neuroelectric activity in vivo, which has the same order of current density and temporal synchrony.     

Impact of intravenous nicotine on BOLD signal response to photic stimulation

Functional magnetic resonance imaging (fMRI) is increasingly being applied in the study of brain effects of nicotine. In addition, because tobacco smoking is common, many subjects studied with fMRI for other reasons may have appreciable levels of nicotine in plasma and brain during scanning. However, there is concern that the vascular effects of nicotine may alter the coupling between blood oxygen level dependent (BOLD) signal and neuronal activity. The objective of this study was to test for evidence of alteration of BOLD signal response of occipital cortex, a region with a relatively low concentration of neuronal nicotine receptors, to photic stimulation during intravenous infusion of nicotine. Nine nicotine dependent healthy smokers were withdrawn from nicotine under controlled conditions and then scanned while receiving photic stimulation and successive intravenous infusions of saline and nicotine. No evidence for an effect of nicotine on BOLD signal response to photic stimulation was detected at the doses studied. This observation suggests that nicotine does not alter the coupling between BOLD signal and neuronal activity in the visual cortex.     

Functional MRI using super-resolved spatiotemporal encoding

Recently, new ultrafast imaging sequences such as rapid acquisition by sequential excitation and refocusing (RASER) and hybrid spatiotemporal encoding (SPEN) magnetic resonance imaging (MRI) have been proposed, in which the phase encoding of conventional echo planar imaging (EPI) is replaced with a SPEN. In contrast to EPI, SPEN provides significantly higher immunity to frequency heterogeneities including those caused by B₀ inhomogeneities and chemical shift offsets. Utilizing the inherent robustness of SPEN, it was previously shown that RASER can be used to successfully perform functional MRI (fMRI) experiments in the orbitofrontal cortex — a task which is challenging using EPI due to strong magnetic susceptibility variation near the air-filled sinuses. Despite this superior performance, systematic analyses have shown that, in its initial implementation, the use of SPEN was penalized by lower signal-to-noise ratio (SNR) and higher radiofrequency power deposition as compared to EPI-based methods. A recently developed reconstruction algorithm based on super-resolution principles is able to alleviate both of these shortcomings; the use of this algorithm is hereby explored within an fMRI context. Specifically, a series of fMRI measurements on the human visual cortex confirmed that the super-resolution algorithm retains the statistical significance of the blood oxygenation level dependent (BOLD) response, while significantly reducing the power deposition associated with SPEN and restoring the SNR to levels that are comparable with those of EPI.     

Relationship between neural and hemodynamic signals during spontaneous activity studied with temporal kernel CCA

Functional magnetic resonance imaging (fMRI) based on the so-called blood oxygen level-dependent (BOLD) contrast is a powerful tool for studying brain function not only locally but also on the large scale. Most studies assume a simple relationship between neural and BOLD activity, in spite of the fact that it is important to elucidate how the “when” and “what” components of neural activity are correlated to the “where” of fMRI data. Here we conducted simultaneous recordings of neural and BOLD signal fluctuations in primary visual (V1) cortex of anesthetized monkeys. We explored the neurovascular relationship during periods of spontaneous activity by using temporal kernel canonical correlation analysis (tkCCA). tkCCA is a multivariate method that can take into account any features in the signals that univariate analysis cannot. The method detects filters in voxel space (for fMRI data) and in frequency–time space (for neural data) that maximize the neurovascular correlation without any assumption of a hemodynamic response function (HRF). Our results showed a positive neurovascular coupling with a lag of 4–5 s and a larger contribution from local field potentials (LFPs) in the γ range than from low-frequency LFPs or spiking activity. The method also detected a higher correlation around the recording site in the concurrent spatial map, even though the pattern covered most of the occipital part of V1. These results are consistent with those of previous studies and represent the first multivariate analysis of intracranial electrophysiology and high-resolution fMRI.     

Comparison between end-tidal CO₂ and respiration volume per time for detecting BOLD signal fluctuations during paced hyperventilation

Respiratory motion and capnometry monitoring were performed during blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) of the brain while a series of paced hyperventilation tasks were performed that caused significant hypocapnia. Respiration volume per time (RVT) and end-tidal carbon dioxide (ETCO(2)) were determined and compared for their ability to explain BOLD contrast changes in the data. A 35% decrease in ETCO(2) was observed along with corresponding changes in RVT. A best-fit ETCO(2) response function, with an average initial peak delay time of 12 s, was empirically determined. ETCO(2) data convolved with this response function was more strongly and prevalently correlated to BOLD signal changes than RVT data convolved with the corresponding respiration response function. The results suggest that ETCO(2) better models BOLD signal fluctuations in fMRI experiments with significant transient hypocapnia. This is due to hysteresis in the ETCO(2) response when moving from hypocapnia to normocapnia, compared to moving from normocapnia to hypocapnia.