Removal of FMRI environment artifacts from EEG data using optimal basis sets

A multi-center, randomized, double-blind, sham-stimulation controlled study of transcranial magnetic stimulation with precision navigation for the treatment of multiple system atrophy

BACKGROUND

Multiple system atrophy (MSA) is recognized as an atypical Parkinsonian syndrome, distinguished by a more rapid progression than that observed in Parkinson's disease. Unfortunately, the prognosis for MSA remains poor, with a notable absence of globally recognized effective treatments. Although preliminary studies suggest that transcranial magnetic stimulation (TMS) could potentially alleviate clinical symptoms in MSA patients, there is a significant gap in the literature regarding the optimal stimulation parameters. Furthermore, the field lacks consensus due to the paucity of robust, large-scale, multicenter trials.

METHODS

This investigation is a multi-center, randomized, double-blind, sham-controlled trial. We aim to enroll 96 individuals diagnosed with MSA, categorized into Parkinsonian type (MSA-P) and cerebellar type (MSA-C) according to their predominant clinical features. Participants will be randomly allocated in a 1:1 ratio to either the TMS or sham stimulation group. Utilizing advanced navigation techniques, we will ensure precise targeting for the intervention, applying theta burst stimulation (TBS). To assess the efficacy of TBS on both motor and non-motor functions, a comprehensive evaluation will be conducted using internationally recognized clinical scales and gait analysis. To objectively assess changes in brain connectivity, functional magnetic resonance imaging (fMRI) and electroencephalography (EEG) will be employed as sensitive indicators before and after the intervention.

DISCUSSION

The primary aim of this study is to ascertain whether TBS can alleviate both motor and non-motor symptoms in patients with MSA. Additionally, a critical component of our research involves elucidating the underlying mechanisms through which TBS exerts its potential therapeutic effects.

ETHICS AND DISSEMINATION

All study protocols have been reviewed and approved by the First Affiliated Medical Ethics Committee of the Air Force Military Medical University (KY20232118-F-1).

TRIAL REGISTRATION

Chinese Clinical Trial Registry ChiCTR2300072658. Registered on 20 June 2023.

State-specific Regulation of Electrical Stimulation in the Intralaminar Thalamus of Macaque Monkeys: Network and Transcriptional Insights into Arousal

Long-range thalamocortical communication is central to anesthesia-induced loss of consciousness and its reversal. However, isolating the specific neural networks connecting thalamic nuclei with various cortical regions for state-specific anesthesia regulation is challenging, with the biological underpinnings still largely unknown. Here, simultaneous electroencephalogram-fuctional magnetic resonance imaging (EEG-fMRI) and deep brain stimulation are applied to the intralaminar thalamus in macaques under finely-tuned propofol anesthesia. This approach led to the identification of an intralaminar-driven network responsible for rapid arousal during slow-wave oscillations. A network-based RNA-sequencing analysis is conducted of region-, layer-, and cell-specific gene expression data from independent transcriptomic atlases and identifies 2489 genes preferentially expressed within this arousal network, notably enriched in potassium channels and excitatory, parvalbumin-expressing neurons, and oligodendrocytes. Comparison with human RNA-sequencing data highlights conserved molecular and cellular architectures that enable the matching of homologous genes, protein interactions, and cell types across primates, providing novel insight into network-focused transcriptional signatures of arousal.

Denoising of the gradient artifact present in simultaneous studies of muscle blood oxygen level dependent (BOLD) signal and electromyography (EMG)

The MR-induced gradient artifact affects EMG recordings during simultaneous muscle BOLD/EMG acquisitions. However, no dedicated hardware can remove the gradient artifact easily, and alternative methods are expensive and time-consuming. This study aimed to develop three denoising methods requiring different processing levels and MR-compatible hardware. At two time points, surface EMG was recorded from the lower leg of 6 participants (50:50 sex ratio, age = 26.24.6 yrs., height = 173.59.2 cm, weight = 71.511.4 kg) using a plantar flexion-based block design consisting of 30s of rest followed by 30s of flexion for 5 min, under three conditions: inside the MRI bore, with and without a BOLD sequence (3 T, BOLD sequence, GRE EPI, 10 slices, 64×64 matrix, 2 mm thickness, and TE/TR/flip = 35/3000 ms/70), and outside the MRI environment. Simultaneous BOLD/EMG recordings were denoised using average artifact subtraction with three methods of artifact template creation, each having varying timing and hardware requirements. Method M1 builds the artifact template by recording the scanner triggers coming from the MRI; M2 creates the artifact template with a constant artifact period computed as TR/[number of slices]; M3 estimates the artifact template by looking at the periodicity of the gradient artifact located in the EMG recordings. Following postprocessing, SNR analysis was performed, comparing rest-to-flexion periods, to assess the efficacy of denoising methods and to compare differences between conditions. Linear mixed-effects models showed no significant differences in the mean SNR between denoising methods (p = 0.656). Furthermore, EMG SNR measurements were significantly affected by the magnetic environment (p < 0.05) but not by muscle fatigue over time (p = 0.975). EMG recordings contaminated with gradient artifacts during simultaneous BOLD/EMG can be efficiently denoised using all proposed methods, with two methods requiring no extra hardware. With minimal post-processing, EMG can easily be performed during muscle BOLD MRI studies.

Brain state dynamics differ between eyes open and eyes closed rest

The human brain exhibits spatio-temporally complex activity even in the absence of external stimuli, cycling through recurring patterns of activity known as brain states. Thus far, brain state analysis has primarily been restricted to unimodal neuroimaging data sets, resulting in a limited definition of state and a poor understanding of the spatial and temporal relationships between states identified from different modalities. Here, we applied hidden Markov model (HMM) to concurrent electroencephalography-functional magnetic resonance imaging (EEG-fMRI) eyes open (EO) and eyes closed (EC) resting-state data, training models on the EEG and fMRI data separately, and evaluated the models' ability to distinguish dynamics between the two rest conditions. Additionally, we employed a general linear model approach to identify the BOLD correlates of the EEG-defined states to investigate whether the fMRI data could be used to improve the spatial definition of the EEG states. Finally, we performed a sliding window-based analysis on the state time courses to identify slower changes in the temporal dynamics, and then correlated these time courses across modalities. We found that both models could identify expected changes during EC rest compared to EO rest, with the fMRI model identifying changes in the activity and functional connectivity of visual and attention resting-state networks, while the EEG model correctly identified the canonical increase in alpha upon eye closure. In addition, by using the fMRI data, it was possible to infer the spatial properties of the EEG states, resulting in BOLD correlation maps resembling canonical alpha-BOLD correlations. Finally, the sliding window analysis revealed unique fractional occupancy dynamics for states from both models, with a selection of states showing strong temporal correlations across modalities. Overall, this study highlights the efficacy of using HMMs for brain state analysis, confirms that multimodal data can be used to provide more in-depth definitions of state and demonstrates that states defined across different modalities show similar temporal dynamics.

Heartbeat-evoked neural response abnormalities in generalized anxiety disorder during peripheral adrenergic stimulation

Hyperarousal symptoms in generalized anxiety disorder (GAD) are often incongruent with the observed physiological state, suggesting that abnormal processing of interoceptive signals is a characteristic feature of the disorder. To examine the neural mechanisms underlying interoceptive dysfunction in GAD, we evaluated whether adrenergic modulation of cardiovascular signaling differentially affects the heartbeat-evoked potential (HEP), an electrophysiological marker of cardiac interoception, during concurrent electroencephalogram and functional magnetic resonance imaging (EEG-fMRI) scanning. Intravenous infusions of the peripheral adrenergic agonist isoproterenol (0.5 and 2.0 micrograms, μg) were administered in a randomized, double-blinded and placebo-controlled fashion to dynamically perturb the cardiovascular system while recording the associated EEG-fMRI responses. During the 0.5 μg isoproterenol infusion, the GAD group (n = 24) exhibited significantly larger changes in HEP amplitude in an opposite direction than the healthy comparison (HC) group (n = 24). In addition, the GAD group showed significantly larger absolute HEP amplitudes than the HC group during saline infusions, when cardiovascular tone did not increase. No significant group differences in HEP amplitude were identified during the 2.0 μg isoproterenol infusion. Using analyzable blood oxygenation level-dependent fMRI data from participants with concurrent EEG-fMRI data (21 GAD and 21 HC), we found that the aforementioned HEP effects were uncorrelated with fMRI signals in the insula, ventromedial prefrontal cortex, dorsal anterior cingulate cortex, amygdala, and somatosensory cortex, brain regions implicated in cardiac signal processing in prior fMRI studies. These findings provide additional evidence of dysfunctional cardiac interoception in GAD and identify neural processes at the electrophysiological level that may be independent from blood oxygen level-dependent responses during peripheral adrenergic stimulation.

Simultaneous EEG-fMRI Investigation of Rhythm-Dependent Thalamo-Cortical Circuits Alteration in Schizophrenia

