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

Effects of isoflurane and urethane anesthetics on glutamate neurotransmission in rat brain using in vivo amperometry

BACKGROUND

Aspects of glutamate neurotransmission implicated in normal and pathological conditions are predominantly evaluated using in vivo recording paradigms in rats anesthetized with isoflurane or urethane. Urethane and isoflurane anesthesia influence glutamate neurotransmission through different mechanisms; however, real-time outcome measures of potassium chloride (KCl)-evoked glutamate overflow and glutamate clearance kinetics have not been compared within and between regions of the brain. In order to maintain rigor and reproducibility within the literature between the two most common methods of anesthetized in vivo recording of glutamate, we compared glutamate signaling as a function of anesthesia and brain region in the rat strain most used in neuroscience.

METHODS

In the following experiments, in vivo amperometric recordings of KCl-evoked glutamate overflow and glutamate clearance kinetics (uptake rate and T80) in the cortex, hippocampus, and thalamus were performed using glutamate-selective microelectrode arrays (MEAs) in young adult male, Sprague-Dawley rats anesthetized with either isoflurane or urethane.

RESULTS

Potassium chloride (KCl)-evoked glutamate overflow was similar under urethane and isoflurane anesthesia in all brain regions studied. Analysis of glutamate clearance determined that the uptake rate was significantly faster (53.2%, p < 0.05) within the thalamus under urethane compared to isoflurane, but no differences were measured in the cortex or hippocampus. Under urethane, glutamate clearance parameters were region-dependent, with significantly faster glutamate clearance in the thalamus compared to the cortex but not the hippocampus (p < 0.05). No region-dependent differences were measured for glutamate overflow using isoflurane.

CONCLUSIONS

These data support that amperometric recordings of KCl-evoked glutamate under isoflurane and urethane anesthesia result in similar and comparable data. However, certain parameters of glutamate clearance can vary based on choice of anesthesia and brain region. In these circumstances, special considerations are needed when comparing previous literature and planning future experiments.

Effects of isoflurane and urethane anesthetics on glutamate neurotransmission in rat brain using in vivo amperometry

Aspects of glutamate neurotransmission implicated in normal and pathological conditions are often evaluated using in vivo recording paradigms in rats anesthetized with isoflurane or urethane. Urethane and isoflurane anesthesia influence glutamate neurotransmission through different mechanisms; however real-time outcome measures of potassium chloride (KCl)-evoked glutamate overflow and glutamate clearance kinetics have not been compared within and between regions of the brain. In the following experiments, in vivo amperometric recordings of KCl-evoked glutamate overflow and glutamate clearance kinetics (uptake rate and T80) in the cortex, hippocampus and thalamus were performed using glutamate-selective microelectrode arrays (MEAs) in young adult male, Sprague-Dawley rats anesthetized with isoflurane or urethane. Potassium chloride (KCl)-evoked glutamate overflow was similar under urethane and isoflurane anesthesia in all brain regions studied. Analysis of glutamate clearance determined that the uptake rate was significantly faster (53.2%, p<0.05) within the thalamus under urethane compared to isoflurane, but no differences were measured in the cortex or hippocampus. Under urethane, glutamate clearance parameters were region dependent, with significantly faster glutamate clearance in the thalamus compared to the cortex but not the hippocampus (p<0.05). No region dependent differences were measured for glutamate overflow using isoflurane. These data support that amperometric recordings of glutamate under isoflurane and urethane anesthesia result in mostly similar and comparable data. However, certain parameters of glutamate uptake vary based on choice of anesthesia and brain region. Special considerations must be given to these areas when considering comparison to previous literature and when planning future experiments.

Sub-millimetre resolution laminar fMRI using Arterial Spin Labelling in humans at 7 T

