A neural mass model for the simulation of cortical activity estimated from high resolution EEG during cognitive or motor tasks

Simulated dynamical transitions in a heterogeneous marmoset pFC cluster

Network analysis of the marmoset cortical connectivity data indicates a significant 3D cluster in and around the pre-frontal cortex. A multi-node, heterogeneous neural mass model of this six-node cluster was constructed. Its parameters were informed by available experimental and simulation data so that each neural mass oscillated in a characteristic frequency band. Nodes were connected with directed, weighted links derived from the marmoset structural connectivity data. Heterogeneity arose from the different link weights and model parameters for each node. Stimulation of the cluster with an incident pulse train modulated in the standard frequency bands induced a variety of dynamical state transitions that lasted in the range of 5-10 s, suggestive of timescales relevant to short-term memory. A short gamma burst rapidly reset the beta-induced transition. The theta-induced transition state showed a spontaneous, delayed reset to the resting state. An additional, continuous gamma wave stimulus induced a new beating oscillatory state. Longer or repeated gamma bursts were phase-aligned with the beta oscillation, delivering increasing energy input and causing shorter transition times. The relevance of these results to working memory is yet to be established, but they suggest interesting opportunities.

Control of epileptic activities in a cortex network of multiple coupled neural populations under electromagnetic induction

Epilepsy is believed to be associated with the abnormal synchronous neuronal activity in the brain, which results from large groups or circuits of neurons. In this paper, we choose to focus on the temporal lobe epilepsy, and establish a cortex network of multiple coupled neural populations to explore the epileptic activities under electromagnetic induction. We demonstrate that the epileptic activities can be controlled and modulated by electromagnetic induction and coupling among regions. In certain regions, these two types of control are observed to show exactly reverse effects. The results show that the strong electromagnetic induction is conducive to eliminating the epileptic seizures. The coupling among regions has a conduction effect that the previous normal background activity of the region gives way to the epileptic discharge, owing to coupling with spike wave discharge regions. Overall, these results highlight the role of electromagnetic induction and coupling among the regions in controlling and modulating epileptic activities, and might provide novel insights into the treatments of epilepsy.

Decreased coherence in the model of the dorsal visual pathway associated with Alzheimer's disease

Decreased coherence in electroencephalogram (EEG) has been reported in Alzheimer's disease (AD) experimentally, which could be considered as a typical electrophysiological characteristic in AD. This work aimed to investigate the effect of AD on coherence in the dorsal visual pathway by the technique of neurocomputation. Firstly, according to the hierarchical organization of the cerebral cortex and the information flows of the dorsal visual pathway, a more physiologically plausible neural mass model including cortical areas v1, v2, and v5 was established in the dorsal visual pathway. The three interconnected cortical areas were connected by ascending and descending projections. Next, the pathological condition of loss of long synaptic projections in AD was simulated by reducing the parameters of long synaptic projections in the model. Then, the loss of long synaptic projections on coherence among different visual cortex areas was explored by means of power spectral analysis and coherence function. The results demonstrate that the coherence between these interconnected cortical areas showed an obvious decline with the gradual decrease of long synaptic projections, i.e. decrease in descending projections from area v2 to v1 and v5 to v2 and ascending projection from area v2 to v5. Hopefully, the results of this study could provide theoretical guidance for understanding the dynamical mechanism of AD.

Effect of acetylcholine deficiency on neural oscillation in a brainstem-thalamus-cortex neurocomputational model related with Alzheimer's disease

Previous works imply that involving brainstem in neuropathological studies of Alzheimer's disease (AD) is of clinically significant. This work constructs a comprehensive neural mass model for cholinergic neuropathogenesis that involves brainstem, thalamus and cortex, wherein how acetylcholine deficiency in AD affects neural oscillation of the model output is systematically explored from the perspective of neurocomputation. By decreasing synapse connectivity parameters in direct cholinergic pathway from brainstem to thalamus or in indirect glutamatergic synapse pathway from cortex to brainstem to mimic the pathological condition of reduced acetylcholine release in patients with AD, the property of neural oscillation in this model is numerically investigated by means of power spectrum in frequency domain and amplitude distribution in time domain. Simulated results demonstrate that decreasing synapse connectivity whether in the direct cholinergic pathway or in the indirect glutamatergic synapse pathway can alter the neural oscillation significantly in three aspects: it induces an obvious decrease of dominant frequency; it leads to a degraded rhythmic activity in the alpha frequency band as well as an enhanced rhythmic activity in the theta frequency band; it results in reduced oscillation amplitude of the model output. These results are agreement with the characteristic of electrophysiological EEG measurement recorded in AD, especially support the hypothesis that cholinergic deficiency is a promising pathophysiological origin of EEG slowing in AD. Our analysis indicates that targeting the cholinergic system may have potential prospects in early diagnosis and treatment of AD.

An asynchronous artifact-enhanced electroencephalogram based control paradigm assisted by slight facial expression

In this study, an asynchronous artifact-enhanced electroencephalogram (EEG)-based control paradigm assisted by slight-facial expressions (sFE-paradigm) was developed. The brain connectivity analysis was conducted to reveal the dynamic directional interactions among brain regions under sFE-paradigm. The component analysis was applied to estimate the dominant components of sFE-EEG and guide the signal processing. Enhanced by the artifact within the detected electroencephalogram (EEG), the sFE-paradigm focused on the mainstream defect as the insufficiency of real-time capability, asynchronous logic, and robustness. The core algorithm contained four steps, including “obvious non-sFE-EEGs exclusion,” “interface ‘ON’ detection,” “sFE-EEGs real-time decoding,” and “validity judgment.” It provided the asynchronous function, decoded eight instructions from the latest 100 ms signal, and greatly reduced the frequent misoperation. In the offline assessment, the sFE-paradigm achieved 96.46% ± 1.07 accuracy for interface “ON” detection and 92.68% ± 1.21 for sFE-EEGs real-time decoding, with the theoretical output timespan less than 200 ms. This sFE-paradigm was applied to two online manipulations for evaluating stability and agility. In “object-moving with a robotic arm,” the averaged intersection-over-union was 60.03 ± 11.53%. In “water-pouring with a prosthetic hand,” the average water volume was 202.5 ± 7.0 ml. During online, the sFE-paradigm performed no significant difference (P = 0.6521 and P = 0.7931) with commercial control methods (i.e., FlexPendant and Joystick), indicating a similar level of controllability and agility. This study demonstrated the capability of sFE-paradigm, enabling a novel solution to the non-invasive EEG-based control in real-world challenges.

Realizing the Application of EEG Modeling in BCI Classification: Based on a Conditional GAN Converter

Electroencephalogram (EEG) modeling in brain-computer interface (BCI) provides a theoretical foundation for its development. However, limited by the lack of guidelines in model parameter selection and the inability to obtain personal tissue information in practice, EEG modeling in BCI is mainly focused on the theoretical qualitative level which shows a gap between the theory and its application. Based on such problems, this work combined the surface EEG simulation with a converter based on the generative adversarial network (GAN), to establish the connection from simulated EEG to its application in BCI classification. For the scalp EEGs modeling, a mathematical model was built according to the physics of surface EEG, which consisted of the parallel 3-population neural mass model, the equivalent dipole, and the forward computation. For application, a converter based on the conditional GAN was designed, to transfer the simulated theoretical-only EEG to its practical version, in the lack of individual bio-information. To verify the feasibility, based on the latest microexpression-assisted BCI paradigm proposed by our group, the converted simulated EEGs were used in the training of BCI classifiers. The results indicated that, compared with training with insufficient real data, by adding the simulated EEGs, the overall performance showed a significant improvement (P = 0.04 < 0.05), and the test performance can be improved by 2.17% ± 4.23, in which the largest increase was up to 12.60% ± 1.81. Through this work, the link from theoretical EEG simulation to BCI classification has been initially established, providing an enhanced novel solution for the application of EEG modeling in BCI.

