Stimulus frequency dependence of the linear relationship between local cerebral blood flow and field potential evoked by activation of rat somatosensory cortex

Correlation of zero echo time functional MRI with neuronal activity in rats

Zero echo time (zero-TE) pulse sequences provide a quiet and artifact-free alternative to conventional functional magnetic resonance imaging (fMRI) pulse sequences. The fast readouts (<1 ms) utilized in zero-TE fMRI produce an image contrast with negligible contributions from blood oxygenation level-dependent (BOLD) mechanisms, yet the zero-TE contrast is highly sensitive to brain function. However, the precise relationship between the zero-TE contrast and neuronal activity has not been determined. Therefore, we aimed to derive a function to model the temporal dynamics of the zero-TE fMRI signal in response to neuronal activity. Furthermore, we examined the correlation of zero-TE fMRI with neuronal activity across stimulation frequencies. To these ends, we performed simultaneous electrophysiological recordings and zero-TE fMRI in rats subjected to whisker stimulation. The presented impulse response function provides a basis for the statistical modeling of neuronal activity-induced changes in the zero-TE fMRI signal. The temporal characteristics of the zero-TE fMRI response were found to be consistent with the previously postulated non-BOLD hemodynamic origin of the functional contrast. The zero-TE fMRI signal was well predicted by electrophysiological recordings, although systematic stimulation-dependent residuals were also observed, suggesting nonlinearities in neurovascular coupling. We conclude that zero-TE fMRI provides a robust proxy for neuronal activity.

Continuous versus discrete robotic feedback for brain-computer interfaces aimed for neurorehabilitation

Introduction

Brain-Computer Interfaces (BCI) can allow control of external devices using motor imagery (MI) decoded from electroencephalography (EEG). Although BCI have a wide range of applications including neurorehabilitation, the low spatial resolution of EEG, coupled to the variability of cortical activations during MI, make control of BCI based on EEG a challenging task.

Methods

An assessment of BCI control with different feedback timing strategies was performed. Two different feedback timing strategies were compared, comprised by passive hand movement provided by a robotic hand orthosis. One of the timing strategies, the continuous, involved the partial movement of the robot immediately after the recognition of each time segment in which hand MI was performed. The other feedback, the discrete, was comprised by the entire movement of the robot after the processing of the complete MI period. Eighteen healthy participants performed two sessions of BCI training and testing, one with each feedback.

Results

Significantly higher BCI performance (65.4 ± 17.9% with the continuous and 62.1 ± 18.6% with the discrete feedback) and pronounced bilateral alpha and ipsilateral beta cortical activations were observed with the continuous feedback.

Discussion

It was hypothesized that these effects, although heterogenous across participants, were caused by the enhancement of attentional and closed-loop somatosensory processes. This is important, since a continuous feedback timing could increase the number of BCI users that can control a MI-based system or enhance cortical activations associated with neuroplasticity, important for neurorehabilitation applications.

Regional variation in neurovascular coupling and why we still lack a Rosetta Stone

Functional magnetic resonance imaging (fMRI) is the dominant tool in cognitive neuroscience although its relation to underlying neural activity, particularly in the human brain, remains largely unknown. A major research goal, therefore, has been to uncover a 'Rosetta Stone' providing direct translation between the blood oxygen level-dependent (BOLD) signal, the local field potential and single-neuron activity. Here, I evaluate the proposal that BOLD signal changes equate to changes in gamma-band activity, which in turn may partially relate to the spiking activity of neurons. While there is some support for this idea in sensory cortices, findings in deeper brain structures like the hippocampus instead suggest both regional and frequency-wise differences. Relatedly, I consider four important factors in linking fMRI to neural activity: interpretation of correlations between these signals, regional variability in local vasculature, distributed neural coding schemes and varying fMRI signal quality. Novel analytic fMRI techniques, such as multivariate pattern analysis (MVPA), employ the distributed patterns of voxels across a brain region to make inferences about information content rather than whether a small number of voxels go up or down relative to baseline in response to a stimulus. Although unlikely to provide a Rosetta Stone, MVPA, therefore, may represent one possible means forward for better linking BOLD signal changes to the information coded by underlying neural activity. This article is part of the theme issue 'Key relationships between non-invasive functional neuroimaging and the underlying neuronal activity'.

Functional Neuroimaging of Nociceptive and Pain-Related Activity in the Spinal Cord and Brain: Insights From Neurovascular Coupling Studies

Spinal cord and brain processes underlie pain perception, which produces systemic cardiovascular changes. In turn, the autonomic nervous system regulates vascular function in the spinal cord and brain in order to adapt to these systemic changes, while neuronal activity induces local vascular changes. Thus, autonomic regulation and pain processes in the brain and spinal cord are tightly linked and interrelated. The objective of this topical review is to discuss work on neurovascular coupling during nociceptive processing in order to highlight supporting evidence and limitations for the use of cerebral and spinal fMRI to investigate pain mechanisms and spinal nociceptive processes. Work on functional neuroimaging of pain is presented and discussed in relation to available neurovascular coupling studies and related issues. Perspectives on future work are also discussed with an emphasis on differences between the brain and the spinal cord and on different approaches that may be useful to improve current methods, data analyses and interpretation. In summary, this review highlights the lack of data on neurovascular coupling during nociceptive stimulation and indicates that hemodynamic and BOLD responses measured with fMRI may be biased by nonspecific vascular changes. Future neuroimaging studies on nociceptive and pain-related processes would gain further understanding of neurovascular coupling in the brain and spinal cord and should take into account the effects of systemic vascular changes that may affect hemodynamic responses. Anat Rec, 301:1585-1595, 2018. © 2018 Wiley Periodicals, Inc.

Regional Changes in Cerebral Oxygenation During Repeated Passive Movement Measured by Functional Near-infrared Spectroscopy

