EEG synchronization measures predict epilepsy-related BOLD-fMRI fluctuations better than commonly used univariate metrics

Extraction of brain function pattern with visual-capture-task fMRI using dynamic time-window method in ADHD children

Attention deficit/Hyperactivity disorder (ADHD) has a great impact on children's development. This paper uses a novel adaptive brain state extraction algorithm to construct a dynamic time-window brain network, which captures the brain function pattern characteristics of ADHD children with higher temporal resolution. The test data were acquired by functional magnetic resonance imaging (fMRI) obtained from 23 children with ADHD during the visual-capture-task [age: (8.27 ± 2.77)]. A spatial standard deviation method is used after the initial data processing, to extract the brain activity pattern state; An improved clustering algorithm is constructed to verify the changes made to the dynamic time-window brain network model. There can be seen clear differences between each state within 0.05 s after the test. The results show that our improved new framework can effectively obtain the characteristics of dynamic brain functional connection strength changes during the task. In addition, the new algorithm is able to capture the dynamic changes of the brain network, with an 80 % improvement compared to traditional methods for the average modularity value Q. This work demonstrates a novel approach to find out the pattern changes between dynamic brain function connections, which can be of great significance for the adjuvant treatment of children with ADHD.

On temporal scale-free non-periodic stimulation and its mechanisms as an infinite improbability drive of the brain's functional connectogram

Rationalized development of electrical stimulation (ES) therapy is of paramount importance. Not only it will foster new techniques and technologies with increased levels of safety, efficacy, and efficiency, but it will also facilitate the translation from basic research to clinical practice. For such endeavor, design of new technologies must dialogue with state-of-the-art neuroscientific knowledge. By its turn, neuroscience is transitioning-a movement started a couple of decades earlier-into adopting a new conceptual framework for brain architecture, in which time and thus temporal patterns plays a central role in the neuronal representation of sampled data from the world. This article discusses how neuroscience has evolved to understand the importance of brain rhythms in the overall functional architecture of the nervous system and, consequently, that neuromodulation research should embrace this new conceptual framework. Based on such support, we revisit the literature on standard (fixed-frequency pulsatile stimuli) and mostly non-standard patterns of ES to put forward our own rationale on how temporally complex stimulation schemes may impact neuromodulation strategies. We then proceed to present a low frequency, on average (thus low energy), scale-free temporally randomized ES pattern for the treatment of experimental epilepsy, devised by our group and termed NPS (Non-periodic Stimulation). The approach has been shown to have robust anticonvulsant effects in different animal models of acute and chronic seizures (displaying dysfunctional hyperexcitable tissue), while also preserving neural function. In our understanding, accumulated mechanistic evidence suggests such a beneficial mechanism of action may be due to the natural-like characteristic of a scale-free temporal pattern that may robustly compete with aberrant epileptiform activity for the recruitment of neural circuits. Delivering temporally patterned or random stimuli within specific phases of the underlying oscillations (i.e., those involved in the communication within and across brain regions) could both potentiate and disrupt the formation of neuronal assemblies with random probability. The usage of infinite improbability drive here is obviously a reference to the "The Hitchhiker's Guide to the Galaxy" comedy science fiction classic, written by Douglas Adams. The parallel is that dynamically driving brain functional connectogram, through neuromodulation, in a manner that would not favor any specific neuronal assembly and/or circuit, could re-stabilize a system that is transitioning to fall under the control of a single attractor. We conclude by discussing future avenues of investigation and their potentially disruptive impact on neurotechnology, with a particular interest in NPS implications in neural plasticity, motor rehabilitation, and its potential for clinical translation.

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.

Synchronization in Hindmarsh-Rose neurons subject to higher-order interactions

