The BOLD plot thickens: sign- and layer-dependent hemodynamic changes with activation

GABAergic inhibition shapes behavior and neural dynamics in human visual working memory

Neuronal inhibition, primarily mediated by GABAergic neurotransmission, is crucial for brain development and healthy cognition. Gamma-aminobutyric acid concentration levels in sensory areas have been shown to correlate with hemodynamic and oscillatory neuronal responses. How these measures relate to one another during working memory, a higher-order cognitive process, is still poorly understood. We address this gap by collecting magnetoencephalography, functional magnetic resonance imaging, and Flumazenil positron emission tomography data within the same subject cohort using an n-back working-memory paradigm. By probing the relationship between GABAA receptor distribution, neural oscillations, and Blood Oxygen Level Dependent (BOLD) modulations, we found that GABAA receptor density in higher-order cortical areas predicted the reaction times on the working-memory task and correlated positively with the peak frequency of gamma power modulations and negatively with BOLD amplitude. These findings support and extend theories linking gamma oscillations and hemodynamic responses to gamma-aminobutyric acid neurotransmission and to the excitation-inhibition balance and cognitive performance in humans. Considering the small sample size of the study, future studies should test whether these findings also hold for other, larger cohorts as well as to examine in detail how the GABAergic system and neural fluctuations jointly support working-memory task performance.

Multi-Echo Investigations of Positive and Negative CBF and Concomitant BOLD Changes: Positive and negative CBF and BOLD changes

Unlike the positive blood oxygenation level-dependent (BOLD) response (PBR), commonly taken as an indication of an 'activated' brain region, the physiological origin of negative BOLD signal changes (i.e. a negative BOLD response, NBR), also referred to as 'deactivation' is still being debated. In this work, an attempt was made to gain a better understanding of the underlying mechanism by obtaining a comprehensive measure of the contributing cerebral blood flow (CBF) and its relationship to the NBR in the human visual cortex, in comparison to a simultaneously induced PBR in surrounding visual regions. To overcome the low signal-to-noise ratio (SNR) of CBF measurements, a newly developed multi-echo version of a center-out echo planar-imaging (EPI) readout was employed with pseudo-continuous arterial spin labeling (pCASL). It achieved very short echo and inter-echo times and facilitated a simultaneous detection of functional CBF and BOLD changes at 3 T with improved sensitivity. Evaluations of the absolute and relative changes of CBF and the effective transverse relaxation rate,R₂^* the coupling ratios, and their dependence on CBF at rest, CBF_rest indicated differences between activated and deactivated regions. Analysis of the shape of the respective functional responses also revealed faster negative responses with more pronounced post-stimulus transients. Resulting differences in the flow-metabolism coupling ratios were further examined for potential distinctions in the underlying neuronal contributions.

Weak Task Synchronization of Default Mode Network in Task Based Paradigms

The brains Default mode network (DMN) is generally characterized by brain areas that gets deactivated upon the presentation of a wide variety of externally focused, attention demanding tasks. These areas also exhibit significant intra-DMN functional connectivity and significant negative functional connectivity with other brain areas, especially with attention networks, in both resting state and task conditions. Therefore, the DMN has been hypothesized to be involved in internally directed cognitive activities such as autobiographical recall of the past, future planning and mind wandering. Recent research has discovered the role of bottom-up attention in modulating the behaviour of DMN. We hypothesize that the de-engagement of the DMN regions upon the presentation of an externally-focused attention-demanding stimulus may not be strictly stimulus locked and may exhibit significant trial-to-trial as well as subject-to-subject variability. Due to the involvement of frontoparietal control network in modulating the anticorrelations between DMN and dorsal attention network (DAN), we expect the DMN regions to have lower inter-trial and inter-subject synchronization in their fMRI BOLD responses as compared to the bottom-up early-sensory task-positive regions. To test this hypothesis, we designed new statistical methods called Inter Trial Temporal Synchronization Analysis (IT-TSA) and Inter Subject TSA (IS-TSA) to analyse variability across trials and subjects respectively. We analysed four publicly available datasets (total 223 subjects) across seven tasks related to different cognitive modalities and found out that there is significantly low stimulus-locked synchronization across trials and subjects in the DMN regions as compared to early sensory task positive regions. Our study challenges the understanding of DMN as a strictly task-negative region and supports the recent findings that DMN acts as an active component associated with intrinsic processing which deactivates differentially and non-linearly across trials and subjects in the presence of extrinsic processes.

