Neurophysiology of the BOLD fMRI signal in awake monkeys

BACKGROUND: Simultaneous intracortical recordings of neural activity and blood-oxygen-level-dependent (BOLD) functional magnetic resonance imaging (fMRI) in primary visual cortex of anesthetized monkeys demonstrated varying degrees of correlation between fMRI signals and the different types of neural activity, such as local field potentials (LFPs), multiple-unit activity (MUA), and single-unit activity (SUA). One important question raised by the aforementioned investigation is whether the reported correlations also apply to alert subjects. RESULTS: Monkeys were trained to perform a fixation task while stimuli within the receptive field of each recording site were used to elicit neural responses followed by a BOLD response. We show -- also in alert behaving monkeys -- that although both LFP and MUA make significant contributions to the BOLD response, LFPs are better and more reliable predictors of the BOLD signal. Moreover, when MUA responses adapt but LFP remains unaffected, the BOLD signal remains unaltered. CONCLUSIONS: The persistent coupling of the BOLD signal to the field potential when LFP and MUA have different time evolutions suggests that BOLD is primarily determined by the local processing of inputs in a given cortical area. In the alert animal the largest portion of the BOLD signal's variance is explained by an LFP range (20-60 Hz) that is most likely related to neuromodulation. Finally, the similarity of the results in alert and anesthetized subjects indicates that at least in V1 anesthesia is not a confounding factor. This enables the comparison of human fMRI results with a plethora of electrophysiological results obtained in alert or anesthetized animals.     

Relations between BOLD fMRI-Derived Resting Brain Activity and Cerebral Blood Flow

Consistent resting brain activity patterns have been repeatedly demonstrated using measures derived from resting BOLD fMRI data. While those metrics are presumed to reflect underlying spontaneous brain activity (SBA), it is challenging to prove that association because resting BOLD fMRI metrics are purely model-free and scale-free variables. Cerebral blood flow (CBF) is typically closely coupled to brain metabolism and is used as a surrogate marker for quantifying regional brain function, including resting function. Assessing the correlations between resting BOLD fMRI measures and CBF correlation should provide a means of linking of those measures to the underlying SBA, and a means to quantify those scale-free measures. The purpose of this paper was to examine the CBF correlations of 3 widely used neuroimaging-based SBA measures, including seed-region based functional connectivity (FC), regional homogeneity (ReHo), and amplitude of low frequency fluctuation (ALFF). Test-retest data were acquired to check the stability of potential correlations across time. Reproducible posterior cingulate cortex (PCC) FC vs regional CBF correlations were found in much of the default mode network and visual cortex. Dorsal anterior cingulate cortex (ACC) FC vs CBF correlations were consistently found in bilateral prefrontal cortex. Both ReHo and ALFF were found to be reliably correlated with CBF in most of brain cortex. None of the assessed SBA measures was correlated with whole brain mean CBF. These findings suggest that resting BOLD fMRI-derived measures are coupled with regional CBF and are therefore linked to regional SBA.     

Optogenetic Activation of CA1 Pyramidal Neurons at the Dorsal and Ventral Hippocampus Evokes Distinct Brain-Wide Responses Revealed by Mouse fMRI

The dorsal and ventral hippocampal regions (dHP and vHP) are proposed to have distinct functions. Electrophysiological studies have revealed intra-hippocampal variances along the dorsoventral axis. Nevertheless, the extra-hippocampal influences of dHP and vHP activities remain unclear. In this study, we compared the spatial distribution of brain-wide responses upon dHP or vHP activation and further estimate connection strengths between the dHP and the vHP with corresponding extra-hippocampal areas. To achieve this, we first investigated responses of local field potential (LFP) and multi unit activities (MUA) upon light stimulation in the hippocampus of an anesthetized transgenic mouse, whose CA1 pyramidal neurons expressed a step-function opsin variant of channelrhodopsin-2 (ChR2). Optogenetic stimulation increased hippocampal LFP power at theta, gamma, and ultra-fast frequency bands, and augmented MUA, indicating light-induced activation of CA1 pyramidal neurons. Brain-wide responses examined using fMRI revealed that optogenetic activation at the dHP or vHP caused blood oxygenation level-dependent (BOLD) fMRI signals in situ. Although activation at the dHP induced BOLD responses at the vHP, the opposite was not observed. Outside the hippocampal formation, activation at the dHP, but not the vHP, evoked BOLD responses at the retrosplenial cortex (RSP), which is in line with anatomical evidence. In contrast, BOLD responses at the lateral septum (LS) were induced only upon vHP activation, even though both dHP and vHP send axonal fibers to the LS. Our findings suggest that the primary targets of dHP and vHP activation are distinct, which concurs with attributed functions of the dHP and RSP in spatial memory, as well as of the vHP and LS in emotional responses.     

