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.     

Vasomotion and neurovascular coupling in the visual thalamus in vivo

Spontaneous contraction and relaxation of arteries (and in some instances venules) has been termed vasomotion and has been observed in an extensive variety of tissues and species. However, its functions and underlying mechanisms are still under discussion. We demonstrate that in vivo spectrophotometry, measured simultaneously with extracellular recordings at the same locations in the visual thalamus of the cat, reveals vasomotion, measured as an oscillation (0.14 hz) in the recorded oxyhemoglobin (OxyHb) signal, which appears spontaneously in the microcirculation and can last for periods of hours. During some non-oscillatory periods, maintained sensory stimulation evokes vasomotion lasting ~30s, resembling an adaptive vascular phenomenon. This oscillation in the oxyhaemoblobin signal is sensitive to pharmacological manipulation: it is inducible by chloralose anaesthesia and it can be temporarily blocked by systemic administration of adrenaline or acetylcholine (ACh). During these oscillatory periods, neurovascular coupling (i.e. the relationship between local neural activity and the rate of blood supply to that location) appears significantly altered. This raises important questions with regard to the interpretation of results from studies currently dependent upon a linear relationship between neural activity and blood flow, such as neuroimaging.     

An artery-specific fluorescent dye for studying neurovascular coupling

We demonstrate that Alexa Fluor 633 hydrazide (Alexa Fluor 633) selectively labels neocortical arteries and arterioles by binding to elastin fibers. We measured sensory stimulus-evoked arteriole dilation dynamics in mouse, rat and cat visual cortex using Alexa Fluor 633 together with neuronal activity using calcium indicators or blood flow using fluorescein dextran. Arteriole dilation decreased fluorescence recorded from immediately underlying neurons, representing a potential artifact during neuronal functional imaging experiments.     

Role of astrocytes in trimethyltin neurotoxicity

Although the neurotoxicity of trimethyltin (TMT) is well known, mechanisms are still not clear. Glia have been proposed to mediate the toxic action of TMT on nerve cells. Accordingly, the effects of TMT were tested in primary neuronal cultures from rat cerebellum and compared to effects in astrocytes and mixed cultures. Neuronal damage observed following TMT exposure was less in the presence of astrocytes and astrocytes alone were resistant to TMT. Thus, astrocytes have a protective effect against TMT-induced neurotoxicity. TMT caused an oxidative stress in granule cell cultures involving a variety of oxidative species (O2)*-, H2O2, NO), but astrocytes were less sensitive to TMT-induced oxidative species generation. Antioxidants, glutathione and 7-nitroindazole attenuated neuronal cell death induced by TMT. It appears that oxidative stress mediates a large part of the destructive action of TMT in neuronal cultures. The presence of astrocytes appears to modulate TMT-induced oxidative stress so that TMT causes only a small increase in lipid peroxidation in mouse brain after systemic administration. Thus, TMT induces a pronounced oxidative stress in cultured neurons, but when astrocytes are present, oxidative species play a lesser role in the neurotoxic action of TMT.     

Models of neurovascular coupling via potassium and EET signalling

Functional hyperemia is an important metabolic autoregulation mechanism by which increased neuronal activity is matched by a rapid and regional increase in blood supply. This mechanism is facilitated by a process known as "neurovascular coupling"--the orchestrated communication system involving neurons, astrocytes and arterioles. Important steps in this process are the production of EETs in the astrocyte and the release of potassium, via two potassium channels (BK and KIR), into the perivascular space. We provide a model which successfully accounts for several observations seen in experiment. The model is capable of simulating the approximate 15% arteriolar dilation caused by a 60-s neuronal activation (modelled as a release of potassium and glutamate into the synaptic cleft). This model also successfully emulates the paradoxical experimental finding that vasoconstriction follows vasodilation when the astrocytic calcium concentration (or perivascular potassium concentration) is increased further. We suggest that the interaction of the changing smooth muscle cell membrane potential and the changing potassium-dependent resting potential of the KIR channel are responsible for this effect. Finally, we demonstrate that a well-controlled mechanism of potassium buffering is potentially important for successful neurovascular coupling.     