Resting state fMRI connectivity is sensitive to laminar connectional architecture in the human brain

Previous invasive studies indicate that human neocortical graymatter contains cytoarchitectonically distinct layers, with notable differences in their structural connectivity with the rest of the brain. Given recent improvements in the spatial resolution of anatomical and functional magnetic resonance imaging (fMRI), we hypothesize that resting state functional connectivity (FC) derived from fMRI is sensitive to layer-specific thalamo-cortical and cortico-cortical microcircuits. Using sub-millimeter resting state fMRI data obtained at 7 T, we found that: (1) FC between the entire thalamus and cortical layers I and VI was significantly stronger than between the thalamus and other layers. Furthermore, FC between somatosensory thalamus (ventral posterolateral nucleus, VPL) and layers IV, VI of the primary somatosensory cortex were stronger than with other layers; (2) Inter-hemispheric cortico-cortical FC between homologous regions in superficial layers (layers I-III) was stronger compared to deep layers (layers V-VI). These findings are in agreement with structural connections inferred from previous invasive studies that showed that: (i) M-type neurons in the entire thalamus project to layer-I; (ii) Pyramidal neurons in layer-VI target all thalamic nuclei, (iii) C-type neurons in the VPL project to layer-IV and receive inputs from layer-VI of the primary somatosensory cortex, and (iv) 80% of collosal projecting neurons between homologous cortical regions connect superficial layers. Our results demonstrate for the first time that resting state fMRI is sensitive to structural connections between cortical layers (previously inferred through invasive studies), specifically in thalamo-cortical and cortico-cortical networks.

Tonic, Burst, and Burst-Cycle Spinal Cord Stimulation Lead to Differential Brain Activation Patterns as Detected by Functional Magnetic Resonance Imaging

OBJECTIVE

The objective of this preclinical study was to examine the responses of the brain to noxious stimulation in the presence and absence of different modes of spinal cord stimulation (SCS) using blood-oxygen-level-dependent functional magnetic resonance imaging (BOLD-fMRI).

MATERIALS AND METHODS

Sprague-Dawley rats were randomized to groups based on the mode of SCS delivered which included tonic stimulation (n = 27), burst stimulation (n = 30), and burst-cycle stimulation (n = 29). The control (sham) group (n = 28) received no SCS. The SCS electrode was inserted between T10 and T12 spinal levels prior to fMRI session. The experimental protocol for fMRI acquisition consisted of an initial noxious stimulation phase, a treatment phase wherein the SCS was turned on concurrently with noxious stimulation, and a residual effect phase wherein the noxious stimulation alone was turned on. The responses were statistically analyzed through paired t-test and the results were presented as z-scores for the quantitative analysis of the fMRI data.

RESULTS

The treatment with different SCS modes attenuated the BOLD brain responses to noxious hindlimb stimulation. The tonic, burst, and burst-cycle SCS treatment attenuated BOLD responses in the caudate putamen (CPu), insula (In), and secondary somatosensory cortex (S2). There was little to no corresponding change in sham control in these three regions. The burst and burst-cycle SCS demonstrated greater attenuation of BOLD signals in CPu, In, and S2 compared to tonic stimulation.

CONCLUSION

The high-resolution fMRI study using a rat model demonstrated the potential of different SCS modes to act on several pain-matrix-related regions of the brain in response to noxious stimulation. The burst and burst-cycle SCS exhibited greater brain activity reduction in response to noxious hindlimb stimulation in the caudate putamen, insula, and secondary somatosensory cortex compared to tonic stimulation.

Arterial blood stealing as a mechanism of negative BOLD response: From the steady-flow with nonlinear phase separation to a windkessel-based model