Controlling Alzheimer's Disease Through the Deep Brain Stimulation to Thalamic Relay Cells

Experimental and clinical studies have shown that the technique of deep brain stimulation (DBS) plays a potential role in the regulation of Alzheimer’s disease (AD), yet it still desires for ongoing studies including clinical trials, theoretical approach and action mechanism. In this work, we develop a modified thalamo-cortico-thalamic (TCT) model associated with AD to explore the therapeutic effects of DBS on AD from the perspective of neurocomputation. First, the neuropathological state of AD resulting from synapse loss is mimicked by decreasing the synaptic connectivity strength from the Inter-Neurons (IN) neuron population to the Thalamic Relay Cells (TRC) neuron population. Under such AD condition, a specific deep brain stimulation voltage is then implanted into the neural nucleus of TRC in this TCT model. The symptom of AD is found significantly relieved by means of power spectrum analysis and nonlinear dynamical analysis. Furthermore, the therapeutic effects of DBS on AD are systematically examined in different parameter space of DBS. The results demonstrate that the controlling effect of DBS on AD can be efficient by appropriately tuning the key parameters of DBS including amplitude A, period P and duration D. This work highlights the critical role of thalamus stimulation for brain disease, and provides a theoretical basis for future experimental and clinical studies in treating AD.

Homology Characteristics of EEG and EMG for Lower Limb Voluntary Movement Intention

In the field of lower limb exoskeletons, besides its electromechanical system design and control, attention has been paid to realizing the linkage of exoskeleton robots to humans via electroencephalography (EEG) and electromyography (EMG). However, even the state of the art performance of lower limb voluntary movement intention decoding still faces many obstacles. In the following work, focusing on the perspective of the inner mechanism, a homology characteristic of EEG and EMG for lower limb voluntary movement intention was conducted. A mathematical model of EEG and EMG was built based on its mechanism, which consists of a neural mass model (NMM), neuromuscular junction model, EMG generation model, decoding model, and musculoskeletal biomechanical model. The mechanism analysis and simulation results demonstrated that EEG and EMG signals were both excited by the same movement intention with a response time difference. To assess the efficiency of the proposed model, a synchronous acquisition system for EEG and EMG was constructed to analyze the homology and response time difference from EEG and EMG signals in the limb movement intention. An effective method of wavelet coherence was used to analyze the internal correlation between EEG and EMG signals in the same limb movement intention. To further prove the effectiveness of the hypothesis in this paper, six subjects were involved in the experiments. The experimental results demonstrated that there was a strong EEG-EMG coherence at 1 Hz around movement onset, and the phase of EEG was leading the EMG. Both the simulation and experimental results revealed that EEG and EMG are homologous, and the response time of the EEG signals are earlier than EMG signals during the limb movement intention. This work can provide a theoretical basis for the feasibility of EEG-based pre-perception and fusion perception of EEG and EMG in human movement detection.

NMMGenerator: an automatic neural mass model generator from population graphs

Neural mass models are among the most popular mathematical models of brain activity, since they enable the rapid simulation of large-scale networks involving different neural types at a spatial scale compatible with electrophysiological experiments (e.g. local field potentials). However, establishing neural mass model (NMM) equations associated with specific neuronal network architectures can be tedious and is an error-prone process, restricting their use to scientists who are familiar with mathematics. In order to overcome this challenge, we have developed a user-friendly software that enables a user to construct rapidly, under the form of a graph, a neuronal network with its populations and connectivity patterns. The resulting graph is then automatically translated into the corresponding set of differential equations, which can be solved and displayed within the same software environment. The software is proposed as open access, and should assist in offering the possibility for a wider audience of scientists to develop NMM corresponding to their specific neuroscience research questions.

A New Neural Mass Model Driven Method and Its Application in Early Epileptic Seizure Detection

OBJECTIVE

Despite numerous neural computational models proposed to explain physiological and pathological mechanisms of brain activity, a large gap remains between theory and application of the models. Building on the successful application of data-driven methods in epileptic seizure detection, we aim to build a bridge between data and models in this paper.

METHODS

We first propose a novel model-driven seizure detection method based on dynamic features in epileptic EEGs, where the rationale for dynamic features in epileptic EEGs can be clarified in theory by characterizing the variation of parameters of the model. Then we apply the proposed D&F-model-driven method to the problem of early epileptic seizure detection, where the evolution of model parameters selected and optimized by the proposed method is measured and used to detect the starting point of the seizure.

RESULTS

Numerical results on two open EEG databases demonstrate that our proposed method does a good job of early epileptic seizure detection. The average detection sensitivity, false positive rate and early detection period attain 100%, 0.1/h, and 7.1 s respectively.

CONCLUSION

This paper provides a strategy to characterize EEG signals using a NMM-related method and the model parameters optimized by real EEG may then serve as features in their own right for early seizure detection.

SIGNIFICANCE

An useful attempt to early detect epileptic seizures by combining the neural mass model with data analysis.

Alpha rhythm slowing in a modified thalamo-cortico-thalamic model related with Alzheimer's disease

A decrease in alpha band power is defined as a hallmark of electroencephalogram (EEG) in Alzheimer’s disease (AD). This study devotes to understanding the neuronal correlates of alpha rhythm slowing associated with AD from the view of neurocomputation. Firstly, a modified computational model of thalamo-cortico-thalamic (TCT) circuitry is constructed by incorporating two important biologically plausible ingredients. One is the disinhibition property between different inhibitory interneurons in the cortical module. The other is the full relay function of thalamic relay nucleus (TCR) to the cortical module. Then, by decreasing synaptic connectivity parameters to mimic the neuropathological condition of synapse loss in AD, the correlation between neuronal synaptic behavior and abnormal alpha rhythm is simulated by means of power spectral analysis. The results indicate that these decreases of synaptic activity, i.e., not only the excitatory synaptic connections from TCR to fast inhibitory interneurons C_fte and from excitatory interneurons to pyramidal neurons C_pxe but also the inhibitory synaptic connections from fast inhibitory interneurons to slow inhibitory interneurons C_lfi and from inhibitory interneurons to TCR C_tii, can significantly diminish the peak power density over the alpha band of the thalamic output, which implies that there is a slowing of alpha band. Furthermore, the underlying mechanism behind the alpha rhythmic changes is analyzed using nonlinear dynamical technique. The results reveal that decreases of C_fte, C_pxe, C_lfi and C_tii can make the thalamic module transfer from a limit cycle mode to a point attractor mode, which may lead to the alpha rhythm slowing in the modified TCT model. We expect this work can be helpful in identifying early biomarkers of AD’s EEG and understanding potential pathogenesis of AD.

Seizure tracking of epileptic EEGs using a model-driven approach

OBJECTIVE

As a chronic neurological disorder, epilepsy is characterized by recurrent and unprovoked epileptic seizures that can disrupt the normal neuro-biologic, cognitive, psychological conditions of patients. Therefore, it is worthwhile to give a detailed account of how the epileptic EEG evolves during a period of seizure so that an effective control can be guided for epileptic patients in clinics.

APPROACH

