Levin Kuhlmann

Patient-specific bivariate-synchrony-based seizure prediction for short prediction horizons

This paper evaluates the patient-specific seizure prediction performance of pre-ictal changes in bivariate-synchrony between pairs of intracranial electroencephalographic (iEEG) signals within 15min of a seizure in patients with pharmacoresistant focal epilepsy. Prediction horizons under 15min reduce the durations of warning times and should provide adequate time for a seizure control device to intervene. Long-term continuous iEEG was obtained from 6 patients. The seizure prediction performance was evaluated for all possible channel pairs and for different prediction methods to find the best performing channel pairs and methods for both pre-ictal decreases and increases in synchrony. The different prediction methods involved changes in window duration, signal filtering, thresholding approach, and prediction horizon durations. Performance for each patient, for all seizures, was first compared with an analytical-Poisson-based random predictor. The performance of the top 5% of channel pairs for each patient closely matched the top 5% of analytical-Poisson-based random predictor performance indicating that patient-specific, bivariate-synchrony-based seizure prediction could be random in general (under the assumption that channel-pair prediction times are statistically independent). Analysis of the spatial patterns of performance showed no clear relationship to the seizure onset zone. For each patient the best channel pair showed better performance than Poisson-based random prediction for a selected subset of prediction thresholds. Given the caveats of comparing with this form of random prediction, alarm time surrogates were employed to assess statistical significance of a four-fold out-of-sample cross-validation analysis applied to the best channel-pairs. The cross-validation analysis obtained reasonable testing performance for most patients when performance was compared to random prediction based on alarm time surrogates. The most significant case was a patient whose testing set sensitivity and false positive rate were 0.67±0.09 and 3.04±0.29h(-1), respectively, for decreases in synchrony, an intervention time of 15min and a seizure onset period of 5min. For each testing set for this patient, performance was better than that obtained by random prediction at the significance level of 0.05 (average sensitivity of 0.47±0.05). Moreover, there were 9 seizures in each testing set which gives greater power to this cross-validation result, although the cross-validation was performed on the best channel pair selected by within-sample optimization for all seizures of the patient. Further validation with larger datasets from individual patients is needed. Improvements in prediction performance should be achievable through investigations of multivariate synchrony combined with non-linear classification methods.

A neural model of 3D shape-from-texture: Multiple-scale filtering, boundary grouping, and surface filling-in

A neural model is presented of how cortical areas V1, V2, and V4 interact to convert a textured 2D image into a representation of curved 3D shape. Two basic problems are solved to achieve this: (1) Patterns of spatially discrete 2D texture elements are transformed into a spatially smooth surface representation of 3D shape. (2) Changes in the statistical properties of texture elements across space induce the perceived 3D shape of this surface representation. This is achieved in the model through multiple-scale filtering of a 2D image, followed by a cooperative-competitive grouping network that coherently binds texture elements into boundary webs at the appropriate depths using a scale-to-depth map and a subsequent depth competition stage. These boundary webs then gate filling-in of surface lightness signals in order to form a smooth 3D surface percept. The model quantitatively simulates challenging psychophysical data about perception of prolate ellipsoids [Todd, J., & Akerstrom, R. (1987). Perception of three-dimensional form from patterns of optical texture. Journal of Experimental Psychology: Human Perception and Performance, 13(2), 242-255]. In particular, the model represents a high degree of 3D curvature for a certain class of images, all of whose texture elements have the same degree of optical compression, in accordance with percepts of human observers. Simulations of 3D percepts of an elliptical cylinder, a slanted plane, and a photo of a golf ball are also presented.

Summation of spatiotemporal input patterns in leaky integrate-and-fire neurons: application to neurons in the cochlear nucleus receiving converging auditory nerve fiber input

