Evidence of a cerebrovascular postarteriole windkessel with delayed compliance

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.

Temporal resolving power of perfusion- and BOLD-based event-related functional MRI

Functional magnetic resonance imaging (fMRI) based on both perfusion and blood oxygenation level-dependent (BOLD) contrasts has been widely applied in spatiotemporal mapping of the human brain function. Temporal resolving power of fMRI is limited by the smoothed hemodynamic response function dispersed from the neuronal activity. In this study, temporal modulation transfer functions were utilized to quantify the resolving powers of perfusion and BOLD fMR signals in time domain. The impulse response function was determined using brief visual stimulations and event-related image acquisition schemes. An important feature of arterial spin labeling techniques is that quantitative perfusion and BOLD signals could be simultaneously acquired. This simultaneous BOLD response may arise from signals that are more proximal to capillary beds, and its temporal resolution may be different from that of the typical BOLD response. Therefore, we assessed and compared the temporal resolving capabilities of perfusion, simultaneous BOLD, and the typical BOLD response obtained from the gradient echo EPI pulse sequence. Full-width-at-half-maximums of perfusion and simultaneous BOLD measurements were significantly smaller than that of BOLD ones (4.3+/- 0.6 s vs 5.5 +/- 0.9 s, p<0.02 and 4.7 +/- 1.3 s vs 5.5 +/- 0.9 s, p<0.01, respectively). The corresponding temporal resolving powers of perfusion and simultaneous BOLD signals were statistically better than that of BOLD signals (0.23 +/- 0.03 Hz vs 0.17 +/- 0.02 Hz, p<0.01 and 0.21 +/- 0.04 Hz vs 0.17 +/- 0.02 Hz, p<0.01, respectively). Our results showed that the typical BOLD response was significantly smoothed from the perfusion response, thus resulting in a degraded temporal resolving power. However, results from the simultaneous BOLD and perfusion measurements were not significantly different. Biophysical implications of the experimental outcomes were further investigated using a computer simulation based on the Balloon model. By fitting the measured data into the model, an apparently longer transit time was obtained for the typical BOLD signal (1.7 s), comparing to that for the simultaneous BOLD one (1.2 s). Therefore, the simultaneous BOLD signal was regarded as less susceptible to the variations from local draining veins. Combining the simulation result with the significantly discrepant resolving powers between the two BOLD signals, we speculated that the blurred effects from large vessels played a predominant role that further reduced the temporal resolution of the BOLD-based fMRI from the perfusion response.

Hemodynamic and metabolic responses to neuronal inhibition

Functional magnetic resonance imaging (fMRI) was used to investigate the changes in blood oxygenation level dependent (BOLD) signal, cerebral blood flow (CBF) and cerebral metabolic rate of oxygen consumption (CMR(O(2))) accompanying neuronal inhibition. Eight healthy volunteers performed a periodic right-hand pinch grip every second using 5% of their maximum voluntary contraction (MVC), a paradigm previously shown to produce robust ipsilateral neuronal inhibition. To simultaneously quantify CBF and BOLD signals, an interleaved multislice pulsed arterial spin labeling (PASL) and T(2)*-weighted gradient echo sequence was employed. The CMR(O(2)) was calculated using the deoxyhemoglobin dilution model, calibrated by data measured during graded hypercapnia. In all subjects, BOLD, CBF and CMR(O(2)) signals increased in the contralateral and decreased in the ipsilateral primary motor (M1) cortex. The relative changes in CMR(O(2)) and CBF were linearly related, with a slope of approximately 0.4. The coupling ratio thus established for both positive and negative CMR(O(2)) and CBF changes is in close agreement with the ones observed by earlier studies investigating M1 perfusion and oxygen consumption increases. These findings characterize the hemodynamic and metabolic downregulation accompanying neuronal inhibition and thereby establish the sustained negative BOLD response as a marker of neuronal deactivation.

Sustained poststimulus elevation in cerebral oxygen utilization after vascular recovery

The brain's response to functional activation is characterized by focal increases in cerebral blood flow. It is generally assumed that this hyperemia is a direct response to the energy demands of activation, the so-called flow-metabolism coupling. Here we report experimental evidence that increases in oxygen metabolism can occur after activation without increases in flow. When using multimodality functional MRI (fMRI) to study visual activation in human brain, we observed a postactivation period of about 30 seconds during which oxygen consumption remained elevated, while blood flow and volume had already returned to baseline levels. The finding of such a prolonged and complete dissociation of vascular response and energy metabolism during the poststimulus period indicates that increased metabolic demand needs not per se cause a concomitant increase in blood flow. The results also show that the postactivation undershoot after the positive blood-oxygen-level-dependent hemodynamic response in fMRI should be reinterpreted as a continued elevation of oxygen metabolism, rather than a delayed blood volume compliance.

Nonlinear coupling of neural activity and CBF in rodent barrel cortex

The relationship between neural activity and accompanying changes in cerebral blood flow (CBF) and oxygenation must be fully understood before data from brain imaging techniques can be correctly interpreted. Whether signals in fMRI reflect the neural input or output of an activated region is still unclear. Similarly, quantitative relationships between neural activity and changes in CBF are not well understood. The present study addresses these issues by using simultaneous laser Doppler flowmetry (LDF) to measure CBF and multichannel electrophysiology to record neural activity in the form of field potentials and multiunit spiking. We demonstrate that CBF-activation coupling is a nonlinear inverse sigmoid function. Comparing the data with previous work suggests that within a cortical model, CBF shows greatest spatial correlation with a current sink 500 microm below the surface corresponding to sensory input. These results show that care must be exercised when interpreting imaging data elicited by particularly strong or weak stimuli and that hemodynamic changes may better reflect the input to a region rather than its spiking output.

Physiological basis and image processing in functional magnetic resonance imaging: neuronal and motor activity in brain