Considering the successful application of the neural mass model (NMM) in exploring the insights into brain activities for epilepsy, in this paper, we aim to construct a model-driven approach to track the development of seizures using epileptic EEGs. We first propose a new time-delay Wendling model with sub-populations (TD-W-SP model) with respect to three aspects of improvements. Then we introduce a model-driven seizure tracking approach, where a model training method is designed based on extracted features from epileptic EEGs and a tracking index is defined as a function of the trained model parameters.

MAIN RESULTS

Numerical results on eight patients on CHB-MIT database demonstrate that our proposed method performs well in simulating epileptic-like EEGs as well as tracking the evolution of three stages (that is, from pre-ictal to ictal and from ictal to post-ictal) during a period of epileptic seizure.

SIGNIFICANCE

A useful attempt to track epileptic seizures by combining the NMM with the data analysis.

A novel neural computational model of generalized periodic discharges in acute hepatic encephalopathy

Acute hepatic encephalopathy (AHE) due to acute liver failure is a common form of delirium, a state of confusion, impaired attention, and decreased arousal. The electroencephalogram (EEG) in AHE often exhibits a striking abnormal pattern of brain activity, which epileptiform discharges repeat in a regular repeating pattern. This pattern is known as generalized periodic discharges, or triphasic-waves (TPWs). While much is known about the neurophysiological mechanisms underlying AHE, how these mechanisms relate to TPWs is poorly understood. In order to develop hypotheses how TPWs arise, our work builds a computational model of AHE (AHE-CM), based on three modifications of the well-studied Liley model which emulate mechanisms believed central to brain dysfunction in AHE: increased neuronal excitability, impaired synaptic transmission, and enhanced postsynaptic inhibition. To relate our AHE-CM to clinical EEG data from patients with AHE, we design a model parameter optimization method based on particle filtering (PF-POM). Based on results from 7 AHE patients, we find that the proposed AHE-CM not only performs well in reproducing important aspects of the EEG, namely the periodicity of triphasic waves (TPWs), but is also helpful in suggesting mechanisms underlying variation in EEG patterns seen in AHE. In particular, our model helps explain what conditions lead to increased frequency of TPWs. In this way, our model represents a starting point for exploring the underlying mechanisms of brain dynamics in delirium by relating microscopic mechanisms to EEG patterns.

Transcallosal Inhibition during Motor Imagery: Analysis of a Neural Mass Model

The EEG rhythmic activities of the somato-sensory cortex reveal event-related desynchronization (ERD) or event-related synchronization (ERS) in beta band (14–30 Hz) as subjects perform certain tasks or react to specific stimuli. Data reported for imagination of movement support the hypothesis that activation of one sensorimotor area (SMA) can be accompanied by deactivation of the other. In order to improve our understanding of beta ERD/ERS generation, two neural mass models (NMM) of a cortical column taken from Wendling et al. (2002) were interconnected to simulate the transmission of information from one cortex to the other. The results show that the excitation of one cortex leads to inhibition of the other and vice versa, enforcing the Theory of Inhibition. This behavior strongly depends on the initial working point (WP) of the neural populations (between the linear and the upper saturation region of a sigmoidal function) and on how the cortical activation or deactivation can move the WP in the upper saturation region ERD or in the linear region ERS, respectively.

Neurophysiological Analysis of the Genesis Mechanism of EEG During the Interictal and Ictal Periods Using a Multiple Neural Masses Model

Electroencephalography (EEG) is an important method to investigate the neurophysiological mechanism underlying epileptogenesis to identify new therapies for the treatment of epilepsy. The neurophysiologically based neural mass model (NMM) can build a bridge between signal processing and neurophysiology, which can be used as a platform to explore the neurophysiological mechanism of epileptogenesis. Most EEG signals cannot be regarded as the outputs of a single NMM with identical model parameters. The outputs of NMM are simple because the diversity of neural signals in the same NMM is ignored. To improve the simulation of EEG signals, a multiple NMM is proposed, the output of which is the linear combination of the outputs of all NMMs. The NMM number is not fixed and is minimized under the premise of guaranteeing the fitting effect. Orthogonal matching pursuit is used to solve a constrained [Formula: see text] norm minimization problem for NMM number and the strength of every NMM. The results showed that the NMM number was significantly lower during the ictal period than during the interictal period, and the strength of major NMMs increased. This indicates that neural masses fuse into fewer larger neural masses with greater strength. The distribution of excitatory and inhibitory strength during the ictal and interictal periods was similar, whereas the excitation/inhibition ratio was higher during the ictal period than during the interictal period.

Topology, Cross-Frequency, and Same-Frequency Band Interactions Shape the Generation of Phase-Amplitude Coupling in a Neural Mass Model of a Cortical Column

Phase-amplitude coupling (PAC), a type of cross-frequency coupling (CFC) where the phase of a low-frequency rhythm modulates the amplitude of a higher frequency, is becoming an important indicator of information transmission in the brain. However, the neurobiological mechanisms underlying its generation remain undetermined. A realistic, yet tractable computational model of the phenomenon is thus needed. Here we analyze a neural mass model of a cortical column, comprising fourteen neuronal populations distributed across four layers (L2/3, L4, L5 and L6). A control analysis showed that the conditional transfer entropy (cTE) measure is able to correctly estimate the flow of information between neuronal populations. Then, we computed cTE from the phases to the amplitudes of the oscillations generated in the cortical column. This approach provides information regarding directionality by distinguishing PAC from APC (amplitude-phase coupling), i.e. the information transfer from amplitudes to phases, and can be used to estimate other types of CFC such as amplitude-amplitude coupling (AAC) and phase-phase coupling (PPC). While experiments often only focus on one or two PAC combinations (e.g., theta-gamma or alpha-gamma), we found that a cortical column can simultaneously generate almost all possible PAC combinations, depending on connectivity parameters, time constants, and external inputs. PAC interactions with and without an anatomical equivalent (direct and indirect interactions, respectively) were analyzed. We found that the strength of PAC between two populations was strongly correlated with the strength of the effective connections between the populations and, on average, did not depend on whether the PAC connection was direct or indirect. When considering a cortical column circuit as a complex network, we found that neuronal populations making indirect PAC connections had, on average, higher local clustering coefficient, efficiency, and betweenness centrality than populations making direct connections and populations not involved in PAC connections. This suggests that their interactions were more effective when transmitting information. Since approximately 60% of the obtained interactions represented indirect connections, our results highlight the importance of the topology of cortical circuits for the generation of the PAC phenomenon. Finally, our results demonstrated that indirect PAC interactions can be explained by a cascade of direct CFC and same-frequency band interactions, suggesting that PAC analysis of experimental data should be accompanied by the estimation of other types of frequency interactions for an integrative understanding of the phenomenon.

Differential frontal theta activity during cognitive and motor tasks

The present study investigated psychophysiological differences between cognitive and motor tasks. Participants were 16 university students (male[Formula: see text][Formula: see text][Formula: see text]10, female[Formula: see text][Formula: see text][Formula: see text]6) aged between 24 to 31 years (mean[Formula: see text][Formula: see text][Formula: see text]27.6 years, SD[Formula: see text][Formula: see text] [Formula: see text] years). They were instructed to perform cognitive and motor tasks while their brain activity was simultaneously recorded using electroencephalography (EEG). A 3(task)[Formula: see text][Formula: see text][Formula: see text]8(area)[Formula: see text][Formula: see text][Formula: see text]4(bands) analysis of variance (ANOVA) was performed on theta, alpha and beta frequency bands. The results showed that the cognitive task resulted in a significantly higher spectral power in the theta band at frontal electrodes than did the motor task. This suggests that the frontal lobe might be engaged in problem solving, resulting in increased theta activity. However, there were no differences in alpha and beta activity between the two tasks. The present study provides psychophysiological evidence for classifying cognitive and motor tasks, which has been a controversial issue for task classifications in motor learning research.

