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.

The resting-state neurovascular coupling relationship: rapid changes in spontaneous neural activity in the somatosensory cortex are associated with haemodynamic fluctuations that resemble stimulus-evoked haemodynamics

Although promise exists for patterns of resting-state blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) brain connectivity to be used as biomarkers of early brain pathology, a full understanding of the nature of the relationship between neural activity and spontaneous fMRI BOLD fluctuations is required before such data can be correctly interpreted. To investigate this issue, we combined electrophysiological recordings of rapid changes in multi-laminar local field potentials from the somatosensory cortex of anaesthetized rats with concurrent two-dimensional optical imaging spectroscopy measurements of resting-state haemodynamics that underlie fluctuations in the BOLD fMRI signal. After neural ‘events’ were identified, their time points served to indicate the start of an epoch in the accompanying haemodynamic fluctuations. Multiple epochs for both neural ‘events’ and the accompanying haemodynamic fluctuations were averaged. We found that the averaged epochs of resting-state haemodynamic fluctuations taken after neural ‘events’ closely resembled the temporal profile of stimulus-evoked cortical haemodynamics. Furthermore, we were able to demonstrate that averaged epochs of resting-state haemodynamic fluctuations resembling the temporal profile of stimulus-evoked haemodynamics could also be found after peaks in neural activity filtered into specific electroencephalographic frequency bands (theta, alpha, beta, and gamma). This technique allows investigation of resting-state neurovascular coupling using methodologies that are directly comparable to that developed for investigating stimulus-evoked neurovascular responses.

Modulation of visually evoked cortical fMRI responses by phase of ongoing occipital alpha oscillations

Using simultaneous electroencephalography as a measure of ongoing activity and functional magnetic resonance imaging (fMRI) as a measure of the stimulus-driven neural response, we examined whether the amplitude and phase of occipital alpha oscillations at the onset of a brief visual stimulus affects the amplitude of the visually evoked fMRI response. When accounting for intrinsic coupling of alpha amplitude and occipital fMRI signal by modeling and subtracting pseudo-trials, no significant effect of prestimulus alpha amplitude on the evoked fMRI response could be demonstrated. Regarding the effect of alpha phase, we found that stimuli arriving at the peak of the alpha cycle yielded a lower blood oxygenation level-dependent (BOLD) fMRI response in early visual cortex (V1/V2) than stimuli presented at the trough of the cycle. Our results therefore show that phase of occipital alpha oscillations impacts the overall strength of a visually evoked response, as indexed by the BOLD signal. This observation complements existing evidence that alpha oscillations reflect periodic variations in cortical excitability and suggests that the phase of oscillations in postsynaptic potentials can serve as a mechanism of gain control for incoming neural activity. Finally, our findings provide a putative neural basis for observations of alpha phase dependence of visual perceptual performance.

The BOLD response during Stroop task-like inhibition paradigms: effects of task difficulty and task-relevant modality

Previous studies of the Stroop task propose two key mediators: the prefrontal and cingulate cortices but hints exist of functional specialization within these regions. This study aimed to examine the effect of task modality upon the prefrontal and cingulate response by examining the response to colour, number, and shape Stroop tasks whilst BOLD fMRI images were acquired on a Siemens 3 T MRI scanner. Behavioural analyses indicated facilitation and interference effects and a noticeable effect of task difficulty. Some modular effects of modality were observed in the prefrontal cortex that survived exclusion of task difficulty related activations. No effect of task-relevant information was observed in the anterior cingulate. Future comparisons of the mediation of selective attention need to consider the effects of task context and task difficulty. (c) 2005 Elsevier Inc. All rights reserved.

A dynamic causal model of the coupling between pulse stimulation and neural activity

We present a dynamic causal model that can explain context-dependent changes in neural responses, in the rat barrel cortex, to an electrical whisker stimulation at different frequencies. Neural responses were measured in terms of local field potentials. These were converted into current source density (CSD) data, and the time series of the CSD sink was extracted to provide a time series response train. The model structure consists of three layers (approximating the responses from the brain stem to the thalamus and then the barrel cortex), and the latter two layers contain nonlinearly coupled modules of linear second-order dynamic systems. The interaction of these modules forms a nonlinear regulatory system that determines the temporal structure of the neural response amplitude for the thalamic and cortical layers. The model is based on the measured population dynamics of neurons rather than the dynamics of a single neuron and was evaluated against CSD data from experiments with varying stimulation frequency (1–40 Hz), random pulse trains, and awake and anesthetized animals. The model parameters obtained by optimization for different physiological conditions (anesthetized or awake) were significantly different. Following Friston, Mechelli, Turner, and Price (2000), this work is part of a formal mathematical system currently being developed (Zheng et al., 2005) that links stimulation to the blood oxygen level dependent (BOLD) functional magnetic resonance imaging (fMRI) signal through neural activity and hemodynamic variables. The importance of the model described here is that it can be used to invert the hemodynamic measurements of changes in blood flow to estimate the underlying neural activity.

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.

The neural basis of maternal responsiveness to infants: an fMRI study

