Nonlinear responses of cerebral blood volume, blood flow and blood oxygenation signals during visual stimulation

Hemodynamic-based functional magnetic resonance imaging (fMRI) techniques provide a great utility for noninvasive functional brain mapping. However, because the hemodynamic signals reflect underlying neural activity indirectly, characterization of these signals following brain activation is essential for experimental design and data interpretation. In this report, the linear (or nonlinear) responses to neuronal activation of three hemodynamic parameters based primarily on changes of cerebral blood volume, blood flow and blood oxygenation were investigated by testing these hemodynamic responses' additivity property. Using a recently developed fMRI technique that acquires vascular space occupancy (VASO), arterial spin labeling (ASL) perfusion and blood oxygenation level-dependent (BOLD) signals simultaneously, the additivity property of the three hemodynamic responses in human visual cortex was assessed using various visual stimulus durations. Experiments on healthy volunteers showed that all three hemodynamic-weighted signals responded nonlinearly to stimulus durations less than 4 s, with the degree of nonlinearity becoming more severe as the stimulus duration decreased. Vascular space occupancy and ASL perfusion signals showed similar nonlinearity properties, whereas the BOLD signal was the most nonlinear. These data suggest that caution should be taken in the interpretation of hemodynamic-based signals in fMRI.     

A nonlinear BOLD model accounting for refractory effect by applying the longitudinal relaxation in NMR to the linear BOLD model

A mathematical model to regress the nonlinear blood oxygen level-dependent (BOLD) fMRI signal has been developed by incorporating the refractory effect into the linear BOLD model of the biphasic gamma variate function. The refractory effect was modeled as a relaxation of two separate BOLD capacities corresponding to the biphasic components of the BOLD signal in analogy with longitudinal relaxation of magnetization in NMR. When tested with the published fMRI data of finger tapping, the nonlinear BOLD model with the refractory effect reproduced the nonlinear BOLD effects such as reduced poststimulus undershoot and saddle pattern in a prolonged stimulation as well as the reduced BOLD signal for repetitive stimulation.     

Steady-state activation in somatosensory cortex after changes in stimulus rate during median nerve stimulation

Passive electrical stimulation activates various human somatosensory cortical systems including the contralateral primary somatosensory area (SI), bilateral secondary somatosensory area (SII) and bilateral insula. The effect of stimulation frequency on blood oxygenation level-dependent (BOLD) activity remains unclear. We acquired 3-T functional magnetic resonance imaging (fMRI) in eight healthy volunteers during electrical median nerve stimulation at frequencies of 1, 3 and 10 Hz. During stimulation BOLD signal changes showed activation in the contralateral SI, bilateral SII and bilateral insula. Results of fMRI analysis showed that these areas were progressively active with the increase of rate of stimulation. As a major finding, the contralateral SI showed an increase of peak of BOLD activation from 1 to 3 Hz but reached a plateau during 10-Hz stimulation. Our finding is of interest for basic research and for clinical applications in subjects unable to perform cognitive tasks in the fMRI scanner.     

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.     

Functional magnetic resonance imaging reference phantom

Functional magnetic resonance imaging (fMRI) is widely used to pinpoint active brain areas. Changes in neuronal activity modulate the local blood oxygenation level, and the associated modulation of the magnetic field homogeneity can be detected with magnetic resonance imaging. Thus, the blood oxygenation level-dependent (BOLD) fMRI indirectly measures neuronal activity. Similar modulation of magnetic field homogeneity was here elicited by other means to generate a BOLD-like change in a new phantom constructed to provide reference activations during fMRI. Magnetic inhomogeneities were produced by applying current to coils located near the phantom containing 1.5 ml of Gd-doped water. The signal-to-noise ratio of the images, produced by gradient-recalled echo-planar imaging, varied between 104 and 107 at a selected voxel when the field was and was not inhomogenized, respectively. The contrast of signals between homogeneous and inhomogeneous conditions was generally stable, except in 3% of time points. During the periods of greatest deviations an observable change would have been detected in a simultaneously measured BOLD signal. Such changes could result from the imaging method or occur through glitches in hardware or alterations in the measurement environment. With identical measurement setups, the phantom could allow comparing intersession or intersubject brain activations.     

Mapping of the cerebral response to acetazolamide using graded asymmetric spin echo EPI

