Cerebral oxygen demand for short-lived and steady-state events

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.

Sparsely Wiring Connectivity in the Upper Beta Band Characterizes the Brains of Top Swimming Athletes

Human brains are extremely energy costly in neural connections and activities. However, it is unknown what is the difference in the brain connectivity between top athletes with long-term professional trainings and age-matched controls. Here we ask whether long-term training can lower brain-wiring cost while have better performance. Since elite swimming requires athletes to move their arms and legs at different tempos in time with high coordination skills, we selected an eye-hand-foot complex reaction (CR) task to examine the relations between the task performance and the brain connections and activities, as well as to explore the energy cost-efficiency of top athletes. Twenty-one master-level professional swimmers and 23 age-matched non-professional swimmers as controls were recruited to perform the CR task with concurrent 8-channel EEG recordings. Reaction time and accuracy of the CR task were recorded. Topological network analysis of various frequency bands was performed using the phase lag index (PLI) technique to avoid volume conduction effects. The wiring number of connections and mean frequency were calculated to reflect the wiring and activity cost, respectively. Results showed that professional athletes demonstrated better eye-hand-foot coordination than controls when performing the CR task, indexing by faster reaction time and higher accuracy. Comparing to controls, athletes' brain demonstrated significantly less connections and weaker correlations in upper beta frequency band between the frontal and parietal regions, while demonstrated stronger connectivity in the low theta frequency band between sites of F3 and Cz/C4. Additionally, athletes showed highly stable and low eye-blinking rates across different reaction performance, while controls had high blinking frequency with high variance. Elite athletes' brain may be characterized with energy efficient sparsely wiring connections in support of superior motor performance and better cognitive performance in the eye-hand-foot complex reaction task.

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.

Resting cerebral oxygen metabolism exhibits archetypal network features

Standard magnetic resonance imaging approaches offer high‐resolution but indirect measures of neural activity, limiting understanding of the physiological processes associated with imaging findings. Here, we used calibrated functional magnetic resonance imaging during the resting state to recover low‐frequency fluctuations of the cerebral metabolic rate of oxygen (CMRO₂). We tested whether functional connections derived from these fluctuations exhibited organization properties similar to those established by previous standard functional and anatomical connectivity studies. Seventeen participants underwent 20 min of resting imaging during dual‐echo, pseudocontinuous arterial spin labeling, and blood‐oxygen‐level dependent (BOLD) signal acquisition. Participants also underwent a 10 min normocapnic and hypercapnic procedure. Brain‐wide, CMRO₂ low‐frequency fluctuations were subjected to graph‐based and voxel‐wise functional connectivity analyses. Results demonstrated that connections derived from resting CMRO₂ fluctuations exhibited complex, small‐world topological properties (i.e., high integration and segregation, cost efficiency) consistent with those observed in previous studies using functional and anatomical connectivity approaches. Voxel‐wise CMRO₂ connectivity also exhibited spatial patterns consistent with four targeted resting‐state subnetworks: two association (i.e., frontoparietal and default mode) and two perceptual (i.e., auditory and occipital‐visual). These are the first findings to support the use of calibration‐derived CMRO₂ low‐frequency fluctuations for detecting brain‐wide organizational properties typical of healthy participants. We discuss interpretations, advantages, and challenges in using calibration‐derived oxygen metabolism signals for examining the intrinsic organization of the human brain.

A cortical rat hemodynamic response function for improved detection of BOLD activation under common experimental conditions

