Global signal modulation of single-trial fMRI response variability: Effect on positive vs negative BOLD response relationship

Can structure predict function at individual level in the human connectome?

Several studies predicting Functional Connectivity (FC) from Structural Connectivity (SC) at individual level have been published in recent years, each promising increased performance and utility. We investigated three of these studies, analyzing whether the results truly represent a meaningful individual-level mapping from SC to FC. Using data from the Human Connectome Project shared accross the three studies, we constructed a predictor by averaging FC of training data and analyzed its performance in the same way. In each case, we found that group average FC is an equivalent or better predictor of individual FC than the predictive models in terms of raw prediction performance. Furthermore, we showed that additional analyses performed by the authors of the three studies, in which they attempt to show that their predicted FC has value beyond raw prediction performance, could also be reproduced using the group average FC predictor. This makes it unclear whether any of the three methods represent a meaningful individual-level predictive model. We conclude that either the methods are not appropriate for the data, that the sample size is too small, or that the data does not contain sufficient information to learn a mapping from SC to FC. We advise future individual-level studies to explicitly report results in comparison to the performance of the group average, and carefully demonstrate that their predictions contain meaningful individual-level information. Finally, we believe that investigating alternatives for the construction of SC and FC may improve the chances of developing a meaningful individual-level mapping from SC to FC.

Altruistic acting caused by a touching hand: neural underpinnings of the Midas touch effect

Giving and receiving touch are some of the most important social stimuli we exchange in daily life. By touching someone, we can communicate various types of information. Previous studies have also demonstrated that interpersonal touch may affect our altruistic behavior. A classic study showed that customers give bigger tips when they are lightly touched by a waitress, which has been called the Midas touch effect. Numerous studies reported similar effects of touch on different kinds of helping or prosocial behaviors. Here, we aim to examine the neural underpinnings of this effect by employing a functional magnetic resonance imaging approach. While lying in the scanner, participants played different rounds of the dictator game, a measure of prosocial behavior. Before each round, participants were touched (or not touched in the control condition) by an experimenter. We found that touching the hand increased the likeliness to behave prosocial (but not the general liking of control stimuli), thereby confirming the Midas touch effect. The effect was predicted by activity in the primary somatosensory cortex, indicating that the somatosensory cortex here plays a causal role in prosocial behavior. We conclude that the tactile modality in social life may be much more important than previously thought.

Pre-processing of Sub-millimeter GE-BOLD fMRI Data for Laminar Applications

Over the past 30 years, brain function has primarily been evaluated non-invasively using functional magnetic resonance imaging (fMRI) with gradient-echo (GE) sequences to measure blood-oxygen-level-dependent (BOLD) signals. Despite the multiple advantages of GE sequences, e.g., higher signal-to-noise ratio, faster acquisitions, etc., their relatively inferior spatial localization compromises the routine use of GE-BOLD in laminar applications. Here, in an attempt to rescue the benefits of GE sequences, we evaluated the effect of existing pre-processing methods on the spatial localization of signals obtained with EPIK, a GE sequence that affords voxel volumes of 0.25 mm3 with near whole-brain coverage. The methods assessed here apply to both task and resting-state fMRI data assuming the availability of reconstructed magnitude and phase images.

Brain-wide ongoing activity is responsible for significant cross-trial BOLD variability

A notorious issue of task-based functional magnetic resonance imaging (fMRI) is its large cross-trial variability. To quantitatively characterize this variability, the blood oxygenation level-dependent (BOLD) signal can be modeled as a linear summation of a stimulation-relevant and an ongoing (i.e. stimulation-irrelevant) component. However, systematic investigation on the spatiotemporal features of the ongoing BOLD component and how these features affect the BOLD response is still lacking. Here we measured fMRI responses to light onsets and light offsets in awake rats. The neuronal response was simultaneously recorded with calcium-based fiber photometry. We established that between-region BOLD signals were highly correlated brain-wide at zero time lag, including regions that did not respond to visual stimulation, suggesting that the ongoing activity co-fluctuates across the brain. Removing this ongoing activity reduced cross-trial variability of the BOLD response by ~30% and increased its coherence with the Ca2+ signal. Additionally, the negative ongoing BOLD activity sometimes dominated over the stimulation-driven response and contributed to the post-stimulation BOLD undershoot. These results suggest that brain-wide ongoing activity is responsible for significant cross-trial BOLD variability, and this component can be reliably quantified and removed to improve the reliability of fMRI response. Importantly, this method can be generalized to virtually all fMRI experiments without changing stimulation paradigms.