Functional magnetic resonance imaging (fMRI) is recently developing as imaging modality used for mapping hemodynamics of neuronal and motor event related tissue blood oxygen level dependence (BOLD) in terms of brain activation. Image processing is performed by segmentation and registration methods. Segmentation algorithms provide brain surface-based analysis, automated anatomical labeling of cortical fields in magnetic resonance data sets based on oxygen metabolic state. Registration algorithms provide geometric features using two or more imaging modalities to assure clinically useful neuronal and motor information of brain activation. This review article summarizes the physiological basis of fMRI signal, its origin, contrast enhancement, physical factors, anatomical labeling by segmentation, registration approaches with examples of visual and motor activity in brain. Latest developments are reviewed for clinical applications of fMRI along with other different neurophysiological and imaging modalities.

A state-space model of the hemodynamic approach: nonlinear filtering of BOLD signals

In this paper, a new procedure is presented which allows the estimation of the states and parameters of the hemodynamic approach from blood oxygenation level dependent (BOLD) responses. The proposed method constitutes an alternative to the recently proposed Friston [Neuroimage 16 (2002) 513] method and has some advantages over it. The procedure is based on recent groundbreaking time series analysis techniques that have been, in this case, adopted to characterize hemodynamic responses in functional magnetic resonance imaging (fMRI). This work represents a fundamental improvement over existing approaches to system identification using nonlinear hemodynamic models and is important for three reasons. First, our model includes physiological noise. Previous models have been based upon ordinary differential equations that only allow for noise or error to enter at the level of observation. Secondly, by using the innovation method and the local linearization filter, not only the parameters, but also the underlying states of the system generating responses can be estimated. These states can include things like a flow-inducing signal triggered by neuronal activation, de-oxyhemoglobine, cerebral blood flow and volume. Finally, radial basis functions have been introduced as a parametric model to represent arbitrary temporal input sequences in the hemodynamic approach, which could be essential to understanding those brain areas indirectly related to the stimulus. Hence, thirdly, by inferring about the radial basis parameters, we are able to perform a blind deconvolution, which permits both the reconstruction of the dynamics of the most likely hemodynamic states and also, to implicitly reconstruct the underlying synaptic dynamics, induced experimentally, which caused these states variations. From this study, we conclude that in spite of the utility of the standard discrete convolution approach used in statistical parametric maps (SPM), nonlinear BOLD phenomena and unspecific input temporal sequences must be included in the fMRI analysis.

A Novel Method for Determining the Nature of Time Series

The delay vector variance (DVV) method, which analyzes the nature of a time series with respect to the prevalence of deterministic or stochastic components, is introduced. Due to the standardization within the DVV method, it is possible both to statistically test for the presence of nonlinearities in a time series, and to visually inspect the results in a DVV scatter diagram. This approach is convenient for interpretation as it conveys information about the linear or nonlinear nature, as well as about the prevalence of deterministic or stochastic components in the time series, thus unifying the existing approaches which deal either with only deterministic versus stochastic, or the linear versus nonlinear aspect. The results on biomedical time series, namely heart rate variability (HRV) and functional Magnetic Resonance Imaging (fMRI) time series, illustrate the applicability of the proposed DVV-method.

Linear and Nonlinear Relationships between Neuronal Activity, Oxygen Metabolism, and Hemodynamic Responses

We investigated the relationship between neuronal activity, oxygen metabolism, and hemodynamic responses in rat somatosensory cortex with simultaneous optical intrinsic signal imaging and spectroscopy, laser Doppler flowmetry, and local field potential recordings. Changes in cerebral oxygen consumption increased linearly with synaptic activity but with a threshold effect consistent with the existence of a tissue oxygen buffer. Modeling analysis demonstrated that the coupling between neuronal activity and hemodynamic response magnitude may appear linear over a narrow range but incorporates nonlinear effects that are better described by a threshold or power law relationship. These results indicate that caution is required in the interpretation of perfusion-based indicators of brain activity, such as functional magnetic resonance imaging (fMRI), and may help to refine quantitative models of neurovascular coupling.

Temporal dynamics and magnitude of the blood flow response at the optic disk in normal subjects during functional retinal flicker-stimulation

Near-infrared laser Doppler flowmetry was applied in 15 normal volunteers to record the time course and magnitude of changes in the velocity (Vel), volume (Vol) and flow (F) of blood and tissue reflectance (R) at the optic disk in response to 40 and 50 s of increased retinal neural activity. This activity was evoked by diffuse luminance flicker of the retinal posterior pole. After 20 s of flicker, the group averages of Vel, Vol, and F were significantly higher than at baseline (pre-flicker) by 12, 24 and 38%. Time constants of the increases in Vel, Vol, and F were 3.4, 12.7 and 9.1 s, respectively. The group average change in R of 1% was not significant. However, in one subject, 15 recordings from the same site of the optic disk showed a significant increase in R of 8%, with a time course similar to that of Vol. Our findings show that, in the human optic nerve, a white matter tissue, the temporal dynamics and magnitude of the response of blood flow to an increase in retinal neural activity are similar to those reported for brain gray matter. Furthermore, although the R-response could be due, in part, to changes in blood volume, other factors, such as activity-evoked tissue scattering changes, may also affect this response.

A dynamic causal modeling study on category effects: bottom-up or top-down mediation?

