Publications

OBJECTIVE: Delineation of the seizure onset zone (SOZ) is required in children with drug resistant epilepsy (DRE) undergoing neurosurgery. Intracranial EEG (icEEG) serves as gold standard but has limitations. Here, we examine the utility of virtual implantation with electrical source imaging (ESI) on ictal scalp EEG for mapping the SOZ and predict surgical outcome.

METHODS: We retrospectively analyzed EEG data from 35 children with DRE who underwent surgery and dichotomized into seizure-free (SF) and non-seizure-free (NSF). We estimated virtual sensors (VSs) at brain locations that matched icEEG implantation and compared ictal patterns at VSs vs icEEG. We calculated the agreement between VSs SOZ and clinically defined SOZ and built receiver operating characteristic (ROC) curves to test whether it predicted outcome.

RESULTS: Twenty-one patients were SF after surgery. Moderate agreement between virtual and icEEG patterns was observed (kappa = 0.45, p < 0.001). Virtual SOZ agreement with clinically defined SOZ was higher in SF vs NSF patients (66.6% vs 41.6%, p = 0.01). Anatomical concordance of virtual SOZ with clinically defined SOZ predicted outcome (AUC = 0.73; 95% CI: 0.57-0.89; sensitivity = 66.7%; specificity = 78.6%; accuracy = 71.4%).

CONCLUSIONS: Virtual implantation on ictal scalp EEG can approximate the SOZ and predict outcome.

SIGNIFICANCE: SOZ mapping with VSs may contribute to tailoring icEEG implantation and predict outcome.

Delineation of resected brain cavities on magnetic resonance images (MRIs) of epilepsy surgery patients is essential for neuroimaging/neurophysiology studies investigating biomarkers of the epileptogenic zone. The gold standard to delineate the resection on MRI remains manual slice-by-slice tracing by experts. Here, we proposed and validated a semiautomated MRI segmentation pipeline, generating an accurate model of the resection and its anatomical labeling, and developed a graphical user interface (GUI) for user-friendly usage. We retrieved pre- and postoperative MRIs from 35 patients who had focal epilepsy surgery, implemented a region-growing algorithm to delineate the resection on postoperative MRIs and tested its performance while varying different tuning parameters. Similarity between our output and hand-drawn gold standards was evaluated via dice similarity coefficient (DSC; range: 0-1). Additionally, the best segmentation pipeline was trained to provide an automated anatomical report of the resection (based on presurgical brain atlas). We found that the best-performing set of parameters presented DSC of 0.83 (0.72-0.85), high robustness to seed-selection variability and anatomical accuracy of 90% to the clinical postoperative MRI report. We presented a novel user-friendly open-source GUI that implements a semiautomated segmentation pipeline specifically optimized to generate resection models and their anatomical reports from epilepsy surgery patients, while minimizing user interaction.

OBJECTIVE: To assess the utility of interictal magnetic and electric source imaging (MSI and ESI) using dipole clustering in magnetic resonance imaging (MRI)-negative patients with drug resistant epilepsy (DRE).

METHODS: We localized spikes in low-density (LD-EEG) and high-density (HD-EEG) electroencephalography as well as magnetoencephalography (MEG) recordings using dipoles from 11 pediatric patients. We computed each dipole's level of clustering and used it to discriminate between clustered and scattered dipoles. For each dipole, we computed the distance from seizure onset zone (SOZ) and irritative zone (IZ) defined by intracranial EEG. Finally, we assessed whether dipoles proximity to resection was predictive of outcome.

RESULTS: LD-EEG had lower clusterness compared to HD-EEG and MEG (p < 0.05). For all modalities, clustered dipoles showed higher proximity to SOZ and IZ than scattered (p < 0.001). Resection percentage was higher in optimal vs. suboptimal outcome patients (p < 0.001); their proximity to resection was correlated to outcome (p < 0.001). No difference in resection percentage was seen for scattered dipoles between groups.

