Regional cerebral blood flow responses to variable frequency whisker stimulation: an autoradiographic analysis

Laminar differences in functional oxygen metabolism in monkey visual cortex measured with calibrated fMRI

Blood-oxygenation-level-dependent functional magnetic resonance imaging (BOLD fMRI) of cortical layers relies on the hemodynamic response and is biased toward large veins on the cortical surface. Functional changes in the cerebral metabolic rate of oxygen (ΔCMRO2) may reflect neural cortical function better than BOLD fMRI, but it is unknown whether the calibrated BOLD model for functional CMRO2 measurement remains valid at high resolution. Here, we measure laminar ΔCMRO2 elicited by visual stimulation in macaque primary visual cortex (V1) and find that ΔCMRO2 peaks in the middle of the cortex, in agreement with autoradiographic measures of metabolism. ΔCMRO2 values in gray matter are similar as found previously. Reductions in CMRO2 are associated with veins at the cortical surface, suggesting that techniques for vein removal may improve the accuracy of the model at very high resolution. However, our results show feasibility of laminar ΔCMRO2 measurement, providing a physiologically meaningful metric of laminar functional metabolism.

Dynamic Contrast Optical Coherence Tomography reveals laminar microvascular hemodynamics in the mouse neocortex in vivo

Studies of flow-metabolism coupling often presume that microvessel architecture is a surrogate for blood flow. To test this assumption, we introduce an in vivo Dynamic Contrast Optical Coherence Tomography (DyC-OCT) method to quantify layer-resolved microvascular blood flow and volume across the full depth of the mouse neocortex, where the angioarchitecture has been previously described. First, we cross-validate average DyC-OCT cortical flow against conventional Doppler OCT flow. Next, with laminar DyC-OCT, we discover that layer 4 consistently exhibits the highest microvascular blood flow, approximately two-fold higher than the outer cortical layers. While flow differences between layers are well-explained by microvascular volume and density, flow differences between subjects are better explained by transit time. Finally, from layer-resolved tracer enhancement, we also infer that microvascular hematocrit increases in deep cortical layers, consistent with predictions of plasma skimming. Altogether, our results show that while the cortical blood supply derives mainly from the pial surface, laminar hemodynamics ensure that the energetic needs of individual cortical layers are met. The laminar trends reported here provide data that links predictions based on the cortical angioarchitecture to cerebrovascular physiology in vivo.

Laminar microvascular transit time distribution in the mouse somatosensory cortex revealed by Dynamic Contrast Optical Coherence Tomography

The transit time distribution of blood through the cerebral microvasculature both constrains oxygen delivery and governs the kinetics of neuroimaging signals such as blood-oxygen-level-dependent functional Magnetic Resonance Imaging (BOLD fMRI). However, in spite of its importance, capillary transit time distribution has been challenging to quantify comprehensively and efficiently at the microscopic level. Here, we introduce a method, called Dynamic Contrast Optical Coherence Tomography (DyC-OCT), based on dynamic cross-sectional OCT imaging of an intravascular tracer as it passes through the field-of-view. Quantitative transit time metrics are derived from temporal analysis of the dynamic scattering signal, closely related to tracer concentration. Since DyC-OCT does not require calibration of the optical focus, quantitative accuracy is achieved even deep in highly scattering brain tissue where the focal spot degrades. After direct validation of DyC-OCT against dilution curves measured using a fluorescent plasma label in surface pial vessels, we used DyC-OCT to investigate the transit time distribution in microvasculature across the entire depth of the mouse somatosensory cortex. Laminar trends were identified, with earlier transit times and less heterogeneity in the middle cortical layers. The early transit times in the middle cortical layers may explain, at least in part, the early BOLD fMRI onset times observed in these layers. The layer-dependencies in heterogeneity may help explain how a single vascular supply manages to deliver oxygen to individual cortical layers with diverse metabolic needs.

Mapping functional connectivity using cerebral blood flow in the mouse brain

Brain function can be assessed from resting-state functional connectivity (rs-fc) maps, most commonly created by analyzing the dynamics of cerebral hemoglobin concentration. Here, we develop the use of Laser Speckle Contrast Imaging (LSCI) for mapping rs-fc using cerebral blood flow (CBF) dynamics. Because LSCI is intrinsically noisy, we used spatial and temporal averaging to sufficiently raise the signal-to-noise ratio for observing robust functional networks. Although CBF-based rs-fc maps in healthy mice are qualitatively similar to simultaneously-acquired [HbO2]-based maps, some quantitative regional differences were observed. These combined flow/concentration maps might help clarify mechanisms involved in network disruption during disease.

Regulation of blood flow in the retinal trilaminar vascular network

