Physiological origin for the BOLD poststimulus undershoot in human brain: vascular compliance versus oxygen metabolism

Enhanced parameter estimation in multiparametric arterial spin labeling using artificial neural networks

PURPOSE

Multiparametric arterial spin labeling (MP-ASL) can quantify cerebral blood flow (CBF) and arterial cerebral blood volume (CBVa). However, its accuracy is compromised owing to its intrinsically low SNR, necessitating complex and time-consuming parameter estimation. Deep neural networks (DNNs) offer a solution to these limitations. Therefore, we aimed to develop simulation-based DNNs for MP-ASL and compared the performance of a supervised DNN (DNNSup), physics-informed unsupervised DNN (DNNUns), and the conventional lookup table method (LUT) using simulation and in vivo data.

METHODS

MP-ASL was performed twice during resting state and once during the breath-holding task. First, the accuracy and noise immunity were evaluated in the first resting state. Second, CBF and CBVa values were statistically compared between the first resting state and the breath-holding task using the Wilcoxon signed-rank test and Cliff's delta. Finally, reproducibility of the two resting states was assessed.

RESULTS

Simulation and first resting-state analyses demonstrated that DNNSup had higher accuracy, noise immunity, and a six-fold faster computation time than LUT. Furthermore, all methods detected task-induced CBF and CBVa elevations, with the effect size being larger with the DNNSup (CBF, p = 0.055, Δ = 0.286; CBVa, p = 0.008, Δ = 0.964) and DNNUns (CBF, p = 0.039, Δ = 0.286; CBVa, p = 0.008, Δ = 1.000) than that with LUT (CBF, p = 0.109, Δ = 0.214; CBVa, p = 0.008, Δ = 0.929). Moreover, all the methods exhibited comparable and satisfactory reproducibility.

CONCLUSION

DNNSup outperforms DNNUns and LUT with respect to estimation performance and computation time.

Scale-free dynamics in the core-periphery topography and task alignment decline from conscious to unconscious states

Scale-free physiological processes are ubiquitous in the human organism. Resting-state functional MRI studies observed the loss of scale-free dynamics under anesthesia. In contrast, the modulation of scale-free dynamics during task-related activity remains an open question. We investigate scale-free dynamics in the cerebral cortex's unimodal periphery and transmodal core topography in rest and task states during three conscious levels (awake, sedation, and anesthesia) complemented by computational modelling (Stuart-Landau model). The empirical findings demonstrate that the loss of the brain's intrinsic scale-free dynamics in the core-periphery topography during anesthesia, where pink noise transforms into white noise, disrupts the brain's neuronal alignment with the task's temporal structure. The computational model shows that the stimuli's scale-free dynamics, namely pink noise distinguishes from brown and white noise, also modulate task-related activity. Together, we provide evidence for two mechanisms of consciousness, temporo-spatial nestedness and alignment, suggested by the Temporo-Spatial Theory of Consciousness (TTC).

Photoacoustic imaging of squirrel monkey cortical and subcortical brain regions during peripheral electrical stimulation

The investigation of neuronal activity in non-human primate models is of critical importance due to their genetic similarity to human brains. In this study, we tested the feasibility of using photoacoustic imaging for the detection of cortical and subcortical responses due to peripheral electrical stimulation in a squirrel monkey model. Photoacoustic computed tomography and photoacoustic microscopy were applied on squirrel monkeys for real-time deep subcortical imaging and optical-resolution cortical imaging, respectively. The electrically evoked hemodynamic changes in primary somatosensory cortex, premotor cortices, primary motor cortex, and underlying subcortical areas were measured. Hemodynamic responses were observed in both cortical and subcortical brain areas at the cortices during external stimulation, demonstrating the feasibility of photoacoustic technique for functional imaging of non-human primate brain.

Multi-vendor and multisite evaluation of cerebrovascular reactivity mapping using hypercapnia challenge

Cerebrovascular reactivity (CVR), which measures the ability of cerebral blood vessels to dilate or constrict in response to vasoactive stimuli such as CO2 inhalation, is an important index of the brain's vascular health. Quantification of CVR using BOLD MRI with hypercapnia challenge has shown great promises in research and clinical studies. However, in order for it to be used as a potential imaging biomarker in large-scale and multi-site studies, the reliability of CO2-CVR quantification across different MRI acquisition platforms and researchers/raters must be examined. The goal of this report from the MarkVCID small vessel disease biomarkers consortium is to evaluate the reliability of CO2-CVR quantification in three studies. First, the inter-rater reliability of CO2-CVR data processing was evaluated by having raters from 5 MarkVCID sites process the same 30 CVR datasets using a cloud-based CVR data processing pipeline. Second, the inter-scanner reproducibility of CO2-CVR quantification was assessed in 10 young subjects across two scanners of different vendors. Third, test-retest repeatability was evaluated in 20 elderly subjects from 4 sites with a scan interval of less than 2 weeks. In all studies, the CO2 CVR measurements were performed using the fixed inspiration method, where the subjects wore a nose clip and a mouthpiece and breathed room air and 5% CO2 air contained in a Douglas bag alternatively through their mouth. The results showed that the inter-rater CoV of CVR processing was 0.08 ± 0.08% for whole-brain CVR values and ranged from 0.16% to 0.88% in major brain regions, with ICC of absolute agreement above 0.9959 for all brain regions. Inter-scanner CoV was found to be 6.90 ± 5.08% for whole-brain CVR values, and ranged from 4.69% to 12.71% in major brain regions, which are comparable to intra-session CoVs obtained from the same scanners on the same day. ICC of consistency between the two scanners was 0.8498 for whole-brain CVR and ranged from 0.8052 to 0.9185 across major brain regions. In the test-retest evaluation, test-retest CoV across different days was found to be 18.29 ± 17.12% for whole-brain CVR values, and ranged from 16.58% to 19.52% in major brain regions, with ICC of absolute agreement ranged from 0.6480 to 0.7785. These results demonstrated good inter-rater, inter-scanner, and test-retest reliability in healthy volunteers, and suggested that CO2-CVR has suitable instrumental properties for use as an imaging biomarker of cerebrovascular function in multi-site and longitudinal observational studies and clinical trials.

Separating spin compartments in arterial spin labeling using delays alternating with nutation for tailored excitation (DANTE) pulse: A validation study using T2 -relaxometry and application to arterial cerebral blood volume imaging

PURPOSE

To clarify the type of spin compartment in arterial spin labeling (ASL) that is eliminated by delays alternating with nutation for tailored excitation (DANTE) pulse using T₂ -relaxometry, and to demonstrate the feasibility of arterial cerebral blood volume (CBV_a ) imaging using DANTE-ASL in combination with a simplified two-compartment model.

METHOD

The DANTE and T₂ -preparation modules were combined into a single ASL sequence. T₂ values under the application of DANTE were determined to evaluate changes in T₂ , along with the post-labeling delay (PLD) and the relationship between transit time without DANTE (TT_noVS ) and T₂ . The reference tissue T₂ (T_{2_ref} ) was also obtained. Subsequently, the DANTE module was embedded into the Hadamard-encoded ASL. Cerebral blood flow (CBF) and CBV_a were computed using two Hadamard-encoding datasets (with and without DANTE) in a rest and breath-holding (BH) task.

RESULTS

While T₂ without DANTE (T_{2_noVS} ) decreased as the PLD increased, T₂ with DANTE (T_{2_DANTE} ) was equivalent to T_{2_ref} and did not change with the PLD. Although there was a significant positive correlation between TT_noVS and T_{2_noVS} with short PLD, T_{2_DANTE} was not correlated with TT_noVS nor PLD. Baseline CBV_a values obtained at rest were 0.64 ± 0.12, 0.64 ± 0.11, and 0.58 ± 0.15 mL/100 g for anterior, middle, and posterior cerebral arteries, respectively. Significant CBF and CBV_a elevations were observed in the BH task.

CONCLUSION

Microvascular compartment signals were eliminated from the total ASL signals by DANTE. CBV_a can be measured using Hadamard-encoded DANTE-ASL in combination with a simplified two-compartment model.

Functional Activities Detected in the Olfactory Bulb and Associated Olfactory Regions in the Human Brain Using T2-Prepared BOLD Functional MRI at 7T

