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 BOLD fMRI and fiber-optic calcium recording in rat neocortex

Functional magnetic resonance imaging (fMRI) based on blood oxygen level-dependent (BOLD) contrast is widely used for probing brain activity, but its relationship to underlying neural activity remains elusive. Here, we combined fMRI with fiber-optic recordings of fluorescent calcium indicator signals to investigate this relationship in rat somatosensory cortex. Electrical forepaw stimulation (1-10 Hz) evoked fast calcium signals of neuronal origin that showed frequency-dependent adaptation. Additionally, slower calcium signals occurred in astrocyte networks, as verified by astrocyte-specific staining and two-photon microscopy. Without apparent glia activation, we could predict BOLD responses well from simultaneously recorded fiber-optic signals, assuming an impulse response function and taking into account neuronal adaptation. In cases with glia activation, we uncovered additional prolonged BOLD signal components. Our findings highlight the complexity of fMRI BOLD signals, involving both neuronal and glial activity. Combined fMRI and fiber-optic recordings should help to clarify cellular mechanisms underlying BOLD signals.     

Repeated BOLD-fMRI Imaging of Deep Brain Stimulation Responses in Rats

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

Optogenetic Functional MRI

The investigation of the functional connectivity of precise neural circuits across the entire intact brain can be achieved through optogenetic functional magnetic resonance imaging (ofMRI), which is a novel technique that combines the relatively high spatial resolution of high-field fMRI with the precision of optogenetic stimulation. Fiber optics that enable delivery of specific wavelengths of light deep into the brain in vivo are implanted into regions of interest in order to specifically stimulate targeted cell types that have been genetically induced to express light-sensitive trans-membrane conductance channels, called opsins. fMRI is used to provide a non-invasive method of determining the brain''s global dynamic response to optogenetic stimulation of specific neural circuits through measurement of the blood-oxygen-level-dependent (BOLD) signal, which provides an indirect measurement of neuronal activity. This protocol describes the construction of fiber optic implants, the implantation surgeries, the imaging with photostimulation and the data analysis required to successfully perform ofMRI. In summary, the precise stimulation and whole-brain monitoring ability of ofMRI are crucial factors in making ofMRI a powerful tool for the study of the connectomics of the brain in both healthy and diseased states.     

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.     

BOLD signals correlate with ensemble unit activities in rat's somatosensory cortex

Evoked neural activity (ensemble single-unit activity and evoked field potential) and functional magnetic resonance imaging (fMRI) changes of the primary somatosensory cortex in response to electrical stimulation of the hind paw were studied in rats under anesthesia. The effects of stimulation frequency (ranging from 0.3 to 10 Hz) and types of anesthetics (alpha-chloralose and sodium pentobarbital) on blood oxygen level dependent (BOLD) activation and neural activation were compared. Both ensemble single-unit activity and BOLD signal changes achieved maximal activation at 3 Hz of stimulation and responses were significantly stronger under alpha-chloralose anesthesia. The maximal activation of the integral evoked potential (sigmaEP), in contrast, was the highest at 10 Hz; and the values were similar for alpha-chloralose and pentobarbital. These analyses revealed that fMRI image changes were better correlated with ensemble single-unit activity than with sigmaEP during somatosensory stimulations.     

Positive Allosteric Modulator of GABA Lowers BOLD Responses in the Cingulate Cortex

Knowledge about the neural underpinnings of the negative blood oxygen level dependent (BOLD) responses in functional magnetic resonance imaging (fMRI) is still limited. We hypothesized that pharmacological GABAergic modulation attenuates BOLD responses, and that blood concentrations of a positive allosteric modulator of GABA correlate inversely with BOLD responses in the cingulate cortex. We investigated whether or not pure task-related negative BOLD responses were co-localized with pharmacologically modulated BOLD responses. Twenty healthy adults received either 5 mg diazepam or placebo in a double blind, randomized design. During fMRI the subjects performed a working memory task. Results showed that BOLD responses in the cingulate cortex were inversely correlated with diazepam blood concentrations; that is, the higher the blood diazepam concentration, the lower the BOLD response. This inverse correlation was most pronounced in the pregenual anterior cingulate cortex and the anterior mid-cingulate cortex. For subjects with diazepam plasma concentration > 0.1 mg/L we observed negative BOLD responses with respect to fixation baseline. There was minor overlap between cingulate regions with task-related negative BOLD responses and regions where the BOLD responses were inversely correlated with diazepam concentration. We interpret that the inverse correlation between the BOLD response and diazepam was caused by GABA-related neural inhibition. Thus, this study supports the hypothesis that GABA attenuates BOLD responses in fMRI. The minimal overlap between task-related negative BOLD responses and responses attenuated by diazepam suggests that these responses might be caused by different mechanisms.     