Laminar fMRI at ultra-high magnetic field strength is typically carried out using the Blood Oxygenation Level-Dependent (BOLD) contrast. Despite its unrivalled sensitivity to detecting activation, the BOLD contrast is limited in its spatial specificity due to signals stemming from intra-cortical ascending and pial veins. Alternatively, regional changes in perfusion (i.e., cerebral blood flow through tissue) are colocalised to neuronal activation, which can be non-invasively measured using Arterial Spin Labelling (ASL) MRI. In addition, ASL provides a quantitative marker of neuronal activation in terms of perfusion signal, which is simultaneously acquired along with the BOLD signal. However, ASL for laminar imaging is challenging due to the lower SNR of the perfusion signal and higher RF power deposition i.e., specific absorption rate (SAR) of ASL sequences. In the present study, we present for the first time in humans, isotropic sub-millimetre spatial resolution functional perfusion images using Flow-sensitive Alternating Inversion Recovery (FAIR) ASL with a 3D-EPI readout at 7 T. We show that robust statistical activation maps can be obtained with perfusion-weighting in a single session. We observed the characteristic BOLD amplitude increase towards the superficial laminae, and, in apparent discrepancy, the relative perfusion profile shows a decrease of the amplitude and the absolute perfusion profile a much smaller increase towards the cortical surface. Considering the draining vein effect on the BOLD signal using model-based spatial “convolution”, we show that the empirically measured perfusion and BOLD profiles are, in fact, consistent with each other. This study demonstrates that laminar perfusion fMRI in humans is feasible at 7 T and that caution must be exercised when interpreting BOLD signal laminar profiles as direct representation of the cortical distribution of neuronal activity.

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.

The effects of capillary transit time heterogeneity on the BOLD signal

Neurovascular coupling mechanisms give rise to vasodilation and functional hyperemia upon neural activation, thereby altering blood oxygenation. This blood oxygenation level dependent (BOLD) contrast allows studies of activation patterns in the working human brain by functional MRI (fMRI). The BOLD-weighted fMRI signal shows characteristic transients in relation to functional activation, such as the so-called initial dip, overshoot, and post-stimulus undershoot. These transients are modulated by other physiological stimuli and in disease, but the underlying physiological mechanisms remain incompletely understood. Capillary transit time heterogeneity (CTH) has been shown to affect oxygen extraction, and hence blood oxygenation. Here, we examine how recently reported redistributions of capillary blood flow during functional activation would be expected to affect BOLD signal transients. We developed a three-compartment (hemoglobin, plasma, and tissue) model to predict the BOLD signal, incorporating the effects of dynamic changes in CTH. Our model predicts that the BOLD signal represents the superposition of a positive component resulting from increases in cerebral blood flow (CBF), and a negative component, resulting from elevated tissue metabolism and homogenization of capillary flows (reduced CTH). The model reproduces salient features of BOLD signal dynamics under conditions such as hypercapnia, hyperoxia, and caffeine intake, where both brain physiology and BOLD characteristics are altered. Neuroglial signaling and metabolism could affect CBF and capillary flow patterns differently. Further studies of neurovascular and neuro-capillary coupling mechanisms may help us relate BOLD signals to the firing of certain neuronal populations based on their respective BOLD "fingerprints."

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.

Laminar fMRI: What can the time domain tell us?

The rapid developments in functional MRI (fMRI) acquisition methods and hardware technologies in recent years, particularly at high field (≥7 T), have enabled unparalleled visualization of functional detail at a laminar or columnar level, bringing fMRI close to the intrinsic resolution of brain function. These advances highlight the potential of high resolution fMRI to be a valuable tool to study the fundamental processing performed in cortical micro-circuits, and their interactions such as feedforward and feedback processes. Notably, because fMRI measures neuronal activity via hemodynamics, the ultimate resolution it affords depends on the spatial specificity of hemodynamics to neuronal activity at a detailed spatial scale, and by the evolution of this specificity over time. Several laminar (≤1 mm spatial resolution) fMRI studies have examined spatial characteristics of the measured hemodynamic signals across cortical depth, in light of understanding or improving the spatial specificity of laminar fMRI. Few studies have examined temporal features of the hemodynamic response across cortical depth. Temporal features of the hemodynamic response offer an additional means to improve the specificity of fMRI, and could help target neuronal processes and neurovascular coupling relationships across laminae, for example by differences in the onset times of the response across cortical depth. In this review, we discuss factors that affect the timing of neuronal and hemodynamic responses across laminae, touching on the neuronal laminar organization, and focusing on the laminar vascular organization. We provide an overview of hemodynamics across the cortical vascular tree based on optical imaging studies, and review temporal aspects of hemodynamics that have been examined across cortical depth in high spatiotemporal resolution fMRI studies. Last, we discuss the limits and potential of high spatiotemporal resolution fMRI to study laminar neurovascular coupling and neuronal processes.

