A critical role for network structure in seizure onset: a computational modeling approach

Treatment effects in epilepsy: a mathematical framework for understanding response over time

Epilepsy is a neurological disorder characterized by recurrent seizures, affecting over 65 million people worldwide. Treatment typically commences with the use of anti-seizure medications, including both mono- and poly-therapy. Should these fail, more invasive therapies such as surgery, electrical stimulation and focal drug delivery are often considered in an attempt to render the person seizure free. Although a significant portion ultimately benefit from these treatment options, treatment responses often fluctuate over time. The physiological mechanisms underlying these temporal variations are poorly understood, making prognosis a significant challenge when treating epilepsy. Here we use a dynamic network model of seizure transition to understand how seizure propensity may vary over time as a consequence of changes in excitability. Through computer simulations, we explore the relationship between the impact of treatment on dynamic network properties and their vulnerability over time that permit a return to states of high seizure propensity. For small networks we show vulnerability can be fully characterised by the size of the first transitive component (FTC). For larger networks, we find measures of network efficiency, incoherence and heterogeneity (degree variance) correlate with robustness of networks to increasing excitability. These results provide a set of potential prognostic markers for therapeutic interventions in epilepsy. Such markers could be used to support the development of personalized treatment strategies, ultimately contributing to understanding of long-term seizure freedom.

Personalized dynamic network models of the human brain as a future tool for planning and optimizing epilepsy therapy

Epilepsy is a common neurological disorder, with one third of patients not responding to currently available antiepileptic drugs. The proportion of pharmacoresistant epilepsies has remained unchanged for many decades. To cure epilepsy and control seizures requires a paradigm shift in the development of new approaches to epilepsy diagnosis and treatment. Contemporary medicine has benefited from the exponential growth of computational modeling, and the application of network dynamics theory to understanding and treating human brain disorders. In epilepsy, the introduction of these approaches has led to personalized epileptic network modeling that can explore the patient's seizure genesis and predict the functional impact of resection on its individual network's propensity to seize. The application of the dynamic systems approach to neurostimulation therapy of epilepsy allows designing stimulation strategies that consider the patient's seizure dynamics and long-term fluctuations in the stability of their epileptic networks. In this article, we review, in a nontechnical fashion suitable for a broad neuroscientific audience, recent progress in personalized dynamic brain network modeling that is shaping the future approach to the diagnosis and treatment of epilepsy.

Development of MAPT S305 mutation models exhibiting elevated 4R tau expression, resulting in altered neuronal and astrocytic function

Due to the importance of 4R tau in the pathogenicity of primary tauopathies, it has been challenging to model these diseases in iPSC-derived neurons, which express very low levels of 4R tau. To address this problem we have developed a panel of isogenic iPSC lines carrying the MAPT splice-site mutations S305S, S305I or S305N, derived from four different donors. All three mutations significantly increased the proportion of 4R tau expression in iPSC-neurons and astrocytes, with up to 80% 4R transcripts in S305N neurons from as early as 4 weeks of differentiation. Transcriptomic and functional analyses of S305 mutant neurons revealed shared disruption in glutamate signaling and synaptic maturity, but divergent effects on mitochondrial bioenergetics. In iPSC-astrocytes, S305 mutations induced lysosomal disruption and inflammation and exacerbated internalization of exogenous tau that may be a precursor to the glial pathologies observed in many tauopathies. In conclusion, we present a novel panel of human iPSC lines that express unprecedented levels of 4R tau in neurons and astrocytes. These lines recapitulate previously characterized tauopathy-relevant phenotypes, but also highlight functional differences between the wild type 4R and mutant 4R proteins. We also highlight the functional importance of MAPT expression in astrocytes. These lines will be highly beneficial to tauopathy researchers enabling a more complete understanding of the pathogenic mechanisms underlying 4R tauopathies across different cell types.

Classification of human chronotype based on fMRI network-based statistics

Chronotype-the relationship between the internal circadian physiology of an individual and the external 24-h light-dark cycle-is increasingly implicated in mental health and cognition. Individuals presenting with a late chronotype have an increased likelihood of developing depression, and can display reduced cognitive performance during the societal 9-5 day. However, the interplay between physiological rhythms and the brain networks that underpin cognition and mental health is not well-understood. To address this issue, we use rs-fMRI collected from 16 people with an early chronotype and 22 people with a late chronotype over three scanning sessions. We develop a classification framework utilizing the Network Based-Statistic methodology, to understand if differentiable information about chronotype is embedded in functional brain networks and how this changes throughout the day. We find evidence of subnetworks throughout the day that differ between extreme chronotypes such that high accuracy can occur, describe rigorous threshold criteria for achieving 97.3% accuracy in the Evening and investigate how the same conditions hinder accuracy for other scanning sessions. Revealing differences in functional brain networks based on extreme chronotype suggests future avenues of research that may ultimately better characterize the relationship between internal physiology, external perturbations, brain networks, and disease.

The distribution and heterogeneity of excitability in focal epileptic network potentially contribute to the seizure propagation

Introduction

Existing dynamical models can explain the transmigration mechanisms involved in seizures but are limited to a single modality. Combining models with networks can reproduce scaled epileptic dynamics. And the structure and coupling interactions of the network, as well as the heterogeneity of both the node and network activities, may influence the final state of the network model.

Methods

We built a fully connected network with focal nodes prominently interacting and established a timescale separated epileptic network model. The factors affecting epileptic network seizure were explored by varying the connectivity patterns of focal network nodes and modulating the distribution of network excitability.

Results

The whole brain network topology as the brain activity foundation affects the consistent delayed clustering seizure propagation. In addition, the network size and distribution heterogeneity of the focal excitatory nodes can influence seizure frequency. With the increasing of the network size and averaged excitability level of focal network, the seizure period decreases. In contrast, the larger heterogeneity of excitability for focal network nodes can lower the functional activity level (average degree) of focal network. There are also subtle effects of focal network topologies (connection patterns of excitatory nodes) that cannot be ignored along with non-focal nodes.

Discussion

Unraveling the role of excitatory factors in seizure onset and propagation can be used to understand the dynamic mechanisms and neuromodulation of epilepsy, with profound implications for the treatment of epilepsy and even for the understanding of the brain.

The role of additive and diffusive coupling on the dynamics of neural populations

Dynamical models consisting of networks of neural masses commonly assume that the interactions between neural populations are via additive or diffusive coupling. When using the additive coupling, a population's activity is affected by the sum of the activities of neighbouring populations. In contrast, when using the diffusive coupling a neural population is affected by the sum of the differences between its activity and the activity of its neighbours. These two coupling functions have been used interchangeably for similar applications. In this study, we show that the choice of coupling can lead to strikingly different brain network dynamics. We focus on a phenomenological model of seizure transitions that has been used both with additive and diffusive coupling in the literature. We consider small networks with two and three nodes, as well as large random and scale-free networks with 64 nodes. We further assess resting-state functional networks inferred from magnetoencephalography (MEG) from people with juvenile myoclonic epilepsy (JME) and healthy controls. To characterize the seizure dynamics on these networks, we use the escape time, the brain network ictogenicity (BNI) and the node ictogenicity (NI), which are measures of the network's global and local ability to generate seizure activity. Our main result is that the level of ictogenicity of a network is strongly dependent on the coupling function. Overall, we show that networks with additive coupling have a higher propensity to generate seizures than those with diffusive coupling. We find that people with JME have higher additive BNI than controls, which is the hypothesized BNI deviation between groups, while the diffusive BNI provides opposite results. Moreover, we find that the nodes that are more likely to drive seizures in the additive coupling case are more likely to prevent seizures in the diffusive coupling case, and that these features correlate to the node's number of connections. Consequently, previous results in the literature involving such models to interrogate functional or structural brain networks could be highly dependent on the choice of coupling. Our results on the MEG functional networks and evidence from the literature suggest that the additive coupling may be a better modeling choice than the diffusive coupling, at least for BNI and NI studies. Thus, we highlight the need to motivate and validate the choice of coupling in future studies involving network models of brain activity.

DNAJC30 disease-causing gene variants in a large Central European cohort of patients with suspected Leber's hereditary optic neuropathy and optic atrophy

BACKGROUND

