Stimulus-dependent BOLD and perfusion dynamics in human V1

Biophysical model for integrating neuronal activity, EEG, fMRI and metabolism

Our goal is to model the coupling between neuronal activity, cerebral metabolic rates of glucose and oxygen consumption, cerebral blood flow (CBF), electroencephalography (EEG) and blood oxygenation level-dependent (BOLD) responses. In order to accomplish this, two previous models are coupled: a metabolic/hemodynamic model (MHM) for a voxel, linking BOLD signals and neuronal activity, and a neural mass model describing the neuronal dynamics within a voxel and its interactions with voxels of the same area (short-range interactions) and other areas (long-range interactions). For coupling both models, we take as the input to the BOLD model, the number of active synapses within the voxel, that is, the average number of synapses that will receive an action potential within the time unit. This is obtained by considering the action potentials transmitted between neuronal populations within the voxel, as well as those arriving from other voxels. Simulations are carried out for testing the integrated model. Results show that realistic evoked potentials (EP) at electrodes on the scalp surface and the corresponding BOLD signals for each voxel are produced by the model. In another simulation, the alpha rhythm was reproduced and reasonable similarities with experimental data were obtained when calculating correlations between BOLD signals and the alpha power curve. The origin of negative BOLD responses and the characteristics of EEG, PET and BOLD signals in Alzheimer's disease were also studied.

Flexible statistical modelling detects clinical functional magnetic resonance imaging activation in partially compliant subjects

Clinical functional magnetic resonance imaging (fMRI) occasionally fails to detect significant activation, often due to variability in task performance. The present study seeks to test whether a more flexible statistical analysis can better detect activation, by accounting for variance associated with variable compliance to the task over time. Experimental results and simulated data both confirm that even at 80% compliance to the task, such a flexible model outperforms standard statistical analysis when assessed using the extent of activation (experimental data), goodness of fit (experimental data), and area under the operator characteristic curve (simulated data). Furthermore, retrospective examination of 14 clinical fMRI examinations reveals that in patients where the standard statistical approach yields activation, there is a measurable gain in model performance in adopting the flexible statistical model, with little or no penalty in lost sensitivity. This indicates that a flexible model should be considered, particularly for clinical patients who may have difficulty complying fully with the study task.

Improved spatial localization of post-stimulus BOLD undershoot relative to positive BOLD

The negative blood oxygenation level-dependent (BOLD) signal following the cessation of stimulation (post-stimulus BOLD undershoot) is observed in functional magnetic resonance imaging (fMRI) studies. However, its spatial characteristics are unknown. To investigate this, gradient-echo BOLD fMRI in response to visual stimulus was obtained in isoflurane-anesthetized cats at 9.4 T. Since the middle cortical layer (layer 4) is known to have the highest metabolic and cerebral blood volume (CBV) responses, images were obtained to view the cortical cross-section. Robust post-stimulus BOLD undershoot was observed in all studies, and lasted longer than 30 s after the cessation of 40-60 s stimulation. The magnitude of post-stimulus BOLD undershoot was linearly dependent on echo time with little intercept when extrapolating to TE = 0, indicating that the T2* change is the major cause of the BOLD undershoot. The post-stimulus BOLD undershoot was observed within the cortex and near the surface of the cortex, while the prolonged CBV elevation was observed only at the middle of the cortex. Within the cortex, the largest post-stimulus undershoot was detected at the middle of the cortex, similar to the CBV increase during the stimulation period. Our findings demonstrate that, even though there is significant contribution from pial vessel signals, the post-stimulus undershoot BOLD signal is useful to improve the spatial localization of fMRI to active cortical sites.

Topography of Cortical Activation Differs for Fundamental and Harmonic Frequencies of the Steady-State Visual-Evoked Responses. An EEG and PET H215O Study

In humans, visual flicker stimuli of graded frequency (2-90 Hz) elicit an electroencephalographic (EEG) steady-state visual-evoked response (SSVER) with the same fundamental frequency as the stimulus and, in addition, a series of harmonic responses. The fundamental component of the SSVER is generated by increased synaptic activity in primary visual cortex (V1). We set out to determine the cortical origin of the harmonic responses in humans. For this purpose, we recorded the SSVERs at 5 different frequencies (5, 10, 15, 25, and 40 Hz) and measured regional cerebral blood flow (rCBF) with positron emission tomography-H(2)(15)O at rest and during visual stimulation at the same frequencies. The rCBF contrast weighted by the amplitude of the SSVERs first harmonics showed activation of a swath of cortex perpendicular to V1, including mostly the inferior half of the parieto-occipital sulcus. This area overlapped minimally with the primary visual cortex activated by the fundamental frequency. A different method, estimating EEG cortical source current density with low-resolution brain electromagnetic tomography, gave the same results. Our finding suggests that the inferior portion of the banks of the parieto-occipital sulci contains association visual cortex involved in the processing of stimuli that can be as simple as a flickering light source.

Detrimental effects of BOLD signal in arterial spin labeling fMRI at high field strength

Arterial spin labeling (ASL) MRI is a useful technique for noninvasive measurement of cerebral blood flow (CBF) in humans. High field strength provides a unique advantage for ASL because of longer blood T(1) relaxation times, making this technique a promising quantitative approach for functional brain mapping. However, higher magnetic field also introduces new challenges. Here it is shown that the CBF response determined using ASL functional MRI (fMRI) at 3.0 T contains significant contamination from blood-oxygenation-level-dependent (BOLD) effects. Due to interleaved acquisitions of label and control images, difference in blood oxygenation status between these two scans can cause incomplete cancellation of the static signal upon image subtraction, resulting in a BOLD-related artifact in the estimated CBF hemodynamics. If not accounted for, such an effect can complicate the interpretation of the ASL results, e.g., causing a delayed onset and offset of the response, or inducing an artifactual poststimulus undershoot. The BOLD contribution also decreases the sensitivity of ASL-based fMRI. Correction methods are proposed to reduce the artifact, giving increased number of activated voxels (18+/-5%, P=0.006) and more accurate estimation of CBF temporal characteristics.

