Hemodynamic and metabolic responses to neuronal inhibition

On the interplay between state-dependent reconfigurations of global signal correlation and BOLD fluctuations: An fMRI study

BACKGROUND

The dynamics of global, state-dependent reconfigurations in brain connectivity are yet unclear. We aimed at assessing reconfigurations of the global signal correlation coefficient (GSCORR), a measure of the connectivity between each voxel timeseries and the global signal, from resting-state to a stop-signal task. The secondary aim was to assess the relationship between GSCORR and blood-oxygen-level-dependent (BOLD) activations or deactivation across three different trial-conditions (GO, STOP-correct, and STOP-incorrect).

METHODS

As primary analysis we computed whole-brain, voxel-wise GSCORR during resting-state (GSCORR-rest) and stop-signal task (GSCORR-task) in 107 healthy subjects aged 21-50, deriving GSCORR-shift as GSCORR-task minus GSCORR-rest. GSCORR-tr and trGSCORR-shift were also computed on the task residual time series to quantify the impact of the task-related activity during the trials. To test the secondary aim, brain regions were firstly divided in one cluster showing significant task-related activation and one showing significant deactivation across the three trial conditions. Then, correlations between GSCORR-rest/task/shift and activation/deactivation in the two clusters were computed. As sensitivity analysis, GSCORR-shift was computed on the same sample after performing a global signal regression and GSCORR-rest/task/shift were correlated with the task performance.

RESULTS

Sensory and temporo-parietal regions exhibited a negative GSCORR-shift. Conversely, associative regions (ie. left lingual gyrus, bilateral dorsal posterior cingulate gyrus, cerebellum areas, thalamus, posterolateral parietal cortex) displayed a positive GSCORR-shift (FDR-corrected p < 0.05). GSCORR-shift showed similar patterns to trGSCORR-shift (magnitude increased) and after global signal regression (magnitude decreased). Concerning BOLD changes, Brodmann area 6 and inferior parietal lobule showed activation, while posterior parietal lobule, cuneus, precuneus, middle frontal gyrus showed deactivation (FDR-corrected p < 0.05). No correlations were found between GSCORR-rest/task/shift and beta-coefficients in the activation cluster, although negative correlations were observed between GSCORR-task and GO/STOP-correct deactivation (Pearson rho=-0.299/-0.273; Bonferroni-p < 0.05). Weak associations between GSCORR and task performance were observed (uncorrected p < 0.05).

CONCLUSION

GSCORR state-dependent reconfiguration indicates a reallocation of functional resources to associative areas during stop-signal task. GSCORR, activation and deactivation may represent distinct proxies of brain states with specific neurofunctional relevance.

Tracing the development of learned song preferences in the female zebra finch brain with functional magnetic resonance imaging

In sexually dimorphic zebra finches (Taeniopygia guttata), only males learn to sing their father's song, whereas females learn to recognize the songs of their father or mate but cannot sing themselves. Memory of learned songs is behaviorally expressed in females by preferring familiar songs over unfamiliar ones. Auditory association regions such as the caudomedial mesopallium (CMM; or caudal mesopallium) have been shown to be key nodes in a network that supports preferences for learned songs in adult females. However, much less is known about how song preferences develop during the sensitive period of learning in juvenile female zebra finches. In this study, we used blood-oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) to trace the development of a memory-based preference for the father's song in female zebra finches. Using BOLD fMRI, we found that only in adult female zebra finches with a preference for learned song over novel conspecific song, neural selectivity for the father's song was localized in the thalamus (dorsolateral nucleus of the medial thalamus; part of the anterior forebrain pathway, AFP) and in CMM. These brain regions also showed a selective response in juvenile female zebra finches, although activation was less prominent. These data reveal that neural responses in CMM, and perhaps also in the AFP, are shaped during development to support behavioral preferences for learned songs.

Rat superior colliculus encodes the transition between static and dynamic vision modes

The visual continuity illusion involves a shift in visual perception from static to dynamic vision modes when the stimuli arrive at high temporal frequency, and is critical for recognizing objects moving in the environment. However, how this illusion is encoded across the visual pathway remains poorly understood, with disparate frequency thresholds at retinal, cortical, and behavioural levels suggesting the involvement of other brain areas. Here, we employ a multimodal approach encompassing behaviour, whole-brain functional MRI, and electrophysiological measurements, for investigating the encoding of the continuity illusion in rats. Behavioural experiments report a frequency threshold of 18±2 Hz. Functional MRI reveal that superior colliculus signals transition from positive to negative at the behaviourally-driven threshold, unlike thalamic and cortical areas. Electrophysiological recordings indicate that these transitions are underpinned by neural activation/suppression. Lesions in the primary visual cortex reveal this effect to be intrinsic to the superior colliculus (under a cortical gain effect). Our findings highlight the superior colliculus' crucial involvement in encoding temporal frequency shifts, especially the change from static to dynamic vision modes.

Thermodynamic limitations on brain oxygen metabolism: physiological implications

Recent thermodynamic modelling indicates that maintaining the brain tissue ratio of O2 to CO2 (abbreviated tissue O2 /CO2 ) is critical for preserving the entropy increase available from oxidative metabolism of glucose, with a fall of that available entropy leading to a reduction of the phosphorylation potential and impairment of brain energy metabolism. This provides a novel perspective for understanding physiological responses under different conditions in terms of preserving tissue O2 /CO2 . To enable estimation of tissue O2 /CO2 in the human brain, a detailed mathematical model of O2 and CO2 transport was developed, and applied to reported physiological responses to different challenges, asking: how well is tissue O2 /CO2 preserved? Reported experimental results for increased neural activity, hypercapnia and hypoxia due to high altitude are consistent with preserving tissue O2 /CO2 . The results highlight two physiological mechanisms that control tissue O2 /CO2 : cerebral blood flow, which modulates tissue O2 ; and ventilation rate, which modulates tissue CO2 . The hypoxia modelling focused on humans at high altitude, including acclimatized lowlanders and Tibetan and Andean adapted populations, with a primary finding that decreasing CO2 by increasing ventilation rate is more effective for preserving tissue O2 /CO2 than increasing blood haemoglobin content to maintain O2 delivery to tissue. This work focused on the function served by particular physiological responses, and the underlying mechanisms require further investigation. The modelling provides a new framework and perspective for understanding how blood flow and other physiological factors support energy metabolism in the brain under a wide range of conditions. KEY POINTS: Thermodynamic modelling indicates that preserving the O2 /CO2 ratio in brain tissue is critical for preserving the entropy change available from oxidative metabolism of glucose and the phosphorylation potential underlying energy metabolism. A detailed model of O2 and CO2 transport was developed to allow estimation of the tissue O2 /CO2 ratio in the human brain in different physiological states. Reported experimental results during hypoxia, hypercapnia and increased oxygen metabolic rate in response to increased neural activity are consistent with maintaining brain tissue O2 /CO2 ratio. The hypoxia modelling of high-altitude acclimatization and adaptation in humans demonstrates the critical role of reducing CO2 with increased ventilation for preserving tissue O2 /CO2 . Preservation of tissue O2 /CO2 provides a novel perspective for understanding the function of observed physiological responses under different conditions in terms of preserving brain energy metabolism, although the mechanisms underlying these functions are not well understood.

Functional and Structural Properties of Interhemispheric Interaction between Bilateral Precentral Hand Motor Regions in a Top Wheelchair Racing Paralympian