Olfaction is a fundamental sense that plays a vital role in daily life in humans, and can be altered in neuropsychiatric and neurodegenerative diseases. Blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging (fMRI) using conventional echo-planar-imaging (EPI) based sequences can be challenging in brain regions important for olfactory processing, such as the olfactory bulb (OB) and orbitofrontal cortex, mainly due to the signal dropout and distortion artifacts caused by large susceptibility effects from the sinonasal cavity and temporal bone. To date, few studies have demonstrated successful fMRI in the OB in humans. T2-prepared (T2prep) BOLD fMRI is an alternative approach developed especially for performing fMRI in regions affected by large susceptibility artifacts. The purpose of this technical study is to evaluate T2prep BOLD fMRI for olfactory functional experiments in humans. Olfactory fMRI scans were performed on 7T in 14 healthy participants. T2prep BOLD showed greater sensitivity than GRE EPI BOLD in the OB, orbitofrontal cortex and the temporal pole. Functional activation was detected using T2prep BOLD in the OB and associated olfactory regions. Habituation effects and a bi-phasic pattern of fMRI signal changes during olfactory stimulation were observed in all regions. Both positively and negatively activated regions were observed during olfactory stimulation. These signal characteristics are generally consistent with literature and showed a good intra-subject reproducibility comparable to previous human BOLD fMRI studies. In conclusion, the methodology demonstrated in this study holds promise for future olfactory fMRI studies in the OB and other brain regions that suffer from large susceptibility artifacts.

Validating layer-specific VASO across species

Cerebral blood volume (CBV) has been shown to be a robust and important physiological parameter for quantitative interpretation of functional (f)MRI, capable of delivering highly localized mapping of neural activity. Indeed, with recent advances in ultra-high-field (≥7T) MRI hardware and associated sequence libraries, it has become possible to capture non-invasive CBV weighted fMRI signals across cortical layers. One of the most widely used approaches to achieve this (in humans) is through vascular-space-occupancy (VASO) fMRI. Unfortunately, the exact contrast mechanisms of layer-dependent VASO fMRI have not been validated for human fMRI and thus interpretation of such data is confounded. Here we validate the signal source of layer-dependent SS-SI VASO fMRI using multi-modal imaging in a rat model in response to neuronal activation (somatosensory cortex) and respiratory challenge (hypercapnia). In particular VASO derived CBV measures are directly compared to concurrent measures of total haemoglobin changes from high resolution intrinsic optical imaging spectroscopy (OIS). Quantified cortical layer profiling is demonstrated to be in agreement between VASO and contrast enhanced fMRI (using monocrystalline iron oxide nanoparticles, MION). Responses show high spatial localisation to layers of cortical processing independent of confounding large draining veins which can hamper BOLD fMRI studies, (depending on slice positioning). Thus, a cross species comparison is enabled using VASO as a common measure. We find increased VASO based CBV reactivity (3.1 ± 1.2 fold increase) in humans compared to rats. Together, our findings confirm that the VASO contrast is indeed a reliable estimate of layer-specific CBV changes. This validation study increases the neuronal interpretability of human layer-dependent VASO fMRI as an appropriate method in neuroscience application studies, in which the presence of large draining intracortical and pial veins limits neuroscientific inference with BOLD fMRI.

Impaired response of cerebral oxygen metabolism to visual stimulation in Huntington's disease

Huntington's disease (HD) is a neurodegenerative disease caused by a CAG triplet repeat expansion in the Huntingtin gene. Metabolic and microvascular abnormalities in the brain may contribute to early physiological changes that subserve the functional impairments in HD. This study is intended to investigate potential abnormality in dynamic changes in cerebral blood volume (CBV) and cerebral blood flow (CBF), and cerebral metabolic rate of oxygen (CMRO2) in the brain in response to functional stimulation in premanifest and early manifest HD patients. A recently developed 3-D-TRiple-acquisition-after-Inversion-Preparation magnetic resonance imaging (MRI) approach was used to measure dynamic responses in CBV, CBF, and CMRO2 during visual stimulation in one single MRI scan. Experiments were conducted in 23 HD patients and 16 healthy controls. Decreased occipital cortex CMRO2 responses were observed in premanifest and early manifest HD patients compared to controls (P < 0.001), correlating with the CAG-Age Product scores in these patients (R2 = 0.4, P = 0.001). The results suggest the potential value of this reduced CMRO2 response during visual stimulation as a biomarker for HD and may illuminate the role of metabolic alterations in the pathophysiology of HD.

Cerebrovascular Reactivity Mapping Without Gas Challenges: A Methodological Guide

Cerebrovascular reactivity (CVR) is defined as the ability of vessels to alter their caliber in response to vasoactive factors, by means of dilating or constricting, in order to increase or decrease regional cerebral blood flow (CBF). Importantly, CVR may provide a sensitive biomarker for pathologies where vasculature is compromised. Furthermore, the spatiotemporal dynamics of CVR observed in healthy subjects, reflecting regional differences in cerebral vascular tone and response, may also be important in functional MRI studies based on neurovascular coupling mechanisms. Assessment of CVR is usually based on the use of a vasoactive stimulus combined with a CBF measurement technique. Although transcranial Doppler ultrasound has been frequently used to obtain global flow velocity measurements, MRI techniques are being increasingly employed for obtaining CBF maps. For the vasoactive stimulus, vasodilatory hypercapnia is usually induced through the manipulation of respiratory gases, including the inhalation of increased concentrations of carbon dioxide. However, most of these methods require an additional apparatus and complex setups, which not only may not be well-tolerated by some populations but are also not widely available. For these reasons, strategies based on voluntary breathing fluctuations without the need for external gas challenges have been proposed. These include the task-based methodologies of breath holding and paced deep breathing, as well as a new generation of methods based on spontaneous breathing fluctuations during resting-state. Despite the multitude of alternatives to gas challenges, existing literature lacks definitive conclusions regarding the best practices for the vasoactive modulation and associated analysis protocols. In this work, we perform an extensive review of CVR mapping techniques based on MRI and CO2 variations without gas challenges, focusing on the methodological aspects of the breathing protocols and corresponding data analysis. Finally, we outline a set of practical guidelines based on generally accepted practices and available data, extending previous reports and encouraging the wider application of CVR mapping methodologies in both clinical and academic MRI settings.

Measurement of CMRO2 and its relationship with CBF in hypoxia with an extended calibrated BOLD method

Cerebral blood flow (CBF) and cerebral metabolic rate of oxygen (CMRO₂) are physiological parameters that not only reflect brain health and disease but also jointly contribute to blood oxygen level-dependent (BOLD) signals. Nevertheless, unsolved issues remain concerning the CBF-CMRO₂ relationship in the working brain under various oxygen conditions. In particular, the CMRO₂ responses to functional tasks in hypoxia are less studied. We extended the calibrated BOLD model to incorporate CMRO₂ measurements in hypoxia. The extended model, which was cross-validated with a multicompartment BOLD model, considers the influences of the reduced arterial saturation level and increased baseline cerebral blood volume (CBV) and deoxyhemoglobin concentration on the changes of BOLD signals in hypoxia. By implementing a pulse sequence to simultaneously acquire the CBV-, CBF- and BOLD-weighted signals, we investigated the effects of mild hypoxia on the CBF and CMRO₂ responses to graded visual stimuli. Compared with normoxia, mild hypoxia caused significant alterations in both the amplitude and the trend of the CMRO₂ responses but did not impact the corresponding CBF responses. Our observations suggested that the flow-metabolism coupling strategies in the brain during mild hypoxia were different from those during normoxia.

Modeling of vascular space occupancy and BOLD functional MRI from first principles using real microvascular angiograms

PURPOSE

The vascular space occupancy (VASO) is a functional MRI technique for probing cerebral blood volume changes noninvasively, including during neuronal activation in humans. An important consideration when implementing VASO is the BOLD effect in the signal. Assessing the physical origin of this BOLD contamination and the capabilities of correction methods could improve the quantification of cerebral blood volume changes with VASO.

METHODS

Given the heterogeneity of cerebral microvascular architecture, the vascular geometry within an MRI voxel can influence both BOLD and VASO signals. To investigate this effect, 3D high-resolution images of mouse cerebral vasculature measured with two-photon microscopy were used to model BOLD and VASO signals from first principles using Monte Carlo diffusion of water protons. Quantitative plots of VASO together with intravascular and extravascular BOLD signals as a function of TE at B₀ fields 1.5 T to 14 T were obtained.

