Epileptic networks studied with EEG-fMRI

Mesoscopic mapping of hemodynamic responses and neuronal activity during pharmacologically induced interictal spikes in awake and anesthetized mice

Imaging hemodynamic responses to interictal spikes holds promise for presurgical epilepsy evaluations. Understanding the hemodynamic response function is crucial for accurate interpretation. Prior interictal neurovascular coupling data primarily come from anesthetized animals, impacting reliability. We simultaneously monitored calcium fluctuations in excitatory neurons, hemodynamics, and local field potentials (LFP) during bicuculline-induced interictal events in both isoflurane-anesthetized and awake mice. Isoflurane significantly affected LFP amplitude but had little impact on the amplitude and area of the calcium signal. Anesthesia also dramatically blunted the amplitude and latency of the hemodynamic response, although not its area of spread. Cerebral blood volume change provided the best spatial estimation of excitatory neuronal activity in both states. Targeted silencing of the thalamus in awake mice failed to recapitulate the impact of anesthesia on hemodynamic responses suggesting that isoflurane's interruption of the thalamocortical loop did not contribute either to the dissociation between the LFP and the calcium signal nor to the alterations in interictal neurovascular coupling. The blood volume increase associated with interictal spikes represents a promising mapping signal in both the awake and anesthetized states.

Electroencephalography-functional magnetic resonance imaging for clinical evaluation in focal epilepsy

OBJECTIVE

We aimed to evaluate the contribution of simultaneous recording of electroencephalography-functional magnetic resonance imaging (EEG-fMRI) in the diagnosis of epilepsy syndrome, localization of the epileptogenic zone (EZ), and decision-making regarding surgical treatment.

METHODS

We performed a retrospective study to evaluate patients with focal epilepsy who underwent EEG-fMRI. Two evaluators assessed epilepsy syndrome, presumed focus, and surgical candidacy and defined confidence levels. They assessed these clinical characteristics first without EEG-fMRI and then including EEG-fMRI to assess how the results of EEG-fMRI changed the evaluations. We also determined how the clinical evaluation was affected by the concordance level between the blood oxygen level-dependent (BOLD) response and the presumed focus location, and by the confidence level of the BOLD response itself based on the t-value of the primary and secondary clusters.

RESULTS

Fifty-one scans from 48 patients were included. The BOLD map affected 66.7% of the evaluations by altering evaluation items (epilepsy syndrome, presumed focus, or surgical candidacy) or their confidence levels. EEG-fMRI results increased the confidence levels of epilepsy syndrome, presumed focus, or surgical candidacy in 47.1% of patients but reduced clinical confidence in these features in 11.8%. More specifically, the confidence levels increased for epilepsy syndrome in 28.5%, identification of presumed focus in 33.9%, and determination of surgical candidacy in 29.4%. The BOLD signal confidence level, whether high or low, did not influence these clinical factors.

SIGNIFICANCE

Previous studies have emphasized the utility of EEG-fMRI for the localization of the epileptogenic zone. This study demonstrated the potential of EEG-fMRI to influence clinical confidence when determining epilepsy syndrome, the presumed epileptic focus, and surgical candidacy.

Altered spread of waves of activities at large scale is influenced by cortical thickness organization in temporal lobe epilepsy: a magnetic resonance imaging-high-density electroencephalography study

Temporal lobe epilepsy is a brain network disorder characterized by alterations at both the structural and the functional levels. It remains unclear how structure and function are related and whether this has any clinical relevance. In the present work, we adopted a novel methodological approach investigating how network structural features influence the large-scale dynamics. The functional network was defined by the spatio-temporal spreading of aperiodic bursts of activations (neuronal avalanches), as observed utilizing high-density electroencephalography in patients with temporal lobe epilepsy. The structural network was modelled as the region-based thickness covariance. Loosely speaking, we quantified the similarity of the cortical thickness of any two brain regions, both across groups and at the individual level, the latter utilizing a novel approach to define the subject-wise structural covariance network. In order to compare the structural and functional networks (at the nodal level), we studied the correlation between the probability that a wave of activity would propagate from a source to a target region and the similarity of the source region thickness as compared with other target brain regions. Building on the recent evidence that large-waves of activities pathologically spread through the epileptogenic network in temporal lobe epilepsy, also during resting state, we hypothesize that the structural cortical organization might influence such altered spatio-temporal dynamics. We observed a stable cluster of structure-function correlation in the bilateral limbic areas across subjects, highlighting group-specific features for left, right and bilateral temporal epilepsy. The involvement of contralateral areas was observed in unilateral temporal lobe epilepsy. We showed that in temporal lobe epilepsy, alterations of structural and functional networks pair in the regions where seizures propagate and are linked to disease severity. In this study, we leveraged on a well-defined model of neurological disease and pushed forward personalization approaches potentially useful in clinical practice. Finally, the methods developed here could be exploited to investigate the relationship between structure-function networks at subject level in other neurological conditions.

Language lateralization mapping (reversibly) masked by non-dominant focal epilepsy: a case report

Language lateralization in patients with focal epilepsy frequently diverges from the left-lateralized pattern that prevails in healthy right-handed people, but the mechanistic explanations are still a matter of debate. Here, we debate the complex interaction between focal epilepsy, language lateralization, and functional neuroimaging techniques by introducing the case of a right-handed patient with unaware focal seizures preceded by aphasia, in whom video-EEG and PET examination suggested the presence of focal cortical dysplasia in the right superior temporal gyrus, despite a normal structural MRI. The functional MRI for language was inconclusive, and the neuropsychological evaluation showed mild deficits in language functions. A bilateral stereo-EEG was proposed confirming the right superior temporal gyrus origin of seizures, revealing how ictal aphasia emerged only once seizures propagated to the left superior temporal gyrus and confirming, by cortical mapping, the left lateralization of the posterior language region. Stereo-EEG-guided radiofrequency thermocoagulations of the (right) focal cortical dysplasia not only reduced seizure frequency but led to the normalization of the neuropsychological assessment and the "restoring" of a classical left-lateralized functional MRI pattern of language. This representative case demonstrates that epileptiform activity in the superior temporal gyrus can interfere with the functioning of the contralateral homologous cortex and its associated network. In the case of presurgical evaluation in patients with epilepsy, this interference effect must be carefully taken into consideration. The multimodal language lateralization assessment reported for this patient further suggests the sensitivity of different explorations to this interference effect. Finally, the neuropsychological and functional MRI changes after thermocoagulations provide unique cues on the network pathophysiology of focal cortical dysplasia and the role of diverse techniques in indexing language lateralization in complex scenarios.

Wearable OPM-MEG: A changing landscape for epilepsy

Magnetoencephalography with optically pumped magnometers (OPM-MEG) is an emerging and novel, cost-effective wearable system that can simultaneously record neuronal activity with high temporal resolution ("when" neuronal activity occurs) and spatial resolution ("where" neuronal activity occurs). This paper will first outline recent methodological advances in OPM-MEG compared to conventional superconducting quantum interference device (SQUID)-MEG before discussing how OPM-MEG can become a valuable and noninvasive clinical support tool in epilepsy surgery evaluation. Although OPM-MEG and SQUID-MEG share similar data features, OPM-MEG is a wearable design that fits children and adults, and it is also robust to head motion within a magnetically shielded room. This means that OPM-MEG can potentially extend the application of MEG into the neurobiology of severe childhood epilepsies with intellectual disabilities (e.g., epileptic encephalopathies) without sedation. It is worth noting that most OPM-MEG sensors are heated, which may become an issue with large OPM sensor arrays (OPM-MEG currently has fewer sensors than SQUID-MEG). Future implementation of triaxial sensors may alleviate the need for large OPM sensor arrays. OPM-MEG designs allowing both awake and sleep recording are essential for potential long-term epilepsy monitoring.

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.

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.