Higher-order interactions might play a significant role in the collective dynamics of the brain. With this motivation, we here consider a simplicial complex of neurons, in particular, studying the effects of pairwise and three-body interactions on the emergence of synchronization. We assume pairwise interactions to be mediated through electrical synapses, while for second-order interactions, we separately study diffusive coupling and nonlinear chemical coupling. For all the considered cases, we derive the necessary conditions for synchronization by means of linear stability analysis, and we compute the synchronization errors numerically. Our research shows that the second-order interactions, even if of weak strength, can lead to synchronization under significantly lower first-order coupling strengths. Moreover, the overall synchronization cost is reduced due to the introduction of three-body interactions if compared to pairwise interactions.

Localizing Epileptic Foci Using Simultaneous EEG-fMRI Recording: Template Component Cross-Correlation

Conventional EEG-fMRI methods have been proven to be of limited use in the sense that they cannot reveal the information existing in between the spikes. To resolve this issue, the current study obtains the epileptic components time series detected on EEG and uses them to fit the Generalized Linear Model (GLM), as a substitution for classical regressors. This approach allows for a more precise localization, and equally importantly, the prediction of the future behavior of the epileptic generators. The proposed method approaches the localization process in the component domain, rather than the electrode domain (EEG), and localizes the generators through investigating the spatial correlation between the candidate components and the spike template, as well as the medical records of the patient. To evaluate the contribution of EEG-fMRI and concordance between fMRI and EEG, this method was applied on the data of 30 patients with refractory epilepsy. The results demonstrated the significant numbers of 29 and 24 for concordance and contribution, respectively, which mark improvement as compared to the existing literature. This study also shows that while conventional methods often fail to properly localize the epileptogenic zones in deep brain structures, the proposed method can be of particular use. For further evaluation, the concordance level between IED-related BOLD clusters and Seizure Onset Zone (SOZ) has been quantitatively investigated by measuring the distance between IED/SOZ locations and the BOLD clusters in all patients. The results showed the superiority of the proposed method in delineating the spike-generating network compared to conventional EEG-fMRI approaches. In all, the proposed method goes beyond the conventional methods by breaking the dependency on spikes and using the outside-the-scanner spike templates and the selected components, achieving an accuracy of 97%. Doing so, this method contributes to improving the yield of EEG-fMRI and creates a more realistic perception of the neural behavior of epileptic generators which is almost without precedent in the literature.

Mapping Epileptic Networks Using Simultaneous Intracranial EEG-fMRI

Background: Potentially curative epilepsy surgery can be offered if a single, discrete epileptogenic zone (EZ) can be identified. For individuals in whom there is no clear concordance between clinical localization, scalp EEG, and imaging data, intracranial EEG (icEEG) may be needed to confirm a predefined hypothesis regarding irritative zone (IZ), seizure onset zone (SOZ), and EZ prior to surgery. However, icEEG has limited spatial sampling and may fail to reveal the full extent of epileptogenic network if predefined hypothesis is not correct. Simultaneous icEEG-fMRI has been safely acquired in humans and allows exploration of neuronal activity at the whole-brain level related to interictal epileptiform discharges (IED) captured intracranially.

Methods: We report icEEG-fMRI in eight patients with refractory focal epilepsy who had resective surgery and good postsurgical outcome. Surgical resection volume in seizure-free patients post-surgically reflects confirmed identification of the EZ. IEDs on icEEG were classified according to their topographic distribution and localization (Focal, Regional, Widespread, and Non-contiguous). We also divided IEDs by their location within the surgical resection volume [primary IZ (IZ1) IED] or outside [secondary IZ (IZ2) IED]. The distribution of fMRI blood oxygen level-dependent (BOLD) changes associated with individual IED classes were assessed over the whole brain using a general linear model. The concordance of resulting BOLD map was evaluated by comparing localization of BOLD clusters with surgical resection volume. Additionally, we compared the concordance of BOLD maps and presence of BOLD clusters in remote brain areas: precuneus, cuneus, cingulate, medial frontal, and thalamus for different IED classes.

