EEG simulation: variation of spectral envelope, pulse synchrony and ≈40 Hz oscillation

 Macroscopic EEG travelling wave phenomena and cortical pulse synchronisation effects are related within a single simple simulation. Non-specific activation acts to control the transfer function of the simulated cortex, and thus determines the relative amplitude of macroscopic EEG waves generated by rhythmic inputs. When concurrent asynchronous excitatory inputs to separate, local, cortical sites are introduced, the simulation reproduces both gamma-band (40 Hz) electrocorticogram (ECoG) activity and synchronous oscillation of action potential pulse density at the separate sites. The gamma-band ECoG and pulse synchrony effects depend on different mechanisms: the former upon local excitatory/inhibitory interactions, and the latter on cortico-cortical interactions. The pattern of synchronous activity depends upon both structural and dynamic aspects of gain, and is sustained by linearised versions of the simulation’s state equations. Dynamic properties of the simulation, which are independent of scale, describe both microscopic and macroscopic phenomena, all in accord with physiological findings. Received: 25 June 1996 / Accepted in revised form: 29 November 1996     

Frequency component selection for an ECoG-based brain-computer interface.

The aim of the present study was to investigate the most significant frequency components in electrocorticogram (ECoG) recordings in order to operate a brain computer interface (BCI). For this purpose the time-frequency ERD/ERS map and the distinction sensitive learning vector quantization (DSLVQ) are applied to ECoG from three subjects, recorded during a self-paced finger movement. The results show that the ERD/ERS pattern found in ECoG generally matches the ERD/ERS pattern found in EEG recordings, but has an increased prevalence of frequency components in the beta range.     

Functional microsurgical partial callosotomy in patients with secondary generalized epilepsies. II. Mesial surface electrocorticography

Eight patients with secondary generalized epilepsy not alleviated by medical treatment underwent partial callosotomy. During the surgical procedure, they had mesial surface ECoG recordings taken from both frontal and parietal lobes, using large flat multilead platinum electrodes, and simultaneously recordings from a number of scalp positions, using needle electrodes. In all cases studied, this approach demonstrated one or, more commonly, several focal areas of epileptiform activity discharging independently over the mesial aspects of one or both hemispheres. The findings were correlated with the pre- and postoperative EEG patterns, in the light of current concepts of generalized epilepsies.     

Nonlinear EEG analysis based on a neural mass model

The well-known neural mass model described by Lopes da Silva et al. (1976) and Zetterberg et al. (1978) is fitted to actual EEG data. This is achieved by reformulating the original set of integral equations as a continuous-discrete state space model. The local linearization approach is then used to discretize the state equation and to construct a nonlinear Kalman filter. On this basis, a maximum likelihood procedure is used for estimating the model parameters for several EEG recordings. The analysis of the noise-free differential equations of the estimated models suggests that there are two different types of alpha rhythms: those with a point attractor and others with a limit cycle attractor. These attractors are also found by means of a nonlinear time series analysis of the EEG recordings. We conclude that the Hopf bifurcation described by Zetterberg et al. (1978) is present in actual brain dynamics. Received: 11 August 1997 / Accepted in revised form: 20 April 1999     

Timings and interactions of skilled musicians

We report here a preliminary study of interactive behavior between two members of a skilled string quartet performing a selected musical passage that required both performers to play several hundred notes in rapid succession at a steady tempo and in synchrony. Bowing movements were recorded using angular velocity sensors attached to their right forearms. The results show a high degree of temporal precision in both players. In addition, both players exhibited embedded rhythmic components in their timekeeping pattern, which arose from the grouping of notes in the musical score: four 16th notes to a beat. Within each group of four notes, we found a consistent timing microstructure: alternate upbows and alternate downbows had different mean durations. Both players’ bowings could be modeled as alternating renewal processes. In addition, we report evidence of interactive coupling between the players as an essential component of their joint performance. The alternating renewal model enables us to propose a note-generation process that has implications for the central generators underlying the observed behavior and their hierarchical organization. We discuss the implications of this model for the organization and execution of more complex motor sequences.     

Behavioral states may be associated with distinct spatial patterns in electrocorticogram

