Stimulus-dependent BOLD and perfusion dynamics in human V1

Physiological modeling of the BOLD signal and implications for effective connectivity: A primer

In this primer, I provide an overview of the physiological processes that contribute to the observed BOLD signal (i.e., the generative biophysical model), including their time course properties within the framework of the physiologically-informed dynamic causal modeling (P-DCM). The BOLD signal is primarily determined by the change in paramagnetic deoxygenated hemoglobin, which results from combination of changes in oxygen metabolism, and cerebral blood flow and volume. Specifically, the physiological origin of the so-called BOLD signal "transients" will be discussed, including the initial overshoot, steady-state activation and the post-stimulus undershoot. I argue that incorrect physiological assumptions in the generative model of the BOLD signal can lead to incorrect inferences pertaining to both local neuronal activity and effective connectivity between brain regions. In addition, I introduce the recent laminar BOLD signal model, which extends P-DCM to cortical depths-resolved BOLD signals, allowing for laminar neuronal activity to be determined using high-resolution fMRI data.

Functional MRS studies of GABA and glutamate/Glx - A systematic review and meta-analysis

Functional magnetic resonance spectroscopy (fMRS) can be used to investigate neurometabolic responses to external stimuli in-vivo, but findings are inconsistent. We performed a systematic review and meta-analysis on fMRS studies of the primary neurotransmitters Glutamate (Glu), Glx (Glutamate + Glutamine), and GABA. Data were extracted, grouped by metabolite, stimulus domain, and brain region, and analysed by determining standardized effect sizes. The quality of individual studies was rated. When results were analysed by metabolite type small to moderate effect sizes of 0.29-0.47 (p < 0.05) were observed for changes in Glu and Glx regardless of stimulus domain and brain region, but no significant effects were observed for GABA. Further analysis suggests that Glu, Glx and GABA responses differ by stimulus domain or task and vary depending on the time course of stimulation and data acquisition. Here, we establish effect sizes and directionality of GABA, Glu and Glx response in fMRS. This work highlights the importance of standardised reporting and minimal best practice for fMRS research.

Separating vascular and neuronal effects of age on fMRI BOLD signals

Accurate identification of brain function is necessary to understand the neurobiology of cognitive ageing, and thereby promote well-being across the lifespan. A common tool used to investigate neurocognitive ageing is functional magnetic resonance imaging (fMRI). However, although fMRI data are often interpreted in terms of neuronal activity, the blood oxygenation level-dependent (BOLD) signal measured by fMRI includes contributions of both vascular and neuronal factors, which change differentially with age. While some studies investigate vascular ageing factors, the results of these studies are not well known within the field of neurocognitive ageing and therefore vascular confounds in neurocognitive fMRI studies are common. Despite over 10 000 BOLD-fMRI papers on ageing, fewer than 20 have applied techniques to correct for vascular effects. However, neurovascular ageing is not only a confound in fMRI, but an important feature in its own right, to be assessed alongside measures of neuronal ageing. We review current approaches to dissociate neuronal and vascular components of BOLD-fMRI of regional activity and functional connectivity. We highlight emerging evidence that vascular mechanisms in the brain do not simply control blood flow to support the metabolic needs of neurons, but form complex neurovascular interactions that influence neuronal function in health and disease. This article is part of the theme issue 'Key relationships between non-invasive functional neuroimaging and the underlying neuronal activity'.

Influence of contrast-reversing frequency on the amplitude and spatial distribution of visual cortex hemodynamic responses

Visual stimulation is one of the most commonly used paradigms for cerebral cortex function investigation. Experiments typically involve presenting to a volunteer a black-and-white checkerboard with contrast-reversing at a frequency of 4 to 16 Hz. The aim of the present study was to investigate the influence of the flickering frequency on the amplitude of changes in the concentration of oxygenated and deoxygenated hemoglobin. The hemoglobin concentrations were assessed with the use of a high resolution diffuse optical tomography method. Spatial distributions of changes in hemoglobin concentrations overlaying the visual cortex are shown for various stimuli frequencies. Moreover, the hemoglobin concentration changes obtained for different source-detector separations (from 1.5 to 5.4 cm) are presented. Our results demonstrate that the flickering frequency had a statistically significant effect on the induced oxyhemoglobin changes (p < 0,001). The amplitude of oxy hemoglobin concentration changes at a frequency of 8 Hz was higher in comparison with that measured at 4 Hz :[median(25th-75thpercentiles) 1.24 (0.94-1.71) vs. 0.92(0.73-1.28)µM, p < 0.001]; 12 Hz:[1.24 (0.94-1.71) vs. 1.04 (0.78-1.32) µM, p < 0.001]; and 16 Hz:[1.24 (0.94-1.71) vs. 1.15(0.87-1.48) µM, p < 0.001]. No significant differences were observed between the size of an area of activation for various frequencies. The demonstrated superiority of 8 Hz over other frequencies can advance understanding of visual stimulations and help guide future fNIRS protocols.

Gradient-echo and spin-echo blood oxygenation level-dependent functional MRI at ultrahigh fields of 9.4 and 15.2 Tesla

Purpose

Sensitivity and specificity of blood oxygenation level–dependent (BOLD) functional MRI (fMRI) is sensitive to magnetic field strength and acquisition methods. We have investigated gradient‐echo (GE)‐ and spin‐echo (SE)‐BOLD fMRI at ultrahigh fields of 9.4 and 15.2  Tesla.

Methods

BOLD fMRI experiments responding to forepaw stimulation were performed with 3 echo times (TE) at each echo type and B₀ in α‐chloralose–anesthetized rats. The contralateral forelimb somatosensory region was selected for quantitative analyses.

Results

At 9.4 T and 15.2 T, average baseline T₂^* (n = 9) was 26.6 and 17.1 msec, whereas baseline T₂ value (n = 9) was 35.7 and 24.5 msec, respectively. Averaged stimulation‐induced ΔR₂^* was –1.72 s^–1 at 9.4 T and –3.09 s^–1 at 15.2 T, whereas ΔR₂ was –1.19 s^–1 at 9.4 T and –1.97 s^–1 at 15.2 T. At the optimal TE of tissue T₂^* or T₂, BOLD percent changes were slightly higher at 15.2 T than at 9.4 T (GE: 7.4% versus 6.4% and SE: 5.7% versus 5.4%). The ΔR₂^* and ΔR₂ ratio of 15.2 T to 9.4 T was 1.8 and 1.66, respectively. The ratio of the macrovessel‐containing superficial to microvessel‐dominant parenchymal BOLD signal was 1.73 to 1.76 for GE‐BOLD versus 1.13 to 1.19 for SE‐BOLD, indicating that the SE‐BOLD contrast is less sensitive to macrovessels than GE‐BOLD.

Conclusion

SE‐BOLD fMRI improves spatial specificity to microvessels compared to GE‐BOLD at both fields. BOLD sensitivity is similar at the both fields and can be improved at ultrahigh fields only for thermal‐noise–dominant ultrahigh‐resolution fMRI.

Response Hand and Motor Set Differentially Modulate the Connectivity of Brain Pathways During Simple Uni-manual Motor Behavior