Foundations of layer-specific fMRI and investigations of neurophysiological activity in the laminarized neocortex and olfactory bulb of animal models

Laminar organization of neuronal circuits is a recurring feature of how the brain processes information. For instance, different layers compartmentalize different cell types, synaptic activities, and have unique intrinsic and extrinsic connections that serve as units for specialized signal processing. Functional MRI is an invaluable tool to investigate laminar processing in the in vivo human brain, but it measures neuronal activity indirectly by way of the hemodynamic response. Therefore, the accuracy of high-resolution laminar fMRI depends on how precisely it can measure localized microvascular changes nearest to the site of evoked activity. To determine the specificity of fMRI responses to the true neurophysiological responses across layers, the flexibility to invasive procedures in animal models has been necessary. In this review, we will examine different fMRI contrasts and their appropriate uses for layer-specific fMRI, and how localized laminar processing was examined in the neocortex and olfactory bulb. Through collective efforts, it was determined that microvessels, including capillaries, are regulated within single layers and that several endogenous and contrast-enhanced fMRI contrast mechanisms can separate these neural-specific vascular changes from the nonspecific, especially cerebral blood volume-weighted fMRI with intravenous contrast agent injection. We will also propose some open questions that are relevant for the successful implementation of layer-specific fMRI and its potential future directions to study laminar processing when combined with optogenetics.

Non-invasive quantification of absolute cerebral blood volume during functional activation applicable to the whole human brain

PURPOSE

Cerebral blood volume (CBV) changes in many diverse pathologic conditions, and in response to functional challenges along with changes in blood flow, blood oxygenation, and the cerebral metabolic rate of oxygen. The feasibility of a new method for non-invasive quantification of absolute cerebral blood volume that can be applicable to the whole human brain was investigated.

METHODS

Multi-slice data were acquired at 3 T using a novel inversion recovery echo planar imaging (IR-EPI) pulse sequence with varying contrast weightings and an efficient rotating slice acquisition order, at rest and during visual activation. A biophysical model was used to estimate absolute cerebral blood volume at rest and during activation, and oxygenation during activation, on data from 13 normal human subjects.

RESULTS

Cerebral blood volume increased by 21.7% from 6.6 ± 0.8 mL/100 mL of brain parenchyma at rest to 8.0 ± 1.3 mL/100 mL of brain parenchyma in the occipital cortex during visual activation, with average blood oxygenation of 84 ± 2.1% during activation, comparing well with literature.

CONCLUSION

The method is feasible, and could foster improved understanding of the fundamental physiological relationship between neuronal activity, hemodynamic changes, and metabolism underlying brain activation; complement existing methods for estimating compartmental changes; and potentially find utility in evaluating vascular health.

Temporal pattern of acoustic imaging noise asymmetrically modulates activation in the auditory cortex

This study investigated the hemisphere-specific effects of the temporal pattern of imaging related acoustic noise on auditory cortex activation. Hemodynamic responses (HDRs) to five temporal patterns of imaging noise corresponding to noise generated by unique combinations of imaging volume and effective repetition time (TR), were obtained using a stroboscopic event-related paradigm with extra-long (≥27.5 s) TR to minimize inter-acquisition effects. In addition to confirmation that fMRI responses in auditory cortex do not behave in a linear manner, temporal patterns of imaging noise were found to modulate both the shape and spatial extent of hemodynamic responses, with classically non-auditory areas exhibiting responses to longer duration noise conditions. Hemispheric analysis revealed the right primary auditory cortex to be more sensitive than the left to the presence of imaging related acoustic noise. Right primary auditory cortex responses were significantly larger during all the conditions. This asymmetry of response to imaging related acoustic noise could lead to different baseline activation levels during acquisition schemes using short TR, inducing an observed asymmetry in the responses to an intended acoustic stimulus through limitations of dynamic range, rather than due to differences in neuronal processing of the stimulus. These results emphasize the importance of accounting for the temporal pattern of the acoustic noise when comparing findings across different fMRI studies, especially those involving acoustic stimulation.

