Biophysical and Physiological Origins of Blood Oxygenation Level-Dependent fMRI Signals

Abstract

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₂ utilization (CMRO₂), (5) dynamic responses of BOLD, CBF, CMRO₂, 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.

Introduction

In 1990, Ogawa and colleagues at AT&T Bell Laboratories studied responses to global physiological stimulus in rat brains at 7 T and reported that functional brain mapping was possible using venous blood oxygenation level-dependent (BOLD) magnetic resonance imaging (MRI) contrast (Ogawa et al, 1990a, 1990b; Ogawa and Lee, 1990). This BOLD contrast relies on changes in deoxyhemoglobin, which acts as an endogenous paramagnetic contrast agent (Pauling and Coryell 1936). Therefore, changes in the local deoxyhemoglobin concentration in the brain lead to alterations in MRI signal intensity (Thulborn et al, 1982; Ogawa et al, 1990a, 1990b; Ogawa and Lee, 1990). Application of BOLD contrast to human functional brain mapping soon followed (Bandettini et al, 1992; Kwong et al, 1992; Ogawa et al, 1992). In recent years, functional MRI (fMRI) has been the tool of choice to visualize neural activity in the human brain, with >3,000 papers with the keyword ‘fMRI’ published annually. Thus, it is critical to understand the origins of BOLD fMRI responses.

Figure 1 depicts the components on cellular (Figure 1A) and vascular (Figure 1B) levels affected by neural activity induced by task execution or stimulus exposure, which give rise to BOLD fMRI signals (Figure 1C). Neural activity can be broadly classified by local field potentials (LFPs) and spiking activity; LFP is generally believed to represent synaptic activity including neural input, while spiking activity in cell bodies represents suprathreshold neural output. Increased neural activity is accompanied by an increased flux of Na⁺, K⁺, and Ca⁺², and increased ATP production via glucose consumption. These changes induce astrocytes (yellow in Figure 1A) and neuronal cells (round green regions in Figure 1A) to send vasoactive signals into nearby arterioles and capillaries (dark bands in Figure 1A), consequently dilating the upstream arterial vessels. The exact signaling mechanism remains an area of intense research, with excellent review articles available (Attwell et al, 2010; Petzold and Murthy, 2011), as this tight neurovascular coupling is a key to understanding vascular-based fMRI approaches. Also relevant (and as yet undetermined) are the contributions of neuronal versus astrocyte signaling, the contributions of excitatory versus inhibitatory neurotransmission, and the spatial extent and dynamic properties of neural versus vascular responses.

Figure 1

Relationship between neurovascular components and BOLD fMRI. (A) Cells are differentiated in this in-vivo rat brain image from two-photon laser scanning microscopy (provided by Alberto Vazquez at the University of Pittsburgh), where neuronal cell bodies and processes appear green (Oregon Green BAPTA-1AM), astrocyte cell bodies and processes are yellow (Sulforhodamine 101), and capillaries appear as dark bands. Increased activity in the neurons and astrocytes near and surrounding the capillaries induces an increase in capillary blood flow and dilation of upstream arterioles. (B) Vascular structure is depicted in this image from human cortex obtained by vascular casting and electron scanning microscopy, as adapted from Figure 4 in Reina-De La Torre et al (1998). The region outlined by the red rectangle approximates the size of the cellular image in (A). A dense network of small-diameter intracortical vessels appears below the dashed horizontal line, while a smaller number of large-diameter pial vessels appear above the dashed horizontal line. Blood flows from the pial arteries to the intracortical arteriolar branches, then into the dense capillary network, draining back through the intracortical venules and pial veins. (C) Whole-brain human BOLD fMRI studies depict functional responses localized to multiple brain regions (Kim and Ugurbil, 1997b). Functional MRI was acquired by the gradient-echo echo planar imaging technique with spatial resolution of 3.1 × 3.1 × 5 mm³ at 4 T. The size of each single fMRI pixel is approximated by the boundaries of the vascular image in (B); thus, there are numerous vessels of different sizes within each fMRI pixel, and contributions from these vessels must be considered to understand BOLD biophysical mechanisms. BOLD, blood oxygenation level dependent; MRI, magnetic resonance imaging; fMRI, functional MRI.

A focal increase in cerebral blood flow (CBF) can be considered to directly relate to neuronal activity because glucose metabolism and CBF changes are closely coupled (Raichle, 1987). However, based on the positron emission tomographic measurements of CBF and cerebral metabolic rate of oxygen utilization (CMRO₂) in humans during somatosensory and visual stimulation (Fox and Raichle, 1986; Fox et al, 1988), the increase in CBF surpasses the increase in CMRO₂, resulting in an increase in capillary and venous oxygenation levels. Regardless of the exact signaling mechanisms for this CBF response, the BOLD signal is related to a mismatch between CBF increase and CMRO₂ change. Thus, the relationship between neural activity and BOLD fMRI signals must be understood to properly interpret the maps. It has been reported that the BOLD fMRI response correlates (1) with underlying LFP rather than spiking activity (Logothetis et al, 2001; Viswanathan and Freeman, 2007), (2) with mostly spiking activity (Heeger and Ress, 2002), or (3) with both LFP and spiking activity (Mukamel et al, 2005). Because spiking activity occurs together with LFP, in most circumstances it is difficult to separate their contributions to the BOLD response (Yen et al, 2011). The exact relationship between BOLD signals and underlying neural activity is still a matter of debate, even after decades of research (Ekstrom, 2010), but excellent review articles on this topic are also available (Logothetis, 2008; Ekstrom, 2010).

In this article, we will focus on the biophysical and physiological sources of BOLD fMRI signals and discuss what is currently known about stimulus-evoked vascular changes. Since BOLD fMRI measures signal changes induced by neural activity, we will evaluate the contribution of various fMRI signal components. We will also discuss nonBOLD fMRI approaches that are sensitive to changes in single parameters such as perfusion or blood volume. Topics also reviewed are the spatial specificity of BOLD fMRI, physiological sources of stimulus-evoked fMRI responses, and resting-state fMRI. Since the BOLD fMRI field is now >20 years old, many excellent review articles are available on human applications and neurophysiology (e.g., Logothetis, 2008; van Eijsden et al, 2009; Glover, 2011). Our review will therefore focus instead on findings based primarily on our animal studies, which afford the systematic assessment of physiological sources and the combination of fMRI with invasive approaches. Topics discussed here, however, should also be relevant to human fMRI studies.

Vascular architecture relevant to hemodynamic-based functional magnetic resonance imaging

Since hemodynamic-based fMRI signals originate from vessels, understanding the vascular structure is highly important. Thus, we review what is known about vascular structure and density within human cortex, based on an influential paper by Duvernoy (Duvernoy et al, 1981). Vessels can be broadly classified as either pial or intracortical (see vessels above and below the dashed red horizontal line in Figure 1B, respectively). Pial vessels are larger and less dense than intracortical vessels, and connect to penetrating arteries and emerging veins. Intracortical vessels can be subdivided into arteries, veins, and the capillary network. Capillaries with ∼5 μm average diameter and ∼100 μm length (Pawlik et al, 1981) are close to active neurons and astrocytes as seen in Figure 1A; the distance between capillaries and neurons is likely <13 μm, since the average intercapillary distance is ∼25 μm (Pawlik et al, 1981). Intracortical arteries and veins (see Figure 2A) can be further classified by the depth of their cortical penetration (Duvernoy et al, 1981); vessels in groups 1 and 2 (arterial diameters of 10 to 25 μm, venous diameters of 20 to 30 μm) penetrate only the upper cortical layers (layers 1 to 3), group 3 vessels (arterial diameters of 15 to 30 μm, venous diameter of 45 μm) penetrate the middle of the cortex (layers 3 to 5), group 4 vessels (arterial diameters of 30 to 40 μm, venous diameter of 65 μm) penetrate lower cortical depths (layer 6), and vessels in group 5 (arterial diameters of 30 to 75 μm, venous diameters of 80 to 125 μm) vascularize both cortex and white matter. Note that group 6 arterial vessels (A6 in Figure 2A) pass through cortex and vascularize only white matter. Venous vasculature can be imaged using BOLD contrast (often referred to as ‘susceptibility-weighted imaging’). Figure 2B illustrates this BOLD contrast, where veins appear as dark lines or dots because T₂^∗ in venous blood is shorter than in tissue and arterial blood, and also because the blood susceptibility effect extension to extravascular tissue enlarges the apparent venous vessel size; groups 3 to 5 intracortical veins (e.g., red arrows) can be clearly visualized with this technique. Typical distances between magnetic resonance (MR)-detectable neighboring intracortical veins are ∼0.5 to 1 mm in rats and cats (Park et al, 2008). In the middle and deep cortical layers, intracortical arteries outnumber intracortical veins by a factor of ∼4 (Duvernoy et al, 1981). The size of a single vascular unit corresponds to the amount of tissue it supplies or drains. These vascular unit diameters are 0.33 to 0.5 mm for groups 2 to 3 arteries, 0.75 to 1 mm for groups 2 to 3 veins, 0.5 to 2 mm for group 5 arteries, and 1 to 4 mm for group 5 veins (Duvernoy et al, 1981; see Figure 2C).

