Negative BOLD fMRI response in the visual cortex carries precise stimulus-specific information

Sustained positive BOLD (blood oxygen level-dependent) activity is employed extensively in functional magnetic resonance imaging (fMRI) studies as evidence for task or stimulus-specific neural responses. However, the presence of sustained negative BOLD activity (i.e., sustained responses that are lower than the fixation baseline) has remained more difficult to interpret. Some studies suggest that it results from local "blood stealing" wherein blood is diverted to neurally active regions without a concomitant change of neural activity in the negative BOLD regions. However, other evidence suggests that negative BOLD is a result of local neural suppression. In both cases, regions of negative BOLD response are usually interpreted as carrying relatively little, if any, stimulus-specific information (hence the predominant reliance on positive BOLD activity in fMRI). Here we show that the negative BOLD response resulting from visual stimulation can carry high information content that is stimulus-specific. Using a general linear model (GLM), we contrasted standard flickering stimuli to a fixation baseline and found regions of the visual cortex that displayed a sustained negative BOLD response, consistent with several previous studies. Within these negative BOLD regions, we compared patterns of fMRI activity generated by flickering Gabors that were systematically shifted in position. As the Gabors were shifted further from each other, the correlation in the spatial pattern of activity across a population of voxels (such as the population of V1 voxels that displayed a negative BOLD response) decreased significantly. Despite the fact that the BOLD signal was significantly negative (lower than fixation baseline), these regions were able to discriminate objects separated by less than 0.5 deg (at approximately 10 deg eccentricity). The results suggest that meaningful, stimulus-specific processing occurs even in regions that display a strong negative BOLD response.     

Sex differences in the cerebral BOLD signal response to painful heat stimuli

There are limited data addressing the question of sex differences in pain-related cerebral processing. This study examined whether pain-related blood oxygenation level-dependent (BOLD) signal change measured with functional magnetic resonance imaging (fMRI) demonstrated sex differences, under conditions of equivalent pain perception. Twenty-eight healthy volunteers (17 women, 11 men) were subject to a fMRI scan while noxious heat stimuli were applied to the dorsum of the left foot. Significant BOLD signal modulation was observed in several nociceptive processing regions of interest (ROIs) in all subjects. There were no sex differences in the spatial extent of BOLD signal change for any ROI, but the signal amplitude was lower for women in most ROIs and significantly so for the primary somatosensory cortex (S1), the midanterior cingulate cortex, and the dorsolateral prefrontal cortex (DLPFC). The BOLD signal response could be positive or negative, and frequently, both polarities were observed within a single ROI. In most ROIs, women show proportionately more voxels with negative signal change than men, and this difference was statistically significant for the S1 and the DLPFC. The time course of the negative signal change was very similar to that of the positive signal change, suggesting that the latter was not "driving" the former. The location of negative and positive clusters formed distinct patterns in several of the ROIs, and these patterns suggest something other than a local "steal" phenomenon as an explanation for the negative signal changes. Sex differences in baseline cerebral blood flow may contribute to the BOLD signal differences observed in this study.     

Positive Allosteric Modulator of GABA Lowers BOLD Responses in the Cingulate Cortex

Knowledge about the neural underpinnings of the negative blood oxygen level dependent (BOLD) responses in functional magnetic resonance imaging (fMRI) is still limited. We hypothesized that pharmacological GABAergic modulation attenuates BOLD responses, and that blood concentrations of a positive allosteric modulator of GABA correlate inversely with BOLD responses in the cingulate cortex. We investigated whether or not pure task-related negative BOLD responses were co-localized with pharmacologically modulated BOLD responses. Twenty healthy adults received either 5 mg diazepam or placebo in a double blind, randomized design. During fMRI the subjects performed a working memory task. Results showed that BOLD responses in the cingulate cortex were inversely correlated with diazepam blood concentrations; that is, the higher the blood diazepam concentration, the lower the BOLD response. This inverse correlation was most pronounced in the pregenual anterior cingulate cortex and the anterior mid-cingulate cortex. For subjects with diazepam plasma concentration > 0.1 mg/L we observed negative BOLD responses with respect to fixation baseline. There was minor overlap between cingulate regions with task-related negative BOLD responses and regions where the BOLD responses were inversely correlated with diazepam concentration. We interpret that the inverse correlation between the BOLD response and diazepam was caused by GABA-related neural inhibition. Thus, this study supports the hypothesis that GABA attenuates BOLD responses in fMRI. The minimal overlap between task-related negative BOLD responses and responses attenuated by diazepam suggests that these responses might be caused by different mechanisms.     

