Slow Spatial Recruitment of Neocortex during Secondarily Generalized Seizures and Its Relation to Surgical Outcome

Overlap of spike and ripple propagation onset predicts surgical outcome in epilepsy

OBJECTIVE

Interictal biomarkers are critical for identifying the epileptogenic focus. However, spikes and ripples lack specificity while fast ripples lack sensitivity. These biomarkers propagate from more epileptogenic onset to areas of spread. The pathophysiological mechanism of these propagations is elusive. Here, we examine zones where spikes and high frequency oscillations co-occur (SHFO), the spatiotemporal propagations of spikes, ripples, and fast ripples, and evaluate the spike-ripple onset overlap (SRO) as an epilepsy biomarker.

METHODS

We retrospectively analyzed intracranial EEG data from 41 patients with drug-resistant epilepsy. We mapped propagations of spikes, ripples, and fast ripples, and identified their onset and spread zones, as well as SHFO and SRO. We then estimated the SRO prognostic value in predicting surgical outcome and compared it to onset and spread zones of spike, ripple, and fast ripple propagations, and SHFO.

RESULTS

We detected spikes and ripples in all patients and fast ripples in 12 patients (29%). We observed spike and ripple propagations in 40 (98%) patients. Spike and ripple onsets overlapped in 35 (85%) patients. In good outcome patients, SRO showed higher specificity and precision (p < 0.05) in predicting resection compared to onset and zones of spikes, ripples, and SHFO. Only SRO resection predicted outcome (p = 0.01) with positive and negative predictive values of 82% and 57%, respectively.

INTERPRETATION

SRO is a specific and precise biomarker of the epileptogenic zone whose removal predicts outcome. SRO is present in most patients with drug-resistant epilepsy. Such a biomarker may reduce prolonged intracranial monitoring and improve outcome.

Focal seizures induce spatiotemporally organized spiking activity in the human cortex

Epileptic seizures are debilitating because of the clinical symptoms they produce. These symptoms, in turn, may stem directly from disruptions in neural coding. Recent evidence has suggested that the specific temporal order, or sequence, of spiking across a population of cortical neurons may encode information. Here, we investigate how seizures disrupt neuronal spiking sequences in the human brain by recording multi-unit activity from the cerebral cortex in five male participants undergoing monitoring for seizures. We find that pathological discharges during seizures are associated with bursts of spiking activity across a population of cortical neurons. These bursts are organized into highly consistent and stereotyped temporal sequences. As the seizure evolves, spiking sequences diverge from the sequences observed at baseline and become more spatially organized. The direction of this spatial organization matches the direction of the ictal discharges, which spread over the cortex as traveling waves. Our data therefore suggest that seizures can entrain cortical spiking sequences by changing the spatial organization of neuronal firing, providing a possible mechanism by which seizures create symptoms.

Intranasal Administration of BDNF Improves Recovery and Promotes Neural Plasticity in a Neonatal Mouse Model of Hypoxic Ischemia

The benefit of intranasal brain derived neurotrophic factor (BDNF) treatment on cognitive function in a neonatal postnatal day 7 (P7) mouse model of hypoxic ischemia (HI) was explored. Intranasal delivery is attractive in that it can promote widespread distribution of BDNF within both the brain and spinal cord. In this study we evaluated the effectiveness of intranasal BDNF to improve cognitive recovery following HI. HI is induced via ligation of the right carotid artery followed by a 45-minute exposure to an 8% oxygen/ 92% nitrogen mixture in an enclosed chamber. Male and female pups were subjected to a 2-hour hypothermia in a temperature-controlled chamber as a standard of care. A solution of saline (control) or recombinant human BDNF (Harlan Laboratories) was administered with a Gilson pipette at the same time each day for 7 days into each nasal cavity in awake mice beginning 24 hours after HI. We evaluated cognitive recovery using the novel object recognition (NOR) and western analysis to analyze neuro-markers and brain health such as synaptophysin and microtubule associated protein -2 (MAP2). The objective of this study was to evaluate the role and therapeutic potential of BDNF in neonatal HI recovery. Our results indicate that intranasal BDNF delivered within 24 hours after HI improved object discrimination at both 28 and 42 days after HI. Our results also demonstrate increased synaptophysin and MAP2 at day 42 in HI animals that received intranasal BDNF treatment compared to HI animals that were administered saline.

Hippocampal barques and their manifestation as 14&6 Hz positive spikes during sleep

OBJECTIVE

This study investigates variations in hippocampal barque occurrence during sleep and compares findings to respective variations of their scalp manifestation as 14&6/sec positive spikes.

METHODS

From 11 epilepsy patients, 12 non-epileptogenic hippocampi with barques were identified for this study. Using the first seizure-free whole-night sleep stereo-encephalography (sEEG) recording, we performed sleep staging and measured the occurrence of barques and 14&6/sec positive spikes variants.

RESULTS