Distinct temporal patterns of electrical stimulation influence neural recruitment during PTZ infusion: an fMRI study

Our working hypothesis is that constant inter-pulse interval (IPI) electrical stimulation (ES) would resonate with endogenous epileptogenic reverberating circuits, favoring seizure, while random inter-interval ES protocol would promote desynchronization of such neural networks, interfering with the abnormal recruitment of neural structures. Male Wistar rats were stereotaxically implanted with a monopolar ES carbon-fiber electrode (minimizing fMRI artifact) in the amygdala. A 7T fMRI scanner was used to evaluate brain activity during ES, fixed four pulses per second ratio, using either a periodic IPI (ES-P) or random IPI (non-periodic ES-NP) stimulation paradigm. Appropriate imaging protocols were used to compare baseline BOLD (blood oxygen level dependent) MRI with scans during ES. A second series of experiments, both without stimuli and under the same ES paradigms, were evaluated during continuous infusion of pentylenetetrazole (PTZ, 4 mg/ml/min) through an i.v. catheter. Our results show that temporal lobe activation during ES-P or ES-NP did not present any statistical differences during ES. However, during PTZ infusion, PTZ-P facilitated recruitment of the temporal lobe ipsilateral to ES while PTZ-NP showed significantly less activation ipsilateral to ES and, in turn, less inter-hemispheric differences. Altogether, our results support the hypothesis of reverberating circuits being synchronized by ES-P and desynchronized by ES-NP. Time-coded low frequency stimulation may be an interesting alternative treatment for patients with refractory epilepsy.     

Functional magnetic resonance imaging: imaging techniques and contrast mechanisms.

Functional magnetic resonance imaging (fMRI) is a widely used technique for generating images or maps of human brain activity. The applications of the technique are widespread in cognitive neuroscience and it is hoped they will eventually extend into clinical practice. The activation signal measured with fMRI is predicated on indirectly measuring changes in the concentration of deoxyhaemoglobin which arise from an increase in blood oxygenation in the vicinity of neuronal firing. The exact mechanisms of this blood oxygenation level dependent (BOLD) contrast are highly complex. The signal measured is dependent on both the underlying physiological events and the imaging physics. BOLD contrast, although sensitive, is not a quantifiable measure of neuronal activity. A number of different imaging techniques and parameters can be used for fMRI, the choice of which depends on the particular requirements of each functional imaging experiment. The high-speed MRI technique, echo-planar imaging provides the basis for most fMRI experiments. The problems inherent to this method and the ways in which these may be overcome are particularly important in the move towards performing functional studies on higher field MRI systems. Future developments in techniques and hardware are also likely to enhance the measurement of brain activity using MRI.     

Functional magnetic resonanse imaging and neuroscience

This review describes functional magnetic resonance imaging (fMRI), the intravital noninvasive dynamic study of the active zones of cerebral structures at the time of their activity, fMR is based on the difference in magnetic properties of oxyhemoglobin--the oxygen carrier, and deoxyhemoglobin--a product produced in the areas of oxygen consumption, the brain parenchyma. This ratio is reflected as the physical phenomenon BOLD (blood oxygenation level dependent) that is a marker of neuronal activity. fMRI has good spatial resolution and the possibility of multiple repetitions of the research. This allows the 3D reconstruction of the sequence of formation and intracerebral "geometry" (stereometry) of newly formed neural ensembles (NE) and/or realization of potentially pre-existing NE. fMRI is an optimal tool for neuronal activity mapping, or more accurate, for functional state of the NE in the reconstruction of neural networks; it should be regarded as the technology for studying the brain of humans and animals, both in terms of natural life activity, and in pathological conditions.     