Near-Infrared Spectroscopy - Electroencephalography-Based Brain-State-Dependent Electrotherapy: A Computational Approach Based on Excitation-Inhibition Balance Hypothesis

Stroke is the leading cause of severe chronic disability and the second cause of death worldwide with 15 million new cases and 50 million stroke survivors. The poststroke chronic disability may be ameliorated with early neuro rehabilitation where non-invasive brain stimulation (NIBS) techniques can be used as an adjuvant treatment to hasten the effects. However, the heterogeneity in the lesioned brain will require individualized NIBS intervention where innovative neuroimaging technologies of portable electroencephalography (EEG) and functional-near-infrared spectroscopy (fNIRS) can be leveraged for Brain State Dependent Electrotherapy (BSDE). In this hypothesis and theory article, we propose a computational approach based on excitation–inhibition (E–I) balance hypothesis to objectively quantify the poststroke individual brain state using online fNIRS–EEG joint imaging. One of the key events that occurs following Stroke is the imbalance in local E–I (that is the ratio of Glutamate/GABA), which may be targeted with NIBS using a computational pipeline that includes individual “forward models” to predict current flow patterns through the lesioned brain or brain target region. The current flow will polarize the neurons, which can be captured with E–I-based brain models. Furthermore, E–I balance hypothesis can be used to find the consequences of cellular polarization on neuronal information processing, which can then be implicated in changes in function. We first review the evidence that shows how this local imbalance between E–I leading to functional dysfunction can be restored in targeted sites with NIBS (motor cortex and somatosensory cortex) resulting in large-scale plastic reorganization over the cortex, and probably facilitating recovery of functions. Second, we show evidence how BSDE based on E–I balance hypothesis may target a specific brain site or network as an adjuvant treatment. Hence, computational neural mass model-based integration of neurostimulation with online neuroimaging systems may provide less ambiguous, robust optimization of NIBS, and its application in neurological conditions and disorders across individual patients.

Particle swarm optimization algorithm based parameters estimation and control of epileptiform spikes in a neural mass model

This paper proposes an epilepsy detection and closed-loop control strategy based on Particle Swarm Optimization (PSO) algorithm. The proposed strategy can effectively suppress the epileptic spikes in neural mass models, where the epileptiform spikes are recognized as the biomarkers of transitions from the normal (interictal) activity to the seizure (ictal) activity. In addition, the PSO algorithm shows capabilities of accurate estimation for the time evolution of key model parameters and practical detection for all the epileptic spikes. The estimation effects of unmeasurable parameters are improved significantly compared with unscented Kalman filter. When the estimated excitatory-inhibitory ratio exceeds a threshold value, the epileptiform spikes can be inhibited immediately by adopting the proportion-integration controller. Besides, numerical simulations are carried out to illustrate the effectiveness of the proposed method as well as the potential value for the model-based early seizure detection and closed-loop control treatment design.

A neural mass model of phase-amplitude coupling

Brain activity shows phase-amplitude coupling between its slow and fast oscillatory components. We study phase-amplitude coupling as recorded at individual sites, using a modified version of the well-known Wendling neural mass model. To the population of fast inhibitory interneurons of this model, we added external modulatory input and dynamic self-feedback. These two modifications together are sufficient to let the inhibitory population serve as a limit-cycle oscillator, with frequency characteristics comparable to the beta and gamma bands. The frequency and power of these oscillations can be tuned through the time constant of the dynamic and modulatory input. Alpha band activity is generated, as is usual in such models, as a result of interactions of pyramidal neurons and a population of slow inhibitory interneurons. The slow inhibitory population activity directly influences the fast oscillations via the synaptic gain between slow and fast inhibitory populations. As a result, the amplitude envelope of the fast oscillation is coupled to the phase of the slow activity; this result is consistent with the notion that phase-amplitude coupling is effectuated by interactions between inhibitory interneurons.

A neurophysiological training evaluation metric for air traffic management

The aim of this work was to analyze the possibility to apply a neuroelectrical cognitive metrics for the evaluation of the training level of subjects during the learning of a task employed by Air Traffic Controllers (ATCos). In particular, the Electroencephalogram (EEG), the Electrocardiogram (ECG) and the Electrooculogram (EOG) signals were gathered from a group of students during the execution of an Air Traffic Management (ATM) task, proposed at three different levels of difficulty. The neuroelectrical results were compared with the subjective perception of the task difficulty obtained by the NASA-TLX questionnaires. From these analyses, we suggest that the integration of information derived from the power spectral density (PSD) of the EEG signals, the heart rate (HR) and the eye-blink rate (EBR) return important quantitative information about the training level of the subjects. In particular, by focusing the analysis on the direct and inverse correlation of the frontal PSD theta (4-7 (Hz)) and HR, and of the parietal PSD alpha (10-12 (Hz)) and EBR, respectively, with the degree of mental and emotive engagement, it is possible to obtain useful information about the training improvement across the training sessions.

Modeling the Generation of Phase-Amplitude Coupling in Cortical Circuits: From Detailed Networks to Neural Mass Models

Phase-amplitude coupling (PAC), the phenomenon where the amplitude of a high frequency oscillation is modulated by the phase of a lower frequency oscillation, is attracting an increasing interest in the neuroscience community due to its potential relevance for understanding healthy and pathological information processing in the brain. PAC is a diverse phenomenon, having been experimentally detected in at least ten combinations of rhythms: delta-theta, delta-alpha, delta-beta, delta-gamma, theta-alpha, theta-beta, theta-gamma, alpha-beta, alpha-gamma, and beta-gamma. However, a complete understanding of the biophysical mechanisms generating this diversity is lacking. Here we review computational models of PAC generation that range from detailed models of neuronal networks, where each cell is described by Hodgkin-Huxley-type equations, to neural mass models (NMMs) where only the average activities of neuronal populations are considered. We argue that NMMs are an appropriate mathematical framework (due to the small number of parameters and variables involved and the richness of the dynamics they can generate) to study the PAC phenomenon.

Bidirectional interactions between neuronal and hemodynamic responses to transcranial direct current stimulation (tDCS): challenges for brain-state dependent tDCS

Transcranial direct current stimulation (tDCS) has been shown to modulate cortical neural activity. During neural activity, the electric currents from excitable membranes of brain tissue superimpose in the extracellular medium and generate a potential at scalp, which is referred as the electroencephalogram (EEG). Respective neural activity (energy demand) has been shown to be closely related, spatially and temporally, to cerebral blood flow (CBF) that supplies glucose (energy supply) via neurovascular coupling. The hemodynamic response can be captured by near-infrared spectroscopy (NIRS), which enables continuous monitoring of cerebral oxygenation and blood volume. This neurovascular coupling phenomenon led to the concept of neurovascular unit (NVU) that consists of the endothelium, glia, neurons, pericytes, and the basal lamina. Here, recent works suggest NVU as an integrated system working in concert using feedback mechanisms to enable proper brain homeostasis and function where the challenge remains in capturing these mostly nonlinear spatiotemporal interactions within NVU for brain-state dependent tDCS. In principal accordance, we propose EEG-NIRS-based whole-head monitoring of tDCS-induced neuronal and hemodynamic alterations during tDCS.

