Nonlinear coupling of neural activity and CBF in rodent barrel cortex

Assessment of vibration modulated regional cerebral blood flow with MRI

Human brain experiences vibration of certain magnitude and frequency during various physical activities such as vehicle transportation and machine operation, which may cause traumatic brain injury or other brain diseases. However, the mechanisms of brain pathogenesis due to vibration are not fully elucidated due to the lack of techniques to study brain functions while applying vibration to the brain at a specific magnitude and frequency. Here, this study reported a custom-built head-worn electromagnetic actuator that applied vibration to the brain in vivo at an accurate frequency inside a magnetic resonance imaging scanner while cerebral blood flow (CBF) was acquired. Using this technique, CBF values from 45 healthy volunteers were quantitatively measured immediately following vibration at 20, 30, 40 Hz, respectively. Results showed increasingly reduced CBF with increasing frequency at multiple regions of the brain, while the size of the regions expanded. Importantly, the vibration-induced CBF reduction regions largely fell inside the brain's default mode network (DMN), with about 58 or 46% overlap at 30 or 40 Hz, respectively. These findings demonstrate that vibration as a mechanical stimulus can change strain conditions, which may induce CBF reduction in the brain with regional differences in a frequency-dependent manner. Furthermore, the overlap between vibration-induced CBF reduction regions and DMN suggested a potential relationship between external mechanical stimuli and cognitive functions.

Neurocapillary-Modulation

OBJECTIVES

We consider two consequences of brain capillary ultrastructure in neuromodulation. First, blood-brain barrier (BBB) polarization as a consequence of current crossing between interstitial space and the blood. Second, interstitial current flow distortion around capillaries impacting neuronal stimulation.

MATERIALS AND METHODS

We developed computational models of BBB ultrastructure morphologies to first assess electric field amplification at the BBB (principle 1) and neuron polarization amplification by the presence of capillaries (principle 2). We adapt neuron cable theory to develop an analytical solution for maximum BBB polarization sensitivity.

RESULTS

Electrical current crosses between the brain parenchyma (interstitial space) and capillaries, producing BBB electric fields (EBBB) that are >400x of the average parenchyma electric field (ĒBRAIN), which in turn modulates transport across the BBB. Specifically, for a BBB space constant (λBBB) and wall thickness (dth-BBB), the analytical solution for maximal BBB electric field (EABBB) is given as: (ĒBRAIN × λBBB)/dth-BBB. Electrical current in the brain parenchyma is distorted around brain capillaries, amplifying neuronal polarization. Specifically, capillary ultrastructure produces ∼50% modulation of the ĒBRAIN over the ∼40 μm inter-capillary distance. The divergence of EBRAIN (Activating function) is thus ∼100 kV/m2 per unit ĒBRAIN.

CONCLUSIONS

BBB stimulation by principle 1 suggests novel therapeutic strategies such as boosting metabolic capacity or interstitial fluid clearance. Whereas the spatial profile of EBRAIN is traditionally assumed to depend only on macroscopic anatomy, principle 2 suggests a central role for local capillary ultrastructure-which impact forms of neuromodulation including deep brain stimulation (DBS), spinal cord stimulation (SCS), transcranial magnetic stimulation (TMS), electroconvulsive therapy (ECT), and transcranial electrical stimulation (tES)/transcranial direct current stimulation (tDCS).

Validating layer-specific VASO across species

Cerebral blood volume (CBV) has been shown to be a robust and important physiological parameter for quantitative interpretation of functional (f)MRI, capable of delivering highly localized mapping of neural activity. Indeed, with recent advances in ultra-high-field (≥7T) MRI hardware and associated sequence libraries, it has become possible to capture non-invasive CBV weighted fMRI signals across cortical layers. One of the most widely used approaches to achieve this (in humans) is through vascular-space-occupancy (VASO) fMRI. Unfortunately, the exact contrast mechanisms of layer-dependent VASO fMRI have not been validated for human fMRI and thus interpretation of such data is confounded. Here we validate the signal source of layer-dependent SS-SI VASO fMRI using multi-modal imaging in a rat model in response to neuronal activation (somatosensory cortex) and respiratory challenge (hypercapnia). In particular VASO derived CBV measures are directly compared to concurrent measures of total haemoglobin changes from high resolution intrinsic optical imaging spectroscopy (OIS). Quantified cortical layer profiling is demonstrated to be in agreement between VASO and contrast enhanced fMRI (using monocrystalline iron oxide nanoparticles, MION). Responses show high spatial localisation to layers of cortical processing independent of confounding large draining veins which can hamper BOLD fMRI studies, (depending on slice positioning). Thus, a cross species comparison is enabled using VASO as a common measure. We find increased VASO based CBV reactivity (3.1 ± 1.2 fold increase) in humans compared to rats. Together, our findings confirm that the VASO contrast is indeed a reliable estimate of layer-specific CBV changes. This validation study increases the neuronal interpretability of human layer-dependent VASO fMRI as an appropriate method in neuroscience application studies, in which the presence of large draining intracortical and pial veins limits neuroscientific inference with BOLD fMRI.

Quantitative relations between BOLD responses, cortical energetics and impulse firing across cortical depth

The blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) signal arises as a consequence of changes in cerebral blood flow (CBF) and cerebral metabolic rate of oxygen ( CMR O 2 ) that in turn are modulated by changes in neural activity. Recent advances in imaging have achieved sub-millimetre resolution and allowed investigation of the BOLD response as a function of cortical depth. Here, we adapt our previous theory relating the BOLD signal to neural activity to produce a quantitative model that incorporates venous blood draining between cortical layers. The adjustable inputs to the model are the neural activity and a parameter governing this blood draining. A three-layer version for transient neural inputs and a multi-layer version for constant or tonic neural inputs are able to account for a variety of experimental results, including negative BOLD signals.