We investigated the flexible modulation of undirected functional connectivity (uFC) of brain pathways during simple uni-manual responding. Two questions were central to our interests: (1) does response hand (dominant vs. non-dominant) differentially modulate connectivity and (2) are these effects related to responding under varying motor sets. fMRI data were acquired in twenty right-handed volunteers who responded with their right (dominant) or left (non-dominant) hand (blocked across acquisitions). Within acquisitions, the task oscillated between periodic responses (promoting the emergence of motor sets) or randomly induced responses (disrupting the emergence of motor sets). Conjunction analyses revealed eight shared nodes across response hand and condition, time series from which were analyzed. For right hand responses connectivity of the M1 ←→ Thalamus and SMA ←→ Parietal pathways was more significantly modulated during periodic responding. By comparison, for left hand responses, connectivity between five network pairs (including M1 and SMA, insula, basal ganglia, premotor cortex, parietal cortex, thalamus) was more significantly modulated during random responding. uFC analyses were complemented by directed FC based on multivariate autoregressive models of times series from the nodes. These results were complementary and highlighted significant modulation of dFC for SMA → Thalamus, SMA → M1, basal ganglia → Insula and basal ganglia → Thalamus. The results demonstrate complex effects of motor organization and task demand and response hand on different connectivity classes of fMRI data. The brain's sub-networks are flexibly modulated by factors related to motor organization and/or task demand, and our results have implications for assessment of medical conditions associated with motor dysfunction.

Determining Excitatory and Inhibitory Neuronal Activity from Multimodal fMRI Data Using a Generative Hemodynamic Model

Hemodynamic responses, in general, and the blood oxygenation level-dependent (BOLD) fMRI signal, in particular, provide an indirect measure of neuronal activity. There is strong evidence that the BOLD response correlates well with post-synaptic changes, induced by changes in the excitatory and inhibitory (E-I) balance between active neuronal populations. Typical BOLD responses exhibit transients, such as the early-overshoot and post-stimulus undershoot, that can be linked to transients in neuronal activity, but they can also result from vascular uncoupling between cerebral blood flow (CBF) and venous cerebral blood volume (venous CBV). Recently, we have proposed a novel generative hemodynamic model of the BOLD signal within the dynamic causal modeling framework, inspired by physiological observations, called P-DCM (Havlicek et al., 2015). We demonstrated the generative model's ability to more accurately model commonly observed neuronal and vascular transients in single regions but also effective connectivity between multiple brain areas (Havlicek et al., 2017b). In this paper, we additionally demonstrate the versatility of the generative model to jointly explain dynamic relationships between neuronal and hemodynamic physiological variables underlying the BOLD signal using multi-modal data. For this purpose, we utilized three distinct data-sets of experimentally induced responses in the primary visual areas measured in human, cat, and monkey brain, respectively: (1) CBF and BOLD responses; (2) CBF, total CBV, and BOLD responses (Jin and Kim, 2008); and (3) positive and negative neuronal and BOLD responses (Shmuel et al., 2006). By fitting the generative model to the three multi-modal experimental data-sets, we showed that the presence or absence of dynamic features in the BOLD signal is not an unambiguous indication of presence or absence of those features on the neuronal level. Nevertheless, the generative model that takes into account the dynamics of the physiological mechanisms underlying the BOLD response allowed dissociating neuronal from vascular transients and deducing excitatory and inhibitory neuronal activity time-courses from BOLD data alone and from multi-modal data.

Alterations in neurovascular coupling following acute traumatic brain injury

Following acute traumatic brain injury (TBI), timely transport to a hospital can significantly improve the prognosis for recovery. There is, however, a dearth of quantitative biomarkers for brain injury that can be rapidly acquired and interpreted in active, field environments in which TBIs are frequently incurred. We explored potential functional indicators for TBI that can be noninvasively obtained through portable detection modalities, namely optical and electrophysiological approaches. By combining diffuse correlation spectroscopy with colocalized electrophysiological measurements in a mouse model of TBI, we observed concomitant alterations in sensory-evoked cerebral blood flow (CBF) and electrical potentials following controlled cortical impact. Injury acutely reduced the peak amplitude of both electrophysiological and CBF responses, which mostly recovered to baseline values within 30 min, and intertrial variability for these parameters was also acutely altered. Notably, the postinjury dynamics of the CBF overshoot and undershoot amplitudes differed significantly; whereas the amplitude of the initial peak of stimulus-evoked CBF recovered relatively rapidly, the ensuing undershoot did not appear to recover within 30 min of injury. Additionally, acute injury induced apparent low-frequency oscillatory behavior in CBF ([Formula: see text]). Histological assessment indicated that these physiological alterations were not associated with any major, persisting anatomical changes. Several time-domain features of the blood flow and electrophysiological responses showed strong correlations in recovery kinetics. Overall, our results reveal an array of stereotyped, injury-induced alterations in electrophysiological and hemodynamic responses that can be rapidly obtained using a combination of portable detection techniques.

Understanding the dynamic relationship between cerebral blood flow and the BOLD signal: Implications for quantitative functional MRI

Calibrated BOLD imaging, in which traditional measurements of the BOLD signal are combined with measurements of cerebral blood flow (CBF) within a BOLD biophysical model to estimate changes in oxygen metabolism (CMRO2), has been a valuable tool for untangling the physiological processes associated with neural stimulus-induced BOLD activation. However, to date this technique has largely been applied to the study of essentially steady-state physiological changes (baseline to activation) associated with block-design stimuli, and it is unclear whether this approach may be directly extended to the study of more dynamic, naturalistic experimental designs. In this study we tested an assumption underlying this technique whose validity is critical to the application of calibrated BOLD to the study of more dynamic stimuli, that information about fluctuations in venous cerebral blood volume (CBVv) can be captured indirectly by measuring fluctuations in CBF, making the independent measurement of CBVv unnecessary. To accomplish this, simultaneous arterial spin labeling and BOLD imaging were used to measure the CBF and BOLD responses to flickering checkerboards with contrasts that oscillated continuously with frequencies of ~0.02-0.16Hz. The measurements were then fit to a dynamic physiological model of the BOLD response in order to explore the range of consistent CMRO2 and CBVv responses. We found that the BOLD and CBF responses were most consistent with relatively tight dynamic coupling between CBF and CMRO2 and a CBVv response that was an order of magnitude slower than either CBF or CMRO2. This finding suggests that the assumption of tight flow-volume coupling may not be strictly valid, complicating the extension of calibrated BOLD to more naturalistic experimental designs.

The coupling of cerebral blood flow and oxygen metabolism with brain activation is similar for simple and complex stimuli in human primary visual cortex

Quantitative functional MRI (fMRI) experiments to measure blood flow and oxygen metabolism coupling in the brain typically rely on simple repetitive stimuli. Here we compared such stimuli with a more naturalistic stimulus. Previous work on the primary visual cortex showed that direct attentional modulation evokes a blood flow (CBF) response with a relatively large oxygen metabolism (CMRO2) response in comparison to an unattended stimulus, which evokes a much smaller metabolic response relative to the flow response. We hypothesized that a similar effect would be associated with a more engaging stimulus, and tested this by measuring the primary human visual cortex response to two contrast levels of a radial flickering checkerboard in comparison to the response to free viewing of brief movie clips. We did not find a significant difference in the blood flow-metabolism coupling (n=%ΔCBF/%ΔCMRO2) between the movie stimulus and the flickering checkerboards employing two different analysis methods: a standard analysis using the Davis model and a new analysis using a heuristic model dependent only on measured quantities. This finding suggests that in the primary visual cortex a naturalistic stimulus (in comparison to a simple repetitive stimulus) is either not sufficient to provoke a change in flow-metabolism coupling by attentional modulation as hypothesized, that the experimental design disrupted the cognitive processes underlying the response to a more natural stimulus, or that the technique used is not sensitive enough to detect a small difference.