RESULTS

The BOLD contamination of the VASO response was on the order of 50% for gradient echo and 5% for spin echo at 7 T and TE = 6 ms and significantly increased with TE and B₀ . Two currently used correction schemes were shown to account for most of this contamination and recover accurate relative signal changes, with optimal correction obtained using TEs as short as possible.

CONCLUSION

These results may provide useful information for optimizing sequence parameters in VASO and BOLD functional MRI, leading the way to a wider application of these techniques in healthy and diseased brain.

Quantitative theory for the transverse relaxation time of blood water

An integrative model is proposed to describe the dependence of the transverse relaxation rate of blood water protons (R_2blood = 1/T_2blood ) on hematocrit fraction and oxygenation fraction (Y). This unified model takes into account (a) the diamagnetic effects of albumin, hemoglobin and the cell membrane; (b) the paramagnetic effect of hemoglobin; (c) the effect of compartmental exchange between plasma and erythrocytes under both fast and slow exchange conditions that vary depending on field strength and compartmental relaxation rates and (d) the effect of diffusion through field gradients near the erythrocyte membrane. To validate the model, whole-blood and lysed-blood R₂ data acquired previously using Carr-Purcell-Meiboom-Gill measurements as a function of inter-echo spacing τ_cp at magnetic fields of 3.0, 7.0, 9.4 and 11.7 T were fitted to determine the lifetimes (field-independent physiological constants) for water diffusion and exchange, as well as several physical constants, some of which are field-independent (magnetic susceptibilities) and some are field-dependent (relaxation rates for water protons in solutions of albumin and oxygenated and deoxygenated hemoglobin, ie, blood plasma and erythrocytes, respectively). This combined exchange-diffusion model allowed excellent fitting of the curve of the τ_cp -dependent relaxation rate dispersion at all four fields using a single average erythrocyte water lifetime, τ_ery = 9.1 ± 1.4 ms, and an averaged diffusional correlation time, τ_D = 3.15 ± 0.43 ms. Using this model and the determined physiological time constants and relaxation parameters, blood T₂ values published by multiple groups based on measurements at magnetic field strengths of 1.5 T and higher could be predicted correctly within error. Establishment of this theory is a fundamental step for quantitative modeling of the BOLD effect underlying functional MRI.

Correcting Task fMRI Signals for Variability in Baseline CBF Improves BOLD-Behavior Relationships: A Feasibility Study in an Aging Model

Blood Oxygen Level Dependent (BOLD) functional MRI is a complex neurovascular signal whose magnitude depends on baseline physiological factors such as cerebral blood flow (CBF). Because baseline CBF varies across the brain and is altered with aging, the interpretation of stand-alone aging-related BOLD changes can be misleading. The primary objective of this study was to develop a methodology that combines task fMRI and arterial spin labeling (ASL) techniques to sensitize task-induced BOLD activity by covarying out the baseline physiology (i.e., CBF) in an aging model. We recruited 11 younger and 13 older healthy participants who underwent ASL and an overt language fMRI task (semantic category member generation). We measured in-scanner language performance to investigate the effect of BOLD sensitization on BOLD-behavior relationships. The results demonstrate that our correction approach is effective at enhancing the specificity and sensitivity of the BOLD signal in both groups. In addition, the correction strengthens the statistical association between task BOLD activity and behavioral performance. Although CBF has inherent age dependence, our results show that retaining the age factor within CBF aides in greater sensitization of task fMRI signals. From a cognitive standpoint, compared to young adults, the older participants showed a delayed domain-general language-related task activity possibly due to compromised vessel compliance. Further, assessment of functional evolution of corrected BOLD activity revealed biphasic BOLD dynamics in both groups where BOLD deactivation may reflect greater semantic demand or increased premium on domain general executive functioning in response to task difficulty. Although it was promising to note that the predictability of behavior using the proposed methodology outperforms other methodologies (i.e., no correction and normalization by division), and provides moderate stability and adequate power, further work with a larger cohort and other task designs is necessary to improve the stability of predicting associated behavior. In summary, we recommend correction of task fMRI signals by covarying out baseline CBF especially when comparing groups with different neurovascular properties. Given that ASL and BOLD fMRI are well established and widely employed techniques, our proposed multi-modal methodology can be readily implemented into data processing pipelines to obtain more accurate BOLD activation maps.

Changes in Hemodynamic Response Function Resulting From Chronic Alcohol Consumption

BACKGROUND

Functional MRI (fMRI) task-related analyses rely on an estimate of the brain's hemodynamic response function (HRF) to model the brain's response to events. Although changes in the HRF have been found after acute alcohol administration, the effects of heavy chronic alcohol consumption on the HRF have not been explored, and the potential benefits or pitfalls of estimating each individual's HRF on fMRI analyses of chronic alcohol use disorder (AUD) are not known.

METHODS

Participants with AUD and controls (CTL) received structural, functional, and vascular scans. During fMRI, participants were cued to tap their fingers, and averaged responses were extracted from the motor cortex. Curve fitting on these HRFs modeled them as a difference between 2 gamma distributions, and the temporal occurrence of the main peak and undershoot of the HRF was computed from the mean of the first and second gamma distributions, respectively.

RESULTS

ANOVA and regression analyses found that the timing of the HRF undershoot increased significantly as a function of total lifetime drinking. Although gray matter volume in the motor cortex decreased with lifetime drinking, this was not sufficient to explain undershoot timing shifts, and vascular factors measured in the motor cortex did not differ among groups. Comparison of random-effects analyses using custom-fitted and canonical HRFs for CTL and AUD groups showed better results throughout the brain for custom-fitted versus canonical HRFs for CTL subjects. For AUD subjects, the same was true except for the basal ganglia.

CONCLUSIONS

These findings suggest that excessive alcohol consumption is associated with changes in the HRF undershoot. HRF changes could provide a possible biomarker for the effects of lifetime drinking on brain function. Changes in HRF topography affect fMRI activation measures, and subject-specific HRFs generally improve fMRI activation results.

Magnetic Resonance Imaging in Aneurysmal Subarachnoid Hemorrhage: Current Evidence and Future Directions

BACKGROUND

Aneurysmal subarachnoid hemorrhage (aSAH) is associated with an unacceptably high mortality and chronic disability in survivors, underscoring a need to validate new approaches for treatment and prognosis. The use of advanced imaging, magnetic resonance imaging (MRI) in particular, could help address this gap given its versatile capacity to quantitatively evaluate and map changes in brain anatomy, physiology and functional activation. Yet there is uncertainty about the real value of brain MRI in the clinical setting of aSAH.

METHODS

In this review, we discuss current and emerging MRI research in aSAH. PubMed was searched from inception to June 2017, and additional studies were then chosen on the basis of relevance to the topics covered in this review.

RESULTS

Available studies suggest that brain MRI is a feasible, safe, and valuable testing modality. MRI detects brain abnormalities associated with neurologic examination, outcomes, and aneurysm treatment and thus has the potential to increase knowledge of aSAH pathophysiology as well as to guide management and outcome prediction. Newer pulse sequences have the potential to reveal structural and physiological changes that could also improve management of aSAH.

CONCLUSION

Research is needed to confirm the value of MRI-based biomarkers in clinical practice and as endpoints in clinical trials, with the goal of improving outcome for patients with aSAH.

More than BOLD: Dual-spin populations create functional contrast

PURPOSE

Functional MRI contrast has generally been associated with changes in transverse relaxivity caused by blood oxygen concentration, the so-called blood oxygen level dependent contrast. However, this interpretation of fMRI contrast has been called into question by several recent experiments at high spatial resolution. Experiments were conducted to examine contrast dependencies that cannot be explained only by differences in relaxivity in a single-spin population.

METHODS

Measurements of functional signal and contrast were obtained in human early visual cortex during a high-contrast visual stimulation over a large range of TEs and for several flip angles. Small voxels (1.5 mm) were used to restrict the measurements to cortical gray matter in early visual areas identified using retinotopic mapping procedures.

RESULTS

Measurements were consistent with models that include 2 spin populations. The dominant population has a relatively short transverse lifetime that is strongly modulated by activation. However, functional contrast is also affected by volume changes between this short-lived population and the longer-lived population.

