Kyle Q. Lepage

Human seizures self-terminate across spatial scales via a critical transition

Why seizures spontaneously terminate remains an unanswered fundamental question of epileptology. Here we present evidence that seizures self-terminate via a discontinuous critical transition or bifurcation. We show that human brain electrical activity at various spatial scales exhibits common dynamical signatures of an impending critical transition—slowing, increased correlation, and flickering—in the approach to seizure termination. In contrast, prolonged seizures (status epilepticus) repeatedly approach, but do not cross, the critical transition. To support these results, we implement a computational model that demonstrates that alternative stable attractors, representing the ictal and postictal states, emulate the observed dynamics. These results suggest that self-terminating seizures end through a common dynamical mechanism. This description constrains the specific biophysical mechanisms underlying seizure termination, suggests a dynamical understanding of status epilepticus, and demonstrates an accessible system for studying critical transitions in nature.

Emergence of Persistent Networks in Long-Term Intracranial EEG Recordings

Over the past two decades, the increased ability to analyze network relationships among neural structures has provided novel insights into brain function. Most network approaches, however, focus on static representations of the brains physical or statistical connectivity. Few studies have examined how brain functional networks evolve spontaneously over long epochs of continuous time. To address this, we examine functional connectivity networks deduced from continuous long-term electrocorticogram recordings. For a population of six human patients, we identify a persistent pattern of connections that form a frequency-band-dependent network template, and a set of core connections that appear frequently and together. These structures are robust, emerging from brief time intervals (∼100 s) regardless of cognitive state. These results suggest that a metastable, frequency-band-dependent scaffold of brain connectivity exists from which transient activity emerges and recedes.

Age-Related Changes in 1/f Neural Electrophysiological Noise

Aging is associated with performance decrements across multiple cognitive domains. The neural noise hypothesis, a dominant view of the basis of this decline, posits that aging is accompanied by an increase in spontaneous, noisy baseline neural activity. Here we analyze data from two different groups of human subjects: intracranial electrocorticography from 15 participants over a 38 year age range (15–53 years) and scalp EEG data from healthy younger (20–30 years) and older (60–70 years) adults to test the neural noise hypothesis from a 1/f noise perspective. Many natural phenomena, including electrophysiology, are characterized by 1/f noise. The defining characteristic of 1/f is that the power of the signal frequency content decreases rapidly as a function of the frequency (f) itself. The slope of this decay, the noise exponent (χ), is often <−1 for electrophysiological data and has been shown to approach white noise (defined as χ = 0) with increasing task difficulty. We observed, in both electrophysiological datasets, that aging is associated with a flatter (more noisy) 1/f power spectral density, even at rest, and that visual cortical 1/f noise statistically mediates age-related impairments in visual working memory. These results provide electrophysiological support for the neural noise hypothesis of aging. SIGNIFICANCE STATEMENT Understanding the neurobiological origins of age-related cognitive decline is of critical scientific, medical, and public health importance, especially considering the rapid aging of the worlds population. We find, in two separate human studies, that 1/f electrophysiological noise increases with aging. In addition, we observe that this age-related 1/f noise statistically mediates age-related working memory decline. These results significantly add to this understanding and contextualize a long-standing problem in cognition by encapsulating age-related cognitive decline within a neurocomputational model of 1/f noise-induced deficits in neural communication.

The dependence of spike field coherence on expected intensity

The coherence between neural spike trains and local-field potential recordings, called spike-field coherence, is of key importance in many neuroscience studies. In this work, aside from questions of estimator performance, we demonstrate that theoretical spike-field coherence for a broad class of spiking models depends on the expected rate of spiking. This rate dependence confounds the phase locking of spike events to field-potential oscillations with overall neuron activity and is demonstrated analytically, for a large class of stochastic models, and in simulation. Finally, the relationship between the spike-field coherence and the intensity field coherence is detailed analytically. This latter quantity is independent of neuron firing rate and, under commonly found conditions, is proportional to the probability that a neuron spikes at a specific phase of field oscillation. Hence, intensity field coherence is a rate-independent measure and a candidate on which to base the appropriate statistical inference of spike field synchrony.

Slow spatial recruitment of neocortex during secondarily generalized seizures and its relation to surgical outcome

Understanding the spatiotemporal dynamics of brain activity is crucial for inferring the underlying synaptic and nonsynaptic mechanisms of brain dysfunction. Focal seizures with secondary generalization are traditionally considered to begin in a limited spatial region and spread to connected areas, which can include both pathological and normal brain tissue. The mechanisms underlying this spread are important to our understanding of seizures and to improve therapies for surgical intervention. Here we study the properties of seizure recruitment—how electrical brain activity transitions to large voltage fluctuations characteristic of spike-and-wave seizures. We do so using invasive subdural electrode arrays from a population of 16 patients with pharmacoresistant epilepsy. We find an average delay of ∼30 s for a broad area of cortex (8 × 8 cm) to be recruited into the seizure, at an estimated speed of ∼4 mm/s. The spatiotemporal characteristics of recruitment reveal two categories of patients: one in which seizure recruitment of neighboring cortical regions follows a spatially organized pattern consistent from seizure to seizure, and a second group without consistent spatial organization of activity during recruitment. The consistent, organized recruitment correlates with a more regular, compared with small-world, connectivity pattern in simulation and successful surgical treatment of epilepsy. We propose that an improved understanding of how the seizure recruits brain regions into large amplitude voltage fluctuations provides novel information to improve surgical treatment of epilepsy and highlights the slow spread of massive local activity across a vast extent of cortex during seizure.

