Oxidative neuroenergetics in event-related paradigms

Thermodynamic limitations on brain oxygen metabolism: physiological implications

Recent thermodynamic modelling indicates that maintaining the brain tissue ratio of O2 to CO2 (abbreviated tissue O2 /CO2 ) is critical for preserving the entropy increase available from oxidative metabolism of glucose, with a fall of that available entropy leading to a reduction of the phosphorylation potential and impairment of brain energy metabolism. This provides a novel perspective for understanding physiological responses under different conditions in terms of preserving tissue O2 /CO2 . To enable estimation of tissue O2 /CO2 in the human brain, a detailed mathematical model of O2 and CO2 transport was developed, and applied to reported physiological responses to different challenges, asking: how well is tissue O2 /CO2 preserved? Reported experimental results for increased neural activity, hypercapnia and hypoxia due to high altitude are consistent with preserving tissue O2 /CO2 . The results highlight two physiological mechanisms that control tissue O2 /CO2 : cerebral blood flow, which modulates tissue O2 ; and ventilation rate, which modulates tissue CO2 . The hypoxia modelling focused on humans at high altitude, including acclimatized lowlanders and Tibetan and Andean adapted populations, with a primary finding that decreasing CO2 by increasing ventilation rate is more effective for preserving tissue O2 /CO2 than increasing blood haemoglobin content to maintain O2 delivery to tissue. This work focused on the function served by particular physiological responses, and the underlying mechanisms require further investigation. The modelling provides a new framework and perspective for understanding how blood flow and other physiological factors support energy metabolism in the brain under a wide range of conditions. KEY POINTS: Thermodynamic modelling indicates that preserving the O2 /CO2 ratio in brain tissue is critical for preserving the entropy change available from oxidative metabolism of glucose and the phosphorylation potential underlying energy metabolism. A detailed model of O2 and CO2 transport was developed to allow estimation of the tissue O2 /CO2 ratio in the human brain in different physiological states. Reported experimental results during hypoxia, hypercapnia and increased oxygen metabolic rate in response to increased neural activity are consistent with maintaining brain tissue O2 /CO2 ratio. The hypoxia modelling of high-altitude acclimatization and adaptation in humans demonstrates the critical role of reducing CO2 with increased ventilation for preserving tissue O2 /CO2 . Preservation of tissue O2 /CO2 provides a novel perspective for understanding the function of observed physiological responses under different conditions in terms of preserving brain energy metabolism, although the mechanisms underlying these functions are not well understood.

Mapping oxidative metabolism in the human brain with calibrated fMRI in health and disease

Conventional functional MRI (fMRI) with blood-oxygenation level dependent (BOLD) contrast is an important tool for mapping human brain activity non-invasively. Recent interest in quantitative fMRI has renewed the importance of oxidative neuroenergetics as reflected by cerebral metabolic rate of oxygen consumption (CMR_O2) to support brain function. Dynamic CMR_O2 mapping by calibrated fMRI require multi-modal measurements of BOLD signal along with cerebral blood flow (CBF) and/or volume (CBV). In human subjects this "calibration" is typically performed using a gas mixture containing small amounts of carbon dioxide and/or oxygen-enriched medical air, which are thought to produce changes in CBF (and CBV) and BOLD signal with minimal or no CMR_O2 changes. However non-human studies have demonstrated that the "calibration" can also be achieved without gases, revealing good agreement between CMR_O2 changes and underlying neuronal activity (e.g., multi-unit activity and local field potential). Given the simpler set-up of gas-free calibrated fMRI, there is evidence of recent clinical applications for this less intrusive direction. This up-to-date review emphasizes technological advances for such translational gas-free calibrated fMRI experiments, also covering historical progression of the calibrated fMRI field that is impacting neurological and neurodegenerative investigations of the human brain.

Thalamic activations in rat brain by fMRI during tactile (forepaw, whisker) and non-tactile (visual, olfactory) sensory stimulations

The thalamus is a crucial subcortical hub that impacts cortical activity. Tracing experiments in animals and post-mortem humans suggest rich morphological specificity of the thalamus. Very few studies reported rodent thalamic activations by functional MRI (fMRI) as compared to cortical activations for different sensory stimuli. Here, we show different portions of the rat thalamus in response to tactile (forepaw, whisker) and non-tactile (visual, olfactory) sensory stimuli with high field fMRI (11.7T) using a custom-build quadrature surface coil to capture high sensitivity signals from superficial and deep brain regions simultaneously. Results demonstrate reproducible thalamic activations during both tactile and non-tactile stimuli. Forepaw and whisker stimuli activated broader regions within the thalamus: ventral posterior lateral (VPL), ventral posterior medial (VPM), lateral posterior mediorostral (LPMR) and posterior medial (POm) thalamic nuclei. Visual stimuli activated dorsal lateral geniculate nucleus (DLG) of the thalamus but also parts of the superior/inferior colliculus, whereas olfactory stimuli activated specifically the mediodorsal nucleus of the thalamus (MDT). BOLD activations in LGN and MDT were much stronger than in VPL, VPM, LPMR and POm. These fMRI-based thalamic activations suggest that forepaw and whisker (i.e., tactile) stimuli engage VPL, VPM, LPMR and POm whereas visual and olfactory (i.e., non-tactile) stimuli, respectively, recruit DLG and MDT exclusively.

High-resolution relaxometry-based calibrated fMRI in murine brain: Metabolic differences between awake and anesthetized states

Functional magnetic resonance imaging (fMRI) techniques using the blood-oxygen level-dependent (BOLD) signal have shown great potential as clinical biomarkers of disease. Thus, using these techniques in preclinical rodent models is an urgent need. Calibrated fMRI is a promising technique that can provide high-resolution mapping of cerebral oxygen metabolism (CMR_O2). However, calibrated fMRI is difficult to use in rodent models for several reasons: rodents are anesthetized, stimulation-induced changes are small, and gas challenges induce noisy CMR_O2 predictions. We used, in mice, a relaxometry-based calibrated fMRI method which uses cerebral blood flow (CBF) and the BOLD-sensitive magnetic relaxation component, R₂', the same parameter derived in the deoxyhemoglobin-dilution model of calibrated fMRI. This method does not use any gas challenges, which we tested on mice in both awake and anesthetized states. As anesthesia induces a whole-brain change, our protocol allowed us to overcome the former limitations of rodent studies using calibrated fMRI. We revealed 1.5-2 times higher CMR_O2, dependent upon brain region, in the awake state versus the anesthetized state. Our results agree with alternative measurements of whole-brain CMR_O2 in the same mice and previous human anesthesia studies. The use of calibrated fMRI in rodents has much potential for preclinical fMRI.

Functional enrichment of alternative splicing events with NEASE reveals insights into tissue identity and diseases

Alternative splicing (AS) is an important aspect of gene regulation. Nevertheless, its role in molecular processes and pathobiology is far from understood. A roadblock is that tools for the functional analysis of AS-set events are lacking. To mitigate this, we developed NEASE, a tool integrating pathways with structural annotations of protein-protein interactions to functionally characterize AS events. We show in four application cases how NEASE can identify pathways contributing to tissue identity and cell type development, and how it highlights splicing-related biomarkers. With a unique view on AS, NEASE generates unique and meaningful biological insights complementary to classical pathways analysis.

Supraspinal sensorimotor and pain-related reorganization after a hemicontusion rat cervical spinal cord injury