Different perception of musical stimuli in patients with monolateral and bilateral cochlear implants

The aim of the present study is to measure the perceived pleasantness during the observation of a musical video clip in a group of cochlear implanted adult patients when compared to a group of normal hearing subjects. This comparison was performed by using the imbalance of the EEG power spectra in alpha band over frontal areas as a metric for the perceived pleasantness. Subjects were asked to watch a musical video clip in three different experimental conditions: with the original audio included (Norm), with a distorted version of the audio (Dist), and without the audio (Mute). The frontal EEG imbalance between the estimated power spectra for the left and right prefrontal areas has been calculated to investigate the differences among the two populations. Results suggested that the perceived pleasantness of the musical video clip in the normal hearing population and in the bilateral cochlear implanted populations has similar range of variation across the different stimulations (Norm, Dist, and Mute), when compared to the range of variation of video clip's pleasantness for the monolateral cochlear implanted population. A similarity exists in the trends of the perceived pleasantness across the different experimental conditions in the mono- and bilaterally cochlear implanted patients.

Neuroelectrical correlates of trustworthiness and dominance judgments related to the observation of political candidates

The present research investigates the neurophysiological activity elicited by fast observations of faces of real candidates during simulated political elections. We used simultaneous recording of electroencephalographic (EEG) signals as well as galvanic skin response (GSR) and heart rate (HR) as measurements of central and autonomic nervous systems. Twenty healthy subjects were asked to give judgments on dominance, trustworthiness, and a preference of vote related to the politicians' faces. We used high-resolution EEG techniques to map statistical differences of power spectral density (PSD) cortical activity onto a realistic head model as well as partial directed coherence (PDC) and graph theory metrics to estimate the functional connectivity networks and investigate the role of cortical regions of interest (ROIs). Behavioral results revealed that judgment of dominance trait is the most predictive of the outcome of the simulated elections. Statistical comparisons related to PSD and PDC values highlighted an asymmetry in the activation of frontal cortical areas associated with the valence of the judged trait as well as to the probability to cast the vote. Overall, our results highlight the existence of cortical EEG features which are correlated with the prediction of vote and with the judgment of trustworthy and dominant faces.

Neurophysiological tools to investigate consumer's gender differences during the observation of TV commercials

Neuromarketing is a multidisciplinary field of research whose aim is to investigate the consumers' reaction to advertisements from a neuroscientific perspective. In particular, the neuroscience field is thought to be able to reveal information about consumer preferences which are unobtainable through conventional methods, including submitting questionnaires to large samples of consumers or performing psychological personal or group interviews. In this scenario, we performed an experiment in order to investigate cognitive and emotional changes of cerebral activity evaluated by neurophysiologic indices during the observation of TV commercials. In particular, we recorded the electroencephalographic (EEG), galvanic skin response (GSR), and heart rate (HR) in a group of 28 healthy subjects during the observation of a series of TV advertisements that have been grouped by commercial categories. Comparisons of cerebral indices have been performed to highlight gender differences between commercial categories and scenes of interest of two specific commercials. Findings show how EEG methodologies, along with the measurements of autonomic variables, could be used to obtain hidden information to marketers not obtainable otherwise. Most importantly, it was suggested how these tools could help to analyse the perception of TV advertisements and differentiate their production according to the consumer's gender.

Engineering a thalamo-cortico-thalamic circuit on SpiNNaker: a preliminary study toward modeling sleep and wakefulness

We present a preliminary study of a thalamo-cortico-thalamic (TCT) implementation on SpiNNaker (Spiking Neural Network architecture), a brain inspired hardware platform designed to incorporate the inherent biological properties of parallelism, fault tolerance and energy efficiency. These attributes make SpiNNaker an ideal platform for simulating biologically plausible computational models. Our focus in this work is to design a TCT framework that can be simulated on SpiNNaker to mimic dynamical behavior similar to Electroencephalogram (EEG) time and power-spectra signatures in sleep-wake transition. The scale of the model is minimized for simplicity in this proof-of-concept study; thus the total number of spiking neurons is ≈1000 and represents a "mini-column" of the thalamocortical tissue. All data on model structure, synaptic layout and parameters is inspired from previous studies and abstracted at a level that is appropriate to the aims of the current study as well as computationally suitable for model simulation on a small 4-chip SpiNNaker system. The initial results from selective deletion of synaptic connectivity parameters in the model show similarity with EEG power spectra characteristics of sleep and wakefulness. These observations provide a positive perspective and a basis for future implementation of a very large scale biologically plausible model of thalamo-cortico-thalamic interactivity-the essential brain circuit that regulates the biological sleep-wake cycle and associated EEG rhythms.

Performance comparison between gPDC and PCMI for measuring directionality of neural information flow

BACKGROUND

General partial directed coherence (gPDC) and permutation conditional mutual information (PCMI) have been widely used to analyze neural activities. These two algorithms are representative of linear and nonlinear methods, respectively. However, there is little known about the difference between their performances in measurements of neural information flow (NIF).

NEW METHOD

Comparison of these two approaches was effectively performed based on the neural mass model (NMM) and real local field potentials.

RESULTS

The results showed that the sensitivity of PCMI was more robust than that of gPDC. The coupling strengths calculated by PCMI were closer to theoretical values in the bidirectional mode of NMM. Furthermore, there was a small Coefficient of Variance (C.V.) for the PCMI results. The gPDC was more sensitive to alterations in the directionality index or the coupling strength of NMM; the gPDC method was more likely to detect a difference between two distinct types of coupling strengths compared to that of PCMI, and gPDC performed well in the identification of the coupling strength in the unidirectional mode.

COMPARISON TO EXISTING METHOD(S)

A comparison between gPDC and PCMI was performed and the advantages of the approaches are discussed.

CONCLUSIONS

The performance of the PCMI is better than that of gPDC in measuring the characteristics of connectivity between neural populations. However, gPDC is recommended to distinguish the differences in connectivity between two states in the same pathway or to detect the coupling strength of the unidirectional mode, such as the hippocampal CA3-CA1 pathway.

EEG oscillatory states: universality, uniqueness and specificity across healthy-normal, altered and pathological brain conditions

For the first time the dynamic repertoires and oscillatory types of local EEG states in 13 diverse conditions (examined over 9 studies) that covered healthy-normal, altered and pathological brain states were quantified within the same methodological and conceptual framework. EEG oscillatory states were assessed by the probability-classification analysis of short-term EEG spectral patterns. The results demonstrated that brain activity consists of a limited repertoire of local EEG states in any of the examined conditions. The size of the state repertoires was associated with changes in cognition and vigilance or neuropsychopathologic conditions. Additionally universal, optional and unique EEG states across 13 diverse conditions were observed. It was demonstrated also that EEG oscillations which constituted EEG states were characteristic for different groups of conditions in accordance to oscillations' functional significance. The results suggested that (a) there is a limit in the number of local states available to the cortex and many ways in which these local states can rearrange themselves and still produce the same global state and (b) EEG individuality is determined by varying proportions of universal, optional and unique oscillatory states. The results enriched our understanding about dynamic microstructure of EEG-signal.

EEG oscillatory states: universality, uniqueness and specificity across healthy-normal, altered and pathological brain conditions