Sustained negative BOLD,blood flow and oxygen consumption response and its coupling to the positive response in the human brain

Most fMRI studies are based on the detection of a positive BOLD response (PBR). Here, we demonstrate and characterize a robust sustained negative BOLD response (NBR) in the human occipital cortex, triggered by stimulating part of the visual field. The NBR was spatially adjacent to but segregated from the PBR. It depended on the stimulus and thus on the pattern of neuronal activity. The time courses of the NBR and PBR were similar, and their amplitudes covaried both with increasing stimulus duration and increasing stimulus contrast. The NBR was associated with reductions in blood flow and with decreases in oxygen consumption. Our findings support the contribution to the NBR of (1) a significant component of reduction in neuronal activity and (2) possibly a component of hemodynamic changes independent of the local changes in neuronal activity.     

The Statistical Analysis of Multi-Voxel Patterns in Functional Imaging

The goal of multi-voxel pattern analysis (MVPA) in BOLD imaging is to determine whether patterns of activation across multiple voxels change with experimental conditions. MVPA is a powerful technique, its use is rapidly growing, but it poses serious statistical challenges. For instance, it is well-known that the slow nature of the BOLD response can lead to greatly exaggerated performance estimates. Methods are available to avoid this overestimation, and we present those here in tutorial fashion. We go on to show that, even with these methods, standard tests of significance such as Students’ T and the binomial tests are invalid in typical MRI experiments. Only a carefully constructed permutation test correctly assesses statistical significance. Furthermore, our simulations show that performance estimates increase with both temporal as well as spatial signal correlations among multiple voxels. This dependence implies that a comparison of MVPA performance between areas, between subjects, or even between BOLD signals that have been preprocessed in different ways needs great care.     

Multi-Regional Investigation of the Relationship between Functional MRI Blood Oxygenation Level Dependent (BOLD) Activation and GABA Concentration

Several recent studies have reported an inter-individual correlation between regional GABA concentration, as measured by MRS, and the amplitude of the functional blood oxygenation level dependent (BOLD) response in the same region. In this study, we set out to investigate whether this coupling generalizes across cortex. In 18 healthy participants, we performed edited MRS measurements of GABA and BOLD-fMRI experiments using regionally related activation paradigms. Regions and tasks were the: occipital cortex with a visual grating stimulus; auditory cortex with a white noise stimulus; sensorimotor cortex with a finger-tapping task; frontal eye field with a saccade task; and dorsolateral prefrontal cortex with a working memory task. In contrast to the prior literature, no correlation between GABA concentration and BOLD activation was detected in any region. The origin of this discrepancy is not clear. Subtle differences in study design or insufficient power may cause differing results; these and other potential reasons for the discrepant results are discussed. This negative result, although it should be interpreted with caution, has a larger sample size than prior positive results, and suggests that the relationship between GABA and the BOLD response may be more complex than previously thought.     

Action Processing and Mirror Neuron Function in Patients with Amyotrophic Lateral Sclerosis: An fMRI Study

Amyotrophic lateral sclerosis (ALS) is a highly debilitating and rapidly fatal neurodegenerative disease. It has been suggested that social cognition may be affected, such as impairment in theory of mind (ToM) ability. Despite these findings, research in this area is scarce and the investigation of neural mechanisms behind such impairment is absent. Nineteen patients with ALS and eighteen healthy controls participated in this study. Because the mirror neuron system (MNS) is thought to be involved in theory of mind, we first implemented a straightforward action-execution and observation task to assess basic MNS function. Second, we examined the social-cognitive ability to understand actions of others, which is a component of ToM. We used fMRI to assess BOLD activity differences between groups during both experiments. Theory of mind was also measured behaviorally using the Reading the Mind in the Eyes test (RME). ALS patients displayed greater BOLD activity during the action-execution and observation task, especially throughout right anterior cortical regions. These areas included the right inferior operculum, premotor and primary motor regions, and left inferior parietal lobe. A conjunction analysis showed significantly more co-activated voxels during both the observation and action-execution conditions in the patient group throughout MNS regions. These results support a compensatory response in the MNS during action processing. In the action understanding experiment, healthy controls performed better behaviorally and subsequently recruited greater regions of activity throughout the prefrontal cortex and middle temporal gyrus. Lastly, action understanding performance was able to cluster patients with ALS into high and lower performing groups, which then differentiated RME performance. Collectively, these data suggest that social cognition, particularly theory of mind, may be affected in a subset of patients with ALS. This impairment may be related to functioning of the MNS and other regions related to action processing and understanding. Implications for future research are discussed.     