Involvement of brain structures in childhood epilepsy with centrotemporal spikes

BACKGROUND

We aimed to investigate electroencephalography (EEG)-functional magnetic resonance imaging (fMRI) findings to elucidate the interictal epileptiform discharge (IED)-related functional alterations in deep brain structures and the neocortex in childhood epilepsy with centrotemporal spikes (CECTS).

METHODS

Ten children with CECTS (median age 8.2 years), referred to our hospital within a year of onset, were eligible for inclusion. They underwent EEG-fMRI recording during sleep. Llongitudinal evaluations, including medical examinations, intelligence tests, and questionnaires about developmental disabilities, were performed. The initial evaluation was performed at the same time as the EEG-fMRI, and the second evaluation was performed over 2 years after the initial evaluation.

RESULTS

Three children were unable to maintain sleep during the EEG-fMRI recording, and the remaining seven children were eligible for further assessment. All patients showed unilateral-dominant centrotemporal spikes during scans. One patient had only positive hemodynamic responses, while the others had both positive and negative hemodynamic responses. All patients showed IED-related hemodynamic responses in the bilateral neocortex. For deep brain structures, IED-related hemodynamic responses were observed in the cingulate gyrus (n = 4), basal ganglia (n = 3), thalamus (n = 2), and default mode network (n = 1). Seizure frequencies at the second evaluation were infrequent or absent, and the longitudinal results of intelligence tests and questionnaires were within normal ranges.

CONCLUSIONS

We demonstrated that IEDs affect broad brain areas, including deep brain structures such as the cingulate gyrus, basal ganglia, and thalamus. Deep brain structures may play an important role in the pathophysiology of CECTS.

New interinstitutional, multimodal presurgical evaluation protocol associated with improved seizure freedom for poorly defined cases of focal epilepsy in children

OBJECTIVE

In an attempt to improve postsurgical seizure outcomes for poorly defined cases (PDCs) of pediatric focal epilepsy (i.e., those that are not visible or well defined on 3T MRI), the authors modified their presurgical evaluation strategy. Instead of relying on concordance between video-electroencephalography and 3T MRI and using functional imaging and intracranial recording in select cases, the authors systematically used a multimodal, 3-tiered investigation protocol that also involved new collaborations between their hospital, the Montreal Children's Hospital, and the Montreal Neurological Institute. In this study, the authors examined how their new strategy has impacted postsurgical outcomes. They hypothesized that it would improve postsurgical seizure outcomes, with the added benefit of identifying a subset of tests contributing the most.

METHODS

Chart review was performed for children with PDCs who underwent resection following the new strategy (i.e., new protocol [NP]), and for the same number who underwent treatment previously (i.e., preprotocol [PP]); ≥ 1-year follow-up was required for inclusion. Well-defined, multifocal, and diffuse hemispheric cases were excluded. Preoperative demographics and clinical characteristics, resection volumes, and pathology, as well as seizure outcomes (Engel class Ia vs > Ia) at 1 year postsurgery and last follow-up were reviewed.

RESULTS

Twenty-two consecutive NP patients were compared with 22 PP patients. There was no difference between the two groups for resection volumes, pathology, or preoperative characteristics, except that the NP group underwent more presurgical evaluation tests (p < 0.001). At 1 year postsurgery, 20 of 22 NP patients and 10 of 22 PP patients were seizure free (OR 11.81, 95% CI 2.00-69.68; p = 0.006). Magnetoencephalography and PET/MRI were associated with improved postsurgical seizure outcomes, but both were highly correlated with the protocol group (i.e., independent test effects could not be demonstrated).

CONCLUSIONS

A new presurgical evaluation strategy for children with PDCs of focal epilepsy led to improved postsurgical seizure freedom. No individual presurgical evaluation test was independently associated with improved outcome, suggesting that it may be the combined systematic protocol and new interinstitutional collaborations that makes the difference rather than any individual test.

Temporal Lobe Spikes Affect Distant Intrinsic Connectivity Networks

Objective: To evaluate local and distant blood oxygen level dependent (BOLD) signal changes related to interictal epileptiform discharges (IED) in drug-resistant temporal lobe epilepsy (TLE).

Methods: Thirty-three TLE patients undergoing EEG–functional Magnetic Resonance Imaging (fMRI) as part of the presurgical workup were consecutively enrolled. First, a single-subject spike-related analysis was performed: (a) to verify the BOLD concordance with the presumed Epileptogenic Zone (EZ); and (b) to investigate the Intrinsic Connectivity Networks (ICN) involvement. Then, a group analysis was performed to search for common BOLD changes in TLE.

Results: Interictal epileptiform discharges were recorded in 25 patients and in 19 (58%), a BOLD response was obtained at the single-subject level. In 42% of the cases, BOLD changes were observed in the temporal lobe, although only one patient had a pure concordant finding, with a single fMRI cluster overlapping (and limited to) the EZ identified by anatomo-electro-clinical correlations. In the remaining 58% of the cases, BOLD responses were localized outside the temporal lobe and the presumed EZ. In every patient, with a spike-related fMRI map, at least one ICN appeared to be involved. Four main ICNs were preferentially involved, namely, motor, visual, auditory/motor speech, and the default mode network. At the single-subject level, EEG–fMRI proved to have high specificity (above 65%) in detecting engagement of an ICN and the corresponding ictal/postictal symptom, and good positive predictive value (above 67%) in all networks except the visual one. Finally, in the group analysis of BOLD changes related to IED revealed common activations at the right precentral gyrus, supplementary motor area, and middle cingulate gyrus.

Significance: Interictal temporal spikes affect several distant extra-temporal areas, and specifically the motor/premotor cortex. EEG–fMRI in patients with TLE eligible for surgery is recommended not for strictly localizing purposes rather it might be useful to investigate ICNs alterations at the single-subject level.

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.

Stationary EEG pattern relates to large-scale resting state networks - An EEG-fMRI study connecting brain networks across time-scales

Relating brain dynamics acting on time scales that differ by at least an order of magnitude is a fundamental issue in brain research. The same is true for the observation of stable dynamical structures in otherwise highly non-stationary signals. The present study addresses both problems by the analysis of simultaneous resting state EEG-fMRI recordings of 53 patients with epilepsy. Confirming previous findings, we observe a generic and temporally stable average correlation pattern in EEG recordings. We design a predictor for the General Linear Model describing fluctuations around the stationary EEG correlation pattern and detect resting state networks in fMRI data. The acquired statistical maps are contrasted to several surrogate tests and compared with maps derived by spatial Independent Component Analysis of the fMRI data. By means of the proposed EEG-predictor we observe core nodes of known fMRI resting state networks with high specificity in the default mode, the executive control and the salience network. Our results suggest that both, the stationary EEG pattern as well as resting state fMRI networks are different expressions of the same brain activity. This activity is interpreted as the dynamics on (or close to) a stable attractor in phase space that is necessary to maintain the brain in an efficient operational mode. We discuss that this interpretation is congruent with the theoretical framework of complex systems as well as with the brain's energy balance.

Functional reorganization of brain regions into a network in childhood absence epilepsy: A magnetoencephalography study

OBJECTIVE