Light stimulation evokes neuronal activity in the retina, resulting in the dilation of retinal blood vessels and increased blood flow. This response, named functional hyperemia, brings oxygen and nutrients to active neurons. However, it remains unclear which vessels mediate functional hyperemia. We have characterized blood flow regulation in the rat retina in vivo by measuring changes in retinal vessel diameter and red blood cell (RBC) flux evoked by a flickering light stimulus. We found that, in first- and second-order arterioles, flicker evoked large (7.5 and 5.0%), rapid (0.73 and 0.70 s), and consistent dilations. Flicker-evoked dilations in capillaries were smaller (2.0%) and tended to have a slower onset (0.97 s), whereas dilations in venules were smaller (1.0%) and slower (1.06 s) still. The proximity of pericyte somata did not predict capillary dilation amplitude. Expression of the contractile protein α-smooth muscle actin was high in arterioles and low in capillaries. Unexpectedly, we found that blood flow in the three vascular layers was differentially regulated. Flicker stimulation evoked far larger dilations and RBC flux increases in the intermediate layer capillaries than in the superficial and deep layer capillaries (2.6 vs 0.9 and 0.7% dilation; 25.7 vs 0.8 and 11.3% RBC flux increase). These results indicate that functional hyperemia in the retina is driven primarily by active dilation of arterioles. The dilation of intermediate layer capillaries is likely mediated by active mechanisms as well. The physiological consequences of differential regulation in the three vascular layers are discussed.

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.

Anesthesia and the quantitative evaluation of neurovascular coupling

Anesthesia has broad actions that include changing neuronal excitability, vascular reactivity, and other baseline physiologies and eventually modifies the neurovascular coupling relationship. Here, we review the effects of anesthesia on the spatial propagation, temporal dynamics, and quantitative relationship between the neural and vascular responses to cortical stimulation. Previous studies have shown that the onset latency of evoked cerebral blood flow (CBF) changes is relatively consistent across anesthesia conditions compared with variations in the time-to-peak. This finding indicates that the mechanism of vasodilation onset is less dependent on anesthesia interference, while vasodilation dynamics are subject to this interference. The quantitative coupling relationship is largely influenced by the type and dosage of anesthesia, including the actions on neural processing, vasoactive signal transmission, and vascular reactivity. The effects of anesthesia on the spatial gap between the neural and vascular response regions are not fully understood and require further attention to elucidate the mechanism of vascular control of CBF supply to the underlying focal and surrounding neural activity. The in-depth understanding of the anesthesia actions on neurovascular elements allows for better decision-making regarding the anesthetics used in specific models for neurovascular experiments and may also help elucidate the signal source issues in hemodynamic-based neuroimaging techniques.

Early and late stimulus-evoked cortical hemodynamic responses provide insight into the neurogenic nature of neurovascular coupling

Understanding neurovascular coupling is a prerequisite for the interpretation of results obtained from modern neuroimaging techniques. This study investigated the hemodynamic and neural responses in rat somatosensory cortex elicited by 16 seconds electrical whisker stimuli. Hemodynamics were measured by optical imaging spectroscopy and neural activity by multichannel electrophysiology. Previous studies have suggested that the whisker-evoked hemodynamic response contains two mechanisms, a transient 'backwards' dilation of the middle cerebral artery, followed by an increase in blood volume localized to the site of neural activity. To distinguish between the mechanisms responsible for these aspects of the response, we presented whisker stimuli during normocapnia ('control'), and during a high level of hypercapnia. Hypercapnia was used to 'predilate' arteries and thus possibly 'inhibit' aspects of the response related to the 'early' mechanism. Indeed, hemodynamic data suggested that the transient stimulus-evoked response was absent under hypercapnia. However, evoked neural responses were also altered during hypercapnia and convolution of the neural responses from both the normocapnic and hypercapnic conditions with a canonical impulse response function, suggested that neurovascular coupling was similar in both conditions. Although data did not clearly dissociate early and late vascular responses, they suggest that the neurovascular coupling relationship is neurogenic in origin.

Microvascular organization of the cat inferior colliculus

Brain neural activity depends critically on the blood supply to a given structure. The blood supply can differ within and between divisions, which may have functional significance. We analyzed the microvascular organization of the cat inferior colliculus (IC) to determine if the capillary distribution is homogenous throughout. The IC consists of the central nucleus (CN), the dorsal cortex (DC), and the lateral cortex (LC), each with different roles in auditory behavior and perception. Plastic-embedded tissue was studied from adult cats in 1-μm thick semi-thin sections stained with toluidine blue; tissue was sampled from the IC in a caudal-rostral series of sections. The architectonic subdivisions were drawn independently based on Golgi impregnations. We used the nearest neighbor distance (NND) method to quantify capillary density between subdivisions. Overall, the distribution of capillary density was non-homogenous across the IC. We found significant capillary NND differences between the CN and LC (Mann-Whitney test; p ≤ 0.05), CN and DC (Mann-Whitney test; p ≤ 0.05), and LC and DC (Mann-Whitney test; p ≤ 0.05). The CN had the lowest NND values among all three divisions, indicating the highest capillary density. NND values changed gradually as analysis moved from the center of the IC towards the periphery. The significantly higher microvascular density in the CN may imply that the lemniscal auditory pathway has higher levels of blood flow and metabolic activity than non-lemniscal areas of the IC. The non-homogenous microvascular organization of the IC supports parcellation schemes that delineate three major subdivisions and confirms that the borders between the three regions are not sharp.