Model-Predicted Balance Between Neural Excitation and Inhibition Was Maintained Despite of Age-Related Decline in Sensory Evoked Local Field Potential in Rat Barrel Cortex

The balance between neural excitation and inhibition has been shown to be crucial for normal brain function. However, it is unclear whether this balance is maintained through healthy aging. This study investigated the effect of aging on the temporal dynamics of the somatosensory evoked local field potential (LFP) in rats and tested the hypothesis that excitatory and inhibitory post-synaptic activities remain balanced during the aging process. The LFP signal was obtained from the barrel cortex of three different age groups of anesthetized rats (pre-adolescence: 4–6 weeks, young adult: 2–3 months, middle-aged adult: 10–20 months) under whisker pad stimulation. To confirm our previous finding that the initial segment of the evoked LFP was solely associated with excitatory post-synaptic activity, we micro-injected gabazine into the barrel cortex to block inhibition while LFP was collected continuously under the same stimulus condition. As expected, the initial slope of the evoked LFP in the granular layer was unaffected by gabazine injection. We subsequently estimated the excitatory and inhibitory post-synaptic activities through a balanced model of the LFP with delayed inhibition as an explicit constraint, and calculated the amplitude ratio of inhibition to excitation. We found an age-dependent slowing of the temporal dynamics in the somatosensory-evoked post-synaptic activity, as well as a significant age-related decrease in the amplitude of the excitatory component and a decreasing trend in the amplitude of the inhibitory component. Furthermore, the delay of inhibition with respect to excitation was significantly increased with age, but the amplitude ratio was maintained. Our findings suggest that aging reduces the amplitude of neural responses, but the balance between sensory evoked excitatory and inhibitory post-synaptic activities is maintained to support normal brain function during healthy aging. Further whole cell patch clamp experiments will be needed to confirm or refute these findings by measuring sensory evoked synaptic excitatory and inhibitory activities in vivo during the normal aging process.

The Relationship Between Local Field Potentials and the Blood-Oxygenation-Level Dependent MRI Signal Can Be Non-linear

Functional magnetic resonance imaging (fMRI) is currently one of the most important neuroimaging methods in neuroscience. The image contrast in fMRI relies on the blood-oxygenation-level dependent (BOLD) signal, which indirectly reflects neural activity through neurovascular coupling. Because the mechanism that links the BOLD signal to neural activities involves multiple complicated processes, where neural activity, regional metabolism, hemodynamics, and the BOLD signal are all inter-connected, understanding the quantitative relationship between the BOLD signal and the underlying neural activities is crucial for interpreting fMRI data. Simultaneous local field potential (LFP) and fMRI recordings provide a method to study neurovascular coupling. There were a few studies that have shown non-linearities in stimulus related responses, but whether there is any non-linearity in LFP—BOLD relationship at rest has not been specifically quantified. In this study, we analyzed the simultaneous LFP and resting state-fMRI data acquired from rodents, and found that the relationship between LFP and BOLD is non-linear under isoflurane (ISO) anesthesia, but linear under dexmedetomidine (DMED) anesthesia. Subsequent analysis suggests that such non-linearity may come from the non-Gaussian distribution of LFP power and switching from LFP power to LFP amplitude can alleviate the problem to a degree. We also confirmed that, despite the non-linearity in the mean LFP—BOLD curve, the Pearson correlation between the two signals is relatively unaffected.

Supplemental Vitamin B-12 Enhances the Neural Response to Sensory Stimulation in the Barrel Cortex of Healthy Rats but Does Not Affect Spontaneous Neural Activity

BACKGROUND

Although vitamin B-12 (B-12) is known to contribute to the structural and functional development of the brain, it is unclear if B-12 supplementation has any beneficial effect in healthy populations in terms of enhanced neurologic status of the brain or improved cognitive function.

OBJECTIVES

We investigated the effect of dietary supplementation of B-12 on the cortical neural activity of well-nourished young adult rats and tested the hypothesis that B-12 supplementation in healthy rats may reduce sensory-evoked neural activity due to enhanced inhibition.

METHODS

Female Lister Hooded rats weighing 190-265 g (2-4 mo old) were included in the study. The experimental group was fed with B-12 (cyanocobalamin)-enriched water at a concentration of 1 mg/L, and the control (CON) group with tap water for 3 wk. Animals were then anesthetized and cortical neural responses to whisker stimulation were recorded in vivo through the use of a multichannel microelectrode, from which local field potentials (LFPs) were extracted.

RESULTS

Somatosensory-evoked LFP was 25% larger in the B-12 group (4.13 ± 0.24 mV) than in the CON group (3.30 ± 0.21 mV) (P = 0.02). Spontaneous neural activity did not differ between groups; frequency spectra at each frequency bin of interest did not pass the cluster-forming threshold at the 5% significance level.

CONCLUSIONS

These findings do not provide evidence supporting the hypothesis of decreased neural activity due to B-12 supplementation. As the spontaneous neural activity was unaffected, the increase in somatosensory-evoked LFP may be due to enhanced afferent signal reaching the barrel cortex from the whisker pad, indicating that B-12-supplemented rats may have enhanced sensitivity to sensory stimulation compared with the CON group. We suggest that this enhancement might be the result of lowered sensory threshold, although the underlying mechanism has yet to be elucidated.