For a reliable estimation of neuronal activation based on BOLD fMRI measurements an accurate model of the hemodynamic response is essential. Since a large part of basic neuroscience research is based on small animal data, it is necessary to characterize a hemodynamic response function (HRF) which is optimized for small animals. Therefore, we have determined and investigated the HRFs of rats obtained under a variety of experimental conditions in the primary somatosensory cortex. Measurements were performed on animals of different sex and strain, under different anesthetics, with and without ventilation and using different stimulation modalities. All modalities of stimulation used in this study induced neuronal activity in the primary somatosensory cortex or in subcortical regions. Since the HRFs of the BOLD responses in the primary somatosensory cortex showed a close concordance for the different conditions, we were able to determine a cortical rat HRF. This HRF is based on 143 BOLD measurements of 76 rats and can be used for statistical parametric mapping. It showed substantially faster progression than the human HRF, with a maximum after 2.8 ± 0.8 s, and a following undershoot after 6.1 ± 3.7 s. If the rat HRF was used statistical analysis of rat data showed a significantly improved detection performance in the somatosensory cortex in comparison to the commonly used HRF based on measurements in humans.

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.

Calibrated imaging reveals altered grey matter metabolism related to white matter microstructure and symptom severity in multiple sclerosis

Multiple sclerosis (MS) involves damage to white matter microstructures. This damage has been related to grey matter function as measured by standard, physiologically-nonspecific neuroimaging indices (i.e., blood-oxygen-level dependent signal [BOLD]). Here, we used calibrated functional magnetic resonance imaging and diffusion tensor imaging to examine the extent to which specific, evoked grey matter physiological processes were associated with white matter diffusion in MS. Evoked changes in BOLD, cerebral blood flow (CBF), and oxygen metabolism (CMRO₂ ) were measured in visual cortex. Individual differences in the diffusion tensor measure, radial diffusivity, within occipital tracts were strongly associated with MS patients' BOLD and CMRO₂ . However, these relationships were in opposite directions, complicating the interpretation of the relationship between BOLD and white matter microstructural damage in MS. CMRO₂ was strongly associated with individual differences in patients' fatigue and neurological disability, suggesting that alterations to evoked oxygen metabolic processes may be taken as a marker for primary symptoms of MS. This work demonstrates the first application of calibrated and diffusion imaging together and details the first application of calibrated functional MRI in a neurological population. Results lend support for neuroenergetic hypotheses of MS pathophysiology and provide an initial demonstration of the utility of evoked oxygen metabolism signals for neurology research. Hum Brain Mapp 38:5375-5390, 2017. © 2017 Wiley Periodicals, Inc.

Evaluation of Visual-Evoked Cerebral Metabolic Rate of Oxygen as a Diagnostic Marker in Multiple Sclerosis

A multiple sclerosis (MS) diagnosis often relies upon clinical presentation and qualitative analysis of standard, magnetic resonance brain images. However, the accuracy of MS diagnoses can be improved by utilizing advanced brain imaging methods. We assessed the accuracy of a new neuroimaging marker, visual-evoked cerebral metabolic rate of oxygen (veCMRO₂), in classifying MS patients and closely age- and sex-matched healthy control (HC) participants. MS patients and HCs underwent calibrated functional magnetic resonance imaging (cfMRI) during a visual stimulation task, diffusion tensor imaging, T₁- and T₂-weighted imaging, neuropsychological testing, and completed self-report questionnaires. Using resampling techniques to avoid bias and increase the generalizability of the results, we assessed the accuracy of veCMRO₂ in classifying MS patients and HCs. veCMRO₂ classification accuracy was also examined in the context of other evoked visuofunctional measures, white matter microstructural integrity, lesion-based measures from T₂-weighted imaging, atrophy measures from T₁-weighted imaging, neuropsychological tests, and self-report assays of clinical symptomology. veCMRO₂ was significant and within the top 16% of measures (43 total) in classifying MS status using both within-sample (82% accuracy) and out-of-sample (77% accuracy) observations. High accuracy of veCMRO₂ in classifying MS demonstrated an encouraging first step toward establishing veCMRO₂ as a neurodiagnostic marker of MS.

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.

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.

Understanding the dynamic relationship between cerebral blood flow and the BOLD signal: Implications for quantitative functional MRI

