Samuel P. Burns

Laminar analysis of visually evoked activity in the primary visual cortex

Studying the laminar pattern of neural activity is crucial for understanding the processing of neural signals in the cerebral cortex. We measured neural population activity [multiunit spike activity (MUA) and local field potential, LFP] in Macaque primary visual cortex (V1) in response to drifting grating stimuli. Sustained visually driven MUA was at an approximately constant level across cortical depth in V1. However, sustained, visually driven, local field potential power, which was concentrated in the γ-band (20–60 Hz), was greatest at the cortical depth corresponding to cortico-cortical output layers 2, 3, and 4B. γ-band power also tends to be more sustained in the output layers. Overall, cortico-cortical output layers accounted for 67% of total γ-band activity in V1, whereas 56% of total spikes evoked by drifting gratings were from layers 2, 3, and 4B. The high-resolution layer specificity of γ-band power, the laminar distribution of MUA and γ-band activity, and their dynamics imply that neural activity in V1 is generated by laminar-specific mechanisms. In particular, visual responses of MUA and γ-band activity in cortico-cortical output layers 2, 3, and 4B seem to be strongly influenced by laminar-specific recurrent circuitry and/or feedback.

Network dynamics of the brain and influence of the epileptic seizure onset zone

Significance In epilepsy, seizures elicit changes in the functional connectivity of the brain that shed insight into the seizures’ nature and onset zone. We investigated the brain connectivity of patients with partial epileptic seizures from continuous multiday recordings and found that (i) the connectivity defines a finite set of brain states, (ii) seizures are characterized by a consistent progression of states, and (iii) the seizure onset zone is isolated from the surrounding regions at seizure onset but becomes most connected toward seizure termination. Our results suggest that a finite-dimensional state space model may characterize the dynamics of the epileptic brain and ultimately help localize the seizure onset zone, which is currently done by clinicians through visual inspection of electrocorticographic recordings. The human brain is a dynamic networked system. Patients with partial epileptic seizures have focal regions that periodically diverge from normal brain network dynamics during seizures. We studied the evolution of brain connectivity before, during, and after seizures with graph-theoretic techniques on continuous electrocorticographic (ECoG) recordings (5.4 ± 1.7 d per patient, mean ± SD) from 12 patients with temporal, occipital, or frontal lobe partial onset seizures. Each electrode was considered a node in a graph, and edges between pairs of nodes were weighted by their coherence within a frequency band. The leading eigenvector of the connectivity matrix, which captures network structure, was tracked over time and clustered to uncover a finite set of brain network states. Across patients, we found that (i) the network connectivity is structured and defines a finite set of brain states, (ii) seizures are characterized by a consistent sequence of states, (iii) a subset of nodes is isolated from the network at seizure onset and becomes more connected with the network toward seizure termination, and (iv) the isolated nodes may identify the seizure onset zone with high specificity and sensitivity. To localize a seizure, clinicians visually inspect seizures recorded from multiple intracranial electrode contacts, a time-consuming process that may not always result in definitive localization. We show that network metrics computed from all ECoG channels capture the dynamics of the seizure onset zone as it diverges from normal overall network structure. This suggests that a state space model can be used to help localize the seizure onset zone in ECoG recordings.

Is Gamma-Band Activity in the Local Field Potential of V1 Cortex a “Clock” or Filtered Noise?

Gamma-band (25–90 Hz) peaks in local field potential (LFP) power spectra are present throughout the cerebral cortex and have been related to perception, attention, memory, and disorders (e.g., schizophrenia and autism). It has been theorized that gamma oscillations provide a “clock” for precise temporal encoding and “binding” of signals about stimulus features across brain regions. For gamma to function as a clock, it must be autocoherent: phase and frequency conserved over a period of time. We computed phase and frequency trajectories of gamma-band bursts, using time-frequency analysis of LFPs recorded in macaque primary visual cortex (V1) during visual stimulation. The data were compared with simulations of random networks and clock signals in noise. Gamma-band bursts in LFP data were statistically indistinguishable from those found in filtered broadband noise. Therefore, V1 LFP data did not contain clock-like gamma-band signals. We consider possible functions for stochastic gamma-band activity, such as a synchronizing pulse signal.

