The spatiotemporal dynamics of human focal seizures

Spontaneous human focal seizures can present with a plethora of behavioral manifestations that vary according to the affected cortical regions; however, several key features have been consistently observed. During my doctoral studies, I applied both theoretical and experimental methods to study mechanisms underpinning these consistently seen dynamics. I first analyzed human intracranial EEG recordings, describing statistical methods for measuring their electrophysiological signatures. I next proposed several neurophysiological hypotheses that could explain seizure dynamics and verified them in rodent seizure models. Finally, a computational model was developed, successfully explaining how the complex spatiotemporal evolution of focal seizures emerges from simple neurophysiological principles. In Chapter 1, the long-standing behavioral manifestations and the most up-to-date electrophysiology findings are reviewed. This section details the inspiration for the studies reported in the subsequent chapters. In Chapter 2, I describe several statistical methods for estimating traveling wave velocities. I show most ictal discharges can be described as traveling waves whose velocities contain rich information about the stages of seizure evolution. I compare performance of various statistical methods and propose a robust approach to boost the quality of each method’s estimation results. In Chapter 3, I show how inhibition modulates seizure propagation patterns. Surround inhibition spatially restrains focal seizures and masks excitatory projections of ictal activities. When compromised, two patterns of seizure propagation emerge according to the position of inhibition defects relative to the ictal focus. I show that two distant ictal foci can communicate via physiological connectivity without any chronic rewiring processes – confirming the existence of long-range propagation pathways that could lead to epileptic network formation. In Chapter 4, I show that thalamic inputs might be necessary for interictal epileptiform discharges (IEDs). The relative positions between IEDs and ictal foci indicate that surround inhibition, shown in the previous chapter, can be exhausted by repetitive exposure to ictal projections. In Chapter 5, I propose a neural network model that can explain both long-standing behavioral observations of seizures and account for the most up-to-date electrophysiological recordings of spontaneous human focal seizures. The model relies on few assumptions, all of which are proved or supported in earlier chapters of this thesis. The model explains phasic evolution of seizure dynamics – how the commonly observed patterns arise from simple neurophysiological principles, as well as seizure onset subtypes, traveling wave directions and speeds. The model also predicts how spontaneous seizures might arise from synaptic plasticity. The chapter ends with a discussion of the model’s implications and future work. The thesis is organized in a way that each chapter can be read independently, with Chapter 5 summarizing the central theory spanning the whole study. Each chapter is also tightly linked to a clinically relevant question. In sum, the dissertation’s goal is to provide an in-principle understanding of focal seizure dynamics. With rapid advancement of clinical and experimental tools, I believe this work provides a roadmap for future therapies for epilepsy patients

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 in-variant across seizures in a same patient, the source of traveling (2-3 Hz) spike-and-wave discharges (SWDs) during seizures can either move with the slower propagating ictal wavefront or remain stationary at the seizure onset area. In addition, although most focal seizures terminate quasi-synchronously across brain areas, some evolve into distinct ictal clusters and terminate asynchronously. To provide a unifying perspective on the observed diversity of spatiotemporal dynamics for seizure spread and termination, we introduce here the Epileptor neural field model. Two mechanisms play an essential role. First, while the slow ictal wavefront propagates as a front in excitable neural media, the faster SWDs propagation results from coupled-oscillator dynamics. Second, multiple time scales interact during seizure spread, allowing for low-voltage fast-activity (>10 Hz) to hamper seizure spread and for SWD propagation to affect the way a seizure terminates. These dynamics, together with variations in short and long-range connectivity strength, play a central role on seizure spread, maintenance and termination. We demonstrate how Epileptor field models incorporating the above mechanisms predict the previously reported diversity in seizure spread patterns. Furthermore, we confirm the predictions for synchronous or asynchronous (clustered) seizure termination in human seizures recorded via stereotactic EEG. Our new insights into seizure spatiotemporal dynamics may also contribute to the development of new closed-loop neuromodulation therapies for focal epilepsy.Comment: 10 pages + 9 pages Supporting Information (SI), 7 figures, 1 SI table, 7 SI figure

Involvement of fast-spiking cells in ictal sequences during spontaneous seizures in rats with chronic temporal lobe epilepsy

Epileptic seizures represent altered neuronal network dynamics, but the temporal evolution and cellular substrates of the neuronal activity patterns associated with spontaneous seizures are not fully understood. We used simultaneous recordings from multiple neurons in the hippocampus and neocortex of rats with chronic temporal lobe epilepsy to demonstrate that subsets of cells discharge in a highly stereotypical sequential pattern during ictal events, and that these stereotypical patterns were reproducible across consecutive seizures. In contrast to the canonical view that principal cell discharges dominate ictal events, the ictal sequences were predominantly composed of fast-spiking, putative inhibitory neurons, which displayed unusually strong coupling to local field potential even before seizures. The temporal evolution of activity was characterized by unique dynamics where the most correlated neuronal pairs before seizure onset displayed the largest increases in correlation strength during the seizures. These results demonstrate the selective involvement of fast spiking interneurons in structured temporal sequences during spontaneous ictal events in hippocampal and neocortical circuits in experimental models of chronic temporal lobe epilepsy

