John Martindale

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.

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.

Nonlinear coupling of neural activity and CBF in rodent barrel cortex

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

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.

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

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

A Model of the Dynamic Relationship Between Blood Flow and Volume Changes During Brain Activation

The temporal relationship between changes in cerebral blood flow (CBF) and cerebral blood volume (CBV) is important in the biophysical modeling and interpretation of the hemodynamic response to activation, particularly in the context of magnetic resonance imaging and the blood oxygen level–dependent signal. Grubb et al. (1974) measured the steady state relationship between changes in CBV and CBF after hypercapnic challenge. The relationship CBVαCBFΦ has been used extensively in the literature. Two similar models, the Balloon (Buxton et al., 1998) and the Windkessel (Mandeville et al., 1999), have been proposed to describe the temporal dynamics of changes in CBV with respect to changes in CBF. In this study, a dynamic model extending the Windkessel model by incorporating delayed compliance is presented. The extended model is better able to capture the dynamics of CBV changes after changes in CBF, particularly in the return-to-baseline stages of the response.

Optical imaging spectroscopy in the unanaesthetised rat

We describe a method for imaging the local cortical haemodynamic response to whisker stimulation in the rat without use of anaesthetic or paralytic agents. Female Hooded Lister rats were anaesthetised and a section of skull overlying somatosensory cortex thinned to translucency. A stainless steel chamber was then secured over the thin cranial window. Following recovery, animals were supported in a harness whilst the head was held by the implanted chamber using a pneumatically driven clamp. Optical imaging and optical imaging spectroscopy (OIS) of somatosensory cortex were performed whilst the contralateral whiskers were stimulated using a computer controlled air-puffer. Imaging sessions lasted approximately 15 min and data were collected for at least three consecutive days. Experiments were then repeated with the animals under urethane anaesthesia. Spectral analysis revealed qualitatively similar haemodynamic response functions across both anaesthetic states. However, our results indicate that the cortical haemodynamic response to somatosensory stimulation is larger by a factor of approximately 5 in the unanaesthetised rat compared with the anaesthetised rat. This preparation may make possible the investigation of the haemodynamic correlates of a broad range of neurological processes in the awake, behaving rodent.

Theory and generalization of Monte Carlo models of the BOLD signal source.

The dependency of the blood oxygenation level dependent (BOLD) signal on underlying hemodynamics is not well understood. Building a forward biophysical model of this relationship is important for the quantitative estimation of the hemodynamic changes and neural activity underlying functional magnetic resonance imaging (fMRI) signals. We have developed a general model of the BOLD signal which can model both intra‐ and extravascular signals for an arbitrary tissue model across a wide range of imaging parameters. The model of the BOLD signal was instantiated as a look‐up‐table (LuT), and was verified against concurrent fMRI and optical imaging measurements of activation induced hemodynamics. Magn Reson Med, 2008.