Optical Coherence Tomography angiography reveals laminar microvascular hemodynamics in the rat somatosensory cortex during activation

The BOLD (blood-oxygen-level dependent) fMRI (functional Magnetic Resonance Imaging) signal is shaped, in part, by changes in red blood cell (RBC) content and flow across vascular compartments over time. These complex dynamics have been challenging to characterize directly due to a lack of appropriate imaging modalities. In this study, making use of infrared light scattering from RBCs, depth-resolved Optical Coherence Tomography (OCT) angiography was applied to image laminar functional hyperemia in the rat somatosensory cortex. After defining and validating depth-specific metrics for changes in RBC content and speed, laminar hemodynamic responses in microvasculature up to cortical depths of >1mm were measured during a forepaw stimulus. The results provide a comprehensive picture of when and where changes in RBC content and speed occur during and immediately following cortical activation. In summary, the earliest and largest microvascular RBC content changes occurred in the middle cortical layers, while post-stimulus undershoots were most prominent superficially. These laminar variations in positive and negative responses paralleled known distributions of excitatory and inhibitory synapses, suggesting neuronal underpinnings. Additionally, the RBC speed response consistently returned to baseline more promptly than RBC content after the stimulus across cortical layers, supporting a "flow-volume mismatch" of hemodynamic origin.

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.

Noninvasive diffusive optical imaging of the auditory response to birdsong in the zebra finch

Songbirds communicate by learned vocalizations with concomitant changes in neurophysiological and genomic activities in discrete parts of the brain. Here, we tested a novel implementation of diffusive optical imaging (also known as diffuse optical imaging, DOI) for monitoring brain physiology associated with vocal signal perception. DOI noninvasively measures brain activity using red and near-infrared light delivered through optic fibers (optodes) resting on the scalp. DOI does not harm subjects, so it raises the possibility of repeatedly measuring brain activity and the effects of accumulated experience in the same subject over an entire life span, all while leaving tissue intact for further study. We developed a custom-made apparatus for interfacing optodes to the zebra finch (Taeniopygia guttata) head using 3D modeling software and rapid prototyping technology, and applied it to record responses to presentations of birdsong in isoflurane-anesthetized zebra finches. We discovered a subtle but significant difference between the hemoglobin spectra of zebra finches and mammals which has a major impact in how hemodynamic responses are interpreted in the zebra finch. Our measured responses to birdsong playback were robust, highly repeatable, and readily observed in single trials. Responses were complex in shape and closely paralleled responses described in mammals. They were localized to the caudal medial portion of the brain, consistent with response localization from prior gene expression, electrophysiological, and functional magnetic resonance imaging studies. These results define an approach for collecting neurophysiological data from songbirds that should be applicable to diverse species and adaptable for studies in awake behaving animals.

High-resolution fMRI reveals laminar differences in neurovascular coupling between positive and negative BOLD responses

The six cortical layers have distinct anatomical and physiological properties, like different energy use and different feedforward and feedback connectivity. It is not known if and how layer-specific neural processes are reflected in the fMRI signal. To address this question we used high-resolution fMRI to measure BOLD, CBV, and CBF responses to stimuli that elicit positive and negative BOLD signals in macaque primary visual cortex. We found that regions with positive BOLD responses had parallel increases in CBV and CBF, whereas areas with negative BOLD responses showed a decrease in CBF but an increase in CBV. For positive BOLD responses, CBF and CBV increased in the center of the cortex, but for negative BOLD responses, CBF decreased superficially while CBV increased in the center. Our findings suggest different mechanisms for neurovascular coupling for BOLD increases and decreases, as well as laminar differences in neurovascular coupling.