The Dynamic Brain: From Spiking Neurons to Neural Masses and Cortical Fields

The cortex is a complex system, characterized by its dynamics and architecture, which underlie many functions such as action, perception, learning, language, and cognition. Its structural architecture has been studied for more than a hundred years; however, its dynamics have been addressed much less thoroughly. In this paper, we review and integrate, in a unifying framework, a variety of computational approaches that have been used to characterize the dynamics of the cortex, as evidenced at different levels of measurement. Computational models at different space–time scales help us understand the fundamental mechanisms that underpin neural processes and relate these processes to neuroscience data. Modeling at the single neuron level is necessary because this is the level at which information is exchanged between the computing elements of the brain; the neurons. Mesoscopic models tell us how neural elements interact to yield emergent behavior at the level of microcolumns and cortical columns. Macroscopic models can inform us about whole brain dynamics and interactions between large-scale neural systems such as cortical regions, the thalamus, and brain stem. Each level of description relates uniquely to neuroscience data, from single-unit recordings, through local field potentials to functional magnetic resonance imaging (fMRI), electroencephalogram (EEG), and magnetoencephalogram (MEG). Models of the cortex can establish which types of large-scale neuronal networks can perform computations and characterize their emergent properties. Mean-field and related formulations of dynamics also play an essential and complementary role as forward models that can be inverted given empirical data. This makes dynamic models critical in integrating theory and experiments. We argue that elaborating principled and informed models is a prerequisite for grounding empirical neuroscience in a cogent theoretical framework, commensurate with the achievements in the physical sciences

Dynamics and precursor signs for phase transitions in neural systems

This thesis investigates neural state transitions associated with sleep, seizure and anaesthesia. The aim is to address the question: How does a brain traverse the critical threshold between distinct cortical states, both healthy and pathological? Specifically we are interested in sub-threshold neural behaviour immediately prior to state transition. We use theoretical neural modelling (single spiking neurons, a network of these, and a mean-field continuum limit) and in vitro experiments to address this question. Dynamically realistic equations of motion for thalamic relay neuron, reticular nuclei, cortical pyramidal and cortical interneuron in different vigilance states are developed, based on the Izhikevich spiking neuron model. A network of cortical neurons is assembled to examine the behaviour of the gamma-producing cortical network and its transition to lower frequencies due to effect of anaesthesia. Then a three-neuron model for the thalamocortical loop for sleep spindles is presented. Numerical simulations of these networks confirms spiking consistent with reported in vivo measurement results, and provides supporting evidence for precursor indicators of imminent phase transition due to occurrence of individual spindles. To complement the spiking neuron networks, we study the Wilson–Cowan neural mass equations describing homogeneous cortical columns and a 1D spatial cluster of such columns. The abstract representation of cortical tissue by a pair of coupled integro-differential equations permits thorough linear stability, phase plane and bifurcation analyses. This model shows a rich set of spatial and temporal bifurcations marking the boundary to state transitions: saddle-node, Hopf, Turing, and mixed Hopf–Turing. Close to state transition, white-noise-induced subthreshold fluctuations show clear signs of critical slowing down with prolongation and strengthening of autocorrelations, both in time and space, irrespective of bifurcation type. Attempts at in vitro capture of these predicted leading indicators form the last part of the thesis. We recorded local field potentials (LFPs) from cortical and hippocampal slices of mouse brain. State transition is marked by the emergence and cessation of spontaneous seizure-like events (SLEs) induced by bathing the slices in an artificial cerebral spinal fluid containing no magnesium ions. Phase-plane analysis of the LFP time-series suggests that distinct bifurcation classes can be responsible for state change to seizure. Increased variance and growth of spectral power at low frequencies (f < 15 Hz) was observed in LFP recordings prior to initiation of some SLEs. In addition we demonstrated prolongation of electrically evoked potentials in cortical tissue, while forwarding the slice to a seizing regime. The results offer the possibility of capturing leading temporal indicators prior to seizure generation, with potential consequences for understanding epileptogenesis. Guided by dynamical systems theory this thesis captures evidence for precursor signs of phase transitions in neural systems using mathematical and computer-based modelling as well as in vitro experiments

Analysis of traveling wave propagation in one-dimensional integrate-and-fire neural networks