From Blood Oxygenation Level Dependent (BOLD) Signals to Brain Temperature Maps

A theoretical framework is presented for converting Blood Oxygenation Level Dependent (BOLD) images to brain temperature maps, based on the idea that disproportional local changes in cerebral blood flow (CBF) as compared with cerebral metabolic rate of oxygen consumption (CMRO ₂) during functional brain activity, lead to both brain temperature changes and the BOLD effect. Using an oxygen limitation model and a BOLD signal model, we obtain a transcendental equation relating CBF and CMRO ₂ changes with the corresponding BOLD signal, which is solved in terms of the Lambert W function. Inserting this result in the dynamic bioheat equation describing the rate of temperature changes in the brain, we obtain a nonautonomous ordinary differential equation that depends on the BOLD response, which is solved numerically for each brain voxel. Temperature maps obtained from a real BOLD dataset registered in an attention to visual motion experiment were calculated, obtaining temperature variations in the range: (−0.15, 0.1) which is consistent with experimental results. The statistical analysis revealed that significant temperature activations have a similar distribution pattern than BOLD activations. An interesting difference was the activation of the precuneus in temperature maps, a region involved in visuospatial processing, an effect that was not observed on BOLD maps. Furthermore, temperature maps were more localized to gray matter regions than the original BOLD maps, showing less activated voxels in white matter and cerebrospinal fluid.     

Reproducibility of swallow-induced cortical BOLD positive and negative fMRI activity

Functional MRI (fMRI) studies have demonstrated that a number of brain regions (cingulate, insula, prefrontal, and sensory/motor cortices) display blood oxygen level-dependent (BOLD) positive activity during swallow. Negative BOLD activations and reproducibility of these activations have not been systematically studied. The aim of our study was to investigate the reproducibility of swallow-related cortical positive and negative BOLD activity across different fMRI sessions. We studied 16 healthy volunteers utilizing an fMRI event-related analysis. Individual analysis using a general linear model was used to remove undesirable signal changes correlated with motion, white matter, and cerebrospinal fluid. The group analysis used a mixed-effects multilevel model to identify active cortical regions. The volume and magnitude of a BOLD signal within each cluster was compared between the two study sessions. All subjects showed significant clustered BOLD activity within the known areas of cortical swallowing network across both sessions. The cross-correlation coefficient of percent fMRI signal change and the number of activated voxels across both positive and negative BOLD networks were similar between the two studies (r ≥ 0.87, P < 0.0001). Swallow-associated negative BOLD activity was comparable to the well-defined "default-mode" network, and positive BOLD activity had noticeable overlap with the previously described "task-positive" network. Swallow activates two parallel cortical networks. These include a positive and a negative BOLD network, respectively, correlated and anticorrelated with swallow stimulus. Group cortical activity maps, as well as extent and amplitude of activity induced by volitional swallowing in the cortical swallowing network, are reproducible between study sessions.     

Trait anxiety and neural efficiency of abstract reasoning: An fMRI investigation

Worries preoccupy the working memory capacity in anxious individuals, thereby affecting their performance during tasks that require efficient attention regulation. According to the attentional control theory (ACT), trait anxiety affects the processing efficiency, i.e. the effort required for task performance, more than the accuracy of task performance. We investigated the relation between trait anxiety and neural response for a reasoning task in healthy subjects. Functional magnetic resonance imaging (fMRI) was carried out on 22 healthy participants and blood oxygenation level dependent (BOLD) contrast estimates were extracted from a priori regions of interest (ROIs) that were earlier implicated in reasoning (i.e., bilaterally caudate head, globus pallidus, thalamus, prefrontal cortex [rostral, dorsal and ventral regions], inferior parietal lobule and middle occipital gyrus). Controlling for the effects of age, gender, state anxiety and depressive symptoms, for equivalent levels of task performance, trait anxiety of the participants was found to be associated with an increase in task related BOLD activation in right globus pallidus, left thalamus and left middle occipital gyrus. Our results suggest a reduced processing efficiency for reasoning in high trait anxiety subjects and provides important brain–behaviour relationships with respect to sub-clinical anxiety.     