Biophysical models of fMRI responses

Functional magnetic resonance imaging (fMRI) is used to investigate where the neural implementation of specific cognitive processes occurs. The standard approach uses linear convolution models that relate experimentally designed inputs, through a haemodynamic response function, to observed blood oxygen level dependent (BOLD) signals. Such models are, however, blind to the causal mechanisms that underlie observed BOLD responses. Recent developments have focused on how BOLD responses are generated and include biophysical input-state-output models with neural and haemodynamic state equations and models of functional integration that explain local dynamics through interactions with remote areas. Forward models with parameters at the neural level, such as dynamic causal modelling, combine both approaches, modelling the whole causal chain from external stimuli, via induced neural dynamics, to observed BOLD responses.     

Negative BOLD fMRI response in the visual cortex carries precise stimulus-specific information

Sustained positive BOLD (blood oxygen level-dependent) activity is employed extensively in functional magnetic resonance imaging (fMRI) studies as evidence for task or stimulus-specific neural responses. However, the presence of sustained negative BOLD activity (i.e., sustained responses that are lower than the fixation baseline) has remained more difficult to interpret. Some studies suggest that it results from local "blood stealing" wherein blood is diverted to neurally active regions without a concomitant change of neural activity in the negative BOLD regions. However, other evidence suggests that negative BOLD is a result of local neural suppression. In both cases, regions of negative BOLD response are usually interpreted as carrying relatively little, if any, stimulus-specific information (hence the predominant reliance on positive BOLD activity in fMRI). Here we show that the negative BOLD response resulting from visual stimulation can carry high information content that is stimulus-specific. Using a general linear model (GLM), we contrasted standard flickering stimuli to a fixation baseline and found regions of the visual cortex that displayed a sustained negative BOLD response, consistent with several previous studies. Within these negative BOLD regions, we compared patterns of fMRI activity generated by flickering Gabors that were systematically shifted in position. As the Gabors were shifted further from each other, the correlation in the spatial pattern of activity across a population of voxels (such as the population of V1 voxels that displayed a negative BOLD response) decreased significantly. Despite the fact that the BOLD signal was significantly negative (lower than fixation baseline), these regions were able to discriminate objects separated by less than 0.5 deg (at approximately 10 deg eccentricity). The results suggest that meaningful, stimulus-specific processing occurs even in regions that display a strong negative BOLD response.     

Functional Imaging of Auditory Cortex in Adult Cats using High-field fMRI

Current knowledge of sensory processing in the mammalian auditory system is mainly derived from electrophysiological studies in a variety of animal models, including monkeys, ferrets, bats, rodents, and cats. In order to draw suitable parallels between human and animal models of auditory function, it is important to establish a bridge between human functional imaging studies and animal electrophysiological studies. Functional magnetic resonance imaging (fMRI) is an established, minimally invasive method of measuring broad patterns of hemodynamic activity across different regions of the cerebral cortex. This technique is widely used to probe sensory function in the human brain, is a useful tool in linking studies of auditory processing in both humans and animals and has been successfully used to investigate auditory function in monkeys and rodents. The following protocol describes an experimental procedure for investigating auditory function in anesthetized adult cats by measuring stimulus-evoked hemodynamic changes in auditory cortex using fMRI. This method facilitates comparison of the hemodynamic responses across different models of auditory function thus leading to a better understanding of species-independent features of the mammalian auditory cortex.     

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.     

Coupling between neuronal firing rate, gamma LFP, and BOLD fMRI is related to interneuronal correlations

BACKGROUND: To what extent is activity of individual neurons coupled to the local field potential (LFP) and to blood-oxygenation-level dependent (BOLD) functional magnetic resonance imaging (fMRI)? This issue is of high significance for understanding brain function and for relating animal studies to fMRI, yet it is still under debate. RESULTS: Here we report data from simultaneous recordings of isolated unit activity and LFP by using multiple electrodes in the human auditory cortex. We found a wide range of coupling levels between the activity of individual neurons and gamma LFP. However, this large variability could be predominantly explained (r = 0.66) by the degree of firing-rate correlations between neighboring neurons. Importantly, this phenomenon occurred during both sensory stimulation and spontaneous activity. Concerning the coupling of neuronal activity to BOLD fMRI, we found that gamma LFP was well coupled to BOLD measured across different individuals (r = 0.62). By contrast, the coupling of single units to BOLD was highly variable and, again, tightly related to interneuronal-firing-rate correlations (r = 0.70). CONCLUSIONS: Our results offer a resolution to a central controversy regarding the coupling between neurons, LFP, and BOLD signals by demonstrating, for the first time, that the coupling of single units to the other measures is variable yet it is tightly related to the degree of interneuronal correlations in the human auditory cortex.     

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

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

