Jason Berwick

Concurrent optical imaging spectroscopy and laser-Doppler flowmetry: the relationship between blood flow, oxygenation, and volume in rodent barrel cortex.

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

Investigating neural-hemodynamic coupling and the hemodynamic response function in the awake rat.

An understanding of the relationship between changes in neural activity and the accompanying hemodynamic response is crucial for accurate interpretation of functional brain imaging data and in particular the blood oxygen level-dependent (BOLD) fMRI signal. Much physiological research investigating this topic uses anesthetized animal preparations, and yet, the effects of anesthesia upon the neural and hemodynamic responses measured in such studies are not well understood. In this study, we electrically stimulated the whisker pad of both awake and urethane anesthetized rats at frequencies of 1-40 Hz. Evoked field potential responses were recorded using electrodes implanted into the contralateral barrel cortex. Changes in hemoglobin oxygenation and concentration were measured using optical imaging spectroscopy, and cerebral blood flow changes were measured using laser Doppler flowmetry. A linear neural-hemodynamic coupling relationship was found in the awake but not the anesthetized animal preparation. Over the range of stimulation conditions studied, hemodynamic response magnitude increased monotonically with summed neural activity in awake, but not in anesthetized, animals. Additionally, the temporal structure of the hemodynamic response function was different in awake compared to anesthetized animals. The responses in each case were well approximated by gamma variates, but these were different in terms of mean latency (approximately 2 s awake; 4 s anesthetized) and width (approximately 0.6 s awake; 2.5 s anesthetized). These findings have important implications for research into the intrinsic signals that underpin BOLD fMRI and for biophysical models of cortical hemodynamics and neural-hemodynamic coupling.

Spectroscopic analysis of neural activity in brain: increased oxygen consumption following activation of barrel cortex.

This research investigates the hemodynamic response to stimulation of the barrel cortex in anaesthetized rats using optical imaging and spectroscopy (Bonhoeffer and Grinvald, 1996; Malonek and Grinvald, 1996; Mayhew et al., 1999). A slit spectrograph was used to collect spectral image data sequences. These were analyzed using an algorithm that corrects for the wavelength dependency in the optical path lengths produced by the light scattering properties of tissue. The analysis produced the changes in the oxy- and deoxygenation of hemoglobin following stimulation. Two methods of stimulation were used. One method mechanically vibrated a single whisker, the other electrically stimulated the whisker pad. The electrical stimulation intensity varied from 0.4 to 1.6 mA. The hemodynamic responses to stimulation increased as a function of intensity. At 0.4 mA they were commensurate with those from the mechanical stimulation; however, the responses at the higher levels were greater by a factor of approximately 10. For both methods of data collection, the results of the spectroscopic analysis showed an early increase in deoxygenated hemoglobin (Hbr) with no evidence for a corresponding decrease in oxygenated hemoglobin (HbO(2)). Evidence for increased oxygen consumption (CMRO(2)) was obtained by converting the fractional changes in blood volume (Hbt) into estimates of changes in blood flow (Grubb et al., 1974) and using the resulting time course to scale the fractional changes in Hbr. The results show an early increase CMRO(2) peaking approximately 2 s after stimulation onset. Using these methods, we find evidence for increased oxygen consumption following increased neural activity even at low levels of stimulation intensity.

Spectroscopic analysis of changes in remitted illumination: the response to increased neural activity in brain.

