Independent component analysis of interictal fMRI in focal epilepsy: Comparison with general linear model-based EEG-correlated fMRI

Association of Normative and Non-Normative Brain Networks With Cognitive Function in Patients With Temporal Lobe Epilepsy

BACKGROUND AND OBJECTIVES

Despite their temporal lobe pathology, a significant subgroup of patients with temporal lobe epilepsy (TLE) is able to maintain normative cognitive functioning. In this study, we identify patients with TLE with intact vs impaired neurocognitive profiles and interrogate for the presence of both normative and highly individual intrinsic connectivity networks (ICNs)-all toward understanding the transition from impaired to intact neurocognitive status.

METHODS

This retrospective cross-sectional study included patients with TLE and matched healthy controls (HCs) from the Thomas Jefferson Comprehensive Epilepsy Center. Functional MRI data were decomposed using independent component analysis to obtain individualized ICNs. In this article, we calculated the degree of match between individualized ICNs and canonical ICNs (e.g., 17 resting-state networks by Yeo et al.) and divided each participant's ICNs into normative or non-normative status based on the degree of match.

RESULTS

100 patients with TLE (mean age 42.0 [SD: 13.7] years, 47 women) and 92 HCs were included in this study. We found that the individualized networks matched to the canonical networks less well in the cognitively impaired (n = 24) compared with the cognitively intact (n = 63) patients with TLE by 2-way mixed-measures analysis of variance (impaired vs intact mean difference [MD] -0.165 [-0.317, -0.013], p = 0.028). The cognitively impaired patients showed significant abnormalities in the profiles of both normative (impaired vs intact MD -0.537 [-0.998, -0.076], p = 0.017, intact vs HC MD -0.221 [-0.536, 0.924], p = 0.220, and impaired vs HC MD -0.759 [-1.200, -0.319], p < 0.001) and non-normative networks (impaired vs intact MD 0.484 [0.030, 0.937], p = 0.033, intact vs HC MD 0.369 [0.059, 0.678], p = 0.014, and impaired vs HC MD 0.853 [0.419, 1.286], p < 0.001) while the intact patients showed abnormalities only in non-normative networks. At the same time, we found that normative networks held a strong, positive association with the neuropsychological measures, with this association negative in non-normative networks.

DISCUSSION

Our data demonstrated that significant cognitive deficits are associated with the status of both canonical and highly individual ICNs, making clear that the transition from intact to impaired cognitive status is not simply the result of disruption to normative brain networks.

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.

EEG-fMRI of epileptiform discharges: non-invasive investigation of the whole brain

Simultaneous EEG-fMRI is a unique and non-invasive method for investigating epileptic activity. Interictal epileptiform discharge-related EEG-fMRI provides cortical and subcortical blood oxygen level-dependent (BOLD) signal changes specific to epileptic discharges. As a result, EEG-fMRI has revealed insights into generators and networks involved in epileptic activity in different types of epilepsy, demonstrating-for instance-the implication of the thalamus in human generalized spike and wave discharges and the role of the Default Mode Network (DMN) in absences and focal epilepsy, and proposed a mechanism for the cortico-subcortical interactions in Lennox-Gastaut syndrome discharges. EEG-fMRI can find deep sources of epileptic activity not available to scalp EEG or MEG and provides critical new information to delineate the epileptic focus when considering surgical treatment or electrode implantation. In recent years, methodological advances, such as artifact removal and automatic detection of events have rendered this method easier to implement, and its clinical potential has since been established by evidence of the impact of BOLD response on clinical decision-making and of the relationship between concordance of BOLD responses with extent of resection and surgical outcome. This review presents the recent developments in EEG-fMRI methodology and EEG-fMRI studies in different types of epileptic disorders as follows: EEG-fMRI acquisition, gradient and pulse artifact removal, statistical analysis, clinical applications, pre-surgical evaluation, altered physiological state in generalized genetic epilepsy, and pediatric EEG-fMRI studies.

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.

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.

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.

Machine learning identifies "rsfMRI epilepsy networks" in temporal lobe epilepsy

OBJECTIVES

Experimental models have provided compelling evidence for the existence of neural networks in temporal lobe epilepsy (TLE). To identify and validate the possible existence of resting-state "epilepsy networks," we used machine learning methods on resting-state functional magnetic resonance imaging (rsfMRI) data from 42 individuals with TLE.

METHODS

Probabilistic independent component analysis (PICA) was applied to rsfMRI data from 132 subjects (42 TLE patients + 90 healthy controls) and 88 independent components (ICs) were obtained following standard procedures. Elastic net-selected features were used as inputs to support vector machine (SVM). The strengths of the top 10 networks were correlated with clinical features to obtain "rsfMRI epilepsy networks."

RESULTS

SVM could classify individuals with epilepsy with 97.5% accuracy (sensitivity = 100%, specificity = 94.4%). Ten networks with the highest ranking were found in the frontal, perisylvian, cingulo-insular, posterior-quadrant, thalamic, cerebello-thalamic, and temporo-thalamic regions. The posterior-quadrant, cerebello-thalamic, thalamic, medial-visual, and perisylvian networks revealed significant correlation (r > 0.40) with age at onset of seizures, the frequency of seizures, duration of illness, and a number of anti-epileptic drugs.

CONCLUSIONS

IC-derived rsfMRI networks contain epilepsy-related networks and machine learning methods are useful in identifying these networks in vivo. Increased network strength with disease progression in these "rsfMRI epilepsy networks" could reflect epileptogenesis in TLE.

KEY POINTS

• ICA of resting-state fMRI carries disease-specific information about epilepsy. • Machine learning can classify these components with 97.5% accuracy. • "Subject-specific epilepsy networks" could quantify "epileptogenesis" in vivo.

Correlating Resting-State Functional Magnetic Resonance Imaging Connectivity by Independent Component Analysis-Based Epileptogenic Zones with Intracranial Electroencephalogram Localized Seizure Onset Zones and Surgical Outcomes in Prospective Pediatric Intractable Epilepsy Study

