Assessment of cerebral oxidative metabolism with breath holding and fMRI

Functional magnetic resonance imaging signal has sub-second temporal accuracy

Neuronal activation sequence information is essential for understanding brain functions. Extracting such timing information from blood-oxygenation-level-dependent functional magnetic resonance imaging (fMRI) signals is confounded by local cerebral vascular reactivity (CVR), which varies across brain locations. Thus, detecting neuronal synchrony as well as inferring inter-regional causal modulation using fMRI signals can be biased. Here we used fast fMRI measurements sampled at 10 Hz to measure the fMRI latency difference between visual and sensorimotor areas when participants engaged in a visuomotor task. The regional fMRI timing was calibrated by subtracting the CVR latency measured by a breath-holding task. After CVR calibration, the fMRI signal at the lateral geniculate nucleus (LGN) preceded that at the visual cortex by 496 ms, followed by the fMRI signal at the sensorimotor cortex with a latency of 464 ms. Sequential LGN, visual, and sensorimotor cortex activations were found in each participant after the CVR calibration. These inter-regional fMRI timing differences across and within participants were more closely related to the reaction time after the CVR calibration. Our results suggested the feasibility of mapping brain activity using fMRI with accuracy in hundreds of milliseconds.

Evaluation of calibrated and uncalibrated optical imaging approaches for relative cerebral oxygen metabolism measurements in awake mice

Objective. The continuous delivery of oxygen is critical to sustain brain function, and therefore, measuring brain oxygen consumption can provide vital physiological insight. In this work, we examine the impact of calibration and cerebral blood flow (CBF) measurements on the computation of the relative changes in the cerebral metabolic rate of oxygen consumption (rCMRO2) from hemoglobin-sensitive intrinsic optical imaging data. Using these data, we calculate rCMRO2, and calibrate the model using an isometabolic stimulus.Approach. We used awake head-fixed rodents to obtain hemoglobin-sensitive optical imaging data to test different calibrated and uncalibrated rCMRO2models. Hypercapnia was used for calibration and whisker stimulation was used to test the impact of calibration.Main results. We found that typical uncalibrated models can provide reasonable estimates of rCMRO2with differences as small as 7%-9% compared to their calibrated models. However, calibrated models showed lower variability and less dependence on baseline hemoglobin concentrations. Lastly, we found that supplying the model with measurements of CBF significantly reduced error and variability in rCMRO2change calculations.Significance. The effect of calibration on rCMRO2calculations remains understudied, and we systematically evaluated different rCMRO2calculation scenarios that consider including different measurement combinations. This study provides a quantitative comparison of these scenarios to evaluate trade-offs that can be vital to the design of blood oxygenation sensitive imaging experiments for rCMRO2calculation.

Resting-state fMRI signals contain spectral signatures of local hemodynamic response timing

Functional magnetic resonance imaging (fMRI) has proven to be a powerful tool for noninvasively measuring human brain activity; yet, thus far, fMRI has been relatively limited in its temporal resolution. A key challenge is understanding the relationship between neural activity and the blood-oxygenation-level-dependent (BOLD) signal obtained from fMRI, generally modeled by the hemodynamic response function (HRF). The timing of the HRF varies across the brain and individuals, confounding our ability to make inferences about the timing of the underlying neural processes. Here, we show that resting-state fMRI signals contain information about HRF temporal dynamics that can be leveraged to understand and characterize variations in HRF timing across both cortical and subcortical regions. We found that the frequency spectrum of resting-state fMRI signals significantly differs between voxels with fast versus slow HRFs in human visual cortex. These spectral differences extended to subcortex as well, revealing significantly faster hemodynamic timing in the lateral geniculate nucleus of the thalamus. Ultimately, our results demonstrate that the temporal properties of the HRF impact the spectral content of resting-state fMRI signals and enable voxel-wise characterization of relative hemodynamic response timing. Furthermore, our results show that caution should be used in studies of resting-state fMRI spectral properties, because differences in fMRI frequency content can arise from purely vascular origins. This finding provides new insight into the temporal properties of fMRI signals across voxels, which is crucial for accurate fMRI analyses, and enhances the ability of fast fMRI to identify and track fast neural dynamics.

Resting-state fMRI signals contain spectral signatures of local hemodynamic response timing

Functional magnetic resonance imaging (fMRI) has proven to be a powerful tool for noninvasively measuring human brain activity; yet, thus far, fMRI has been relatively limited in its temporal resolution. A key challenge is understanding the relationship between neural activity and the blood-oxygenation-level-dependent (BOLD) signal obtained from fMRI, generally modeled by the hemodynamic response function (HRF). The timing of the HRF varies across the brain and individuals, confounding our ability to make inferences about the timing of the underlying neural processes. Here we show that resting-state fMRI signals contain information about HRF temporal dynamics that can be leveraged to understand and characterize variations in HRF timing across both cortical and subcortical regions. We found that the frequency spectrum of resting-state fMRI signals significantly differs between voxels with fast versus slow HRFs in human visual cortex. These spectral differences extended to subcortex as well, revealing significantly faster hemodynamic timing in the lateral geniculate nucleus of the thalamus. Ultimately, our results demonstrate that the temporal properties of the HRF impact the spectral content of resting-state fMRI signals and enable voxel-wise characterization of relative hemodynamic response timing. Furthermore, our results show that caution should be used in studies of resting-state fMRI spectral properties, as differences can arise from purely vascular origins. This finding provides new insight into the temporal properties of fMRI signals across voxels, which is crucial for accurate fMRI analyses, and enhances the ability of fast fMRI to identify and track fast neural dynamics.

fMRI studies evaluating central respiratory control in humans

A plethora of neural centers in the central nervous system control the fundamental respiratory pattern. This control is ensured by neurons that act as pacemakers, modulating activity through chemical control driven by changes in the O2/CO2 balance. Most of the respiratory neural centers are located in the brainstem, but difficult to localize on magnetic resonance imaging (MRI) due to their small size, lack of visually-detectable borders with neighboring areas, and significant physiological noise hampering detection of its activity with functional MRI (fMRI). Yet, several approaches make it possible to study the normal response to different abnormal stimuli or conditions such as CO2 inhalation, induced hypercapnia, volitional apnea, induced hypoxia etc. This review provides a comprehensive overview of the majority of available studies on central respiratory control in humans.

Real-Time Functional MRI in the Treatment of Mental Health Disorders

Multiple mental disorders have been associated with dysregulation of precise brain processes. However, few therapeutic approaches can correct such specific patterns of brain activity. Since the late 1960s and early 1970s, many researchers have hoped that this feat could be achieved by closed-loop brain imaging approaches, such as neurofeedback, that aim to modulate brain activity directly. However, neurofeedback never gained mainstream acceptance in mental health, in part due to methodological considerations. In this review, we argue that, when contemporary methodological guidelines are followed, neurofeedback is one of the few intervention methods in psychology that can be assessed in double-blind placebo-controlled trials. Furthermore, using new advances in machine learning and statistics, it is now possible to target very precise patterns of brain activity for therapeutic purposes. We review the recent literature in functional magnetic resonance imaging neurofeedback and discuss current and future applications to mental health.

Expected final online publication date for the Annual Review of Clinical Psychology, Volume 18 is May 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Separating vascular and neuronal effects of age on fMRI BOLD signals

