Detection of epileptic activity in fMRI without recording the EEG

Resting-State Functional Magnetic Resonance Imaging Network Association With Mortality, Epilepsy, Cognition, and Motor Two-Year Outcomes in Suspected Severe Neonatal Acute Brain Injury

BACKGROUND AND OBJECTIVES

In acute brain injury of neonates, resting-state functional magnetic resonance imaging (MRI) (RS) showed incremental association with consciousness, mortality, cognitive and motor development, and epilepsy, with correction for multiple comparisons, at six months postgestation in neonates with suspected acute brain injury (ABI). However, there are relatively few developmental milestones at six months to benchmark against, thus, we extended this cohort study to evaluate two-year outcomes.

METHODS

In 40 consecutive neonates with ABI and RS, ordinal scores of resting-state networks; MRI, magnetic resonance spectroscopy, and electroencephalography; and up to 42-month outcomes of mortality, general and motor development, Pediatric Cerebral Performance Category Scale (PCPC), and epilepsy informed associations between tests and outcomes.

RESULTS

Mean gestational age was 37.8 weeks, 68% were male, and 60% had hypoxic-ischemic encephalopathy. Three died in-hospital, four at six to 42 months, and five were lost to follow-up. Associations included basal ganglia network with PCPC (P = 0.0003), all-mortality (P = 0.005), and motor (P = 0.0004); language/frontoparietal network with developmental delay (P = 0.009), PCPC (P = 0.006), and all-mortality (P = 0.01); default mode network with developmental delay (P = 0.003), PCPC (P = 0.004), neonatal intensive care unit mortality (P = 0.01), and motor (P = 0.009); RS seizure onset zone with epilepsy (P = 0.01); and anatomic MRI with epilepsy (P = 0.01).

CONCLUSION

For the first time, at any age, resting state functional MRI in ABI is associated with long-term epilepsy and RSNs predicted mortality in neonates. Severity of RSN abnormality was associated with incrementally worsened neurodevelopment including cognition, language, and motor function over two years.

The expert's knowledge combined with AI outperforms AI alone in seizure onset zone localization using resting state fMRI

We evaluated whether integration of expert guidance on seizure onset zone (SOZ) identification from resting state functional MRI (rs-fMRI) connectomics combined with deep learning (DL) techniques enhances the SOZ delineation in patients with refractory epilepsy (RE), compared to utilizing DL alone. Rs-fMRI was collected from 52 children with RE who had subsequently undergone ic-EEG and then, if indicated, surgery for seizure control (n = 25). The resting state functional connectomics data were previously independently classified by two expert epileptologists, as indicative of measurement noise, typical resting state network connectivity, or SOZ. An expert knowledge integrated deep network was trained on functional connectomics data to identify SOZ. Expert knowledge integrated with DL showed a SOZ localization accuracy of 84.8 ± 4.5% and F1 score, harmonic mean of positive predictive value and sensitivity, of 91.7 ± 2.6%. Conversely, a DL only model yielded an accuracy of <50% (F1 score 63%). Activations that initiate in gray matter, extend through white matter, and end in vascular regions are seen as the most discriminative expert-identified SOZ characteristics. Integration of expert knowledge of functional connectomics can not only enhance the performance of DL in localizing SOZ in RE but also lead toward potentially useful explanations of prevalent co-activation patterns in SOZ. RE with surgical outcomes and preoperative rs-fMRI studies can yield expert knowledge most salient for SOZ identification.

Whole-brain multivariate hemodynamic deconvolution for functional MRI with stability selection