Cortical functional reorganization in response to intact forelimb stimulation from acute to chronic stage in rodent amputation model

Brain plasticity after amputation is related to the short-term unmasking of latent synapses as well as the long-term reorganization due to the sprouting new synaptic connections. The cortical functional reorganization has been reported along the intact somatosensory pathway after unilateral deafferentation. Cerebral blood flow (CBF) change serves as an important biomarker of the functional reorganization of brain. Using laser speckle contrast imaging (LSCI) technology, we performed a longitudinal study to unveil the cortical functional reorganization after forelimb amputation in rodent model, particularly along the intact somatosensory pathway. Our results showed that the CBF response to electrical stimulation of the intact forepaw increased significantly at 9 hours after amputation in acute stage. While in chronic stage (>14 days), the CBF response showed a pattern similar to the control group. The results showed the dynamic brain functional response along the intact somatosensory pathway at different stages after amputation and indicated that cortical functional reorganization occurred within the acute stage. Our work provided additional insights in understanding the inter-hemispheric functional changes from acute to chronic stages of amputation.

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.

Wide-field optical mapping of neural activity and brain haemodynamics: considerations and novel approaches

Although modern techniques such as two-photon microscopy can now provide cellular-level three-dimensional imaging of the intact living brain, the speed and fields of view of these techniques remain limited. Conversely, two-dimensional wide-field optical mapping (WFOM), a simpler technique that uses a camera to observe large areas of the exposed cortex under visible light, can detect changes in both neural activity and haemodynamics at very high speeds. Although WFOM may not provide single-neuron or capillary-level resolution, it is an attractive and accessible approach to imaging large areas of the brain in awake, behaving mammals at speeds fast enough to observe widespread neural firing events, as well as their dynamic coupling to haemodynamics. Although such wide-field optical imaging techniques have a long history, the advent of genetically encoded fluorophores that can report neural activity with high sensitivity, as well as modern technologies such as light emitting diodes and sensitive and high-speed digital cameras have driven renewed interest in WFOM. To facilitate the wider adoption and standardization of WFOM approaches for neuroscience and neurovascular coupling research, we provide here an overview of the basic principles of WFOM, considerations for implementation of wide-field fluorescence imaging of neural activity, spectroscopic analysis and interpretation of results.

This article is part of the themed issue ‘Interpreting BOLD: a dialogue between cognitive and cellular neuroscience’.

Impact of Altered Cholinergic Tones on the Neurovascular Coupling Response to Whisker Stimulation

Brain imaging techniques that use vascular signals to map changes in neuronal activity rely on the coupling between electrophysiology and hemodynamics, a phenomenon referred to as "neurovascular coupling" (NVC). It is unknown whether this relationship remains reliable under altered brain states associated with acetylcholine (ACh) levels, such as attention and arousal and in pathological conditions such as Alzheimer's disease. We therefore assessed the effects of varying ACh tone on whisker-evoked NVC responses in rat barrel cortex, measured by cerebral blood flow (CBF) and neurophysiological recordings (local field potentials, LFPs). We found that acutely enhanced ACh tone significantly potentiated whisker-evoked CBF responses through muscarinic ACh receptors and concurrently facilitated neuronal responses, as illustrated by increases in the amplitude and power in high frequencies of the evoked LFPs. However, the cellular identity of the activated neuronal network within the responsive barrel was unchanged, as characterized by c-Fos upregulation in pyramidal cells and GABA interneurons coexpressing vasoactive intestinal polypeptide. In contrast, chronic ACh deprivation hindered whisker-evoked CBF responses and the amplitude and power in most frequency bands of the evoked LFPs and reduced the rostrocaudal extent and area of the activated barrel without altering its identity. Correlations between LFP power and CBF, used to estimate NVC, were enhanced under high ACh tone and disturbed significantly by ACh depletion. We conclude that ACh is not only a facilitator but also a prerequisite for the full expression of sensory-evoked NVC responses, indicating that ACh may alter the fidelity of hemodynamic signals in assessing changes in evoked neuronal activity.SIGNIFICANCE STATEMENT Neurovascular coupling, defined as the tight relationship between activated neurons and hemodynamic responses, is a fundamental brain function that underlies hemodynamic-based functional brain imaging techniques. However, the impact of altered brain states on this relationship is largely unknown. We therefore investigated how acetylcholine (ACh), known to drive brain states of attention and arousal and to be deficient in pathologies such as Alzheimer's disease, would alter neurovascular coupling responses to sensory stimulation. Whereas acutely increased ACh enhanced neuronal responses and the resulting hemodynamic signals, chronic loss of cholinergic input resulted in dramatic impairments in both types of sensory-evoked signals. We conclude that ACh is not only a potent modulator but also a requirement for the full expression of sensory-evoked neurovascular coupling responses.

Resting-state hemodynamics are spatiotemporally coupled to synchronized and symmetric neural activity in excitatory neurons