Visual Interactions Conform to Pattern Decorrelation in Multiple Cortical Areas

Neural responses to visual stimuli are strongest in the classical receptive field, but they are also modulated by stimuli in a much wider region. In the primary visual cortex, physiological data and models suggest that such contextual modulation is mediated by recurrent interactions between cortical areas. Outside the primary visual cortex, imaging data has shown qualitatively similar interactions. However, whether the mechanisms underlying these effects are similar in different areas has remained unclear. Here, we found that the blood oxygenation level dependent (BOLD) signal spreads over considerable cortical distances in the primary visual cortex, further than the classical receptive field. This indicates that the synaptic activity induced by a given stimulus occurs in a surprisingly extensive network. Correspondingly, we found suppressive and facilitative interactions far from the maximum retinotopic response. Next, we characterized the relationship between contextual modulation and correlation between two spatial activation patterns. Regardless of the functional area or retinotopic eccentricity, higher correlation between the center and surround response patterns was associated with stronger suppressive interaction. In individual voxels, suppressive interaction was predominant when the center and surround stimuli produced BOLD signals with the same sign. Facilitative interaction dominated in the voxels with opposite BOLD signal signs. Our data was in unison with recently published cortical decorrelation model, and was validated against alternative models, separately in different eccentricities and functional areas. Our study provides evidence that spatial interactions among neural populations involve decorrelation of macroscopic neural activation patterns, and suggests that the basic design of the cerebral cortex houses a robust decorrelation mechanism for afferent synaptic input.     

Measuring Granger Causality between cortical regions from voxelwise fMRI BOLD signals with LASSO

Functional brain network studies using the Blood Oxygen-Level Dependent (BOLD) signal from functional Magnetic Resonance Imaging (fMRI) are becoming increasingly prevalent in research on the neural basis of human cognition. An important problem in functional brain network analysis is to understand directed functional interactions between brain regions during cognitive performance. This problem has important implications for understanding top-down influences from frontal and parietal control regions to visual occipital cortex in visuospatial attention, the goal motivating the present study. A common approach to measuring directed functional interactions between two brain regions is to first create nodal signals by averaging the BOLD signals of all the voxels in each region, and to then measure directed functional interactions between the nodal signals. Another approach, that avoids averaging, is to measure directed functional interactions between all pairwise combinations of voxels in the two regions. Here we employ an alternative approach that avoids the drawbacks of both averaging and pairwise voxel measures. In this approach, we first use the Least Absolute Shrinkage Selection Operator (LASSO) to pre-select voxels for analysis, then compute a Multivariate Vector AutoRegressive (MVAR) model from the time series of the selected voxels, and finally compute summary Granger Causality (GC) statistics from the model to represent directed interregional interactions. We demonstrate the effectiveness of this approach on both simulated and empirical fMRI data. We also show that averaging regional BOLD activity to create a nodal signal may lead to biased GC estimation of directed interregional interactions. The approach presented here makes it feasible to compute GC between brain regions without the need for averaging. Our results suggest that in the analysis of functional brain networks, careful consideration must be given to the way that network nodes and edges are defined because those definitions may have important implications for the validity of the analysis.     

Observations on the pyramidal tract of the slow loris

The pyramidal system of the slow loris was investigated physiologically and histologically. The cerebral response to medullary pyramid stimulation was mapped and found to distribute in two separate foci, rostral to the hand and foot somatosensory foci. The tract was found to contain about 204,000 fibers at the midolivary level on one side, with about 85% of the fibers less than 2 Μm and 1% greater than 4 Μm. Calculation of the amplitude and configuration of the true antidromic response, based on the histological findings, disagreed with the observed response. The latter was found to be analogous to the rwave of the American opossum, which appears to result from intracortical pyramidal tract collateral activation of the cerebral neurons. This response is not present in the macaque monkey.     