Calibrated BOLD imaging, in which traditional measurements of the BOLD signal are combined with measurements of cerebral blood flow (CBF) within a BOLD biophysical model to estimate changes in oxygen metabolism (CMRO2), has been a valuable tool for untangling the physiological processes associated with neural stimulus-induced BOLD activation. However, to date this technique has largely been applied to the study of essentially steady-state physiological changes (baseline to activation) associated with block-design stimuli, and it is unclear whether this approach may be directly extended to the study of more dynamic, naturalistic experimental designs. In this study we tested an assumption underlying this technique whose validity is critical to the application of calibrated BOLD to the study of more dynamic stimuli, that information about fluctuations in venous cerebral blood volume (CBVv) can be captured indirectly by measuring fluctuations in CBF, making the independent measurement of CBVv unnecessary. To accomplish this, simultaneous arterial spin labeling and BOLD imaging were used to measure the CBF and BOLD responses to flickering checkerboards with contrasts that oscillated continuously with frequencies of ~0.02-0.16Hz. The measurements were then fit to a dynamic physiological model of the BOLD response in order to explore the range of consistent CMRO2 and CBVv responses. We found that the BOLD and CBF responses were most consistent with relatively tight dynamic coupling between CBF and CMRO2 and a CBVv response that was an order of magnitude slower than either CBF or CMRO2. This finding suggests that the assumption of tight flow-volume coupling may not be strictly valid, complicating the extension of calibrated BOLD to more naturalistic experimental designs.

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.

Are glutamate and lactate increases ubiquitous to physiological activation? A (1)H functional MR spectroscopy study during motor activation in human brain at 7Tesla

Recent studies at high field (7Tesla) have reported small metabolite changes, in particular lactate and glutamate (below 0.3μmol/g) during visual stimulation. These studies have been limited to the visual cortex because of its high energy metabolism and good magnetic resonance spectroscopy (MRS) sensitivity using surface coil. The aim of this study was to extend functional MRS (fMRS) to investigate for the first time the metabolite changes during motor activation at 7T. Small but sustained increases in lactate (0.17μmol/g±0.05μmol/g, p<0.001) and glutamate (0.17μmol/g±0.09μmol/g, p<0.005) were detected during motor activation followed by a return to the baseline after the end of activation. The present study demonstrates that increases in lactate and glutamate during motor stimulation are small, but similar to those observed during visual stimulation. From the observed glutamate and lactate increase, we inferred that these metabolite changes may be a general manifestation of the increased neuronal activity. In addition, we propose that the measured metabolite concentration increases imply an increase in ΔCMRO2 that is transiently below that of ΔCMRGlc during the first 1 to 2min of the stimulation.

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.

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.

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.

Evolution of the dynamic changes in functional cerebral oxidative metabolism from tissue mitochondria to blood oxygen

The dynamic properties of the cerebral metabolic rate of oxygen consumption (CMR(O2)) during changes in brain activity remain unclear. Therefore, the spatial and temporal evolution of functional increases in CMR(O2) was investigated in the rat somato-sensory cortex during forelimb stimulation under a suppressed blood flow response condition. Temporally, stimulation elicited a fast increase in tissue mitochondria CMR(O2) described by a time constant of ~1 second measured using flavoprotein autofluorescence imaging. CMR(O2)-driven changes in the tissue oxygen tension measured using an oxygen electrode and blood oxygenation measured using optical imaging of intrinsic signal followed; however, these changes were slow with time constants of ~5 and ~10 seconds, respectively. This slow change in CMR(O2)-driven blood oxygenation partly explains the commonly observed post-stimulus blood oxygen level-dependent (BOLD) undershoot. Spatially, the changes in mitochondria CMR(O2) were similar to the changes in blood oxygenation. Finally, the increases in CMR(O2) were well correlated with the evoked multi-unit spiking activity. These findings show that dynamic CMR(O2) calculations made using only blood oxygenation data (e.g., BOLD functional magnetic resonance imaging (fMRI)) do not directly reflect the temporal changes in the tissue's mitochondria metabolic rate; however, the findings presented can bridge the gap between the changes in cellular oxidative rate and blood oxygenation.

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.

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.

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.

Interpreting oxygenation-based neuroimaging signals: the importance and the challenge of understanding brain oxygen metabolism

