Experimental observation of a theoretically predicted nonlinear sleep spindle harmonic in human EEG

Relationships between lognormal distributions of neural properties, activity, criticality, and connectivity

Relationships between convergence of inputs onto neurons, divergence of outputs from them, synaptic strengths, nonlinear firing response properties, and randomness of axonal ranges are systematically explored by interrelating means and variances of synaptic strengths, firing rates, and soma voltages. When self-consistency is imposed, it is found that broad distributions of synaptic strength are a necessary concomitant of the known massive convergence of inputs to individual neurons, and observed widths of lognormal distributions of synaptic strength and firing rate are explained provided the brain is in a near-critical state, consistent with independent observations. The strongest individual synapses are shown to have an effect on soma voltage comparable to the effect of all others combined, which supports suggestions that they may have a key role in neural communication. Remarkably, inclusion of moderate randomness in characteristic axonal ranges is shown to account for the observed [Formula: see text]-fold variability in two-point connectivity at a given separation and [Formula: see text]-fold overall when the known mean exponential fall-off is included, consistent with observed near-lognormal distributions. Inferred axonal deviations from straight-line paths are also consistent with independent estimates.

Coupling of gamma band activity to sleep spindle oscillations - a combined EEG/MEG study

Sleep spindles are crucial to memory consolidation. Cortical gamma oscillations (30-100 Hz) are considered to reflect processing of memory in local cortical networks. The temporal and regulatory relationship between spindles and gamma activity might therefore provide clues into how sleep strengthens cortical memory representations. Here, combining EEG with MEG recordings during sleep in healthy humans (n = 12), we investigated the temporal relationships of cortical gamma band activity, always measured by MEG, during fast (12-16 Hz) and slow (8-12 Hz) sleep spindles detected in the EEG or MEG. Time-frequency distributions did not show a consistent coupling of gamma to the spindle oscillation, although activity in the low gamma (30-40 Hz) and neighboring beta range (<30 Hz) was generally increased during spindles. However, more fine-grained analyses of cross-frequency interactions revealed that both low and high gamma power (30-100 Hz) was coupled to the phase of slow and fast EEG spindles, importantly, with this coupling at a fixed phase only for the oscillations within an individual spindle, but with variable phase across spindles. We did not observe any coupling of gamma activity for spindles detected solely in the MEG and not in parallel EEG recordings, raising the possibility that these are more local spindles of different quality. Similar to fast spindle activity, low gamma band power followed a ~0.025 Hz infraslow rhythm during sleep whose frequency, however, was significantly faster than that of spindle activity. Our findings suggest a general function of fast and slow spindles that by spanning larger cortical networks might serve to synchronize gamma band activity occurring in more local but distributed networks. Thereby, spindles might help linking local memory processing between distributed networks.

Neural field theory of evoked response potentials in a spherical brain geometry

Evoked response potentials (ERPs) are calculated in spherical and planar geometries using neural field theory of the corticothalamic system. The ERP is modeled as an impulse response and the resulting modal effects of spherical corticothalamic dynamics are explored, showing that results for spherical and planar geometries converge in the limit of large brain size. Cortical modal effects can lead to a double-peak structure in the ERP time series. It is found that the main difference between infinite planar geometry and spherical geometry is that the ERP peak is sharper and stronger in the spherical geometry. It is also found that the magnitude of the response decreases with increasing spatial width of the stimulus at the cortex. The peak is slightly delayed at large angles from the stimulus point, corresponding to group velocities of 6-10 m s^{-1}. Strong modal effects are found in the spherical geometry, with the lowest few modes sufficing to describe the main features of ERPs, except very near to spatially narrow stimuli.

Sleep Modelling across Physiological Levels

Sleep and circadian rhythms are regulated across multiple functional, spatial and temporal levels: from genes to networks of coupled neurons and glial cells, to large scale brain dynamics and behaviour. The dynamics at each of these levels are complex and the interaction between the levels is even more so, so research have mostly focused on interactions within the levels to understand the underlying mechanisms—the so-called reductionist approach. Mathematical models were developed to test theories of sleep regulation and guide new experiments at each of these levels and have become an integral part of the field. The advantage of modelling, however, is that it allows us to simulate and test the dynamics of complex biological systems and thus provides a tool to investigate the connections between the different levels and study the system as a whole. In this paper I review key models of sleep developed at different physiological levels and discuss the potential for an integrated systems biology approach for sleep regulation across these levels. I also highlight the necessity of building mechanistic connections between models of sleep and circadian rhythms across these levels.