Leber's hereditary optic neuropathy (LHON) has been considered a prototypical mitochondriopathy and a textbook example for maternal inheritance linked to certain disease-causing variants in the mitochondrial genome. Recently, an autosomal recessive form of LHON (arLHON) has been described, caused by disease-causing variants in the nuclear encoded gene DNAJC30.

METHODS AND RESULTS

In this study, we screened the DNAJC30 gene in a large Central European cohort of patients with a clinical diagnosis of LHON or other autosomal inherited optic atrophies (OA). We identified likely pathogenic variants in 35/1202 patients, corresponding to a detection rate of 2.9%. The previously described missense variant c.152A>G;p.(Tyr51Cys) accounts for 90% of disease-associated alleles in our cohort and we confirmed a strong founder effect. Furthermore, we identified two novel pathogenic variants in DNAJC30: the nonsense variant c.610G>T;p.(Glu204*) and the in-frame deletion c.230_232del;p.(His77del). Clinical investigation of the patients with arLHON revealed a younger age of onset, a more frequent bilateral onset and an increased clinically relevant recovery compared with LHON associated with disease-causing variants in the mitochondrial DNA.

CONCLUSION

This study expands previous findings on arLHON and emphasises the importance of DNAJC30 in the genetic diagnostics of LHON and OA in European patients.

Circle of Willis variation and outcome after intra-arterial treatment

Background

Intra-arterial treatment (IAT) improves outcomes in acute ischaemic stroke. Presence of collaterals increases likelihood of good outcome. We investigated whether variations in the circle of Willis (CoW) and contributing carotid arteries influence outcome in patients who had a stroke treated with IAT.

Methods

CT angiography data on patients who had an acute stroke treated with IAT were retrospectively collected. CoW was regarded complete if the contralateral A1 segment, anterior communicating artery and ipsilateral posterior communicating artery were fully developed, and the P1 segment was visible. Carotid artery contribution was studied with a self-developed carotid artery score ranging from 0 to 2 depending on the number of arteries supplying the occluded side of the CoW. Good clinical outcome was defined as modified Rankin Score ≤2 and measured at discharge and 3 months. We calculated risk ratios for the relation between completeness of the CoW, carotid score and good outcome, and performed a trend analysis for good outcome according to the carotid score.

Results

126 patients were included for analysis. Patients with a complete and incomplete CoW had a comparable risk for good outcome at discharge and 3 months. A higher carotid score was associated with a higher likelihood of good clinical outcome (p for trend 0.24 at discharge and 0.05 at 3 months).

Conclusion

In patients with acute ischaemic stroke treated with IAT, chances of good clinical outcome tended to improve with number of carotid arteries supplying the cerebral circulation. Completeness of the CoW was not related to clinical outcome.

Heterogeneity of resting-state EEG features in juvenile myoclonic epilepsy and controls