Stochastic Generation of Gamma-Band Activity in Primary Visual Cortex of Awake and Anesthetized Monkeys

Oscillatory neural activity within the gamma band (25–90 Hz) is generally thought to be able to provide a timing signal for harmonizing neural computations across different brain regions. Using time–frequency analyses of the dynamics of gamma-band activity in the local field potentials recorded from monkey primary visual cortex, we found identical temporal characteristics of gamma activity in both awake and anesthetized brain states, including large variability of peak frequency, brief oscillatory epochs (<100 ms on average), and stochastic statistics of the incidence and duration of oscillatory events. These findings indicate that gamma-band activity is temporally unstructured and is inherently a stochastic signal generated by neural networks. This idea was corroborated further by our neural-network simulations. Our results suggest that gamma-band activity is too random to serve as a clock signal for synchronizing neuronal responses in awake as in anesthetized monkeys. Instead, gamma-band activity is more likely to be filtered neuronal network noise. Its mean frequency changes with global state and is reduced under anesthesia.

Comparisons of the Dynamics of Local Field Potential and Multiunit Activity Signals in Macaque Visual Cortex

The local field potential (LFP) and multiunit activity (MUA) are extracellularly recorded signals that describe local neuronal network dynamics. In our experiments, the LFP and MUA, recorded from the same electrode in macaque primary visual cortex V1 in response to drifting grating visual stimuli, were evaluated on coarse timescales (∼1–5 s) and fine timescales (<0.1 s). On coarse timescales, MUA and the LFP both produced sustained visual responses to optimal and non-optimal oriented visual stimuli. The sustainedness of the two signals across the population of recording sites was correlated (correlation coefficient, ∼0.4). At most recording sites, the MUA was at least as sustained as the LFP and significantly more sustained for optimal orientations. In previous literature, the blood oxygen level-dependent (BOLD) signal of functional magnetic resonance imaging studies was found to be more strongly correlated with the LFP than with the MUA as a result of the lack of sustained response in the MUA signal. Because we found that MUA was as sustained as the LFP, MUA may also be correlated with BOLD. On fine timescales, we computed the coherence between the LFP and MUA over the frequency range 10–150 Hz. The LFP and MUA were weakly but significantly coherent (∼0.14) in the gamma band (20–90 Hz). The amount of gamma-band coherence was correlated with the power in the gamma band of the LFP. The data were consistent with the proposal that the LFP and MUA are generated in a noisy, resonant cortical network.

Quickest detection of drug-resistant seizures: an optimal control approach.

Epilepsy affects 50 million people worldwide, and seizures in 30% of the cases remain drug resistant. This has increased interest in responsive neurostimulation, which is most effective when administered during seizure onset. We propose a novel framework for seizure onset detection that involves (i) constructing statistics from multichannel intracranial EEG (iEEG) to distinguish nonictal versus ictal states; (ii) modeling the dynamics of these statistics in each state and the state transitions; you can remove this word if there is no room. (iii) developing an optimal control-based “quickest detection” (QD) strategy to estimate the transition times from nonictal to ictal states from sequential iEEG measurements. The QD strategy minimizes a cost function of detection delay and false positive probability. The solution is a threshold that non-monotonically decreases over time and avoids responding to rare events that normally trigger false positives. We applied QD to four drug resistant epileptic patients (168 hour continuous recordings, 26–44 electrodes, 33 seizures) and achieved 100% sensitivity with low false positive rates (0.16 false positive/hour). This article is part of a Supplemental Special Issue entitled The Future of Automated Seizure Detection and Prediction.

Searching for Autocoherence in the Cortical Network with a Time-Frequency Analysis of the Local Field Potential