For the first time the dynamic repertoires and oscillatory types of local EEG states in 13 diverse conditions (examined over 9 studies) that covered healthy-normal, altered and pathological brain states were quantified within the same methodological and conceptual framework. EEG oscillatory states were assessed by the probability-classification analysis of short-term EEG spectral patterns. The results demonstrated that brain activity consists of a limited repertoire of local EEG states in any of the examined conditions. The size of the state repertoires was associated with changes in cognition and vigilance or neuropsychopathologic conditions. Additionally universal, optional and unique EEG states across 13 diverse conditions were observed. It was demonstrated also that EEG oscillations which constituted EEG states were characteristic for different groups of conditions in accordance to oscillations' functional significance. The results suggested that (a) there is a limit in the number of local states available to the cortex and many ways in which these local states can rearrange themselves and still produce the same global state and (b) EEG individuality is determined by varying proportions of universal, optional and unique oscillatory states. The results enriched our understanding about dynamic microstructure of EEG-signal.

A thalamo-cortical neural mass model for the simulation of brain rhythms during sleep

Cortico-thalamic interactions are known to play a pivotal role in many brain phenomena, including sleep, attention, memory consolidation and rhythm generation. Hence, simple mathematical models that can simulate the dialogue between the cortex and the thalamus, at a mesoscopic level, have a great cognitive value. In the present work we describe a neural mass model of a cortico-thalamic module, based on neurophysiological mechanisms. The model includes two thalamic populations (a thalamo-cortical relay cell population, TCR, and its related thalamic reticular nucleus, TRN), and a cortical column consisting of four connected populations (pyramidal neurons, excitatory interneurons, inhibitory interneurons with slow and fast kinetics). Moreover, thalamic neurons exhibit two firing modes: bursting and tonic. Finally, cortical synapses among pyramidal neurons incorporate a disfacilitation mechanism following prolonged activity. Simulations show that the model is able to mimic the different patterns of rhythmic activity in cortical and thalamic neurons (beta and alpha waves, spindles, delta waves, K-complexes, slow sleep waves) and their progressive changes from wakefulness to deep sleep, by just acting on modulatory inputs. Moreover, simulations performed by providing short sensory inputs to the TCR show that brain rhythms during sleep preserve the cortex from external perturbations, still allowing a high cortical activity necessary to drive synaptic plasticity and memory consolidation. In perspective, the present model may be used within larger cortico-thalamic networks, to gain a deeper understanding of mechanisms beneath synaptic changes during sleep, to investigate the specific role of brain rhythms, and to explore cortical synchronization achieved via thalamic influences.

Transcranial current brain stimulation (tCS): models and technologies

In this paper, we provide a broad overview of models and technologies pertaining to transcranial current brain stimulation (tCS), a family of related noninvasive techniques including direct current (tDCS), alternating current (tACS), and random noise current stimulation (tRNS). These techniques are based on the delivery of weak currents through the scalp (with electrode current intensity to area ratios of about 0.3-5 A/m2) at low frequencies (typically < 1 kHz) resulting in weak electric fields in the brain (with amplitudes of about 0.2-2 V/m). Here we review the biophysics and simulation of noninvasive, current-controlled generation of electric fields in the human brain and the models for the interaction of these electric fields with neurons, including a survey of in vitro and in vivo related studies. Finally, we outline directions for future fundamental and technological research.

Binding and segmentation via a neural mass model trained with Hebbian and anti-Hebbian mechanisms

Synchronization of neural activity in the gamma band, modulated by a slower theta rhythm, is assumed to play a significant role in binding and segmentation of multiple objects. In the present work, a recent neural mass model of a single cortical column is used to analyze the synaptic mechanisms which can warrant synchronization and desynchronization of cortical columns, during an autoassociation memory task. The model considers two distinct layers communicating via feedforward connections. The first layer receives the external input and works as an autoassociative network in the theta band, to recover a previously memorized object from incomplete information. The second realizes segmentation of different objects in the gamma band. To this end, units within both layers are connected with synapses trained on the basis of previous experience to store objects. The main model assumptions are: (i) recovery of incomplete objects is realized by excitatory synapses from pyramidal to pyramidal neurons in the same object; (ii) binding in the gamma range is realized by excitatory synapses from pyramidal neurons to fast inhibitory interneurons in the same object. These synapses (both at points i and ii) have a few ms dynamics and are trained with a Hebbian mechanism. (iii) Segmentation is realized with faster AMPA synapses, with rise times smaller than 1 ms, trained with an anti-Hebbian mechanism. Results show that the model, with the previous assumptions, can correctly reconstruct and segment three simultaneous objects, starting from incomplete knowledge. Segmentation of more objects is possible but requires an increased ratio between the theta and gamma periods.

A realistic neural mass model of the cortex with laminar-specific connections and synaptic plasticity - evaluation with auditory habituation

In this work we propose a biologically realistic local cortical circuit model (LCCM), based on neural masses, that incorporates important aspects of the functional organization of the brain that have not been covered by previous models: (1) activity dependent plasticity of excitatory synaptic couplings via depleting and recycling of neurotransmitters and (2) realistic inter-laminar dynamics via laminar-specific distribution of and connections between neural populations. The potential of the LCCM was demonstrated by accounting for the process of auditory habituation. The model parameters were specified using Bayesian inference. It was found that: (1) besides the major serial excitatory information pathway (layer 4 to layer 2/3 to layer 5/6), there exists a parallel "short-cut" pathway (layer 4 to layer 5/6), (2) the excitatory signal flow from the pyramidal cells to the inhibitory interneurons seems to be more intra-laminar while, in contrast, the inhibitory signal flow from inhibitory interneurons to the pyramidal cells seems to be both intra- and inter-laminar, and (3) the habituation rates of the connections are unsymmetrical: forward connections (from layer 4 to layer 2/3) are more strongly habituated than backward connections (from Layer 5/6 to layer 4). Our evaluation demonstrates that the novel features of the LCCM are of crucial importance for mechanistic explanations of brain function. The incorporation of these features into a mass model makes them applicable to modeling based on macroscopic data (like EEG or MEG), which are usually available in human experiments. Our LCCM is therefore a valuable building block for future realistic models of human cognitive function.

Modulation of epileptic activity by deep brain stimulation: a model-based study of frequency-dependent effects

A number of studies showed that deep brain stimulation (DBS) can modulate the activity in the epileptic brain and that a decrease of seizures can be achieved in “responding” patients. In most of these studies, the choice of stimulation parameters is critical to obtain desired clinical effects. In particular, the stimulation frequency is a key parameter that is difficult to tune. A reason is that our knowledge about the frequency-dependant mechanisms according to which DBS indirectly impacts the dynamics of pathological neuronal systems located in the neocortex is still limited. We address this issue using both computational modeling and intracerebral EEG (iEEG) data. We developed a macroscopic (neural mass) model of the thalamocortical network. In line with already-existing models, it includes interconnected neocortical pyramidal cells and interneurons, thalamocortical cells and reticular neurons. The novelty was to introduce, in the thalamic compartment, the biophysical effects of direct stimulation. Regarding clinical data, we used a quite unique data set recorded in a patient (drug-resistant epilepsy) with a focal cortical dysplasia (FCD). In this patient, DBS strongly reduced the sustained epileptic activity of the FCD for low-frequency (LFS, < 2 Hz) and high-frequency stimulation (HFS, > 70 Hz) while intermediate-frequency stimulation (IFS, around 50 Hz) had no effect. Signal processing, clustering, and optimization techniques allowed us to identify the necessary conditions for reproducing, in the model, the observed frequency-dependent stimulation effects. Key elements which explain the suppression of epileptic activity in the FCD include: (a) feed-forward inhibition and synaptic short-term depression of thalamocortical connections at LFS, and (b) inhibition of the thalamic output at HFS. Conversely, modeling results indicate that IFS favors thalamic oscillations and entrains epileptic dynamics.