Hippocampal barques (total count: 9,183; mean count per record: 765.2 ± 251.2) occurred predominantly during non-rapid eye movement (NREM) II sleep (total: 5,744; mean: 478.6 ± 176.1; 62.2 ± 6.0%) and slow-wave sleep (SWS) (total: 2,950; mean: 245.83 ± 92.9; 32.0 ± 6.2%), with rare to occasional occurrence in NREM I (total: 85; mean: 7.0 ± 2.8; 0.9 ± 0.4%), rapid eye movement (REM) (total: 153; mean: 12.75 ± 4.0; 1.7 ± 0.6) and wakefulness (total: 251; mean: 20.9 ± 6.3; 2.9 ± 0.9%). Barque rate increased during SWS (mean: 2.7 ± 1.0 per min) compared to NREM II (2.2 ± 1.0 per min) and other states (wakefulness: 0.1 ± 0.0 per min; NREM I: 0.3 ± 0.1 per min; REM: 0.1 ± 0.0 per min). The 14&6/sec positive spikes variant (total count: 2,406; mean: 343.7 ± 106.7) was present in NREM II (total: 2,059; mean: 249.1 ± 100.2, 84.9 ± 3.6%) and SWS (total: 347; mean: 49.5 ± 12.8, 15.0 ± 3.6%) stages, and absent from the rest of sleep and wakefulness. While all 14&6/sec positive spikes correlated with barques, only 44.7 ± 6.1% of barques manifested as 14&6/sec positive spikes.

CONCLUSIONS

Hippocampal barques are predominant in NREM II and SWS, and tend to increase their presence during SWS. Their scalp manifestation as 14&6/sec positive spikes is confounded by wakefulness, REM and NREM I stages, and "masked" by the co-occurrence of NREM II and SWS slow waves, and overlapping reactive micro-arousal elements.

SIGNIFICANCE

Our study highlighted the overnight profile of hippocampal barques, in relation to the respective profile of their scalp manifestation, the 14&6/sec positive spikes variant.

Human anterior thalamic stimulation evoked cortical potentials align with intrinsic functional connectivity

Characterizing human thalamocortical network is fundamental for understanding a vast array of human behaviors since the thalamus plays a central role in cortico-subcortical communication. Over the past few decades, advances in functional magnetic resonance imaging have allowed for spatial mapping of intrinsic resting-state functional connectivity (RSFC) between both cortical regions and in cortico-subcortical networks. Despite these advances, identifying the electrophysiological basis of human thalamocortical network architecture remains challenging. By leveraging stereoelectroencephalography electrodes temporarily implanted into distributed cortical regions and the anterior nucleus of the thalamus (ANT) of 10 patients with refractory focal epilepsy, we tested whether ANT stimulation evoked cortical potentials align with RSFC from the stimulation site, derived from a normative functional connectome (n=1000). Our study identifies spatial convergence of ANT stimulation evoked cortical potentials and normative RSFC. Other than connections to the Papez circuit, the ANT was found to be closely connected to several distinct higher-order association cortices, including the precuneus, angular gyrus, dorsal lateral prefrontal cortex, and anterior insula. Remarkably, we found that the spatial distribution and magnitude of cortical-evoked responses to single-pulse electrical stimulation of the ANT aligned with the spatial pattern and strength of normative RSFC of the stimulation site. The present study provides electrophysiological evidence that stimulation evoked electrical activity flows along intrinsic brain networks connected on a thalamocortical level.

Evaluation of Pharmacology and Pathophysiology Knowledge of Epilepsy among Senior Pharmacy Students: A Single Center Experience

Background and Objectives: Epilepsy is a chronic disease that causes substantial morbidity and mortality. Pharmacists represent an integral role in managing patients with epilepsy. The aim of this study was to evaluate the level of knowledge about the pharmacology and pathophysiology of epilepsy among senior pharmacy students. Materials and Methods: Cross-sectional study using a designed questionnaire to measure the pharmacological and physiological knowledge of senior pharmacy students regarding epilepsy who are studying at Umm Al-Qura University, Makkah, Saudi Arabia, from August to October 2022. Results: A total of 211 senior clinical pharmacy students responded to the questionnaire. The majority of the respondents were 4th year pharmacy students. The numbers of female and male participants were equal (106 and 105 students, respectively). The participants represented an acceptable level of knowledge about the pathophysiology aspects of epilepsy, with a mean total score of 6.22 ± 1.9 out of a maximum score of 10. The respondents reported that epilepsy could be due to genetic predisposition combined with environmental conditions (80.1%) or brain stroke (17.1%). Regarding the respondent knowledge about the pharmacology of epilepsy, the total score was 4.6 ± 2.1 (maximum attainable score: 9). Conclusions: The majority of pharmacy students had knowledge about the pathophysiology concept of the disease; however, low knowledge was shown by the respondents regarding the pharmacology of epilepsy. Thus, there is a need to identify better strategies to improve students' education.

Localized chemogenetic silencing of inhibitory neurons: a novel mouse model of focal cortical epileptic activity

Focal cortical epilepsies are frequently refractory to available anticonvulsant drug therapies. One key factor contributing to this state is the limited availability of animal models that allow to reliably study focal cortical seizures and how they recruit surrounding brain areas in vivo. In this study, we selectively expressed the inhibitory chemogenetic receptor, hM4D, in GABAergic neurons in focal cortical areas using viral gene transfer. GABAergic silencing using Clozapine-N-Oxide (CNO) demonstrated reliable induction of local epileptiform events in the electroencephalogram signal of awake freely moving mice. Anesthetized mice experiments showed consistent induction of focal epileptiform-events in both the barrel cortex (BC) and the medial prefrontal cortex (mPFC), accompanied by high-frequency oscillations, a known characteristic of human seizures. Epileptiform-events showed propagation indication with favored propagation pathways: from the BC on 1 hemisphere to its counterpart and from the BC to the mPFC, but not vice-versa. Lastly, sensory whisker-pad stimulation evoked BC epileptiform events post-CNO, highlighting the potential use of this model in studying sensory-evoked seizures. Combined, our results show that targeted chemogenetic inhibition of GABAergic neurons using hM4D can serve as a novel, versatile, and reliable model of focal cortical epileptic activity suitable for systematically studying cortical ictogenesis in different cortical areas.