Results: A total of 38 different topographic IED classes were identified across the 8 patients: Focal (22) and non-focal (16, Regional = 9, Widespread = 2, Non-contiguous = 5). Twenty-nine IEDs originated from IZ1 and 9 from IZ2. All IED classes were associated with BOLD changes. BOLD maps were concordant with the surgical resection volume for 27/38 (71%) IED classes, showing statistical global maximum BOLD cluster or another cluster in the surgical resection volume. The concordance of BOLD maps with surgical resection volume was greater (p < 0.05) for non-focal (87.5%, 14/16) as compared to Focal (59%, 13/22) IED classes. Additionally, BOLD clusters in remote cortical and deep brain areas were present in 84% (32/38) of BOLD maps, more commonly (15/16; 93%) for non-focal IED-related BOLD maps.

Conclusions: Simultaneous icEEG-fMRI can reveal BOLD changes at the whole-brain level for a wide range of IEDs on icEEG. BOLD clusters within surgical resection volume and remote brain areas were more commonly seen for non-focal IED classes, suggesting that a wider hemodynamic network is at play.

Nonperiodic stimulation for the treatment of refractory epilepsy: Applications, mechanisms, and novel insights

Electrical stimulation of the central nervous system is a promising alternative for the treatment of pharmacoresistant epilepsy. Successful clinical and experimental stimulation is most usually carried out as continuous trains of current or voltage pulses fired at rates of 100 Hz or above, since lower frequencies yield controversial results. On the other hand, stimulation frequency should be as low as possible, in order to maximize implant safety and battery efficiency. Moreover, the development of stimulation approaches has been largely empirical in general, while they should be engineered with the neurobiology of epilepsy in mind if a more robust, efficient, efficacious, and safe application is intended. In an attempt to reconcile evidence of therapeutic effect with the understanding of the underpinnings of epilepsy, our group has developed a nonstandard form of low-frequency stimulation with randomized interpulse intervals termed nonperiodic stimulation (NPS). The rationale was that an irregular temporal pattern would impair neural hypersynchronization, which is a hallmark of epilepsy. In this review, we start by briefly revisiting the literature on the molecular, cellular, and network level mechanisms of epileptic phenomena in order to highlight this often-overlooked emergent property of cardinal importance in the pathophysiology of the disease. We then review our own studies on the efficacy of NPS against acute and chronic experimental seizures and also on the anatomical and physiological mechanism of the method, paying special attention to the hypothesis that the lack of temporal regularity induces desynchronization. We also put forward a novel insight regarding the temporal structure of NPS that may better encompass the set of findings published by the group: the fact that intervals between stimulation pulses have a distribution that follows a power law and thus may induce natural-like activity that would compete with epileptiform discharge for the recruitment of networks. We end our discussion by mentioning ongoing research and future projects of our lab.

Localization of Epileptic Foci Based on Simultaneous EEG-fMRI Data

Combining functional magnetic resonance imaging (fMRI) and electroencephalography (EEG) enables a non-invasive investigation of the human brain function and evaluation of the correlation of these two important modalities of brain activity. This paper explores recent reports on using advanced simultaneous EEG–fMRI methods proposed to map the regions and networks involved in focal epileptic seizure generation. One of the applications of EEG and fMRI combination as a valuable clinical approach is the pre-surgical evaluation of patients with epilepsy to map and localize the precise brain regions associated with epileptiform activity. In the process of conventional analysis using EEG–fMRI data, the interictal epileptiform discharges (IEDs) are visually extracted from the EEG data to be convolved as binary events with a predefined hemodynamic response function (HRF) to provide a model of epileptiform BOLD activity and use as a regressor for general linear model (GLM) analysis of the fMRI data. This review examines the methodologies involved in performing such studies, including techniques used for the recording of EEG inside the scanner, artifact removal, and statistical analysis of the fMRI signal. It then discusses the results reported for patients with primary generalized epilepsy and patients with different types of focal epileptic disorders. An important matter that these results have brought to light is that the brain regions affected by interictal epileptic discharges might not be limited to the ones where they have been generated. The developed methods can help reveal the regions involved in or affected by a seizure onset zone (SOZ). As confirmed by the reviewed literature, EEG–fMRI provides information that comes particularly useful when evaluating patients with refractory epilepsy for surgery.