The local and non-local components of the local field potential in awake primate visual cortex

The Local Field Potential (LFP) is the analog signal recorded from a microelectrode inserted into cortex, typically in the frequency band of approximately 1 to 200 Hz. Here visual stimuli were flashed on in the receptive fields of primary visual cortical neurons in awake behaving macaques, and both isolated single units (neurons) and the LFP signal were recorded from the same unipolar microelectrode. The fall-off of single unit activity as a visual stimulus was moved from near the center to near the edge of the receptive field paralleled the fall-off of the stimulus-locked (evoked) LFP response. This suggests that the evoked LFP strongly reflects local neuronal activity. However, the evoked LFP could be significant even when the visual stimulus was completely outside the receptive field and the single unit response had fallen to zero, although this phenomenon was variable. Some of the non-local components of the LFP may be related to the slow distributed, or non-retinotopic, LFP signal previously observed in anesthetized animals. The induced (not time-locked to stimulus onset) component of the LFP showed significant increases only for stimuli within the receptive field of the single units. While the LFP primarily reflects local neuronal activity, it can also reflect neuronal activity at more distant sites, although these non-local components are typically more variable, slower, and weaker than the local components.     

Detecting Static and Dynamic Differences between Eyes-Closed and Eyes-Open Resting States Using ASL and BOLD fMRI

Resting-state fMRI studies have increasingly focused on multi-contrast techniques, such as BOLD and ASL imaging. However, these techniques may reveal different aspects of brain activity (e.g., static vs. dynamic), and little is known about the similarity or disparity of these techniques in detecting resting-state brain activity. It is therefore important to assess the static and dynamic characteristics of these fMRI techniques to guide future applications. Here we acquired fMRI data while subjects were in eyes-closed (EC) and eyes-open (EO) states, using both ASL and BOLD techniques, at two research centers (NIDA and HNU). Static brain activity was calculated as voxel-wise mean cerebral blood flow (CBF) using ASL, i.e., CBF-mean, while dynamic activity was measured by the amplitude of low frequency fluctuations (ALFF) of BOLD, i.e., BOLD-ALFF, at both NIDA and HNU, and CBF, i.e., CBF-ALFF, at NIDA. We showed that mean CBF was lower under EC than EO in the primary visual cortex, while BOLD-ALFF was higher under EC in the primary somatosensory cortices extending to the primary auditory cortices and lower in the lateral occipital area. Interestingly, mean CBF and BOLD-ALFF results overlapped at the visual cortex to a very small degree. Importantly, these findings were largely replicated by the HNU dataset. State differences found by CBF-ALFF were located in the primary auditory cortices, which were generally a subset of BOLD-ALFF and showed no spatial overlap with CBF-mean. In conclusion, static brain activity measured by mean CBF and dynamic brain activity measured by BOLD- and CBF-ALFF may reflect different aspects of resting-state brain activity and a combination of ASL and BOLD may provide complementary information on the biophysical and physiological processes of the brain.     

The neural basis of the blood-oxygen-level-dependent functional magnetic resonance imaging signal

Magnetic resonance imaging (MRI) has rapidly become an important tool in clinical medicine and biological research. Its functional variant (functional magnetic resonance imaging; fMRI) is currently the most widely used method for brain mapping and studying the neural basis of human cognition. While the method is widespread, there is insufficient knowledge of the physiological basis of the fMRI signal to interpret the data confidently with respect to neural activity. This paper reviews the basic principles of MRI and fMRI, and subsequently discusses in some detail the relationship between the blood-oxygen-level-dependent (BOLD) fMRI signal and the neural activity elicited during sensory stimulation. To examine this relationship, we conducted the first simultaneous intracortical recordings of neural signals and BOLD responses. Depending on the temporal characteristics of the stimulus, a moderate to strong correlation was found between the neural activity measured with microelectrodes and the BOLD signal averaged over a small area around the microelectrode tips. However, the BOLD signal had significantly higher variability than the neural activity, indicating that human fMRI combined with traditional statistical methods underestimates the reliability of the neuronal activity. To understand the relative contribution of several types of neuronal signals to the haemodynamic response, we compared local field potentials (LFPs), single- and multi-unit activity (MUA) with high spatio-temporal fMRI responses recorded simultaneously in monkey visual cortex. At recording sites characterized by transient responses, only the LFP signal was significantly correlated with the haemodynamic response. Furthermore, the LFPs had the largest magnitude signal and linear systems analysis showed that the LFPs were better than the MUAs at predicting the fMRI responses. These findings, together with an analysis of the neural signals, indicate that the BOLD signal primarily measures the input and processing of neuronal information within a region and not the output signal transmitted to other brain regions.     