Using fMRI, we examined the neural correlates of maternal responsiveness. Ten healthy mothers viewed alternating blocks of video: (i) 40 s of their own infant; (ii) 20 s of a neutral video; (iii) 40 s of an unknown infant and (iv) 20 s of neutral video, repeated 4 times. Predominant BOLD signal change to the contrast of infants minus neutral stimulus occurred in bilateral visual processing regions BA minus neutral stimulus occurred in bilateral visual processing regions (BA 38), left amygdala and visual cortex (BA 19), and to the unknown infant minus own infant contrast in bilateral orbitofrontal cortex (BA 10,47) and medial prefrontal cortex (BA 8). These findings suggest that amygdala and temporal pole may be key sites in mediating a mother's response to her infant and reaffirms their importance in face emotion processing and social behaviour.

Connecting mean field models of neural activity to EEG and fMRI data

Progress in functional neuroimaging of the brain increasingly relies on the integration of data from complementary imaging modalities in order to improve spatiotemporal resolution and interpretability. However, the usefulness of merely statistical combinations is limited, since neural signal sources differ between modalities and are related non-trivially. We demonstrate here that a mean field model of brain activity can simultaneously predict EEG and fMRI BOLD with proper signal generation and expression. Simulations are shown using a realistic head model based on structural MRI, which includes both dense short-range background connectivity and long-range specific connectivity between brain regions. The distribution of modeled neural masses is comparable to the spatial resolution of fMRI BOLD, and the temporal resolution of the modeled dynamics, importantly including activity conduction, matches the fastest known EEG phenomena. The creation of a cortical mean field model with anatomically sound geometry, extensive connectivity, and proper signal expression is an important first step towards the model-based integration of multimodal neuroimages.

Realistic mean field forward predictions for the integration of co-registered EEG/fMRI

Brain activity can be measured non-invasively with functional imaging techniques. Each pixel in such an image represents a neural mass of about 105 to 107 neurons. Mean field models (MFMs) approximate their activity by averaging out neural variability while retaining salient underlying features, like neurotransmitter kinetics. However, MFMs incorporating the regional variability, realistic geometry and connectivity of cortex have so far appeared intractable. This lack of biological realism has led to a focus on gross temporal features of the EEG. We address these impediments and showcase a "proof of principle" forward prediction of co-registered EEG/fMRI for a full-size human cortex in a realistic head model with anatomical connectivity, see figure 1. MFMs usually assume homogeneous neural masses, isotropic long-range connectivity and simplistic signal expression to allow rapid computation with partial differential equations. But these approximations are insufficient in particular for the high spatial resolution obtained with fMRI, since different cortical areas vary in their architectonic and dynamical properties, have complex connectivity, and can contribute non-trivially to the measured signal. Our code instead supports the local variation of model parameters and freely chosen connectivity for many thousand triangulation nodes spanning a cortical surface extracted from structural MRI. This allows the introduction of realistic anatomical and physiological parameters for cortical areas and their connectivity, including both intra- and inter-area connections. Proper cortical folding and conduction through a realistic head model is then added to obtain accurate signal expression for a comparison to experimental data. To showcase the synergy of these computational developments, we predict simultaneously EEG and fMRI BOLD responses by adding an established model for neurovascular coupling and convolving "Balloon-Windkessel" hemodynamics. We also incorporate regional connectivity extracted from the CoCoMac database [1]. Importantly, these extensions can be easily adapted according to future insights and data. Furthermore, while our own simulation is based on one specific MFM [2], the computational framework is general and can be applied to models favored by the user. Finally, we provide a brief outlook on improving the integration of multi-modal imaging data through iterative fits of a single underlying MFM in this realistic simulation framework.

Satiation attenuates BOLD activity in brain regions involved in reward and increases activity in dorsolateral prefrontal cortex: an fMRI study in healthy volunteers

BACKGROUND: Neural responses to rewarding food cues are significantly different in the fed vs. fasted (>8 h food-deprived) state. However, the effect of eating to satiety after a shorter (more natural) intermeal interval on neural responses to both rewarding and aversive cues has not been examined. OBJECTIVE: With the use of a novel functional magnetic resonance imaging (fMRI) task, we investigated the effect of satiation on neural responses to both rewarding and aversive food tastes and pictures. DESIGN: Sixteen healthy participants (8 men, 8 women) were scanned on 2 separate test days, before and after eating a meal to satiation or after not eating for 4 h (satiated vs. premeal). fMRI blood oxygen level-dependent (BOLD) signals to the sight and/or taste of the stimuli were recorded. RESULTS: A whole-brain cluster-corrected analysis (P < 0.05) showed that satiation attenuated the BOLD response to both stimulus types in the ventromedial prefrontal cortex (vmPFC), orbitofrontal cortex, nucleus accumbens, hypothalamus, and insula but increased BOLD activity in the dorsolateral prefrontal cortex (dlPFC; local maxima corrected to P ≤ 0.001). A psychophysiological interaction analysis showed that the vmPFC was more highly connected to the dlPFC when individuals were exposed to food stimuli when satiated than when not satiated. CONCLUSIONS: These results suggest that natural satiation attenuates activity in reward-related brain regions and increases activity in the dlPFC, which may reflect a "top down" cognitive influence on satiation. This trial was registered at clinicaltrials.gov as NCT02298049.

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.