CONCLUSION

Some of the previously observed "nonclassical" behaviors of functional contrast can be explained by these interacting dual-spin populations.

Cerebral blood volume changes during the BOLD post-stimulus undershoot measured with a combined normoxia/hyperoxia method

Cerebral blood flow (CBF) and blood oxygenation level dependent (BOLD) signal measurements make it possible to estimate steady-state changes in the cerebral metabolic rate of oxygen (CMRO2) with a calibrated BOLD method. However, extending this approach to measure the dynamics of CMRO2 requires an additional assumption: that deoxygenated cerebral blood volume (CBVdHb) follows CBF in a predictable way. A test-case for this assumption is the BOLD post-stimulus undershoot, for which one proposed explanation is a strong uncoupling of flow and blood volume with an elevated level of CBVdHb during the post-stimulus period compared to baseline due to slow blood volume recovery (Balloon Model). A challenge in testing this model is that CBVdHb differs from total blood volume, which can be measured with other techniques. In this study, the basic hypothesis of elevated CBVdHb during the undershoot was tested, based on the idea that the BOLD signal change when a subject switches from breathing a normoxic gas to breathing a hyperoxic gas is proportional to the absolute CBVdHb. In 19 subjects (8F), dual-echo BOLD responses were measured in primary visual cortex during a flickering radial checkerboard stimulus in normoxia, and the identical experiment was repeated in hyperoxia (50% O2/balance N2). The BOLD signal differences between normoxia and hyperoxia for the pre-stimulus baseline, stimulus, and post-stimulus periods were compared using an equivalent BOLD signal calculated from measured R2* changes to eliminate signal drifts. Relative to the pre-stimulus baseline, the average BOLD signal change from normoxia to hyperoxia was negative during the undershoot period (p = 0.0251), consistent with a reduction of CBVdHb and contrary to the prediction of the Balloon Model. Based on these results, the BOLD post-stimulus undershoot does not represent a case of strong uncoupling of CBVdHb and CBF, supporting the extension of current calibrated BOLD methods to estimate the dynamics of CMRO2.

The effects of capillary transit time heterogeneity on the BOLD signal

Neurovascular coupling mechanisms give rise to vasodilation and functional hyperemia upon neural activation, thereby altering blood oxygenation. This blood oxygenation level dependent (BOLD) contrast allows studies of activation patterns in the working human brain by functional MRI (fMRI). The BOLD-weighted fMRI signal shows characteristic transients in relation to functional activation, such as the so-called initial dip, overshoot, and post-stimulus undershoot. These transients are modulated by other physiological stimuli and in disease, but the underlying physiological mechanisms remain incompletely understood. Capillary transit time heterogeneity (CTH) has been shown to affect oxygen extraction, and hence blood oxygenation. Here, we examine how recently reported redistributions of capillary blood flow during functional activation would be expected to affect BOLD signal transients. We developed a three-compartment (hemoglobin, plasma, and tissue) model to predict the BOLD signal, incorporating the effects of dynamic changes in CTH. Our model predicts that the BOLD signal represents the superposition of a positive component resulting from increases in cerebral blood flow (CBF), and a negative component, resulting from elevated tissue metabolism and homogenization of capillary flows (reduced CTH). The model reproduces salient features of BOLD signal dynamics under conditions such as hypercapnia, hyperoxia, and caffeine intake, where both brain physiology and BOLD characteristics are altered. Neuroglial signaling and metabolism could affect CBF and capillary flow patterns differently. Further studies of neurovascular and neuro-capillary coupling mechanisms may help us relate BOLD signals to the firing of certain neuronal populations based on their respective BOLD "fingerprints."

Is There a Nonadditive Interaction Between Spontaneous and Evoked Activity? Phase-Dependence and Its Relation to the Temporal Structure of Scale-Free Brain Activity

The aim of our study was to use functional magnetic resonance imaging to investigate how spontaneous activity interacts with evoked activity, as well as how the temporal structure of spontaneous activity, that is, long-range temporal correlations, relate to this interaction. Using an extremely sparse event-related design (intertrial intervals: 52-60 s), a novel blood oxygen level-dependent signal correction approach (accounting for spontaneous fluctuations using pseudotrials) and phase analysis, we provided direct evidence for a nonadditive interaction between spontaneous and evoked activity. We demonstrated the discrepancy between the present and previous observations on why a linear superposition between spontaneous and evoked activity can be seen by using co-occurring signals from homologous brain regions. Importantly, we further demonstrated that the nonadditive interaction can be characterized by phase-dependent effects of spontaneous activity, which is closely related to the degree of long-range temporal correlations in spontaneous activity as indexed by both power-law exponent and phase-amplitude coupling. Our findings not only contribute to the understanding of spontaneous brain activity and its scale-free properties, but also bear important implications for our understanding of neural activity in general.

MRI techniques to measure arterial and venous cerebral blood volume

The measurement of cerebral blood volume (CBV) has been the topic of numerous neuroimaging studies. To date, however, most in vivo imaging approaches can only measure CBV summed over all types of blood vessels, including arterial, capillary and venous vessels in the microvasculature (i.e. total CBV or CBVtot). As different types of blood vessels have intrinsically different anatomy, function and physiology, the ability to quantify CBV in different segments of the microvascular tree may furnish information that is not obtainable from CBVtot, and may provide a more sensitive and specific measure for the underlying physiology. This review attempts to summarize major efforts in the development of MRI techniques to measure arterial (CBVa) and venous CBV (CBVv) separately. Advantages and disadvantages of each type of method are discussed. Applications of some of the methods in the investigation of flow-volume coupling in healthy brains, and in the detection of pathophysiological abnormalities in brain diseases such as arterial steno-occlusive disease, brain tumors, schizophrenia, Huntington's disease, Alzheimer's disease, and hypertension are demonstrated. We believe that the continual development of MRI approaches for the measurement of compartment-specific CBV will likely provide essential imaging tools for the advancement and refinement of our knowledge on the exquisite details of the microvasculature in healthy and diseased brains.

Determining Excitatory and Inhibitory Neuronal Activity from Multimodal fMRI Data Using a Generative Hemodynamic Model

Hemodynamic responses, in general, and the blood oxygenation level-dependent (BOLD) fMRI signal, in particular, provide an indirect measure of neuronal activity. There is strong evidence that the BOLD response correlates well with post-synaptic changes, induced by changes in the excitatory and inhibitory (E-I) balance between active neuronal populations. Typical BOLD responses exhibit transients, such as the early-overshoot and post-stimulus undershoot, that can be linked to transients in neuronal activity, but they can also result from vascular uncoupling between cerebral blood flow (CBF) and venous cerebral blood volume (venous CBV). Recently, we have proposed a novel generative hemodynamic model of the BOLD signal within the dynamic causal modeling framework, inspired by physiological observations, called P-DCM (Havlicek et al., 2015). We demonstrated the generative model's ability to more accurately model commonly observed neuronal and vascular transients in single regions but also effective connectivity between multiple brain areas (Havlicek et al., 2017b). In this paper, we additionally demonstrate the versatility of the generative model to jointly explain dynamic relationships between neuronal and hemodynamic physiological variables underlying the BOLD signal using multi-modal data. For this purpose, we utilized three distinct data-sets of experimentally induced responses in the primary visual areas measured in human, cat, and monkey brain, respectively: (1) CBF and BOLD responses; (2) CBF, total CBV, and BOLD responses (Jin and Kim, 2008); and (3) positive and negative neuronal and BOLD responses (Shmuel et al., 2006). By fitting the generative model to the three multi-modal experimental data-sets, we showed that the presence or absence of dynamic features in the BOLD signal is not an unambiguous indication of presence or absence of those features on the neuronal level. Nevertheless, the generative model that takes into account the dynamics of the physiological mechanisms underlying the BOLD response allowed dissociating neuronal from vascular transients and deducing excitatory and inhibitory neuronal activity time-courses from BOLD data alone and from multi-modal data.

Accounting for the role of hematocrit in between-subject variations of MRI-derived baseline cerebral hemodynamic parameters and functional BOLD responses