Cerebral blood volume changes during brain activation

Cerebral blood volume (CBV) changes significantly with brain activation, whether measured using positron emission tomography, functional magnetic resonance imaging (fMRI), or optical microscopy. If cerebral vessels are considered to be impermeable, the contents of the skull incompressible, and the skull itself inextensible, task- and hypercapnia-related changes of CBV could produce intolerable changes of intracranial pressure. Because it is becoming clear that CBV may be useful as a well-localized marker of neural activity changes, a resolution of this apparent paradox is needed. We have explored the idea that much of the change in CBV is facilitated by exchange of water between capillaries and surrounding tissue. To this end, we developed a novel hemodynamic boundary-value model and found approximate solutions using a numerical algorithm. We also constructed a macroscopic experimental model of a single capillary to provide biophysical insight. Both experiment and theory model capillary membranes as elastic and permeable. For a realistic change of input pressure, a relative pipe volume change of 21±5% was observed when using the experimental setup, compared with the value of approximately 17±1% when this quantity was calculated from the mathematical model. Volume, axial flow, and pressure changes are in the expected range.

Is optical imaging spectroscopy a viable measurement technique for the investigation of the negative BOLD phenomenon? A concurrent optical imaging spectroscopy and fMRI study at high field (7 T)

Traditionally functional magnetic resonance imaging (fMRI) has been used to map activity in the human brain by measuring increases in the Blood Oxygenation Level Dependent (BOLD) signal. Often accompanying positive BOLD fMRI signal changes are sustained negative signal changes. Previous studies investigating the neurovascular coupling mechanisms of the negative BOLD phenomenon have used concurrent 2D-optical imaging spectroscopy (2D-OIS) and electrophysiology (Boorman et al., 2010). These experiments suggested that the negative BOLD signal in response to whisker stimulation was a result of an increase in deoxy-haemoglobin and reduced multi-unit activity in the deep cortical layers. However, Boorman et al. (2010) did not measure the BOLD and haemodynamic response concurrently and so could not quantitatively compare either the spatial maps or the 2D-OIS and fMRI time series directly. Furthermore their study utilised a homogeneous tissue model in which is predominantly sensitive to haemodynamic changes in more superficial layers.

Here we test whether the 2D-OIS technique is appropriate for studies of negative BOLD. We used concurrent fMRI with 2D-OIS techniques for the investigation of the haemodynamics underlying the negative BOLD at 7 Tesla. We investigated whether optical methods could be used to accurately map and measure the negative BOLD phenomenon by using 2D-OIS haemodynamic data to derive predictions from a biophysical model of BOLD signal changes. We showed that despite the deep cortical origin of the negative BOLD response, if an appropriate heterogeneous tissue model is used in the spectroscopic analysis then 2D-OIS can be used to investigate the negative BOLD phenomenon.

MRI in animal models of psychiatric disorders

Here we describe MRI and (1)H MRS protocols for the investigation of animal models (mainly mice and rats) of psychiatric disorders. The introduction provides general findings from brain imaging studies in patients with psychiatric diseases and refers to general rules regarding the use of animal models in research. The methods section includes a selection of basic 9.4 T MRI and MRS protocols applicable for the investigation of animal models of psychiatric disorders (T1W, T2W, FLAIR, (1)H MRS). The notes section discusses in detail a series of factors that can influence the outcome of the experiment: from animal handling, stress-triggering aspects, and experimental design-related factors to technical aspects that affect T (1) and T (2) measurements.

Dynamic models of BOLD contrast