Imaging the spatio-temporal dynamics of supragranular activity in the rat somatosensory cortex in response to stimulation of the paws

We employed voltage-sensitive dye (VSD) imaging to investigate the spatio-temporal dynamics of the responses of the supragranular somatosensory cortex to stimulation of the four paws in urethane-anesthetized rats. We obtained the following main results. (1) Stimulation of the contralateral forepaw evoked VSD responses with greater amplitude and smaller latency than stimulation of the contralateral hindpaw, and ipsilateral VSD responses had a lower amplitude and greater latency than contralateral responses. (2) While the contralateral stimulation initially activated only one focus, the ipsilateral stimulation initially activated two foci: one focus was typically medial to the focus activated by contralateral stimulation and was stereotaxically localized in the motor cortex; the other focus was typically posterior to the focus activated by contralateral stimulation and was stereotaxically localized in the somatosensory cortex. (3) Forepaw and hindpaw somatosensory stimuli activated large areas of the sensorimotor cortex, well beyond the forepaw and hindpaw somatosensory areas of classical somatotopic maps, and forepaw stimuli activated larger cortical areas with greater activation velocity than hindpaw stimuli. (4) Stimulation of the forepaw and hindpaw evoked different cortical activation dynamics: forepaw responses displayed a clear medial directionality, whereas hindpaw responses were much more uniform in all directions. In conclusion, this work offers a complete spatio-temporal map of the supragranular VSD cortical activation in response to stimulation of the paws, showing important somatotopic differences between contralateral and ipsilateral maps as well as differences in the spatio-temporal activation dynamics in response to forepaw and hindpaw stimuli.     

Amphetamine Sensitization Alters Reward Processing in the Human Striatum and Amygdala

Dysregulation of mesolimbic dopamine transmission is implicated in a number of psychiatric illnesses characterised by disruption of reward processing and goal-directed behaviour, including schizophrenia, drug addiction and impulse control disorders associated with chronic use of dopamine agonists. Amphetamine sensitization (AS) has been proposed to model the development of this aberrant dopamine signalling and the subsequent dysregulation of incentive motivational processes. However, in humans the effects of AS on the dopamine-sensitive neural circuitry associated with reward processing remains unclear. Here we describe the effects of acute amphetamine administration, following a sensitising dosage regime, on blood oxygen level dependent (BOLD) signal in dopaminoceptive brain regions during a rewarded gambling task performed by healthy volunteers. Using a randomised, double-blind, parallel-groups design, we found clear evidence for sensitization to the subjective effects of the drug, while rewarded reaction times were unchanged. Repeated amphetamine exposure was associated with reduced dorsal striatal BOLD signal during decision making, but enhanced ventromedial caudate activity during reward anticipation. The amygdala BOLD response to reward outcomes was blunted following repeated amphetamine exposure. Positive correlations between subjective sensitization and changes in anticipation- and outcome-related BOLD signal were seen for the caudate nucleus and amygdala, respectively. These data show for the first time in humans that AS changes the functional impact of acute stimulant exposure on the processing of reward-related information within dopaminoceptive regions. Our findings accord with pathophysiological models which implicate aberrant dopaminergic modulation of striatal and amygdala activity in psychosis and drug-related compulsive disorders.     

Cognitive Control in Auditory Working Memory Is Enhanced in Musicians

Musical competence may confer cognitive advantages that extend beyond processing of familiar musical sounds. Behavioural evidence indicates a general enhancement of both working memory and attention in musicians. It is possible that musicians, due to their training, are better able to maintain focus on task-relevant stimuli, a skill which is crucial to working memory. We measured the blood oxygenation-level dependent (BOLD) activation signal in musicians and non-musicians during working memory of musical sounds to determine the relation among performance, musical competence and generally enhanced cognition. All participants easily distinguished the stimuli. We tested the hypothesis that musicians nonetheless would perform better, and that differential brain activity would mainly be present in cortical areas involved in cognitive control such as the lateral prefrontal cortex. The musicians performed better as reflected in reaction times and error rates. Musicians also had larger BOLD responses than non-musicians in neuronal networks that sustain attention and cognitive control, including regions of the lateral prefrontal cortex, lateral parietal cortex, insula, and putamen in the right hemisphere, and bilaterally in the posterior dorsal prefrontal cortex and anterior cingulate gyrus. The relationship between the task performance and the magnitude of the BOLD response was more positive in musicians than in non-musicians, particularly during the most difficult working memory task. The results confirm previous findings that neural activity increases during enhanced working memory performance. The results also suggest that superior working memory task performance in musicians rely on an enhanced ability to exert sustained cognitive control. This cognitive benefit in musicians may be a consequence of focused musical training.     

