Theory and generalization of Monte Carlo models of the BOLD signal source

The dependency of the blood oxygenation level dependent (BOLD) signal on underlying hemodynamics is not well understood. Building a forward biophysical model of this relationship is important for the quantitative estimation of the hemodynamic changes and neural activity underlying functional magnetic resonance imaging (fMRI) signals. We have developed a general model of the BOLD signal which can model both intra- and extravascular signals for an arbitrary tissue model across a wide range of imaging parameters. The model of the BOLD signal was instantiated as a look-up-table (LuT), and was verified against concurrent fMRI and optical imaging measurements of activation induced hemodynamics. Magn Reson Med, 2008. © 2008 Wiley-Liss, Inc.

Long-latency reductions in gamma power predict hemodynamic changes that underlie the negative BOLD signal

Studiesthat use prolonged periods of sensory stimulation report associations between regional reductions in neural activity and negative blood oxygenation level-dependent (BOLD) signaling. However, the neural generators of the negative BOLD response remain to be characterized. Here, we use single-impulse electrical stimulation of the whisker pad in the anesthetized rat to identify components of the neural response that are related to “negative” hemodynamic changes in the brain. Laminar multiunit activity and local field potential recordings of neural activity were performed concurrently withtwo-dimensional optical imaging spectroscopy measuring hemodynamic changes. Repeated measurements over multiple stimulation trials revealed significant variations in neural responses across session and animal datasets. Within this variation, we found robust long-latency decreases (300 and 2000 ms after stimulus presentation) in gammaband power (30 – 80 Hz) in the middle-superficial cortical layers in regions surrounding the activated whisker barrel cortex. This reduction in gamma frequency activity was associated with corresponding decreases in the hemodynamic responses that drive the negative BOLD signal. These findings suggest a close relationship between BOLD responses and neural events that operate over time scales that outlast the initiating sensory stimulus, and provide important insights into the neurophysiological basis of negative neuroimaging signals.

Nonlinear estimation of neural processing time from BOLD signal with application to decision-making

The extraction of information about neural activity timing from BOLD signal is a challenging task as the shape of the BOLD curve does not directly reflect the temporal characteristics of electrical activity of neurons. In this work, we introduce the concept of neural processing time (NPT) as a parameter of the biophysical model of the hemodynamic response function (HRF). Through this new concept we aim to infer more accurately the duration of neuronal response from the highly nonlinear BOLD effect. The face validity and applicability of the concept of NPT are evaluated through simulations and analysis of experimental time series. The results of both simulation and application were compared with summary measures of HRF shape. The experiment that was analyzed consisted of a decision-making paradigm with simultaneous emotional distracters. We hypothesize that the NPT in primary sensory areas, like the fusiform gyrus, is approximately the stimulus presentation duration. On the other hand, in areas related to processing of an emotional distracter, the NPT should depend on the experimental condition. As predicted, the NPT in fusiform gyrus is close to the stimulus duration and the NPT in dorsal anterior cingulate gyrus depends on the presence of an emotional distracter. Interestingly, the NPT in right but not left dorsal lateral prefrontal cortex depends on the stimulus emotional content. The summary measures of HRF obtained by a standard approach did not detect the variations observed in the NPT. Hum Brain Mapp, 2012. (C) 2010 Wiley Periodicals, Inc.

Developmental changes of BOLD signal correlations with global human EEG power and synchronization during working memory

In humans, theta band (5-7 Hz) power typically increases when performing cognitively demanding working memory (WM) tasks, and simultaneous EEG-fMRI recordings have revealed an inverse relationship between theta power and the BOLD (blood oxygen level dependent) signal in the default mode network during WM. However, synchronization also plays a fundamental role in cognitive processing, and the level of theta and higher frequency band synchronization is modulated during WM. Yet, little is known about the link between BOLD, EEG power, and EEG synchronization during WM, and how these measures develop with human brain maturation or relate to behavioral changes. We examined EEG-BOLD signal correlations from 18 young adults and 15 school-aged children for age-dependent effects during a load-modulated Sternberg WM task. Frontal load (in-)dependent EEG theta power was significantly enhanced in children compared to adults, while adults showed stronger fMRI load effects. Children demonstrated a stronger negative correlation between global theta power and the BOLD signal in the default mode network relative to adults. Therefore, we conclude that theta power mediates the suppression of a task-irrelevant network. We further conclude that children suppress this network even more than adults, probably from an increased level of task-preparedness to compensate for not fully mature cognitive functions, reflected in lower response accuracy and increased reaction time. In contrast to power, correlations between instantaneous theta global field synchronization and the BOLD signal were exclusively positive in both age groups but only significant in adults in the frontal-parietal and posterior cingulate cortices. Furthermore, theta synchronization was weaker in children and was--in contrast to EEG power--positively correlated with response accuracy in both age groups. In summary we conclude that theta EEG-BOLD signal correlations differ between spectral power and synchronization and that these opposite correlations with different distributions undergo similar and significant neuronal developments with brain maturation.