Since the presence of pain impedes motor recovery in individuals with spinal cord injury (SCI), it is necessary to understand their supraspinal substrates in translational animal models. Using functional magnetic resonance imaging (fMRI) in a rat model of hemicontusion cervical SCI, supraspinal changes were mapped and correlated with sensorimotor behavioral outcomes. Female adult rats underwent sham or SCI using a 2.5 mm impactor and 150 kDyne force. SCI permanently impaired motor activity in only the ipsilesional forelimb along with thermal hyperalgesia at 5 and 6 wks. Spinal MRI at 8 wks after SCI showed ipsilateral T1 and T2 lesions with no discernable lesions across shams. fMRI mapping during electrical forepaw stimulation indicated SCI-induced sensorimotor reorganization with an expansion of the contralesional forelimb representation. Resting state fMRI based functional connectivity density (FCD), a marker of regional neuronal hubs increased or decreased across brain regions involved in nociception. FCD increases after SCI were in the primary and secondary somatosensory cortices (S1 and S2), anterior cingulate cortex (ACC), insula and the prefrontal cortex (PFC) and decreases were across the hippocampus, thalamus, hypothalamus and amygdala in SCI. Resting state functional connectivity (RSFC) assessments from the FCD altered regions of interest indicated cortico-cortical RSFC increases and cortico-insular, cortico-thalamic and cortico-hypothalamic RSFC decreases after SCI. Hippocampus, amygdala and thalamus showed decreased RSFC with most cortical regions and between themselves except the hippocampus-amygdala network, which showed increased RSFC after SCI. While select nociceptive region's intrinsic activity associated strongly with evoked pain behaviors after SCI (eg., PFC, ACC, hippocampus, thalamus, hypothalamus, M1 and S1BF) other nociceptive regions had weaker associations (eg., amygdala, insula, auditory cortex, S1FL, S1HL, S2 and M2), but differed significantly in their intrinsic activities between sham and SCI. The weaker associated nociceptive regions may possibly encode both the evoked and affective components of pain.

The neural-vascular basis of age-related processing speed decline

Most studies examining neurocognitive aging are based on the blood-oxygen level-dependent signal obtained during functional magnetic resonance imaging (fMRI). The physiological basis of this signal is neural-vascular coupling, the process by which neurons signal cerebrovasculature to dilate in response to an increase in active neural metabolism due to stimulation. These fMRI studies of aging rely on the hemodynamic equivalence assumption that this process is not disrupted by physiologic deterioration associated with aging. Studies of neural-vascular coupling challenge this assumption and show that neural-vascular coupling is closely related to cognition. In this review, we put forward a theory of processing speed decline in aging and how it is related to age-related neural-vascular coupling changes based on the results of studies elucidating the relationships between cognition, cerebrovascular dynamics, and aging.

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

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

Local oxygen homeostasis during various neuronal network activity states in the mouse hippocampus

Cortical information processing comprises various activity states emerging from timed synaptic excitation and inhibition. However, the underlying energy metabolism is widely unknown. We determined the cerebral metabolic rate of oxygen (CMRO2) along a tissue depth of <0.3 mm in the hippocampal CA3 region during various network activities, including gamma oscillations and sharp wave-ripples that occur during wakefulness and sleep. These physiological states associate with sensory perception and memory formation, and critically depend on perisomatic GABA inhibition. Moreover, we modelled vascular oxygen delivery based on quantitative microvasculature analysis. (1) Local CMRO2 was highest during gamma oscillations (3.4 mM/min), medium during sharp wave-ripples, asynchronous activity and isoflurane application (2.0-1.6 mM/min), and lowest during tetrodotoxin application (1.4 mM/min). (2) Energy expenditure of axonal and synaptic signaling accounted for >50% during gamma oscillations. (3) CMRO2 positively correlated with number and synchronisation of activated synapses, and neural multi-unit activity. (4) The median capillary distance was 44 µm. (5) The vascular oxygen partial pressure of 33 mmHg was needed to sustain oxidative phosphorylation during gamma oscillations. We conclude that gamma oscillations featuring high energetics require a hemodynamic response to match oxygen consumption of respiring mitochondria, and that perisomatic inhibition significantly contributes to the brain energy budget.

Effective Connectivity Measured Using Optogenetically Evoked Hemodynamic Signals Exhibits Topography Distinct from Resting State Functional Connectivity in the Mouse

Brain connectomics has expanded from histological assessment of axonal projection connectivity (APC) to encompass resting state functional connectivity (RS-FC). RS-FC analyses are efficient for whole-brain mapping, but attempts to explain aspects of RS-FC (e.g., interhemispheric RS-FC) based on APC have been only partially successful. Neuroimaging with hemoglobin alone lacks specificity for determining how activity in a population of cells contributes to RS-FC. Wide-field mapping of optogenetically defined connectivity could provide insights into the brain's structure-function relationship. We combined optogenetics with optical intrinsic signal imaging to create an efficient, optogenetic effective connectivity (Opto-EC) mapping assay. We examined EC patterns of excitatory neurons in awake, Thy1-ChR2 transgenic mice. These Thy1-based EC (Thy1-EC) patterns were evaluated against RS-FC over the cortex. Compared to RS-FC, Thy1-EC exhibited increased spatial specificity, reduced interhemispheric connectivity in regions with strong RS-FC, and appreciable connection strength asymmetry. Comparing the topography of Thy1-EC and RS-FC patterns to maps of APC revealed that Thy1-EC more closely resembled APC than did RS-FC. The more general method of Opto-EC mapping with hemoglobin can be determined for 100 sites in single animals in under an hour, and is amenable to other neuroimaging modalities. Opto-EC mapping represents a powerful strategy for examining evolving connectivity-related circuit plasticity.

Evaluating the gray and white matter energy budgets of human brain function

The insatiable appetite for energy to support human brain function is mainly supplied by glucose oxidation (CMR_glc(ox)). But how much energy is consumed for signaling and nonsignaling processes in gray/white matter is highly debated. We examined this issue by combining metabolic measurements of gray/white matter and a theoretical calculation of bottom-up energy budget using biophysical properties of neuronal/glial cells in conjunction with species-exclusive electrophysiological and morphological data. We calculated a CMR_glc(ox)-derived budget and confirmed it with experimental results measured by PET, autoradiography, ¹³C-MRS, and electrophysiology. Several conserved principles were observed regarding the energy costs for brain's signaling and nonsignaling components in both human and rat. The awake resting cortical signaling processes and mass-dependent nonsignaling processes, respectively, demand ∼70% and ∼30% of CMR_glc(ox). Inhibitory neurons and glia need 15-20% of CMR_glc(ox), with the rest demanded by excitatory neurons. Nonsignaling demands dominate in white matter, in near opposite contrast to gray matter demands. Comparison between ¹³C-MRS data and calculations suggests ∼1.2 Hz glutamatergic signaling rate in the awake human cortex, which is ∼4 times lower than signaling in the rat cortex. Top-down validated bottom-up budgets could allow computation of anatomy-based CMR_glc(ox) maps and accurate cellular level interpretation of brain metabolic imaging.

Advances in Imaging Brain Metabolism

Metabolism is central to neuroimaging because it can reveal pathways by which neuronal and glial cells use nutrients to fuel their growth and function. We focus on advanced magnetic resonance imaging (MRI) and magnetic resonance spectroscopy (MRS) methods used in brain metabolic studies. ¹⁷O-MRS and ³¹P-MRS, respectively, provide rates of oxygen use and ATP synthesis inside mitochondria, whereas ¹⁹F-MRS enables measurement of cytosolic glucose metabolism. Calibrated functional MRI (fMRI), an advanced form of fMRI that uses contrast generated by deoxyhemoglobin, provides maps of oxygen use that track neuronal firing across brain regions. ¹³C-MRS is the only noninvasive method of measuring both glutamatergic neurotransmission and cell-specific energetics with signaling and nonsignaling purposes. Novel MRI contrasts, arising from endogenous diamagnetic agents and exogenous paramagnetic agents, permit pH imaging of glioma. Overall, these magnetic resonance methods for imaging brain metabolism demonstrate translational potential to better understand brain disorders and guide diagnosis and treatment.