Brain hemodynamics serve as a proxy for neural activity in a range of noninvasive neuroimaging techniques including functional magnetic resonance imaging (fMRI). In resting-state fMRI, hemodynamic fluctuations have been found to exhibit patterns of bilateral synchrony, with correlated regions inferred to have functional connectivity. However, the relationship between resting-state hemodynamics and underlying neural activity has not been well established, making the neural underpinnings of functional connectivity networks unclear. In this study, neural activity and hemodynamics were recorded simultaneously over the bilateral cortex of awake and anesthetized Thy1-GCaMP mice using wide-field optical mapping. Neural activity was visualized via selective expression of the calcium-sensitive fluorophore GCaMP in layer 2/3 and 5 excitatory neurons. Characteristic patterns of resting-state hemodynamics were accompanied by more rapidly changing bilateral patterns of resting-state neural activity. Spatiotemporal hemodynamics could be modeled by convolving this neural activity with hemodynamic response functions derived through both deconvolution and gamma-variate fitting. Simultaneous imaging and electrophysiology confirmed that Thy1-GCaMP signals are well-predicted by multiunit activity. Neurovascular coupling between resting-state neural activity and hemodynamics was robust and fast in awake animals, whereas coupling in urethane-anesthetized animals was slower, and in some cases included lower-frequency (<0.04 Hz) hemodynamic fluctuations that were not well-predicted by local Thy1-GCaMP recordings. These results support that resting-state hemodynamics in the awake and anesthetized brain are coupled to underlying patterns of excitatory neural activity. The patterns of bilaterally-symmetric spontaneous neural activity revealed by wide-field Thy1-GCaMP imaging may depict the neural foundation of functional connectivity networks detected in resting-state fMRI.

The neurogenesis of P1 and N1: A concurrent EEG/LFP study

It is generally recognised that event related potentials (ERPs) of electroencephalogram (EEG) primarily reflect summed post-synaptic activity of the local pyramidal neural population(s). However, it is still not understood how the positive and negative deflections (e.g. P1, N1 etc) observed in ERP recordings are related to the underlying excitatory and inhibitory post-synaptic activity. We investigated the neurogenesis of P1 and N1 in ERPs by pharmacologically manipulating inhibitory post-synaptic activity in the somatosensory cortex of rodent, and concurrently recording EEG and local field potentials (LFPs). We found that the P1 wave in the ERP and LFP of the supragranular layers is determined solely by the excitatory post-synaptic activity of the local pyramidal neural population, as is the initial segment of the N1 wave across cortical depth. The later part of the N1 wave was modulated by inhibitory post-synaptic activity, with its peak and the pulse width increasing as inhibition was reduced. These findings suggest that the temporal delay of inhibition with respect to excitation observed in intracellular recordings is also reflected in extracellular field potentials (FPs), resulting in a temporal window during which only excitatory post-synaptic activity and leak channel activity are recorded in the ERP and evoked LFP time series. Based on these findings, we provide clarification on the interpretation of P1 and N1 in terms of the excitatory and inhibitory post-synaptic activities of the local pyramidal neural population(s).

Encoding of graded changes in spatial specificity of prior cues in human visual cortex

Prior information about the relevance of spatial locations can vary in specificity; a single location, a subset of locations, or all locations may be of potential importance. Using a contrast-discrimination task with four possible targets, we asked whether performance benefits are graded with the spatial specificity of a prior cue and whether we could quantitatively account for behavioral performance with cortical activity changes measured by blood oxygenation level-dependent (BOLD) imaging. Thus we changed the prior probability that each location contained the target from 100 to 50 to 25% by cueing in advance 1, 2, or 4 of the possible locations. We found that behavioral performance (discrimination thresholds) improved in a graded fashion with spatial specificity. However, concurrently measured cortical responses from retinotopically defined visual areas were not strictly graded; response magnitude decreased when all 4 locations were cued (25% prior probability) relative to the 100 and 50% prior probability conditions, but no significant difference in response magnitude was found between the 100 and 50% prior probability conditions for either cued or uncued locations. Also, although cueing locations increased responses relative to noncueing, this cue sensitivity was not graded with prior probability. Furthermore, contrast sensitivity of cortical responses, which could improve contrast discrimination performance, was not graded. Instead, an efficient-selection model showed that even if sensory responses do not strictly scale with prior probability, selection of sensory responses by weighting larger responses more can result in graded behavioral performance benefits with increasing spatial specificity of prior information.

Contributions and complexities from the use of in vivo animal models to improve understanding of human neuroimaging signals

Many of the major advances in our understanding of how functional brain imaging signals relate to neuronal activity over the previous two decades have arisen from physiological research studies involving experimental animal models. This approach has been successful partly because it provides opportunities to measure both the hemodynamic changes that underpin many human functional brain imaging techniques and the neuronal activity about which we wish to make inferences. Although research into the coupling of neuronal and hemodynamic responses using animal models has provided a general validation of the correspondence of neuroimaging signals to specific types of neuronal activity, it is also highlighting the key complexities and uncertainties in estimating neural signals from hemodynamic markers. This review will detail how research in animal models is contributing to our rapidly evolving understanding of what human neuroimaging techniques tell us about neuronal activity. It will highlight emerging issues in the interpretation of neuroimaging data that arise from in vivo research studies, for example spatial and temporal constraints to neuroimaging signal interpretation, or the effects of disease and modulatory neurotransmitters upon neurovascular coupling. We will also give critical consideration to the limitations and possible complexities of translating data acquired in the typical animals models used in this area to the arena of human fMRI. These include the commonplace use of anesthesia in animal research studies and the fact that many neuropsychological questions that are being actively explored in humans have limited homologs within current animal models for neuroimaging research. Finally we will highlighting approaches, both in experimental animals models (e.g. imaging in conscious, behaving animals) and human studies (e.g. combined fMRI-EEG), that mitigate against these challenges.

Phase-amplitude coupling and infraslow (<1 Hz) frequencies in the rat brain: relationship to resting state fMRI

