Publications

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.

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.

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.

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.

OBJECTIVE: Intracranial electroencephalographic (icEEG) studies show that interictal ripples propagate across the brain of children with medically refractory epilepsy (MRE), and the onset of this propagation (ripple onset zone [ROZ]) estimates the epileptogenic zone. It is still unknown whether we can map this propagation noninvasively. The goal of this study is to map ripples (ripple zone [RZ]) and their propagation onset (ROZ) using high-density EEG (HD-EEG) and magnetoencephalography (MEG), and to estimate their prognostic value in pediatric epilepsy surgery.

METHODS: We retrospectively analyzed simultaneous HD-EEG and MEG data from 28 children with MRE who underwent icEEG and epilepsy surgery. Using electric and magnetic source imaging, we estimated virtual sensors (VSs) at brain locations that matched the icEEG implantation. We detected ripples on VSs, defined the virtual RZ and virtual ROZ, and estimated their distance from icEEG. We assessed the predictive value of resecting virtual RZ and virtual ROZ for postsurgical outcome. Interictal spike localization on HD-EEG and MEG was also performed and compared with ripples.

RESULTS: We mapped ripple propagation in all patients with HD-EEG and in 27 (96%) patients with MEG. The distance from icEEG did not differ between HD-EEG and MEG when mapping the RZ (26-27mm, p = 0.6) or ROZ (22-24mm, p = 0.4). Resecting the virtual ROZ, but not virtual RZ or the sources of spikes, was associated with good outcome for HD-EEG (p = 0.016) and MEG (p = 0.047).

INTERPRETATION: HD-EEG and MEG can map interictal ripples and their propagation onset (virtual ROZ). Noninvasively mapping the ripple onset may augment epilepsy surgery planning and improve surgical outcome of children with MRE. ANN NEUROL 2021;89:911-925.

A term neonate is born with the ability to suck; this neuronal network is already formed and functional by 28 weeks gestational age and continues to evolve into adulthood. Because of the necessity of acquiring nutrition, the complexity of the neuronal network needed to suck, and neuroplasticity in infancy, the skill of sucking has the unique ability to give insight into areas of the brain that may be damaged either during or before birth. Interpretation of the behaviors during sucking shows promise in guiding therapies and how to potentially repair the damage early in life, when neuroplasticity is high. Sucking requires coordinated suck-swallow-breathe actions and is classified into two basic types, nutritive and non-nutritive. Each type of suck has particular characteristics that can be measured and used to learn about the infant's neuronal circuitry. Basic sucking and swallowing are present in embryos and further develop to incorporate breathing ex utero. Due to the rhythmic nature of the suck-swallow-breathe process, these motor functions are controlled by central pattern generators. The coordination of swallowing, breathing, and sucking is an enormously complex sensorimotor process. Because of this complexity, brain injury before birth can have an effect on these sucking patterns. Clinical assessments allow evaluators to score the oral-motor pattern, however, they remain ultimately subjective. Thus, clinicians are in need of objective measures to identify the specific area of deficit in the sucking pattern of each infant to tailor therapies to their specific needs. Therapeutic approaches involve pacifiers, cheek/chin support, tactile, oral kinesthetic, auditory, vestibular, and/or visual sensorimotor inputs. These therapies are performed to train the infant to suck appropriately using these subjective assessments along with the experience of the therapist (usually a speech therapist), but newer, more objective measures are coming along. Recent studies have correlated pathological sucking patterns with neuroimaging data to get a map of the affected brain regions to better inform therapies. The purpose of this review is to provide a broad scope synopsis of the research field of infant nutritive and non-nutritive feeding, their underlying neurophysiology, and relationship of abnormal activity with brain injury in preterm and term infants.