Unified neurophysical model of EEG spectra and evoked potentials

Modulation of brain function and antidepressant effects by transcranial alternating current stimulation in patients with major depressive disorder: Evidence from ERP

Transcranial alternating current stimulation (tACS) is an emerging non-invasive neuromodulation treatment for major depressive disorder (MDD), but its mechanism remains unclear. Therefore, we evaluated the effects of tACS on event-related potentials (ERP) based on a randomized controlled study. All patients were divided into two groups to receive either 20 sessions 77.5Hz-tACS or 20 sessions of sham stimulation during 4 weeks. The Hamilton Depression Rating Scale for Depression -17 item (HAMD-17) and ERP during face-word Stroop task were recorded before and after the treatment (the fourth weekend). Our findings indicate a significant alleviation of depressive symptoms after tACS. For the behavioral performance, sham group showed a significant decrease in reaction time to the sad incongruent condition and an increase in accuracy to the happy condition. The active group showed an increase in accuracy to the incongruent condition. ERP analysis revealed that tACS significantly shortened the latency of P2 to incongruent condition, decreased the amplitude and prolonged the latency of N2 to negative condition. These ERP alterations suggest a potential rectification of negative bias and enhancement of cognitive functioning in patients with MDD, offering insights into the antidepressant mechanisms of tACS.

Recruiting neural field theory for data augmentation in a motor imagery brain-computer interface

We introduce a novel approach to training data augmentation in brain-computer interfaces (BCIs) using neural field theory (NFT) applied to EEG data from motor imagery tasks. BCIs often suffer from limited accuracy due to a limited amount of training data. To address this, we leveraged a corticothalamic NFT model to generate artificial EEG time series as supplemental training data. We employed the BCI competition IV '2a' dataset to evaluate this augmentation technique. For each individual, we fitted the model to common spatial patterns of each motor imagery class, jittered the fitted parameters, and generated time series for data augmentation. Our method led to significant accuracy improvements of over 2% in classifying the "total power" feature, but not in the case of the "Higuchi fractal dimension" feature. This suggests that the fit NFT model may more favorably represent one feature than the other. These findings pave the way for further exploration of NFT-based data augmentation, highlighting the benefits of biophysically accurate artificial data.

Neural field theory of adaptive effects on auditory evoked responses and mismatch negativity in multifrequency stimulus sequences

Physiologically based neural field theory (NFT) of the corticothalamic system, including adaptation, is used to calculate the responses evoked by trains of auditory stimuli that differ in frequency. In oddball paradigms, fully distinguishable frequencies lead to different standard (common stimulus) and deviant (rare stimulus) responses; the signal obtained by subtracting the standard response from the deviant is termed the mismatch negativity (MMN). In this analysis, deviant responses are found to correspond to unadapted cortex, whereas the part of auditory cortex that processes the standard stimuli adapts over several stimulus presentations until the final standard response form is achieved. No higher-order memory processes are invoked. In multifrequency experiments, the deviant response approaches the standard one as the deviant frequency approaches that of the standard and analytic criteria for this effect to be obtained. It is shown that these criteria can also be used to understand adaptation in random tone sequences. A method of probing MMNs and adaptation in random tone sequences is suggested to makes more use of such data.

Method for quantifying arousal and consciousness in healthy states and severe brain injury via EEG-based measures of corticothalamic physiology

BACKGROUND

Characterization of normal arousal states has been achieved by fitting predictions of corticothalamic neural field theory (NFT) to electroencephalographic (EEG) spectra to yield relevant physiological parameters.

NEW METHOD

A prior fitting method is extended to distinguish conscious and unconscious states in healthy and brain injured subjects by identifying additional parameters and clusters in parameter space.

RESULTS

Fits of NFT predictions to EEG spectra are used to estimate neurophysiological parameters in healthy and brain injured subjects. Spectra are used from healthy subjects in wake and sleep and from patients with unresponsive wakefulness syndrome, in a minimally conscious state (MCS), and emerged from MCS. Subjects cluster into three groups in parameter space: conscious healthy (wake and REM), sleep, and brain injured. These are distinguished by the difference X-Y between corticocortical (X) and corticothalamic (Y) feedbacks, and by mean neural response rates α and β to incoming spikes. X-Y tracks consciousness in healthy individuals, with smaller values in wake/REM than sleep, but cannot distinguish between brain injuries. Parameters α and β differentiate deep sleep from wake/REM and brain injury.

COMPARISON WITH EXISTING METHODS

Other methods typically rely on laborious clinical assessment, manual EEG scoring, or evaluation of measures like Φ from integrated information theory, for which no efficient method exists. In contrast, the present method can be automated on a personal computer.

CONCLUSION

The method provides a means to quantify consciousness and arousal in healthy and brain injured subjects, but does not distinguish subtypes of brain injury.

The mesoanatomy of the cortex, minimization of free energy, and generative cognition