Volumetric Spatial Correlations of Neurovascular Coupling Studied using Single Pulse Opto-fMRI

Neurovascular coupling describes the link between neuronal activity and cerebral blood flow. This relationship has been the subject of intense scrutiny, with most previous work seeking to understand temporal correlations that describe neurovascular coupling. However, to date, the study of spatial correlations has been limited to two-dimensional mapping of neuronal or vascular derived signals emanating from the brain’s surface, using optical imaging techniques. Here, we investigate spatial correlations of neurovascular coupling in three dimensions, by applying a single 10 ms pulse of light to trigger optogenetic activation of cortical neurons transduced to express channelrhodopsin2, with concurrent fMRI. We estimated the spatial extent of increased neuronal activity using a model that takes into the account the scattering and absorption of blue light in brain tissue together with the relative density of channelrhodopsin2 expression across cortical layers. This method allows precise modulation of the volume of activated tissue in the cerebral cortex with concurrent three-dimensional mapping of functional hyperemia. Single pulse opto-fMRI minimizes adaptation, avoids heating artefacts and enables confined recruitment of the neuronal activity. Using this novel method, we present evidence for direct proportionality of volumetric spatial neurovascular coupling in the cerebral cortex.

Mitochondrial calcium homeostasis: Implications for neurovascular and neurometabolic coupling

Mitochondrial function is critical to maintain high rates of oxidative metabolism supporting energy demands of both spontaneous and evoked neuronal activity in the brain. Mitochondria not only regulate energy metabolism, but also influence neuronal signaling. Regulation of "energy metabolism" and "neuronal signaling" (i.e. neurometabolic coupling), which are coupled rather than independent can be understood through mitochondria's integrative functions of calcium ion (Ca²⁺) uptake and cycling. While mitochondrial Ca²⁺ do not affect hemodynamics directly, neuronal activity changes are mechanistically linked to functional hyperemic responses (i.e. neurovascular coupling). Early in vitro studies lay the foundation of mitochondrial Ca²⁺ homeostasis and its functional roles within cells. However, recent in vivo approaches indicate mitochondrial Ca²⁺ homeostasis as maintained by the role of mitochondrial Ca²⁺ uniporter (mCU) influences system-level brain activity as measured by a variety of techniques. Based on earlier evidence of subcellular cytoplasmic Ca²⁺ microdomains and cellular bioenergetic states, a mechanistic model of Ca²⁺ mobilization is presented to understand systems-level neurovascular and neurometabolic coupling. This integrated view from molecular and cellular to the systems level, where mCU plays a major role in mitochondrial and cellular Ca²⁺ homeostasis, may explain the wide range of activation-induced coupling across neuronal activity, hemodynamic, and metabolic responses.

Metabolic demands of neural-hemodynamic associated and disassociated areas in brain

Interpretation of regional blood oxygenation level-dependent (BOLD) responses in functional magnetic resonance imaging (fMRI) is contingent on whether local field potential (LFP) and multi-unit activity (MUA) is either dissociated or associated. To examine whether neural-hemodynamic associated and dissociated areas have different metabolic demands, we recorded sensory-evoked responses of BOLD signal, blood flow (CBF), and blood volume (CBV), which with calibrated fMRI provided oxidative metabolism (CMRO2) from rat's ventral posterolateral thalamic nucleus (VPL) and somatosensory forelimb cortex (S1FL) and compared these neuroimaging signals to neurophysiological recordings. MUA faithfully recorded evoked latency differences between VPL and S1FL because evoked MUA in these regions were similar in magnitude. Since evoked LFP was significantly attenuated in VPL, we extracted the time courses of the weaker thalamic LFP to compare with the stronger cortical LFP using wavelet transform. BOLD and CBV responses were greater in S1FL than in VPL, similar to LFP regional differences. CBF and CMRO2 responses were both comparably larger in S1FL and VPL. Despite different levels of CBF-CMRO2 and LFP-MUA couplings in VPL and S1FL, the CMRO2 was well matched with MUA in both regions. These results suggest that neural-hemodynamic associated and dissociated areas in VPL and S1FL can have similar metabolic demands.

New horizons in neurometabolic and neurovascular coupling from calibrated fMRI

Neurovascular coupling relates changes in neuronal activity to constriction/dilation of microvessels. However neurometabolic coupling, which is less well known, relates alterations in neuronal activity with metabolic demands. The link between the blood oxygenation level dependent (BOLD) signal and neural activity opened doors for functional MRI (fMRI) to be a powerful neuroimaging tool in the neurosciences. But due to the complex makeup of BOLD contrast, researchers began to investigate the relationship between BOLD signal and blood flow and/or volume changes during functional brain activation, which together provided the tools to measure oxygen consumption on the basis of the biophysical model of BOLD. This field is called calibrated fMRI, thereby allowed probing of both neurometabolic and neurovascular couplings for a variety of health conditions in animals and humans. Calibrated fMRI may provide brain disorder biomarkers that could be used for monitoring effective therapies.

Comparison of glomerular activity patterns by fMRI and wide-field calcium imaging: Implications for principles underlying odor mapping

Functional imaging signals arise from distinct metabolic and hemodynamic events at the neuropil, but how these processes are influenced by pre- and post-synaptic activities need to be understood for quantitative interpretation of stimulus-evoked mapping data. The olfactory bulb (OB) glomeruli, spherical neuropil regions with well-defined neuronal circuitry, can provide insights into this issue. Optical calcium-sensitive fluorescent dye imaging (OICa(2+)) reflects dynamics of pre-synaptic input to glomeruli, whereas high-resolution functional magnetic resonance imaging (fMRI) using deoxyhemoglobin contrast reveals neuropil function within the glomerular layer where both pre- and post-synaptic activities contribute. We imaged odor-specific activity patterns of the dorsal OB in the same anesthetized rats with fMRI and OICa(2+) and then co-registered the respective maps to compare patterns in the same space. Maps by each modality were very reproducible as trial-to-trial patterns for a given odor, overlapping by ~80%. Maps evoked by ethyl butyrate and methyl valerate for a given modality overlapped by ~80%, suggesting activation of similar dorsal glomerular networks by these odors. Comparison of maps generated by both methods for a given odor showed ~70% overlap, indicating similar odor-specific maps by each method. These results suggest that odor-specific glomerular patterns by high-resolution fMRI primarily tracks pre-synaptic input to the OB. Thus combining OICa(2+) and fMRI lays the framework for studies of OB processing over a range of spatiotemporal scales, where OICa(2+) can feature the fast dynamics of dorsal glomerular clusters and fMRI can map the entire glomerular sheet in the OB.

Role of mitochondrial calcium uptake homeostasis in resting state fMRI brain networks

Mitochondrial Ca(2+) uptake influences both brain energy metabolism and neural signaling. Given that brain mitochondrial organelles are distributed in relation to vascular density, which varies considerably across brain regions, we hypothesized different physiological impacts of mitochondrial Ca(2+) uptake across brain regions. We tested the hypothesis by monitoring brain "intrinsic activity" derived from the resting state functional MRI (fMRI) blood oxygen level dependent (BOLD) fluctuations in different functional networks spanning the somatosensory cortex, caudate putamen, hippocampus and thalamus, in normal and perturbed mitochondrial Ca(2+) uptake states. In anesthetized rats at 11.7 T, mitochondrial Ca(2+) uptake was inhibited or enhanced respectively by treatments with Ru360 or kaempferol. Surprisingly, mitochondrial Ca(2+) uptake inhibition by Ru360 and enhancement by kaempferol led to similar dose-dependent decreases in brain-wide intrinsic activities in both the frequency domain (spectral amplitude) and temporal domain (resting state functional connectivity; RSFC). The fact that there were similar dose-dependent decreases in the frequency and temporal domains of the resting state fMRI-BOLD fluctuations during mitochondrial Ca(2+) uptake inhibition or enhancement indicated that mitochondrial Ca(2+) uptake and its homeostasis may strongly influence the brain's functional organization at rest. Interestingly, the resting state fMRI-derived intrinsic activities in the caudate putamen and thalamic regions saturated much faster with increasing dosage of either drug treatment than the drug-induced trends observed in cortical and hippocampal regions. Regional differences in how the spectral amplitude and RSFC changed with treatment indicate distinct mitochondrion-mediated spontaneous neuronal activity coupling within the various RSFC networks determined by resting state fMRI.