Multiple mechanisms shape the relationship between pathway and duration of focal seizures

A seizure's electrographic dynamics are characterized by its spatiotemporal evolution, also termed dynamical 'pathway', and the time it takes to complete that pathway, which results in the seizure's duration. Both seizure pathways and durations have been shown to vary within the same patient. However, it is unclear whether seizures following the same pathway will have the same duration or if these features can vary independently. We compared within-subject variability in these seizure features using (i) epilepsy monitoring unit intracranial EEG (iEEG) recordings of 31 patients (mean: 6.7 days, 16.5 seizures/subject), (ii) NeuroVista chronic iEEG recordings of 10 patients (mean: 521.2 days, 252.6 seizures/subject) and (iii) chronic iEEG recordings of three dogs with focal-onset seizures (mean: 324.4 days, 62.3 seizures/subject). While the strength of the relationship between seizure pathways and durations was highly subject-specific, in most subjects, changes in seizure pathways were only weakly to moderately associated with differences in seizure durations. The relationship between seizure pathways and durations was strengthened by seizures that were 'truncated' versions, both in pathway and duration, of other seizures. However, the relationship was weakened by seizures that had a common pathway, but different durations ('elasticity'), or had similar durations, but followed different pathways ('semblance'). Even in subjects with distinct populations of short and long seizures, seizure durations were not a reliable indicator of different seizure pathways. These findings suggest that seizure pathways and durations are modulated by multiple different mechanisms. Uncovering such mechanisms may reveal novel therapeutic targets for reducing seizure duration and severity.

Human interictal epileptiform discharges are bidirectional traveling waves echoing ictal discharges

Interictal epileptiform discharges (IEDs), also known as interictal spikes, are large intermittent electrophysiological events observed between seizures in patients with epilepsy. Although they occur far more often than seizures, IEDs are less studied, and their relationship to seizures remains unclear. To better understand this relationship, we examined multi-day recordings of microelectrode arrays implanted in human epilepsy patients, allowing us to precisely observe the spatiotemporal propagation of IEDs, spontaneous seizures, and how they relate. These recordings showed that the majority of IEDs are traveling waves, traversing the same path as ictal discharges during seizures, and with a fixed direction relative to seizure propagation. Moreover, the majority of IEDs, like ictal discharges, were bidirectional, with one predominant and a second, less frequent antipodal direction. These results reveal a fundamental spatiotemporal similarity between IEDs and ictal discharges. These results also imply that most IEDs arise in brain tissue outside the site of seizure onset and propagate toward it, indicating that the propagation of IEDs provides useful information for localizing the seizure focus.

Ictal wavefront propagation in slices and simulations with conductance-based refractory density model

The mechanisms determining ictal discharge (ID) propagation are still not clear. In the present study, we aimed to examine these mechanisms in animal and mathematical models of epileptiform activity. Using double-patch and extracellular potassium ion concentration recordings in rat hippocampal-cortical slices, we observed that IDs moved at a speed of about 1 mm/s or less. The mechanisms of such slow propagation have been studied with a mathematical, conductance-based refractory density (CBRD) model that describes the GABA- and glutamatergic neuronal populations’ interactions and ion dynamics in brain tissue. The modeling study reveals two main factors triggerring IDs: (i) increased interneuronal activity leading to chloride ion accumulation and a consequent depolarizing GABAergic effect and (ii) the elevation of extracellular potassium ion concentration. The local synaptic transmission followed by local potassium ion extrusion and GABA receptor-mediated chloride ion accumulation underlies the ID wavefront’s propagation. In contrast, potassium ion diffusion in the extracellular space is slower and does not affect ID’s speed. The short discharges, constituting the ID, propagate much faster than the ID front. The accumulation of sodium ions inside neurons due to their hyperactivity and glutamatergic currents boosts the Na⁺/K⁺ pump, which terminates the ID. Knowledge of the mechanism of ID generation and propagation contributes to the development of new treatments against epilepsy.

During an epileptic seizure, neuronal excitation spreads across the brain tissue and is accompanied by significant changes in ionic concentrations. Ictal discharge front spreads at low speeds, less than 1 mm/s. Mechanisms underlying this phenomenon are not yet well understood. We study these mechanisms using electrophysiological recordings in brain slices and computer simulations. Our detailed biophysical model describing neuronal populations’ interaction, spatial propagation, and ionic dynamics reproduces the generation and propagation of spontaneously repeating ictal discharges. The simulations are consistent with our recordings of the electrical activity and the extracellular potassium ion concentration. We distinguished between the two alternative mechanisms of the ictal wavefront propagation: (i) the diffusion of potassium ions released from excited neurons, which depolarizes distant neurons and thus supports excitation, and (ii) the axonal spread of excitation followed by the local extracellular potassium ion accumulation that supports the excitation. Our simulations provide evidence in favor of the latter mechanism. Our experiment-based modeling contributes to a mathematical description of brain tissue functioning and potentially contributes to developing new treatments against epilepsy.

Simultaneous two-photon imaging of action potentials and subthreshold inputs in vivo