Separating vascular and neuronal effects of age on fMRI BOLD signals

Accurate identification of brain function is necessary to understand the neurobiology of cognitive ageing, and thereby promote well-being across the lifespan. A common tool used to investigate neurocognitive ageing is functional magnetic resonance imaging (fMRI). However, although fMRI data are often interpreted in terms of neuronal activity, the blood oxygenation level-dependent (BOLD) signal measured by fMRI includes contributions of both vascular and neuronal factors, which change differentially with age. While some studies investigate vascular ageing factors, the results of these studies are not well known within the field of neurocognitive ageing and therefore vascular confounds in neurocognitive fMRI studies are common. Despite over 10 000 BOLD-fMRI papers on ageing, fewer than 20 have applied techniques to correct for vascular effects. However, neurovascular ageing is not only a confound in fMRI, but an important feature in its own right, to be assessed alongside measures of neuronal ageing. We review current approaches to dissociate neuronal and vascular components of BOLD-fMRI of regional activity and functional connectivity. We highlight emerging evidence that vascular mechanisms in the brain do not simply control blood flow to support the metabolic needs of neurons, but form complex neurovascular interactions that influence neuronal function in health and disease. This article is part of the theme issue 'Key relationships between non-invasive functional neuroimaging and the underlying neuronal activity'.

Cortical Processing of Multimodal Sensory Learning in Human Neonates

Following birth, infants must immediately process and rapidly adapt to the array of unknown sensory experiences associated with their new ex-utero environment. However, although it is known that unimodal stimuli induce activity in the corresponding primary sensory cortices of the newborn brain, it is unclear how multimodal stimuli are processed and integrated across modalities. The latter is essential for learning and understanding environmental contingencies through encoding relationships between sensory experiences; and ultimately likely subserves development of life-long skills such as speech and language. Here, for the first time, we map the intracerebral processing which underlies auditory-sensorimotor classical conditioning in a group of 13 neonates (median gestational age at birth: 38 weeks + 4 days, range: 32 weeks + 2 days to 41 weeks + 6 days; median postmenstrual age at scan: 40 weeks + 5 days, range: 38 weeks + 3 days to 42 weeks + 1 days) with blood-oxygen-level-dependent (BOLD) functional magnetic resonance imaging (MRI) and magnetic resonance (MR) compatible robotics. We demonstrate that classical conditioning can induce crossmodal changes within putative unimodal sensory cortex even in the absence of its archetypal substrate. Our results also suggest that multimodal learning is associated with network wide activity within the conditioned neural system. These findings suggest that in early life, external multimodal sensory stimulation and integration shapes activity in the developing cortex and may influence its associated functional network architecture.

Functional brain changes in the elderly for the perception of hand movements: A greater impairment occurs in proprioception than touch

Unlike age-related brain changes linked to motor activity, neural alterations related to self-motion perception remain unknown. Using fMRI data, we investigated age-related changes in the central processing of somatosensory information by inducing illusions of right-hand rotations with specific proprioceptive and tactile stimulation. Functional connectivity during resting-state (rs-FC) was also compared between younger and older participants. Results showed common sensorimotor activations in younger and older adults during proprioceptive and tactile illusions, but less deactivation in various right frontal regions and the precuneus were found in the elderly. Older participants exhibited a less-lateralized pattern of activity across the primary sensorimotor cortices (SM1) in the proprioceptive condition only. This alteration of the interhemispheric balance correlated with declining individual performance in illusion velocity perception from a proprioceptive, but not a tactile, origin. By combining task-related data, rs-FC and behavioral performance, this study provided consistent results showing that hand movement perception was altered in the elderly, with a more pronounced deterioration of the proprioceptive system, likely due to the breakdown of inhibitory processes with aging. Nevertheless, older people could benefit from an increase in internetwork connectivity to overcome this kinesthetic decline.