Brain region and activity-dependent properties of M for calibrated fMRI

Calibrated fMRI extracts changes in oxidative energy demanded by neural activity based on hemodynamic and metabolic dependencies of the blood oxygenation level-dependent (BOLD) response. This procedure requires the parameter M, which is determined from the dynamic range of the BOLD signal between deoxyhemoglobin (paramagnetic) and oxyhemoglobin (diamagnetic). Since it is unclear if the range of M-values in human calibrated fMRI is due to regional/state differences, we conducted a 9.4T study to measure M-values across brain regions in deep (α-chloralose) and light (medetomidine) anesthetized rats, as verified by electrophysiology. Because BOLD signal is captured differentially by gradient-echo (R2*) and spin-echo (R2) relaxation rates, we measured M-values by the product of the fMRI echo time and R2' (i.e., the reversible magnetic susceptibility component), which is given by the absolute difference between R2* and R2. While R2' mapping was shown to be dependent on the k-space sampling method used, at nominal spatial resolutions achieved at high magnetic field of 9.4T the M-values were quite homogenous across cortical gray matter. However cortical M-values varied in relation to neural activity between brain states. The findings from this study could improve precision of future calibrated fMRI studies by focusing on the global uniformity of M-values in gray matter across different resting activity levels.

Neural and hemodynamic responses to optogenetic and sensory stimulation in the rat somatosensory cortex

Introducing optogenetics into neurovascular research can provide novel insights into the cell-specific control of the hemodynamic response. To generalize findings from molecular approaches, it is crucial to determine whether light-activated circuits have the same effect on the vasculature as sensory-activated ones. For that purpose, rats expressing channelrhodopsin (ChR2) specific to excitatory glutamatergic neurons were used to measure neural activity, blood flow, hemoglobin-based optical intrinsic signal, and blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging (fMRI) during optogenetic and sensory stimulation. The magnitude of the evoked hemodynamic responses was monotonically correlated with optogenetic stimulus strength. The BOLD hemodynamic response function was consistent for optogenetic and sensory stimuli. The relationship between electrical activities and hemodynamic responses was comparable for optogenetic and sensory stimuli, and better explained by the local field potential (LFP) than the firing rate. The LFP was well correlated with cerebral blood flow, moderately with cerebral blood volume, and less with deoxyhemoglobin (dHb) level. The presynaptic firing rate had little impact on evoking vascular response. Contribution of the postsynaptic LFP to the blood flow response induced by optogenetic stimulus was further confirmed by the application of glutamate receptor antagonists. Overall, neurovascular coupling during optogenetic control of glutamatergic neurons largely conforms to that of a sensory stimulus.

Imaging of prolonged BOLD response in the somatosensory cortex of the rat

Blood oxygenation level-dependent (BOLD) functional MRI is a widely employed methodology in experimental and clinical neuroscience, although its nature is not fully understood. To gain insights into BOLD mechanisms and take advantage of the new functional methods, it is of interest to investigate prolonged paradigms of activation suitable for long experimental protocols and to observe any long-term modifications induced by these functional challenges. While different types of sustained stimulation paradigm have been explored in human studies, the BOLD response is typically limited to a few minutes in animal models, due to fatigue, anesthesia effects and physiological instability. In the present study, the rat forepaw was electrically stimulated for 2 h, which resulted in a prolonged and localized cortical BOLD response over that period. The stimulation paradigm, including an inter-stimulus interval (ISI) of 10 s, that is 25% of the total time, was applied at constant or variable frequency over 2 h. The steady-state level of the BOLD response was reached after 15-20 min of stimulation and was maintained until the end of the stimulation. On average, no substantial loss in activated volume was observed at the end of the stimulation, but less variability in the fraction of remaining activated volume and higher steady-state BOLD amplitude were observed when stimulation frequency was varied between 2 and 3 Hz every 5 min. We conclude that the combination of ISI and variable stimulus frequency reproducibly results in robust, prolonged and localized BOLD activation.

Insights from neuroenergetics into the interpretation of functional neuroimaging: an alternative empirical model for studying the brain's support of behavior

Functional neuroimaging measures quantitative changes in neurophysiological parameters coupled to neuronal activity during observable behavior. These results have usually been interpreted by assuming that mental causation of behavior arises from the simultaneous actions of distinct psychological mechanisms or modules. However, reproducible localization of these modules in the brain using functional magnetic resonance imaging (MRI) and positron emission tomography (PET) imaging has been elusive other than for sensory systems. In this paper, we show that neuroenergetic studies using PET, calibrated functional magnetic resonance imaging (fMRI), (13)C magnetic resonance spectroscopy, and electrical recordings do not support the standard approach, which identifies the location of mental modules from changes in brain activity. Of importance in reaching this conclusion is that changes in neuronal activities underlying the fMRI signal are many times smaller than the high ubiquitous, baseline neuronal activity, or energy in resting, awake humans. Furthermore, the incremental signal depends on the baseline activity contradicting theoretical assumptions about linearity and insertion of mental modules. To avoid these problems, while making use of these valuable results, we propose that neuroimaging should be used to identify observable brain activities that are necessary for a person's observable behavior rather than being used to seek hypothesized mental processes.

Purinergic Modulation of Spinal Neuroglial Maladaptive Plasticity Following Peripheral Nerve Injury

Modulation of spinal reactive gliosis following peripheral nerve injury (PNI) is a promising strategy to restore synaptic homeostasis. Oxidized ATP (OxATP), a nonselective antagonist of purinergic P2X receptors, was found to recover a neuropathic behavior following PNI. We investigated the role of intraperitoneal (i.p.) OxATP treatment in restoring the expression of neuronal and glial markers in the mouse spinal cord after sciatic spared nerve injury (SNI). Using in vivo two-photon microscopy, we imaged Ca(2+) transients in neurons and astrocytes of the dorsal horn of spinal cord at rest and upon right hind paw electrical stimulation in sham, SNI, and OxATP-treated mice. Neuropathic behavior was investigated by von Frey and thermal plantar test. Glial [glial fibrillary acidic protein (GFAP), ionized calcium-binding adaptor molecule 1 (Iba1)] and GABAergic [vesicular GABA transporter (vGAT) and glutamic acid decarboxylase 65/76 (GAD65/67)] markers and glial [glutamate transporter (GLT1) and GLAST] and neuronal amino acid [EAAC1, vesicular glutamate transporter 1 (vGLUT1)] transporters have been evaluated. In SNI mice, we found (i) increased glial response, (ii) decreased glial amino acid transporters, and (iii) increased levels of neuronal amino acid transporters, and (iv) in vivo analysis of spinal neurons and astrocytes showed a persistent increase of Ca(2+) levels. OxATP administration reduced glial activation, modulated the expression of glial and neuronal glutamate/GABA transporters, restored neuronal and astrocytic Ca(2+) levels, and prevented neuropathic behavior. In vitro studies validated that OxATP (i) reduced levels of reactive oxygen species (ROS), (ii) reduced astrocytic proliferation, (iii) increase vGLUT expression. All together, these data support the correlation between reactive gliosis and perturbation of the spinal synaptic homeostasis and the role played by the purinergic system in modulating spinal plasticity following PNI.