The aim of this study is to investigate the influence of passive movement repetition frequency at 1.5-Hz and 1-Hz on changes in cerebral oxygenation and assess the temporal properties of these changes using functional near-infrared spectroscopy (fNIRS). No significant differences in systemic hemodynamics were observed between resting and passive movement phases for either 1.5-Hz or 1-Hz trial. Changes in cortical oxygenation as measured by fNIRS in bilateral supplementary motor cortex (SMC), left primary motor cortex (M1), left primary somatosensory cortex (S1), and left posterior association area (PAA) during passive movement of the right index finger revealed greater cortical activity at only 1.5-Hz movement frequency. However, there were no significant differences in the time for peak oxyhemoglobin (oxyHb) among regions (bilateral SMC, 206.4 ± 14.4 s; left M1, 199.1 ± 14.8 s; left S1, 207.3 ± 9.4 s; left PAA, 219.1 ± 10.2 s). Therefore, our results that passive movement above a specific frequency may be required to elicit a changed in cerebral oxygenation, and the times of peak ΔoxyHb did not differ significantly among measured regions.

Electrophysiological correlates of the BOLD signal for EEG-informed fMRI

Electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) are important tools in cognitive and clinical neuroscience. Combined EEG-fMRI has been shown to help to characterise brain networks involved in epileptic activity, as well as in different sensory, motor and cognitive functions. A good understanding of the electrophysiological correlates of the blood oxygen level-dependent (BOLD) signal is necessary to interpret fMRI maps, particularly when obtained in combination with EEG. We review the current understanding of electrophysiological-haemodynamic correlates, during different types of brain activity. We start by describing the basic mechanisms underlying EEG and BOLD signals and proceed by reviewing EEG-informed fMRI studies using fMRI to map specific EEG phenomena over the entire brain (EEG-fMRI mapping), or exploring a range of EEG-derived quantities to determine which best explain colocalised BOLD fluctuations (local EEG-fMRI coupling). While reviewing studies of different forms of brain activity (epileptic and nonepileptic spontaneous activity; cognitive, sensory and motor functions), a significant attention is given to epilepsy because the investigation of its haemodynamic correlates is the most common application of EEG-informed fMRI. Our review is focused on EEG-informed fMRI, an asymmetric approach of data integration. We give special attention to the invasiveness of electrophysiological measurements and the simultaneity of multimodal acquisitions because these methodological aspects determine the nature of the conclusions that can be drawn from EEG-informed fMRI studies. We emphasise the advantages of, and need for, simultaneous intracranial EEG-fMRI studies in humans, which recently became available and hold great potential to improve our understanding of the electrophysiological correlates of BOLD fluctuations.

Repeated BOLD-fMRI imaging of deep brain stimulation responses in rats

Functional magnetic resonance imaging (fMRI) provides a picture of the global spatial activation pattern of the brain. Interest is growing regarding the application of fMRI to rodent models to investigate adult brain plasticity. To date, most rodent studies used an electrical forepaw stimulation model to acquire fMRI data, with α-chloralose as the anesthetic. However, α-chloralose is harmful to animals, and not suitable for longitudinal studies. Moreover, peripheral stimulation models enable only a limited number of brain regions to be studied. Processing between peripheral regions and the brain is multisynaptic, and renders interpretation difficult and uncertain. In the present study, we combined the medetomidine-based fMRI protocol (a noninvasive rodent fMRI protocol) with chronic implantation of an MRI-compatible stimulation electrode in the ventroposterior (VP) thalamus to repetitively sample thalamocortical responses in the rat brain. Using this model, we scanned the forebrain responses evoked by the VP stimulation repeatedly of individual rats over 1 week. Cortical BOLD responses were compared between the 2 profiles obtained at day1 and day8. We discovered reproducible frequency- and amplitude-dependent BOLD responses in the ipsilateral somatosensory cortex (S1). The S1 BOLD responses during the 2 sessions were conserved in maximal response amplitude, area size (size ratio from 0.88 to 0.91), and location (overlap ratio from 0.61 to 0.67). The present study provides a long-term chronic brain stimulation protocol for studying the plasticity of specific neural circuits in the rodent brain by BOLD-fMRI.

Facial vibrotactile stimulation activates the parasympathetic nervous system: study of salivary secretion, heart rate, pupillary reflex, and functional near-infrared spectroscopy activity

We previously found that the greatest salivation response in healthy human subjects is produced by facial vibrotactile stimulation of 89 Hz frequency with 1.9 μm amplitude (89 Hz-S), as reported by Hiraba et al. (2012, 20011, and 2008). We assessed relationships between the blood flow to brain via functional near-infrared spectroscopy (fNIRS) in the frontal cortex and autonomic parameters. We used the heart rate (HRV: heart rate variability analysis in RR intervals), pupil reflex, and salivation as parameters, but the interrelation between each parameter and fNIRS measures remains unknown. We were to investigate the relationship in response to established paradigms using simultaneously each parameter-fNIRS recording in healthy human subjects. Analysis of fNIRS was examined by a comparison of various values between before and after various stimuli (89 Hz-S, 114 Hz-S, listen to classic music, and “Ahh” vocalization). We confirmed that vibrotactile stimulation (89 Hz) of the parotid glands led to the greatest salivation, greatest increase in heart rate variability, and the most constricted pupils. Furthermore, there were almost no detectable differences between fNIRS during 89 Hz-S and fNIRS during listening to classical music of fans. Thus, vibrotactile stimulation of 89 Hz seems to evoke parasympathetic activity.

BOLD consistently matches electrophysiology in human sensorimotor cortex at increasing movement rates: a combined 7T fMRI and ECoG study on neurovascular coupling

Blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging (fMRI) is widely used to measure human brain function and relies on the assumption that hemodynamic changes mirror the underlying neuronal activity. However, an often reported saturation of the BOLD response at high movement rates has led to the notion of a mismatch in neurovascular coupling. We combined BOLD fMRI at 7T and intracranial electrocorticography (ECoG) to assess the relationship between BOLD and neuronal population activity in human sensorimotor cortex using a motor task with increasing movement rates. Though linear models failed to predict BOLD responses from the task, the measured BOLD and ECoG responses from the same tissue were in good agreement. Electrocorticography explained almost 80% of the mismatch between measured- and model-predicted BOLD responses, indicating that in human sensorimotor cortex, a large portion of the BOLD nonlinearity with respect to behavior (movement rate) is well predicted by electrophysiology. The results further suggest that other reported examples of BOLD mismatch may be related to neuronal processes, rather than to neurovascular uncoupling.

Neurovascular saturation thresholds under high intensity auditory stimulation during wake

Coupling between neural activity and hemodynamic responses is important in understanding brain function, interpreting brain-imaging signals, and assessing pathological conditions. Tissue state is a major factor in neurovascular coupling and may alter the relationship between neural and hemodynamic activity. However, most neurovascular-coupling studies are performed under anesthetized or sedated states which may have severe consequences on coupling mechanisms. Our previous studies showed that following prolonged periods of sleep deprivation, evoked hemodynamic responses were muted despite consistent electrical responses, suggesting that sustained neural activity may decrease vascular compliance and limit blood perfusion. To investigate potential perfusion limitations during natural waking conditions, we simultaneously measured evoked response potentials (ERPs) and evoked hemodynamic responses using optical-imaging techniques to increase intensity auditory stimulation. The relationship between evoked hemodynamic responses and integrated ERPs followed a sigmoid relationship where the hemodynamic response approached saturation at lower stimulus intensities than the ERP. If limits in blood perfusion are caused by stretching of the vessel wall, then these results suggest there may be decreased vascular compliance due to sustained neural activity during wake, which could limit vascular responsiveness and local blood perfusion. Conditions that stress cerebral vasculature, such as sleep deprivation and some pathologies (e.g., epilepsy), may further decrease vascular compliance, limit metabolic delivery, and cause tissue trauma. While ERPs and evoked hemodynamic responses provide an indication of the correlated neural activity and metabolic demand, the relationship between these two responses is complex and the different measurement techniques are not directly correlated. Future studies are required to verify these findings and further explore neurovascular coupling during wake by assessing local field potentials, vascular expansion, hemodynamic response localization.

Analysis of time and space invariance of BOLD responses in the rat visual system

Neuroimaging studies of functional magnetic resonance imaging (fMRI) and electrophysiology provide the linkage between neural activity and the blood oxygenation level-dependent (BOLD) response. Here, BOLD responses to light flashes were imaged at 11.7T and compared with neural recordings from superior colliculus (SC) and primary visual cortex (V1) in rat brain--regions with different basal blood flow and energy demand. Our goal was to assess neurovascular coupling in V1 and SC as reflected by temporal/spatial variances of impulse response functions (IRFs) and assess, if any, implications for general linear modeling (GLM) of BOLD responses. Light flashes induced high magnitude neural/BOLD responses reproducibly from both regions. However, neural/BOLD responses from SC and V1 were markedly different. SC signals followed the boxcar shape of the stimulation paradigm at all flash rates, whereas V1 signals were characterized by onset/offset transients that exhibited different flash rate dependencies. We find that IRF(SC) is generally time-invariant across wider flash rate range compared with IRF(V1), whereas IRF(SC) and IRF(V1) are both space invariant. These results illustrate the importance of measured neural signals for interpretation of fMRI by showing that GLM of BOLD responses may lead to misinterpretation of neural activity in some cases.

Diversity of neural-hemodynamic relationships associated with differences in cortical processing during bilateral somatosensory activation in rats

The neural-hemodynamic relationships may vary depending on cortical processing patterns. To investigate how cortical hemodynamics reflects neural activity involving different cortical processing patterns, we delivered electrical stimulation pulses to rat hindpaws, unilaterally or bilaterally, and simultaneously measured electrophysiological (local field potential, LFP < 100 Hz; multiunit activity, MUA>300 Hz) and optical intrinsic signals associated with changes in cerebral blood volume (CBV). Unilateral stimulation evoked neural and optical signals in bilateral primary somatosensory cortices. Ipsilateral optical responses indicating an increased CBV exhibited a peak magnitude of ~30% and mediocaudal shifts relative to contralateral responses. Correlation analyses revealed different scale factors between contralateral and ipsilateral responses in LFP-MUA and LFP-CBV relationships. Bilateral stimulation at varying time intervals evoked hemodynamic responses that were strongly suppressed at 40-ms intervals. This suppression quantitatively reflected suppressed LFP responses to contralateral testing stimulation and not linear summation, with slowly fluctuating LFP responses to ipsilateral conditioning stimulation. Consequently, in the overall responses to bilateral stimulation, CBV-related responses were more linearly correlated with MUA than with LFPs. When extracting high-frequency components (>30 Hz) from LFPs, we found similar scale factors between contralateral and ipsilateral responses in LFP-MUA and LFP-CBV relationships, resulting in significant linear relationships among these components, MUA, and cortical hemodynamics in overall responses to bilateral stimulation. The dependence of LFP-MUA-hemodynamic relationships on cortical processing patterns and the LFP temporal/spectral structure is important for interpreting hemodynamic signals in complex functional paradigms driving diverse cortical processing.

EEG-fMRI integration: a critical review of biophysical modeling and data analysis approaches

The diverse nature of cerebral activity, as measured using neuroimaging techniques, has been recognised long ago. It seems obvious that using single modality recordings can be limited when it comes to capturing its complex nature. Thus, it has been argued that moving to a multimodal approach will allow neuroscientists to better understand the dynamics and structure of this activity. This means that integrating information from different techniques, such as electroencephalography (EEG) and the blood oxygenated level dependent (BOLD) signal recorded with functional magnetic resonance imaging (fMRI), represents an important methodological challenge. In this work, we review the work that has been done thus far to derive EEG/fMRI integration approaches. This leads us to inspect the conditions under which such an integration approach could work or fail, and to disclose the types of scientific questions one could (and could not) hope to answer with it.