The purpose of this study was to prospectively investigate the agreement between the epileptogenic zone(s) (EZ) localization by resting-state functional magnetic resonance imaging (rs-fMRI) and the seizure onset zone(s) (SOZ) identified by intracranial electroencephalogram (ic-EEG) using novel differentiating and ranking criteria of rs-fMRI abnormal independent components (ICs) in a large consecutive heterogeneous pediatric intractable epilepsy population without an a priori alternate modality informing EZ localization or prior declaration of total SOZ number. The EZ determination criteria were developed by using independent component analysis (ICA) on rs-fMRI in an initial cohort of 350 pediatric patients evaluated for epilepsy surgery over a 3-year period. Subsequently, these rs-fMRI EZ criteria were applied prospectively to an evaluation cohort of 40 patients who underwent ic-EEG for SOZ identification. Thirty-seven of these patients had surgical resection/disconnection of the area believed to be the primary source of seizures. One-year seizure frequency rate was collected postoperatively. Among the total 40 patients evaluated, agreement between rs-fMRI EZ and ic-EEG SOZ was 90% (36/40; 95% confidence interval [CI], 0.76-0.97). Of the 37 patients who had surgical destruction of the area believed to be the primary source of seizures, 27 (73%) rs-fMRI EZ could be classified as true positives, 7 (18%) false positives, and 2 (5%) false negatives. Sensitivity of rs-fMRI EZ was 93% (95% CI 78-98%) with a positive predictive value of 79% (95% CI, 63-89%). In those with cryptogenic localization-related epilepsy, agreement between rs-fMRI EZ and ic-EEG SOZ was 89% (8/9; 95% CI, 0.52-99), with no statistically significant difference between the agreement in the cryptogenic and symptomatic localization-related epilepsy subgroups. Two children with negative ic-EEG had removal of the rs-fMRI EZ and were seizure free 1 year postoperatively. Of the 33 patients where at least 1 rs-fMRI EZ agreed with the ic-EEG SOZ, 24% had at least 1 additional rs-fMRI EZ outside the resection area. Of these patients with un-resected rs-fMRI EZ, 75% continued to have seizures 1 year later. Conversely, among 75% of patients in whom rs-fMRI agreed with ic-EEG SOZ and had no anatomically separate rs-fMRI EZ, only 24% continued to have seizures 1 year later. This relationship between extraneous rs-fMRI EZ and seizure outcome was statistically significant (p = 0.01). rs-fMRI EZ surgical destruction showed significant association with postoperative seizure outcome. The pediatric population with intractable epilepsy studied prospectively provides evidence for use of resting-state ICA ranking criteria, to identify rs-fMRI EZ, as developed by the lead author (V.L.B.). This is a high yield test in this population, because no seizure nor particular interictal epilepiform activity needs to occur during the study. Thus, rs-fMRI EZ detected by this technique are potentially informative for epilepsy surgery evaluation and planning in this population. Independent of other brain function testing modalities, such as simultaneous EEG-fMRI or electrical source imaging, contextual ranking of abnormal ICs of rs-fMRI localized EZs correlated with the gold standard of SOZ localization, ic-EEG, across the broad range of pediatric epilepsy surgery candidates, including those with cryptogenic epilepsy.

ICN_Atlas: Automated description and quantification of functional MRI activation patterns in the framework of intrinsic connectivity networks

Generally, the interpretation of functional MRI (fMRI) activation maps continues to rely on assessing their relationship to anatomical structures, mostly in a qualitative and often subjective way. Recently, the existence of persistent and stable brain networks of functional nature has been revealed; in particular these so-called intrinsic connectivity networks (ICNs) appear to link patterns of resting state and task-related state connectivity. These networks provide an opportunity of functionally-derived description and interpretation of fMRI maps, that may be especially important in cases where the maps are predominantly task-unrelated, such as studies of spontaneous brain activity e.g. in the case of seizure-related fMRI maps in epilepsy patients or sleep states. Here we present a new toolbox (ICN_Atlas) aimed at facilitating the interpretation of fMRI data in the context of ICN. More specifically, the new methodology was designed to describe fMRI maps in function-oriented, objective and quantitative way using a set of 15 metrics conceived to quantify the degree of 'engagement' of ICNs for any given fMRI-derived statistical map of interest. We demonstrate that the proposed framework provides a highly reliable quantification of fMRI activation maps using a publicly available longitudinal (test-retest) resting-state fMRI dataset. The utility of the ICN_Atlas is also illustrated on a parametric task-modulation fMRI dataset, and on a dataset of a patient who had repeated seizures during resting-state fMRI, confirmed on simultaneously recorded EEG. The proposed ICN_Atlas toolbox is freely available for download at http://icnatlas.com and at http://www.nitrc.org for researchers to use in their fMRI investigations.

Two-Dimensional Temporal Clustering Analysis for Patients with Epilepsy: Detecting Epilepsy-Related Information in EEG-fMRI Concordant, Discordant and Spike-Less Patients

EEG acquired simultaneously with fMRI (EEG-fMRI) is a multimodal method that has shown promise in mapping the seizure onset zone in patients with focal epilepsy. However, there are many instances when this method is unsuccessful or not applicable, and other data driven fMRI methods may be utilized. One such method is the two-dimensional temporal clustering analysis (2dTCA). In this study we compared the classic EEG-fMRI and 2dTCA performance in mapping regions related to the seizure onset region in 18 focal epilepsy patients (12 presenting interictal epileptiform discharges (IEDs), during EEG-fMRI acquisition) with Engel I or II surgical outcome. Activation maps of both 2dTCA timing outputs (positive and negative histograms) and EEG detected IEDs were computed and compared to the region of epilepsy surgical resection. Patients were evaluated in three categories based on frequency of EEG detected spiking during the MRI. EEG-fMRI maps were concordant to the epilepsy region in 5/12 subjects, four with frequent IEDs on EEG. The 2dTCA was successful in mapping 13/18 patients including 3/6 with no IEDs detected (10/12 with IEDs detected). The epilepsy-related activities were successfully mapped by both methods in only 4/12 patients. This work suggests that the epilepsy-related information detected by each method may be different: while EEG-fMRI is more accurate in patients with high rather than lower numbers of EEG detected IEDs; 2dTCA can be useful in evaluating patients even when no concurrent EEG spikes are detected or EEG-fMRI is not effective. Therefore, our results support that 2dTCA might be an alternative for mapping epilepsy-related BOLD activity in negative EEG-fMRI (6/7 patients) and spike-less patients.

Altered brain function in new onset childhood acute lymphoblastic leukemia before chemotherapy: A resting-state fMRI study

OBJECTIVE

Cognitive impairments had been reported in childhood acute lymphoblastic leukemia, what caused the impairments needed to be demonstrated, chemotherapy-related or the disease itself. The primary aim of this exploratory investigation was to determine if there were changes in brain function of children with acute lymphoblastic leukemia before chemotherapy.

METHODS

In this study, we advanced a measure named regional homogeneity to evaluate the resting-state brain activities, intelligence quotient test was performed at same time. Using regional homogeneity, we first investigated the resting state brain function in patients with new onset childhood acute lymphoblastic leukemia before chemotherapy, healthy children as control.

RESULTS