Quantitative basis for neuroimaging of cortical laminae with calibrated functional MRI

Layer-specific neurophysiologic, hemodynamic, and metabolic measurements are needed to interpret high-resolution functional magnetic resonance imaging (fMRI) data in the cerebral cortex. We examined how neurovascular and neurometabolic couplings vary vertically in the rat's somatosensory cortex. During sensory stimulation we measured dynamic layer-specific responses of local field potential (LFP) and multiunit activity (MUA) as well as blood oxygenation level-dependent (BOLD) signal and cerebral blood volume (CBV) and blood flow (CBF), which in turn were used to calculate changes in oxidative metabolism (CMR(O2)) with calibrated fMRI. Both BOLD signal and CBV decreased from superficial to deep laminae, but these responses were not well correlated with either layer-specific LFP or MUA. However, CBF changes were quite stable across laminae, similar to LFP. However, changes in CMR(O2) and MUA varied across cortex in a correlated manner and both were reduced in superficial lamina. These results lay the framework for quantitative neuroimaging across cortical laminae with calibrated fMRI methods.

Mitochondrial calcium uptake capacity modulates neocortical excitability

Local calcium (Ca(2+)) changes regulate central nervous system metabolism and communication integrated by subcellular processes including mitochondrial Ca(2+) uptake. Mitochondria take up Ca(2+) through the calcium uniporter (mCU) aided by cytoplasmic microdomains of high Ca(2+). Known only in vitro, the in vivo impact of mCU activity may reveal Ca(2+)-mediated roles of mitochondria in brain signaling and metabolism. From in vitro studies of mitochondrial Ca(2+) sequestration and cycling in various cell types of the central nervous system, we evaluated ranges of spontaneous and activity-induced Ca(2+) distributions in multiple subcellular compartments in vivo. We hypothesized that inhibiting (or enhancing) mCU activity would attenuate (or augment) cortical neuronal activity as well as activity-induced hemodynamic responses in an overall cytoplasmic and mitochondrial Ca(2+)-dependent manner. Spontaneous and sensory-evoked cortical activities were measured by extracellular electrophysiology complemented with dynamic mapping of blood oxygen level dependence and cerebral blood flow. Calcium uniporter activity was inhibited and enhanced pharmacologically, and its impact on the multimodal measures were analyzed in an integrated manner. Ru360, an mCU inhibitor, reduced all stimulus-evoked responses, whereas Kaempferol, an mCU enhancer, augmented all evoked responses. Collectively, the results confirm aforementioned hypotheses and support the Ca(2+) uptake-mediated integrative role of in vivo mitochondria on neocortical activity.

Neural and hemodynamic responses elicited by forelimb- and photo-stimulation in channelrhodopsin-2 mice: insights into the hemodynamic point spread function

Hemodynamic responses are commonly used to map brain activity; however, their spatial limits have remained unclear because of the lack of a well-defined and malleable spatial stimulus. To examine the properties of neural activity and hemodynamic responses, multiunit activity, local field potential, cerebral blood volume (CBV)-sensitive optical imaging, and laser Doppler flowmetry were measured from the somatosensory cortex of transgenic mice expressing Channelrhodopsin-2 in cortex Layer 5 pyramidal neurons. The magnitude and extent of neural and hemodynamic responses were modulated using different photo-stimulation parameters and compared with those induced by somatosensory stimulation. Photo-stimulation-evoked spiking activity across cortical layers was similar to forelimb stimulation, although their activity originated in different layers. Hemodynamic responses induced by forelimb- and photo-stimulation were similar in magnitude and shape, although the former were slightly larger in amplitude and wider in extent. Altogether, the neurovascular relationship differed between these 2 stimulation pathways, but photo-stimulation-evoked changes in neural and hemodynamic activities were linearly correlated. Hemodynamic point spread functions were estimated from the photo-stimulation data and its full-width at half-maximum ranged between 103 and 175 µm. Therefore, submillimeter functional structures separated by a few hundred micrometers may be resolved using hemodynamic methods, such as optical imaging and functional magnetic resonance imaging.

Adaptation of cerebral oxygen metabolism and blood flow and modulation of neurovascular coupling with prolonged stimulation in human visual cortex

Prolonged visual stimulation results in neurophysiologic and hemodynamic adaptation. However, the hemodynamic adaptation appears to be small compared to neural adaptation. It is not clear how the cerebral metabolic rate of oxygen (CMRO2) is affected by adaptation. We measured cerebral blood flow (CBF) and CMRO2 change in responses to peripheral stimulation either continuously, or intermittently (on/off cycles). A linear system's response to the continuous input should be equal to the sum of the original response to the intermittent input and a version of that response shifted by half a cycle. The CMRO2 response showed a large non-linearity consistent with adaptation, the CBF response adapted to a lesser degree, and the blood oxygenation level dependent (BOLD) response was nearly linear. The metabolic response was coupled with a larger flow in the continuous condition than in the intermittent condition. Our results suggest that contrast adaptation improves energy economy of visual processing. However BOLD modulations may not accurately represent the underlying metabolic nonlinearity due to modulation of the coupling of blood flow and oxygen metabolism changes.

Functional MRI and neural responses in a rat model of Alzheimer's disease

Based on the hypothesis that brain plaques and tangles can affect cortical function in Alzheimer's disease (AD), we investigated functional responses in an AD rat model (called the Samaritan Alzheimer's rat achieved by ventricular infusion of amyloid peptide) and age-matched healthy control. High-field functional magnetic resonance imaging (fMRI) and extracellular neural activity measurements were applied to characterize sensory-evoked responses. Electrical stimulation of the forepaw led to BOLD and neural responses in the contralateral somatosensory cortex and thalamus. In AD brain we noted much smaller BOLD activation patterns in the somatosensory cortex (i.e., about 50% less activated voxels compared to normal brain). While magnitudes of BOLD and neural responses in the cerebral cortex were markedly attenuated in AD rats compared to normal rats (by about 50%), the dynamic coupling between the BOLD and neural responses in the cerebral cortex, as assessed by transfer function analysis, remained unaltered between the groups. However thalamic BOLD and neural responses were unaltered in AD brain compared to controls. Thus cortical responses in the AD model were indeed diminished compared to controls, but the thalamic responses in the AD and control rats were quite similar. Therefore these results suggest that Alzheimer's disease may affect cortical function more than subcortical function, which may have implications for interpreting altered human brain functional responses in fMRI studies of Alzheimer's disease.

Mitochondrial functional state impacts spontaneous neocortical activity and resting state FMRI

Mitochondrial Ca²⁺ uptake, central to neural metabolism and function, is diminished in aging whereas enhanced after acute/sub-acute traumatic brain injury. To develop relevant translational models for these neuropathologies, we determined the impact of perturbed mitochondrial Ca²⁺ uptake capacities on intrinsic brain activity using clinically relevant markers. From a multi-compartment estimate of probable baseline Ca²⁺ ranges in the brain, we hypothesized that reduced or enhanced mitochondrial Ca²⁺ uptake capacity would decrease or increase spontaneous neuronal activity respectively. As resting state fMRI-BOLD fluctuations and stimulus-evoked BOLD responses have similar physiological origins [1] and stimulus-evoked neuronal and hemodynamic responses are modulated by mitochondrial Ca²⁺ uptake capacity [2], [3] respectively, we tested our hypothesis by measuring hemodynamic fluctuations and spontaneous neuronal activities during normal and altered mitochondrial functional states. Mitochondrial Ca²⁺ uptake capacity was perturbed by pharmacologically inhibiting or enhancing the mitochondrial Ca²⁺ uniporter (mCU) activity. Neuronal electrical activity and cerebral blood flow (CBF) fluctuations were measured simultaneously and integrated with fMRI-BOLD fluctuations at 11.7T. mCU inhibition reduced spontaneous neuronal activity and the resting state functional connectivity (RSFC), whereas mCU enhancement increased spontaneous neuronal activity but reduced RSFC. We conclude that increased or decreased mitochondrial Ca²⁺ uptake capacities lead to diminished resting state modes of brain functional connectivity.