Layer-specific BOLD activation in human V1

The neocortex is known to have a distinct laminar structure which has previously been probed in animals using high-resolution fMRI. Detection of layer-specific activation in humans has however to date proven elusive. In this study we demonstrate for the first time such layer-specific activation, specifically at a depth corresponding to layer IV of human primary visual cortex (V1). We used a gradient-echo (GE) sequence at 3T with an isotropic resolution of 0.75 mm, in which a stria at the depth of layer IV was visible in the averaged time series, and could be used as an anatomical landmark. Upon visual stimulation (7.5 Hz flickering checkerboard) the signal increase of 3% in layer IV was significantly higher than in the neighboring laminae. The width of this activation peak was 0.8-1 mm. Based on this result and known laminar organization of the intracortical vasculature we conclude that in the direction perpendicular to the cortical surface the intrinsic spatial resolution of the GE-BOLD fMRI signal is in the submillimetre range. Human laminar fMRI is a significant development which may improve our understanding of intracortical activation patterns and of the way in which different cortical regions interact.

Study of neurovascular coupling by modulating neuronal activity with GABA

Fundamental to the interpretation of neurovascular coupling is determining the neuronal activity that accounts for functional hyperemia. Recently, synaptic and not spiking activity has been found to be responsible for the hemodynamic response. Using pharmacological manipulation in rats, we want to further determine whether the cortical synaptic activity generated by the thalamic input or the subsequent synaptic activity related to secondary cortical processing is driving the hemodynamic response. In this study, we topically applied γ-aminobutyric acid (GABA) in the somatosensory cortex and used electrical forepaw stimulation to evoke neural and vascular activity. In a group of 8 animals, using laminar electrophysiology, we verified that topical application of GABA for 20min does not affect layer IV synaptic activity but reduces subsequent activity in the supragranular and infragranular layers. In another group of 8 animals, we simultaneously measured the electrical and vascular responses with scalp electroencephalography (EEG) and diffuse optical imaging (DOI), respectively. We decomposed somatosensory evoked potentials (SEP) into three major components: P1, N1, and P2, where P1 represents the thalamic input activity originating in layer IV and N1 and P2 represent the subsequent cortical transmissions. We verified that GABA infusion in the somatosensory cortex does not significantly reduce the P1 SEP component but strongly reduces the N1 and P2 components. We found that GABA also elicits a large reduction in the hemodynamic responses, which correlate with the reduction in N1 and P2 components. These results suggest that the hemodynamic response is predominantly driven by cortico-cortical interactions and not by the initial thalamocortical activity in layer IV.

Cerebral blood flow modeling in primate cortex

We report new results on blood flow modeling over large volumes of cortical gray matter of primate brain. We propose a network method for computing the blood flow, which handles realistic boundary conditions, complex vessel shapes, and complex nonlinear blood rheology. From a detailed comparison of the available models for the blood flow rheology and the phase separation effect, we are able to derive important new results on the impact of network structure on blood pressure, hematocrit, and flow distributions. Our findings show that the network geometry (vessel shapes and diameters), the boundary conditions associated with the arterial inputs and venous outputs, and the effective viscosity of the blood are essential components in the flow distribution. In contrast, we show that the phase separation effect has a minor function in the global microvascular hemodynamic behavior. The behavior of the pressure, hematocrit, and blood flow distributions within the network are described through the depth of the primate cerebral cortex and are discussed.

A dynamic model of neurovascular coupling: implications for blood vessel dilation and constriction

Neurovascular coupling in response to stimulation of the rat barrel cortex was investigated using concurrent multichannel electrophysiology and laser Doppler flowmetry. The data were used to build a linear dynamic model relating neural activity to blood flow. Local field potential time series were subject to current source density analysis, and the time series of a layer IV sink of the barrel cortex was used as the input to the model. The model output was the time series of the changes in regional cerebral blood flow (CBF). We show that this model can provide excellent fit of the CBF responses for stimulus durations of up to 16 s. The structure of the model consisted of two coupled components representing vascular dilation and constriction. The complex temporal characteristics of the CBF time series were reproduced by the relatively simple balance of these two components. We show that the impulse response obtained under the 16-s duration stimulation condition generalised to provide a good prediction to the data from the shorter duration stimulation conditions. Furthermore, by optimising three out of the total of nine model parameters, the variability in the data can be well accounted for over a wide range of stimulus conditions. By establishing linearity, classic system analysis methods can be used to generate and explore a range of equivalent model structures (e.g., feed-forward or feedback) to guide the experimental investigation of the control of vascular dilation and constriction following stimulation.