Capacity for generativity and unlimited association is the defining characteristic of sentience, and this capacity somehow arises from neuronal self-organization in the cortex. We have previously argued that, consistent with the free energy principle, cortical development is driven by synaptic and cellular selection maximizing synchrony, with effects manifesting in a wide range of features of mesoscopic cortical anatomy. Here, we further argue that in the postnatal stage, as more structured inputs reach the cortex, the same principles of self-organization continue to operate at multitudes of local cortical sites. The unitary ultra-small world structures that emerged antenatally can represent sequences of spatiotemporal images. Local shifts of presynapses from excitatory to inhibitory cells result in the local coupling of spatial eigenmodes and the development of Markov blankets, minimizing prediction errors in each unit's interactions with surrounding neurons. In response to the superposition of inputs exchanged between cortical areas, more complicated, potentially cognitive structures are competitively selected by the merging of units and the elimination of redundant connections that result from the minimization of variational free energy and the elimination of redundant degrees of freedom. The trajectory along which free energy is minimized is shaped by interaction with sensorimotor, limbic, and brainstem mechanisms, providing a basis for creative and unlimited associative learning.

Preservation of thalamocortical circuitry is essential for good recovery after cardiac arrest

Continuous electroencephalographam (EEG) monitoring contributes to prediction of neurological outcome in comatose cardiac arrest survivors. While the phenomenology of EEG abnormalities in postanoxic encephalopathy is well known, the pathophysiology, especially the presumed role of selective synaptic failure, is less understood. To further this understanding, we estimate biophysical model parameters from the EEG power spectra from individual patients with a good or poor recovery from a postanoxic encephalopathy. This biophysical model includes intracortical, intrathalamic, and corticothalamic synaptic strengths, as well as synaptic time constants and axonal conduction delays. We used continuous EEG measurements from hundred comatose patients recorded during the first 48 h postcardiac arrest, 50 with a poor neurological outcome [cerebral performance category ( CPC = 5 ) ] and 50 with a good neurological outcome ( CPC = 1 ). We only included patients that developed (dis-)continuous EEG activity within 48 h postcardiac arrest. For patients with a good outcome, we observed an initial relative excitation in the corticothalamic loop and corticothalamic propagation that subsequently evolved towards values observed in healthy controls. For patients with a poor outcome, we observed an initial increase in the cortical excitation-inhibition ratio, increased relative inhibition in the corticothalamic loop, delayed corticothalamic propagation of neuronal activity, and severely prolonged synaptic time constants that did not return to physiological values. We conclude that the abnormal EEG evolution in patients with a poor neurological recovery after cardiac arrest may result from persistent and selective synaptic failure that includes corticothalamic circuitry and also delayed corticothalamic propagation.

Unification of free energy minimization, spatiotemporal energy, and dimension reduction models of V1 organization: Postnatal learning on an antenatal scaffold

Developmental selection of neurons and synapses so as to maximize pulse synchrony has recently been used to explain antenatal cortical development. Consequences of the same selection process—an application of the Free Energy Principle—are here followed into the postnatal phase in V1, and the implications for cognitive function are considered. Structured inputs transformed via lag relay in superficial patch connections lead to the generation of circumferential synaptic connectivity superimposed upon the antenatal, radial, “like-to-like” connectivity surrounding each singularity. The spatiotemporal energy and dimension reduction models of cortical feature preferences are accounted for and unified within the expanded model, and relationships of orientation preference (OP), space frequency preference (SFP), and temporal frequency preference (TFP) are resolved. The emergent anatomy provides a basis for “active inference” that includes interpolative modification of synapses so as to anticipate future inputs, as well as learn directly from present stimuli. Neurodynamic properties are those of heteroclinic networks with coupled spatial eigenmodes.

Neural Field Theory of Evoked Response Sequences and Mismatch Negativity With Adaptation

Physiologically based neural field theory of the corticothalamic system is used to calculate the responses evoked by trains of auditory stimuli that correspond to different cortical locations via the tonotopic map. The results are shown to account for standard and deviant evoked responses to frequent and rare stimuli, respectively, in the auditory oddball paradigms widely used in human cognitive studies, and the so-called mismatch negativity between them. It also reproduces a wide range of other effects and variants, including the mechanism by which a change in standard responses relative to deviants can develop through adaptation, different responses when two deviants are presented in a row or a standard is presented after two deviants, relaxation of standard responses back to deviant form after a stimulus-free period, and more complex sequences. Some cases are identified in which adaptation does not account for the whole difference between standard and deviant responses. The results thus provide a systematic means to determine how much of the response is due to adaptation in the system comprising the primary auditory cortex and medial geniculate nucleus, and how much requires involvement of higher-level processing.

Determination of Dynamic Brain Connectivity via Spectral Analysis