In this study, we combined functional magnetic resonance imaging (fMRI) and dynamic causal modeling (DCM) to investigate whether object category effects in the occipital anti temporal cortex are mediated by inputs from early visual cortex or parietal regions. Resolving this issue may provide anatomical constraints on theories of category specificity--which make different assumptions about the underlying neurophysiology. The data were acquired by Ishai, Ungerleider, Martin, Schouten, and Haxby (1999, 2000) and provided by the National fiNRI Data Center (http://www.fmridcc.org). The original authors used a conventional analysis to estimate differential effects in the occipital anti temporal cortex in response to pictures of chairs, faces, and houses. We extended this approach by estimating neuronal interactions that mediate category effects using DCM. DCM uses a Bayesian framework to estimate and make inferences about the influence that one region exerts over another and how this is affected by experimental changes. DCM differs from previous approaches to brain connectivity, such as multivariate autoregressive models and structural equation modeling, as it assumes that the observed hemodynamic responses are driven by experimental changes rather than endogenous noise. DCM therefore brings the analysis of brain connectivity much closer to the analysis of regionally specific effects usually applied to functional imaging data. We used DCM to estimate the influence that V3 and the superior/inferior parietal cortex exerted over category-responsive regions and how this was affected by the presentation of houses, faces, and chairs. We found that category effects in occipital and temporal cortex were mediated by inputs from early visual cortex. In contrast,the connectivity from the superior/inferior parietal area to the category-responsive areas was unaffected by the presentation of chairs, faces, or houses. These findings indicate that category effects in the occipital and temporal cortex can be mediated by bottom-up mechanisms-a finding that needs to be embraced by models of category specificity.

Modeling the hemodynamic response to brain activation

Neural activity in the brain is accompanied by changes in cerebral blood flow (CBF) and blood oxygenation that are detectable with functional magnetic resonance imaging (fMRI) techniques. In this paper, recent mathematical models of this hemodynamic response are reviewed and integrated. Models are described for: (1) the blood oxygenation level dependent (BOLD) signal as a function of changes in cerebral oxygen extraction fraction (E) and cerebral blood volume (CBV); (2) the balloon model, proposed to describe the transient dynamics of CBV and deoxy-hemoglobin (Hb) and how they affect the BOLD signal; (3) neurovascular coupling, relating the responses in CBF and cerebral metabolic rate of oxygen (CMRO(2)) to the neural activity response; and (4) a simple model for the temporal nonlinearity of the neural response itself. These models are integrated into a mathematical framework describing the steps linking a stimulus to the measured BOLD and CBF responses. Experimental results examining transient features of the BOLD response (post-stimulus undershoot and initial dip), nonlinearities of the hemodynamic response, and the role of the physiologic baseline state in altering the BOLD signal are discussed in the context of the proposed models. Quantitative modeling of the hemodynamic response, when combined with experimental data measuring both the BOLD and CBF responses, makes possible a more specific and quantitative assessment of brain physiology than is possible with standard BOLD imaging alone. This approach has the potential to enhance numerous studies of brain function in development, health, and disease.

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.

Differences in the hemodynamic response to event-related motor and visual paradigms as measured by near-infrared spectroscopy

Several current brain imaging techniques rest on the assumption of a tight coupling between neural activity and hemodynamic response. The nature of this neurovascular coupling, however, is not completely understood. There is some evidence for a decoupling of these processes at the onset of neural activity, which manifests itself as a momentary increase in the relative concentration of deoxyhemoglobin (HbR). The existence of this early component of the hemodynamic response function, however, is controversial, as it is inconsistently found. Near infrared spectroscopy (NIRS) allows quantification of levels of oxyhemoglobin (HbO(2)) and HbR during task performance in humans. We acquired NIRS data during performance of simple motor and visual tasks, using rapid-presentation event-related paradigms. Our results demonstrate that rapid, event-related NIRS can provide robust estimates of the hemodynamic response without artifacts due to low-frequency signal components, unlike data from blocked designs. In both the motor and visual data the onset of the increase in HbO(2) occurs before HbR decreases, and there is a poststimulus undershoot. Our results also show that total blood volume (HbT) drops before HbO(2) and undershoots baseline, raising a new issue for neurovascular models. We did not find early deoxygenation in the motor data using physiologically plausible values for the differential pathlength factor, but did find one in the visual data. We suggest that this difference, which is consistent with functional magnetic resonance imaging (fMRI) data, may be attributable to different capillary transit times in these cortices.

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.

Sustained blood oxygenation and volume response to repetition rate-modulated sound in human auditory cortex

The blood oxygen level-dependent (BOLD) signal time course in the auditory cortex is characterized by two components, an initial transient peak and a subsequent sustained plateau with smaller amplitude. Because the T(2)(*) signal detected by functional magnetic resonance imaging (fMRI) depends on at least two counteracting factors, blood oxygenation and volume, we examined whether the reduction in the sustained BOLD signal results from decreased levels of oxygenation or from increased levels of blood volume. We used conventional fMRI to quantify the BOLD signal and fMRI in combination with superparamagnetic contrast agent to quantify blood volume and employed repetition rate-modulated sounds in a silent background to manipulate the response amplitude in the auditory cortex. In the BOLD signal, the initial peak reached 3.3% with pulsed sound and 1.9% with continuous sound, whereas the sustained BOLD signal fell to 2.2% with pulsed sound and to 0.5% with continuous sound, respectively. The repetition rate-dependent reduction in the sustained BOLD amplitude was accompanied by concordant changes in sustained blood volume levels, which, compared to silence, increased by approximately 30% with pulsed and by approximately 10% with continuous sound. Thus, our data suggest that the reduced amplitude of the sustained BOLD signal reflects stimulus-dependent modulation of blood oxygenation rather than blood volume-related effects.

Studying spontaneous EEG activity with fMRI

The multifaceted technological challenge of acquiring simultaneous EEG-correlated fMRI data has now been met and the potential exists for mapping electrophysiological activity with unprecedented spatio-temporal resolution. Work has already begun on studying a host of spontaneous EEG phenomena ranging from alpha rhythm and sleep patterns to epileptiform discharges and seizures, with far reaching clinical implications. However, the transformation of EEG data into linear models suitable for voxel-based statistical hypothesis testing is central to the endeavour. This in turn is predicated upon a number of assumptions regarding the manner in which the generators of EEG phenomena may engender changes in the blood oxygen level dependent (BOLD) signal. Furthermore, important limitations are posed by a set of considerations quite unique to 'paradigmless fMRI'. Here, these issues are assembled and explored to provide an overview of progress made and unresolved questions, with an emphasis on applications in epilepsy.

Can the cerebral metabolic rate of oxygen be estimated with near-infrared spectroscopy?

We have measured the changes in oxy-haemoglobin and deoxy-haemoglobin in the adult human brain during a brief finger tapping exercise using near-infrared spectroscopy (NIRS). The cerebral metabolic rate of oxygen (CMRO2) can be estimated from these NIRS data provided certain model assumptions. The change in CMRO2 is related to changes in the total haemoglobin concentration, deoxy-haemoglobin concentration and blood flow. As NIRS does not provide a measure of dynamic changes in blood flow during brain activation, we relied on a Windkessel model that relates dynamic blood volume and flow changes, which has been used previously for estimating CMRO2 from functional magnetic resonance imaging (fMRI) data. Because of the partial volume effect we are unable to quantify the absolute changes in the local brain haemoglobin concentrations with NIRS and thus are unable to obtain an estimate of the absolute CMRO2 change. An absolute estimate is also confounded by uncertainty in the flow-volume relationship. However, the ratio of the flow change to the CMRO2 change is relatively insensitive to these uncertainties. For the linger tapping task, we estimate a most probable flow-consumption ratio ranging from 1.5 to 3 in agreement with previous findings presented in the literature, although we cannot exclude the possibility that there is no CMRO2 change. The large range in the ratio arises from the large number of model parameters that must be estimated from the data. A more precise estimate of the flow-consumption ratio will require better estimates of the model parameters or flow information, as can be provided by combining NIRS with fMRI.

Functional magnetic resonance imaging based on changes in vascular space occupancy

During brain activation, local control of oxygen delivery is facilitated through microvascular dilatation and constriction. A new functional MRI (fMRI) methodology is reported that is sensitive to these microvascular adjustments. This contrast is accomplished by eliminating the blood signal in a manner that is independent of blood oxygenation and flow. As a consequence, changes in cerebral blood volume (CBV) can be assessed through changes in the remaining extravascular water signal (i.e., that of parenchymal tissue) without need for exogenous contrast agents or any other invasive procedures. The feasibility of this vascular space occupancy (VASO)-dependent functional MRI (fMRI) approach is demonstrated for visual stimulation, breath-hold (hypercapnia), and hyperventilation (hypocapnia). During visual stimulation and breath-hold, the VASO signal shows an inverse correlation with the stimulus paradigm, consistent with local vasodilatation. This effect is reversed during hyperventilation. Comparison of the hemodynamic responses of VASO-fMRI, cerebral blood flow (CBF)-based fMRI, and blood oxygenation level-dependent (BOLD) fMRI indicates both arteriolar and venular temporal characteristics in VASO. The effect of changes in water exchange rate and partial volume contamination with CSF were calculated to be negligible. At the commonly-used fMRI resolution of 3.75 x 3.75 x 5 mm(3), the contrast-to-noise-ratio (CNR) of VASO-fMRI was comparable to that of CBF-based fMRI, but a factor of 3 lower than for BOLD-fMRI. Arguments supporting a better gray matter localization for the VASO-fMRI approach compared to BOLD are provided.

Neuronal activity-induced changes of local cerebral microvascular blood oxygenation in the rat: effect of systemic hyperoxia or hypoxia

To investigate the effect of resting blood oxygen concentration on the hemodynamic response to functional brain activation, we compared activation-induced changes in hemoglobin oxygenation during normoxia with systemic hyperoxia or mild hypoxia. Hemoglobin oxygenation changes were measured by microfiber optical spectroscopy (500-590 nm) in response to physiological whisker barrel cortex activation by whole whisker pad deflection (4 s, 4 Hz) in alpha-chloralose/urethane anesthetised male Wistar rats. During systemic hyperoxia (n=6), the stimulation-induced hyperoxygenation response was decreased and prolonged, whereas during mild systemic hypoxia (n=7) the peak response was significantly increased followed by a faster return to baseline. During normoxia, a poststimulation under- (oxy-hemoglobin) and overshoot (deoxy-hemoglobin) was observed, which disappeared during systemic hyperoxia and was pronounced during systemic hypoxia. Although averaging out below statistical significance when combining all animals, during mild systemic hypoxia a very small early increase of deoxy-Hb at the beginning of the stimulation period was conspicuous more often than during normoxic or even hyperoxic conditions. This small early increase of deoxy-Hb never preceded the onset of the oxy-Hb response, and was not accompanied by a concomitant decrease in oxy-Hb. Hyperoxia or hypoxia did not affect the induced neuronal responses. Our findings support the concept that the hemodynamic response is regulated according to the metabolic demand of oxygen within the activated tissue.

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.

Relation between cerebral blood flow and metabolism explained by a model of oxygen exchange

The cerebral blood flow (CBF) and cerebral metabolic rate of oxygen (CMRo(2)) are major determinants of the contrast in functional magnetic resonance imaging and optical imaging. However, the coupling between CBF and CMRo(2) during cerebral activation remains controversial. Whereas most of the previous models tend to show a nonlinear coupling, experimental studies have led to conflicting conclusions. A physiologic model was developed of oxygen transport through the blood-brain barrier (BBB) for dynamic and stationary states. Common model simplifications are proposed and their implications for the CBF/CMRo(2) relation are studied. The tissue oxygen pool, the BBB permeability, and the hemoglobin dissociation curve are physiologic parameters directly involved in the CBF/CMRo(2) relation. We have been shown that the hypothesis of a negligible tissue oxygen pool, which was admitted by most of the previous models, implies a tight coupling between CBF and CMRo(2). By relaxing this hypothesis, a real uncoupling was allowed that gives a more coherent view of the CBF/CMRo(2) relation, in better agreement with the hypercapnia data and with the variability reported in experimental works for the relative changes of those two variables. This also allows a temporal mismatch between CBF and CMRo(2), which influences the temporal shape of oxygenation at the capillary end.

Imaging methods for evaluating brain function in man

Imaging of brain function and neurotransmission is an important bridge between basic and clinical research. Regional cerebral energy metabolism and blood flow are normally coupled to regional cerebral function. Positron tomography (PET) studies of cerebral glucose metabolism and blood flow, single photon tomography (SPECT) and MRI studies of cerebral perfusion, have been used to image cerebral development and aging in man. The sensitivity, temporal resolution, spatial resolution and lack of radiation have led to the widespread utilization of blood oxygen level dependent (BOLD) and MRI perfusion techniques. PET and SPECT methods for studying cerebral neurotransmission include studies of dopaminergic, serotonergic, cholinergic, opiate and GABAergic neurotransmission in man. Studies of cerebral neurotransmission in man have helped to delineate the mechanisms of action of antipsychotic and antidepressant drugs, the diagnosis and progression of Parkinson's disease, and to evaluate neuroprotective drugs. The strengths, limitations, and application of these modalities are reviewed. The application of these methods to cerebral development and aging are briefly discussed.

The hemodynamic impulse response to a single neural event

This article investigates the relation between stimulus-evoked neural activity and cerebral hemodynamics. Specifically, the hypothesis is tested that hemodynamic responses can be modeled as a linear convolution of experimentally obtained measures of neural activity with a suitable hemodynamic impulse response function. To obtain a range of neural and hemodynamic responses, rat whisker pad was stimulated using brief (</=2 seconds) electrical stimuli consisting of single pulses (0.3 millisecond, 1.2 mA) combined both at different frequencies and in a paired-pulse design. Hemodynamic responses were measured using concurrent optical imaging spectroscopy and laser Doppler flowmetry, whereas neural responses were assessed through current source density analysis of multielectrode recordings from a single barrel. General linear modeling was used to deconvolve the hemodynamic impulse response to a single "neural event" from the hemodynamic and neural responses to stimulation. The model provided an excellent fit to the empirical data. The implications of these results for modeling schemes and for physiologic systems coupling neural and hemodynamic activity are discussed.

Simultaneous imaging of total cerebral hemoglobin concentration, oxygenation, and blood flow during functional activation

A simple instrument is demonstrated for high-resolution simultaneous imaging of total hemoglobin concentration and oxygenation and blood flow in the brain by combining rapid multiwavelength imaging with laser speckle contrast imaging. The instrument was used to image changes in oxyhemoglobin and deoxyhemoglobin and blood flow during cortical spreading depression and single whisker stimulation in rats through a thinned skull. The ability to image blood flow and hemoglobin concentration changes simultaneously with high resolution will permit detailed quantitative analysis of the spatiotemporal hemodynamics of functional brain activation, including imaging of oxygen metabolism. This is of significance to the neuroscience community and will lead to a better understanding of the interrelationship of neural, metabolic, and hemodynamic processes in normal and diseased brains.

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.

Cerebral venous and arterial blood volumes can be estimated separately in humans using magnetic resonance imaging

Approaches to obtain quantitative, noninvasive estimates of total cerebral blood volume (tCBV) and cerebral venous blood volume (vCBV) separately in humans are proposed. Two sequences were utilized, including a 3D high-resolution gradient-echo (GE) sequence and a 2D multi-echo GE/spin-echo (MEGESE) sequence. Images acquired by the former sequence provided an estimate of background magnetic field variations (DeltaB), while images obtained by the latter sequence were utilized to obtain separate measures of tCBV and vCBV with and without contrast agent. Prior to the calculation of vCBV and tCBV, the acquired images were corrected for signal loss induced by the presence of DeltaB. vCBV and tCBV were estimated to be 2.46% +/- 0.28% and 3.20% +/- 0.41%, respectively, after the DeltaB correction, which in turn provided a vCBV/tCBV ratio of 0.77 +/- 0.04, in excellent agreement with results reported in the literature. Our results demonstrate that quantitative estimates of vCBV and tCBV can be obtained in vivo.

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.

Effect of basal conditions on the magnitude and dynamics of the blood oxygenation level-dependent fMRI response

The effect of the basal cerebral blood flow (CBF) on both the magnitude and dynamics of the functional hemodynamic response in humans has not been fully investigated. Thus, the hemodynamic response to visual stimulation was measured using blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging (fMRI) in human subjects in a 7-T magnetic field under different basal conditions: hypocapnia, normocapnia, and hypercapnia. Hypercapnia was induced by inhalation of a 5% carbon dioxide gas mixture and hypocapnia was produced by hyperventilation. As the fMRI baseline signal increased linearly with expired CO2 from hypocapnic to hypercapnic levels, the magnitude of the BOLD response to visual stimulation decreased linearly. Measures of the dynamics of the visually evoked BOLD response (onset time, full-width-at-half-maximum, and time-to-peak) increased linearly with the basal fMRI signal and the end-tidal CO2 level. The basal CBF level, modulated by the arterial partial pressure of CO2, significantly affects both the magnitude and dynamics of the BOLD response induced by neural activity. These results suggest that caution should be exercised when comparing stimulus-induced fMRI responses under different physiologic or pharmacologic states.

The neural basis of the blood-oxygen-level-dependent functional magnetic resonance imaging signal

Magnetic resonance imaging (MRI) has rapidly become an important tool in clinical medicine and biological research. Its functional variant (functional magnetic resonance imaging; fMRI) is currently the most widely used method for brain mapping and studying the neural basis of human cognition. While the method is widespread, there is insufficient knowledge of the physiological basis of the fMRI signal to interpret the data confidently with respect to neural activity. This paper reviews the basic principles of MRI and fMRI, and subsequently discusses in some detail the relationship between the blood-oxygen-level-dependent (BOLD) fMRI signal and the neural activity elicited during sensory stimulation. To examine this relationship, we conducted the first simultaneous intracortical recordings of neural signals and BOLD responses. Depending on the temporal characteristics of the stimulus, a moderate to strong correlation was found between the neural activity measured with microelectrodes and the BOLD signal averaged over a small area around the microelectrode tips. However, the BOLD signal had significantly higher variability than the neural activity, indicating that human fMRI combined with traditional statistical methods underestimates the reliability of the neuronal activity. To understand the relative contribution of several types of neuronal signals to the haemodynamic response, we compared local field potentials (LFPs), single- and multi-unit activity (MUA) with high spatio-temporal fMRI responses recorded simultaneously in monkey visual cortex. At recording sites characterized by transient responses, only the LFP signal was significantly correlated with the haemodynamic response. Furthermore, the LFPs had the largest magnitude signal and linear systems analysis showed that the LFPs were better than the MUAs at predicting the fMRI responses. These findings, together with an analysis of the neural signals, indicate that the BOLD signal primarily measures the input and processing of neuronal information within a region and not the output signal transmitted to other brain regions.

A model of the hemodynamic response and oxygen delivery to brain

A recent nonlinear system by Friston et al. (2000. NeuroImage 12: 466-477) links the changes in BOLD response to changes in neural activity. The system consists of five subsystems, linking: (1) neural activity to flow changes; (2) flow changes to oxygen delivery to tissue; (3) flow changes to changes in blood volume and venous outflow; (4) changes in flow, volume, and oxygen extraction fraction to deoxyhemoglobin changes; and finally (5) volume and deoxyhemoglobin changes to the BOLD response. Friston et al. exploit, in subsystem 2, a model by Buxton and Frank coupling flow changes to changes in oxygen metabolism which assumes tissue oxygen concentration to be close to zero. We describe below a model of the coupling between flow and oxygen delivery which takes into account the modulatory effect of changes in tissue oxygen concentration. The major development has been to extend the original Buxton and Frank model for oxygen transport to a full dynamic capillary model making the model applicable to both transient and steady state conditions. Furthermore our modification enables us to determine the time series of CMRO(2) changes under different conditions, including CO(2) challenges. We compare the differences in the performance of the "Friston system" using the original model of Buxton and Frank and that of our model. We also compare the data predicted by our model (with appropriate parameters) to data from a series of OIS studies. The qualitative differences in the behaviour of the models are exposed by different experimental simulations and by comparison with the results of OIS data from brief and extended stimulation protocols and from experiments using hypercapnia.

Bayesian estimation of dynamical systems: an application to fMRI

This paper presents a method for estimating the conditional or posterior distribution of the parameters of deterministic dynamical systems. The procedure conforms to an EM implementation of a Gauss-Newton search for the maximum of the conditional or posterior density. The inclusion of priors in the estimation procedure ensures robust and rapid convergence and the resulting conditional densities enable Bayesian inference about the model parameters. The method is demonstrated using an input-state-output model of the hemodynamic coupling between experimentally designed causes or factors in fMRI studies and the ensuing BOLD response. This example represents a generalization of current fMRI analysis models that accommodates nonlinearities and in which the parameters have an explicit physical interpretation. Second, the approach extends classical inference, based on the likelihood of the data given a null hypothesis about the parameters, to more plausible inferences about the parameters of the model given the data. This inference provides for confidence intervals based on the conditional density.

Hemodynamic response in the unanesthetized rat: intrinsic optical imaging and spectroscopy of the barrel cortex

Optical imaging spectroscopy was used to measure the hemodynamic response of somatosensory cortex to stimulation of the whiskers. Responses to brief puffs of air were compared in anesthetized and unanesthetized rats. The hemodynamic response was approximately four times larger in the unanesthetized animal than the corresponding anesthetized animal. In unanesthetized animals, a short-latency (approximately 400 milliseconds) short-duration (approximately 300 milliseconds) hemodynamic startle response was observed. General linear model analysis was used to extract this component from the time series, and revealed an underlying short-latency increase in deoxygenated hemoglobin in response to somatosensory stimulation. It is proposed that anesthesia can have a marked affect on the relation between changes in blood volume and blood flow. This work represents a step in the development of an experimental model that can be used to investigate fundamental neurologic processes in the awake-behaving rodent.

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).

