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

Neuronal networks and mediators of cortical neurovascular coupling responses in normal and altered brain states

Brain imaging techniques that use vascular signals to map changes in neuronal activity, such as blood oxygenation level-dependent functional magnetic resonance imaging, rely on the spatial and temporal coupling between changes in neurophysiology and haemodynamics, known as 'neurovascular coupling (NVC)'. Accordingly, NVC responses, mapped by changes in brain haemodynamics, have been validated for different stimuli under physiological conditions. In the cerebral cortex, the networks of excitatory pyramidal cells and inhibitory interneurons generating the changes in neural activity and the key mediators that signal to the vascular unit have been identified for some incoming afferent pathways. The neural circuits recruited by whisker glutamatergic-, basal forebrain cholinergic- or locus coeruleus noradrenergic pathway stimulation were found to be highly specific and discriminative, particularly when comparing the two modulatory systems to the sensory response. However, it is largely unknown whether or not NVC is still reliable when brain states are altered or in disease conditions. This lack of knowledge is surprising since brain imaging is broadly used in humans and, ultimately, in conditions that deviate from baseline brain function. Using the whisker-to-barrel pathway as a model of NVC, we can interrogate the reliability of NVC under enhanced cholinergic or noradrenergic modulation of cortical circuits that alters brain states.This article is part of the themed issue 'Interpreting BOLD: a dialogue between cognitive and cellular neuroscience'.

Coupling mechanism and significance of the BOLD signal: a status report

Functional magnetic resonance imaging (fMRI) provides a unique view of the working human mind. The blood-oxygen-level-dependent (BOLD) signal, detected in fMRI, reflects changes in deoxyhemoglobin driven by localized changes in brain blood flow and blood oxygenation, which are coupled to underlying neuronal activity by a process termed neurovascular coupling. Over the past 10 years, a range of cellular mechanisms, including astrocytes, pericytes, and interneurons, have been proposed to play a role in functional neurovascular coupling. However, the field remains conflicted over the relative importance of each process, while key spatiotemporal features of BOLD response remain unexplained. Here, we review current candidate neurovascular coupling mechanisms and propose that previously overlooked involvement of the vascular endothelium may provide a more complete picture of how blood flow is controlled in the brain. We also explore the possibility and consequences of conditions in which neurovascular coupling may be altered, including during postnatal development, pathological states, and aging, noting relevance to both stimulus-evoked and resting-state fMRI studies.

Comprehensive chronic laminar single-unit, multi-unit, and local field potential recording performance with planar single shank electrode arrays

BACKGROUND

Intracortical electrode arrays that can record extracellular action potentials from small, targeted groups of neurons are critical for basic neuroscience research and emerging clinical applications. In general, these electrode devices suffer from reliability and variability issues, which have led to comparative studies of existing and emerging electrode designs to optimize performance. Comparisons of different chronic recording devices have been limited to single-unit (SU) activity and employed a bulk averaging approach treating brain architecture as homogeneous with respect to electrode distribution.

NEW METHOD

In this study, we optimize the methods and parameters to quantify evoked multi-unit (MU) and local field potential (LFP) recordings in eight mice visual cortices.

RESULTS

These findings quantify the large recording differences stemming from anatomical differences in depth and the layer dependent relative changes to SU and MU recording performance over 6-months. For example, performance metrics in Layer V and stratum pyramidale were initially higher than Layer II/III, but decrease more rapidly. On the other hand, Layer II/III maintained recording metrics longer. In addition, chronic changes at the level of layer IV are evaluated using visually evoked current source density.

COMPARISON WITH EXISTING METHOD(S)

The use of MU and LFP activity for evaluation and tracking biological depth provides a more comprehensive characterization of the electrophysiological performance landscape of microelectrodes.

CONCLUSIONS

A more extensive spatial and temporal insight into the chronic electrophysiological performance over time will help uncover the biological and mechanical failure mechanisms of the neural electrodes and direct future research toward the elucidation of design optimization for specific applications.

GABAA Receptor-Mediated Bidirectional Control of Synaptic Activity, Intracellular Ca2+, Cerebral Blood Flow, and Oxygen Consumption in Mouse Somatosensory Cortex In Vivo

Neural activity regulates local increases in cerebral blood flow (ΔCBF) and the cortical metabolic rate of oxygen (ΔCMRO2) that constitutes the basis of BOLD functional neuroimaging signals. Glutamate signaling plays a key role in brain vascular and metabolic control; however, the modulatory effect of GABA is incompletely understood. Here we performed in vivo studies in mice to investigate how THIP (which tonically activates extrasynaptic GABAARs) and Zolpidem (a positive allosteric modulator of synaptic GABAARs) impact stimulation-induced ΔCBF, ΔCMRO2, local field potentials (LFPs), and fluorescent cytosolic Ca(2+) transients in neurons and astrocytes. Low concentrations of THIP increased ΔCBF and ΔCMRO2 at low stimulation frequencies. These responses were coupled to increased synaptic activity as indicated by LFP responses, and to Ca(2+) activities in neurons and astrocytes. Intermediate and high concentrations of THIP suppressed ΔCBF and ΔCMRO2 at high stimulation frequencies. Zolpidem had similar but less-pronounced effects, with similar dependence on drug concentration and stimulation frequency. Our present findings suggest that slight increases in both synaptic and extrasynaptic GABAAR activity might selectively gate and amplify transient low-frequency somatosensory inputs, filter out high-frequency inputs, and enhance vascular and metabolic responses that are likely to be reflected in BOLD functional neuroimaging signals.

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.

