Fluctuations of the EEG-fMRI correlation reflect intrinsic strength of functional connectivity in default mode network

Multi-scale analysis of acupuncture mechanisms for motor and sensory cortex activity based on SEEG data

Acupuncture, a traditional Chinese therapy, is gaining attention for its impact on the brain. While existing electroencephalogram and functional magnetic resonance image research has made significant contributions, this paper utilizes stereo-electroencephalography data for a comprehensive exploration of neurophysiological effects. Employing a multi-scale approach, channel-level analysis reveals notable $\delta $-band activity changes during acupuncture. At the brain region level, acupuncture modulated connectivity between the paracentral lobule and the precentral gyrus. Whole-brain analysis indicates acupuncture's influence on network organization, and enhancing $E_{glob}$ and increased interaction between the motor and sensory cortex. Brain functional reorganization is an important basis for functional recovery or compensation after central nervous system injury. The use of acupuncture to stimulate peripheral nerve trunks, muscle motor points, acupoints, etc., in clinical practice may contribute to the reorganization of brain function. This multi-scale perspective provides diverse insights into acupuncture's effects. Remarkably, this paper pioneers the introduction of stereo-electroencephalography data, advancing our understanding of acupuncture's mechanisms and potential therapeutic benefits in clinical settings.

Neurovascular coupling in eye-open-eye-close task and resting state: Spectral correspondence between concurrent EEG and fMRI

Neurovascular coupling serves as an essential neurophysiological mechanism in functional neuroimaging, which is generally presumed to be robust and invariant across different physiological states, encompassing both task engagement and resting state. Nevertheless, emerging evidence suggests that neurovascular coupling may exhibit state dependency, even in normal human participants. To investigate this premise, we analyzed the cross-frequency spectral correspondence between concurrently recorded electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) data, utilizing them as proxies for neurovascular coupling during the two conditions: an eye-open-eye-close (EOEC) task and a resting state. We hypothesized that given the state dependency of neurovascular coupling, EEG-fMRI spectral correspondences would change between the two conditions in the visual system. During the EOEC task, we observed a negative phase-amplitude-coupling (PAC) between EEG alpha-band and fMRI visual activity. Conversely, in the resting state, a pronounced amplitude-amplitude-coupling (AAC) emerged between EEG and fMRI signals, as evidenced by the spectral correspondence between the EEG gamma-band of the midline occipital channel (Oz) and the high-frequency fMRI signals (0.15-0.25 Hz) in the visual network. This study reveals distinct scenarios of EEG-fMRI spectral correspondence in healthy participants, corroborating the state-dependent nature of neurovascular coupling.

State-transition dynamics of resting-state functional magnetic resonance imaging data: model comparison and test-to-retest analysis

Electroencephalogram (EEG) microstate analysis entails finding dynamics of quasi-stable and generally recurrent discrete states in multichannel EEG time series data and relating properties of the estimated state-transition dynamics to observables such as cognition and behavior. While microstate analysis has been widely employed to analyze EEG data, its use remains less prevalent in functional magnetic resonance imaging (fMRI) data, largely due to the slower timescale of such data. In the present study, we extend various data clustering methods used in EEG microstate analysis to resting-state fMRI data from healthy humans to extract their state-transition dynamics. We show that the quality of clustering is on par with that for various microstate analyses of EEG data. We then develop a method for examining test-retest reliability of the discrete-state transition dynamics between fMRI sessions and show that the within-participant test-retest reliability is higher than between-participant test-retest reliability for different indices of state-transition dynamics, different networks, and different data sets. This result suggests that state-transition dynamics analysis of fMRI data could discriminate between different individuals and is a promising tool for performing fingerprinting analysis of individuals.

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.

Infra-slow fluctuations in cortical potentials and respiration drive fast cortical EEG rhythms in sleeping and waking states

OBJECTIVE