Laminar circuit organization and response modulation in mouse visual cortex

The mouse has become an increasingly important animal model for visual system studies, but few studies have investigated local functional circuit organization of mouse visual cortex. Here we used our newly developed mapping technique combining laser scanning photostimulation (LSPS) with fast voltage-sensitive dye (VSD) imaging to examine the spatial organization and temporal dynamics of laminar circuit responses in living slice preparations of mouse primary visual cortex (V1). During experiments, LSPS using caged glutamate provided spatially restricted neuronal activation in a specific cortical layer, and evoked responses from the stimulated layer to its functionally connected regions were detected by VSD imaging. In this study, we first provided a detailed analysis of spatiotemporal activation patterns at specific V1 laminar locations and measured local circuit connectivity. Then we examined the role of cortical inhibition in the propagation of evoked cortical responses by comparing circuit activity patterns in control and in the presence of GABAa receptor antagonists. We found that GABAergic inhibition was critical in restricting layer-specific excitatory activity spread and maintaining topographical projections. In addition, we investigated how AMPA and NMDA receptors influenced cortical responses and found that blocking AMPA receptors abolished interlaminar functional projections, and the NMDA receptor activity was important in controlling visual cortical circuit excitability and modulating activity propagation. The NMDA receptor antagonist reduced neuronal population activity in time-dependent and laminar-specific manners. Finally, we used the quantitative information derived from the mapping experiments and presented computational modeling analysis of V1 circuit organization. Taken together, the present study has provided important new information about mouse V1 circuit organization and response modulation.

Critical fluctuations in cortical models near instability

Computational studies often proceed from the premise that cortical dynamics operate in a linearly stable domain, where fluctuations dissipate quickly and show only short memory. Studies of human electroencephalography (EEG), however, have shown significant autocorrelation at time lags on the scale of minutes, indicating the need to consider regimes where non-linearities influence the dynamics. Statistical properties such as increased autocorrelation length, increased variance, power law scaling, and bistable switching have been suggested as generic indicators of the approach to bifurcation in non-linear dynamical systems. We study temporal fluctuations in a widely-employed computational model (the Jansen–Rit model) of cortical activity, examining the statistical signatures that accompany bifurcations. Approaching supercritical Hopf bifurcations through tuning of the background excitatory input, we find a dramatic increase in the autocorrelation length that depends sensitively on the direction in phase space of the input fluctuations and hence on which neuronal subpopulation is stochastically perturbed. Similar dependence on the input direction is found in the distribution of fluctuation size and duration, which show power law scaling that extends over four orders of magnitude at the Hopf bifurcation. We conjecture that the alignment in phase space between the input noise vector and the center manifold of the Hopf bifurcation is directly linked to these changes. These results are consistent with the possibility of statistical indicators of linear instability being detectable in real EEG time series. However, even in a simple cortical model, we find that these indicators may not necessarily be visible even when bifurcations are present because their expression can depend sensitively on the neuronal pathway of incoming fluctuations.

Influence of connection type on phase synchrony: analysis of a neural mass model

Empirical studies have demonstrated synchronized frontal and parietal electrophysiological signals at 22-34 Hz during a conjunctive visual search task and at 36-56 Hz during a pop-out visual search task. Bidirectional (conjunctive) versus unidirectional (pop-out) information transfer between neuronal populations is hypothesized to underly this difference in synchronization frequency. This study modeled the influence of connection type (i.e., unidirectional vs. bidirectional) on phase synchrony between two neural populations using a neural mass model. Phase-locking values (PLVs) were used as the measure of synchrony between populations. Consistent with the connectivity hypothesis, the model revealed greater PLVs at 22-34 Hz when the two populations were connected bidirectionally than unidirectionally, but greater PLVs at 34-52 Hz when connected unidirectionally than bidirectionally. The model suggests that inter-population connectivity also changes with bottom-up versus top-down control of attention.

Towards a model-based integration of co-registered electroencephalography/functional magnetic resonance imaging data with realistic neural population meshes

Brain activity can be measured with several non-invasive neuroimaging modalities, but each modality has inherent limitations with respect to resolution, contrast and interpretability. It is hoped that multimodal integration will address these limitations by using the complementary features of already available data. However, purely statistical integration can prove problematic owing to the disparate signal sources. As an alternative, we propose here an advanced neural population model implemented on an anatomically sound cortical mesh with freely adjustable connectivity, which features proper signal expression through a realistic head model for the electroencephalogram (EEG), as well as a haemodynamic model for functional magnetic resonance imaging based on blood oxygen level dependent contrast (fMRI BOLD). It hence allows simultaneous and realistic predictions of EEG and fMRI BOLD from the same underlying model of neural activity. As proof of principle, we investigate here the influence on simulated brain activity of strengthening visual connectivity. In the future we plan to fit multimodal data with this neural population model. This promises novel, model-based insights into the brain's activity in sleep, rest and task conditions.

On the use of EEG or MEG brain imaging tools in neuromarketing research

Here we present an overview of some published papers of interest for the marketing research employing electroencephalogram (EEG) and magnetoencephalogram (MEG) methods. The interest for these methodologies relies in their high-temporal resolution as opposed to the investigation of such problem with the functional Magnetic Resonance Imaging (fMRI) methodology, also largely used in the marketing research. In addition, EEG and MEG technologies have greatly improved their spatial resolution in the last decades with the introduction of advanced signal processing methodologies. By presenting data gathered through MEG and high resolution EEG we will show which kind of information it is possible to gather with these methodologies while the persons are watching marketing relevant stimuli. Such information will be related to the memorization and pleasantness related to such stimuli. We noted that temporal and frequency patterns of brain signals are able to provide possible descriptors conveying information about the cognitive and emotional processes in subjects observing commercial advertisements. These information could be unobtainable through common tools used in standard marketing research. We also show an example of how an EEG methodology could be used to analyze cultural differences between fruition of video commercials of carbonated beverages in Western and Eastern countries.

Dynamic causal models and physiological inference: a validation study using isoflurane anaesthesia in rodents

Generative models of neuroimaging and electrophysiological data present new opportunities for accessing hidden or latent brain states. Dynamic causal modeling (DCM) uses Bayesian model inversion and selection to infer the synaptic mechanisms underlying empirically observed brain responses. DCM for electrophysiological data, in particular, aims to estimate the relative strength of synaptic transmission at different cell types and via specific neurotransmitters. Here, we report a DCM validation study concerning inference on excitatory and inhibitory synaptic transmission, using different doses of a volatile anaesthetic agent (isoflurane) to parametrically modify excitatory and inhibitory synaptic processing while recording local field potentials (LFPs) from primary auditory cortex (A1) and the posterior auditory field (PAF) in the auditory belt region in rodents. We test whether DCM can infer, from the LFP measurements, the expected drug-induced changes in synaptic transmission mediated via fast ionotropic receptors; i.e., excitatory (glutamatergic) AMPA and inhibitory GABA(A) receptors. Cross- and auto-spectra from the two regions were used to optimise three DCMs based on biologically plausible neural mass models and specific network architectures. Consistent with known extrinsic connectivity patterns in sensory hierarchies, we found that a model comprising forward connections from A1 to PAF and backward connections from PAF to A1 outperformed a model with forward connections from PAF to A1 and backward connections from A1 to PAF and a model with reciprocal lateral connections. The parameter estimates from the most plausible model indicated that the amplitude of fast glutamatergic excitatory postsynaptic potentials (EPSPs) and inhibitory postsynaptic potentials (IPSPs) behaved as predicted by previous neurophysiological studies. Specifically, with increasing levels of anaesthesia, glutamatergic EPSPs decreased linearly, whereas fast GABAergic IPSPs displayed a nonlinear (saturating) increase. The consistency of our model-based in vivo results with experimental in vitro results lends further validity to the capacity of DCM to infer on synaptic processes using macroscopic neurophysiological data.