Spectral analysis based on neural field theory is used to analyze dynamic connectivity via methods based on the physical eigenmodes that are the building blocks of brain dynamics. These approaches integrate over space instead of averaging over time and thereby greatly reduce or remove the temporal averaging effects, windowing artifacts, and noise at fine spatial scales that have bedeviled the analysis of dynamical functional connectivity (FC). The dependences of FC on dynamics at various timescales, and on windowing, are clarified and the results are demonstrated on simple test cases, demonstrating how modes provide directly interpretable insights that can be related to brain structure and function. It is shown that FC is dynamic even when the brain structure and effective connectivity are fixed, and that the observed patterns of FC are dominated by relatively few eigenmodes. Common artifacts introduced by statistical analyses that do not incorporate the physical nature of the brain are discussed and it is shown that these are avoided by spectral analysis using eigenmodes. Unlike most published artificially discretized “resting state networks” and other statistically-derived patterns, eigenmodes overlap, with every mode extending across the whole brain and every region participating in every mode—just like the vibrations that give rise to notes of a musical instrument. Despite this, modes are independent and do not interact in the linear limit. It is argued that for many purposes the intrinsic limitations of covariance-based FC instead favor the alternative of tracking eigenmode coefficients vs. time, which provide a compact representation that is directly related to biophysical brain dynamics.

Neural field theory of neural avalanche exponents

The power-law exponents of observed size and lifetime distributions of near-critical neural avalanches are calculated from neural field theory using diagrammatic methods. This brings neural avalanches within the ambit of neural field theory, which has also previously explained near-critical 1/f spectra and many other observed features of neural activity. This strengthens the case for near-criticality of the brain and opens the way for these other phenomena to be interrelated with avalanches and their dynamics.

Analysis of the Hindriks and van Putten model for propofol anesthesia: Limitations and extensions

We present a detailed analysis of the Hindriks and van Putten thalamocortical mean-field model for propofol anesthesia [NeuroImage 60(23), 2012]. The Hindriks and van Putten (HvP) model predicts increases in delta and alpha power for moderate (up to 130%) prolongation of GABA_A inhibitory response, corresponding to light anesthetic sedation. Our analysis reveals that, for deeper anesthetic effect, the model exhibits an unexpected abrupt jump in cortical activity from a low-firing state to an extremely high-firing stable state (∼250 spikes/s), and remains locked there even at GABA_A prolongations as high as 300% which would be expected to induce full comatose suppression of all firing activity. We demonstrate that this unphysiological behavior can be completely suppressed with appropriate tuning of the parameters controlling the sigmoidal functions that map soma voltage to firing rate for the excitatory and inhibitory neural populations, coupled with elimination of the putative population-dependent anesthetic efficacies introduced in the HvP model. The modifications reported here constrain the anesthetized brain activity into a biologically plausible range in which the cortex now has access to a moderate-firing state ("awake") and a low-firing ("anesthetized") state such that the brain can transition from "awake" to "anesthetized" states at a critical level of drug concentration. The modified HvP model predicts a drug-effect hysteresis in which the drug concentration required for induction is larger than that at emergence. In addition, the revised model shows a decrease in the intensity and frequency of alpha-band fluctuations, transitioning to delta-band dominance, with deepening anesthesia. These predicted drug concentration-dependent changes in EEG dynamics are consistent with clinical reports.

Evoked response activity eigenmode analysis in a convoluted cortex via neural field theory

Neural field theory of the corticothalamic system is used to explore evoked response potentials (ERPs) caused by spatially localized impulse stimuli on the convoluted cortex and on a spherical cortex. Eigenfunctions are calculated analytically on the spherical cortex and numerically on the convoluted cortex via eigenfunction expansions. Eigenmodes on a convoluted cortex are similar to those of the spherical cortex, and a few such modes are found to be sufficient to reproduce the main ERP features. It is found that the ERP peak is stronger in spherical cortex than convoluted cortex, but in both cases the peak decreases monotonically with increasing distance from the stimulus point. In the convoluted case, cortical folding causes ERPs to differ between locations at the same distance from the stimulus point and spherical symmetries are only approximately preserved.

Nonlinear wave-wave interactions in the brain

Neural field theory of the corticothalamic system is used to analyze nonlinear wave-wave interactions in steady state visual evoked potential responses. The nonlinear power spectrum is analytically calculated by convolving the linear power spectrum with itself and other factors. Periodic sine and square wave stimuli are used to generate steady state visual evoked potential responses and to study stimulus-driven nonlinear corticothalamic dynamic interactions. Moreover, we use dual sine drives to analyze the driven dynamics. Numerical analysis shows that the nonlinear power spectrum embodies key nonlinear features, including harmonic and subharmonic generation, entrainment of the alpha rhythm to periodic stimuli at the drive frequency, sum and difference frequencies due to wave-wave coalescence and decay. Further, the scaling properties of the key phenomena observed in nonlinear interactions are studied, verifying some of the theoretical predictions for these being generated by three-wave processes.

Neural Field Theory of Evoked Response Potentials With Attentional Gain Dynamics