Abnormal EEG features are a hallmark of epilepsy, and abnormal frequency and network features are apparent in EEGs from people with idiopathic generalized epilepsy in both ictal and interictal states. Here, we characterize differences in the resting-state EEG of individuals with juvenile myoclonic epilepsy and assess factors influencing the heterogeneity of EEG features. We collected EEG data from 147 participants with juvenile myoclonic epilepsy through the Biology of Juvenile Myoclonic Epilepsy study. Ninety-five control EEGs were acquired from two independent studies [Chowdhury et al. (2014) and EU-AIMS Longitudinal European Autism Project]. We extracted frequency and functional network-based features from 10 to 20 s epochs of resting-state EEG, including relative power spectral density, peak alpha frequency, network topology measures and brain network ictogenicity: a computational measure of the propensity of networks to generate seizure dynamics. We tested for differences between epilepsy and control EEGs using univariate, multivariable and receiver operating curve analysis. In addition, we explored the heterogeneity of EEG features within and between cohorts by testing for associations with potentially influential factors such as age, sex, epoch length and time, as well as testing for associations with clinical phenotypes including anti-seizure medication, and seizure characteristics in the epilepsy cohort. P-values were corrected for multiple comparisons. Univariate analysis showed significant differences in power spectral density in delta (2-5 Hz) (P = 0.0007, hedges' g = 0.55) and low-alpha (6-9 Hz) (P = 2.9 × 10-8, g = 0.80) frequency bands, peak alpha frequency (P = 0.000007, g = 0.66), functional network mean degree (P = 0.0006, g = 0.48) and brain network ictogenicity (P = 0.00006, g = 0.56) between epilepsy and controls. Since age (P = 0.009) and epoch length (P = 1.7 × 10-8) differed between the two groups and were potential confounders, we controlled for these covariates in multivariable analysis where disparities in EEG features between epilepsy and controls remained. Receiver operating curve analysis showed low-alpha power spectral density was optimal at distinguishing epilepsy from controls, with an area under the curve of 0.72. Lower average normalized clustering coefficient and shorter average normalized path length were associated with poorer seizure control in epilepsy patients. To conclude, individuals with juvenile myoclonic epilepsy have increased power of neural oscillatory activity at low-alpha frequencies, and increased brain network ictogenicity compared with controls, supporting evidence from studies in other epilepsies with considerable external validity. In addition, the impact of confounders on different frequency-based and network-based EEG features observed in this study highlights the need for careful consideration and control of these factors in future EEG research in idiopathic generalized epilepsy particularly for their use as biomarkers.

In vitro ictogenesis is stochastic at the single neuron level

Seizure initiation is the least understood and most disabling element of epilepsy. Studies of ictogenesis require high speed recordings at cellular resolution in the area of seizure onset. However, in vivo seizure onset areas cannot be determined at the level of resolution necessary to enable such studies. To circumvent these challenges, we used novel GCaMP7-based calcium imaging in the organotypic hippocampal slice culture model of post-traumatic epilepsy in mice. Organotypic hippocampal slice cultures generate spontaneous, recurrent seizures in a preparation in which it is feasible to image the activity of the entire network (with no unseen inputs existing). Chronic calcium imaging of the entire hippocampal network, with paired electrophysiology, revealed three patterns of seizure onset: (i) low amplitude fast activity; (ii) sentinel spike; and (iii) spike burst and low amplitude fast activity onset. These patterns recapitulate common features of human seizure onset, including low voltage fast activity and spike discharges. Weeks-long imaging of seizure activity showed a characteristic evolution in onset type and a refinement of the seizure onset zone. Longitudinal tracking of individual neurons revealed that seizure onset is stochastic at the single neuron level, suggesting that seizure initiation activates neurons in non-stereotyped sequences seizure to seizure. This study demonstrates for the first time that transitions to seizure are not initiated by a small number of neuronal 'bad actors' (such as overly connected hub cells), but rather by network changes which enable the onset of pathology among large populations of neurons.

A large-scale brain network mechanism for increased seizure propensity in Alzheimer's disease

People with Alzheimer’s disease (AD) are 6-10 times more likely to develop seizures than the healthy aging population. Leading hypotheses largely consider hyperexcitability of local cortical tissue as primarily responsible for increased seizure prevalence in AD. However, in the general population of people with epilepsy, large-scale brain network organization additionally plays a role in determining seizure likelihood and phenotype. Here, we propose that alterations to large-scale brain network organization seen in AD may contribute to increased seizure likelihood. To test this hypothesis, we combine computational modelling with electrophysiological data using an approach that has proved informative in clinical epilepsy cohorts without AD. EEG was recorded from 21 people with probable AD and 26 healthy controls. At the time of EEG acquisition, all participants were free from seizures. Whole brain functional connectivity derived from source-reconstructed EEG recordings was used to build subject-specific brain network models of seizure transitions. As cortical tissue excitability was increased in the simulations, AD simulations were more likely to transition into seizures than simulations from healthy controls, suggesting an increased group-level probability of developing seizures at a future time for AD participants. We subsequently used the model to assess seizure propensity of different regions across the cortex. We found the most important regions for seizure generation were those typically burdened by amyloid-beta at the early stages of AD, as previously reported by in-vivo and post-mortem staging of amyloid plaques. Analysis of these spatial distributions also give potential insight into mechanisms of increased susceptibility to generalized (as opposed to focal) seizures in AD vs controls. This research suggests avenues for future studies testing patients with seizures, e.g. co-morbid AD/epilepsy patients, and comparisons with PET and MRI scans to relate regional seizure propensity with AD pathologies.

People with Alzheimer’s disease (AD) are more likely to develop seizures than cognitively healthy people. In this study, we aimed to understand whether whole-brain network structure is related to this increased seizure likelihood. We used electroencephalography (EEG) to estimate brain networks from people with AD and healthy controls. We subsequently inserted these networks into a model brain and simulated disease progression by increasing the excitability of brain tissue. We found the simulated AD brains were more likely to develop seizures than the simulated control brains. No participants had seizures when we collected data, so our results suggest an increased probability of developing seizures at a future time for AD participants. Therefore functional brain network structure may play a role in increased seizure likelihood in AD. We also used the model to examine which brain regions were most important for generating seizures, and found that the seizure-generating regions corresponded to those typically affected in early AD. Our results also provide a potential explanation for why people with AD are more likely to have generalized seizures (i.e. seizures involving the whole brain, as opposed to ‘focal’ seizures which only involve certain areas) than the general population with epilepsy.

Ictal gamma-band interactions localize ictogenic nodes of the epileptic network in focal cortical dysplasia

OBJECTIVE

Epilepsy surgery fails in > 30% of patients with focal cortical dysplasia (FCD). The seizure persistence after surgery can be attributed to the inability to precisely localize the tissue with an endogenous potential to generate seizures. In this study, we aimed to identify the critical components of the epileptic network that were actively involved in seizure genesis.

METHODS

The directed transfer function was applied to intracranial EEG recordings and the effective connectivity was determined with a high temporal and frequency resolution. Pre-ictal network properties were compared with ictal epochs to identify regions actively generating ictal activity and discriminate them from the areas of propagation.

RESULTS

Analysis of 276 seizures from 30 patients revealed the existence of a seizure-related network reconfiguration in the gamma-band (25-170 Hz; p < 0.005) - ictogenic nodes. Unlike seizure onset zone, resecting the majority of ictogenic nodes correlated with favorable outcomes (p < 0.012).

CONCLUSION

The prerequisite to successful epilepsy surgery is the accurate identification of brain areas from which seizures arise. We show that in FCD-related epilepsy, gamma-band network markers can reliably identify and distinguish ictogenic areas in macroelectrode recordings, improve intracranial EEG interpretation and better delineate the epileptogenic zone.

SIGNIFICANCE

Ictogenic nodes localize the critical parts of the epileptogenic tissue and increase the diagnostic yield of intracranial evaluation.

Combined external and endoscopic transnasal approach with use of a diamond burr in the removal of a fishing harpoon hook traversing bilateral sphenoid sinuses in a 22-year-old man

A 22-year old construction worker was shot with a fishing harpoon gun on the left side of his face. He consulted at the emergency room 12 days postinjury, stable but with blurring of vision on the right. The shaft of the harpoon was protruding at the left preauricular area; the tip was neither visible nor palpable. Craniofacial CT scan and skull anteroposterolateral radiographs revealed the tip of the harpoon to be at the right orbital apex. A hook attached 1 cm from the tip was lodged in the sphenoid sinus. The hook was dismantled from the shaft via a combined external and endoscopic transnasal approach, enabling the shaft to be gently pulled. The hook, together with the tip, were removed endoscopically. The patient's visual acuity improved. He was discharged after 2 days on oral antibiotics with no deficits on follow-up.

Remote and Long-Term Self-Monitoring of Electroencephalographic and Noninvasive Measurable Variables at Home in Patients With Epilepsy (EEG@HOME): Protocol for an Observational Study

Background

Epileptic seizures are spontaneous events that severely affect the lives of patients due to their recurrence and unpredictability. The integration of new wearable and mobile technologies to collect electroencephalographic (EEG) and extracerebral signals in a portable system might be the solution to prospectively identify times of seizure occurrence or propensity. The performances of several seizure detection devices have been assessed by validated studies, and patient perspectives on wearables have been explored to better match their needs. Despite this, there is a major gap in the literature on long-term, real-life acceptability and performance of mobile technology essential to managing chronic disorders such as epilepsy.

Objective

EEG@HOME is an observational, nonrandomized, noninterventional study that aims to develop a new feasible procedure that allows people with epilepsy to independently, continuously, and safely acquire noninvasive variables at home. The data collected will be analyzed to develop a general model to predict periods of increased seizure risk.

Methods

A total of 12 adults with a diagnosis of pharmaco-resistant epilepsy and at least 20 seizures per year will be recruited at King’s College Hospital, London. Participants will be asked to self-apply an easy and portable EEG recording system (ANT Neuro) to record scalp EEG at home twice daily. From each serial EEG recording, brain network ictogenicity (BNI), a new biomarker of the propensity of the brain to develop seizures, will be extracted. A noninvasive wrist-worn device (Fitbit Charge 3; Fitbit Inc) will be used to collect non-EEG biosignals (heart rate, sleep quality index, and steps), and a smartphone app (Seer app; Seer Medical) will be used to collect data related to seizure occurrence, medication taken, sleep quality, stress, and mood. All data will be collected continuously for 6 months. Standardized questionnaires (the Post-Study System Usability Questionnaire and System Usability Scale) will be completed to assess the acceptability and feasibility of the procedure. BNI, continuous wrist-worn sensor biosignals, and electronic survey data will be correlated with seizure occurrence as reported in the diary to investigate their potential values as biomarkers of seizure risk.

Results

The EEG@HOME project received funding from Epilepsy Research UK in 2018 and was approved by the Bromley Research Ethics Committee in March 2020. The first participants were enrolled in October 2020, and we expect to publish the first results by the end of 2022.

Conclusions

With the EEG@HOME study, we aim to take advantage of new advances in remote monitoring technology, including self-applied EEG, to investigate the feasibility of long-term disease self-monitoring. Further, we hope our study will bring new insights into noninvasively collected personalized risk factors of seizure occurrence and seizure propensity that may help to mitigate one of the most difficult aspects of refractory epilepsy: the unpredictability of seizure occurrence.

International Registered Report Identifier (IRRID)

PRR1-10.2196/25309

In-Motion-App for remote General Movement Assessment: a multi-site observational study

OBJECTIVES

To determine whether videos taken by parents of their infants' spontaneous movements were in accordance with required standards in the In-Motion-App, and whether the videos could be remotely scored by a trained General Movement Assessment (GMA) observer. Additionally, to assess the feasibility of using home-based video recordings for automated tracking of spontaneous movements, and to examine parents' perceptions and experiences of taking videos in their homes.

DESIGN

The study was a multi-centre prospective observational study.

SETTING

Parents/families of high-risk infants in tertiary care follow-up programmes in Norway, Denmark and Belgium.

METHODS

Parents/families were asked to video record their baby in accordance with the In-Motion standards which were based on published GMA criteria and criteria covering lighting and stability of smartphone. Videos were evaluated as GMA 'scorable' or 'non-scorable' based on predefined criteria. The accuracy of a 7-point body tracker software was compared with manually annotated body key points. Parents were surveyed about the In-Motion-App information and clarity.

PARTICIPANTS

The sample comprised 86 parents/families of high-risk infants.

RESULTS

The 86 parent/families returned 130 videos, and 121 (96%) of them were in accordance with the requirements for GMA assessment. The 7-point body tracker software detected more than 80% of body key point positions correctly. Most families found the instructions for filming their baby easy to follow, and more than 90% reported that they did not become more worried about their child's development through using the instructions.

CONCLUSIONS

This study reveals that a short instructional video enabled parents to video record their infant's spontaneous movements in compliance with the standards required for remote GMA. Further, an accurate automated body point software detecting infant body landmarks in smartphone videos will facilitate clinical and research use soon. Home-based video recordings could be performed without worrying parents about their child's development.

TRIALS REGISTRATION NUMBER

NCT03409978.

A computational biomarker of juvenile myoclonic epilepsy from resting-state MEG

•

Computational modelling is combined with MEG to differentiate people with juvenile myoclonic epilepsy from healthy controls.

•

Brain network ictogenicity (BNI) was found higher in people with juvenile myoclonic epilepsy relative to healthy controls.

•

BNI’s classification accuracy in our cohort was 73%.

Objective

For people with idiopathic generalized epilepsy, functional networks derived from their resting-state scalp electrophysiological recordings have shown an inherent higher propensity to generate seizures than those from healthy controls when assessed using the concept of brain network ictogenicity (BNI). Herein we tested whether the BNI framework is applicable to resting-state magnetoencephalography (MEG) from people with juvenile myoclonic epilepsy (JME).

Methods

The BNI framework consists in deriving a functional network from apparently normal brain activity, placing a mathematical model of ictogenicity into the network and then computing how often such network generates seizures in silico. We considered data from 26 people with JME and 26 healthy controls.

Results

We found that resting-state MEG functional networks from people with JME are characterized by a higher propensity to generate seizures (i.e., higher BNI) than those from healthy controls. We found a classification accuracy of 73%.

Conclusions

The BNI framework is applicable to MEG and was capable of differentiating people with epilepsy from healthy controls.

Significance

The BNI framework may be applied to resting-state MEG to aid in epilepsy diagnosis.

FitzHugh-Nagumo oscillators on complex networks mimic epileptic-seizure-related synchronization phenomena

We study patterns of partial synchronization in a network of FitzHugh-Nagumo oscillators with empirical structural connectivity measured in human subjects. We report the spontaneous occurrence of synchronization phenomena that closely resemble the ones seen during epileptic seizures in humans. In order to obtain deeper insights into the interplay between dynamics and network topology, we perform long-term simulations of oscillatory dynamics on different paradigmatic network structures: random networks, regular nonlocally coupled ring networks, ring networks with fractal connectivities, and small-world networks with various rewiring probability. Among these networks, a small-world network with intermediate rewiring probability best mimics the findings achieved with the simulations using the empirical structural connectivity. For the other network topologies, either no spontaneously occurring epileptic-seizure-related synchronization phenomena can be observed in the simulated dynamics, or the overall degree of synchronization remains high throughout the simulation. This indicates that a topology with some balance between regularity and randomness favors the self-initiation and self-termination of episodes of seizure-like strong synchronization.

Epilepsy surgery: Evaluating robustness using dynamic network models

Epilepsy is one of the most common neurological conditions affecting over 65 million people worldwide. Over one third of people with epilepsy are considered refractory: they do not respond to drug treatments. For this significant cohort of people, surgery is a potentially transformative treatment. However, only a small minority of people with refractory epilepsy are considered suitable for surgery, and long-term seizure freedom is only achieved in half the cases. Recently, several computational approaches have been proposed to support presurgical planning. Typically, these approaches use a dynamic network model to explore the potential impact of surgical resection in silico. The network component of the model is informed by clinical imaging data and is considered static thereafter. This assumption critically overlooks the plasticity of the brain and, therefore, how continued evolution of the brain network post-surgery may impact upon the success of a resection in the longer term. In this work, we use a simplified dynamic network model, which describes transitions to seizures, to systematically explore how the network structure influences seizure propensity, both before and after virtual resections. We illustrate key results in small networks, before extending our findings to larger networks. We demonstrate how the evolution of brain networks post resection can result in a return to increased seizure propensity. Our results effectively determine the robustness of a given resection to possible network reconfigurations and so provide a potential strategy for optimizing long-term seizure freedom.

The value of intra-operative electrographic biomarkers for tailoring during epilepsy surgery: from group-level to patient-level analysis

Signal analysis biomarkers, in an intra-operative setting, may be complementary tools to guide and tailor the resection in drug-resistant focal epilepsy patients. Effective assessment of biomarker performances are needed to evaluate their clinical usefulness and translation. We defined a realistic ground-truth scenario and compared the effectiveness of different biomarkers alone and combined to localize epileptogenic tissue during surgery. We investigated the performances of univariate, bivariate and multivariate signal biomarkers applied to 1 min inter-ictal intra-operative electrocorticography to discriminate between epileptogenic and non-epileptogenic locations in 47 drug-resistant people with epilepsy (temporal and extra-temporal) who had been seizure-free one year after the operation. The best result using a single biomarker was obtained using the phase-amplitude coupling measure for which the epileptogenic tissue was localized in 17 out of 47 patients. Combining the whole set of biomarkers provided an improvement of the performances: 27 out of 47 patients. Repeating the analysis only on the temporal-lobe resections we detected the epileptogenic tissue in 29 out of 30 combining all the biomarkers. We suggest that the assessment of biomarker performances on a ground-truth scenario is required to have a proper estimate on how biomarkers translate into clinical use. Phase-amplitude coupling seems the best performing single biomarker and combining biomarkers improves localization of epileptogenic tissue. Performance achieved is not adequate as a tool in the operation theater yet, but it can improve the understanding of pathophysiological process.

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.

The Bayesian Virtual Epileptic Patient: A probabilistic framework designed to infer the spatial map of epileptogenicity in a personalized large-scale brain model of epilepsy spread

Despite the importance and frequent use of Bayesian frameworks in brain network modeling for parameter inference and model prediction, the advanced sampling algorithms implemented in probabilistic programming languages to overcome the inference difficulties have received relatively little attention in this context. In this technical note, we propose a probabilistic framework, namely the Bayesian Virtual Epileptic Patient (BVEP), which relies on the fusion of structural data of individuals to infer the spatial map of epileptogenicity in a personalized large-scale brain model of epilepsy spread. To invert the individualized whole-brain model employed in this study, we use the recently developed algorithms known as No-U-Turn Sampler (NUTS) as well as Automatic Differentiation Variational Inference (ADVI). Our results indicate that NUTS and ADVI accurately estimate the degree of epileptogenicity of brain regions, therefore, the hypothetical brain areas responsible for the seizure initiation and propagation, while the convergence diagnostics and posterior behavior analysis validate the reliability of the estimations. Moreover, we illustrate the efficiency of the transformed non-centered parameters in comparison to centered form of parameterization. The Bayesian framework used in this work proposes an appropriate patient-specific strategy for estimating the epileptogenicity of the brain regions to improve outcome after epilepsy surgery.

The Role of Excitability and Network Structure in the Emergence of Focal and Generalized Seizures

Epileptic seizures are generally classified as either focal or generalized. It had been traditionally assumed that focal seizures imply localized brain abnormalities, whereas generalized seizures involve widespread brain pathologies. However, recent evidence suggests that large-scale brain networks are involved in the generation of focal seizures, and generalized seizures can originate in localized brain regions. Herein we study how network structure and tissue heterogeneities underpin the emergence of focal and widespread seizure dynamics. Mathematical modeling of seizure emergence in brain networks enables the clarification of the characteristics responsible for focal and generalized seizures. We consider neural mass network dynamics of seizure generation in exemplar synthetic networks and we measure the variance in ictogenicity across the network. Ictogenicity is defined as the involvement of network nodes in seizure activity, and its variance is used to quantify whether seizure patterns are focal or widespread across the network. We address both the influence of network structure and different excitability distributions across the network on the ictogenic variance. We find that this variance depends on both network structure and excitability distribution. High variance, i.e., localized seizure activity, is observed in networks highly heterogeneous with regard to the distribution of connections or excitabilities. However, networks that are both heterogeneous in their structure and excitability can underlie the emergence of generalized seizures, depending on the interplay between structure and excitability. Thus, our results imply that the emergence of focal and generalized seizures is underpinned by an interplay between network structure and excitability distribution.

Computational modelling in source space from scalp EEG to inform presurgical evaluation of epilepsy

OBJECTIVE

The effectiveness of intracranial electroencephalography (iEEG) to inform epilepsy surgery depends on where iEEG electrodes are implanted. This decision is informed by noninvasive recording modalities such as scalp EEG. Herein we propose a framework to interrogate scalp EEG and determine epilepsy lateralization to aid in electrode implantation.

METHODS

We use eLORETA to map source activities from seizure epochs recorded from scalp EEG and consider 15 regions of interest (ROIs). Functional networks are then constructed using the phase-locking value and studied using a mathematical model. By removing different ROIs from the network and simulating their impact on the network's ability to generate seizures in silico, the framework provides predictions of epilepsy lateralization. We consider 15 individuals from the EPILEPSIAE database and study a total of 62 seizures. Results were assessed by taking into account actual intracranial implantations and surgical outcome.

RESULTS

The framework provided potentially useful information regarding epilepsy lateralization in 12 out of the 15 individuals (p=0.02, binomial test).

CONCLUSIONS

Our results show promise for the use of this framework to better interrogate scalp EEG to determine epilepsy lateralization.

SIGNIFICANCE

The framework may aid clinicians in the decision process to define where to implant electrodes for intracranial monitoring.

Quantification and Selection of Ictogenic Zones in Epilepsy Surgery

Network models of brain dynamics provide valuable insight into the healthy functioning of the brain and how this breaks down in disease. A pertinent example is the use of network models to understand seizure generation (ictogenesis) in epilepsy. Recently, computational models have emerged to aid our understanding of seizures and to predict the outcome of surgical perturbations to brain networks. Such approaches provide the opportunity to quantify the effect of removing regions of tissue from brain networks and thereby search for the optimal resection strategy. Here, we use computational models to elucidate how sets of nodes contribute to the ictogenicity of networks. In small networks we fully elucidate the ictogenicity of all possible sets of nodes and demonstrate that the distribution of ictogenicity across sets depends on network topology. However, the full elucidation is a combinatorial problem that becomes intractable for large networks. Therefore, we combine computational models with a genetic algorithm to search for minimal sets of nodes that contribute significantly to ictogenesis. We demonstrate the potential applicability of these methods in practice by identifying optimal sets of nodes to resect in networks derived from 20 individuals who underwent resective surgery for epilepsy. We show that they have the potential to aid epilepsy surgery by suggesting alternative resection sites as well as facilitating the avoidance of brain regions that should not be resected.

Seizure prediction - ready for a new era

Epilepsy is a common disorder characterized by recurrent seizures. An overwhelming majority of people with epilepsy regard the unpredictability of seizures as a major issue. More than 30 years of international effort have been devoted to the prediction of seizures, aiming to remove the burden of unpredictability and to couple novel, time-specific treatment to seizure prediction technology. A highly influential review published in 2007 concluded that insufficient evidence indicated that seizures could be predicted. Since then, several advances have been made, including successful prospective seizure prediction using intracranial EEG in a small number of people in a trial of a real-time seizure prediction device. In this Review, we examine advances in the field, including EEG databases, seizure prediction competitions, the prospective trial mentioned and advances in our understanding of the mechanisms of seizures. We argue that these advances, together with statistical evaluations, set the stage for a resurgence in efforts towards the development of seizure prediction methodologies. We propose new avenues of investigation involving a synergy between mechanisms, models, data, devices and algorithms and refine the existing guidelines for the development of seizure prediction technology to instigate development of a solution that removes the burden of the unpredictability of seizures.

Intracortical neural activity distal to seizure-onset-areas predicts human focal seizures

The apparent unpredictability of epileptic seizures has a major impact in the quality of life of people with pharmacologically resistant seizures. Here, we present initial results and a proof-of-concept of how focal seizures can be predicted early in advance based on intracortical signals recorded from small neocortical patches away from identified seizure onset areas. We show that machine learning algorithms can discriminate between interictal and preictal periods based on multiunit activity (i.e. thresholded action potential counts) and multi-frequency band local field potentials recorded via 4 X 4 mm² microelectrode arrays. Microelectrode arrays were implanted in 5 patients undergoing neuromonitoring for resective surgery. Post-implant analysis revealed arrays were outside the seizure onset areas. Preictal periods were defined as the 1-hour period leading to a seizure. A 5-minute gap between the preictal period and the putative seizure onset was enforced to account for potential errors in the determination of actual seizure onset times. We used extreme gradient boosting and long short-term memory networks for prediction. Prediction accuracy based on the area under the receiver operating characteristic curves reached 90% for at least one feature type in each patient. Importantly, successful prediction could be achieved based exclusively on multiunit activity. This result indicates that preictal activity in the recorded neocortical patches involved not only subthreshold postsynaptic potentials, perhaps driven by the distal seizure onset areas, but also neuronal spiking in distal recurrent neocortical networks. Beyond the commonly identified seizure onset areas, our findings point to the engagement of large-scale neuronal networks in the neural dynamics building up toward a seizure. Our initial results obtained on currently available human intracortical microelectrode array recordings warrant new studies on larger datasets, and open new perspectives for seizure prediction and control by emphasizing the contribution of multiscale neural signals in large-scale neuronal networks.

Revealing epilepsy type using a computational analysis of interictal EEG

Seizure onset in epilepsy can usually be classified as focal or generalized, based on a combination of clinical phenomenology of the seizures, EEG recordings and MRI. This classification may be challenging when seizures and interictal epileptiform discharges are infrequent or discordant, and MRI does not reveal any apparent abnormalities. To address this challenge, we introduce the concept of Ictogenic Spread (IS) as a prediction of how pathological electrical activity associated with seizures will propagate throughout a brain network. This measure is defined using a person-specific computer representation of the functional network of the brain, constructed from interictal EEG, combined with a computer model of the transition from background to seizure-like activity within nodes of a distributed network. Applying this method to a dataset comprising scalp EEG from 38 people with epilepsy (17 with genetic generalized epilepsy (GGE), 21 with mesial temporal lobe epilepsy (mTLE)), we find that people with GGE display a higher IS in comparison to those with mTLE. We propose IS as a candidate computational biomarker to classify focal and generalized epilepsy using interictal EEG.

Dynamic brain network states in human generalized spike-wave discharges

Generalized spike-wave discharges in idiopathic generalized epilepsy are conventionally assumed to have abrupt onset and offset. However, in rodent models, discharges emerge during a dynamic evolution of brain network states, extending several seconds before and after the discharge. In human idiopathic generalized epilepsy, simultaneous EEG and functional MRI shows cortical regions may be active before discharges, and network connectivity around discharges may not be normal. Here, in human idiopathic generalized epilepsy, we investigated whether generalized spike-wave discharges emerge during a dynamic evolution of brain network states. Using EEG-functional MRI, we studied 43 patients and 34 healthy control subjects. We obtained 95 discharges from 20 patients. We compared data from patients with discharges with data from patients without discharges and healthy controls. Changes in MRI (blood oxygenation level-dependent) signal amplitude in discharge epochs were observed only at and after EEG onset, involving a sequence of parietal and frontal cortical regions then thalamus (P < 0.01, across all regions and measurement time points). Examining MRI signal phase synchrony as a measure of functional connectivity between each pair of 90 brain regions, we found significant connections (P < 0.01, across all connections and measurement time points) involving frontal, parietal and occipital cortex during discharges, and for 20 s after EEG offset. This network prominent during discharges showed significantly low synchrony (below 99% confidence interval for synchrony in this network in non-discharge epochs in patients) from 16 s to 10 s before discharges, then ramped up steeply to a significantly high level of synchrony 2 s before discharge onset. Significant connections were seen in a sensorimotor network in the minute before discharge onset. This network also showed elevated synchrony in patients without discharges compared to healthy controls (P = 0.004). During 6 s prior to discharges, additional significant connections to this sensorimotor network were observed, involving prefrontal and precuneus regions. In healthy subjects, significant connections involved a posterior cortical network. In patients with discharges, this posterior network showed significantly low synchrony during the minute prior to discharge onset. In patients without discharges, this network showed the same level of synchrony as in healthy controls. Our findings suggest persistently high sensorimotor network synchrony, coupled with transiently (at least 1 min) low posterior network synchrony, may be a state predisposing to generalized spike-wave discharge onset. Our findings also show that EEG onset and associated MRI signal amplitude change is embedded in a considerably longer period of evolving brain network states before and after discharge events.

The role that choice of model plays in predictions for epilepsy surgery

Mathematical modelling has been widely used to predict the effects of perturbations to brain networks. An important example is epilepsy surgery, where the perturbation in question is the removal of brain tissue in order to render the patient free of seizures. Different dynamical models have been proposed to represent transitions to ictal states in this context. However, our choice of which mathematical model to use to address this question relies on making assumptions regarding the mechanism that defines the transition from background to the seizure state. Since these mechanisms are unknown, it is important to understand how predictions from alternative dynamical descriptions compare. Herein we evaluate to what extent three different dynamical models provide consistent predictions for the effect of removing nodes from networks. We show that for small, directed, connected networks the three considered models provide consistent predictions. For larger networks, predictions are shown to be less consistent. However consistency is higher in networks that have sufficiently large differences in ictogenicity between nodes. We further demonstrate that heterogeneity in ictogenicity across nodes correlates with variability in the number of connections for each node.

A Model-Based Assessment of the Seizure Onset Zone Predictive Power to Inform the Epileptogenic Zone

Epilepsy surgery is a clinical procedure that aims to remove the brain tissue responsible for the emergence of seizures, the epileptogenic zone (EZ). It is preceded by an evaluation to determine the brain tissue that must be resected. The identification of the seizure onset zone (SOZ) from intracranial EEG recordings stands as one of the key proxies for the EZ. In this study we used computational models of epilepsy to assess to what extent the SOZ may or may not represent the EZ. We considered a set of different synthetic networks (e.g., regular, small-world, random, and scale-free networks) to represent large-scale brain networks and a phenomenological network model of seizure generation. In the model, the SOZ was inferred from the seizure likelihood (SL), a measure of the propensity of single nodes to produce epileptiform dynamics, whilst a surgery corresponded to the removal of nodes and connections from the network. We used the concept of node ictogenicity (NI) to quantify the effectiveness of each node removal on reducing the network's propensity to generate seizures. This framework enabled us to systematically compare the SOZ and the seizure control achieved by each considered surgery. Specifically, we compared the distributions of SL and NI across different networks. We found that SL and NI were concordant when all nodes were similarly ictogenic, whereas when there was a small fraction of nodes with high NI, the SL was not specific at identifying these nodes. We further considered networks with heterogeneous node excitabilities, i.e., nodes with different susceptibilities of being engaged in seizure activity, to understand how such heterogeneity may affect the relationship between SL and NI. We found that while SL and NI are concordant when there is a small fraction of hyper-excitable nodes in a network that is otherwise homogeneous, they do diverge if the network is heterogeneous, such as in scale-free networks. We observe that SL is highly dependent on node excitabilities, whilst the effect of surgical resections as revealed by NI is mostly determined by network structure. Together our results suggest that the SOZ is not always a good marker of the EZ.

Sensorimotor network hypersynchrony as an endophenotype in families with genetic generalized epilepsy: A resting-state functional magnetic resonance imaging study

Recent evidence suggests that three specific brain networks show state-dependent levels of synchronization before, during, and after episodes of generalized spike-wave discharges (GSW) in patients with genetic generalized epilepsy (GGE). Here, we investigate whether synchronization in these networks differs between patients with GGE (n = 13), their unaffected first-degree relatives (n = 17), and healthy controls (n = 18). All subjects underwent two 10-minute simultaneous electroencephalographic-functional magnetic resonance imaging (fMRI) recordings without GSW. Whole-brain data were divided into 90 regions, and blood oxygen level-dependent (BOLD) phase synchrony in a 0.04-0.07-Hz band was estimated between all pairs of regions. Three networks were defined: (1) the network with highest synchrony during GSW events, (2) a sensorimotor network, and (3) an occipital network. Average synchrony (mean node degree) was inferred across each network over time. Notably, synchrony was significantly higher in the sensorimotor network in patients and in unaffected relatives, compared to controls. There was a trend toward higher synchrony in the GSW network in patients and in unaffected relatives. There was no difference between groups for the occipital network. Our findings provide evidence that elevated fMRI BOLD synchrony in a sensorimotor network is a state-independent endophenotype of GGE, present in patients in the absence of GSW, and present in unaffected relatives.

Slower alpha rhythm associates with poorer seizure control in epilepsy

Objective

Slowing and frontal spread of the alpha rhythm have been reported in multiple epilepsy syndromes. We investigated whether these phenomena are associated with seizure control.

Methods

We prospectively acquired resting-state electroencephalogram (EEG) in 63 patients with focal and idiopathic generalized epilepsy (FE and IGE) and 39 age- and gender-matched healthy subjects (HS). Patients were divided into good and poor (≥4 seizures/12 months) seizure control groups based on self-reports and clinical records. We computed spectral power from 20-sec EEG segments during eyes-closed wakefulness, free of interictal abnormalities, and quantified power in high- and low-alpha bands. Analysis of covariance and post hoc t-tests were used to assess group differences in alpha-power shift across all EEG channels. Permutation-based statistics were used to assess the topography of this shift across the whole scalp.

Results

Compared to HS, patients showed a statistically significant shift of spectral power from high- to low-alpha frequencies (effect size g = 0.78 [95% confidence interval 0.43, 1.20]). This alpha-power shift was driven by patients with poor seizure control in both FE and IGE (g = 1.14, [0.65, 1.74]), and occurred over midline frontal and bilateral occipital regions. IGE exhibited less alpha power shift compared to FE over bilateral frontal regions (g = -1.16 [-0.68, -1.74]). There was no interaction between syndrome and seizure control. Effects were independent of antiepileptic drug load, time of day, or subgroup definitions.

Interpretation

Alpha slowing and anteriorization are a robust finding in patients with epilepsy and might represent a generic indicator of seizure liability.

Rationale and design of a single-blind, randomised controlled trial of exercise training for managing learning and memory impairment in persons with multiple sclerosis

INTRODUCTION

This randomised controlled trial (RCT) examines treadmill walking exercise training effects on learning and memory performance, hippocampal volume, and hippocampal resting-state functional connectivity in persons with multiple sclerosis (MS) who have objective impairments in learning new information.

METHODS AND ANALYSIS

Forty fully ambulatory persons with MS who demonstrate objective learning and memory impairments will be randomly assigned into either the intervention or active control study conditions. The intervention condition involves supervised, progressive treadmill walking exercise training three times per week for a 3-month period. The active control condition involves supervised, progressive low-intensity resistive exercise that will be delivered at the same frequency as the intervention condition. The primary outcome will involve composite performance on neuropsychological learning and memory tests, and the secondary outcomes involve MRI measures of hippocampal volume and resting-state functional connectivity administered before and after the 3-month study period. Outcomes will be administered by treatment-blinded assessors using alternate test forms to minimise practice effects, and MRI data processing will be performed by blinded data analysts.

ETHICS AND DISSEMINATION

This study has been approved by a university institutional review board. The primary results will be disseminated via peer-reviewed publications and the final data will be made available to third parties in applicable data repositories. If successful, the results from this study will eventually inform subsequent RCTs for developing physical rehabilitation interventions (ie, treadmill walking exercise training) for improving learning and memory and its relationship with hippocampal outcomes in larger samples of cognitively impaired persons with MS. The results from this early-phase RCT will further lay preliminary groundwork for ultimately providing clinicians and patients with guidelines for better using chronic treadmill walking exercise for improving cognition and brain health. This approach is paramount as learning and memory impairment is common, burdensome and poorly managed in MS.

TRIAL REGISTRATION NUMBER

NCT03319771; Pre-results.

Characterizing EEG Cortical Dynamics and Connectivity with Responses to Single Pulse Electrical Stimulation (SPES)

Objectives: To model cortical connections in order to characterize their oscillatory behavior and role in the generation of spontaneous electroencephalogram (EEG). Methods: We studied averaged responses to single pulse electrical stimulation (SPES) from the non-epileptogenic hemisphere of five patients assessed with intracranial EEG who became seizure free after contralateral temporal lobectomy. Second-order control system equations were modified to characterize the systems generating a given response. SPES responses were modeled as responses to a unit step input. EEG power spectrum was calculated on the 20[Formula: see text]s preceding SPES. Results: 121 channels showed responses to 32 stimulation sites. A single system could model the response in 41.3% and two systems were required in 58.7%. Peaks in the frequency response of the models tended to occur within the frequency range of most activity on the spontaneous EEG. Discrepancies were noted between activity predicted by models and activity recorded in the spontaneous EEG. These discrepancies could be explained by the existence of alpha rhythm or interictal epileptiform discharges. Conclusions: Cortical interactions shown by SPES can be described as control systems which can predict cortical oscillatory behavior. The method is unique as it describes connectivity as well as dynamic interactions.

Progressive topological disorganization of brain network in focal epilepsy

OBJECTIVE

Increasing evidence has suggested that epilepsy is a network disease. Graph theory is a mathematical tool that allows for the analysis and quantification of the brain network. We aimed to evaluate the influences of duration of epilepsy on the topological organization of brain network in focal epilepsy patients with normal MRI using the graph theoretical analysis based on diffusion tenor imaging.

METHODS

We prospectively enrolled 66 patients with focal epilepsy (18/66 patients were newly diagnosed) and 84 healthy subjects. All of the patients with epilepsy had normal MRI on visual inspection. All of the subjects underwent diffusion tensor imaging that was analyzed using graph theory to obtain network measures.

RESULTS

The measures of characteristic path length and small-worldness in the patients with focal epilepsy were significantly decreased, even after multiple corrections (P < .01). Moreover, the measures including mean clustering coefficient and global efficiency in the patients with epilepsy had strong tendency to decrease compared to those in healthy subjects (P = .0153 and P = .0138, respectively). When comparing the measures among the patients with newly diagnosed/chronic epilepsy and healthy subjects using ANOVA, the characteristic path length (P = .006), small-worldness (P = .032), and global efficiency (P = .004) were significantly different. In addition, the duration of epilepsy was negatively correlated with global efficiency (r = -.249, P = .0454).

CONCLUSIONS

We newly found a progressive topological disorganization of the brain network in focal epilepsy. In addition, we demonstrated disrupted topological organization in focal epilepsy, shifting toward a more random state.

Chimera states in brain networks: Empirical neural vs. modular fractal connectivity

Complex spatiotemporal patterns, called chimera states, consist of coexisting coherent and incoherent domains and can be observed in networks of coupled oscillators. The interplay of synchrony and asynchrony in complex brain networks is an important aspect in studies of both the brain function and disease. We analyse the collective dynamics of FitzHugh-Nagumo neurons in complex networks motivated by its potential application to epileptology and epilepsy surgery. We compare two topologies: an empirical structural neural connectivity derived from diffusion-weighted magnetic resonance imaging and a mathematically constructed network with modular fractal connectivity. We analyse the properties of chimeras and partially synchronized states and obtain regions of their stability in the parameter planes. Furthermore, we qualitatively simulate the dynamics of epileptic seizures and study the influence of the removal of nodes on the network synchronizability, which can be useful for applications to epileptic surgery.

Elevated Ictal Brain Network Ictogenicity Enables Prediction of Optimal Seizure Control

Recent studies have shown that mathematical models can be used to analyze brain networks by quantifying how likely they are to generate seizures. In particular, we have introduced the quantity termed brain network ictogenicity (BNI), which was demonstrated to have the capability of differentiating between functional connectivity (FC) of healthy individuals and those with epilepsy. Furthermore, BNI has also been used to quantify and predict the outcome of epilepsy surgery based on FC extracted from pre-operative ictal intracranial electroencephalography (iEEG). This modeling framework is based on the assumption that the inferred FC provides an appropriate representation of an ictogenic network, i.e., a brain network responsible for the generation of seizures. However, FC networks have been shown to change their topology depending on the state of the brain. For example, topologies during seizure are different to those pre- and post-seizure. We therefore sought to understand how these changes affect BNI. We studied peri-ictal iEEG recordings from a cohort of 16 epilepsy patients who underwent surgery and found that, on average, ictal FC yield higher BNI relative to pre- and post-ictal FC. However, elevated ictal BNI was not observed in every individual, rather it was typically observed in those who had good post-operative seizure control. We therefore hypothesize that elevated ictal BNI is indicative of an ictogenic network being appropriately represented in the FC. We evidence this by demonstrating superior model predictions for post-operative seizure control in patients with elevated ictal BNI.

From Structure to Activity: Using Centrality Measures to Predict Neuronal Activity

It is clear that the topological structure of a neural network somehow determines the activity of the neurons within it. In the present work, we ask to what extent it is possible to examine the structural features of a network and learn something about its activity? Specifically, we consider how the centrality (the importance of a node in a network) of a neuron correlates with its firing rate. To investigate, we apply an array of centrality measures, including In-Degree, Closeness, Betweenness, Eigenvector, Katz, PageRank, Hyperlink-Induced Topic Search (HITS) and NeuronRank to Leaky-Integrate and Fire neural networks with different connectivity schemes. We find that Katz centrality is the best predictor of firing rate given the network structure, with almost perfect correlation in all cases studied, which include purely excitatory and excitatory–inhibitory networks, with either homogeneous connections or a small-world structure. We identify the properties of a network which will cause this correlation to hold. We argue that the reason Katz centrality correlates so highly with neuronal activity compared to other centrality measures is because it nicely captures disinhibition in neural networks. In addition, we argue that these theoretical findings are applicable to neuroscientists who apply centrality measures to functional brain networks, as well as offer a neurophysiological justification to high level cognitive models which use certain centrality measures.

An optimal strategy for epilepsy surgery: Disruption of the rich-club?

Surgery is a therapeutic option for people with epilepsy whose seizures are not controlled by anti-epilepsy drugs. In pre-surgical planning, an array of data modalities, often including intra-cranial EEG, is used in an attempt to map regions of the brain thought to be crucial for the generation of seizures. These regions are then resected with the hope that the individual is rendered seizure free as a consequence. However, post-operative seizure freedom is currently sub-optimal, suggesting that the pre-surgical assessment may be improved by taking advantage of a mechanistic understanding of seizure generation in large brain networks. Herein we use mathematical models to uncover the relative contribution of regions of the brain to seizure generation and consequently which brain regions should be considered for resection. A critical advantage of this modeling approach is that the effect of different surgical strategies can be predicted and quantitatively compared in advance of surgery. Herein we seek to understand seizure generation in networks with different topologies and study how the removal of different nodes in these networks reduces the occurrence of seizures. Since this a computationally demanding problem, a first step for this aim is to facilitate tractability of this approach for large networks. To do this, we demonstrate that predictions arising from a neural mass model are preserved in a lower dimensional, canonical model that is quicker to simulate. We then use this simpler model to study the emergence of seizures in artificial networks with different topologies, and calculate which nodes should be removed to render the network seizure free. We find that for scale-free and rich-club networks there exist specific nodes that are critical for seizure generation and should therefore be removed, whereas for small-world networks the strategy should instead focus on removing sufficient brain tissue. We demonstrate the validity of our approach by analysing intra-cranial EEG recordings from a database comprising 16 patients who have undergone epilepsy surgery, revealing rich-club structures within the obtained functional networks. We show that the postsurgical outcome for these patients was better when a greater proportion of the rich club was removed, in agreement with our theoretical predictions.

Temporally correlated fluctuations drive epileptiform dynamics

Macroscopic models of brain networks typically incorporate assumptions regarding the characteristics of afferent noise, which is used to represent input from distal brain regions or ongoing fluctuations in non-modelled parts of the brain. Such inputs are often modelled by Gaussian white noise which has a flat power spectrum. In contrast, macroscopic fluctuations in the brain typically follow a 1/fb spectrum. It is therefore important to understand the effect on brain dynamics of deviations from the assumption of white noise. In particular, we wish to understand the role that noise might play in eliciting aberrant rhythms in the epileptic brain. To address this question we study the response of a neural mass model to driving by stochastic, temporally correlated input. We characterise the model in terms of whether it generates "healthy" or "epileptiform" dynamics and observe which of these dynamics predominate under different choices of temporal correlation and amplitude of an Ornstein-Uhlenbeck process. We find that certain temporal correlations are prone to eliciting epileptiform dynamics, and that these correlations produce noise with maximal power in the δ and θ bands. Crucially, these are rhythms that are found to be enhanced prior to seizures in humans and animal models of epilepsy. In order to understand why these rhythms can generate epileptiform dynamics, we analyse the response of the model to sinusoidal driving and explain how the bifurcation structure of the model gives rise to these findings. Our results provide insight into how ongoing fluctuations in brain dynamics can facilitate the onset and propagation of epileptiform rhythms in brain networks. Furthermore, we highlight the need to combine large-scale models with noise of a variety of different types in order to understand brain (dys-)function.

Bursts of seizures in long-term recordings of human focal epilepsy

OBJECTIVE

We report on temporally clustered seizures detected from continuous long-term ambulatory human electroencephalographic data. The objective was to investigate short-term seizure clustering, which we have termed bursting, and consider implications for patient care, seizure prediction, and evaluating therapies.

METHODS

Chronic ambulatory intracranial electroencephalography (EEG) data collected for the purpose of seizure prediction were annotated to identify seizure events. A detection algorithm was used to identify bursts of events. Burst events were compared to nonburst events to evaluate event dispersion, duration and dynamics.

RESULTS

Bursts of seizures were present in 6 of 15 subjects, and detections were consistent over long-term monitoring (>2 years). Subjects with bursts of seizures had highly overdispersed seizure rates, compared to other subjects. There was a complicated relationship between bursts and clinical seizures, although bursts were associated with multimodal distributions of seizure duration, and poorer predictive outcomes. For three subjects, bursts demonstrated distinctive preictal dynamics compared to clinical seizures.

SIGNIFICANCE

We have previously hypothesized that there are distinct physiologic pathways underlying short- and long-duration seizures. Herein we show that burst seizures fall almost exclusively within the short population of seizure durations; however, a short duration event was not sufficient to induce or imply bursting. We can therefore conclude that in addition to distinct mechanisms underlying seizure duration, there are separate factors regulating bursts of seizures. We show that bursts were a robust phenomenon in our patient cohort, which were consistent with overdispersed seizure rates, suggesting long-memory dynamics.

Recurring Functional Interactions Predict Network Architecture of Interictal and Ictal States in Neocortical Epilepsy

Human epilepsy patients suffer from spontaneous seizures, which originate in brain regions that also subserve normal function. Prior studies demonstrate focal, neocortical epilepsy is associated with dysfunction, several hours before seizures. How does the epileptic network perpetuate dysfunction during baseline periods? To address this question, we developed an unsupervised machine learning technique to disentangle patterns of functional interactions between brain regions, or subgraphs, from dynamic functional networks constructed from approximately 100 h of intracranial recordings in each of 22 neocortical epilepsy patients. Using this approach, we found: (1) subgraphs from ictal (seizure) and interictal (baseline) epochs are topologically similar, (2) interictal subgraph topology and dynamics can predict brain regions that generate seizures, and (3) subgraphs undergo slower and more coordinated fluctuations during ictal epochs compared to interictal epochs. Our observations suggest that seizures mark a critical shift away from interictal states that is driven by changes in the dynamical expression of strongly interacting components of the epileptic network.

Predicting neurosurgical outcomes in focal epilepsy patients using computational modelling

See Eissa and Schevon (doi:10.1093/aww332) for a scientific commentary on this article.

Surgery can be a last resort for patients with intractable, medically refractory epilepsy. For many of these patients, however, there is substantial risk that the surgery will be ineffective. The prediction of who is likely to benefit from a surgical approach is crucial for being able to inform patients better, conduct principled prospective clinical trials, and ultimately tailor therapeutic approaches to these patients more effectively. Dynamical computational models, informed with patient data, can be used to make predictions and give mechanistic insight. In this study, we develop patient-specific dynamical network models of epileptogenic cortex. We infer the network connectivity matrix from non-seizure electrographic recordings of patients and use these connectivity matrices as the network structure in our model. The model simulates the dynamics of a bi-stable switch at every node in this network, meaning that every node starts in a background state, but has the ability to transit to a co-existing seizure state. Whether a transition happens in a node is partly determined by the stochastic nature of the input to the node, but also by the input the node receives from other connected nodes in the network. By conducting simulations with such a model, we can detect the average transition time for nodes in a given network, and therefore define nodes with a short transition time as highly epileptogenic. In a retrospective study, we found that in some patients the regions with high epileptogenicity in the model overlap with those identified clinically as the seizure onset zone. Moreover, it was found that the resection of these regions in the model reduces the overall likelihood of a seizure. Following removal of these regions in the model, we predicted surgical outcomes and compared these to actual patient outcomes. Our predictions were found to be 81.3% accurate on a dataset of 16 patients with intractable epilepsy. Intriguingly, in patients with unsuccessful outcomes, the proposed computational approach is able to suggest alternative resection sites. The model presented here gives mechanistic insight as to why surgery may be unsuccessful in some patients. This may aid clinicians in presurgical evaluation by providing a tool to explore various surgical options, offering complementary information to existing clinical techniques.

Estimation of brain network ictogenicity predicts outcome from epilepsy surgery

Surgery is a valuable option for pharmacologically intractable epilepsy. However, significant post-operative improvements are not always attained. This is due in part to our incomplete understanding of the seizure generating (ictogenic) capabilities of brain networks. Here we introduce an in silico, model-based framework to study the effects of surgery within ictogenic brain networks. We find that factors conventionally determining the region of tissue to resect, such as the location of focal brain lesions or the presence of epileptiform rhythms, do not necessarily predict the best resection strategy. We validate our framework by analysing electrocorticogram (ECoG) recordings from patients who have undergone epilepsy surgery. We find that when post-operative outcome is good, model predictions for optimal strategies align better with the actual surgery undertaken than when post-operative outcome is poor. Crucially, this allows the prediction of optimal surgical strategies and the provision of quantitative prognoses for patients undergoing epilepsy surgery.

Epileptogenic zone localization and seizure control in coupled neural mass models

Exact localization of the epileptogenic zone (EZ) is the first priority for ensuring epilepsy treatments and reducing side effects. The results of traditional visual methods for localizing the origin of seizures are far from satisfactory in some cases. Signal processing methods could extract substantial information that may complement visual inspection of EEG signals. In this study, EZ localization is changed into a driver identification problem, and a nonlinear interdependence measure, the weighted rank interdependence, is proposed and used as a driver indicator because it can detect coupling information, especially directionality, from EEG signals. A proportional integral derivative (PID) controller is then explored, using simulations, to establish its suitability for seizure control. The seizure control we propose rests on identifying the EZ using nonlinear interdependence measures of directed functional connectivity. Two directionally coupled neural mass models are employed for simulation investigation. Two parameters can adjust the sensitivity and completeness of the weighted rank interdependence for different applications, and their effect is discussed in the context of neural mass models. Simulation results demonstrate that use of the weighted rank interdependence for EZ identification can be applied to different EZ types, and the approach achieves an overall identification rate of 98.84 % for several EZ types. Simulations also indicate that PID control can effectively regulate synchronization between neural masses.

Reflex seizures, traits, and epilepsies: from physiology to pathology

Epileptic seizures are generally unpredictable and arise spontaneously. Patients often report non-specific triggers such as stress or sleep deprivation, but only rarely do seizures occur as a reflex event, in which they are objectively and consistently modulated, precipitated, or inhibited by external sensory stimuli or specific cognitive processes. The seizures triggered by such stimuli and processes in susceptible individuals can have different latencies. Once seizure-suppressing mechanisms fail and a critical mass (the so-called tipping point) of cortical activation is reached, reflex seizures stereotypically manifest with common motor features independent of the physiological network involved. The complexity of stimuli increases from simple sensory to complex cognitive-emotional with increasing age of onset. The topography of physiological networks involved follows the posterior-to-anterior trajectory of brain development, reflecting age-related changes in brain excitability. Reflex seizures and traits probably represent the extremes of a continuum, and understanding of their underlying mechanisms might help to elucidate the transition of normal physiological function to paroxysmal epileptic activity.

White Matter Connectivity of the Thalamus Delineates the Functional Architecture of Competing Thalamocortical Systems

There is an increasing awareness of the involvement of thalamic connectivity on higher level cortical functioning in the human brain. This is reflected by the influence of thalamic stimulation on cortical activity and behavior as well as apparently cortical lesion syndromes occurring as a function of small thalamic insults. Here, we attempt to noninvasively test the correspondence of structural and functional connectivity of the human thalamus using diffusion-weighted and resting-state functional MRI. Using a large sample of 102 adults, we apply tensor independent component analysis to diffusion MRI tractography data to blindly parcellate bilateral thalamus according to diffusion tractography-defined structural connectivity. Using resting-state functional MRI collected in the same subjects, we show that the resulting structurally defined thalamic regions map to spatially distinct, and anatomically predictable, whole-brain functional networks in the same subjects. Although there was significant variability in the functional connectivity patterns, the resulting 51 structural and functional patterns could broadly be reduced to a subset of 7 similar core network types. These networks were distinct from typical cortical resting-state networks. Importantly, these networks were distributed across the brain and, in a subset, map extremely well to known thalamocortico-basal-ganglial loops.