Long-term motor training can cause functional and structural changes in the human brain. Assessing how the training of specific movements affects specific parts of the neural circuitry is essential to understand better the underlying mechanisms of motor training-induced plasticity in the human brain. We report a single-case neuroimaging study that investigated functional and structural properties in a professional athlete of wheelchair racing. As wheelchair racing requires bilateral synchronization of upper limb movements, we hypothesized that functional and structural properties of interhemispheric interactions in the central motor system might differ between the professional athlete and controls. Functional and diffusion magnetic resonance imaging (fMRI and dMRI) data were obtained from a top Paralympian (P1) in wheelchair racing. With 23 years of wheelchair racing training starting at age eight, she holds an exceptional competitive record. Furthermore, fMRI and dMRI data were collected from three other paraplegic participants (P2-P4) with long-term wheelchair sports training other than wheelchair racing and 37 able-bodied control volunteers. Based on the fMRI data analyses, P1 showed activation in the bilateral precentral hand sections and greater functional connectivity between these sections during a right-hand unimanual task. In contrast, other paraplegic participants and controls showed activation in the contralateral hemisphere and deactivation in the ipsilateral hemisphere. Moreover, dMRI data analysis revealed that P1 exhibited significantly lower mean diffusivity along the transcallosal pathway connecting the bilateral precentral motor regions than control participants, which was not observed in the other paraplegic participants. These results suggest that long-term training with bilaterally synchronized upper-limb movements may promote bilateral recruitment of the precentral hand sections. Such recruitment may affect the structural circuitry involved in the interhemispheric interaction between the bilateral precentral regions. This study provides valuable evidence of the extreme adaptability of the human brain.

Correlation Between Cerebral Tissue Oxygen Saturation and Oxygen Extraction Fraction During Anesthesia: Monitoring Cerebral Metabolic Demand-supply Balance During Vasopressor Administration

BACKGROUND

The speculation that cerebral tissue oxygen saturation (SctO 2 ) measured using tissue near-infrared spectroscopy reflects the balance between cerebral metabolic rate of oxygen and cerebral oxygen delivery has not been validated. Our objective was to correlate SctO 2 with cerebral oxygen extraction fraction (OEF) measured using positron emission tomography; OEF is the ratio between cerebral metabolic rate of oxygen and cerebral oxygen delivery and reflects the balance between these 2 variables.

MATERIALS AND METHODS

This cohort study was based on data collected in a previously published trial assessing phenylephrine versus ephedrine treatment in anesthetized patients undergoing brain tumor surgery. The variables of interest were measured twice over the healthy hemisphere before surgery: the first measurement performed after anesthesia induction and the second measurement performed after induction of a ∼20% increase in blood pressure using either phenylephrine or ephedrine.

RESULTS

Data from 24 patients were analyzed. The overall vasopressor-induced relative changes in SctO 2 (ΔSctO 2 ) and OEF (ΔOEF) were 3.16% [interquartile range, -0.73% to 6.04%] and -12.5% [interquartile range, -24.0% to -6.19%], respectively. ΔSctO 2 negatively correlated with ΔOEF after phenylephrine treatment (Spearman rank correlation coefficient [ rs ]=-0.76; P =0.007), ephedrine treatment ( rs =-0.76; P =0.006), and any treatment ( rs =-0.79; P <0.001). ΔSctO 2 significantly associated with ΔOEF based on multivariable analysis with ΔOEF, relative changes in mean arterial pressure, arterial blood oxygen tension, and the bispectral index as covariates ( P =0.036).

CONCLUSIONS

The negative correlation between changes in SctO 2 and OEF suggests that SctO 2 may reflect the cerebral metabolic demand-supply balance during vasopressor treatment. The generalizability of our findings in other clinical scenarios remains to be determined.

Tracing development of song memory with fMRI in zebra finches after a second tutoring experience

Sensory experiences in early development shape higher cognitive functions such as language acquisition in humans and song learning in birds. Zebra finches (Taeniopygia guttata) sequentially exposed to two different song 'tutors' during the sensitive period in development are able to learn from their second tutor and eventually imitate aspects of his song, but the neural substrate involved in learning a second song is unknown. We used fMRI to examine neural activity associated with learning two songs sequentially. We found that acquisition of a second song changes lateralization of the auditory midbrain. Interestingly, activity in the caudolateral Nidopallium (NCL), a region adjacent to the secondary auditory cortex, was related to the fidelity of second-song imitation. These findings demonstrate that experience with a second tutor can permanently alter neural activity in brain regions involved in auditory perception and song learning.

Ultra-high-field pharmacological functional MRI of dopamine D1 receptor-related interventions in anesthetized rats

The dopamine D1 receptor (D1R) is associated with schizophrenia, Parkinson's disease, and attention deficit hyperactivity disorder. Although the receptor is considered a therapeutic target for these diseases, its neurophysiological function has not been fully elucidated. Pharmacological functional MRI (phfMRI) has been used to evaluate regional brain hemodynamic changes induced by neurovascular coupling resulting from pharmacological interventions, thus phfMRI studies can be used to help understand the neurophysiological function of specific receptors. Herein, the blood oxygenation level-dependent (BOLD) signal changes associated with D1R action in anesthetized rats was investigated by using a preclinical ultra-high-field 11.7-T MRI scanner. PhfMRI was performed before and after administration of the D1-like receptor agonist (SKF82958), antagonist (SCH39166), or physiological saline subcutaneously. Compared to saline, the D1-agonist induced a BOLD signal increase in the striatum, thalamus, prefrontal cortex, and cerebellum. At the same time, the D1-antagonist reduced the BOLD signal in the striatum, thalamus, and cerebellum by evaluating temporal profiles. PhfMRI detected D1R-related BOLD signal changes in the brain regions associated with high expression of D1R. We also measured the early expression of c-fos at the mRNA level to evaluate the effects of SKF82958 and isoflurane anesthesia on neuronal activity. Regardless of the presence of isoflurane anesthesia, c-fos expression level was increased in the region where positive BOLD responses were observed with administration of SKF82958. These findings demonstrated that phfMRI could be used to identify the effects of direct D1 blockade on physiological brain functions and also for neurophysiological assessment of dopamine receptor functions in living animals.

Dynamic Effective Connectivity using Physiologically informed Dynamic Causal Model with Recurrent Units: A functional Magnetic Resonance Imaging simulation study

Functional MRI (fMRI) is an indirect reflection of neuronal activity. Using generative biophysical model of fMRI data such as Dynamic Causal Model (DCM), the underlying neuronal activities of different brain areas and their causal interactions (i.e., effective connectivity) can be calculated. Most DCM studies typically consider the effective connectivity to be static for a cognitive task within an experimental run. However, changes in experimental conditions during complex tasks such as movie-watching might result in temporal variations in the connectivity strengths. In this fMRI simulation study, we leverage state-of-the-art Physiologically informed DCM (P-DCM) along with a recurrent window approach and discretization of the equations to infer the underlying neuronal dynamics and concurrently the dynamic (time-varying) effective connectivities between various brain regions for task-based fMRI. Results from simulation studies on 3- and 10-region models showed that functional magnetic resonance imaging (fMRI) blood oxygenation level-dependent (BOLD) responses and effective connectivity time-courses can be accurately predicted and distinguished from faulty graphical connectivity models representing cognitive hypotheses. In summary, we propose and validate a novel approach to determine dynamic effective connectivity between brain areas during complex cognitive tasks by combining P-DCM with recurrent units.

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.

Cerebral blood flow predicts multiple demand network activity and fluid intelligence across the adult lifespan

The preservation of cognitive function in old age is a public health priority. Cerebral hypoperfusion is a hallmark of dementia but its impact on maintaining cognitive ability across the lifespan is less clear. We investigated the relationship between baseline cerebral blood flow (CBF) and blood oxygenation level-dependent (BOLD) response during a fluid reasoning task in a population-based adult lifespan cohort. As age differences in CBF could lead to non-neuronal contributions to the BOLD signal, we introduced commonality analysis to neuroimaging to dissociate performance-related CBF effects from the physiological confounding effects of CBF on the BOLD response. Accounting for CBF, we confirmed that performance- and age-related differences in BOLD responses in the multiple-demand network were implicated in fluid reasoning. Age differences in CBF explained not only performance-related BOLD responses but also performance-independent BOLD responses. Our results suggest that CBF is important for maintaining cognitive function, while its non-neuronal contributions to BOLD signals reflect an age-related confound. Maintaining perfusion into old age may serve to support brain function and preserve cognitive performance.