Baseline hematocrit fraction (Hct) is a determinant for baseline cerebral blood flow (CBF) and between-subject variation of Hct thus causes variation in task-based BOLD fMRI signal changes. We first verified in healthy volunteers (n = 12) that Hct values can be derived reliably from venous blood T₁ values by comparison with the conventional lab test. Together with CBF measured using phase-contrast MRI, this noninvasive estimation of Hct, instead of using a population-averaged Hct value, enabled more individual determination of oxygen delivery (DO₂ ), oxygen extraction fraction (OEF), and cerebral metabolic rate of oxygen (CMRO₂ ). The inverse correlation of CBF and Hct explained about 80% of between-subject variation of CBF in this relatively uniform cohort of subjects, as expected based on the regulation of DO₂ to maintain constant CMRO₂ . Furthermore, we compared the relationships of visual task-evoked BOLD response with Hct and CBF. We showed that Hct and CBF contributed 22%-33% of variance in BOLD signal and removing the positive correlation with Hct and negative correlation with CBF allowed normalization of BOLD signal with 16%-22% lower variability. The results of this study suggest that adjustment for Hct effects is useful for studies of MRI perfusion and BOLD fMRI. Hum Brain Mapp 39:344-353, 2018. © 2017 Wiley Periodicals, Inc.

Post-stimulus fMRI and EEG responses: Evidence for a neuronal origin hypothesised to be inhibitory

Post-stimulus undershoots, negative responses following cessation of stimulation, are widely observed in functional magnetic resonance (fMRI) blood oxygenation level dependent (BOLD) data. However, the debate surrounding whether the origin of this response phase is neuronal or vascular, and whether it provides functionally relevant information, that is additional to what is contained in the primary response, means that undershoots are widely overlooked. We simultaneously recorded electroencephalography (EEG), BOLD and cerebral blood-flow (CBF) [obtained from arterial spin labelled (ASL) fMRI] fMRI responses to hemifield checkerboard stimulation to test the potential neural origin of the fMRI post-stimulus undershoot. The post-stimulus BOLD and CBF signal amplitudes in both contralateral and ipsilateral visual cortex depended on the post-stimulus power of the occipital 8-13Hz (alpha) EEG neuronal activity, such that trials with highest EEG power showed largest fMRI undershoots in contralateral visual cortex. This correlation in post-stimulus EEG-fMRI responses was not predicted by the primary response amplitude. In the contralateral visual cortex we observed a decrease in both cerebral rate of oxygen metabolism (CMRO₂) and CBF during the post-stimulus phase. In addition, the coupling ratio (n) between CMRO₂ and CBF was significantly lower during the positive contralateral primary response phase compared with the post-stimulus phase and we propose that this reflects an altered balance of excitatory and inhibitory neuronal activity. Together our data provide strong evidence that the post-stimulus phase of the BOLD response has a neural origin which reflects, at least partially, an uncoupling of the neuronal responses driving the primary and post-stimulus responses, explaining the uncoupling of the signals measured in the two response phases. We suggest our results are consistent with inhibitory processes driving the post-stimulus EEG and fMRI responses. We therefore propose that new methods are required to model the post-stimulus and primary responses independently, enabling separate investigation of response phases in cognitive function and neurological disease.

Foundations of layer-specific fMRI and investigations of neurophysiological activity in the laminarized neocortex and olfactory bulb of animal models

Laminar organization of neuronal circuits is a recurring feature of how the brain processes information. For instance, different layers compartmentalize different cell types, synaptic activities, and have unique intrinsic and extrinsic connections that serve as units for specialized signal processing. Functional MRI is an invaluable tool to investigate laminar processing in the in vivo human brain, but it measures neuronal activity indirectly by way of the hemodynamic response. Therefore, the accuracy of high-resolution laminar fMRI depends on how precisely it can measure localized microvascular changes nearest to the site of evoked activity. To determine the specificity of fMRI responses to the true neurophysiological responses across layers, the flexibility to invasive procedures in animal models has been necessary. In this review, we will examine different fMRI contrasts and their appropriate uses for layer-specific fMRI, and how localized laminar processing was examined in the neocortex and olfactory bulb. Through collective efforts, it was determined that microvessels, including capillaries, are regulated within single layers and that several endogenous and contrast-enhanced fMRI contrast mechanisms can separate these neural-specific vascular changes from the nonspecific, especially cerebral blood volume-weighted fMRI with intravenous contrast agent injection. We will also propose some open questions that are relevant for the successful implementation of layer-specific fMRI and its potential future directions to study laminar processing when combined with optogenetics.

Cerebrovascular reactivity measured by functional magnetic resonance imaging during breath-hold challenge: A systematic review

Cerebrovascular reactivity (CVR) is the cerebral hemodynamic response to a vasoactive substance. Breath-hold (BH) induced CVR has the advantage of being non-invasive and easy to implement during magnetic resonance imaging (MRI). We systematically reviewed the literature regarding MRI measurement of BH induced CVR. The literature was searched using MEDLINE with the search terms breath-hold; and MRI or cerebrovascular reactivity. The search yielded 2244 results and 54 articles were included. Between-group comparisons have found that CVR was higher among healthy controls than patients with various pathologies (e.g. sleep apnea, diabetes, hypertension etc.). However, counter-intuitive findings have also been reported, including higher CVR among smokers, sedentary individuals, and patients with schizophrenia vs.

CONTROLS

Methodological studies have highlighted important measurement characteristics (e.g. normalizing signal to end-tidal CO₂), and comparisons of BH induced CVR to non-BH methods. Future studies are warranted to address questions about group differences, treatment response, disease progression, and other salient clinical themes. Standardization of CVR and BH designs is needed to fully exploit the potential of this practical non-invasive method.

Noise concerns and post-processing procedures in cerebral blood flow (CBF) and cerebral blood volume (CBV) functional magnetic resonance imaging

Functional neuroimaging with blood oxygenation level-dependent (BOLD) contrast has emerged as the most popular method for evaluating qualitative changes in brain function in humans. At typical human field strengths (1.5-3.0T), BOLD contrast provides a measure of changes in transverse water relaxation rates in and around capillary and venous blood, and as such provides only a surrogate marker of brain function that depends on dynamic changes in hemodynamics (e.g., cerebral blood flow and volume) and metabolism (e.g., oxygen extraction fraction and the cerebral metabolic rate of oxygen consumption). Alternative functional neuroimaging methods that are specifically sensitive to these constituents of the BOLD signal are being developed and applied in a growing number of clinical and neuroscience applications of quantitative cerebral physiology. These methods require additional considerations for interpreting and quantifying their contrast responsibly. Here, an overview of two popular methods, arterial spin labeling and vascular space occupancy, is presented specifically in the context of functional neuroimaging. Appropriate post-processing and experimental acquisition strategies are summarized with the motivation of reducing sensitivity to noise and unintended signal sources and improving quantitative accuracy of cerebral hemodynamics.

Hemodynamic signals in fNIRS

Near-infrared spectroscopy (NIRS) was originally designed for clinical monitoring of tissue oxygenation, and it has also been developed into a useful tool in neuroimaging studies, with the so-called functional NIRS (fNIRS). With NIRS, cerebral activation is detected by measuring the cerebral hemoglobin (Hb), where however, the precise correlation between NIRS signal and neural activity remains to be fully understood. This can in part be attributed to the situation that NIRS signals are inherently subject to contamination by signals arising from extracerebral tissue. In recent years, several approaches have been investigated to distinguish between NIRS signals originating in cerebral tissue and signals originating in extracerebral tissue. Selective measurements of cerebral Hb will enable a further evolution of fNIRS. This chapter is divided into six sections: first a summary of the basic theory of NIRS, NIRS signals arising in the activated areas, correlations between NIRS signals and fMRI signals, correlations between NIRS signals and neural activities, and the influence of a variety of extracerebral tissue on NIRS signals and approaches to this issue are reviewed. Finally, future prospects of fNIRS are described.

Layer-Specific fMRI Responses to Excitatory and Inhibitory Neuronal Activities in the Olfactory Bulb