Infra-slow fluctuations (ISF, 0.008-0.1 Hz) characterize hemodynamic and electric potential signals of human brain. ISFs correlate with the amplitude dynamics of fast (>1 Hz) neuronal oscillations, and may arise from permeability fluctuations of the blood-brain barrier (BBB). It is unclear if physiological rhythms like respiration drive or track fast cortical oscillations, and the role of sleep in this coupling is unknown.

METHODS

We used high-density full-band electroencephalography (EEG) in healthy human volunteers (N = 21) to measure concurrently the ISFs, respiratory pulsations, and fast neuronal oscillations during periods of wakefulness and sleep, and to assess the strength and direction of their phase-amplitude coupling.

RESULTS

The phases of ISFs and respiration were both coupled with the amplitude of fast neuronal oscillations, with stronger ISF coupling being evident during sleep. Phases of ISF and respiration drove the amplitude dynamics of fast oscillations in sleeping and waking states, with different contributions.

CONCLUSIONS

ISFs in slow cortical potentials and respiration together significantly determine the dynamics of fast cortical oscillations.

SIGNIFICANCE

We propose that these slow physiological phases play a significant role in coordinating cortical excitability, which is a fundamental aspect of brain function.

Aging affects the phase coherence between spontaneous oscillations in brain oxygenation and neural activity

The risk of neurodegenerative disorders increases with age, due to reduced vascular nutrition and impaired neural function. However, the interactions between cardiovascular dynamics and neural activity, and how these interactions evolve in healthy aging, are not well understood. Here, the interactions are studied by assessment of the phase coherence between spontaneous oscillations in cerebral oxygenation measured by fNIRS, the electrical activity of the brain measured by EEG, and cardiovascular functions extracted from ECG and respiration effort, all simultaneously recorded. Signals measured at rest in 21 younger participants (31.1 ± 6.9 years) and 24 older participants (64.9 ± 6.9 years) were analysed by wavelet transform, wavelet phase coherence and ridge extraction for frequencies between 0.007 and 4 Hz. Coherence between the neural and oxygenation oscillations at ∼ 0.1 Hz is significantly reduced in the older adults in 46/176 fNIRS-EEG probe combinations. This reduction in coherence cannot be accounted for in terms of reduced power, thus indicating that neurovascular interactions change with age. The approach presented promises a noninvasive means of evaluating the efficiency of the neurovascular unit in aging and disease.

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.

Transcranial Direct Current Stimulation Modulates EEG Microstates in Low-Functioning Autism: A Pilot Study

Autism spectrum disorder (ASD) is a heterogeneous disorder that affects several behavioral domains of neurodevelopment. Transcranial direct current stimulation (tDCS) is a new method that modulates motor and cognitive function and may have potential applications in ASD treatment. To identify its potential effects on ASD, differences in electroencephalogram (EEG) microstates were compared between children with typical development (n = 26) and those with ASD (n = 26). Furthermore, children with ASD were divided into a tDCS (experimental) and sham stimulation (control) group, and EEG microstates and Autism Behavior Checklist (ABC) scores before and after tDCS were compared. Microstates A, B, and D differed significantly between children with TD and those with ASD. In the experimental group, the scores of microstates A and C and ABC before tDCS differed from those after tDCS. Conversely, in the control group, neither the EEG microstates nor the ABC scores before the treatment period (sham stimulation) differed from those after the treatment period. This study indicates that tDCS may become a viable treatment for ASD.

Multimodal resting-state connectivity predicts affective neurofeedback performance