Blood Oxygen Level Dependent (BOLD) signal indirectly characterizes neuronal activity by measuring hemodynamic and metabolic changes in the nearby microvasculature. A deeper understanding of how localized changes in electrical, metabolic and hemodynamic factors translate into a BOLD signal is crucial for the interpretation of functional brain imaging techniques. While positive BOLD responses (PBR) are widely considered to be linked with neuronal activation, the origins of negative BOLD responses (NBR) have remained largely unknown. As NBRs are sometimes observed in close proximity of regions with PBR, a blood "stealing" effect, i.e., redirection of blood from a passive periphery to the area with high neuronal activity, has been postulated. In this study, we used the Hagen-Poiseuille equation to model hemodynamics in an idealized microvascular network that account for the particulate nature of blood and nonlinearities arising from the red blood cell (RBC) distribution (i.e., the Fåhraeus, Fåhraeus-Lindqvist and the phase separation effects). Using this detailed model, we evaluate determinants driving this "stealing" effect in a microvascular network with geometric parameters within physiological ranges. Model simulations predict that during localized cerebral blood flow (CBF) increases due to neuronal activation-hyperemic response, blood from surrounding vessels is reallocated towards the activated region. This stealing effect depended on the resistance of the microvasculature and the uneven distribution of RBCs at vessel bifurcations. A parsimonious model consisting of two-connected windkessel regions sharing a supplying artery was proposed to simulate the stealing effect with a minimum number of parameters. Comparison with the detailed model showed that the parsimonious model can reproduce the observed response for hematocrit values within the physiological range for different species. Our novel parsimonious model promise to be of use for statistical inference (top-down analysis) from direct blood flow measurements (e.g., arterial spin labeling and laser Doppler/Speckle flowmetry), and when combined with theoretical models for oxygen extraction/diffusion will help account for some types of NBRs.

Identification of Negative BOLD Responses in Epilepsy Using Windkessel Models

Alongside positive blood oxygenation level–dependent (BOLD) responses associated with interictal epileptic discharges, a variety of negative BOLD responses (NBRs) are typically found in epileptic patients. Previous studies suggest that, in general, up to four mechanisms might underlie the genesis of NBRs in the brain: (i) neuronal disruption of network activity, (ii) altered balance of neurometabolic/vascular couplings, (iii) arterial blood stealing, and (iv) enhanced cortical inhibition. Detecting and classifying these mechanisms from BOLD signals are pivotal for the improvement of the specificity of the electroencephalography–functional magnetic resonance imaging (EEG-fMRI) image modality to identify the seizure-onset zones in refractory local epilepsy. This requires models with physiological interpretation that furnish the understanding of how these mechanisms are fingerprinted by their BOLD responses. Here, we used a Windkessel model with viscoelastic compliance/inductance in combination with dynamic models of both neuronal population activity and tissue/blood O₂ to classify the hemodynamic response functions (HRFs) linked to the above mechanisms in the irritative zones of epileptic patients. First, we evaluated the most relevant imprints on the BOLD response caused by variations of key model parameters. Second, we demonstrated that a general linear model is enough to accurately represent the four different types of NBRs. Third, we tested the ability of a machine learning classifier, built from a simulated ensemble of HRFs, to predict the mechanism underlying the BOLD signal from irritative zones. Cross-validation indicates that these four mechanisms can be classified from realistic fMRI BOLD signals. To demonstrate proof of concept, we applied our methodology to EEG-fMRI data from five epileptic patients undergoing neurosurgery, suggesting the presence of some of these mechanisms. We concluded that a proper identification and interpretation of NBR mechanisms in epilepsy can be performed by combining general linear models and biophysically inspired models.

Moderating Effect of Cortical Thickness on BOLD Signal Variability Age-Related Changes

The time course of neuroanatomical structural and functional measures across the lifespan is commonly reported in association with aging. Blood oxygen-level dependent signal variability, estimated using the standard deviation of the signal, or "BOLDSD," is an emerging metric of variability in neural processing, and has been shown to be positively correlated with cognitive flexibility. Generally, BOLDSD is reported to decrease with aging, and is thought to reflect age-related cognitive decline. Additionally, it is well established that normative aging is associated with structural changes in brain regions, and that these predict functional decline in various cognitive domains. Nevertheless, the interaction between alterations in cortical morphology and BOLDSD changes has not been modeled quantitatively. The objective of the current study was to investigate the influence of cortical morphology metrics [i.e., cortical thickness (CT), gray matter (GM) volume, and cortical area (CA)] on age-related BOLDSD changes by treating these cortical morphology metrics as possible physiological confounds using linear mixed models. We studied these metrics in 28 healthy older subjects scanned twice at approximately 2.5 years interval. Results show that BOLDSD is confounded by cortical morphology metrics. Respectively, changes in CT but not GM volume nor CA, show a significant interaction with BOLDSD alterations. Our study highlights that CT changes should be considered when evaluating BOLDSD alternations in the lifespan.