To determine if behavioral states are associated with unique spatial electrocorticographic (ECoG) patterns, we obtained recordings with a microgrid electrode array applied to the cortical surface of a human subject. The array was constructed with the intent of extracting maximal spatial information by optimizing interelectrode distances. A 34-year-old patient with intractable epilepsy underwent intracranial ECoG monitoring after standard methods failed to reveal localization of seizures. During the 8-day period of invasive recording, in addition to standard clinical electrodes a square 1 × 1 cm microgrid array with 64 electrodes (1.25 mm separation) was placed on the right inferior temporal gyrus. Careful review of video recordings identified four extended naturalistic behaviors: reading, conversing on the telephone, looking at photographs, and face-to-face interactions. ECoG activity recorded with the microgrid that corresponded to these behaviors was collected and ECoG spatial patterns were analyzed. During periods of ECoG selected for analysis, no electrographic seizures or epileptiform patterns were present. Moments of maximal spatial variance are shown to cluster by behavior. Comparisons between conditions using a permutation test reveal significantly different spatial patterns for each behavior. We conclude that ECoG recordings obtained on the cortical surface with optimal high spatial frequency resolution reveal distinct local spatial patterns that reflect different behavioral states, and we predict that similar patterns will be found in many if not most cortical areas on which a microgrid is placed.     

Divergent Cortical Generators of MEG and EEG during Human Sleep Spindles Suggested by Distributed Source Modeling

Background

Sleep spindles are ∼1-second bursts of 10–15 Hz activity, occurring during normal stage 2 sleep. In animals, sleep spindles can be synchronous across multiple cortical and thalamic locations, suggesting a distributed stable phase-locked generating system. The high synchrony of spindles across scalp EEG sites suggests that this may also be true in humans. However, prior MEG studies suggest multiple and varying generators.

Methodology/Principal Findings

We recorded 306 channels of MEG simultaneously with 60 channels of EEG during naturally occurring spindles of stage 2 sleep in 7 healthy subjects. High-resolution structural MRI was obtained in each subject, to define the shells for a boundary element forward solution and to reconstruct the cortex providing the solution space for a noise-normalized minimum norm source estimation procedure. Integrated across the entire duration of all spindles, sources estimated from EEG and MEG are similar, diffuse and widespread, including all lobes from both hemispheres. However, the locations, phase and amplitude of sources simultaneously estimated from MEG versus EEG are highly distinct during the same spindles. Specifically, the sources estimated from EEG are highly synchronous across the cortex, whereas those from MEG rapidly shift in phase, hemisphere, and the location within the hemisphere.

Conclusions/Significance

The heterogeneity of MEG sources implies that multiple generators are active during human sleep spindles. If the source modeling is correct, then EEG spindles are generated by a different, diffusely synchronous system. Animal studies have identified two thalamo-cortical systems, core and matrix, that produce focal or diffuse activation and thus could underlie MEG and EEG spindles, respectively. Alternatively, EEG spindles could reflect overlap at the sensors of the same sources as are seen from the MEG. Although our results generally match human intracranial recordings, additional improvements are possible and simultaneous intra- and extra-cranial measures are needed to test their accuracy.     

Increase trend of correlation and phase synchrony of microwire iEEG before macroseizure onset

Micro/macrowire intracranial EEG (iEEG) signals recorded from implanted micro/macroelectrodes in epileptic patients have received great attention and are considered to include much information of neuron activities in seizure transition compared to scalp EEG from cortical electrodes. Microelectrode is contacted more close to neurons than macroelectrode and it is more sensitive to neuron activity changes than macroelectrode. Microwire iEEG recordings are inevitably advantageous over macrowire iEEG recordings to reveal neuronal mechanisms contributing to the generation of seizures. In this study, we investigate the seizure generation from microwire iEEG recordings and discuss synchronization of microwire iEEGs in four frequency bands: alpha (1−30 Hz), gamma (30−80 Hz), ripple (80–250 Hz), and fast ripple (>250 Hz) via two measures: correlation and phase synchrony. We find that an increase trend of correlation or phase synchrony exists before the macroseizure onset mostly in gamma and ripple bands where the duration of the preictal states varied in different seizures ranging up to a few seconds (minutes). This finding is contrast to the well-known result that a decrease of synchronization in macro domains exists before the macroseizure onset. The finding demonstrates that it is only when the seizure has recruited enough surrounding brain tissue does the signal become strong enough to be observed on the clinical macroelectrode and as a result support the hypothesis of progressive coalescence of microseizure domains. The potential ramifications of such an early detection of microscale seizure activity may open a new window on treatment by making possible disruption of seizure activity before it becomes fully established.     