This personal recollection looks at the evolution of ideas about the dynamics of the blood oxygenation level dependent (BOLD) signal, with an emphasis on the balloon model. From the first detection of the BOLD response it has been clear that the signal exhibits interesting dynamics, such as a pronounced and long-lasting post-stimulus undershoot. The BOLD response, reflecting a change in local deoxyhemoglobin, is a combination of a hemodynamic response, related to changes in blood flow and venous blood volume, and a metabolic response related to oxygen metabolism. Modeling is potentially a way to understand the complex path from changes in neural activity to the BOLD signal. In the early days of fMRI it was hoped that the hemodynamic/metabolic response could be modeled in a unitary way, with blood flow, oxygen metabolism, and venous blood volume-the physiological factors that affect local deoxyhemoglobin-all tightly linked. The balloon model was an attempt to do this, based on the physiological ideas of limited oxygen delivery at baseline and a slow recovery of venous blood volume after the stimulus (the balloon effect), and this simple model of the physiology worked well to simulate the BOLD response. However, subsequent experiments suggest a more complicated picture of the underlying physiology, with blood flow and oxygen metabolism driven in parallel, possibly by different aspects of neural activity. In addition, it is still not clear whether the post-stimulus undershoot is a hemodynamic or a metabolic phenomenon, although the original venous balloon effect is unlikely to be the full explanation, and a flow undershoot is likely to be important. Although our understanding of the physics of the BOLD response is now reasonably solid, our understanding of the underlying physiological relationships is still relatively poor, and this is the primary hurdle for future models of BOLD dynamics.

fMRI of pain processing in the brain: a within-animal comparative study of BOLD vs. CBV and noxious electrical vs. noxious mechanical stimulation in rat

This study aims to identify fMRI signatures of nociceptive processing in whole brain of anesthetized rats during noxious electrical stimulation (NES) and noxious mechanical stimulation (NMS) of paw. Activation patterns for NES were mapped with blood oxygen level dependent (BOLD) and cerebral blood volume (CBV) fMRI, respectively, to investigate the spatially-dependent hemodynamic responses during nociception processing. A systematic evaluation of fMRI responses to varying frequencies of electrical stimulus was carried out to optimize the NES protocol. Both BOLD and CBV fMRI showed widespread activations, but with different spatial characteristics. While BOLD and CBV showed well-localized activations in ipsilateral dorsal column nucleus, contralateral primary somatosensory cortex (S1), and bilateral caudate putamen (CPu), CBV fMRI showed additional bilateral activations in the regions of pons, midbrain and thalamus compared to BOLD fMRI. CBV fMRI that offers higher sensitivity compared to BOLD was then used to compare the nociception processing during NES and NMS in the same animal. The activations in most regions were similar. In the medulla, however, NES induced a robust activation in the ipsilateral dorsal column nucleus while NMS showed no activation. This study demonstrates that (1) the hemodynamic response to nociception is spatial-dependent; (2) the widespread activations during nociception in CBV fMRI are similar to what have been observed in (14)C-2-deoxyglucose (2DG) autoradiography and PET; (3) the bilateral activations in the brain originate from the divergence of neural responses at supraspinal level; and (4) the similarity of activation patterns suggests that nociceptive processing in rats is similar during NES and NMS.

Physiological origin for the BOLD poststimulus undershoot in human brain: vascular compliance versus oxygen metabolism

The poststimulus blood oxygenation level-dependent (BOLD) undershoot has been attributed to two main plausible origins: delayed vascular compliance based on delayed cerebral blood volume (CBV) recovery and a sustained increased oxygen metabolism after stimulus cessation. To investigate these contributions, multimodal functional magnetic resonance imaging was employed to monitor responses of BOLD, cerebral blood flow (CBF), total CBV, and arterial CBV (CBV(a)) in human visual cortex after brief breath hold and visual stimulation. In visual experiments, after stimulus cessation, CBV(a) was restored to baseline in 7.9±3.4 seconds, and CBF and CBV in 14.8±5.0 seconds and 16.1±5.8 seconds, respectively, all significantly faster than BOLD signal recovery after undershoot (28.1±5.5 seconds). During the BOLD undershoot, postarterial CBV (CBV(pa), capillaries and venules) was slightly elevated (2.4±1.8%), and cerebral metabolic rate of oxygen (CMRO(2)) was above baseline (10.6±7.4%). Following breath hold, however, CBF, CBV, CBV(a) and BOLD signals all returned to baseline in ∼20 seconds. No significant BOLD undershoot, and residual CBV(pa) dilation were observed, and CMRO(2) did not substantially differ from baseline. These data suggest that both delayed CBV(pa) recovery and enduring increased oxidative metabolism impact the BOLD undershoot. Using a biophysical model, their relative contributions were estimated to be 19.7±15.9% and 78.7±18.6%, respectively.