Role of astrocytes in cerebrovascular regulation.

Astrocytes send processes to synapses and blood vessels, communicate with other astrocytes through gap junctions and by release of ATP, and thus are an integral component of the neurovascular unit. Electrical field stimulations in brain slices demonstrate an increase in intracellular calcium in astrocyte cell bodies transmitted to perivascular end-feet, followed by a decrease in vascular smooth muscle calcium oscillations and arteriolar dilation. The increase in astrocyte calcium after neuronal activation is mediated, in part, by activation of metabotropic glutamate receptors. Calcium signaling in vitro can also be influenced by adenosine acting on A2B receptors and by epoxyeicosatrienoic acids (EETs) shown to be synthesized in astrocytes. Prostaglandins, EETs, arachidonic acid, and potassium ions are candidate mediators of communication between astrocyte end-feet and vascular smooth muscle. In vivo evidence supports a role for cyclooxygenase-2 metabolites, EETs, adenosine, and neuronally derived nitric oxide in the coupling of increased blood flow to increased neuronal activity. Combined inhibition of the EETs, nitric oxide, and adenosine pathways indicates that signaling is not by parallel, independent pathways. Indirect pharmacological results are consistent with astrocytes acting as intermediaries in neurovascular signaling within the neurovascular unit. For specific stimuli, astrocytes are also capable of transmitting signals to pial arterioles on the brain surface for ensuring adequate inflow pressure to parenchymal feeding arterioles. Therefore, evidence from brain slices and indirect evidence in vivo with pharmacological approaches suggest that astrocytes play a pivotal role in regulating the fundamental physiological response coupling dynamic changes in cerebral blood flow to neuronal synaptic activity. Future work using in vivo imaging and genetic manipulation will be required to provide more direct evidence for a role of astrocytes in neurovascular coupling.     

The role of astrocytic calcium and TRPV4 channels in neurovascular coupling

Neuronal activity evokes a localised change in cerebral blood flow in a response known as neurovascular coupling (NVC). Although NVC has been widely studied the exact mechanisms that mediate this response remain unclear; in particular the role of astrocytic calcium is controversial. Mathematical modelling can be a useful tool for investigating the contribution of various signalling pathways towards NVC and for analysing the underlying cellular mechanisms. The lumped parameter model of a neurovascular unit with both potassium and nitric oxide (NO) signalling pathways and comprised of neurons, astrocytes, and vascular cells has been extended to include the glutamate induced astrocytic calcium pathway with epoxyeicosatrienoic acid (EET) signalling and the stretch dependent TRPV4 calcium channel on the astrocytic endfoot. Results show that the potassium pathway governs the fast onset of vasodilation while the NO pathway has a delayed response, maintaining dilation longer following neuronal stimulation. Increases in astrocytic calcium concentration via the calcium signalling pathway and/or TRPV4 channel to levels consistent with experimental data are insufficient for inducing either vasodilation or constriction, in contrast to a number of experimental results. It is shown that the astrocyte must depolarise in order to produce a significant potassium flux through the astrocytic BK channel. However astrocytic calcium is shown to strengthen potassium induced NVC by opening the BK channel further, consequently allowing more potassium into the perivascular space. The overall effect is vasodilation with a higher maximal vessel radius.     

A model of neurovascular coupling and the BOLD response PART II

A mathematical model is developed which describes a signalling mechanism of neurovascular coupling with a model of a pyramidal neuron and its corresponding fMRI BOLD response. In the first part of two papers (Part I) we described the integration of the neurovascular coupling unit extended to include a complex neuron model, which includes the important Na/K ATPase pump, with a model that provides a BOLD signal taking its input from the cerebral blood flow and the metabolic rate of oxygen consumption. We showed that this produced a viable signal in terms of initial dip, positive and negative BOLD signals. In this paper (PART II) our model predicts the variations of the BOLD response due to variations in neuronal activity and indicates that the BOLD signal could be used as an initial biomarker for neuronal dysfunction or variations in the perfusion of blood to the cerebral tissue. We have compared the simulated hypoxic BOLD response to experimental BOLD signals observed in the hippocampus during hypoxia showing good agreement. This approach of combined quantitative modelling of neurovascular coupling response and its BOLD response will enable more specific assessment of a brain region.     