Decoding Individual Finger Movements from One Hand Using Human EEG Signals

Brain computer interface (BCI) is an assistive technology, which decodes neurophysiological signals generated by the human brain and translates them into control signals to control external devices, e.g., wheelchairs. One problem challenging noninvasive BCI technologies is the limited control dimensions from decoding movements of, mainly, large body parts, e.g., upper and lower limbs. It has been reported that complicated dexterous functions, i.e., finger movements, can be decoded in electrocorticography (ECoG) signals, while it remains unclear whether noninvasive electroencephalography (EEG) signals also have sufficient information to decode the same type of movements. Phenomena of broadband power increase and low-frequency-band power decrease were observed in EEG in the present study, when EEG power spectra were decomposed by a principal component analysis (PCA). These movement-related spectral structures and their changes caused by finger movements in EEG are consistent with observations in previous ECoG study, as well as the results from ECoG data in the present study. The average decoding accuracy of 77.11% over all subjects was obtained in classifying each pair of fingers from one hand using movement-related spectral changes as features to be decoded using a support vector machine (SVM) classifier. The average decoding accuracy in three epilepsy patients using ECoG data was 91.28% with the similarly obtained features and same classifier. Both decoding accuracies of EEG and ECoG are significantly higher than the empirical guessing level (51.26%) in all subjects (p<0.05). The present study suggests the similar movement-related spectral changes in EEG as in ECoG, and demonstrates the feasibility of discriminating finger movements from one hand using EEG. These findings are promising to facilitate the development of BCIs with rich control signals using noninvasive technologies.     

Role of the Brainstem Reticular Formation in the Mechanisms of Cortical Electrogenesis

In rats, we studied peculiarities of the electrocorticogram (ECoG) and its segments characterized by synchronization and desynchronization phenomena; these periods were differentiated using a segmentation procedure. The electrocorticogram was recorded under conditions of the intact brain (IB) and in animals with the isolated forebrain (IFB) using intercollicular transection. When ECoG was recorded from IFB preparations, it differed from ECoG of the IB in decreased amplitudes of the rhythms and their reorganization with shifts of the frequency characteristics toward lower values, as well as in greater amplitudes of ECoG rhythms in the left hemisphere (lateralization of the ECoG activity toward this hemisphere). Significant differences of the indices of functional asymmetry were observed in the two examined experimental situations. Using calculations of multiple linear regression and subsequent visualization of the results using polycyclic multigraphs, we found a clear specificity in the formation of correlation links between the amplitudes of different ECoG rhythms under conditions of the IB and IFB preparation. We discuss the role of the brainstem structures, in particular that of the reticular formation, in the organization of integral ECoG activity and the possible importance of mutual influences between hypothetical generators of the ECoG rhythmus for this organization.Neirofiziologiya/Neurophysiology, Vol. 37, No. 1, pp. 39–51, January–February, 2005.     

A synchronization-desynchronization code for natural communication signals

Synchronous spiking of neural populations is hypothesized to play important computational roles in forming neural assemblies and solving the binding problem. Although the opposite phenomenon of desynchronization is well known from EEG studies, it is largely neglected on the neuronal level. We here provide an example of in vivo recordings from weakly electric fish demonstrating that, depending on the social context, different types of natural communication signals elicit transient desynchronization as well as synchronization of the electroreceptor population without changing the mean firing rate. We conclude that, in general, both positive and negative changes in the degree of synchrony can be the relevant signals for neural information processing.     

Cortical Source Analysis of High-Density EEG Recordings in Children

EEG is traditionally described as a neuroimaging technique with high temporal and low spatial resolution. Recent advances in biophysical modelling and signal processing make it possible to exploit information from other imaging modalities like structural MRI that provide high spatial resolution to overcome this constraint¹. This is especially useful for investigations that require high resolution in the temporal as well as spatial domain. In addition, due to the easy application and low cost of EEG recordings, EEG is often the method of choice when working with populations, such as young children, that do not tolerate functional MRI scans well. However, in order to investigate which neural substrates are involved, anatomical information from structural MRI is still needed. Most EEG analysis packages work with standard head models that are based on adult anatomy. The accuracy of these models when used for children is limited², because the composition and spatial configuration of head tissues changes dramatically over development³. In the present paper, we provide an overview of our recent work in utilizing head models based on individual structural MRI scans or age specific head models to reconstruct the cortical generators of high density EEG. This article describes how EEG recordings are acquired, processed, and analyzed with pediatric populations at the London Baby Lab, including laboratory setup, task design, EEG preprocessing, MRI processing, and EEG channel level and source analysis.      