A thalamo-cortico-thalamic neural mass model to study alpha rhythms in Alzheimer's disease

We present a lumped computational model of the thalamo-cortico-thalamic circuitry. The model essentially consists of two modules: a thalamic module and a cortical module. The thalamic module circuitry is a modified version of a classic neural mass computational model of the thalamic circuitry to simulate cortical alpha rhythms and which we have used in previous research to study EEG abnormality associated with Alzheimer's Disease (AD). Here, we introduce a modified synaptic structure representing a neuronal population in the thalamic model. Furthermore, the synaptic organisation and connectivity parameter values in the model are based on experimental data reported from the dorsal Lateral Geniculate Nucleus of different species. The cortical module circuitry is based on a recent work studying cortical brain rhythms. We vary the synaptic connectivity parameters in the thalamic module of the model to simulate the effects of AD on brain synaptic circuitry and study power within the alpha frequency bands. The power and dominant frequencies of the model output are studied in three sub-bands within the alpha band: lower alpha (7-9 Hz), middle alpha (9-11 Hz) and upper alpha (11-13 Hz). Such an analytical method conforms to recent comparative EEG studies on young adults, healthy aged adults and MCI or early stage AD patients. The results show a remarkable role of the synaptic connectivities in the inhibitory thalamic cell populations on the alpha band power and frequency. Furthermore, the total number of active synapses in the thalamic cell populations produces the slowing of alpha rhythms and a simultaneous decrease of alpha band power in the brain as a result of AD.

Intermittent spike-wave dynamics in a heterogeneous, spatially extended neural mass model

Generalised epileptic seizures are frequently accompanied by sudden, reversible transitions from low amplitude, irregular background activity to high amplitude, regular spike-wave discharges (SWD) in the EEG. The underlying mechanisms responsible for SWD generation and for the apparently spontaneous transitions to SWD and back again are still not fully understood. Specifically, the role of spatial cortico-cortical interactions in ictogenesis is not well studied. We present a macroscopic, neural mass model of a cortical column which includes two distinct time scales of inhibition. This model can produce both an oscillatory background and a pathological SWD rhythm. We demonstrate that coupling two of these cortical columns can lead to a bistability between out-of-phase, low amplitude background dynamics and in-phase, high amplitude SWD activity. Stimuli can cause state-dependent transitions from background into SWD. In an extended local area of cortex, spatial heterogeneities in a model parameter can lead to spontaneous reversible transitions from a desynchronised background to synchronous SWD due to intermittency. The deterministic model is therefore capable of producing absence seizure-like events without any time dependent adjustment of model parameters. The emergence of such mechanisms due to spatial coupling demonstrates the importance of spatial interactions in modelling ictal dynamics, and in the study of ictogenesis.

Dynamical principles of emotion-cognition interaction: mathematical images of mental disorders

The key contribution of this work is to introduce a mathematical framework to understand self-organized dynamics in the brain that can explain certain aspects of itinerant behavior. Specifically, we introduce a model based upon the coupling of generalized Lotka-Volterra systems. This coupling is based upon competition for common resources. The system can be regarded as a normal or canonical form for any distributed system that shows self-organized dynamics that entail winnerless competition. Crucially, we will show that some of the fundamental instabilities that arise in these coupled systems are remarkably similar to endogenous activity seen in the brain (using EEG and fMRI). Furthermore, by changing a small subset of the system's parameters we can produce bifurcations and metastable sequential dynamics changing, which bear a remarkable similarity to pathological brain states seen in psychiatry. In what follows, we will consider the coupling of two macroscopic modes of brain activity, which, in a purely descriptive fashion, we will label as cognitive and emotional modes. Our aim is to examine the dynamical structures that emerge when coupling these two modes and relate them tentatively to brain activity in normal and non-normal states.

Connecting mean field models of neural activity to EEG and fMRI data

Progress in functional neuroimaging of the brain increasingly relies on the integration of data from complementary imaging modalities in order to improve spatiotemporal resolution and interpretability. However, the usefulness of merely statistical combinations is limited, since neural signal sources differ between modalities and are related non-trivially. We demonstrate here that a mean field model of brain activity can simultaneously predict EEG and fMRI BOLD with proper signal generation and expression. Simulations are shown using a realistic head model based on structural MRI, which includes both dense short-range background connectivity and long-range specific connectivity between brain regions. The distribution of modeled neural masses is comparable to the spatial resolution of fMRI BOLD, and the temporal resolution of the modeled dynamics, importantly including activity conduction, matches the fastest known EEG phenomena. The creation of a cortical mean field model with anatomically sound geometry, extensive connectivity, and proper signal expression is an important first step towards the model-based integration of multimodal neuroimages.

The generation of rhythms within a cortical region: analysis of a neural mass model

Rhythms in brain electrical activity are assumed to play a significant role in many cognitive and perceptual processes. It is thus of great value to analyze these rhythms and their mutual relationships in large scale models of cortical regions. In the present work, we modified the neural mass model by Wendling et al. (Eur. J. Neurosci. 15 (2002) 1499-1508) by including a new inhibitory self-loop among GABAA,fast interneurons. A theoretical analysis was performed to demonstrate that, thanks to this loop, GABAA,fast interneurons can produce a gamma rhythm in the power spectral density (PSD) even without the participation of the other neural populations. Then, the model of a whole cortical region, built upon four interconnected neural populations (pyramidal cells, excitatory, GABAA,slow and GABAA,fast interneurons) was investigated by changing the internal connectivity parameters. Results show that different rhythm combinations (beta and gamma, alpha and gamma, or a wide spectrum) can be obtained within the same region by simply altering connectivity values, without the need to change synaptic kinetics. Finally, two or three cortical regions were connected by using different topologies of long range connections. Results show that long-range connections directed from pyramidal neurons to GABAA,fast interneurons are the most efficient to transmit rhythms from one region to another. In this way, PSD with three or four peaks can be obtained using simple connectivity patterns. The model can be of value to gain a deeper insight into the mechanisms involved in the generation of gamma rhythms and provide a better understanding of cortical EEG spectra.

A neural mass model to simulate different rhythms in a cortical region

An original neural mass model of a cortical region has been used to investigate the origin of EEG rhythms. The model consists of four interconnected neural populations: pyramidal cells, excitatory interneurons and inhibitory interneurons with slow and fast synaptic kinetics, GABA(A, slow) and GABA(A,fast) respectively. A new aspect, not present in previous versions, consists in the inclusion of a self-loop among GABA(A,fast) interneurons. The connectivity parameters among neural populations have been changed in order to reproduce different EEG rhythms. Moreover, two cortical regions have been connected by using different typologies of long range connections. Results show that the model of a single cortical region is able to simulate the occurrence of multiple power spectral density (PSD) peaks; in particular the new inhibitory loop seems to have a critical role in the activation in gamma (gamma) band, in agreement with experimental studies. Moreover the effect of different kinds of connections between two regions has been investigated, suggesting that long range connections toward GABA(A,fast) interneurons have a major impact than connections toward pyramidal cells. The model can be of value to gain a deeper insight into mechanisms involved in the generation of gamma rhythms and to provide better understanding of cortical EEG spectra.