Accurate identification of brain function is necessary to understand the neurobiology of cognitive ageing, and thereby promote well-being across the lifespan. A common tool used to investigate neurocognitive ageing is functional magnetic resonance imaging (fMRI). However, although fMRI data are often interpreted in terms of neuronal activity, the blood oxygenation level-dependent (BOLD) signal measured by fMRI includes contributions of both vascular and neuronal factors, which change differentially with age. While some studies investigate vascular ageing factors, the results of these studies are not well known within the field of neurocognitive ageing and therefore vascular confounds in neurocognitive fMRI studies are common. Despite over 10 000 BOLD-fMRI papers on ageing, fewer than 20 have applied techniques to correct for vascular effects. However, neurovascular ageing is not only a confound in fMRI, but an important feature in its own right, to be assessed alongside measures of neuronal ageing. We review current approaches to dissociate neuronal and vascular components of BOLD-fMRI of regional activity and functional connectivity. We highlight emerging evidence that vascular mechanisms in the brain do not simply control blood flow to support the metabolic needs of neurons, but form complex neurovascular interactions that influence neuronal function in health and disease. This article is part of the theme issue 'Key relationships between non-invasive functional neuroimaging and the underlying neuronal activity'.

Toward a comprehensive understanding of the neural mechanisms of decoded neurofeedback

Real-time functional magnetic resonance imaging (fMRI) neurofeedback is an experimental framework in which fMRI signals are presented to participants in a real-time manner to change their behaviors. Changes in behaviors after real-time fMRI neurofeedback are postulated to be caused by neural plasticity driven by the induction of specific targeted activities at the neuronal level (targeted neural plasticity model). However, some research groups argued that behavioral changes in conventional real-time fMRI neurofeedback studies are explained by alternative accounts, including the placebo effect and physiological artifacts. Recently, decoded neurofeedback (DecNef) has been developed as a result of adapting new technological advancements, including implicit neurofeedback and fMRI multivariate analyses. DecNef provides strong evidence for the targeted neural plasticity model while refuting the abovementioned alternative accounts. In this review, we first discuss how DecNef refutes the alternative accounts. Second, we propose a model that shows how targeted neural plasticity occurs at the neuronal level during DecNef training. Finally, we discuss computational and empirical evidence that supports the model. Clarification of the neural mechanisms of DecNef would lead to the development of more advanced fMRI neurofeedback methods that may serve as powerful tools for both basic and clinical research.

The Impact of Echo Time Shifts and Temporal Signal Fluctuations on BOLD Sensitivity in Presurgical Planning at 7 T

OBJECTIVES

Gradients in the static magnetic field caused by tissues with differing magnetic susceptibilities lead to regional variations in the effective echo time, which modifies both image signal and BOLD sensitivity. Local echo time changes are not considered in the most commonly used metric for BOLD sensitivity, temporal signal-to-noise ratio (tSNR), but may be significant, particularly at ultrahigh field close to air cavities (such as the sinuses and ear canals) and near gross brain pathologies and postoperative sites.

MATERIALS AND METHODS

We have studied the effect of local variations in echo time and tSNR on BOLD sensitivity in 3 healthy volunteers and 11 patients with tumors, postoperative cavities, and venous malformations at 7 T. Temporal signal-to-noise ratio was estimated from a 5-minute run of resting state echo planar imaging with a nominal echo time of 22 milliseconds. Maps of local echo time were derived from the phase of a multiecho GE scan. One healthy volunteer performed 10 runs of a breath-hold task. The t-map from this experiment served as a criterion standard BOLD sensitivity measure. Two runs of a less demanding breath-hold paradigm were used for patients.

RESULTS

In all subjects, a strong reduction in the echo time (from 22 milliseconds to around 11 milliseconds) was found close to the ear canals and sinuses. These regions were characterized by high tSNR but low t-values in breath-hold t-maps. In some patients, regions of particular interest in presurgical planning were affected by reductions in the echo time to approximately 13-15 milliseconds. These included the primary motor cortex, Broca's area, and auditory cortex. These regions were characterized by high tSNR values (70 and above). Breath-hold results were corrupted by strong motion artifacts in all patients.

CONCLUSIONS

Criterion standard BOLD sensitivity estimation using hypercapnic experiments is challenging, especially in patient populations. Taking into consideration the tSNR, commonly used for BOLD sensitivity estimation, but ignoring local reductions in the echo time (eg, from 22 to 11 milliseconds), would erroneously suggest functional sensitivity sufficient to map BOLD signal changes. It is therefore important to consider both local variations in the echo time and temporal variations in signal, using the product metric of these two indices for instance. This should ensure a reliable estimation of BOLD sensitivity and to facilitate the identification of potential false-negative results. This is particularly true at high fields, such as 7 T and in patients with large pathologies and postoperative cavities.

Reduced Dynamic Coupling Between Spontaneous BOLD-CBF Fluctuations in Older Adults: A Dual-Echo pCASL Study

Measurement of the dynamic coupling between spontaneous Blood Oxygenation Level Dependent (BOLD) and cerebral blood flow (CBF) fluctuations has been recently proposed as a method to probe resting-state brain physiology. Here we investigated how the dynamic BOLD-CBF coupling during resting-state is affected by aging. Fifteen young subjects and 17 healthy elderlies were studied using a dual-echo pCASL sequence. We found that the dynamic BOLD-CBF coupling was markedly reduced in elderlies, in particular in the left supramarginal gyrus, an area known to be involved in verbal working memory and episodic memory. Moreover, correcting for temporal shift between BOLD and CBF timecourses resulted in an increased correlation of the two signals for both groups, but with a larger increase for elderlies. However, even after temporal shift correction, a significantly decreased correlation was still observed for elderlies in the left supramarginal gyrus, indicating that the age-related dynamic BOLD-CBF uncoupling in this region is more pronounced and can be only partially explained with a simple time-shift between the two signals. Interestingly, these results were observed in a group of elderlies with normal cognitive functions, suggesting that the study of dynamic BOLD-CBF coupling during resting-state is a promising technique, potentially able to provide early biomarkers of functional changes in the aging brain.

Prefrontal hemodynamic after-effects caused by rebreathing may predict affective states - A multimodal functional near-infrared spectroscopy study

Brain activity has been shown to be influenced by respiratory behavior. Here, we evaluated whether respiration-induced hypo- or hypercapnia may support differentiation between physiological versus pathological respiratory behavior. In particular, we investigated whether systemic physiological measures could predict the brain's time-frequency hemodynamics after three respiratory challenges (i.e., breath-holding, rebreathing, and hyperventilation) compared to resting-state. Prefrontal hemodynamics were assessed in healthy subjects (N = 27) using functional near-infrared spectroscopy (fNIRS). Systemic physiological measures were assessed in form of heart rate, partial end-tidal carbon dioxide, respiration rate, and saturation of peripheral oxygen. Time-frequency dynamics were quantified using the wavelet transform coherence (i.e., defined here as cortical-systemic coherence). We found that the three respiratory challenges modulated cortical-systemic coherence differently: (1) After rebreathing, cortical-systemic coherence could be predicted from the amplitude of the heart rate (strong negative correlation). (2) After breath-holding, the same observation was made (moderate negative correlation). (3) After hyperventilation, no significant effect was observed. (4) These effects were found only in the frequency range of very low-frequency oscillations. The presented findings highlight a distinct role of rebreathing in predicting cortical-systemic coupling based on heart rate changes, which may represents a measure of affective states in the brain. The applied multimodal assessment of hemodynamic and systemic physiological measures during respiratory challenges may therefore have potential applications in the differentiation between physiological and pathological respiratory behavior.