Paradoxical correlation between signal in functional magnetic resonance imaging and deoxygenated haemoglobin content in capillaries: a new theoretical explanation

Signal increases in functional magnetic resonance imaging (fMRI) are believed to be a result of decreased paramagnetic deoxygenated haemoglobin (deoxyHb) content in the neural activation area. However, discrepancies in this canonical blood oxygenation level dependent (BOLD) theory have been pointed out in studies using optical techniques, which directly measure haemoglobin changes. To explain the discrepancies, we developed a new theory bridging magnetic resonance (MR) signal and haemoglobin changes. We focused on capillary influences, which have been neglected in most previous fMRI studies and performed a combined fMRI and near-infrared spectroscopy (NIRS) study using a language task. Paradoxically, both the MR signal and deoxyHb content increased in Broca's area. On the other hand, fMRI activation in the auditory area near large veins correlated with a mirror-image decrease in deoxyHb and increase in oxygenated haemoglobin (oxyHb), in agreement with canonical BOLD theory. All fMRI signal changes correlated consistently with changes in oxyHb, the diamagnetism of which is insensitive to MR. We concluded that the discrepancy with the canonical BOLD theory is caused by the fact that the BOLD theory ignores the effect of the capillaries. Our theory explains the paradoxical phenomena of the oxyHb and deoxyHb contributions to the MR signal and gives a new insight into the precise haemodynamics of activation by analysing fMRI and NIRS data.