External drivers of BOLD signal's non-stationarity

A fundamental challenge in neuroscience is to uncover the principles governing how the brain interacts with the external environment. However, assumptions about external stimuli fundamentally constrain current computational models. We show in silico that unknown external stimulation can produce error in the estimated linear time-invariant dynamical system. To address these limitations, we propose an approach to retrieve the external (unknown) input parameters and demonstrate that the estimated system parameters during external input quiescence uncover spatiotemporal profiles of external inputs over external stimulation periods more accurately. Finally, we unveil the expected (and unexpected) sensory and task-related extra-cortical input profiles using functional magnetic resonance imaging data acquired from 96 subjects (Human Connectome Project) during the resting-state and task scans. This dynamical systems model of the brain offers information on the structure and dimensionality of the BOLD signal's external drivers and shines a light on the likely external sources contributing to the BOLD signal's non-stationarity. Our findings show the role of exogenous inputs in the BOLD dynamics and highlight the importance of accounting for external inputs to unravel the brain's time-varying functional dynamics.

Temporal stability of the hemodynamic response function across the majority of human cerebral cortex

The hemodynamic response function (HRF) measured with functional magnetic resonance imaging is generated by vascular and metabolic responses evoked by brief (<4 s) stimuli. It is known that the human HRF varies across cortex, between subjects, with stimulus paradigms, and even between different measurements in the same cortical location. However, our results demonstrate that strong HRFs are remarkably repeatable across sessions separated by time intervals up to 3 months. In this study, a multisensory stimulus was used to activate and measure the HRF across the majority of cortex (>70%, with lesser reliability observed in some areas of prefrontal cortex). HRFs were measured with high spatial resolution (2‐mm voxels) in central gray matter to minimize variations caused by partial‐volume effects. HRF amplitudes and temporal dynamics were highly repeatable across four sessions in 20 subjects. Positive and negative HRFs were consistently observed across sessions and subjects. Negative HRFs were generally weaker and, thus, more variable than positive HRFs. Statistical measurements showed that across‐session variability is highly correlated to the variability across events within a session; these measurements also indicated a normal distribution of variability across cortex. The overall repeatability of the HRFs over long time scales generally supports the long‐term use of event‐related functional magnetic resonance imaging protocols.

Evidence suggesting common mechanisms underlie contralateral and ipsilateral negative BOLD responses in human visual cortex

The task-evoked positive BOLD response (PBR) to a unilateral visual hemi-field stimulation is often accompanied by robust and sustained contralateral as well as ipsilateral negative BOLD responses (NBRs) in the visual cortex. The signal characteristics and the neural and/or vascular mechanisms that underlie these two types of NBRs are not completely understood. In this paper, we investigated the properties of these two types of NBRs. We first demonstrated the linearity of both NBRs with respect to stimulus duration. Next, we showed that the hemodynamic response functions (HRFs) of the two NBRs were similar to each other, but significantly different from that of the PBR. Moreover, the subject-wise expressions of the two NBRs were tightly coupled to the degree that the correlation between the two NBRs was significantly higher than the correlation between each NBR and the PBR. However, the activation patterns of the two NBRs did not show a high level of interhemispheric spatial similarity, and the functional connectivity between them was not different than the interhemispheric functional connectivity between the NBRs and PBR. Finally, while attention did modulate both NBRs, the attention-related changes in their HRFs were similar. Our findings suggest that the two NBRs might be generated through common neural and/or vascular mechanisms involving distal/deep brain regions that project to the two hemispheres.

A Systematic Review of the Impact of Physical Exercise-Induced Increased Resting Cerebral Blood Flow on Cognitive Functions

Brain perfusion declines with aging. Physical exercise represents a low-cost accessible form of intervention to increase cerebral blood flow; however, it remains unclear if exercise-induced amelioration of brain perfusion has any impact on cognition. We aimed to provide a state-of-the art review on this subject. A comprehensive search of the PubMed (MEDLINE) database was performed. On the basis of the inclusion and exclusion criteria, 14 studies were included in the analysis. Eleven of the studies conducted well-controlled exercise programs that lasted 12–19 weeks for 10–40 participants and two studies were conducted in much larger groups of subjects for more than 5 years, but the exercise loads were indirectly measured, and three of them were focused on acute exercise. Literature review does not show a direct link between exercise-induced augmentation of brain perfusion and better cognitive functioning. However, in none of the reviewed studies was such an association the primary study endpoint. Carefully designed clinical studies with focus on cognitive and perfusion variables are needed to provide a response to the question whether exercise-induced cerebral perfusion augmentation is of clinical importance.

Characterization of cortical hemodynamic changes following sensory, visual, and speech activation by intraoperative optical imaging utilizing phase-based evaluation methods

Alterations within cerebral hemodynamics are the intrinsic signal source for a wide variety of neuroimaging techniques. Stimulation of specific functions leads due to neurovascular coupling, to changes in regional cerebral blood flow, oxygenation and volume. In this study, we investigated the temporal characteristics of cortical hemodynamic responses following electrical, tactile, visual, and speech activation for different stimulation paradigms using Intraoperative Optical Imaging (IOI). Image datasets from a total of 22 patients that underwent surgical resection of brain tumors were evaluated. The measured reflectance changes at different light wavelength bands, representing alterations in regional cortical blood volume (CBV), and deoxyhemoglobin (HbR) concentration, were assessed by using Fourier-based evaluation methods. We found a decrease of CBV connected to an increase of HbR within the contralateral primary sensory cortex (SI) in patients that were prolonged (30 s/15 s) electrically stimulated. Additionally, we found differences in amplitude as well as localization of activated areas for different stimulation patterns. Contrary to electrical stimulation, prolonged tactile as well as prolonged visual stimulation are provoking increases in CBV within the corresponding activated areas (SI, visual cortex). The processing of the acquired data from awake patients performing speech tasks reveals areas with increased, as well as areas with decreased CBV. The results lead us to the conclusion, that the CBV decreases in connection with HbR increases in SI are associated to processing of nociceptive stimuli and that stimulation type, as well as paradigm have a nonnegligible impact on the temporal characteristics of the following hemodynamic response.

Increased migraine-free intervals with multifocal repetitive transcranial magnetic stimulation

INTRODUCTION

Episodic migraine is a debilitating condition associated with vast impairments of health, daily living, and life quality. Several prophylactic treatments exist, having a moderate ratio of action related to side effects and therapy costs. Repetitive transcranial magnetic stimulation (rTMS) is an evidence based therapy in several neuropsychiatric conditions, showing robust efficacy in alleviating specific symptoms. However, its efficacy in migraine disorders is unequivocal and might be tightly linked to the applied rTMS protocol. We hypothesized that multifocal rTMS paradigm could improve clinical outcomes in patients with episodic migraine by reducing the number of migraine days, frequency and intensity of migraine attacks, and improve the quality of life.

METHODS

We conducted an experimental, double-blind, randomized controlled study by applying a multifocal rTMS paradigm. Patients with episodic migraine with or without aura were enrolled in two centers from August 2018, to December 2019, and randomized to receive either real (n = 37) or sham (sham coil stimulation, n = 28) multifocal rTMS for six sessions over two weeks. Patients, physicians, and raters were blinded to the applied protocol. The experimental multifocal rTMS protocol included two components; first, swipe stimulation of 13 trains of 140 pulses/train, 67 Hz, 60% of RMT, and 2s intertrain interval and second, spot burst stimulation of 33 trains of 15 pulses/train, 67 Hz, 85% of RMT, and 8s intertrain interval. Reduction >50% from the baseline in migraine days (as primary outcome) and frequency and intensity of migraine attacks (as key secondary outcomes) over a 12-week period were assessed. To balance the baseline variables between the treatment arms, we applied the propensity score matching through the logistic regression.

RESULTS

