Simultaneous Blood Oxygen Level-dependent Research Articles

Modelling the depth-dependent VASO and BOLD responses in human primary visual cortex.

Functional magnetic resonance imaging (fMRI) using a blood-oxygenation-level-dependent (BOLD) contrast is a common method for studying human brain function noninvasively. Gradient-echo (GRE) BOLD is highly sensitive to the blood oxygenation change in blood vessels; however, the spatial signal specificity can be degraded due to signal leakage from activated lower layers to superficial layers in depth-dependent (also called laminar or layer-specific) fMRI. Alternatively, physiological variables such as cerebral blood volume using the VAscular-Space-Occupancy (VASO) contrast have shown higher spatial specificity compared to BOLD. To better understand the physiological mechanisms such as blood volume and oxygenation changes and to interpret the measured depth-dependent responses, models are needed which reflect vascular properties at this scale. For this purpose, we extended and modified the "cortical vascular model" previously developed to predict layer-specific BOLD signal changes in human primary visual cortex to also predict a layer-specific VASO response. To evaluate the model, we compared the predictions with experimental results of simultaneous VASO and BOLD measurements in a group of healthy participants. Fitting the model to our experimental data provided an estimate of CBV change in different vascular compartments upon neural activity. We found that stimulus-evoked CBV change mainly occurs in small arterioles, capillaries, and intracortical arteries and that the contribution from venules and ICVs is smaller. Our results confirm that VASO is less susceptible to large vessel effects compared to BOLD, as blood volume changes in intracortical arteries did not substantially affect the resulting depth-dependent VASO profiles, whereas depth-dependent BOLD profiles showed a bias towards signal contributions from intracortical veins.

Imaging brain function with simultaneous BOLD and viscoelasticity contrast: fMRI/fMRE

Magnetic resonance elastography (MRE) is emerging as a new tool for studying viscoelastic changes in the brain resulting from functional processes. Here, we demonstrate a novel time series method to generate robust functional magnetic resonance elastography (fMRE) activation maps in response to a visual task with a flashing checkerboard stimulus. Using a single-shot spin-echo (SS-SE) pulse sequence, the underlying raw images inherently contain blood-oxygen-level dependent (BOLD) contrast, allowing simultaneous generation of functional magnetic resonance imaging (fMRI) activation maps from the magnitude and functional magnetic resonance elastography (fMRE) maps from the phase. This allows an accurate comparison of the spatially localized stiffness (fMRE) and BOLD (fMRI) changes within a single scan, eliminating confounds inherent in separately acquired scans. Results indicate that tissue stiffness within the visual cortex increases 6–11% with visual stimuli, whereas the BOLD signal change was 1–2%. Furthermore, the fMRE and fMRI activation maps have strong spatial overlap within the visual cortex, providing convincing evidence that fMRE is possible in the brain. However, the fMRE temporal SNR (tSNRfMRE) maps are heterogeneous across the brain. Using a dictionary matching approach to characterize the time series, the viscoelastic changes are consistent with a viscoelastic response function (VRF) time constant of 12.1 ​s ± 3.0 ​s for a first-order exponential decay, or a shape parameter of 8.1 ​s ± 1.4 ​s for a gamma-variate.

Fiber-optic implant for simultaneous fluorescence-based calcium recordings and BOLD fMRI in mice.

Despite the growing popularity of blood oxygen level-dependent (BOLD) functional MRI (fMRI), understanding of its underlying principles is still limited. This protocol describes a technique for simultaneous measurement of neural activity using fluorescent calcium indicators together with the corresponding hemodynamic BOLD fMRI response in the mouse brain. Our early work using small-molecule fluorophores in rats gave encouraging results but was limited to acute measurements using synthetic dyes. Our latest procedure combines fMRI with optical detection of cell-type-specific virally delivered GCaMP6, a genetically encoded calcium indicator (GECI). GCaMP6 fluorescence, which increases upon calcium binding, is collected by a chronically implanted optical fiber, allowing longitudinal studies in mice. The chronic implant, placed horizontally on the skull, has an angulated tip that reflects light into the brain and is connected via fiber optics to a remote optical setup. The technique allows access to the neocortex and does not require adaptations of commercial MRI hardware. The hybrid approach permits fiber-optic calcium recordings with simultaneous artifact-free BOLD fMRI with full brain coverage and 1-s temporal resolution using standard gradient-echo echo-planar imaging (GE-EPI) sequences. The method provides robust, cell-type-specific readouts to link neural activity to BOLD signals, as emonstrated for task-free ('resting-state') conditions and in response to hind-paw stimulation. These results highlight the power of fiber photometry combined with fMRI, which we aim to further advance in this protocol. The approach can be easily adapted to study other molecular processes using suitable fluorescent indicators.