Hemodynamic Correlates of Electrophysiological Activity in the Default Mode Network

Hemodynamic fluctuations in the default mode network (DMN), observed through functional magnetic resonance imaging (fMRI), have been linked to electrophysiological oscillations detected by electroencephalography (EEG). It has been reported that, among the electrophysiological oscillations, those in the alpha frequency range (8–13 Hz) are the most dominant during resting state. We hypothesized that DMN spatial configuration closely depends on the specific neuronal oscillations considered, and that alpha oscillations would mainly correlate with increased blood oxygen-level dependent (BOLD) signal in the DMN. To test this hypothesis, we used high-density EEG (hdEEG) data simultaneously collected with fMRI scanning in 20 healthy volunteers at rest. We first detected the DMN from source reconstructed hdEEG data for multiple frequency bands, and we then mapped the correlation between temporal profile of hdEEG-derived DMN activity and fMRI–BOLD signals on a voxel-by-voxel basis. In line with our hypothesis, we found that the correlation map associated with alpha oscillations, more than with any other frequency bands, displayed a larger overlap with DMN regions. Overall, our study provided further evidence for a primary role of alpha oscillations in supporting DMN functioning. We suggest that simultaneous EEG–fMRI may represent a powerful tool to investigate the neurophysiological basis of human brain networks.

Simultaneous EEG-fMRI: Evaluating the Effect of the EEG Cap-Cabling Configuration on the Gradient Artifact

Electroencephalography (EEG) data recorded during simultaneous EEG-fMRI experiments are contaminated by large gradient artifacts (GA). The amplitude of the GA depends on the area of the wire loops formed by the EEG leads, as well as on the rate of switching of the magnetic field gradients, which are essential for MR imaging. Average artifact subtraction (AAS), the most commonly used method for GA correction, relies on the EEG amplifier having a large enough dynamic range to characterize the artifact voltages. Low-pass filtering (250 Hz cut-off) is generally used to attenuate the high-frequency voltage fluctuations of the GA, but even with this precaution channel saturation can occur, particularly during acquisition of high spatial resolution MRI data. Previous work has shown that the ribbon cable, used to connect the EEG cap and amplifier, makes a significant contribution to the GA, since the cable geometry produces large effective wire-loop areas. However, by appropriately connecting the wires of the ribbon cable to the EEG cap it should be possible to minimize the overall range and root mean square (RMS) amplitude of the GA by producing partial cancelation of the cap and cable contributions. Here by modifying the connections of the EEG cap to a 1 m ribbon cable we were able to reduce the range of the GA for a high-resolution coronal echo planar Imaging (EPI) acquisition by a factor of ∼ 1.6 and by a factor of ∼ 1.15 for a standard axial EPI acquisition. These changes could potentially be translated into a reduction in the required dynamic range, an increase in the EEG bandwidth or an increase in the achievable image resolution without saturation, all of which could be beneficially exploited in EEG-fMRI studies. The re-wiring could also prevent the system from saturating when small subject movements occur using the standard recording bandwidth.

The relationship between negative BOLD responses and ERS and ERD of alpha/beta oscillations in visual and motor cortex