Simultaneous EEG-fMRI during a Working Memory Task: Modulations in Low and High Frequency Bands

Background

EEG studies of working memory (WM) have demonstrated load dependent frequency band modulations. FMRI studies have localized load modulated activity to the dorsolateral prefrontal cortex (DLPFC), medial prefrontal cortex (MPFC), and posterior parietal cortex (PPC). Recently, an EEG-fMRI study found that low frequency band (theta and alpha) activity negatively correlated with the BOLD signal during the retention phase of a WM task. However, the coupling of higher (beta and gamma) frequencies with the BOLD signal during WM is unknown.

Methodology

In 16 healthy adult subjects, we first investigated EEG-BOLD signal correlations for theta (5–7 Hz), alpha1 (8–10), alpha2 (10–12 Hz), beta1 (13–20), beta2 (20–30 Hz), and gamma (30–40 Hz) during the retention period of a WM task with set size 2 and 5. Secondly, we investigated whether load sensitive brain regions are characterised by effects that relate frequency bands to BOLD signals effects.

Principal Findings

We found negative theta-BOLD signal correlations in the MPFC, PPC, and cingulate cortex (ACC and PCC). For alpha1 positive correlations with the BOLD signal were found in ACC, MPFC, and PCC; negative correlations were observed in DLPFC, PPC, and inferior frontal gyrus (IFG). Negative alpha2-BOLD signal correlations were observed in parieto-occipital regions. Beta1-BOLD signal correlations were positive in ACC and negative in precentral and superior temporal gyrus. Beta2 and gamma showed only positive correlations with BOLD, e.g., in DLPFC, MPFC (gamma) and IFG (beta2/gamma). The load analysis revealed that theta and—with one exception—beta and gamma demonstrated exclusively positive load effects, while alpha1 showed only negative effects.

Conclusions

We conclude that the directions of EEG-BOLD signal correlations vary across brain regions and EEG frequency bands. In addition, some brain regions show both load sensitive BOLD and frequency band effects. Our data indicate that lower as well as higher frequency brain oscillations are linked to neurovascular processes during WM.     

Multi-electrode Array Recordings of Neuronal Avalanches in Organotypic Cultures

The cortex is spontaneously active, even in the absence of any particular input or motor output. During development, this activity is important for the migration and differentiation of cortex cell types and the formation of neuronal connections¹. In the mature animal, ongoing activity reflects the past and the present state of an animal into which sensory stimuli are seamlessly integrated to compute future actions. Thus, a clear understanding of the organization of ongoing i.e. spontaneous activity is a prerequisite to understand cortex function. Numerous recording techniques revealed that ongoing activity in cortex is comprised of many neurons whose individual activities transiently sum to larger events that can be detected in the local field potential (LFP) with extracellular microelectrodes, or in the electroencephalogram (EEG), the magnetoencephalogram (MEG), and the BOLD signal from functional magnetic resonance imaging (fMRI). The LFP is currently the method of choice when studying neuronal population activity with high temporal and spatial resolution at the mesoscopic scale (several thousands of neurons). At the extracellular microelectrode, locally synchronized activities of spatially neighbored neurons result in rapid deflections in the LFP up to several hundreds of microvolts. When using an array of microelectrodes, the organizations of such deflections can be conveniently monitored in space and time. Neuronal avalanches describe the scale-invariant spatiotemporal organization of ongoing neuronal activity in the brain^2,3. They are specific to the superficial layers of cortex as established in vitro^4,5, in vivo in the anesthetized rat ⁶, and in the awake monkey⁷. Importantly, both theoretical and empirical studies^2,8-10 suggest that neuronal avalanches indicate an exquisitely balanced critical state dynamics of cortex that optimizes information transfer and information processing.In order to study the mechanisms of neuronal avalanche development, maintenance, and regulation, in vitro preparations are highly beneficial, as they allow for stable recordings of avalanche activity under precisely controlled conditions. The current protocol describes how to study neuronal avalanches in vitro by taking advantage of superficial layer development in organotypic cortex cultures, i.e. slice cultures, grown on planar, integrated microelectrode arrays (MEA; see also ^11-14).     