Artifact Reduction in Simultaneous EEG-fMRI: A Systematic Review of Methods and Contemporary Usage

Simultaneous electroencephalography-functional MRI (EEG-fMRI) is a technique that combines temporal (largely from EEG) and spatial (largely from fMRI) indicators of brain dynamics. It is useful for understanding neuronal activity during many different event types, including spontaneous epileptic discharges, the activity of sleep stages, and activity evoked by external stimuli and decision-making tasks. However, EEG recorded during fMRI is subject to imaging, pulse, environment and motion artifact, causing noise many times greater than the neuronal signals of interest. Therefore, artifact removal methods are essential to ensure that artifacts are accurately removed, and EEG of interest is retained. This paper presents a systematic review of methods for artifact reduction in simultaneous EEG-fMRI from literature published since 1998, and an additional systematic review of EEG-fMRI studies published since 2016. The aim of the first review is to distill the literature into clear guidelines for use of simultaneous EEG-fMRI artifact reduction methods, and the aim of the second review is to determine the prevalence of artifact reduction method use in contemporary studies. We find that there are many published artifact reduction techniques available, including hardware, model based, and data-driven methods, but there are few studies published that adequately compare these methods. In contrast, recent EEG-fMRI studies show overwhelming use of just one or two artifact reduction methods based on literature published 15–20 years ago, with newer methods rarely gaining use outside the group that developed them. Surprisingly, almost 15% of EEG-fMRI studies published since 2016 fail to adequately describe the methods of artifact reduction utilized. We recommend minimum standards for reporting artifact reduction techniques in simultaneous EEG-fMRI studies and suggest that more needs to be done to make new artifact reduction techniques more accessible for the researchers and clinicians using simultaneous EEG-fMRI.

Augmenting interictal mapping with neurovascular coupling biomarkers by structured factorization of epileptic EEG and fMRI data

EEG-correlated fMRI analysis is widely used to detect regional BOLD fluctuations that are synchronized to interictal epileptic discharges, which can provide evidence for localizing the ictal onset zone. However, the typical, asymmetrical and mass-univariate approach cannot capture the inherent, higher order structure in the EEG data, nor multivariate relations in the fMRI data, and it is nontrivial to accurately handle varying neurovascular coupling over patients and brain regions. We aim to overcome these drawbacks in a data-driven manner by means of a novel structured matrix-tensor factorization: the single-subject EEG data (represented as a third-order spectrogram tensor) and fMRI data (represented as a spatiotemporal BOLD signal matrix) are jointly decomposed into a superposition of several sources, characterized by space-time-frequency profiles. In the shared temporal mode, Toeplitz-structured factors account for a spatially specific, neurovascular 'bridge' between the EEG and fMRI temporal fluctuations, capturing the hemodynamic response's variability over brain regions. By analyzing interictal data from twelve patients, we show that the extracted source signatures provide a sensitive localization of the ictal onset zone (10/12). Moreover, complementary parts of the IOZ can be uncovered by inspecting those regions with the most deviant neurovascular coupling, as quantified by two entropy-like metrics of the hemodynamic response function waveforms (9/12). Hence, this multivariate, multimodal factorization provides two useful sets of EEG-fMRI biomarkers, which can assist the presurgical evaluation of epilepsy. We make all code required to perform the computations available at https://github.com/svaneynd/structured-cmtf.

Ballistocardiogram Artifact Reduction in Simultaneous EEG-fMRI Using Deep Learning

OBJECTIVE