Cortical depth dependence of the BOLD initial dip and poststimulus undershoot in human visual cortex at 7 Tesla

PURPOSE

Owing to variability in vascular dynamics across cerebral cortex, blood-oxygenation-level-dependent (BOLD) spatial and temporal characteristics should vary as a function of cortical-depth. Here, the positive response, initial dip (ID), and post-stimulus undershoot (PSU) of the BOLD response in human visual cortex are investigated as a function of cortical depth and stimulus duration at 7 Tesla (T).

METHODS

Gradient-echo echo-planar-imaging BOLD fMRI with high spatial and temporal resolution was performed in 7 healthy volunteers and measurements of the ID, PSU, and positive BOLD response were made as a function of cortical depth and stimulus duration (0.5-8 s). Exploratory analyses were applied to understand whether functional mapping could be achieved using the ID, rather than positive, BOLD signal characteristics

RESULTS

The ID was largest in outer cortical layers, consistent with previously reported upstream propagation of vasodilation along the diving arterioles in animals. The positive BOLD signal and PSU showed different relationships across the cortical depth with respect to stimulus duration.

CONCLUSION

The ID and PSU were measured in humans at 7T and exhibited similar trends to those recently reported in animals. Furthermore, while evidence is provided for the ID being a potentially useful feature for better understanding BOLD signal dynamics, such as laminar neurovascular coupling, functional mapping based on the ID is extremely difficult.

Measuring cerebrovascular reactivity: what stimulus to use?

Cerebrovascular reactivity is the change in cerebral blood flow in response to a vasodilatory or vasoconstrictive stimulus. Measuring variations of cerebrovascular reactivity between different regions of the brain has the potential to not only advance understanding of how the cerebral vasculature controls the distribution of blood flow but also to detect cerebrovascular pathophysiology. While there are standardized and repeatable methods for estimating the changes in cerebral blood flow in response to a vasoactive stimulus, the same cannot be said for the stimulus itself. Indeed, the wide variety of vasoactive challenges currently employed in these studies impedes comparisons between them. This review therefore critically examines the vasoactive stimuli in current use for their ability to provide a standard repeatable challenge and for the practicality of their implementation. Such challenges include induced reductions in systemic blood pressure, and the administration of vasoactive substances such as acetazolamide and carbon dioxide. We conclude that many of the stimuli in current use do not provide a standard stimulus comparable between individuals and in the same individual over time. We suggest that carbon dioxide is the most suitable vasoactive stimulus. We describe recently developed computer-controlled MRI compatible gas delivery systems which are capable of administering reliable and repeatable vasoactive CO2 stimuli.

Sensitivity of BOLD response to increasing visual contrast: spin echo versus gradient echo EPI

Previous evidence showed that spin-echo (SE) BOLD signals offer an increased linearity and promptness with respect to gradient-echo (GE) acquisition, possibly providing a more accurate estimate of the amplitude of neuronal activity. However there is no evidence that the two sequences differ in representing different activation levels due to changes in stimulus intensity. In this study at 3T we compared GE and SE BOLD responses to visual stimuli at increasing contrast levels (5%, 20%, 60%, and 100%). Both sequences showed a monotonic increase of the BOLD response with stimulus contrast. While the larger sensitivity of GE yielded overall larger signal changes, step-wise increase in activation for GE was significant only when comparing 20% with 5% contrast, whereas for SE a significant increase was observed also when comparing 60% with 20% contrast. Moreover, BOLD responses normalized to the lowest contrast showed that relative increases of SE fMRI signal with increasing stimulus strength are larger with respect to the corresponding values of GE signal. This difference was observed also when excluding voxels attributed to large vessels, suggesting a non negligible role of the extravascular contribution to the modulation of SE fMRI signal with stimulus intensity. These results are shown to be in agreement with theoretical predictions that we derived from a recently proposed model of GE and SE functional signals. The present findings suggest that, despite the limited increase in functional localization accuracy at low field, SE fMRI might offer a potential advantage in distinguishing different levels of stimulus-evoked brain activity.

Coupling of cerebral blood flow and oxygen metabolism is conserved for chromatic and luminance stimuli in human visual cortex

The ratio of the changes in cerebral blood flow (CBF) and cerebral metabolic rate of oxygen (CMRO(2)) during brain activation is a critical determinant of the magnitude of the blood oxygenation level dependent (BOLD) response measured with functional magnetic resonance imaging (fMRI). Cytochrome oxidase (CO), a key component of oxidative metabolism in the mitochondria, is non-uniformly distributed in visual area V1 in distinct blob and interblob regions, suggesting significant spatial variation in the capacity for oxygen metabolism. The goal of this study was to test whether CBF/CMRO(2) coupling differed when these subpopulations of neurons were preferentially stimulated, using chromatic and luminance stimuli to preferentially stimulate either the blob or interblob regions. A dual-echo spiral arterial spin labeling (ASL) technique was used to measure CBF and BOLD responses simultaneously in 7 healthy human subjects. When the stimulus contrast levels were adjusted to evoke similar CBF responses (mean 65.4% ± 19.0% and 64.6% ± 19.9%, respectively for chromatic and luminance contrast), the BOLD responses were remarkably similar (1.57% ± 0.39% and 1.59% ± 0.35%) for both types of stimuli. We conclude that CBF-CMRO(2) coupling is conserved for the chromatic and luminance stimuli used, suggesting a consistent coupling for blob and inter-blob neuronal populations despite the difference in CO concentration.

Luminance contrast of a visual stimulus modulates the BOLD response more than the cerebral blood flow response in the human brain

The blood oxygenation level dependent (BOLD) response measured with functional magnetic resonance imaging (fMRI) depends on the evoked changes in cerebral blood flow (CBF) and cerebral metabolic rate of oxygen (CMRO(2)) in response to changes in neural activity. This response is strongly modulated by the CBF/CMRO(2) coupling relationship with activation, defined as n, the ratio of the fractional changes. The reliability of the BOLD signal as a quantitative reflection of underlying physiological changes depends on the stability of n in response to different stimuli. The effect of visual stimulus contrast on this coupling ratio was tested in 9 healthy human subjects, measuring CBF and BOLD responses to a flickering checkerboard at four visual contrast levels. The theory of the BOLD effect makes a robust prediction-independent of details of the model-that if the CBF/CMRO(2) coupling ratio n remains constant, then the response ratio between the lowest and highest contrast levels should be higher for the BOLD response than the CBF response because of the ceiling effect on the BOLD response. Instead, this response ratio was significantly lower for the BOLD response (BOLD response: 0.23 ± 0.13, mean ± SD; CBF response: 0.42 ± 0.18; p=0.0054). This data is consistent with a reduced dynamic range (strongest/weakest response ratio) of the CMRO(2) response (~1.7-fold) compared to that of the CBF response (~2.4-fold) as luminance contrast increases, corresponding to an increase of n from 1.7 at the lowest contrast level to 2.3 at the highest contrast level. The implication of these results for fMRI studies is that the magnitude of the BOLD response does not accurately reflect the magnitude of underlying physiological processes.

Robust quantification of microvascular transit times via linear dynamical systems using two-photon fluorescence microscopy data