Blood-oxygen-level-dependent (BOLD) signal changes in total cortex and subcortex, pain, reward, and vigilance regions during mechanical pressure pain after morphine administration

Morphine is the most common clinical choice in the management of severe pain. Although the molecular mechanisms of morphine have already been characterized, the cerebral circuits by which it attenuates the sensation of pain have not yet been studied in humans. The objective of this two-arm (morphine versus placebo), between-subjects study was to examine whether morphine affects pain via pain-related cortical circuits, but also via reward regions that relate to the motivational state, as well as prefrontal regions that relate to vigilance as a result of morphine's sedative effects. Cortical activity was measured by the blood-oxygen-level-dependent (BOLD) signal changes using functional magnetic resonance imaging (fMRI). ^ The novelty of this study is at three levels: (i) to develop a methodology that will assess the average BOLD signal across subjects for the pain, reward, and vigilance cortical systems; (ii) to examine whether the reward and/or sedative effects of morphine are contributing factors to cortical regions associated with the motivational state and vigilance; and (iii) to propose a neuroanatomical model related to the opioid-sensitive effects of reward and sedation as a function of cortical activity related to pain in an effort to assess future analgesics. ^ Consistent with our hypotheses, our findings showed that the decrease in total pain-related volume activated between the post- and the pre-treatment morphine group was about 78%, while the post-treatment placebo group displayed only a 5% decrease when compared to pre-treatment levels of activation. The volume increase in reward regions was 451% in the post-treatment compared to the pre-treatment morphine condition. Finally, the volumetric decrease in vigilance regions was 63% in the posttreatment compared to the pre-treatment morphine condition. ^ These findings imply that changes in the blood flow of the reward and vigilance regions may be contributing factors in producing the analgesic effect under morphine administration. Future studies need to replicate this study in a higher resolution fMRI environment and to assess the proposed neuroanatomical model in patient populations. The necessity of pain research is apparent, since pain cuts across different diseases especially chronic ones, and thus, is recognized as a vital public health developing area. ^

The Stability of the Blood Oxygenation Level-Dependent Functional MRI Response to Motor Tasks Is Altered in Patients With Chronic Ischemic Stroke

Background and Purpose-Functional MRI is a powerful tool to investigate recovery of brain function in patients with stroke. An inherent assumption in functional MRI data analysis is that the blood oxygenation level-dependent (BOLD) signal is stable over the course of the examination. In this study, we evaluated the validity of such assumption in patients with chronic stroke. Methods-Fifteen patients performed a simple motor task with repeated epochs using the paretic and the unaffected hand in separate runs. The corresponding BOLD signal time courses were extracted from the primary and supplementary motor areas of both hemispheres. Statistical maps were obtained by the conventional General Linear Model and by a parametric General Linear Model. Results-Stable BOLD amplitude was observed when the task was executed with the unaffected hand. Conversely, the BOLD signal amplitude in both primary and supplementary motor areas was progressively attenuated in every patient when the task was executed with the paretic hand. The conventional General Linear Model analysis failed to detect brain activation during movement of the paretic hand. However, the proposed parametric General Linear Model corrected the misdetection problem and showed robust activation in both primary and supplementary motor areas. Conclusions-The use of data analysis tools that are built on the premise of a stable BOLD signal may lead to misdetection of functional regions and underestimation of brain activity in patients with stroke. The present data urge the use of caution when relying on the BOLD response as a marker of brain reorganization in patients with stroke. (Stroke. 2010; 41:1921-1926.)

From EEG to BOLD: Brain mapping and estimating transfer functions in simultaneous EEG-fMRI acquisitions