Among 65 randomized patients, sixty (age 39.7 ± 11.6; 52 females; real rTMS n = 33 and sham rTMS n = 27) completed the trial and five patients dropped out. Over 12 weeks, the responder's rate in the number of migraine days was significantly higher in the real rTMS compared to the sham group (42% vs. 26%, p < 0.05). The mean migraine days per month decreased from 7.6 to 4.3 days in the real rTMS group and from 6.2 to 4.3 days in the sham rTMS group, resulting in a difference with real vs. sham rTMS of -3.2 days (p < 0.05). Similarly, over the 12-week period, the responder's rate in the reduction of migraine attacks frequency was higher in the real rTMS compared to the sham group (42% vs 33%, p < 0.05). No serious adverse events were observed.

CONCLUSION

Our pilot study shows compelling evidence in a double placebo-controlled trial that multifocal rTMS is an effective and well-tolerated preventive treatment in patients with episodic migraine.

Acupuncture With deqi Modulates the Hemodynamic Response and Functional Connectivity of the Prefrontal-Motor Cortical Network

As a world intangible cultural heritage, acupuncture is considered an essential modality of complementary and alternative therapy to Western medicine. Despite acupuncture’s long history and public acceptance, how the cortical network is modulated by acupuncture remains largely unclear. Moreover, as the basic acupuncture unit for regulating the central nervous system, how the cortical network is modulated during acupuncture at the Hegu acupoint is mostly unclear. Here, multi-channel functional near-infrared spectroscopy (fNIRS) data were recorded from twenty healthy subjects for acupuncture manipulation, pre- and post-manipulation tactile controls, and pre- and post-acupuncture rest controls. Results showed that: (1) acupuncture manipulation caused significantly increased acupuncture behavioral deqi performance compared with tactile controls. (2) The bilateral prefrontal cortex (PFC) and motor cortex were significantly inhibited during acupuncture manipulation than controls, which was evidenced by the decreased power of oxygenated hemoglobin (HbO) concentration. (3) The bilateral PFC’s hemodynamic responses showed a positive correlation trend with acupuncture behavioral performance. (4) The network connections with bilateral PFC as nodes showed significantly increased functional connectivity during acupuncture manipulation compared with controls. (5) Meanwhile, the network’s efficiency was improved by acupuncture manipulation, evidenced by the increased global efficiency and decreased shortest path length. Taken together, these results reveal that a cooperative PFC-Motor functional network could be modulated by acupuncture manipulation at the Hegu acupoint. This study provides neuroimaging evidence that explains acupuncture’s neuromodulation effects on the cortical network.

Existence of Interhemispheric Inhibition between Foot Sections of Human Primary Motor Cortices: Evidence from Negative Blood Oxygenation-Level Dependent Signal

Interhemispheric inhibition (IHI) between the left and right primary motor cortices (M1) plays an important role when people perform an isolated unilateral limb movement. Moreover, negative blood oxygenation-level dependent signal (deactivation) obtained from the M1 ipsilateral to the limb could be a surrogate IHI marker. Studies have reported deactivation in the hand section of the ipsilateral M1 during simple unilateral hand movement. However, deactivation in the foot section during unilateral foot movement has not been reported. Therefore, IHI between the foot sections of the bilateral M1s has been considered very weak or absent. Thirty-seven healthy adults performed active control of the right foot and also passively received vibration to the tendon of the tibialis anterior muscle of the right foot, which activates the foot section of the contralateral M1, with brain activity being examined through functional magnetic resonance imaging. The vibration and active tasks significantly and non-significantly, respectively, deactivated the foot section of the ipsilateral M1, with a corresponding 86% and 60% of the participants showing decreased activity. Thus, there could be IHI between the foot sections of the bilateral M1s. Further, our findings demonstrate between-task differences and similarities in cross-somatotopic deactivation.

FMRI and intra-cranial electrocorticography recordings in the same human subjects reveals negative BOLD signal coupled with silenced neuronal activity

Positive blood oxygenation level-dependent (BOLD) responses (PBR), as measured by functional Magnetic Resonance Imaging (fMRI), are the most utilized measurements to non-invasively map activity in the brain. Recent studies have consistently shown that BOLD responses are not exclusively positive. Negative BOLD responses (NBR) have been reported in response to specific sensory stimulations and tasks. However, the exact relationship between NBR and the underlying metabolic and neuronal demand is still under debate. In this study, we investigated the neurophysiological basis of negative BOLD using fMRI and intra-cranial electrophysiology (electrocorticography, ECoG) measurements from the same human participants. We show that, for those electrodes that responded to visual stimulation, PBR are correlated with high-frequency band (HFB) responses. Crucially, NBR were associated with an absence of HFB power responses and an unpredicted decrease in the alpha power responses.

Hemodynamic Responses Link Individual Differences in Informational Masking to the Vicinity of Superior Temporal Gyrus

Suppressing unwanted background sound is crucial for aural communication. A particularly disruptive type of background sound, informational masking (IM), often interferes in social settings. However, IM mechanisms are incompletely understood. At present, IM is identified operationally: when a target should be audible, based on suprathreshold target/masker energy ratios, yet cannot be heard because target-like background sound interferes. We here confirm that speech identification thresholds differ dramatically between low- vs. high-IM background sound. However, speech detection thresholds are comparable across the two conditions. Moreover, functional near infrared spectroscopy recordings show that task-evoked blood oxygenation changes near the superior temporal gyrus (STG) covary with behavioral speech detection performance for high-IM but not low-IM background sound, suggesting that the STG is part of an IM-dependent network. Moreover, listeners who are more vulnerable to IM show increased hemodynamic recruitment near STG, an effect that cannot be explained based on differences in task difficulty across low- vs. high-IM. In contrast, task-evoked responses near another auditory region of cortex, the caudal inferior frontal sulcus (cIFS), do not predict behavioral sensitivity, suggesting that the cIFS belongs to an IM-independent network. Results are consistent with the idea that cortical gating shapes individual vulnerability to IM.

The neural-vascular basis of age-related processing speed decline

Most studies examining neurocognitive aging are based on the blood-oxygen level-dependent signal obtained during functional magnetic resonance imaging (fMRI). The physiological basis of this signal is neural-vascular coupling, the process by which neurons signal cerebrovasculature to dilate in response to an increase in active neural metabolism due to stimulation. These fMRI studies of aging rely on the hemodynamic equivalence assumption that this process is not disrupted by physiologic deterioration associated with aging. Studies of neural-vascular coupling challenge this assumption and show that neural-vascular coupling is closely related to cognition. In this review, we put forward a theory of processing speed decline in aging and how it is related to age-related neural-vascular coupling changes based on the results of studies elucidating the relationships between cognition, cerebrovascular dynamics, and aging.

Quantitative relations between BOLD responses, cortical energetics and impulse firing across cortical depth

The blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) signal arises as a consequence of changes in cerebral blood flow (CBF) and cerebral metabolic rate of oxygen ( CMR O 2 ) that in turn are modulated by changes in neural activity. Recent advances in imaging have achieved sub-millimetre resolution and allowed investigation of the BOLD response as a function of cortical depth. Here, we adapt our previous theory relating the BOLD signal to neural activity to produce a quantitative model that incorporates venous blood draining between cortical layers. The adjustable inputs to the model are the neural activity and a parameter governing this blood draining. A three-layer version for transient neural inputs and a multi-layer version for constant or tonic neural inputs are able to account for a variety of experimental results, including negative BOLD signals.

Opposed hemodynamic responses following increased excitation and parvalbumin-based inhibition