Figure 2

Vascular architecture relevant to hemodynamic-based fMRI. (A) This schematic illustrates a section from the cortical surface (top) to white matter (bottom) with the various types of arterial (A) and venous (V) intracortical vessels as classified by penetration depth (adapted from Figure 27 in Duvernoy et al, 1981). Group 6 (arteries only) extends to white matter with no cortical branches, group 5 also extends to white matter, but has cortical branches, group 4 extends to cortical layer 6 (highlighted in green), group 3 extends to layers 3 to 5 (layer 4 highlighted in yellow), and groups 1 and 2 extend only to the upper cortical layers. (B) This venogram is a 3-D T₂^∗-weighted coronal MR image of cat brain obtained at 9.4 T with 78 mm isotropic resolution, 1.25 mm-thick slab selection, and minimum intensity projection performed to enhance the contrast of venous vessels; data acquisition and processing methods are reported elsewhere (Park et al, 2008). The bottom image is an expansion of signal within the dashed yellow box. White matter (outlined by dashed blue lines) is clearly differentiated from gray matter. Venous vessels draining from white matter (indicated by red arrows) belong to group 5. (C) This schematic of a tangential section (parallel to the cortical surface) shows venous vascular units (blue dots with black traces) with their arterial rings (red dots) (adapted from Figure 63 in Duvernoy et al, 1981), where the size of a single venous vascular unit (typically 1 to 2 mm in diameter) corresponds to the amount of tissue it drains. MR, magnetic resonance; MRI, MR imaging; fMRI, functional MRI.

Contribution of physiological parameters to baseline magnetic resonance imaging signals

Examining individual contributions to MRI signals is an important step toward understanding BOLD fMRI mechanisms. Pixel sizes in typical fMRI studies are a few millimeters; each pixel may therefore include blood, extravascular tissue, and cerebrospinal fluid (CSF) components. Since arterial blood and venous blood have different T₂ values, these will be considered separately, with the assumption that capillary content is partly arterial blood and partly venous blood. Thus, MR signal intensity from a pixel is the sum of signals originating from multiple compartments with different spin density and relaxation parameters. The MRI intensity, S, can be described as

$$S=\sum {\rho }_i\times V_i\times M_{\mathrm{s}\mathrm{s},i}\times e^{-\mathrm{T}\mathrm{E}/T_{2,i}^{\ast }}$$

(1)

where subscripts, i, indicate the compartment; ρ is the water proton spin density; V is the volume fraction; M_SS is the steady-state magnetization, which is $(1-e^{-\text{TR}/T_1^{\ast }})\mathit{\text{sin}}\theta /(1-\mathit{\text{cos}}\theta \times e^{-\text{TR}/T_1^{\ast }})$ , where TR is the repetition time, T₁^∗ is the apparent longitudinal relaxation time in the presence of inflow, and θ is the flip angle; TE is the echo time; and T₂^∗ is the apparent transverse relaxation time including the BOLD-contrast susceptibility effect. Therefore, fMRI signal changes are dependent not only on imaging parameters (including TR and θ), but also on biophysical responses that significantly change any of the above parameters. It is therefore critical to understand all potentially significant contributions to MRI signals, including those that are not BOLD effects.

Water Proton Spin Density (ρ) and Volume Fractions (V)

Water proton spin density is directly related to a water content, which is 0.88 g/mL for extravascular tissue (obtained from the gray matter value in Table 1 in Herscovitch and Raichle, 1985) and 0.87 g/mL for arterial and venous blood with a hematocrit level of 35% (see further details in Table 2 in Herscovitch and Raichle, 1985). Arterial and venous blood compartment values for V are ∼1% (Ito et al, 2001) and 2.5% to 3% (An and Lin, 2002a, 2002b), respectively, in human brain.

Inflow and Cerebral Blood Flow (T₁^∗)

Inflow and CBF can act to decrease the T₁^∗ values of blood and extravascular tissue components. Flowing blood moves spins from outside the imaging plane into the slice pixels. When spins in the imaging plane are saturated, signal from the unsaturated inflowing blood is enhanced relative to the surrounding stationary spins; this effect is the basis of time-of-flight angiography, which visualizes arterial blood vessels. The magnitude of this inflow contribution in vessels depends on the MRI parameters, TR and θ. When the TR allows arterial blood spins from outside the imaging slice to flow into capillaries and exchange with extravascular tissue water is sufficiently long (e.g., >1 second), spatially specific perfusion contrast appears. In the extravascular tissue pool, 1/T₁^∗=1/T₁+f/λ where f is CBF (mL blood/g tissue per second) and λ is the blood-to-tissue partition coefficient (mL blood/g tissue) (Detre et al, 1992). Both inflow and CBF effects are negligible only when TR>T₁^∗. At 3 T, gray matter T₁ is ∼1.5 seconds, while human gray matter f is ∼0.01 mL/g tissue per second and λ is ∼0.96 mL/g (Herscovitch and Raichle, 1985), giving a typical baseline T₁^∗ value of ∼1.48 seconds.

Blood Oxygenation Level-Dependent Relaxation Parameters (T₂ and T₂^∗)

R₂^∗ (=1/T₂^∗) is the sum of R₂ (=1/T₂) and susceptibility-induced R₂′ (=1/T₂′), where R₂^∗ and R₂ can be measured by gradient-echo (GE) and spin-echo (SE) MR techniques, respectively. Blood water R₂ and R₂^∗ values are directly related to paramagnetic deoxyhemoglobin content (Thulborn et al, 1982). Water rapidly exchanges between deoxyhemoglobin-containing red blood cells and plasma, and also diffuses in the presence of the magnetic field gradients generated by the deoxyhemoglobin inside the vessel. These exchange and diffusion processes result in a loss of phase coherence, referred to as ‘dynamic’ (time irreversible) averaging, and water molecules inside a given vessel experience these processes similarly. In addition to the deoxyhemoglobin content-related R₂ change, a frequency change is also observed, which is dependent on magnetic field, oxygen saturation level, and the angle between vessel direction and B₀ (Ogawa et al, 1993b). Since multiple vessels at different orientations typically exist within any given pixel, their multiple frequency shifts cause a phase dispersion (rather than a net phase shift), and a reduction in blood T₂^∗.