Methods for cleaning the BOLD fMRI signal

Blood oxygen-level-dependent functional magnetic resonance imaging (BOLD fMRI) has rapidly become a popular technique for the investigation of brain function in healthy individuals, patients as well as in animal studies. However, the BOLD signal arises from a complex mixture of neuronal, metabolic and vascular processes, being therefore an indirect measure of neuronal activity, which is further severely corrupted by multiple non-neuronal fluctuations of instrumental, physiological or subject-specific origin. This review aims to provide a comprehensive summary of existing methods for cleaning the BOLD fMRI signal. The description is given from a methodological point of view, focusing on the operation of the different techniques in addition to pointing out the advantages and limitations in their application. Since motion-related and physiological noise fluctuations are two of the main noise components of the signal, techniques targeting their removal are primarily addressed, including both data-driven approaches and using external recordings. Data-driven approaches, which are less specific in the assumed model and can simultaneously reduce multiple noise fluctuations, are mainly based on data decomposition techniques such as principal and independent component analysis. Importantly, the usefulness of strategies that benefit from the information available in the phase component of the signal, or in multiple signal echoes is also highlighted. The use of global signal regression for denoising is also addressed. Finally, practical recommendations regarding the optimization of the preprocessing pipeline for the purpose of denoising and future venues of research are indicated. Through the review, we summarize the importance of signal denoising as an essential step in the analysis pipeline of task-based and resting state fMRI studies.

Identifying topological motif patterns of human brain functional networks

Recent imaging connectome studies demonstrated that the human functional brain network follows an efficient small-world topology with cohesive functional modules and highly connected hubs. However, the functional motif patterns that represent the underlying information flow remain largely unknown. Here, we investigated motif patterns within directed human functional brain networks, which were derived from resting-state functional magnetic resonance imaging data with controlled confounding hemodynamic latencies. We found several significantly recurring motifs within the network, including the two-node reciprocal motif and five classes of three-node motifs. These recurring motifs were distributed in distinct patterns to support intra- and inter-module functional connectivity, which also promoted integration and segregation in network organization. Moreover, the significant participation of several functional hubs in the recurring motifs exhibited their critical role in global integration. Collectively, our findings highlight the basic architecture governing brain network organization and provide insight into the information flow mechanism underlying intrinsic brain activities. Hum Brain Mapp 38:2734-2750, 2017. © 2017 Wiley Periodicals, Inc.

Multiplexed MRI methods for rapid estimation of global cerebral metabolic rate of oxygen consumption

The global cerebral metabolic rate of oxygen (CMRO₂), which reflects metabolic activity of the brain under various physiologic conditions, can be quantified using a method, referred to as 'OxFlow', which simultaneously measures hemoglobin oxygen saturation in a draining vein (Y_v) and total cerebral blood flow (tCBF). Conventional OxFlow (Conv-OxFlow) entails four interleaves incorporated in a single pulse sequence - two for phase-contrast based measurement of tCBF in the supplying arteries of the neck, and two to measure the intra- to extravascular phase difference in the superior sagittal sinus to derive Y_v [Jain et al., JCBFM 2010]. However, this approach limits achievable temporal resolution thus precluding capture of rapid changes of brain metabolic states such as the response to apneic stimuli. Here, we developed a time-efficient, multiplexed OxFlow method and evaluated its potential for measuring dynamic alterations in global CMRO₂ during a breath-hold challenge. Two different implementations of multiplexed OxFlow were investigated: 1) simultaneous-echo-refocusing based OxFlow (SER-OxFlow) and 2) simultaneous-multi-slice imaging-based dual-band OxFlow (DB-OxFlow). The two sequences were implemented on 3T scanners (Siemens TIM Trio and Prisma) and their performance was evaluated in comparison to Conv-OxFlow in ten healthy subjects for baseline CMRO₂ quantification. Comparison of measured parameters (Y_v, tCBF, CMRO₂) revealed no significant bias of SER-OxFlow and DB-OxFlow, with respect to the reference Conv-OxFlow while improving temporal resolution two-fold (12.5 versus 25s). Further acceleration shortened scan time to 8 and 6s for SER and DB-OxFlow, respectively, for time-resolved CMRO₂ measurement. The two sequences were able of capturing smooth transitions of Y_v, tCBF, and CMRO₂ over the time course consisting of 30s of normal breathing, 30s of volitional apnea, and 90s of recovery. While both SER- and DB-OxFlow techniques provide significantly improved temporal resolution (by a factor of 3 - 4 relative to Conv-OxFlow), DB-OxFlow was found to be superior for the study of short physiologic stimuli.

Imaging Posture Veils Neural Signals

Whereas modern brain imaging often demands holding body positions incongruent with everyday life, posture governs both neural activity and cognitive performance. Humans commonly perform while upright; yet, many neuroimaging methodologies require participants to remain motionless and adhere to non-ecological comportments within a confined space. This inconsistency between ecological postures and imaging constraints undermines the transferability and generalizability of many a neuroimaging assay. Here we highlight the influence of posture on brain function and behavior. Specifically, we challenge the tacit assumption that brain processes and cognitive performance are comparable across a spectrum of positions. We provide an integrative synthesis regarding the increasingly prominent influence of imaging postures on autonomic function, mental capacity, sensory thresholds, and neural activity. Arguing that neuroimagers and cognitive scientists could benefit from considering the influence posture wields on both general functioning and brain activity, we examine existing imaging technologies and the potential of portable and versatile imaging devices (e.g., functional near infrared spectroscopy). Finally, we discuss ways that accounting for posture may help unveil the complex brain processes of everyday cognition.

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.

Non-invasive quantification of absolute cerebral blood volume during functional activation applicable to the whole human brain

PURPOSE

Cerebral blood volume (CBV) changes in many diverse pathologic conditions, and in response to functional challenges along with changes in blood flow, blood oxygenation, and the cerebral metabolic rate of oxygen. The feasibility of a new method for non-invasive quantification of absolute cerebral blood volume that can be applicable to the whole human brain was investigated.

METHODS

Multi-slice data were acquired at 3 T using a novel inversion recovery echo planar imaging (IR-EPI) pulse sequence with varying contrast weightings and an efficient rotating slice acquisition order, at rest and during visual activation. A biophysical model was used to estimate absolute cerebral blood volume at rest and during activation, and oxygenation during activation, on data from 13 normal human subjects.

RESULTS

Cerebral blood volume increased by 21.7% from 6.6 ± 0.8 mL/100 mL of brain parenchyma at rest to 8.0 ± 1.3 mL/100 mL of brain parenchyma in the occipital cortex during visual activation, with average blood oxygenation of 84 ± 2.1% during activation, comparing well with literature.

CONCLUSION