The decreased ReHo values were mainly founded in the default mode network and left frontal lobe, bilateral inferior parietal lobule, bilateral temporal lobe, bilateral occipital lobe, precentral gyrus, bilateral cerebellum in the newly diagnosed acute lymphoblastic leukemia patients compared with the healthy control. While in contrast, increased ReHo values were mainly shown in the right frontal lobe (language area), superior frontal gyrus-R, middle frontal gyrus-R and inferior parietal lobule-R for acute lymphoblastic leukemia patients group. There were no significant differences for intelligence quotient measurements between the acute lymphoblastic leukemia patient group and the healthy control in performance intelligence quotient, verbal intelligence quotient, total intelligence quotient.

CONCLUSION

The altered brain functions are associated with cognitive change and language, it is suggested that there may be cognition impairment before the chemotherapy. Regional homogeneity by functional magnetic resonance image is a sensitive way for early detection on brain damage in childhood acute lymphoblastic leukemia.

Hand classification of fMRI ICA noise components

We present a practical "how-to" guide to help determine whether single-subject fMRI independent components (ICs) characterise structured noise or not. Manual identification of signal and noise after ICA decomposition is required for efficient data denoising: to train supervised algorithms, to check the results of unsupervised ones or to manually clean the data. In this paper we describe the main spatial and temporal features of ICs and provide general guidelines on how to evaluate these. Examples of signal and noise components are provided from a wide range of datasets (3T data, including examples from the UK Biobank and the Human Connectome Project, and 7T data), together with practical guidelines for their identification. Finally, we discuss how the data quality, data type and preprocessing can influence the characteristics of the ICs and present examples of particularly challenging datasets.

Delineating potential epileptogenic areas utilizing resting functional magnetic resonance imaging (fMRI) in epilepsy patients

Seizure localization includes neuroimaging like electroencephalogram, and magnetic resonance imaging (MRI) with limited ability to characterize the epileptogenic network. Temporal clustering analysis (TCA) characterizes epileptogenic network congruent with interictal epileptiform discharges by clustering together voxels with transient signals. We generated epileptogenic areas for 12 of 13 epilepsy patients with TCA, congruent with different areas of seizure onset. Resting functional MRI (fMRI) scans are noninvasive, and can be acquired quickly, in patients with different levels of severity and function. Analyzing resting fMRI data using TCA is quick and can complement clinical methods to characterize the epileptogenic network.

Epileptic networks are strongly connected with and without the effects of interictal discharges

OBJECTIVE

Epilepsy is increasingly considered as the dysfunction of a pathologic neuronal network (epileptic network) rather than a single focal source. We aimed to assess the interactions between the regions that comprise the epileptic network and to investigate their dependence on the occurrence of interictal epileptiform discharges (IEDs).

METHODS

We analyzed resting state simultaneous electroencephalography-functional magnetic resonance imaging (EEG-fMRI) recordings in 10 patients with drug-resistant focal epilepsy with multifocal IED-related blood oxygen level-dependent (BOLD) responses and a maximum t-value in the IED field. We computed functional connectivity (FC) maps of the epileptic network using two types of seed: (1) a 10-mm diameter sphere centered in the global maximum of IED-related BOLD map, and (2) the independent component with highest correlation to the IED-related BOLD map, named epileptic component. For both approaches, we compared FC maps before and after regressing out the effect of IEDs in terms of maximum and mean t-values and percentage of map overlap.

RESULTS

Maximum and mean FC maps t-values were significantly lower after regressing out IEDs at the group level (p < 0.01). Overlap extent was 85% ± 12% and 87% ± 12% when the seed was the 10-mm diameter sphere and the epileptic component, respectively.

SIGNIFICANCE

Regions involved in a specific epileptic network show coherent BOLD fluctuations independent of scalp EEG IEDs. FC topography and strength is largely preserved by removing the IED effect. This could represent a signature of a sustained pathologic network with contribution from epileptic activity invisible to the scalp EEG.

Independent Component Analysis of EEG-fMRI data for studying epilepsy and epileptic seizures

Here we present a method for classifying fMRI independent components (ICs) by using an optimized algorithm for the individuation of noisy signals from sources of interest. The method was applied to estimate brain activations from combined EEG-fMRI data for the exploration of epilepsy. Spatial ICA was performed using the above-mentioned optimized algorithm and other three popular algorithms. ICs were sorted considering the value: of the coefficients of determination R2, obtained from the multiple regression analysis with morphometric maps of cerebral matter; of the kurtosis, which features the signal energy. The validation of the method was performed comparing the brain activations obtained with those resulted using the General Linear Model (GLM). The ICA-derived activations in different datasets comprised subareas of the GLM-revealed activations, even if the volume and the shape of activated areas do not correspond exactly. The method proposed also detects additional negative regions implicated in a default mode of brain activity, and not clearly identified by GLM. Compared with a traditional GLM approach, the ICA one provides a flexible way to analyze fMRI data that reduces the assumptions placed upon the hemodynamic response of the brain and the temporal constrains.

Network analysis of EEG related functional MRI changes due to medication withdrawal in focal epilepsy

Anti-epileptic drugs (AEDs) have a global effect on the neurophysiology of the brain which is most likely reflected in functional brain activity recorded with EEG and fMRI. These effects may cause substantial inter-subject variability in studies where EEG correlated functional MRI (EEG-fMRI) is used to determine the epileptogenic zone in patients who are candidate for epilepsy surgery. In the present study the effects on resting state fMRI are quantified in conditions with AED administration and after withdrawal of AEDs. EEG-fMRI data were obtained from 10 patients in the condition that the patient was on the steady-state maintenance doses of AEDs as prescribed (condition A) and after withdrawal of AEDs (condition B), at the end of a clinically standard pre-surgical long term video-EEG monitoring session. Resting state networks (RSN) were extracted from fMRI. The epileptic component (ICE) was identified by selecting the RSN component with the largest overlap with the EEG-fMRI correlation pattern. Changes in RSN functional connectivity between conditions A and B were quantified. EEG-fMRI correlation analysis was successful in 30% and 100% of the cases in conditions A and B, respectively. Spatial patterns of ICEs are comparable in conditions A and B, except for one patient for whom it was not possible to identify the ICE in condition A. However, the resting state functional connectivity is significantly increased in the condition after withdrawal of AEDs (condition B), which makes resting state fMRI potentially a new tool to study AED effects. The difference in sensitivity of EEG-fMRI in conditions A and B, which is not related to the number of epileptic EEG events occurring during scanning, could be related to the increased functional connectivity in condition B.

A Comparison of Independent Component Analysis (ICA) of fMRI and Electrical Source Imaging (ESI) in Focal Epilepsy Reveals Misclassification Using a Classifier