Schizophrenia is accompanied by aberrant interactions of intrinsic brain networks. However, the modulatory effect of electroencephalography (EEG) rhythms on the functional connectivity (FC) in schizophrenia remains unclear. This study aims to provide new insight into network communication in schizophrenia by integrating FC and EEG rhythm information. After collecting simultaneous resting-state EEG-functional magnetic resonance imaging data, the effect of rhythm modulations on FC was explored using what we term "dynamic rhythm information." We also investigated the synergistic relationships among three networks under rhythm modulation conditions, where this relationship presents the coupling between two brain networks with other networks as the center by the rhythm modulation. This study found FC between the thalamus and cortical network regions was rhythm-specific. Further, the effects of the thalamus on the default mode network (DMN) and salience network (SN) were less similar under alpha rhythm modulation in schizophrenia patients than in controls ([Formula: see text]). However, the similarity between the effects of the central executive network (CEN) on the DMN and SN under gamma modulation was greater ([Formula: see text]), and the degree of coupling was negatively correlated with the duration of disease ([Formula: see text], [Formula: see text]). Moreover, schizophrenia patients exhibited less coupling with the thalamus as the center and greater coupling with the CEN as the center. These results indicate that modulations in dynamic rhythms might contribute to the disordered functional interactions seen in schizophrenia.

Altered correlation of concurrently recorded EEG-fMRI connectomes in temporal lobe epilepsy

Whole-brain functional connectivity networks (connectomes) have been characterized at different scales in humans using EEG and fMRI. Multimodal epileptic networks have also been investigated, but the relationship between EEG and fMRI defined networks on a whole-brain scale is unclear. A unified multimodal connectome description, mapping healthy and pathological networks would close this knowledge gap. Here, we characterize the spatial correlation between the EEG and fMRI connectomes in right and left temporal lobe epilepsy (rTLE/lTLE). From two centers, we acquired resting-state concurrent EEG-fMRI of 35 healthy controls and 34 TLE patients. EEG-fMRI data was projected into the Desikan brain atlas, and functional connectomes from both modalities were correlated. EEG and fMRI connectomes were moderately correlated. This correlation was increased in rTLE when compared to controls for EEG-delta/theta/alpha/beta. Conversely, multimodal correlation in lTLE was decreased in respect to controls for EEG-beta. While the alteration was global in rTLE, in lTLE it was locally linked to the default mode network. The increased multimodal correlation in rTLE and decreased correlation in lTLE suggests a modality-specific lateralized differential reorganization in TLE, which needs to be considered when comparing results from different modalities. Each modality provides distinct information, highlighting the benefit of multimodal assessment in epilepsy.

Reliability of task-based fMRI in the dorsal horn of the human spinal cord

The application of functional magnetic resonance imaging (fMRI) to the human spinal cord is still a relatively small field of research and faces many challenges. Here we aimed to probe the limitations of task-based spinal fMRI at 3T by investigating the reliability of spinal cord blood oxygen level dependent (BOLD) responses to repeated nociceptive stimulation across two consecutive days in 40 healthy volunteers. We assessed the test-retest reliability of subjective ratings, autonomic responses, and spinal cord BOLD responses to short heat pain stimuli (1s duration) using the intraclass correlation coefficient (ICC). At the group level, we observed robust autonomic responses as well as spatially specific spinal cord BOLD responses at the expected location, but no spatial overlap in BOLD response patterns across days. While autonomic indicators of pain processing showed good-to-excellent reliability, both β-estimates and z-scores of task-related BOLD responses showed poor reliability across days in the target region (gray matter of the ipsilateral dorsal horn). When taking into account the sensitivity of gradient-echo echo planar imaging (GE-EPI) to draining vein signals by including the venous plexus in the analysis, we observed BOLD responses with fair reliability across days. Taken together, these results demonstrate that heat pain stimuli as short as one second are able to evoke a robust and spatially specific BOLD response, which is however strongly variable within participants across time, resulting in low reliability in the dorsal horn gray matter. Further improvements in data acquisition and analysis techniques are thus necessary before event-related spinal cord fMRI as used here can be reliably employed in longitudinal designs or clinical settings.

EEG decoding with spatiotemporal convolutional neural network for visualization and closed-loop control of sensorimotor activities: A simultaneous EEG-fMRI study

Closed-loop neurofeedback training utilizes neural signals such as scalp electroencephalograms (EEG) to manipulate specific neural activities and the associated behavioral performance. A spatiotemporal filter for high-density whole-head scalp EEG using a convolutional neural network can overcome the ambiguity of the signaling source because each EEG signal includes information on the remote regions. We simultaneously acquired EEG and functional magnetic resonance images in humans during the brain-computer interface (BCI) based neurofeedback training and compared the reconstructed and modeled hemodynamic responses of the sensorimotor network. Filters constructed with a convolutional neural network captured activities in the targeted network with spatial precision and specificity superior to those of the EEG signals preprocessed with standard pipelines used in BCI-based neurofeedback paradigms. The middle layers of the trained model were examined to characterize the neuronal oscillatory features that contributed to the reconstruction. Analysis of the layers for spatial convolution revealed the contribution of distributed cortical circuitries to reconstruction, including the frontoparietal and sensorimotor areas, and those of temporal convolution layers that successfully reconstructed the hemodynamic response function. Employing a spatiotemporal filter and leveraging the electrophysiological signatures of the sensorimotor excitability identified in our middle layer analysis would contribute to the development of a further effective neurofeedback intervention.

Uncovering Hemispheric Asymmetry and Directed Oscillatory Brain-Heart Interplay in Anxiety Processing: An fMRI Study

Brain-heart interactions (BHI) are critical for generating and processing emotions, including anxiety. Understanding specific neural correlates would be instrumental for greater comprehension and potential therapeutic interventions of anxiety disorders. While prior work has implicated the pontine structure as a central processor in cardiac regulation in anxiety, the distributed nature of anxiety processing across the cortex remains elusive. To address this, we performed a whole-brain-heart analysis using the full frequency directed transfer function to study resting-state spectral differences in BHI between high and low anxiety groups undergoing fMRI scans. Our findings revealed a hemispheric asymmetry in low-frequency interplay (0.05 Hz - 0.15 Hz) characterized by ascending BHI to the left insula and descending BHI from the right insula. Furthermore, we provide evidence supporting the "pacemaker hypothesis", highlighting the pons' function in regulating cardiac activity. Higher frequency interplay (0.2 Hz - 0.4Hz) demonstrate a preference for ascending interactions, particularly towards ventral prefrontal cortical activity in high anxiety groups, suggesting the heart's role in triggering a cognitive response to regulate anxiety. These findings highlight the impact of anxiety on BHI, contributing to a better understanding of its effect on the resting-state fMRI signal, with further implications for potential therapeutic interventions in treating anxiety disorders.

Striatum- and Cerebellum-Modulated Epileptic Networks Varying Across States with and without Interictal Epileptic Discharges

Idiopathic generalized epilepsy (IGE) is characterized by cryptogenic etiology and the striatum and cerebellum are recognized as modulators of epileptic network. We collected simultaneous electroencephalogram and functional magnetic resonance imaging data from 145 patients with IGE, 34 of whom recorded interictal epileptic discharges (IEDs) during scanning. In states without IEDs, hierarchical connectivity was performed to search core cortical regions which might be potentially modulated by striatum and cerebellum. Node-node and edge-edge moderation models were constructed to depict direct and indirect moderation effects in states with and without IEDs. Patients showed increased hierarchical connectivity with sensorimotor cortices (SMC) and decreased connectivity with regions in the default mode network (DMN). In the state without IEDs, striatum, cerebellum, and thalamus were linked to weaken the interactions of regions in the salience network (SN) with DMN and SMC. In periods with IEDs, overall increased moderation effects on the interaction between regions in SN and DMN, and between regions in DMN and SMC were observed. The thalamus and striatum were implicated in weakening interactions between regions in SN and SMC. The striatum and cerebellum moderated the cortical interaction among DMN, SN, and SMC in alliance with the thalamus, contributing to the dysfunction in states with and without IEDs in IGE. The current work revealed state-specific modulation effects of striatum and cerebellum on thalamocortical circuits and uncovered the potential core cortical targets which might contribute to develop new clinical neuromodulation techniques.

Multi-Algorithm Artifact Correction (MAAC) procedure part one: Algorithm and example

The Multi-Algorithm Artifact Correction (MAAC) procedure is presented for electroencephalographic (EEG) data, as made freely available in the open-source EP Toolkit (Dien, 2010). First the major EEG artifact correction methods (regression, spatial filters, principal components analysis, and independent components analysis) are reviewed. Contrary to the dominant approach of picking one method that is thought to be most effective, this review concludes that none are globally superior, but rather each has strengths and weaknesses. Then each of the major artifact types are reviewed (Blink, Corneo-Retinal Dipole, Saccadic Spike Potential, and Movement). For each one, it is proposed that one of the major correction methods is best matched to address it, resulting in the MAAC procedure. The MAAC itself is then presented, as implemented in the EP Toolkit, in order to provide a sense of the user experience. The primary goal of this present paper is to make the conceptual argument for the MAAC approach.

EEG as a potential ground truth for the assessment of cognitive state in software development activities: A multimodal imaging study