Developmental Switch in Neurovascular Coupling in the Immature Rodent Barrel Cortex

Neurovascular coupling (NVC) in the adult central nervous system (CNS) is a mechanism that provides regions of the brain with more oxygen and glucose upon increased levels of neural activation. Hemodynamic changes that go along with neural activation evoke a blood oxygen level-dependent (BOLD) signal in functional magnetic resonance imaging (fMRI) that can be used to study brain activity non-invasively. A correct correlation of the BOLD signal to neural activity is pivotal to understand this signal in neuronal development, health and disease. However, the function of NVC during development is largely unknown. The rodent whisker-to-barrel cortex is an experimentally well established model to study neurovascular interdependences. Using extracellular multi-electrode recordings and laser-Doppler-flowmetry (LDF) we show in the murine barrel cortex of postnatal day 7 (P7) and P30 mice in vivo that NVC undergoes a physiological shift during the first month of life. In the mature CNS it is well accepted that cortical sensory processing results in a rise in regional cerebral blood flow (rCBF). We show in P7 animals that rCBF decreases during prolonged multi-whisker stimulation and goes along with multi unit activity (MUA) fatigue. In contrast at P30, MUA remains stable during repetitive stimulation and is associated with an increase in rCBF. Further we characterize in both age groups the responses in NVC to single sensory stimuli. We suggest that the observed shift in NVC is an important process in cortical development that may be of high relevance for the correct interpretation of brain activity e.g. in fMRI studies of the immature central nervous system (CNS).     

Self-regulation of local brain activity using real-time functional magnetic resonance imaging (fMRI)

Functional magnetic resonance imaging (fMRI) measures the blood oxygen level-dependent (BOLD) signal related to neuronal activity. So far, this technique has been limited by time-consuming data analysis impeding on-line analysis. In particular, no brain-computer interface (BCI) was available which provided on-line feedback to learn physiological self-regulation of the BOLD signal. Recently, studies have shown that fMRI feedback is feasible and facilitates voluntary control of brain activity. Here we review these studies to make the fMRI feedback methodology accessible to a broader scientific community such as researchers concerned with functional brain imaging and the neurobiology of learning. Methodological and conceptual limitations were substantially reduced by artefact control, sensitivity improvements, real-time algorithms, and adapted experimental designs. Physiological self-regulation of the local BOLD response is a new paradigm for cognitive neuroscience to study brain plasticity and the functional relevance of regulated brain areas by modification of behaviour. Voluntary control of abnormal activity in circumscribed brain areas may even be applied as psychophysiological treatment.     

Speech facilitation by left inferior frontal cortex stimulation

Electrophysiological studies in humans and animals suggest that noninvasive neurostimulation methods such as transcranial direct current stimulation (tDCS) can elicit long-lasting [1], polarity-dependent [2] changes in neocortical excitability. Application of tDCS can have significant and selective behavioral consequences that are associated with the cortical location of the stimulation electrodes and the task engaged during stimulation [3-8]. However, the mechanism by which tDCS affects human behavior is unclear. Recently, functional magnetic resonance imaging (fMRI) has been used to determine the spatial topography of tDCS effects [9-13], but no behavioral data were collected during stimulation. The present study is unique in this regard, in that both neural and behavioral responses were recorded using a novel combination of left frontal anodal tDCS during an overt picture-naming fMRI study. We found that tDCS had significant behavioral and regionally specific neural facilitation effects. Furthermore, faster naming responses correlated with decreased blood oxygen level-dependent (BOLD) signal in Broca's area. Our data support the importance of Broca's area within the normal naming network and as such indicate that Broca's area may be a suitable candidate site for tDCS in neurorehabilitation of anomic patients, whose brain damage spares this region.