To better understand the input-output computations of neuronal populations, we developed ArcLight-ST, a genetically-encoded voltage indicator, to specifically measure subthreshold membrane potentials. We combined two-photon imaging of voltage and calcium, and successfully discriminated subthreshold inputs and spikes with cellular resolution in vivo. We demonstrate the utility of the method by mapping epileptic seizures progression through cortical circuits, revealing divergent sub- and suprathreshold dynamics within compartmentalized epileptic micronetworks. Two-photon, two-color imaging of calcium and voltage enables mapping of inputs and outputs in neuronal populations in living animals.

Rapidly spreading seizures arise from large-scale functional brain networks in focal epilepsy

It remains unclear whether epileptogenic networks in focal epilepsy develop on physiological networks. This work aimed to explore the association between the rapid spread of ictal fast activity (IFA), a proposed biomarker for epileptogenic networks, and the functional connectivity or networks of healthy subjects. We reviewed 45 patients with focal epilepsy who underwent electrocorticographic (ECoG) recordings to identify the patients showing the rapid spread of IFA. IFA power was quantified as normalized beta-gamma band power. Using published resting-state functional magnetic resonance imaging databases, we estimated resting-state functional connectivity of healthy subjects (RSFC-HS) and resting-state networks of healthy subjects (RSNs-HS) at the locations corresponding to the patients' electrodes. We predicted the IFA power of each electrode based on RSFC-HS between electrode locations (RSFC-HS-based prediction) using a recently developed method, termed activity flow mapping. RSNs-HS were identified using seed-based and atlas-based methods. We compared IFA power with RSFC-HS-based prediction or RSNs-HS using non-parametric correlation coefficients. RSFC and seed-based RSNs of each patient (RSFC-PT and seed-based RSNs-PT) were also estimated using interictal ECoG data and compared with IFA power in the same way as RSFC-HS and seed-based RSNs-HS. Spatial autocorrelation-preserving randomization tests were performed for significance testing. Nine patients met the inclusion criteria. None of the patients had reflex seizures. Six patients showed pathological evidence of a structural etiology. In total, we analyzed 49 seizures (2-13 seizures per patient). We observed significant correlations between IFA power and RSFC-HS-based prediction, seed-based RSNs-HS, or atlas-based RSNs-HS in 28 (57.1%), 21 (42.9%), and 28 (57.1%) seizures, respectively. Thirty-two (65.3%) seizures showed a significant correlation with either seed-based or atlas-based RSNs-HS, but this ratio varied across patients: 27 (93.1%) of 29 seizures in six patients correlated with either of them. Among atlas-based RSNs-HS, correlated RSNs-HS with IFA power included the default mode, control, dorsal attention, somatomotor, and temporal-parietal networks. We could not obtain RSFC-PT and RSNs-PT in one patient due to frequent interictal epileptiform discharges. In the remaining eight patients, most of the seizures showed significant correlations between IFA power and RSFC-PT-based prediction or seed-based RSNs-PT. Our study provides evidence that the rapid spread of IFA in focal epilepsy can arise from physiological RSNs. This finding suggests an overlap between epileptogenic and functional networks, which may explain why functional networks in patients with focal epilepsy frequently disrupt.

Travelling waves reveal a dynamic seizure source in human focal epilepsy

Treatment of patients with drug-resistant focal epilepsy relies upon accurate seizure localization. Ictal activity captured by intracranial EEG has traditionally been interpreted to suggest that the underlying cortex is actively involved in seizures. Here, we hypothesize that such activity instead reflects propagated activity from a relatively focal seizure source, even during later time points when ictal activity is more widespread. We used the time differences observed between ictal discharges in adjacent electrodes to estimate the location of the hypothesized focal source and demonstrated that the seizure source, localized in this manner, closely matches the clinically and neurophysiologically determined brain region giving rise to seizures. Moreover, we determined this focal source to be a dynamic entity that moves and evolves over the time course of a seizure. Our results offer an interpretation of ictal activity observed by intracranial EEG that challenges the traditional conceptualization of the seizure source.

Evidence for spreading seizure as a cause of theta-alpha activity electrographic pattern in stereo-EEG seizure recordings

Intracranial electroencephalography is a standard tool in clinical evaluation of patients with focal epilepsy. Various early electrographic seizure patterns differing in frequency, amplitude, and waveform of the oscillations are observed. The pattern most common in the areas of seizure propagation is the so-called theta-alpha activity (TAA), whose defining features are oscillations in the θ - α range and gradually increasing amplitude. A deeper understanding of the mechanism underlying the generation of the TAA pattern is however lacking. In this work we evaluate the hypothesis that the TAA patterns are caused by seizures spreading across the cortex. To do so, we perform simulations of seizure dynamics on detailed patient-derived cortical surfaces using the spreading seizure model as well as reference models with one or two homogeneous sources. We then detect the occurrences of the TAA patterns both in the simulated stereo-electroencephalographic signals and in the signals of recorded epileptic seizures from a cohort of fifty patients, and we compare the features of the groups of detected TAA patterns to assess the plausibility of the different models. Our results show that spreading seizure hypothesis is qualitatively consistent with the evidence available in the seizure recordings, and it can explain the features of the detected TAA groups best among the examined models.

Multimodal in vivo recording using transparent graphene microelectrodes illuminates spatiotemporal seizure dynamics at the microscale