Simultaneous acquisition of electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) aims to disentangle the description of brain processes by exploiting the advantages of each technique. Most studies in this field focus on exploring the relationships between fMRI signals and the power spectrum at some specific frequency bands (alpha, beta, etc.). On the other hand, brain mapping of EEG signals (e.g., interictal spikes in epileptic patients) usually assumes an haemodynamic response function for a parametric analysis applying the GLM, as a rough approximation. The integration of the information provided by the high spatial resolution of MR images and the high temporal resolution of EEG may be improved by referencing them by transfer functions, which allows the identification of neural driven areas without strong assumptions about haemodynamic response shapes or brain haemodynamic`s homogeneity. The difference on sampling rate is the first obstacle for a full integration of EEG and fMRI information. Moreover, a parametric specification of a function representing the commonalities of both signals is not established. In this study, we introduce a new data-driven method for estimating the transfer function from EEG signal to fMRI signal at EEG sampling rate. This approach avoids EEG subsampling to fMRI time resolution and naturally provides a test for EEG predictive power over BOLD signal fluctuations, in a well-established statistical framework. We illustrate this concept in resting state (eyes closed) and visual simultaneous fMRI-EEG experiments. The results point out that it is possible to predict the BOLD fluctuations in occipital cortex by using EEG measurements. (C) 2010 Elsevier Inc. All rights reserved.

Is optical imaging spectroscopy a viable measurement technique for the investigation of the negative BOLD phenomenon? A concurrent optical imaging spectroscopy and fMRI study at high field (7T)

Traditionally functional magnetic resonance imaging (fMRI) has been used to map activity in the human brain by measuring increases in the Blood Oxygenation Level Dependent (BOLD) signal. Often accompanying positive BOLD fMRI signal changes are sustained negative signal changes. Previous studies investigating the neurovascular coupling mechanisms of the negative BOLD phenomenon have used concurrent 2D-optical imaging spectroscopy (2D-OIS) and electrophysiology (Boorman et al., 2010). These experiments suggested that the negative BOLD signal in response to whisker stimulation was a result of an increase in deoxy-haemoglobin and reduced multi-unit activity in the deep cortical layers. However, Boorman et al. (2010) did not measure the BOLD and haemodynamic response concurrently and so could not quantitatively compare either the spatial maps or the 2D-OIS and fMRI time series directly. Furthermore their study utilised a homogeneous tissue model in which is predominantly sensitive to haemodynamic changes in more superficial layers. Here we test whether the 2D-OIS technique is appropriate for studies of negative BOLD. We used concurrent fMRI with 2D-OIS techniques for the investigation of the haemodynamics underlying the negative BOLD at 7 Tesla. We investigated whether optical methods could be used to accurately map and measure the negative BOLD phenomenon by using 2D-OIS haemodynamic data to derive predictions from a biophysical model of BOLD signal changes. We showed that despite the deep cortical origin of the negative BOLD response, if an appropriate heterogeneous tissue model is used in the spectroscopic analysis then 2D-OIS can be used to investigate the negative BOLD phenomenon.

Time course of blood oxygenation level-dependent signal response after theta burst transcranial magnetic stimulation of the frontal eye field

The aim of the current study was to examine the effect of theta burst repetitive transcranial magnetic stimulation (rTMS) on the blood oxygenation level-dependent (BOLD) activation during repeated functional magnetic resonance imaging (fMRI) measurements. Theta burst rTMS was applied over the right frontal eye field in seven healthy subjects. Subsequently, repeated fMRI measurements were performed during a saccade-fixation task (block design) 5, 20, 35, and 60 min after stimulation. We found that theta burst rTMS induced a strong and long-lasting decrease of the BOLD signal response of the stimulated frontal eye field at 20 and 35 min. Furthermore, less pronounced alterations of the BOLD signal response with different dynamics were found for remote oculomotor areas such as the left frontal eye field, the pre-supplementary eye field, the supplementary eye field, and both parietal eye fields. Recovery of the BOLD signal changes in the anterior remote areas started earlier than in the posterior remote areas. These results show that a) the major inhibitory impact of theta burst rTMS occurs directly in the stimulated area itself, and that b) a lower effect on remote, oculomotor areas can be induced.

BOLD correlates of continuously fluctuating epileptic activity isolated by independent component analysis

Combined EEG/fMRI recordings offer a promising opportunity to detect brain areas with altered BOLD signal during interictal epileptic discharges (IEDs). These areas are likely to represent the irritative zone, which is itself a reflection of the epileptogenic zone. This paper reports on the imaging findings using independent component analysis (ICA) to continuously quantify epileptiform activity in simultaneously acquired EEG and fMRI. Using ICA derived factors coding for the epileptic activity takes into account that epileptic activity is continuously fluctuating with each spike differing in amplitude, duration and maybe topography, including subthreshold epileptic activity besides clear IEDs and may thus increase the sensitivity and statistical power of combined EEG/fMRI in epilepsy. Twenty patients with different types of focal and generalized epilepsy syndromes were investigated. ICA separated epileptiform activity from normal physiological brain activity and artifacts. In 16/20 patients, BOLD correlates of epileptic activity matched the EEG sources, the clinical semiology, and, if present, the structural lesions. In clinically equivocal cases, the BOLD correlates aided to attribute proper diagnosis of the underlying epilepsy syndrome. Furthermore, in one patient with temporal lobe epilepsy, BOLD correlates of rhythmic delta activity could be employed to delineate the affected hippocampus. Compared to BOLD correlates of manually identified IEDs, the sensitivity was improved from 50% (10/20) to 80%. The ICA EEG/fMRI approach is a safe, non-invasive and easily applicable technique, which can be used to identify regions with altered hemodynamic effects related to IEDs as well as intermittent rhythmic discharges in different types of epilepsy.

BOLD correlates of EEG alpha phase-locking and the fMRI default mode network

Phase locking or synchronization of brain areas is a key concept of information processing in the brain. Synchronous oscillations have been observed and investigated extensively in EEG during the past decades. EEG oscillations occur over a wide frequency range. In EEG, a prominent type of oscillations is alpha-band activity, present typically when a subject is awake, but at rest with closed eyes. The spectral power of alpha rhythms has recently been investigated in simultaneous EEG/fMRI recordings, establishing a wide-range cortico-thalamic network. However, spectral power and synchronization are different measures and little is known about the correlations between BOLD effects and EEG synchronization. Interestingly, the fMRI BOLD signal also displays synchronous oscillations across different brain regions. These oscillations delineate so-called resting state networks (RSNs) that resemble the correlation patterns of simultaneous EEG/fMRI recordings. However, the nature of these BOLD oscillations and their relations to EEG activity is still poorly understood. One hypothesis is that the subunits constituting a specific RSN may be coordinated by different EEG rhythms. In this study we report on evidence for this hypothesis. The BOLD correlates of global EEG synchronization (GFS) in the alpha frequency band are located in brain areas involved in specific RSNs, e.g. the 'default mode network'. Furthermore, our results confirm the hypothesis that specific RSNs are organized by long-range synchronization at least in the alpha frequency band. Finally, we could localize specific areas where the GFS BOLD correlates and the associated RSN overlap. Thus, we claim that not only the spectral dynamics of EEG are important, but also their spatio-temporal organization.

Functional connectivity in BOLD and CBF data: similarity and reliability of resting brain networks

Resting-state functional connectivity (FC) fMRI (rs-fcMRI) offers an appealing approach to mapping the brain's intrinsic functional organization. Blood oxygen level dependent (BOLD) and arterial spin labeling (ASL) are the two main rs-fcMRI approaches to assess alterations in brain networks associated with individual differences, behavior and psychopathology. While the BOLD signal is stronger with a higher temporal resolution, ASL provides quantitative, direct measures of the physiology and metabolism of specific networks. This study systematically investigated the similarity and reliability of resting brain networks (RBNs) in BOLD and ASL. A 2×2×2 factorial design was employed where each subject underwent repeated BOLD and ASL rs-fcMRI scans on two occasions on two MRI scanners respectively. Both independent and joint FC analyses revealed common RBNs in ASL and BOLD rs-fcMRI with a moderate to high level of spatial overlap, verified by Dice Similarity Coefficients. Test-retest analyses indicated more reliable spatial network patterns in BOLD (average modal Intraclass Correlation Coefficients: 0.905±0.033 between-sessions; 0.885±0.052 between-scanners) than ASL (0.545±0.048; 0.575±0.059). Nevertheless, ASL provided highly reproducible (0.955±0.021; 0.970±0.011) network-specific CBF measurements. Moreover, we observed positive correlations between regional CBF and FC in core areas of all RBNs indicating a relationship between network connectivity and its baseline metabolism. Taken together, the combination of ASL and BOLD rs-fcMRI provides a powerful tool for characterizing the spatiotemporal and quantitative properties of RBNs. These findings pave the way for future BOLD and ASL rs-fcMRI studies in clinical populations that are carried out across time and scanners.

Engagement of multifocal neural circuits during recall of autobiographical happy events

Happy emotional states have not been extensively explored in functional magnetic resonance imaging studies using autobiographic recall paradigms. We investigated the brain circuitry engaged during induction of happiness by standardized script-driven autobiographical recall in 11 healthy subjects (6 males), aged 32.4 ± 7.2 years, without physical or psychiatric disorders, selected according to their ability to vividly recall personal experiences. Blood oxygen level-dependent (BOLD) changes were recorded during auditory presentation of personal scripts of happiness, neutral content and negative emotional content (irritability). The same uniform structure was used for the cueing narratives of both emotionally salient and neutral conditions, in order to decrease the variability of findings. In the happiness relative to the neutral condition, there was an increased BOLD signal in the left dorsal prefrontal cortex and anterior insula, thalamus bilaterally, left hypothalamus, left anterior cingulate gyrus, and midportions of the left middle temporal gyrus (P < 0.05, corrected for multiple comparisons). Relative to the irritability condition, the happiness condition showed increased activity in the left insula, thalamus and hypothalamus, and in anterior and midportions of the inferior and middle temporal gyri bilaterally (P < 0.05, corrected), varying in size between 13 and 64 voxels. Findings of happiness-related increased activity in prefrontal and subcortical regions extend the results of previous functional imaging studies of autobiographical recall. The BOLD signal changes identified reflect general aspects of emotional processing, emotional control, and the processing of sensory and bodily signals associated with internally generated feelings of happiness. These results reinforce the notion that happiness induction engages a wide network of brain regions.

Wavelet correlation between subjects: A time-scale data driven analysis for brain mapping using fMRI

Functional magnetic resonance imaging (fMRI) based on BOLD signal has been used to indirectly measure the local neural activity induced by cognitive tasks or stimulation. Most fMRI data analysis is carried out using the general linear model (GLM), a statistical approach which predicts the changes in the observed BOLD response based on an expected hemodynamic response function (HRF). In cases when the task is cognitively complex or in cases of diseases, variations in shape and/or delay may reduce the reliability of results. A novel exploratory method using fMRI data, which attempts to discriminate between neurophysiological signals induced by the stimulation protocol from artifacts or other confounding factors, is introduced in this paper. This new method is based on the fusion between correlation analysis and the discrete wavelet transform, to identify similarities in the time course of the BOLD signal in a group of volunteers. We illustrate the usefulness of this approach by analyzing fMRI data from normal subjects presented with standardized human face pictures expressing different degrees of sadness. The results show that the proposed wavelet correlation analysis has greater statistical power than conventional GLM or time domain intersubject correlation analysis. (C) 2010 Elsevier B.V. All rights reserved.

White matter microstructure underlying default mode network connectivity in the human brain

Resting state functional magnetic resonance imaging (fMRI) reveals a distinct network of correlated brain function representing a default mode state of the human brain The underlying structural basis of this functional connectivity pattern is still widely unexplored We combined fractional anisotropy measures of fiber tract integrity derived from diffusion tensor imaging (DTI) and resting state fMRI data obtained at 3 Tesla from 20 healthy elderly subjects (56 to 83 years of age) to determine white matter microstructure e 7 underlying default mode connectivity We hypothesized that the functional connectivity between the posterior cingulate and hippocampus from resting state fMRI data Would be associated with the white matter microstructure in the cingulate bundle and fiber tracts connecting posterior cingulate gyrus With lateral temporal lobes, medial temporal lobes, and precuneus This was demonstrated at the p<0001 level using a voxel-based multivariate analysis of covariance (MANCOVA) approach In addition, we used a data-driven technique of joint independent component analysis (ICA) that uncovers spatial pattern that are linked across modalities. It revealed a pattern of white matter tracts including cingulate bundle and associated fiber tracts resembling the findings from the hypothesis-driven analysis and was linked to the pattern of default mode network (DMN) connectivity in the resting state fMRI data Out findings support the notion that the functional connectivity between the posterior cingulate and hippocampus and the functional connectivity across the entire DMN is based oil distinct pattern of anatomical connectivity within the cerebral white matter (C) 2009 Elsevier Inc All rights reserved