The concurrent recording of electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) is a technique that has received much attention due to its potential for combined high temporal and spatial resolution. However, the ballistocardiogram (BCG), a large-amplitude artifact caused by cardiac induced movement contaminates the EEG during EEG-fMRI recordings. Removal of BCG in software has generally made use of linear decompositions of the corrupted EEG. This is not ideal as the BCG signal propagates in a manner which is non-linearly dependent on the electrocardiogram (ECG). In this paper, we present a novel method for BCG artifact suppression using recurrent neural networks (RNNs).

METHODS

EEG signals were recovered by training RNNs on the nonlinear mappings between ECG and the BCG corrupted EEG. We evaluated our model's performance against the commonly used Optimal Basis Set (OBS) method at the level of individual subjects, and investigated generalization across subjects.

RESULTS

We show that our algorithm can generate larger average power reduction of the BCG at critical frequencies, while simultaneously improving task relevant EEG based classification.

CONCLUSION

The presented deep learning architecture can be used to reduce BCG related artifacts in EEG-fMRI recordings.

SIGNIFICANCE

We present a deep learning approach that can be used to suppress the BCG artifact in EEG-fMRI without the use of additional hardware. This method may have scope to be combined with current hardware methods, operate in real-time and be used for direct modeling of the BCG.

Modulation of epileptic networks by transient interictal epileptic activity: A dynamic approach to simultaneous EEG-fMRI

•

EEG-fMRI has been instrumental in characterizing brain networks in epilepsy.

•

Its value is documented in the pre-surgical assessment of drug-resistant epilepsy.

•

The delineation of brain areas to resect is fundamental for the post-surgical outcome.

•

Standard EEG-fMRI in epilepsy assesses static functional connectivity of the network.

•

EEG-fMRI dynamic connectivity identifies transitory features of specific connections.

•

We integrate dynamic fMRI connectivity and dynamic patterns of simultaneous scalp EEG.

•

This allows to better characterize the spatiotemporal aspects of epileptic networks.

•

This may help in more efficiently target the surgical intervention.

Epileptic networks, defined as brain regions involved in epileptic brain activity, have been mapped by functional connectivity in simultaneous electroencephalography and functional magnetic resonance imaging (EEG-fMRI) recordings. This technique allows to define brain hemodynamic changes, measured by the Blood Oxygen Level Dependent (BOLD) signal, associated to the interictal epileptic discharges (IED), which together with ictal events constitute a signature of epileptic disease. Given the highly time-varying nature of epileptic activity, a dynamic functional connectivity (dFC) analysis of EEG-fMRI data appears particularly suitable, having the potential to identify transitory features of specific connections in epileptic networks. In the present study, we propose a novel method, defined dFC-EEG, that integrates dFC assessed by fMRI with the information recorded by simultaneous scalp EEG, in order to identify the connections characterised by a dynamic profile correlated with the occurrence of IED, forming the dynamic epileptic subnetwork. Ten patients with drug-resistant focal epilepsy were included, with different aetiology and showing a widespread (or multilobar) BOLD activation, defined as involving at least two distinct clusters, located in two different lobes and/or extended to the hemisphere contralateral to the epileptic focus. The epileptic focus was defined from the IED-related BOLD map. Regions involved in the occurrence of interictal epileptic activity; i.e., forming the epileptic network, were identified by a general linear model considering the timecourse of the fMRI-defined focus as main regressor. dFC between these regions was assessed with a sliding-window approach. dFC timecourses were then correlated with the sliding-window variance of the IED signal (VarIED), to identify connections whose dynamics related to the epileptic activity; i.e., the dynamic epileptic subnetwork. As expected, given the very different clinical picture of each individual, the extent of this subnetwork was highly variable across patients, but was but was reduced of at least 30% with respect to the initially identified epileptic network in 9/10 patients. The connections of the dynamic subnetwork were most commonly close to the epileptic focus, as reflected by the laterality index of the subnetwork connections, reported higher than the one within the original epileptic network. Moreover, the correlation between dFC timecourses and VarIED was predominantly positive, suggesting a strengthening of the dynamic subnetwork associated to the occurrence of IED. The integration of dFC and scalp IED offers a more specific description of the epileptic network, identifying connections strongly influenced by IED. These findings could be relevant in the pre-surgical evaluation for the resection or disconnection of the epileptogenic zone and help in reaching a better post-surgical outcome. This would be particularly important for patients characterised by a widespread pathological brain activity which challenges the surgical intervention.