Epilepsy is considered as a network disorder. However, it is unknown how normal brain activity develops into the highly synchronized discharging activity seen in disordered networks. This study aimed to explore the epilepsy brain network and the significant re-combined brain areas in childhood absence epilepsy (CAE).

METHODS

Twenty-two children with CAE were recruited to study the neural source activity during ictal-onset and interictal periods at frequency bands of 1-30 Hz and 30-80 Hz with magnetoencephalography (MEG) scanning. Accumulated source imaging (ASI) was used to analyze the locations of neural source activity and peak source strength.

RESULTS

Most of the participants had more active source activity locations in the ictal-onset period rather than in the interictal period, both at 1-30 Hz and 30-80 Hz. The frontal lobe (FL), the temporo-parietal junction (T-P), and the parietal lobe (PL) became the main active areas of source activity during the ictal period, while the precuneus (PC), cuneus, and thalamus were relatively inactive.

CONCLUSIONS

Some brain areas become more excited and have increased source activity during seizures. These significant brain regions might be re-combined to form an epilepsy network that regulates the process of absence seizures.

SIGNIFICANCE

The study confirmed that important brain regions are reorganized in an epilepsy network, which provides a basis for exploring the network mechanism of CAE development. Imaging findings may provide a reference for clinical characteristics.

Investigation of dynamic functional connectivity of the source reconstructed epileptiform discharges in focal epilepsy: A graph theory approach

OBJECTIVE

The aim of the present study is to investigate with noninvasive methods the modulation of dynamic functional connectivity during interictal epileptiform discharge (IED).

METHOD

We reconstructed the cortical source of the EEG recorded IED of 17 patients with focal epilepsy. We then computed dynamic connectivity using the time resolved phase locking value (PLV). We derived graph theory indices (i.e. degree, strength, local efficiency, clustering coefficient and global efficiency). Finally, we selected the atlas node with the maximum activation as the IED cortical source investigating the graph indices dynamics in theta, alpha, beta and gamma frequency bands.

RESULTS

We observed IED-locked modulations of the graph indexes depending on the frequency bands. We detected a modulation of the strength, clustering coefficient, local and global efficiency both in theta and in alpha bands, which also displayed modulations of the degree index. In the beta band only the global efficiency was modulated by the IED, while no effects were detected in the gamma band. Finally, we found a correlation between alpha and theta local efficiency, as well as alpha global efficiency, and the epilepsy duration.

SIGNIFICANCE

Our findings suggest that the neural synchronization is not limited to the IED cortical source, but implies a phase synchronization across multiple brain areas. We hypothesize that the aberrant electrical activity originating from the IED locus is spread amongst the other network nodes throughout the low frequency bands (i.e. theta and alpha). Moreover, IED-dependent increase in the global efficiency indicates that the IED interfere with the whole network functioning. We finally discussed possible application of this methodology for future investigation.

Different Modes of Low-Frequency Focused Ultrasound-Mediated Attenuation of Epilepsy Based on the Topological Theory

Epilepsy is common brain dysfunction, where abnormal synchronized activities can be observed across multiple brain regions. Low-frequency focused pulsed ultrasound has been proven to modulate the epileptic brain network. In this study, we used two modes of low-intensity focused ultrasound (pulsed-wave and continuous-wave) to sonicate the brains of KA-induced epileptic rats, analyzed the EEG functional brain connections to explore their respective effect on the epileptic brain network, and discuss the mechanism of ultrasound neuromodulation. By comparing the brain network characteristics before and after sonication, we found that two modes of ultrasound both significantly affected the functional brain network, especially in the low-frequency band below 12 Hz. After two modes of sonication, the power spectral density of the EEG signals and the connection strength of the brain network were significantly reduced, but there was no significant difference between the two modes. Our results indicated that the ultrasound neuromodulation could effectively regulate the epileptic brain connections. The ultrasound-mediated attenuation of epilepsy was independent of modes of ultrasound.

Executive functioning in children with self-limited epilepsy with centrotemporal spikes: a systematic review and meta-analysis

Self-limited epilepsy with centrotemporal spikes (ECTS) is one of the most frequent focal epilepsies amongst children. Because remission usually occurs before 16 years old and patients present infrequent clinical manifestation, ECTS was considered benign for a long time. Despite the reports on cognitive deficits associated with ECTS in the last years, knowledge about the condition's specific executive function domains (inhibitory control, working memory, cognitive flexibility, verbal fluency, and higher-order executive functions) is still lacking. The following systematic review was conducted according to PRISMA guidelines. The PubMed and Scopus databases and gray literature were searched according to the following eligibility criteria: (1) original articles published in peer-review journals; (2) studies that present assessment of children with ECTS; and (3) studies with an available assessment of the executive function of the participants. A total of 43 studies (1179 patients and 1086 healthy controls) met the inclusion criteria. Data from 19 studies were extracted, and meta-analysis methods were used to compare results in the three main executive function domains and verbal fluency. The study quality was measured through the Newcastle-Ottawa Scale (NOS) and the evidence quality with the GRADEpro tool. Results and conclusions: The present systematic review is the first to gather information about executive functioning in children with ECTS. According to the meta-analyses, children with ECTS show weaker performances when compared with a control group in inhibitory control, cognitive flexibility, and verbal fluency. However, because the quality of evidence was classified as very low, caution is needed when interpreting the strength of the results.

Safety and data quality of EEG recorded simultaneously with multi-band fMRI

PURPOSE

Simultaneously recorded electroencephalography and functional magnetic resonance imaging (EEG-fMRI) is highly informative yet technically challenging. Until recently, there has been little information about EEG data quality and safety when used with newer multi-band (MB) fMRI sequences. Here, we measure the relative heating of a MB protocol compared with a standard single-band (SB) protocol considered to be safe. We also evaluated EEG quality recorded concurrently with the MB protocol on humans.

MATERIALS AND METHODS

We compared radiofrequency (RF)-related heating at multiple electrodes and magnetic field magnitude, B1+RMS, of a MB fMRI sequence with whole-brain coverage (TR = 440 ms, MB factor = 4) against a previously recommended, safe SB sequence using a phantom outfitted with a 64-channel EEG cap. Next, 9 human subjects underwent eyes-closed resting state EEG-fMRI using the MB sequence. Additionally, in three of the subjects resting state EEG was recorded also during the SB sequence and in an fMRI-free condition to directly compare EEG data quality across scanning conditions. EEG data quality was assessed by the ability to remove gradient and cardioballistic artifacts along with a clean spectrogram.

RESULTS

The heating induced by the MB sequence was lower than that of the SB sequence by a factor of 0.73 ± 0.38. This is consistent with an expected heating ratio of 0.64, calculated from the square of the ratio of B1+RMS values of the sequences. In the resting state EEG data, gradient and cardioballistic artifacts were successfully removed using traditional template subtraction. All subjects showed an individual alpha peak in the spectrogram with a posterior topography characteristic of eyes-closed EEG. The success of artifact rejection for the MB sequence was comparable to that in traditional SB sequences.

CONCLUSIONS

Our study shows that B1+RMS is a useful indication of the relative heating of fMRI protocols. This observation indicates that simultaneous EEG-fMRI recordings using this MB sequence can be safe in terms of RF-related heating, and that EEG data recorded using this sequence is of acceptable quality after traditional artifact removal techniques.

Different Functional Network Connectivity Patterns in Epilepsy: A Rest-State fMRI Study on Mesial Temporal Lobe Epilepsy and Benign Epilepsy With Centrotemporal Spike