A generalized neural field model of large-scale activity in the corticothalamic system is used to predict standard evoked potentials. This model embodies local feedbacks that modulate the gains of neural activity as part of the response to incoming stimuli and thus enables both activity changes and effective connectivity changes to be calculated as parts of a generalized evoked response, and their relative contributions to be determined. The results show that incorporation of gain modulations enables a compact and physically justifiable description of the differences in gain between background-EEG and standard-ERP conditions, with the latter able to be initiated from the background state, rather than requiring distinct parameters as in earlier work. In particular, top-down gains are found to be reduced during an ERP, consistent with recent theoretical suggestions that the role of internal models is diminished in favor of external inputs when the latter change suddenly. The static-gain and modulated-gain system transfer functions are analyzed via control theory in terms of system resonances that were recently shown to implement data filtering whose gain adjustments can be interpreted as attention. These filters are shown to govern early and late features in standard evoked responses and their gain parameters are shown to be dynamically adjusted in a way that implements a form of attention. The results show that dynamically modulated resonant filters responsible for the low-frequency oscillations in an evoked potential response have different parameters than those responsible for low-frequency resting EEG responses, while both responses share similar mid- and high-frequency resonant filters. These results provide a biophysical mechanism by which oscillatory activity in the theta, alpha, and beta frequency ranges of an evoked response are modulated as reflections of attention; notably theta is enhanced and alpha suppressed during the latter parts of the ERP. Furthermore, the model enables the part of the ERP response induced by gain modulations to be estimated and interpreted in terms of attention.

The growth of cognition: Free energy minimization and the embryogenesis of cortical computation

The assumption that during cortical embryogenesis neurons and synaptic connections are selected to form an ensemble maximising synchronous oscillation explains mesoscopic cortical development, and a mechanism for cortical information processing is implied by consistency with the Free Energy Principle and Dynamic Logic. A heteroclinic network emerges, with stable and unstable fixed points of oscillation corresponding to activity in symmetrically connected, versus asymmetrically connected, sets of neurons. Simulations of growth explain a wide range of anatomical observations for columnar and non-columnar cortex, superficial patch connections, and the organization and dynamic interactions of neurone response properties. An antenatal scaffold is created, upon which postnatal learning can establish continuously ordered neuronal representations, permitting matching of co-synchronous fields in multiple cortical areas to solve optimization problems as in Dynamic Logic. Fast synaptic competition partitions equilibria, minimizing "the curse of dimensionality", while perturbations between imperfectly partitioned synchronous fields, under internal reinforcement, enable the cortex to become adaptively self-directed. As learning progresses variational free energy is minimized and entropy bounded.

Effects of physiological parameter evolution on the dynamics of tonic-clonic seizures

The temporal and spectral characteristics of tonic-clonic seizures are investigated using a neural field model of the corticothalamic system in the presence of a temporally varying connection strength between the cerebral cortex and thalamus. Increasing connection strength drives the system into ∼ 10 Hz seizure oscillations once a threshold is passed and a subcritical Hopf bifurcation occurs. In this study, the spectral and temporal characteristics of tonic-clonic seizures are explored as functions of the relevant properties of physiological connection strengths, such as maximum strength, time above threshold, and the ramp rate at which the strength increases or decreases. Analysis shows that the seizure onset time decreases with the maximum connection strength and time above threshold, but increases with the ramp rate. Seizure duration and offset time increase with maximum connection strength, time above threshold, and rate of change. Spectral analysis reveals that the power of nonlinear harmonics and the duration of the oscillations increase as the maximum connection strength and the time above threshold increase. A secondary limit cycle at ∼ 18 Hz, termed a saddle-cycle, is also seen during seizure onset and becomes more prominent and robust with increasing ramp rate. If the time above the threshold is too small, the system does not reach the 10 Hz limit cycle, and only exhibits 18 Hz saddle-cycle oscillations. It is also seen that the time to reach the saturated large amplitude limit-cycle seizure oscillation from both the instability threshold and from the end of the saddle-cycle oscillations is inversely proportional to the square root of the ramp rate.

The impact of driver sleepiness on fixation-related brain potentials

The effects of driver sleepiness are often quantified as deteriorated driving performance, increased blink durations and high levels of subjective sleepiness. Driver sleepiness has also been associated with increasing levels of electroencephalogram (EEG) power, especially in the alpha range. The present exploratory study investigated a new measure of driver sleepiness, the EEG fixation-related lambda response. Thirty young male drivers (23.6 ± 1.7 years old) participated in a driving simulator experiment in which they drove on rural and suburban roads in simulated daylight versus darkness during both the daytime (full sleep) and night-time (sleep deprived). The results show lower lambda responses during night driving and with longer time on task, indicating that sleep deprivation and time on task cause a general decrement in cortical responsiveness to incoming visual stimuli. Levels of subjective sleepiness and line crossings were higher under the same conditions. Furthermore, results of a linear mixed-effects model showed that low lambda responses are associated with high subjective sleepiness and more line crossings. We suggest that the fixation-related lambda response can be used to investigate driving impairment induced by sleep deprivation while driving and that, after further refinement, it may be useful as an objective measure of driver sleepiness.

Neural Field Theory of Corticothalamic Attention With Control System Analysis

Neural field theory is used to analyze attention by extending an existing model of the large-scale activity in the corticothalamic system to incorporate local feedbacks that modulate the gains of neural connectivity as part of the response to incoming stimuli. Treatment of both activity changes and connectivity changes as part of a generalized response enables generalized linear transfer functions of the combined response to be derived. These are then analyzed and interpreted via control theory in terms of stimulus-driven changes in system resonances that were recently shown to implement data filtering and prediction of the inputs. Using simple visual stimuli as a test case, it is shown that the gain response can implement attention by evaluating two main features of the stimuli: the magnitude and the rate of change, by increasing the weight placed on the rate of change in response to sudden changes, while reducing the contribution of stimuli value in tandem. These changes of filter parameters are shown to improve the prediction of the upcoming stimuli based on its recent time course. This outcome is analogous to controller-parameter tuning for performance enhancement in engineering control theory.