Neurological disorders such as epilepsy arise from disrupted brain networks. Our capacity to treat these disorders is limited by our inability to map these networks at sufficient temporal and spatial scales to target interventions. Current best techniques either sample broad areas at low temporal resolution (e.g. calcium imaging) or record from discrete regions at high temporal resolution (e.g. electrophysiology). This limitation hampers our ability to understand and intervene in aberrations of network dynamics. Here we present a technique to map the onset and spatiotemporal spread of acute epileptic seizures in vivo by simultaneously recording high bandwidth microelectrocorticography and calcium fluorescence using transparent graphene microelectrode arrays. We integrate dynamic data features from both modalities using non-negative matrix factorization to identify sequential spatiotemporal patterns of seizure onset and evolution, revealing how the temporal progression of ictal electrophysiology is linked to the spatial evolution of the recruited seizure core. This integrated analysis of multimodal data reveals otherwise hidden state transitions in the spatial and temporal progression of acute seizures. The techniques demonstrated here may enable future targeted therapeutic interventions and novel spatially embedded models of local circuit dynamics during seizure onset and evolution.

Early seizure spread and epilepsy surgery: A systematic review

OBJECTIVE

A fundamental question in epilepsy surgery is how to delineate the margins of cortex that must be resected to result in seizure freedom. Whether and which areas showing seizure activity early in ictus must be removed to avoid postoperative recurrence of seizures is an area of ongoing research. Seizure spread dynamics in the initial seconds of ictus are often correlated with postoperative outcome; there is neither a consensus definition of early spread nor a concise summary of the existing literature linking seizure spread to postsurgical seizure outcomes. The present study is intended to summarize the literature that links seizure spread to postoperative seizure outcome and to provide a framework for quantitative assessment of early seizure spread.

METHODS

A systematic review was carried out according to Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines. A Medline search identified clinical studies reporting data on seizure spread measured by intracranial electrodes, having at least 10 subjects and reporting at least 1-year postoperative outcome in the English literature from 1990 to 2019. Studies were evaluated regarding support for a primary hypothesis: Areas of early seizure spread represent cortex with seizure-generating potential.

RESULTS

The search yielded 4562 studies: 15 studies met inclusion criteria and 7 studies supported the primary hypothesis. The methods and metrics used to describe seizure spread were heterogenous. The timeframe of seizure spread associated with seizure outcome ranged from 1-14 seconds, with large, well-designed, retrospective studies pointing to 3-10 seconds as most likely to provide meaningful correlates of postoperative seizure freedom.

SIGNIFICANCE

The complex correlation between electrophysiologic seizure spread and the potential for seizure generation needs further elucidation. Prospective cohort studies or trials are needed to evaluate epilepsy surgery targeting cortex involved in the first 3-10 seconds of ictus.

Domino-like transient dynamics at seizure onset in epilepsy

The International League Against Epilepsy (ILAE) groups seizures into "focal", "generalized" and "unknown" based on whether the seizure onset is confined to a brain region in one hemisphere, arises in several brain region simultaneously, or is not known, respectively. This separation fails to account for the rich diversity of clinically and experimentally observed spatiotemporal patterns of seizure onset and even less so for the properties of the brain networks generating them. We consider three different patterns of domino-like seizure onset in Idiopathic Generalized Epilepsy (IGE) and present a novel approach to classification of seizures. To understand how these patterns are generated on networks requires understanding of the relationship between intrinsic node dynamics and coupling between nodes in the presence of noise, which currently is unknown. We investigate this interplay here in the framework of domino-like recruitment across a network. In particular, we use a phenomenological model of seizure onset with heterogeneous coupling and node properties, and show that in combination they generate a range of domino-like onset patterns observed in the IGE seizures. We further explore the individual contribution of heterogeneous node dynamics and coupling by interpreting in-vitro experimental data in which the speed of onset can be chemically modulated. This work contributes to a better understanding of possible drivers for the spatiotemporal patterns observed at seizure onset and may ultimately contribute to a more personalized approach to classification of seizure types in clinical practice.

Feasibility of focal brain cooling for partial epilepsy with secondary generalization: A computational study. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2020;2020:2532-2535. [PMID: 33018522 DOI: 10.1109/embc44109.2020.9175157] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Experiments with animal models of epilepsy have consistently shown that focal cooling of epilepsy-induced brain region reversibly suppresses or terminates epileptic discharge activity. Recently, we formulated a physiologically plausible temperature dependence in a neural mass model that can reproduce the effect of focal cooling on epileptic discharge activity. This can be used to implement a temperature control in an implantable cooling device for thermal neuromodulation of the epileptogenic zone in patients with partial epilepsy when seizure activity is detected. However, there have been no experiments that looked into the effect of focal cooling in animal models of epilepsy with secondary generalization in which the seizure activity spreads from the pathologic region to other regions of the brain. Using the temperature-dependent neural mass model and a physiological coupling model, we show that focal cooling stops the propagation of low-frequency discharge activity; on the other hand, it increases the amount of coupling required to propagate high-frequency discharge activity. Moreover, discharge activities that are propagated with cooling are lower in both magnitude and frequency compared to those propagated without cooling. These results suggest the feasibility of focal cooling as an effective alternative therapeutic treatment for medically intractable partial epilepsy even with secondary generalization.Clinical Relevance- The computational study establishes focal cooling of the brain region with partial epilepsy not only suppresses epileptic discharges but can also prevent its generalization to other brain regions.

Robust dynamic community detection with applications to human brain functional networks