Cerebral vascular reactivity in different regions of the rat brain was quantitatively characterized by spatial and temporal measurements of blood oxygenation level-dependent (BOLD)-fMRI signals following intravenous administration of the carbonic anhydrase inhibitor acetazolamide: this causes cerebral vasodilatation through a cerebral extracellular acidosis that spares neuronal metabolism and vascular smooth muscle function, thus separating vascular and cerebral metabolic events. An asymmetric spin echo-echo planar imaging (ASE-EPI) pulse sequence sensitised images selectively to oxygenation changes in the microvasculature; use of a surface coil receiver enhanced image signal-to-noise ratios (SNRs). Image SNRs and hardware integrity were verified by incorporating quality assurance procedures; cardiorespiratory stability in the physiological preparations were monitored and maintained through the duration of the experiments. These conditions made it possible to apply BOLD contrast fMRI to map regional changes in cerebral perfusion in response to acetazolamide administration. Thus, fMRI findings demonstrated cerebral responses to acetazolamide that directly paralleled the known physiological actions of acetazolamide and whose time courses were similar through all regions of interest, consistent with acetazolamide's initial distribution in brain plasma, where it affects cerebral haemodynamics by acting at cerebral capillary endothelial cells. However, marked variations in the magnitude of the responses suggested relative perfusion deficits in the hippocampus and white matter regions correlating well with their relatively low vascularity and the known vulnerability of the hippocampus to ischaemic damage.     

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.     

Neurophysiological and functional MRI evidence of reorganization of cortical motor areas in cerebral arteriovenous malformation

Functional magnetic resonance imaging (fMRI) research has shown that brain arteriovenous malformations (AVMs) lead to reorganization of cortical motor areas. Since it is known that blood oxygenation level-dependent signal in fMRI may be influenced by the hemodynamic perturbation associated with the presence of the AVM, in the present study, a combined exploration with fMRI and transcranial magnetic stimulation was performed in a patient with a right rolandic AVM in order to explore the relationship between neuronal and hemodynamic activity. The combined protocol of investigation adopted in this study was able to provide significant information regarding neuronal activity of the different cortical areas that partake to post-lesional reorganization.     

BOLD signal and vessel dynamics: a hierarchical cluster analysis

The aim of the present study was to analyze blood oxygenation level-dependent (BOLD) signal variation during an apnea-based task in order to assess the capability of a functional magnetic resonance imaging (fMRI) procedure to estimate cerebral vascular dynamic effects. We measured BOLD contrast by hierarchical cluster analysis in healthy subjects undergoing an fMRI experiment, in which the task paradigm was one phase of inspirational apnea (IA). By processing the time courses of the fMRI experiment, analysis was performed only on a subclass of all the possible signal patterns; basically, root mean square and absolute variation differences have been calculated. Considering the baseline value obtained by computing the mean value of the initial rest period as reference, particular voxels showed relative important variations during the IA task and during the recovery phase following the IA. We focused our interest on the signal response of voxels that would correspond mainly to white and gray matter regions and that also may be affected by the proximity of large venous vessels. The results are presented as maps of space-temporal distribution of time series variations with two levels of hierarchical clustering among voxels with low to high initial response. Furthermore, we have presented a clustering of the signal response delay, conducting to a partition and identification of specified brain sites.     

On the origin of respiratory artifacts in BOLD-EPI of the human brain

BOLD-based functional MRI (fMRI) can be used to explicitly measure hemodynamic aspects and functions of human neuro-physiology. As fMRI measures changes in regional cerebral blood flow and volume as well as blood oxygenation, rather than neuronal brain activity directly, other processes that may change the above parameters have to be examined closely to assess sensitivity and specificity of fMRI results. Physiological processes that can cause artifacts include cardiac action, breathing and vasomotion. Although there has been substantial research on physiological artifacts and appropriate compensation methods, controversy still remains on the mechanisms that cause the fMRI signal fluctuations. Respiratory-correlated fluctuations may either be induced by changes of the magnetic field homogeneity due to moving organs, intra-thoracic pressure differences, respiration-dependent vasodilation or oxygenation differences. The aim of this study was to characterize the impact of different breathing patterns by varying respiration frequency and/or tidal volume on EPI time courses of the resting human brain. The amount of respiration-related oscillations during three respiration patterns was quantified, and statistically significant differences were obtained in white matter only: p < 0.03 between 6 vs. 12 ml/kg body weight end tidal volume at a respiration frequency of 15/min, p < 0.03 between 12 vs. 6 ml/kg body weight and 15 vs. 10 respiration cycles/min. There was no significant difference between 15 vs. 10 respiration cycles/min at an end tidal volume of 6 ml/kg body weight (p = 0.917). In addition, the respiration-affected brain regions were very similar with EPI readout in the a-p and l-r direction. Based on our results and published literature we hypothesize that venous oxygenation oscillations due to changing intra-thoracic pressure represent a major factor for respiration-related signal fluctuations and increase significantly with increasing end tidal volume in white matter only.     

EEG-fMRI mapping of asymmetrical delta activity in a patient with refractory epilepsy is concordant with the epileptogenic region determined by intracranial EEG