Previous work has investigated the electrophysiological origins of the intra-modal (within the stimulated sensory cortex) negative BOLD fMRI response (NBR, decrease from baseline) but little attention has been paid to the origin of cross-modal NBRs, those in a different sensory cortex. In the current study we use simultaneous EEG-fMRI recordings to assess the neural correlates of both intra- and cross-modal responses to left-hemifield visual stimuli and right-hand motor tasks, and evaluate the balance of activation and deactivation between the visual and motor systems. Within- and between-subject covariations of EEG and fMRI responses to both tasks are assessed to determine how patterns of event-related desynchronization/synchronisation (ERD/ERS) of alpha/beta frequency oscillations relate to the NBR in the two sensory cortices. We show that both visual and motor tasks induce intra-modal NBR and cross-modal NBR (e.g. visual stimuli evoked NBRs in both visual and motor cortices). In the EEG data, bilateral intra-modal alpha/beta ERD were consistently observed to both tasks, whilst the cross-modal EEG response varied across subjects between alpha/beta ERD and ERS. Both the mean cross-modal EEG and fMRI response amplitudes showed a small increase in magnitude with increasing task intensity. In response to the visual stimuli, subjects displaying cross-modal ERS of motor beta power displayed a significantly larger magnitude of cross-modal NBR in motor cortex. However, in contrast to the motor stimuli, larger cross-modal ERD of visual alpha power was associated with larger cross-modal visual NBR. Single-trial correlation analysis provided further evidence of relationship between EEG signals and the NBR, motor cortex beta responses to motor tasks were significantly negatively correlated with cross-modal visual cortex NBR amplitude, and positively correlated with intra-modal motor cortex PBR. This study provides a new body of evidence that the coupling between BOLD and low-frequency (alpha/beta) sensory cortex EEG responses extends to cross-modal NBR.

Intersubject Synchronisation Analysis of Brain Activity Associated with the Instant Effects of Acupuncture: An Fmri Study

Objective

To use a promising analytical method, namely intersubject synchronisation (ISS), to evaluate the brain activity associated with the instant effects of acupuncture and compare the findings with traditional general linear model (GLM) methods.

Methods

30 healthy volunteers were recruited for this study. Block-designed manual acupuncture stimuli were delivered at SP6, and de qi sensations were measured after acupuncture stimulation. All subjects underwent functional MRI (fMRI) scanning during the acupuncture stimuli. The fMRI data were separately analysed by ISS and traditional GLM methods.

Results

All subjects experienced de qi sensations. ISS analysis showed that the regions activated during acupuncture stimulation at SP6 were mainly divided into five clusters based on the time courses. The time courses of clusters 1 and 2 were in line with the acupuncture stimulation pattern, and the active regions were mainly involved in the sensorimotor system and salience network. Clusters 3, 4 and 5 displayed an almost contrary time course relative to the stimulation pattern. The brain regions activated included the default mode network, descending pain modulation pathway and visual cortices. GLM analysis indicated that the brain responses associated with the instant effects of acupuncture were largely implicated in sensory and motor processing and sensory integration.

Conclusion

The ISS analysis considered the sustained effect of acupuncture and uncovered additional information not shown by GLM analysis. We suggest that ISS may be a suitable approach to investigate the brain responses associated with the instant effects of acupuncture.

Neuronal dynamics enable the functional differentiation of resting state networks in the human brain

Intrinsic brain activity is organized in spatial-temporal patterns, called resting-state networks (RSNs), exhibiting specific structural-functional architecture. These networks presumably reflect complex neurophysiological processes and have a central role in distinct perceptual and cognitive functions. In this work, we propose an innovative approach for characterizing RSNs according to their underlying neural oscillations. We investigated specific electrophysiological properties, including spectral features, fractal dimension, and entropy, associated with eight core RSNs derived from high-density electroencephalography (EEG) source-reconstructed signals. Specifically, we found higher synchronization of the gamma-band activity and higher fractal dimension values in perceptual (PNs) compared with higher cognitive (HCNs) networks. The inspection of this underlying rapid activity becomes of utmost importance for assessing possible alterations related to specific brain disorders. The disruption of the coordinated activity of RSNs may result in altered behavioral and perceptual states. Thus, this approach could potentially be used for the early detection and treatment of neurological disorders.

Reprint of 'Noise contributions to the fMRI signal: An Overview'