Dynamic causal modeling applied to fMRI data shows high reliability

Sensitivity, specificity, and reproducibility are vital to interpret neuroscientific results from functional magnetic resonance imaging (fMRI) experiments. Here we examine the scan–rescan reliability of the percent signal change (PSC) and parameters estimated using Dynamic Causal Modeling (DCM) in scans taken in the same scan session, less than 5 min apart. We find fair to good reliability of PSC in regions that are involved with the task, and fair to excellent reliability with DCM. Also, the DCM analysis uncovers group differences that were not present in the analysis of PSC, which implies that DCM may be more sensitive to the nuances of signal changes in fMRI data.

fMRI delineation of working memory for emotional prosody in the brain: commonalities with the lexico-semantic emotion network

Decoding emotional prosody is crucial for successful social interactions, and continuous monitoring of emotional intent via prosody requires working memory. It has been proposed by Ross and others that emotional prosody cognitions in the right hemisphere are organized in an analogous fashion to propositional language functions in the left hemisphere. This study aimed to test the applicability of this model in the context of prefrontal cortex working memory functions. BOLD response data were therefore collected during performance of two emotional working memory tasks by participants undergoing fMRI. In the prosody task, participants identified the emotion conveyed in pre-recorded sentences, and working memory load was manipulated in the style of an N-back task. In the matched lexico-semantic task, participants identified the emotion conveyed by sentence content. Block-design neuroimaging data were analyzed parametrically with SPM5. At first, working memory for emotional prosody appeared to be right-lateralized in the PFC, however, further analyses revealed that it shared much bilateral prefrontal functional neuroanatomy with working memory for lexico-semantic emotion. Supplementary separate analyses of males and females suggested that these language functions were less bilateral in females, but their inclusion did not alter the direction of laterality. It is concluded that Ross et al.'s model is not applicable to prefrontal cortex working memory functions, that evidence that working memory cannot be subdivided in prefrontal cortex according to material type is increased, and that incidental working memory demands may explain the frontal lobe involvement in emotional prosody comprehension as revealed by neuroimaging studies. (c) 2007 Elsevier Inc. All rights reserved.

The contribution of primary and secondary somatosensory cortices to the representation of body parts and body sides: an fMRI adaptation study

Although the somatosensory homunculus is a classically used description of the way somatosensory inputs are processed in the brain, the actual contributions of primary (SI) and secondary (SII) somatosensory cortices to the spatial coding of touch remain poorly understood. We studied adaptation of the fMRI BOLD response in the somatosensory cortex by delivering pairs of vibrotactile stimuli to the finger tips of the index and middle fingers. The first stimulus (adaptor) was delivered either to the index or to the middle finger of the right or left hand, whereas the second stimulus (test) was always administered to the left index finger. The overall BOLD response evoked by the stimulation was primarily contralateral in SI and was more bilateral in SII. However, our fMRI adaptation approach also revealed that both somatosensory cortices were sensitive to ipsilateral as well as to contralateral inputs. SI and SII adapted more after subsequent stimulation of homologous as compared with nonhomologous fingers, showing a distinction between different fingers. Most importantly, for both somatosensory cortices, this finger-specific adaptation occurred irrespective of whether the tactile stimulus was delivered to the same or to different hands. This result implies integration of contralateral and ipsilateral somatosensory inputs in SI as well as in SII. Our findings suggest that SI is more than a simple relay for sensory information and that both SI and SII contribute to the spatial coding of touch by discriminating between body parts (fingers) and by integrating the somatosensory input from the two sides of the body (hands).

Epoch-specific functional networks involved in working memory

Working memory (WM) is not a unitary construct. There are distinct processes involved in encoding information, maintaining it on-line, and using it to guide responses. The anatomical configurations of these processes are more accurately analyzed as functionally connected networks than collections of individual regions. In the current study we analyzed event-related functional magnetic resonance imaging (fMRI) data from a Sternberg Item Recognition Paradigm WM task using a multivariate analysis method that allowed the linking of functional networks to temporally-separated WM epochs. The length of the delay epochs was varied to optimize isolation of the hemodynamic response (HDR) for each task epoch. All extracted functional networks displayed statistically significant sensitivity to delay length. Novel information extracted from these networks that was not apparent in the univariate analysis of these data included involvement of the hippocampus in encoding/probe, and decreases in BOLD signal in the superior temporal gyrus (STG), along with default-mode regions, during encoding/delay. The bilateral hippocampal activity during encoding/delay fits with theoretical models of WM in which memoranda held across the short term are activated long-term memory representations. The BOLD signal decreases in the STG were unexpected, and may reflect repetition suppression effects invoked by internal repetition of letter stimuli. Thus, analysis methods focusing on how network dynamics relate to experimental conditions allowed extraction of novel information not apparent in univariate analyses, and are particularly recommended for WM experiments for which task epochs cannot be randomized.