Vascular transit time is an important indicator of microcirculatory health. We present a second-order-plus-dead-time (SOPDT) model for robust estimation of kinetic parameters characterizing microvascular bolus passage using two-photon fluorescence microscopy (2PFM) in anesthetized rats receiving somatosensory stimulation. This methodology enables quantification of transit time, time-to-peak, overshoot, and rate of bolus passage through the microvascular network. The overall transit time during stimulation, of 2.2 ± 0.1 seconds, was shorter (P ~ 0.0008) than that at rest (2.7 ± 0.2 seconds). When compared with conventional γ-variate modeling, the SOPDT modeling yielded better quality of fit both at rest (P<0.0001) and on activation (P<0.001).

Biophysical and physiological origins of blood oxygenation level-dependent fMRI signals

After its discovery in 1990, blood oxygenation level-dependent (BOLD) contrast in functional magnetic resonance imaging (fMRI) has been widely used to map brain activation in humans and animals. Since fMRI relies on signal changes induced by neural activity, its signal source can be complex and is also dependent on imaging parameters and techniques. In this review, we identify and describe the origins of BOLD fMRI signals, including the topics of (1) effects of spin density, volume fraction, inflow, perfusion, and susceptibility as potential contributors to BOLD fMRI, (2) intravascular and extravascular contributions to conventional gradient-echo and spin-echo BOLD fMRI, (3) spatial specificity of hemodynamic-based fMRI related to vascular architecture and intrinsic hemodynamic responses, (4) BOLD signal contributions from functional changes in cerebral blood flow (CBF), cerebral blood volume (CBV), and cerebral metabolic rate of O(2) utilization (CMRO(2)), (5) dynamic responses of BOLD, CBF, CMRO(2), and arterial and venous CBV, (6) potential sources of initial BOLD dips, poststimulus BOLD undershoots, and prolonged negative BOLD fMRI signals, (7) dependence of stimulus-evoked BOLD signals on baseline physiology, and (8) basis of resting-state BOLD fluctuations. These discussions are highly relevant to interpreting BOLD fMRI signals as physiological means.

Cellular Links between Neuronal Activity and Energy Homeostasis

Neuronal activity, astrocytic responses to this activity, and energy homeostasis are linked together during baseline, conscious conditions, and short-term rapid activation (as occurs with sensory or motor function). Nervous system energy homeostasis also varies during long-term physiological conditions (i.e., development and aging) and with adaptation to pathological conditions, such as ischemia or low glucose. Neuronal activation requires increased metabolism (i.e., ATP generation) which leads initially to substrate depletion, induction of a variety of signals for enhanced astrocytic function, and increased local blood flow and substrate delivery. Energy generation (particularly in mitochondria) and use during ATP hydrolysis also lead to considerable heat generation. The local increases in blood flow noted following neuronal activation can both enhance local substrate delivery but also provides a heat sink to help cool the brain and removal of waste by-products. In this review we highlight the interactions between short-term neuronal activity and energy metabolism with an emphasis on signals and factors regulating astrocyte function and substrate supply.

Contributions of dynamic venous blood volume versus oxygenation level changes to BOLD fMRI

Blood-oxygenation-level-dependent (BOLD) fMRI has contributions from venous oxygenation and venous cerebral blood volume (CBV) changes. To examine the relative contribution of venous CBV change (ΔCBV(v)) to BOLD fMRI, BOLD and arterial CBV changes (ΔCBV(a)) to a 40-s forepaw stimulation in six α-chloralose anesthetized rats were measured using a magnetization transfer-varied fMRI technique, while total CBV change (ΔCBV(t)) was measured with injection of iron oxide nanoparticles. ΔCBV(v) was obtained by subtracting ΔCBV(a) from ΔCBV(t). We observed a fast ΔCBV(a) response with a time constant of 2.9 ± 2.3s and a slower ΔCBV(v) response with a time constant of 13.5 ± 5.7s and an onset delay of 6.1 ± 3.3s. These results are consistent with earlier studies under different anesthesia and stimulus, supporting that fast CBV(a) and slow CBV(v) responses are generalizable. Assuming the observed post-stimulus BOLD undershoot is at least partly explained by the ΔCBV(v) contribution, the relative contribution of the ΔCBV(v)- and oxygenation-change-related components to the BOLD response was estimated. The relative ΔCBV(v) contribution increases with time during stimulation; whereby it is <0.14 during initial 10s and reaches a maximum possible value of ~0.45 relative to the oxygenation contribution during the 30-40s period after stimulus onset. Our data indicates that the contribution of venous oxygenation change to BOLD fMRI is dominant for short stimulations.

IRON fMRI measurements of CBV and implications for BOLD signal

Changes in cerebral blood volume (CBV) and blood magnetization each induce changes in the transverse relaxation rate of MRI signal that are associated with changes in cerebral activity. BOLD signal, the preeminent method for non-invasive localization of task-induced brain function in human subjects, reflects a combination of changes in CBV and blood magnetization. Intravenous injection of paramagnetic contrast media, usually iron oxide particles surrounded by larger macromolecules, can overwhelm the BOLD response and sensitize signal to blood plasma volume, a method we have deemed "IRON" fMRI. The practical advantage of this technique is the ability to optimize blood magnetization at any echo time, enabling high detection power and the use of short echo times; for these reasons, IRON fMRI has become a valuable imaging tool in animal models. The temporal response of blood plasma volume is quite different from blood flow and BOLD signal; thus, CBV has been identified as a prominent source of transient features of the BOLD response. This article reviews the methodological advantages of the IRON method and how CBV measurements have informed our understanding of the BOLD response.

The BOLD post-stimulus undershoot, one of the most debated issues in fMRI

This paper provides a brief overview of how we got involved in fMRI work and of our efforts to elucidate the mechanisms underlying BOLD signal changes. The phenomenon discussed here in particular is the post-stimulus undershoot (PSU), the interpretation of which has captivated many fMRI scientists and is still under debate to date. This controversy is caused both by the convoluted physiological origin of the BOLD effect, which allows many possible explanations, and the lack of comprehensive data in the early years. BOLD effects reflect changes in cerebral blood flow (CBF), volume (CBV), metabolic rate of oxygen (CMRO(2)), and hematocrit fraction (Hct). However, the size of such effects is modulated by vascular origin such as intravascular, extravascular, macro and microvascular, venular and capillary, the relative contributions of which depend not only on the spatial resolution of the measurements, but also on stimulus duration, on magnetic field strength and on whether spin echo (SE) or gradient echo (GRE) detection is used. The two most dominant explanations of the PSU have been delayed vascular compliance (first venular, later arteriolar, and recently capillary) and sustained increases in CMRO(2), while post-activation reduction in CBF is a distant third. MRI has the capability to independently measure CBF and arteriolar, venous, and total CBV contributions in humans and animals, which has been of great assistance in improving the understanding of BOLD phenomena. Using currently available MRI and optical data, we conclude that the predominant PSU origin is a sustained increase in CMRO(2). However, some contributions from delayed vascular compliance are likely, and small CBF undershoot contributions that are difficult to detect with current arterial spin labeling technology can also not be excluded. The relative contribution of these different processes, which are not mutually exclusive and can act together, is likely to vary with stimulus duration and type.

Quantitative functional MRI: concepts, issues and future challenges