Understanding cellular contributions to hemodynamic activity is essential for interpreting blood-based brain mapping signals. Optogenetic studies examining cell-specific influences on local hemodynamics have reported that excitatory activity results in cerebral perfusion and blood volume increase, while inhibitory activity contributes to both vasodilation and vasoconstriction. How specific subpopulations of interneurons regulate the brain's blood supply is less examined. Parvalbumin interneurons are the largest subpopulation of GABAergic neurons in the brain, critical for brain development, plasticity, and long-distance excitatory neurotransmission. Despite their essential role in brain function, the contribution of parvalbumin neurons to neurovascular coupling has been relatively unexamined. Using optical intrinsic signal imaging and laser speckle contrast imaging, we photostimulated awake and anesthetized transgenic mice expressing channelrhodopsin under a parvalbumin promoter. Increased parvalbumin activity reduced local oxygenation, cerebral blood volume, and cerebral blood flow. These "negative" hemodynamic responses were consistent within and across mice and reproducible across a broad range of photostimulus parameters. However, the sign and magnitude of the hemodynamic response resulting from increased parvalbumin activity depended on the type and level of anesthesia used. Opposed hemodynamic responses following increased excitation or parvalbumin-based inhibition suggest unique contributions from different cell populations to neurovascular coupling.

Retinotopic variations of the negative blood-oxygen-level dependent hemodynamic response function in human primary visual cortex

Functional magnetic resonance imaging (fMRI) measures blood-oxygen-level-dependent (BOLD) contrast that is generally assumed to be linearly related to excitatory neural activity. The positive hemodynamic response function (pHRF) is the positive BOLD response (PBR) evoked by a brief neural stimulation; the pHRF is often used as the impulse response for linear analysis of neural excitation. Many fMRI studies have observed a negative BOLD response (NBR) that is often associated with neural suppression. However, the temporal dynamics of the NBR evoked by a brief stimulus, the negative HRF (nHRF), remains unclear. Here, a unilateral visual stimulus was presented in a slow event-related design to elicit both pHRFs in the stimulus representation (SR), and nHRFs elsewhere. The observed nHRFs were not inverted versions of the pHRF previously reported. They were characterized by a stronger initial negative response followed by a significantly later positive peak. In contralateral primary visual cortex (V1), these differences varied with eccentricity from the SR. Similar nHRFs were observed in ipsilateral V1 with less eccentricity variation. Experiments with the blocked version of the same stimulus confirmed that brain regions presenting the unexpected nHRF dynamics correspond to those presenting a strong NBR. These data demonstrated that shift-invariant temporal linearity did not hold for the NBR while confirming that the PBR maintained rough linearity. Modeling indicated that the observed nHRFs can be created by suppression of both blood flow and oxygen metabolism. Critically, the nHRF can be misinterpreted as a pHRF due to their similarity, which could confound linear analysis for event-related fMRI experiments.NEW & NOTEWORTHY We investigate dynamics of the negative hemodynamic response function (nHRF), the negative blood-oxygen-level-dependent (BOLD) response (NBR) evoked by a brief stimulus, in human early visual cortex. Here, we show that the nHRFs are not inverted versions of the corresponding pHRFs. The nHRF has complex dynamics that varied significantly with eccentricity. The results also show shift-invariant temporal linearity does not hold for the NBR.

Altered task-induced cerebral blood flow and oxygen metabolism underlies motor impairment in multiple sclerosis