Assessment of temporal state-dependent interactions between auditory fMRI responses to desired and undesired acoustic sources

A confounding factor in auditory functional magnetic resonance imaging (fMRI) experiments is the presence of the acoustic noise inherently associated with the echo planar imaging acquisition technique. Previous studies have demonstrated that this noise can induce unwanted neuronal responses that can mask stimulus-induced responses. Similarly, activation accumulated over multiple stimuli has been demonstrated to elevate the baseline, thus reducing the dynamic range available for subsequent responses. To best evaluate responses to auditory stimuli, it is necessary to account for the presence of all recent acoustic stimulation, beginning with an understanding of the attenuating effects brought about by interaction between and among induced unwanted neuronal responses, and responses to desired auditory stimuli. This study focuses on the characterization of the duration of this temporal memory and qualitative assessment of the associated response attenuation. Two experimental parameters--inter-stimulus interval (ISI) and repetition time (TR)--were varied during an fMRI experiment in which participants were asked to passively attend to an auditory stimulus. Results present evidence of a state-dependent interaction between induced responses. As expected, attenuating effects of these interactions become less significant as TR and ISI increase and in contrast to previous work, persist up to 18s after a stimulus presentation.

Spatially extended FMRI signal response to stimulus in non-functionally relevant regions of the human brain: preliminary results

The blood-oxygenation level dependent (BOLD) haemodynamic response function (HDR) in functional magnetic resonance imaging (fMRI) is a delayed and indirect marker of brain activity. In this single case study a small BOLD response synchronised with the stimulus paradigm is found globally, i.e. in all areas outside those of expected activation in a single subject study. The nature of the global response has similar shape properties to the archetypal BOLD HDR, with an early positive signal and a late negative response typical of the negative overshoot. Fitting Poisson curves to these responses showed that voxels were potentially split into two sets: one with dominantly positive signal and the other predominantly negative. A description, quantification and mapping of the global BOLD response is provided along with a 2 × 2 classification table test to demonstrate existence with very high statistical confidence. Potential explanations of the global response are proposed in terms of 1) global HDR balancing; 2) resting state network modulation; and 3) biological systems synchronised with the stimulus cycle. Whilst these widespread and low-level patterns seem unlikely to provide additional information for determining activation in functional neuroimaging studies as conceived in the last 15 years, knowledge of their properties may assist more comprehensive accounts of brain connectivity in the future.

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.

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.

Interpreting oxygenation-based neuroimaging signals: the importance and the challenge of understanding brain oxygen metabolism

Functional magnetic resonance imaging is widely used to map patterns of brain activation based on blood oxygenation level dependent (BOLD) signal changes associated with changes in neural activity. However, because oxygenation changes depend on the relative changes in cerebral blood flow (CBF) and cerebral metabolic rate of oxygen (CMRO(2)), a quantitative interpretation of BOLD signals, and also other functional neuroimaging signals related to blood or tissue oxygenation, is fundamentally limited until we better understand brain oxygen metabolism and how it is related to blood flow. However, the positive side of the complexity of oxygenation signals is that when combined with dynamic CBF measurements they potentially provide the best tool currently available for investigating the dynamics of CMRO(2). This review focuses on the problem of interpreting oxygenation-based signals, the challenges involved in measuring CMRO(2) in general, and what is needed to put oxygenation-based estimates of CMRO(2) on a firm foundation. The importance of developing a solid theoretical framework is emphasized, both as an essential tool for analyzing oxygenation-based multimodal measurements, and also potentially as a way to better understand the physiological phenomena themselves. The existing data, integrated within a simple theoretical framework of O(2) transport, suggests the hypothesis that an important functional role of the mismatch of CBF and CMRO(2) changes with neural activation is to prevent a fall of tissue pO(2). Future directions for better understanding brain oxygen metabolism are discussed.

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

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.