Cognitive human error and recent cognitive taxonomy on human error causes of software defects support the intuitive idea that, for instance, mental overload, attention slips, and working memory overload are important human causes for software bugs. In this paper, we approach the EEG as a reliable surrogate to MRI-based reference of the programmer's cognitive state to be used in situations where heavy imaging techniques are infeasible. The idea is to use EEG biomarkers to validate other less intrusive physiological measures, that can be easily recorded by wearable devices and useful in the assessment of the developer's cognitive state during software development tasks. Herein, our EEG study, with the support of fMRI, presents an extensive and systematic analysis by inspecting metrics and extracting relevant information about the most robust features, best EEG channels and the best hemodynamic time delay in the context of software development tasks. From the EEG-fMRI similarity analysis performed, we found significant correlations between a subset of EEG features and the Insula region of the brain, which has been reported as a region highly related to high cognitive tasks, such as software development tasks. We concluded that despite a clear inter-subject variability of the best EEG features and hemodynamic time delay used, the most robust and predominant EEG features, across all the subjects, are related to the Hjorth parameter Activity and Total Power features, from the EEG channels F4, FC4 and C4, and considering in most of the cases a hemodynamic time delay of 4 seconds used on the hemodynamic response function. These findings should be taken into account in future EEG-fMRI studies in the context of software debugging.

The control patterns of affective processing and cognitive reappraisal: insights from brain controllability analysis

Perceiving and modulating emotions is vital for cognitive function and is often impaired in neuropsychiatric conditions. Current tools for evaluating emotional dysregulation suffer from subjectivity and lack of precision, especially when it comes to understanding emotion from a regulatory or control-based perspective. To address these limitations, this study leverages an advanced methodology known as functional brain controllability analysis. We simultaneously recorded electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) data from 17 healthy subjects engaged in emotion processing and regulation tasks. We then employed a novel EEG/fMRI integration technique to reconstruct cortical activity in a high spatiotemporal resolution manner. Subsequently, we conducted functional brain controllability analysis to explore the neural network control patterns underlying different emotion conditions. Our findings demonstrated that the dorsolateral and ventrolateral prefrontal cortex exhibited increased controllability during the processing and regulation of negative emotions compared to processing of neutral emotion. Besides, the anterior cingulate cortex was notably more active in managing negative emotion than in either controlling neutral emotion or regulating negative emotion. Finally, the posterior parietal cortex emerged as a central network controller for the regulation of negative emotion. This study offers valuable insights into the cortical control mechanisms that support emotion perception and regulation.

Exploring changes in brain function in IBD patients using SPCCA: a study of simultaneous EEG-fMRI

Research on functional changes in the brain of inflammatory bowel disease (IBD) patients is emerging around the world, which brings new perspectives to medical research. In this paper, the methods of canonical correlation analysis (CCA), kernel canonical correlation analysis (KCCA), and sparsity preserving canonical correlation analysis (SPCCA) were applied to the fusion of simultaneous EEG-fMRI data from 25 IBD patients and 15 healthy individuals. The CCA, KCCA and SPCCA fusion methods were used for data processing to compare the results obtained by the three methods. The results clearly show that there is a significant difference in the activation intensity between IBD and healthy control (HC), not only in the frontal lobe (p < 0.01) and temporal lobe (p < 0.01) regions, but also in the posterior cingulate gyrus (p < 0.01), gyrus rectus (p < 0.01), and amygdala (p < 0.01) regions, which are usually neglected. The mean difference in the SPCCA activation intensity was 60.1. However, the mean difference in activation intensity was only 36.9 and 49.8 by using CCA and KCCA. In addition, the correlation of the relevant components selected during the SPCCA calculation was high, with correlation components of up to 0.955; alternatively, the correlations obtained from CCA and KCCA calculations were only 0.917 and 0.926, respectively. It can be seen that SPCCA is indeed superior to CCA and KCCA in processing high-dimensional multimodal data. This work reveals the process of analyzing the brain activation state in IBD disease, provides a further perspective for the study of brain function, and opens up a new avenue for studying the SPCCA method and the change in the intensity of brain activation in IBD disease.

Temporal cluster-based organization of sleep spindles underlies motor memory consolidation

Sleep benefits motor memory consolidation, which is mediated by sleep spindle activity and associated memory reactivations during non-rapid eye movement (NREM) sleep. However, the particular role of NREM2 and NREM3 sleep spindles and the mechanisms triggering this memory consolidation process remain unclear. Here, simultaneous electroencephalographic and functional magnetic resonance imaging (EEG-fMRI) recordings were collected during night-time sleep following the learning of a motor sequence task. Adopting a time-based clustering approach, we provide evidence that spindles iteratively occur within clustered and temporally organized patterns during both NREM2 and NREM3 sleep. However, the clustering of spindles in trains is related to motor memory consolidation during NREM2 sleep only. Altogether, our findings suggest that spindles' clustering and rhythmic occurrence during NREM2 sleep may serve as an intrinsic rhythmic sleep mechanism for the timed reactivation and subsequent consolidation of motor memories, through synchronized oscillatory activity within a subcortical-cortical network involved during learning.

Maturation-dependent changes in cortical and thalamic activity during sleep slow waves: Insights from a combined EEG-fMRI study

INTRODUCTION

Studies using scalp EEG have shown that slow waves (0.5-4 Hz), the most prominent hallmark of NREM sleep, undergo relevant changes from childhood to adulthood, mirroring brain structural modifications and the acquisition of cognitive skills. Here we used simultaneous EEG-fMRI to investigate the cortical and subcortical correlates of slow waves in school-age children and determine their relative developmental changes.

METHODS

We analyzed data from 14 school-age children with self-limited focal epilepsy of childhood who fell asleep during EEG-fMRI recordings. Brain regions associated with slow-wave occurrence were identified using a voxel-wise regression that also modelled interictal epileptic discharges and sleep spindles. At the group level, a mixed-effects linear model was used. The results were qualitatively compared with those obtained from 2 adolescents with epilepsy and 17 healthy adults.

RESULTS

Slow waves were associated with hemodynamic-signal decreases in bilateral somatomotor areas. Such changes extended more posteriorly relative to those in adults. Moreover, the involvement of areas belonging to the default mode network changes as a function of age. No significant hemodynamic responses were observed in subcortical structures. However, we identified a significant correlation between age and thalamic hemodynamic changes.

CONCLUSIONS

Present findings indicate that the somatomotor cortex may have a key role in slow-wave expression throughout the lifespan. At the same time, they are consistent with a posterior-to-anterior shift in slow-wave distribution mirroring brain maturational changes. Finally, our results suggest that slow-wave changes may not reflect only neocortical modifications but also the maturation of subcortical structures, including the thalamus.

Heartbeat-evoked neural response abnormalities in generalized anxiety disorder during peripheral adrenergic stimulation

Hyperarousal symptoms in generalized anxiety disorder (GAD) are often incongruent with the observed physiological state, suggesting that abnormal processing of interoceptive signals is a characteristic feature of the disorder. To examine the neural mechanisms underlying interoceptive dysfunction in GAD, we evaluated whether adrenergic modulation of cardiovascular signaling differentially affects the heartbeat evoked potential (HEP), an electrophysiological marker of cardiac interoception, during concurrent electroencephalogram and functional magnetic resonance imaging (EEG-fMRI) scanning. Intravenous infusions of the peripheral adrenergic agonist isoproterenol (0.5 and 2.0 micrograms, μg) were administered in a randomized, double-blinded and placebo-controlled fashion to dynamically perturb the cardiovascular system while recording the associated EEG-fMRI responses. During the 0.5 μg isoproterenol infusion, the GAD group (n=24) exhibited significantly larger changes in HEP amplitude in an opposite direction than the HC group (n=24). In addition, the GAD group showed significantly larger absolute HEP amplitudes than HC during saline infusions, when cardiovascular tone did not increase. No significant group differences in HEP amplitude were identified during the 2.0 μg isoproterenol infusion. Using analyzable blood oxygenation level dependent fMRI data from participants with concurrent EEG-fMRI data (21 GAD and 21 HC), we found that the aforementioned HEP effects were uncorrelated with fMRI signals in the insula, ventromedial prefrontal cortex, dorsal anterior cingulate cortex, amygdala, and somatosensory cortex, brain regions implicated in cardiac signal processing according to prior fMRI studies. These findings provide additional evidence of dysfunctional cardiac interoception in GAD and identify neural processes at the electrophysiological level that may be independent from blood oxygen level-dependent responses during peripheral adrenergic stimulation.

Alterations of sleep initiation in NREM parasomnia after sleep deprivation - A multimodal pilot study

Objectives

NREM parasomnias also known as disorders of arousal (DOA) are characterised by abnormal motor and autonomic activation during arousals primarily from slow wave sleep. Dissociative state between sleep and wake is likely responsible for clinical symptoms of DOA. We therefore investigated potential dissociation outside of parasomnic events by using simultaneous 256-channel EEG (hdEEG) and functional magnetic resonance imaging (fMRI).

Methods

Eight DOA patients (3 women, mean age = 27.8; SD = 4.2) and 8 gender and age matched healthy volunteers (3 women, mean age = 26,5; SD = 4.0) were included into the study. They underwent 30-32 h of sleep deprivation followed by hdEEG and fMRI recording. We determined 2 conditions: falling asleep (FA) and arousal (A), that occurred outside of deep sleep and/or parasomnic event. We used multimodal approach using data obtained from EEG, fMRI and EEG-fMRI integration approach.