We studied a patient with refractory focal epilepsy using continuous EEG-correlated fMRI. Seizures were characterized by head turning to the left and clonic jerking of the left arm, suggesting a right frontal epileptogenic region. Interictal EEG showed occasional runs of independent nonlateralized slow activity in the delta band with right frontocentral dominance and had no lateralizing value. Ictal scalp EEG had no lateralizing value. Ictal scalp EEG suggested right-sided central slow activity preceding some seizures. Structural 3-T MRI showed no abnormality. There was no clear epileptiform abnormality during simultaneous EEG-fMRI. We therefore modeled asymmetrical EEG delta activity at 1-3 Hz near frontocentral electrode positions. Significant blood oxygen level-dependent (BOLD) signal changes in the right superior frontal gyrus correlated with right frontal oscillations at 1-3 Hz but not at 4-7 Hz and with neither of the two frequency bands when derived from contralateral or posterior electrode positions, which served as controls. Motor fMRI activations with a finger-tapping paradigm were asymmetrical: they were more anterior for the left hand compared with the right and were near the aforementioned EEG-correlated signal changes. A right frontocentral perirolandic seizure onset was identified with a subdural grid recording, and electric stimulation of the adjacent contact produced motor responses in the left arm and after discharges. The fMRI localization of the left hand motor and the detected BOLD activation associated with modeled slow activity suggest a role for localization of the epileptogenic region with EEG-fMRI even in the absence of clear interictal discharges.     

Fast spectroscopic imaging strategies for potential applications in fMRI

Technical aspects of two general fast spectroscopic imaging (SI) strategies, one based on gradient echo trains and the other on spin echo trains, are reviewed within the context of potential applications in the field of functional magnetic resonance imaging (fMRI). Fast spectroscopic imaging of water may prove useful for identifying mechanisms underlying the blood oxygenation level dependence (BOLD) of the water signal during brain activation studies. Reasonably rapid mapping of changes in proton signals from brain metabolites, like lactate, creatine or even neurotransmitter associated metabolites like GABA, is substantially more challenging but technically feasible particularly as higher field strengths become available. Fast spectroscopic methods directed towards the ³¹P signals from phosphocreatine (PCr) and adenosine tri-phosphates (ATP) are also technically feasible and may prove useful for studying cerebral energetics within fMRI contexts.     

Visual cortex reactivity in sedated children examined with perfusion MRI (FAIR)

Sleeping and sedated children can respond to visual stimulation with a decrease in blood oxygenation level dependent (BOLD) functional MRI signal response. The contribution of metabolic and hemodynamic parameters to this inverse signal response is incompletely understood. It has been hypothesized that it is caused by a relatively greater increase of oxygen consumption compared to rCBF (regional cerebral blood flow) increase. We studied the rCBF changes during visual stimulation in four sedated children, aged 4-71 months, and four alert adults, with an arterial water spin labeling technique (FAIR) and BOLD fMRI in a 1.5T MR scanner. In the children, FAIR signal decreased by a mean of 0.96% (range 0.77-1.05) of the baseline periods of the non-selective images, while BOLD signal decreased by 2.03% (range 1.99-2.93). In the adults, FAIR and BOLD signal increased by 0.88% (range 0.8-0.99) and 2.63% (range 1.99-2.93), respectively. Thus, in the children, an rCBF increase could not be detected by perfusion MRI, but indications of a FAIR signal decrease were found. An rCBF decrease in the primary visual cortex during stimulation has not been reported previously, but it is a possible explanation for the negative BOLD response. Future studies will have to address if this response pattern is a consequence of age or sleep/sedation.     

Component structure of event-related fMRI responses in the different neurovascular compartments

In most functional magnetic resonance imaging (fMRI) studies, brain activity is localized by observing changes in the blood oxygenation level-dependent (BOLD) signal that are believed to arise from capillaries, venules and veins in and around the active neuronal population. However, the contribution from veins can be relatively far downstream from active neurons, thereby limiting the ability of BOLD imaging methods to precisely pinpoint neural generators. Hemodynamic measures based on apparent diffusion coefficients (ADCs) have recently been used to identify more upstream functional blood flow changes in the capillaries, arterioles and arteries. In particular, we recently showed that, due to the complementary vascular sensitivities of ADC and BOLD signals, the voxels conjointly activated by both measures may identify the capillary networks of the active neuronal areas. In this study, we first used simultaneously acquired ADC and BOLD functional imaging signals to identify brain voxels activated by ADC only, by both ADC and BOLD and by BOLD only, thereby delineating voxels relatively dominated by the arterial, capillary, and draining venous neurovascular compartments, respectively. We then examined the event-related fMRI BOLD responses in each of these delineated neurovascular compartments, hypothesizing that their event-related responses would show different temporal componentries. In the regions activated by both the BOLD and ADC contrasts, but not in the BOLD-only areas, we observed an initial transient signal reduction (an initial dip), consistent with the local production of deoxyhemoglobin by the active neuronal population. In addition, the BOLD-ADC overlap areas and the BOLD-only areas showed a clear poststimulus undershoot, whereas the compartment activated by only ADC did not show this component. These results indicate that using ADC contrast in conjunction with BOLD imaging can help delineate the various neurovascular compartments, improve the localization of active neural populations, and provide insight into the physiological mechanisms underlying the hemodynamic signals.     