Unified dynamics of interictal events and absence seizures

The dynamics of interictal events between absence seizures and their relationship to seizures themselves are investigated by employing a neural field model of the corticothalamic system. Interictal events are modeled as being due to transient parameter excursions beyond the seizure threshold, in the present case by sufficiently temporally varying the connection strength between the cerebral cortex and the thalamus. Increasing connection strength drives the system into ∼3-Hz seizure oscillations via a supercritical Hopf bifurcation once the linear instability threshold is passed. Depending on the time course of the excursion above threshold, different interictal activity event dynamics are seen in the time series of corticothalamic fields. These resemble experimental interictal time series observed via electroencephalography. It is found that the morphology of these events depends on the magnitude and duration of the excursion above threshold. For a large-amplitude excursion of short duration, events resemble interictal spikes, where one large spike is seen, followed by small damped oscillations. For a short excursion with long duration, events like observed interictal periodic sharp waves are seen. When both amplitude and duration above threshold are large, seizure oscillations are seen. Using these outcomes, proximity to seizure can be estimated and tracked.

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.

Optimal Model Parameter Estimation from EEG Power Spectrum Features Observed during General Anesthesia

Mathematical modeling is a powerful tool that enables researchers to describe the experimentally observed dynamics of complex systems. Starting with a robust model including model parameters, it is necessary to choose an appropriate set of model parameters to reproduce experimental data. However, estimating an optimal solution of the inverse problem, i.e., finding a set of model parameters that yields the best possible fit to the experimental data, is a very challenging problem. In the present work, we use different optimization algorithms based on a frequentist approach, as well as Monte Carlo Markov Chain methods based on Bayesian inference techniques to solve the considered inverse problems. We first probe two case studies with synthetic data and study models described by a stochastic non-delayed linear second-order differential equation and a stochastic linear delay differential equation. In a third case study, a thalamo-cortical neural mass model is fitted to the EEG spectral power measured during general anesthesia induced by anesthetics propofol and desflurane. We show that the proposed neural mass model fits very well to the observed EEG power spectra, particularly to the power spectral peaks within δ - (0 - 4 Hz) and α - (8 - 13 Hz) frequency ranges. Furthermore, for each case study, we perform a practical identifiability analysis by estimating the confidence regions of the parameter estimates and interpret the corresponding correlation and sensitivity matrices. Our results indicate that estimating the model parameters from analytically computed spectral power, we are able to accurately estimate the unknown parameters while avoiding the computational costs due to numerical integration of the model equations.

Abnormalities in P300 components in depression: an ERP-sLORETA study

BACKGROUND

Alterations in P300 components occur in depressed patients, but the brain regions contributing to these changes remain unclear.

AIMS

Thus, the aim of the present study was to examine the underlying neural activation of P300 components in patients with depression to explore brain regions related to depression.

METHODS

P300 components were evoked by an oddball auditory paradigm and recorded from 30 patients with current depression, as well as 30 age-, gender-, and education level-matched healthy controls. The standardized Low-Resolution Brain Electromagnetic Tomography (sLORETA) method was used to explore the source activation of P300 components.

RESULTS

Compared with healthy controls, depressed patients tended to exhibit lower P200 and P300 amplitudes and prolonged P300 latency. In depressed patients, P200 source activations were reduced in the right insula, right precentral gyrus, left anterior cingulate, medial frontal gyrus, superior frontal gyrus, and middle frontal gyrus. Decreased source activations of P300 were identified in the right insula, postcentral gyrus, superior temporal gyrus, inferior parietal lobule, transverse temporal gyrus, cingulate gyrus, precentral gyrus, middle frontal gyrus, superior frontal gyrus, medial frontal gyrus, and paracentral gyrus.

CONCLUSIONS

Extensive dysfunction over the right hemisphere and bilateral prefrontal dysfunction may be involved in the pathophysiology of depression.

Suppression of Parkinsonian Beta Oscillations by Deep Brain Stimulation: Determination of Effective Protocols

A neural field model of the corticothalamic-basal ganglia system is developed that describes enhanced beta activity within subthalamic and pallidal circuits in Parkinson's disease (PD) via system resonances. A model of deep brain stimulation (DBS) of typical clinical targets, the subthalamic nucleus (STN) and globus pallidus internus (GPi), is added and studied for several distinct stimulation protocols that are used for treatment of the motor symptoms of PD and that reduce pathological beta band activity (13-30 Hz) in the corticothalamic-basal ganglia network. The resulting impact of DBS on enhanced beta activity in the STN and GPi, as well as cortico-subthalamic and cortico-pallidal coherence, are studied. Both STN-DBS and GPi-DBS are found to be effective for suppressing peak STN and GPi power in the beta band, with GPi-DBS being slightly more effective in both the STN and the GPi for all stimulus protocols tested. The largest decrease in cortico-STN coherence is observed during STN-DBS, whereas GPi-DBS is most effective for reducing cortico-GPi coherence. A reduction of the pathologically large STN connection strengths that define the parkinsonian state results in enhanced 6 Hz activity and could thus represent a compensatory mechanism that has the side effect of driving parkinsonian tremor-like oscillations. This model provides a method for systematically testing effective DBS protocols that agrees with experimental and clinical findings. Furthermore, the model suggests GPi-DBS and STN-DBS have distinct impacts on elevated synchronization between the basal ganglia and motor cortex in PD.