Interictal epileptiform discharges (IEDs) can produce haemodynamic responses that can be detected by electroencephalography-functional magnetic resonance imaging (EEG-fMRI) using different analysis methods such as the general linear model (GLM) of IEDs or independent component analysis (ICA). The IEDs can also be mapped by electrical source imaging (ESI) which has been demonstrated to be useful in presurgical evaluation in a high proportion of cases with focal IEDs. ICA advantageously does not require IEDs or a model of haemodynamic responses but its use in EEG-fMRI of epilepsy has been limited by its ability to separate and select epileptic components. Here, we evaluated the performance of a classifier that aims to filter all non-BOLD responses and we compared the spatial and temporal features of the selected independent components (ICs). The components selected by the classifier were compared to those components selected by a strong spatial correlation with ESI maps of IED sources. Both sets of ICs were subsequently compared to a temporal model derived from the convolution of the IEDs (derived from the simultaneously acquired EEG) with a standard haemodynamic response. Selected ICs were compared to the patients' clinical information in 13 patients with focal epilepsy. We found that the misclassified ICs clearly related to IED in 16/25 cases. We also found that the classifier failed predominantly due to the increased spectral range of fMRIs temporal responses to IEDs. In conclusion, we show that ICA can be an efficient approach to separate responses related to epilepsy but that contemporary classifiers need to be retrained for epilepsy data. Our findings indicate that, for ICA to contribute to the analysis of data without IEDs to improve its sensitivity, classification strategies based on data features other than IC time course frequency is required.

Methods and utility of EEG-fMRI in epilepsy

Brain activity data in general and more specifically in epilepsy can be represented as a matrix that includes measures of electrophysiology, anatomy and behaviour. Each of these sub-matrices has a complex interaction depending upon the brain state i.e., rest, cognition, seizures and interictal periods. This interaction presents significant challenges for interpretation but also potential for developing further insights into individual event types. Successful treatments in epilepsy hinge on unravelling these complexities, and also on the sensitivity and specificity of methods that characterize the nature and localization of underlying physiological and pathological networks. Limitations of pharmacological and surgical treatments call for refinement and elaboration of methods to improve our capability to localise the generators of seizure activity and our understanding of the neurobiology of epilepsy. Simultaneous electroencephalography and functional magnetic resonance imaging (EEG-fMRI), by potentially circumventing some of the limitations of EEG in terms of sensitivity, can allow the mapping of haemodynamic networks over the entire brain related to specific spontaneous and triggered epileptic events in humans, and thereby provide new localising information. In this work we review the published literature, and discuss the methods and utility of EEG-fMRI in localising the generators of epileptic activity. We draw on our experience and that of other groups, to summarise the spectrum of information provided by an increasing number of EEG-fMRI case-series, case studies and group studies in patients with epilepsy, for its potential role to elucidate epileptic generators and networks. We conclude that EEG-fMRI provides a multidimensional view that contributes valuable clinical information to localize the epileptic focus with potential important implications for the surgical treatment of some patients with drug-resistant epilepsy, and insights into the resting state and cognitive network dynamics.

Using functional MRI alone for localization in focal epilepsy

In the present study, we developed a method for the purpose of localizing epilepsy related hemodynamic foci for patients suffering intractable focal epilepsy using resting state fMRI alone. We studied two groups of subjects: five patients with intractable focal epilepsy, and ten healthy volunteers performing motor tasks. Spatial independent component analysis (ICA) was performed on the fMRI alone data and a set of independent component (IC) selection criteria was developed to identify epilepsy related ICs. The method was then evaluated in the healthy group with motor tasks. In all five surgery patients, there was at least one identified IC concordant with surgical resection. In the motor task study of healthy subjects, our method revealed components with concordant spatial and temporal features as expected from the unilateral motor tasks. These results suggest the lateralization and localization value of fMRI alone in presurgical evaluation for patients with intractable unilateral focal epilepsy. The proposed method is noninvasive in nature and easy to implement. It has the potential to be incorporated in current presurgical workup for the diagnosis of intractable focal epilepsy patients.

Automatic selection of epileptic independent fMRI components

EEG-correlated fMRI analysis has proven to be useful in localizing regions of BOLD activation related to epileptic activity. However, as EEG does not always provide reliable information, purely fMRI-based data-driven techniques are invaluable. Recently, we have shown that independent component analysis (ICA) can extract sources related to the epileptic network even in such EEG-negative cases [1]. Moreover, these sources were shown to be informative with respect to the seizure onset zone (SOZ). In order to utilize this concept in clinical practice in a prospective manner, this work aims at developing an automatic technique for selecting the epileptic sources. The proposed approach applies a cascade of two classifiers. In the first step artifact related sources are discarded. In the second step the sources are characterized by four discriminative features and epileptic sources are selected from among other BOLD-related components. Our technique reaches a promising 77% specificity and provides concordant sources with the EEG-correlated fMRI activation maps or with the SOZ in 71% of the cases.

A reliable and time-saving semiautomatic spike-template-based analysis of interictal EEG-fMRI

OBJECTIVE

A prerequisite for the implementation of interictal electroencephalography-correlated functional magnetic resonance imaging (EEG-fMRI) in the presurgical work-up for epilepsy surgery is straightforward processing. We propose a new semi-automatic method as alternative for the challenging and time-consuming visual spike identification.

METHODS

Our method starts from a patient-specific spike-template, built by averaging spikes recorded on the EEG outside the scanner. Spatiotemporal cross-correlations between the template and the EEG measured during fMRI were calculated. To minimize false-positive detections, this time course of cross-correlations was binarized by means of a spike-template-specific threshold determined in healthy controls. To inform our model for statistical parametric mapping, this binarized regressor was convolved with the canonical hemodynamic response function. We validated our "template-based" method in 21 adult patients with refractory focal epilepsy with a well-defined epileptogenic zone and interictal spikes during EEG-fMRI. Sensitivity and specificity for detecting the epileptogenic zone were calculated and represented in receiver operating characteristic (ROC) curves. Our approach was compared with a previously proposed semiautomatic "topography-based" method that used the topographic amplitude distribution of spikes as a starting point for correlation-based fitting.

RESULTS

Good diagnostic performance could be reached with our template-based method. The optimal area under the ROC curve was 0.77. Diagnostic performance of the topography-based method was overall low.

SIGNIFICANCE

Our new template-based method is more standardized and time-saving than visual spike identification on intra-scanner EEG recordings, and preserves good diagnostic performance for detecting the epileptogenic zone.

Are Epilepsy-Related fMRI Components Dependent on the Presence of Interictal Epileptic Discharges in Scalp EEG?