Dynamic magnetic resonance imaging of cerebral blood flow using arterial spin labeling

Modern functional neuroimaging techniques, including positron emission tomography, optical imaging of intrinsic signals, and magnetic resonance imaging (MRI) rely on a tight coupling between neural activity and cerebral blood flow (CBF) to visualize brain activity using CBF as a surrogate marker. Because the spatial and temporal resolution of neuroimaging modalities is ultimately determined by the spatial and temporal specificity of the underlying hemodynamic signals, characterization of the spatial and temporal profiles of the hemodynamic response to focal brain stimulation is of paramount importance for the correct interpretation and quantification of functional data. The ability to properly measure and quantify CBF with MRI is a major determinant of progress in our understanding of brain function. We review the dynamic arterial spin labeling (DASL) method to measure CBF and the CBF functional response with high temporal resolution.

BOLD fMRI and somatosensory evoked potentials are well correlated over a broad range of frequency content of somatosensory stimulation of the rat forepaw

Electrical stimulation of the rat paw is commonly used to study the hemodynamic, metabolic and neuronal mechanisms of functional MRI (fMRI) responses in somatosensory cortex. Several groups have reported good correlation between the blood oxygenation level-dependent (BOLD) fMRI signal and somatosensory evoked potentials (SEPs) using short, typically 300 micros, square stimulation pulses. The spectral power of these short pulses is evenly distributed over a wide range of frequencies and thus the effects of the frequency content of the stimulation pulse on fMRI responses have not been previously described. Here, the effects that different stimulation pulse waveforms with a range of frequency content have on neuronal activity, as measured by SEPs, and on the amplitude of the BOLD fMRI signal in rat somatosensory cortex are investigated. The peak-to-peak SEP amplitudes increased as the power in the high frequency harmonics of the different pulse waveforms increased, using either triangular or sinusoidal stimuli waveforms from 9 Hz to 180 Hz. Similarly, BOLD fMRI response increased with increased high frequency content of the stimulation pulse. There was a linear correlation between SEPs and BOLD fMRI over the full range of frequency content in the stimulations.

The hemo-neural hypothesis: on the role of blood flow in information processing

Brain vasculature is a complex and interconnected network under tight regulatory control that exists in intimate communication with neurons and glia. Typically, hemodynamics are considered to exclusively serve as a metabolic support system. In contrast to this canonical view, we propose that hemodynamics also play a role in information processing through modulation of neural activity. Functional hyperemia, the basis of the functional MRI (fMRI) BOLD signal, is a localized influx of blood correlated with neural activity levels. Functional hyperemia is considered by many to be excessive from a metabolic standpoint, but may be appropriate if interpreted as having an activity-dependent neuro-modulatory function. Hemodynamics may impact neural activity through direct and indirect mechanisms. Direct mechanisms include delivery of diffusible blood-borne messengers and mechanical and thermal modulation of neural activity. Indirect mechanisms are proposed to act through hemodynamic modulation of astrocytes, which can in turn regulate neural activity. These hemo-neural mechanisms should alter the information processing capacity of active local neural networks. Here, we focus on analysis of neocortical sensory processing. We predict that hemodynamics alter the gain of local cortical circuits, modulating the detection and discrimination of sensory stimuli. This novel view of information processing-that includes hemodynamics as an active and significant participant-has implications for understanding neural representation and the construction of accurate brain models. There are also potential medical benefits of an improved understanding of the role of hemodynamics in neural processing, as it directly bears on interpretation of and potential treatment for stroke, dementia, and epilepsy.

Simultaneous laser Doppler flowmetry and arterial spin labeling MRI for measurement of functional perfusion changes in the cortex

This study compares laser Doppler flowmetry (LDF) and arterial spin labeling (ASL) for the measurement of functional changes in cerebral blood flow (CBF). The two methods were applied concurrently in a paradigm of electrical whisker stimulation in the anaesthetised rat. Multi-channel LDF was used, with each channel corresponding to different fiber separation (and thus measurement depth). Continuous ASL was applied using separate imaging and labeling coils at 3 T. Careful experimental set up ensured that both techniques recorded from spatially concordant regions of the barrel cortex, where functional responses were maximal. Strong correlations were demonstrated between CBF changes measured by each LDF channel and ASL in terms of maximum response magnitude and response time-course within a 6-s-long temporal resolution imposed by ASL. Quantitatively, the measurements of the most superficial LDF channels agreed strongly with those of ASL, whereas the deeper LDF channels underestimated consistently the ASL measurement. It was thus confirmed that LDF quantifies CBF changes consistently at a superficial level, and for this case the two methods provided concordant measures of functional CBF changes, despite their essentially different physical principles and spatiotemporal characteristics.

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.