Echo-time and field strength dependence of BOLD reactivity in veins and parenchyma using flow-normalized hypercapnic manipulation

While the BOLD (Blood Oxygenation Level Dependent) contrast mechanism has demonstrated excellent sensitivity to neuronal activation, its specificity with regards to differentiating vascular and parenchymal responses has been an area of ongoing concern. By inducing a global increase in Cerebral Blood Flow (CBF), we examined the effect of magnetic field strength and echo-time (TE) on the gradient-echo BOLD response in areas of cortical gray matter and in resolvable veins. In order to define a quantitative index of BOLD reactivity, we measured the percent BOLD response per unit fractional change in global gray matter CBF induced by inhaling carbon dioxide (CO(2)). By normalizing the BOLD response to the underlying CBF change and determining the BOLD response as a function of TE, we calculated the change in R(2)(*) (ΔR(2)(*)) per unit fractional flow change; the Flow Relaxation Coefficient, (FRC) for 3T and 1.5T in parenchymal and large vein compartments. The FRC in parenchymal voxels was 1.76±0.54 fold higher at 3T than at 1.5T and was 2.96±0.66 and 3.12±0.76 fold higher for veins than parenchyma at 1.5T and 3T respectively, showing a quantitative measure of the increase in specificity to parenchymal sources at 3T compared to 1.5T. Additionally, the results allow optimization of the TE to prioritize either maximum parenchymal BOLD response or maximum parenchymal specificity. Parenchymal signals peaked at TE values of 62.0±11.5 ms and 41.5±7.5 ms for 1.5T and 3T, respectively, while the response in the major veins peaked at shorter TE values; 41.0±6.9 ms and 21.5±1.0 ms for 1.5T and 3T. These experiments showed that at 3T, the BOLD CNR in parenchymal voxels exceeded that of 1.5T by a factor of 1.9±0.4 at the optimal TE for each field.     

BOLD signals correlate with ensemble unit activities in rat's somatosensory cortex

Evoked neural activity (ensemble single-unit activity and evoked field potential) and functional magnetic resonance imaging (fMRI) changes of the primary somatosensory cortex in response to electrical stimulation of the hind paw were studied in rats under anesthesia. The effects of stimulation frequency (ranging from 0.3 to 10 Hz) and types of anesthetics (alpha-chloralose and sodium pentobarbital) on blood oxygen level dependent (BOLD) activation and neural activation were compared. Both ensemble single-unit activity and BOLD signal changes achieved maximal activation at 3 Hz of stimulation and responses were significantly stronger under alpha-chloralose anesthesia. The maximal activation of the integral evoked potential (sigmaEP), in contrast, was the highest at 10 Hz; and the values were similar for alpha-chloralose and pentobarbital. These analyses revealed that fMRI image changes were better correlated with ensemble single-unit activity than with sigmaEP during somatosensory stimulations.     

Physiological Changes during Recovery from a Primate Dorsal Column Lesion

The evoked potential (EP) over primary somatosensory cortex (SI) was monitored before and after a complete lesion of the primate dorsal column (DC) pathway on one side. The EP was elicited by electrocutaneous or mechanical stimulation of either foot, and was recorded from the contralateral cortical surface for periods of up to 3 months after the lesion. The amplitudes of the three major peaks (P20, N50, and P90) of the cortical somatosensory EP were significantly reduced following interruption of the contralateral DC. Over weeks following the lesion, there was a significant increase in amplitude of the P90 component of the EP that was not evident in the other peaks. The postlesion increases in P90 amplitude were correlated with improved performance on a task that required grasping with either foot, suggesting that behavioral recovery from a DC lesion results in part from neural plasticity, as opposed to a simple relearning of the task.