Functional MRI today

Most brain imaging researchers would agree with the assertion that functional MRI (fMRI) is progressing. Since fMRI began in 1991, the number of people, papers, and abstracts related to fMRI has been increasing; the technology and methodology has shown advances in robustness and sophistication; the physiology of the signal is better understood; and, even though it hasn't yet made significant headway into the clinical setting, applications are widening. Questions that stem from this optimistic and perhaps overly general set of observations include those that ask what the ultimate theoretical and practical limits of fMRI are and how close are we to approaching these limits. In this commentary, I attempt to provide a snapshot of fMRI as it exists at the end of 2005, and to give a clear impression that not only are we progressing by "dotting the i's and crossing the t's" but that fundamental changes in fMRI methodology and processing are being put forth as the field matures.

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.

The spatial dependence of the poststimulus undershoot as revealed by high-resolution BOLD- and CBV-weighted fMRI

The hemodynamic response to neural activity consists of changes in blood flow, blood volume and oxygen metabolism. Changes in the vascular state after sensory stimulation have different spatial and temporal characteristics in the brain. This has been shown using imaging techniques, such as BOLD functional magnetic resonance imaging (fMRI), which monitor vascular changes once the stimulus is turned on, and the eventual return to baseline levels, once the stimulus is turned off. The BOLD fMRI signal during sensory stimulation has been well characterized and modeled in terms of the spatial and temporal characteristics of the vascular response. However, the return of the signals to baseline levels after sensory stimulation is not as well characterized. During this period, a poststimulus undershoot in the BOLD signal is observed. This poststimulus undershoot has been modeled and investigated to characterize the physiological mechanisms (cerebral blood flow (CBF), cerebral blood volume (CBV), and cerebral oxygen consumption) associated with the response. However, the data in the literature, which lack any spatially dependent information, appear to be contradictory in terms of the mechanisms associated with this poststimulus response. With a high spatial resolution cat model at 9.4 T, we show that CBV changes in the tissue persist once the stimulus is turned off, while CBV changes in the surface vessels quickly return to baseline levels, despite a concurrent undershoot in the BOLD signal in both the tissue and surface vessel areas. In addition, the BOLD data alone indicate that different physiological mechanisms regulate the poststimulus response in the tissue versus the surface vessel regions.

Comparison of spatial and temporal pattern for fMRI obtained with BOLD and arterial spin labeling

Functional magnetic resonance imaging (fMRI) is presently either performed using blood oxygenation level-dependent (BOLD) contrast or using cerebral blood flow (CBF), measured with arterial spin labeling (ASL) technique. The present fMRI study aimed to provide practical hints to favour one method over the other. It involved three different acquisition methods during visual checkerboard stimulation on nine healthy subjects: 1) CBF contrast obtained from ASL, 2) BOLD contrast extracted from ASL and 3) BOLD contrast from Echo planar imaging. Previous findings were replicated; i) no differences between the three measurements were found in the location of the activated region; ii) differences were found in the temporal characteristics of the signals and iii) BOLD has significantly higher sensitivity than ASL perfusion. ASL fMRI was favoured when the investigation demands for perfusion and task related signal changes. BOLD fMRI is more suitable in conjunction with fast event related design.

A temporal comparison of BOLD, ASL, and NIRS hemodynamic responses to motor stimuli in adult humans

In this study, we have preformed simultaneous near-infrared spectroscopy (NIRS) along with BOLD (blood oxygen level dependent) and ASL (arterial spin labeling)-based fMRI during an event-related motor activity in human subjects in order to compare the temporal dynamics of the hemodynamic responses recorded in each method. These measurements have allowed us to examine the validity of the biophysical models underlying each modality and, as a result, gain greater insight into the hemodynamic responses to neuronal activation. Although prior studies have examined the relationships between these two methodologies through similar experiments, they have produced conflicting results in the literature for a variety of reasons. Here, by employing a short-duration, event-related motor task, we have been able to emphasize the subtle temporal differences between the hemodynamic parameters with a high contrast-to-noise ratio. As a result of this improved experimental design, we are able to report that the fMRI measured BOLD response is more correlated with the NIRS measure of deoxy-hemoglobin (R = 0.98; P < 10(-20)) than with oxy-hemoglobin (R = 0.71), or total hemoglobin (R = 0.53). This result was predicted from the theoretical grounds of the BOLD response and is in agreement with several previous works [Toronov, V.A.W., Choi, J.H., Wolf, M., Michalos, A., Gratton, E., Hueber, D., 2001. "Investigation of human brain hemodynamics by simultaneous near-infrared spectroscopy and functional magnetic resonance imaging." Med. Phys. 28 (4) 521-527.; MacIntosh, B.J., Klassen, L.M., Menon, R.S., 2003. "Transient hemodynamics during a breath hold challenge in a two part functional imaging study with simultaneous near-infrared spectroscopy in adult humans". NeuroImage 20 1246-1252.; Toronov, V.A.W., Walker, S., Gupta, R., Choi, J.H., Gratton, E., Hueber, D., Webb, A., 2003. "The roles of changes in deoxyhemoglobin concentration and regional cerebral blood volume in the fMRI BOLD signal" Neuroimage 19 (4) 1521-1531]. These data have also allowed us to examine more detailed measurement models of the fMRI signal and comment on the roles of the oxygen saturation and blood volume contributions to the BOLD response. In addition, we found high correlation between the NIRS measured total hemoglobin and ASL measured cerebral blood flow (R = 0.91; P < 10(-10)) and oxy-hemoglobin with flow (R = 0.83; P < 10(-05)) as predicted by the biophysical models. Finally, we note a significant amount of cross-modality, correlated, inter-subject variability in amplitude change and time-to-peak of the hemodynamic response. The observed co-variance in these parameters between subjects is in agreement with hemodynamic models and provides further support that fMRI and NIRS have similar vascular sensitivity.

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.