Neurofeedback has been suggested as a potential complementary therapy to different psychiatric disorders. Of interest for this approach is the prediction of individual performance and outcomes. In this study, we applied functional connectivity-based modeling using electroencephalography (EEG) and functional near-infrared spectroscopy (fNIRS) modalities to (i) investigate whether resting-state connectivity predicts performance during an affective neurofeedback task and (ii) evaluate the extent to which predictive connectivity profiles are correlated across EEG and fNIRS techniques. The fNIRS oxyhemoglobin and deoxyhemoglobin concentrations and the EEG beta and gamma bands modulated by the alpha frequency band (beta-m-alpha and gamma-m-alpha, respectively) recorded over the frontal cortex of healthy subjects were used to estimate functional connectivity from each neuroimaging modality. For each connectivity matrix, relevant edges were selected in a leave-one-subject-out procedure, summed into "connectivity summary scores" (CSS), and submitted as inputs to a support vector regressor (SVR). Then, the performance of the left-out-subject was predicted using the trained SVR model. Linear relationships between the CSS across both modalities were evaluated using Pearson's correlation. The predictive model showed a mean absolute error smaller than 20%, and the fNIRS oxyhemoglobin CSS was significantly correlated with the EEG gamma-m-alpha CSS (r = -0.456, p = 0.030). These results support that pre-task electrophysiological and hemodynamic resting-state connectivity are potential predictors of neurofeedback performance and are meaningfully coupled. This investigation motivates the use of joint EEG-fNIRS connectivity as outcome predictors, as well as a tool for functional connectivity coupling investigation.

Simultaneous electroencephalography-functional magnetic resonance imaging for assessment of human brain function

Electroencephalography (EEG) and functional Magnetic Resonance Imaging (MRI) have long been used as tools to examine brain activity. Since both methods are very sensitive to changes of synaptic activity, simultaneous recording of EEG and fMRI can provide both high temporal and spatial resolution. Therefore, the two modalities are now integrated into a hybrid tool, EEG-fMRI, which encapsulates the useful properties of the two. Among other benefits, EEG-fMRI can contribute to a better understanding of brain connectivity and networks. This review lays its focus on the methodologies applied in performing EEG-fMRI studies, namely techniques used for the recording of EEG inside the scanner, artifact removal, and statistical analysis of the fMRI signal. We will investigate simultaneous resting-state and task-based EEG-fMRI studies and discuss their clinical and technological perspectives. Moreover, it is established that the brain regions affected by a task-based neural activity might not be limited to the regions in which they have been initiated. Advanced methods can help reveal the regions responsible for or affected by a developed neural network. Therefore, we have also looked into studies related to characterization of structure and dynamics of brain networks. The reviewed literature suggests that EEG-fMRI can provide valuable complementary information about brain neural networks and functions.

Source localized infraslow neurofeedback training in people with chronic painful knee osteoarthritis: A randomized, double-blind, sham-controlled feasibility clinical trial

Persistent pain is a key symptom in people living with knee osteoarthritis (KOA). Infra-slow Neurofeedback (ISF-NF) training is a recent development focusing on modulating cortical slow-wave activity to improve pain outcomes. A parallel, two-armed double-blinded, randomized sham-controlled, feasibility clinical trial aimed to determine the feasibility and safety of a novel electroencephalography-based infraslow fluctuation neurofeedback (EEG ISF-NF) training in people with KOA and determine the variability of clinical outcomes and EEG changes following NF training. Eligible participants attended nine 30-min ISF-NF training sessions involving three cortical regions linked to pain. Feasibility measures were monitored during the trial period. Pain and functional outcomes were measured at baseline, post-intervention, and follow-up after 2 weeks. Resting-state EEG was recorded at baseline and immediate post-intervention. Participants were middle-aged (61.7 ± 7.6 years), New Zealand European (90.5%), and mostly females (62%) with an average knee pain duration of 4 ± 3.4 years. The study achieved a retention rate of 91%, with 20/22 participants completing all the sessions. Participants rated high levels of acceptance and “moderate to high levels of perceived effectiveness of the training.” No serious adverse events were reported during the trial. Mean difference (95% CI) for clinical pain and function measures are as follows for pain severity [active: 0.89 ± 1.7 (−0.27 to 2.0); sham: 0.98 ± 1.1 (0.22–1.7)], pain interference [active: 0.75 ± 2.3 (−0.82 to 2.3); Sham: 0.89 ± 2.1 (−0.60 to 2.4)], pain unpleasantness [active: 2.6 ± 3.7 (0.17–5.1); sham: 2.8 ± 3 (0.62–5.0)] and physical function [active: 6.2 ± 13 (−2.6 to 15); sham: 1.6 ± 12 (−6.8 to 10)]. EEG sources demonstrated frequency-specific neuronal activity, functional connectivity, and ISF ratio changes following NF training. The findings of the study indicated that the ISF-NF training is a feasible, safe, and acceptable intervention for pain management in people with KOA, with high levels of perceived effectiveness. The study also reports the variability in clinical, brain activity, and connectivity changes following training.