Functional connectivity with cortical depth assessed by resting state fMRI of subregions of S1 in squirrel monkeys

Whereas resting state blood oxygenation-level dependent (BOLD) functional MRI has been widely used to assess functional connectivity between cortical regions, the laminar specificity of such measures is poorly understood. This study aims to determine: (a) whether the resting state functional connectivity (rsFC) between two functionally related cortical regions varies with cortical depth, (b) the relationship between layer-resolved tactile stimulus-evoked activation pattern and interlayer rsFC pattern between two functionally distinct but related somatosensory areas 3b and 1, and (c) the effects of spatial resolution on rsFC measures. We examined the interlayer rsFC between areas 3b and 1 of squirrel monkeys under anesthesia using tactile stimulus-driven and resting state BOLD acquisitions at submillimeter resolution. Consistent with previous observations in the areas 3b and 1, we detected robust stimulus-evoked BOLD activations with foci were confined mainly to the upper layers (centered at 21% of the cortical depth). By carefully placing seeds in upper, middle, and lower layers of areas 3b and 1, we observed strong rsFC between upper and middle layers of these two areas. The layer-resolved activation patterns in areas 3b and 1 agree with their interlayer rsFC patterns, and are consistent with the known anatomical connections between layers. In summary, using BOLD rsFC pattern, we identified an interlayer interareal microcircuit that shows strong intrinsic functional connections between upper and middle layer areas 3b and 1. RsFC can be used as a robust invasive tool to probe interlayer corticocortical microcircuits.

Investigation of the neurovascular coupling in positive and negative BOLD responses in human brain at 7 T

Decreases in stimulus-dependent blood oxygenation level dependent (BOLD) signal and their underlying neurovascular origins have recently gained considerable interest. In this study a multi-echo, BOLD-corrected vascular space occupancy (VASO) functional magnetic resonance imaging (fMRI) technique was used to investigate neurovascular responses during stimuli that elicit positive and negative BOLD responses in human brain at 7 T. Stimulus-induced BOLD, cerebral blood volume (CBV), and cerebral blood flow (CBF) changes were measured and analyzed in 'arterial' and 'venous' blood compartments in macro- and microvasculature. We found that the overall interplay of mean CBV, CBF and BOLD responses is similar for tasks inducing positive and negative BOLD responses. Some aspects of the neurovascular coupling however, such as the temporal response, cortical depth dependence, and the weighting between 'arterial' and 'venous' contributions, are significantly different for the different task conditions. Namely, while for excitatory tasks the BOLD response peaks at the cortical surface, and the CBV change is similar in cortex and pial vasculature, inhibitory tasks are associated with a maximum negative BOLD response in deeper layers, with CBV showing strong constriction of surface arteries and a faster return to baseline. The different interplays of CBV, CBF and BOLD during excitatory and inhibitory responses suggests different underlying hemodynamic mechanisms.

Weighing brain activity with the balance: Angelo Mosso's original manuscripts come to light

Neuroimaging techniques, such as positron emission tomography and functional magnetic resonance imaging are essential tools for the analysis of organized neural systems in working and resting states, both in physiological and pathological conditions. They provide evidence of coupled metabolic and cerebral local blood flow changes that strictly depend upon cellular activity. In 1890, Charles Smart Roy and Charles Scott Sherrington suggested a link between brain circulation and metabolism. In the same year William James, in his introduction of the concept of brain blood flow variations during mental activities, briefly reported the studies of the Italian physiologist Angelo Mosso, a multifaceted researcher interested in the human circulatory system. James focused on Mosso's recordings of brain pulsations in patients with skull breaches, and in the process only briefly referred to another invention of Mosso's, the 'human circulation balance', which could non-invasively measure the redistribution of blood during emotional and intellectual activity. However, the details and precise workings of this instrument and the experiments Mosso performed with it have remained largely unknown. Having found Mosso's original manuscripts in the archives, we remind the scientific community of his experiments with the 'human circulation balance' and of his establishment of the conceptual basis of non-invasive functional neuroimaging techniques. Mosso unearthed and investigated several critical variables that are still relevant in modern neuroimaging such as the 'signal-to-noise ratio', the appropriate choice of the experimental paradigm and the need for the simultaneous recording of differing physiological parameters.