Results

DOA patients showed increase in delta and beta activity over postcentral gyrus and cuneus during awakening period. This group expressed increased connectivity between motor cortex and cingulate during arousals unrelated to parasomnic events in the beta frequency band. They also showed lower connectivity between different portions of cingulum. In contrast, the greater connectivity was found between thalamus and some cortical areas, such as occipital cortex.

Conclusion

Our findings suggest a complex alteration in falling asleep and arousal mechanisms at both subcortical and cortical levels in response to sleep deprivation. As this alteration is present also outside of slow wave sleep and/or parasomnic episodes we believe this could be a trait factor of DOA.

Broadly applicable methods for the detection of artefacts in electroencephalography acquired simultaneously with hemodynamic recordings

Electroencephalography (EEG) data, acquired simultaneously with magnetic resonance imaging (MRI), must be corrected for artefacts related to MR gradient switches (GS) and the cardioballistic (CB) effect. Canonical approaches require additional signal acquisition for artefact detection (e.g., MR volume onsets, ECG), without which the EEG data would be rendered uncleanable from these artefacts.•We present two broadly applicable methods for artefact detection based on peak detection combined with temporal constraints with respect to periodicity directly from the EEG data itself; no additional signals are required. We validated the performance of our methods versus the two canonical approaches for detection of GS/CB artefact, respectively, on 26 healthy human EEG-functional MRI resting-state datasets. Utilising various performance metrics, we found our methods to perform as well as - and sometimes better than - the canonical standard approaches. With as little as one EEG channel recording, our methods can be applied to detect GS/CB artefacts in EEG data acquired simultaneously with MRI in the absence of MR volume onsets and/or an ECG recording. The detected artefact onsets can then be fed into the standard artefact correction software.

NeuXus open-source tool for real-time artifact reduction in simultaneous EEG-fMRI