The field gradient generated by deoxyhemoglobin decreases by (r/a)², where r is the distance from vessel to the region of interest and a is the vessel radius. Thus, the susceptibility effect extends to extravascular tissue and nearby CSF. Additionally, the susceptibility effect depends on vessel orientation; there is no extravascular susceptibility effect around vessels running parallel to B₀, while the maximum effect occurs around vessels orthogonal to B₀ (Ogawa and Lee, 1990). Water proton spins dephase depending on its local susceptibility effect during TE; spins nearby a vein dephase larger than those far away. Thus, multiple frequency shifts within a given pixel result in phase dispersion and consequently MRI signal loss (see dark vessels in Figure 2B). In a pixel with numerous randomly oriented venous vessels, R₂′ in extravascular tissue is closely related to the amount of deoxyhemoglobin (Ogawa et al, 1993b; Kennan et al, 1994; Weisskoff et al, 1994; Yablonskiy and Haacke, 1994; Boxerman et al, 1995). A simple biophysical model of susceptibility-induced R₂′ is described here for ease of understanding, but a more complete theoretical description can be found in recent articles (Frohlich et al, 2005; Dickson et al, 2011). Based on Monte Carlo simulation (Ogawa et al, 1993b; Boxerman et al, 1995), the dephasing effect within a pixel can be approximated as R₂′, which is A·CBV_v{Δχ₀ω₀(1−Y)} ^β where CBV_v is venous cerebral blood volume, Δχ₀ is the susceptibility difference between fully deoxygenated blood and fully oxygenated blood, ω₀ is the water resonance frequency (=γB₀, where γ is the gyromagnetic ratio for the proton nucleus), and A and β are constants. The constant, A, is dependent on vessel size and orientation (note that α was used in Ogawa et al, 1993b), while the constant, β, has values between 1 (for large vessels) and 2 (for capillaries). Susceptibility effects (and therefore R₂^∗ and R₂′) will increase whenever there is (1) an increase in venous CBV and consequent increase in deoxyhemoglobin content in the voxel, (2) a decrease in venous oxygen saturation level, or (3) an increase in magnetic field (ω₀).

Non-blood oxygenation level-dependent contribution to functional magnetic resonance imaging signals

The relative contributions from spin density/volume fraction, inflow/CBF (T₁^∗), and BOLD effects (T₂ and T₂^∗) are all dependent on the imaging parameters. Figure 3 includes a schematic of vascular responses to increased neural activity, and summarizes the expected contribution of these changes to MR-related parameters.

Figure 3

Vascular responses to neural activity and expected changes to MR-related parameters. (A) This schematic adapted from Kim and Fukuda (2008) shows oxyhemoglobin (red dots) and deoxyhemoglobin (blue dots) in blood flowing through arteries, arterioles, capillaries, venules, and finally to veins. Oxygen delivered via oxyhemoglobin diffuses into extravascular tissue, and where it is used as a metabolic substrate. At prestimulus baseline conditions, blood oxygen saturation is ∼100% in arteries, while it is ∼60% in veins, although actual values vary with physiological condition. Increases in neural activity trigger an increase in blood velocity (indicated by the size of arrows) and dilation of vessels. The resulting increase in perfusion exceeds what is required by the increase in oxygen consumption rate, and the consequence is higher oxyhemoglobin levels and lower deoxyhemoglobin levels, such that capillary and venous oxygenation levels actually increase. (B) This table graphically illustrates the expected changes in the MR-related parameters of equation (1) for different compartments within each pixel due to vascular responses to neural activity. For simplicity, it is assumed that capillaries are composed of both arterial and venous blood. Upward arrows indicate an increase during stimulus, downward arrows indicate a decrease, and horizontal arrows indicate no change. BOLD, blood oxygenation level dependent; MR, magnetic resonance.

Effect of Water Proton Spin Density (ρ) and Volume Fractions (V)

Stimulus-induced dilation of vessels increases CBV, thus increasing the volume fractions of arterial and venous blood. Under the constraint of constant intracranial brain volume dictated by the Monro–Kellie doctrine, there must be a corresponding reduction of the volume fractions of extravascular tissue and CSF. Indeed, the vascular space occupancy technique (Lu et al, 2003) showed that increases in CBV reduce the tissue volume fraction in pixels without CSF. Cortical surface boundaries may then expand, reducing the CSF volume fraction (Jin and Kim, 2010).