Necessity of the sleep-wake cycle for synaptic homeostasis: system-level analysis of plasticity in the corticothalamic system

Neural field theory is used to study the system-level effects of plasticity in the corticothalamic system, where arousal states are represented parametrically by the connection strengths of the system, among other physiologically based parameters. It is found that the plasticity dynamics have no fixed points or closed cycles in the parameter space of the connection strengths, but parameter subregions exist where flows have opposite signs. Remarkably, these subregions coincide with previously identified regions that correspond to wake and slow-wave sleep, thus demonstrating state dependence of the sign of synaptic modification. We then show that a closed cycle in the parameter space is possible when the plasticity dynamics are driven by the ascending arousal system, which cycles the brain between sleep and wake to complete a closed loop that includes arcs through the opposite-flow subregions. Thus, it is concluded that both wake and sleep are necessary, and together are able to stabilize connection weights in the brain over the daily cycle, thereby providing quantitative realization of the synaptic homeostasis hypothesis.

Physiology-based ERPs in normal and abnormal states

Evoked response potentials (ERPs) and other transients are modeled as impulse responses using physiology-based neural field theory (NFT) of the corticothalamic system of neural activity in the human brain that incorporates synaptic and dendritic dynamics, firing response, axonal propagation, and corticocortical and corticothalamic pathways. The properties of model-predicted ERPs are explored throughout the stability zone of the corticothalamic system, and predicted time series and wavelet spectra are also analyzed. This provides a unified treatment of predicted ERPs for both normal and abnormal states within the brain's stability zone, including likely parameters to represent abnormal states of reduced arousal.

A Blind Module Identification Approach for Predicting Effective Connectivity Within Brain Dynamical Subnetworks

Model-based network discovery measures, such as the brain effective connectivity, require fitting of generative process models to measurements obtained from key areas across the network. For distributed dynamic phenomena, such as generalized seizures and slow-wave sleep, studying effective connectivity from real-time recordings is significantly complicated since (i) outputs from only a subnetwork can be practically measured, and (ii) exogenous subnetwork inputs are unobservable. Model fitting, therefore, constitutes a challenging blind module identification or model inversion problem for finding both the parameters and the many unknown inputs of the subnetwork. We herein propose a novel estimation framework for identifying nonlinear dynamic subnetworks in the case of slowly-varying, otherwise unknown local inputs. Starting with approximate predictions obtained using Cubature Kalman filtering, residuals of local output predictions are utilized to improve upon local input estimates. The algorithm performance is tested on both simulated and clinical EEG of induced seizures under electroconvulsive therapy (ECT). For the simulated network, the algorithm significantly boosted the estimation accuracy for inputs and connections from noisy EEG. For the clinical data, the algorithm predicted increased subnetwork inputs during the pre-stimulus anesthesia condition. Importantly, it predicted an increased frontocentral connectivity during the generalized seizure that is commensurate with electrode placement and that corroborates the clinical hypothesis of increased frontal focality of therapeutic ECT seizures. The proposed framework can be extended to account for several input configurations and can in principle be applied to study effective connectivity within brain subnetworks defined at the microscale (cortical lamina interaction) or at the macroscale (sensory integration).

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.

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.

Cortical geometry as a determinant of brain activity eigenmodes: Neural field analysis

Perturbation analysis of neural field theory is used to derive eigenmodes of neural activity on a cortical hemisphere, which have previously been calculated numerically and found to be close analogs of spherical harmonics, despite heavy cortical folding. The present perturbation method treats cortical folding as a first-order perturbation from a spherical geometry. The first nine spatial eigenmodes on a population-averaged cortical hemisphere are derived and compared with previous numerical solutions. These eigenmodes contribute most to brain activity patterns such as those seen in electroencephalography and functional magnetic resonance imaging. The eigenvalues of these eigenmodes are found to agree with the previous numerical solutions to within their uncertainties. Also in agreement with the previous numerics, all eigenmodes are found to closely resemble spherical harmonics. The first seven eigenmodes exhibit a one-to-one correspondence with their numerical counterparts, with overlaps that are close to unity. The next two eigenmodes overlap the corresponding pair of numerical eigenmodes, having been rotated within the subspace spanned by that pair, likely due to second-order effects. The spatial orientations of the eigenmodes are found to be fixed by gross cortical shape rather than finer-scale cortical properties, which is consistent with the observed intersubject consistency of functional connectivity patterns. However, the eigenvalues depend more sensitively on finer-scale cortical structure, implying that the eigenfrequencies and consequent dynamical properties of functional connectivity depend more strongly on details of individual cortical folding. Overall, these results imply that well-established tools from perturbation theory and spherical harmonic analysis can be used to calculate the main properties and dynamics of low-order brain eigenmodes.