Gamma-band peaks in the power spectrum of local field potentials (LFP) are found in multiple brain regions. It has been theorized that gamma oscillations may serve as a ‘clock’ signal for the purposes of precise temporal encoding of information and ‘binding’ of stimulus features across regions of the brain. Neurons in model networks may exhibit periodic spike firing or synchronized membrane potentials that give rise to a gamma-band oscillation that could operate as a ‘clock’. The phase of the oscillation in such models is conserved over the length of the stimulus. We define these types of oscillations to be ‘autocoherent’. We investigated the hypothesis that autocoherent oscillations are the basis of the experimentally observed gamma-band peaks: the autocoherent oscillator (ACO) hypothesis. To test the ACO hypothesis, we developed a new technique to analyze the autocoherence of a time-varying signal. This analysis used the continuous Gabor transform to examine the time evolution of the phase of each frequency component in the power spectrum. Using this analysis method, we formulated a statistical test to compare the ACO hypothesis with measurements of the LFP in macaque primary visual cortex, V1. The experimental data were not consistent with the ACO hypothesis. Gamma-band activity recorded in V1 did not have the properties of a ‘clock’ signal during visual stimulation. We propose instead that the source of the gamma-band spectral peak is the resonant V1 network driven by random inputs.

Brain state evolution during seizure and under anesthesia: A network-based analysis of stereotaxic eeg activity in drug-resistant epilepsy patients

Epilepsy is a neurological condition with a prevalence of 1%, and 14-34% have medically refractory epilepsy (MRE). Seizures in focal MRE are generated by a single epileptogenic zone (or focus), thus there is potentially a curative procedure - surgical resection. This procedure depends significantly on correct identification of the focus, which is often uncertain in clinical practice. In this study, we analyzed intracranial stereotaxic EEG (sEEG) data recorded in two human patients with drug-resistant epilepsy prior to undergoing resection surgery. We view the sEEG data as samples from the brain network and hypothesize that seizure foci can be identified based on their network connectivity during seizure. Specifically, we computed a time sequence of connectivity matrices from EEG recordings that represent network structure over time. For each patient, connectivity between electrodes was measured using the coherence in a given frequency band. Matrix structure was analyzed using singular value decomposition and the leading singular vector was used to estimate each electrodes time dependent centrality (importance to the networks connectivity). Our preliminary study suggests that seizure foci may be the most weakly connected regions in the brain during the beginning of a seizure and the most strongly connected regions towards the end of a seizure. Additionally, in one of the patients analyzed, the network connectivity under anesthesia highlights seizure foci. Ultimately, network centrality computed from sEEG activity may be used to develop an automated, reliable, and computationally efficient algorithm for identifying seizure foci.

Multivariate analysis of SEEG signals during seizure

Epilepsy is a neurological disorder that affects tens of millions of people every year and is characterized by sudden-onset seizures which are often associated with physical convulsions. Effective treatment and management of epilepsy would be greatly improved if convulsions could be caught quickly through early seizure detection. However, this is still a largely open problem due to the challenge of finding a robust statistic from the neural measurements. This paper suggests a new multivariate statistic by combining spectral techniques with matrix theory. Specifically, stereoelectroencephalography (SEEG) data was used to generate a series of coherence connectivity matrices which were then examined using singular value decomposition. Tracking the relative angles of the first singular vectors generated from this data provides an effective way of defining the most dominant characteristics of the SEEG during the normal, the pre-ictal, and the ictal states. This paper indicates that the first singular vector has a characteristic direction indicative of the seizure state and illustrates a data analysis method that incorporates all neural data as opposed to a small selection of channels.

State Dynamics of the Epileptic Brain

Communication between specialized regions of the brain is a dynamic process allowing for different connections to accomplish different tasks. While the content of interregional communication is complex, the pattern of connectivity (i.e., which regions communicate) may lie in a lower dimensional state-space. In epilepsy, seizures elicit changes in connectivity, whose patterns shed insight into the nature of seizures and the seizure focus. We investigated connectivity in 3 patients by applying network-based analysis on multi-day subdural electrocorticographic recordings (ECoG). We found that (i) the network connectivity defines a finite set of brain states, (ii) seizures are characterized by a consistent progression of states, and (iii) the focus is isolated from surrounding regions at the seizure onset and becomes most connected in the network towards seizure termination. Our results suggest that a finite-dimensional state-space model may characterize the dynamics of the epileptic brain, and may ultimately be used to localize seizure foci.© 2013 ASME