An empirical EEG analysis in brain death diagnosis for adults

Electroencephalogram (EEG) is often used in the confirmatory test for brain death diagnosis in clinical practice. Because EEG recording and monitoring is relatively safe for the patients in deep coma, it is believed to be valuable for either reducing the risk of brain death diagnosis (while comparing other tests such as the apnea) or preventing mistaken diagnosis. The objective of this paper is to study several statistical methods for quantitative EEG analysis in order to help bedside or ambulatory monitoring or diagnosis. We apply signal processing and quantitative statistical analysis for the EEG recordings of 32 adult patients. For EEG signal processing, independent component analysis (ICA) was applied to separate the independent source components, followed by Fourier and time-frequency analysis. For quantitative EEG analysis, we apply several statistical complexity measures to the EEG signals and evaluate the differences between two groups of patients: the subjects in deep coma, and the subjects who were categorized as brain death. We report statistically significant differences of quantitative statistics with real-life EEG recordings in such a clinical study, and we also present interpretation and discussions on the preliminary experimental results.

Phase Synchronization in Electroencephalographic Recordings Prognosticates Outcome in Paediatric Coma

Brain injury from trauma, cardiac arrest or stroke is the most important cause of death and acquired disability in the paediatric population. Due to the lifetime impact of brain injury, there is a need for methods to stratify patient risk and ultimately predict outcome. Early prognosis is fundamental to the implementation of interventions to improve recovery, but no clinical model as yet exists. Healthy physiology is associated with a relative high variability of physiologic signals in organ systems. This was first evaluated in heart rate variability research. Brain variability can be quantified through electroencephalographic (EEG) phase synchrony. We hypothesised that variability in brain signals from EEG recordings would correlate with patient outcome after brain injury. Lower variability in EEG phase synchronization, would be associated with poor patient prognosis. A retrospective study, spanning 10 years (2000–2010) analysed the scalp EEGs of children aged 1 month to 17 years in coma (Glasgow Coma Scale, GCS, <8) admitted to the paediatric critical care unit (PCCU) following brain injury from TBI, cardiac arrest or stroke. Phase synchrony of the EEGs was evaluated using the Hilbert transform and the variability of the phase synchrony calculated. Outcome was evaluated using the 6 point Paediatric Performance Category Score (PCPC) based on chart review at the time of hospital discharge. Outcome was dichotomized to good outcome (PCPC score 1 to 3) and poor outcome (PCPC score 4 to 6). Children who had a poor outcome following brain injury secondary to cardiac arrest, TBI or stroke, had a higher magnitude of synchrony (R index), a lower spatial complexity of the synchrony patterns and a lower temporal variability of the synchrony index values at 15 Hz when compared to those patients with a good outcome.     

Establishing a Statistical Link between Network Oscillations and Neural Synchrony

Pairs of active neurons frequently fire action potentials or “spikes” nearly synchronously (i.e., within 5 ms of each other). This spike synchrony may occur by chance, based solely on the neurons’ fluctuating firing patterns, or it may occur too frequently to be explicable by chance alone. When spike synchrony above chances levels is present, it may subserve computation for a specific cognitive process, or it could be an irrelevant byproduct of such computation. Either way, spike synchrony is a feature of neural data that should be explained. A point process regression framework has been developed previously for this purpose, using generalized linear models (GLMs). In this framework, the observed number of synchronous spikes is compared to the number predicted by chance under varying assumptions about the factors that affect each of the individual neuron’s firing-rate functions. An important possible source of spike synchrony is network-wide oscillations, which may provide an essential mechanism of network information flow. To establish the statistical link between spike synchrony and network-wide oscillations, we have integrated oscillatory field potentials into our point process regression framework. We first extended a previously-published model of spike-field association and showed that we could recover phase relationships between oscillatory field potentials and firing rates. We then used this new framework to demonstrate the statistical relationship between oscillatory field potentials and spike synchrony in: 1) simulated neurons, 2) in vitro recordings of hippocampal CA1 pyramidal cells, and 3) in vivo recordings of neocortical V4 neurons. Our results provide a rigorous method for establishing a statistical link between network oscillations and neural synchrony.     