Anesthetic action on the transmission delay between cortex and thalamus explains the beta-buzz observed under propofol anesthesia

In recent years, more and more surgeries under general anesthesia have been performed with the assistance of electroencephalogram (EEG) monitors. An increase in anesthetic concentration leads to characteristic changes in the power spectra of the EEG. Although tracking the anesthetic-induced changes in EEG rhythms can be employed to estimate the depth of anesthesia, their precise underlying mechanisms are still unknown. A prominent feature in the EEG of some patients is the emergence of a strong power peak in the β-frequency band, which moves to the α-frequency band while increasing the anesthetic concentration. This feature is called the beta-buzz. In the present study, we use a thalamo-cortical neural population feedback model to reproduce observed characteristic features in frontal EEG power obtained experimentally during propofol general anesthesia, such as this beta-buzz. First, we find that the spectral power peak in the α- and δ-frequency ranges depend on the decay rate constant of excitatory and inhibitory synapses, but the anesthetic action on synapses does not explain the beta-buzz. Moreover, considering the action of propofol on the transmission delay between cortex and thalamus, the model reveals that the beta-buzz may result from a prolongation of the transmission delay by increasing propofol concentration. A corresponding relationship between transmission delay and anesthetic blood concentration is derived. Finally, an analytical stability study demonstrates that increasing propofol concentration moves the systems resting state towards its stability threshold.

Predicting epileptic seizures from scalp EEG based on attractor state analysis

BACKGROUND AND OBJECTIVE

Epilepsy is the second most common disease of the brain. Epilepsy makes it difficult for patients to live a normal life because it is difficult to predict when seizures will occur. In this regard, if seizures could be predicted a reasonable period of time before their occurrence, epilepsy patients could take precautions against them and improve their safety and quality of life. In this paper, we investigate a novel seizure precursor based on attractor state analysis for seizure prediction.

METHODS

We analyze the transition process from normal to seizure attractor state and investigate a precursor phenomenon seen before reaching the seizure attractor state. From the result of an analysis, we define a quantified spectral measure in scalp EEG for seizure prediction. From scalp EEG recordings, the Fourier coefficients of six EEG frequency bands are extracted, and the defined spectral measure is computed based on the coefficients for each half-overlapped 20-second-long window. The computed spectral measure is applied to seizure prediction using a low-complexity methodology.

RESULTS

Within scalp EEG, we identified an early-warning indicator before an epileptic seizure occurs. Getting closer to the bifurcation point that triggers the transition from normal to seizure state, the power spectral density of low frequency bands of the perturbation of an attractor in the EEG, showed a relative increase. A low-complexity seizure prediction algorithm using this feature was evaluated, using ∼583h of scalp EEG in which 143 seizures in 16 patients were recorded. With the test dataset, the proposed method showed high sensitivity (86.67%) with a false prediction rate of 0.367h^-1 and average prediction time of 45.3min.

CONCLUSIONS

A novel seizure prediction method using scalp EEG, based on attractor state analysis, shows potential for application with real epilepsy patients. This is the first study in which the seizure-precursor phenomenon of an epileptic seizure is investigated based on attractor-based analysis of the macroscopic dynamics of the brain. With the scalp EEG, we first propose use of a spectral feature identified for seizure prediction, in which the dynamics of an attractor are excluded, and only the perturbation dynamics from the attractor are considered.

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.

Mapping how local perturbations influence systems-level brain dynamics

The human brain exhibits a distinct spatiotemporal organization that supports brain function and can be manipulated via local brain stimulation. Such perturbations to local cortical dynamics are globally integrated by distinct neural systems. However, it remains unclear how local changes in neural activity affect large-scale system dynamics. Here, we briefly review empirical and computational studies addressing how localized perturbations affect brain activity. We then systematically analyze a model of large-scale brain dynamics, assessing how localized changes in brain activity at the different sites affect whole-brain dynamics. We find that local stimulation induces changes in brain activity that can be summarized by relatively smooth tuning curves, which relate a region's effectiveness as a stimulation site to its position within the cortical hierarchy. Our results also support the notion that brain hubs, operating in a slower regime, are more resilient to focal perturbations and critically contribute to maintain stability in global brain dynamics. In contrast, perturbations of peripheral regions, characterized by faster activity, have greater impact on functional connectivity. As a parallel with this region-level result, we also find that peripheral systems such as the visual and sensorimotor networks were more affected by local perturbations than high-level systems such as the cingulo-opercular network. Our findings highlight the importance of a periphery-to-core hierarchy to determine the effect of local stimulation on the brain network. This study also provides novel resources to orient empirical work aiming at manipulating functional connectivity using non-invasive brain stimulation.

Neurophysiological Mechanisms of Auditory Information Processing in Adolescence: A Study on Sex Differences