Resting state functional magnetic resonance imaging (fMRI) can identify network alterations that occur in complex psychiatric diseases and behaviors, but its interpretation is difficult because the neural basis of the infraslow BOLD fluctuations is poorly understood. Previous results link dynamic activity during the resting state to both infraslow frequencies in local field potentials (LFP) (<1 Hz) and band-limited power in higher frequency LFP (>1 Hz). To investigate the relationship between these frequencies, LFPs were recorded from rats under two anesthetics: isoflurane and dexmedetomidine. Signal phases were calculated from low-frequency LFP and compared to signal amplitudes from high-frequency LFP to determine if modulation existed between the two frequency bands (phase-amplitude coupling). Isoflurane showed significant, consistent phase-amplitude coupling at nearly all pairs of frequencies, likely due to the burst-suppression pattern of activity that it induces. However, no consistent phase-amplitude coupling was observed in rats that were anesthetized with dexmedetomidine. fMRI-LFP correlations under isoflurane using high frequency LFP were reduced when the low frequency LFP's influence was accounted for, but not vice-versa, or in any condition under dexmedetomidine. The lack of consistent phase-amplitude coupling under dexmedetomidine and lack of shared variance between high frequency and low frequency LFP as it relates to fMRI suggests that high and low frequency neural electrical signals may contribute differently, possibly even independently, to resting state fMRI. This finding suggests that researchers take care in interpreting the neural basis of resting state fMRI, as multiple dynamic factors in the underlying electrophysiology could be driving any particular observation.

Neurovascular coupling in normal aging: a combined optical, ERP and fMRI study

Brain aging is characterized by changes in both hemodynamic and neuronal responses, which may be influenced by the cardiorespiratory fitness of the individual. To investigate the relationship between neuronal and hemodynamic changes, we studied the brain activity elicited by visual stimulation (checkerboard reversals at different frequencies) in younger adults and in older adults varying in physical fitness. Four functional brain measures were used to compare neuronal and hemodynamic responses obtained from BA17: two reflecting neuronal activity (the event-related optical signal, EROS, and the C1 response of the ERP), and two reflecting functional hemodynamic changes (functional magnetic resonance imaging, fMRI, and near-infrared spectroscopy, NIRS). The results indicated that both younger and older adults exhibited a quadratic relationship between neuronal and hemodynamic effects, with reduced increases of the hemodynamic response at high levels of neuronal activity. Although older adults showed reduced activation, similar neurovascular coupling functions were observed in the two age groups when fMRI and deoxy-hemoglobin measures were used. However, the coupling between oxy- and deoxy-hemoglobin changes decreased with age and increased with increasing fitness. These data indicate that departures from linearity in neurovascular coupling may be present when using hemodynamic measures to study neuronal function.

The effects of focal epileptic activity on regional sensory-evoked neurovascular coupling and postictal modulation of bilateral sensory processing

While it is known that cortical sensory dysfunction may occur in focal neocortical epilepsy, it is unknown whether sensory-evoked neurovascular coupling is also disrupted during epileptiform activity. Addressing this open question may help to elucidate both the effects of focal neocortical epilepsy on sensory responses and the neurovascular characteristics of epileptogenic regions in sensory cortex. We therefore examined bilateral sensory-evoked neurovascular responses before, during, and after 4-aminopyridine (4-AP, 15 mmol/L, 1 μL) induced focal neocortical seizures in right vibrissal cortex of the rat. Stimulation consisted of electrical pulse trains (16 seconds, 5 Hz, 1.2 mA) presented to the mystacial pad. Consequent current-source density neural responses and epileptic activity in both cortices and across laminae were recorded via two 16-channel microelectrodes bilaterally implanted in vibrissal cortices. Concurrent two-dimensional optical imaging spectroscopy was used to produce spatiotemporal maps of total, oxy-, and deoxy-hemoglobin concentration. Compared with control, sensory-evoked neurovascular coupling was altered during ictal activity, but conserved postictally in both ipsilateral and contralateral vibrissal cortices, despite neurovascular responses being significantly reduced in the former, and enhanced in the latter. Our results provide insights into sensory-evoked neurovascular dynamics and coupling in epilepsy, and may have implications for the localization of epileptogenic foci and neighboring eloquent cortex.

Hemodynamic and electrophysiological spontaneous low-frequency oscillations in the cortex: directional influences revealed by Granger causality

We used a combined electrophysiological/hemodynamic system to examine low-frequency oscillations (LFOs) in spontaneous neuronal activities (spike trains and local field potentials) and hemodynamic signals (cerebral blood flow) recorded from the anesthetized rat somatosensory and visual cortices. The laser Doppler flowmetry (LDF) probe was tilted slightly to approach the area in which a microelectrode array (MEA) was implanted for simultaneous recordings. Spike trains (STs) were converted into continuous-time rate functions (CRFs) using the ST instantaneous firing rates. LFOs were detected for all three of the components using the multi-taper method (MTM). The frequencies of these LFOs ranged from 0.052 to 0.167 Hz (mean±SD, 0.10±0.026 Hz) for cerebral blood flow (CBF), from 0.027 to 0.26 Hz (mean±SD, 0.12±0.041 Hz) for the CRFs of the STs and from 0.04 to 0.19 Hz (mean±SD, 0.11±0.035 Hz) for local field potentials (LFPs). We evaluated the Granger causal relationships of spontaneous LFOs among CBF, LFPs and CRFs using Granger causality (GC) analysis. Significant Granger causal relationships were observed from LFPs to CBF, from STs to CBF and from LFPs to STs at approximately 0.1 Hz. The present results indicate that spontaneous LFOs exist not only in hemodynamic components but also in neuronal activities of the rat cortex. To the best of our knowledge, the present study is the first to identify Granger causal influences among CBF, LFPs and STs and show that spontaneous LFOs carry important Granger causal influences from neural activities to hemodynamic signals.