Blind Visualization of Task-Related Networks From Visual Oddball Simultaneous EEG-fMRI Data: Spectral or Spatiospectral Model?

Various disease conditions can alter EEG event-related responses and fMRI-BOLD signals. We hypothesized that event-related responses and their clinical alterations are imprinted in the EEG spectral domain as event-related (spatio)spectral patterns (ERSPat). We tested four EEG-fMRI fusion models utilizing EEG power spectra fluctuations (i.e., absolute spectral model - ASM; relative spectral model - RSM; absolute spatiospectral model - ASSM; and relative spatiospectral model - RSSM) for fully automated and blind visualization of task-related neural networks. Two (spatio)spectral patterns (high δ 4 band and low β 1 band) demonstrated significant negative linear relationship (p FWE < 0.05) to the frequent stimulus and three patterns (two low δ 2 and δ 3 bands, and narrow θ 1 band) demonstrated significant positive relationship (p < 0.05) to the target stimulus. These patterns were identified as ERSPats. EEG-fMRI F-map of each δ 4 model showed strong engagement of insula, cuneus, precuneus, basal ganglia, sensory-motor, motor and dorsal part of fronto-parietal control (FPCN) networks with fast HRF peak and noticeable trough. ASM and RSSM emphasized spatial statistics, and the relative power amplified the relationship to the frequent stimulus. For the δ 4 model, we detected a reduced HRF peak amplitude and a magnified HRF trough amplitude in the frontal part of the FPCN, default mode network (DMN) and in the frontal white matter. The frequent-related β 1 patterns visualized less significant and distinct suprathreshold spatial associations. Each θ 1 model showed strong involvement of lateralized left-sided sensory-motor and motor networks with simultaneous basal ganglia co-activations and reduced HRF peak and amplified HRF trough in the frontal part of the FPCN and DMN. The ASM θ 1 model preserved target-related EEG-fMRI associations in the dorsal part of the FPCN. For δ 4, β 1, and θ 1 bands, all models provided high local F-statistics in expected regions. The most robust EEG-fMRI associations were observed for ASM and RSSM.

Topography-Related EEG-fMRI in Surgically Confirmed Epileptic Foci: A Comparison to Spike-Related EEG-fMRI in Clinical Practice

EEG-fMRI has gained increasing importance in epilepsy pre-surgical diagnosis. However, 40-70% of EEG-fMRI recordings in patients lack interictal epileptiform discharges (IEDs) during the scan, which could be overcome by detecting matching topography maps. We tried to validate this method in clinical settings taking various electroclinical factors into consideration. Eleven patients who had undergone EEG-fMRI during pre-surgical evaluation for drug-resistant epilepsy and who had had clinical long-term video-EEG were studied. Spike-related blood oxygen level-dependent (BOLD) maps were created using IEDs occurring during the EEG-fMRI scan. Separate maps were then generated from IEDs marked on the clinical long-term EEG recordings, which were averaged to produce topographical IED maps and correlated with the EEGs recorded inside the scanner yielding a correlation coefficient time course. Epileptogenic zones were defined by an expert panel during pre-surgical evaluation and validated by an epilepsy surgery resulting in a good outcome. Both techniques' performance was evaluated according to factors including arousal during IED recording, IED topography and lateralization, lesion type, and localization. Topography-related EEG-fMRI yielded more specific results compared to the spike-related method. Superficial lesion location and ipsilateral IED seem to result in a higher concordance of BOLD maps. The polarity of BOLD responses may be lesion-dependent, and both positive and negative BOLD changes may be associated with the irritative zone. Topography-related EEG-fMRI may show improved specificity especially for superficial lesions producing ipsilateral spikes. This method can be used as an alternative either in the absence of spikes during the simultaneous EEG-fMRI acquisition or to sharpen a diffusely activated BOLD-map.