Functional magnetic resonance imaging is widely used to map patterns of brain activation based on blood oxygenation level dependent (BOLD) signal changes associated with changes in neural activity. However, because oxygenation changes depend on the relative changes in cerebral blood flow (CBF) and cerebral metabolic rate of oxygen (CMRO(2)), a quantitative interpretation of BOLD signals, and also other functional neuroimaging signals related to blood or tissue oxygenation, is fundamentally limited until we better understand brain oxygen metabolism and how it is related to blood flow. However, the positive side of the complexity of oxygenation signals is that when combined with dynamic CBF measurements they potentially provide the best tool currently available for investigating the dynamics of CMRO(2). This review focuses on the problem of interpreting oxygenation-based signals, the challenges involved in measuring CMRO(2) in general, and what is needed to put oxygenation-based estimates of CMRO(2) on a firm foundation. The importance of developing a solid theoretical framework is emphasized, both as an essential tool for analyzing oxygenation-based multimodal measurements, and also potentially as a way to better understand the physiological phenomena themselves. The existing data, integrated within a simple theoretical framework of O(2) transport, suggests the hypothesis that an important functional role of the mismatch of CBF and CMRO(2) changes with neural activation is to prevent a fall of tissue pO(2). Future directions for better understanding brain oxygen metabolism are discussed.

Linear and nonlinear relationships between visual stimuli, EEG and BOLD fMRI signals

In the present study, the cascaded interactions between stimuli and neural and hemodynamic responses were modeled using linear systems. These models provided the theoretical hypotheses that were tested against the electroencephalography (EEG) and blood oxygen level dependent (BOLD) functional magnetic resonance imaging (fMRI) data recorded from human subjects during prolonged periods of repeated visual stimuli with a variable setting of the inter-stimulus interval (ISI) and visual contrast. Our results suggest that (1) neural response is nonlinear only when ISI<0.2 s, (2) BOLD response is nonlinear with an exclusively vascular origin when 0.25<ISI<4.2 s, (3) vascular response nonlinearity reflects a refractory effect, rather than a ceiling effect, and (4) there is a strong linear relationship between the BOLD effect size and the integrated power of event-related synaptic current activity, after modeling and taking into account the vascular refractory effect. These conclusions offer important insights into the origins of BOLD nonlinearity and the nature of neurovascular coupling, and suggest an effective means to quantitatively interpret the BOLD signal in terms of neural activity. The validated cross-modal relationship between fMRI and EEG may provide a theoretical basis for the integration of these two modalities.

Oxidative neuroenergetics in event-related paradigms

Energetic basis of neural activity provides a solid foundation for noninvasive neuroimaging with calibrated functional magnetic resonance imaging (fMRI). Calculating dynamic changes in cerebral oxidative energy utilization (CMR(O(2))) is limited by uncertainties about whether or not the conventional blood oxygenation level-dependent (BOLD) model can be applied transiently using multimodal measurements of blood flow (CBF) and volume (CBV) that affect the BOLD signal. A prerequisite for dynamic calibrated fMRI is testing the linearity of multimodal signals within a temporal regimen, as assessed by signal strength (i.e., both intensity and width). If each hyperemic component (BOLD, CBV, CBF) is demonstrated to be linear with neural activity under various experimental conditions, then the respective transfer functions generated by deconvolution with neural activity should be time invariant and thus could potentially be used for calculating CMR(O(2)) transients. Hyperemic components were investigated at 11.7 T in alpha-chloralose-anesthetized rats and combined with electrophysiological recordings of local field potential (LFP) and multiunit activity (MUA) from the cortex during forepaw stimulation, in which stimulus number and frequency were varied. Although relationships between neural activity and stimulus features ranged from linear to nonlinear, associations between hyperemic components and neural activity were linear. Specific to each hyperemic component, a universal transfer function (with LFP or MUA) yielded predictions in agreement with experimental measurements. The results identified a component of the BOLD signal that can be attributed to significant changes in CMR(O(2)), even for temporal events separated by <200 ms.