Conventionally, analysis of functional MRI (fMRI) data relies on available information about the experimental paradigm to establish hypothesized models of brain activity. However, this information can be inaccurate, incomplete or unavailable in multiple scenarios such as resting-state, naturalistic paradigms or clinical conditions. In these cases, blind estimates of neuronal-related activity can be obtained with paradigm-free analysis methods such as hemodynamic deconvolution. Yet, current formulations of the hemodynamic deconvolution problem have three important limitations: (1) their efficacy strongly depends on the appropriate selection of regularization parameters, (2) being univariate, they do not take advantage of the information present across the brain, and (3) they do not provide any measure of statistical certainty associated with each detected event. Here we propose a novel approach that addresses all these limitations. Specifically, we introduce multivariate sparse paradigm free mapping (Mv-SPFM), a novel hemodynamic deconvolution algorithm that operates at the whole brain level and adds spatial information via a mixed-norm regularization term over all voxels. Additionally, Mv-SPFM employs a stability selection procedure that removes the need to select regularization parameters and also lets us obtain an estimate of the true probability of having a neuronal-related BOLD event at each voxel and time-point based on the area under the curve (AUC) of the stability paths. Besides, we present a formulation tailored for multi-echo fMRI acquisitions (MvME-SPFM), which allows us to better isolate fluctuations of BOLD origin on the basis of their linear dependence with the echo time (TE) and to assign physiologically interpretable units (i.e., changes in the apparent transverse relaxation ΔR2∗) to the resulting deconvolved events. Remarkably, we demonstrate that Mv-SPFM achieves comparable performance even when using a single-echo formulation. We demonstrate that this algorithm outperforms existing state-of-the-art deconvolution approaches, and shows higher spatial and temporal agreement with the activation maps and BOLD signals obtained with a standard model-based linear regression approach, even at the level of individual neuronal events. Furthermore, we show that by employing stability selection, the performance of the algorithm depends less on the selection of temporal and spatial regularization parameters λ and ρ. Consequently, the proposed algorithm provides more reliable estimates of neuronal-related activity, here in terms of ΔR2∗, for the study of the dynamics of brain activity when no information about the timings of the BOLD events is available. This algorithm will be made publicly available as part of the splora Python package.

The confound of hemodynamic response function variability in human resting-state functional MRI studies

Functional magnetic resonance imaging (fMRI) is an indirect measure of neural activity with the hemodynamic response function (HRF) coupling it with unmeasured neural activity. The HRF, modulated by several non-neural factors, is variable across brain regions, individuals and populations. Yet, a majority of human resting-state fMRI connectivity studies continue to assume a non-variable HRF. In this article, with supportive prior evidence, we argue that HRF variability cannot be ignored as it substantially confounds within-subject connectivity estimates and between-subjects connectivity group differences. We also discuss its clinical relevance with connectivity impairments confounded by HRF aberrations in several disorders. We present limited data on HRF differences between women and men, which resulted in a 15.4% median error in functional connectivity estimates in a group-level comparison. We also discuss the implications of HRF variability for fMRI studies in the spinal cord. There is a need for more dialogue within the community on the HRF confound, and we hope that our article is a catalyst in the process.

Automated seizure onset zone locator from resting-state functional MRI in drug-resistant epilepsy

OBJECTIVE

Accurate localization of a seizure onset zone (SOZ) from independent components (IC) of resting-state functional magnetic resonance imaging (rs-fMRI) improves surgical outcomes in children with drug-resistant epilepsy (DRE). Automated IC sorting has limited success in identifying SOZ localizing ICs in adult normal rs-fMRI or uncategorized epilepsy. Children face unique challenges due to the developing brain and its associated surgical risks. This study proposes a novel SOZ localization algorithm (EPIK) for children with DRE.

METHODS

EPIK is developed in a phased approach, where fMRI noise-related biomarkers are used through high-fidelity image processing techniques to eliminate noise ICs. Then, the SOZ markers are used through a maximum likelihood-based classifier to determine SOZ localizing ICs. The performance of EPIK was evaluated on a unique pediatric DRE dataset (n = 52). A total of 24 children underwent surgical resection or ablation of an rs-fMRI identified SOZ, concurrently evaluated with an EEG and anatomical MRI. Two state-of-art techniques were used for comparison: (a) least squares support-vector machine and (b) convolutional neural networks. The performance was benchmarked against expert IC sorting and Engel outcomes for surgical SOZ resection or ablation. The analysis was stratified across age and sex.

RESULTS

EPIK outperformed state-of-art techniques for SOZ localizing IC identification with a mean accuracy of 84.7% (4% higher), a precision of 74.1% (22% higher), a specificity of 81.9% (3.2% higher), and a sensitivity of 88.6% (16.5% higher). EPIK showed consistent performance across age and sex with the best performance in those < 5 years of age. It helped achieve a ~5-fold reduction in the number of ICs to be potentially analyzed during pre-surgical screening.

SIGNIFICANCE

Automated SOZ localization from rs-fMRI, validated against surgical outcomes, indicates the potential for clinical feasibility. It eliminates the need for expert sorting, outperforms prior automated methods, and is consistent across age and sex.

Association of network connectivity via resting state functional MRI with consciousness, mortality, and outcomes in neonatal acute brain injury

•

Basal ganglia and seizure onset zone networks were associated with motor outcomes.

•