The stark discrepancy in the prognosis of epilepsy is closely related to brain damage features and underlying mechanisms, which have not yet been unraveled. In this study, differences in the epileptic brain functional connectivity states were explored through a network-based connectivity analysis between intractable mesial temporal lobe epilepsy (MTLE) patients and benign epilepsy with centrotemporal spikes (BECT). Resting state fMRI imaging data were collected for 14 MTLE patients, 12 BECT patients and 16 healthy controls (HCs). Independent component analysis (ICA) was performed to identify the cortical functional networks. Subcortical nuclei of interest were extracted from the Harvard-Oxford probability atlas. Network-based statistics were used to detect functional connectivity (FC) alterations across intranetworks and internetworks, including the connectivity between cortical networks and subcortical nuclei. Compared with HCs, MTLE patients showed significant lower activity between the connectivity of cortical networks and subcortical nuclei (especially hippocampus) and lower internetwork FC involving the lateral temporal lobe; BECT patients showed normal cortical-subcortical FC with hyperconnectivity between cortical networks. Together, cortical-subcortical hypoconnectivity in MTLE suggested a low efficiency and collaborative network pattern, and this might be relevant to the final decompensatory state and the intractable prognosis. Conversely, cortical-subcortical region with normal connectivity remained well in global cooperativity, and compensatory internetwork hyperconnectivity caused by widespread cortical abnormal discharge, which might account for the self-limited clinical outcome in BECT. Based on the fMRI functional network study, different brain network patterns might provide a better explanation of mechanisms in different types of epilepsy.

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.

Connective profiles and antagonism between dynamic and static connectivity underlying generalized epilepsy

This study aims to characterize the connective profiles and the coupling relationship between dynamic and static functional connectivity (dFC and sFC) in large-scale brain networks in patients with generalized epilepsy (GE). Functional, structural and diffuse MRI data were collected from 83 patients with GE and 106 matched healthy controls (HC). Resting-state BOLD time course was deconvolved to neural time course using a blind hemodynamic deconvolution method. Then, five connective profiles, including the structural connectivity (SC) and BOLD/neural time course-derived sFC/dFC networks, were constructed based on the proposed whole brain atlas. Network-level weighted correlation probability (NWCP) were proposed to evaluate the association between dFC and sFC. Both the BOLD signal and neural time course showed highly concordant findings and the present study emphasized the consistent findings between two functional approaches. The patients with GE showed hypervariability and enhancement of FC, and notably decreased SC in the subcortical network. Besides, increased dFC, weaker anatomic links, and complex alterations of sFC were observed in the default mode network of GE. Moreover, significantly increased SC and predominantly increased sFC were found in the frontoparietal network. Remarkably, antagonism between dFC and sFC was observed in large-scale networks in HC, while patients with GE showed significantly decreased antagonism in core epileptic networks. In sum, our study revealed distinct connective profiles in different epileptic networks and provided new clues to the brain network mechanism of epilepsy from the perspective of antagonism between dynamic and static functional connectivity.

Shared and distinct global signal topography disturbances in subcortical and cortical networks in human epilepsy

Epilepsy is a common brain network disorder associated with disrupted large-scale excitatory and inhibitory neural interactions. Recent resting-state fMRI evidence indicates that global signal (GS) fluctuations that have commonly been ignored are linked to neural activity. However, the mechanisms underlying the altered global pattern of fMRI spontaneous fluctuations in epilepsy remain unclear. Here, we quantified GS topography using beta weights obtained from a multiple regression model in a large group of epilepsy with different subtypes (98 focal temporal epilepsy; 116 generalized epilepsy) and healthy population (n = 151). We revealed that the nonuniformly distributed GS topography across association and sensory areas in healthy controls was significantly shifted in patients. Particularly, such shifts of GS topography disturbances were more widespread and bilaterally distributed in the midbrain, cerebellum, visual cortex, and medial and orbital cortex in generalized epilepsy, whereas in focal temporal epilepsy, these networks spread beyond the temporal areas but mainly remain lateralized. Moreover, we found that these abnormal GS topography patterns were likely to evolve over the course of a longer epilepsy disease. Our study demonstrates that epileptic processes can potentially affect global excitation/inhibition balance and shift the normal GS topological distribution. These progressive topographical GS disturbances in subcortical-cortical networks may underlie pathophysiological mechanisms of global fluctuations in human epilepsy.

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.

Ecstatic and gelastic seizures related to the hypothalamus

Ecstatic seizures constitute a rare form of epilepsy, and the semiology is diverse. Previously, brain areas including the temporal lobe and the insula have been identified to be involved in clinical expression. The aim of this report is to review changes in ecstatic seizures in a patient before and after operation for a hypothalamic hamartoma, and to scrutinize the relation to gelastic seizures. In this case, the ecstatic seizures disappeared after surgery of the hamartoma but reappeared eleven years later. Clinical information was retrospectively obtained from medical records, interviews, and a questionnaire covering seizure semiology that pertained to ecstatic and gelastic seizures. Our findings imply a possible connection between gelastic and ecstatic seizures, originating from a hypothalamic hamartoma. To our knowledge, this location has not previously been described in ecstatic seizures. Gelastic seizures may in this case be associated with ecstatic seizures. We speculate that patients with ecstatic seizures may have an ictal activation of neuronal networks that involve the insula. Our case may add information to the knowledge concerning ecstatic seizures.

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.

Dynamic analysis of fMRI activation during epileptic spikes can help identify the seizure origin

OBJECTIVE

We use the dynamic electroencephalography-functional magnetic resonance imaging (EEG-fMRI) method to incorporate variability in the amplitude and field of the interictal epileptic discharges (IEDs) into the fMRI analysis. We ask whether IED variability analysis can (a) identify additional activated brain regions during the course of IEDs, not seen in standard analysis; and (b) demonstrate the origin and spread of epileptic activity. We explore whether these functional changes recapitulate the structural connections and propagation of epileptic activity during seizures.

METHODS

Seventeen patients with focal epilepsy and at least 30 IEDs of a single type during simultaneous EEG-fMRI were studied. IED variability and EEG source imaging (ESI) analysis extracted time-varying dynamic changes. General linear modeling (GLM) generated static functional maps. Dynamic maps were compared to static functional maps. The dynamic sequence from IED variability was compared to the ESI results. In a subset of patients, we investigated structural connections between active brain regions using diffusion-based fiber tractography.

RESULTS

IED variability distinguished the origin of epileptic activity from its propagation in 15 of 17 (88%) patients. This included two cases where no result was obtained from the standard GLM analysis. In both of these cases, IED variability revealed activation in line with the presumed epileptic focus. Two cases showed no result from either method. Both had very high spike rates associated with dysplasia in the postcentral gyrus. In all 15 cases with dynamic activation, the observed dynamics were concordant with ESI. Fiber tractography identified specific white matter pathways between brain regions that were active at IED onset and propagation.

SIGNIFICANCE

Dynamic techniques involving IED variability can provide additional power for EEG-fMRI analysis, compared to standard analysis, revealing additional biologically plausible information in cases with no result from the standard analysis and gives insight into the origin and spread of IEDs.

Respiratory-related brain pulsations are increased in epilepsy-a two-centre functional MRI study