The simultaneous acquisition of electroencephalography and functional magnetic resonance imaging (EEG-fMRI) allows the complementary study of the brain's electrophysiology and hemodynamics with high temporal and spatial resolution. One application with great potential is neurofeedback training of targeted brain activity, based on the real-time analysis of the EEG and/or fMRI signals. This depends on the ability to reduce in real time the severe artifacts affecting the EEG signal acquired with fMRI, mainly the gradient and pulse artifacts. A few methods have been proposed for this purpose, but they are either slow, hardware-dependent, publicly unavailable, or proprietary software. Here, we present a fully open-source and publicly available tool for real-time EEG artifact reduction in simultaneous EEG-fMRI recordings that is fast and applicable to any hardware. Our tool is integrated in the Python toolbox NeuXus for real-time EEG processing and adapts to a real-time scenario well-established artifact average subtraction methods combined with a long short-term memory network for R peak detection. We benchmarked NeuXus on three different datasets, in terms of artifact power reduction and background signal preservation in resting state, alpha-band power reactivity to eyes closure, and event-related desynchronization during motor imagery. We showed that NeuXus performed at least as well as the only available real-time tool for conventional hardware setups (BrainVision's RecView) and a well-established offline tool (EEGLAB's FMRIB plugin). We also demonstrated NeuXus' real-time ability by reporting execution times under 250 ms. In conclusion, we present and validate the first fully open-source and hardware-independent solution for real-time artifact reduction in simultaneous EEG-fMRI studies.

Neural modulation enhancement using connectivity-based EEG neurofeedback with simultaneous fMRI for emotion regulation

Emotion regulation plays a key role in human behavior and overall well-being. Neurofeedback is a non-invasive self-brain training technique used for emotion regulation to enhance brain function and treatment of mental disorders through behavioral changes. Previous neurofeedback research often focused on using activity from a single brain region as measured by fMRI or power from one or two EEG electrodes. In a new study, we employed connectivity-based EEG neurofeedback through recalling positive autobiographical memories and simultaneous fMRI to upregulate positive emotion. In our novel approach, the feedback was determined by the coherence of EEG electrodes rather than the power of one or two electrodes. We compared the efficiency of this connectivity-based neurofeedback to traditional activity-based neurofeedback through multiple experiments. The results showed that connectivity-based neurofeedback effectively improved BOLD signal change and connectivity in key emotion regulation regions such as the amygdala, thalamus, and insula, and increased EEG frontal asymmetry, which is a biomarker for emotion regulation and treatment of mental disorders such as PTSD, anxiety, and depression and coherence among EEG channels. The psychometric evaluations conducted both before and after the neurofeedback experiments revealed that participants demonstrated improvements in enhancing positive emotions and reducing negative emotions when utilizing connectivity-based neurofeedback, as compared to traditional activity-based and sham neurofeedback approaches. These findings suggest that connectivity-based neurofeedback may be a superior method for regulating emotions and could be a useful alternative therapy for mental disorders, providing individuals with greater control over their brain and mental functions.

Comparison of Resting-State EEG Network Analyses With and Without Parallel MRI in Genetic Generalized Epilepsy

Genetic generalized epilepsy (GGE) is conceptualized as a brain disorder involving distributed bilateral networks. To study these networks, simultaneous EEG-fMRI measurements can be used. However, inside-MRI EEG suffers from strong MR-related artifacts; it is not established whether EEG-based metrics in EEG-fMRI resting-state measurements are suitable for the analysis of group differences at source-level. We evaluated the impact of the inside-MR measurement condition on statistical group comparisons of EEG on source-level power and functional connectivity in patients with GGE versus healthy controls. We studied the cross-modal spatial relation of statistical group differences in seed-based FC derived from EEG and parallel fMRI. We found a significant increase in power and a frequency-specific change in functional connectivity for the inside MR-scanner compared to the outside MR-scanner condition. For power, we found reduced group difference between GGE and controls both in terms of statistical significance as well as effect size. Group differences for ImCoh remained similar both in terms of statistical significance as well as effect size. We found increased seed-based FC for GGE patients from the thalamus to the precuneus cortex region in fMRI, and in the theta band of simultaneous EEG. Our findings suggest that the analysis of EEG functional connectivity based on ImCoh is suitable for MR-EEG, and that relative group difference in a comparison of patients with GGE against controls are preserved. Spatial correspondence of seed-based FC group differences between the two modalities was found for the thalamus.

An open-access dataset of naturalistic viewing using simultaneous EEG-fMRI

In this work, we present a dataset that combines functional magnetic imaging (fMRI) and electroencephalography (EEG) to use as a resource for understanding human brain function in these two imaging modalities. The dataset can also be used for optimizing preprocessing methods for simultaneously collected imaging data. The dataset includes simultaneously collected recordings from 22 individuals (ages: 23-51) across various visual and naturalistic stimuli. In addition, physiological, eye tracking, electrocardiography, and cognitive and behavioral data were collected along with this neuroimaging data. Visual tasks include a flickering checkerboard collected outside and inside the MRI scanner (EEG-only) and simultaneous EEG-fMRI recordings. Simultaneous recordings include rest, the visual paradigm Inscapes, and several short video movies representing naturalistic stimuli. Raw and preprocessed data are openly available to download. We present this dataset as part of an effort to provide open-access data to increase the opportunity for discoveries and understanding of the human brain and evaluate the correlation between electrical brain activity and blood oxygen level-dependent (BOLD) signals.

Enhancing precision in human neuroscience

Human neuroscience has always been pushing the boundary of what is measurable. During the last decade, concerns about statistical power and replicability - in science in general, but also specifically in human neuroscience - have fueled an extensive debate. One important insight from this discourse is the need for larger samples, which naturally increases statistical power. An alternative is to increase the precision of measurements, which is the focus of this review. This option is often overlooked, even though statistical power benefits from increasing precision as much as from increasing sample size. Nonetheless, precision has always been at the heart of good scientific practice in human neuroscience, with researchers relying on lab traditions or rules of thumb to ensure sufficient precision for their studies. In this review, we encourage a more systematic approach to precision. We start by introducing measurement precision and its importance for well-powered studies in human neuroscience. Then, determinants for precision in a range of neuroscientific methods (MRI, M/EEG, EDA, Eye-Tracking, and Endocrinology) are elaborated. We end by discussing how a more systematic evaluation of precision and the application of respective insights can lead to an increase in reproducibility in human neuroscience.

Reliability of resting-state functional connectivity in the human spinal cord: Assessing the impact of distinct noise sources

The investigation of spontaneous fluctuations of the blood-oxygen-level-dependent (BOLD) signal has recently been extended from the brain to the spinal cord, where it has stimulated interest from a clinical perspective. A number of resting-state functional magnetic resonance imaging (fMRI) studies have demonstrated robust functional connectivity between the time series of BOLD fluctuations in bilateral dorsal horns and between those in bilateral ventral horns, in line with the functional neuroanatomy of the spinal cord. A necessary step prior to extension to clinical studies is assessing the reliability of such resting-state signals, which we aimed to do here in a group of 45 healthy young adults at the clinically prevalent field strength of 3T. When investigating connectivity in the entire cervical spinal cord, we observed fair to good reliability for dorsal-dorsal and ventral-ventral connectivity, whereas reliability was poor for within- and between-hemicord dorsal-ventral connectivity. Considering how prone spinal cord fMRI is to noise, we extensively investigated the impact of distinct noise sources and made two crucial observations: removal of physiological noise led to a reduction in functional connectivity strength and reliability - due to the removal of stable and participant-specific noise patterns - whereas removal of thermal noise considerably increased the detectability of functional connectivity without a clear influence on reliability. Finally, we also assessed connectivity within spinal cord segments and observed that while the pattern of connectivity was similar to that of whole cervical cord, reliability at the level of single segments was consistently poor. Taken together, our results demonstrate the presence of reliable resting-state functional connectivity in the human spinal cord even after thoroughly accounting for physiological and thermal noise, but at the same time urge caution if focal changes in connectivity (e.g. due to segmental lesions) are to be studied, especially in a longitudinal manner.

EEG-LLAMAS: A low-latency neurofeedback platform for artifact reduction in EEG-fMRI

Simultaneous EEG-fMRI is a powerful multimodal technique for imaging the brain, but its use in neurofeedback experiments has been limited by EEG noise caused by the MRI environment. Neurofeedback studies typically require analysis of EEG in real time, but EEG acquired inside the scanner is heavily contaminated with ballistocardiogram (BCG) artifact, a high-amplitude artifact locked to the cardiac cycle. Although techniques for removing BCG artifacts do exist, they are either not suited to real-time, low-latency applications, such as neurofeedback, or have limited efficacy. We propose and validate a new open-source artifact removal software called EEG-LLAMAS (Low Latency Artifact Mitigation Acquisition Software), which adapts and advances existing artifact removal techniques for low-latency experiments. We first used simulations to validate LLAMAS in data with known ground truth. We found that LLAMAS performed better than the best publicly-available real-time BCG removal technique, optimal basis sets (OBS), in terms of its ability to recover EEG waveforms, power spectra, and slow wave phase. To determine whether LLAMAS would be effective in practice, we then used it to conduct real-time EEG-fMRI recordings in healthy adults, using a steady state visual evoked potential (SSVEP) task. We found that LLAMAS was able to recover the SSVEP in real time, and recovered the power spectra collected outside the scanner better than OBS. We also measured the latency of LLAMAS during live recordings, and found that it introduced a lag of less than 50 ms on average. The low latency of LLAMAS, coupled with its improved artifact reduction, can thus be effectively used for EEG-fMRI neurofeedback. A limitation of the method is its use of a reference layer, a piece of EEG equipment which is not commercially available, but can be assembled in-house. This platform enables closed-loop experiments which previously would have been prohibitively difficult, such as those that target short-duration EEG events, and is shared openly with the neuroscience community.

Analysis of Respiratory Sinus Arrhythmia and Directed Information Flow between Brain and Body Indicate Different Management Strategies of fMRI-Related Anxiety

Background: Respiratory sinus arrhythmia (RSA) denotes decrease of cardiac beat-to-beat intervals (RRI) during inspiration and RRI increase during expiration, but an inverse pattern (termed negative RSA) was also found in healthy humans with elevated anxiety. It was detected using wave-by-wave analysis of cardiorespiratory rhythms and was considered to reflect a strategy of anxiety management involving the activation of a neural pacemaker. Results were consistent with slow breathing, but contained uncertainty at normal breathing rates (0.2–0.4 Hz). Objectives and methods: We combined wave-by-wave analysis and directed information flow analysis to obtain information on anxiety management at higher breathing rates. We analyzed cardiorespiratory rhythms and blood oxygen level-dependent (BOLD) signals from the brainstem and cortex in 10 healthy fMRI participants with elevated anxiety. Results: Three subjects with slow respiratory, RRI, and neural BOLD oscillations showed 57 ± 26% negative RSA and significant anxiety reduction by 54 ± 9%. Six participants with breathing rate of ~0.3 Hz showed 41 ± 16% negative RSA and weaker anxiety reduction. They presented significant information flow from RRI to respiration and from the middle frontal cortex to the brainstem, which may result from respiration-entrained brain oscillations, indicating another anxiety management strategy. Conclusion: The two analytical approaches applied here indicate at least two different anxiety management strategies in healthy subjects.

Exploring Neural Mechanisms of Reward Processing Using Coupled Matrix Tensor Factorization: A Simultaneous EEG-fMRI Investigation

BACKGROUND

It is crucial to understand the neural feedback mechanisms and the cognitive decision-making of the brain during the processing of rewards. Here, we report the first attempt for a simultaneous electroencephalography (EEG)-functional magnetic resonance imaging (fMRI) study in a gambling task by utilizing tensor decomposition.

METHODS

First, the single-subject EEG data are represented as a third-order spectrogram tensor to extract frequency features. Next, the EEG and fMRI data are jointly decomposed into a superposition of multiple sources characterized by space-time-frequency profiles using coupled matrix tensor factorization (CMTF). Finally, graph-structured clustering is used to select the most appropriate model according to four quantitative indices.

RESULTS

The results clearly show that not only are the regions of interest (ROIs) found in other literature activated, but also the olfactory cortex and fusiform gyrus which are usually ignored. It is found that regions including the orbitofrontal cortex and insula are activated for both winning and losing stimuli. Meanwhile, regions such as the superior orbital frontal gyrus and anterior cingulate cortex are activated upon winning stimuli, whereas the inferior frontal gyrus, cingulate cortex, and medial superior frontal gyrus are activated upon losing stimuli.

CONCLUSION

This work sheds light on the reward-processing progress, provides a deeper understanding of brain function, and opens a new avenue in the investigation of neurovascular coupling via CMTF.

Unbalance between working memory task-activation and task-deactivation networks in epilepsy: Simultaneous EEG-fMRI study

Working memory (WM) is a cognitive function involving emergent properties of theta oscillations and large-scale network interactions. The synchronization of WM task-related networks in the brain enhanced WM performance. However, how these networks regulate WM processing is not well known, and the alteration of the interaction among these networks may play an important role in patients with cognitive dysfunction. In this study, we used simultaneous EEG-fMRI to examine the features of theta oscillations and the functional interactions among activation/deactivation networks during the n-back WM task in patients with idiopathic generalized epilepsy (IGE). The results showed that there was more enhancement of frontal theta power along with WM load increase in IGE, and the theta power was positively correlated with the accuracy of the WM tasks. Moreover, fMRI activations/deactivations correlated with n-back tasks were estimated, and we found that the IGE group had increased and widespread activations in high-load WM tasks, including the frontoparietal activation network and task-related deactivation areas, such as the default mode network and primary visual and auditory networks. In addition, the network connectivity results demonstrated decreased counteraction between the activation network and deactivation network, and the counteraction was correlated with the higher theta power in IGE. These results indicated the important role of the interactions between activation and deactivation networks during the WM process, and the unbalance among them may indicate the pathophysiological mechanism of cognitive dysfunction in generalized epilepsy.

Respiration-entrained brain oscillations in healthy fMRI participants with high anxiety

Brain-body interactions can be studied by using directed coupling measurements of fMRI oscillations in the low (0.1-0.2 Hz) and high frequency bands (HF; 0.2-0.4 Hz). Recently, a preponderance of oscillations in the information flow between the brainstem and the prefrontal cortex at around 0.15/0.16 Hz was shown. The goal of this study was to investigate the information flow between BOLD-, respiratory-, and heart beat-to-beat interval (RRI) signals in the HF band in healthy subjects with high anxiety during fMRI examinations. A multivariate autoregressive model was concurrently applied to the BOLD signals from the middle frontal gyrus (MFG), precentral gyrus and the brainstem, as well as to respiratory and RRI signals. Causal coupling between all signals was determined using the Directed Transfer Function (DTF). We found a salience of fast respiratory waves with a period of 3.1 s (corresponding to ~ 0.32 Hz) and a highly significant (p < 0.001) top-down information-flow from BOLD oscillations in the MFG to the brainstem. Additionally, there was a significant (p < 0.01) information flow from RRI to respiratory oscillations. We speculate that brain oscillations around 0.32 Hz, triggered by nasal breathing, are projected downwards to the brainstem. Particularly interesting is the driving force of cardiac to respiratory waves with a ratio of 1:1 or 1:2. These results support the binary hierarchy model with preferred respiratory frequencies at 0.32 Hz and 0.16 Hz.

Probing fMRI brain connectivity and activity changes during emotion regulation by EEG neurofeedback

Despite the existence of several emotion regulation studies using neurofeedback, interactions among a small number of regions were evaluated, and therefore, further investigation is needed to understand the interactions of the brain regions involved in emotion regulation. We implemented electroencephalography (EEG) neurofeedback with simultaneous functional magnetic resonance imaging (fMRI) using a modified happiness-inducing task through autobiographical memories to upregulate positive emotion. Then, an explorative analysis of whole brain regions was done to understand the effect of neurofeedback on brain activity and the interaction of whole brain regions involved in emotion regulation. The participants in the control and experimental groups were asked to do emotion regulation while viewing positive images of autobiographical memories and getting sham or real (based on alpha asymmetry) EEG neurofeedback, respectively. The proposed multimodal approach quantified the effects of EEG neurofeedback in changing EEG alpha power, fMRI blood oxygenation level-dependent (BOLD) activity of prefrontal, occipital, parietal, and limbic regions (up to 1.9% increase), and functional connectivity in/between prefrontal, parietal, limbic system, and insula in the experimental group. New connectivity links were identified by comparing the brain functional connectivity between experimental conditions (Upregulation and View blocks) and also by comparing the brain connectivity of the experimental and control groups. Psychometric assessments confirmed significant changes in positive and negative mood states in the experimental group by neurofeedback. Based on the exploratory analysis of activity and connectivity among all brain regions involved in emotion regions, we found significant BOLD and functional connectivity increases due to EEG neurofeedback in the experimental group, but no learning effect was observed in the control group. The results reveal several new connections among brain regions as a result of EEG neurofeedback which can be justified according to emotion regulation models and the role of those regions in emotion regulation and recalling positive autobiographical memories.

Preservation of EEG spectral power features during simultaneous EEG-fMRI

Introduction

Electroencephalographic (EEG) data quality is severely compromised when recorded inside the magnetic resonance (MR) environment. Here we characterized the impact of the ballistocardiographic (BCG) artifact on resting-state EEG spectral properties and compared the effectiveness of seven common BCG correction methods to preserve EEG spectral features. We also assessed if these methods retained posterior alpha power reactivity to an eyes closure-opening (EC-EO) task and compared the results from EEG-informed fMRI analysis using different BCG correction approaches.

Method

Electroencephalographic data from 20 healthy young adults were recorded outside the MR environment and during simultaneous fMRI acquisition. The gradient artifact was effectively removed from EEG-fMRI acquisitions using Average Artifact Subtraction (AAS). The BCG artifact was corrected with seven methods: AAS, Optimal Basis Set (OBS), Independent Component Analysis (ICA), OBS followed by ICA, AAS followed by ICA, PROJIC-AAS and PROJIC-OBS. EEG signal preservation was assessed by comparing the spectral power of traditional frequency bands from the corrected rs-EEG-fMRI data with the data recorded outside the scanner. We then assessed the preservation of posterior alpha functional reactivity by computing the ratio between the EC and EO conditions during the EC-EO task. EEG-informed fMRI analysis of the EC-EO task was performed using alpha power-derived BOLD signal predictors obtained from the EEG signals corrected with different methods.

Results

The BCG artifact caused significant distortions (increased absolute power, altered relative power) across all frequency bands. Artifact residuals/signal losses were present after applying all correction methods. The EEG reactivity to the EC-EO task was better preserved with ICA-based correction approaches, particularly when using ICA feature extraction to isolate alpha power fluctuations, which allowed to accurately predict hemodynamic signal fluctuations during the EEG-informed fMRI analysis.

Discussion

Current software solutions for the BCG artifact problem offer limited efficiency to preserve the EEG spectral power properties using this particular EEG setup. The state-of-the-art approaches tested here can be further refined and should be combined with hardware implementations to better preserve EEG signal properties during simultaneous EEG-fMRI. Existing and novel BCG artifact correction methods should be validated by evaluating signal preservation of both ERPs and spontaneous EEG spectral power.

Sleep, learning, and memory in human research using noninvasive neuroimaging techniques

An accumulating body of evidence indicates that sleep is beneficial for learning and memory. Task performance improves significantly after a period that includes sleep, whereas a lack of sleep nullifies or impairs such improvements. Our current knowledge about sleep's role in learning and memory has been obtained based on studies that were conducted in both animal models and human subjects. Nevertheless, how sleep promotes learning and memory in humans is not fully understood. In this review, we overview our current understating of how sleep may contribute to learning and memory, covering different roles of non-rapid eye movement and rapid eye movement sleep. We then discuss cutting-edge advanced techniques that are currently available, including simultaneous functional magnetic resonance imaging (fMRI) and electroencephalogram (EEG) and simultaneous functional magnetic resonance spectroscopy (fMRS) and EEG measurements, and evaluate how these may contribute to advance the understanding of the role of sleep in human cognition. We also highlight the current limitations and challenges using these methods and discuss ways that may allow us to overcome these limitations.

Individual slow wave events give rise to macroscopic fMRI signatures and drive the strength of the BOLD signal in human resting-state EEG-fMRI recordings

The slow wave state is a general state of quiescence interrupted by sudden bursts of activity or so-called slow wave events (SWEs). Recently, the relationship between SWEs and blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) signals was assessed in rodent models which revealed cortex-wide BOLD activation. However, it remains unclear which macroscopic signature corresponds to these specific neurophysiological events in the human brain. Therefore, we analyzed simultaneous electroencephalographic (EEG)-fMRI data during human non-REM sleep. SWEs individually detected in the EEG data were used as predictors in event-related fMRI analyses to examine the relationship between SWEs and fMRI signals. For all 10 subjects we identified significant changes in BOLD activity associated with SWEs covering substantial parts of the gray matter. As demonstrated in rodents, we observed a direct relation of a neurophysiological event to specific BOLD activation patterns. We found a correlation between the number of SWEs and the spatial extent of these BOLD activation patterns and discovered that the amplitude of the BOLD response strongly depends on the SWE amplitude. As altered SWE propagation has recently been found in neuropsychiatric diseases, it is critical to reveal the brain's physiological slow wave state networks to potentially establish early imaging biomarkers for various diseases long before disease onset.