Broad language/cognitive region networks were associated with developmental delay.

•

Discharge with mortality was linked to default mode and language/cognitive networks.

•

Exams were not linked to networks after multiple testing corrections.

•

Lack of detection of all studied networks only occurred in those who did not survive.

Background

An accurate and comprehensive test of integrated brain network function is needed for neonates during the acute brain injury period to inform on morbidity. This retrospective cohort study assessed whether integrated brain network function acquired by resting state functional MRI during the acute period in neonates with brain injury, is associated with acute exam, neonatal mortality, and 6-month outcomes.

Methods

Study subjects included 40 consecutive neonates with resting state functional MRI acquired within 31 days after suspected brain insult from March 2018 to July 2019 at Phoenix Children’s Hospital. Acute-period exam and test results were assigned ordinal scores based on severity as documented by respective treating specialists. Analyses (Fisher exact, Wilcoxon-rank sum test, ordinal/multinomial logistic regression) examined association of resting state networks with demographics, presentation, neurological exam, electroencephalogram, anatomical MRI, magnetic resonance spectroscopy, passive task functional MRI, and outcomes of discharge condition, outpatient development, motor tone, seizure, and mortality.

Results

Subjects had a mean (standard deviation) gestational age of 37.8 (2.6) weeks, a majority were male (63%), with a diagnosis of hypoxic ischemic encephalopathy (68%). Findings at birth included mild distress (48%), moderately abnormal neurological exam (33%), and consciousness characterized as awake but irritable (40%). Significant associations after multiple testing corrections were detected for resting state networks: basal ganglia with outpatient developmental delay (odds ratio [OR], 14.5; 99.4% confidence interval [CI], 2.00–105; P < .001) and motor tone/weakness (OR, 9.98; 99.4% CI, 1.72–57.9; P < .001); language/frontoparietal network with discharge condition (OR, 5.13; 99.4% CI, 1.22–21.5; P = .002) and outpatient developmental delay (OR, 4.77; 99.4% CI, 1.21–18.7; P=.002); default mode network with discharge condition (OR, 3.72; 99.4% CI, 1.01–13.78; P=.006) and neurological exam (P = .002 (FE); OR, 11.8; 99.4% CI, 0.73–191; P = .01 (OLR)); and seizure onset zone with motor tone/weakness (OR, 3.31; 99.4% CI, 1.08–10.1; P=.003). Resting state networks were not detected in three neonates, who died prior to discharge.

Conclusions

This study provides level 3 evidence (OCEBM Levels of Evidence Working Group) demonstrating that in neonatal acute brain injury, the degree of abnormality of resting state networks is associated with acute exam and outcomes. Total lack of brain network detection was only found in patients who did not survive.

Hemodynamic matrix factorization for functional magnetic resonance imaging

The General Linear Model (GLM) used in task-fMRI relates activated brain areas to extrinsic task conditions. The translation of resulting neural activation into a hemodynamic response is commonly approximated with a linear convolution model using a hemodynamic response function (HRF). There are two major limitations in GLM analysis. Firstly, the GLM assumes that neural activation is either on or off and matches the exact stimulus duration in the corresponding task timings. Secondly, brain networks observed in resting-state fMRI experiments present also during task experiments, but the GLM approach models these task-unrelated brain activity as noise. A novel kernel matrix factorization approach, called hemodynamic matrix factorization (HMF), is therefore proposed that addresses both limitations by assuming that task-related and task-unrelated brain activity can be modeled with the same convolution model as in GLM analysis. By contrast to the GLM, the proposed HMF is a blind source separation (BSS) technique, which decomposes fMRI data into modes. Each mode comprises of a neural activation time course and a spatial mapping. Two versions of HMF are proposed in which the neural activation time course of each mode is convolved with either the canonical HRF or predetermined subject-specific HRFs. Firstly, HMF with the canonical HRF is applied to two open-source cohorts. These cohorts comprise of several task experiments including motor, incidental memory, spatial coherence discrimination, verbal discrimination task and a very short localization task, engaging multiple parts of the eloquent cortex. HMF modes were obtained whose neural activation time course followed original task timings and whose corresponding spatial map matched cortical areas known to be involved in the respective task processing. Secondly, the alignment of these neural activation time courses to task timings were further improved by replacing the canonical HRF with subject-specific HRFs during HMF mode computation. In addition to task-related modes, HMF also produced seemingly task-unrelated modes whose spatial maps matched known resting-state networks. The validity of a fMRI task experiment relies on the assumption that the exposure to a stimulus for a given time causes an imminent increase in neural activation of equal duration. The proposed HMF is an attempt to falsify this assumption and allows to identify subject task participation that does not comply with the experiment instructions.

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.