Neurovascular coupling is brain region-dependent

Despite recent advances in alternative brain imaging technologies, functional magnetic resonance imaging (fMRI) remains the workhorse for both medical diagnosis and primary research. Indeed, the number of research articles that utilise fMRI have continued to rise unabated since its conception in 1991, despite the limitation that recorded signals originate from the cerebral vasculature rather than neural tissue. Consequently, understanding the relationship between brain activity and the resultant changes in metabolism and blood flow (neurovascular coupling) remains a vital area of research. In the past, technical constraints have restricted investigations of neurovascular coupling to cortical sites and have led to the assumption that coupling in non-cortical structures is the same as in the cortex, despite the lack of any evidence. The current study investigated neurovascular coupling in the rat using whole-brain blood oxygenation level-dependent (BOLD) fMRI and multi-channel electrophysiological recordings and measured the response to a sensory stimulus as it proceeded through brainstem, thalamic and cortical processing sites - the so-called whisker-to-barrel pathway. We found marked regional differences in the amplitude of BOLD activation in the pathway and non-linear neurovascular coupling relationships in non-cortical sites. The findings have important implications for studies that use functional brain imaging to investigate sub-cortical function and caution against the use of simple, linear mapping of imaging signals onto neural activity.

Mapping cortical representations of the rodent forepaw and hindpaw with BOLD fMRI reveals two spatial boundaries

Electrical stimulation of the rat forepaw and hindpaw was employed to study the spatial distribution of BOLD fMRI. Averaging of multiple fMRI sessions significantly improved the spatial stability of the BOLD signal and enabled quantitative determination of the boundaries of the BOLD fMRI maps. The averaged BOLD fMRI signal was distributed unevenly over the extent of the map and the data at the boundaries could be modeled with major and minor spatial components. Comparison of three-dimensional echo-planar imaging (EPI) fMRI at isotropic 300 μm resolution demonstrated that the border locations of the major spatial component of BOLD signal did not overlap between the forepaw and hindpaw maps. Interestingly, the border positions of the minor BOLD fMRI spatial components extended significantly into neighboring representations. Similar results were found for cerebral blood volume (CBV) weighted fMRI obtained using iron oxide particles, suggesting that the minor spatial components may not be due to vascular mislocalization typically associated with BOLD fMRI. Comparison of the BOLD fMRI maps of the forepaw and hindpaw to histological determination of these representations using cytochrome oxidase (CO) staining demonstrated that the major spatial component of the BOLD fMRI activation maps accurately localizes the borders. Finally, 2-3 weeks following peripheral nerve denervation, cortical reorganization/plasticity at the boundaries of somatosensory limb representations in adult rat brain was studied. Denervation of the hindpaw caused a growth in the major component of forepaw representation into the adjacent border of hindpaw representation, such that fitting to two components no longer led to a better fit as compared to using one major component. The border of the representation after plasticity was the same as the border of its minor component in the absence of any plasticity. It is possible that the minor components represent either vascular effects that extend from the real neuronal representations or the neuronal communication between neighboring regions. Either way the results will be useful for studying mechanisms of plasticity that cause alterations in the boundaries of neuronal representations.

Evoked local field potentials can explain temporal variation in blood oxygenation level-dependent responses in rat somatosensory cortex

The aim of this study was to explain the temporal variations between subjects in the blood oxygenation level-dependent (BOLD) response. Somatosensory responses were elicited with the electrical forepaw stimulus at a frequency of 10 Hz in urethane-anesthetized rats, and functional magnetic resonance imaging (fMRI) with BOLD contrast and local field potential (LFP) measurements were performed simultaneously. BOLD fMRI activation was evaluated by two different models, one based on the stimulus paradigm (the block model) and the other on the simultaneously measured evoked LFP responses. In the initial analysis, the LFP model captured the BOLD activation in the primary somatosensory cortex in all cases, and the block model in 10 of 12 rats. A statistical comparison of the two models revealed that the LFP-derived model was able to explain additional BOLD variation over the block model in the somatosensory cortex in nine of 12 rats. These results suggest that there is more information regarding neuronal activity in the BOLD signal than can be exploited using the block model alone. Furthermore, the hemodynamic coupling remains unchanged in the case of temporally variable BOLD signals.

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.

Neurovascular and Neurometabolic Couplings in Dynamic Calibrated fMRI: Transient Oxidative Neuroenergetics for Block-Design and Event-Related Paradigms

Functional magnetic resonance imaging (fMRI) with blood-oxygenation level dependent (BOLD) contrast is an important tool for mapping brain activity. Interest in quantitative fMRI has renewed awareness in importance of oxidative neuroenergetics, as reflected by cerebral metabolic rate of oxygen consumption(CMRO2), for supporting brain function. Relationships between BOLD signal and the underlying neurophysiological parameters have been elucidated to allow determination of dynamic changes inCMRO2 by "calibrated fMRI," which require multi-modal measurements of BOLD signal along with cerebral blood flow (CBF) and volume (CBV). But how doCMRO2 changes, steady-state or transient, derived from calibrated fMRI compare with neural activity recordings of local field potential (LFP) and/or multi-unit activity (MUA)? Here we discuss recent findings primarily from animal studies which allow high magnetic fields studies for superior BOLD sensitivity as well as multi-modal CBV and CBF measurements in conjunction with LFP and MUA recordings from activated sites. A key observation is that while relationships between neural activity and sensory stimulus features range from linear to non-linear, associations between hyperemic components (BOLD, CBF, CBV) and neural activity (LFP, MUA) are almost always linear. More importantly, the results demonstrate good agreement between the changes inCMRO2 and independent measures of LFP or MUA. The tight neurovascular and neurometabolic couplings, observed from steady-state conditions to events separated by <200 ms, suggest rapid oxygen equilibration between blood and tissue pools and thus calibrated fMRI at high magnetic fields can provide high spatiotemporal mapping ofCMRO2 changes.