Complexity of Wake Electroencephalography Correlates With Slow Wave Activity After Sleep Onset

Sleep electroencephalography (EEG) provides an opportunity to study sleep scientifically, whose chaotic, dynamic, complex, and dissipative nature implies that non-linear approaches could uncover some mechanism of sleep. Based on well-established complexity theories, one hypothesis in sleep medicine is that lower complexity of brain waves at pre-sleep state can facilitate sleep initiation and further improve sleep quality. However, this has never been studied with solid data. In this study, EEG collected from healthy subjects was used to investigate the association between pre-sleep EEG complexity and sleep quality. Multiscale entropy analysis (MSE) was applied to pre-sleep EEG signals recorded immediately after light-off (while subjects were awake) for measuring the complexities of brain dynamics by a proposed index, CI₁₋₃₀. Slow wave activity (SWA) in sleep, which is commonly used as an indicator of sleep depth or sleep intensity, was quantified based on two methods, traditional Fast Fourier transform (FFT) and ensemble empirical mode decomposition (EEMD). The associations between wake EEG complexity, sleep latency, and SWA in sleep were evaluated. Our results demonstrated that lower complexity before sleep onset is associated with decreased sleep latency, indicating a potential facilitating role of reduced pre-sleep complexity in the wake-sleep transition. In addition, the proposed EEMD-based method revealed an association between wake complexity and quantified SWA in the beginning of sleep (90 min after sleep onset). Complexity metric could thus be considered as a potential indicator for sleep interventions, and further studies are encouraged to examine the application of EEG complexity before sleep onset in populations with difficulty in sleep initiation. Further studies may also examine the mechanisms of the causal relationships between pre-sleep brain complexity and SWA, or conduct comparisons between normal and pathological conditions.

NFTsim: Theory and Simulation of Multiscale Neural Field Dynamics

A user ready, portable, documented software package, NFTsim, is presented to facilitate numerical simulations of a wide range of brain systems using continuum neural field modeling. NFTsim enables users to simulate key aspects of brain activity at multiple scales. At the microscopic scale, it incorporates characteristics of local interactions between cells, neurotransmitter effects, synaptodendritic delays and feedbacks. At the mesoscopic scale, it incorporates information about medium to large scale axonal ranges of fibers, which are essential to model dissipative wave transmission and to produce synchronous oscillations and associated cross-correlation patterns as observed in local field potential recordings of active tissue. At the scale of the whole brain, NFTsim allows for the inclusion of long range pathways, such as thalamocortical projections, when generating macroscopic activity fields. The multiscale nature of the neural activity produced by NFTsim has the potential to enable the modeling of resulting quantities measurable via various neuroimaging techniques. In this work, we give a comprehensive description of the design and implementation of the software. Due to its modularity and flexibility, NFTsim enables the systematic study of an unlimited number of neural systems with multiple neural populations under a unified framework and allows for direct comparison with analytic and experimental predictions. The code is written in C++ and bundled with Matlab routines for a rapid quantitative analysis and visualization of the outputs. The output of NFTsim is stored in plain text file enabling users to select from a broad range of tools for offline analysis. This software enables a wide and convenient use of powerful physiologically-based neural field approaches to brain modeling. NFTsim is distributed under the Apache 2.0 license.

Dependence of absence seizure dynamics on physiological parameter evolution

A neural field model of the corticothalamic system is applied to investigate the temporal and spectral characteristics of absence seizures in the presence of a temporally varying connection strength between the cerebral cortex and thalamus. Increasing connection strength drives the system into an absence seizure-like state once a threshold is passed and a supercritical Hopf bifurcation occurs. The dynamics and spectral characteristics of the resulting model seizures are explored as functions of maximum connection strength, time above threshold, and the rate at which the connection strength increases (ramp rate). Our results enable spectral and temporal characteristics of seizures to be related to changes in the underlying physiological evolution of connections via nonlinear dynamics and neural field theory. Spectral analysis reveals that the power of the harmonics and the duration of the oscillations increase as the maximum connection strength and the time above threshold increase. It is also found that the time to reach the stable limit-cycle seizure oscillation from the instability threshold decreases with the square root of the ramp rate.