Abnormal Profiles of Local Functional Connectivity Proximal to Focal Cortical Dysplasias

Introduction

Focal cortical dysplasia (FCD) is a congenital malformation of cortical development that often leads to medically refractory epilepsy. Focal resection can be an effective treatment, but is challenging as the surgically relevant abnormality may exceed the MR-visible lesion. The aim of the current study is to develop methodology to characterize the profile of functional connectivity around FCDs using resting-state functional MRI and in the individual patient. The detection of aberrant connectivity may provide a means to more completely delineate the clinically relevant lesion.

Materials and Methods

Fifteen FCD patients (age, mean±SD: 31±11 years; 11 males) and 16 matched healthy controls (35±9 years; 7 males) underwent structural and functional imaging at 3 Tesla. The cortical surface was reconstructed from the T1-weighted scan and the registered functional MRI data was spatially normalized to a common anatomical standard space employing the gyral pattern. Seed-based functional connectivity was determined in all subjects for all dysplasia locations. A single patient was excluded based on an aberrant FCD seed time series. Functional connectivity as a function of geodesic distance (along the cortical surface) was compared between the individual patients and the homotopic normative connectivity profiles derived from the controls.

Results

In 12/14 patients, aberrant profiles of functional connectivity were found, which demonstrated both hyper- and hypoconnectivity as well as combinations. Abnormal functional connectivity was typically found (also) beyond the lesion visible on structural MRI, while functional connectivity profiles not related to a lesion appeared normal in patients.

Conclusion

This novel functional MRI technique has potential for delineating functionally aberrant from normal cortex beyond the structural lesion in FCD, which remains to be confirmed in future research.

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.

Automatic epileptic seizure detection using scalp EEG and advanced artificial intelligence techniques

The epilepsies are a heterogeneous group of neurological disorders and syndromes characterised by recurrent, involuntary, paroxysmal seizure activity, which is often associated with a clinicoelectrical correlate on the electroencephalogram. The diagnosis of epilepsy is usually made by a neurologist but can be difficult to be made in the early stages. Supporting paraclinical evidence obtained from magnetic resonance imaging and electroencephalography may enable clinicians to make a diagnosis of epilepsy and investigate treatment earlier. However, electroencephalogram capture and interpretation are time consuming and can be expensive due to the need for trained specialists to perform the interpretation. Automated detection of correlates of seizure activity may be a solution. In this paper, we present a supervised machine learning approach that classifies seizure and nonseizure records using an open dataset containing 342 records. Our results show an improvement on existing studies by as much as 10% in most cases with a sensitivity of 93%, specificity of 94%, and area under the curve of 98% with a 6% global error using a k-class nearest neighbour classifier. We propose that such an approach could have clinical applications in the investigation of patients with suspected seizure disorders.

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.

Study on the Relationships between Intrinsic Functional Connectivity of the Default Mode Network and Transient Epileptic Activity

Rationale: Simultaneous recording of electroencephalogram and functional MRI (EEG–fMRI) is a powerful tool for localizing epileptic networks via the detection of hemodynamic changes correlated with interictal epileptic discharges (IEDs). fMRI can be used to study the long-lasting effect of epileptic activity by assessing stationary functional connectivity during the resting-state period [especially, the connectivity of the default mode network (DMN)]. Temporal lobe epilepsy (TLE) and idiopathic generalized epilepsy (IGE) are associated with low responsiveness and disruption of DMN activity. A dynamic functional connectivity approach might enable us to determine the effect of IEDs on DMN connectivity and to better understand the correlation between DMN connectivity changes and altered consciousness.

Method: We studied dynamic changes in DMN intrinsic connectivity and their relation to IEDs. Six IGE patients (with generalized spike and slow-waves) and 6 TLE patients (with unilateral left temporal spikes) were included. Functional connectivity before, during, and after IEDs was estimated using a sliding window approach and compared with the baseline period.

Results: No dependence on window size was observed. The baseline DMN connectivity was decreased in the left hemisphere (ipsilateral to the epileptic focus) in TLEs and was less strong but remained bilateral in IGEs. We observed an overall increase in DMN intrinsic connectivity prior to the onset of IEDs in both IGEs and TLEs. After IEDs in TLEs, we found that DMN connectivity increased before it returned to baseline values. Most of the DMN regions with increased connectivity before and after IEDs were lateralized to the left hemisphere in TLE (i.e., ipsilateral to the epileptic focus).