One-dimensional neural networks comprised of large numbers of Integrate-and-Fire neurons have been widely used to model electrical activity propagation in neural slices. Despite these efforts, the vast majority of these computational models have no analytical solutions. Consequently, my Ph.D. research focuses on a specific class of homogeneous Integrate-and-Fire neural network, for which analytical solutions of network dynamics can be derived. One crucial analytical finding is that the traveling wave acceleration quadratically depends on the instantaneous speed of the activity propagation, which means that two speed solutions exist in the activities of wave propagation: one is fast-stable and the other is slow-unstable. Furthermore, via this property, we analytically compute temporal-spatial spiking dynamics to help gain insights into the stability mechanisms of traveling wave propagation. Indeed, the analytical solutions are in perfect agreement with the numerical solutions. This analytical method also can be applied to determine the effects induced by a non-conductive gap of brain tissue and extended to more general synaptic connectivity functions, by converting the evolution equations for network dynamics into a low-dimensional system of ordinary differential equations. Building upon these results, we investigate how periodic inhomogeneities affect the dynamics of activity propagation. In particular, two types of periodic inhomogeneities are studied: alternating regions of additional fixed excitation and inhibition, and cosine form inhomogeneity. Of special interest are the conditions leading to propagation failure. With similar analytical procedures, explicit expressions for critical speeds of activity propagation are obtained under the influence of additional inhibition and excitation. However, an explicit formula for speed modulations is difficult to determine in the case of cosine form inhomogeneity. Instead of exact solutions from the system of equations, a series of speed approximations are constructed, rendering a higher accuracy with a higher order approximation of speed

Capture of fixation by rotational flow; a deterministic hypothesis regarding scaling and stochasticity in fixational eye movements.

Visual scan paths exhibit complex, stochastic dynamics. Even during visual fixation, the eye is in constant motion. Fixational drift and tremor are thought to reflect fluctuations in the persistent neural activity of neural integrators in the oculomotor brainstem, which integrate sequences of transient saccadic velocity signals into a short term memory of eye position. Despite intensive research and much progress, the precise mechanisms by which oculomotor posture is maintained remain elusive. Drift exhibits a stochastic statistical profile which has been modeled using random walk formalisms. Tremor is widely dismissed as noise. Here we focus on the dynamical profile of fixational tremor, and argue that tremor may be a signal which usefully reflects the workings of oculomotor postural control. We identify signatures reminiscent of a certain flavor of transient neurodynamics; toric traveling waves which rotate around a central phase singularity. Spiral waves play an organizational role in dynamical systems at many scales throughout nature, though their potential functional role in brain activity remains a matter of educated speculation. Spiral waves have a repertoire of functionally interesting dynamical properties, including persistence, which suggest that they could in theory contribute to persistent neural activity in the oculomotor postural control system. Whilst speculative, the singularity hypothesis of oculomotor postural control implies testable predictions, and could provide the beginnings of an integrated dynamical framework for eye movements across scales

Spatio-temporal modelling and analysis of epileptiform EEG

In this thesis we investigate the mechanisms underlying the generation of abnormal EEG rhythms in epilepsy, which is a crucial step towards better treatment of this disorder in the future. To this end, macroscopic scale mathematical models of the interactions between neuronal populations are examined. In particular, the role of interactions between neural masses that are spatially distributed in cortical networks are explored. In addition, two other important aspects of the modelling process are addressed, namely the conversion of macroscopic model variables into EEG output and the comparison of multivariate, spatio-temporal data. For the latter, we adopt a vectorisation of the correlation matrix of windowed data and subsequent comparison of data by vector distance measures. Our modelling studies indicate that excitatory connectivity between neural masses facilitates self-organised dynamics. In particular, we report for the first time the production of complex rhythmic transients and the generation of intermittent periods of 'abnormal' rhythmic activity in two different models of epileptogenic tissue. These models therefore provide novel accounts of the spontaneous, intermittent transition between normal and pathological rhythms in primarily generalised epilepsies and the evocation of complex, self-terminating, spatio-temporal dynamics by brief stimulation in focal epilepsies. Two key properties of these models are excitability at the macroscopic level and the presence of spatial heterogeneities. The identification of neural mass excitability as an important processes in spatially extended brain networks is a step towards uncovering the multi-scale nature of the pathological mechanisms of epilepsy. A direct consequence of this work is therefore that novel experimental investigations are proposed, which in itself is a validation of our modelling approach. In addition, new considerations regarding the nature of dynamical systems as applied to problems of transitions between rhythmic states are proposed and will prompt future investigations of complex transients in spatio-temporal excitable systems.EThOS - Electronic Theses Online ServiceGBUnited Kingdo