Binocular activation elicits differences in neurometabolic coupling in visual cortex

Non-invasive brain imaging requires comprehensive interpretation of hemodynamic signals. In functional magnetic resonance imaging, blood oxygen level dependent (BOLD) signals are used to infer neural processes. This necessitates a clear understanding of how BOLD signals and neural activity are related. One fundamental question concerns the relative importance of synaptic activity and spiking discharge. Although these two components are related, most previous work shows that synaptic activity is better reflected in the BOLD signal. However, the mechanisms of this relationship are not clear. The BOLD signal depends on relative changes in cerebral blood flow and cerebral metabolic rate of oxygen. Oxygen metabolism changes are difficult to measure with current imaging techniques, but it is possible to obtain direct quantitative simultaneous in vivo measurement of tissue oxygen and co-localized underlying neural activity. Here, we use this approach with a specific binocular stimulus protocol in order to activate inhibitory and excitatory neuronal pathways in the visual cortex. During excitatory binocular interaction, we find that metabolic, spiking, and local field potential responses are correlated. However, during suppressive binocular interaction, spiking activity and local field potentials (LFP) are dissociated while only the latter is coupled with metabolic response. These results suggest that inhibitory connections may be a key factor in the dissociation between LFP and spiking activity, which may contribute substantially to the close coupling between the BOLD signal and synchronized synaptic activity in the brain.

Ultra high-resolution fMRI and electrophysiology of the rat primary somatosensory cortex

High-resolution functional-magnetic-resonance-imaging (fMRI) has been used to study brain functions at increasingly finer scale, but whether fMRI can accurately reflect layer-specific neuronal activities is less well understood. The present study investigated layer-specific cerebral-blood-volume (CBV) fMRI and electrophysiological responses in the rat cortex. CBV fMRI at 40×40 μm in-plane resolution was performed on an 11.7-T scanner. Electrophysiology used a 32-channel electrode array that spanned the entire cortical depth. Graded electrical stimulation was used to study activations in different cortical layers, exploiting the notion that most of the sensory-specific neurons are in layers II-V and most of the nociceptive-specific neurons are in layers V-VI. CBV response was strongest in layer IV of all stimulus amplitudes. Current source density analysis showed strong sink currents at cortical layers IV and VI. Multi-unit activities mainly appeared at layers IV-VI and peaked at layer V. Although our measures showed scaled activation profiles during modulation of stimulus amplitude and failed to detect specific recruitment at layers V and VI during noxious electrical stimuli, there appears to be discordance between CBV fMRI and electrophysiological peak responses, suggesting neurovascular uncoupling at laminar resolution. The technique implemented in the present study offers a means to investigate intracortical neurovascular function in the normal and diseased animal models at laminar resolution.

Inter-trial variability in sensory-evoked cortical hemodynamic responses: the role of the magnitude of pre-stimulus fluctuations

Brain imaging techniques utilize hemodynamic changes that accompany brain activation. However, stimulus-evoked hemodynamic responses display considerable inter-trial variability and the sources of this variability are poorly understood. One of the sources of this response variation could be ongoing spontaneous hemodynamic fluctuations. We recently investigated this issue by measuring cortical hemodynamics in response to sensory stimuli in anesthetized rodents using 2-dimensional optical imaging spectroscopy. We suggested that sensory-evoked cortical hemodynamics displayed distinctive response characteristics and magnitudes depending on the phase of ongoing fluctuations at stimulus onset due to a linear superposition of evoked and ongoing hemodynamics (Saka et al., 2010). However, the previous analysis neglected to examine the possible influence of variability of the size of ongoing fluctuations. Consequently, data were further analyzed to examine whether the size of pre-stimulus hemodynamic fluctuations also influenced the magnitude of subsequent stimulus-evoked responses. Indeed, in the case of all individual trials, a moderate correlation between the size of the pre-stimulus fluctuations and the magnitudes of the subsequent sensory-evoked responses were observed. However, different correlations between the size of the pre-stimulus fluctuations and magnitudes of the subsequent sensory-evoked cortical hemodynamic responses could be observed depending on their phase at stimulus onset. These analyses suggest that both the size and phase of pre-stimulus fluctuations in cortical hemodynamics contribute to inter-trial variability in sensory-evoked responses.

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.

Calibrating the BOLD signal during a motor task using an extended fusion model incorporating DOT, BOLD and ASL data

Multimodal imaging improves the accuracy of the localization and the quantification of brain activation when measuring different manifestations of the hemodynamic response associated with cerebral activity. In this study, we incorporated cerebral blood flow (CBF) changes measured with arterial spin labeling (ASL), Diffuse Optical Tomography (DOT) and blood oxygen level-dependent (BOLD) recordings to reconstruct changes in oxy- (ΔHbO(2)) and deoxyhemoglobin (ΔHbR). Using the Grubb relation between relative changes in CBF and cerebral blood volume (CBV), we incorporated the ASL measurement as a prior to the total hemoglobin concentration change (ΔHbT). We applied this ASL fusion model to both synthetic data and experimental multimodal recordings during a 2-s finger-tapping task. Our results show that the new approach is very powerful in estimating ΔHbO(2) and ΔHbR with high spatial and quantitative accuracy. Moreover, our approach allows the computation of baseline total hemoglobin concentration (HbT(0)) as well as of the BOLD calibration factor M on a single subject basis. We obtained an average HbT(0) of 71 μM, an average M value of 0.18 and an average increase of 13% in cerebral metabolic rate of oxygen (CMRO(2)), all of which are in agreement with values previously reported in the literature. Our method yields an independent measurement of M, which provides an alternative measurement to validate the hypercapnic calibration of the BOLD signal.

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.