Resting-state functional MRI has shown potential for detecting changes in cerebral blood oxygen level-dependent signal in patients with epilepsy, even in the absence of epileptiform activity. Furthermore, it has been suggested that coefficient of variation mapping of fast functional MRI signal may provide a powerful tool for the identification of intrinsic brain pulsations in neurological diseases such as dementia, stroke and epilepsy. In this study, we used fast functional MRI sequence (magnetic resonance encephalography) to acquire ten whole-brain images per second. We used the functional MRI data to compare physiological brain pulsations between healthy controls (n = 102) and patients with epilepsy (n = 33) and furthermore to drug-naive seizure patients (n = 9). Analyses were performed by calculating coefficient of variation and spectral power in full band and filtered sub-bands. Brain pulsations in the respiratory-related frequency sub-band (0.11-0.51 Hz) were significantly (P < 0.05) increased in patients with epilepsy, with an increase in both signal variance and power. At the individual level, over 80% of medicated and drug-naive seizure patients exhibited areas of abnormal brain signal power that correlated well with the known clinical diagnosis, while none of the controls showed signs of abnormality with the same threshold. The differences were most apparent in the basal brain structures, respiratory centres of brain stem, midbrain and temporal lobes. Notably, full-band, very low frequency (0.01-0.1 Hz) and cardiovascular (0.8-1.76 Hz) brain pulses showed no differences between groups. This study extends and confirms our previous results of abnormal fast functional MRI signal variance in epilepsy patients. Only respiratory-related brain pulsations were clearly increased with no changes in either physiological cardiorespiratory rates or head motion between the subjects. The regional alterations in brain pulsations suggest that mechanisms driving the cerebrospinal fluid homeostasis may be altered in epilepsy. Magnetic resonance encephalography has both increased sensitivity and high specificity for detecting the increased brain pulsations, particularly in times when other tools for locating epileptogenic areas remain inconclusive.

Ecstatic and gelastic seizures relate to the hypothalamus

Ecstatic seizures constitute a rare form of epilepsy, and the semiology is diverse. Previously, brain areas including the temporal lobe and the insula have been identified to be involved in clinical expression. The aim of this report is to review changes in ecstatic seizures in a patient before and after operation of a hypothalamic hamartoma, and to scrutinize the relation to gelastic seizures. In this case, the ecstatic seizures disappeared after surgery of the hamartoma but reappeared eleven years later. Clinical information was retrospectively obtained from medical records, interviews, and a questionnaire covering seizure semiology that pertained to ecstatic and gelastic seizures. Our findings imply a possible connection between gelastic and ecstatic seizures, originating from a hypothalamic hamartoma. To our knowledge, this location has not previously been described in ecstatic seizures. Gelastic seizures may in this case were associated with ecstatic seizures. We speclate patients with ecstatic seizures may have an ictal activation of neuronal networks that involves the insula. Our case may add information to the growing knowledge concerning ecstatic seizures.

Subsecond EEG-fMRI analysis for presurgical evaluation in focal epilepsy

OBJECTIVE

The authors recently reported a novel subsecond analysis method of analyzing EEG-functional MRI (fMRI) to improve the detection rate of epileptic focus. This study aims to validate the utility of this method for presurgical evaluation in pharmacoresistant focal epilepsy.

METHODS

Among 13 patients with focal epilepsy undergoing presurgical examinations including simultaneous EEG-fMRI at 3T, 11 patients had interictal epileptiform discharges (IEDs) during fMRI. The authors used the sequence of topographic maps during the IEDs as a reference to obtain subsecond fMRI activation maps with the same temporal resolution as the EEG data, and constructed "spike-and-slow-wave-activation-summary" (SSWAS) maps that showed the activation frequency of voxels during IEDs. Clusters were defined by thresholding the SSWAS maps (voxel value > 10), and those containing voxels with the top 3 highest activation frequencies were considered significant. Significant hemodynamic responses using conventional event-related (ER) analysis and SSWAS maps were compared with the resection areas and surgical outcomes at 1 year after surgery.

RESULTS

Using ER analysis, 4 (36%) of 11 patients had significant hemodynamic responses. One of 4 patients had significant hemodynamic responses in the resection area and good surgical outcome. Using SSWAS maps, 10 (91%) of 11 patients had significant hemodynamic responses. Six of 10 patients had significant hemodynamic responses in the resection area, and 5 of the 6 patients achieved good surgical outcomes. The remaining 4 patients had significant hemodynamic responses distant from the resection area, and only 1 of the 4 patients achieved good surgical outcomes. The sensitivity, specificity, positive predictive value, and negative predictive value of SSWAS maps were 83.3%, 75.0%, 83.3%, and 75.0%, respectively.

CONCLUSIONS

This study demonstrated the clinical utility of SSWAS maps for presurgical evaluation of pharmacoresistant focal epilepsy. The findings indicated that subsecond EEG-fMRI analysis may help surgeons choose the resection areas that could lead to good surgical outcomes.

Beyond the Epileptic Focus: Functional Epileptic Networks in Focal Epilepsy

Focal epilepsy can be conceptualized as a network disorder, and the functional epileptic network can be described as a complex system of multiple brain areas that interact dynamically to generate epileptic activity. However, we still do not fully understand the functional architecture of epileptic networks. We studied a cohort of 21 patients with extratemporal focal epilepsy. We used independent component analysis of functional magnetic resonance imaging (fMRI) data. In order to identify the epilepsy-related components, we examined the general linear model-derived electroencephalography-fMRI (EEG–fMRI) time courses associated with interictal epileptic activity as intrinsic hemodynamic epileptic biomarkers. Independent component analysis revealed components related to the epileptic time courses in all 21 patients. Each epilepsy-related component described a network of spatially distributed brain areas that corresponded to the specific epileptic network in each patient. We also provided evidence for the interaction between the epileptic activity generated at the epileptic network and the physiological resting state networks. Our findings suggest that independent component analysis, guided by EEG–fMRI epileptic time courses, have the potential to define the functional architecture of the epileptic network in a noninvasive way. These data could be useful in planning invasive EEG electrode placement, guiding surgical resections, and more effective therapeutic interventions.

Hydrogel-Based Organic Subdural Electrode with High Conformability to Brain Surface

A totally soft organic subdural electrode has been developed by embedding an array of poly(3,4-ethylenedioxythiophene)-modified carbon fabric (PEDOT-CF) into the polyvinyl alcohol (PVA) hydrogel substrate. The mesh structure of the stretchable PEDOT-CF allowed stable structural integration with the PVA substrate. The electrode performance for monitoring electrocorticography (ECoG) was evaluated in saline solution, on ex vivo brains, and in vivo animal experiments using rats and porcines. It was demonstrated that the large double-layer capacitance of the PEDOT-CF brings low impedance at the frequency of brain wave including epileptic seizures, and PVA hydrogel substrate minimized the contact impedance on the brain. The most important unique feature of the hydrogel-based ECoG electrode was its shape conformability to enable tight adhesion even to curved, grooved surface of brains by just being placed. In addition, since the hydrogel-based electrode is totally organic, the simultaneous ECoG-fMRI measurements could be conducted without image artifacts, avoiding problems induced by conventional metallic electrodes.

Background EEG Connectivity Captures the Time-Course of Epileptogenesis in a Mouse Model of Epilepsy

Large-scale brain networks are increasingly recognized as important for the generation of seizures in epilepsy. However, how a network evolves from a healthy state through the process of epileptogenesis remains unclear. To address this question, here, we study longitudinal epicranial background EEG recordings (30 electrodes, EEG free from epileptiform activity) of a mouse model of mesial temporal lobe epilepsy. We analyze functional connectivity networks and observe that over the time course of epileptogenesis the networks become increasingly asymmetric. Furthermore, computational modelling reveals that a set of nodes, located outside of the region of initial insult, emerges as particularly important for the network dynamics. These findings are consistent with experimental observations, thus demonstrating that ictogenic mechanisms can be revealed on the EEG, that computational models can be used to monitor unfolding epileptogenesis and that both the primary focus and epileptic network play a role in epileptogenesis.