Cellular-Scale Microelectrode Arrays to Monitor Movement-Related Neuron Activities in the Epileptic Hippocampus of Awake Mice

OBJECTIVE

Epilepsy affects 50 million people worldwide and its pathogenesis is still unknown. In particular, the movement-related neural activities involving glutamate (Glu) and electrophysiological signals at cellular level remains unclear.

METHODS

A cellular-scale implantable microelectrode array (MEA) was fabricated to detect the movement-related neural activities involving Glu concentration and electrophysiological signals. Platinum and reduced graphene oxide nanocomposites were deposited to enhance the surface area. Glu oxidase (Gluox) were coated to effectively recognize Glu molecule.

RESULTS

Neural activities in the hippocampus of normal and epileptic mice is different, and the changes are closely connected with movement. Glu concentration and spike firing rate in the epileptic mice were much higher than those in the normal ones. And the neural activities with significant synchronization were detected in the epileptic mice even without seizure occurrence. Meanwhile, the spikes fire more intensively and Glu level became much higher during the movement of the mice compared to the stationary state.

CONCLUSION

The existing abnormality of neural activities in the epileptic mice are potential factors to induce a seizure. Movement may impact the neural activities and the duration of seizure.

SIGNIFICANCE

The MEA can monitor changes of movement, Glu and neuron discharges synchronously and provides us an effective technology to understand the neuronal disease.

Comparison of Visual and Auditory Modalities for Upper-Alpha EEG-Neurofeedback

Electroencephalography (EEG) neurofeedback (NF) training has been shown to produce long-lasting effects on the improvement of cognitive function as well as the normalization of aberrant brain activity in disease. However, the impact of the sensory modality used as the NF reinforcement signal on training effectiveness has not been systematically investigated. In this work, an EEG-based NF-training system was developed targeting the individual upper-alpha (UA) band and using either a visual or an auditory reinforcement signal, so as to compare the effects of the two sensory modalities. Sixteen healthy volunteers were randomly assigned to the Visual or Auditory group, where a radius-varying sphere or a volume-varying sound, respectively, reflected the relative amplitude of UA measured at EEG electrode Cz. Each participant underwent a total of four NF sessions, of approximately 40 min each, on consecutive days. Both groups showed significant increases in UA at Cz within sessions, and also across sessions. Effects subsequent to NF training were also found beyond the target frequency UA and scalp location Cz, namely in the lower-alpha and theta bands and in posterior brain regions, respectively. Only small differences were found on the EEG between the Visual and Auditory groups, suggesting that auditory reinforcement signals may be as effective as the more commonly used visual signals. The use of auditory NF may potentiate training protocols conducted under mobile conditions, which are now possible due to the increasing availability of wireless EEG systems.

Semi-automated EEG Enhancement Improves Localization of Ictal Onset Zone With EEG-Correlated fMRI

Objective: To improve the accuracy of detecting the ictal onset zone, we propose to enhance the epilepsy-related activity present in the EEG signals, before mapping their BOLD correlates through EEG-correlated fMRI analysis. Methods: Based solely on a segmentation of interictal epileptic discharges (IEDs) on the EEG, we train multi-channel Wiener filters (MWF) which enhance IED-like waveforms, and suppress background activity and noisy influences. Subsequently, we use EEG-correlated fMRI to find the brain regions in which the BOLD signal fluctuation corresponds to the filtered signals' time-varying power (after convolving with the hemodynamic response function), and validate the identified regions by quantitatively comparing them to ground-truth maps of the (resected or hypothesized) ictal onset zone. We validate the performance of this novel predictor vs. that of commonly used unitary or power-weighted predictors and a recently introduced connectivity-based metric, on a cohort of 12 patients with refractory epilepsy. Results: The novel predictor, derived from the filtered EEG signals, allowed the detection of the ictal onset zone in a larger percentage of epileptic patients (92% vs. at most 83% for the other predictors), and with higher statistical significance, compared to existing predictors. At the same time, the new method maintains maximal specificity by not producing false positive activations in healthy controls. Significance: The findings of this study advocate for the use of the MWF to maximize the signal-to-noise ratio of IED-like events in the interictal EEG, and subsequently use time-varying power as a sensitive predictor of the BOLD signal, to localize the ictal onset zone.