The ability to discriminate signal from noise plays a key role in the analysis and interpretation of functional magnetic resonance imaging (fMRI) measures of brain activity. Over the past two decades, a number of major sources of noise have been identified, including system-related instabilities, subject motion, and physiological fluctuations. This article reviews the characteristics of the various noise sources as well as the mechanisms through which they affect the fMRI signal. Approaches for distinguishing signal from noise and the associated challenges are also reviewed. These challenges reflect the fact that some noise sources, such as respiratory activity, are generated by the same underlying brain networks that give rise to functional signals that are of interest.

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.

Both activated and less-activated regions identified by functional MRI reconfigure to support task executions

INTRODUCTION

Functional magnetic resonance imaging (fMRI) has become very important for noninvasively characterizing BOLD signal fluctuations, which reflect the changes in neuronal firings in the brain. Unlike the activation detection strategy utilized with fMRI, which only emphasizes the synchronicity between the functional nodes (activated regions) and the task design, brain connectivity and network theory are able to decipher the interactive structure across the entire brain. However, little is known about whether and how the activated/less-activated interactions are associated with the functional changes that occur when the brain changes from the resting state to a task state. What are the key networks that play important roles in the brain state changes?

METHODS

We used the fMRI data from the Human Connectome Project S500 release to examine the changes of network efficiency, interaction strength, and fractional modularity contributions of both the local and global networks, when the subjects change from the resting state to seven different task states.

RESULTS

We found that, from the resting state to each of seven task states, both the activated and less-activated regions had significantly changed to be in line with, and comparably contributed to, a global network reconfiguration. We also found that three networks, the default mode network, frontoparietal network, and salience network, dominated the flexible reconfiguration of the brain.

CONCLUSIONS

This study shows quantitatively that contributions from both activated and less-activated regions enable the global functional network to respond when the brain switches from the resting state to a task state and suggests the necessity of considering large-scale networks (rather than only activated regions) when investigating brain functions in imaging cognitive neuroscience.

Functional Alignment with Anatomical Networks is Associated with Cognitive Flexibility

Cognitive flexibility describes the human ability to switch between modes of mental function to achieve goals. Mental switching is accompanied by transient changes in brain activity, which must occur atop an anatomical architecture that bridges disparate cortical and subcortical regions by underlying white matter tracts. However, an integrated perspective regarding how white matter networks might constrain brain dynamics during cognitive processes requiring flexibility has remained elusive. To address this challenge, we applied emerging tools from graph signal processing to examine whether BOLD signals measured at each point in time correspond to complex underlying anatomical networks in 28 individuals performing a perceptual task that probed cognitive flexibility. We found that the alignment between functional signals and the architecture of the underlying white matter network was associated with greater cognitive flexibility across subjects. By computing a concise measure using multi-modal neuroimaging data, we uncovered an integrated structure-function correlate of human behavior.

The global signal in fMRI: Nuisance or Information?

The global signal is widely used as a regressor or normalization factor for removing the effects of global variations in the analysis of functional magnetic resonance imaging (fMRI) studies. However, there is considerable controversy over its use because of the potential bias that can be introduced when it is applied to the analysis of both task-related and resting-state fMRI studies. In this paper we take a closer look at the global signal, examining in detail the various sources that can contribute to the signal. For the most part, the global signal has been treated as a nuisance term, but there is growing evidence that it may also contain valuable information. We also examine the various ways that the global signal has been used in the analysis of fMRI data, including global signal regression, global signal subtraction, and global signal normalization. Furthermore, we describe new ways for understanding the effects of global signal regression and its relation to the other approaches.

Noise contributions to the fMRI signal: An overview

The ability to discriminate signal from noise plays a key role in the analysis and interpretation of functional magnetic resonance imaging (fMRI) measures of brain activity. Over the past two decades, a number of major sources of noise have been identified, including system-related instabilities, subject motion, and physiological fluctuations. This article reviews the characteristics of the various noise sources as well as the mechanisms through which they affect the fMRI signal. Approaches for distinguishing signal from noise and the associated challenges are also reviewed. These challenges reflect the fact that some noise sources, such as respiratory activity, are generated by the same underlying brain networks that give rise to functional signals that are of interest.