CONCLUSION: MSI and ESI dipole clustering helps to localize the SOZ and IZ and facilitate the prognostic assessment of MRI-negative patients with DRE.

SIGNIFICANCE: Assessing the MSI and ESI clustering allows recognizing epileptogenic areas whose removal is associated with optimal outcome.

Children with medically refractory epilepsy (MRE) require resective neurosurgery to achieve seizure freedom, whose success depends on accurate delineation of the epileptogenic zone (EZ). Functional connectivity (FC) can assess the extent of epileptic brain networks since intracranial EEG (icEEG) studies have shown its link to the EZ and predictive value for surgical outcome in these patients. Here, we propose a new noninvasive method based on magnetoencephalography (MEG) and high-density (HD-EEG) data that estimates FC metrics at the source level through an "implantation" of virtual sensors (VSs). We analyzed MEG, HD-EEG, and icEEG data from eight children with MRE who underwent surgery having good outcome and performed source localization (beamformer) on noninvasive data to build VSs at the icEEG electrode locations. We analyzed data with and without Interictal Epileptiform Discharges (IEDs) in different frequency bands, and computed the following FC matrices: Amplitude Envelope Correlation (AEC), Correlation (CORR), and Phase Locking Value (PLV). Each matrix was used to generate a graph using Minimum Spanning Tree (MST), and for each node (i.e., each sensor) we computed four centrality measures: betweenness, closeness, degree, and eigenvector. We tested the reliability of VSs measures with respect to icEEG (regarded as benchmark) via linear correlation, and compared FC values inside vs. outside resection. We observed higher FC inside than outside resection (p<0.05) for AEC [alpha (8-12 Hz), beta (12-30 Hz), and broadband (1-50 Hz)] on data with IEDs and AEC theta (4-8 Hz) on data without IEDs for icEEG, AEC broadband (1-50 Hz) on data without IEDs for MEG-VSs, as well as for all centrality measures of icEEG and MEG/HD-EEG-VSs. Additionally, icEEG and VSs metrics presented high correlation (0.6-0.9, p<0.05). Our data support the notion that the proposed method can potentially replicate the icEEG ability to map the epileptogenic network in children with MRE.Clinical Relevance - The estimation of FC with noninvasive techniques, such as MEG and HD-EEG, via VSs is a promising tool that would help the presurgical evaluation by delineating the EZ without waiting for a seizure to occur, and potentially improve the surgical outcome of patients with MRE undergoing surgery.

Studies on intracranial electroencephalography (icEEG) recordings of patients with drug resistant epilepsy (DRE) show that epilepsy biomarkers propagate in time across brain areas. Here, we propose a novel method that estimates critical features of these propagations for different epilepsy biomarkers (spikes, ripples, and fast ripples), and assess their common onset as a reliable biomarker of the epileptogenic zone (EZ). For each biomarker, an automatic algorithm ranked the icEEG electrodes according to their timing occurrence in propagations and finally dichotomized them as onset or spread. We validated our algorithm on icEEG recordings of 8 children with DRE having a good surgical outcome (Engel score = 1). We estimated the overlap of the onset, spread, and entire zone of propagation with the resection (RZ) and the seizure onset zone (SOZ). Spike and ripple propagations were seen in all patients, whereas fast ripple propagations were seen in 6 patients. Spike, ripple, and fast ripple propagations had a mean duration of 28.3 ± 24.3 ms, 38.7 ± 37 ms, and 25 ± 14 ms respectively. Onset electrodes predicted the RZ and SOZ with higher specificity compared to the entire zone for all three biomarkers (p<0.05). Overlap of spike and ripple onsets presented a higher specificity than each separate biomarker onset: for the SOZ, the onsets overlap was more specific (97.89 ± 2.36) than the ripple onset (p=0.031); for the RZ, the onsets overlap was more specific (98.48 ± 1.5) than the spike onset (p=0.016). Yet, the entire zone for spike and ripple propagations predicted the RZ with higher sensitivity compared to each biomarker's onset (or spread) (p<0.05). We present, for the first time, preliminary evidence from icEEG data that fast ripples propagate in time across large areas of the brain. The onsets overlap of spike and ripple propagations constitutes an extremely specific (but not sensitive) biomarker of the EZ, compared to areas of spread (and entire areas) in propagation.