Changes in blood flow, oxygenation, and volume following extended stimulation of rodent barrel cortex

Simultaneous optical imaging spectroscopy and laser-Doppler flowmetry were used in rodent barrel cortex to examine the hemodynamic response to extended electrical stimulation (20 s, 5 Hz) of the whisker pad. Stimulation results in a fast early increase in deoxyhemoglobin concentration (Hbr) followed by a later decrease to a "plateau" phase approximately 4 s after stimulation onset. There was no corresponding decrease in oxyhemoglobin (HbO(2)), which simply increased after stimulation, reaching a plateau at approximately 5 s. The time series of flow and volume had similar onset times and did not differ significantly throughout the presentation of the stimuli. Following stimulation cessation all aspects of the hemodynamic response returned to baseline with a long decay constant (>20 s), CBV doing so at a slower rate than CBF. The time courses of CBF, CBV, Hbr, and HbO(2) were very similar to that produced by a brief stimulus up to peak. The relationship between the flow and the volume changes is well approximated by the expression CBV = CBF(phi). We find phi to be slightly lower under stimulation (0.26 +/- 0.0152) than during hypercapnia (0.32 +/- 0.0172). Saturation and flow data were used to estimate changes in CMRO(2) for a range of baseline oxygen extraction fractions (OEF). In the case of hypercapnia CMRO(2) was biphasic, increasing after onset and sharply decreasing below baseline following cessation. If it is assumed that there is no "net" increase in CMRO(2) (i.e., SigmaDeltaCMRO(2) = 0) following the onset and offset of hypercapnia, then the corresponding estimate of baseline OEF is 0.45. Evidence for increased oxygen consumption was obtained for all stimulation intensities assuming a baseline OEF of 0.45.

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.