Investigaton of the neuronal efficacy and EEG source power under steady-state visual stimulation

Understanding the nature of the link between neuronal activity and BOLD signal plays a crucial role i) for improving the interpretability of BOLD images and ii) on the design of more realistic models for the integration of EEG and fMRI. The aim of this study is to investigate the neural mechanism underlying hemodynamic behavior in a series of visual stimulation frequencies and explore possible implications for the neurovascular coupling. We studied the relationship between electrophysiological and hemodynamic measures by performing simultaneous steady state electroencephalography (EEG) and fMRI recordings in a healthy human subject during a series of visual stimulation frequencies (6 Hz, 8 Hz, 10 Hz, 12 Hz). BOLD amplitudes were computed for voxels within an anatomical mask which was obtained by mapping the significantly active voxels using general linear modelling (GLM) on fMRI data. On the same anatomical map, EEG power time series belonging to the fundamental frequency and its harmonics due to the stimulation are estimated using a distributed source imaging technique. The neuronal efficacies which represent the vascular inputs driving the BOLD response are estimated by use of an extended version of Balloon model. A nonlinear relationship is demonstrated between the mean EEG source powers and the neuronal efficacies driving the BOLD response. The result suggests that BOLD signal which is an indicator of the metabolic demand of both synchronized and non-synchronized neuronal activities; changes independent of EEG activity which is a measure sensitive to the synchronicity of neuronal activity.

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.

Degeneration of cortical function in the Royal College of Surgeons rat

The purpose of the current study was to determine the progress of cortical functional degeneration in the Royal College of Surgeons (RCS) rat. Cortical responses were measured with optical imaging of intrinsic signals using gratings of various spatial frequencies. Subsequently, electrophysiological recordings were also taken across cortical layers in response to a pulse of broad-spectrum light. We found significant degeneration in the cortical processing of visual information as early as 4 weeks of age. These results show that degeneration in the cortical response of the RCS rat starts before development has been properly completed.

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.

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.

High-speed vascular dynamics of the hemodynamic response

While a range of cellular mechanisms have been proposed to underlie control of neurovascular coupling, a comprehensive, reconciliatory model has yet to be determined. To fit with such a model, it is essential that candidate mechanisms exhibit reaction times, spatial ranges, and speeds of propagation that are consistent with the vascular manifestations of the 'hemodynamic response'. Understanding these vascular dynamics is therefore a critical step towards developing a robust model of neurovascular coupling. In this study, we utilize high-speed optical imaging of exposed rodent somatosensory cortex to explore and characterize the spatiotemporal dynamics of surface vessels during functional hyperemia. Our high-speed, high-resolution optical imaging approach allows us to study the hemodynamic response independently in individual vessels, and in discrete regions of the parenchyma with enough resolution to precisely characterize subtle spatial and temporal features of the response. Specifically, we explore when and where the first hemodynamic changes occur in response to stimuli, the direction and speed at which these changes propagate in arterioles and regions of the parenchyma, and the relative timing at which each of these compartments returns to its original baseline state. From these results, we are able to conclude that the hemodynamic response appears to initiate in the parenchyma and then spreads rapidly to surface arterioles. Following the initial onset we find evidence that the response spreads spatially outwards via the dilation of targeted arterioles. This propagation of vasodilation is independent of the direction of blood flow within each arteriole. We also find evidence of a decay phase that acts with a more uniform spatial dependence, rather than along targeted vessels, causing the periphery of the responding region to return to baseline first. We hypothesize that different underlying cellular mechanisms/signaling pathways are responsible for the response initiation and the response decay. Our results advance the fundamental understanding of the hemodynamic response, as well as our ability to evaluate potential cellular mechanisms for their involvement in neurovascular coupling.

Linear superposition of sensory-evoked and ongoing cortical hemodynamics

Modern non-invasive brain imaging techniques utilize changes in cerebral blood flow, volume and oxygenation that accompany brain activation. However, stimulus-evoked hemodynamic responses display considerable inter-trial variability even when identical stimuli are presented and the sources of this variability are poorly understood. One of the sources of this response variation could be ongoing spontaneous hemodynamic fluctuations. To investigate this issue, 2-dimensional optical imaging spectroscopy was used to measure cortical hemodynamics in response to sensory stimuli in anesthetized rodents. Pre-stimulus cortical hemodynamics displayed spontaneous periodic fluctuations and as such, data from individual stimulus presentation trials were assigned to one of four groups depending on the phase angle of pre-stimulus hemodynamic fluctuations and averaged. This analysis revealed that sensory evoked cortical hemodynamics displayed distinctive response characteristics and magnitudes depending on the phase angle of ongoing fluctuations at stimulus onset. To investigate the origin of this phenomenon, "null-trials" were collected without stimulus presentation. Subtraction of phase averaged "null trials" from their phase averaged stimulus-evoked counterparts resulted in four similar time series that resembled the mean stimulus-evoked response. These analyses suggest that linear superposition of evoked and ongoing cortical hemodynamic changes may be a property of the structure of inter-trial variability.

Does neural input or processing play a greater role in the magnitude of neuroimaging signals?