Neural field theory of perceptual echo and implications for estimating brain connectivity

Neural field theory is used to predict and analyze the phenomenon of perceptual echo in which random input stimuli at one location are correlated with electroencephalographic responses at other locations. It is shown that this echo correlation (EC) yields an estimate of the transfer function from the stimulated point to other locations. Modal analysis then explains the observed spatiotemporal structure of visually driven EC and the dominance of the alpha frequency; two eigenmodes of similar amplitude dominate the response, leading to temporal beating and a line of low correlation that runs from the crown of the head toward the ears. These effects result from mode splitting and symmetry breaking caused by interhemispheric coupling and cortical folding. It is shown how eigenmodes obtained from functional magnetic resonance imaging experiments can be combined with temporal dynamics from EC or other evoked responses to estimate the spatiotemporal transfer function between any two points and hence their effective connectivity.

Theory of corticothalamic brain activity in a spherical geometry: Spectra, coherence, and correlation

Corticothalamic neural field theory is applied to a spherical geometry to better model neural activity in the human brain and is also compared with planar approximations. The frequency power spectrum, correlation, and coherence functions are computed analytically and numerically. The effects of cortical boundary conditions and resulting modal aspects of spherical corticothalamic dynamics are explored, showing that the results of spherical and finite planar geometries converge to those for the infinite planar geometry in the limit of large brain size. Estimates are made of the point at which modal series can be truncated and it is found that for physiologically plausible parameters only the lowest few spatial eigenmodes are needed for an accurate representation of macroscopic brain activity. A difference between the geometries is that there is a low-frequency 1/f spectrum in the infinite planar geometry, whereas in the spherical geometry it is 1/f^{2}. Another difference is that the alpha peak in the spherical geometry is sharper and stronger than in the planar geometry. Cortical modal effects can lead to a double alpha peak structure in the power spectrum, although the main determinant of the alpha peak is corticothalamic feedback. In the spherical geometry, the cross spectrum between two points is found to only depend on their relative distance apart. At small spatial separations the low-frequency cross spectrum is stronger than for an infinite planar geometry and the alpha peak is sharper and stronger due to the partitioning of the energy into discrete modes. In the spherical geometry, the coherence function between points decays monotonically as their separation increases at a fixed frequency, but persists further at resonant frequencies. The correlation between two points is found to be positive, regardless of the time lag and spatial separation, but decays monotonically as the separation increases at fixed time lag. At fixed distance the correlation has peaks at multiples of the period of the dominant frequency of system activity.

Spectral signatures of activity-dependent neural feedback in the corticothalamic system

The modulation of neural quantities by presynaptic and postsynaptic activities via local feedback processes is investigated by incorporating nonlinear phenomena such as relative refractory period, synaptic enhancement, synaptic depression, and habituation. This is done by introducing susceptibilities, which quantify the response in either firing threshold or synaptic strength to unit change in either presynaptic or postsynaptic activity. Effects on the power spectra are then analyzed for a realistic corticothalamic model to determine the spectral signatures of various nonlinear processes and to what extent these are distinct. Depending on the feedback processes, there can be enhancements or reductions in low-frequency and/or alpha power, splitting of the alpha resonance, and/or appearance of new resonances at high frequencies. These features in the power spectra allow processes to be fully distinguished where they are unique, or partly distinguished if they are common to only a subset of feedbacks, and can potentially be used to constrain the types, strengths, and dynamics of feedbacks present.

K-complexes, spindles, and ERPs as impulse responses: unification via neural field theory

To interrelate K-complexes, spindles, evoked response potentials (ERPs), and spontaneous electroencephalography (EEG) using neural field theory (NFT), physiology-based NFT of the corticothalamic system is used to model cortical excitatory and inhibitory populations and thalamic relay and reticular nuclei. The impulse response function of the model is used to predict the responses to impulses, which are compared with transient waveforms in sleep studies. Fits to empirical data then allow underlying brain physiology to be inferred and compared with other waves. Spontaneous K-complexes, spindles, and other transient waveforms can be reproduced using NFT by treating them as evoked responses to impulsive stimuli with brain parameters appropriate to spontaneous EEG in sleep stage 2. Using this approach, spontaneous K-complexes and sleep spindles can be analyzed using the same single theory as previously been used to account for waking ERPs and other EEG phenomena. As a result, NFT can explain a wide variety of transient waveforms that have only been phenomenologically classified to date. This enables noninvasive fitting to be used to infer underlying physiological parameters. This physiology-based model reproduces the time series of different transient EEG waveforms; it has previously reproduced experimental EEG spectra, and waking ERPs, and many other observations, thereby unifying transient sleep waveforms with these phenomena.