Relationship of spikes, synaptic activity, and local changes of cerebral blood flow

The coupling of electrical activity in the brain to changes in cerebral blood flow (CBF) is of interest because hemodynamic changes are used to track brain function. Recent studies, especially those investigating the cerebellar cortex, have shown that the spike rate in the principal target cell of a brain region (i.e. the efferent cell) does not affect vascular response amplitude. Subthreshold integrative synaptic processes trigger changes in the local microcirculation and local glucose consumption. The spatial specificity of the vascular response on the brain surface is limited because of the functional anatomy of the pial vessels. Within the cortex there is a characteristic laminar flow distribution, the largest changes of which are observed at the depth of maximal synaptic activity (i.e. layer IV) for an afferent input system. Under most conditions, increases in CBF are explained by activity in postsynaptic neurons, but presynaptic elements can contribute. Neurotransmitters do not mediate increases in CBF that are triggered by the concerted action of several second messenger molecules. It is important to distinguish between effective synaptic inhibition and deactivation that increase and decrease CBF and glucose consumption, respectively. In summary, hemodynamic changes evoked by neuronal activity depend on the afferent input function (i.e. all aspects of presynaptic and postsynaptic processing), but are totally independent of the efferent function (i.e., the spike rate of the same region). Thus, it is not possible to conclude whether the output level of activity of a region is increased based on brain maps that use blood-flow changes as markers.