Modulation of somatosensory-evoked cortical blood flow changes by GABAergic inhibition of the nucleus basalis of Meynert in urethane-anaesthetized rats

Vascular changes associated with brain functions are thought to be tightly coupled with neuronal activity through neuronal glucose consumption or the local release of vasoactive agents. In contrast, another view suggests that cortical blood flow is strongly regulated by the nucleus basalis of Meynert (NBM), independently of regional metabolism. Thus, although cortical regional cerebral blood flow (rCBF) variations induced by somatosensory stimulation are strongly linked to neuronal activity, they may also be partly controlled by the NBM. In the present study, cortical rCBF alterations in response to innocuous brushing of the hindlimb (HL) were investigated by laser speckle contrast imaging. The contribution of NBM to these changes was examined after injection of the GABAergic agonist muscimol into the right NBM, allowing comparison of somatosensory-evoked cortical rCBF modifications before and after NBM inactivation. As expected, HL brushing elicited a robust rCBF increase in the contralateral parietal cortex (PC), over the representation of the HL. However, these alterations were decreased, by approximately 40%, in the hemisphere ipsilateral to muscimol inactivation of NBM, whereas vehicle injection did not produce any significant variation. The results demonstrate that cortical rCBF changes induced by somatosensory stimulation are partly regulated by NBM.

Frequency-dependent neural activity, CBF, and BOLD fMRI to somatosensory stimuli in isoflurane-anesthetized rats

Inhalation anesthetics (e.g. isoflurane) are preferable for longitudinal fMRI experiments in the same animals. We previously implemented isoflurane anesthesia for rodent forepaw stimulation studies, and optimized the stimulus parameters with short stimuli (1-3-s long stimulation with ten electric pulses). These parameters, however, may not be applicable for long periods of stimulation because repetitive stimuli induce neural adaptation. Here we evaluated frequency-dependent responses (pulse width of 1.0 ms and current of 1.5 mA) for 30-s long stimulation under 1.3-1.5% isoflurane anesthesia. The cerebral blood flow (CBF) response (using laser Doppler flowmetry: CBF(LDF)) and field potential (FP) changes were simultaneously measured for nine stimulus frequencies (1-24 Hz). CBF (using arterial spin labeling: CBF(ASL)) and blood oxygenation level dependent (BOLD) fMRI responses were measured at 9.4 T for four stimulus frequencies (1.5-12 Hz). Higher stimulus frequencies (12-24 Hz) produced a larger FP per unit time initially, but decreased more rapidly later due to neural adaptation effects. On the other hand, lower stimulus frequencies (1-3 Hz) induced smaller, but sustained FP activities over the entire stimulus period. Similar frequency-dependencies were observed in CBF(LDF), CBF(ASL) and BOLD responses. A linear relationship between FP and CBF(LDF) was observed for all stimulus frequencies. Stimulation frequency for the maximal cumulative neural and hemodynamic changes is dependent on stimulus duration; 8-12 Hz for short stimulus durations (<10s) and 6-8 Hz for 30-s stimulation. Our findings suggest that neural adaptation should be considered in determining the somatosensory stimulation frequency and duration under isoflurane anesthesia.

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.

How and when the fMRI BOLD signal relates to underlying neural activity: the danger in dissociation

Functional magnetic resonance imaging (fMRI) has become the dominant means of measuring behavior-related neural activity in the human brain. Yet the relation between the blood oxygen-level dependent (BOLD) signal and underlying neural activity remains an open and actively researched question. A widely accepted model, established for sensory neo-cortex, suggests that the BOLD signal reflects peri-synaptic activity in the form of the local field potential rather than the spiking rate of individual neurons. Several recent experimental results, however, suggest situations in which BOLD, spiking, and the local field potential dissociate. Two different models are discussed, based on the literature reviewed to account for this dissociation, a circuitry-based and vascular-based explanation. Both models are found to account for existing data under some testing situations and in certain brain regions. Because both the vascular and local circuitry-based explanations challenge the BOLD-LFP coupling model, these models provide guidance in predicting when BOLD can be expected to reflect neural processing and when the underlying relation with BOLD may be more complex than a direct correspondence.

Cerebral control of face, jaw, and tongue movements in awake cats: Changes in regional cerebral blood flow during lateral feeding

Mastication is achieved by cooperation among facial, masticatory, and lingual muscles. However, cortical control in cats for the masticatory performance is processed by two systems: facial movement processed by facial SI (the first somatosensory cortex), area C, and area M (motor areas), and jaw and tongue movements performed by intraoral SI, masticatory area, and area P (motor area). In particular, outputs from area P organized in the corticobulbar tract are projected bilaterally in the brainstem. In this present study, the aim is to explore changes in the regional cerebral blood flow (rCBF) in the facial SI, area M, and area P during trained lateral feeding (licking or chewing from the right or left side) of milk, fish paste, and small dry fish. The rCBF in area M showed contralateral dominance, and rCBF in area P during chewing or licking from the right or left side was almost the same value. Furthermore, activities of genioglossus and masseter muscles in the left side showed almost the same values during licking of milk and of fish paste, and chewing of small dry fish during lateral feeding. These findings suggest that the cortical process for facial, jaw, and tongue movements may be regulated by the contralateral dominance of area M and the bilateral one of area P.

Cerebral oxygen demand for short-lived and steady-state events