Multiband multi-echo simultaneous ASL/BOLD for task-induced functional MRI.

Typical simultaneous blood oxygenation-level dependent (BOLD) and arterial spin labeling (ASL) sequences acquire two echoes, one perfusion-sensitive and one BOLD-sensitive. However, for ASL, spatial resolution and brain coverage are limited due to the T1 decay of the labeled blood. This study applies a sequence combining a multiband acquisition with four echoes for simultaneous BOLD and pseudo-continuous ASL (pCASL) echo planar imaging (MBME ASL/BOLD) for block-design task-fMRI. A multiband acceleration of four was employed to increase brain coverage and reduce slice-timing effects on the ASL signal. Multi-echo independent component analysis (MEICA) was implemented to automatically denoise the BOLD signal by regressing non-BOLD components. This technique led to increased temporal signal-to-noise ratio (tSNR) and BOLD sensitivity. The MEICA technique was also modified to denoise the ASL signal by regressing artifact and BOLD signals from the first echo time-series. The MBME ASL/BOLD sequence was applied to a finger-tapping task functional MRI (fMRI) experiment. Signal characteristics and activation were evaluated using single echo BOLD, combined ME BOLD, combined ME BOLD after MEICA denoising, perfusion-weighted (PW), and perfusion-weighted after MEICA denoising time-series. The PW data was extracted using both surround subtraction and high-pass filtering followed by demodulation. In addition, the CBF/BOLD response ratio and CBF/BOLD coupling were analyzed. Results showed that the MEICA denoising procedure significantly improved the BOLD signal, leading to increased BOLD sensitivity, tSNR, and activation statistics compared to conventional single echo BOLD data. At the same time, the denoised PW data showed increased tSNR and activation statistics compared to the non-denoised PW data. CBF/BOLD coupling was also increased using the denoised ASL and BOLD data. Our preliminary data suggest that the MBME ASL/BOLD sequence can be employed to collect whole-brain task-fMRI with improved data quality for both BOLD and PW time series, thus improving the results of block-design task fMRI.

Simultaneous BOLD detection and incomplete fMRI data reconstruction.

The problem of simultaneous blood oxygenation level dependent (BOLD) detection and data completion is addressed in this paper. It is assumed that a set of fMRI data with significant number of missing samples is available and the aim is to recover those samples with least possible quality degradation. At the same time, BOLD should be detected. We propose a new cost function comprising both BOLD detection and data reconstruction terms. A solution based on singular value thresholding and sparsity-inducing approach is proposed. Due to the low-rank nature of the fMRI data, it is expected that the related techniques to be very effective for reconstruction. Extensive experiments are conducted on different datasets in noisy conditions. The achieved results, both in terms of data quality and data analysis accuracy, are promising and confirm that the proposed method can be effective for recovery of compressed/incomplete fMRI data. Several state-of-the art image reconstruction techniques are compared with the proposed method. In addition, the results of applying general linear model (GLM) using statistical parameter mapping (SPM) toolbox are compared with those of the proposed method.

Coherence of BOLD signal and electrical activity in the human brain during deep sevoflurane anesthesia.

IntroductionChanges in neural activity induce changes in functional magnetic resonance (fMRI) blood oxygenation level dependent (BOLD) signal. Commonly, increases in BOLD signal are ascribed to cellular excitation.ObjectiveThe relationship between electrical activity and BOLD signal in the human brain was probed on the basis of burst suppression EEG. This condition includes two distinct states of high and low electrical activity.MethodsResting‐state simultaneous EEG and BOLD measurements were acquired during deep sevoflurane anesthesia with burst suppression EEG in nineteen healthy volunteers. Afterwards, fMRI volumes were assigned to one of the two states (burst or suppression) as defined by the EEG.ResultsIn the frontal, parietal and temporal lobes as well as in the basal ganglia, BOLD signal increased after burst onset in the EEG and decreased after onset of EEG suppression. In contrast, BOLD signal in the occipital lobe was anticorrelated to electrical activity. This finding was obtained consistently in a general linear model and in raw data.ConclusionsIn human brains exhibiting burst suppression EEG induced by sevoflurane, the positive correlation between BOLD signal and electrical brain activity could be confirmed in most gray matter. The exceptional behavior of the occipital lobe with an anticorrelation of BOLD signal and electrical activity might be due to specific neurovascular coupling mechanisms that are pronounced in the deeply anesthetized brain.