Since its inception 20 years ago, functional magnetic resonance imaging (fMRI) of the human brain based on the blood oxygenation level dependent (BOLD) contrast phenomenon has proliferated and matured. Today it is the predominant functional brain imaging modality with the majority of applications being in basic cognitive neuroscience where it has primarily been used as a tool to localize brain activity. While the magnitude of the BOLD response is often used in these studies as a surrogate for the level of neuronal activity, the link between the two is, in fact, quite indirect. The BOLD response is dependent upon hemodynamic (blood flow and volume) and metabolic (oxygen consumption) responses as well as acquisition details. Furthermore, the relationship between neuronal activity and the hemodynamic response, termed neurovascular coupling, is itself complex and incompletely understood. Quantitative fMRI techniques have therefore been developed to measure the hemodynamic and metabolic responses to modulations in brain activity. These methods have not only helped clarify the behaviour and origins of the BOLD signal under normal physiological conditions but they have also provided a potentially valuable set of tools for exploring pathophysiological conditions. Such quantitative methods will be critical to realize the potential of fMRI in a clinical context, where simple BOLD measurements cannot be uniquely interpreted, and to enhance the power of fMRI in basic neuroscience research. In this article, recent advances in human quantitative fMRI methods are reviewed, outstanding issues discussed and future challenges and opportunities highlighted.

Linear coupling of undershoot with BOLD response in ER-fMRI and nonlinear BOLD response in rapid-presentation ER-fMRI

In event-related (ER) BOLD-fMRI brain activation studies, understanding the relationship between the elicited BOLD signal and its underlying neuronal activity is essential for any quantitative interpretation of the neural events from the BOLD measurements. This requires a better understanding of the dynamic BOLD response. Besides the neuronal activity-induced positive BOLD response, the dynamic response is also characterized by a profound post-stimulus undershoot. The relationship between the positive response and the post-stimulus undershoot, however, remains poorly understood. Earlier studies using block-design paradigms with long stimulation durations (>10s) do not suggest a quantitative relationship. Using an ER paradigm, this study revealed a linear coupling between the positive BOLD response and the post-stimulus undershoot across the human visual cortex. The voxelwise linear coupling across the visual cortex strongly supports a homogeneous hemodynamic response in ER paradigms, though the BOLD response magnitude varies substantially over a wide range across the visual cortex. Although underlying neuronal activity is responsible for a BOLD response, the blood volume fraction affects the magnitude of the BOLD response; the larger the blood volume fraction, the larger the magnitude. This effect needs to be accounted for in any quantitative interpretation of the BOLD measurements. In the absence of nonlinear neuronal activities, the nonlinear vascular response renders the estimated BOLD responses smaller in rapid presentation (RP) ER paradigms compared to that in ER paradigms, and this reduction effect also needs to be considered when interpreting the estimated BOLD responses in RP-ER paradigms. Interestingly, this nonlinear effect might be simply accounted for by a scaling factor across the visual cortex.

Exploring the post-stimulus undershoot with spin-echo fMRI: implications for models of neurovascular response

As a consequence of neural stimulation the blood oxygenation-level dependent (BOLD) contrast in gradient-echo echo-planar imaging (GE-EPI) based functional MRI (fMRI) leads to an increased MR signal in activated brain regions. Following this, a BOLD signal undershoot below baseline is generally observed with GE-EPI. The origin of this undershoot has been the focus of many investigations using fMRI and optical modalities, but the underlying mechanisms remain disputed. Here, we investigate the BOLD undershoot following visual stimulation by using a purely T₂-weighted fMRI sequence at 1.5 and 3 T. By taking advantage of the field strength dependency of the T₂ BOLD contrast and complete absence of static dephasing effects due to the pure spin echoes, one can draw conclusions about the origin of the BOLD undershoot and test the hypotheses in the literature. We observe a significant undershoot at both field strengths, with constant undershoot-to-main response ratio. This provides strong evidence that the undershoot is caused by BOLD changes due to elevated post-stimulus deoxyhaemoglobin concentration in the small vessels. 'Delayed vascular compliance' as suggested by the well-known Balloon and Windkessel models does not appear capable of explaining the undershoot. Our results also suggest that blood volume changes in arterioles and capillaries, for which there is consistent evidence from optical imaging studies, cannot alone cause the undershoot. This has important implications for models of neurovascular response and provides further support for the decoupling of changes in the rate of oxygen metabolism and blood flow. In addition, we found that an 'arteriolar balloon' (delayed arterial compliance) may provide a plausible explanation for the temporal characteristics of the BOLD undershoot.

Basal cerebral blood volume during the poststimulation undershoot in BOLD MRI of the human brain

One of the characteristics of the blood oxygenation level-dependent (BOLD) magnetic resonance imaging (MRI) response to functional challenges of the brain is the poststimulation undershoot, which has been suggested to originate from a delayed recovery of either cerebral blood volume (CBV) or cerebral metabolic rate of oxygen to baseline. Using bolus-tracking MRI in humans, we recently showed that relative CBV rapidly normalizes after the end of stimulation. As this observation contradicts at least part of the blood-pool contrast agent studies performed in animals, we reinvestigated the CBV contribution by dynamic T1-weighted three-dimensional MRI (8 seconds temporal resolution) and Vasovist at 3 T (12 subjects). Initially, we determined the time constants of individual BOLD responses. After injection of Vasovist, CBV-related T1-weighted signal changes revealed a signal increase during visual stimulation (1.7% ± 0.4%), but no change relative to baseline in the poststimulation phase (0.2% ± 0.3%). This finding renders the specific nature of the contrast agent unlikely to be responsible for the discrepancy between human and animal studies. With the assumption of normalized cerebral blood flow after stimulus cessation, a normalized CBV lends support to the idea that the BOLD MRI undershoot reflects a prolonged elevation of oxidative metabolism.

Improved fMRI calibration: precisely controlled hyperoxic versus hypercapnic stimuli

The calibration of functional magnetic resonance imaging (fMRI) for the estimation of neuronal activation-induced changes in cerebral metabolic rate of oxygen (CMRO(2)) has been achieved through hypercapnic-induced iso-metabolic increases in cerebral blood flow (CBF). Hypercapnia (HC) has been traditionally implemented through alterations in the fixed inspired fractional concentrations of carbon dioxide (F(I)CO(2)) without otherwise controlling end-tidal partial pressures of carbon dioxide (P(ET)CO(2)) or oxygen (P(ET)O(2)). There are several shortcomings to the use of this manual HC method that may be improved by using precise targeting of P(ET)CO(2) while maintaining iso-oxia. Similarly, precise control of blood gases can be used to induce isocapnic hyperoxia (HO) to reduce venous deoxyhaemoglobin (dHb) and thus increase BOLD signals, without appreciably altering CMRO(2) or CBF. The aim of our study was to use precise end-tidal targeting to compare the calibration of BOLD signals under an isocapnic hyperoxic protocol (HOP) (rises in P(ET)O(2) to 140, 240 and 340 mm Hg from baseline) to that of an iso-oxic hypercapnic protocol (HCP) (rises in P(ET)CO(2) of 3, 5, 7 and 9 mm Hg from baseline). Nine healthy volunteers were imaged at 3T while monitoring end-tidal gas concentrations and simultaneously measuring BOLD and CBF signals, via arterial spin labeling (ASL), during graded HCP and HOP, alternating with normocapnic states in a blocked experimental design. The variability of the calibration constant obtained under HOP (M(HOP)) was 0.3-0.5 that of the HCP one (M(HCP)). In addition, M-variances with precise gas targeting (M(HCP) and M(HOP)) were less than those reported in studies using traditional F(I)CO(2) and F(I)O(2) methods (M(HC) and M(HO), respectively). We conclude that precise controlled gas delivery markedly improves BOLD-calibration for fMRI studies of oxygen metabolism with both the HCP and the more precise HOP-alternative.