Imaging of neural activation has been used to produce maps of functional architecture and metabolic activity. There is some uncertainty associated with the sources underlying the intrinsic signals. It has been reported that following increased neural activity there was little increased oxygen consumption ( approximately 5%), although glucose consumption increased by approximately 50%. The research we describe uses a modification of the Beer-Lambert Law called path-length scaling analysis (PLSA) to analyze the spectra of the hemodynamic and metabolic responses to vibrissal stimulation in rat somatosensory cortex. The results of the PLSA algorithm were compared with those obtained using a linear spectrographic analysis method (we refer to this as LMCA). There are differences in the results of the analysis depending on which of the two algorithms (PLSA or LMCA) is used. Using the LMCA algorithm, we obtain results showing an increase in the volume of Hbr at approximately 2 s, following onset of stimulation but no complementary decrease in oxygenated haemoglobin (HbO(2)). These results are similar to a previous report. In contrast, after using the PLSA algorithm, the time series of the chromophore changes shows no evidence for an increase in the volume of deoxygenated haemoglobin (Hbr). However, after further analysis of the time series from the PLSA using general linear models (GLM) to remove contributions from low frequency baseline oscillations, both the HbO(2) and Hbr times series of the response to stimulation were found to be biphasic with an early decrease in saturation peaking approximately 1 s after onset of stimulation followed by a larger increase in saturation peaking at approximately 3 s. Finally, following the PLSA-then-GLM analysis procedure, we do not find convincing evidence for an increase in cytochrome oxidation following stimulation, though we demonstrate the PLSA algorithm to be capable of disassociating changes in cytochrome oxidation state from changes in hemoglobin oxygenation.

The hemodynamic impulse response to a single neural event.

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

Negative Blood Oxygen Level Dependence in the Rat:A Model for Investigating the Role of Suppression in Neurovascular Coupling

Modern neuroimaging techniques rely on neurovascular coupling to show regions of increased brain activation. However, little is known of the neurovascular coupling relationships that exist for inhibitory signals. To address this issue directly we developed a preparation to investigate the signal sources of one of these proposed inhibitory neurovascular signals, the negative blood oxygen level-dependent (BOLD) response (NBR), in rat somatosensory cortex. We found a reliable NBR measured in rat somatosensory cortex in response to unilateral electrical whisker stimulation, which was located in deeper cortical layers relative to the positive BOLD response. Separate optical measurements (two-dimensional optical imaging spectroscopy and laser Doppler flowmetry) revealed that the NBR was a result of decreased blood volume and flow and increased levels of deoxyhemoglobin. Neural activity in the NBR region, measured by multichannel electrodes, varied considerably as a function of cortical depth. There was a decrease in neuronal activity in deep cortical laminae. After cessation of whisker stimulation there was a large increase in neural activity above baseline. Both the decrease in neuronal activity and increase above baseline after stimulation cessation correlated well with the simultaneous measurement of blood flow suggesting that the NBR is related to decreases in neural activity in deep cortical layers. Interestingly, the magnitude of the neural decrease was largest in regions showing stimulus-evoked positive BOLD responses. Since a similar type of neural suppression in surround regions was associated with a negative BOLD signal, the increased levels of suppression in positive BOLD regions could importantly moderate the size of the observed BOLD response.

A three-compartment model of the hemodynamic response and oxygen delivery to brain

We describe a mathematical model linking changes in cerebral blood flow, blood volume and the blood oxygenation state in response to stimulation. The model has three compartments to take into account the fact that the cerebral blood flow and volume as measured concurrently using laser Doppler flowmetry and optical imaging spectroscopy have contributions from the arterial, capillary as well as the venous compartments of the vasculature. It is an extension to previous one-compartment hemodynamic models which assume that the measured blood volume changes are from the venous compartment only. An important assumption of the model is that the tissue oxygen concentration is a time varying state variable of the system and is driven by the changes in metabolic demand resulting from changes in neural activity. The model takes into account the pre-capillary oxygen diffusion by flexibly allowing the saturation of the arterial compartment to be less than unity. Simulations are used to explore the sensitivity of the model and to optimise the parameters for experimental data. We conclude that the three-compartment model was better than the one-compartment model at capturing the hemodynamics of the response to changes in neural activation following stimulation.

Neurovascular coupling investigated with two-dimensional optical imaging spectroscopy in rat whisker barrel cortex