Resting-state EEG Dynamics Reveals Differences in Network Organization and its Fluctuation between Frequency Bands

Functional connectivity in EEG resting-state is not stable but fluctuates considerably. The aim of this study was to investigate how efficient information flows through a network, i.e. how resting-state EEG networks are organized and whether this organization it also subject to fluctuations. Differences of the network organization (small-worldness), degree of clustered connectivity, and path length as an indicator of how information is integrated into the network across time was compared between theta, alpha and beta bands. We show robust differences in network organization (small-worldness) between frequency bands. Fluctuations in network organization were larger in the theta, compared to the alpha and beta frequency. Variation in network organization and not the frequency of fluctuations differs between frequency bands. Furthermore, the degree of clustered connectivity and its modulation across time is the same across frequency bands, but the path length revealed the same modulatory pattern as the small-world metric. It is therefore the interplay of local processing efficiency and global information processing efficiency in the brain that fluctuates in a frequency-specific way. Properties of how information can be integrated is subject to fluctuations in a frequency-specific way in the resting-state. The possible relevance of these resting-state EEG properties is discussed including its clinical relevance.

15 Years MR-encephalography

Objective

This review article gives an account of the development of the MR-encephalography (MREG) method, which started as a mere ‘Gedankenexperiment’ in 2005 and gradually developed into a method for ultrafast measurement of physiological activities in the brain. After going through different approaches covering k-space with radial, rosette, and concentric shell trajectories we have settled on a stack-of-spiral trajectory, which allows full brain coverage with (nominal) 3 mm isotropic resolution in 100 ms. The very high acceleration factor is facilitated by the near-isotropic k-space coverage, which allows high acceleration in all three spatial dimensions.

Methods

The methodological section covers the basic sequence design as well as recent advances in image reconstruction including the targeted reconstruction, which allows real-time feedback applications, and—most recently—the time-domain principal component reconstruction (tPCR), which applies a principal component analysis of the acquired time domain data as a sparsifying transformation to improve reconstruction speed as well as quality.

Applications

Although the BOLD-response is rather slow, the high speed acquisition of MREG allows separation of BOLD-effects from cardiac and breathing related pulsatility. The increased sensitivity enables direct detection of the dynamic variability of resting state networks as well as localization of single interictal events in epilepsy patients. A separate and highly intriguing application is aimed at the investigation of the glymphatic system by assessment of the spatiotemporal patterns of cardiac and breathing related pulsatility.

Discussion

MREG has been developed to push the speed limits of fMRI. Compared to multiband-EPI this allows considerably faster acquisition at the cost of reduced image quality and spatial resolution.

Modeling regional changes in dynamic stability during sleep and wakefulness

Global brain states are frequently placed within a unidimensional continuum by correlational studies, ranging from states of deep unconsciousness to ordinary wakefulness. An alternative is their multidimensional and mechanistic characterization in terms of different cognitive capacities, using computational models to reproduce the underlying neural dynamics. We explore this alternative by introducing a semi-empirical model linking regional activation and long-range functional connectivity in the different brain states visited during the natural wake-sleep cycle. Our model combines functional magnetic resonance imaging (fMRI) data, in vivo estimates of structural connectivity, and anatomically-informed priors to constrain the independent variation of regional activation. The best fit to empirical data was achieved using priors based on functionally coherent networks, with the resulting model parameters dividing the cortex into regions presenting opposite dynamical behavior. Frontoparietal regions approached a bifurcation from dynamics at a fixed point governed by noise, while sensorimotor regions approached a bifurcation from oscillatory dynamics. In agreement with human electrophysiological experiments, sleep onset induced subcortical deactivation with low correlation, which was subsequently reversed for deeper stages. Finally, we introduced periodic forcing of variable intensity to simulate external perturbations, and identified the key regions relevant for the recovery of wakefulness from deep sleep. Our model represents sleep as a state with diminished perceptual gating and the latent capacity for global accessibility that is required for rapid arousals. To the extent that the qualitative characterization of local dynamics is exhausted by the dichotomy between unstable and stable behavior, our work highlights how expanding the model parameter space can describe states of consciousness in terms of multiple dimensions with interpretations given by the choice of anatomically-informed priors.