Lactate preserves neuronal metabolism and function following antecedent recurrent hypoglycemia

Hypoglycemia occurs frequently during intensive insulin therapy in patients with both type 1 and type 2 diabetes and remains the single most important obstacle in achieving tight glycemic control. Using a rodent model of hypoglycemia, we demonstrated that exposure to antecedent recurrent hypoglycemia leads to adaptations of brain metabolism so that modest increments in circulating lactate allow the brain to function normally under acute hypoglycemic conditions. We characterized 3 major factors underlying this effect. First, we measured enhanced transport of lactate both into as well as out of the brain that resulted in only a small increase of its contribution to total brain oxidative capacity, suggesting that it was not the major fuel. Second, we observed a doubling of the glucose contribution to brain metabolism under hypoglycemic conditions that restored metabolic activity to levels otherwise only observed at euglycemia. Third, we determined that elevated lactate is critical for maintaining glucose metabolism under hypoglycemia, which preserves neuronal function. These unexpected findings suggest that while lactate uptake was enhanced, it is insufficient to support metabolism as an alternate substrate to replace glucose. Lactate is, however, able to modulate metabolic and neuronal activity, serving as a "metabolic regulator" instead.

Role of ongoing, intrinsic activity of neuronal populations for quantitative neuroimaging of functional magnetic resonance imaging-based networks

A primary objective in neuroscience is to determine how neuronal populations process information within networks. In humans and animal models, functional magnetic resonance imaging (fMRI) is gaining increasing popularity for network mapping. Although neuroimaging with fMRI-conducted with or without tasks-is actively discovering new brain networks, current fMRI data analysis schemes disregard the importance of the total neuronal activity in a region. In task fMRI experiments, the baseline is differenced away to disclose areas of small evoked changes in the blood oxygenation level-dependent (BOLD) signal. In resting-state fMRI experiments, the spotlight is on regions revealed by correlations of tiny fluctuations in the baseline (or spontaneous) BOLD signal. Interpretation of fMRI-based networks is obscured further, because the BOLD signal indirectly reflects neuronal activity, and difference/correlation maps are thresholded. Since the small changes of BOLD signal typically observed in cognitive fMRI experiments represent a minimal fraction of the total energy/activity in a given area, the relevance of fMRI-based networks is uncertain, because the majority of neuronal energy/activity is ignored. Thus, another alternative for quantitative neuroimaging of fMRI-based networks is a perspective in which the activity of a neuronal population is accounted for by the demanded oxidative energy (CMR(O2)). In this article, we argue that network mapping can be improved by including neuronal energy/activity of both the information about baseline and small differences/fluctuations of BOLD signal. Thus, total energy/activity information can be obtained through use of calibrated fMRI to quantify differences of ΔCMR(O2) and through resting-state positron emission tomography/magnetic resonance spectroscopy measurements for average CMR(O2).

Cortical energy demands of signaling and nonsignaling components in brain are conserved across mammalian species and activity levels

The continuous need for ion gradient restoration across the cell membrane, a prerequisite for synaptic transmission and conduction, is believed to be a major factor for brain's high oxidative demand. However, do energy requirements of signaling and nonsignaling components of cortical neurons and astrocytes vary with activity levels and across species? We derived oxidative ATP demand associated with signaling (P(s)) and nonsignaling (P(ns)) components in the cerebral cortex using species-specific physiologic and anatomic data. In rat, we calculated glucose oxidation rates from layer-specific neuronal activity measured across different states, spanning from isoelectricity to awake and sensory stimulation. We then compared these calculated glucose oxidation rates with measured glucose metabolic data for the same states as reported by 2-deoxy-glucose autoradiography. Fixed values for P(s) and P(ns) were able to predict the entire range of states in the rat. We then calculated glucose oxidation rates from human EEG data acquired under various conditions using fixed P(s) and P(ns) values derived for the rat. These calculated metabolic data in human cerebral cortex compared well with glucose metabolism measured by PET. Independent of species, linear relationship was established between neuronal activity and neuronal oxidative demand beyond isoelectricity. Cortical signaling requirements dominated energy demand in the awake state, whereas nonsignaling requirements were ∼20% of awake value. These predictions are supported by (13)C magnetic resonance spectroscopy results. We conclude that mitochondrial energy support for signaling and nonsignaling components in cerebral cortex are conserved across activity levels in mammalian species.

Calcium imaging of living astrocytes in the mouse spinal cord following sensory stimulation

Astrocytic Ca²⁺ dynamics have been extensively studied in ex vivo models; however, the recent development of two-photon microscopy and astrocyte-specific labeling has allowed the study of Ca²⁺ signaling in living central nervous system. Ca²⁺ waves in astrocytes have been described in cultured cells and slice preparations, but evidence for astrocytic activation during sensory activity is lacking. There are currently few methods to image living spinal cord: breathing and heart-beating artifacts have impeded the widespread application of this technique. We here imaged the living spinal cord by two-photon microscopy in C57BL6/J mice. Through pressurized injection, we specifically loaded spinal astrocytes using the red fluorescent dye sulforhodamine 101 (SR101) and imaged astrocytic Ca²⁺ levels with Oregon-Green BAPTA-1 (OGB). Then, we studied astrocytic Ca²⁺ levels at rest and after right electrical hind paw stimulation. Sensory stimulation significantly increased astrocytic Ca²⁺ levels within the superficial dorsal horn of the spinal cord compared to rest. In conclusion, in vivo morphofunctional imaging of living astrocytes in spinal cord revealed that astrocytes actively participate to sensory stimulation.

Quantitative fMRI and oxidative neuroenergetics

The discovery of functional magnetic resonance imaging (fMRI) has greatly impacted neuroscience. The blood oxygenation level-dependent (BOLD) signal, using deoxyhemoglobin as an endogenous paramagnetic contrast agent, exposes regions of interest in task-based and resting-state paradigms. However the BOLD contrast is at best a partial measure of neuronal activity, because the functional maps obtained by differencing or correlations ignore the total neuronal activity in the baseline state. Here we describe how studies of brain energy metabolism at Yale, especially with (13)C magnetic resonance spectroscopy and related techniques, contributed to development of quantitative functional brain imaging with fMRI by providing a reliable measurement of baseline energy. This narrative takes us on a journey, from molecules to mind, with illuminating insights about neuronal-glial activities in relation to energy demand of synaptic activity. These results, along with key contributions from laboratories worldwide, comprise the energetic basis for quantitative interpretation of fMRI data.

Prediction of Task-Related BOLD fMRI with Amplitude Signatures of Resting-State fMRI

Blood oxygen contrast-functional magnetic resonance imaging (fMRI) signals are a convolution of neural and vascular components. Several studies indicate that task-related (T-fMRI) or resting-state (R-fMRI) responses linearly relate to hypercapnic task responses. Based on the linearity of R-fMRI and T-fMRI with hypercapnia demonstrated by different groups using different study designs, we hypothesized that R-fMRI and T-fMRI signals are governed by a common physiological mechanism and that resting-state fluctuation of amplitude (RSFA) should be linearly related to T-fMRI responses. We tested this prediction in a group of healthy younger humans where R-fMRI, T-fMRI, and hypercapnic (breath hold, BH) task measures were obtained form the same scan session during resting state and during performance of motor and BH tasks. Within individual subjects, significant linear correlations were observed between motor and BH task responses across voxels. When averaged over the whole brain, the subject-wise correlation between the motor and BH tasks showed a similar linear relationship within the group. Likewise, a significant linear correlation was observed between motor-task activity and RSFA across voxels and subjects. The linear rest-task (R-T) relationship between motor activity and RSFA suggested that R-fMRI and T-fMRI responses are governed by similar physiological mechanisms. A practical use of the R-T relationship is its potential to estimate T-fMRI responses in special populations unable to perform tasks during fMRI scanning. Using the R-T relationship determined from the first group of 12 healthy subjects, we predicted the T-fMRI responses in a second group of 7 healthy subjects. RSFA in both the lower and higher frequency ranges robustly predicted the magnitude of T-fMRI responses at the subject and voxel levels. We propose that T-fMRI responses are reliably predictable to the voxel level in situations where only R-fMRI measures are possible, and may be useful for assessing neural activity in task non-compliant clinical populations.