While current technology permits inference of dynamic brain networks over long time periods at high temporal resolution, the detailed structure of dynamic network communities during human seizures remains poorly understood. We introduce a new methodology that addresses critical aspects unique to the analysis of dynamic functional networks inferred from noisy data. We propose a dynamic plex percolation method (DPPM) that is robust to edge noise, and yields well-defined spatiotemporal communities that span forward and backwards in time. We show in simulation that DPPM outperforms existing methods in accurately capturing certain stereotypical dynamic community behaviors in noisy situations. We then illustrate the ability of this method to track dynamic community organization during human seizures, using invasive brain voltage recordings at seizure onset. We conjecture that application of this method will yield new targets for surgical treatment of epilepsy, and more generally could provide new insights in other network neuroscience applications.

Understanding how brain networks evolve in time remains a challenge, with the potential for significant impact to human health and disease. Here, the authors introduce a new methodology to track dynamic functional networks that is robust to edge noise, and yields well-defined spatiotemporal communities that span forward and backwards in time.

Seizure pathways change on circadian and slower timescales in individual patients with focal epilepsy

Personalized medicine requires that treatments adapt to not only the patient but also changing factors within each individual. Although epilepsy is a dynamic disorder characterized by pathological fluctuations in brain state, surprisingly little is known about whether and how seizures vary in the same patient. We quantitatively compared within-patient seizure network evolutions using intracranial electroencephalographic (iEEG) recordings of over 500 seizures from 31 patients with focal epilepsy (mean 16.5 seizures per patient). In all patients, we found variability in seizure paths through the space of possible network dynamics. Seizures with similar pathways tended to occur closer together in time, and a simple model suggested that seizure pathways change on circadian and/or slower timescales in the majority of patients. These temporal relationships occurred independent of whether the patient underwent antiepileptic medication reduction. Our results suggest that various modulatory processes, operating at different timescales, shape within-patient seizure evolutions, leading to variable seizure pathways that may require tailored treatment approaches.

Association of Seizure Spread With Surgical Failure in Epilepsy

Importance

Seizures recur in as many as half of patients who undergo surgery for drug-resistant temporal lobe epilepsy (TLE). Understanding why TLE is resistant to surgery in some patients may reveal insights into epileptogenic networks and direct new therapies to improve outcomes.

Objective

To characterize features of surgically refractory TLE.

Design, Setting, and Participants

Medical records from a comprehensive epilepsy center were retrospectively reviewed for 131 patients who received a standard anteromedial temporal resection by a single surgeon from January 1, 2000, to December 31, 2015. Thirteen patients were excluded for having less than 1 year of follow-up. Patients at the highest risk for seizure recurrence were identified. Intracranial electroencephalogram (iEEG) analyses generated 3-dimensional seizure spread representations and quantified rapid seizure spread. The final analyses of seizure outcome and follow-up data were performed in June 2017.

Main Outcomes and Measures

The Engel class seizure outcome following surgery was evaluated for all patients, defining seizure recurrence as Engel class II or greater. Intracranial recordings of neocortical grids/strips and depth electrodes were analyzed visually for seizure spread. Fast β power was projected onto reconstructions of patients' brain magnetic resonance imaging scans to visualize spread patterns and was quantified to compare power within vs outside resective margins.

Results

Of 118 patients with 1 year of follow-up or more (mean [SD], 6.5 [4.6] years), 66 (55.9%) were women and 52 (44.1%) were men (median age, 39 years [range, 4-66 years]). The cumulative probability of continuous Engel class I seizure freedom since surgery at postoperative year 10 and afterward was 65.6%, with 92% of recurrences in years 1 to 3. Multivariable statistical analyses found that the selection for iEEG study was the most reliable predictor of seizure recurrence, with a mixed-effects model estimating that the Engel score in the iEEG cohort was higher by a mean (SD) of 1.1 (0.33) (P = .001). In patients with iEEG results, rapid seizure spread in less than 10 seconds was associated with recurrence (hazard ratio, 5.99; 95% CI, 1.7-21.1; P < .01). In the first 10 seconds of seizures, fast β power activity outside the resective margins in the lateral temporal cortex was significantly greater in patients whose seizures recurred compared with patients who were seizure-free (mean [SEM], 137.5% [16.8%] vs 93.4% [4.6%]; P < .05).

Conclusions and significance

Rapid seizure spread outside anteromedial temporal resection resective margins plays a significant role in the surgical failure of drug-resistant TLE. Seizure control after epilepsy surgery might be improved by investigating areas of early spread as candidates for resection or neuromodulation.

Unified analysis of global and focal aspects of absence epilepsy via neural field theory of the corticothalamic system

Absence epilepsy is characterized by a sudden paroxysmal loss of consciousness accompanied by oscillatory activity propagating over many brain areas. Although primary generalized absence seizures are supported by the global corticothalamic system, converging experimental evidence supports a focal theory of absence epilepsy. Here a physiology-based corticothalamic model is investigated with spatial heterogeneity due to focal epilepsy to unify global and focal aspects of absence epilepsy. Numeric and analytic calculations are employed to investigate the emergent spatiotemporal dynamics as well as their underlying dynamical mechanisms. They can be categorized into three scenarios: suppressed epilepsy, focal seizures, or generalized seizures, as summarized from a phase diagram vs focal width and characteristic axon range. The corresponding temporal frequencies and spatial extents of cortical waves in generalized seizures and focal seizures agree well with experimental observations of global and focal aspects of absence epilepsy, respectively. The emergence of the spatiotemporal dynamics corresponding to focal seizures provides a biophysical explanation of the temporally higher frequency but spatially more localized cortical waves observed in genetic rat models that display characteristics of human absence epilepsy. Predictions are also presented for further experimental test.