How local is the local field potential?

Local field potentials (LFPs) are of growing importance in neurophysiological investigations. LFPs supplement action potential recordings by indexing activity relevant to EEG, magnetoencephalographic, and hemodynamic (fMRI) signals. Recent reports suggest that LFPs reflect activity within very small domains of several hundred micrometers. We examined this conclusion by comparing LFP, current source density (CSD), and multiunit activity (MUA) signals in macaque auditory cortex. Estimated by frequency tuning bandwidths, these signals' "listening areas" differ systematically with an order of MUA?< CSD?< LFP. Computational analyses confirm that observed LFPs receive local contributions. Direct measurements indicate passive spread of LFPs to sites more than a centimeter from their origins. These findings appear to be independent of the frequency content of the LFP. Our results challenge the idea that LFP recordings typically integrate over extremely circumscribed local domains. Rather, LFPs appear as a mixture of local potentials with "volume conducted" potentials from distant sites.     

Mechanistic Mathematical Modeling Tests Hypotheses of the Neurovascular Coupling in fMRI

Functional magnetic resonance imaging (fMRI) measures brain activity by detecting the blood-oxygen-level dependent (BOLD) response to neural activity. The BOLD response depends on the neurovascular coupling, which connects cerebral blood flow, cerebral blood volume, and deoxyhemoglobin level to neuronal activity. The exact mechanisms behind this neurovascular coupling are not yet fully investigated. There are at least three different ways in which these mechanisms are being discussed. Firstly, mathematical models involving the so-called Balloon model describes the relation between oxygen metabolism, cerebral blood volume, and cerebral blood flow. However, the Balloon model does not describe cellular and biochemical mechanisms. Secondly, the metabolic feedback hypothesis, which is based on experimental findings on metabolism associated with brain activation, and thirdly, the neurotransmitter feed-forward hypothesis which describes intracellular pathways leading to vasoactive substance release. Both the metabolic feedback and the neurotransmitter feed-forward hypotheses have been extensively studied, but only experimentally. These two hypotheses have never been implemented as mathematical models. Here we investigate these two hypotheses by mechanistic mathematical modeling using a systems biology approach; these methods have been used in biological research for many years but never been applied to the BOLD response in fMRI. In the current work, model structures describing the metabolic feedback and the neurotransmitter feed-forward hypotheses were applied to measured BOLD responses in the visual cortex of 12 healthy volunteers. Evaluating each hypothesis separately shows that neither hypothesis alone can describe the data in a biologically plausible way. However, by adding metabolism to the neurotransmitter feed-forward model structure, we obtained a new model structure which is able to fit the estimation data and successfully predict new, independent validation data. These results open the door to a new type of fMRI analysis that more accurately reflects the true neuronal activity.     

Comparison of LFP-based and spike-based spectro-temporal receptive fields and cross-correlation in cat primary auditory cortex

Multi-electrode array recordings of spike and local field potential (LFP) activity were made from primary auditory cortex of 12 normal hearing, ketamine-anesthetized cats. We evaluated 259 spectro-temporal receptive fields (STRFs) and 492 frequency-tuning curves (FTCs) based on LFPs and spikes simultaneously recorded on the same electrode. We compared their characteristic frequency (CF) gradients and their cross-correlation distances. The CF gradient for spike-based FTCs was about twice that for 2-40 Hz-filtered LFP-based FTCs, indicating greatly reduced frequency selectivity for LFPs. We also present comparisons for LFPs band-pass filtered between 4-8 Hz, 8-16 Hz and 16-40 Hz, with spike-based STRFs, on the basis of their marginal frequency distributions. We find on average a significantly larger correlation between the spike based marginal frequency distributions and those based on the 16-40 Hz filtered LFP, compared to those based on the 4-8 Hz, 8-16 Hz and 2-40 Hz filtered LFP. This suggests greater frequency specificity for the 16-40 Hz LFPs compared to those of lower frequency content. For spontaneous LFP and spike activity we evaluated 1373 pair correlations for pairs with >200 spikes in 900 s per electrode. Peak correlation-coefficient space constants were similar for the 2-40 Hz filtered LFP (5.5 mm) and the 16-40 Hz LFP (7.4 mm), whereas for spike-pair correlations it was about half that, at 3.2 mm. Comparing spike-pairs with 2-40 Hz (and 16-40 Hz) LFP-pair correlations showed that about 16% (9%) of the variance in the spike-pair correlations could be explained from LFP-pair correlations recorded on the same electrodes within the same electrode array. This larger correlation distance combined with the reduced CF gradient and much broader frequency selectivity suggests that LFPs are not a substitute for spike activity in primary auditory cortex.     