Tracking the effects of crusher gradients on gradient-echo BOLD signal in space and time during rat sensory stimulation

A unique method to map the effect of crusher gradients in space and time on the gradient echo blood oxygen level dependent (BOLD) signal is introduced. Using the Radial Correlation Contrast (RCC) analysis method, amplitude-RCC maps at different time segments and different gradient strengths were obtained. The ratio of amplitude-RCC cluster volumes, with and without crusher gradients, showed a temporal dependency with stronger volume reduction for stimulation-onset versus stimulation-decline. Aside from signal-to-noise ratio reduction in diffusion weighted images, the average temporal patterns were equal. Comparison of the data with and without crushers showed a stronger reduction in local coherence for stimulation-onset times. We hypothesize that the stimulation decline was weighted by extravascular effects originating in expanded veins due to their larger volume and long range susceptibility which couples neighboring voxels. The ratio of amplitude-RCC with and without crushers calculated for each voxel at each time segment yielded a spatial-temporal mapping of the crusher effect. These maps suggest that early stimulation-onset ( approximately 9 s) is weighted by flow; later a dynamic steady-state between intra- and extravascular effects is obtained. Stimulation-decline was dominated by extravascular effects, and at late stimulation decline as well as at early stimulation onset, clusters were small and localized to expected site of neuronal activity.

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.

Coupling between neuronal activity and microcirculation: implications for functional brain imaging

In the neocortex, neurons with similar response properties are often clustered together in column-like structures, giving rise to what has become known as functional architecture-the mapping of various stimulus feature dimensions onto the cortical sheet. At least partially, we owe this finding to the availability of several functional brain imaging techniques, both post-mortem and in-vivo, which have become available over the last two generations, revolutionizing neuroscience by yielding information about the spatial organization of active neurons in the brain. Here, we focus on how our understanding of such functional architecture is linked to the development of those functional imaging methodologies, especially to those that image neuronal activity indirectly, through metabolic or haemodynamic signals, rather than directly through measurement of electrical activity. Some of those approaches allow exploring functional architecture at higher spatial resolution than others. In particular, optical imaging of intrinsic signals reaches the striking detail of approximately 50 mum, and, together with other methodologies, it has allowed characterizing the metabolic and haemodynamic responses induced by sensory-evoked neuronal activity. Here, we review those findings about the spatio-temporal characteristics of neurovascular coupling and discuss their implications for functional brain imaging, including position emission tomography, and non-invasive neuroimaging techniques, such as funtional magnetic resonance imaging, applicable also to the human brain.

BOLD study of stimulation-induced neural activity and resting-state connectivity in medetomidine-sedated rat

Functional magnetic resonance imaging (fMRI) in anesthetized-animals is critical in studying the mechanisms of fMRI and investigating animal models of various diseases. Medetomidine was recently introduced for independent anesthesia for longitudinal (survival) fMRI studies in rats. Since stimulation-induced fMRI signal is anesthesia-dependent and its characteristics in rats under medetomidine are not fully elucidated, the blood oxygenation level dependent (BOLD) fMRI response to electrical forepaw stimulation under medetomidine was systematically investigated at 9.4 T. Robust activations in contralateral primary somatosensory cortex (SI) and thalamus were observed and peaked at the stimulus frequency of 9 Hz. The response in SI saturates at the stimulus strength of 4 mA while that in thalamus monotonically increases. In addition to fMRI data acquired with the forepaw stimulation, data were also acquired during the resting-state to investigate the synchronization of low frequency fluctuations (LFF) in the BOLD signal (<0.08 Hz) in different brain regions. LFF during resting-state have been observed to be synchronized between functionally related brain regions in human subjects while its origin is not fully understood. LFF have not been extensively studied or widely reported in anesthetized-animals. In our data, synchronized LFF of BOLD signals are found in clustered, bilaterally symmetric regions, including SI and caudate-putamen and the magnitude of the LFF is approximately 1.5%, comparable to the stimulation-induced BOLD signals. Similar to resting-state data reported in human subjects, LFF in rats under medetomidine likely reflect functional connectivity of these brain regions.