Reliable and Elastic Propagation of Cortical Seizures In Vivo

Mapping the fine-scale neural activity that underlies epilepsy is key to identifying potential control targets of this frequently intractable disease. Yet, the detailed in vivo dynamics of seizure progression in cortical microcircuits remain poorly understood. We combine fast (30-Hz) two-photon calcium imaging with local field potential (LFP) recordings to map, cell by cell, the spread of locally induced (4-AP or picrotoxin) seizures in anesthetized and awake mice. Using single-layer and microprism-assisted multilayer imaging in different cortical areas, we uncover reliable recruitment of local neural populations within and across cortical layers, and we find layer-specific temporal delays, suggesting an initial supra-granular invasion followed by deep-layer recruitment during lateral seizure spread. Intriguingly, despite consistent progression pathways, successive seizures show pronounced temporal variability that critically depends on GABAergic inhibition. We propose an epilepsy circuit model resembling an elastic meshwork, wherein ictal progression faithfully follows preexistent pathways but varies flexibly in time, depending on the local inhibitory restraint.

Predicting the spatiotemporal diversity of seizure propagation and termination in human focal epilepsy

Recent studies have shown that seizures can spread and terminate across brain areas via a rich diversity of spatiotemporal patterns. In particular, while the location of the seizure onset area is usually invariant across seizures in an individual patient, the source of traveling (2–3 Hz) spike-and-wave discharges during seizures can either move with the slower propagating ictal wavefront or remain stationary at the seizure onset area. Furthermore, although many focal seizures terminate synchronously across brain areas, some evolve into distinct ictal clusters and terminate asynchronously. Here, we introduce a unifying perspective based on a new neural field model of epileptic seizure dynamics. Two main mechanisms, the co-existence of wave propagation in excitable media and coupled-oscillator dynamics, together with the interaction of multiple time scales, account for the reported diversity. We confirm our predictions in seizures and tractography data obtained from patients with pharmacologically resistant epilepsy. Our results contribute toward patient-specific seizure modeling.

A major goal of epilepsy research is understanding the spatiotemporal dynamics of seizure. Here, the authors extend the Epileptor neural mass model into a neural field model, in order to provide a unified and patient-specific model of seizure initiation, propagation, and termination.

Slow moving neural source in the epileptic hippocampus can mimic progression of human seizures

Fast and slow neural waves have been observed to propagate in the human brain during seizures. Yet the nature of these waves is difficult to study in a surgical setting. Here, we report an observation of two different traveling waves propagating in the in-vitro epileptic hippocampus at speeds similar to those in the human brain. A fast traveling spike and a slow moving wave were recorded simultaneously with a genetically encoded voltage sensitive fluorescent protein (VSFP Butterfly 1.2) and a high speed camera. The results of this study indicate that the fast traveling spike is NMDA-sensitive but the slow moving wave is not. Image analysis and model simulation demonstrate that the slow moving wave is moving slowly, generating the fast traveling spike and is, therefore, a moving source of the epileptiform activity. This slow moving wave is associated with a propagating neural calcium wave detected with calcium dye (OGB-1) but is independent of NMDA receptors, not related to ATP release, and much faster than those previously recorded potassium waves. Computer modeling suggests that the slow moving wave can propagate by the ephaptic effect like epileptiform activity. These findings provide an alternative explanation for slow propagation seizure wavefronts associated with fast propagating spikes.

Individual brain structure and modelling predict seizure propagation

See Lytton (doi:10.1093/awx018) for a scientific commentary on this article.Neural network oscillations are a fundamental mechanism for cognition, perception and consciousness. Consequently, perturbations of network activity play an important role in the pathophysiology of brain disorders. When structural information from non-invasive brain imaging is merged with mathematical modelling, then generative brain network models constitute personalized in silico platforms for the exploration of causal mechanisms of brain function and clinical hypothesis testing. We here demonstrate with the example of drug-resistant epilepsy that patient-specific virtual brain models derived from diffusion magnetic resonance imaging have sufficient predictive power to improve diagnosis and surgery outcome. In partial epilepsy, seizures originate in a local network, the so-called epileptogenic zone, before recruiting other close or distant brain regions. We create personalized large-scale brain networks for 15 patients and simulate the individual seizure propagation patterns. Model validation is performed against the presurgical stereotactic electroencephalography data and the standard-of-care clinical evaluation. We demonstrate that the individual brain models account for the patient seizure propagation patterns, explain the variability in postsurgical success, but do not reliably augment with the use of patient-specific connectivity. Our results show that connectome-based brain network models have the capacity to explain changes in the organization of brain activity as observed in some brain disorders, thus opening up avenues towards discovery of novel clinical interventions.

Management of diffuse low-grade gliomas in adults - use of molecular diagnostics