Hemodynamic functional connectivity optimization of frequency EEG microstates enables attention LSTM framework to classify distinct temporal cortical communications of different cognitive tasks

Temporal analysis of global cortical communication of cognitive tasks in coarse EEG information is still challenging due to the underlying complex neural mechanisms. This study proposes an attention-based time-series deep learning framework that processes fMRI functional connectivity optimized quasi-stable frequency microstates for classifying distinct temporal cortical communications of the cognitive task. Seventy volunteers were subjected to visual target detection tasks, and their electroencephalogram (EEG) and functional MRI (fMRI) were acquired simultaneously. At first, the acquired EEG information was preprocessed and bandpass to delta, theta, alpha, beta, and gamma bands and then subjected to quasi-stable frequency-microstate estimation. Subsequently, time-series elicitation of each frequency microstates is optimized with graph theory measures of simultaneously eliciting fMRI functional connectivity between frontal, parietal, and temporal cortices. The distinct neural mechanisms associated with each optimized frequency-microstate were analyzed using microstate-informed fMRI. Finally, these optimized, quasi-stable frequency microstates were employed to train and validate the attention-based Long Short-Term Memory (LSTM) time-series architecture for classifying distinct temporal cortical communications of the target from other cognitive tasks. The temporal, sliding input sampling windows were chosen between 180 to 750 ms/segment based on the stability of transition probabilities of the optimized microstates. The results revealed 12 distinct frequency microstates capable of deciphering target detections' temporal cortical communications from other task engagements. Particularly, fMRI functional connectivity measures of target engagement were observed significantly correlated with the right-diagonal delta (r = 0.31), anterior-posterior theta (r = 0.35), left-right theta (r = - 0.32), alpha (r = - 0.31) microstates. Further, neuro-vascular information of microstate-informed fMRI analysis revealed the association of delta/theta and alpha/beta microstates with cortical communications and local neural processing, respectively. The classification accuracies of the attention-based LSTM were higher than the traditional LSTM architectures, particularly the frameworks that sampled the EEG data with a temporal width of 300 ms/segment. In conclusion, the study demonstrates reliable temporal classifications of global cortical communication of distinct tasks using an attention-based LSTM utilizing fMRI functional connectivity optimized quasi-stable frequency microstates.

Involvement of the Thalamus, Hippocampus, and Brainstem in Hypsarrhythmia of West Syndrome: Simultaneous Recordings of Electroencephalography and fMRI Study

BACKGROUND AND PURPOSE

West syndrome is a developmental and epileptic encephalopathy characterized by epileptic spasms, neurodevelopmental regression, and a specific EEG pattern called hypsarrhythmia. Our aim was to investigate the brain activities related to hypsarrhythmia at onset and focal epileptiform discharges in the remote period in children with West syndrome using simultaneous electroencephalography and fMRI recordings.

MATERIALS AND METHODS

Fourteen children with West syndrome underwent simultaneous electroencephalography and fMRI at the onset of West syndrome. Statistically significant blood oxygen level-dependent responses related to hypsarrhythmia were analyzed using an event-related design of 4 hemodynamic response functions with peaks at 3, 5, 7, and 9 seconds after the onset of each event. Six of 14 children had focal epileptiform discharges after treatment and underwent simultaneous electroencephalography and fMRI from 12 to 25 months of age.

RESULTS

At onset, positive blood oxygen level-dependent responses were seen in the brainstem (14/14 patients), thalami (13/14), basal ganglia (13/14), and hippocampi (13/14), in addition to multiple cerebral cortices. Group analysis using hemodynamic response functions with peaks at 3, 5, and 7 seconds showed positive blood oxygen level-dependent responses in the brainstem, thalamus, and hippocampus, while positive blood oxygen level-dependent responses in multiple cerebral cortices were seen using hemodynamic response functions with peaks at 5 and 7 seconds. In the remote period, 3 of 6 children had focal epileptiform discharge-related positive blood oxygen level-dependent responses in the thalamus, hippocampus, and brainstem.