The response of leaky integrate-and-fire neurons is analyzed for periodic inputs whose phases vary with their spatial location. The model gives the relationship between the spatial summation distance and the degree of phase locking of the output spikes (i.e., locking to the periodic stochastic inputs, measured by the synchronization index). The synaptic inputs are modeled as an inhomogeneous Poisson process, and the analysis is carried out in the Gaussian approximation. The model has been applied to globular bushy cells of the cochlear nucleus, which receive converging inputs from auditory nerve fibers that originate at neighboring sites in the cochlea. The model elucidates the roles played by spatial summation and coincidence detection, showing how synchronization decreases with an increase in both frequency and spatial spread of inputs. It also shows under what conditions an enhancement of synchronization of the output relative to the input takes place.

Seizure Prediction: Science Fiction or Soon to Become Reality?

This review highlights recent developments in the field of epileptic seizure prediction. We argue that seizure prediction is possible; however, most previous attempts have used data with an insufficient amount of information to solve the problem. The review discusses four methods for gaining more information above standard clinical electrophysiological recordings. We first discuss developments in obtaining long-term data that enables better characterisation of signal features and trends. Then, we discuss the usage of electrical stimulation to probe neural circuits to obtain robust information regarding excitability. Following this, we present a review of developments in high-resolution micro-electrode technologies that enable neuroimaging across spatial scales. Finally, we present recent results from data-driven model-based analyses, which enable imaging of seizure generating mechanisms from clinical electrophysiological measurements. It is foreseeable that the field of seizure prediction will shift focus to a more probabilistic forecasting approach leading to improvements in the quality of life for the millions of people who suffer uncontrolled seizures. However, a missing piece of the puzzle is devices to acquire long-term high-quality data. When this void is filled, seizure prediction will become a reality.

Role of multiple-scale modeling of epilepsy in seizure forecasting

Summary: Over the past three decades, a number of seizure prediction, or forecasting, methods have been developed. Although major achievements were accomplished regarding the statistical evaluation of proposed algorithms, it is recognized that further progress is still necessary for clinical application in patients. The lack of physiological motivation can partly explain this limitation. Therefore, a natural question is raised: can computational models of epilepsy be used to improve these methods? Here, we review the literature on the multiple-scale neural modeling of epilepsy and the use of such models to infer physiologic changes underlying epilepsy and epileptic seizures. The authors argue how these methods can be applied to advance the state-of-the-art in seizure forecasting.

A Computational Study of How Orientation Bias in the Lateral Geniculate Nucleus Can Give Rise to Orientation Selectivity in Primary Visual Cortex

Controversy remains about how orientation selectivity emerges in simple cells of the mammalian primary visual cortex. In this paper, we present a computational model of how the orientation-biased responses of cells in lateral geniculate nucleus (LGN) can contribute to the orientation selectivity in simple cells in cats. We propose that simple cells are excited by lateral geniculate fields with an orientation-bias and disynaptically inhibited by unoriented lateral geniculate fields (or biased fields pooled across orientations), both at approximately the same retinotopic co-ordinates. This interaction, combined with recurrent cortical excitation and inhibition, helps to create the sharp orientation tuning seen in simple cell responses. Along with describing orientation selectivity, the model also accounts for the spatial frequency and length–response functions in simple cells, in normal conditions as well as under the influence of the GABAA antagonist, bicuculline. In addition, the model captures the response properties of LGN and simple cells to simultaneous visual stimulation and electrical stimulation of the LGN. We show that the sharp selectivity for stimulus orientation seen in primary visual cortical cells can be achieved without the excitatory convergence of the LGN input cells with receptive fields along a line in visual space, which has been a core assumption in classical models of visual cortex. We have also simulated how the full range of orientations seen in the cortex can emerge from the activity among broadly tuned channels tuned to a limited number of optimum orientations, just as in the classical case of coding for color in trichromatic primates.

Mechanisms of seizure propagation in 2-dimensional centre-surround recurrent networks.