Optical imaging of visually guided reaching in macaque posterior parietal cortex

Sensorimotor transformation for reaching movements in primates requires a large network of visual, parietal, and frontal cortical areas. We performed intrinsic optical imaging over posterior parietal cortex including areas 7a and the dorsal perilunate in macaque monkeys during visually guided hand movements. Reaching was performed while foveating one of nine static reach targets; thus eye-position-varied concurrently with reach position. The hemodynamic reflectance signal was analyzed during specific phases of the task including pre-reach, reach, and touch epochs. The eye position maps changed substantially as the task progressed: First, direction of spatial tuning shifted from a weak preference close to the center to the lower eye positions in both cortical areas. Overall tuning strength was greater in area 7a. Second, strength of spatial tuning increased from the early pre-reach to the later touch epoch. These consistent temporal changes suggest that dynamic properties of the reflectance signal were modulated by task parameters. The peak amplitude and peak delay of the reflectance signal showed considerable differences between eye position but were similar between areas. Compared with a detection task using a lever response, the reach task yielded higher amplitudes and longer delays. These findings demonstrate a spatially tuned topographical representation for reaching in both areas and suggest a strong synergistic combination of various feedback signals that result in a spatially tuned amplification of the hemodynamic response in posterior parietal cortex.

Anesthesia and the quantitative evaluation of neurovascular coupling

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

Is optical imaging spectroscopy a viable measurement technique for the investigation of the negative BOLD phenomenon? A concurrent optical imaging spectroscopy and fMRI study at high field (7 T)

Traditionally functional magnetic resonance imaging (fMRI) has been used to map activity in the human brain by measuring increases in the Blood Oxygenation Level Dependent (BOLD) signal. Often accompanying positive BOLD fMRI signal changes are sustained negative signal changes. Previous studies investigating the neurovascular coupling mechanisms of the negative BOLD phenomenon have used concurrent 2D-optical imaging spectroscopy (2D-OIS) and electrophysiology (Boorman et al., 2010). These experiments suggested that the negative BOLD signal in response to whisker stimulation was a result of an increase in deoxy-haemoglobin and reduced multi-unit activity in the deep cortical layers. However, Boorman et al. (2010) did not measure the BOLD and haemodynamic response concurrently and so could not quantitatively compare either the spatial maps or the 2D-OIS and fMRI time series directly. Furthermore their study utilised a homogeneous tissue model in which is predominantly sensitive to haemodynamic changes in more superficial layers.

Here we test whether the 2D-OIS technique is appropriate for studies of negative BOLD. We used concurrent fMRI with 2D-OIS techniques for the investigation of the haemodynamics underlying the negative BOLD at 7 Tesla. We investigated whether optical methods could be used to accurately map and measure the negative BOLD phenomenon by using 2D-OIS haemodynamic data to derive predictions from a biophysical model of BOLD signal changes. We showed that despite the deep cortical origin of the negative BOLD response, if an appropriate heterogeneous tissue model is used in the spectroscopic analysis then 2D-OIS can be used to investigate the negative BOLD phenomenon.

NMDA-dependent mechanisms only affect the BOLD response in the rat dentate gyrus by modifying local signal processing

The role of N-methyl-D-aspartate (NMDA) receptor-mediated mechanisms in the formation of a blood oxygen level-dependent (BOLD) response was studied using electrical stimulation of the right perforant pathway. Stimulation of this fiber bundle triggered BOLD responses in the right hippocampal formation and in the left entorhinal cortex. The perforant pathway projects to and activates the dentate gyrus monosynaptically, activation in the contralateral entorhinal cortex is multisynaptic and requires forwarding and processing of signals. Application of the NMDA receptor antagonist MK801 during stimulation had no effect on BOLD responses in the right dentate gyrus, but reduced the BOLD responses in the left entorhinal cortex. In contrast, application of MK801 before the first stimulation train reduced the BOLD response in both regions. Electrophysiological recordings revealed that the initial stimulation trains changed the local processing of the incoming signals in the dentate gyrus. This altered electrophysiological response was not further changed by a subsequent application of MK801, which is in agreement with an unchanged BOLD response. When MK801 was present during the first stimulation train, a dissimilar electrophysiological response pattern was observed and corresponds to an altered BOLD response, indicating that NMDA-dependent mechanisms indirectly affect the BOLD response, mainly via modifying local signal processing and subsequent propagation.

Low frequency stimulation of the perforant pathway generates anesthesia-specific variations in neural activity and BOLD responses in the rat dentate gyrus