Spatial independent component analysis (ICA) is increasingly being used to extract resting-state networks from fMRI data. Previous studies showed that ICA also reveals independent components (ICs) related to the seizure onset zone. However, it is currently unknown how these epileptic ICs depend on the presence of interictal epileptic discharges (IEDs) in the EEG. The goal of this study was to explore the relation between ICs obtained from fMRI epochs during the occurrence of IEDs in the EEG and those without IEDs. fMRI data sets with co-registered EEG were retrospectively selected of patients from whom the location of the epileptogenic zone was confirmed by outcome of surgery (n = 8). The fMRI data were split into two epochs: one with IEDs visible in scalp EEG and one without. Spatial ICA was applied to the fMRI data of each part separately. The maps of all resulting components were compared to the resection area and the EEG-fMRI correlation pattern by computing a spatial correlation coefficient to detect the epilepsy-related component. For all patients, except one, there was a remarkable resemblance between the epilepsy-related components selected during epochs with IEDs and those without IEDs. These findings suggest that epilepsy-related ICs are not dependent on the presence of IEDs in scalp EEG. Since these epileptic ICs showed partial overlap with resting-state networks of healthy volunteers (n = 10), our study supports the need for new ways to classify epileptic ICs.

De-noising with a SOCK can improve the performance of event-related ICA

Event-related ICA (eICA) is a partially data-driven analysis method for event-related fMRI that is particularly suited to analysis of simultaneous EEG-fMRI of patients with epilepsy. EEG-fMRI studies in epileptic patients are typically analyzed using the general linear model (GLM), often with assumption that the onset and offset of neuronal activity match EEG event onset and offset, the neuronal activation is sustained at a constant level throughout the epileptiform event and that associated fMRI signal changes follow the canonical HRF. The eICA method allows for less constrained analyses capable of detecting early, non-canonical responses. A key step of eICA is the initial deconvolution which can be confounded by various sources of structured noise present in the fMRI signal. To help overcome this, we have extend the eICA procedure by utilizing a fully standalone and automated fMRI de-noising procedure to process the fMRI data from an EEG-fMRI acquisition prior to running eICA. Specifically we first apply ICA to the entire fMRI time-series and use a classifier to remove noise-related components. The automated objective de-noiser, “Spatially Organized Component Klassificator” (SOCK) is used; it has previously been shown to distinguish a substantial fraction of noise from true activation, without rejecting the latter, in resting-state fMRI. A second ICA is then performed, this time on the event-related response estimates derived from the denoised data (according to the usual eICA procedure). We hypothesize that SOCK + eICA has the potential to be more sensitive than eICA alone. We test the effectiveness of SOCK by comparing activation obtained in an eICA analysis of EEG-fMRI data with and without the use of SOCK for 14 patients with rolandic epilepsy who exhibited stereotypical IEDs arising from a focus in the rolandic fissure.

Network Connectivity in Epilepsy: Resting State fMRI and EEG-fMRI Contributions

There is a growing body of evidence pointing toward large-scale networks underlying the core phenomena in epilepsy, from seizure generation to cognitive dysfunction or response to treatment. The investigation of networks in epilepsy has become a key concept to unlock a deeper understanding of the disease. Functional imaging can provide valuable information to characterize network dysfunction; in particular resting state fMRI (RS-fMRI), which is increasingly being applied to study brain networks in a number of diseases. In patients with epilepsy, network connectivity derived from RS-fMRI has found connectivity abnormalities in a number of networks; these include the epileptogenic, cognitive and sensory processing networks. However, in majority of these studies, the effect of epileptic transients in the connectivity of networks has been neglected. EEG–fMRI has frequently shown networks related to epileptic transients that in many cases are concordant with the abnormalities shown in RS studies. This points toward a relevant role of epileptic transients in the network abnormalities detected in RS-fMRI studies. In this review, we summarize the network abnormalities reported by these two techniques side by side, provide evidence of their overlapping findings, and discuss their significance in the context of the methodology of each technique. A number of clinically relevant factors that have been associated with connectivity changes are in turn associated with changes in the frequency of epileptic transients. These factors include different aspects of epilepsy ranging from treatment effects, cognitive processes, or transition between different alertness states (i.e., awake–sleep transition). For RS-fMRI to become a more effective tool to investigate clinically relevant aspects of epilepsy it is necessary to understand connectivity changes associated with epileptic transients, those associated with other clinically relevant factors and the interaction between them, which represents a gap in the current literature. We propose a framework for the investigation of network connectivity in patients with epilepsy that can integrate epileptic processes that occur across different time scales such as epileptic transients and disease duration and the implications of this approach are discussed.

Lateralization and localization of epilepsy related hemodynamic foci using presurgical fMRI

OBJECTIVE

The aim was to develop a method for the purpose of localizing epilepsy related hemodynamic foci for patients suffering intractable focal epilepsy using task-free fMRI alone.

METHODS

We studied three groups of subjects: patients with intractable focal epilepsy, healthy volunteers performing motor tasks, and healthy volunteers in resting state. We performed spatial independent component analysis (ICA) on the fMRI alone data and developed a set of IC selection criteria to identify epilepsy related ICs. The method was then tested in the two healthy groups.

RESULTS

In seven out of the nine surgery patients, identified ICs were concordant with surgical resection. Our results were also consistent with presurgical evaluation of the remaining one patient without surgery and may explain why she was not suitable for resection treatment. In the motor task study of ten healthy subjects, our method revealed components with concordant spatial and temporal features as expected from the unilateral motor tasks. In the resting state study of seven healthy subjects, the method successfully rejected all components in four out of seven subjects as non-epilepsy related components.

CONCLUSION

These results suggest the lateralization and localization value of fMRI alone in presurgical evaluation for patients with intractable unilateral focal epilepsy.

SIGNIFICANCE

The proposed method is noninvasive in nature and easy to implement. It has the potential to be incorporated in current presurgical workup for treating intractable focal epilepsy patients.

The role of functional neuroimaging in pre-surgical epilepsy evaluation

The prevalence of epilepsy is about 1% and one-third of cases do not respond to medical treatment. In an eligible subset of patients with drug-resistant epilepsy, surgical resection of the epileptogenic zone is the only treatment that can possibly cure the disease. Non-invasive techniques provide information for the localization of the epileptic focus in the majority of cases, whereas in others invasive procedures are required. In the last years, non-invasive neuroimaging techniques, such as simultaneous recording of functional magnetic resonance imaging and electroencephalogram (EEG-fMRI), positron emission tomography (PET), single photon emission computed tomography (SPECT), electric and magnetic source imaging (MSI, ESI), spectroscopy (MRS), have proved their usefulness in defining the epileptic focus. The combination of these functional techniques can yield complementary information and their concordance is crucial for guiding clinical decision, namely the planning of invasive EEG recordings or respective surgery. The aim of this review is to present these non-invasive neuroimaging techniques, their potential combination, and their role in the pre-surgical evaluation of patients with pharmaco-resistant epilepsy.