Optical imaging slit spectroscopy is a powerful method for estimating quantitative changes in cerebral haemodynamics, such as deoxyhaemoglobin, oxyhaemoglobin and blood volume (Hbr, HbO2 and Hbt, respectively). Its disadvantage is that there is a large loss of spatial data as one image dimension is used to encode spectral wavelength information. Single wavelength optical imaging, on the other hand, produces high‐resolution spatiotemporal maps of brain activity, but yields only indirect measures of Hbr, HbO2 and Hbt. In this study we perform two‐dimensional optical imaging spectroscopy (2D‐OIS) in rat barrel cortex during contralateral whisker stimulation to obtain two‐dimensional maps over time of Hbr, HbO2 and Hbt. The 2D‐OIS was performed by illuminating the cortex with four wavelengths of light (575, 559, 495 and 587 nm), which were presented sequentially at a high frame rate (32 Hz). The contralateral whisker pad was stimulated using two different durations: 1 and 16 s (5 Hz, 1.2 mA). Control experiments used a hypercapnic (5% CO2) challenge to manipulate baseline blood flow and volume in the absence of corresponding neural activation. The 2D‐OIS method allowed separation of artery, vein and parenchyma regions. The magnitude of the haemodynamic response elicited varied considerably between different vascular compartments; the largest responses in Hbt were in the arteries and the smallest in the veins. Phase lags in the HbO2 response between arteries and veins suggest that a process of upstream signalling maybe responsible for dilating the arteries. There was also a consistent increase in Hbr from arterial regions after whisker stimulation.

Concurrent fMRI and optical measures for the investigation of the hemodynamic response function

Functional magnetic resonance imaging (fMRI) signal variations are based on a combination of changes in cerebral blood flow (CBF) and volume (CBV), and blood oxygenation. We investigated the relationship between these hemodynamic parameters in the rodent barrel cortex by performing fMRI concurrently with laser Doppler flowmetry (LDF) or optical imaging spectroscopy (OIS), following whisker stimulation and hypercapnic challenge. A difference between the positions of the maximum blood oxygenation level‐dependent (BOLD) and CBV changes was observed in coronal fMRI maps, with the BOLD region being more superficial. A 6.5% baseline blood volume fraction in this superficial region dropped to 4% in deeper cortical layers (corresponding to total hemoglobin baseline volumes Hbt0 = 110 μM and 67 μM, respectively), as inferred from maps of ΔR  2* . Baseline volume profiles were used to parameterize the Monte Carlo simulations (MCS) to interpret the 2D OIS. From this it was found that the optical blood volume measurements (i.e., changes in total hemoglobin) equated with CBV‐MRI measurements when the MRI data were taken from superficial cortical layers. Optical measures of activation showed a good spatial overlap with fMRI measurements taken in the same plane (covering the right hemisphere surface). Changes in CBV and CBF followed the scaling relationship CBV = CBFα, with mean α = 0.38 ± 0.06. Magn Reson Med 54:354–365, 2005.

Further nonlinearities in neurovascular coupling in rodent barrel cortex

An essential prerequisite for the accurate interpretation of noninvasive functional brain imaging techniques, such as blood oxygen level dependent (BOLD) fMRI, is a thorough understanding of the coupling relationship between neural activity and the haemodynamic response. The current study investigates this relationship using rat barrel cortex as a model. Neural input was measured by applying current source density (CSD) analysis to multi-laminar field potentials to remove ambiguities regarding the origin of the signal inherent in single electrode recordings. Changes in cerebral blood flow (CBF) were recorded with a laser Doppler flowmetry probe. The magnitude of neural and CBF responses were modulated over a large range by altering both the intensity and frequency of electrical whisker pad stimulation. Consistent with previous findings [Devor, A., et al., 2003. Neuron 39, 353-359; Sheth, S.A., et al., 2004. Neuron 42, 347-355] a power law function well described the relationship between neural activity and haemodynamics. Despite the nonlinearity of the coupling over the whole data set, the relationship was very well approximated by a linear function over mid-range stimuli. Altering the frequency of stimulation at 1.2 mA shifted the neural activity and corresponding haemodynamic response along this linear region, reconciling recent reports of a nonlinear relationship [Devor, A., et al., 2003. Neuron 39, 353-359; Jones, M., et al., 2004. NeuroImage 22, 956-965; Sheth, S.A., et al., 2004. Neuron 42, 347-355] with previous work that found a linear coupling relationship when altering stimulation frequency [Martindale, J., et al., 2003. J. Cereb. Blood Flow Metab. 23, 546-555; Ngai, A.C., et al., 1999. Brain Res. 837, 221-228; Sheth, S., et al., 2003. NeuroImage 19, 884-894]. Using stimuli within this linear range in imaging studies would simplify the interpretation of findings.