The method is feasible, and could foster improved understanding of the fundamental physiological relationship between neuronal activity, hemodynamic changes, and metabolism underlying brain activation; complement existing methods for estimating compartmental changes; and potentially find utility in evaluating vascular health.

Vascular autorescaling of fMRI (VasA fMRI) improves sensitivity of population studies: A pilot study

The blood oxygenation level-dependent (BOLD) signal is widely used for functional magnetic resonance imaging (fMRI) of brain function in health and disease. The statistical power of fMRI group studies is significantly hampered by high inter-subject variance due to differences in baseline vascular physiology. Several methods have been proposed to account for physiological vascularization differences between subjects and hence improve the sensitivity in group studies. However, these methods require the acquisition of additional reference scans (such as a full resting-state fMRI session or ASL-based calibrated BOLD). We present a vascular autorescaling (VasA) method, which does not require any additional reference scans. VasA is based on the observation that slow oscillations (<0.1Hz) in arterial blood CO2 levels occur naturally due to changes in respiration patterns. These oscillations yield fMRI signal changes whose amplitudes reflect the blood oxygenation levels and underlying local vascularization and vascular responsivity. VasA estimates proxies of the amplitude of these CO2-driven oscillations directly from the residuals of task-related fMRI data without the need for reference scans. The estimates are used to scale the amplitude of task-related fMRI responses, to account for vascular differences. The VasA maps compared well to cerebrovascular reactivity (CVR) maps and cerebral blood volume maps based on vascular space occupancy (VASO) measurements in four volunteers, speaking to the physiological vascular basis of VasA. VasA was validated in a wide variety of tasks in 138 volunteers. VasA increased t-scores by up to 30% in specific brain areas such as the visual cortex. The number of activated voxels was increased by up to 200% in brain areas such as the orbital frontal cortex while still controlling the nominal false-positive rate. VasA fMRI outperformed previously proposed rescaling approaches based on resting-state fMRI data and can be readily applied to any task-related fMRI data set, even retrospectively.

Quantitative β mapping for calibrated fMRI

The metabolic and hemodynamic dependencies of the blood oxygenation level-dependent (BOLD) signal form the basis for calibrated fMRI, where the focus is on oxidative energy demanded by neural activity. An important part of calibrated fMRI is the power-law relationship between the BOLD signal and the deoxyhemoglobin concentration, which in turn is related to the ratio between oxidative demand (CMRO2) and blood flow (CBF). The power-law dependence between BOLD signal and deoxyhemoglobin concentration is signified by a scaling exponent β. Until recently most studies assumed a β value of 1.5, which is based on numerical simulations of the extravascular BOLD component. Since the basal value of CMRO2 and CBF can vary from subject-to-subject and/or region-to-region, a method to independently measure β in vivo should improve the accuracy of calibrated fMRI results. We describe a new method for β mapping through characterizing R2' - the most sensitive relaxation component of BOLD signal (i.e., the reversible magnetic susceptibility component that is predominantly of extravascular origin at high magnetic field) - as a function of intravascular magnetic susceptibility induced by an FDA-approved superparamagnetic contrast agent. In α-chloralose anesthetized rat brain, at 9.4 T, we measured β values of ~0.8 uniformly across large neocortical swathes, with lower magnitude and more heterogeneity in subcortical areas. Comparison of β maps in rats anesthetized with medetomidine and α-chloralose revealed that β is independent of neural activity levels at these resting states. We anticipate that this method for β mapping can help facilitate calibrated fMRI for clinical studies.

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.

Functional Magnetic Resonance Imaging Methods

Since its inception in 1992, Functional Magnetic Resonance Imaging (fMRI) has become an indispensible tool for studying cognition in both the healthy and dysfunctional brain. FMRI monitors changes in the oxygenation of brain tissue resulting from altered metabolism consequent to a task-based evoked neural response or from spontaneous fluctuations in neural activity in the absence of conscious mentation (the "resting state"). Task-based studies have revealed neural correlates of a large number of important cognitive processes, while fMRI studies performed in the resting state have demonstrated brain-wide networks that result from brain regions with synchronized, apparently spontaneous activity. In this article, we review the methods used to acquire and analyze fMRI signals.

Time-frequency dynamics of the sum of intra- and extracerebral hemodynamic functional connectivity during resting-state and respiratory challenges assessed by multimodal functional near-infrared spectroscopy

Monitoring respiratory processes is important for evaluating neuroimaging data, given their influence on time-frequency dynamics of intra- and extracerebral hemodynamics. Here we investigated the time-frequency dynamics of the sum of intra- and extracerebral hemodynamic functional connectivity states during hypo- and hypercapnia by using three different respiratory challenge tasks (i.e., hyperventilation, breath-holding, and rebreathing) compared to resting-state. The sum of intra- and extracerebral hemodynamic responses were assessed using functional near-infrared spectroscopy (fNIRS) within two regions of interest (i.e., the dorsolateral and the medial prefrontal cortex). Time-frequency fNIRS analysis was performed based on wavelet transform coherence to quantify functional connectivity in terms of positive and negative phase-coupling within each region of interest. Physiological measures were assessed in the form of partial end-tidal carbon dioxide, heart rate, arterial tissue oxygen saturation, and respiration rate. We found that the three respiration challenges modulated time-frequency dynamics differently with respect to resting-state: 1) Hyperventilation and breath-holding exhibited inverse patterns of positive and negative phase-coupling. 2) In contrast, rebreathing had no significant effect. 3) Low-frequency oscillations contributed to a greater extent to time-frequency dynamics compared to high-frequency oscillations. The results highlight that there exist distinct differences in time-frequency dynamics of the sum of intra- and extracerebral functional connectivity not only between hypo- (hyperventilation) and hypercapnia but also between different states of hypercapnia (breath-holding versus rebreathing). This suggests that a multimodal assessment of intra-/extracerebral and systemic physiological changes during respiratory challenges compared to resting-state may have potential use in the differentiation between physiological and pathological respiratory behavior accompanied by the psycho-physiological state of a human.

Validation of the hypercapnic calibrated fMRI method using DOT-fMRI fusion imaging

Calibrated functional magnetic resonance imaging (fMRI) is a widely used method to investigate brain function in terms of physiological quantities such as the cerebral metabolic rate of oxygen (CMRO2). The first and one of the most common methods of fMRI calibration is hypercapnic calibration. This is achieved via simultaneous measures of the blood-oxygenation-level dependent (BOLD) and the arterial spin labeling (ASL) signals during a functional task that evokes regional changes in CMRO2. A subsequent acquisition is then required during which the subject inhales carbon dioxide for short periods of time. A calibration constant, typically labeled M, is then estimated from the hypercapnic data and is subsequently used together with the BOLD-ASL recordings to compute evoked changes in CMRO2 during the functional task. The computation of M assumes a constant CMRO2 during the CO2 inhalation, an assumption that has been questioned since the origin of calibrated fMRI. In this study we used diffuse optical tomography (DOT) together with BOLD and ASL--an alternative calibration method that does not require any gas manipulation and therefore no constant CMRO2 assumption--to cross-validate the estimation of M obtained from a traditional hypercapnic calibration. We found a high correlation between the M values (R=0.87, p<0.01) estimated using these two approaches. The findings serve to validate the hypercapnic fMRI calibration technique and suggest that the inter-subject variability routinely obtained for M is reproducible with an alternative method and might therefore reflect inter-subject physiological variability.