Dual echo vessel-encoded ASL for simultaneous BOLD and CBF reactivity assessment in patients with ischemic cerebrovascular disease.

Blood oxygenation level-dependent (BOLD)-weighted and vessel-encoded arterial spin labeling (VE-ASL) MRI provide complementary information and can be used in sequence to gauge hemodynamic contributions to cerebrovascular reactivity. Here, cerebrovascular reactivity is assessed using dual echo VE-ASL MRI to understand how VE labeling preparations influence BOLD and ASL contrast in flow-limited and healthy perfusion territories. Patients (n = 12; age = 55 +/- 14 years; 6F/6M) presenting with ischemic steno-occlusive cerebrovascular disease underwent 3.0T angiographic imaging, T1 -weighted structural, and planning-free dual echo hypercarbic hyperoxic (i.e., carbogen) VE-ASL MRI. Vasculopathy extent, timecourses, and cerebrovascular reactivity (signal change and Z-statistic) for different VE-ASL images were contrasted across flow territories and Bonferroni-corrected P-values reported. BOLD cerebrovascular reactivity (i.e., long-TE VE-ASL) Z-statistics were similarly sensitive to asymmetric disease (P ≤ 0.002) regardless of labeling scenario. Cerebral blood flow reactivity correlated significantly with BOLD reactivity (Z-statistic). However, BOLD signal changes did not differ significantly between labeling scenarios (P > 0.003) or across territories (P > 0.002), indicating BOLD signal changes in response to carbogen offer low sensitivity to lateralizing disease. Dual echo VE-ASL can provide simultaneous cerebral blood flow and qualitative BOLD contrast consistent with lateralizing disease severity in patients with asymmetric steno-occlusive disease. The methodological strengths and limitations of composite BOLD and VE-ASL measurements in the clinic are discussed.

Localization of the hand motor area by arterial spin labeling and blood oxygen level‐dependent functional magnetic resonance imaging

The new clinically available arterial spin labeling (ASL) perfusion imaging sequences present some advantages relatively to the commonly used blood oxygen level-dependent (BOLD) method for functional brain studies using magnetic resonance imaging (MRI). In particular, regional cerebral blood flow (CBF) changes are thought to be more directly related with neuronal activation. In this study, we aimed to investigate the accuracy of the functional localization of the hand motor area obtained by simultaneous CBF and BOLD contrasts provided by ASL functional MRI (fMRI) and compare it with a standard BOLD fMRI protocol. For this purpose, we measured the distance between the center of gravity of the activation clusters obtained with each contrast (CBF, BOLD(ASL), and Standard BOLD) and 11 positions defined on a well-established anatomical landmark of the hand motor area (the omega in the axial plane of the precentral gyrus). We found that CBF measurements were significantly closer to the anatomical landmark than the ones obtained using either simultaneous BOLD(ASL) or standard BOLD contrasts. Moreover, we also observed reduced intersubject variability of the functional localization, as well as percent signal change, for CBF relative to both BOLD contrast measurements. In conclusion, our results add further evidence in support to the notion that CBF provides a more accurate localization of motor activation than BOLD contrast, indicating that ASL may be an appropriate technique for clinical fMRI studies.

Simultaneous Dynamic Blood Oxygen Level-Dependent Magnetic Resonance Imaging of Foot and Calf Muscles