Investigating the post-stimulus undershoot of the BOLD signal--a simultaneous fMRI and fNIRS study

Measuring the hemodynamic response with functional magnetic resonance imaging (fMRI) together with functional near-infrared spectroscopy (fNIRS) may overcome limitations of single-method approaches. Accordingly, we measured the event-related hemodynamic response with both imaging methods simultaneously in young subjects during visual stimulation. An intertrial interval of 60 s was chosen to include the prolonged post-stimulus undershoot of the blood oxygenation level dependent (BOLD) signal. During visual stimulation, the BOLD signal, oxy-, and total hemoglobin (Hb) increased, whereas deoxy-Hb decreased. The post-stimulus period was characterized by an undershoot of the BOLD signal, oxy-Hb, and an overshoot of deoxy-Hb. Total Hb as measured by fNIRS returned to baseline immediately after the end of stimulation. Results suggest that the post-stimulus events as measured by fNIRS are dominated by a prolonged high-level oxygen consumption in the microvasculature. The contribution of a delayed return of blood volume to the BOLD post-stimulus undershoot in post-capillary veins as suggested by the Balloon and Windkessel models remains ambiguous. Temporal changes in the BOLD signal were highly correlated with deoxy-Hb, with lower correlation values for oxy- and total Hb. Furthermore, data show that fNIRS covers the outer 1 cm of the brain cortex. These results were confirmed by simultaneous fMRI/fNIRS measurements during rest. In conclusion, multimodal imaging approaches may contribute to the understanding of neurovascular coupling.

Simultaneous perfusion and blood-oxygenation-level-dependent measurements using single-shot interleaved z-shim echo-planar imaging

Single-shot interleaved z-shim EPI (SSIZS-EPI) was extended to a simultaneous perfusion and blood-oxygenation-level-dependent (BOLD) imaging technique that reduces susceptibility-induced signal loss while preserving rapid image acquisition. Experiments on human brains showed that images acquired with this technique had improved signal-to-noise ratio in the inferior prefrontal, meso-, and lateral-temporal lobes compared with a conventional EPI. Perfusion maps obtained from the SSIZS-EPI images at resting state illustrated substantial signal recovery in these brain areas. Perfusion and BOLD images collected with a sensorimotor paradigm demonstrated the feasibility of the technique to simultaneously measure cerebral blood flow and blood oxygenation signals associated with brain activation. Functional experiments with a neuropsychiatric paradigm showed increased brain activities in the periamygdalar regions in both perfusion and BOLD maps, consistent with a previous H(2) (15)O PET study. The proposed technique, with its advantages of reducing susceptibility artifacts and fast scanning speed, would be useful for obtaining more reliable measurements of functional signals, particularly in the brain regions with field inhomogeneities.

The BOLD onset transient: identification of novel functional differences in schizophrenia

Blood oxygen level dependent (BOLD) signals characteristically exhibit an overshoot (transient signal increase) at the beginning of fMRI task blocks. This onset transient has often been overlooked as an independent measure of neuronal activity, but it may represent unique functional processes. We examined onset transient responses in normal subjects and individuals with schizophrenia performing three cognitive tasks. These analyses revealed a regionally specific and task specific attenuation of the onset transient in individuals with schizophrenia during performance of a working memory task. Furthermore, this attenuation was often not accompanied by a corresponding population difference in the sustained response, and is missed through conventional fMRI analysis techniques. Relevance of these findings to both an interpretation of the onset transient and the pathology of schizophrenia are discussed.

Transient BOLD responses at block transitions

Block-design fMRI responses include sustained components present for the duration of each task block as well as transient components at the beginning and end of each block. Almost all prior block-design fMRI studies have focused on the sustained response components while the transient responses at block transitions have been largely ignored. These transients, therefore, remain poorly characterized. We here present a systematic study of block-transition transient responses obtained using four widely divergent tasks. We characterize transient response topography and examine the extent to which these responses vary across different tasks and between block onset and offset. Our analysis reveals that certain regions show transient responses regardless of task or transition type. However, our analysis also shows that specific task state transitions give rise to transient responses with unique spatial profiles. Relevance of the current findings to studies of exogenous attention, task shifting, and the BOLD overshoot is discussed.

A signal processing model for arterial spin labeling functional MRI

A model of the signal path in arterial spin labeling (ASL)-based functional magnetic resonance imaging (fMRI) is presented. Three subtraction-based methods for forming a perfusion estimate are considered and shown to be specific cases of a generalized estimate consisting of a modulator followed by a low pass filter. The performance of the methods is evaluated using the signal model. Contamination of the perfusion estimate by blood oxygenation level dependent contrast (BOLD) is minimized by using either sinc subtraction or surround subtraction for block design experiments and by using pair-wise subtraction for randomized event-related experiments. The subtraction methods all tend to decorrelate the 1/f type low frequency noise often observed in fMRI experiments. Sinc subtraction provides the flattest noise power spectrum at low frequencies, while pair-wise subtraction yields the narrowest autocorrelation function. The formation of BOLD estimates from the ASL data is also considered and perfusion weighting of the estimates is examined using the signal model.