The role of functional neuroimaging in pre-surgical epilepsy evaluation

The prevalence of epilepsy is about 1% and one-third of cases do not respond to medical treatment. In an eligible subset of patients with drug-resistant epilepsy, surgical resection of the epileptogenic zone is the only treatment that can possibly cure the disease. Non-invasive techniques provide information for the localization of the epileptic focus in the majority of cases, whereas in others invasive procedures are required. In the last years, non-invasive neuroimaging techniques, such as simultaneous recording of functional magnetic resonance imaging and electroencephalogram (EEG-fMRI), positron emission tomography (PET), single photon emission computed tomography (SPECT), electric and magnetic source imaging (MSI, ESI), spectroscopy (MRS), have proved their usefulness in defining the epileptic focus. The combination of these functional techniques can yield complementary information and their concordance is crucial for guiding clinical decision, namely the planning of invasive EEG recordings or respective surgery. The aim of this review is to present these non-invasive neuroimaging techniques, their potential combination, and their role in the pre-surgical evaluation of patients with pharmaco-resistant epilepsy.

ICA extracts epileptic sources from fMRI in EEG-negative patients: a retrospective validation study

Simultaneous EEG-fMRI has proven to be useful in localizing interictal epileptic activity. However, the applicability of traditional GLM-based analysis is limited as interictal spikes are often not seen on the EEG inside the scanner. Therefore, we aim at extracting epileptic activity purely from the fMRI time series using independent component analysis (ICA). To our knowledge, we show for the first time that ICA can find sources related to epileptic activity in patients where no interictal spikes were recorded in the EEG. The epileptic components were identified retrospectively based on the known localization of the ictal onset zone (IOZ). We demonstrate that the selected components truly correspond to epileptic activity, as sources extracted from patients resemble significantly better the IOZ than sources found in healthy controls. Furthermore, we show that the epileptic components in patients with and without spikes recorded inside the scanner resemble the IOZ in the same degree. We conclude that ICA of fMRI has the potential to extend the applicability of EEG-fMRI for presurgical evaluation in epilepsy.

Epileptic discharges specifically affect intrinsic connectivity networks during absence seizures

Intrinsic connectivity network (ICN) technique provides a feasible way for evaluating cognitive impairments in epilepsy. This EEG-fMRI study aims to comprehensively assess the alterations of ICNs affected by generalized spike-and-wave discharge (GSWD) during absence seizure (AS). Twelve fMRI sessions with GSWD, and individually paired non-GSWD sessions were acquired from 16 patients with AS. Ten ICNs corresponding to seizure origination and cognitive processes were extracted using independent component analysis. Intra- and inter-network connectivity alterations of the ICNs were observed through comparisons between GSWD and non-GSWD sessions. Sequential correlation analysis between GSWD and the ICN time courses addressed the immediate effects of GSWD on ICNs during AS. GSWD-related increase of intra-network connectivity was found only in the thalamus, and extensive decreases were found in the ICNs corresponding to higher-order cognitive processes including the default-mode network, dorsal attention network, central executive network and salience network. The perceptive networks and motor network were less affected by GSWD. Sequential correlation analysis further demonstrated different responses of the ICNs to GSWD. In addition to GSWD-related functional excitation in the thalamus and functional suspension in the default-mode network, this study revealed extensive inhibitions in the other ICNs corresponding to higher-order cognitive processes, and spared perceptive and motor processes in AS. GSWD elevated synchronization of brain network activity and sequentially affected the ICNs.

An automated method for identifying artifact in independent component analysis of resting-state FMRI

An enduring issue with data-driven analysis and filtering methods is the interpretation of results. To assist, we present an automatic method for identification of artifact in independent components (ICs) derived from functional MRI (fMRI). The method was designed with the following features: does not require temporal information about an fMRI paradigm; does not require the user to train the algorithm; requires only the fMRI images (additional acquisition of anatomical imaging not required); is able to identify a high proportion of artifact-related ICs without removing components that are likely to be of neuronal origin; can be applied to resting-state fMRI; is automated, requiring minimal or no human intervention. We applied the method to a MELODIC probabilistic ICA of resting-state functional connectivity data acquired in 50 healthy control subjects, and compared the results to a blinded expert manual classification. The method identified between 26 and 72% of the components as artifact (mean 55%). About 0.3% of components identified as artifact were discordant with the manual classification; retrospective examination of these ICs suggested the automated method had correctly identified these as artifact. We have developed an effective automated method which removes a substantial number of unwanted noisy components in ICA analyses of resting-state fMRI data. Source code of our implementation of the method is available.

Group independent component analysis and functional MRI examination of changes in language areas associated with brain tumors at different locations

Object

This study investigates the effect of tumor location on alterations of language network by brain tumors at different locations using blood oxygenation level dependent (BOLD) fMRI and group independent component analysis (ICA).

Subjects and Methods

BOLD fMRI data were obtained from 43 right handed brain tumor patients. Presurgical mapping of language areas was performed on all 43 patients with a picture naming task. All data were retrospectively analyzed using group ICA. Patents were divided into three groups based on tumor locations, i.e., left frontal region, left temporal region or right hemisphere. Laterality index (LI) was used to assess language lateralization in each group.

Results

The results from BOLD fMRI and ICA revealed the different language activation patterns in patients with brain tumors located in different brain regions. Language areas, such as Broca’s and Wernicke’s areas, were intact in patients with tumors in the right hemisphere. Significant functional changes were observed in patients with tumor in the left frontal and temporal areas. More specifically, the tumors in the left frontal region affect both Broca’s and Wernicke’s areas, while tumors in the left temporal lobe affect mainly Wernicke’s area. The compensated activation increase was observed in the right frontal areas in patients with left hemisphere tumors.

Conclusion

Group ICA provides a model free alternative approach for mapping functional networks in brain tumor patients. Altered language activation by different tumor locations suggested reorganization of language functions in brain tumor patients and may help better understanding of the language plasticity.

Mapping brain activity using event-related independent components analysis (eICA): specific advantages for EEG-fMRI