Interictal epileptiform discharges (IEDs) serve as sensitive but not specific biomarkers of epilepsy that can delineate the epileptogenic zone (EZ) in patients with drug resistant epilepsy (DRE) undergoing surgery. Intracranial EEG (icEEG) studies have shown that IEDs propagate in time across large areas of the brain. The onset of this propagation is regarded as a more specific biomarker of epilepsy than areas of spread. Yet, the limited spatial resolution of icEEG does not allow to identify the onset of this activity with high precision. Here, we propose a new method of mapping the spatiotemporal propagation of IEDs (and identify its onset) by using Electrical Source Imaging (ESI) on icEEG bypassing the spatial limitations of icEEG. We validated our method on icEEG recordings from 8 children with DRE who underwent surgery with good outcome (Engel score =1). On each icEEG channel, we detected IEDs and identified the propagation onset using an automated algorithm. We localized the propagation of IEDs with dynamic Statistical Parametric Mapping (dSPM) using a time-sliding window approach. We defined two brain regions: the ESI-onset and ESI-spread zone. We estimated the overlap of these regions with resection volume (in percentage), which served as the gold-standard of the EZ. We also estimated the mean distance of these regions from resection and clinically defined seizure onset zone (SOZ). We observed spatio-temporal propagation of IEDs in all patients across several channels (98 [85-102]) with a mean duration of 155 ms [96-186 ms]. A higher overlap with resection was seen for the ESI-onset zone compared to spread (73.3 % [ 47.4-100 %], 36.5 % [20.3-59.9 %], p = 0.008). The distance of the ESI-onset from resection was shorter compared to the ESI-spread zone (4.3 mm [3.4-5.5 mm], 7.4 mm [6.0-20.6 mm], p = 0.008) and the same trend was observed for the distance from the SOZ (11.9 mm [7.2-15.1 mm], 20.6 mm [15.4-27.2 mm], p = 0.02). These findings show that our method can map the spatiotemporal propagation of IEDs and de-lineate its onset, which is a reliable and focal biomarker of the EZ in children with DRE.Clinical Relevance - ESI on icEEG recordings of children with DRE can localize the spikes propagation phenomenon and help in the delineation of the EZ.

INTRODUCTION: Generalized tonic-clonic seizures (GTCS) are associated with elevated electrodermal activity (EDA) and postictal generalized electroencephalographic suppression (PGES), markers that may indicate sudden unexpected death in epilepsy (SUDEP) risk. This study investigated the association of GTCS semiology, EDA, and PGES in children with epilepsy.

METHODS: Patients admitted to the Boston Children's Hospital long-term video-EEG monitoring unit wore a sensor that records EDA. We selected patients with at least one GTCS and reviewed video-EEGs for semiology, tonic and clonic phase duration, total clinical seizure duration, electrographic onset, offset, and PGES. We grouped patients into three semiology classes: GTCS 1: bilateral symmetric tonic arm extension, GTCS 2: no specific tonic arm extension or flexion, GTCS 3: unilateral or asymmetrical arm extension, tonic arm flexion or posturing that does not fit into GTCS 1 or 2. We analyzed the correlation between semiology, EDA, and PGES, and measured the area under the curve (AUC) of the ictal EDA (seizure onset to one hour after), subtracting baseline EDA (one-hour seizure-free before seizure onset). Using generalized estimating equation (GEE) and linear regression, we analyzed all seizures and single episodes per patient.