Similar Spectral Power Densities Within the Schumann Resonance and a Large Population of Quantitative Electroencephalographic Profiles: Supportive Evidence for Koenig and Pobachenko

In 1954 and 1960 Koenig and his colleagues described the remarkable similarities of spectral power density profiles and patterns between the earth-ionosphere resonance and human brain activity which also share magnitudes for both electric field (mV/m) and magnetic field (pT) components. In 2006 Pobachenko and colleagues reported real time coherence between variations in the Schumann and brain activity spectra within the 6-16 Hz band for a small sample. We examined the ratios of the average potential differences (~3 μV) obtained by whole brain quantitative electroencephalography (QEEG) between rostral-caudal and left-right (hemispheric) comparisons of 238 measurements from 184 individuals over a 3.5 year period. Spectral densities for the rostral-caudal axis revealed a powerful peak at 10.25 Hz while the left-right peak was 1.95 Hz with beat-differences of ~7.5 to 8 Hz. When global cerebral measures were employed, the first (7-8 Hz), second (13-14 Hz) and third (19-20 Hz) harmonics of the Schumann resonances were discernable in averaged QEEG profiles in some but not all participants. The intensity of the endogenous Schumann resonance was related to the 'best-of-fitness' of the traditional 4-class microstate model. Additional measurements demonstrated real-time coherence for durations approximating microstates in spectral power density variations between Schumann frequencies measured in Sudbury, Canada and Cumiana, Italy with the QEEGs of local subjects. Our results confirm the measurements reported by earlier researchers that demonstrated unexpected similarities in the spectral patterns and strengths of electromagnetic fields generated by the human brain and the earth-ionospheric cavity.

Real-time automated EEG tracking of brain states using neural field theory

A real-time fitting system is developed and used to fit the predictions of an established physiologically-based neural field model to electroencephalographic spectra, yielding a trajectory in a physiological parameter space that parametrizes intracortical, intrathalamic, and corticothalamic feedbacks as the arousal state evolves continuously over time. This avoids traditional sleep/wake staging (e.g., using Rechtschaffen-Kales stages), which is fundamentally limited because it forces classification of continuous dynamics into a few discrete categories that are neither physiologically informative nor individualized. The classification is also subject to substantial interobserver disagreement because traditional staging relies in part on subjective evaluations. The fitting routine objectively and robustly tracks arousal parameters over the course of a full night of sleep, and runs in real-time on a desktop computer. The system developed here supersedes discrete staging systems by representing arousal states in terms of physiology, and provides an objective measure of arousal state which solves the problem of interobserver disagreement. Discrete stages from traditional schemes can be expressed in terms of model parameters for backward compatibility with prior studies.

Physiologically based arousal state estimation and dynamics

A neural field model of the brain is used to represent brain states using physiologically based parameters rather than arbitrary, discrete sleep stages. Each brain state is represented as a point in a physiologically parametrized space. Over time, changes in brain state cause these points to trace continuous trajectories, unlike the artificial discrete jumps in sleep stage that occur with traditional sleep staging. The discrete Rechtschaffen and Kales sleep stages are associated with regions in the physiological parameter space based on their electroencephalographic features, which enables interpretation of traditional sleep stages in terms of physiological trajectories. Wake states are found to be associated with strong positive corticothalamic feedback compared to sleep. The existence of physiologically valid trajectories between brain states in the model is demonstrated. Actual trajectories for an individual can be determined by fitting the model using EEG alone, and enable analysis of the physiological differences between subjects.

Neural field theory of nonlinear wave-wave and wave-neuron processes

Systematic expansion of neural field theory equations in terms of nonlinear response functions is carried out to enable a wide variety of nonlinear wave-wave and wave-neuron processes to be treated systematically in systems involving multiple neural populations. The results are illustrated by analyzing second-harmonic generation, and they can also be applied to wave-wave coalescence, multiharmonic generation, facilitation, depression, refractoriness, and other nonlinear processes.