Because of the importance of oxidative energetics for cerebral function, extraction of oxygen consumption (CMR(O2)) from blood oxygenation level-dependent (BOLD) signal using multi-modal measurements of blood flow (CBF) and volume (CBV) has become an accepted functional magnetic resonance imaging (fMRI) technique. This approach, termed calibrated fMRI, is based on a biophysical model which describes tissue oxygen extraction at steady-state. A problem encountered for calculating dynamic CMR(O2) relates to concerns whether the conventional BOLD model can be applied transiently. In particular, it is unclear whether calculation of CMR(O2) differs between short and long stimuli. Linearity was experimentally demonstrated between BOLD-related components and neural activity, thereby making it possible to use calibrated fMRI in a dynamic manner. We used multi-modal fMRI and electrophysiology, in alpha-chloralose anesthetized rats during forepaw stimulation to show that respective transfer functions (of BOLD, CBV, CBF) generated by deconvolution with neural activity are time invariant, for events in the millisecond to minute range. These results allowed extraction of a significant component of the BOLD signal that can be ascribed to CMR(O2) transients. We discuss the importance of minimizing residual signal, represented by the difference between modeled and raw signals, in convolution analysis of multi-modal signals.

Study of neurovascular coupling in humans via simultaneous magnetoencephalography and diffuse optical imaging acquisition

By combining diffuse optical imaging (DOI) and magnetoencephalography (MEG) we investigate neurovascular coupling non-invasively in human subjects using median-nerve stimulation. Previous fMRI studies have shown a habituation effect in the hemodynamic blood oxygen level-dependent (BOLD) response for stimulation periods longer than 2 s. With DOI and MEG we can test whether this effect in hemodynamic response can be accounted for by a habituation effect in the neural response. Our experimental results show that the habituation effect in the hemodynamic response is stronger than that in the earliest cortical neural response (N20). Using a linear convolution model to predict hemodynamic responses we found that including late neural components (> or = 30 ms) improves the prediction of the hemoglobin response. This finding suggests that in addition to the initial evoked-response deflections related to the talamic afferent input, later cortical activity is needed to predict the hemodynamic response.

Oxidative neuroenergetics in event-related paradigms

Energetic basis of neural activity provides a solid foundation for noninvasive neuroimaging with calibrated functional magnetic resonance imaging (fMRI). Calculating dynamic changes in cerebral oxidative energy utilization (CMR(O(2))) is limited by uncertainties about whether or not the conventional blood oxygenation level-dependent (BOLD) model can be applied transiently using multimodal measurements of blood flow (CBF) and volume (CBV) that affect the BOLD signal. A prerequisite for dynamic calibrated fMRI is testing the linearity of multimodal signals within a temporal regimen, as assessed by signal strength (i.e., both intensity and width). If each hyperemic component (BOLD, CBV, CBF) is demonstrated to be linear with neural activity under various experimental conditions, then the respective transfer functions generated by deconvolution with neural activity should be time invariant and thus could potentially be used for calculating CMR(O(2)) transients. Hyperemic components were investigated at 11.7 T in alpha-chloralose-anesthetized rats and combined with electrophysiological recordings of local field potential (LFP) and multiunit activity (MUA) from the cortex during forepaw stimulation, in which stimulus number and frequency were varied. Although relationships between neural activity and stimulus features ranged from linear to nonlinear, associations between hyperemic components and neural activity were linear. Specific to each hyperemic component, a universal transfer function (with LFP or MUA) yielded predictions in agreement with experimental measurements. The results identified a component of the BOLD signal that can be attributed to significant changes in CMR(O(2)), even for temporal events separated by <200 ms.

Frequency-dependent tactile responses in rat brain measured by functional MRI

We measured frequency-dependent functional MRI (fMRI) activations (at 11.7 T) in the somatosensory cortex with whisker and forepaw stimuli in the same alpha-chloralose anesthetized rats. Whisker and forepaw stimuli were attained by computer-controlled pulses of air puffs and electrical currents, respectively. Air puffs deflected (+/-2 mm) the chosen whisker(s) in the right snout in the rostral to caudal direction, and electrical currents (2 mA amplitude, 0.3 ms duration) stimulated the left forepaw with subcutaneous copper electrodes placed between the second and fourth digits. In the same subject, unimodal stimulation of whisker and forepaw gave rise to significant blood oxygen level-dependent (BOLD) signal increases in corresponding contralateral somatosensory areas of whisker barrel field (S1BF) and forelimb (S1FL), respectively, with no significant spatial overlap between these regions. The BOLD responses in S1(BF) and S1(FL) regions were found to be differentially variable with frequency of each stimulus type. In the S1BF, a linear increase in the BOLD response was observed with whisker stimulation frequency of up to approximately 12 Hz, beyond which the response seemed to saturate (and/or slightly attenuate) up to the maximum frequency studied (i.e. 30 Hz). In the S1FL, the magnitude of the BOLD response was largest at forepaw stimulation frequency between 1.5 and 3 Hz, beyond which the response diminished with little or no activity at frequencies higher than 20 Hz. The volume of tissue activated by each stimulus type followed a similar pattern to that of the stimulation frequency dependence. These results of bimodal whisker and forepaw stimuli in the same subject may provide a framework to study interactions of different tactile modules, with both fMRI and neurophysiology (i.e. inside and outside the magnet).

Coupling between somatosensory evoked potentials and hemodynamic response in the rat

We studied the relationship between somatosensory evoked potentials (SEP) recorded with scalp electroencephalography (EEG) and hemoglobin responses recorded non-invasively with diffuse optical imaging (DOI) during parametrically varied electrical forepaw stimulation in rats. Using these macroscopic techniques we verified that the hemodynamic response is not linearly coupled to the somatosensory evoked potentials, and that a power or threshold law best describes the coupling between SEP and the hemoglobin response, in agreement with the results of most invasive studies. We decompose the SEP response in three components (P1, N1, and P2) to determine which best predicts the hemoglobin response. We found that N1 and P2 predict the hemoglobin response significantly better than P1 and the input stimuli (S). Previous electrophysiology studies reported in the literature show that P1 originates in layer IV directly from thalamocortical afferents, while N1 and P2 originate in layers I and II and reflect the majority of local cortico-cortical interactions. Our results suggest that the evoked hemoglobin response is driven by the cortical synaptic activity and not by direct thalamic input. The N1 and P2 components, and not P1, need to be considered to correctly interpret neurovascular coupling.

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.