RESULTS: We included 30 patients (median age 13.8 ± 3.6 years, 46.7% females) and 53 seizures. With GEE, GTCS 1 was associated with longer PGES duration compared to GTCS 2 (Estimate (β) = -26.32 s, 95% Confidence Interval (CI): -36.46 to -16.18, p < 0.001), and the presence of PGES was associated with greater EDA change (β = 429604 μS, 95% CI: 3550.96 to 855657.04, p = 0.048). With single-episode analysis, GTCS 1 had greater EDA change than GTCS 2 ((β = -601339 μS, 95% CI: -1167016.56 to -35661.44, p = 0.047). EDA increased with PGES presence (β = 637500 μS, 95% CI: 183571.84 to 1091428.16, p = 0.01) and duration (β = 16794 μS, 95% CI: 5729.8 to 27858.2, p = 0.006). Patients with GTCS 1 had longer PGES duration compared to GTCS 2 (β = -30.53 s, 95% CI: -44.6 to -16.46, p < 0.001) and GTCS 3 (β = -22.07 s, 95% CI: -38.95 to -5.19, p = 0.016).

CONCLUSION: In children with epilepsy, PGES correlates with greater ictal EDA. GTCS 1 correlated with longer PGES duration and may indirectly correlate with greater ictal EDA. Our study suggests potential applications in monitoring and preventing SUDEP in these patients.

About 30% of children with drug-resistant epilepsy (DRE) continue to have seizures after epilepsy surgery. Since epilepsy is increasingly conceptualized as a network disorder, understanding how brain regions interact may be critical for planning re-operation in these patients. We aimed to estimate functional brain connectivity using scalp EEG and its evolution over time in patients who had repeated surgery (RS-group, n = 9) and patients who had one successful surgery (seizure-free, SF-group, n = 12). We analyzed EEGs without epileptiform activity at varying time points (before and after each surgery). We estimated functional connectivity between cortical regions and their relative centrality within the network. We compared the pre- and post-surgical centrality of all the non-resected (untouched) regions (far or adjacent to resection) for each group (using the Wilcoxon signed rank test). In alpha, theta, and beta frequency bands, the post-surgical centrality of the untouched cortical regions increased in the SF group (p < 0.001) whereas they decreased (p < 0.05) or did not change (p > 0.05) in the RS group after failed surgeries; when re-operation was successful, the post-surgical centrality of far regions increased (p < 0.05). Our data suggest that removal of the epileptogenic focus in children with DRE leads to a gain in the network centrality of the untouched areas. In contrast, unaltered or decreased connectivity is seen when seizures persist after surgery.

OBJECTIVE: To assess whether ictal electric source imaging (ESI) on low-density scalp EEG can approximate the seizure onset zone (SOZ) location and predict surgical outcome in children with refractory epilepsy undergoing surgery.

METHODS: We examined 35 children with refractory epilepsy. We dichotomized surgical outcome into seizure- and non-seizure-free. We identified ictal onsets recorded with scalp and intracranial EEG and localized them using equivalent current dipoles and standardized low-resolution magnetic tomography (sLORETA). We estimated the localization accuracy of scalp EEG as distance of scalp dipoles from intracranial dipoles. We also calculated the distances of scalp dipoles from resection, as well as their resection percentage and compared between seizure-free and non-seizure-free patients. We built receiver operating characteristic curves to test whether resection percentage predicted outcome.

RESULTS: Resection distance was lower in seizure-free patients for both dipoles (p = 0.006) and sLORETA (p = 0.04). Resection percentage predicted outcome with a sensitivity of 57.1% (95% CI, 34-78.2%), a specificity of 85.7% (95% CI, 57.2-98.2%) and an accuracy of 68.6% (95% CI, 50.7-83.5%) (p = 0.01).

CONCLUSION: Ictal ESI performed on low-density scalp EEG can delineate the SOZ and predict outcome.

SIGNIFICANCE: Such an application may increase the number of children who are referred for epilepsy surgery and improve their outcome.