Complex, multifocal, individual-specific attention-related cortical functional circuits

Recent studies focusing on the analysis of individual patterns of non-sensory-motor CNS activity may significantly alter our view of CNS functional mapping. We have recently provided evidence for highly variable attention-related Slow Potential (SP) generating cortical areas across individuals (Basile et al., 2003, 2006). In this work, we present new evidence, searching for other physiological indexes of attention by a new use of a well established method, for individual-specific sets of cortical areas active during expecting attention. We applied latency corrected peak averaging to oscillatory bursts, from 124-channel EEG recordings, and modeled their generators by current density reconstruction. We first computed event-related total power, and averaging was based on individual patterns of narrow task-induced band-power. This method is sensitive to activity out of synchrony with stimuli, and may detect task-related changes missed by regular Event-Related Potential (ERP) averaging. We additionally analyzed overall inter-electrode phase-coherence. The main results were (1) the detection of two bands of attention-induced beta range oscillations (around 25 and 21 Hz), whose scalp topography and current density cortical distribution were complex multi-focal, and highly variable across subjects, including prefrontal and posterior cortical areas. Most important, however, was the observation that (2) the generators of task-induced oscillations are largely the same individual-specific sets of cortical areas active during the resting, baseline state. We concluded that attention-related electrical cortical activity is highly individual-specific (significantly different from sensory-related visual evoked potentials or delta and theta induced band-power), and to a great extent already established during mere wakefulness. We discuss the critical implications of those results, in combination with other studies presenting individual data, to functional mapping: the need to abandon group averaging of task-related cortical activity and to revise studies on group averaged data, since the assumption of universal function to each cortical area appears deeply challenged. Clinical implications regard the interpretation of focal lesion consequences, functional reorganization, and neurosurgical planning.     

Nonlinear autoregressive analysis of the 3/s ictal electroencephalogram: implications for underlying dynamics

 In a previous study, nonlinear autoregressive (NLAR) models applied to ictal electroencephalogram (EEG) recordings in six patients revealed nonlinear signal interactions that correlated with seizure type and clinical diagnosis. Here we interpret these models from a theoretical viewpoint. Extended models with multiple nonlinear terms are employed to demonstrate the independence of nonlinear dynamical interactions identified in the ‘NLAR fingerprint’ of patients with 3/s seizure discharges. Analysis of the role of periodicity in the EEG signal reveals that the fingerprints reflect the dynamics not only of the periodic discharge itself, but also of the fluctuations of each cycle about an average waveform. A stability analysis is used to make qualitative inferences concerning the network properties of the ictal generators. Finally, the NLAR fingerprint is analyzed in the context of Volterra-Weiner theory. Received: 6 April 1994/Accepted in revised form: 18 November 1994     

Synchrony measures for biological neural networks

 Synchronous firing of a population of neurons has been observed in many experimental preparations; in addition, various mathematical neural network models have been shown, analytically or numerically, to contain stable synchronous solutions. In order to assess the level of synchrony of a particular network over some time interval, quantitative measures of synchrony are needed. We develop here various synchrony measures which utilize only the spike times of the neurons; these measures are applicable in both experimental situations and in computer models. Using a mathematical model of the CA3 region of the hippocampus, we evaluate these synchrony measures and compare them with pictorial representations of network activity. We illustrate how synchrony is lost and synchrony measures change as heterogeneity amongst cells increases. Theoretical expected values of the synchrony measures for different categories of network solutions are derived and compared with results of simulations. Received: 6 June 1994/Accepted in revised form: 13 January 1995     

An autoregressive method for the measurement of synchronization of interictal and ictal EEG signals

We propose a new measure of synchronization of multichannel ictal and interictal EEG signals. The measure is based on the residual covariance matrix of a multichannel autoregressive model. A major advantage of this measure is its ability to be interpreted both in the framework of stochastic and deterministic models. A preliminary analysis of EEG data from three patients using this measure documents the expected increased synchronization during ictal periods but also reveals that increased synchrony persists for prolonged periods (up to 2 h or more) in the postictal period. Received: 20 July 1997 / Accepted in revised form: 26 January 1999