To study how various anesthetics affect the relationship between stimulus frequency and generated functional magnetic resonance imaging (fMRI) signals in the rat dentate gyrus, the perforant pathway was electrically stimulated with repetitive low frequency (i.e., 0.625, 1.25, 2.5, 5, and 10 Hz) stimulation trains under isoflurane/N(2)O, isoflurane, medetomidine, and α-chloralose. During stimulation, the blood oxygen level-dependent signal intensity (BOLD response) and local field potentials in the dentate gyrus were simultaneously recorded to prove whether the present anesthetic controls the generation of a BOLD response via targeting general hemodynamic parameters, by affecting mechanisms of neurovascular coupling, or by disrupting local signal processing. Using this combined electrophysiological/fMRI approach, we found that the threshold frequency (i.e., the minimal frequency required to trigger significant BOLD responses), the optimal frequency (i.e., the frequency that elicit the strongest BOLD response), and the spatial distribution of generated BOLD responses are specific for each anesthetic used. Concurrent with anesthetic-dependent characteristics of the BOLD response, we found the pattern of stimulus-induced neuronal activity in the dentate gyrus is also specific for each anesthetic. Consequently, the anesthetic-specific influence on local signaling processes is the underlying cause for the observation that an identical stimulus elicits different BOLD responses under various anesthetics.

Bayesian comparison of neurovascular coupling models using EEG-fMRI

Functional magnetic resonance imaging (fMRI), with blood oxygenation level-dependent (BOLD) contrast, is a widely used technique for studying the human brain. However, it is an indirect measure of underlying neuronal activity and the processes that link this activity to BOLD signals are still a topic of much debate. In order to relate findings from fMRI research to other measures of neuronal activity it is vital to understand the underlying neurovascular coupling mechanism. Currently, there is no consensus on the relative roles of synaptic and spiking activity in the generation of the BOLD response. Here we designed a modelling framework to investigate different neurovascular coupling mechanisms. We use Electroencephalographic (EEG) and fMRI data from a visual stimulation task together with biophysically informed mathematical models describing how neuronal activity generates the BOLD signals. These models allow us to non-invasively infer the degree of local synaptic and spiking activity in the healthy human brain. In addition, we use Bayesian model comparison to decide between neurovascular coupling mechanisms. We show that the BOLD signal is dependent upon both the synaptic and spiking activity but that the relative contributions of these two inputs are dependent upon the underlying neuronal firing rate. When the underlying neuronal firing is low then the BOLD response is best explained by synaptic activity. However, when the neuronal firing rate is high then both synaptic and spiking activity are required to explain the BOLD signal.

Functional magnetic resonance imaging (fMRI), with blood oxygenation level-dependent (BOLD) contrast, is a widely used technique for studying the human brain. However, the relationship between neuronal activity and blood flow, the basis of fMRI, is still under much debate. A growing body of evidence from animal studies suggests that fMRI signals are more closely coupled to synaptic input activity than to the spiking output of a neuronal population. However, data from neurosurgical patients does not seem to support this view and this hypothesis hasn't yet been tested in the healthy human brain. Here we design a powerful and efficient modelling framework that can be used to non-invasively compare different biologically plausible hypotheses of neurovascular coupling. We use this framework to explore the contribution of these two aspects of neuronal activity (synaptic and spiking) to the generation of hemodynamic signals in human visual cortex, with Electroencephalographic (EEG)-fMRI data. Our results provide preliminary evidence that depending on the frequency of the visual stimulus and underlying firing rate, fMRI relates closer to synaptic activity (low-frequencies) or to both synaptic and spiking activities (high-frequencies).

Frequency-specific electrocorticographic correlates of working memory delay period fMRI activity

Electrocorticography (ECoG) and functional MRI (BOLD-fMRI) have been used previously to measure brain activity during working memory delay periods. These studies have separately reported oscillation changes in the theta (4-8 Hz) band and BOLD-fMRI increases during delay periods when information is maintained in memory. However, it is not known how intracranial cortical field potential (CFP) changes relate to BOLD-fMRI responses during delay periods. To answer this question, fMRI was obtained from six epilepsy patients during a visual working memory task. Then, following subdural macroelectrode implant, continuous ECoG was used to record CFPs during the same task. Time-frequency analyses showed delay period gamma band oscillation amplitude increases on electrodes located near fMRI activity, while in the theta band changes were higher for electrodes located away from fMRI activation. The amplitude of the ECoG gamma band response was significantly positively correlated with the fMRI response, while a negative correlation was found for the theta band. The findings are consistent with previous reports of local field potential (LFP) coupling in the gamma band with BOLD-fMRI responses during visual stimulation in monkeys, but are novel in that the relationship reported here persists after the disappearance of visual stimuli while information is being maintained in memory. We conclude that there is a relationship between BOLD-fMRI increases and human working memory delay period gamma oscillation increases and theta decreases. The spectral profile change provides a basis for comparison of working memory delay period BOLD-fMRI with field potential recordings in animals and other human intracranial EEG studies.