Coupling of cerebral blood flow and oxygen consumption during physiological activation and deactivation measured with fMRI

The physiological basis of the blood oxygenation level dependent (BOLD) signal and its dependence on baseline cerebral blood flow (CBF) were investigated by comparing responses to a visual stimulus after physiological changes of the baseline. Eight human subjects were imaged with 3 and 4 T MRI scanners, and both BOLD signal and CBF were simultaneously measured. Subjects viewed a flickering radial checkerboard in a block design experiment, alternating between eyes open or closed during the off periods. Compared to a baseline state with eyes open in a darkened room, substantial deactivation (average change: 2.9 +/- 0.3% BOLD, 22 +/- 2.1% CBF) in the occipital cortex was observed when the eyes were closed. The absolute response during stimulation (average change: 4.4 +/- 0.4% BOLD, 36.3 +/- 3.1% CBF) was independent of the preceding resting condition. We estimated the fractional change in CBF to be approximately 2.2 +/- 0.15 times greater than the fractional change in metabolic rate of oxygen (CMRO2). The changes in CBF and CMRO2 were consistently linearly coupled during activation and deactivation with CBF changes being between approximately 60% and 150% compared to baseline with eyes open. Relative to an assumed baseline oxygen extraction fraction (OEF) of 40%, the estimated OEF decreased to 33 +/- 1.4% during activation and increased to 46 +/- 1.2% during rest with eyes closed. In conclusion, we found that simply closing the eyes creates a large physiological deactivation in the visual cortex, and provides a robust paradigm for studying baseline effects in fMRI. In addition, we propose a feed-forward model for neurovascular coupling which accounts for the changes in OEF seen following baseline changes, including both the current physiological perturbations as well as previously reported pharmacologically induced changes.

Perfusion imaging using arterial spin labeling

Arterial spin labeling is a magnetic resonance method for the measurement of cerebral blood flow. In its simplest form, the perfusion contrast in the images gathered by this technique comes from the subtraction of two successively acquired images: one with, and one without, proximal labeling of arterial water spins after a small delay time. Over the last decade, the method has moved from the experimental laboratory to the clinical environment. Furthermore, numerous improvements, ranging from new pulse sequence implementations to extensive theoretical studies, have broadened its reach and extended its potential applications. In this review, the multiple facets of this powerful yet difficult technique are discussed. Different implementations are compared, the theoretical background is summarized, and potential applications of various implementations in research as well as in the daily clinical routine are proposed. Finally, a summary of the new developments and emerging techniques in this field is provided.

The effect of task block arrangement on the detectability of activation in fMRI

The effect of task block arrangements on the detection of brain activation was investigated. Sessions of functional magnetic resonance imaging (fMRI) including the same number of two different task conditions but with different arrangements were compared. The two task conditions were, A) Ellipse-shaped black and white checkerboard flicker stimulation at 4.2 Hz covering the bilateral visual field, and B) the same flicker stimuli covering only the left visual field. In the rest blocks (0), the subjects looked at a fixation point. Four different task block arrangements were compared, 1) A0 (0A0A0A0) and B0 (0B0B0B0), 2) A0B0 (0A0B0A0B0A0B0), 3) AB0 (0AB0AB0AB0) and 4) AB (0ABABAB). Bilateral V1, V2, V3 and the left V5 were activated by condition A, and the right V1 and V2 by B. The activation in the left visual field by A0 was larger than in the other three conditions. In a differential analysis between conditions A and B, activation in the left V3 and V5 was declined by AB0 or AB. When rest blocks were located in the post-stimulus undershoot phase, the % signal change of the BOLD signal was emphasized, which caused augmented significance in the detection of the activity. It was indicated that the outcome of the activation map was influenced by the arrangement of task blocks, even though the same number of task blocks were repeated within the sessions. In fMRI studies, task conditions should be carefully compared within or across sessions considering the characteristics of hemodynamic response functions.

Detection and quantification of a wide range of fMRI temporal responses using a physiologically-motivated basis set

The temporal dynamics of fMRI responses can span a broad range, indicating a rich underlying physiology, but also posing a significant challenge for detection. For instance, in human auditory cortex, prolonged sound stimuli ( approximately 30 sec) can evoke responses ranging from sustained to highly phasic (i.e., characterized by prominent peaks just after sound onset and offset). In the present study, we developed a method capable of detecting a wide variety of responses, while simultaneously extracting information about individual response components, which may have different neurophysiological underpinnings. Specifically, we implemented the general linear model using a novel set of basis functions chosen to reflect temporal features of cortical fMRI responses. This physiologically-motivated basis set (the "OSORU" basis set) was tested against (1) the commonly employed "sustained-only" basis "set" (i.e., a single smoothed "boxcar" function), and (2) a sinusoidal basis set, which is capable of detecting a broad range of responses, but lacks a direct relationship to individual response components. On data that included many different temporal responses, the OSORU basis set performed far better overall than the sustained-only set, and as well or better than the sinusoidal basis set. The OSORU basis set also proved effective in exploring brain physiology. As an example, we demonstrate that the OSORU basis functions can be used to spatially map the relative amount of transient vs. sustained activity within auditory cortex. The OSORU basis set provides a powerful means for response detection and quantification that should be broadly applicable to any brain system and to both human and non-human species.

Discrepancies between BOLD and flow dynamics in primary and supplementary motor areas: application of the balloon model to the interpretation of BOLD transients