Modulatory role of cyclooxygenase-2 in cerebrovascular coupling

To investigate the role of cyclooxygenase-2 (COX-2) in the cerebrovascular coupling, hemodynamic and neuronal responses to forepaw stimulation were measured in alpha-chloralose-anesthetized rats (N = 18) before and after intravenous administration of Meloxicam (MEL), a preferential COX-2 inhibitor, and following a bolus of prostaglandin E(2) (PGE(2)), a prominent vasodilatatory product of COX-2 catalyzed metabolism of arachidonic acid. The cerebral blood flow (CBF) and blood-oxygenation-level-dependent (BOLD) response was quantified using continuous arterial spin labeling magnetic resonance imaging. Neuronal activity was measured by recording somatosensory-evoked potentials (SEPs) via intracranial electrodes. Both MEL and PGE(2) had a significant effect on the activation-elicited CBF (P < 10(-6)) and BOLD (P < 10(-6)) responses, without affecting the baseline perfusion. Meloxicam decreased brain COX enzymatic activity by 57 +/- 14% and decreased the stimulation-induced CBF response to 32 +/- 2% and BOLD to 46 +/- 1% of their respective pre-drug amplitudes. In turn, PGE(2) bolus resulted in a partial recovery of functional hyperemia, with the CBF response recovering to 52 +/- 3% and the BOLD response to 56 +/- 2% of their values prior to MEL administration. There was no concomitant decrease in either amplitudes or latencies of SEP components. These findings suggest a modulatory role of COX-2 products in the cerebrovascular coupling and provide evidence for existence of a functional metabolic buffer.

Spatial flow-volume dissociation of the cerebral microcirculatory response to mild hypercapnia

The spatial and temporal response of the cerebral microcirculation to mild hypercapnia was investigated via two-photon laser-scanning microscopy. Cortical vessels, traversing the top 200 microm of somatosensory cortex, were visualized in alpha-chloralose-anesthetized Sprague-Dawley rats equipped with a cranial window. Intraluminal vessel diameters, transit times of fluorescent dextrans and red blood cells (RBC) velocities in individual capillaries were measured under normocapnic (PaCO2= 32.6 +/- 2.6 mm Hg) and slightly hypercapnic (PaCO2= 45 +/- 7 mm Hg) conditions. This gentle increase in PaCO2 was sufficient to produce robust and significant increases in both arterial and venous vessel diameters, concomitant to decreases in transit times of a bolus of dye from artery to venule (14%, P < 0.05) and from artery to vein (27%, P < 0.05). On the whole, capillaries exhibited a significant increase in diameter (16 +/- 33%, P < 0.001, n = 393) and a substantial increase in RBC velocities (75 +/- 114%, P < 0.001, n = 46) with hypercapnia. However, the response of the cerebral microvasculature to modest increases in PaCO2 was spatially heterogeneous. The maximal relative dilatation (range: 5-77%; mean +/- SD: 25 +/- 34%, P < 0.001, n = 271) occurred in the smallest capillaries (1.6 microm-4.0 microm resting diameter), while medium and larger capillaries (4.4 microm-6.8 microm resting diameter) showed no significant changes in diameter (P > 0.08, n = 122). In contrast, on average, RBC velocities increased less in the smaller capillaries (39 +/- 5%, P < 0.002, n = 22) than in the medium and larger capillaries (107 +/- 142%, P < 0.003, n = 24). Thus, the changes in capillary RBC velocities were spatially distinct from the observed volumetric changes and occurred to homogenize cerebral blood flow along capillaries of all diameters.

Perfusion-based fMRI: insights from animal models

Modern functional neuroimaging techniques, including functional magnetic resonance imaging (fMRI), positron emission tomography (PET), and optical imaging of intrinsic signals (OIS), rely on a tight coupling between neural activity and cerebral blood flow (CBF) to visualize brain activity using CBF as a surrogate marker. Because CBF is a uniquely defined physiological parameter, fMRI techniques based on CBF contrast have the advantage of being specific to tissue signal change, and the potential to provide more direct and quantitative measures of brain activation than blood oxygenation level-dependent (BOLD)- or cerebral blood volume (CBV)-based techniques. The changes in CBF elicited by increased neural activity are an excellent index of the magnitude of electrical activity. Increases in CBF are more closely localized to the foci of increased electrical activity, and occur more promptly to the stimulus than BOLD- or CBV-based contrast. In addition, CBF-based fMRI is less affected by confounds from venous drainage common to BOLD. Animal studies of brain activation have yielded considerable insights into the advantages of CBF-based fMRI. Based on results provided by animal studies, CBF fMRI may offer a means of better assessing the magnitude, spatial extent, and temporal response of neural activity, and may be more specific to tissue state. These properties are expected to be particularly useful for longitudinal and quantitative fMRI studies.

Spatial extent of oxygen metabolism and hemodynamic changes during functional activation of the rat somatosensory cortex