To demonstrate the feasibility of simultaneous blood oxygen level-dependent (BOLD) magnetic resonance imaging of calf and foot muscles and investigate age-related changes of BOLD signal changes during ischemia and postocclusive hyperemia in healthy volunteers. In this study, 15 healthy elderly volunteers (mean age: 69.0 +/- 7.4 years) and 15 healthy young volunteers (mean age: 26.1 +/- 3.9 years) were enrolled. In both legs, simultaneous BOLD imaging of calf and foot muscles was performed at 1.5 Tesla. Short-term ischemia and consecutive reactive hyperemia were provoked by a cuff-compression paradigm. T2*-weighted signal time courses were obtained from foot and calf muscles simultaneously. Ischemia was assessed by T2* minimum ischemic value (MIV) and the time to half ischemic minimum (THIM). Reperfusion was assessed by the time to reach T2* half hyperemia peak (THHP). Reactive hyperemia was characterized by hyperemia peak value (HPV), time to peak (TTP), and relative T2* change from end of ischemia to HPV (deltaS). Parameter differences were assessed using a 2-sided Student t test. Dynamic BOLD measurement of foot and calf muscles was techniqually feasible and successful in all volunteers.In comparison, THIM was significantly longer in elderly than in young volunteers for calf (P < 0.01) muscles (young: 28.9 +/- 3.7 seconds; elderly: 57.8 +/- 31.4 seconds) and foot (P = 0.01) muscles (young: 36.8 +/- 25.5 seconds; elderly: 56.6 +/- 31.7 seconds). MIV relative to baseline T2*-signal was significantly (P < 0.01) lower in the elderly for calf (young: 96.0% +/- 2.6%; elderly: 91.3% +/- 4.4%) and foot (young: 95.8% +/- 2.5%; elderly: 91.1% +/- 8.2%) muscles. TTP was significantly (calf: P = 0.01; foot: P = 0.02) delayed in the elderly (elderly calf: 103.0 +/- 92.7 seconds and foot: 157.1 +/- 109.9 seconds vs. young calf: 54.8 +/- 42.1 seconds and foot: 95.1 +/- 77.6 seconds). HPV was significantly (calf: P < 0.01 and foot: P = 0.03) higher in (young calf: 114.1% +/- 7.2% and foot: 105.8% +/- 3.3% vs. elderly calf: 104.0% +/- 2.1% and foot: 103.9% +/- 3.2% seconds) young volunteers.In a muscle-group comparison, no significant differences in THIM and MIV were observed between calf and foot.THHP was significantly (P = 0.02) longer in foot muscles (foot young: 32.5 +/- 29.8 seconds and elderly: 34.1 +/- 25.0 seconds vs. calf young: 16.8 +/- 14.1 seconds and elderly: 23.6 +/- 14.1 seconds) of both age groups. TTP was significantly (P = 0.01 and 0.02) longer in foot muscles (foot young: 95.1 +/- 77.6 seconds and elderly: 157.1 +/- 109.9 seconds vs. calf young: 54.8 +/- 42.1 seconds and elderly: 103.0 +/- 92.7 seconds) of both age groups. HPV was lower (P < 0.01) in foot muscles of the young (calf: 114.1% +/- 7.2% vs. foot: 105.8 +/- 3.3%). Simultaneous BOLD-imaging of calf and foot muscles is feasible and reveals statistically significant age-related differences during ischemia and postocclusive hyperemia in healthy volunteers.

Quantitative fMRI using hyperoxia calibration: Reproducibility during a cognitive Stroop task

Arterial spin labelling allows simultaneous measurement of both the blood-oxygenation-level-dependent (BOLD) and the cerebral blood flow (CBF) response to changes in neural activity. The addition of a hypercapnia or hyperoxia calibration allows additional quantification of changes in the cerebral metabolic rate of oxygen (CMRO2). In this study we test the reproducibility of measurements derived using the hyperoxia approach, during a cognitive Stroop task. A QUIPSSII sequence is used at 3 T to collect simultaneous CBF and BOLD signal during two 3 min periods of hyperoxia and an 8 min Stroop task. Hyperoxia was administered via an open system and end-tidal values were sampled via a nasal cannula; average end-tidal values of 60% were reached. This procedure is repeated to allow the reproducibility of the estimated parameters to be tested. The use of a cognitive Stroop task allows testing of the measurements in frontal and parietal regions as well as sensorimotor areas in which previous studies have been focussed. We find reduced reproducibility of the calculated parameters compared to the hypercapnia approach, thought to be attributable to lower absolute BOLD and CBF responses. In particular we do not find ‘n’ to have improved reproducibility compared to other parameters, as has been found in previous work using the hypercapnia approach. Across all brain areas we report a value of ΔCMRO2 of 12% and neurovascular coupling constant n of 2.5. Interestingly we find n to be higher in parietal and frontal areas in comparison to the primary motor cortex.

Brain metabolic fluctuations in absence of stimuli