High-resolution functional magnetic resonance imaging (fMRI) detects localized neuronal activity via the hemodynamic response, but it is unclear whether it accurately identifies neuronal activity specific to individual layers. To address this issue, we preferentially evoked neuronal activity in superficial, middle, and deep layers of the rat olfactory bulb: the glomerular layer by odor (5% amyl acetate), the external plexiform layer by electrical stimulation of the lateral olfactory tract (LOT), and the granule cell layer by electrical stimulation of the anterior commissure (AC), respectively. Electrophysiology, laser-Doppler flowmetry of cerebral blood flow (CBF), and blood oxygenation level-dependent (BOLD) and cerebral blood volume-weighted (CBV) fMRI at 9.4 T were performed independently. We found that excitation of inhibitory granule cells by stimulating LOT and AC decreased the spontaneous multi-unit activities of excitatory mitral cells and subsequently increased CBF, CBV, and BOLD signals. Odor stimulation also increased the hemodynamic responses. Furthermore, the greatest CBV fMRI responses were discretely separated into the same layers as the evoked neuronal activities for all three stimuli, whereas BOLD was poorly localized with some exception to the poststimulus undershoot. In addition, the temporal dynamics of the fMRI responses varied depending on the stimulation pathway, even within the same layer. These results indicate that the vasculature is regulated within individual layers and CBV fMRI has a higher fidelity to the evoked neuronal activity compared with BOLD. Our findings are significant for understanding the neuronal origin and spatial specificity of hemodynamic responses, especially for the interpretation of laminar-resolution fMRI.

SIGNIFICANCE STATEMENT

Functional magnetic resonance imaging (fMRI) is a noninvasive, in vivo technique widely used to map function of the entire brain, including deep structures, in animals and humans. However, it measures neuronal activity indirectly by way of the vascular response. It is currently unclear how finely the hemodynamic response is regulated within single cortical layers and whether increased inhibitory neuronal activities affect fMRI signal changes. Both laminar specificity and the neural origins of fMRI are important to interpret functional maps properly, which we investigated by activating discrete rat olfactory bulb circuits.

Hemodynamic response function abnormalities in schizophrenia during a multisensory detection task

Functional magnetic resonance imaging (fMRI) of the blood oxygen level dependent (BOLD) response has commonly been used to investigate the neuropathology underlying cognitive and sensory deficits in patients with schizophrenia (SP) by examining the positive phase of the BOLD response, assuming a fixed shape for the hemodynamic response function (HRF). However, the individual phases (positive and post-stimulus undershoot (PSU)) of the HRF may be differentially affected by a variety of underlying pathologies. The current experiment used a multisensory detection task with a rapid event-related fMRI paradigm to investigate both the positive and PSU phases of the HRF in SP and healthy controls (HC). Behavioral results indicated no significant group differences during task performance. Analyses that examined the shape of the HRF indicated two distinct group differences. First, SP exhibited a reduced and/or prolonged PSU following normal task-related positive BOLD activation in secondary auditory and visual sensory areas relative to HC. Second, SP did not show task-induced deactivation in the anterior node of the default-mode network (aDMN) relative to HC. In contrast, when performing traditional analyses that focus on the positive phase, there were no group differences. Interestingly, the magnitude of the PSU in secondary auditory and visual areas was positively associated with the magnitude of task-induced deactivation within the aDMN, suggesting a possible common neural mechanism underlying both of these abnormalities (failure in neural inhibition). Results are consistent with recent views that separate neural processes underlie the two phases of the HRF and that they are differentially affected in SP. Hum Brain Mapp 37:745-755, 2016. © 2015 Wiley Periodicals, Inc.

Temporal pattern of acoustic imaging noise asymmetrically modulates activation in the auditory cortex

This study investigated the hemisphere-specific effects of the temporal pattern of imaging related acoustic noise on auditory cortex activation. Hemodynamic responses (HDRs) to five temporal patterns of imaging noise corresponding to noise generated by unique combinations of imaging volume and effective repetition time (TR), were obtained using a stroboscopic event-related paradigm with extra-long (≥27.5 s) TR to minimize inter-acquisition effects. In addition to confirmation that fMRI responses in auditory cortex do not behave in a linear manner, temporal patterns of imaging noise were found to modulate both the shape and spatial extent of hemodynamic responses, with classically non-auditory areas exhibiting responses to longer duration noise conditions. Hemispheric analysis revealed the right primary auditory cortex to be more sensitive than the left to the presence of imaging related acoustic noise. Right primary auditory cortex responses were significantly larger during all the conditions. This asymmetry of response to imaging related acoustic noise could lead to different baseline activation levels during acquisition schemes using short TR, inducing an observed asymmetry in the responses to an intended acoustic stimulus through limitations of dynamic range, rather than due to differences in neuronal processing of the stimulus. These results emphasize the importance of accounting for the temporal pattern of the acoustic noise when comparing findings across different fMRI studies, especially those involving acoustic stimulation.

Venous cerebral blood volume increase during voluntary locomotion reflects cardiovascular changes

Understanding how changes in the cardiovascular system contribute to cerebral blood flow (CBF) and volume (CBV) increases is critical for interpreting hemodynamic signals. Here we investigated how systemic cardiovascular changes affect the cortical hemodynamic response during voluntary locomotion. In the mouse, voluntary locomotion drives an increase in cortical CBF and arterial CBV that is localized to the forelimb/hindlimb representation in the somatosensory cortex, as well as a diffuse venous CBV increase. To determine if the heart rate increases that accompany locomotion contribute to locomotion-induced CBV and CBF increases, we occluded heart rate increases with the muscarinic cholinergic receptor antagonist glycopyrrolate, and reduced heart rate with the β1-adrenergic receptor antagonist atenolol. We quantified the effects of these cardiovascular manipulations on CBV and CBF dynamics by comparing the hemodynamic response functions (HRF) to locomotion across these conditions. Neither the CBF HRF nor the arterial component of the CBV HRF was significantly affected by pharmacological disruption of the heart rate. In contrast, the amplitude and spatial extent of the venous component of the CBV HRF were decreased by atenolol. These results suggest that the increase in venous CBV during locomotion was partially driven by peripheral cardiovascular changes, whereas CBF and arterial CBV increases associated with locomotion reflect central processes.

Measurement of parenchymal extravascular R2* and tissue oxygen extraction fraction using multi-echo vascular space occupancy MRI at 7 T

Parenchymal extravascular R2* is an important parameter for quantitative blood oxygenation level-dependent (BOLD) studies. Total and intravascular R2* values and changes in R2* values during functional stimulations have been reported in a number of studies. The purpose of this study was to measure absolute extravascular R2* values in human visual cortex and to estimate the intra- and extravascular contributions to the BOLD effect at 7 T. Vascular space occupancy (VASO) MRI was employed to separate out the extravascular tissue signal. Multi-echo VASO and BOLD functional MRI (fMRI) with visual stimulation were performed at 7 T for R2* measurement at a spatial resolution of 2.5 × 2.5 × 2.5 mm(3) in healthy volunteers (n = 6). The ratio of changes in extravascular and total R2* (ΔR2*) was used to estimate the extravascular fraction of the BOLD effect. Extravascular R2* values were found to be 44.66 ± 1.55 and 43.38 ± 1.51 s(-1) (mean ± standard error of the mean, n = 6) at rest and activation, respectively, in human visual cortex at 7 T. The extravascular BOLD fraction was estimated to be 91 ± 3%. The parenchymal oxygen extraction fraction (OEF) during activation was estimated to be 0.24 ± 0.01 based on the R2* measurements, indicating an approximately 37% decrease compared with OEF at rest.

Three-dimensional acquisition of cerebral blood volume and flow responses during functional stimulation in a single scan

In addition to the BOLD scan, quantitative functional MRI studies require measurement of both cerebral blood volume (CBV) and flow (CBF) dynamics. The ability to detect CBV and CBF responses in a single additional scan would shorten the total scan time and reduce temporal variations. Several approaches for simultaneous CBV and CBF measurement during functional MRI experiments have been proposed in two-dimensional (2D) mode covering one to three slices in one repetition time (TR). Here, we extended the principles from previous work and present a three-dimensional (3D) whole-brain MRI approach that combines the vascular-space-occupancy (VASO) and flow-sensitive alternating inversion recovery (FAIR) arterial spin labeling (ASL) techniques, allowing the measurement of CBV and CBF dynamics, respectively, in a single scan. 3D acquisitions are complicated for such a scan combination as the time to null blood signal during a steady state needs to be known. We estimated this using Bloch simulations and demonstrate that the resulting 3D acquisition can detect activation patterns and relative signal changes of quality comparable to that of the original separate scans. The same was found for temporal signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR). This approach provides improved acquisition efficiency when both CBV and CBF responses need to be monitored during a functional task.