Spike-related haemodynamic responses overlap with high frequency oscillations in patients with focal epilepsy

Simultaneous scalp EEG/functional MRI measures non-invasively haemodynamic responses to interictal epileptic discharges, which are related to the epileptogenic zone. High frequency oscillations are also an excellent indicator of this zone, but are primarily recorded from intracerebral EEG. We studied the spatial overlap of these two important markers in patients with drug-resistant epilepsy to assess if their combination could help better define the extent of the epileptogenic zone. We included patients who underwent EEG-functional MRI and later intracerebral EEG. Based on intracerebral EEG findings, we separated patients with unifocal seizures from patients with multifocal or unknown onset seizures. Haemodynamic t-maps were coregistered with the intracerebral electrode positions. Each EEG channel was classified as pertaining to one of the following categories: primary haemodynamic cluster (maximum t-value), secondary cluster (t-value > 90% of the primary cluster) or outside the primary and secondary clusters. We marked high frequency oscillations (ripples: 80-250 Hz; fast ripples: 250-500 Hz) during 1 h of slow wave sleep, and compared their rates in each haemodynamic category. After classifying channels as high- or low-rate, the proportion of high-rate channels within the primary or primary plus secondary clusters was compared to the proportion expected by chance. Twenty-five patients, 11 with unifocal and 14 with multifocal/unknown seizure onsets, were studied. We found a significantly higher median high frequency oscillation rate in the primary cluster compared to secondary cluster and outside these two clusters for the unifocal group (P < 0.0001), but not for the multifocal/unknown group. For the unifocal group, the number of high-rate channels within the primary or primary plus secondary clusters was significantly higher than expected by chance. This held only for the high-ripple-rate channels in the multifocal/unknown group. At the patient level, most patients (18/25, or 72%) had at least one high-rate channel within a primary cluster. In patients with unifocal epilepsy, the maximum haemodynamic response (primary cluster) related to scalp interictal discharges overlaps with the tissue generating high frequency oscillations at high rates. If intracranial EEG is warranted, this response should be explored. As a tentative clinical use of the combination of these techniques we propose that higher high frequency oscillation rates inside than outside the maximum response indicates that the patient has indeed a focal epileptogenic zone demarcated by this response, whereas similar rates inside and outside may indicate a widespread epileptogenic zone or an epileptogenic zone not covered by the implantation.

Component-related BOLD response to localize epileptic focus using simultaneous EEG-fMRI recordings at 3T

BACKGROUND

Simultaneous EEG-fMRI experiments record spatiotemporal dynamics of epileptic activity. A shortcoming of spike-based EEG-fMRI studies is their inability to provide information about behavior of epileptic generators when no spikes are visible.

NEW METHOD

We extract time series of epileptic components identified on EEG and fit them with Generalized Linear Model (GLM) model. This allows a precise and reliable localization of epileptic foci in addition to predicting generator's behavior. The proposed method works in the source domain and delineates generators considering spatial correlation between spike template and candidate components in addition to patient's medical records.

RESULTS

The proposed method was applied on 20 patients with refractory epilepsy and 20 age- and gender-matched healthy controls. The identified components were examined statistically and threshold of localization accuracy was determined as 86% based on Receiver Operating Characteristic (ROC) curve analysis. Accuracy, sensitivity, and specificity were found to be 88%, 85%, and 95%, respectively. Contribution of EEG-fMRI and concordance between EEG and fMRI were also evaluated. Concordance was found in 19 patients and contribution in 17.

COMPARISON WITH EXISTING METHODS

We compared the proposed method with conventional methods. Our comparisons showed superiority of the proposed method. In particular, when epileptogenic zone was located deep in the brain, the method outperformed existing methods.

CONCLUSIONS

This study contributes substantially to increasing the yield of EEG-fMRI and presents a realistic estimate of the neural behavior of epileptic generators, to the best of our knowledge, for the first time in the literature.

Spectral entropy indicates electrophysiological and hemodynamic changes in drug-resistant epilepsy - A multimodal MREG study

Objective

Epilepsy causes measurable irregularity over a range of brain signal frequencies, as well as autonomic nervous system functions that modulate heart and respiratory rate variability. Imaging dynamic neuronal signals utilizing simultaneously acquired ultra-fast 10 Hz magnetic resonance encephalography (MREG), direct current electroencephalography (DC-EEG), and near-infrared spectroscopy (NIRS) can provide a more comprehensive picture of human brain function. Spectral entropy (SE) is a nonlinear method to summarize signal power irregularity over measured frequencies. SE was used as a joint measure to study whether spectral signal irregularity over a range of brain signal frequencies based on synchronous multimodal brain signals could provide new insights in the neural underpinnings of epileptiform activity.

Methods

Ten patients with focal drug-resistant epilepsy (DRE) and ten healthy controls (HC) were scanned with 10 Hz MREG sequence in combination with EEG, NIRS (measuring oxygenated, deoxygenated, and total hemoglobin: HbO, Hb, and HbT, respectively), and cardiorespiratory signals. After pre-processing, voxelwise SE_MREG was estimated from MREG data. Different neurophysiological and physiological subfrequency band signals were further estimated from MREG, DC-EEG, and NIRS: fullband (0–5 Hz, FB), near FB (0.08–5 Hz, NFB), brain pulsations in very-low (0.009–0.08 Hz, VLFP), respiratory (0.12–0.4 Hz, RFP), and cardiac (0.7–1.6 Hz, CFP) frequency bands. Global dynamic fluctuations in MREG and NIRS were analyzed in windows of 2 min with 50% overlap.

Results

Right thalamus, cingulate gyrus, inferior frontal gyrus, and frontal pole showed significantly higher SE_MREG in DRE patients compared to HC. In DRE patients, SE of cortical Hb was significantly reduced in FB (p = .045), NFB (p = .017), and CFP (p = .038), while both HbO and HbT were significantly reduced in RFP (p = .038, p = .045, respectively). Dynamic SE of HbT was reduced in DRE patients in RFP during minutes 2 to 6. Fitting to the frontal MREG and NIRS results, DRE patients showed a significant increase in SE_EEG in FB in fronto-central and parieto-occipital regions, in VLFP in parieto-central region, accompanied with a significant decrease in RFP in frontal pole and parietal and occipital (O2, Oz) regions.

Conclusion

This is the first study to show altered spectral entropy from synchronous MREG, EEG, and NIRS in DRE patients. Higher SE_MREG in DRE patients in anterior cingulate gyrus together with SE_EEG and SE_NIRS results in 0.12–0.4 Hz can be linked to altered parasympathetic function and respiratory pulsations in the brain. Higher SE_MREG in thalamus in DRE patients is connected to disturbances in anatomical and functional connections in epilepsy. Findings suggest that spectral irregularity of both electrophysiological and hemodynamic signals are altered in specific way depending on the physiological frequency range.

•

Simultaneous imaging methods indicate spectral irregularity in neurovascular and electrophysiological brain pulsations in DRE.

•

Altered spectral entropy in EEG, NIRS and BOLD indicate dysfunctional brain pulsations in respiratory frequency in epilepsy.

•

Spectral irregularity (0-5 Hz) of BOLD in right thalamus supports previous structural and functional findings in epilepsy.

Conductive gel bridge sensor for motion tracking in simultaneous EEG-fMRI recordings