Mapping of brain function with blood oxygen level dependent (BOLD) magnetic resonance imaging (MRI) BOLD-MRI relies on the detection of focal changes in cerebral blood flow (CBF) in response to conditioned stimuli. However, in the absence of stimuli, the brain continues to show temporal fluctuations in blood flow as measured with BOLD and perfusion-based MRI (Biswal, NMR. Biomed. 1997). These fluctuations are substantially synchronized across functional regions that have an apparent functional relationship, and therefore might allow mapping and classification of the networks that underlie human brain function without the need for carefully conditioned stimuli. Despite their potential functional siginificance, the origin of the blood flow fluctuations has not been established. Potential sources include cognitive processes, fluctuations in vigilance and conscious awareness, and homeostatic (restorative) processes. Alternatively, the fluctuations might be caused by purely vascular events (e.g. vasomotor effect), without a substantial metabolic or neuronal component. In this study we used MRI methods to establish whether blood flow fluctuations in absence of stimuli subserve a metabolic process. For this purpose, we performed simultaneous BOLD and perfusion MRI at 3.0 T on 7 normal volunteers during a visual task (5 min checkerboard-grey paradigm) and an extended (25 min) rest period. The ratio between BOLD and perfusion changes was used as an indicator of metabolism (Hoge, Proc. Nat. Acad. Sci. 1999). As reference, a 10-minute breath-holding paradigm was inserted between the two conditions. The paradigm consisted of 5 stages of 40 s breath-hold, 80 s normal breathing. Physiologic monitoring included respiratory and cardiac cycles, as well as scalp potentials using a 64-channel EEG system. During rest conditions, we found substantial BOLD and perfusion signal fluctuations in most of neocortex, with strong correlations between functionally related regions. No significant correlations were found with cardiac and respiratory cycles. Frequency analysis of the fluctuations suggested that most energy was concentrated in the 0.01–0.1 Hz band. The strongest fluctuations were observed during early (non-REM) sleep, as identified from EEG. Within visual cortex (VC), the amplitude of the fluctuations reached a level similar to that evoked with the task. BOLD and perfusion signals were highly correlated over most of cortex. The ratio of their fluctuation amplitude in VC averaged 1.66(0.41), 2.82(0.70), and 1.80(0.57) for task, breath holding paradigm, and rest (sleep) respectively. The difference between task and rest was insignificant (p=0.56), while the ratios during task and rest were significantly (p=0.0003/0.017) lower than during the breath-holding paradigm. The lower ratios are indicative of a relatively high oxygen extraction, suggesting involvement of a metabolic process during task and rest, and less so during breathholding. The findings indicate that spatially correlated fluctuations in CBF in absence of stimuli have a metabolic component and continue during early sleep. This has important implication for the design of brain mapping experiments and could shed light on the physiologic processes present during of rest and sleep. Further experiments are planned to investigate a potential role of processes facilitating synaptic plasticity (Huber, Nature 2004).

Temporal resolving power of perfusion- and BOLD-based event-related functional MRI.

Functional magnetic resonance imaging (fMRI) based on both perfusion and blood oxygenation level-dependent (BOLD) contrasts has been widely applied in spatiotemporal mapping of the human brain function. Temporal resolving power of fMRI is limited by the smoothed hemodynamic response function dispersed from the neuronal activity. In this study, temporal modulation transfer functions were utilized to quantify the resolving powers of perfusion and BOLD fMR signals in time domain. The impulse response function was determined using brief visual stimulations and event-related image acquisition schemes. An important feature of arterial spin labeling techniques is that quantitative perfusion and BOLD signals could be simultaneously acquired. This simultaneous BOLD response may arise from signals that are more proximal to capillary beds, and its temporal resolution may be different from that of the typical BOLD response. Therefore, we assessed and compared the temporal resolving capabilities of perfusion, simultaneous BOLD, and the typical BOLD response obtained from the gradient echo EPI pulse sequence. Full-width-at-half-maximums of perfusion and simultaneous BOLD measurements were significantly smaller than that of BOLD ones (4.3+/- 0.6 s vs 5.5 +/- 0.9 s, p<0.02 and 4.7 +/- 1.3 s vs 5.5 +/- 0.9 s, p<0.01, respectively). The corresponding temporal resolving powers of perfusion and simultaneous BOLD signals were statistically better than that of BOLD signals (0.23 +/- 0.03 Hz vs 0.17 +/- 0.02 Hz, p<0.01 and 0.21 +/- 0.04 Hz vs 0.17 +/- 0.02 Hz, p<0.01, respectively). Our results showed that the typical BOLD response was significantly smoothed from the perfusion response, thus resulting in a degraded temporal resolving power. However, results from the simultaneous BOLD and perfusion measurements were not significantly different. Biophysical implications of the experimental outcomes were further investigated using a computer simulation based on the Balloon model. By fitting the measured data into the model, an apparently longer transit time was obtained for the typical BOLD signal (1.7 s), comparing to that for the simultaneous BOLD one (1.2 s). Therefore, the simultaneous BOLD signal was regarded as less susceptible to the variations from local draining veins. Combining the simulation result with the significantly discrepant resolving powers between the two BOLD signals, we speculated that the blurred effects from large vessels played a predominant role that further reduced the temporal resolution of the BOLD-based fMRI from the perfusion response.