Fine detail of neurovascular coupling revealed by spatiotemporal analysis of the hemodynamic response to single whisker stimulation in rat barrel cortex

The spatial resolution of hemodynamic-based neuroimaging techniques, including functional magnetic resonance imaging, is limited by the degree to which neurons regulate their blood supply on a fine scale. Here we investigated the spatial detail of neurovascular events with a combination of high spatiotemporal resolution two-dimensional spectroscopic optical imaging, multichannel electrode recordings and cytochrome oxidase histology in the rodent whisker barrel field. After mechanical stimulation of a single whisker, we found two spatially distinct cortical hemodynamic responses: a transient response in the "upstream" branches of surface arteries and a later highly localized increase in blood volume centered on the activated cortical column. Although the spatial representation of this localized response exceeded that of a single "barrel," the spread of hemodynamic activity accurately reflected the neural response in neighboring columns rather than being due to a passive "overspill." These data confirm hemodynamics are capable of providing accurate "single-condition" maps of neural activity.

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.

Apparent diffusion time of oxygen from blood to tissue in rat cerebral cortex: implication for tissue oxygen dynamics during brain functions

To investigate the dynamics of tissue oxygen demand and supply during brain functions, we simultaneously recorded Po(2) and local cerebral blood flow (LCBF) with an oxygen microelectrode and laser Doppler flowmetry, respectively, in rat somatosensory cortex. Electrical hindlimb stimuli were applied for 1, 2, and 5 s to vary the duration of evoked cerebral metabolic rate of oxygen (CMR(O(2))). The electrical stimulation induced a robust increase in Po(2) (4-9 Torr at peak) after an increase in LCBF (14-26% at peak). A consistent lag of approximately 1.2 s (0.6-2.3 s for individual animals) in the Po(2) relative to LCBF was found, irrespective of stimulus length. It is argued that the lag in Po(2) was predominantly caused by the time required for oxygen to diffuse through tissue. During brain functions, the supply of fresh oxygen further lagged because of the latency of LCBF onset ( approximately 0.4 s). The results indicate that the tissue oxygen supports excess demand until the arrival of fresh oxygen. However, a large drop in Po(2) was not observed, indicating that the evoked neural activity demands little extra oxygen or that the time course of excess demand is as slow as the increase in supply. Thus the dynamics of Po(2) during brain functions predominantly depend on the time course of LCBF. Possible factors influencing the lag between demand and supply are discussed, including vascular spacing, reactivity of the vessels, and diffusivity of oxygen.

Functional uncoupling of hemodynamic from neuronal response by inhibition of neuronal nitric oxide synthase

The cerebrovascular coupling under neuronal nitric oxide synthase (nNOS) inhibition was investigated in alpha-chloralose anesthetized rats. Cerebral blood flow (CBF), cerebral blood volume (CBV), and blood oxygenation level dependent (BOLD) responses to electrical stimulation of the forepaw were measured before and after an intraperitoneal bolus of 7-nitroindazole (7-NI), an in vivo inhibitor of the neuronal isoform of nitric oxide synthase. Neuronal activity was measured by recording somatosensory-evoked potentials (SEPs) via intracranial electrodes. 7-Nitroindazole produced a significant attenuation of the activation-elicited CBF (P<10(-6)), CBV (P<10(-6)), and BOLD responses (P<10(-6)), without affecting the baseline perfusion level. The average DeltaCBF was nulled, while DeltaBOLD and DeltaCBV decreased to approximately 30% of their respective amplitudes before 7-NI administration. The average SEP amplitude decreased (P<10(-5)) to approximately 60% of its pretreatment value. These data describe a pharmacologically induced uncoupling between neuronal and hemodynamic responses to functional activation, and provide further support for the critical role of neuronally produced NO in the cerebrovascular coupling.

Nonlinear neurovascular coupling in rat sensory cortex by activation of transcallosal fibers

Functional neuroimaging and normal brain function rely on the robust coupling between neural activity and cerebral blood flow (CBF), that is neurovascular coupling. We examined neurovascular coupling in rat sensory cortex in response to direct stimulation of transcallosal pathways, which allows examination of brain regions inaccessible to peripheral stimulation techniques. Using laser-Doppler flowmetry to record CBF and electrophysiologic recordings of local field potentials (LFPs), we show an exponential relation between CBF responses and summed LFP amplitudes. Hemodynamic responses were dependent on glutamate receptor activation. CNQX, an AMPA receptor blocker, strongly attenuated evoked CBF responses and LFP amplitudes at all stimulation frequencies. In comparison, N-methyl D-aspartate (NMDA) receptor blockade by MK801 attenuated CBF responses at high (>7 Hz) but not low (<7 Hz) stimulation frequencies, without affecting evoked LFP amplitudes. This shows the limitation of using LFP amplitudes as indicators of synaptic activity. 7-Nitroindazole, a neuronal nitric oxide synthase inhibitor, and indomethacin, a nonspecific cyclooxygenase inhibitor, attenuated the hemodynamic responses by 50%+/-1% and 48%+/-1%, respectively, without affecting LFP amplitudes. The data suggest that preserved activity of both AMPA and NMDA receptors is necessary for the full CBF response evoked by stimulation of rodent interhemispheric connections. AMPA receptor activation gives rise to a measurable LFP, but NMDA receptor activation does not. The lack of a measurable LFP from neural processes that contribute importantly to CBF may explain some of the difficulties in transforming extracellular current or voltage measurements to a hemodynamic response.

Modelling the role of excitatory and inhibitory neuronal activity in the generation of the BOLD signal