Conclusion: Results suggest that DMN connectivity may facilitate IED generation and may be affected at the time of the IED. However, these results need to be confirmed in a larger independent cohort.

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.

[MRI in epilepsy and migraine]

In vivo observation of the ictal events are thought to help understanding the etiology and pathology in both epilepsy and migraine. While simultaneous recording of EEG and fMRI is actively conducted for the former, in some cases, epileptogenic activity is undetectable by EEG. Attempts to detect such abnormal brain activity by using fMRI are underway. Analysis methods for resting-state fMRI can be applicable for such purposes. For migraine, fMRI is also highly valuable in detecting series of ictal events. However, since the disease is suspected to involve abnormal neuro-vascular coupling, it is not always straightforward how to interpret the observation by vascular-dependent methods. Therefore development of non-vascular methods is critical for future advances.

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.

Electroencephalography/functional magnetic resonance imaging responses help predict surgical outcome in focal epilepsy

PURPOSE

Simultaneous electroencephalography/functional magnetic resonance imaging (EEG/fMRI) recording can noninvasively map in the whole brain the hemodynamic response following an interictal epileptic discharge. EEG/fMRI is gaining interest as a presurgical evaluation tool. This study aims to determine how hemodynamic responses related to epileptic activity can help predict surgical outcome in patients considered for epilepsy surgery.

METHODS

Thirty-five consecutive patients with focal epilepsy who had significant hemodynamic responses and eventually surgical resection, were studied. The statistical map of hemodynamic responses were generated and co-registered to postoperative anatomic imaging. Patients were classified into four groups defined by the relative relationship between the location of the maximum hemodynamic response and the resection: group 1, fully concordant; group 2, partially concordant; group 3, partially discordant; and group 4, fully discordant. These findings were correlated with surgical outcome with at least 12-month follow-up.

KEY FINDINGS

Ten patients in group 1 had the maximum t value (t-max) inside the resection; nine in group 2 had the t-max outside but close to the resection and the cluster with t-max overlapped the resection; five in group 3 had the t-max remote from resection, but with another less significant cluster in the resection; and 11 in group 4 had no response in the resection. The degree of concordance correlated largely with surgical outcome: a good surgical outcome (Engel's class I) was found in 7 of 10 patients of group 1, 4 of 9 of group 2, 3 of 5 of group 3, and only 1 of 11 of group 4. These results indicate that the partially concordant and partially discordant groups are best considered as inconclusive. In contrast, in the fully concordant and fully discordant groups, the sensitivity, specificity, positive predictive value, and negative predictive value were high, 87.5%, 76.9%, 70%, and 90.9%, respectively.

SIGNIFICANCE

This study demonstrates that hemodynamic responses related to epileptic activity can help delineate the epileptogenic region. Full concordance between maximum response and surgical resection is indicative of seizure freedom, whereas a resection leaving the maximum response intact is likely to lead to a poor outcome. EEG/fMRI is noninvasive but is limited to patients in whom interictal epileptic discharges can be recorded during the 60-90 min scan.

On the definition of signal-to-noise ratio and contrast-to-noise ratio for FMRI data

Signal-to-noise ratio, the ratio between signal and noise, is a quantity that has been well established for MRI data but is still subject of ongoing debate and confusion when it comes to fMRI data. fMRI data are characterised by small activation fluctuations in a background of noise. Depending on how the signal of interest and the noise are identified, signal-to-noise ratio for fMRI data is reported by using many different definitions. Since each definition comes with a different scale, interpreting and comparing signal-to-noise ratio values for fMRI data can be a very challenging job. In this paper, we provide an overview of existing definitions. Further, the relationship with activation detection power is investigated. Reference tables and conversion formulae are provided to facilitate comparability between fMRI studies.

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.

Resting-state fMRI studies in epilepsy

Epilepsy is a disease characterized by abnormal spontaneous activity in the brain. Resting-state functional magnetic resonance imaging (RS-fMRI) is a powerful technique for exploring this activity. With good spatial and temporal resolution, RS-fMRI is a promising approach for accurate localization of the focus of seizure activity. Although simultaneous electroencephalogram-fMRI has been performed with patients in the resting state, most studies focused on activation. This mini-review focuses on RS-fMRI alone, including its computational methods and its application to epilepsy.