To date, little is known about sex differences in the neurophysiological correlates underlying auditory information processing. In the present study, auditory evoked potentials were evoked in typically developing male (n = 15) and female (n = 14) adolescents (13-18 years) during an auditory oddball task. Girls compared to boys displayed lower N100 and P300 amplitudes to targets. Larger N100 amplitudes in adolescent boys might indicate higher neural sensitivity to changes of incoming auditory information. The P300 findings point toward sex differences in auditory working memory and might suggest that adolescent boys might allocate more attentional resources when processing relevant auditory stimuli than adolescent girls.

Eigenmodes of brain activity: Neural field theory predictions and comparison with experiment

Neural field theory of the corticothalamic system is applied to predict and analyze the activity eigenmodes of the bihemispheric brain, focusing particularly on their spatial structure. The eigenmodes of a single brain hemisphere are found to be close analogs of spherical harmonics, which are the natural modes of the sphere. Instead of multiple eigenvalues being equal, as in the spherical case, cortical folding splits them to have distinct values. Inclusion of interhemispheric connections between homologous regions via the corpus callosum leads to further splitting that depends on symmetry or antisymmetry of activity between brain hemispheres, and the strength and sign of the interhemispheric connections. Symmetry properties of the lowest observed eigenmodes strongly constrain the interhemispheric connectivity strengths and unihemispheric mode spectra, and it is predicted that most spontaneous brain activity will be symmetric between hemispheres, consistent with observations. Comparison with the eigenmodes of an experimental anatomical connectivity matrix confirms these results, permits the relative strengths of intrahemispheric and interhemispheric connectivities to be approximately inferred from their eigenvalues, and lays the foundation for further experimental tests. The results are consistent with brain activity being in corticothalamic eigenmodes, rather than discrete "networks" and open the way to new approaches to brain analysis.

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.

Generalized seizures in a neural field model with bursting dynamics

The mechanisms underlying generalized seizures are explored with neural field theory. A corticothalamic neural field model that has accounted for multiple brain activity phenomena and states is used to explore changes leading to pathological seizure states. It is found that absence seizures arise from instabilities in the system and replicate experimental studies in numerous animal models and clinical studies.

How the cortico-thalamic feedback affects the EEG power spectrum over frontal and occipital regions during propofol-induced sedation

Increasing concentrations of the anaesthetic agent propofol initially induces sedation before achieving full general anaesthesia. During this state of anaesthesia, the observed specific changes in electroencephalographic (EEG) rhythms comprise increased activity in the δ- (0.5-4 Hz) and α- (8-13 Hz) frequency bands over the frontal region, but increased δ- and decreased α-activity over the occipital region. It is known that the cortex, the thalamus, and the thalamo-cortical feedback loop contribute to some degree to the propofol-induced changes in the EEG power spectrum. However the precise role of each structure to the dynamics of the EEG is unknown. In this paper we apply a thalamo-cortical neuronal population model to reproduce the power spectrum changes in EEG during propofol-induced anaesthesia sedation. The model reproduces the power spectrum features observed experimentally both in frontal and occipital electrodes. Moreover, a detailed analysis of the model indicates the importance of multiple resting states in brain activity. The work suggests that the α-activity originates from the cortico-thalamic relay interaction, whereas the emergence of δ-activity results from the full cortico-reticular-relay-cortical feedback loop with a prominent enforced thalamic reticular-relay interaction. This model suggests an important role for synaptic GABAergic receptors at relay neurons and, more generally, for the thalamus in the generation of both the δ- and the α- EEG patterns that are seen during propofol anaesthesia sedation.

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.

Mean-field thalamocortical modeling of longitudinal EEG acquired during intensive meditation training

Meditation training has been shown to enhance attention and improve emotion regulation. However, the brain processes associated with such training are poorly understood and a computational modeling framework is lacking. Modeling approaches that can realistically simulate neurophysiological data while conforming to basic anatomical and physiological constraints can provide a unique opportunity to generate concrete and testable hypotheses about the mechanisms supporting complex cognitive tasks such as meditation. Here we applied the mean-field computational modeling approach using the scalp-recorded electroencephalogram (EEG) collected at three assessment points from meditating participants during two separate 3-month-long shamatha meditation retreats. We modeled cortical, corticothalamic, and intrathalamic interactions to generate a simulation of EEG signals recorded across the scalp. We also present two novel extensions to the mean-field approach that allow for: (a) non-parametric analysis of changes in model parameter values across all channels and assessments; and (b) examination of variation in modeled thalamic reticular nucleus (TRN) connectivity over the retreat period. After successfully fitting whole-brain EEG data across three assessment points within each retreat, two model parameters were found to replicably change across both meditation retreats. First, after training, we observed an increased temporal delay between modeled cortical and thalamic cells. This increase provides a putative neural mechanism for a previously observed reduction in individual alpha frequency in these same participants. Second, we found decreased inhibitory connection strength between the TRN and secondary relay nuclei (SRN) of the modeled thalamus after training. This reduction in inhibitory strength was found to be associated with increased dynamical stability of the model. Altogether, this paper presents the first computational approach, taking core aspects of physiology and anatomy into account, to formally model brain processes associated with intensive meditation training. The observed changes in model parameters inform theoretical accounts of attention training through meditation, and may motivate future study on the use of meditation in a variety of clinical populations.