Visual motion processing investigated using contrast agent-enhanced fMRI in awake behaving monkeys

To reduce the information gap between human neuroimaging and macaque physiology and anatomy, we mapped fMRI signals produced by moving and stationary stimuli (random dots or lines) in fixating monkeys. Functional sensitivity was increased by a factor of approximately 5 relative to the BOLD technique by injecting a contrast agent (monocrystalline iron oxide nanoparticle [MION]). Areas identified as motion sensitive included V2, V3, MT/V5, vMST, FST, VIP, and FEF (with moving dots), as well as V4, TE, LIP, and PIP (with random lines). These regions sensitive for moving dots are largely in agreement with monkey single unit data and (except for V3A) with human fMRI results. Moving lines activate some regions that have not been previously implicated in motion processing. Overall, the results clarify the relationship between the motion pathway and the dorsal stream in primates.

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.

Nonlinear coupling between evoked rCBF and BOLD signals: a simulation study of hemodynamic responses

The aim of this work was to investigate the dependence of BOLD responses on different patterns of stimulus input/neuronal changes. In an earlier report, we described an input-state-output model that combined (i) the Balloon/Windkessel model of nonlinear coupling between rCBF and BOLD signals, and (ii) a linear model of how regional flow changes with synaptic activity. In the present investigation, the input-state-output model was used to explore the dependence of simulated PET (rCBF) and fMRI (BOLD) signals on various parameters pertaining to experimental design. Biophysical simulations were used to estimate rCBF and BOLD responses as functions of (a) a prior stimulus, (b) epoch length (for a fixed SOA), (c) SOA (for a fixed number of events), and (d) stimulus amplitude. We also addressed the notion that a single neuronal response may differ, in terms of the relative contributions of early and late neural components, and investigated the effect of (e) the relative size of the late or "endogenous" neural component. We were interested in the estimated average rCBF and BOLD responses per stimulus or event, not in the statistical efficiency with which these responses are detected. The BOLD response was underestimated relative to rCBF with a preceding stimulus, increasing epoch length, and increasing SOA. Furthermore, the BOLD response showed some highly nonlinear behaviour when varying stimulus amplitude, suggesting some form of hemodynamic "rectification." Finally, the BOLD response was underestimated in the context of large late neuronal components. The difference between rCBF and BOLD is attributed to the nonlinear transduction of rCBF to BOLD signal. Our simulations support the idea that varying parameters that specify the experimental design may have differential effects in PET and fMRI. Moreover, they show that fMRI can be asymmetric in its ability to detect deactivations relative to activations when an absolute baseline is stipulated. Finally, our simulations suggest that relative insensitivity to BOLD signal in specific regions, such as the temporal lobe, may be partly explained by higher cognitive functions eliciting a relatively large late endogenous neuronal component.

Spatial heterogeneity of the nonlinear dynamics in the FMRI BOLD response

Recent studies of blood oxygenation level dependent (BOLD) signal responses averaged over a region of interest have demonstrated that the response is nonlinear with respect to stimulus duration. Specifically, shorter duration stimuli produce signal changes larger than expected from a linear system. The focus of this study is to characterize the spatial heterogeneity of this nonlinear effect. A series of MR images of the visual and motor cortexes were acquired during visual stimulation and finger tapping, respectively, at five different stimulus durations (SD). The nonlinearity was assessed by fitting ideal linear responses to the responses at each SD. This amplitude, which is constant for different SD in a linear system, was normalized by the amplitude of the response to a blocked design, thus describing the amount by which the stimulus is larger than predicted from a linear extrapolation of the response to the long duration stimulus. The amplitude of the BOLD response showed a nonlinear behavior that varied considerably and consistently over space, ranging from almost linear to 10 times larger than a linear prediction at short SD. In the motor cortex different nonlinear behavior was found in the primary and supplementary motor cortexes.

Coupling of neural activity and BOLD fMRI response: new insights by combination of fMRI and VEP experiments in transition from single events to continuous stimulation