Cortical depth dependence of the BOLD initial dip and poststimulus undershoot in human visual cortex at 7 Tesla

PURPOSE

Owing to variability in vascular dynamics across cerebral cortex, blood-oxygenation-level-dependent (BOLD) spatial and temporal characteristics should vary as a function of cortical-depth. Here, the positive response, initial dip (ID), and post-stimulus undershoot (PSU) of the BOLD response in human visual cortex are investigated as a function of cortical depth and stimulus duration at 7 Tesla (T).

METHODS

Gradient-echo echo-planar-imaging BOLD fMRI with high spatial and temporal resolution was performed in 7 healthy volunteers and measurements of the ID, PSU, and positive BOLD response were made as a function of cortical depth and stimulus duration (0.5-8 s). Exploratory analyses were applied to understand whether functional mapping could be achieved using the ID, rather than positive, BOLD signal characteristics

RESULTS

The ID was largest in outer cortical layers, consistent with previously reported upstream propagation of vasodilation along the diving arterioles in animals. The positive BOLD signal and PSU showed different relationships across the cortical depth with respect to stimulus duration.

CONCLUSION

The ID and PSU were measured in humans at 7T and exhibited similar trends to those recently reported in animals. Furthermore, while evidence is provided for the ID being a potentially useful feature for better understanding BOLD signal dynamics, such as laminar neurovascular coupling, functional mapping based on the ID is extremely difficult.

Temporal profile of fronto-striatal-limbic activity during implicit decisions in drug dependence

BACKGROUND

Substance dependence is associated with impaired decision-making and altered fronto-striatal-limbic activity. Both greater and lesser brain activity have been reported in drug users compared to controls during decision-making. Inconsistent results might be explained by group differences in the temporal profile of the functional magnetic resonance imaging (fMRI) response. While most previous studies model a canonical hemodynamic response, a finite impulse response (FIR) model measures fMRI signal at discrete time points without assuming a temporal profile. We compared brain activity during decision-making and feedback in substance users and controls using two models: a canonical hemodynamic response function (HRF) and a FIR model.

METHODS

37 substance-dependent individuals (SDI) and 43 controls performed event-related decision-making during fMRI scanning. Brain activity was compared across group using canonical HRF and FIR models.

RESULTS

Compared to controls, SDI were impaired at decision-making. The canonical HRF model showed that SDI had significantly greater fronto-striatal-limbic activity during decisions and less activity during feedback than controls. The FIR model confirmed greater activity in SDI during decisions. However, lower activity in SDI during feedback corresponded to a lower post-stimulus undershoot of the hemodynamic response.

CONCLUSIONS

Greater activity in fronto-striatal-limbic pathways in SDI compared to controls is consistent with prior work, further supporting the hypothesis that abnormalities in these circuits underlie impaired decision-making. We demonstrate for the first time using FIR analysis that lower activity during feedback may simply reflect the tail end of the hemodynamic response to decision, the post-stimulus undershoot, rather than an actual difference in feedback response.

Time-resolved MRI oximetry for quantifying CMRO(2) and vascular reactivity

This brief review of magnetic resonance susceptometry summarizes the methods conceived in the authors' laboratory during the past several years. This article shows how venous oxygen saturation is quantified in large draining veins by field mapping and how this information, in concert with simultaneous measurement of cerebral blood flow, yields cerebral metabolic rate of oxygen, the brain's rate of oxygen consumption. The accuracy of this model-based approach in which the blood vessel is approximated as a long, straight cylinder, for which an analytical solution for the induced field exists, is discussed. It is shown that the approach is remarkably robust, allowing for time-resolved quantification of whole-brain metabolism at rest and in response to stimuli, thereby providing detailed information on cerebral physiology in health and disease not previously amenable by noninvasive methods.

Kids, candy, brain and behavior: age differences in responses to candy gains and losses

•

Neural responses to candy loss feedback vary greatly between children and adults.

•

Children and adults show largely similar neural responses to candy gain feedback.

•

Striatal responses to candy gains/losses are influenced by age and by task behavior.

•

Insular responses to loss relate to age, even when controlling for task behavior.

The development of reward-related neural systems, from adolescence through adulthood, has received much recent attention in the developmental neuroimaging literature. However, few studies have investigated behavioral and neural responses to both gains and losses in pre-pubertal child populations. To address this gap in the literature, in the present study healthy children aged 7–11 years and young-adults completed an fMRI card-guessing game using candy pieces delivered post-scan as an incentive. Age differences in behavioral and neural responses to candy gains/losses were investigated. Adults and children displayed similar responses to gains, but robust age differences were observed following candy losses within the caudate, thalamus, insula, and hippocampus. Interestingly, when task behavior was included as a factor in post hoc mediation analyses, activation following loss within the caudate/thalamus related to task behavior and relationships with age were no longer significant. Conversely, relationships between response to loss and age within the hippocampus and insula remained significant even when controlling for behavior, with children showing heightened loss responses within the dorsal/posterior insula. These results suggest that both age and task behavior influence responses within the extended reward circuitry, and that children seem to be more sensitive than adults to loss feedback particularly within the dorsal/posterior insula.

Investigating the properties of the hemodynamic response function after mild traumatic brain injury

Abstract Although several functional magnetic resonance imaging (fMRI) studies have been conducted in human models of mild traumatic brain injury (mTBI), to date no studies have explicitly examined how injury may differentially affect both the positive phase of the hemodynamic response function (HRF) as well as the post-stimulus undershoot (PSU). Animal models suggest that the acute and semi-acute stages of mTBI are associated with significant disruptions in metabolism and to the microvasculature, both of which could impact on the HRF. Therefore, fMRI data were collected on a cohort of 30 semi-acute patients with mTBI (16 males; 27.83±9.97 years old; 13.00±2.18 years of education) and 30 carefully matched healthy controls (HC; 16 males; 27.17±10.08 years old; 13.37±2.31 years of education) during a simple sensory-motor task. Patients reported increased cognitive, somatic, and emotional symptoms relative to controls, although no group differences were detected on traditional neuropsychological examination. There were also no differences between patients with mTBI and controls on fMRI data using standard analytic techniques, although mTBI exhibited a greater volume of activation during the task qualitatively. A significant Group×Time interaction was observed in the right supramarginal gyrus, bilateral primary and secondary visual cortex, and the right parahippocampal gyrus. The interaction was the result of an earlier time-to-peak and positive magnitude shift throughout the estimated HRF in patients with mTBI relative to HC. This difference in HRF shape combined with the greater volume of activated tissue may be indicative of a potential compensatory mechanism to injury. The current study demonstrates that direct examination and modeling of HRF characteristics beyond magnitude may provide additional information about underlying neuropathology that is not available with more standard fMRI analyses.

High temporal resolution MRI quantification of global cerebral metabolic rate of oxygen consumption in response to apneic challenge

We present a technique for quantifying global cerebral metabolic rate of oxygen consumption (CMRO2) in absolute physiologic units at 3-second temporal resolution and apply the technique to quantify the dynamic CMRO2 response to volitional apnea. Temporal resolution of 3 seconds was achieved via a combination of view sharing and superior sagittal sinus-based estimation of total cerebral blood flow (tCBF) rather than tCBF measurement in the neck arteries. These modifications were first validated in three healthy adults and demonstrated to produce minimal errors in image-derived blood flow and venous oxygen saturation (SvO2) values. The technique was then applied in 10 healthy adults during an apnea paradigm of three repeated 30-second breath-holds. Subject-averaged baseline tCBF, arteriovenous oxygen difference (AVO2D), and CMRO2 were 48.6 ± 7.0 mL/100 g per minute, 29.4 ± 3.4 %HbO2, and 125.1 ± 11.4 μmol/100 g per minute, respectively. Subject-averaged maximum changes in tCBF and AVO2D were 43.5 ± 9.4% and -32.1 ± 5.7%, respectively, resulting in a small (6.0 ± 3.5%) but statistically significant (P=0.00044, two-tailed t-test) increase in average end-apneic CMRO2. This method could be used to investigate neurometabolic-hemodynamic relationships in normal physiology, to better define the biophysical origins of the BOLD signal, and to quantify neurometabolic responsiveness in diseases of altered neurovascular reactivity.