NADPH-Diaphorase reactivity and neurovascular coupling in the basal forebrain and motor cortex

Using NADPH-diaphorase histochemistry, distribution of reactive neurons in the forebrain structures and motor cortex of rats was studied. Some reactive (NO-generating) neurons are associated with microvessels and, thus, can be involved in the regulation of regional blood flow. Neirofiziologiya/Neurophysiology, Vol. 39, Nos. 4/5, pp. 405–407, July–October, 2007.     

Antioxidant and bioenergetic coupling between neurons and astrocytes

Oxidative and nitrosative stress underlie the pathogenesis of a broad range of human diseases, in particular neurodegenerative disorders. Within the brain, neurons are the cells most vulnerable to excess reactive oxygen and nitrogen species; their survival relies on the antioxidant protection promoted by neighbouring astrocytes. However, neurons are also intrinsically equipped with a biochemical mechanism that links glucose metabolism to antioxidant defence. Neurons actively metabolize glucose through the pentose phosphate pathway, which maintains the antioxidant glutathione in its reduced state, hence exerting neuroprotection. This process is tightly controlled by a key glycolysis-promoting enzyme and is dependent on an appropriate supply of energy substrates from astrocytes. Thus brain bioenergetic and antioxidant defence is coupled between neurons and astrocytes. A better understanding of the regulation of this intercellular coupling should be important for identifying novel targets for future therapeutic interventions.     

Electrophysiological properties and gap junction coupling of striatal astrocytes

The striatum is the biggest nucleus of the basal ganglia and receives input from almost all cortical regions, substantia nigra and the thalamus. Striatal neuronal circuitry is well characterized, but less is known about glial physiology. To this end, we evaluated astrocyte electrophysiological properties using whole-cell patch-clamp recording in dorsal striatal brain slices from P15 to P21 rat. The majority of cells (95%) were passive astrocytes that do not express any detectable voltage-gated channels. Passive astrocytes were subcategorized into three groups based on time-dependent current properties. The observed proportion of the different astrocyte subtypes did not change within the age range evaluated here, but was modulated during reduction of specific conductances and gap junction coupling. Striatal astrocytes were extensively interconnected and closure of gap junctions with octanol (1 mM), carbenoxolone (100 μM) or increased intracellular calcium (2 mM), significantly altered intrinsic properties. When simultaneously blocking potassium channels and gap junction coupling almost no passive conductance was detected, implying that the major currents in striatal astrocytes derive from potassium and gap junction conductance. Uncoupling of the syncytium reduced currents activated in response to a hyperpolarizing pulse, suggesting that changes in gap junction coupling alters astrocyte electrophysiological responses. Our findings indicate that the prevalent gap junction coupling is vital for astrocyte function in the striatum, and that whole-cell recordings will be distorted by currents activated in neighboring cells.     

Somatostatin alters beta-adrenergic receptor-effector coupling in cultured rat astrocytes.

The neuropeptide somatostatin potentiates beta-adrenergic receptor-mediated cAMP formation in astrocytes derived from neonatal rat cortex but does not affect cAMP levels by itself. beta-Adrenergic receptors in these cells can be specifically labeled with the high affinity antagonist [125I] cyanopindolol ([125I]CYP). In addition, astrocytes display both high and low affinity binding sites for the agonist isoproterenol, which are thought to represent receptors which are coupled or uncoupled, respectively, to the guanine nucleotide regulatory protein. We find that somatostatin does not modify beta-receptor density, nor receptor affinity for either the antagonist ([125I]CYP) or for the agonist isoproterenol. In the presence of the guanine nucleotide analogue, Gpp(NH)p, only low affinity (uncoupled) displacement of [125I]CYP binding by isoproterenol is observed. However, somatostatin (1 microM), when added to the cells together with Gpp(NH)p, prevents the nucleotide-induced loss of the high affinity (coupled) component of agonist displacement. This result suggests that somatostatin increases noradrenaline-induced cAMP production by enhancing coupling between the beta-receptor and the stimulatory guanine nucleotide regulatory protein.     

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