In fMRI studies, it may be not possible to separate changes in ρ from changes in V, since both components are independent of TE and TR. There is no strong evidence that ρ increases with stimulation, but any change in V could be interpreted as a change in ρ. Stroman et al (2003) proposed a nonBOLD contrast mechanism for fMRI based on the changes in extravascular spin density, and called it signal enhancement by extravascular water protons. In fast-spin-echo BOLD fMRI studies of human spinal cord, TE-dependent signal changes were nonlinear, indicating the existence of more than one compartment (Stroman et al, 2005). Spin-echo BOLD fMRI studies of the cerebral cortex in animals and humans also show nonlinear TE-dependent signal changes (similar to Stroman's findings). However, when vascular signal contributions were removed with flow-crushing gradients (Yacoub et al, 2003; Jochimsen et al, 2005; Jin et al, 2006), the TE-dependent signal changes were linear, indicating that the nonBOLD contrast observed in vivo is not due to a change in spin density, but rather to changes in vascular contributions.

Effect of Inflow and Cerebral Blood Flow (T₁^∗)

Stimulus-induced increases in blood velocity/inflow in pixels containing large vessels act to reduce blood T₁^∗, thereby increasing the fMRI response (Segebarth et al, 1994). This large-vessel inflow contribution is undesirable in fMRI studies, as it decreases spatial specificity. If the CBF contribution is significant in MRI signals, then any stimulus-induced increases in CBF will decrease T₁^∗ in extravascular tissue and increase fMRI signals. However, in typical fMRI studies with ∼3 × 3 × 3 mm³ resolution and 2 to 3 seconds temporal resolution, stimulus-induced nonBOLD effects are not likely significant contributions to fMRI signal changes.

Blood oxygenation level-dependent contributions to gradient-echo and spin-echo functional magnetic resonance imaging signals

To understand the biophysical basis of BOLD-based fMRI, it is important to examine intravascular and extravascular contributions to the BOLD signal. Note that more complete BOLD modeling of intravascular and extravascular contributions has been reported recently (Uludag et al, 2009; Griffeth and Buxton, 2011), which is beyond the scope of this review.

Intravascular Blood Oxygenation Level-Dependent Components

The intravascular BOLD signal almost linearly increases when TE<blood T₂ for SE studies or TE<blood T₂^∗ for GE studies becomes maximal at TE=T₂ or T₂^∗ of blood, and decreases thereafter (Duong et al, 2003; Jin et al, 2006; Uludag et al, 2009). This intravascular BOLD contribution can be estimated from blood relaxation times for different experimental conditions. Blood R₂ (=1/T₂) is quadratically dependent on the deoxygenation level (1−Y, where Y is blood oxygen saturation fraction) and magnetic field (B₀) (Wright et al, 1991). T₂ values of blood water for Y=0.6 are ∼127 ms at 1.5 T (Wright et al, 1991), ∼31 ms at 3 T (Zhao et al, 2007), and ∼12 to 15 ms at 7 T (Ogawa et al, 1993a), while T₂ values of gray matter water are 90 ms at 1.5 T (Breger et al, 1989), 71 ms at 3 T (Gelman et al, 1999), and 55 ms at 7 T (Yacoub et al, 2001). Note that in-vivo blood T₂ and T₂^∗ values are typically shorter than the values measured in vitro due to the additional inflow effects during measurements. As the magnetic field increases, venous blood T₂ and T₂^∗ values decrease faster than corresponding tissue T₂ and T₂^∗ values; therefore at ultrahigh fields, the intravascular component can be significantly reduced.

Extravascular Blood Oxygenation Level-Dependent Components

During the TE typically chosen for fMRI studies (∼50 ms), water molecules diffuse 10 to 20 μm. Thus, phase coherence from extravascular tissue water near capillaries will be dynamically averaged over the steep field gradients produced by blood deoxyhemoglobin, while near large veins it will be locally averaged across smaller field gradients (see schematic in Kim and Ugurbil, 2003). The dephasing effect around large veins can be refocused by a 180° radiofrequency pulse. In GE BOLD fMRI, the extravascular effect around small vessels is supralinearly dependent on B₀, while the extravascular effect around large vessels is linearly dependent on B₀ (Ogawa et al, 1993b); thus, at higher B₀, the extravascular contribution of small vessel-related signals increases; however, even at 9.4 T, the extravascular effect around large vessels is still significant in GE BOLD fMRI studies (Zhao et al, 2004). These extravascular contributions from large vessels can be reduced with SE techniques. Both SE and GE extravascular BOLD fMRI signals are almost linearly dependent on TE within the commonly used TE ranges (e.g., Yacoub et al, 2003; Jin et al, 2006).

Separation of Intravascular and Extravascular Blood Oxygenation Level-Dependent Components

The relative contribution of intravascular versus extravascular signals to BOLD fMRI responses can be examined using small diffusion gradients (e.g., b value <200 to 300 s/mm²; Song et al, 1996). Diffusion-weighted gradients induce velocity-dependent phase shifts in the presence of flow, and consequently reduce signals from blood due to the inhomogeneous velocities within a vessel and the presence of blood vessels at different orientations within a pixel. Increasing the magnitude of diffusion weighting progressively attenuates the intravascular signals from faster flowing to slower flowing blood. Table 1 summarizes intravascular fractions (=1 −BOLD with diffusion gradients/BOLD without diffusion gradients) from the literature. It should be noted that small diffusion gradients may not totally remove the intravascular contributions from small vessels (including capillaries), and thus the literature intravascular fraction values are likely to be slightly underestimated. In summary, the intravascular contribution decreases as magnetic field strength increases (1.5 T versus 7 and 9.4 T) and as TE increases (e.g., at 4 T, the intravascular fractions are 0.75 (TE=32 ms) and 0.2 (TE=64 ms)). As expected, pial vessel regions have higher intravascular contributions than cortical regions (see Figure 3D in Jin et al, 2006); even at 9.4 T, when TE=16 ms, the intravascular fraction from pial vessel regions can be up to 60%.

Table 1

Summary of intravascular attenuation by diffusion gradients in BOLD studiesa

Magnetic field (T)	MR technique, echo time (ms)	b value (s/mm 2 )	Subject, brain region	Intravenous fraction	Reference
1.5	ASE, 165	10–690	Human, V1	0.5–0.7	Boxerman et al (1995)
	SE, 125	200–700	Human, V1	0.6–1.0	Zhong et al (1998) (Figure 4B)
	GE, 82	165–660	Human, M1	0.6–0.9	Song et al (1996)
	GE, 33–93	100	Human, V1	0.46–0.62b	Donahue et al (2011) c
3.0	SE, 96.5	14–454	Human, V1	∼0.5	Jochimsen et al (2004)
	GE, 32.7–70.7	100	Human, V1	0.22–0.4b	Donahue et al (2011) c
4.0	SE, 32 and 65	40–400	Human, V1	0.2–0.75b	Duong et al (2003)
	GE, 80	20–500	Human, V1	∼0	Duong et al (2003)
7.0	SE, 32 and 55	40–400	Human, V1	0.1–0.2b	Duong et al (2003)
	GE, 55	20–500	Human, V1	∼0	Duong et al (2003)
	GE, 28.6–46.6	100	Human, V1	0.08–0.16b	Donahue et al (2011) c
9.4	SE, 40	30–500	Rat, S1	∼0	Lee et al (1999)
	SE, 16–70	200	Cat, V1	0–0.6b	Jin et al (2006) d

ASE, asymmetric spin echo; SE, spin echo; GE, gradient echo; V1, visual cortex; M1, motor cortex; S1, somatosensory cortex; TE, echo time; BOLD, blood oxygenation level-dependent; MR, magnetic resonance.

a

Intravenous fraction is determined as (1–the percent BOLD signal with diffusion gradients normalized by the percent BOLD signal without diffusion gradients).

b

Higher intravenous fraction values were obtained from the lowest TE values.

c

Region of interest (ROI) data from Table 2 of Donahue et al (2011).

d

Although two ROIs were selected (a region with pial veins at the cortical surface and a cortical region containing venous vessels, see Figure 3 in Jin et al, 2006), only data from cortical surface ROIs are reported here to include a larger intravenous contribution.

Spatial specificity and sensitivity of blood oxygenation level-dependent functional magnetic resonance imaging

Spatial resolution of BOLD fMRI is determined not only by intrinsic hemodynamic responses, but also by the vascular structure. One important issue for extremely high-resolution studies is that pixels with different baseline CBV_v values will have different BOLD fMRI signal changes for similar oxygenation level changes. To compare BOLD signals across pixels, it is therefore important that pixels have similar vascular properties (both density and diameters).

In BOLD fMRI studies, draining of deoxyhemoglobin from the capillaries (spatially specific to regions of neural activation) to intracortical veins and then to pial veins is one major concern of spatial specificity. As the diameters of these downstream vessels progressively increase, this deoxygenated blood combines with blood draining from other sources. The deoxyhemoglobin change in these draining veins is therefore dependent on their distance from the activated area (Turner, 2002), and on the spatial extent of neural response. When the area of activation is small, this deoxygenated blood is quickly diluted with blood from inactivated regions, effectively reducing the oxygenation level change in adjacent regions. Thus, GE BOLD fMRI is capable of differentiating small functional modules such as single whisker barrels (Yang et al, 1996). However, when a large brain area responds to stimulus, much of the blood draining into the large downstream vessels also has origins from activated regions, so in addition to the spatially specific contributions, there will also be nonspecific contributions from large vessels. This is the case in the GE BOLD fMRI example of cat visual cortex in Figure 4, where nonspecific contributions are seen in pial vessels at the cortical surface and in CSF regions (yellow pixels in green-outlined regions). These nonspecific signal changes can be significantly reduced with SE BOLD fMRI (Figure 4B), improving the spatial specificity to small-sized vessels and nearby tissue, whose location correlates well with the dilations of cortical vessels as seen in the results of CBV fMRI using exogenous intravascular contrast agent with high susceptibility (Kennan et al, 1998; Mandeville et al, 1998; van Bruggen et al, 1998; see Figure 4C).

Figure 4

High-resolution fMRI of cat visual cortex. GE BOLD, SE BOLD, and CBV fMRI studies were performed at 9.4 T with 156 μm in-plane resolution and 2 mm slice thickness (Zhao et al, 2006); CBV fMRI was obtained after intravascular injection of long half-life iron oxide nanoparticles. During full-field binocular visual stimulus, it is expected that layer 4 will have the highest metabolic and CBF responses (Lowel et al, 1987), as well as the highest synaptic density and highest cytochrome oxidase activity (Woolsey et al, 1996). If an fMRI technique is highly specific to metabolic response and neural activity, then the middle of the cortex should therefore shows the highest signal change. In these studies, the highest CBV change (C) indeed occurred at the middle of the visual cortex in layer 4. GE BOLD fMRI (A) has the highest percent signal change at the cortical surface, where large pial vessels are located (green contours), while large vessel contributions are suppressed in SE BOLD (B). BOLD, blood oxygenation level dependent; fMRI, functional magnetic resonance imaging; GE, gradient echo; SE, spin echo; CBV, cerebral blood volume; TE, echo time; CBF, cerebral blood flow.

The SE BOLD fMRI technique at high magnetic field is thus preferable for mapping submillimeter functional structures if the nonspecific large vessel contribution is significant. Note that SE BOLD fMRI with small diffusion gradients improves spatial specificity (similarly to the effect of acquiring SE BOLD fMRI at high magnetic field). But since SE BOLD signal sensitivity is significantly reduced due to refocusing of the dephasing effect around large vessels, and since the additional refocusing radiofrequency pulse increases power deposition to tissue, the GE BOLD technique is often the tool of choice. Large vessel contributions can also be reduced by postprocessing approaches as described below. The location of large pial venous vessels can also be determined from venographic images obtained with high-resolution T₂^∗-weighted MR techniques (see Figure 2B), or can be inferred from structures adjacent to the cortex such as sulci and CSF determined from anatomical images. Large venous vessel areas tend to show large BOLD percentage changes (see Figure 4A), delayed responses (Lee et al, 1995), significant phase changes (Menon, 2002), and large baseline fluctuations (Kim et al, 1994). Importantly, although relative signal changes are high in large vessel regions (Figure 4A), absolute signal changes (and statistical values) are reduced at high magnetic fields due to reduced baseline signal intensities within and near large veins, as discussed in the prior section ‘Intravascular BOLD components’.

The ultimate spatial specificity of hemodynamic-based fMRI ultimately depends on how finely the blood flow is regulated. If intracortical arteries are independently regulated, then spatial resolution can be as fine as 0.33 to 2 mm for arterial blood-based fMRI measurements (Duvernoy et al, 1981). If precapillary arterioles regulate the hemodynamic changes, then achievable spatial resolution is even better. Studies by fMRI suggest that intrinsic CBF and CBV fMRI responses are reasonably specific to submillimeter functional domains (Duong et al, 2001; Zhao et al, 2005), which in cat cortex are on the order of 0.5 to 0.7 mm in diameter (roughly corresponding to the diameter of the arterial unit for groups 2 to 3 arteries). Recent papers indicate that precapillary arterioles indeed dilate during stimulus via astrocyte-capillary signaling (see review articles, Attwell et al, 2010; Petzold and Murthy, 2011). In fact, capillary network responses precisely outline regions of neural activity in the rat olfactory bulb (Chaigneau et al, 2003), suggesting that spatial resolution of hemodynamic-based fMRI can be on the order of ∼100 μm if only capillary signals are detected. Intrinsic spatial resolution for the venous vessel-based BOLD response is ∼1 mm (the diameter of the venous unit for groups 2 to 3 veins). If only small-sized venous vessels are detected, then BOLD spatial resolution improves; however, its spatial specificity is still poorer than techniques based on capillary changes.

Physiological sources of blood oxygenation level-dependent functional magnetic resonance imaging signals

The functional change in venous blood oxygenation level (Y) is dependent on the mismatch between CBF and CMRO₂ changes. Assuming arterial Y=1.0 and constant hematocrit, the relative change in venous Y can be determined from relative changes to both CBF and CMRO₂ as

$$\frac{\mathrm{\Delta }Y}{1-Y}=\frac{\mathrm{\Delta }\mathrm{C}\mathrm{B}\mathrm{F}/\mathrm{C}\mathrm{B}\mathrm{F}-\mathrm{\Delta }{\mathrm{C}\mathrm{M}\mathrm{R}\mathrm{O}}_2/{\mathrm{C}\mathrm{M}\mathrm{R}\mathrm{O}}_2}{(\mathrm{\Delta }\mathrm{C}\mathrm{B}\mathrm{F}/\mathrm{C}\mathrm{B}\mathrm{F}+1)}$$

(2)

where Δ refers to the stimulus-induced change. An increase in venous oxygenation level decreases the blood OEF (oxygen extraction fraction) from blood to tissue (1−Y=OEF). When relative changes in CBF and CMRO₂ are similar, ΔY will be close to zero, but when the relative change in CBF is much larger than the relative change in CMRO₂, ΔY will be highly correlated with the CBF change.

Assuming that hematocrit level in venous blood does not change with stimulus, the BOLD signal change in extravascular tissue can be approximated (Davis et al, 1998) as

$$\mathrm{\%}\mathrm{B}\mathrm{O}\mathrm{L}\mathrm{D}=M(1-{\left(\frac{(1+\mathrm{\Delta }{\mathrm{C}\mathrm{M}\mathrm{R}\mathrm{O}}_2/{\mathrm{C}\mathrm{M}\mathrm{R}\mathrm{O}}_2)}{(1+\mathrm{\Delta }\mathrm{C}\mathrm{B}\mathrm{F}/\mathrm{C}\mathrm{B}\mathrm{F})}\right)}^{\beta }(1+\frac{\mathrm{\Delta }\mathrm{C}\mathrm{B}\mathrm{V}}{\mathrm{C}\mathrm{B}\mathrm{V}}))$$

(3)

where M is a constant, related to baseline physiological, vascular, and imaging parameters. Value for β of 1.5 is commonly assumed (Davis et al, 1998). The first physiological term of equation (3) relates to the mismatch between relative oxygen consumption and CBF change, while the second physiological term relates to the relative CBV change, calculated from Grubb's formula as (1+ΔCBV/CBV)=(1+ΔCBF/CBF) ^α , where α is the flow-volume converting power term (Grubb et al, 1974). An α value of 0.38 was obtained in anesthetized monkeys under global challenge (Grubb et al, 1974) and has been used for most human BOLD modeling studies (e.g., Davis et al, 1998).

Calibration constant M is highly dependent on imaging sequences (SE versus GE) and parameters (TE and B₀), thus should not directly compare across literature values. M can be simulated assuming normal physiological and vascular parameter values (Ogawa et al, 1993b; Kim and Ugurbil, 1997a; Hyder et al, 2001), or can be experimentally determined by inducing hemodynamic changes without metabolic changes (Davis et al, 1998). Calibration studies with hypercapnia (Davis et al, 1998; Hoge et al, 1999b; Kim et al, 1999), hyperoxia (Chiarelli et al, 2007), or breath hold (Kastrup et al, 1999) measured CBF and BOLD to determine M from equation (3) with the assumption that ΔCMRO₂=0. The fMRI community commonly determines oxygen metabolic changes from calibration with equation (3). The choice of initial assumed β value does not change the calculated ΔCMRO₂/CMRO₂ because the determined M factor is also adjusted accordingly (Davis et al, 1998). However, one important underlying assumption is that the total CBV change obtained from the CBF change with Grubb's equation is the same as the BOLD-induced venous CBV (CBV_v) change (see additional assumptions in Kim et al, 1999). This assumption may be invalid (see later), and if so, then reported literature values for quantitative functional CMRO₂ changes obtained from BOLD fMRI data may contain significant errors.

To easily conceptualize equation (3), a first-order linear approximation can be performed (Kim et al, 1999) as

$$\begin{array}{rl}\mathrm{\%}\mathrm{B}\mathrm{O}\mathrm{L}\mathrm{D}& =M(\frac{\mathrm{\Delta }\mathrm{C}\mathrm{B}\mathrm{F}/\mathrm{C}\mathrm{B}\mathrm{F}-\mathrm{\Delta }{\mathrm{C}\mathrm{M}\mathrm{R}\mathrm{O}}_2/{\mathrm{C}\mathrm{M}\mathrm{R}\mathrm{O}}_2}{\mathrm{\Delta }\mathrm{C}\mathrm{B}\mathrm{F}/\mathrm{C}\mathrm{B}\mathrm{F}+1}-\frac{1}{\beta }\cdot \frac{\mathrm{\Delta }{\mathrm{C}\mathrm{B}\mathrm{V}}_v}{{\mathrm{C}\mathrm{B}\mathrm{V}}_v})\\ & =M(\frac{\mathrm{\Delta }Y}{1-Y}-\frac{1}{\beta }\cdot \frac{\mathrm{\Delta }{\mathrm{C}\mathrm{B}\mathrm{V}}_v}{{\mathrm{C}\mathrm{B}\mathrm{V}}_v})\end{array}$$

(4)

This approximation shows that an increase in Y (a decrease in OEF) will increase fMRI signal, while an increase in CBV_v decreases fMRI signal. The term ΔY/(1−Y) is larger than the term ΔCBV_v/(β·CBV_v), and when venous CBV changes are negligible in equation (4), the BOLD signal change is simply related to the venous oxygenation level change weighted by baseline venous blood volume.

Since CBF, CMRO₂, and CBV_v responses contribute to BOLD fMRI dynamics (see equation (4)), these physiological contributions can be explained with biophysical models (e.g., Ogawa et al, 1993b; Boxerman et al, 1995; Ogawa et al, 1998). Quantitative contribution of physiological changes to BOLD fMRI responses requires knowledge of constants (e.g., M and α in equation (2)) not easily determined. Therefore, we describe only the properties of different physiological parameters and their impact on dynamic BOLD fMRI signals.

Dynamic Cerebral Blood Flow Responses

The CBF changes can be measured by perfusion-based fMRI techniques using arterial blood water as an endogenous flow tracer. Numerous CBF measurements have been performed with various stimuli (e.g., Kwong et al, 1992; Edelman et al, 1994; Kim, 1995). Perfusion-based fMRI methods have poor temporal resolution due to the requirement that spins from arterial blood outside the imaging slice must have sufficient time to move into capillaries and exchange with extravascular tissue water, and because control images must also be acquired. Usually, stimulus onset causes CBF to initially rise to an ‘overshoot’ level, after which the increase in CBF decreases to a slightly lower ‘plateau’ level during stimulus continuance. These increased CBF levels will increase venous oxygen saturation, which is directly related to a positive BOLD response.

Dynamic CBF responses at high temporal resolution can be obtained with a conventional laser Doppler flowmeter or from fMRI (e.g., Silva and Kim, 1999). Since blood flow dictates that arterial blood responses precede those of venous vessels, the CBF response is ∼0.5 to 1 second earlier than the BOLD response. To show the difference in hemodynamic onset times, CBF and BOLD responses in the rat somatosensory cortex to forepaw stimulus are plotted in Figure 5A (Silva et al, 2000); the CBF increase occurs ∼0.5 second after stimulus onset, while the BOLD response occurs ∼0.5 second later.

Figure 5

Dynamic changes in stimulus-induced CBF, BOLD, CBV, and CMRO₂-related fMRI signals. Stimulus starts at time=0 second and stimulus periods are indicated by black horizontal bars. (A) Dynamic CBF and BOLD responses to forepaw stimulus were obtained in the rat somatosensory cortex (Silva et al, 2000). These time courses for regions of interest from the cortical ‘surface’ and the middle of the cortex (‘deep’) show that the CBF response is earlier than the BOLD response. (B) Dynamic BOLD, CBV, and CMRO₂-related responses were obtained in cat visual cortex (Nagaoka et al, 2006). BOLD and CBV data were obtained at normal physiological conditions, while the CMRO₂-related change was acquired with infusion of a vasodilator to suppress stimulus-induced CBF and CBV changes. The CBV response was obtained after intravascular injection of iron oxide nanoparticles. Peak intensity from the time course of each response type was used for normalization. The BOLD response starts slightly later than CBV_t and CMRO₂ responses, but BOLD and CBV_t peak times are similar during 10 seconds stimulation, with a lag in CMRO₂ peak time. Note that a small BOLD dip was observed. (C) The amplitude and dynamics of BOLD responses in human visual cortex (Cohen et al, 2002) are closely dependent on baseline CBF levels. When baseline CBF is low due to hypocapnia, the BOLD response is of higher amplitude and faster versus normocapnia. Conversely, when baseline CBF is high due to hypercapnia, the BOLD response is of lower amplitude and sluggish versus normocapnia. BOLD, blood oxygenation level dependent; fMRI, functional magnetic resonance imaging; CBV, cerebral blood volume; CBF, cerebral blood flow; CMRO₂, cerebral metabolic rate of oxygen utilization; CBV_t, total CBV.

Dynamic Cerebral Metabolic Rate of Oxygen Utilization-Related Changes

It is well accepted that the CBF increase with stimulus more than compensates for the functional CMRO₂ increase, and it was also generally assumed that CMRO₂ increases rapidly and peaks quickly after stimulus onset. Most in-vivo assessments of dynamic changes in CMRO₂ rely on blood oxygenation measurements. The aforementioned calibration studies have suggested that functional increases in CMRO₂ are concomitant with increases in CBF (e.g., Davis et al, 1998). This observation is largely due to the overwhelming contribution of stimulus-induced CBF increases to blood oxygenation changes (Chiarelli et al, 2007). In addition, the calibration method assumes that the oxygen extracted from blood is instantaneously consumed by tissue cells, and the observed CMRO₂ is same as the actual cellular response.

The above confounds are removed when a vasodilator is infused to suppress CBF and CBV functional responses (Fukuda et al, 2006; Nagaoka et al, 2006), such that stimulus-induced dynamic changes in blood oxygenation levels measured by fMRI indicate solely CMRO₂ responses. This approach shows a functional CMRO₂ response that begins earlier, but surprisingly peaks later than the BOLD response (Nagaoka et al, 2006) (see blue versus red time courses in Figure 5B). A similar slow functional CMRO₂ response compared with CBF response was also observed in tissue oxygen tension data (Masamoto et al, 2008). The delayed peak in functional CMRO₂ response (versus BOLD and CBF) may be explained by the astrocyte-neuron-lactate-shuttle hypothesis (reviews in Pellerin et al, 2007), in which lactate is produced by increased glycolysis in astrocytes and then becomes the primary substrate for oxidative metabolism in neurons.

An alternative reason for observing delayed peaks in CMRO₂ responses by fMRI may involve temporal latencies between mitochondrial oxygen consumption and changes in venous blood oxygenation levels. This issue was recently examined by measuring dynamic cellular CMRO₂ changes with flavoprotein autofluorescence imaging (Vazquez et al, 2012), where it was found that the tissue mitochondrial CMRO₂ responds relatively fast (time constant ∼1 second) compared with the venous oxygenation level change (time constant ∼10 seconds), showing that dynamic CMRO₂ responses obtained from BOLD fMRI do not instantaneously reflect temporal changes in tissue mitochondrial metabolic rates.

Dynamic Cerebral Blood Volume Responses

Dynamic CBF and CBV changes are intercorrelated, since CBF is dependent on CBV and velocity changes. The total CBV (CBV_t) change can be estimated from the CBF change (Grubb et al, 1974), or can be measured with intravascular injection of contrast agent (Kennan et al, 1998; Mandeville et al, 1998; van Bruggen et al, 1998) (see Figures 4C and 5B) or by the vascular space occupancy method (Lu et al, 2003). In some early studies, the BOLD contribution to CBV-weighted fMRI was not corrected, leading to errors in quantifying dynamic CBV changes. More recently, it has been consistently observed that the CBV_t response occurs earlier than the BOLD response (see green versus red time courses in Figure 5B).

Total CBV can be subdivided into arterial and venous components. Since ∼60% to 80% of baseline CBV_t is from CBV_v (Ito et al, 2001; Lee et al, 2001; An and Lin, 2002b), it is often assumed that stimulus-induced CBV_t changes are dominated by CBV_v changes (Buxton et al, 1998; Mandeville et al, 1999). This assumption has been used to determine changes in venous oxygenation level and CMRO₂ from BOLD fMRI data. But determination of CMRO₂ from CBF and BOLD fMRI data depends on understanding venous—not total—CBV changes (see equation (4)).

Compartment-specific CBV studies indicate that the CBV_v change is actually much less than the CBV_t change. Within a short stimulus period (<15 seconds), rat fMRI studies showed that arterial CBV (CBV_a) changes are dominant, while CBV_v changes are minimal (Kim et al, 2007). During a long stimulus period, CBV_v increases slowly in cat studies (Figure 6A), eventually reaching a magnitude similar to the CBV_a change (Kim and Kim, 2011). At a steady state, the relative change in CBV_v is about half that of CBV_t in rat studies (Lee et al, 2001). Similar observations were also detected in optical studies in rodents, where arterial dilation was dominant with short stimulus duration (Hillman et al, 2007; Vazquez et al, 2010; Drew et al, 2011; see also Figure 6D), and venous dilation occurred slowly with extended stimulus (Drew et al, 2011). After stimulus offset, cat fMRI studies show that CBV_a rapidly decreases to baseline values or exhibits a small prolonged poststimulus undershoot, while CBV_v is slow to return to baseline values (Figure 6A). Although low sensitivity and large signal fluctuations indeed complicate identification of poststimulus CBV_a undershoots, it should be noted that similar poststimulus undershoots were observed in CBF fMRI data from the same animal model (see Figure 1D in Jin and Kim, 2008) and in surface arteriolar vessel diameters measured by two-photon microscopy in mice after termination of a 30-second vibrissae stimulus (see Figure 2 in Drew et al, 2011). In humans, positron emission tomographic studies showed that the CBV_a change is dominant during hypocapnia and hypercapnia, while the CBV_v change is minimal (Ito et al, 2005). Venous CBV can be noninvasively determined with the MRI technique named VERVE (venous refocusing for volume estimation) (Stefanovic and Pike, 2005), which relies on multispin-echo measurements with different interpulse delays. Recent fMRI measurements with VERVE in humans show that functional CBV_v responses are significant during stimulation and the power exponent, α, is 0.23 during neural stimulus (Chen and Pike, 2009a) and 0.18 during steady-state global hypocapnia and hypercapnia challenges (Chen and Pike, 2010).

Figure 6

Compartment-specific CBV changes during sensory stimulus. (A) Corresponding total, arterial, and venous CBV changes were measured by fMRI in the cat visual cortex (Kim and Kim, 2011). The arterial CBV change was obtained by fMRI using the varied magnetization transfer effect (Kim et al, 2008), while total CBV was determined with intravascular infusion of iron oxide nanoparticles. The venous CBV change was then estimated by the difference in total CBV and arterial CBV changes. The arterial CBV response is rapid like the total CBV response, while the venous CBV change is relatively slow. (B–D) Imaging of the rat somatosensory cortex followed isoflurane anesthesia and craniotomy (provided by Alberto Vazquez at the University of Pittsburgh). Intrinsic optical image with signal weighted to red blood cells shows pial vessels as hypointense structures (B). Baseline fluorescent microscopic images (C) enhanced by intravascular injection of Texas red dye show arterial and venous vessels (indicated by ‘A’ and ‘V’, respectively) as hyperintense structures. In the functional change map (D) determined by the difference between baseline images and images acquired during 4 seconds of forepaw stimulus, the highest signal intensities indicate stimulus-induced increases in dye content due to dilation. Only arterial vessels dilate, while venous vessels do not. fMRI, functional magnetic resonance imaging; CBV, cerebral blood volume.

The recent CBV_v studies have three important implications for BOLD quantification: (1) since during stimulus the relative CBV_v change is less than relative CBV_t change, relative CMRO₂ changes are underestimated when obtained with equation (3) and an α value of 0.38, as in most fMRI calibration studies (e.g., Davis et al, 1998), even at a steady state; (2) a decrease in BOLD signals with constant CBF during the stimulus period (Frahm et al, 1996) can be explained by a slow increase in CBV_v; and (3) the time-independent BOLD linearity may not hold due to the time-dependent CBV_v contribution when short versus long stimulus data are compared. If functional CBV_v response is not experimentally measured, then the CBV_v response determined from the relationship between CBF and CBV_v should be used for the determination of ΔY and ΔCMRO₂ from BOLD data. One simple approach for determining CMRO₂ change from the BOLD fMRI signal is to minimize the neural stimulus-induced CBV_v response with a relatively short stimulus (such as <15 seconds), and then to consider only venous oxygenation changes in the BOLD biophysical model (see equation (4)).

Dynamic Blood Oxygenation Level-Dependent Responses

The BOLD response function is determined by CBF, CMRO₂, and CBV_v changes. Since the BOLD response is sensitive to the draining effects from capillaries to pial veins, pixel-to-pixel differences in vascular content mean that even if stimulus-induced changes in CBF, CMRO₂, and CBV_v were identical across pixels, BOLD dynamic responses (onset time, time-to-peak, etc.) would still be expected to vary across pixels (Lee et al, 1995). This draining of blood into pial veins makes the BOLD response in the middle of the cortex start ∼0.5 second before the BOLD response at the cortical surface in rat somatosensory cortex (Silva et al, 2000; Tian et al, 2010), cat visual cortex (Jin and Kim, 2008), and human visual cortex (Siero et al, 2011). If instead fMRI signals from all pixels within an image originate mostly from microvessels including capillaries, then onset times should correspond to the capillary response times, which are likely to be similar across regions. If the microvessel contribution is dominant, then the difference in hemodynamic responses between regions may then be used to infer sequential neural activity (Buckner et al, 1996; Richter et al, 1997; Menon et al, 1998; Ogawa et al, 2000).

Sources of initial blood oxygenation level-dependent dips, poststimulus blood oxygenation level-dependent undershoots, and prolonged negative blood oxygenation level-dependent functional magnetic resonance imaging signals

Prolonged positive responses (where stimulus-induced signals are higher than baseline signals) are not the only type of response observed with conventional BOLD fMRI. Understanding signal sources of negative BOLD fMRI signals (signals lower than baseline), such as initial dips, poststimulus undershoots, and prolonged negative signals is of great interest.

Initial Blood Oxygenation Level-Dependent Dips

Optical spectroscopy studies by Malonek and Grinvald (1996) showed a delayed hemodynamic response that was diffuse with large spatial extent and also an initial CMRO₂-based deoxyhemoglobin signal increase (termed the ‘optical dip’) that was highly localized to the site of neuronal activity. Although the exact origin of the deoxyhemoglobin-weighted optical dip is still controversial (Malonek and Grinvald, 1996; Sirotin et al, 2009), it has been commonly observed. Similarly to the optical dip, a small initial BOLD dip has been detected in humans (Ernst and Hennig, 1994; Menon et al, 1995; Hu et al, 1997) and in animals (Kim et al, 2000) (see also red BOLD time course in Figure 5B), though existence of this response is still debated (Buxton, 2001). The initial BOLD dip may reflect (1) an increase in deoxyhemoglobin concentration due to CMRO₂ increases before the CBF response and (2) an increase in CBV_v preceding the CBF response. Recent measurements of vessel diameters (Drew et al, 2011) and CBV fMRI (Kim and Kim, 2011) during stimulus show a slow increase in CBV_v (see Figure 6A), with no evidence of significant venous vessel dilation during the initial period, suggesting that the initial BOLD dip is due to an early mismatch between CMRO₂ and CBF changes. However, detection of a highly localized initial BOLD dip requires both high spatial resolution and high temporal resolution, and the resulting reduction in sensitivity may not be practical for high-resolution fMRI studies.

Poststimulus Blood Oxygenation Level-Dependent Undershoots

One or more physiological factors may contribute to a poststimulus BOLD undershoot, including a poststimulus CBF undershoot (Hoge et al, 1999a; Jin and Kim, 2008; Chen and Pike, 2009b), a slow poststimulus return of CBV_v to baseline (Mandeville et al, 1998; Chen and Pike, 2009b; Hua et al, 2011; Kim and Kim, 2011), and a slow poststimulus return of CMRO₂-related oxygenation values to baseline (Frahm et al, 1996; Yacoub et al, 2006; Hua et al, 2011). One should be careful when interpreting early findings, because (1) CBV_t responses are often measured instead of CBV_v, (2) at high magnetic field studies, both CBF and CBV responses often contaminate the BOLD signal, and (3) there is no direct means of dynamically measuring CMRO₂. In recent BOLD, CBF, and CBV_v fMRI studies to explain the poststimulus BOLD fMRI undershoot, one report (Hua et al, 2011) suggested that ∼80% of the BOLD fMRI undershoot is from a slow return of CMRO₂-related contributions, with the remainder related to a slow return of CBV_v contributions, and no CBF undershoot contributions, while another paper (Chen and Pike, 2009b) reported that both a CBF undershoot and a slow return of CBV_v to baseline are contributors. Exact contributors of the poststimulus undershoot may also depend on stimulus type and duration.

Even if CMRO₂-related factors are dominant contributors to the poststimulus BOLD undershoot, this is still not an indication of an increase in cellular CMRO₂ during the poststimulus period, as discussed previously in ‘Dynamic CMRO₂-related changes’. More importantly, is the poststimulus BOLD undershoot useful for high-resolution fMRI? The answer is the same as for initial BOLD dips: although any CMRO₂-related BOLD signal would improve spatial specificity to areas of neuronal activity, both high spatial resolution and high temporal resolution would be required to localize deoxyhemoglobin in capillaries before it drains into large vessels, with the resulting loss of sensitivity. Therefore, the usefulness of the poststimulus BOLD undershoot for high-resolution functional imaging studies is dependent on its signal sources.

Prolonged Negative Blood Oxygenation Level-Dependent Functional Magnetic Resonance Imaging Signals

Negative BOLD signals occurring during a stimulus period may be explained by (1) a decrease in CBF due to neural inhibition (Shmuel et al, 2002, 2006; Boorman et al, 2010; Wade and Rowland, 2010), (2) a decrease in CBF due to redistribution of CBF (the steal effect) into nearby active regions (Harel et al, 2002), (3) an increase in CMRO₂ without concomitant CBF increase (Schridde et al, 2008), and (4) an increase in vasoconstrictive neurotransmitter with increased neural activity in subcortical areas (Shih et al, 2009). In areas where vascular reactivity is hampered, a negative BOLD response is likely due to increased CMRO₂ without a concomitant increase in CBF. The physiological basis of negative BOLD signals therefore depends on stimulus type and brain region (e.g., cortical versus subcortical areas).

Influence of global hemodynamic baseline conditions on stimulus-evoked blood oxygenation level-dependent responses

It is well known that global CBF can be modulated by changes in arterial partial pressures of carbon dioxide and oxygen. It has also been shown or suggested that global perturbation of CBF can result from the intake of commonly used substances (e.g., caffeine, nicotine, and alcohol), changes in the concentration of endogenous substances (e.g., estrogen and adrenaline), or experimental administration of various drugs (e.g., cocaine or acetazolamide). These global CBF changes can affect the dynamics and magnitude of BOLD responses (Hoge et al, 1999a; Kemna and Posse, 2001; Cohen et al, 2002). Typically, when baseline CBF (and consequently baseline venous oxygenation) is low, the BOLD response is fast and intense (Cohen et al, 2002; Lu et al, 2008). For example, a faster and larger stimulus-evoked BOLD response occurs under hypocapnic conditions (black trace in Figure 5C) or with intake of the vasoconstrictor caffeine (Mulderink et al, 2002), while a slower and smaller response occurs with alcohol intake (Levin et al, 1998). Baseline-condition dependence of BOLD responses may be due to (1) a nonlinear relationships between BOLD responses versus baseline-dependent values of Y or (2) baseline-condition dependence of CBF responses. At higher baseline venous oxygenation levels, the maximal allowable value for ΔY decreases; and as ΔY approaches zero, the BOLD response no longer linearly correlates with CBF changes, especially for extremely large CBF responses (Lee et al, 2002). However, within the range of normal physiological conditions, it is most likely that baseline-condition dependence of BOLD responses is due to baseline-condition dependence of CBF responses.

Additionally, the hemodynamic response is sensitive to age (Yamada et al, 2000; Desposito et al, 2003), patterns of resting-state fluctuations, and resting-state neurotransmitter levels (Muthukumaraswamy et al, 2009). When the stimulus-induced fMRI response is compared under different physiological or pharmacological states, investigators must distinguish between effects resulting from altered neural activity and those resulting merely from global hemodynamic influences.

Basis of resting-state blood oxygenation level-dependent fluctuations

A seminal report by Biswal et al (1995) showed that in the resting-state, slow temporal fluctuations (<0.1 Hz) in BOLD responses were correlated between functionally related areas, as defined by conventional stimulus-evoked fMRI, suggesting that large-scale neural networks could be mapped without task performance. Despite some concern that some contributions to these fluctuations could be due to vasomotor oscillations (Mitra et al, 1997), there is sufficient evidence of neural contributions to resting-state BOLD fluctuations. The field of resting-state fMRI expanded greatly after large-scale ‘default-mode networks’ was discovered in humans (Raichle et al, 2001) and in anesthetized monkeys (Vincent et al, 2007). Resting-state BOLD fMRI has become an important tool for empirically mapping various interhemispheric and intrahemispheric local and global networks and assessing connectivity in normal and diseased brains. Still, the relationship between resting-state BOLD fMRI and neural activity or neurovascular coupling is an area of intense research (see review article, Fox and Raichle, 2007). In monkey studies, Leopold and Logothetis (Leopold et al, 2003) measured fluctuations in band-limited power of LFP signals and found high coherence over cortical distances >1 cm, indicating that resting-state BOLD fluctuations reflect ongoing (spontaneous and coherent) neural activity. A question that needs to be answered by systematic investigation is whether the sources of resting-state BOLD responses are similar to those of stimulus-induced responses.

Concluding remarks

The BOLD effect in fMRI is very complex, and thus is still an area of intense research. However, the following conclusions can be drawn from our experiences and review of the existing literature.

Contrast in BOLD fMRI is determined by both intravascular and extravascular signal changes. With higher magnetic fields and longer TEs, the extravascular contribution increases and the intravascular contribution decreases; and therefore, the spatial specificity to parenchyma improves.

The magnitude of a BOLD fMRI response is weighted by baseline CBV_v and stimulus-induced venous oxygenation level changes. This means that if one pixel has a higher value for baseline CBV_v than another pixel, that pixel would show a higher BOLD response if all other factors were equal (including oxygenation level changes). Since vascular density and structure are unlikely to be similar across pixels in extremely high-resolution images, BOLD fMRI response magnitudes should therefore not compared across pixels.

Spatial resolution in BOLD fMRI is determined by the size of each venous vascular unit (∼1 mm diameter when signals from intracortical emerging veins contribute). Spatial resolution is also dependent on the area of activation, and can be improved by acquisition with SE sequences at high magnetic fields.

Steady-state CMRO₂ changes are often determined from BOLD and CBF data by estimating the CBV_t change from the CBF response with Grubb's formula and a volume-flow power term of 0.38. However, it is the CBV_v response (not CBV_t) that is relevant for BOLD fMRI, and at steady state, the CBV_v responses are only about half of CBV_t changes; therefore, functional CMRO₂ changes reported in the literature are significantly underestimated when calculated from CBV_t.

Functional CBF, CBV_v, and CMRO₂ responses contribute to the dynamics of BOLD fMRI responses. The CBV_v response is sluggish, and negligible for short stimulus periods. This lack of significant CBV_v contributions to short stimulus-induced BOLD responses simplifies the BOLD biophysical model.

Dynamic CMRO₂ changes obtained from the BOLD biophysical model do not reflect dynamic cellular CMRO₂ changes in mitochondria; and therefore, dynamic CMRO₂ values obtained from BOLD data should not be interpreted as dynamic cellular responses. This pertains, for example, to interpretation of CMRO₂-related contributions to the poststimulus BOLD undershoot.

The initial BOLD dip originates from an early increase in CMRO₂, but the source of poststimulus BOLD undershoot is still controversial. Investigating the sources of these signals will definitely aid our understanding of neurovascular coupling, but both high spatial and temporal resolutions would be required to use these signals for fMRI, and thus sensitivity would likely be low.

Dynamic properties and magnitudes of BOLD functional responses are dependent on many physiological parameters as well as baseline conditions. In patients with neurovascular disorders, the BOLD response could be sluggish, or even decreased relative to baseline. This should not be interpreted simply as a decrease in neural activity, because neurovascular coupling may be hampered.

Resting-state fMRI studies are widely performed, but its physiological source needs to be systematically investigated.