Effects of L-dopa during Auditory Instrumental Learning in Humans

The dopaminergic neurotransmitter system is critically involved in promoting plasticity in auditory cortex. We combined functional magnetic resonance imaging (fMRI) and a pharmacological manipulation to investigate dopaminergic modulation of neural activity in auditory cortex during instrumental learning. Volunteers either received 100 mg L-dopa (Madopar) or placebo in an appetitive, differential instrumental conditioning paradigm, which involved learning that a specific category of frequency modulated tones predicts a monetary reward when fast responses were made in a subsequent reaction time task. The other category of frequency modulated tones was not related to a reward. Our behavioral data provides evidence that dopaminergic stimulation differentially impacts on the speed of instrumental responding in rewarded and unrewarded trials. L-dopa increased neural BOLD activity in left auditory cortex to tones in rewarded and unrewarded trials. This increase was related to plasma L-dopa levels and learning rate. Our data thus provides evidence for dopaminergic modulation of neural activity in auditory cortex, which occurs for both auditory stimuli related to a later reward and those not related to a reward.     

Modulation of activity in human visual area V1 during memory masking

Neurons in the primary visual cortex, V1, are specialized for the processing of elemental features of the visual stimulus, such as orientation and spatial frequency. Recent fMRI evidence suggest that V1 neurons are also recruited in visual perceptual memory; a number of studies using multi-voxel pattern analysis have successfully decoded stimulus-specific information from V1 activity patterns during the delay phase in memory tasks. However, consistent fMRI signal modulations reflecting the memory process have not yet been demonstrated. Here, we report evidence, from three subjects, that the low V1 BOLD activity during retention of low-level visual features is caused by competing interactions between neural populations coding for different values along the spectrum of the dimension remembered. We applied a memory masking paradigm in which the memory representation of a masker stimulus interferes with a delayed spatial frequency discrimination task when its frequency differs from the discriminanda with ±1 octave and found that impaired behavioral performance due to masking is reflected in weaker V1 BOLD signals. This cross-channel inhibition in V1 only occurs with retinotopic overlap between the masker and the sample stimulus of the discrimination task. The results suggest that memory for spatial frequency is a local process in the retinotopically organized visual cortex.     

Conversion of Phase Information into a Spike-Count Code by Bursting Neurons

Single neurons in the cerebral cortex are immersed in a fluctuating electric field, the local field potential (LFP), which mainly originates from synchronous synaptic input into the local neural neighborhood. As shown by recent studies in visual and auditory cortices, the angular phase of the LFP at the time of spike generation adds significant extra information about the external world, beyond the one contained in the firing rate alone. However, no biologically plausible mechanism has yet been suggested that allows downstream neurons to infer the phase of the LFP at the soma of their pre-synaptic afferents. Therefore, so far there is no evidence that the nervous system can process phase information. Here we study a model of a bursting pyramidal neuron, driven by a time-dependent stimulus. We show that the number of spikes per burst varies systematically with the phase of the fluctuating input at the time of burst onset. The mapping between input phase and number of spikes per burst is a robust response feature for a broad range of stimulus statistics. Our results suggest that cortical bursting neurons could play a crucial role in translating LFP phase information into an easily decodable spike count code.     

Functional imaging reveals numerous fields in the monkey auditory cortex

Anatomical studies propose that the primate auditory cortex contains more fields than have actually been functionally confirmed or described. Spatially resolved functional magnetic resonance imaging (fMRI) with carefully designed acoustical stimulation could be ideally suited to extend our understanding of the processing within these fields. However, after numerous experiments in humans, many auditory fields remain poorly characterized. Imaging the macaque monkey is of particular interest as these species have a richer set of anatomical and neurophysiological data to clarify the source of the imaged activity. We functionally mapped the auditory cortex of behaving and of anesthetized macaque monkeys with high resolution fMRI. By optimizing our imaging and stimulation procedures, we obtained robust activity throughout auditory cortex using tonal and band-passed noise sounds. Then, by varying the frequency content of the sounds, spatially specific activity patterns were observed over this region. As a result, the activity patterns could be assigned to many auditory cortical fields, including those whose functional properties were previously undescribed. The results provide an extensive functional tessellation of the macaque auditory cortex and suggest that 11 fields contain neurons tuned for the frequency of sounds. This study provides functional support for a model where three fields in primary auditory cortex are surrounded by eight neighboring “belt” fields in non-primary auditory cortex. The findings can now guide neurophysiological recordings in the monkey to expand our understanding of the processing within these fields. Additionally, this work will improve fMRI investigations of the human auditory cortex.     