The blood-oxygen-level-dependent (BOLD) signal measured in the brain with functional magnetic resonance imaging (fMRI) during an activation experiment often exhibits pronounced transients at the beginning and end of the stimulus. Such transients could be a reflection of transients in the underlying neural activity, or they could result from transients in cerebral blood flow (CBF), cerebral metabolic rate of oxygen (CMRO2), or cerebral blood volume (CBV). These transients were investigated using an arterial spin labeling (ASL) method that allows simultaneous measurements of BOLD and CBF responses. Responses to a finger-tapping task (40-s stimulus, 80-s rest) were measured in primary motor area (M1) and supplementary motor area (SMA) in five healthy volunteers. In SMA, the average BOLD response was pronounced near the beginning and end of the stimulus, while in M1, the BOLD response was nearly flat. However, CBF responses in the two regions were rather similar, and did not exhibit the same transient features as the BOLD response in SMA. Because this suggests a hemodynamic rather than a neural origin for the transients of the BOLD response in SMA, we used a generalization of the balloon model to test the degree of hemodynamic transients required to produce the measured curves. Both data sets could be approximated with modest differences in the shapes of the CMRO2 and CBV responses. This study illustrates the utility and the limitations of using theoretical models combined with ASL techniques to understand the dynamics of the BOLD response.

Neuroimaging of direction-selective mechanisms for second-order motion

Psychophysical findings have revealed a functional segregation of processing for 1st-order motion (movement of luminance modulation) and 2nd-order motion (e.g., movement of contrast modulation). However neural correlates of this psychophysical distinction remain controversial. To test for a corresponding anatomical segregation, we conducted a new functional magnetic resonance imaging (fMRI) study to localize direction-selective cortical mechanisms for 1st- and 2nd-order motion stimuli, by measuring direction-contingent response changes induced by motion adaptation, with deliberate control of attention. The 2nd-order motion stimulus generated direction-selective adaptation in a wide range of visual cortical areas, including areas V1, V2, V3, VP, V3A, V4v, and MT+. Moreover, the pattern of activity was similar to that obtained with 1st-order motion stimuli. Contrary to expectations from psychophysics, these results suggest that in the human visual cortex, the direction of 2nd-order motion is represented as early as V1. In addition, we found no obvious anatomical segregation in the neural substrates for 1st- and 2nd-order motion processing that can be resolved using standard fMRI.

Four facets of a single brain: behaviour, cerebral blood flow/metabolism, neuronal activity and neurotransmitter dynamics

Is functional neuroimaging a royal way to understand brain function or is it a new phrenology without an exact understanding what we measure? After two decades of imaging revolution, more and more authors ask this question. Brain functions are multidimensional, which can be approached from the point of (1) behavioural measures, (2) brain activation as reflected by blood flow and metabolic changes, (3) electrical activity of cells and cell-populations, and (4) neurotransmitter dynamics (release, receptor binding and reuptake). Using imaging techniques, we must take into consideration that even during the simplest task all of these processes operate in a closely interacting manner. Therefore, before drawing final conclusions about brain functions on the basis of a single aspect of these mechanisms, we must clarify the exact relationship among them. In this paper, we address this issue in order to draw attention to a number of uncertainties and controversies in the relationship of the four facets of brain functions.

Temporal comparison of functional brain imaging with diffuse optical tomography and fMRI during rat forepaw stimulation

The time courses of oxyhaemoglobin ([HbO2]), deoxyhaemoglobin ([HbR]) and total haemoglobin ([HbT]) concentration changes following cortical activation in rats by electrical forepaw stimulation were measured using diffuse optical tomography (DOT) and compared to similar measurements performed previously with fMRI at 2.0 T and 4.7 T. We also explored the qualitative effects of varying stimulus parameters on the temporal evolution of the hemodynamic response. DOT images were reconstructed at a depth of 1.5 mm over a 1 cm square area from 2 mm anterior to bregma to 8 mm posterior to bregma. The measurement set included 9 sources and 16 detectors with an imaging frame rate of 10 Hz. Both DOT [HbR] and [HbO2] time courses were compared to the fMRI BOLD time course during stimulation, and the DOT [HbT] time course was compared to the fMRI cerebral plasma volume (CPV) time course. We believe that DOT and fMRI can provide similar temporal information for both blood volume and deoxyhaemoglobin changes, which helps to cross-validate these two techniques and to demonstrate that DOT can be useful as a complementary modality to fMRI for investigating the hemodynamic response to neuronal activity.

Signal nonlinearity in fMRI: a comparison between BOLD and MION

In this paper, we introduce a methodology for comparing the nonlinearities present in sets of time series using four different nonlinearity measures, one of which, the "delay vector variance" method, is a novel approach to the characterization of a time series. It is then applied to examine the difference in nonlinearity between functional magnetic resonance imaging (fMRI) signals that have been recorded using different contrast agents. Recently, an exogenous contrast agent, monocrystalline iron oxide particle (MION), has been introduced for fMRI, which has been shown to increase the functional sensitivity compared with the traditional blood oxygen level dependent (BOLD) technique. The resulting fMRI signals are influenced by cerebral blood volume, whereas the more traditionally recorded BOLD signals are influenced not only by cerebral blood volume, but also by the cerebral blood flow and the metabolic rate of oxygen. The proposed methodology is applied to address the question whether this difference in the number of physiological variables is reflected in a difference in the degree of nonlinearity. We therefore analyze two sets of fMRI signals, one from a BOLD and the other from a MION monkey study with similar experimental designs. In the neuroimaging context, the proposed nonlinearity analyses are different from those described in the literature, since no a priori model is assumed: rather than pinpointing the source(s) of nonlinearity, nonparametric analyses are performed on BOLD and MION fMRI signals. Furthermore, we introduce a strategy for analyzing a population of fMRI signals, rather than focusing the analysis on one signal, as is traditionally done in the domain of nonlinear signal processing. Our results show that, overall, the BOLD signals are more nonlinear in nature than the MION ones, which is in agreement with current hypotheses.