Inferring evoked brain connectivity through adaptive perturbation

Inference of functional networks—representing the statistical associations between time series recorded from multiple sensors—has found important applications in neuroscience. However, networksexhibiting time-locked activity between physically independent elements can bias functional connectivity estimates employing passive measurements. Here, a perturbative and adaptive method of inferring network connectivity based on measurement and stimulation—so called “evoked network connectivity” is introduced. This procedure, employing a recursive Bayesian update scheme, allows principled network stimulation given a current network estimate inferred from all previous stimulations and recordings. The method decouples stimulus and detector design from network inference and can be suitably applied to a wide range of clinical and basic neuroscience related problems. The proposed method demonstrates improved accuracy compared to network inference based on passive observation of node dynamics and an increased rate of convergence relative to network estimation employing a naïve stimulation strategy.

A statistically robust EEG re-referencing procedure to mitigate reference effect

BACKGROUND The electroencephalogram (EEG) remains the primary tool for diagnosis of abnormal brain activity in clinical neurology and for in vivo recordings of human neurophysiology in neuroscience research. In EEG data acquisition, voltage is measured at positions on the scalp with respect to a reference electrode. When this reference electrode responds to electrical activity or artifact all electrodes are affected. Successful analysis of EEG data often involves re-referencing procedures that modify the recorded traces and seek to minimize the impact of reference electrode activity upon functions of the original EEG recordings. NEW METHOD We provide a novel, statistically robust procedure that adapts a robust maximum-likelihood type estimator to the problem of reference estimation, reduces the influence of neural activity from the re-referencing operation, and maintains good performance in a wide variety of empirical scenarios. RESULTS The performance of the proposed and existing re-referencing procedures are validated in simulation and with examples of EEG recordings. To facilitate this comparison, channel-to-channel correlations are investigated theoretically and in simulation. COMPARISON WITH EXISTING METHODS The proposed procedure avoids using data contaminated by neural signal and remains unbiased in recording scenarios where physical references, the common average reference (CAR) and the reference estimation standardization technique (REST) are not optimal. CONCLUSION The proposed procedure is simple, fast, and avoids the potential for substantial bias when analyzing low-density EEG data.

Frontal beta-theta network during REM sleep

We lack detailed knowledge about the spatio-temporal physiological signatures of REM sleep, especially in humans. By analyzing intracranial electrode data from humans, we demonstrate for the first time that there are prominent beta (15–35 Hz) and theta (4–8 Hz) oscillations in both the anterior cingulate cortex (ACC) and the DLPFC during REM sleep. We further show that these theta and beta activities in the ACC and the DLPFC, two relatively distant but reciprocally connected regions, are coherent. These findings suggest that, counter to current prevailing thought, the DLPFC is active during REM sleep and likely interacting with other areas. Since the DLPFC and the ACC are implicated in memory and emotional regulation, and the ACC has motor areas and is thought to be important for error detection, the dialogue between these two areas could play a role in the regulation of emotions and in procedural motor and emotional memory consolidation. DOI: http://dx.doi.org/10.7554/eLife.18894.001

An Approach to Time-Frequency Analysis With Ridges of the Continuous Chirplet Transform

We propose a time-frequency representation based on the ridges of the continuous chirplet transform to identify both fast transients and components with well-defined instantaneous frequency in noisy data. At each chirplet modulation rate, every ridge corresponds to a territory in the time-frequency plane such that the territories form a partition of the time-frequency plane. For many signals containing sparse signal components, ridge length and ridge stability are maximized when the analysis kernel is adapted to the signal content. These properties provide opportunities for enhancing signal in noise and for elementary stream segregation by combining information across multiple analysis chirp rates.

Some sampling properties of common phase estimators

The instantaneous phase of neural rhythms is important to many neuroscience-related studies. In this letter, we show that the statistical sampling properties of three instantaneous phase estimators commonly employed to analyze neuroscience data share common features, allowing an analytical investigation into their behavior. These three phase estimators—the Hilbert, complex Morlet, and discrete Fourier transform—are each shown to maximize the likelihood of the data, assuming the observation of different neural signals. This connection, explored with the use of a geometric argument, is used to describe the bias and variance properties of each of the phase estimators, their temporal dependence, and the effect of model misspecification. This analysis suggests how prior knowledge about a rhythmic signal can be used to improve the accuracy of phase estimates.