Event-related analyses of functional MRI (fMRI) typically assume that the onset and offset of neuronal activity match stimuli onset and offset, and that evoked fMRI signal changes follow the canonical haemodynamic response function (HRF). Some event types, however, may be unsuited to this approach: brief stimuli might elicit an extended neuronal response; anticipatory effects might result in activity preceding the event; or altered neurovascular coupling may result in a non-canonical HRF. An example is interictal epileptiform discharges (IEDs), which may show a non-canonical HRF and fMRI signal changes preceding their onset as detected on EEG. In such cases, less constrained analyses - capable of detecting early, non-canonical responses - may be necessary. A consequence of less constrained analyses, however, is that artefactual sources of signal change - motion or physiological noise for example - may also be detected and mixed with the neuronally-generated signals. In this paper, to address this issue, we describe an event-related independent components analysis (eICA) that identifies different sources of event-related signal change that can then be separately assessed to identify likely artefacts and separate primary from propagated activity. We also describe a group analysis that identifies eICA components that are spatially and temporally consistent across subjects and provides an objective approach for selecting group-specific components likely to be of neural origin. We apply eICA to patients with rolandic epilepsy - with stereotypical IEDs arising from a focus in the rolandic fissure - and demonstrate that a single event-related component, concordant with this source location, is detected.

Mapping interictal epileptic discharges using mutual information between concurrent EEG and fMRI

OBJECTIVE

The mapping of haemodynamic changes related to interictal epileptic discharges (IED) in simultaneous electroencephalography (EEG) and functional MRI (fMRI) studies is usually carried out by means of EEG-correlated fMRI analyses where the EEG information specifies the model to test on the fMRI signal. The sensitivity and specificity critically depend on the accuracy of EEG detection and the validity of the haemodynamic model. In this study we investigated whether an information theoretic analysis based on the mutual information (MI) between the presence of epileptic activity on EEG and the fMRI data can provide further insights into the haemodynamic changes related to interictal epileptic activity. The important features of MI are that: 1) both recording modalities are treated symmetrically; 2) no requirement for a-priori models for the haemodynamic response function, or assumption of a linear relationship between the spiking activity and BOLD responses, and 3) no parametric model for the type of noise or its probability distribution is necessary for the computation of MI.

METHODS

Fourteen patients with pharmaco-resistant focal epilepsy underwent EEG-fMRI and intracranial EEG and/or surgical resection with positive postoperative outcome (seizure freedom or considerable reduction in seizure frequency) was available in 7/14 patients. We used nonparametric statistical assessment of the MI maps based on a four-dimensional wavelet packet resampling method. The results of MI were compared to the statistical parametric maps obtained with two conventional General Linear Model (GLM) analyses based on the informed basis set (canonical HRF and its temporal and dispersion derivatives) and the Finite Impulse Response (FIR) models.

RESULTS

The MI results were concordant with the electro-clinically or surgically defined epileptogenic area in 8/14 patients and showed the same degree of concordance as the results obtained with the GLM-based methods in 12 patients (7 concordant and 5 discordant). In one patient, the information theoretic analysis improved the delineation of the irritative zone compared with the GLM-based methods.

DISCUSSION

Our findings suggest that an information theoretic analysis can provide clinically relevant information about the BOLD signal changes associated with the generation and propagation of interictal epileptic discharges. The concordance between the MI, GLM and FIR maps support the validity of the assumptions adopted in GLM-based analyses of interictal epileptic activity with EEG-fMRI in such a manner that they do not significantly constrain the localization of the epileptogenic zone.

Spontaneous EEG-Functional MRI in Mesial Temporal Lobe Epilepsy: Implications for the Neural Correlates of Consciousness

The combination of electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) has been shown to have great potential for providing a greater understanding of normal and diseased states in both human and animal studies. Simultaneous EEG-fMRI is particularly well suited for the study of epilepsy in that it may reveal the neurobiology of ictal and interictal epileptiform discharges and noninvasively localize epileptogenic foci. Spontaneous, coherent fluctuations of neuronal activity and the coupled hemodynamic responses have also been shown to provide diagnostic markers of disease, extending our understanding of intrinsically structured ongoing brain activity. Following a short summary of the hardware and software development of simultaneous EEG-fMRI, this paper reviews a unified framework of integrating neuronal and hemodynamic processes during epileptic seizures and discusses the role and impact of spontaneous activity in the mesial temporal lobe epilepsies with particular emphasis on the neural and physiological correlates of consciousness.

The impact of EEG/MEG signal processing and modeling in the diagnostic and management of epilepsy

This overview covers recent advances in the field of EEG/MEG signal processing and modeling in epilepsy regarding both interictal and ictal phenomena. In the first part, the main methods used in the analysis of interictal EEG/MEG epileptiform spikes are presented and discussed. Source and volume conductor models are passed in review, namely the equivalent dipole source concept, the requirements for adequate time and spatial sampling, the question of how to validate source solutions, particularly by comparing solutions obtained using scalp and intracranial EEG signals, EEG & MEG data, or EEG simultaneously recorded with fMRI (BOLD signals). In the second part, methods used for the characterization of seizures are considered, namely dipolar modeling of spikes at seizure onset, decomposition of seizure EEG signals into sets of orthogonal spatio-temporal components, and also methods (linear and nonlinear) of estimating seizure propagation. In the third part, the crucial issue of how the transition between interictal and seizure activity takes place is examined. In particular the vicissitudes of the efforts along the road to seizure prediction are shortly reviewed. It is argued that this question can be reduced to the problem of estimating the excitability state of neuronal populations in the course of time as a seizure approaches. The value of active probing methods in contrast with passive analytical methods is emphasized. In the fourth part modeling aspects are considered in the light of two special kinds of epilepsies, absences characterized by spike-and-wave discharges and mesial temporal lobe epilepsy. These two types correspond to different scenarios regarding the transition to epileptic seizures, namely the former is a case of a jump transition and the latter is a typical case of gradual transition. In conclusion, the necessity of developing comprehensive computational models of epileptic seizures is emphasized.

A multimodal imaging approach to the evaluation of post-traumatic epilepsy

Object

Electroencephalography-functional magnetic resonance imaging (EEG-fMRI) coregistration and high-density EEG (hdEEG) can be combined to map noninvasively abnormal brain activation elicited by epileptic processes. By combining noninvasive imaging techniques in a multimodal approach, we sought to investigate pathophysiological mechanisms underlying epileptic activity in seven patients with severe traumatic brain injury.

Materials and methods

Standard EEG and fMRI data were acquired during a single scanning session. The EEG-fMRI data were analyzed using the general linear model and independent component analysis. Source localization of interictal epileptiform discharges (IEDs) was performed using 256-channel hdEEG. Blood oxygenation level dependent (BOLD) localizations were then compared to EEG source reconstruction.

Results

On hdEEG, focal source localization was detected in all seven patients; in six out of seven it was concordant with the expected epileptic activity as defined by EEG data and clinical evaluation; and in four out of seven in whom IEDs were recorded, BOLD signal changes were observed. These activities were partially concordant with the source localization.