A model of the coupling between brain electrical activity, metabolism, and hemodynamics: application to the interpretation of functional neuroimaging

In order to improve the interpretation of functional neuroimaging data, we implemented a mathematical model of the coupling between membrane ionic currents, energy metabolism (i.e., ATP regeneration via phosphocreatine buffer effect, glycolysis, and mitochondrial respiration), blood-brain barrier exchanges, and hemodynamics. Various hypotheses were tested for the variation of the cerebral metabolic rate of oxygen (CMRO(2)): (H1) the CMRO(2) remains at its baseline level; (H2) the CMRO(2) is enhanced as soon as the cerebral blood flow (CBF) increases; (H3) the CMRO(2) increase depends on intracellular oxygen and pyruvate concentrations, and intracellular ATP/ADP ratio; (H4) in addition to hypothesis H3, the CMRO(2) progressively increases, due to the action of a second messenger. A good agreement with experimental data from magnetic resonance imaging and spectroscopy (MRI and MRS) was obtained when we simulated sustained and repetitive activation protocols using hypotheses (H3) or (H4), rather than hypotheses (H1) or (H2). Furthermore, by studying the effect of the variation of some physiologically important parameters on the time course of the modeled blood-oxygenation-level-dependent (BOLD) signal, we were able to formulate hypotheses about the physiological or biochemical significance of functional magnetic resonance data, especially the poststimulus undershoot and the baseline drift.

Human brain activity during illusory visual jitter as revealed by functional magnetic resonance imaging

One central problem in vision is how to compensate for retinal slip. A novel illusion (visual jitter) suggests the compensation mechanism is based solely on retinal motion. Adaptation to visual noise attenuates the motion signals used by the compensation stage, producing illusory jitter due to the undercompensation of retinal slip. Here, we investigated the neural substrate of retinal slip compensation during this illusion using high-field fMRI and retinotopic mapping in flattened cortical format. When jitter perception occurred, MR signal decreased in lower stages of the visual system but increased prominently in area MT+. In conclusion, visual areas as early as V1 are responsible for the adaptation stage, and MT+ is involved in the compensation stage. The present finding suggests the pathway from V1 to MT+ has an important role in stabilizing the visual world.

BOLD fMRI identifies limbic, paralimbic, and cerebellar activation during air hunger

Air hunger (uncomfortable urge to breathe) is a component of dyspnea (shortness of breath). Three human H(2)(15)O positron emission tomography (PET) studies have identified activation of phylogenetically ancient structures in limbic and paralimbic regions during dyspnea. Other studies have shown activation of these structures during other sensations that alert the organism to urgent homeostatic imbalance: pain, thirst, and hunger for food. We employed blood oxygen level dependent (BOLD) functional magnetic resonance imaging (fMRI) to examine activation during air hunger. fMRI conferred several advantages over PET: enhanced signal-to-noise, greater spatial resolution, and lack of ionizing radiation, enabling a greater number of trials in each subject. Six healthy men and women were mechanically ventilated at 12-14 breaths/min. The primary experiment was conducted at mean end-tidal PCO(2) of 41 Torr. Moderate to severe air hunger was evoked during 42-s epochs of lower tidal volume (mean = 0.75 L). Subjects described the sensation as "like breath-hold," "urge to breathe," and "starved for air." In the baseline condition, air hunger was consistently relieved by epochs of higher tidal volume (mean = 1.47 L). A control experiment in the same subjects under a background of mild hypocapnia (mean end-tidal PCO(2) = 33 Torr) employed similar tidal volumes but did not evoke air hunger, controlling for stimulus variables not related to dyspnea. During each experiment, we maintained constant end-tidal PCO(2) and PO(2) to avoid systematic changes in global cerebral blood flow. Whole-brain images were acquired every 5 s (T2*, 56 slices, voxel resolution 3 x 3 x 3 mm). Activations associated with air hunger were determined using voxel-based interaction analysis of covariance that compared data between primary and control experiments (SPM99). We detected activations not seen in the earlier PET study using a similar air hunger stimulus (Banzett et al. 2000). Limbic and paralimbic loci activated in the present study were within anterior insula (seen in all 3 published studies of dyspnea), anterior cingulate, operculum, cerebellum, amygdala, thalamus, and basal ganglia. Elements of frontoparietal attentional networks were also identified. The consistency of anterior insular activation across subjects in this study and across published studies suggests that the insula is essential to dyspnea perception, although present data suggest that the insula acts in concert with a larger neural network.

Sound repetition rate in the human auditory pathway: representations in the waveshape and amplitude of fMRI activation