BOLD-fMRI activity informed by network variation of scalp EEG in juvenile myoclonic epilepsy

Epilepsy is marked by hypersynchronous bursts of neuronal activity, and seizures can propagate variably to any and all areas, leading to brain network dynamic organization. However, the relationship between the network characteristics of scalp EEG and blood oxygenation level-dependent (BOLD) responses in epilepsy patients is still not well known. In this study, simultaneous EEG and fMRI data were acquired in 18 juvenile myoclonic epilepsy (JME) patients. Then, the adapted directed transfer function (ADTF) values between EEG electrodes were calculated to define the time-varying network. The variation of network information flow within sliding windows was used as a temporal regressor in fMRI analysis to predict the BOLD response. To investigate the EEG-dependent functional coupling among the responding regions, modulatory interactions were analyzed for network variation of scalp EEG and BOLD time courses. The results showed that BOLD activations associated with high network variation were mainly located in the thalamus, cerebellum, precuneus, inferior temporal lobe and sensorimotor-related areas, including the middle cingulate cortex (MCC), supplemental motor area (SMA), and paracentral lobule. BOLD deactivations associated with medium network variation were found in the frontal, parietal, and occipital areas. In addition, modulatory interaction analysis demonstrated predominantly directional negative modulation effects among the thalamus, cerebellum, frontal and sensorimotor-related areas. This study described a novel method to link BOLD response with simultaneous functional network organization of scalp EEG. These findings suggested the validity of predicting epileptic activity using functional connectivity variation between electrodes. The functional coupling among the thalamus, frontal regions, cerebellum and sensorimotor-related regions may be characteristically involved in epilepsy generation and propagation, which provides new insight into the pathophysiological mechanisms and intervene targets for JME.

Widespread interictal epileptic discharge more likely than focal discharges to unveil the seizure onset zone in EEG-fMRI

OBJECTIVE

We hypothesized that the number of interictal epileptic discharges (IEDs) during scan and their spatial extent are contributing factors in obtaining appropriate activations that reveal the seizure onset zone (SOZ) in EEG-fMRI.

METHODS

157 IED types, each corresponding to one EEG scalp distribution, in 64 consecutive EEG-fMRI studies from 64 patients with refractory localization-related epilepsy were reviewed. To determine reliable activation, we used the threshold corresponding to corrected whole-brain topological false discovery rate (FDR). The location with maximum activation was compared to the presumed SOZ as defined by a comprehensive evaluation for each patient.

RESULTS

The number of IEDs was significantly higher in the types with t-value above FDR than with t-value below FDR. The presumed SOZ could be delineated in 30 of the 64 patients. Among these patients, the types of IED concordant with the SOZ had significantly larger extent on scalp EEG than the IED types discordant with the SOZ.

CONCLUSIONS

The number of IEDs is important factor in obtaining reliable activations in EEG-fMRI. IEDs with larger spatial extent are more likely to reveal, on maximum BOLD, accurate location of the SOZ.

SIGNIFICANCE

Widespread discharges are more likely to yield a reliable activation for SOZ in EEG-fMRI.

Identification of epileptic brain states by dynamic functional connectivity analysis of simultaneous EEG-fMRI: a dictionary learning approach