Noninvasive MRI measurement of the absolute cerebral blood volume-cerebral blood flow relationship during visual stimulation in healthy humans

PURPOSE

The relationship between cerebral blood volume (CBV) and cerebral blood flow (CBF) underlies blood oxygenation level-dependent functional MRI signal. This study investigates the potential for improved characterization of the CBV-CBF relationship in humans, and examines sex effects as well as spatial variations in the CBV-CBF relationship.

METHODS

Healthy subjects were imaged noninvasively at rest and during visual stimulation, constituting the first MRI measurement of the absolute CBV-CBF relationship in humans with complete coverage of the functional areas of interest.

RESULTS

CBV and CBF estimates were consistent with the literature, and their relationship varied both spatially and with sex. In a region of interest with stimulus-induced activation in CBV and CBF at a significance level of the P < 0.05, a power function fit resulted in CBV = 2.1 CBF(0.32) across all subjects, CBV = 0.8 CBF(0.51) in females and CBV = 4.4 CBF(0.15) in males. Exponents decreased in both sexes as ROIs were expanded to include less significantly activated regions.

CONCLUSION

Consideration for potential sex-related differences, as well as regional variations under a range of physiological states, may reconcile some of the variation across literature and advance our understanding of the underlying cerebrovascular physiology.

Functional imaging of cerebral perfusion

The functional imaging of perfusion enables the study of its properties such as the vasoreactivity to circulating gases, the autoregulation and the neurovascular coupling. Downstream from arterial stenosis, this imaging can estimate the vascular reserve and the risk of ischemia in order to adapt the therapeutic strategy. This method reveals the hemodynamic disorders in patients suffering from Alzheimer's disease or with arteriovenous malformations revealed by epilepsy. Functional MRI of the vasoreactivity also helps to better interpret the functional MRI activation in practice and in clinical research.

A review of calibrated blood oxygenation level-dependent (BOLD) methods for the measurement of task-induced changes in brain oxygen metabolism

The dynamics of the blood oxygenation level-dependent (BOLD) response are dependent on changes in cerebral blood flow, cerebral blood volume and the cerebral metabolic rate of oxygen consumption. Furthermore, the amplitude of the response is dependent on the baseline physiological state, defined by the haematocrit, oxygen extraction fraction and cerebral blood volume. As a result of this complex dependence, the accurate interpretation of BOLD data and robust intersubject comparisons when the baseline physiology is varied are difficult. The calibrated BOLD technique was developed to address these issues. However, the methodology is complex and its full promise has not yet been realised. In this review, the theoretical underpinnings of calibrated BOLD, and issues regarding this theory that are still to be resolved, are discussed. Important aspects of practical implementation are reviewed and reported applications of this methodology are presented.

End-tidal CO2: an important parameter for a correct interpretation in functional brain studies using speech tasks