The spatial extent of the changes in oxy-hemoglobin (HbO), deoxy-hemoglobin (HbR), total hemoglobin concentration (HbT), cerebral blood flow (CBF), and the cerebral metabolic rate of oxygen (CMRO(2)) in response to forepaw and whisker stimulation were compared in the rat somatosensory cortex using a combination of multi-wavelength reflectance imaging and laser speckle contrast imaging of cerebral blood flow. The spatial extents of the response of each hemodynamic parameter and CMRO(2) were found to be comparable at the time of peak response, and at early times following stimulation onset, the spatial extent of the change in HbR was smaller than that of HbO, HbT, CBF, and CMRO(2). In addition, a slight spatial dependence was found in the power law coefficient relating changes in CBF and HbT. Although the CMRO(2) response is a metabolic measure and thus expected to have a more localized response than the hemodynamic parameters, the results presented here suggest that this may not be the case in general, possibly due to the increased sensitivity of optical imaging techniques to superficial cortical layers where the lateral extent of the metabolic and neuronal activation is larger compared to that in layer IV. In addition, we found that the measured spatial extent of the CMRO(2) changes was insensitive to assumptions made in the calculation of the CMRO(2) changes such as baseline hemoglobin concentrations, vascular weighting constants, and wavelength dependence of tissue scattering. Multi-parameter full field imaging of the functional response provides a more complete picture of the hemodynamic response to functional activation including the spatial and temporal estimation of CMRO(2) changes.

Long duration stimuli and nonlinearities in the neural-haemodynamic coupling

Recent studies have shown that the haemodynamic responses to brief (<2 secs) stimuli can be well characterised as a linear convolution of neural activity with a suitable haemodynamic impulse response. In this paper, we show that the linear convolution model cannot predict measurements of blood flow responses to stimuli of longer duration (>2 secs), regardless of the impulse response function chosen. Modifying the linear convolution scheme to a nonlinear convolution scheme was found to provide a good prediction of the observed data. Whereas several studies have found a nonlinear coupling between stimulus input and blood flow responses, the current modelling scheme uses neural activity as an input, and thus implies nonlinearity in the coupling between neural activity and blood flow responses. Neural activity was assessed by current source density analysis of depth-resolved evoked field potentials, while blood flow responses were measured using laser Doppler flowmetry. All measurements were made in rat whisker barrel cortex after electrical stimulation of the whisker pad for 1 to 16 secs at 5 Hz and 1.2 mA (individual pulse width 0.3 ms).

Spatial correlations of laminar BOLD and CBV responses to rat whisker stimulation with neuronal activity localized by Fos expression

The spatial relationship between a measured fMRI signal and its underlying neuronal activity remains unclear. One obstacle is the localization of neuronal activity; another is the spatial resolution of fMRI. In the present study, high-resolution BOLD and CBV fMRI experiments (voxel size: 156 x 156 x 2000 microm3) were conducted in the rat whisker barrel cortex at 3 T; neuronal activity across cortical layers was mapped using the Fos expression technique. Results show that BOLD response is weighted by blood volume and that pixels with high BOLD response can be located at the cortical surface or in deep layers, depending on local vasculature. In contrast to BOLD response, the pixels with high CBV response were consistently clustered in the deep cortical layers. Percentage-CBV change in cortical layers IV-V was 7.3 +/- 1.5%, which was significantly higher than in layers I-III (4.1 +/- 0.9%) and VI (4.3 +/- 0.7%) (mean +/- SEM). The laminar distribution of CBV response correlates well with neuronal activity localized by Fos expression. We conclude that neuronal activity can be inferred from CBV fMRI data with high spatial accuracy. The data indicate that both intracolumn functional connectivity and neurovascular coupling can be studied using CBV fMRI.

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.

Successive depth variations in microvascular distribution of rat somatosensory cortex

Although hemodynamic-based functional brain imaging techniques are powerful tools to explore the brain functions noninvasively, hemodynamic-based signal is strongly affected by spatial configuration of microvessels. Understanding the quantitative relations between microvascular structure and functional activity is therefore significant to make a valid signal interpretation for the imaging techniques. In the present study, we evaluated depth profiles of microvascular distributions in rat somatosensory subfields (barrel field, forelimb region, trunk region and hindlimb region) and characterized depth variations in microvascular structures, such as locations, lengths and directions of microvessels, throughout the cortical layers (I-VI). To obtain the accurate microvascular structure, we made a customized casting method by using confocal laser scanning microscope. We observed that microvascular distribution successively varied throughout the cortical layers (I-VI) and that the maximum number density of microvessels was consistently found in middle layers (III-V). In addition, superficial layers had relatively long microvessels, almost perpendicular to the cortical surface, whereas middle layers had short microvessels propagating in all directions. These regional differences in microvascular structures were closely related to the somatosensory subfields, e.g., barrel field was the greatest number density of microvessels among the investigated subfields. Based on these observations, we compared microvascular profiles with previously reported distribution patterns of tissue partial pressure of oxygen (pO2). The results showed that tissue pO2 was correlated with microvascular distribution in some but not all of the subfields. This finding shows that detailed microvascular profiles are helpful to investigate causal relationships between microvascular structure and functional activities in cerebral cortex.