Model of the transient neurovascular response based on prompt arterial dilation

Brief neural stimulation results in a stereotypical pattern of vascular and metabolic response that is the basis for popular brain-imaging methods such as functional magnetic resonance imagine. However, the mechanisms of transient oxygen transport and its coupling to cerebral blood flow (CBF) and oxygen metabolism (CMRO2) are poorly understood. Recent experiments show that brief stimulation produces prompt arterial vasodilation rather than venous vasodilation. This work provides a neurovascular response model for brief stimulation based on transient arterial effects using one-dimensional convection-diffusion transport. Hemoglobin oxygen dissociation is included to enable predictions of absolute oxygen concentrations. Arterial CBF response is modeled using a lumped linear flow model, and CMRO2 response is modeled using a gamma function. Using six parameters, the model successfully fit 161/166 measured extravascular oxygen time courses obtained for brief visual stimulation in cat cerebral cortex. Results show how CBF and CMRO2 responses compete to produce the observed features of the hemodynamic response: initial dip, hyperoxic peak, undershoot, and ringing. Predicted CBF and CMRO2 response amplitudes are consistent with experimental measurements. This model provides a powerful framework to quantitatively interpret oxygen transport in the brain; in particular, its intravascular oxygen concentration predictions provide a new model for fMRI responses.

The physics of functional magnetic resonance imaging (fMRI)

Functional magnetic resonance imaging (fMRI) is a methodology for detecting dynamic patterns of activity in the working human brain. Although the initial discoveries that led to fMRI are only about 20 years old, this new field has revolutionized the study of brain function. The ability to detect changes in brain activity has a biophysical basis in the magnetic properties of deoxyhemoglobin, and a physiological basis in the way blood flow increases more than oxygen metabolism when local neural activity increases. These effects translate to a subtle increase in the local magnetic resonance signal, the blood oxygenation level dependent (BOLD) effect, when neural activity increases. With current techniques, this pattern of activation can be measured with resolution approaching 1 mm(3) spatially and 1 s temporally. This review focuses on the physical basis of the BOLD effect, the imaging methods used to measure it, the possible origins of the physiological effects that produce a mismatch of blood flow and oxygen metabolism during neural activation, and the mathematical models that have been developed to understand the measured signals. An overarching theme is the growing field of quantitative fMRI, in which other MRI methods are combined with BOLD methods and analyzed within a theoretical modeling framework to derive quantitative estimates of oxygen metabolism and other physiological variables. That goal is the current challenge for fMRI: to move fMRI from a mapping tool to a quantitative probe of brain physiology.

Noninvasive functional imaging of cerebral blood volume with vascular-space-occupancy (VASO) MRI

Functional MRI (fMRI) based on changes in cerebral blood volume (CBV) can probe directly vasodilatation and vasoconstriction during brain activation or physiologic challenges, and can provide important insights into the mechanism of blood oxygenation level-dependent (BOLD) signal changes. At present, the most widely used CBV fMRI technique in humans is called vascular-space-occupancy (VASO) MRI, and this article provides a technical review of this method. VASO MRI utilizes T1 differences between blood and tissue to distinguish between these two compartments within a voxel, and employs a blood-nulling inversion recovery sequence to yield an MR signal proportional to 1 - CBV. As such, vasodilatation will result in a VASO signal decrease and vasoconstriction will have the reverse effect. The VASO technique can be performed dynamically with a temporal resolution comparable with several other fMRI methods, such as BOLD or arterial spin labeling (ASL), and is particularly powerful when conducted in conjunction with these complementary techniques. The pulse sequence and imaging parameters of VASO can be optimized such that the signal change is predominantly of CBV origin, but careful considerations should be taken to minimize other contributions, such as those from the BOLD effect, cerebral blood flow (CBF) and cerebrospinal fluid (CSF). The sensitivity of the VASO technique is the primary disadvantage when compared with BOLD, but this technique is increasingly demonstrating its utility in neuroscientific and clinical applications.

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.

Elevated arteriolar cerebral blood volume in prodromal Huntington's disease

BACKGROUND

Neurovascular alterations have been implicated in the pathophysiology of Huntington's disease (HD). Because arterioles are most responsive to metabolic alterations, arteriolar cerebral blood volume (CBVa) is an important indicator of cerebrovascular regulation. The objective of this pilot study was to investigate potential neurovascular (CBVa ) abnormality in prodromal-HD patients and compare it with the widely used imaging marker: brain atrophy.

METHODS

CBVa and brain volumes were measured with ultra-high-field (7.0-Telsa) magnetic resonance imaging in seven prodromal-HD patients and nine age-matched controls.

RESULTS

Cortical CBVa was elevated significantly in prodromal-HD patients compared with controls (relative difference, 38.5%; effect size, 1.48). Significant correlations were found between CBVa in the frontal cortex and genetic measures. By contrast, no significant brain atrophy was detected in the prodromal-HD patients.

CONCLUSIONS

CBVa may be abnormal in prodromal-HD, even before substantial brain atrophy occurs. Further investigation with a larger cohort and longitudinal follow-up is merited to determine whether CBVa could be used as a potential biomarker for clinical trials.

Physiologic underpinnings of negative BOLD cerebrovascular reactivity in brain ventricles

With a growing need for specific biomarkers in vascular diseases, there has been a surging interest in mapping cerebrovascular reactivity (CVR) of the brain. This index can be measured by conducting a hypercapnia challenge while acquiring blood-oxygenation-level-dependent (BOLD) signals. A BOLD signal increase with hypercapnia is the expected outcome and represents the majority of literature reports; in this work we report an intriguing observation of an apparently negative BOLD CVR response at 3T, during inhalation of 5% CO2 with balance medical air. These "negative-CVR" clusters were specifically located in the ventricular regions of the brain, where CSF is abundant and results in an intense baseline signal. The amplitude of the CVR response was -0.51±0.44% (N=14, age 26±4 years). We hypothesized that this observation might not be due to a decrease in oxygenation but rather a volume effect in which bright CSF signal is replaced by a less intensive blood signal as a result of vasodilation. To test this, we performed an inversion-recovery (IR) experiment to suppress the CSF signal (N=10, age 27±5 years). This maneuver in imaging sequence reversed the sign of the signal response (to 0.66±0.25%), suggesting that the volume change was the predominant reason for the apparently negative CVR in the BOLD experiment. Further support of this hypothesis was provided by a BOLD hyperoxia experiment, in which no voxels showed a negative response, presumably because vasodilation is not usually associated with this challenge. Absolute CBF response to hypercapnia was measured in a new group of subjects (N=8, age 29±7 years) and it was found that CBF in ventricular regions increased by 48% upon CO2 inhalation, suggesting that blood oxygenation most likely increased rather than decreased. The findings from this study suggest that CO2 inhalation results in the dilation of ventricular vessels accompanied by shrinkage in CSF space, which is responsible for the apparently negative CVR in brain ventricles.

Sensitivity of BOLD response to increasing visual contrast: spin echo versus gradient echo EPI

Previous evidence showed that spin-echo (SE) BOLD signals offer an increased linearity and promptness with respect to gradient-echo (GE) acquisition, possibly providing a more accurate estimate of the amplitude of neuronal activity. However there is no evidence that the two sequences differ in representing different activation levels due to changes in stimulus intensity. In this study at 3T we compared GE and SE BOLD responses to visual stimuli at increasing contrast levels (5%, 20%, 60%, and 100%). Both sequences showed a monotonic increase of the BOLD response with stimulus contrast. While the larger sensitivity of GE yielded overall larger signal changes, step-wise increase in activation for GE was significant only when comparing 20% with 5% contrast, whereas for SE a significant increase was observed also when comparing 60% with 20% contrast. Moreover, BOLD responses normalized to the lowest contrast showed that relative increases of SE fMRI signal with increasing stimulus strength are larger with respect to the corresponding values of GE signal. This difference was observed also when excluding voxels attributed to large vessels, suggesting a non negligible role of the extravascular contribution to the modulation of SE fMRI signal with stimulus intensity. These results are shown to be in agreement with theoretical predictions that we derived from a recently proposed model of GE and SE functional signals. The present findings suggest that, despite the limited increase in functional localization accuracy at low field, SE fMRI might offer a potential advantage in distinguishing different levels of stimulus-evoked brain activity.