Functional magnetic resonance imaging (fMRI) measures the correlation between the fMRI response and stimulus properties. A linear relationship between neural activity and fMRI response is commonly assumed. However, the response to repetitive stimulation cannot be explained by a simple superposition of single-event responses. This might be due to neural adaptation or the hemodynamic changes underlying the fMRI BOLD response. To assess the influence of adaptation, the BOLD responses and visual evoked potentials (VEPs) to identical stimuli were recorded. To achieve different adaptation levels, 2-s stimulus epochs alternated with different interstimulus intervals (ISI = 0.0, 0.4, 0.8, 2.0, and 12 s) were presented. Neural adaptation during the checkerboard reversal paradigm used for fMRI measurements is demonstrated. Even if the measured VEP amplitude is used as the weighting function for a linear model, the measured BOLD fMRI signal time-course is not adequately predicted.

In vivo cerebrovascular measurement combining diffuse near-infrared absorption and correlation spectroscopies

We combine two near-infrared diffuse optical techniques to study variations of blood flow, haemoglobin concentration, and blood oxygen saturation in the functioning rat brain. Diffuse correlation spectroscopy (or flowmetry) monitors changes in the cerebral blood flow, without the use of the principles of tracer clearance, by measuring the optical phase-shifts caused by moving blood cells. Near-infrared absorption spectroscopy concurrently measures tissue absorption at two wavelengths to determine haemoglobin concentration and blood oxygen saturation in this same tissue volume. This optical probe is non-invasive and was employed through the intact skull. The utility of the technique is demonstrated in vivo by measuring the temporal changes in the regional vascular dynamics of rat brain during hypercapnia. Temporal and spatial variations of cerebral blood flow, haemoglobin concentration and blood oxygen saturation during hypercapnia are compared with other measurements in the literature, and a quantitative analysis demonstrating the self-consistency of our combined observations of vascular response is presented.

Increased Oxygen Consumption Following Activation of Brain: Theoretical Footnotes Using Spectroscopic Data from Barrel Cortex

Optical imaging spectroscopy (OIS) and laser Doppler flowmetry (LDF) data sequences from anesthetized rats were used to determine the relationship between changes in oxy-and deoxygenated hemoglobin concentration and changes in blood volume and flow in the presence and absence of stimulation. The data from Jones et al. (accompanying paper) were used to explore the differences between two theoretical models of flow activation coupling. The essential difference between the two models is the extension of the model of Buxton and Frank by Hyder et al. (1998, J. Appl. Physiol. 85: 554--564) to incorporate change in capillary diffusivity coupled to flow. In both models activation-increased flow changes increase oxygen transport from the capillary; however, in Hyder et al.'s model the diffusivity of the capillary itself is increased. Hyder et al. proposed a parameter (Omega), a scaling "constant" linking increased blood flow and oxygen "diffusivity" in the capillary bed. Thus, in Buxton and Frank's theory, Omega = 0; i.e., there are no changes in diffusivity. In Hyder et al.'s theory, 0 < Omega < 1, and changes in diffusivity are assumed to be linearly related to flow changes. We elaborate the theoretical position of both models to show that, in principle, the different predictions from the two theories can be evaluated using optical imaging spectroscopy data. We find that both theoretical positions have limitations when applied to data from brief stimulation and when applied to data from mild hypercapnia. In summary, the analysis showed that although Hyder et al.'s proposal that diffusivity increased during activation did occur; it was shown to arise from an implementation of Buxton and Frank's theory under episodes of brief stimulation. The results also showed that the scaling parameter Omega is not a constant as the Hyder et al. model entails but in fact varies over the time course of the flow changes. Data from experiments in which mild hypercapnia was administered also indicated changes in the diffusivity of the capillary bed, but in this case the changes were negative; i.e., oxygen transport from the capillary decreased relative to baseline under hypercapnia. Neither of the models could account for the differences between the hypercapnia and activation data when matched for equivalent flow changes. A modification to the models to allow non-null tissue oxygen concentrations that can be moderated by changes due to increased metabolic demand following increased neural activity is proposed. This modification would allow modulation of oxygen transport from the capillary bed (e.g., changes in diffusivity) by tissue oxygen tension and would allow a degree of decoupling of flow and oxygen delivery, which can encompass both the data from stimulation and from hypercapnia.

Concurrent Optical Imaging Spectroscopy and Laser-Doppler Flowmetry: The Relationship between Blood Flow, Oxygenation, and Volume in Rodent Barrel Cortex

Functional magnetic resonance imaging (fMRI) is based on the coupling between neural activity and changes in the concentration of the endogenous paramagnetic contrast agent deoxygenated hemoglobin. Changes in the blood oxygen level-dependent (BOLD) signal result from a complex interplay of blood volume, flow, and oxygen consumption. Optical imaging spectroscopy (OIS) has been used to measure changes in blood volume and saturation in response to increased neural activity, while laser Doppler Flowmetry (LDF) can be used to measure flow changes and is now commonplace in neurovascular research. Here, we use concurrent OIS and LDF to examine the hemodynamic response in rodent barrel cortex using electrical stimulation of the whisker pad at varying intensities. Spectroscopic analysis showed that stimulation produced a biphasic early increase in deoxygenated hemoglobin (Hbr), followed by a decrease below baseline, reaching minima at approximately 3.7 s. There was no evidence for a corresponding early decrease in oxygenated hemoglobin (HbO(2)), which simply increased after stimulation, reaching maximum at approximately 3.2 s. The time courses of changes in blood volume (CBV) and blood flow (CBF) were similar. Both increased within a second of stimulation onset and peaked at approximately 2.7 s, after which CBV returned to baseline at a slower rate than CBF. The changes in Hbr, Hbt, and CBF were used to estimate changes in oxygen consumption (CMRO(2)), which increased within a second of stimulation and peaked approximately 2.2 s after stimulus onset. Analysis of the relative magnitudes of CBV and CBF indicates that the fractional changes of CBV could be simply scaled to match those of CBF. We found the relationship to be well approximated by CBV = CBF(0.29). A similar relationship was found using the response to elevated fraction of inspired carbon dioxide (FICO(2)).