Cerebral blood flow measurement using fMRI and PET: a cross-validation study

An important aspect of functional magnetic resonance imaging (fMRI) is the study of brain hemodynamics, and MR arterial spin labeling (ASL) perfusion imaging has gained wide acceptance as a robust and noninvasive technique. However, the cerebral blood flow (CBF) measurements obtained with ASL fMRI have not been fully validated, particularly during global CBF modulations. We present a comparison of cerebral blood flow changes (ΔCBF) measured using a flow-sensitive alternating inversion recovery (FAIR) ASL perfusion method to those obtained using H₂¹⁵O PET, which is the current gold standard for in vivo imaging of CBF. To study regional and global CBF changes, a group of 10 healthy volunteers were imaged under identical experimental conditions during presentation of 5 levels of visual stimulation and one level of hypercapnia. The CBF changes were compared using 3 types of region-of-interest (ROI) masks. FAIR measurements of CBF changes were found to be slightly lower than those measured with PET (average ΔCBF of 21.5 ± 8.2% for FAIR versus 28.2 ± 12.8% for PET at maximum stimulation intensity). Nonetheless, there was a strong correlation between measurements of the two modalities. Finally, a t-test comparison of the slopes of the linear fits of PET versus ASL ΔCBF for all 3 ROI types indicated no significant difference from unity (P > .05).

Quantification of CMRO(2) without hypercapnia using simultaneous near-infrared spectroscopy and fMRI measurements

Estimation of the cerebral metabolic rate of oxygen (CMRO(2)) and cerebral blood flow (CBF) is important to investigate the neurovascular coupling and physiological components in blood oxygenation level-dependent (BOLD) signals quantitatively. Although there are methods that can determine CMRO(2) changes using functional MRI (fMRI) or near-infrared spectroscopy (NIRS), current approaches require a separate hypercapnia calibration process and have the potential to incur bias in many assumed model parameters. In this paper, a novel method to estimate CMRO(2) without hypercapnia is described using simultaneous measurements of NIRS and fMRI. Specifically, an optimization framework is proposed that minimizes the differences between the two forms of the relative CMRO(2)-CBF coupling ratio from BOLD and NIRS biophysical models, from which hypercapnia calibration and model parameters are readily estimated. Based on the new methods, we found that group average CBF, CMRO(2), cerebral blood volume (CBV), and BOLD changes within activation of the primary motor cortex during a finger tapping task increased by 39.5 +/- 21.4%, 18.4 +/- 8.7%, 12.9 +/- 6.7%, and 0.5 +/- 0.2%, respectively. The group average estimated flow-metabolism coupling ratio was 2.38 +/- 0.65 and the hypercapnia parameter was 7.7 +/- 1.7%. These values are within the range of values reported from other literatures. Furthermore, the activation maps from CBF and CMRO(2) were well localized on the primary motor cortex, which is the main target region of the finger tapping task.

Neural activity-induced modulation of BOLD poststimulus undershoot independent of the positive signal

Despite intense research on the blood oxygenation level-dependent (BOLD) signal underlying functional magnetic resonance imaging, our understanding of its physiological basis is far from complete. In this study, it was investigated whether the so-called poststimulus BOLD signal undershoot is solely a passive vascular effect or actively induced by neural responses. Prolonged static and flickering black-white checkerboard stimulation with isoluminant grey screen as baseline condition were employed on eight human subjects. Within the same region of interest, the positive BOLD time courses for static and flickering stimuli were identical over the entire stimulus duration. In contrast, the static stimuli exhibited no poststimulus BOLD signal undershoot, whereas the flickering stimuli caused a strong BOLD poststimulus undershoot. To ease the interpretation, we performed an additional study measuring both BOLD signal and cerebral blood flow (CBF) using arterial spin labeling. Also for CBF, a difference in the poststimulus period was found for the two stimuli. Thus, a passive blood volume effect as the only contributor to the poststimulus undershoot comes short in explaining the BOLD poststimulus undershoot phenomenon for this particular experiment. Rather, an additional active neuronal activation or deactivation can strongly modulate the BOLD poststimulus behavior. In summary, the poststimulus time course of BOLD signal could potentially be used to differentiate neuronal activity patterns that are otherwise indistinguishable using the positive evoked response.

BOLD-specific cerebral blood volume and blood flow changes during neuronal activation in humans

To understand and predict the blood-oxygenation level-dependent (BOLD) fMRI signal, an accurate knowledge of the relationship between cerebral blood flow (DeltaCBF) and volume (DeltaCBV) changes is critical. Currently, this relationship is widely assumed to be characterized by Grubb's power-law, derived from primate data, where the power coefficient (alpha) was found to be 0.38. The validity of this general formulation has been examined previously, and an alpha of 0.38 has been frequently cited when calculating the cerebral oxygen metabolism change (DeltaCMRo(2)) using calibrated BOLD. However, the direct use of this relationship has been the subject of some debate, since it is well established that the BOLD signal is primarily modulated by changes in 'venous' CBV (DeltaCBV(v), comprising deoxygenated blood in the capillary, venular, and to a lesser extent, in the arteriolar compartments) instead of total CBV, and yet DeltaCBV(v) measurements in humans have been extremely scarce. In this work, we demonstrate reproducible DeltaCBV(v) measurements at 3 T using venous refocusing for the volume estimation (VERVE) technique, and report on steady-state DeltaCBV(v) and DeltaCBF measurements in human subjects undergoing graded visual and sensorimotor stimulation. We found that: (1) a BOLD-specific flow-volume power-law relationship is described by alpha = 0.23 +/- 0.05, significantly lower than Grubb's constant of 0.38 for total CBV; (2) this power-law constant was not found to vary significantly between the visual and sensorimotor areas; and (3) the use of Grubb's value of 0.38 in gradient-echo BOLD modeling results in an underestimation of DeltaCMRo(2).

The Possible Role of CO(2) in Producing A Post-Stimulus CBF and BOLD Undershoot

Comprehending the underlying mechanisms of neurovascular coupling is important for understanding the pathogenesis of neurodegenerative diseases related to uncoupling. Moreover, it elucidates the casual relation between the neural signaling and the hemodynamic responses measured with various imaging modalities such as functional magnetic resonance imaging (fMRI). There are mainly two hypotheses concerning this mechanism: a metabolic hypothesis and a neurogenic hypothesis. We have modified recent models of neurovascular coupling adding the effects of both NO (nitric oxide) kinetics, which is a well-known neurogenic vasodilator, and CO₂ kinetics as a metabolic vasodilator. We have also added the Hodgkin–Huxley equations relating the membrane potentials to sodium influx through the membrane. Our results show that the dominant factor in the hemodynamic response is NO, however CO₂ is important in producing a brief post-stimulus undershoot in the blood flow response that in turn modifies the fMRI blood oxygenation level-dependent post-stimulus undershoot. Our results suggest that increased cerebral blood flow during stimulation causes CO₂ washout which then results in a post-stimulus hypocapnia induced vasoconstrictive effect.

Linearity of the fMRI response in category-selective regions of human visual cortex