Conclusion

Multimodal integration of EEG-fMRI and hdEEG combining two different methods to localize the same epileptic foci appears to be a promising tool to noninvasively map abnormal brain activation in patients with post-traumatic brain injury.

Mapping hemodynamic correlates of seizures using fMRI: A review

Functional magnetic resonance imaging (fMRI) is able to detect changes in blood oxygenation level associated with neuronal activity throughout the brain. For more than a decade, fMRI alone or in combination with simultaneous EEG recording (EEG-fMRI) has been used to investigate the hemodynamic changes associated with interictal and ictal epileptic discharges. This is the first literature review to focus on the various fMRI acquisition and data analysis methods applied to map epileptic seizure-related hemodynamic changes from the first report of an fMRI scan of a seizure to the present day. Two types of data analysis approaches, based on temporal correlation and data driven, are explained and contrasted. The spatial and temporal relationship between the observed hemodynamic changes using fMRI and other non-invasive and invasive electrophysiological and imaging data is considered. We then describe the role of fMRI in localizing and exploring the networks involved in spontaneous and triggered seizure onset and propagation. We also discuss that fMRI alone and combined with EEG hold great promise in the investigation of seizure-related hemodynamic changes non-invasively in humans. We think that this will lead to significant improvements in our understanding of seizures with important consequences for the treatment of epilepsy.

With or without spikes: localization of focal epileptic activity by simultaneous electroencephalography and functional magnetic resonance imaging

In patients with medically refractory focal epilepsy who are candidates for epilepsy surgery, concordant non-invasive neuroimaging data are useful to guide invasive electroencephalographic recordings or surgical resection. Simultaneous electroencephalography and functional magnetic resonance imaging recordings can reveal regions of haemodynamic fluctuations related to epileptic activity and help localize its generators. However, many of these studies (40-70%) remain inconclusive, principally due to the absence of interictal epileptiform discharges during simultaneous recordings, or lack of haemodynamic changes correlated to interictal epileptiform discharges. We investigated whether the presence of epilepsy-specific voltage maps on scalp electroencephalography correlated with haemodynamic changes and could help localize the epileptic focus. In 23 patients with focal epilepsy, we built epilepsy-specific electroencephalographic voltage maps using averaged interictal epileptiform discharges recorded during long-term clinical monitoring outside the scanner and computed the correlation of this map with the electroencephalographic recordings in the scanner for each time frame. The time course of this correlation coefficient was used as a regressor for functional magnetic resonance imaging analysis to map haemodynamic changes related to these epilepsy-specific maps (topography-related haemodynamic changes). The method was first validated in five patients with significant haemodynamic changes correlated to interictal epileptiform discharges on conventional analysis. We then applied the method to 18 patients who had inconclusive simultaneous electroencephalography and functional magnetic resonance imaging studies due to the absence of interictal epileptiform discharges or absence of significant correlated haemodynamic changes. The concordance of the results with subsequent intracranial electroencephalography and/or resection area in patients who were seizure free after surgery was assessed. In the validation group, haemodynamic changes correlated to voltage maps were similar to those obtained with conventional analysis in 5/5 patients. In 14/18 patients (78%) with previously inconclusive studies, scalp maps related to epileptic activity had haemodynamic correlates even when no interictal epileptiform discharges were detected during simultaneous recordings. Haemodynamic changes correlated to voltage maps were spatially concordant with intracranial electroencephalography or with the resection area. We found better concordance in patients with lateral temporal and extratemporal neocortical epilepsy compared to medial/polar temporal lobe epilepsy, probably due to the fact that electroencephalographic voltage maps specific to lateral temporal and extratemporal epileptic activity are more dissimilar to maps of physiological activity. Our approach significantly increases the yield of simultaneous electroencephalography and functional magnetic resonance imaging to localize the epileptic focus non-invasively, allowing better targeting for surgical resection or implantation of intracranial electrode arrays.

Independent component analysis (ICA) of generalized spike wave discharges in fMRI: comparison with general linear model-based EEG-fMRI

Most EEG-fMRI studies in epileptic patients are analyzed using the general linear model (GLM), which assumes a known hemodynamic response function (HRF) to epileptic spikes. In contrast, independent component analysis (ICA) can extract blood-oxygenation level dependent (BOLD) responses without imposing constraints on the HRF. ICA might therefore detect responses that vary in time and shape, and that are not detected in the GLM analysis. In this study, we compared the findings of ICA and GLM analyses in 12 patients with idiopathic generalized epilepsy. Spatial ICA was used to extract independent components from the functional magnetic resonance imaging (fMRI) data. A deconvolution method identified component time courses significantly related to the generalized EEG discharges, without constraining the shape of the HRF. The results from the ICA analysis were compared to those from the GLM analysis. GLM maps and ICA maps showed significant correlation and revealed BOLD responses in the thalamus, caudate nucleus, and default mode areas. In patients with a low rate of discharges per minute, the GLM analysis detected BOLD signal changes within the thalamus and the caudate nucleus that were not revealed by the ICA. In conclusion, ICA is a viable alternative technique to GLM analyses in EEG-fMRI studies related to generalized discharges. This study demonstrated that the BOLD response largely resembles the standard HRF and that GLM analysis is adequate. However, ICA is more dependent on a sufficient number of events than GLM analysis.

Imaging foci of epileptic discharges from simultaneous EEG and fMRI using the canonical HRF

PURPOSE

Simultaneous electroencephalography and functional magnetic resonance imaging (EEG-fMRI) is considered as a powerful and non-invasive method that allows definition of the irritative zone. However, the complex interictal epileptic discharge (IED) may be present in some patients, and sometimes no active foci can be localized using General Linear Model (GLM) which is a widely adopted tool in EEG-fMRI study. The purpose of this study is to develop a new scheme to improve the detectability and localize the canonical HRF localizable foci.

METHOD

Various IEDs are classified using a combination of an independent component analysis (ICA) and a temporal correlation analysis between the independent components and the raw EEG channel; and the classified IEDs are then separately used for foci localization. This scheme is tested by ten patients with variable IEDs, including two patients whose activity could not be identified by common method.

RESULT

Applying this scheme to the two patients, some foci consistent with electroclinical data were localized. When it was applied to the remaining eight patients with positive results using common method, 2-4 types of IEDs were classified, and the activity could be identified from at least one type of IED. The results were similar to that received from common method.

CONCLUSION

These results indicate that the proposed scheme could enhance the imaging of the localizable foci by isolating its IEDs. This scheme is potentially a useful tool for epilepsy clinic.