CONCLUSIONS

Positive blood oxygen level-dependent responses with hypsarrhythmia appeared in the brainstem, thalamus, and hippocampus on earlier hemodynamic response functions than the cerebral cortices, suggesting the propagation of epileptogenic activities from the deep brain structures to the neocortices. Activation of the hippocampus, thalamus, and brainstem was still seen in half of the patients with focal epileptiform discharges after adrenocorticotropic hormone therapy.

Brain-wide functional connectivity of face patch neurons during rest

In resting state functional magnetic resonance imaging (fMRI), areas showing coherent hemodynamic fluctuations across the brain are operationally defined to be functionally connected. However, it is unknown how the activity of single units residing within a voxel contributes to this network structure. Here we demonstrate a shared but restricted pattern of functional connectivity among neighboring neurons residing in functionally defined face patches. Unexpectedly, such neurons also exhibited a prominent inverse correlation with thalamic structures and brainstem neuromodulatory centers. Single unit maps differed from analogous maps obtained with local field potentials and seed-based fMRI. These findings suggest that during rest, individual cortical neurons have a restricted set of functional connections, which is governed in part by anatomical projections and in part by neuromodulation.

The brain is a highly organized, dynamic system whose network architecture is often assessed through resting functional magnetic resonance imaging (fMRI) functional connectivity. The functional interactions between brain areas, including those observed during rest, are assumed to stem from the collective influence of action potentials carried by long-range neural projections. However, the contribution of individual neurons to brain-wide functional connectivity has not been systematically assessed. Here we developed a method to concurrently measure and compare the spiking activity of local neurons with fMRI signals measured across the brain during rest. We recorded spontaneous activity from neural populations in cortical face patches in the macaque during fMRI scanning sessions. Individual cells exhibited prominent, bilateral coupling with fMRI fluctuations in a restricted set of cortical areas inside and outside the face patch network, partially matching the pattern of known anatomical projections. Within each face patch population, a subset of neurons was positively coupled with the face patch network and another was negatively coupled. The same cells showed inverse correlations with distinct subcortical structures, most notably the lateral geniculate nucleus and brainstem neuromodulatory centers. Corresponding connectivity maps derived from fMRI seeds and local field potentials differed from the single unit maps, particularly in subcortical areas. Together, the results demonstrate that the spiking fluctuations of neurons are selectively coupled with discrete brain regions, with the coupling governed in part by anatomical network connections and in part by indirect neuromodulatory pathways.

Clinical Usefulness of Simultaneous Electroencephalography and Functional Magnetic Resonance Imaging in Children With Focal Epilepsy

Background and Purpose

The current study analyzed the interictal epileptiform discharge (IED)-related hemodynamic response and aimed to determine the clinical usefulness of simultaneous electroencephalography and functional magnetic resonance imaging (EEG-fMRI) in defining the epileptogenic zone (EZ) in children with focal epilepsy.

Methods

Patients with focal epilepsy showing IEDs on conventional EEG were evaluated using EEG-fMRI. Statistical analyses were performed using the times of spike as events modeled with multiple hemodynamic response functions. The area showing the most significant t-value for blood-oxygen-level-dependent (BOLD) changes was compared with the presumed EZ. Moreover, BOLD responses between -9 and +9 s around the spike times were analyzed to track the hemodynamic response patterns over time.

Results

Half (n=13) of 26 EEG-fMRI investigations of 19 patients were successful. Two patients showed 2 different types of spikes, resulting in 15 analyses. The maximum BOLD response was concordant with the EZ in 11 (73.3%) of the 15 analyses. In 10 (66.7%) analyses, the BOLD response localized the EZs more specifically. Focal BOLD responses in the EZs occurred before IEDs in 11 analyses and were often widespread after IEDs. Hemodynamic response patterns were consistent in the same epilepsy syndrome or when repeating the investigation in the same patients.

Conclusions

EEG-fMRI can provide additional information for localizing the EZ in children with focal epilepsy, and also reveal the pathogenesis of pediatric epilepsy by evaluating the patterns in the hemodynamic response across time windows of IEDs.

How Reliable Are Ultra-Short-Term HRV Measurements during Cognitively Demanding Tasks?

Ultra-short-term HRV features assess minor autonomous nervous system variations such as variations resulting from cognitive stress peaks during demanding tasks. Several studies compare ultra-short-term and short-term HRV measurements to investigate their reliability. However, existing experiments are conducted in low cognitively demanding environments. In this paper, we propose to evaluate these measurements' reliability under cognitively demanding tasks using a near real-life setting. For this purpose, we selected 31 HRV features, extracted from data collected from 21 programmers performing code comprehension, and compared them across 18 different time frames, ranging from 3 min to 10 s. Statistical significance and correlation tests were performed between the features extracted using the larger window (3 min) and the same features extracted with the other 17 time frames. We paired these analyses with Bland-Altman plots to inspect how the extraction window size affects the HRV features. The main results show 13 features that presented at least 50% correlation when using 60-second windows. The HF and mNN features achieved around 50% correlation using a 30-second window. The 30-second window was the smallest time frame considered to have reliable measurements. Furthermore, the mNN feature proved to be quite robust to the shortening of the time resolution.

Evaluation and comparison of most prevalent artifact reduction methods for EEG acquired simultaneously with fMRI

Multimodal neuroimaging using EEG and fMRI provides deeper insights into brain function by improving the spatial and temporal resolution of the acquired data. However, simultaneous EEG-fMRI inevitably compromises the quality of the EEG and fMRI signals due to the high degree of interaction between the two systems. Fluctuations in the magnetic flux flowing through the participant and the EEG system, whether due to movement within the magnetic field of the scanner or to changes in magnetic field strength, induce electrical potentials in the EEG recordings that mask the much weaker electrical activity of the neuronal populations. A number of different methods have been proposed to reduce MR artifacts. We present an overview of the most commonly used methods and an evaluation of the methods using three sets of diverse EEG data. We limited the evaluation to open-access and easy-to-use methods and a reference signal regression method using a set of six carbon-wire loops (CWL), which allowed evaluation of their added value. The evaluation was performed by comparing EEG signals recorded outside the MRI scanner with artifact-corrected EEG signals recorded simultaneously with fMRI. To quantify and evaluate the quality of artifact reduction methods in terms of the spectral content of the signal, we analyzed changes in oscillatory activity during a resting-state and a finger tapping motor task. The quality of artifact reduction in the time domain was assessed using data collected during a visual stimulation task. In the study we utilized hierarchical Bayesian probabilistic modeling for statistical inference and observed significant differences between the evaluated methods in the success of artifact reduction and associated signal quality in both the frequency and time domains. In particular, the CWL system proved superior to the other methods evaluated in improving spectral contrast in the alpha and beta bands and in recovering visual evoked responses. Based on the results of the evaluation study, we proposed guidelines for selecting the optimal method for MR artifact reduction.

Source localization and functional network analysis in emotion cognitive reappraisal with EEG-fMRI integration

Background

The neural activity and functional networks of emotion-based cognitive reappraisal have been widely investigated using electroencephalography (EEG) and functional magnetic resonance imaging (fMRI). However, single-mode neuroimaging techniques are limited in exploring the regulation process with high temporal and spatial resolution.

Objectives

We proposed a source localization method with multimodal integration of EEG and fMRI and tested it in the source-level functional network analysis of emotion cognitive reappraisal.

Methods

EEG and fMRI data were simultaneously recorded when 15 subjects were performing the emotional cognitive reappraisal task. Fused priori weighted minimum norm estimation (FWMNE) with sliding windows was proposed to trace the dynamics of EEG source activities, and the phase lag index (PLI) was used to construct the functional brain network associated with the process of downregulating negative affect using the reappraisal strategy.

Results

The functional networks were constructed with the measure of PLI, in which the important regions were indicated. In the gamma band source-level network analysis, the cuneus, the lateral orbitofrontal cortex, the superior parietal cortex, the postcentral gyrus, and the pars opercularis were identified as important regions in reappraisal with high betweenness centrality.

Conclusion

The proposed multimodal integration method for source localization identified the key cortices involved in emotion regulation, and the network analysis demonstrated the important brain regions involved in the cognitive control of reappraisal. It shows promise in the utility in the clinical setting for affective disorders.

A-PASS: An automated pipeline to analyze simultaneously acquired EEG-fMRI data for studying brain activities during sleep

OBJECTIVE

Concurrent electroencephalography and functional magnetic resonance imaging (EEG-fMRI) signals can be used to uncover the nature of brain activities during sleep. However, analyzing simultaneously acquired EEG-fMRI data is extremely time consuming and experience dependent. Thus, we developed a pipeline, which we named A-PASS, to automatically analyze simultaneously acquired EEG-fMRI data for studying brain activities during sleep.

APPROACH

A deep learning model was trained on a sleep EEG-fMRI dataset from 45 subjects and used to perform sleep stage scoring. Various fMRI indices can be calculated with A-PASS to depict the neurophysiological characteristics across different sleep stages. We tested the performance of A-PASS on an independent sleep EEG-fMRI dataset from 28 subjects. Statistical maps regarding the main effect of sleep stages and differences between each pair of stages of fMRI indices were generated and compared using both A-PASS and manual processing methods.