Time-dependent spatial specificity of high-resolution fMRI: insights into mesoscopic neurovascular coupling

High-resolution functional magnetic resonance imaging (fMRI) is becoming increasingly popular because of the growing availability of ultra-high magnetic fields which are capable of improving sensitivity and spatial resolution. However, it is debatable whether increased spatial resolutions for haemodynamic-based techniques, like fMRI, can accurately detect the true location of neuronal activity. We have addressed this issue in functional columns and layers of animals with haemoglobin-based optical imaging and different fMRI contrasts, such as blood oxygenation level-dependent, cerebral blood flow and cerebral blood volume fMRI. In this review, we describe empirical evidence primarily from our own studies on how well these fMRI signals are spatially specific to the neuronally active site and discuss insights into neurovascular coupling at the mesoscale.This article is part of the theme issue ‘Key relationships between non-invasive functional neuroimaging and the underlying neuronal activity’.     

A model of neurovascular coupling and the BOLD response: PART I

The mechanisms with which neurons communicate with the vasculature to increase blood flow, termed neurovascular coupling is still unclear primarily due to the complex interactions between many parameters and the difficulty in accessing, monitoring and measuring them in the highly heterogeneous brain. Hence a solid theoretical framework based on existing experimental knowledge is necessary to study the relation between neural activity, the associated vasoactive factors released and their effects on the vasculature. Such a framework should also be related to experimental data so that it can be validated against repetitive experiments and generate verifiable hypothesis. We have developed a mathematical model which describes a signaling mechanism of neurovascular coupling with a model of pyramidal neuron and its corresponding fMRI BOLD response. In the first part of two papers we describe the integration of the neurovascular coupling unit extended to include a complex neuron model, which includes the important Na/K ATPase pump, with a model that provides a BOLD signal taking its input from the cerebral blood flow and the metabolic rate of oxygen consumption. We show that this produces a viable signal in terms of initial dip, positive and negative BOLD signals.     

A quantitative model for human neurovascular coupling with translated mechanisms from animals

Neurons regulate the activity of blood vessels through the neurovascular coupling (NVC). A detailed understanding of the NVC is critical for understanding data from functional imaging techniques of the brain. Many aspects of the NVC have been studied both experimentally and using mathematical models; various combinations of blood volume and flow, local field potential (LFP), hemoglobin level, blood oxygenation level-dependent response (BOLD), and optogenetics have been measured and modeled in rodents, primates, or humans. However, these data have not been brought together into a unified quantitative model. We now present a mathematical model that describes all such data types and that preserves mechanistic behaviors between experiments. For instance, from modeling of optogenetics and microscopy data in mice, we learn cell-specific contributions; the first rapid dilation in the vascular response is caused by NO-interneurons, the main part of the dilation during longer stimuli is caused by pyramidal neurons, and the post-peak undershoot is caused by NPY-interneurons. These insights are translated and preserved in all subsequent analyses, together with other insights regarding hemoglobin dynamics and the LFP/BOLD-interplay, obtained from other experiments on rodents and primates. The model can predict independent validation-data not used for training. By bringing together data with complementary information from different species, we both understand each dataset better, and have a basis for a new type of integrative analysis of human data.     

星形胶质细胞和神经元之间谷氨酸-谷氨酰胺的代谢偶联

谷氨酸-谷氨酰胺循环是星形胶质细胞和神经元代谢偶联最重要的途径之一。在中枢神经系统中葡萄糖经糖酵解和三羧酸循环,合成三羧酸循环的中间产物。神经元因缺乏丙酮酸羧化酶,不能由葡萄糖直接合成谷氨酸,而必须依赖于星形胶质细胞的三羧酸循环来产生作为谷氨酸前体的三羧酸循环中间代谢产物。星形胶质细胞的谷氨酸载体从突触间隙摄取谷氨酸,在星形胶质细胞中转变成谷氨酰胺并释放到细胞外,然后重新被神经元摄取,转变成谷氨酸进入新一轮的循环。本文介绍了该循环,以及星形胶质细胞谷氨酸载体的功能、特性及调控。