Conserved fMRI and LFP signals during new associative learning in the human and macaque monkey medial temporal lobe

We measured local field potential (LFP) and blood-oxygen-level-dependent (BOLD) functional magnetic resonance imaging (fMRI) in the medial temporal lobes of monkeys and humans, respectively, as they performed the same conditional motor associative learning task. Parallel analyses were used to examine both data sets. Despite significantly faster learning in humans relative to monkeys, we found equivalent neural signals differentiating new versus highly familiar stimuli, first stimulus presentation, trial outcome, and learning strength in the entorhinal cortex and hippocampus of both species. Thus, the use of parallel behavioral tasks and analyses in monkeys and humans revealed conserved patterns of neural activity across the medial temporal lobe during an associative learning task.     

Covariation between Spike and LFP Modulations Revealed with Focal and Asynchronous Stimulation of Receptive Field Surround in Monkey Primary Visual Cortex

A focal visual stimulus outside the classical receptive field (RF) of a V1 neuron does not evoke a spike response by itself, and yet evokes robust changes in the local field potential (LFP). This subthreshold LFP provides a unique opportunity to investigate how changes induced by surround stimulation leads to modulation of spike activity. In the current study, two identical Gabor stimuli were sequentially presented with a variable stimulus onset asynchrony (SOA) ranging from 0 to 100 ms: the first (S1) outside the RF and the second (S2) over the RF of primary visual cortex neurons, while trained monkeys performed a fixation task. This focal and asynchronous stimulation of the RF surround enabled us to analyze the modulation of S2-evoked spike activity and covariation between spike and LFP modulation across SOA. In this condition, the modulation of S2-evoked spike response was dominantly facilitative and was correlated with the change in LFP amplitude, which was pronounced for the cells recorded in the upper cortical layers. The time course of covariation between the SOA-dependent spike modulation and LFP amplitude suggested that the subthreshold LFP evoked by the S1 can predict the magnitude of upcoming spike modulation.     

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.     

Cortical Neurovascular Coupling Driven by Stimulation of Channelrhodopsin-2

While functional imaging is widely used in studies of the brain, how well the hemodynamic signal represents the underlying neural activity is still unclear. And there is a debate on whether hemodynamic signal is more tightly related to synaptic activity or action potentials. This study intends to address these questions by examining neurovascular coupling driven by pyramidal cells in the motor cortex of rats. Pyramidal cells in the motor cortex of rats were selectively transduced with the light sensitive cation channel channelrhodopsin-2 (ChR2). Electrophysiological recordings and optical intrinsic signal imaging were performed simultaneously and synchronously to capture the neural activity and hemodynamics induced by optical stimulation of ChR2-expressing pyramidal cells. Our results indicate that both synaptic activity (local field potential, LFP) and action potentials (multi-unit activity, MUA) are tightly related to hemodynamic signals. While LFPs in γ band are better in predicting hemodynamic signals elicited by short stimuli, MUA has better predictions to hemodynamic signals elicited by long stimuli. Our results also indicate that strong nonlinearity exists in neurovascular coupling.     

Methods for predicting cortical UP and DOWN states from the phase of deep layer local field potentials

During anesthesia, slow-wave sleep and quiet wakefulness, neuronal membrane potentials collectively switch between de- and hyperpolarized levels, the cortical UP and DOWN states. Previous studies have shown that these cortical UP/DOWN states affect the excitability of individual neurons in response to sensory stimuli, indicating that a significant amount of the trial-to-trial variability in neuronal responses can be attributed to ongoing fluctuations in network activity. However, as intracellular recordings are frequently not available, it is important to be able to estimate their occurrence purely from extracellular data. Here, we combine in vivo whole cell recordings from single neurons with multi-site extracellular microelectrode recordings, to quantify the performance of various approaches to predicting UP/DOWN states from the deep-layer local field potential (LFP). We find that UP/DOWN states in deep cortical layers of rat primary auditory cortex (A1) are predictable from the phase of LFP at low frequencies (< 4 Hz), and that the likelihood of a given state varies sinusoidally with the phase of LFP at these frequencies. We introduce a novel method of detecting cortical state by combining information concerning the phase of the LFP and ongoing multi-unit activity.