A biophysical model of the coupling between neuronal activity and the BOLD signal that allows for explicitly evaluating the role of both excitatory and inhibitory activity is formulated in this paper. The model is based on several physiological assumptions. Firstly, in addition to glycolysis, the "glycogen shunt" is assumed for excitatory synapses as a mechanism for energy production in the astrocytes. As a result, oxygen-to-glucose index (OGI) is not constant but varies with excitatory neuronal activity. In contrast, a constant OGI=6 (glycolysis) is assumed for inhibitory synapses. Finally we assume that cerebral blood flow is not directly controlled by energy usage, but it is only related to excitatory activity. Simulations' results show that increases in excitatory activity amplify the oscillations associated with the transient BOLD response, by increasing the initial dip, the maximum, and the post-stimulus undershoot of the signal. In contrast, increasing the inhibitory activity evoked an overall decrease of the BOLD signal along the whole time interval of the response. Simultaneous increase of both types of activity is then expected to reinforce the initial dip and the post-stimulus undershoot, while respective effects on the maximum tend to counteract each other. Two mechanisms for negative BOLD response (NBS) generation were predicted by the model: (i) when inhibition was present alone or together with low activation levels and (ii) when deactivation occurred independently of the accompanying inhibition level. Interestingly, NBS was associated with negative oxygen consumption changes only for the case of mechanism (ii).

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.

Rate code and temporal code for frequency of whisker stimulation in rat primary and secondary somatic sensory cortex

We recorded responses to frequencies of whisker stimulation from 479 neurons in primary (S1) and secondary (S2) somatic sensory cortex of 26 urethane-anesthetized rats. Five whiskers on the right side of the snout were deflected with air puffs at seven frequencies between 1 and 18/s. In left S1 (barrels and septa) and S2, subsets of neurons (5%) responded to whisker stimulation across the entire range of frequencies with > or = 1 electrical discharges/ten stimuli (full responders). In contrast, 60% of the recorded cells responded above threshold only at stimulus frequencies below 6/s and 35% remained subthreshold at all frequencies tested. Thus, the full responders are unique in that they were always responsive and appeared particularly suited to facilitate a dynamic, broadband processing of stimulus frequency. Full responders were most responsive at 1 stimulus/s, and showed greatest synchrony with whisker motion at 18 stimuli/s. The barrel cells responded with the greatest temporal accuracy between 3 and 15 stimuli/s. The septum cells responded less accurately, but maintained their accuracy at all frequencies. Only septum cells continued to increase their discharge rate with increasing stimulus frequency. The S2 cells discharged with lowest temporal accuracy modulated only by stimulus frequencies < or = 6/s and exhibited the steepest decrease in discharge/stimulus with increasing stimulus frequency. Our observations suggest that full responders in the septa are well suited to encode high frequencies of whisker stimulation in timing and rate of discharge. The barrel cells, in contrast, showed the strongest temporal coding at stimulus frequencies in the middle range, and S2 cells were most sensitive to differences in low frequencies. The ubiquitous decline in discharge/stimulus in S1 and S2 may explain the decrease in blood flow observed at increasing stimulus frequency with functional imaging.

Nonlinear correlation between field potential and local cerebral blood flow in rat somatosensory cortex evoked by changing the stimulus current

The relationship between local cerebral blood flow (LCBF) and field potential (FP) evoked by hindpaw stimulation in rat somatosensory cortex has been investigated while changing stimulus current. The change in LCBF was measured using laser-Doppler flowmetry and the field potential was acquired using a tungsten electrode inserted into the cortex of alpha-chloralose-anesthetized rats. The cortex was activated by electrical stimulation of the hind paw with 5 Hz pulses (0.1 ms) applied at currents of 1.0, 1.5, 2.0 and 2.5 mA for 5 s. It was found that the summed FP is nonlinear with respect to stimulus current, whereas the integrated LCBF response is linear across the range of currents used in the experiment. This means that the relationship between the summed FP and integrated LCBF is nonlinear as a function of stimulus current.

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

Reading vascular changes in brain imaging: is dendritic calcium the key?

A key goal in functional neuroimaging is to use signals that are related to local changes in metabolism and blood flow to track the neuronal correlates of mental activity. Recent findings indicate that the dendritic processing of excitatory synaptic inputs correlates more closely than the generation of spikes with brain imaging signals. The correlation is often nonlinear and context-sensitive, and cannot be generalized for every condition or brain region. The vascular signals are mainly produced by increases in intracellular calcium in neurons and possibly astrocytes, which activate important enzymes that produce vasodilators to generate increments in flow and the positive blood oxygen level dependent signal. Our understanding of the cellular mechanisms of functional imaging signals places constraints on the interpretation of the data.

Principal neuron spiking: neither necessary nor sufficient for cerebral blood flow in rat cerebellum

Neuronal activity, cerebral blood flow, and metabolic responses are all strongly coupled, although the mechanisms behind the coupling remain unclear. One of the key questions is whether or not increases in spiking activity in the stimulated neurons are sufficient to drive the activity-dependent rises in cerebral blood flow (CBF) that form the basis of the signals used in functional neuroimaging such as the blood oxygen level-dependent (BOLD) signal. To this end the present study examined the effect of enhanced spike activity per se on CBF in rat cerebellar cortex under conditions of disinhibition, achieved by blocking GABA(A) receptors using either bicuculline or picrotoxin. Purkinje cell spiking activity and local field potentials were recorded by glass microelectrodes, and laser Doppler flowmetry was used to monitor CBF. Disinhibition increased Purkinje cell spiking rate to 200-300% of control without incurring any increase in basal CBF. This demonstrates that increased spike activity per se is not sufficient to affect basal CBF. The neurovascular coupling between excitatory synaptic activity and CBF responses evoked by inferior olive (climbing fibre) stimulation was preserved during disinhibition. Thus, the unchanged basal CBF in the presence of the dramatic rise in Purkinje cell spiking rate was not explained by impaired synaptic activity-CBF coupling. On the basis of our previous and the present studies, we conclude that increased spiking activity of principal neurons is neither sufficient nor necessary to elicit CBF responses and in turn BOLD signals, and that activation-dependent vascular signals reflect excitatory synaptic activity.