Dual responses of tissue partial pressure of oxygen after functional stimulation in rat somatosensory cortex

To compare the spatial heterogeneity of brain tissue partial pressure of oxygen (pO(2)) among local brain regions, we focused on functional and anatomical variations in rat somatosensory cortex. Tissue pO(2) was measured by using an oxygen microelectrode with high spatio-temporal resolution, and investigated in three somatosensory areas including hindlimb (HL), forelimb (FL), and trunk region (Tr). Their anatomical structures were determined with histological techniques (Nissl stain). In addition to the measurement of baseline tissue pO(2), we examined temporal shifts in tissue pO(2) distribution elicited by functional stimulation using the brushing stimulation to the hindlimb, forelimb, and trunk regions of the body. We observed that average tissue pO(2) in the Tr (14+/-10 Torr) was significantly lower than those in the HL (25+/-13 Torr) and FL (24+/-13 Torr). Such regional differences in tissue pO(2) were closely related to the cytoarchitectonic variations among these three areas. In addition, the functional stimulation enlarged the regional differences in the pO(2) depending on each somatosensory area; the pO(2) in the HL increased by 3.6+/-2.9% after the stimulation to hindlimb, whereas that in the Tr decreased by -2.9+/-2.5% after the stimulation to trunk region. Such dual responses of tissue pO(2) (i.e. increase or decrease) after the functional stimulation to the corresponding body regions may provide a criterion to clinically predict regions susceptible to tissue hypoxia, because considerable decrease in tissue pO(2) occurred in the Tr showing the lowest baseline pO(2).

Laminar specificity of functional MRI onset times during somatosensory stimulation in rat

The blood oxygenation level-dependent (BOLD) response to somatosensory stimulation was measured in alpha-chloralose-anesthetized rats. BOLD fMRI was obtained at 40-ms temporal resolution and spatial resolution of 200 x 200 x 2,000 microm(3) by using a gated activation paradigm in an 11.7 T MRI. Results show a consistent heterogeneity of fMRI onset times and amplitudes. The earliest onset time (0.59 +/- 0.17 s, n = 9) corresponded anatomically to layer IV, with superficial and deeper layers starting significantly later (1.27 +/- 0.43 s in layers I-III, and 1.11 +/- 0.45 s in layer VI). The amplitude of BOLD signal changes also varied with the cortical depth from the pial surface. Changes in the supragranular layers (8.3%) were 44% bigger than changes in the intermediate layers (5.5%), located only approximately 700 microm below, and 144% larger than the bottom layer (3.5%), located approximately 1.4 mm below the pial surface. The data presented demonstrate that BOLD signal changes have distinct amplitude and temporal characteristics, which vary spatially across cortical layers.

Arteriolar flow recruitment with vitronectin receptor stimulation linked to remote wall shear stress

Our purpose was to investigate the hemodynamic mechanism of flow recruitment within arteriolar networks as stimulated by the vitronectin receptor. Since remote stimulation of the vitronectin receptor initiates a flow-mediated dilation, we tested whether the response could be mimicked by exogenous NO. Arteriolar diameter, red blood cell flux, and velocity were measured in terminal branch arterioles arising from a central feed (8-10 microm) arteriole (together previously defined as an arteriolar network) in the cheek pouch tissue of anesthetized hamsters (n = 29, pentobarbital 70 mg/kg). Terminal branch arterioles were stimulated (micropipette, 10 s) with 10 microg/ml LM609 (agonist, vitronectin receptor), 10(-4) M nitroprusside (SNP), or control. The same amount of flow was recruited into the network by applying LM609 or SNP to only one branch of the network; however, flow distribution differed. With LM609, all recruited flow went directly to the branch that was stimulated. With SNP, recruited flow was distributed evenly throughout the network, no matter where the stimulus was applied. Thus, flow-mediated dilation induced by vitronectin receptor stimulation initiates a robust flow recruitment response, directing flow to the stimulus flow path. In contrast, NO alone initiates flow recruitment to the entire network, along multiple parallel flow paths.

Age-related and laminar-specific dendritic changes in the medial frontal cortex of the rat