Sound repetition rate plays an important role in stream segregation, temporal pattern recognition, and the perception of successive sounds as either distinct or fused. This study was aimed at elucidating the neural coding of repetition rate and its perceptual correlates. We investigated the representations of rate in the auditory pathway of human listeners using functional magnetic resonance imaging (fMRI), an indicator of population neural activity. Stimuli were trains of noise bursts presented at rates ranging from low (1-2/s; each burst is perceptually distinct) to high (35/s; individual bursts are not distinguishable). There was a systematic change in the form of fMRI response rate-dependencies from midbrain to thalamus to cortex. In the inferior colliculus, response amplitude increased with increasing rate while response waveshape remained unchanged and sustained. In the medial geniculate body, increasing rate produced an increase in amplitude and a moderate change in waveshape at higher rates (from sustained to one showing a moderate peak just after train onset). In auditory cortex (Heschl's gyrus and the superior temporal gyrus), amplitude changed somewhat with rate, but a far more striking change occurred in response waveshape-low rates elicited a sustained response, whereas high rates elicited an unusual phasic response that included prominent peaks just after train onset and offset. The shift in cortical response waveshape from sustained to phasic with increasing rate corresponds to a perceptual shift from individually resolved bursts to fused bursts forming a continuous (but modulated) percept. Thus at high rates, a train forms a single perceptual "event," the onset and offset of which are delimited by the on and off peaks of phasic cortical responses. While auditory cortex showed a clear, qualitative correlation between perception and response waveshape, the medial geniculate body showed less correlation (since there was less change in waveshape with rate), and the inferior colliculus showed no correlation at all. Overall, our results suggest a population neural representation of the beginning and the end of distinct perceptual events that is weak or absent in the inferior colliculus, begins to emerge in the medial geniculate body, and is robust in auditory cortex.

Repeated fMRI using iron oxide contrast agent in awake, behaving macaques at 3 Tesla

Iron oxide contrast agents have been employed extensively in anesthetized rodents to enhance fMRI sensitivity and to study the physiology of cerebral blood volume (CBV) in relation to blood oxygen level-dependent (BOLD) signal following neuronal activation. This study quantified the advantages of exogenous agent for repeated neuroimaging in awake, nonhuman primates using a clinical 3 Tesla scanner. A monocrystalline iron oxide nanoparticle (MION) solution was injected at iron doses of 8 to 10 mg/kg in two macaque monkeys. Adverse behavioral effects due to contrast agent were not observed in either monkey using cumulative doses in excess of 60 mg/kg. Relative to BOLD imaging at 3 Tesla, MION increased functional sensitivity by an average factor of 3 across the brain for a stimulus of long duration. Rapid stimulus presentation attenuated MION signal changes more than BOLD signal changes, due to the slower time constant of the blood volume response relative to BOLD signal. Overall, the contrast agent produced a dramatic improvement in functional brain imaging results in the awake, behaving primate at this field strength. (c) 2002 Elsevier Science (USA).

What does fMRI tell us about neuronal activity?

In recent years, cognitive neuroscientists have taken great advantage of functional magnetic resonance imaging (fMRI) as a non-invasive method of measuring neuronal activity in the human brain. But what exactly does fMRI tell us? We know that its signals arise from changes in local haemodynamics that, in turn, result from alterations in neuronal activity, but exactly how neuronal activity, haemodynamics and fMRI signals are related is unclear. It has been assumed that the fMRI signal is proportional to the local average neuronal activity, but many factors can influence the relationship between the two. A clearer understanding of how neuronal activity influences the fMRI signal is needed if we are correctly to interpret functional imaging data.

Simultaneous BOLD/perfusion measurement using dual-echo FAIR and UNFAIR: sequence comparison at 1.5T and 3.0T

Functional MRI (fMRI) studies designed for simultaneously measuring Blood Oxygenation Level Dependent (BOLD) and Cerebral Blood Flow (CBF) signal often employ the standard Flow Alternating Inversion Recovery (FAIR) technique. However, some sensitivity is lost in the BOLD data due to inherent T1 relaxation. We sought to minimize the preceding problem by employing a modified UN-inverted FAIR (UNFAIR) technique, which (in theory) should provide identical CBF signal as FAIR with minimal degradation of the BOLD signal. UNFAIR BOLD maps acquired from human subjects (n = 8) showed significantly higher mean z-score of approximately 17% (p < 0.001), and number of activated voxels at 1.5T. On the other hand, the corresponding FAIR perfusion maps were superior to the UNFAIR perfusion maps as reflected in a higher mean z-score of approximately 8% (p = 0.013), and number of activated voxels. The reduction in UNFAIR sensitivity for perfusion is attributed to increased motion sensitivity related to its higher background signal, and, T2 related losses from the use of an extra inversion pulse. Data acquired at 3.0T demonstrating similar trends are also presented.

A qualitative test of the balloon model for BOLD-based MR signal changes at 3T

The aim of this study was to adapt the balloon model for BOLD-based MR signal changes to a magnetic field strength of 3T and to examine its validity. The simultaneous measurement of BOLD and diffusion-weighted BOLD responses was performed. The amplitude of the BOLD peak was found to be similar for all subjects when a short visual stimulus of 6 sec was used. The rise-time to the BOLD peak and the shape and depth of the poststimulus undershoot varied significantly. A fit of the experimental BOLD responses was found to be possible by use of parameters within a reasonable physiological range. The relations between these parameters and their influence on the modeled BOLD responses is discussed. A prediction of the balloon model is the occurrence of a BOLD overshoot, i.e., a lag between the changes of the blood volume and the blood flow after the start of the stimulation. Experimental evidence for the existence of a BOLD overshoot is presented.

The effects of graded hypercapnia on the activation flow coupling response due to forepaw stimulation in alpha-chloralose anesthetized rats

Activation flow coupling (AFC), changes in cerebral blood flow (CBF) due to changes in neural activity with functional stimulation, provides the physiological basis of many neuroimaging techniques. Hypercapnia leads to an increase in CBF while neural activity remains unaffected. Laser Doppler (LD) flowmetry was used to measure CBF changes (LD(CBF)) in the somatosensory cortex due to periodic electrical forepaw stimulation (4 s in duration) before and during graded hypercapnia (3% CO(2), 5% CO(2) and 10% CO(2)). With increasing CO(2) concentrations, the baseline LD(CBF) progressively increased. The peak height (PH) of the LD(CBF) response, expressed as a percent change from the observed baseline for each hypercapnic state, significantly decreased (P<0.05) with increasing CO(2) concentrations. However, the absolute magnitude of the LD(CBF) change was independent of CO(2) concentration. The temporal dynamics of the LD(CBF) response during hypercapnia were significantly prolonged compared to baseline conditions (P<0.05).

Nonlinear temporal dynamics of the cerebral blood flow response

The linearity of the cerebral perfusion response relative to stimulus duration is an important consideration in the characterization of the relationship between regional cerebral blood flow (CBF), cerebral metabolism, and the blood oxygenation level dependent (BOLD) signal. It is also a critical component in the design and analysis of functional neuroimaging studies. To study the linearity of the CBF response to different duration stimuli, the perfusion response in primary motor and visual cortices was measured during stimulation using an arterial spin labeling technique with magnetic resonance imaging (MRI) that allows simultaneous measurement of CBF and BOLD changes. In each study, the perfusion response was measured for stimuli lasting 2, 6, and 18 sec. The CBF response was found in general to be nonlinearly related to stimulus duration, although the strength of nonlinearity varied between the motor and visual cortices. In contrast, the BOLD response was found to be strongly nonlinear in both regions studied, in agreement with previous findings. The observed nonlinearities are consistent with a model with a nonlinear step from stimulus to neural activity, a linear step from neural activity to CBF change, and a nonlinear step from CBF change to BOLD signal change.

Post-stimulus response in hemodynamics observed by functional magnetic resonance imaging--difference between the primary sensorimotor area and the supplementary motor area

The blood oxygen level dependency (BOLD) contrast is a useful tool for functional neuroimaging based on the hemodynamic response to neuronal activation. We observed different hemodynamic responses in the BOLD signal between the primary sensorimotor area (SM1) and the supplementary motor area (SMA) in the sequential finger movement task. In the SMA, a stronger initial overshoot and a post-stimulus overshoot were observed. It was hypothesized from the time course analysis that the stronger initial overshoot reflected the activation of the SMA for motor control programming in the initial phase. Although the post-stimulus overshoot may be partially explained by cerebral blood flow (CBF) cerebral blood volume (CBV) uncoupling, its mechanism remained unknown. In the SM1, only the initial overshoot was observed and the level of BOLD signal was almost constant after the initial overshoot during the task period. These observations suggested that the BOLD signal is characterized by both CBF-CBV uncoupling and the neuronal activation characteristics in each region.

Comparison of the hemodynamic response to different visual stimuli in single-event and block stimulation fMRI experiments

Experiments with three different types of basic visual stimulation were performed to compare cortical activation in single-event and block trials. Independent of the stimulation paradigm, the single-event presentation leads to highly consistent signal responses regarding both the activated cortical areas and the dynamics of the signal time course. In contrast, signal time courses during block paradigms depend on the stimulus applied and are a complex and nonlinear function of the single-event responses. Additionally, the initial dip during the first 2 seconds after stimulus onset is consistently observed. However, the small amplitude change (-0.1% to -0.3%) requires signal averaging to establish statistical significance of the effect. Furthermore, different patterns of activation were observed within the primary visual cortex. In an anterior part of the primary visual cortex, activation was only observed at the onset and at the cessation of stimulation involving luminance changes.

Investigation of BOLD signal dependence on cerebral blood flow and oxygen consumption: the deoxyhemoglobin dilution model

The relationship between blood oxygenation level-dependent (BOLD) MRI signals, cerebral blood flow (CBF), and oxygen consumption (CMR(O2)) in the physiological steady state was investigated. A quantitative model, based on flow-dependent dilution of metabolically generated deoxyhemoglobin, was validated by measuring BOLD signals and relative CBF simultaneously in the primary visual cortex (V1) of human subjects (N = 12) during graded hypercapnia at different levels of visual stimulation. BOLD and CBF responses to specific conditions were averaged across subjects and plotted as points in the BOLD-CBF plane, tracing out lines of constant CMR(O2). The quantitative deoxyhemoglobin dilution model could be fit to these measured iso-CMR(O2) contours without significant (P </= 0.05) residual error and yielded MRI-based CMR(O2) measurements that were in agreement with PET results for equivalent stimuli. BOLD and CBF data acquired during graded visual stimulation were then substituted into the model with constant parameters varied over plausible ranges. Relative changes in CBF and CMR(O2) appeared to be coupled in an approximate ratio of approximately 2:1 for all realistic parameter settings. Magn Reson Med 42:849-863, 1999.

Linear coupling between cerebral blood flow and oxygen consumption in activated human cortex

The aim of this study was to test the hypothesis that, within a specific cortical unit, fractional changes in cerebral blood flow (CBF) and cerebral metabolic rate of oxygen consumption (CMR(O(2))) are coupled through an invariant relationship during physiological stimulation. This aim was achieved by simultaneously measuring relative changes in these quantities in human primary visual cortex (V1) during graded stimulation with patterns designed to selectively activate different populations of V1 neurons. Primary visual cortex was delineated individually in each subject by using phase-encoded retinotopic mapping. Flow-sensitive alternating inversion recovery MRI, in conjunction with blood oxygenation-sensitive MRI and hypercapnic calibration, was used to monitor CBF and CMR(O(2)). The stimuli used included (i) diffuse isoluminant chromatic displays; (ii) high spatial-frequency achromatic luminance gratings; and (iii) radial checkerboard patterns containing both color and luminance contrast modulated at different temporal rates. Perfusion responses to each pattern were graded by varying luminance and/or color modulation amplitudes. For all stimulus types, fractional changes in blood flow and oxygen uptake were found to be linearly coupled in a consistent ratio of approximately 2:1. The most potent stimulus produced CBF and CMR(O(2)) increases of 48 +/- 5% and 25 +/- 4%, respectively, with no evidence of a plateau for oxygen consumption. Estimation of aerobic ATP yields from the observed CMR(O(2)) increases and comparison with the maximum possible anaerobic ATP contribution indicate that elevated energy demands during brain activation are met largely through oxidative metabolism.