MAIN RESULTS

The deep learning model implemented in A-PASS achieved both an accuracy and F1-score higher than 70% for sleep stage classification on EEG data acquired during fMRI scanning. The statistical maps generated from A-PASS largely resembled those produced from manually scored stages plus a combination of multiple software programs.

SIGNIFICANCE

A-PASS allowed efficient EEG-fMRI data processing without manual operation and could serve as a reliable and powerful tool for simultaneous EEG-fMRI studies on sleep.

Toward Completely Sampled Extracellular Neural Recording During fMRI

This work aims to estimate severe fMRI scanning artifacts in extracellular neural recordings made at ultrahigh magnetic field strengths in order to remove the artifact interferences and uncover the complete neural electrophysiology signal. We build on previous work that used PCA to denoise EEG recorded during fMRI, adapting it to cover the much larger frequency range (1-6000 Hz) of the extracellular field potentials (EFPs) observed by extracellular neural recordings. We examine the singular value decomposition (SVD)-PCA singular value shrinkage (SVS) and compare two shrinkage rules and a sliding template subtraction approach. Additionally, we present a new technique for estimating the singular value upper bounds in spontaneous neural activity recorded in the isoflurane anesthetized rat that uses the temporal first difference of the neural signal. The approaches are tested on artificial datasets to examine their efficacy in detecting extracellular action potentials (EAPs: 300-6000 Hz) recorded during fMRI gradient interferences. Our results indicate that it is possible to uncover the EAPs recorded during gradient interferences. The methods are then tested on natural (non-artificial) datasets recorded from the cortex of isoflurane anesthetized rats, where both local field potential (LFP: 1-300 Hz) and EAP signals are analyzed. The SVS methods are shown to be advantageous compared to sliding template subtraction, especially in the high frequency range corresponding to EAPs. Our novel approach moves us towards simultaneous fMRI and completely sampled neural recording (1-6000Hz with no temporal gaps), providing the opportunity for further study of spontaneous brain function and neurovascular coupling at ultrahigh field in the isoflurane anesthetized rat.

Classification-guided Neural Network-based Correction of Magnetic Resonance-related Gradient Artifact Residuals in Simultaneously Recorded Surface Electromyography

Spontaneous muscular activities can be studied by simultaneous recordings of surface electromyography (sEMG) and diffusion-weighted magnetic resonance imaging (DW-MRI). For reliable assessment of the spontaneous activity rate in sEMG data during active MR imaging, it is necessary to have a decent gradient artifact (GA) correction algorithm enabling the detection of small spontaneous activities with an amplitude of few microvolts. In this work, a neural network with weak label annotations during the training process is utilized for enhanced correction of GA residuals in the sEMG recordings. Based on sEMG signal decomposition and class-activation maps from the neural network classification, the amount of GA residuals is iteratively decreased in the sEMG signal. This leads to a reduction of the false-positive rate in automated spontaneous activity detection. Quality of GA residual correction is therefore estimated by using a specialized second neural network model. Clinical relevance- This work establishes an improved GA residual correction for simultaneously recorded sEMG data during MRI to enhance the ability for small spontaneous activity detection.

Detection of Transcranial Alternating Current Stimulation Aftereffects Is Improved by Considering the Individual Electric Field Strength and Self-Rated Sleepiness

Non-invasive electrical stimulation methods, such as transcranial alternating current stimulation (tACS), are increasingly used in human neuroscience research and offer potential new avenues to treat neurological and psychiatric disorders. However, their often variable effects have also raised concerns in the scientific and clinical communities. This study aims to investigate the influence of subject-specific factors on the alpha tACS-induced aftereffect on the alpha amplitude (measured with electroencephalography, EEG) as well as on the connectivity strength between nodes of the default mode network (DMN) [measured with functional magnetic resonance imaging (fMRI)]. As subject-specific factors we considered the individual electrical field (EFIELD) strength at target regions in the brain, the frequency mismatch between applied stimulation and individual alpha frequency (IAF) and as a covariate, subject’s changes in mental state, i.e., sleepiness. Eighteen subjects participated in a tACS and a sham session conducted on different days. Each session consisted of three runs (pre/stimulation/). tACS was applied during the second run at each subject’s individual alpha frequency (IAF), applying 1 mA peak-to-peak intensity for 7 min, using an occipital bihemispheric montage. In every run, subjects watched a video designed to increase in-scanner compliance. To investigate the aftereffect of tACS on EEG alpha amplitude and on DMN connectivity strength, EEG data were recorded simultaneously with fMRI data. Self-rated sleepiness was documented using a questionnaire. Conventional statistics (ANOVA) did not show a significant aftereffect of tACS on the alpha amplitude compared to sham stimulation. Including individual EFIELD strengths and self-rated sleepiness scores in a multiple linear regression model, significant tACS-induced aftereffects were observed. However, the subject-wise mismatch between tACS frequency and IAF had no contribution to our model. Neither standard nor extended statistical methods confirmed a tACS-induced aftereffect on DMN functional connectivity. Our results show that it is possible and necessary to disentangle alpha amplitude changes due to intrinsic mechanisms and to external manipulation using tACS on the alpha amplitude that might otherwise be overlooked. Our results suggest that EFIELD is really the most significant factor that explains the alpha amplitude modulation during a tACS session. This knowledge helps to understand the variability of the tACS-induced aftereffects.

Ballistocardiogram suppression in concurrent EEG-MRI by dynamic modeling of heartbeats

The ballistocardiogram (BCG), the induced electric potentials by the head motion originating from heartbeats, is a prominent source of noise in electroencephalography (EEG) data during magnetic resonance imaging (MRI). Although methods have been proposed to suppress the BCG artifact, more work considering the variability of cardiac cycles and head motion across time and subjects is needed to provide highly robust correction. Here, a method called "dynamic modeling of heartbeats" (DMH) is proposed to reduce BCG artifacts in EEG data recorded inside an MRI system. The DMH method models BCG artifacts by combining EEG points at time instants with similar dynamics. The modeled BCG artifact is then subtracted from the EEG recording to suppress the BCG artifact. Performance of DMH was tested and specifically compared with the Optimal Basis Set (OBS) method on EEG data recorded inside a 3T MRI system with either no MRI acquisition (Inside-MRI), echo-planar imaging (EPI-EEG), or fast MRI acquisition using simultaneous multi-slice and inverse imaging methods (SMS-InI-EEG). In a steady-state visual evoked response (SSVEP) paradigm, the 15-Hz oscillatory neuronal activity at the visual cortex after DMH processing was about 130% of that achieved by OBS processing for Inside-MRI, SMS-InI-EEG, and EPI-EEG conditions. The DMH method is computationally efficient for suppressing BCG artifacts and in the future may help to improve the quality of EEG data recorded in high-field MRI systems for neuroscientific and clinical applications.

Processing of fMRI-related anxiety and information flow between brain and body revealed a preponderance of oscillations at 0.15/0.16 Hz

Slow oscillations of different center frequencies and their coupling play an important role in brain-body interactions. The crucial question analyzed by us is, whether the low frequency (LF) band (0.05-0.15 Hz) or the intermediate frequency (IMF) band (0.1-0.2 Hz) is more eminent in respect of the information flow between body (heart rate and respiration) and BOLD signals in cortex and brainstem. A recently published study with the LF band in fMRI-naïve subjects revealed an intensive information flow from the cortex to the brainstem and a weaker flow from the brainstem to the cortex. The comparison of both bands revealed a significant information flow from the middle frontal gyrus (MFG) to the precentral gyrus (PCG) and from brainstem to PCG only in the IMF band. This pattern of directed coupling between slow oscillations in the cortex and brainstem not only supports the existence of a pacemaker-like structure in brainstem, but provides first evidence that oscillations centered at 0.15/0.16 Hz can also emerge in brain networks. BOLD oscillations in resting states are dominating at ~ 0.08 Hz and respiratory rates at ~ 0.32 Hz. Therefore, the frequency component at ~ 0.16 Hz (doubling-halving 0.08 Hz or 0.32 Hz) is of special interest, because phase coupled oscillations can reduce the energy demand.

Spontaneous alpha-band amplitude predicts subjective visibility but not discrimination accuracy during high-level perception

Near-threshold perception is a paradigm case of awareness diverging from reality - the perception of an unchanging stimulus can vacillate from undetected to clearly perceived. The amplitude of low-frequency brain oscillations - particularly in the alpha-band (8-13 Hz) - has emerged as a reliable predictor of trial-to-trial variability in perceptual decisions based on simple, low-level stimuli. Here, we addressed the question of how spontaneous oscillatory amplitude impacts subjective and objective aspects of perception using high-level visual stimuli. Human observers completed a near-threshold face/house discrimination task with subjective visibility ratings while electroencephalograms (EEG) were recorded. Using single-trial multiple regression analysis, we found that spontaneous fluctuations in prestimulus alpha-band amplitude were negatively related to visibility judgments but did not predict trial-by-trial accuracy. These results extend previous findings that indicate that strong prestimulus alpha diminishes subjective perception without affecting the accuracy or sensitivity (d') of perceptual decisions into the domain of high-level perception.