The goal of this study was to determine the linearity of the blood oxygen level-dependent (BOLD) response, as measured by functional magnetic resonance imaging (fMRI), in category-selective regions of human visual cortex. We defined regions of the temporal lobe that were selective to faces (fusiform face area, FFA) and places (parahippocampal place area, PPA). We then determined the linearity of the BOLD response in these regions to their preferred and nonpreferred stimuli. First, we tested the principle of scaling. As we increased the visibility of the stimulus, there was a corresponding linear increase in the fMRI signal in the FFA and PPA to their preferred stimulus (face and place, respectively). In contrast, responses in the FFA and PPA to the nonpreferred stimulus did not conform to the principle of scaling. Next, we asked whether the fMRI response in these regions of visual cortex conformed to the principle of additivity. To assess this, we determined whether the response to a long stimulus block could be predicted by adding the response to multiple shorter duration blocks. Although the fMRI response in the FFA and PPA was generally linear to the preferred stimulus, a more nonlinear response was apparent to the nonpreferred stimulus. In conclusion, the linearity of the BOLD response in the human ventral visual pathway varied across cortical region and stimulus category. This suggests that measures of linearity may provide a useful indication of neural selectivity in the brain.

Estimating cerebral oxygen metabolism from fMRI with a dynamic multicompartment Windkessel model

Stimulus evoked changes in cerebral blood flow, volume, and oxygenation arise from responses to underlying neuronally mediated changes in vascular tone and cerebral oxygen metabolism. There is increasing evidence that the magnitude and temporal characteristics of these evoked hemodynamic changes are additionally influenced by the local properties of the vasculature including the levels of baseline cerebral blood flow, volume, and blood oxygenation. In this work, we utilize a physiologically motivated vascular model to describe the temporal characteristics of evoked hemodynamic responses and their expected relationships to the structural and biomechanical properties of the underlying vasculature. We use this model in a temporal curve-fitting analysis of the high-temporal resolution functional MRI data to estimate the underlying cerebral vascular and metabolic responses in the brain. We present evidence for the feasibility of our model-based analysis to estimate transient changes in the cerebral metabolic rate of oxygen (CMRO(2)) in the human motor cortex from combined pulsed arterial spin labeling (ASL) and blood oxygen level dependent (BOLD) MRI. We examine both the numerical characteristics of this model and present experimental evidence to support this model by examining concurrently measured ASL, BOLD, and near-infrared spectroscopy to validate the calculated changes in underlying CMRO(2).

Origins of the BOLD post-stimulus undershoot

The interpretation of the blood-oxygenation level-dependent (BOLD) post-stimulus undershoot has been a topic of considerable interest, as the mechanisms behind this prominent BOLD transient may provide valuable clues on the neurovascular response process and energy supply routes of the brain. Biomechanical theories explain the origin of the BOLD undershoot through the passive ballooning of post-capillary vessels which leads to an increase in venous blood volume (CBV(v), comprising deoxygenated blood in capillary, venular and arteriolar compartments), resulting in susceptibility-induced signal decrease. While there has been substantial evidence supporting a role for venous ballooning, there have also been reports arguing for a prolonged post-stimulus elevation in cerebral oxygenation consumption (CMRo(2)) as the primary cause. Furthermore, a contribution of post-stimulus cerebral blood flow (CBF) undershoots has also been demonstrated. To clarify the role of the venous compartment in causing the BOLD undershoot, we performed in vivo fMRI measurements of the transient DeltaCBV(v), DeltaCBF and DeltaBOLD responses in healthy humans. We observed a slow post-stimulus return to baseline in venous CBV which supports the existence of a passive "balloon" effect, implying that previous observations of a quicker recovery of the total CBV response may be dominated by arterial CBV change. Our findings also support a significant contribution from the CBF undershoots, which, combined with a slow venous CBV response, would account for much of the BOLD undershoot.

Caffeine's effects on cerebrovascular reactivity and coupling between cerebral blood flow and oxygen metabolism

The blood-oxygenation-level-dependent (BOLD) signal is dependent on multiple physiological factors such as cerebral blood flow (CBF), local oxygen metabolism (CMRO(2)) and cerebral blood volume (CBV). Since caffeine affects both CBF and neural activity, its effects on BOLD remain controversial. The calibrated BOLD approach is an excellent tool to study caffeine because it combines CBF and BOLD measures to estimate changes in CMRO(2). The present study used the calibrated BOLD approach with 5% CO(2) to determine if a 2.5 mg/kg intravenous injection of caffeine changes the coupling between CBF and CMRO(2) during motor and visual tasks. The results show that caffeine decreases n, the CBF:CMRO(2) coupling ratio, from 2.58 to 2.33 in motor (p=0.006) and from 2.45 to 2.23 in visual (p=0.002) areas respectively. The current study also demonstrated that caffeine does not alter cerebrovascular reactivity to CO(2). These results highlight the importance of the calibrated BOLD approach in improving interpretation of the BOLD signal in the presence of substances like caffeine.

Identification and comparison of stochastic metabolic/hemodynamic models (sMHM) for the generation of the BOLD signal

This paper extends a previously formulated deterministic metabolic/hemodynamic model for the generation of blood oxygenated level dependent (BOLD) responses to include both physiological and observation stochastic components (sMHM). This adds a degree of flexibility when fitting the model to actual data by accounting for un-modelled activity. We then show how the innovation method can be used to estimate unobserved metabolic/hemodynamic as well as vascular variables of the sMHM, from simulated and actual BOLD data. The proposed estimation method allowed for doing model comparison by calculating the model's AIC and BIC. This methodology was then used to select between different neurovascular coupling assumptions underlying sMHM. The proposed framework was first validated on simulations and then applied to BOLD data from a motor task experiment. The models under comparison in the analysis of the actual data considered that vascular response was coupled to: (I) inhibition, (II) excitation, (III) both excitation and inhibition. Data was best described by model II, although model III was also supported.

Differentiating sensitivity of post-stimulus undershoot under diffusion weighting: implication of vascular and neuronal hierarchy

The widely used blood oxygenation level dependent (BOLD) signal during brain activation, as measured in typical fMRI methods, is composed of several distinct phases, the last of which, and perhaps the least understood, is the post-stimulus undershoot. Although this undershoot has been consistently observed, its hemodynamic and metabolic sources are still under debate, as evidences for sustained blood volume increases and metabolic activities have been presented. In order to help differentiate the origins of the undershoot from vascular and neuronal perspectives, we applied progressing diffusion weighting gradients to investigate the BOLD signals during visual stimulation. Three distinct regions were established and found to have fundamentally different properties in post-stimulus signal undershoot. The first region, with a small but focal spatial extent, shows a clear undershoot with decreasing magnitude under increasing diffusion weighting, which is inferred to represent intravascular signal from larger vessels with large apparent diffusion coefficients (ADC), or high mobility. The second region, with a large continuous spatial extent in which some surrounds the first region while some spreads beyond, also shows a clear undershoot but no change in undershoot amplitude with progressing diffusion weighting. This would indicate a source based on extravascular and small vessel signal with smaller ADC, or lower mobility. The third region shows no significant undershoot, and is largely confined to higher order visual areas. Given their intermediate ADC, it would likely include both large and small vessels. Thus the consistent observation of this third region would argue against a vascular origin but support a metabolic basis for the post-stimulus undershoot, and would appear to indicate a lack of sustained metabolic rate likely due to a lower oxygen metabolism in these higher visual areas. Our results are the first, to our knowledge, to suggest that the post-stimulus undershoots have a spatial dependence on the vascular and neuronal hierarchy, and that progressing flow-sensitized diffusion weighting can help delineate these dependences.