An important constraint on how hemodynamic neuroimaging signals such as fMRI can be interpreted in terms of the underlying evoked activity is an understanding of neurovascular coupling mechanisms that actually generate hemodynamic responses. The predominant view at present is that the hemodynamic response is most correlated with synaptic input and subsequent neural processing rather than spiking output. It is still not clear whether input or processing is more important in the generation of hemodynamics responses. In order to investigate this we measured the hemodynamic and neural responses to electrical whisker pad stimuli in rat whisker barrel somatosensory cortex both before and after the local cortical injections of the GABA(A) agonist muscimol. Muscimol would not be expected to affect the thalamocortical input into the cortex but would inhibit subsequent intra-cortical processing. Pre-muscimol infusion whisker stimuli elicited the expected neural and accompanying hemodynamic responses to that reported previously. Following infusion of muscimol, although the temporal profile of neural responses to each pulse of the stimulus train was similar, the average response was reduced in magnitude by approximately 79% compared to that elicited pre-infusion. The whisker-evoked hemodynamic responses were reduced by a commensurate magnitude suggesting that, although the neurovascular coupling relationships were similar for synaptic input as well as for cortical processing, the magnitude of the overall response is dominated by processing rather than from that produced from the thalamocortical input alone.

Investigation of Linear Coupling Between Single-Event Blood Flow Responses and Interictal Discharges in a Model of Experimental Epilepsy

A successful outcome of epilepsy neurosurgery relies on an accurate delineation of the epileptogenic region to be resected. Functional magnetic resonance imaging (fMRI) would allow doing this noninvasively at high spatial resolution. However, a clear, quantitative description of the relationship between hemodynamic changes and the underlying epileptiform neuronal activity is still missing, thereby preventing the systematic use of fMRI for routine epilepsy surgery planning. To this aim, we used a local epilepsy model to record simultaneously cerebral blood flow (CBF) with laser Doppler (LD) and local field potentials (LFP) in rat frontal cortex. CBF responses to individual interictal-like spikes were large and robust. Their amplitude correlated linearly with spike amplitude. Moreover, the CBF response added linearly in time over a large range of spiking rates. CBF responses could thus be predicted by a linear model of the kind currently used for the interpretation of fMRI data, but including also the spikes’ amplitudes as additional information. Predicted and measured CBF responses matched accurately. For high spiking frequencies (above ∼0.2 Hz), the responses saturated but could eventually recover, indicating the presence of multiple neurovascular coupling mechanisms, which might act at different spatiotemporal scales. Spatially, CBF responses peaked at the center of epileptic activity and displayed a spatial specificity at least as good as the millimeter. These results suggest that simultaneous electroencephalographic and blood flow-based fMRI recordings should be suitable for the noninvasive precise localization of hyperexcitable regions in epileptic patients candidate for neurosurgery.

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.

A Dynamic Causal Model of the Coupling Between Pulse Stimulation and Neural Activity

We present a dynamic causal model that can explain context-dependent changes in neural responses, in the rat barrel cortex, to an electrical whisker stimulation at different frequencies. Neural responses were measured in terms of local field potentials. These were converted into current source density (CSD) data, and the time series of the CSD sink was extracted to provide a time series response train. The model structure consists of three layers (approximating the responses from the brain stem to the thalamus and then the barrel cortex), and the latter two layers contain nonlinearly coupled modules of linear second-order dynamic systems. The interaction of these modules forms a nonlinear regulatory system that determines the temporal structure of the neural response amplitude for the thalamic and cortical layers.

The model is based on the measured population dynamics of neurons rather than the dynamics of a single neuron and was evaluated against CSD data from experiments with varying stimulation frequency (1–40 Hz), random pulse trains, and awake and anesthetized animals. The model parameters obtained by optimization for different physiological conditions (anesthetized or awake) were significantly different.

Following Friston, Mechelli, Turner, and Price ( 2000 ), this work is part of a formal mathematical system currently being developed (Zheng et al., 2005 ) that links stimulation to the blood oxygen level dependent (BOLD) functional magnetic resonance imaging (fMRI) signal through neural activity and hemodynamic variables. The importance of the model described here is that it can be used to invert the hemodynamic measurements of changes in blood flow to estimate the underlying neural activity.

Current trends in intraoperative optical imaging for functional brain mapping and delineation of lesions of language cortex

Resection of a cerebral arteriovenous malformation (AVM), epileptic focus, or glioma, ideally has a prerequisite of microscopic delineation of the lesion borders in relation to the normal gray and white matter that mediate critical functions. Currently, Wada testing and functional magnetic resonance imaging (fMRI) are used for preoperative mapping of critical function, whereas electrical stimulation mapping (ESM) is used for intraoperative mapping. For lesion delineation, MRI and positron emission tomography (PET) are used preoperatively, whereas microscopy and histological sectioning are used intraoperatively. However, for lesions near eloquent cortex, these imaging techniques may lack sufficient resolution to define the relationship between the lesion and language function, and thus not accurately determine which patients will benefit from neurosurgical resection of the lesion without iatrogenic aphasia. Optical techniques such as intraoperative optical imaging of intrinsic signals (iOIS) show great promise for the precise functional mapping of cortices, as well as delineation of the borders of AVMs, epileptic foci, and gliomas. Here we first review the physiology of neuroimaging, and then progress towards the validation and justification of using intraoperative optical techniques, especially in relation to neurosurgical planning of resection AVMs, epileptic foci, and gliomas near or in eloquent cortex. We conclude with a short description of potential novel intraoperative optical techniques.