S phase entry of neural progenitor cells correlates with increased blood flow in the young subventricular zone

The postnatal subventricular zone (SVZ) contains proliferating neural progenitor cells in close proximity to blood vessels. Insults and drug treatments acutely stimulate cell proliferation in the SVZ, which was assessed by labeling cells entering S phase. Although G1-to-S progression is metabolically demanding on a minute-to-hour time scale, it remains unknown whether increased SVZ cell proliferation is accompanied by a local hemodynamic response. This neurovascular coupling provides energy substrates to active neuronal assemblies. Transcardial dye perfusion revealed the presence of capillaries throughout the SVZ that constrict upon applications of the thromboxane A₂ receptor agonist U-46119 in acute brain slice preparations. We then monitored in vivo blood flow using laser Doppler flowmetry via a microprobe located either in the SVZ or a mature network. U-46119 injections into the lateral ventricle decreased blood flow in the SVZ and the striatum, which are near the ventricle. A 1-hour ventricular injection of epidermal and basic fibroblast growth factor (EGF and bFGF) significantly increased the percentage of Sox2 transcription factor-positive cells in S phase 1.5 hours post-injection. This increase was accompanied by a sustained rise in blood flow in the SVZ but not in the striatum. Direct growth factor injections into the cortex did not alter local blood flow, ruling out direct effects on capillaries. These findings suggest that an acute increase in the number of G1-to-S cycling SVZ cells is accompanied by neurometabolic-vascular coupling, which may provide energy and nutrient for cell cycle progression.

Analysis of time and space invariance of BOLD responses in the rat visual system

Neuroimaging studies of functional magnetic resonance imaging (fMRI) and electrophysiology provide the linkage between neural activity and the blood oxygenation level-dependent (BOLD) response. Here, BOLD responses to light flashes were imaged at 11.7T and compared with neural recordings from superior colliculus (SC) and primary visual cortex (V1) in rat brain--regions with different basal blood flow and energy demand. Our goal was to assess neurovascular coupling in V1 and SC as reflected by temporal/spatial variances of impulse response functions (IRFs) and assess, if any, implications for general linear modeling (GLM) of BOLD responses. Light flashes induced high magnitude neural/BOLD responses reproducibly from both regions. However, neural/BOLD responses from SC and V1 were markedly different. SC signals followed the boxcar shape of the stimulation paradigm at all flash rates, whereas V1 signals were characterized by onset/offset transients that exhibited different flash rate dependencies. We find that IRF(SC) is generally time-invariant across wider flash rate range compared with IRF(V1), whereas IRF(SC) and IRF(V1) are both space invariant. These results illustrate the importance of measured neural signals for interpretation of fMRI by showing that GLM of BOLD responses may lead to misinterpretation of neural activity in some cases.

Where fMRI and electrophysiology agree to disagree: corticothalamic and striatal activity patterns in the WAG/Rij rat

The relationship between neuronal activity and hemodynamic changes plays a central role in functional neuroimaging. Under normal conditions and in neurological disorders such as epilepsy, it is commonly assumed that increased functional magnetic resonance imaging (fMRI) signals reflect increased neuronal activity and that fMRI decreases represent neuronal activity decreases. Recent work suggests that these assumptions usually hold true in the cerebral cortex. However, less is known about the basis of fMRI signals from subcortical structures such as the thalamus and basal ganglia. We used WAG/Rij rats (Wistar albino Glaxo rats of Rijswijk), an established animal model of human absence epilepsy, to perform fMRI studies with blood oxygen level-dependent and cerebral blood volume (CBV) contrasts at 9.4 tesla, as well as laser Doppler cerebral blood flow (CBF), local field potential (LFP), and multiunit activity (MUA) recordings. We found that, during spike-wave discharges, the somatosensory cortex and thalamus showed increased fMRI, CBV, CBF, LFP, and MUA signals. However, the caudate-putamen showed fMRI, CBV, and CBF decreases despite increases in LFP and MUA signals. Similarly, during normal whisker stimulation, the cortex and thalamus showed increases in CBF and MUA, whereas the caudate-putamen showed decreased CBF with increased MUA. These findings suggest that neuroimaging-related signals and electrophysiology tend to agree in the cortex and thalamus but disagree in the caudate-putamen. These opposite changes in vascular and electrical activity indicate that caution should be applied when interpreting fMRI signals in both health and disease from the caudate-putamen, as well as possibly from other subcortical structures.

Fractal analysis of spontaneous fluctuations of the BOLD signal in rat brain

Analysis of task-evoked fMRI data ignores low frequency fluctuations (LFF) of the resting-state the BOLD signal, yet LFF of the spontaneous BOLD signal is crucial for analysis of resting-state connectivity maps. We characterized the LFF of resting-state BOLD signal at 11.7T in α-chloralose and domitor anesthetized rat brain and modeled the spontaneous signal as a scale-free (i.e., fractal) distribution of amplitude power (|A|²) across a frequency range (f) compatible with an |A(f)|² ∝ 1/f(β) model where β is the scaling exponent (or spectral index). We compared β values from somatosensory forelimb area (S1FL), cingulate cortex (CG), and caudate putamen (CPu). With α-chloralose, S1FL and CG β values dropped from ~0.7 at in vivo to ~0.1 at post mortem (p<0.0002), whereas CPu β values dropped from ~0.3 at in vivo to ~0.1 at post mortem (p<0.002). With domitor, cortical (S1FL, CG) β values were slightly higher than with α-chloralose, while subcortical (CPu) β values were similar with α-chloralose. Although cortical and subcortical β values with both anesthetics were significantly different in vivo (p<0.002), at post mortem β values in these regions were not significantly different and approached zero (i.e., range of -0.1 to 0.2). Since a water phantom devoid of susceptibility gradients had a β value of zero (i.e., random), we conclude that deoxyhemoglobin present in voxels post-sacrifice still impacts tissue water diffusion. These results suggest that in the anesthetized rat brain the LFF of BOLD signal at 11.7T follow a general 1/f(β) model of fractality where β is a variable responding to physiology. We describe typical experimental pitfalls which may elude detection of fractality in the resting-state BOLD signal.

Functional micro-ultrasound imaging of rodent cerebral hemodynamics

Healthy cerebral microcirculation is crucial to neuronal functioning. We present a new method to investigate microvascular hemodynamics in living rodent brain through a focal cranial window based on high-frequency ultrasound imaging. The method has a temporal resolution of 40ms, and a 100μm in-plane and 600μm through-plane spatial resolution. We use a commercially available high-frequency ultrasound imaging system to quantify changes in the relative cerebral blood volume (CBV) by measuring the scattered signal intensity from an ultrasound contrast agent circulating in the vasculature. Generalized linear model analysis is then used to produce effect size and significance maps of changes in cerebral blood volume upon electrical stimulation of the forepaw. We observe larger CBV increases in the forelimb representation of the primary somatosensory cortex than in the deep gray matter with stimuli as short as 2s (5.1 ± 1.3% vs. 3.3 ± 0.6%). We also investigate the temporal evolution of the blood volume changes in cortical and subcortical gray matter, pial vessels and subcortical major vessels, and show shorter response onset times in the parenchymal regions than in the neighboring large vessels (1.6 ± 1.0s vs. 2.6 ± 1.3s in the cortex for a 10 second stimulus protocol). This method, which we termed functional micro-ultrasound imaging or fMUS, is a novel, highly accessible, and cost-effective way of imaging rodent brain microvascular topology and hemodynamics in vivo at 100micron resolution over a 1-by-1cm field of view with 10s-100s frames per second that opens up a new set of questions regarding brain function in preclinical models of health and disease.