Cortical layer-dependent dynamic blood oxygenation, cerebral blood flow and cerebral blood volume responses during visual stimulation

The spatiotemporal characteristics of cerebral blood volume (CBV) and flow (CBF) responses are important for understanding neurovascular coupling mechanisms and blood oxygenation level-dependent (BOLD) signals. For this, cortical layer-dependent BOLD, CBV and CBF responses were measured at the cat visual cortex using fMRI. Major findings are: (i) the time-dependent fMRI cortical profile is dependent on imaging modality. Overall, the peak across the cortex occurs at the cortical surface for BOLD, but at the middle cortical layer for CBV and CBF. Compared to an initial stimulation period (4-10 s), the spatial specificity of CBV to the middle cortical layer increases significantly at a later time, while the specificity of BOLD and CBF slightly changes. (ii) The CBV response at the upper cortical area containing large pial vessels has a faster onset time and time to peak than the BOLD response at the same area, and a faster time to peak than CBV at the middle cortical area with microvessels. This suggests that the dilation of microvessels at the middle cortical area follows arterial volume increase at the surface of the cortex. (iii) For all three modalities, the post-stimulus undershoot was observed with the 60-s stimulation paradigm, indicating that the post-stimulus BOLD undershoot cannot be explained by the delayed venous CBV recovery theory under our experimental conditions. (iv) The relationship between CBV and CBF responses is both spatially and temporally dependent. Thus, a single power-law scaling constant (gamma value) may not be applicable for high-resolution study.

Assessment of human brain temperature by 1H MRS during visual stimulation and hypercapnia

Brain temperature is determined by the interplay between the cerebral metabolic rate of oxygen (CMRO2) and cerebral blood flow (CBF). In this study, single-voxel 1H nuclear MRS, with an accuracy of +/-0.2 degrees C for temperature determination, was used at 3 T to measure human brain temperature during visual stimulation (which increases both CBF and CMRO2) and hypercapnia (which increases CBF only). Visual stimulation had no detectable effect on brain temperature in the parenchyma showing blood oxygenation level dependent activation. Hypercapnia, leading to an increase in the end tidal CO2 by 8 +/- 2 mm Hg above the baseline, caused a short-lasting decrease in brain temperature of 0.30 +/- 0.33 degrees C. These results indicate that increased CBF may be a key factor, bringing about a small decrease in brain temperature during brain activation. However, the increase in CBF is not sufficient to lower brain temperature in the presence of a concomitant increase in endogenous heat production.

Neural substrates of narrative comprehension and memory

When reading a narrative, comprehension and retention of information benefit considerably from the use of situation models--coherent representations of the characters, locations, and activities described in the text. Here we used functional magnetic resonance imaging (fMRI) to explore the neural mechanisms supporting situation model processing. Participants read blocks of sentences that were either unrelated to one another or formed coherent narratives. A timecourse-based approach was used to identify regions that differentiated narrative-level comprehension from sentence-level comprehension. Most brain regions that showed modulation of activation during narrative-level comprehension were also modulated to a lesser extent during sentence-level comprehension, suggesting a shared reliance on general coherence-building mechanisms. However, tentative evidence was found for narrative-specific activation in dorsomedial prefrontal cortex. Additional analyses identified spatiotemporally distinct neural contributions to situation model processing, with posterior parietal regions supporting situation model construction and frontotemporal regions supporting situation model maintenance. Finally, a set of subsequent memory analyses demonstrated that the boost in comprehension and memory performance observed for coherent materials was attributable to the use of integrative situation models rather than lower-level differences in sentence-level or word-level encoding. These results clarify the functional contributions of distinct brain systems to situation model processing and their mapping onto existing psychological models of narrative comprehension.

The post-stimulation undershoot in BOLD fMRI of human brain is not caused by elevated cerebral blood volume

Functional magnetic resonance imaging (fMRI) based on blood oxygenation level dependent (BOLD) contrast is the most widely used technique for imaging human brain function. However, the dynamic interplay of altered cerebral blood flow (CBF), cerebral blood volume (CBV), and oxidative metabolism (CMRO2) is not yet fully understood. One of the characteristics of the BOLD response is the post-stimulation undershoot, that is increased deoxyhemoglobin, which has been suggested to originate from a delayed recovery of elevated CBV or CMRO2 to baseline. To investigate the CBV contribution to the post-stimulation BOLD undershoot, we performed bolus-tracking experiments using a paramagnetic contrast agent in eight healthy subjects at 3 T. In an initial BOLD experiment without contrast agent, we determined the individual hemodynamic responsiveness. In two separate experiments, we then evaluated the relative CBV (rCBV) during visual stimulation and the post-stimulation undershoot, respectively. The results confirm a pronounced rCBV increase during stimulation (31.4+/-8.6%), but reveal no change in rCBV relative to baseline in the post-stimulation phase (0.7+/-7.2%). This finding renders a CBV contribution to the BOLD MRI undershoot unlikely and--in conjunction with a rapid post-stimulation return of CBF to baseline--supports the idea of a prolonged elevation of oxidative metabolism.

Vascular space occupancy weighted imaging with control of residual blood signal and higher contrast-to-noise ratio

It has been recently proposed that the local cerebral blood volume change during brain activation can be measured by a series of images whose contrast is dependent on vascular space occupancy (VASO). VASO takes advantage of the inversion recovery sequence to acquire images when the longitudinal magnetization (Mz) of blood is relaxing through zero. The degree of blood suppression, however, is not always well controlled as a consequence of spatial variations in inversion efficiency and blood T1. Furthermore, while blood is eliminated, the Mz of other tissues is also small, which makes the contrast-to-noise ratio inherently low in VASO. In this paper, diffusion gradients were applied to demonstrate residual intravascular signal in the original VASO. An alternative VASO-weighted imaging was then proposed using a longer inversion time at which the Mz difference between blood and gray matter was optimized. A global saturation immediately after image acquisition was employed to eliminate the Mz disparity between inflowing blood and the residual in-plane blood from previous acquisition. Feasibility was evaluated by numerical simulation and functional experiments. In human visual cortex, the fractional VASO signal and cerebral blood volume changes were found to be -0.6% and 44%, respectively (voxel size = 3.4 x 3.4 x 5.0 mm3). As compared to the original VASO, the presented method provided a largely comparable activation map and hemodynamic curve but was not confounded by the existence of blood. Results also demonstrated its advantages of 1.6-fold higher CNR and insensitivity to variant tissue/blood T1 as well as inversion efficiency.

A calibration method for quantitative BOLD fMRI based on hyperoxia

The estimation of changes in CMR(O2) using functional MRI involves an essential calibration step using a vasoactive agent to induce an isometabolic change in CBF. This calibration procedure is performed most commonly using hypercapnia as the isometabolic stimulus. However, hypercapnia possesses a number of detrimental side effects. Here, a new method is presented using hyperoxia to perform the same calibration step. This procedure requires independent measurement of Pa(O2), the BOLD signal, and CBF. We demonstrate that this method yields results that are comparable to those derived using other methods. Further, the hyperoxia technique is able to provide an estimate of the calibration constant that has lower overall intersubject and intersession variability compared to the hypercapnia approach.