Early hypotheses that normal brain aging involves widespread loss of neurons have been revised in light of accumulating evidence that, in most regions of the brain, the number of neurons is stable throughout adulthood and senescence. It is not clear, however, that all aspects of neuronal structure are similarly maintained, and anatomical changes are likely to contribute to age-related declines in cognitive function. The extent and pattern of dendritic branches is one likely target for age-dependent regulation since dendrites remain plastic into adulthood and since dendrites, as the site of most synapses, critically regulate neuronal function. This study quantified the dendritic extent and geometry of superficial and deep pyramidal neurons in the medial frontal cortex of Brown Norway rats from young adulthood through senescence. This region of cortex is of specific interest given its involvement in a variety of cognitive functions that change with age. In the present study, age-related changes in dendritic extent were found to occur with remarkable specificity. Superficial, but not deep, pyramidal neurons exhibited ongoing dendritic growth after 2 months-of-age and then dendritic regression after 18 months-of-age. Apical and basal dendrites were similarly regulated; in each arbor adult growth and regression were limited to terminal dendritic segments. The focal specificity of age-related changes suggests several possible regulatory mechanisms, including regional changes in trophic support and in neuronal activity. Although restricted to specific neuronal populations, dendritic regression in aged animals is likely to contribute to cognitive changes associated with senescence.

Relationship of spikes, synaptic activity, and local changes of cerebral blood flow

The coupling of electrical activity in the brain to changes in cerebral blood flow (CBF) is of interest because hemodynamic changes are used to track brain function. Recent studies, especially those investigating the cerebellar cortex, have shown that the spike rate in the principal target cell of a brain region (i.e. the efferent cell) does not affect vascular response amplitude. Subthreshold integrative synaptic processes trigger changes in the local microcirculation and local glucose consumption. The spatial specificity of the vascular response on the brain surface is limited because of the functional anatomy of the pial vessels. Within the cortex there is a characteristic laminar flow distribution, the largest changes of which are observed at the depth of maximal synaptic activity (i.e. layer IV) for an afferent input system. Under most conditions, increases in CBF are explained by activity in postsynaptic neurons, but presynaptic elements can contribute. Neurotransmitters do not mediate increases in CBF that are triggered by the concerted action of several second messenger molecules. It is important to distinguish between effective synaptic inhibition and deactivation that increase and decrease CBF and glucose consumption, respectively. In summary, hemodynamic changes evoked by neuronal activity depend on the afferent input function (i.e. all aspects of presynaptic and postsynaptic processing), but are totally independent of the efferent function (i.e., the spike rate of the same region). Thus, it is not possible to conclude whether the output level of activity of a region is increased based on brain maps that use blood-flow changes as markers.

Coupling and uncoupling of activity-dependent increases of neuronal activity and blood flow in rat somatosensory cortex

1. Electrical stimulation of the infraorbital nerve was used to examine the coupling between neuronal activity and cerebral blood flow (CBF) in rat somatosensory cortex by laser Doppler flowmetry and extracellular recordings of field potentials. 2. The relationship between field potential (FP) and CBF amplitudes was examined as a function of the stimulus intensity (0--2.0 mA) at fixed frequency (3 Hz). FP amplitudes up to 2.0-2.5 mV were unaccompanied by increases of CBF. Above this threshold, CBF and FP amplitudes increased proportionally. 3. At fixed stimulus intensity of 1.5 mA, CBF increases were highly correlated to FP amplitudes at low frequencies of stimulation (< 2 Hz), but uncoupling was observed at stimulation frequencies of 2--5 Hz. The evoked responses were independent of stimulus duration (8--32 s). 4. The first rise in CBF occurred within the first 0.2 s after onset of stimulation in the upper 0--250 microm of the cortex. Latencies were longer (1.0--1.2 s) in lower cortical layers in which CBF and FP amplitudes were larger. 5. Local AMPA receptor blockade attenuated CBF and FP amplitudes proportionally. 6. This study showed that activity-dependent increases in neuronal activity and CBF were linearly coupled under defined conditions, but neuronal activity was well developed before CBF started to increase. Consequently, the absence of a rise in CBF does not exclude the presence of significant neuronal activity. The CBF increase in upper cortical layers preceded the rise in lower layers suggesting that vessels close to or at the brain surface are the first to react to neuronal activity. The activity-dependent rise in CBF was explained by postsynaptic activity in glutamatergic neurons.

Nerve fibers innervating the cranial and spinal meninges: morphology of nerve fiber terminals and their structural integration

Pachymeninx and leptomeninx of cranial cavity and spine are considerably different in their collagenous fiber texture, cellular composition, vascularization, and innervation. The majority of meningeal nerve fibers terminate as free nerve endings whereas encapsulated and lamellated nerve terminals additionally occur in higher vertebrates including man. With respect to nerve fiber classification, arborization pattern, topography, and organization of the microenvironment at the termination site afferent and efferent nerve terminals are differentiated. Only the dura mater and the pial subcompartment of the leptomeninx possess the morphological prerequisites for neurogenic inflammation. In the current review, the results of morphological studies regarding the meningeal innervation including the sites of CSF (cerebrospinal fluid) production and absorption are discussed with emphasis on their structure-function relationships.