Understanding how seizures spread throughout the brain is an important problem in the treatment of epilepsy, especially for implantable devices that aim to avert focal seizures before they spread to, and overwhelm, the rest of the brain. This paper presents an analysis of the speed of propagation in a computational model of seizure-like activity in a 2-dimensional recurrent network of integrate-and-fire neurons containing both excitatory and inhibitory populations and having a difference of Gaussians connectivity structure, an approximation to that observed in cerebral cortex. In the same computational model network, alternative mechanisms are explored in order to simulate the range of seizure-like activity propagation speeds (0.1–100 mm/s) observed in two animal-slice-based models of epilepsy: (1) low extracellular , which creates excess excitation and (2) introduction of gamma-aminobutyric acid (GABA) antagonists, which reduce inhibition. Moreover, two alternative connection topologies are considered: excitation broader than inhibition, and inhibition broader than excitation. It was found that the empirically observed range of propagation velocities can be obtained for both connection topologies. For the case of the GABA antagonist model simulation, consistent with other studies, it was found that there is an effective threshold in the degree of inhibition below which waves begin to propagate. For the case of the low extracellular model simulation, it was found that activity-dependent reductions in inhibition provide a potential explanation for the emergence of slowly propagating waves. This was simulated as a depression of inhibitory synapses, but it may also be achieved by other mechanisms. This work provides a localised network understanding of the propagation of seizures in 2-dimensional centre-surround networks that can be tested empirically.

A robust circle criterion observer with application to neural mass models

A robust circle criterion observer is designed and applied to neural mass models. At present, no existing circle criterion observers apply to the considered models, i.e. the required linear matrix inequality is infeasible. Therefore, we generalise available results to derive a suitable estimation algorithm. Additionally, the design also takes into account input uncertainty and measurement noise. We show how to apply the observer to estimate the mean membrane potential of neuronal populations of a popular single cortical column model from the literature.

Modulation of functional EEG networks by the NMDA antagonist nitrous oxide

Parietal networks are hypothesised to play a central role in the cortical information synthesis that supports conscious experience and behavior. Significant reductions in parietal level functional connectivity have been shown to occur during general anesthesia with propofol and a range of other GABAergic general anesthetic agents. Using two analysis approaches (1) a graph theoretic analysis based on surrogate-corrected zero-lag correlations of scalp EEG, and (2) a global coherence analysis based on the EEG cross-spectrum, we reveal that sedation with the NMDA receptor antagonist nitrous oxide (N2O), an agent that has quite different electroencephalographic effects compared to the inductive general anesthetics, also causes significant alterations in parietal level functional networks, as well as changes in full brain and frontal level networks. A total of 20 subjects underwent N2O inhalation at either 20%, 40% or 60% peak N2O/O2 gas concentration levels. N2O-induced reductions in parietal network level functional connectivity (on the order of 50%) were exclusively detected by utilising a surface Laplacian derivation, suggesting that superficial, smaller spatial scale, cortical networks were most affected. In contrast reductions in frontal network functional connectivity were optimally discriminated using a common-reference derivation (reductions on the order of 10%), indicating that the NMDA antagonist N2O induces spatially coherent and widespread perturbations in frontal activity. Our findings not only give important weight to the idea of agent invariant final network changes underlying drug-induced reductions in consciousness, but also provide significant impetus for the application and development of multiscale functional analyses to systematically characterise the network level cortical effects of NMDA receptor related hypofunction. Future work at the source space level will be needed to verify the consistency between cortical network changes seen at the source level and those presented here at the EEG sensor space level.

Parameter and State Estimation of Nonlinear Systems Using a Multi-Observer Under the Supervisory Framework

We present a hybrid scheme for the parameter and state estimation of nonlinear continuous-time systems, which is inspired by the supervisory setup used for control. State observers are synthesized for some nominal parameter values and a criterion is designed to select one of these observers at any given time instant, which provides state and parameter estimates. Assuming that a persistency of excitation condition holds, the convergence of the parameter and state estimation errors to zero is ensured up to a margin, which can be made as small as desired by increasing the number of observers. To reduce the potential computational complexity of the scheme, we explain how the sampling of the parameter set can be dynamically updated using a zoom-in procedure. This strategy typically requires a fewer number of observers for a given estimation error margin compared to the static sampling policy. The results are shown to be applicable to linear systems and to a class of nonlinear systems. We illustrate the applicability of the approach by estimating the synaptic gains and the mean membrane potentials of a neural mass model.