EEG-fMRI allows the localization of the hemodynamic correlates of neural activity and has been shown to be useful as a diagnostic tool in pre-surgical evaluation of refractory epilepsy. However, EEG recordings may be highly contaminated by artifacts induced by movements inside the magnetic field thus rendering the scan difficult for interpretation. Existing methods for motion correction require additional equipment or hardware modification. We introduce a simple method for motion artifact detection, the conductive gel bridge sensor (CGBS), easily applicable using the standard setup. We report examples of CGBS use in two patients with epilepsy and demonstrate the method's ability to successfully differentiate between epochs of brain activity and those of movement.

Altered physiological brain variation in drug-resistant epilepsy

INTRODUCTION

Functional magnetic resonance imaging (fMRI) combined with simultaneous electroencephalography (EEG-fMRI) has become a major tool in mapping epilepsy sources. In the absence of detectable epileptiform activity, the resting state fMRI may still detect changes in the blood oxygen level-dependent signal, suggesting intrinsic alterations in the underlying brain physiology.

METHODS

In this study, we used coefficient of variation (CV) of critically sampled 10 Hz ultra-fast fMRI (magnetoencephalography, MREG) signal to compare physiological variance between healthy controls (n = 10) and patients (n = 10) with drug-resistant epilepsy (DRE).

RESULTS

We showed highly significant voxel-level (p < 0.01, TFCE-corrected) increase in the physiological variance in DRE patients. At individual level, the elevations range over three standard deviations (σ) above the control mean (μ) CV_MREG values solely in DRE patients, enabling patient-specific mapping of elevated physiological variance. The most apparent differences in group-level analysis are found on white matter, brainstem, and cerebellum. Respiratory (0.12-0.4 Hz) and very-low-frequency (VLF = 0.009-0.1 Hz) signal variances were most affected.

CONCLUSIONS

The CV_MREG increase was not explained by head motion or physiological cardiorespiratory activity, that is, it seems to be linked to intrinsic physiological pulsations. We suggest that intrinsic brain pulsations play a role in DRE and that critically sampled fMRI may provide a powerful tool for their identification.

Abnormal Blood Oxygen Level-Dependent Fluctuations in Focal Cortical Dysplasia and the Perilesional Zone: Initial Findings

BACKGROUND AND PURPOSE

Focal cortical dysplasia is a common cause of intractable epilepsy for which neurosurgery is an option. Delineations of a focal cortical dysplasia lesion on structural brain images may not necessarily reflect the functional borders of normal tissue. Our objective was to determine whether abnormalities in spontaneous blood oxygen level-dependent fluctuations arise in focal cortical dysplasia lesions and proximal regions.

MATERIALS AND METHODS

Fourteen patients with focal cortical dysplasia-related epilepsy and 16 healthy controls underwent structural and resting-state functional MR imaging. Three known blood oxygen level-dependent measures were determined, including the fractional amplitude of low-frequency fluctuations, regional homogeneity, and wavelet entropy. These measures were evaluated in the lesion and perilesional zone and normalized to the contralateral cortex of patients with focal cortical dysplasia and healthy controls.

RESULTS

Patients showed significantly decreased fractional amplitude of low-frequency fluctuations and increased wavelet entropy in the focal cortical dysplasia lesion and the perilesional zone (≤2 cm) relative to the contralateral homotopic cortex and the same regions in healthy controls. Regional homogeneity was significantly increased in the focal cortical dysplasia lesion compared with the contralateral homotopic cortex and healthy controls.

CONCLUSIONS

Abnormalities in spontaneous blood oxygen level-dependent fluctuations were seen up to 2 cm distant from the radiologically visible boundary. It was demonstrated that functional boundaries go beyond structural boundaries of focal cortical dysplasia lesions. Validation is required to reveal whether this information is valuable for surgical planning and outcome evaluation of focal cortical dysplasia lesions and comparing current results with electrophysiologic analysis.

Multimodal, noninvasive seizure network mapping software: A novel tool for preoperative epilepsy evaluation

BACKGROUND

Despite rigorous preoperative evaluation, epilepsy surgery achieves seizure freedom in only two-thirds of cases. Current preoperative evaluation does not include a detailed network analysis despite the association of network-level changes with epilepsy.

OBJECTIVE

We sought to create a software algorithm to map individualized epilepsy networks by combining noninvasive electroencephalography (EEG) source localization and nonconcurrent resting state functional magnetic resonance imaging (rsfMRI).

METHODS

Scalp EEG and rsfMRI data were acquired for three sample cases: one healthy control case, one case of right temporal lobe epilepsy, and one case of bitemporal seizure onset. Data from rsfMRI were preprocessed, and a time-series function was extracted. Connection coefficients were used to threshold out spurious connections and model global functional networks in a 3D map. Epileptic discharges were localized using a forward model of cortical mesh dipoles followed by an empirical Bayesian approach of inverse source reconstruction and co-registered with rsfMRI. Co-activating brain regions were mapped.

RESULTS

Three illustrative sample cases are presented. In the healthy control case, the software showed symmetrical global connectivity. In the right temporal lobe epilepsy case, asymmetry was found in the global connectivity metrics with a paucity of connectivity ipsilateral to the epileptogenic cortex. The superior longitudinal fasciculus, uncinate fasciculus, and commissural fibers connecting disparate and discontinuous cortical regions involved in the epilepsy network were visualized. In the case with bitemporal lobe epilepsy, global connectivity was symmetric. It showed a network of correlating cortical activity local to epileptogenic tissue in both temporal lobes. The network involved white matter tracks in a similar pattern to those seen in the right temporal case.

CONCLUSIONS

This modeling algorithm allows better definition of the global brain network and potentially demonstrates differences in connectivity between an epileptic and a non-epileptic brain. This finding may be useful for mapping cortico-cortical connections representing the putative epilepsy networks. With this methodology, we localized the epileptogenic brain and showed network asymmetry and long-distance cortical co-activation. This software tool is the first to use a multimodal, nonconcurrent, and noninvasive approach to model and visualize the epilepsy network.

DeepIED: An epileptic discharge detector for EEG-fMRI based on deep learning

Presurgical evaluation that can precisely delineate the epileptogenic zone (EZ) is one important step for successful surgical resection treatment of refractory epilepsy patients. The noninvasive EEG-fMRI recording technique combined with general linear model (GLM) analysis is considered an important tool for estimating the EZ. However, the manual marking of interictal epileptic discharges (IEDs) needed in this analysis is challenging and time-consuming because the quality of the EEG recorded inside the scanner is greatly deteriorated compared to the usual EEG obtained outside the scanner. This is one of main impediments to the widespread use of EEG-fMRI in epilepsy. We propose a deep learning based semi-automatic IED detector that can find the candidate IEDs in the EEG recorded inside the scanner which resemble sample IEDs marked in the EEG recorded outside the scanner. The manual marking burden is greatly reduced as the expert need only edit candidate IEDs. The model is trained on data from 30 patients. Validation of IEDs detection accuracy on another 37 consecutive patients shows our method can improve the median sensitivity from 50.0% for the previously proposed template-based method to 84.2%, with false positive rate as 5 events/min. Reproducibility validation on 15 patients is applied to evaluate if our method can produce similar hemodynamic response maps compared with the manual marking ground truth results. We explore the concordance between the maximum hemodynamic response and the intracerebral EEG defined EZ and find that both methods produce similar percentage of concordance (76.9%, 10 out of 13 patients, electrode was absent in the maximum hemodynamic response in two patients). This tool will make EEG-fMRI analysis more practical for clinical usage.

•

A deep learning based epileptic discharge detector for EEG-fMRI is proposed.

•

The burden of manually marking epileptic discharges is greatly reduced.