The neural mechanisms underlying motor impairment in multiple sclerosis (MS) remain unknown. Motor cortex dysfunction is implicated in blood-oxygen-level-dependent (BOLD) functional magnetic resonance imaging (fMRI) studies, but the role of neural-vascular coupling underlying BOLD changes remains unknown. We sought to independently measure the physiologic factors (i.e., cerebral blood flow (ΔCBF), cerebral metabolic rate of oxygen (ΔCMRO2), and flow-metabolism coupling (ΔCBF/ΔCMRO2), utilizing dual-echo calibrated fMRI (cfMRI) during a bilateral finger-tapping task. We utilized cfMRI to measure physiologic responses in 17 healthy volunteers and 32 MS patients (MSP) with and without motor impairment during a thumb-button-press task in thumb-related (task-central) and surrounding primary motor cortex (task-surround) regions of interest (ROIs). We observed significant ΔCBF and ΔCMRO2 increases in all MSP compared to healthy volunteers in the task-central ROI and increased flow-metabolism coupling (ΔCBF/ΔCMRO2) in the MSP without motor impairment. In the task-surround ROI, we observed decreases in ΔCBF and ΔCMRO2 in MSP with motor impairment. Additionally, ΔCBF and ΔCMRO2 responses in the task-surround ROI were associated with motor function and white matter damage in MSP. These results suggest an important role for task-surround recruitment in the primary motor cortex to maintain motor dexterity and its dependence on intact white matter microstructure and neural-vascular coupling.

The thermodynamics of thinking: connections between neural activity, energy metabolism and blood flow

Several current functional neuroimaging methods are sensitive to cerebral metabolism and cerebral blood flow (CBF) rather than the underlying neural activity itself. Empirically, the connections between metabolism, flow and neural activity are complex and somewhat counterintuitive: CBF and glycolysis increase more than seems to be needed to provide oxygen and pyruvate for oxidative metabolism, and the oxygen extraction fraction is relatively low in the brain and decreases when oxygen metabolism increases. This work lays a foundation for the idea that this unexpected pattern of physiological changes is consistent with basic thermodynamic considerations related to metabolism. In the context of this thermodynamic framework, the apparent mismatches in metabolic rates and CBF are related to preserving the entropy change of oxidative metabolism, specifically the O₂/CO₂ ratio in the mitochondria. However, the mechanism supporting this CBF response is likely not owing to feedback from a hypothetical O₂ sensor in tissue, but rather is consistent with feed-forward control by signals from both excitatory and inhibitory neural activity. Quantitative predictions of the thermodynamic framework, based on models of O₂ and CO₂ transport and possible neural drivers of CBF control, are in good agreement with a wide range of experimental data, including responses to neural activation, hypercapnia, hypoxia and high-altitude acclimatization. This article is part of the theme issue 'Key relationships between non-invasive functional neuroimaging and the underlying neuronal activity'.

Voxel-based morphometry and task functional magnetic resonance imaging in essential tremor: evidence for a disrupted brain network

The pathophysiology of essential tremor (ET) is controversial and might be further elucidated by advanced neuroimaging. Focusing on homogenous ET patients diagnosed according to the 2018 consensus criteria, this study aimed to: (1) investigate whether task functional MRI (fMRI) can identify networks of activated and deactivated brain areas, (2) characterize morphometric and functional modulations, relative to healthy controls (HC). Ten ET patients and ten HC underwent fMRI while performing two motor tasks with their upper limb: (1) maintaining a posture (both groups); (2) simulating tremor (HC only). Activations/deactivations were obtained from General Linear Model and compared across groups/tasks. Voxel-based morphometry and linear regressions between clinical and fMRI data were also performed. Few cerebellar clusters of gray matter loss were found in ET. Conversely, widespread fMRI alterations were shown. Tremor in ET (task 1) was associated with extensive deactivations mainly involving the cerebellum, sensory-motor cortex, and basal ganglia compared to both tasks in HC, and was negatively correlated with clinical tremor scales. Homogeneous ET patients demonstrated deactivation patterns during tasks triggering tremor, encompassing a network of cortical and subcortical regions. Our results point towards a marked cerebellar involvement in ET pathophysiology and the presence of an impaired cerebello-thalamo-cortical tremor network.

Evidence for asymmetric inhibitory activity during motor planning phases of sensorimotor synchronization

Sensorimotor synchronization (SMS) is frequently dependent on coordination of excitatory and inhibitory activity across hemispheres, as well as the cognitive control over environmental distractors. However, the timing (motor planning versus execution) and cortical regions involved in these processes remain actively debated. Functional magnetic resonance imaging data were therefore analyzed from 34 strongly right-handed healthy adults performing a cued (to initiate motor planning) SMS task with either their right or left hand (motor execution phase) based on spatially congruent or incongruent visual stimuli. Behavioral effects of incongruent stimuli were limited to the first stimulus. Functionally, greater activation was observed in left sensorimotor cortex (SMC) and right cerebellar Lobule V for congruent versus incongruent stimuli. A negative blood-oxygen level dependent response, a putative marker of neural inhibition, was present in bilateral SMC, right supplemental motor area (SMA) and bilateral cerebellar Lobule V during the motor planning, but not execution phase. The magnitude of the inhibitory response was greater in right cortical regions and cerebellar Lobule V. Homologue connectivity was associated with inhibitory activity in the right SMA, suggesting that individual differences in intrinsic connectivity may mediate transcallosal inhibition. In summary, results suggest increased inhibition (i.e., greater negative BOLD response) within the right relative to left hemisphere, which was released once motor programs were executed. Both task and intrinsic functional connectivity results highlight a critical role of the left SMA in interhemispheric inhibition and motor planning.

Key Aspects of Neurovascular Control Mediated by Specific Populations of Inhibitory Cortical Interneurons

Inhibitory interneurons can evoke vasodilation and vasoconstriction, making them potential cellular drivers of neurovascular coupling. However, the specific regulatory roles played by particular interneuron subpopulations remain unclear. Our purpose was therefore to adopt a cell-specific optogenetic approach to investigate how somatostatin (SST) and neuronal nitric oxide synthase (nNOS)-expressing interneurons might influence the neurovascular relationship. In mice, specific activation of SST- or nNOS-interneurons was sufficient to evoke hemodynamic changes. In the case of nNOS-interneurons, robust hemodynamic changes occurred with minimal changes in neural activity, suggesting that the ability of blood oxygen level dependent functional magnetic resonance imaging (BOLD fMRI) to reliably reflect changes in neuronal activity may be dependent on type of neuron recruited. Conversely, activation of SST-interneurons produced robust changes in evoked neural activity with shallow cortical excitation and pronounced deep layer cortical inhibition. Prolonged activation of SST-interneurons often resulted in an increase in blood volume in the centrally activated area with an accompanying decrease in blood volume in the surrounding brain regions, analogous to the negative BOLD signal. These results demonstrate the role of specific populations of cortical interneurons in the active control of neurovascular function.

The neurovascular basis of processing speed differences in humans: A model-systems approach using multiple sclerosis

Behavioral studies investigating fundamental cognitive abilities provide evidence that processing speed accounts for large proportions of performance variability between individuals. Processing speed decline is a hallmark feature of the cognitive disruption observed in healthy aging and in demyelinating diseases such as multiple sclerosis (MS), neuromyelitis optica, and Wilson's disease. Despite the wealth of evidence suggesting a central role for processing speed in cognitive decline, the neural mechanisms of this fundamental ability remain unknown. Intact neurovascular coupling, acute localized blood flow increases following neural activity, is essential for optimal neural function. We hypothesized that efficient coupling forms the neural basis of processing speed. Because MS features neural-glial-vascular system disruption, we used it as a model to test this hypothesis. To assess the integrity of the coupling system, we measured blood-oxygen-level-dependent (BOLD) signal in healthy controls (HCs) and MS patients using a 3T MRI scanner while they viewed radial checkerboards that flickered periodically at 8 ​Hz. To assess processing speed and cognitive function, we administered a battery of neuropsychological tests. While MS patients exhibited reduced ΔBOLD with reductions in processing speed, no such relationships were observed in HCs. To further investigate the mechanisms that underlie ΔBOLD-processing speed relationships, we assessed the physiologic components that constitute ΔBOLD signal (i.e., cerebral blood flow, ΔCBF; cerebral metabolic rate of oxygen, ΔCMRO₂; neurovascular coupling ratio) in speed-preserved and -impaired MS patients. While ΔCBF and ΔCMRO₂ showed no group-differences, the neurovascular coupling ratio was significantly reduced in speed-impaired MS patients compared to speed-preserved MS patients. Together, these results suggest that neurovascular uncoupling might underlie cognitive slowing in MS and might be the central pathogenic mechanism governing processing speed decline.

Metabolite concentration changes associated with positive and negative BOLD responses in the human visual cortex: A functional MRS study at 7 Tesla

Negative blood oxygenation-level dependent (BOLD) signal observed during task execution in functional magnetic resonance imaging (fMRI) can be caused by different mechanisms, such as a blood-stealing effect or neuronal deactivation. Electrophysiological recordings showed that neuronal deactivation underlies the negative BOLD observed in the occipital lobe during visual stimulation. In this study, the metabolic demand of such a response was studied by measuring local metabolite concentration changes during a visual checkerboard stimulation using functional magnetic resonance spectroscopy (fMRS) at 7 Tesla. The results showed increases of glutamate and lactate concentrations during the positive BOLD response, consistent with previous fMRS studies. In contrast, during the negative BOLD response, decreasing concentrations of glutamate, lactate and gamma-aminobutyric acid (GABA) were found, suggesting a reduction of glycolytic and oxidative metabolic demand below the baseline. Additionally, the respective changes of the BOLD signal, glutamate and lactate concentrations of both groups suggest that a local increase of inhibitory activity might occur during the negative BOLD response.

Inhibitory Neuron Activity Contributions to Hemodynamic Responses and Metabolic Load Examined Using an Inhibitory Optogenetic Mouse Model

Hemodynamic signals are routinely used to noninvasively assess brain function in humans and animals. This work examined the contribution of inhibitory neuron activity on hemodynamic responses captured by changes in blood flow, volume and oxygenation in the cortex of lightly anesthetized mice. Because cortical activity is not commonly initiated by inhibitory neurons, experiments were conducted to examine the neuronal activity properties elicited by photo-stimulation. We observed comparable increases in neuronal activity evoked by forelimb and photo-stimulation; however, significantly larger increases in blood flow and volume were produced by photo-stimulation of inhibitory neurons compared with forelimb stimulation. Following blockade of glutamate and GABA-A receptors to reduce postsynaptic activity contributions, neuronal activity was reliably modulated and hemodynamic changes persisted, though slightly reduced. More importantly, photo-stimulation-evoked changes in blood flow and volume were suppressed by 75–80% with the administration of a nitric oxide synthase inhibitor, suggesting that inhibitory neurons regulate blood flow mostly via nitric oxide. Lastly, forelimb and photo-stimulation of excitatory neurons produced local decreases in blood oxygenation, while large increases were generated by photo-stimulation of inhibitory neurons. Estimates of oxygen metabolism suggest that inhibitory neuron activity has a small impact on tissue metabolic load, indicating a mismatch between the metabolic demand and blood flow regulation properties of inhibitory and excitatory neurons.

Dynamic metabolic changes in human visual cortex in regions with positive and negative blood oxygenation level-dependent response

Dynamic metabolic changes were investigated by functional magnetic resonance spectroscopy (fMRS) during sustained stimulation of human primary visual cortex. Two established paradigms, consisting of either a full-field or a small-circle flickering checkerboard, were employed to generate wide-spread areas of positive or negative blood oxygenation level-dependent (BOLD) responses, respectively. Compared to baseline, the glutamate concentration increased by 5.3% (p = 0.007) during activation and decreased by -3.8% (p = 0.017) during deactivation. These changes were positively correlated with the amplitude of the BOLD response (R = 0.60, p = 0.002) and probably reflect changes of tricarboxylic acid cycle activity. During deactivation, the glucose concentration decreased by -7.9% (p = 0.025) presumably suggesting increased consumption or reduced glucose supply. Other findings included an increased concentration of glutathione (4.2%, p = 0.023) during deactivation and a negative correlation of glutathione and BOLD signal changes (R = -0.49, p = 0.012) as well as positive correlations of aspartate (R = 0.44, p = 0.035) and N-acetylaspartylglutamate (R = 0.42, p = 0.035) baseline concentrations with the BOLD response. It remains to be shown in future work if the observed effects on glutamate and glucose levels deviate from the assumption of a direct link between glucose utilization and regulation of blood flow or support previous suggestions that the hemodynamic response is mainly driven by feedforward release of vasoactive messengers.

Negative BOLD responses during hand and foot movements: An fMRI study

The present study employed functional magnetic resonance imaging (fMRI) to examine the characteristics of negative blood oxygen level-dependent (Negative BOLD) signals during motor execution. Subjects repeated extension and flexion of one of the following: the right hand, left hand, right ankle, or left ankle. Negative BOLD responses during hand movements were observed in the ipsilateral hemisphere of the hand primary sensorimotor area (SMI), medial frontal gyrus (MeFG), middle frontal gyrus (MFG), and superior frontal gyrus (SFG). Negative BOLD responses during foot movements were also noted in the bilateral hand SMI, MeFG, MFG, SFG, inferior frontal gyrus, middle temporal gyrus, parahippocampal gyrus, anterior cingulate cortex, cingulate gyrus (CG), fusiform gyrus, and precuneus. A conjunction analysis showed that portions of the MeFG and CG involving similar regions to those of the default mode network were commonly deactivated during voluntary movements of the right/left hand or foot. The present results suggest that three mechanisms are involved in the Negative BOLD responses observed during voluntary movements: (1) transcallosal inhibition from the contralateral to ipsilateral hemisphere in the SMI, (2) the deactivated neural network with several brain regions, and (3) the default mode network in the MeFG and CG.

Preserved canonicality of the BOLD hemodynamic response reflects healthy cognition: Insights into the healthy brain through the window of Multiple Sclerosis

The hemodynamic response function (HRF), a model of brain blood-flow changes in response to neural activity, reflects communication between neurons and the vasculature that supplies these neurons in part by means of glial cell intermediaries (e.g., astrocytes). Intact neural-vascular communication might play a central role in optimal cognitive performance. This hypothesis can be tested by comparing healthy individuals to those with known white-matter damage and impaired performance, as seen in Multiple Sclerosis (MS). Glial cell intermediaries facilitate the ability of neurons to adequately convey metabolic needs to cerebral vasculature for sufficient oxygen and nutrient perfusion. In this study, we isolated measurements of the HRF that could quantify the extent to which white-matter affects neural-vascular coupling and cognitive performance. HRFs were modeled from multiple brain regions during multiple cognitive tasks using piecewise cubic spline functions, an approach that minimized assumptions regarding HRF shape that may not be valid for diseased populations, and were characterized using two shape metrics (peak amplitude and time-to-peak). Peak amplitude was reduced, and time-to-peak was longer, in MS patients relative to healthy controls. Faster time-to-peak was predicted by faster reaction time, suggesting an important role for vasodilatory speed in the physiology underlying processing speed. These results support the hypothesis that intact neural-glial-vascular communication underlies optimal neural and cognitive functioning.

Feasibility of functional MRI at ultralow magnetic field via changes in cerebral blood volume

We investigate the feasibility of performing functional MRI (fMRI) at ultralow field (ULF) with a Superconducting QUantum Interference Device (SQUID), as used for detecting magnetoencephalography (MEG) signals from the human head. While there is negligible magnetic susceptibility variation to produce blood oxygenation level-dependent (BOLD) contrast at ULF, changes in cerebral blood volume (CBV) may be a sensitive mechanism for fMRI given the five-fold spread in spin-lattice relaxation time (T₁) values across the constituents of the human brain. We undertook simulations of functional signal strength for a simplified brain model involving activation of a primary cortical region in a manner consistent with a blocked task experiment. Our simulations involve measured values of T₁ at ULF and experimental parameters for the performance of an upgraded ULFMRI scanner. Under ideal experimental conditions we predict a functional signal-to-noise ratio of between 3.1 and 7.1 for an imaging time of 30 min, or between 1.5 and 3.5 for a blocked task experiment lasting 7.5 min. Our simulations suggest it may be feasible to perform fMRI using a ULFMRI system designed to perform MRI and MEG in situ.

An fMRI study of finger tapping in children and adults

Functional brain imaging studies have characterized the neural bases of voluntary movement for finger tapping in adults, but equivalent information for children is lacking. When contrasted to adults, one would expect children to have relatively greater activation, reflecting compensation for an underdeveloped motor system combined with less experience in the execution of voluntary movement. To test this hypothesis, we acquired functional magnetic resonance imaging (fMRI) data on 17 healthy right-handed children (7.48 ± 0.66 years) and 15 adults (24.9 ± 2.9 years) while they performed an irregularly paced finger-tapping task in response to a visual cue (left- and right-hand examined separately). Whole-brain within-group analyses revealed that finger tapping in either age group and for either hand activated contralateral SM1, SMA, ipsilateral anterior cerebellum, and occipital cortices. We used an ANOVA factorial design to test for main effects of Age Group (children vs adults), Hand (left vs. right), and their interactions. For main effects of Age Group, children showed relatively greater activity in left SM1 (extending into bilateral SMA), and, surprisingly, adults exhibited relatively greater activity in right pre-SMA/SMA (extending into left pre-SMA/SMA), right lateral globus pallidus, left putamen, and right anterior cerebellum. The interaction of Age Group × Hand revealed that while both groups activated right SM1 during left finger tapping and exhibited signal decreases (i.e., below fixation baseline) during right finger tapping, both these responses were attenuated in children relative to adults. These data provide an important foundation by which to study children with motor disorders.

Characterization of the hemodynamic response function across the majority of human cerebral cortex

A brief (<4 s) period of neural activation evokes a stereotypical sequence of vascular and metabolic events to create the hemodynamic response function (HRF) measured using functional magnetic resonance imaging (fMRI). Linear analysis of fMRI data requires that the HRF be treated as an impulse response, so the character and temporal stability of the HRF are critical issues. Here, a simple audiovisual stimulus combined with a fast-paced task was used to evoke a strong HRF across a majority, ∼77%, of cortex during a single scanning session. High spatiotemporal resolution (2-mm voxels, 1.25-s acquisition time) was used to focus HRF measurements specifically on the gray matter for whole brain. The majority of activated cortex responds with positive HRFs, while ∼27% responds with negative (inverted) HRFs. Spatial patterns of the HRF response amplitudes were found to be similar across subjects. Timing of the initial positive lobe of the HRF was relatively stable across the cortical surface with a mean of 6.1 ± 0.6 s across subjects, yet small but significant timing variations were also evident in specific regions of cortex. The results provide guidance for linear analysis of fMRI data. More importantly, this method provides a means to quantify neurovascular function across most of the brain, with potential clinical utility for the diagnosis of brain pathologies such as traumatic brain injury.

Quantitative relations between BOLD responses, cortical energetics, and impulse firing

The blood oxygen level-dependent (BOLD) functional magnetic resonance imaging signal arises as a consequence of changes in blood flow and oxygen usage that in turn are modulated by changes in neural activity. Much attention has been given to both theoretical and experimental aspects of the energetics but not to the neural activity. Here we identify the best energetic theory for the steady-state BOLD signal on the basis of correct predictions of experimental observations. This theory is then used, together with the recently determined relationship between energetics and neural activity, to predict how the BOLD signal changes with activity. Unlike existing treatments, this new theory incorporates a nonzero baseline activity in a completely consistent way and is thus able to account for both sustained positive and negative BOLD signals. We also show that the increase in BOLD signal for a given increase in activity is significantly smaller the larger the baseline activity, as is experimentally observed. Furthermore, the decline of the positive BOLD signal arising from deeper cortical laminae in response to an increase in neural firing is shown to arise as a consequence of the larger baseline activity in deeper laminae. Finally, we provide quantitative relations integrating BOLD responses, energetics, and impulse firing, which among other predictions give the same results as existing theories when the baseline activity is zero. NEW & NOTEWORTHY We use a recently established relation between energetics and neural activity to give a quantitative account of BOLD dependence on neural activity. The incorporation of a nonzero baseline neural activity accounts for positive and negative BOLD signals, shows that changes in neural activity give BOLD changes that are smaller the larger the baseline, and provides a basis for the observed inverse relation between BOLD responses and the depth of cortical laminae giving rise to them.

Mapping and characterization of positive and negative BOLD responses to visual stimulation in multiple brain regions at 7T

External stimuli and tasks often elicit negative BOLD responses in various brain regions, and growing experimental evidence supports that these phenomena are functionally meaningful. In this work, the high sensitivity available at 7T was explored to map and characterize both positive (PBRs) and negative BOLD responses (NBRs) to visual checkerboard stimulation, occurring in various brain regions within and beyond the visual cortex. Recently-proposed accelerated fMRI techniques were employed for data acquisition, and procedures for exclusion of large draining vein contributions, together with ICA-assisted denoising, were included in the analysis to improve response estimation. Besides the visual cortex, significant PBRs were found in the lateral geniculate nucleus and superior colliculus, as well as the pre-central sulcus; in these regions, response durations increased monotonically with stimulus duration, in tight covariation with the visual PBR duration. Significant NBRs were found in the visual cortex, auditory cortex, default-mode network (DMN) and superior parietal lobule; NBR durations also tended to increase with stimulus duration, but were significantly less sustained than the visual PBR, especially for the DMN and superior parietal lobule. Responses in visual and auditory cortex were further studied for checkerboard contrast dependence, and their amplitudes were found to increase monotonically with contrast, linearly correlated with the visual PBR amplitude. Overall, these findings suggest the presence of dynamic neuronal interactions across multiple brain regions, sensitive to stimulus intensity and duration, and demonstrate the richness of information obtainable when jointly mapping positive and negative BOLD responses at a whole-brain scale, with ultra-high field fMRI.

Bilateral versus ipsilesional cortico-subcortical activity patterns in stroke show hemispheric dependence

Background Understanding of interhemispheric interactions in stroke patients during motor control is an important clinical neuroscience quest that may provide important clues for neurorehabilitation. In stroke patients bilateral overactivation in both hemispheres has been interpreted as a poor prognostic indicator of functional recovery. In contrast, ipsilesional patterns have been linked with better motor outcomes. Aim We investigated the pathophysiology of hemispheric interactions during limb movement without and with contralateral restraint, to mimic the effects of constraint-induced movement therapy. We used neuroimaging to probe brain activity with such a movement-dependent interhemispheric modulation paradigm. Methods We used a functional magnetic resonance imaging block design during which the plegic/paretic upper limb was recruited/mobilized to perform unilateral arm elevation, as a function of presence versus absence of contralateral limb restriction (n = 20, with balanced left/right lesion sites). Results Analysis of 10 right hemispheric stroke participants yielded bilateral sensorimotor cortex activation in all movement phases in contrast with the unilateral dominance seen in the 10 left hemispheric stroke participants. Superimposition of contralateral restriction led to a prominent shift from activation to deactivation response patterns, in particular in cortical and basal ganglia motor areas in right hemispheric stroke. Left hemispheric stroke was, in general, characterized by reduced activation patterns, even in the absence of restriction, which induced additional cortical silencing. Conclusion The observed hemispheric-dependent activation/deactivation shifts is novel and these pathophysiological observations suggest short-term neuroplasticity that may be useful for hemisphere-tailored neurorehabilitation.

High-Resolution CBV-fMRI Allows Mapping of Laminar Activity and Connectivity of Cortical Input and Output in Human M1

Layer-dependent fMRI allows measurements of information flow in cortical circuits, as afferent and efferent connections terminate in different cortical layers. However, it is unknown to what level human fMRI is specific and sensitive enough to reveal directional functional activity across layers. To answer this question, we developed acquisition and analysis methods for blood-oxygen-level-dependent (BOLD) and cerebral-blood-volume (CBV)-based laminar fMRI and used these to discriminate four different tasks in the human motor cortex (M1). In agreement with anatomical data from animal studies, we found evidence for somatosensory and premotor input in superficial layers of M1 and for cortico-spinal motor output in deep layers. Laminar resting-state fMRI showed directional functional connectivity of M1 with somatosensory and premotor areas. Our findings demonstrate that CBV-fMRI can be used to investigate cortical activity in humans with unprecedented detail, allowing investigations of information flow between brain regions and outperforming conventional BOLD results that are often buried under vascular biases.

Cortical contributions to sensory gating in the ipsilateral somatosensory cortex during voluntary activity

KEY POINTS

It has long been known that the somatosensory cortex gates sensory inputs from the contralateral side of the body. Here, we examined the contribution of the ipsilateral somatosensory cortex (iS1) to sensory gating during index finger voluntary activity. The amplitude of the P25/N33, but not other somatosensory evoked potential (SSEP) components, was reduced during voluntary activity compared with rest. Interhemispheric inhibition between S1s and intracortical inhibition in the S1 modulated the amplitude of the P25/N33. Note that changes in interhemispheric inhibition between S1s correlated with changes in cortical circuits in the ipsilateral motor cortex. Our findings suggest that cortical circuits, probably from somatosensory and motor cortex, contribute to sensory gating in the iS1 during voluntary activity in humans.

ABSTRACT

An important principle in the organization of the somatosensory cortex is that it processes afferent information from the contralateral side of the body. The role of the ipsilateral somatosensory cortex (iS1) in sensory gating in humans remains largely unknown. Using electroencephalographic (EEG) recordings over the iS1 and electrical stimulation of the ulnar nerve at the wrist, we examined somatosensory evoked potentials (SSEPs; P14/N20, N20/P25 and P25/N33 components) and paired-pulse SSEPs between S1s (interhemispheric inhibition) and within (intracortical inhibition) the iS1 at rest and during tonic index finger voluntary activity. We found that the amplitude of the P25/N33, but not other SSEP components, was reduced during voluntary activity compared with rest. Interhemispheric inhibition increased the amplitude of the P25/N33 and intracortical inhibition reduced the amplitude of the P25/N33, suggesting a cortical origin for this effect. The P25/N33 receives inputs from the motor cortex, so we also examined the contribution of distinct sets of cortical interneurons by testing the effect of ulnar nerve stimulation on motor-evoked potentials (MEPs) elicited by transcranial magnetic stimulation over the ipsilateral motor cortex with the coil in the posterior-anterior (PA) and anterior-posterior (AP) orientation. Afferent input attenuated PA, but not AP, MEPs during voluntary activity compared with rest. Notably, changes in interhemispheric inhibition correlated with changes in PA MEPs. Our novel findings suggest that interhemispheric projections between S1s and intracortical circuits, probably from somatosensory and motor cortex, contribute to sensory gating in the iS1 during voluntary activity in humans.

Compression at Myofascial Trigger Point on Chronic Neck Pain Provides Pain Relief through the Prefrontal Cortex and Autonomic Nervous System: A Pilot Study

Compression at myofascial trigger points (MTrPs), known as “ischemic compression,” has been reported to provide immediate relief of musculoskeletal pain and reduce the sympathetic activity that exacerbates chronic pain. We conducted a pilot study to investigate the possible involvement of the prefrontal cortex in pain relief obtained by MTrP compression in the present study, and analyzed the relationships among prefrontal hemodynamic activity, activity of the autonomic nervous system, and subjective pain in patients with chronic neck pain, with and without MTrP compression. Twenty-one female subjects with chronic neck pain were randomly assigned to two groups: MTrP compression (n = 11) or Non-MTrP compression (n = 10). Compression for 30 s was conducted 4 times. During the experiment, prefrontal hemodynamic activity [changes in Oxy-hemoglobin (Hb), Deoxy-Hb, and Total-Hb concentrations] and autonomic activity based on heart rate variability (HRV) were monitored by using near infrared spectroscopy (NIRS) and electrocardiography (ECG), respectively. The results indicated that MTrP compression significantly reduced subjective pain compared with Non-MTrP compression. The spectral frequency-domain analyses of HRV indicated that a low frequency (LF) component of HRV was decreased, and a high frequency (HF) component of HRV was increased during MTrP compression, while LF/HF ratio was decreased during MTrP compression. In addition, prefrontal hemodynamic activity was significantly decreased during MTrP compression compared with Non-MTrP compression. Furthermore, changes in autonomic activity were significantly correlated with changes in subjective pain and prefrontal hemodynamic activity. Along with previous studies indicating a role for sympathetic activity in the exacerbation of chronic pain, the present results suggest that MTrP compression in the neck region alters the activity of the autonomic nervous system via the prefrontal cortex to reduce subjective pain.