Neural and BOLD responses across the brain

Functional Magnetic Resonance Imaging (fMRI) has quickly grown into one of the most important tools for studying brain function, especially in humans. Despite its prevalence, we still do not have a clear picture of what exactly the blood oxygenation level dependent (BOLD) signal represents or how it compares to the signals obtained with other methods (e.g., electrophysiology). We particularly refer to single neuron recordings and electroencephalography when we mention 'electrophysiological methods', given that these methods have been used for more than 50 years, and have formed the basis of much of our current understanding of brain function. Brain function involves the coordinated activity of many different areas and many different cell types that can participate in an enormous variety of processes (neural firing, inhibitory and excitatory synaptic activity, neuromodulation, oscillatory activity, etc.). Of these cells and processes, only a subset is sampled with electrophysiological techniques, and their contribution to the recorded signals is not exactly known. Functional imaging signals are driven by the metabolic needs of the active cells, and are most likely also biased toward certain cell types and certain neural processes, although we know even less about which processes actually drive the hemodynamic response. This article discusses the current status on the interpretation of the BOLD signal and how it relates to neural activity measured with electrophysiological techniques. WIREs Cogn Sci 2012, 3:75-86. doi: 10.1002/wcs.153 For further resources related to this article, please visit the WIREs website.

Spatiotemporal evolution of the functional magnetic resonance imaging response to ultrashort stimuli

The specificity of the hemodynamic response function (HRF) is determined spatially by the vascular architecture and temporally by the evolution of hemodynamic changes. The stimulus duration has additional influence on the spatiotemporal evolution of the HRF, as brief stimuli elicit responses that engage only the local vasculature, whereas long stimuli lead to the involvement of remote vascular supply and drainage. Here, we used functional magnetic resonance imaging to investigate the spatiotemporal evolution of the blood oxygenation level-dependent (BOLD), cerebral blood flow (CBF), and cerebral blood volume (CBV) HRF to ultrashort forelimb stimulation in an anesthetized rodent model. The HRFs to a single 333-μs-long stimulus were robustly detected and consisted of a rapid response in both CBF and CBV, with an onset time (OT) of 350 ms and a full width at half-maximum of 1 s. In contrast, longer stimuli elicited a dispersive transit of oxygenated blood across the cortical microvasculature that significantly prolonged the evolution of the CBV HRF, but not the CBF. The CBF and CBV OTs suggest that vasoactive messengers are synthesized, released, and effective within 350 ms. However, the difference between the BOLD and CBV OT (∼100 ms) was significantly smaller than the arteriolar-venular transit time (∼500 ms), indicating an arterial contribution to the BOLD HRF. Finally, the rapid rate of growth of the active region with stimulus elongation suggests that functional hyperemia is an integrative process that involves the entire functional cortical depth. These findings offer a new view into the spatiotemporal dynamics of functional hemodynamic regulation in the brain.

Neurovascular and Neurometabolic Couplings in Dynamic Calibrated fMRI: Transient Oxidative Neuroenergetics for Block-Design and Event-Related Paradigms

Functional magnetic resonance imaging (fMRI) with blood-oxygenation level dependent (BOLD) contrast is an important tool for mapping brain activity. Interest in quantitative fMRI has renewed awareness in importance of oxidative neuroenergetics, as reflected by cerebral metabolic rate of oxygen consumption(CMRO2), for supporting brain function. Relationships between BOLD signal and the underlying neurophysiological parameters have been elucidated to allow determination of dynamic changes inCMRO2 by "calibrated fMRI," which require multi-modal measurements of BOLD signal along with cerebral blood flow (CBF) and volume (CBV). But how doCMRO2 changes, steady-state or transient, derived from calibrated fMRI compare with neural activity recordings of local field potential (LFP) and/or multi-unit activity (MUA)? Here we discuss recent findings primarily from animal studies which allow high magnetic fields studies for superior BOLD sensitivity as well as multi-modal CBV and CBF measurements in conjunction with LFP and MUA recordings from activated sites. A key observation is that while relationships between neural activity and sensory stimulus features range from linear to non-linear, associations between hyperemic components (BOLD, CBF, CBV) and neural activity (LFP, MUA) are almost always linear. More importantly, the results demonstrate good agreement between the changes inCMRO2 and independent measures of LFP or MUA. The tight neurovascular and neurometabolic couplings, observed from steady-state conditions to events separated by <200 ms, suggest rapid oxygen equilibration between blood and tissue pools and thus calibrated fMRI at high magnetic fields can provide high spatiotemporal mapping ofCMRO2 changes.

No increase of the blood oxygenation level-dependent functional magnetic resonance imaging signal with higher field strength: implications for brain activation studies

Experimental data up to 7.0 T show that the blood oxygenation level-dependent (BOLD) signal of functional magnetic resonance imaging (fMRI) increases with higher magnetic field strength. Although several studies at 11.7 T report higher BOLD signal compared with studies at 7.0 T, no direct comparison at these two field strengths has been performed under the exact same conditions. It therefore remains unclear whether the expected increase of BOLD effect with field strength will still continue to hold for fields >7.0 T. To examine this issue, we compared the BOLD activation signal at 7.0 and 11.7 T with the two common sequences, spin-echo (SE) and gradient-echo (GE) echo planar imaging (EPI). We chose the physiologically well controlled rat model of electrical forepaw stimulation under medetomidine sedation. While a linear to superlinear increase in activation with field strengths up to 7.0 T was reported in the literature, we observed no significant activation difference between 7.0 and 11.7 T with either SE or GE. Discussing the results in light of the four-component model of the BOLD signal, we showed that at high field only two extravascular contributions remain relevant, while both intravascular components vanish. Constancy of the BOLD effect is discussed due to motional narrowing, i.e., susceptibility gradients become so strong that phase variance of diffusing spins decreases and therefore the BOLD signal also decreases. This finding will be of high significance for the planning of future human and animal fMRI studies at high fields and their quantitative analysis.

Quantitative functional magnetic resonance imaging of brain activity using bolus-tracking arterial spin labeling

Blood oxygen level dependent (BOLD) functional magnetic resonance imaging (fMRI) is the most widely used method for mapping neural activity in the brain. The interpretation of altered BOLD signals is problematic when cerebral blood flow (CBF) or cerebral blood volume change because of aging and/or neurodegenerative diseases. In this study, a recently developed quantitative arterial spin labeling (ASL) approach, bolus-tracking ASL (btASL), was applied to an fMRI experiment in the rat brain. The mean transit time (MTT), capillary transit time (CTT), relative cerebral blood volume of labeled water (rCBV(lw)), relative cerebral blood flow (rCBF), and perfusion coefficient in the forelimb region of the somatosensory cortex were quantified during neuronal activation and in the resting state. The average MTT and CTT were 1.939+/-0.175 and 1.606+/-0.106 secs, respectively, in the resting state. Both times decreased significantly to 1.616+/-0.207 and 1.305+/-0.201 secs, respectively, during activation. The rCBV(lw), rCBF, and perfusion coefficient increased on average by a factor of 1.123+/-0.006, 1.353+/-0.078, and 1.479+/-0.148, respectively, during activation. In contrast to BOLD techniques, btASL yields physiologically relevant indices of the functional hyperemia that accompanies neuronal activation.