The aim was to investigate the effect of different speech tasks, i.e. recitation of prose (PR), alliteration (AR) and hexameter (HR) verses and a control task (mental arithmetic (MA) with voicing of the result on end-tidal CO2 (PETCO2), cerebral hemodynamics and oxygenation. CO2 levels in the blood are known to strongly affect cerebral blood flow. Speech changes breathing pattern and may affect CO2 levels. Measurements were performed on 24 healthy adult volunteers during the performance of the 4 tasks. Tissue oxygen saturation (StO2) and absolute concentrations of oxyhemoglobin ([O2Hb]), deoxyhemoglobin ([HHb]) and total hemoglobin ([tHb]) were measured by functional near-infrared spectroscopy (fNIRS) and PETCO2 by a gas analyzer. Statistical analysis was applied to the difference between baseline before the task, 2 recitation and 5 baseline periods after the task. The 2 brain hemispheres and 4 tasks were tested separately. A significant decrease in PETCO2 was found during all 4 tasks with the smallest decrease during the MA task. During the recitation tasks (PR, AR and HR) a statistically significant (p<0.05) decrease occurred for StO2 during PR and AR in the right prefrontal cortex (PFC) and during AR and HR in the left PFC. [O2Hb] decreased significantly during PR, AR and HR in both hemispheres. [HHb] increased significantly during the AR task in the right PFC. [tHb] decreased significantly during HR in the right PFC and during PR, AR and HR in the left PFC. During the MA task, StO2 increased and [HHb] decreased significantly during the MA task. We conclude that changes in breathing (hyperventilation) during the tasks led to lower CO2 pressure in the blood (hypocapnia), predominantly responsible for the measured changes in cerebral hemodynamics and oxygenation. In conclusion, our findings demonstrate that PETCO2 should be monitored during functional brain studies investigating speech using neuroimaging modalities, such as fNIRS, fMRI to ensure a correct interpretation of changes in hemodynamics and oxygenation.

Biophysical and physiological origins of blood oxygenation level-dependent fMRI signals

After its discovery in 1990, blood oxygenation level-dependent (BOLD) contrast in functional magnetic resonance imaging (fMRI) has been widely used to map brain activation in humans and animals. Since fMRI relies on signal changes induced by neural activity, its signal source can be complex and is also dependent on imaging parameters and techniques. In this review, we identify and describe the origins of BOLD fMRI signals, including the topics of (1) effects of spin density, volume fraction, inflow, perfusion, and susceptibility as potential contributors to BOLD fMRI, (2) intravascular and extravascular contributions to conventional gradient-echo and spin-echo BOLD fMRI, (3) spatial specificity of hemodynamic-based fMRI related to vascular architecture and intrinsic hemodynamic responses, (4) BOLD signal contributions from functional changes in cerebral blood flow (CBF), cerebral blood volume (CBV), and cerebral metabolic rate of O(2) utilization (CMRO(2)), (5) dynamic responses of BOLD, CBF, CMRO(2), and arterial and venous CBV, (6) potential sources of initial BOLD dips, poststimulus BOLD undershoots, and prolonged negative BOLD fMRI signals, (7) dependence of stimulus-evoked BOLD signals on baseline physiology, and (8) basis of resting-state BOLD fluctuations. These discussions are highly relevant to interpreting BOLD fMRI signals as physiological means.

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.

Quantitative measurement of cerebral physiology using respiratory-calibrated MRI

Functional magnetic resonance imaging typically measures signal increases arising from changes in the transverse relaxation rate over small regions of the brain and associates these with local changes in cerebral blood flow, blood volume and oxygen metabolism. Recent developments in pulse sequences and image analysis methods have improved the specificity of the measurements by focussing on changes in blood flow or changes in blood volume alone. However, FMRI is still unable to match the physiological information obtainable from positron emission tomography (PET), which is capable of quantitative measurements of blood flow and volume, and can indirectly measure resting metabolism. The disadvantages of PET are its cost, its availability, its poor spatial resolution and its use of ionising radiation. The MRI techniques introduced here address some of these limitations and provide physiological data comparable with PET measurements. We present an 18-minute MRI protocol that produces multi-slice whole-brain coverage and yields quantitative images of resting cerebral blood flow, cerebral blood volume, oxygen extraction fraction, CMRO(2), arterial arrival time and cerebrovascular reactivity of the human brain in the absence of any specific functional task. The technique uses a combined hyperoxia and hypercapnia paradigm with a modified arterial spin labelling sequence.

Functional connectivity MRI in infants: exploration of the functional organization of the developing brain

Advanced neuroimaging techniques have been increasingly applied to the study of preterm and term infants in an effort to further define the functional cerebral architecture of the developing brain. Despite improved understanding of the complex relationship between structure and function obtained through these investigations, significant questions remain regarding the nature, location, and timing of the maturational changes which occur during early development. Functional connectivity magnetic resonance imaging (fcMRI) utilizes spontaneous, low frequency (< 0.1 Hz), coherent fluctuations in blood oxygen level dependent (BOLD) signal to identify networks of functional cerebral connections. Due to the intrinsic characteristics of its image acquisition and analysis, fcMRI offers a novel neuroimaging approach well suited to investigation of infants. Recently, this methodology has been successfully applied to examine neonatal populations, defining normative patterns of large-scale neural network development in the maturing brain. The resting-state networks (RSNs) identified in these studies reflect the evolving cerebral structural architecture, presumably driven by varied genetic and environmental influences. Principal features of these investigations and their role in characterization of the tenets of neural network development during this critical developmental period are highlighted in this review. Despite these successes, optimal methods for fcMRI data acquisition and analysis for this population have not yet been defined. Further, appropriate schemes for interpretation and translation of fcMRI results remain unknown, a matter of increasing importance as functional neuroimaging findings are progressively applied in the clinical arena. Notwithstanding these concerns, fcMRI provides insight into the earliest forms of cerebral connectivity and therefore holds great promise for future neurodevelopmental investigations.

Increasing measurement accuracy of age-related BOLD signal change: minimizing vascular contributions by resting-state-fluctuation-of-amplitude scaling

In this report we demonstrate a hemodynamic scaling method with resting-state fluctuation of amplitude (RSFA) in healthy adult younger and older subject groups. We show that RSFA correlated with breath hold (BH) responses throughout the brain in groups of younger and older subjects which RSFA and BH performed comparably in accounting for age-related hemodynamic coupling changes, and yielded more veridical estimates of age-related differences in task-related neural activity. BOLD data from younger and older adults performing motor and cognitive tasks were scaled using RSFA and BH related signal changes. Scaling with RSFA and BH reduced the skew of the BOLD response amplitude distribution in each subject and reduced mean BOLD amplitude and variability in both age groups. Statistically significant differences in intrasubject amplitude variation across regions of activated cortex, and intersubject amplitude variation in regions of activated cortex were observed between younger and older subject groups. Intra- and intersubject variability differences were mitigated after scaling. RSFA, though similar to BH in minimizing skew in the unscaled BOLD amplitude distribution, attenuated the neural activity-related BOLD amplitude significantly less than BH. The amplitude and spatial extent of group activation were lower in the older than in the younger group before and after scaling. After accounting for vascular variability differences through scaling, age-related decreases in activation volume were observed during the motor and cognitive tasks. The results suggest that RSFA-scaled data yield age-related neural activity differences during task performance with negligible effects from non-neural (i.e., vascular) sources.

Quantification of CMRO(2) without hypercapnia using simultaneous near-infrared spectroscopy and fMRI measurements

Estimation of the cerebral metabolic rate of oxygen (CMRO(2)) and cerebral blood flow (CBF) is important to investigate the neurovascular coupling and physiological components in blood oxygenation level-dependent (BOLD) signals quantitatively. Although there are methods that can determine CMRO(2) changes using functional MRI (fMRI) or near-infrared spectroscopy (NIRS), current approaches require a separate hypercapnia calibration process and have the potential to incur bias in many assumed model parameters. In this paper, a novel method to estimate CMRO(2) without hypercapnia is described using simultaneous measurements of NIRS and fMRI. Specifically, an optimization framework is proposed that minimizes the differences between the two forms of the relative CMRO(2)-CBF coupling ratio from BOLD and NIRS biophysical models, from which hypercapnia calibration and model parameters are readily estimated. Based on the new methods, we found that group average CBF, CMRO(2), cerebral blood volume (CBV), and BOLD changes within activation of the primary motor cortex during a finger tapping task increased by 39.5 +/- 21.4%, 18.4 +/- 8.7%, 12.9 +/- 6.7%, and 0.5 +/- 0.2%, respectively. The group average estimated flow-metabolism coupling ratio was 2.38 +/- 0.65 and the hypercapnia parameter was 7.7 +/- 1.7%. These values are within the range of values reported from other literatures. Furthermore, the activation maps from CBF and CMRO(2) were well localized on the primary motor cortex, which is the main target region of the finger tapping task.

Neural and vascular variability and the fMRI-BOLD response in normal aging

Neural, vascular and structural variables contributing to the blood oxygen level-dependent (BOLD) signal response variability were investigated in younger and older humans. Twelve younger healthy human subjects (six male and six female; mean age: 24 years; range: 19-27 years) and 12 older healthy subjects (five male and seven female; mean age: 58 years; range: 55-71 years) with no history of head trauma and neurological disease were scanned. Functional magnetic resonance imaging measurements using the BOLD contrast were made when participants performed a motor, cognitive or a breath hold (BH) task. Activation volume and the BOLD response amplitude were estimated for the younger and older at both group and subject levels. Mean activation volume was reduced by 45%, 40% and 38% in the elderly group during the motor, cognitive and BH tasks, respectively, compared to the younger. Reduction in activation volume was substantially higher compared to the reduction in the gray matter volume of 14% in the older compared to the younger. A significantly larger variability in the intersubject BOLD signal change occurred during the motor task, compared to the cognitive task. BH-induced BOLD signal change between subjects was significantly less-variable in the motor task-activated areas in the younger compared to older whereas such a difference between age groups was not observed during the cognitive task. Hemodynamic scaling using the BH signal substantially reduced the BOLD signal variability during the motor task compared to the cognitive task. The results indicate that the origin of the BOLD signal variability between subjects was predominantly vascular during the motor task while being principally a consequence of neural variability during the cognitive task. Thus, in addition to gray matter differences, the type of task performed can have different vascular variability weighting that can influence age-related differences in brain functional response.

Effects of aging on cerebral blood flow, oxygen metabolism, and blood oxygenation level dependent responses to visual stimulation

Calibrated functional magnetic resonance imaging (fMRI) provides a noninvasive technique to assess functional metabolic changes associated with normal aging. We simultaneously measured both the magnitude of the blood oxygenation level dependent (BOLD) and cerebral blood flow (CBF) responses in the visual cortex for separate conditions of mild hypercapnia (5% CO(2)) and a simple checkerboard stimulus in healthy younger (n = 10, mean: 28-years-old) and older (n = 10, mean: 53-years-old) adults. From these data we derived baseline CBF, the BOLD scaling parameter M, the fractional change in the cerebral metabolic rate of oxygen consumption (CMRO(2)) with activation, and the coupling ratio n of the fractional changes in CBF and CMRO(2). For the functional activation paradigm, the magnitude of the BOLD response was significantly lower for the older group (0.57 +/- 0.07%) compared to the younger group (0.95 +/- 0.14%), despite the finding that the fractional CBF and CMRO(2) changes were similar for both groups. The weaker BOLD response for the older group was due to a reduction in the parameter M, which was significantly lower for older (4.6 +/- 0.4%) than younger subjects (6.5 +/- 0.8%), most likely reflecting a reduction in baseline CBF for older (41.7 +/- 4.8 mL/100 mL/min) compared to younger (59.6 +/- 9.1 mL/100 mL/min) subjects. In addition to these primary responses, for both groups the BOLD response exhibited a post-stimulus undershoot with no significant difference in this magnitude. However, the post-undershoot period of the CBF response was significantly greater for older compared to younger subjects. We conclude that when comparing two populations, the BOLD response can provide misleading reflections of underlying physiological changes. A calibrated approach provides a more quantitative reflection of underlying metabolic changes than the BOLD response alone.

Comparison of hypercapnia-based calibration techniques for measurement of cerebral oxygen metabolism with MRI

MRI may be used to measure fractional changes in cerebral oxygen metabolism via a metabolic model. One step commonly used in this measurement is calibration with image data acquired during hypercapnia, which is a state of increased CO2 content of the blood. In this study some commonly used hypercapnia-inducing stimuli were compared to assess their suitability for the calibration step. The following stimuli were investigated: (a) inspiration of a mixture of 4% CO2, 21% O2 and balance N2; (b) 30-s breath holding; and (c) inspiration of a mixture of 4% CO2 and 96% O2 (i.e., carbogen). Measurements of BOLD and cerebral blood flow made on nine subjects during the different hypercapnia-inducing stimuli showed that each stimulus leads to a different calibration of the model. We argue that of the aforementioned stimuli, inspiration of 4% CO2, 21% O2 and balance N2 should be preferred for the calibration as the other stimuli produce responses that violate assumptions of the metabolic model.

Influence of heart rate on the BOLD signal: the cardiac response function

It has previously been shown that low-frequency fluctuations in both respiratory volume and cardiac rate can induce changes in the blood-oxygen level dependent (BOLD) signal. Such physiological noise can obscure the detection of neural activation using fMRI, and it is therefore important to model and remove the effects of this noise. While a hemodynamic response function relating respiratory variation (RV) and the BOLD signal has been described [Birn, R.M., Smith, M.A., Jones, T.B., Bandettini, P.A., 2008b. The respiration response function: The temporal dynamics of fMRI signal fluctuations related to changes in respiration. Neuroimage 40, 644-654.], no such mapping for heart rate (HR) has been proposed. In the current study, the effects of RV and HR are simultaneously deconvolved from resting state fMRI. It is demonstrated that a convolution model including RV and HR can explain significantly more variance in gray matter BOLD signal than a model that includes RV alone, and an average HR response function is proposed that well characterizes our subject population. It is observed that the voxel-wise morphology of the deconvolved RV responses is preserved when HR is included in the model, and that its form is adequately modeled by Birn et al.'s previously-described respiration response function. Furthermore, it is shown that modeling out RV and HR can significantly alter functional connectivity maps of the default-mode network.

Mapping and correction of vascular hemodynamic latency in the BOLD signal

Correlation and causality metrics can be applied to blood-oxygen level-dependent (BOLD) signal time series in order to infer neural synchrony and directions of information flow from fMRI data. However, the BOLD signal reflects both the underlying neural activity and the vascular response, the latter of which is governed by local vasomotor physiology. The presence of potential vascular latency differences thus poses a confound in the detection of neural synchrony as well as inferences about the causality of neural processes. In the present study, we investigate the use of a breath holding (BH) task for characterizing and correcting for voxel-wise neurovascular latency differences across the whole brain. We demonstrate that BH yields reliable measurements of relative timing differences between voxels, and further show that a BH-derived latency correction can impact both functional connectivity maps of the resting-state default-mode network and activation maps of an event-related working memory (WM) task.

Detection and scaling of task-induced fMRI-BOLD response using resting state fluctuations

This study evaluated a calibration technique for scaling the fMRI-BOLD response during a simple motor task. A novel scaling parameter, resting state physiological fluctuation amplitude (RSFA), was tested using previously established scaling factors such as breath hold or 5% CO(2). RSFA was hypothesized to contain vascular reactivity information present in the resting state fMRI signal. Subjects were scanned under various stimulus conditions: (a) rest while breathing room air, (b) bilateral fingertapping, (c) breath holding and (d) moderate hypercapnia (breathing 5% CO(2)+air). In all subjects who breathed 5% CO(2), RSFA correlated highly with the BOLD response amplitude during 5% CO(2) inhalation. Also, RSFA correlated highly with the amplitude of the BOLD response elicited by breath hold. RSFA was therefore used as a hemodynamic scaling factor to calibrate both the amplitude and spatial extent of the fMRI-BOLD response during the motor task (fingertapping). Results revealed that amplitude scaling using RSFA was similar to that using breath hold or 5% CO(2), where the spatial extent of activation diminished by 20-37% over all subjects. Spatial extent of activation changed significantly after scaling and only 30-40% of the activated area overlapped with the unscaled activation. RSFA-scaled task-induced fMRI-BOLD response in both amplitude and spatial extent was comparable to that obtained using breath hold or 5% CO(2). We conclude that RSFA may be used to hemodynamically scale the fMRI-BOLD response and does not require the use of a hypercapnic challenge (which may not be purely non-neural), which can be difficult to implement in special populations.

Multimodal investigation of fMRI and fNIRS derived breath hold BOLD signals with an expanded balloon model

Multimodal investigation of blood oxygenation level-dependent (BOLD) signals, using both functional near-infrared spectroscopy (fNIRS) and functional magnetic resonance imaging (fMRI), may give further insight to the underlying physiological principles and the detailed transient dynamics of the vascular response. Utilizing a breath hold task (BHT), we measured deoxy-hemoglobin (HbR) and oxy-hemoglobin (HbO) changes via fNIRS and blood oxygen level dependent (BOLD) changes by fMRI. Measurements were taken in four volunteers asynchronously and carefully aligned for comparative analysis. In order to describe the main stimulus in BHT, partial pressure of carbon dioxide (PaCO(2)) parameter was integrated into the balloon model as the driving function of cerebral blood flow (CBF) which led to the development of an expanded balloon model (EBM). During BHT, the increase in HbR was observed later than the BOLD peak and coincided temporally with its post-stimulus undershoot. Further investigation of these transients with a PaCO(2) integrated balloon model suggests that post-stimulus undershoot measured by fMRI is dominated by slow return of cerebral blood volume (CBV). This was confirmed by fNIRS measurements. In addition, the BOLD signal decreased with the increase of the initial level of PaCO(2) derived from EBM, indicating an effect of basal CBF level on the BOLD signal. In conclusion, a multimodal approach with an appropriate biophysical model gave a comprehensive description of the hemodynamic response during BHT.

Regional differences in the coupling of cerebral blood flow and oxygen metabolism changes in response to activation: implications for BOLD-fMRI

Functional magnetic resonance imaging (fMRI) based on blood oxygenation level dependent (BOLD) signal changes is a sensitive tool for mapping brain activation, but quantitative interpretation of the BOLD response is problematic. The BOLD response is primarily driven by cerebral blood flow (CBF) changes, but is moderated by M, a scaling parameter reflecting baseline deoxyhemoglobin, and n, the ratio of fractional changes in CBF to cerebral metabolic rate of oxygen consumption (CMRO(2)). We compared M and n between cortical (visual cortex, VC) and subcortical (lentiform nuclei, LN) regions using a quantitative approach based on calibrating the BOLD response with a hypercapnia experiment. Although M was similar in both regions (~5.8%), differences in n (2.21+/-0.03 in VC and 1.58+/-0.03 in LN; Cohen d=1.71) produced substantially weaker (~3.7x) subcortical than cortical BOLD responses relative to CMRO(2) changes. Because of this strong sensitivity to n, BOLD response amplitudes cannot be interpreted as a quantitative reflection of underlying metabolic changes, particularly when comparing cortical and subcortical regions.

A calibration method for quantitative BOLD fMRI based on hyperoxia

The estimation of changes in CMR(O2) using functional MRI involves an essential calibration step using a vasoactive agent to induce an isometabolic change in CBF. This calibration procedure is performed most commonly using hypercapnia as the isometabolic stimulus. However, hypercapnia possesses a number of detrimental side effects. Here, a new method is presented using hyperoxia to perform the same calibration step. This procedure requires independent measurement of Pa(O2), the BOLD signal, and CBF. We demonstrate that this method yields results that are comparable to those derived using other methods. Further, the hyperoxia technique is able to provide an estimate of the calibration constant that has lower overall intersubject and intersession variability compared to the hypercapnia approach.

Differential effects of the D- and L- isomers of amphetamine on pharmacological MRI BOLD contrast in the rat

RATIONALE

The D - and L-amphetamine sulphate isomers are used in the formulation of Adderall XR(R), which is effective in the treatment of attention-deficit hyperactivity disorder (ADHD). The effects of these isomers on brain activity has not been examined using neuroimaging.

OBJECTIVES

This study determines the pharmacological magnetic resonance imaging blood-oxygenation-level-dependent (BOLD) response in rat brain regions after administration of each isomer.

MATERIALS AND METHODS

Rats were individually placed into a 2.35 T Bruker magnet for 60 min to achieve basal recording of variation in signal intensity. Either saline (n = 9), D-amphetamine sulphate (2 mg/kg, i.p.; n = 9) or L: -amphetamine sulphate (4 mg/kg, i.p.; n = 9) were administered, and recording continued for a further 90 min. Data were analysed for BOLD effects using statistical parametric maps. Blood pressure, blood gases and respiratory rate were monitored during scanning.

RESULTS

The isomers show overlapping effects on the BOLD responses in areas including nucleus accumbens, medial entorhinal cortex, colliculi, field CA1 of hippocampus and thalamic nuclei. The L-isomer produced greater global changes in the positive BOLD response than the D-isomer, including the somatosensory and motor cortices and frontal brain regions such as the orbitofrontal cortices, prelimbic and infralimbic cortex which were not observed with the D-isomer.

CONCLUSIONS

The amphetamine isomers produce different BOLD responses in brain areas related to cognition, pleasure, pain processing and motor control probably because of variations on brain amine systems such as dopamine and noradrenaline. The isomers may, therefore, have distinct actions on brain regions affected in ADHD patients.

Hemodynamic scaling of fMRI-BOLD signal: validation of low-frequency spectral amplitude as a scalability factor

Functional magnetic resonance imaging blood-oxygenation-level-dependent (fMRI-BOLD) signal representing neural activity may be optimized by discriminating MR signal components related to neural activity and those related to intrinsic properties of the cortical vasculature. The objective of this study was to reduce the hemodynamic change independent of neural activity to obtain a scaled fMRI-BOLD response using two factors, namely, low-frequency spectral amplitude (LFSA) and breath-hold amplitude (BHA). Ten subjects (age range, 22-38 years) were scanned during four task conditions: (a) rest while breathing room air, (b) bilateral finger tapping while breathing room air, (c) rest during a partial inspirational breath-hold, and (d) rest during moderate hypercapnia (breathing 5% CO2, 20% O2 and 75% N2). In all subjects who breathed 5% CO2, regions with significant BOLD response during breath-hold correlated significantly with the percent signal increase during 5% CO2 inhalation. Finger-tapping-induced responses in the motor cortex were diminished to a similar extent after scaling using either LFSA or BHA. Inter- and intrasubject variation in the amplitude of the BOLD signal response reduced after hemodynamic scaling using LFSA or BHA. The results validated the hemodynamic amplitude scaling using LFSA with the earlier established BHA. LFSA free from motor-task contamination can be used to calibrate the fMRI-BOLD response in lieu of BHA or hypercapnia to minimize intra- and intersubject variation arising from vascular anatomy and vasodilative capacity.

Calibration of BOLD fMRI using breath holding reduces group variance during a cognitive task

The proportionality of blood oxygen level-dependent (BOLD) response during a cognitive task and that from a hypercapnic challenge was investigated in cortical structures involved in working memory (WM). Breath holding (BH) following inspiration was used to induce a BOLD response characteristic of regional vasomotor reactivity but devoid of metabolic changes. BOLD effects measured during BH were used to normalize individual subject activations during WM, which effectively reduced the confounding influence of individual- and region-specific differences in hemodynamic responsivity common to both tasks. In a study of seven subjects, the BH calibration reduced intersubject variability in WM effect amplitude by 24.8% (P < 0.03). Reduced intersubject variability resulted in a 23.7% increase in group WM activation voxel extent significant at P < 0.001, with further increases at more stringent thresholds. Because the BH task does not require CO(2) inhalation or other invasive manipulations and is broadly applicable across cortical regions, the proposed approach is simple to implement and may be beneficial for use not only in quantitative group fMRI analyses, but also for multicenter and longitudinal studies.

Increases in oxygen consumption without cerebral blood volume change during visual stimulation under hypotension condition

The magnitude of the blood oxygenation level-dependent (BOLD) signal depends on cerebral blood flow (CBF), cerebral blood volume (CBV) and cerebral metabolic rate of oxygen (CMRO2). Thus, it is difficult to separate CMRO2 changes from CBF and CBV changes. To detect the BOLD signal changes induced only by CMRO2 responses without significant evoked CBF and CBV changes, BOLD and CBV functional magnetic resonance imaging (fMRI) responses to visual stimulation were measured under normal and hypotension conditions in isoflurane-anesthetized cats at 4.7 T. When the mean arterial blood pressure (MABP) decreased from 89+/-10 to 50+/-1 mm Hg (mean+/-standard deviation, n=5) by infusion of vasodilator sodium nitroprusside, baseline CBV in the visual cortex increased by 28.4%+/-8.3%. The neural activity-evoked CBV increase in the visual cortex was 10.8%+/-3.9% at normal MABP, but was negligible at hypotension. Positive BOLD changes of +1.8%+/-0.5% (gradient echo time=25 ms) at normal MABP condition became prolonged negative changes of -1.2%+/-0.3% at hypotension. The negative BOLD response at hypotension starts approximately 1 sec earlier than positive BOLD response, but similar to CBV change at normal MABP condition. Our finding shows that the negative BOLD signals in an absence of CBV changes are indicative of an increase in CMRO2. The vasodilator-induced hypotension model simplifies the physiological source of the BOLD fMRI signals, providing an insight into spatial and temporal CMRO2 changes.