Most fMRI studies of the brain's intrinsic functional connectivity (FC) have assumed that this is static; however, it is now clear that it changes over time. This is particularly relevant in epilepsy, which is characterized by a continuous interchange between epileptic and normal brain states associated with the occurrence of epileptic activity. Interestingly, recurrent states of dynamic FC (dFC) have been found in fMRI data using unsupervised learning techniques, assuming either their sparse or non-sparse combination. Here, we propose an l₁-norm regularized dictionary learning (l₁-DL) approach for dFC state estimation, which allows an intermediate and flexible degree of sparsity in time, and demonstrate its application in the identification of epilepsy-related dFC states using simultaneous EEG-fMRI data. With this l₁-DL approach, we aim to accommodate a potentially varying degree of sparsity upon the interchange between epileptic and non-epileptic dFC states. The simultaneous recording of the EEG is used to extract time courses representative of epileptic activity, which are incorporated into the fMRI dFC state analysis to inform the selection of epilepsy-related dFC states. We found that the proposed l₁-DL method performed best at identifying epilepsy-related dFC states, when compared with two alternative methods of extreme sparsity (k-means clustering, maximum; and principal component analysis, minimum), as well as an l₀-norm regularization framework (l₀-DL), with a fixed amount of temporal sparsity. We further showed that epilepsy-related dFC states provide novel insights into the dynamics of epileptic networks, which go beyond the information provided by more conventional EEG-correlated fMRI analysis, and which were concordant with the clinical profile of each patient. In addition to its application in epilepsy, our study provides a new dFC state identification method of potential relevance for studying brain functional connectivity dynamics in general.

Adaptive optimal basis set for BCG artifact removal in simultaneous EEG-fMRI

Electroencephalography (EEG) signals recorded during simultaneous functional magnetic resonance imaging (fMRI) are contaminated by strong artifacts. Among these, the ballistocardiographic (BCG) artifact is the most challenging, due to its complex spatio-temporal dynamics associated with ongoing cardiac activity. The presence of BCG residuals in EEG data may hide true, or generate spurious correlations between EEG and fMRI time-courses. Here, we propose an adaptive Optimal Basis Set (aOBS) method for BCG artifact removal. Our method is adaptive, as it can estimate the delay between cardiac activity and BCG occurrence on a beat-to-beat basis. The effective creation of an optimal basis set by principal component analysis (PCA) is therefore ensured by a more accurate alignment of BCG occurrences. Furthermore, aOBS can automatically estimate which components produced by PCA are likely to be BCG artifact-related and therefore need to be removed. The aOBS performance was evaluated on high-density EEG data acquired with simultaneous fMRI in healthy subjects during visual stimulation. As aOBS enables effective reduction of BCG residuals while preserving brain signals, we suggest it may find wide application in simultaneous EEG-fMRI studies.

EEG-Informed fMRI: A Review of Data Analysis Methods

The simultaneous acquisition of electroencephalography (EEG) with functional magnetic resonance imaging (fMRI) is a very promising non-invasive technique for the study of human brain function. Despite continuous improvements, it remains a challenging technique, and a standard methodology for data analysis is yet to be established. Here we review the methodologies that are currently available to address the challenges at each step of the data analysis pipeline. We start by surveying methods for pre-processing both EEG and fMRI data. On the EEG side, we focus on the correction for several MR-induced artifacts, particularly the gradient and pulse artifacts, as well as other sources of EEG artifacts. On the fMRI side, we consider image artifacts induced by the presence of EEG hardware inside the MR scanner, and the contamination of the fMRI signal by physiological noise of non-neuronal origin, including a review of several approaches to model and remove it. We then provide an overview of the approaches specifically employed for the integration of EEG and fMRI when using EEG to predict the blood oxygenation level dependent (BOLD) fMRI signal, the so-called EEG-informed fMRI integration strategy, the most commonly used strategy in EEG-fMRI research. Finally, we systematically review methods used for the extraction of EEG features reflecting neuronal phenomena of interest.