Intrinsic connectome organization across temporal scales: New insights from cross-modal approaches

The discovery of a stable, whole-brain functional connectivity organization that is largely independent of external events has drastically extended our view of human brain function. However, this discovery has been primarily based on functional magnetic resonance imaging (fMRI). The role of this whole-brain organization in fast oscillation-based connectivity as measured, for example, by electroencephalography (EEG) and magnetoencephalography (MEG) is only beginning to emerge. Here, we review studies of intrinsic connectivity and its whole-brain organization in EEG, MEG, and intracranial electrophysiology with a particular focus on direct comparisons to connectome studies in fMRI. Synthesizing this literature, we conclude that irrespective of temporal scale over four orders of magnitude, intrinsic neurophysiological connectivity shows spatial similarity to the connectivity organization commonly observed in fMRI. A shared structural connectivity basis and cross-frequency coupling are possible mechanisms contributing to this similarity. Acknowledging that a stable whole-brain organization governs long-range coupling across all timescales of neural processing motivates researchers to take "baseline" intrinsic connectivity into account when investigating brain-behavior associations, and further encourages more widespread exploration of functional connectomics approaches beyond fMRI by using EEG and MEG modalities.

Simultaneous EEG-fMRI for Functional Neurological Assessment

The increasing incidence of neurodegenerative and psychiatric diseases requires increasingly sophisticated tools for their diagnosis and monitoring. Clinical assessment takes advantage of objective parameters extracted by electroencephalogram and magnetic resonance imaging (MRI) among others, to support clinical management of neurological diseases. The complementarity of these two tools can be now emphasized by the possibility of integrating the two technologies in a hybrid solution, allowing simultaneous acquisition of the two signals by the novel EEG-fMRI technology. This review will focus on simultaneous EEG-fMRI technology and related early studies, dealing about issues related to the acquisition and processing of simultaneous signals, and including critical discussion about clinical and technological perspectives.

EEG spatiospectral patterns and their link to fMRI BOLD signal via variable hemodynamic response functions

BACKGROUND

Spatial and temporal resolution of brain network activity can be improved by combining different modalities. Functional Magnetic Resonance Imaging (fMRI) provides full brain coverage with limited temporal resolution, while electroencephalography (EEG), estimates cortical activity with high temporal resolution. Combining them may provide improved network characterization.

NEW METHOD

We examined relationships between EEG spatiospectral pattern timecourses and concurrent fMRI BOLD signals using canonical hemodynamic response function (HRF) with its 1^st and 2^nd temporal derivatives in voxel-wise general linear models (GLM). HRF shapes were derived from EEG-fMRI time courses during "resting-state", visual oddball and semantic decision paradigms.

RESULTS

The resulting GLM F-maps self-organized into several different large-scale brain networks (LSBNs) often with different timing between EEG and fMRI revealed through differences in GLM-derived HRF shapes (e.g., with a lower time to peak than the canonical HRF). We demonstrate that some EEG spatiospectral patterns (related to concurrent fMRI) are weakly task-modulated.

COMPARISON WITH EXISTING METHOD(S)

Previously, we demonstrated 14 independent EEG spatiospectral patterns within this EEG dataset, stable across the resting-state, visual oddball and semantic decision paradigms. Here, we demonstrate that their time courses are significantly correlated with fMRI dynamics organized into LSBN structures. EEG-fMRI derived HRF peak appears earlier than the canonical HRF peak, which suggests limitations when assuming a canonical HRF shape in EEG-fMRI.

CONCLUSIONS

This is the first study examining EEG-fMRI relationships among independent EEG spatiospectral patterns over different paradigms. The findings highlight the importance of considering different HRF shapes when spatiotemporally characterizing brain networks using EEG and fMRI.