Percept-related activity in the human somatosensory system: functional magnetic resonance imaging studies

In this paper, we review blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging (fMRI) studies addressing the neural correlates of touch, thermosensation, pain and the mechanisms of their cognitive modulation in healthy human subjects. There is evidence that fMRI signal changes can be elicited in the parietal cortex by stimulation of single mechanoceptive afferent fibers at suprathreshold intensities for conscious perception. Positive linear relationships between the amplitude or the spatial extents of BOLD fMRI signal changes, stimulus intensity and the perceived touch or pain intensity have been described in different brain areas. Some recent fMRI studies addressed the role of cortical areas in somatosensory perception by comparing the time course of cortical activity evoked by different kinds of stimuli with the temporal features of touch, heat or pain perception. Moreover, parametric single-trial functional MRI designs have been adopted in order to disentangle subprocesses within the nociceptive system.

Available evidence suggest that studies that combine fMRI with psychophysical methods may provide a valuable approach for understanding complex perceptual mechanisms and top-down modulation of the somatosensory system by cognitive factors specifically related to selective attention and to anticipation. The brain networks underlying somatosensory perception are complex and highly distributed. A deeper understanding of perceptual-related brain mechanisms therefore requires new approaches suited to investigate the spatial and temporal dynamics of activation in different brain regions and their functional interaction.     

Functional magnetic resonance imaging of tonic pain and vasopressor effects in rats

In functional magnetic resonance imaging (fMRI) studies, an elevation in blood pressure (BP) in individuals with a poor autoregulatory response may increase cerebral blood flow, potentially enhancing the blood oxygenation level dependent response. To investigate the role of BP changes, the cerebral activation to either tonic pain or the infusion of the vasopressor norepinephrine was correlated with the accompanying BP changes in alpha-chloralose anesthetized rats. Immediately after formalin (2%) injection into the forepaw, fMRI detected an activation that was correlated with the BP increase and additional activations that were independent of blood pressure changes 5–40 minutes later. The activation detected with the administration of the vasopressor norepinephrine, which does not cross the blood-brain barrier was correlated to both the amount and rate of increase in BP. The response ranged from being sparse, localized within cortex or widespread during modest, moderate or severe elevations in BP, respectively. The cerebral circulatory effects of hypertension should be considered as contributing to changes in cerebral blood oxygenation in fMRI studies involving increases in BP.     

BOLD responses to trigeminal nerve stimulation

The current study investigates a new model of barrel cortex activation using stimulation of the infraorbital branch of the trigeminal nerve. A robust and reproducible activation of the rat barrel cortex was obtained following trigeminal nerve stimulation. Blood oxygen level-dependent (BOLD) effects were obtained in the primary somatosensory barrel cortex (S1BF), the secondary somatosensory cortex (S2) and the motor cortex. These cortical areas were reached from afferent pathways from the trigeminal ganglion, the trigeminal nuclei and thalamic nuclei from which neurons project their axons upon whisker stimulation. The maximum BOLD responses were obtained for a stimulus frequency of 1 Hz, a stimulus pulse width of 100 μs and for current intensities between 1.5 and 3 mA. The BOLD response was nonlinear as a function of frequency and current intensity. Additionally, modeling BOLD responses in the rat barrel cortex from separate cerebral blood flow (CBF) and cerebral metabolic rate of oxygen (CMRO₂) measurements showed good agreement with the shape and amplitude of measured BOLD responses as a function of stimulus frequency and will potentially allow to identify the sources of BOLD nonlinearities. Activation of the rat barrel cortex using trigeminal nerve stimulation will contribute to the interpretation of the BOLD signals from functional magnetic resonance imaging studies.     

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.     

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.     

Nonlinear blood oxygen level-dependent responses for transient activations and deactivations in V1 - insights into the hemodynamic response function with the balloon model

We report studies of the nonlinear nature of blood oxygen level-dependent (BOLD) responses to short transient deactivations in human visual cortex. Both functional magnetic resonance imaging (fMRI) and near-infrared spectroscopy (NIRS) have been used to compare and contrast the hemodynamic response functions (HRFs) associated with transient activation and deactivation in primary visual cortex. We show that signal decreases for short duration deactivations are smaller than corresponding signal increases in activation studies. Moreover, the standard balloon model of BOLD effects may be modified to account for the observed nonlinear nature of deactivations by appropriate changes to simple hemodynamic parameters without recourse to new assumptions about the nature of the coupling between activity and oxygen use.