•

Our method can produce similar EEG-fMRI results compared with traditional method.

Dynamic coupling between fMRI local connectivity and interictal EEG in focal epilepsy: A wavelet analysis approach

Simultaneous scalp EEG-fMRI recording is a noninvasive neuroimaging technique for combining electrophysiological and hemodynamic aspects of brain function. Despite the time-varying nature of both measurements, their relationship is usually considered as time-invariant. The aim of this study was to detect direct associations between scalp-recorded EEG and regional changes of hemodynamic brain connectivity in focal epilepsy through a time-frequency paradigm. To do so, we developed a voxel-wise framework that analyses wavelet coherence between dynamic regional phase synchrony (DRePS, calculated from fMRI) and band amplitude fluctuation (BAF) of a target EEG electrode with dominant interictal epileptiform discharges (IEDs). As a proof of concept, we applied this framework to seven patients with focal epilepsy. The analysis produced patient-specific spatial maps of DRePS-BAF coupling, which highlight regions with a strong link between EEG power and local fMRI connectivity. Although we observed DRePS-BAF coupling proximate to the suspected seizure onset zone in some patients, our results suggest that DRePS-BAF is more likely to identify wider 'epileptic networks'. We also compared DRePS-BAF with standard EEG-fMRI analysis based on general linear modelling (GLM). There was, in general, little overlap between the DRePS-BAF maps and GLM maps. However, in some subjects the spatial clusters revealed by these two analyses appeared to be adjacent, particularly in medial posterior cortices. Our findings suggest that (1) there is a strong time-varying relationship between local fMRI connectivity and interictal EEG power in focal epilepsy, and (2) that DRePS-BAF reflect different aspects of epileptic network activity than standard EEG-fMRI analysis. These two techniques, therefore, appear to be complementary. Hum Brain Mapp 38:5356-5374, 2017. © 2017 Wiley Periodicals, Inc.

The Advents of Hybrid Imaging Modalities: A New Era in Neuroimaging Applications

Hybrid Imaging modalities have shown great potential in medical imaging and diagnosis. A more comprehensive and targeted view of neurological disorders can be achieved by blending the anatomical and functional perspectives through hybridization. With consistently improving technologies, there have been many developments in fused imaging techniques over the past few decades. This article provides an overview of various bimodal and trimodal hybrid imaging techniques being developed and explored for neuroimaging applications. Recent advancements and potentials are discussed for single photon emission computed tomography-computed tomography (SPECT-CT), positron emission tomography-CT (PET-CT), PET-magnetic resonance imaging (PET-MRI), electroencephalography-functional magnetic resonance imaging (EEG-fMRI), magnetoencephalography-fMRI (MEG-fMRI), EEG-near-infrared spectroscopy (EEG-NIRS), magnetic resonance-PET-EEG (MR-PET-EEG) and MR-PET-CT in the perspective of neuroimaging. A comparison of these hybrid approaches is provided on a single platform to analyze their performance on the basis of several common factors essential for imaging and analyzing neurological disorders and in vivo molecular processes. This article also provides an overview of recently developed advanced imaging technologies that are being hybridized with other imaging modalities and being explored as potential techniques for neuroscience. Novel approaches and clinical applications of hybrid neuroimaging are anticipated with inclusion of new technologies, better sensing capabilities, multimodal probes, and improved hybridization.

Contributions of EEG-fMRI to Assessing the Epileptogenicity of Focal Cortical Dysplasia

Purpose: To examine the ability of the BOLD response to EEG spikes to assess the epileptogenicity of the lesion in patients with focal cortical dysplasia (FCD).

Method: Patients with focal epilepsy and FCD who underwent 3T EEG-fMRI from 2006 to 2010 were included. Diagnosis of FCD was based on neuroradiology (MRI+), or histopathology in MRI-negative cases (MRI−). Patients underwent 120 min EEG-fMRI recording session. Spikes similar to those recorded outside the scanner were marked in the filtered EEG. The lesion (in MRI+) or the removed cortex (in MRI−) was marked on the anatomical T1 sequence, blindly to the BOLD response, after reviewing the FLAIR images. For each BOLD response we assessed the concordance with the spike field and with the lesion in MRI+ or the removed cortex in MRI−. BOLD responses were considered “concordant” if the maximal t-value was inside the marking. Follow-up after resection was used as gold-standard.

Results: Twenty patients were included (13 MRI+, 7 MRI−), but in seven the EEG was not active or there were artifacts during acquisition. In all 13 studied patients, at least one BOLD response was concordant with the spike field; in 9/13 (69%) at least one BOLD response was concordant with the lesion: in 6/7 (86%) MRI+ and in 3/6 (50%) MRI− patients.

Conclusions: Our study shows a high level of concordance between FCD and BOLD response. This data could provide useful information especially for MRI negative patients. Moreover, it shows in almost all FCD patients, a metabolic involvement of remote cortical or subcortical structures, corroborating the concept of epileptic network.

Multichannel continuous electroencephalography-functional near-infrared spectroscopy recording of focal seizures and interictal epileptiform discharges in human epilepsy: a review

Functional near-infrared spectroscopy (fNIRS) has emerged as a promising neuroimaging technique as it allows noninvasive and long-term monitoring of cortical hemodynamics. Recent work by our group and others has revealed the potential of fNIRS, combined with electroencephalography (EEG), in the context of human epilepsy. Hemodynamic brain responses attributed to epileptic events, such as seizures and interictal epileptiform discharges (IEDs), are routinely observed with a good degree of statistical significance and in concordance with clinical presentation. Recording done with over 100 channels allows sufficiently large coverage of the epileptic focus and other areas. Three types of seizures have been documented: frontal lobe seizures, temporal lobe seizures, and posterior seizures. Increased oxygenation was observed in the epileptic focus in most cases, while rapid but similar hemodynamic variations were identified in the contralateral homologous region. While investigating IEDs, it was shown that their hemodynamic effect is observable with fNIRS, that their response is associated with significant (inhibitive) nonlinearities, and that the sensitivity and specificity of fNIRS to localize the epileptic focus can be estimated in a sample of 40 patients. This paper first reviews recent EEG-fNIRS developments in epilepsy research and then describes applications to the study of focal seizures and IEDs.

Marker-based ballistocardiographic artifact correction improves spike identification in EEG-fMRI of focal epilepsy patients

OBJECTIVES

Ballistocardiographic (BCG) artifacts resemble interictal epileptic discharges (IEDs) and can lead to incorrect IED identification in EEG-fMRI. This study investigates IEDs marked in EEGs corrected using information from a moiré phase tracking (MPT) marker.

METHODS

EEG-fMRI from 18 patients was processed with conventional methods for BCG removal, while 9 patients used a MPT marker. IEDs were marked first without ECG information. In a second review, suspicious IEDs synchronous with the BCG were discarded. After each review, an event-related fMRI analysis was performed on the marked IEDs.

RESULTS

No difference was found in the proportion of suspicious IEDs in the 2 patient groups. However, the distribution of IED timings was significantly related to the cardiac cycle in 11 of 18 patients recorded without MPT marker, but only 2 of 9 patients with marker. In patients recorded without marker, failing to discard suspicious IEDs led to more inaccurate fMRI maps and more distant activations.

CONCLUSIONS

BCG artifact correction based on MPT recordings allowed a more straightforward identification of IEDs that did not require ECG information in the large majority of patients.

SIGNIFICANCE

Marker-based ballistocardiographic artifact correction greatly facilitates the study of the generators of interictal discharges with EEG-fMRI.