The long-term impact of cerebellar tumor resection on executive functioning, anxiety, and fear of pain: A mixed methodology pilot study

This pilot study investigated the long-term impact of a surgery-only treatment (no exposure to other treatments, such as chemotherapy and radiation) for pediatric cerebellar low-grade gliomas on executive function, anxiety, and fear of pain (FOP) beliefs. Twelve patients who underwent surgical glioma resection during childhood (surgery age was 4-16 years, study visit age was 10-28 years), and 12 pain-free controls matched for age, sex, race, and handedness were tested. The spatial extent of resection was precisely mapped using magnetic resonance imaging (MRI). Executive function, anxiety, and FOP were assessed using validated self-report age-appropriate questionnaires for children and adults. Structured clinical interviews at a post-surgery follow-up visit were completed (average: 89 months, range: 20-99). No significant differences in FOP (FOPQ-C t[14 = 1.81, p = 0.09; FOPQ-III t[4] = 0.29, p = 0.79), executive function scores (BRIEF t[20] = 0.30, p = 0.28), or anxiety scores (MASC t[16] = 0.19, p = 0.85; MAQ t[4] = 1.80, p = 0.15) were found in pediatric or adult patients compared to pain-free controls. Clinical interviews mainly categorized pediatric patients as not anxious. One participant reported mild/subclinical anxiety, and one had moderate clinical anxiety. Neither psychologists nor patients endorsed impairments to executive functioning, anxiety, or FOP. Our pilot results suggest that pediatric cerebellar tumor survivors treated with surgery-only have favorable long-term functioning related to these themes. While these results are promising, they will need to be replicated in a larger patient sample.

Interictal EEG source connectivity to localize the epileptogenic zone in patients with drug-resistant epilepsy: A machine learning approach

OBJECTIVE

To deconstruct the epileptogenic networks of patients with drug-resistant epilepsy (DRE) using source functional connectivity (FC) analysis; unveil the FC biomarkers of the epileptogenic zone (EZ); and develop machine learning (ML) models to estimate the EZ using brief interictal electroencephalography (EEG) data.

METHODS

We analyzed scalp EEG from 50 patients with DRE who had surgery. We reconstructed the activity (electrical source imaging [ESI]) of virtual sensors (VSs) across the whole cortex and computed FC separately for epileptiform and non-epileptiform EEG epochs (with or without spikes). In patients with good outcome (Engel 1a), four cortical regions were defined: EZ (resection) and three non-epileptogenic zones (NEZs) in the same and opposite hemispheres. Region-specific FC features in six frequency bands and three spatial ranges (long, short, inner) were compared between regions (Wilcoxon sign-rank). We developed ML classifiers to identify the VSs in the EZ using VS-specific FC features. Cross-validation was performed using good outcome data. Performance was compared with poor outcomes and interictal spike localization.

RESULTS

FC differed between EZ and NEZs (p < .05) during non-epileptiform and epileptiform epochs, showing higher FC in the EZ than its homotopic contralateral NEZ. During epileptiform epochs, the NEZ in the epileptogenic hemisphere showed higher FC than its contralateral NEZ. In good outcome patients, the ML classifiers reached 75% accuracy to the resection (91% sensitivity; 74% specificity; distance from EZ: 38 mm) using epileptiform epochs (gamma and beta frequency bands) and 62% accuracy using broadband non-epileptiform epochs, both outperforming spike localization (accuracy = 47%; p < .05; distance from EZ: 57 mm). Lower performance was seen in poor outcomes.

SIGNIFICANCE

We present an FC approach to extract EZ biomarkers from brief EEG data. Increased FC in various frequencies characterized the EZ during epileptiform and non-epileptiform epochs. FC-based ML models identified the resection better in good than poor outcome patients, demonstrating their potential for presurgical use in pediatric DRE.

The Human Affectome

Over the last decades, theoretical perspectives in the interdisciplinary field of the affective sciences have proliferated rather than converged due to differing assumptions about what human affective phenomena are and how they work. These metaphysical and mechanistic assumptions, shaped by academic context and values, have dictated affective constructs and operationalizations. However, an assumption about the purpose of affective phenomena can guide us to a common set of metaphysical and mechanistic assumptions. In this capstone paper, we home in on a nested teleological principle for human affective phenomena in order to synthesize metaphysical and mechanistic assumptions. Under this framework, human affective phenomena can collectively be considered algorithms that either adjust based on the human comfort zone (affective concerns) or monitor those adaptive processes (affective features). This teleologically-grounded framework offers a principled agenda and launchpad for both organizing existing perspectives and generating new ones. Ultimately, we hope the Human Affectome brings us a step closer to not only an integrated understanding of human affective phenomena, but an integrated field for affective research.

Sleep Spindle Generation Before and After Epilepsy Surgery: A Source Imaging Study in Children with Drug-Resistant Epilepsy

INTRODUCTION

Literature lacks studies investigating the cortical generation of sleep spindles in drug-resistant epilepsy (DRE) and how they evolve after resection of the epileptogenic zone (EZ). Here, we examined sleep EEGs of children with focal DRE who became seizure-free after focal epilepsy surgery, and aimed to investigate the changes in the spindle generation before and after the surgery using low-density scalp EEG and electrical source imaging (ESI).

METHODS

We analyzed N2-sleep EEGs from 19 children with DRE before and after surgery. We identified slow (8-12 Hz) and fast spindles (13-16 Hz), computed their spectral features and cortical generators through ESI and computed their distance from the EZ and irritative zone (IZ). We performed two-way ANOVA testing the effect of spindle type (slow vs. fast) and surgical phase (pre-surgery vs. post-surgery) on each feature.

RESULTS

Power, frequency and cortical activation of slow spindles increased after surgery (p < 0.005), while this was not seen for fast spindles. Before surgery, the cortical generators of slow spindles were closer to the EZ (57.3 vs. 66.2 mm, p = 0.007) and IZ (41.3 vs. 55.5 mm, p = 0.02) than fast spindle generators.

CONCLUSIONS

Our data indicate alterations in the EEG slow spindles after resective epilepsy surgery. Fast spindle generation on the contrary did not change after surgery. Although the study is limited by its retrospective nature, lack of healthy controls, and reduced cortical spatial sampling, our findings suggest a spatial relationship between the slow spindles and the epileptogenic generators.

Reduced evoked cortical beta and gamma activity and neuronal synchronization in succinic semialdehyde dehydrogenase deficiency, a disorder of γ-aminobutyric acid metabolism

Succinic semialdehyde dehydrogenase deficiency is a rare autosomal recessively inherited metabolic disorder of γ-aminobutyric acid catabolism manifested by intellectual disability, expressive aphasia, movement disorders, psychiatric ailments and epilepsy. Subjects with succinic semialdehyde dehydrogenase deficiency are characterized by elevated γ-aminobutyric acid and related metabolites, such as γ-guanidinobutyric acid, and an age-dependent downregulation of cerebral γ-aminobutyric acid receptors. These findings indicate impaired γ-aminobutyric acid and γ-aminobutyric acid sub-type A (GABAA) receptor signalling as major factors underlying the pathophysiology of this neurometabolic disorder. We studied the cortical oscillation patterns and their relationship with γ-aminobutyric acid metabolism in 18 children affected by this condition and 10 healthy controls. Using high-density EEG, we recorded somatosensory cortical responses and resting-state activity. Using electrical source imaging, we estimated the relative power changes (compared with baseline) in both stimulus-evoked and stimulus-induced responses for physiologically relevant frequency bands and resting-state power. Stimulus-evoked oscillations are phase locked to the stimulus, whereas induced oscillations are not. Power changes for both evoked and induced responses as well as resting-state power were correlated with plasma γ-aminobutyric acid and γ-guanidinobutyric acid concentrations and with cortical γ-aminobutyric acid measured by proton magnetic resonance spectroscopy. Plasma γ-aminobutyric acid, γ-guanidinobutyric acid and cortical γ-aminobutyric acid were higher in patients than in controls (P < 0.001 for both). Beta and gamma relative power were suppressed for evoked responses in patients versus controls (P < 0.01). No group differences were observed for induced activity (P > 0.05). The mean gamma frequency of evoked responses was lower in patients versus controls (P = 0.002). Resting-state activity was suppressed in patients for theta (P = 0.011) and gamma (P < 0.001) bands. Evoked power changes were inversely correlated with plasma γ-aminobutyric acid and with γ-guanidinobutyric acid for beta (P < 0.001) and gamma (P < 0.001) bands. Similar relationships were observed between the evoked power changes and cortical γ-aminobutyric acid for all tested areas in the beta band (P < 0.001) and for the posterior cingulate gyrus in the gamma band (P < 0.001). We also observed a negative correlation between resting-state activity and plasma γ-aminobutyric acid and γ-guanidinobutyric acid for theta (P < 0.001; P = 0.003), alpha (P = 0.003; P = 0.02) and gamma (P = 0.02; P = 0.01) bands. Our findings indicate that increased γ-aminobutyric acid concentration is associated with reduced sensory-evoked beta and gamma activity and impaired neuronal synchronization in patients with succinic semialdehyde dehydrogenase deficiency. This further elucidates the pathophysiology of this neurometabolic disorder and serves as a potential biomarker for therapeutic trials.

Spike propagation mapping reveals effective connectivity and predicts surgical outcome in epilepsy

Neurosurgical intervention is the best available treatment for selected patients with drug resistant epilepsy. For these patients, surgical planning requires biomarkers that delineate the epileptogenic zone, the brain area that is indispensable for the generation of seizures. Interictal spikes recorded with electrophysiological techniques are considered key biomarkers of epilepsy. Yet, they lack specificity, mostly because they propagate across brain areas forming networks. Understanding the relationship between interictal spike propagation and functional connections among the involved brain areas may help develop novel biomarkers that can delineate the epileptogenic zone with high precision. Here, we reveal the relationship between spike propagation and effective connectivity among onset and areas of spread and assess the prognostic value of resecting these areas. We analysed intracranial EEG data from 43 children with drug resistant epilepsy who underwent invasive monitoring for neurosurgical planning. Using electric source imaging, we mapped spike propagation in the source domain and identified three zones: onset, early-spread and late-spread. For each zone, we calculated the overlap and distance from surgical resection. We then estimated a virtual sensor for each zone and the direction of information flow among them via Granger causality. Finally, we compared the prognostic value of resecting these zones, the clinically-defined seizure onset zone and the spike onset on intracranial EEG channels by estimating their overlap with resection. We observed a spike propagation in source space for 37 patients with a median duration of 95 ms (interquartile range: 34-206), a spatial displacement of 14 cm (7.5-22 cm) and a velocity of 0.5 m/s (0.3-0.8 m/s). In patients with good surgical outcome (25 patients, Engel I), the onset had higher overlap with resection [96% (40-100%)] than early-spread [86% (34-100%), P = 0.01] and late-spread [59% (12-100%), P = 0.002], and it was also closer to resection than late-spread [5 mm versus 9 mm, P = 0.007]. We found an information flow from onset to early-spread in 66% of patients with good outcomes, and from early-spread to onset in 50% of patients with poor outcome. Finally, resection of spike onset, but not area of spike spread or the seizure onset zone, predicted outcome with positive predictive value of 79% and negative predictive value of 56% (P = 0.04). Spatiotemporal mapping of spike propagation reveals information flow from onset to areas of spread in epilepsy brain. Surgical resection of the spike onset disrupts the epileptogenic network and may render patients with drug resistant epilepsy seizure-free without having to wait for a seizure to occur during intracranial monitoring.

Electromagnetic source imaging predicts surgical outcome in children with focal cortical dysplasia

OBJECTIVE

To evaluate the diagnostic accuracy of electromagnetic source imaging (EMSI) in localizing spikes and predict surgical outcome in children with drug resistant epilepsy (DRE) due to focal cortical dysplasia (FCD).

METHODS

We retrospectively analyzed magnetoencephalography (MEG) and high-density (HD-EEG) data from 23 children with FCD-associated DRE who underwent intracranial EEG and surgery. We localized spikes using equivalent current dipole (ECD) fitting, dipole clustering, and dynamical statistical parametric mapping (dSPM) on EMSI, electric source imaging (ESI), and magnetic source imaging (MSI). We calculated the distance from the seizure onset zone (DSOZ) and resection (DRES). We estimated receiver operating characteristic (ROC) curves with Youden's index (J) to predict outcome.

RESULTS

EMSI presented shorter DSOZ (15.18 ± 9.06 mm) and DRES (8.56 ± 6.24 mm) compared to ESI (DSOZ: 25.04 ± 16.20 mm, p < 0.009; DRES: 18.88 ± 17.30 mm, p < 0.03) and MSI (DSOZ: 23.37 ± 8.98 mm, p < 0.03; DRES: 15.51 ± 10.11 mm, p < 0.02) for clustering in patients with good outcome. Clustering showed shorter DSOZ and DRES compared to ECD fitting and dSPM (p < 0.05). EMSI had higher performance as outcome predictor (J = 70.63%) compared to ESI (J = 41.27%) and MSI (J = 33.33%) for clustering.

CONCLUSIONS

EMSI provides superior localization and improved predictive performance than individual modalities.

SIGNIFICANCE

EMSI can help the surgical planning and facilitate the localization of epileptogenic foci.

Development and characterization of transfontanelle photoacoustic imaging system for detection of intracranial hemorrhages and measurement of brain oxygenation: Ex-vivo

We have developed and optimized an imaging system to study and improve the detection of brain hemorrhage and to quantify oxygenation. Since this system is intended to be used for brain imaging in neonates through the skull opening, i.e., fontanelle, we called it, Transfontanelle Photoacoustic Imaging (TFPAI) system. The system is optimized in terms of optical and acoustic designs, thermal safety, and mechanical stability. The lower limit of quantification of TFPAI to detect the location of hemorrhage and its size is evaluated using in-vitro and ex-vivo experiments. The capability of TFPAI in measuring the tissue oxygenation and detection of vasogenic edema due to brain blood barrier disruption are demonstrated. The results obtained from our experimental evaluations strongly suggest the potential utility of TFPAI, as a portable imaging modality in the neonatal intensive care unit. Confirmation of these findings in-vivo could facilitate the translation of this promising technology to the clinic.

Functional connectivity discriminates epileptogenic states and predicts surgical outcome in children with drug resistant epilepsy

Normal brain functioning emerges from a complex interplay among regions forming networks. In epilepsy, these networks are disrupted causing seizures. Highly connected nodes in these networks are epilepsy surgery targets. Here, we assess whether functional connectivity (FC) using intracranial electroencephalography can quantify brain regions epileptogenicity and predict surgical outcome in children with drug resistant epilepsy (DRE). We computed FC between electrodes on different states (i.e. interictal without spikes, interictal with spikes, pre-ictal, ictal, and post-ictal) and frequency bands. We then estimated the electrodes' nodal strength. We compared nodal strength between states, inside and outside resection for good- (n = 22, Engel I) and poor-outcome (n = 9, Engel II-IV) patients, respectively, and tested their utility to predict the epileptogenic zone and outcome. We observed a hierarchical epileptogenic organization among states for nodal strength: lower FC during interictal and pre-ictal states followed by higher FC during ictal and post-ictal states (p < 0.05). We further observed higher FC inside resection (p < 0.05) for good-outcome patients on different states and bands, and no differences for poor-outcome patients. Resection of nodes with high FC was predictive of outcome (positive and negative predictive values: 47-100%). Our findings suggest that FC can discriminate epileptogenic states and predict outcome in patients with DRE.

Non-invasive mapping of epileptogenic networks predicts surgical outcome

Epilepsy is increasingly considered a disorder of brain networks. Studying these networks with functional connectivity can help identify hubs that facilitate the spread of epileptiform activity. Surgical resection of these hubs may lead patients who suffer from drug-resistant epilepsy to seizure freedom. Here, we aim to map non-invasively epileptogenic networks, through the virtual implantation of sensors estimated with electric and magnetic source imaging, in patients with drug-resistant epilepsy. We hypothesize that highly connected hubs identified non-invasively with source imaging can predict the epileptogenic zone and the surgical outcome better than spikes localized with conventional source localization methods (dipoles). We retrospectively analysed simultaneous high-density electroencephalography (EEG) and magnetoencephalography data recorded from 37 children and young adults with drug-resistant epilepsy who underwent neurosurgery. Using source imaging, we estimated virtual sensors at locations where intracranial EEG contacts were placed. On data with and without spikes, we computed undirected functional connectivity between sensors/contacts using amplitude envelope correlation and phase locking value for physiologically relevant frequency bands. From each functional connectivity matrix, we generated an undirected network containing the strongest connections within sensors/contacts using the minimum spanning tree. For each sensor/contact, we computed graph centrality measures. We compared functional connectivity and their derived graph centrality of sensors/contacts inside resection for good (n = 22, ILAE I) and poor (n = 15, ILAE II-VI) outcome patients, tested their ability to predict the epileptogenic zone in good-outcome patients, examined the association between highly connected hubs removal and surgical outcome and performed leave-one-out cross-validation to support their prognostic value. We also compared the predictive values of functional connectivity with those of dipoles. Finally, we tested the reliability of virtual sensor measures via Spearman's correlation with intracranial EEG at population- and patient-level. We observed higher functional connectivity inside than outside resection (P < 0.05, Wilcoxon signed-rank test) for good-outcome patients, on data with and without spikes across different bands for intracranial EEG and electric/magnetic source imaging and few differences for poor-outcome patients. These functional connectivity measures were predictive of both the epileptogenic zone and outcome (positive and negative predictive values ≥55%, validated using leave-one-out cross-validation) outperforming dipoles on spikes. Significant correlations were found between source imaging and intracranial EEG measures (0.4 ≤ rho ≤ 0.9, P < 0.05). Our findings suggest that virtual implantation of sensors through source imaging can non-invasively identify highly connected hubs in patients with drug-resistant epilepsy, even in the absence of frank epileptiform activity. Surgical resection of these hubs predicts outcome better than dipoles.

An approach for reliably identifying high-frequency oscillations and reducing false-positive detections

OBJECTIVE

Aiming to improve the feasibility and reliability of using high-frequency oscillations (HFOs) for translational studies of epilepsy, we present a pipeline with features specifically designed to reject false positives for HFOs to improve the automatic HFO detector.

METHODS

We presented an integrated, multi-layered procedure capable of automatically rejecting HFOs from a variety of common false positives, such as motion, background signals, and sharp transients. This method utilizes a time-frequency contour approach that embeds three different layers including peak constraints, power thresholds, and morphological identification to discard false positives. Four experts were involved in rating detected HFO events that were randomly selected from different posttraumatic epilepsy (PTE) animals for a comprehensive evaluation.

RESULTS

The algorithm was run on 768-h recordings of intracranial electrodes in 48 PTE animals. A total of 453 917 HFOs were identified by initial HFO detection, of which 450 917 were implemented for HFO refinement and 203 531 events were retained. Random sampling was used to evaluate the performance of the detector. The HFO detection yielded an overall accuracy of 0.95 ± 0.03 , with precision, recall, and F1 scores of 0.92 ± 0.05 , 0.99 ± 0.01 , and 0.94 ± 0.03 , respectively. For the HFO classification, our algorithm obtained an accuracy of 0.97 ± 0.02 . For the inter-rater reliability of algorithm evaluation, the agreement among four experts was 0.94 ± 0.03 for HFO detection and 0.85 ± 0.04 for HFO classification.

SIGNIFICANCE

Our approach shows that a segregated pipeline design with a focus on false-positive rejection can improve the detection efficiency and provide reliable results. This pipeline does not require customization and uses fixed parameters, making it highly feasible and translatable for basic and clinical applications of epilepsy.

Presurgical accuracy of dipole clustering in MRI-negative pediatric patients with epilepsy: Validation against intracranial EEG and resection

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.

Virtual implantation using conventional scalp EEG delineates seizure onset and predicts surgical outcome in children with epilepsy

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.

Protocol for multicentre comparison of interictal high-frequency oscillations as a predictor of seizure freedom

In drug-resistant focal epilepsy, interictal high-frequency oscillations (HFOs) recorded from intracranial EEG (iEEG) may provide clinical information for delineating epileptogenic brain tissue. The iEEG electrode contacts that contain HFO are hypothesized to delineate the epileptogenic zone; their resection should then lead to postsurgical seizure freedom. We test whether our prospective definition of clinically relevant HFO is in agreement with postsurgical seizure outcome. The algorithm is fully automated and is equally applied to all data sets. The aim is to assess the reliability of the proposed detector and analysis approach. We use an automated data-independent prospective definition of clinically relevant HFO that has been validated in data from two independent epilepsy centres. In this study, we combine retrospectively collected data sets from nine independent epilepsy centres. The analysis is blinded to clinical outcome. We use iEEG recordings during NREM sleep with a minimum of 12 epochs of 5 min of NREM sleep. We automatically detect HFO in the ripple (80-250 Hz) and in the fast ripple (250-500 Hz) band. There is no manual rejection of events in this fully automated algorithm. The type of HFO that we consider clinically relevant is defined as the simultaneous occurrence of a fast ripple and a ripple. We calculate the temporal consistency of each patient's HFO rates over several data epochs within and between nights. Patients with temporal consistency <50% are excluded from further analysis. We determine whether all electrode contacts with high HFO rate are included in the resection volume and whether seizure freedom (ILAE 1) was achieved at ≥2 years follow-up. Applying a previously validated algorithm to a large cohort from several independent epilepsy centres may advance the clinical relevance and the generalizability of HFO analysis as essential next step for use of HFO in clinical practice.

Mapping the brain basis of feelings, emotions and much more: A special issue focused on 'The Human Affectome'

The Human Affectome Project was launched by the non-profit organization Neuroqualia (www.neuroqualia.org) in 2015 with the seemingly impossible goal: To map a psychological process and form possible definitions and working models for affective states and related emotions. Twelve reviews based on emotions, feelings and motivation were written dedicated to mapping the brain basis of affect. A capstone piece 'The Human Affectome' provides a foundation for the special issue by giving detailed up-to-date definitions for key terms including feeling, affect, emotion and mood. Critically, the piece offers an overall model synthesizing three main features of affect: valence, motivation, and arousal. Affect itself is explored as the main umbrella function capturing all feeling states and related processes. Overall, the project and the special issue has been a highly successful interdisciplinary effort producing a novel approach that can be used to understand, guide and revise contemporary research on the brain basis of feeling and how diverse feeling states interact with each other in typical and atypical fashions.

Novel User-Friendly Application for MRI Segmentation of Brain Resection following Epilepsy Surgery

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.

Case Report: Laser Ablation Guided by State of the Art Source Imaging Ends an Adolescent's 16-Year Quest for Seizure Freedom

Epilepsy surgery is the most effective therapeutic approach for children with drug resistant epilepsy (DRE). Recent advances in neurosurgery, such as the Laser Interstitial Thermal Therapy (LITT), improved the safety and non-invasiveness of this method. Electric and magnetic source imaging (ESI/MSI) plays critical role in the delineation of the epileptogenic focus during the presurgical evaluation of children with DRE. Yet, they are currently underutilized even in tertiary epilepsy centers. Here, we present a case of an adolescent who suffered from DRE for 16 years and underwent surgery at Cook Children's Medical Center (CCMC). The patient was previously evaluated in a level 4 epilepsy center and treated with multiple antiseizure medications for several years. Presurgical evaluation at CCMC included long-term video electroencephalography (EEG), magnetoencephalography (MEG) with simultaneous conventional EEG (19 channels) and high-density EEG (256 channels) in two consecutive sessions, MRI, and fluorodeoxyglucose - positron emission tomography (FDG-PET). Video long-term EEG captured nine focal-onset clinical seizures with a maximal evolution over the right frontal/frontal midline areas. MRI was initially interpreted as non-lesional. FDG-PET revealed a small region of hypometabolism at the anterior right superior temporal gyrus. ESI and MSI performed with dipole clustering showed a tight cluster of dipoles in the right anterior insula. The patient underwent intracranial EEG which indicated the right anterior insular as seizure onset zone. Eventually LITT rendered the patient seizure free (Engel 1; 12 months after surgery). Retrospective analysis of ESI and MSI clustered dipoles found a mean distance of dipoles from the ablated volume ranging from 10 to 25 mm. Our findings highlight the importance of recent technological advances in the presurgical evaluation and surgical treatment of children with DRE, and the underutilization of epilepsy surgery in children with DRE.

Electric Source Imaging on Intracranial EEG Localizes Spatiotemporal Propagation of Interictal Spikes in Children with Epilepsy

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.

Mapping Functional Connectivity of Epileptogenic Networks through Virtual Implantation

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.

Mapping Propagation of Interictal Spikes, Ripples, and Fast Ripples in Intracranial EEG of Children with Refractory Epilepsy

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.

Localizing the Epileptogenic Zone with Novel Biomarkers

Several noninvasive methods, such as high-density EEG or magnetoencephalography, are currently used to delineate the epileptogenic zone (EZ) during the presurgical evaluation of patients with drug resistant epilepsy (DRE). Yet, none of these methods can reliably identify the EZ by their own. In most cases a multimodal approach is needed. Challenging cases often require the implantation of intracranial electrodes, either through stereo-taxic EEG or electro-corticography. Recently, a growing body of literature introduces novel biomarkers of epilepsy that can be used for analyzing both invasive as well as noninvasive electrophysiological data. Some of these biomarkers are able to delineate the EZ with high precision, augment the presurgical evaluation, and predict the surgical outcome of patients with DRE undergoing surgery. However, the use of these epilepsy biomarkers in clinical practice is limited. Here, we summarize and discuss the latest technological advances in the presurgical neurophysiological evaluation of children with DRE with emphasis on electric and magnetic source imaging, high frequency oscillations, and functional connectivity.

Association between semiologic, autonomic, and electrographic seizure characteristics in children with generalized tonic-clonic seizures

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.

Aberrant somatosensory phase synchronization in children with hemiplegic cerebral palsy

Children with hemiplegic cerebral palsy (HCP) often show disturbances of somatosensation. Despite extensive evidence of somatosensory deficits, neurophysiological alterations associated with somatosensory deficits in children with HCP have not been elucidated. Here, we aim to assess phase synchrony within and between contralateral primary (S1) and secondary (S2) somatosensory areas in children with HCP. Intra-regional and inter-regional phase synchronizations within and between S1 and S2 were estimated from somatosensory evoked fields (SEFs) in response to passive pneumatic stimulation of contralateral upper extremities and recorded with pediatric magnetoencephalography (MEG) in children with HCP and typically developing (TD) children. We found aberrant phase synchronizations within S1 and between S1 and S2 in both hemispheres in children with HCP. Specifically, the less-affected (LA) hemisphere demonstrated diminished phase synchronizations after the stimulus onset up to ~120 ms compared to the more-affected (MA) hemisphere and the dominant hemisphere of TD children, while the MA hemisphere showed enhanced phase synchronizations after ~100 ms compared to the LA hemisphere and the TD dominant hemisphere. Our findings indicate abnormal somatosensory functional connectivity in both hemispheres of children with HCP.

Noninvasive Mapping of Ripple Onset Predicts Outcome in Epilepsy Surgery

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.

Source imaging of seizure onset predicts surgical outcome in pediatric epilepsy

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.

Atypical spatiotemporal activation of cerebellar lobules during emotional face processing in adolescents with autism

Autism spectrum disorder (ASD) is characterized by social deficits and atypical facial processing of emotional expressions. The underlying neuropathology of these abnormalities is still unclear. Recent studies implicate cerebellum in emotional processing; other studies show cerebellar abnormalities in ASD. Here, we elucidate the spatiotemporal activation of cerebellar lobules in ASD during emotional processing of happy and angry faces in adolescents with ASD and typically developing (TD) controls. Using magnetoencephalography, we calculated dynamic statistical parametric maps across a period of 500 ms after emotional stimuli onset and determined differences between group activity to happy and angry emotions. Following happy face presentation, adolescents with ASD exhibited only left‐hemispheric cerebellar activation in a cluster extending from lobule VI to lobule V (compared to TD controls). Following angry face presentation, adolescents with ASD exhibited only midline cerebellar activation (posterior IX vermis). Our findings indicate an early (125–175 ms) overactivation in cerebellar activity only for happy faces and a later overactivation for both happy (250–450 ms) and angry (250–350 ms) faces in adolescents with ASD. The prioritized hemispheric activity (happy faces) could reflect the promotion of a more flexible and adaptive social behavior, while the latter midline activity (angry faces) may guide conforming behavior.

The neuroscience of positive emotions and affect: Implications for cultivating happiness and wellbeing

This review paper provides an integrative account regarding neurophysiological correlates of positive emotions and affect that cumulatively contribute to the scaffolding for happiness and wellbeing in humans and other animals. This paper reviews the associations among neurotransmitters, hormones, brain networks, and cognitive functions in the context of positive emotions and affect. Consideration of lifespan developmental perspectives are incorporated, and we also examine the impact of healthy social relationships and environmental contexts on the modulation of positive emotions and affect. The neurophysiological processes that implement positive emotions are dynamic and modifiable, and meditative practices as well as flow states that change patterns of brain function and ultimately support wellbeing are also discussed. This review is part of "The Human Affectome Project" (http://neuroqualia.org/background.php), and in order to advance a primary aim of the Human Affectome Project, we also reviewed relevant linguistic dimensions and terminology that characterizes positive emotions and wellbeing. These linguistic dimensions are discussed within the context of the neuroscience literature with the overarching goal of generating novel recommendations for advancing neuroscience research on positive emotions and wellbeing.

Abnormal Nutritive Sucking as an Indicator of Neonatal Brain Injury

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.

Succinic Semialdehyde Dehydrogenase Deficiency: Review of the Natural History Study

OBJECTIVE

The SSADHD Natural History Study was initiated in 2019 to define the natural course and identify biomarkers correlating with severity.

METHODS

The study is conducted by 4 institutions: BCH (US clinical), WSU (bioanalytical core), USF (biostatistical core), and Heidelberg (iNTD), with support from the family advocacy group (SSADH Association). Recruitment goals were to study 20 patients on-site at BCH, 10 with iNTD, and 25 as a standard-of care cohort.

RESULTS

At this half-way point of this longitudinal study, 28 subjects have been recruited (57% female, mean 9 years, range 18 months-40 years). Epilepsy is present in half and increases in incidence and severity, as do psychiatric symptoms, in adolescence and adulthood. The average Full Scale IQ (FSIQ) was 53 (Verbal score of 56, Non Verbal score of 49), and half scored as having ASD. Although there was no correlation between gene variant and phenotypic severity, there were extreme cases of lowest functioning in one individual and highest in another that may have genotype-phenotype correlation. The most common EEG finding was mild background slowing with rare epileptiform activity, whereas high-density EEG and magnetoencephalography showed reduction in the gamma frequency band consistent with GABAergic dysfunction. MR spectroscopy showed elevations in the GABA/NAA ratio in all regions studied with no crossover between subjects and controls.

CONCLUSIONS

The SSADH Natural History Study is providing a unique opportunity to study the complex pathophysiology longitudinally and derive electrophysiologic, neuroimaging, and laboratory data for correlation and to serve as biomarkers for clinical trials and prognostic assessments in this ultra-rare inherited disorder of GABA metabolism.

Clinical Correlation Advised: Measuring Functional Connectivity in PNES

Activation of Functional Brain Networks in Children With Psychogenic Non-Epileptic Seizures

Radmanesh M, Jalili M, Kozlowska K. Front Hum Neurosci. 2020. doi:10.3389/fnhum.2020.00339.

Objectives:

Psychogenic nonepileptic seizures (PNES) have been hypothesized to emerge in the context of neural networks instability. To explore this hypothesis in children, we applied a graph theory approach to examine connectivity in neural networks in the resting-state electroencephalogram in 35 children with PNES, 31 children with other functional neurological symptoms (but no PNES), and 75 healthy controls.

Methods:

The networks were extracted from Laplacian-transformed time series by a coherence connectivity estimation method.

Results:

Children with PNES (vs controls) showed widespread changes in network metrics: increased global efficiency (α and β bands), increased local efficiency (γ band), and increased modularity (γ and α bands). Compared to controls, they also had higher levels of autonomic arousal (eg, lower heart variability); more anxiety, depression, and stress on the Depression Anxiety and Stress Scales; and more adverse childhood experiences on the Early Life Stress Questionnaire. Increases in network metrics correlated with arousal. Children with other functional neurological symptoms (but no PNES) showed scattered and less pronounced changes in network metrics.

Conclusion:

The results indicate that children with PNES present with increased activation of neural networks coupled with increased physiological arousal. Although this shift in functional organization may confer a short-term adaptive advantage—one that facilitates neural communication and the child’s capacity to respond self-protectively in the face of stressful life events—it may also have a significant biological cost. It may predispose the child’s neural networks to periods of instability—presenting clinically as PNES—when the neural networks are faced with perturbations in energy flow or with additional demands.

Photoplethysmographic evaluation of generalized tonic-clonic seizures

OBJECTIVE

Photoplethysmography (PPG) is an optical technique measuring variations of blood perfusion in peripheral tissues. We evaluated alterations in PPG signals in relationship to the occurrence of generalized tonic-clonic seizures (GTCSs) in patients with epilepsy to evaluate the feasibility of seizure detection.

METHODS

During electroencephalographic (EEG) long-term monitoring, patients wore portable wristband sensor(s) on their wrists or ankles recording PPG signals. We analyzed PPG signals during three time periods, which were defined with respect to seizures detected on EEG: (1) baseline (>30 minutes prior to seizure), (2) preseizure period, and (3) postseizure period. Furthermore, we selected five random control segments during seizure-free periods. PPG features, including frequency, amplitude, duration, slope, smoothness, and area under the curve, were automatically calculated. We used a linear mixed-effect model to evaluate changes in PPG features between different time periods in an attempt to identify signal changes that detect seizures.

RESULTS

We prospectively enrolled 174 patients from the epilepsy monitoring unit at Boston Children's Hospital. Twenty-five GTCSs were recorded from 13 patients. Data from the first recorded GTCS of each patient were included in the analysis. We observed an increase in PPG frequency during pre- and postseizure periods that was higher than the changes during seizure-free periods (frequency increase: preseizure = 0.22 Hz, postseizure = 0.58 Hz vs changes during seizure-free period = 0.05 Hz). The PPG slope decreased significantly by 56.71 nW/s during preseizure periods compared to seizure-free periods. Additionally, the smoothness increased significantly by 0.22 nW/s during the postseizure period compared to seizure-free periods.

SIGNIFICANCE

Monitoring of PPG signals may assist in the detection of GTCSs in patients with epilepsy. PPG may serve as a promising biomarker for future seizure detection systems and may contribute to future seizure prediction systems.

The Emergence of Hierarchical Somatosensory Processing in Late Prematurity

The somatosensory system has a hierarchical organization. Information processing increases in complexity from the contralateral primary sensory cortex to bilateral association cortices and this is represented by a sequence of somatosensory-evoked potentials recorded with scalp electroencephalographies. The mammalian somatosensory system matures over the early postnatal period in a rostro-caudal progression, but little is known about the development of hierarchical information processing in the human infant brain. To investigate the normal human development of the somatosensory hierarchy, we recorded potentials evoked by mechanical stimulation of hands and feet in 34 infants between 34 and 42 weeks corrected gestational age, with median postnatal age of 3 days. We show that the shortest latency potential was evoked for both hands and feet at all ages with a contralateral somatotopic source in the primary somatosensory cortex (SI). However, the longer latency responses, localized in SI and beyond, matured with age. They gradually emerged for the foot and, although always present for the hand, showed a shift from purely contralateral to bilateral hemispheric activation. These results demonstrate the rostro-caudal development of human somatosensory hierarchy and suggest that the development of its higher tiers is complete only just before the time of normal birth.

Noninvasive Localization of High-Frequency Oscillations in Children with Epilepsy: Validation against Intracranial Gold-Standard

INTRODUCTION

Patients with medically refractory epilepsy (MRE) need surgical resection of the epileptogenic zone (EZ) to gain seizure-freedom. High-frequency oscillations (HFOs, > 80 Hz) are promising biomarkers of the EZ that are typically localized using intracranial electroencephalography (icEEG). The goal of this study was to localize the cortical generators of HFOs non-invasively using high-density (HD) EEG and magnetoencephalography (MEG) and validate the localization against the gold-standard given by the icEEGdefined HFO-zone.

METHODS

We analyzed simultaneous HDEEG and MEG data from six children with MRE who underwent icEEG and surgery. We detected interictal HFOs (80-160 Hz) on HD-EEG and MEG separately, using an inhouse automatic detector followed by visual human review, and distinguished between HFOs with and without spikes. We localized the cortical generators of each HFO on HD-EEG or MEG using the wavelet Maximum Entropy on the Mean (wMEM). For the HFOs localized in the brain area covered by icEEG, we estimated the localization error (E_loc) with respect to the gold-standard, and classified them as either concordant (E_loc≤15mm) or not.

RESULTS

We found that: (i) HD-EEG presented a higher rate of HFOs than MEG (1 vs 0.5 HFOs/min, p=0.031); (ii) HFOs without spikes were more likely to be localized outside the brain regions of interest (i.e. covered by icEEG) than HFOs with spikes; and (iii) both HD-EEG and MEG showed high precision to the gold-standard (92% and 96%).

CONCLUSION

We reported quantitative evidence that HDEEG and MEG can localize the HFO cortical generators with high precision to the icEEG gold-standard in children with MRE, suggesting that they may possibly limit the need for icEEG prior to surgery. We also showed that HFOs with spikes on HD-EEG/MEG are more likely to be epileptogenic than those independent from spikes, which may represent physiological events from normal brain.

Scalp ripples as prognostic biomarkers of epileptogenicity in pediatric surgery

Objective

To assess the ability of high‐density Electroencephalography (HD‐EEG) and magnetoencephalography (MEG) to localize interictal ripples, distinguish between ripples co‐occurring with spikes (ripples‐on‐spike) and independent from spikes (ripples‐alone), and evaluate their localizing value as biomarkers of epileptogenicity in children with medically refractory epilepsy.

Methods

We retrospectively studied 20 children who underwent epilepsy surgery. We identified ripples on HD‐EEG and MEG data, localized their generators, and compared them with intracranial EEG (icEEG) ripples. When ripples and spikes co‐occurred, we performed source imaging distinctly on the data above 80 Hz (to localize ripples) and below 70 Hz (to localize spikes). We assessed whether missed resection of ripple sources predicted poor outcome, separately for ripples‐on‐spikes and ripples‐alone. Similarly, predictive value of spikes was calculated.

Results

We observed scalp ripples in 16 patients (10 good outcome). Ripple sources were highly concordant to the icEEG ripples (HD‐EEG concordance: 79%; MEG: 83%). When ripples and spikes co‐occurred, their sources were spatially distinct in 83‐84% of the cases. Removing the sources of ripples‐on‐spikes predicted good outcome with 90% accuracy for HD‐EEG (P = 0.008) and 86% for MEG (P = 0.044). Conversely, removing ripples‐alone did not predict outcome. Resection of spike sources (generated at the same time as ripples) predicted good outcome for HD‐EEG (P = 0.036; accuracy = 87%), while did not reach significance for MEG (P = 0.1; accuracy = 80%).

Interpretation

HD‐EEG and MEG localize interictal ripples with high precision in children with refractory epilepsy. Scalp ripples‐on‐spikes are prognostic, noninvasive biomarkers of epileptogenicity, since removing their cortical generators predicts good outcome. Conversely, scalp ripples‐alone are most likely generated by non‐epileptogenic areas.

Ictal and interictal source imaging on intracranial EEG predicts epilepsy surgery outcome in children with focal cortical dysplasia

OBJECTIVE

To localize the seizure onset zone (SOZ) and irritative zone (IZ) using electric source imaging (ESI) on intracranial EEG (iEEG) and assess their clinical value in predicting epilepsy surgery outcome in children with focal cortical dysplasia (FCD).

METHODS

We analyzed iEEG data from 25 children with FCD-associated medically refractory epilepsy (MRE) who underwent surgery. We performed ESI on ictal onset to localize SOZ (ESI-SOZ) and on interictal discharges to localize IZ (ESI-IZ). We tested whether resection of ESI-SOZ and ESI-IZ predicted good surgical outcome (Engel 1). We further compared the prediction performance of ESI-SOZ and ESI-IZ to those of SOZ and IZ defined using conventional methods, i.e. by identifying iEEG-contacts showing ictal onsets (conventional-SOZ) or being the most interictally active (conventional-IZ).

RESULTS

The proximity of ESI-SOZ (p = 0.043, odds-ratio = 3.9) and ESI-IZ (p = 0.011, odds-ratio = 7.04) to resection has higher effect on patients' outcome than proximity of conventional-SOZ (p = 0.17, odds-ratio = 1.7) and conventional-IZ (p = 0.038, odds-ratio = 2.6). Resection of ESI-SOZ and ESI-IZ presented higher discriminative power in predicting outcome (68% and 60%) than conventional-SOZ and conventional-IZ (48% and 53%).

CONCLUSIONS

Localizing SOZ and IZ via ESI on iEEG offers higher predictive value compared to conventional-iEEG interpretation.

SIGNIFICANCE

iEEG-ESI may help surgical planning and facilitate prognostic assessment of children with FCD-associated MRE.

Nutritive sucking abnormalities and brain microstructural abnormalities in infants with established brain injury: a pilot study

OBJECTIVE

To determine the relationship between nutritive sucking and microstructural integrity of sensorimotor tracts in newborns with brain injury.

STUDY DESIGN

Diffusion imaging was performed in ten newborns with brain injury. Nutritive sucking was assessed using Nfant®. The motor, sensory, and corpus callosum tracts were reconstructed via tractography. Fractional anisotropy, radial, axial, and mean diffusivity were estimated for these tracts. Multiple regression models were developed to test the association between sucking features and diffusion parameters.

RESULTS

Low-sucking smoothness correlated with low-fractional anisotropy of motor tracts (p = 0.0096). High-sucking irregularity correlated with high-mean diffusivity of motor (p = 0.030) and corpus callosum tracts (p = 0.032). For sensory tracts, high-sucking irregularity (p = 0.018) and low-smoothness variability (p = 0.002) correlated with high-mean diffusivity.

INTERPRETATION

We show a correlation between neuroimaging-demonstrated microstructural brain abnormalities and variations in sucking patterns of newborns. The consistency of this relationship should be shown on larger cohorts.

Maturation of Corticospinal Tracts in Children With Hemiplegic Cerebral Palsy Assessed by Diffusion Tensor Imaging and Transcranial Magnetic Stimulation

Aim: To assess changes in the developmental trajectory of corticospinal tracts (CST) maturation in children with hemiplegic cerebral palsy (HCP).

Methods: Neuroimaging data were obtained from 36 children with HCP for both the more affected (MA) and less affected (LA) hemispheres, and, for purposes of direct comparison, between groups, 15 typically developing (TD) children. With diffusion tensor imaging (DTI), we estimated the mean fractional anisotropy (FA), axial diffusivity (AD), mean diffusivity (MD), and radial diffusivity (RD) of the corticospinal tract, parameters indicative of factors including myelination and axon density. Transcranial magnetic stimulation (TMS) was performed as a neurophysiologic measure of corticospinal tract integrity and organization. Resting motor threshold (rMT) was obtained per hemisphere, per patient.

Results: We observed a significant AD and MD developmental trajectory, both of which were inversely related to age (decrease in AD and diffusivity corresponding to increased age) in both hemispheres of TD children (p < 0.001). This maturation process was absent in both MA and LA hemispheres of children with HCP. Additionally, the TMS-derived previously established rMT developmental trajectory was preserved in the LA hemisphere of children with HCP (n = 26; p < 0.0001) but this trajectory was absent in the MA hemisphere.

Conclusions: Corticospinal tract maturation arrests in both hemispheres of children with HCP, possibly reflecting perinatal disruption of corticospinal tract myelination and axonal integrity.

Altered White Matter Connectivity Associated with Intergyral Brain Disorganization in Hemiplegic Cerebral Palsy

Despite extensive literature showing damages in the sensorimotor projection fibers of children with hemiplegic cerebral palsy (HCP), little is known about how these damages affect the global brain network. In this study, we assess the relationship between the structural integrity of sensorimotor projection fibers and the integrity of intergyral association white matter connections in children with HCP. Diffusion tensor imaging was performed in 10 children with HCP and 16 typically developing children. We estimated the regional and global white-matter connectivity using a region-of-interest (ROI)-based approach and a whole-brain gyrus-based parcellation method. Using the ROI-based approach, we tracked the spinothalamic (STh), thalamocortical (ThC), corticospinal (CST), and sensorimotor U- (SMU) fibers. Using the whole-brain parcellation method, we tracked the short-, middle-, and long-range association fibers. We observed for the more affected hemisphere of children with HCP: (i) an increase in axial diffusivity (AD), mean diffusivity (MD), and radial diffusivity (RD) for the STh and ThC fibers; (ii) a decrease in fractional anisotropy (FA) and an increase in MD and RD for the CST and SMU fibers; in (iii) a decrease in FA and an increase in AD, MD, and RD for the middle- and long-range association fibers; and (iv) an association between the integrity of sensorimotor projection and intergyral association fibers. Our findings indicate that altered structural integrity of the sensorimotor projection fibers disorganizes the intergyral association white matter connections among local and distant regions in children with HCP.

Assessing the localization accuracy and clinical utility of electric and magnetic source imaging in children with epilepsy

OBJECTIVE

To evaluate the accuracy and clinical utility of conventional 21-channel EEG (conv-EEG), 72-channel high-density EEG (HD-EEG) and 306-channel MEG in localizing interictal epileptiform discharges (IEDs).

METHODS

Twenty-four children who underwent epilepsy surgery were studied. IEDs on conv-EEG, HD-EEG, MEG and intracranial EEG (iEEG) were localized using equivalent current dipoles and dynamical statistical parametric mapping (dSPM). We compared the localization error (E_Loc) with respect to the ground-truth Irritative Zone (IZ), defined by iEEG sources, between non-invasive modalities and the distance from resection (D_res) between good- (Engel 1) and poor-outcomes. For each patient, we estimated the resection percentage of IED sources and tested whether it predicted outcome.

RESULTS

MEG presented lower E_Loc than HD-EEG and conv-EEG. For all modalities, D_res was shorter in good-outcome than poor-outcome patients, but only the resection percentage of the ground-truth IZ and MEG-IZ predicted surgical outcome.

CONCLUSIONS

MEG localizes the IZ more accurately than conv-EEG and HD-EEG. MSI may help the presurgical evaluation in terms of patient's outcome prediction. The promising clinical value of ESI for both conv-EEG and HD-EEG prompts the use of higher-density EEG-systems to possibly achieve MEG performance.

SIGNIFICANCE

Localizing the IZ non-invasively with MSI/ESI facilitates presurgical evaluation and surgical prognosis assessment.

Surgical resection of ripple onset predicts outcome in pediatric epilepsy

OBJECTIVE

In patients with medically refractory epilepsy (MRE), interictal ripples (80-250Hz) are observed in large brain areas whose resection may be unnecessary for seizure freedom. This limits their utility as epilepsy biomarkers for surgery. We assessed the spatiotemporal propagation of interictal ripples on intracranial electroencephalography (iEEG) in children with MRE, compared it with the propagation of spikes, identified ripples that initiated propagation (onset-ripples), and evaluated their clinical value as epilepsy biomarkers.

METHODS

Twenty-seven children who underwent epilepsy surgery were studied. We identified propagation sequences of ripples and spikes across multiple iEEG contacts and calculated each ripple or spike latency from the propagation onset. We classified ripples and spikes into categories (ie, onset, spread, and isolated) based on their spatiotemporal characteristics and correlated their mean rate inside and outside resection with outcome (good outcome, Engel 1 versus poor outcome, Engel≥2). We determined, as onset-zone, spread-zone, and isolated-zone, the areas generating the corresponding ripple or spike category and evaluated the predictive value of their resection.

RESULTS

We observed ripple propagation in all patients and spike propagation in 25 patients. Mean rate of onset-ripples inside resection predicted the outcome (odds ratio = 5.37; p = 0.02) and correlated with Engel class (rho = -0.55; p = 0.003). Resection of the onset-ripple-zone was associated with good outcome (p = 0.047). No association was found for the spread-ripple-zone, isolated-ripple-zone, or any spike-zone.

INTERPRETATION

Interictal ripples propagate across iEEG contacts in children with MRE. The association between the onset-ripple-zone resection and good outcome indicates that onset-ripples are promising epilepsy biomarkers, which estimate the epileptogenic tissue better than spread-ripples or onset-spikes. Ann Neurol 2018;84:331-346.

Mapping the Spatiotemporal Evolution of Emotional Processing: An MEG Study Across Arousal and Valence Dimensions

Electrophysiological and functional neuroimaging findings indicate that the neural mechanisms underlying the processing of emotional dimensions (i.e., valence, arousal) constitute a spatially and temporally distributed emotional network, modulated by the arousal and/or valence of the emotional stimuli. We examined the time course and source distribution of gamma time-locked magnetoencephalographic activity in response to a series of emotional stimuli viewed by healthy adults. We used a beamformer and a sliding window analysis to generate a succession of spatial maps of event-related brain responses across distinct levels of valence (pleasant/unpleasant) and arousal (high/low) in 30–100 Hz. Our results show parallel emotion-related responses along specific temporal windows involving mainly dissociable neural pathways for valence and arousal during emotional picture processing. Pleasant valence was localized in the left inferior frontal gyrus, while unpleasant valence in the right occipital gyrus, the precuneus, and the left caudate nucleus. High arousal was processed by the left orbitofrontal cortex, amygdala, and inferior frontal gyrus, as well as the right middle temporal gyrus, inferior parietal lobule, and occipital gyrus. Pleasant by high arousal interaction was localized in the left inferior and superior frontal gyrus, as well as the right caudate nucleus, putamen, and gyrus rectus. Unpleasant by high arousal interaction was processed by the right superior parietal gyrus. Valence was prioritized (onset at ∼60 ms) to all other effects, while pleasant valence was short lived in comparison to unpleasant valence (offsets at ∼110 and ∼320 ms, respectively). Both arousal and valence × arousal interactions emerged relatively early (onset at ∼150 ms, and ∼170 ms, respectively). Our findings support the notion that brain regions differentiate between valence and arousal, and demonstrate, for the first time, that these brain regions may also respond to distinct combinations of these two dimensions within specific time windows.

Neonatal brain resting-state functional connectivity imaging modalities

Infancy is the most critical period in human brain development. Studies demonstrate that subtle brain abnormalities during this state of life may greatly affect the developmental processes of the newborn infants. One of the rapidly developing methods for early characterization of abnormal brain development is functional connectivity of the brain at rest. While the majority of resting-state studies have been conducted using magnetic resonance imaging (MRI), there is clear evidence that resting-state functional connectivity (rs-FC) can also be evaluated using other imaging modalities. The aim of this review is to compare the advantages and limitations of different modalities used for the mapping of infants' brain functional connectivity at rest. In addition, we introduce photoacoustic tomography, a novel functional neuroimaging modality, as a complementary modality for functional mapping of infants' brain.

Reorganization of the somatosensory cortex in hemiplegic cerebral palsy associated with impaired sensory tracts

Functional neuroimaging studies argue that sensory deficits in hemiplegic cerebral palsy (HCP) are related to deviant somatosensory processing in the ipsilesional primary somatosensory cortex (S1). A separate body of structural neuroimaging literature argues that these deficits are due to structural damage of the ascending sensory tracts (AST). The relationship between the functional and structural integrity of the somatosensory system and the sensory performance is largely unknown in HCP. To address this relationship, we combined findings from magnetoencephalography (MEG) and probabilistic diffusion tractography (PDT) in 10 children with HCP and 13 typically developing (TD) children. With MEG, we mapped the functionally active regions in the contralateral S1 during tactile stimulation of the thumb, middle, and little fingers of both hands. Using these MEG-defined functional active regions as regions of interest for PDT, we estimated the diffusion parameters of the AST. Somatosensory function was assessed via two-point discrimination tests. Our MEG data showed: (i) an abnormal somatotopic organization in all children with HCP in either one or both of their hemispheres; (ii) longer Euclidean distances between the digit maps in the S1 of children with HCP compared to TD children; (iii) suppressed gamma responses at early latencies for both hemispheres of children with HCP; and (iv) a positive correlation between the Euclidean distances and the sensory tests for the more affected hemisphere of children with HCP. Our MEG-guided PDT data showed: (i) higher mean and radian diffusivity of the AST in children with HCP; (ii) a positive correlation between the axial diffusivity of the AST with the sensory tests for the more affected hemisphere; and (iii) a negative correlation between the gamma power change and the AD of the AST for the MA hemisphere. Our findings associate for the first time bilateral cortical functional reorganization in the S1 of HCP children with abnormalities in the structural integrity of the AST, and correlate these abnormalities with behaviorally-assessed sensory deficits.

•

Hemiplegic CP exhibits abnormal somatotopy in one or both hemispheres.

•

Hemiplegic CP shows longer distances between digit somatosensory representations.

•

Suppressed gamma somatosensory responses are seen in hemiplegic CP.

•

Altered somatosensory organization is associated with deficient sensory fibers.

•

Injury of sensory fibers correlates with the severity of sensory deficits.

Consensus Paper: Cerebellum and Emotion

Over the past three decades, insights into the role of the cerebellum in emotional processing have substantially increased. Indeed, methodological refinements in cerebellar lesion studies and major technological advancements in the field of neuroscience are in particular responsible to an exponential growth of knowledge on the topic. It is timely to review the available data and to critically evaluate the current status of the role of the cerebellum in emotion and related domains. The main aim of this article is to present an overview of current facts and ongoing debates relating to clinical, neuroimaging, and neurophysiological findings on the role of the cerebellum in key aspects of emotion. Experts in the field of cerebellar research discuss the range of cerebellar contributions to emotion in nine topics. Topics include the role of the cerebellum in perception and recognition, forwarding and encoding of emotional information, and the experience and regulation of emotional states in relation to motor, cognitive, and social behaviors. In addition, perspectives including cerebellar involvement in emotional learning, pain, emotional aspects of speech, and neuropsychiatric aspects of the cerebellum in mood disorders are briefly discussed. Results of this consensus paper illustrate how theory and empirical research have converged to produce a composite picture of brain topography, physiology, and function that establishes the role of the cerebellum in many aspects of emotional processing.

Quantifying Neonatal Sucking Performance: Promise of New Methods

Neonatal feeding has been traditionally understudied so guidelines and evidence-based support for common feeding practices are limited. A major contributing factor to the paucity of evidence-based practice in this area has been the lack of simple-to-use, low-cost tools for monitoring sucking performance. We describe new methods for quantifying neonatal sucking performance that hold significant clinical and research promise. We present early results from an ongoing study investigating neonatal sucking as a marker of risk for adverse neurodevelopmental outcomes. We include quantitative measures of sucking performance to better understand how movement variability evolves during skill acquisition. Results showed the coefficient of variation of suck duration was significantly different between preterm neonates at high risk for developmental concerns (HRPT) and preterm neonates at low risk for developmental concerns (LRPT). For HRPT, results indicated the coefficient of variation of suck smoothness increased from initial feeding to discharge and remained significantly greater than healthy full-term newborns (FT) at discharge. There was no significant difference in our measures between FT and LRPT at discharge. Our findings highlight the need to include neonatal sucking assessment as part of routine clinical care in order to capture the relative risk of adverse neurodevelopmental outcomes at discharge.

Current and Emerging Potential of Magnetoencephalography in the Detection and Localization of High-Frequency Oscillations in Epilepsy

Up to one-third of patients with epilepsy are medically intractable and need resective surgery. To be successful, epilepsy surgery requires a comprehensive preoperative evaluation to define the epileptogenic zone (EZ), the brain area that should be resected to achieve seizure freedom. Due to lack of tools and methods that measure the EZ directly, this area is defined indirectly based on concordant data from a multitude of presurgical non-invasive tests and intracranial recordings. However, the results of these tests are often insufficiently concordant or inconclusive. Thus, the presurgical evaluation of surgical candidates is frequently challenging or unsuccessful. To improve the efficacy of the surgical treatment, there is an overriding need for reliable biomarkers that can delineate the EZ. High-frequency oscillations (HFOs) have emerged over the last decade as new potential biomarkers for the delineation of the EZ. Multiple studies have shown that HFOs are spatially associated with the EZ. Despite the encouraging findings, there are still significant challenges for the translation of HFOs as epileptogenic biomarkers to the clinical practice. One of the major barriers is the difficulty to detect and localize them with non-invasive techniques, such as magnetoencephalography (MEG) or scalp electroencephalography (EEG). Although most literature has studied HFOs using invasive recordings, recent studies have reported the detection and localization of HFOs using MEG or scalp EEG. MEG seems to be particularly advantageous compared to scalp EEG due to its inherent advantages of being less affected by skull conductivity and less susceptible to contamination from muscular activity. The detection and localization of HFOs with MEG would largely expand the clinical utility of these new promising biomarkers to an earlier stage in the diagnostic process and to a wider range of patients with epilepsy. Here, we conduct a thorough critical review of the recent MEG literature that investigates HFOs in patients with epilepsy, summarizing the different methodological approaches and the main findings. Our goal is to highlight the emerging potential of MEG in the non-invasive detection and localization of HFOs for the presurgical evaluation of patients with medically refractory epilepsy (MRE).