Diffuse WHO grade II gliomas are histologically and genetically heterogeneous. The 2016 WHO classification redefines grade II gliomas with respect to morphological and molecular tumour alterations: grade II oligodendrogliomas are defined by the presence of whole-arm codeletion in chromosomal arms 1p/19q, whereas isocitrate dehydrogenase (IDH) mutations define subclasses of astrocytoma. Although histological grade remains useful, the prognoses of patients with glioma are more tightly associated with molecular alterations than with grade, and chromosomal and gene array technologies are becoming increasingly beneficial in understanding tumour genetic heterogeneity. The indolent nature of the disease often creates subtle neurological symptoms that can be overlooked or misunderstood, resulting in delayed diagnosis. Seizures often herald the diagnosis, especially in patients who have IDH mutations, which are associated with an increased production of 2-hydroxyglutarate. Treatment paradigms have shifted, owing to new diagnostic criteria and new clinical trial evidence. Patients benefit more from chemoradiation than radiation alone, especially those with tumour IDH1 Arg132His mutations; gross total resection of the tumour, including tumours with IDH mutations, is associated with prolonged survival. Initial observation remains appropriate in patients whose rate of disease growth is not yet completely defined; such patients could include those with completely resected disease and those with 1p/19q codeleted tumours.

Mechanisms underlying different onset patterns of focal seizures

Focal seizures are episodes of pathological brain activity that appear to arise from a localised area of the brain. The onset patterns of focal seizure activity have been studied intensively, and they have largely been distinguished into two types-low amplitude fast oscillations (LAF), or high amplitude spikes (HAS). Here we explore whether these two patterns arise from fundamentally different mechanisms. Here, we use a previously established computational model of neocortical tissue, and validate it as an adequate model using clinical recordings of focal seizures. We then reproduce the two onset patterns in their most defining properties and investigate the possible mechanisms underlying the different focal seizure onset patterns in the model. We show that the two patterns are associated with different mechanisms at the spatial scale of a single ECoG electrode. The LAF onset is initiated by independent patches of localised activity, which slowly invade the surrounding tissue and coalesce over time. In contrast, the HAS onset is a global, systemic transition to a coexisting seizure state triggered by a local event. We find that such a global transition is enabled by an increase in the excitability of the "healthy" surrounding tissue, which by itself does not generate seizures, but can support seizure activity when incited. In our simulations, the difference in surrounding tissue excitability also offers a simple explanation of the clinically reported difference in surgical outcomes. Finally, we demonstrate in the model how changes in tissue excitability could be elucidated, in principle, using active stimulation. Taken together, our modelling results suggest that the excitability of the tissue surrounding the seizure core may play a determining role in the seizure onset pattern, as well as in the surgical outcome.

Human seizures couple across spatial scales through travelling wave dynamics

Epilepsy—the propensity toward recurrent, unprovoked seizures—is a devastating disease affecting 65 million people worldwide. Understanding and treating this disease remains a challenge, as seizures manifest through mechanisms and features that span spatial and temporal scales. Here we address this challenge through the analysis and modelling of human brain voltage activity recorded simultaneously across microscopic and macroscopic spatial scales. We show that during seizure large-scale neural populations spanning centimetres of cortex coordinate with small neural groups spanning cortical columns, and provide evidence that rapidly propagating waves of activity underlie this increased inter-scale coupling. We develop a corresponding computational model to propose specific mechanisms—namely, the effects of an increased extracellular potassium concentration diffusing in space—that support the observed spatiotemporal dynamics. Understanding the multi-scale, spatiotemporal dynamics of human seizures—and connecting these dynamics to specific biological mechanisms—promises new insights to treat this devastating disease.

The authors record both local and long-range neural activity during human epileptic seizures to study the underlying multi-scale dynamics. They find that coupling of activity across spatial scales increases during seizures through propagating waves that are fit by a model that combines neural activity and potassium concentration dynamics.

Spatiotemporal Mapping of Interictal Spike Propagation: A Novel Methodology Applied to Pediatric Intracranial EEG Recordings

Synchronized cortical activity is implicated in both normative cognitive functioning and many neurologic disorders. For epilepsy patients with intractable seizures, irregular synchronization within the epileptogenic zone (EZ) is believed to provide the network substrate through which seizures initiate and propagate. Mapping the EZ prior to epilepsy surgery is critical for detecting seizure networks in order to achieve postsurgical seizure control. However, automated techniques for characterizing epileptic networks have yet to gain traction in the clinical setting. Recent advances in signal processing and spike detection have made it possible to examine the spatiotemporal propagation of interictal spike discharges across the epileptic cortex. In this study, we present a novel methodology for detecting, extracting, and visualizing spike propagation and demonstrate its potential utility as a biomarker for the EZ. Eighteen presurgical intracranial EEG recordings were obtained from pediatric patients ultimately experiencing favorable (i.e., seizure-free, n = 9) or unfavorable (i.e., seizure-persistent, n = 9) surgical outcomes. Novel algorithms were applied to extract multichannel spike discharges and visualize their spatiotemporal propagation. Quantitative analysis of spike propagation was performed using trajectory clustering and spatial autocorrelation techniques. Comparison of interictal propagation patterns revealed an increase in trajectory organization (i.e., spatial autocorrelation) among Sz-Free patients compared with Sz-Persist patients. The pathophysiological basis and clinical implications of these findings are considered.

Secondary generalization of focal-onset seizures: examining the relationship between seizure propagation and epilepsy surgery outcome

Surgical intervention often fails to achieve seizure-free results in patients with intractable epilepsy. Identifying features of the epileptic brain that dispose certain patients to unfavorable outcomes is critical for improving surgical candidacy assessments. Recent research by Martinet, Ahmad, Lepage, Cash, and Kramer (J Neurosci 35: 9477-9490, 2015) suggests that pathways of secondary seizure generalization distinguish patients with favorable (i.e., seizure free) vs. unfavorable (i.e., seizure persistent) surgical outcomes, lending insights into the network mechanisms of epilepsy surgery failure.