COALIA: A Computational Model of Human EEG for Consciousness Research

Transient alteration of Awareness triggered by direct electrical stimulation of the brain

BACKGROUND

Awareness is a state of consciousness that enables a subject to interact with the environment. Transient alteration of awareness (AA) is a disabling sign of many types of epileptic seizures. The brain mechanisms of awareness and its alteration are not well known.

OBJECTIVE/HYPOTHESIS

Transient and isolated AA induced by electrical brain stimulation during a stereoelectroencephalography (SEEG) recording represents an ideal model for studying the associated modifications of functional connectivity and locating the hubs of awareness networks.

METHODS

We investigated the SEEG signals-based brain functional connectivity (FC) changes vs background occurring during AA triggered by three thalamic and two insular stimulations in three patients explored by SEEG in the frame of presurgical evaluation for focal drug-resistant epilepsy. The results were compared to the stimulations of the same sites that did not induce clinical changes (negative stimulations).

RESULTS

We observed decreased node strength in the pulvinar, insula, and parietal associative cortices during the thalamic and insular stimulations that induced AA. The link strengths characterizing functional coupling between the thalamus and the insular, prefrontal, temporal, or parietal associative cortices were also decreased. In contrast, there was an increased synchronization between the precuneus and the temporal lateral cortex. These FC changes were absent during the negative stimulations.

CONCLUSION

Our study highlights the role of the pulvinar, insular, and parietal hubs in maintaining the awareness networks and paves the way for invasive or non-invasive neuromodulation protocols to reduce AA manifestations during epileptic seizures.

Whole-brain simulation of interictal epileptic discharges for patient-specific interpretation of interictal SEEG data

In patients with refractory epilepsy, the clinical interpretation of stereoelectroencephalographic (SEEG) signals is crucial to delineate the epileptogenic network that should be targeted by surgery. We propose a pipeline of patient-specific computational modeling of interictal epileptic activity to improve the definition of regions of interest. Comparison between the computationally defined regions of interest and the resected region confirmed the efficiency of the pipeline. This result suggests that computational modeling can be used to reconstruct signals and aid clinical interpretation.

Modeling the contribution of theta-gamma coupling to sequential memory, imagination, and dreaming

Gamma oscillations nested in a theta rhythm are observed in the hippocampus, where are assumed to play a role in sequential episodic memory, i.e., memorization and retrieval of events that unfold in time. In this work, we present an original neurocomputational model based on neural masses, which simulates the encoding of sequences of events in the hippocampus and subsequent retrieval by exploiting the theta-gamma code. The model is based on a three-layer structure in which individual Units oscillate with a gamma rhythm and code for individual features of an episode. The first layer (working memory in the prefrontal cortex) maintains a cue in memory until a new signal is presented. The second layer (CA3 cells) implements an auto-associative memory, exploiting excitatory and inhibitory plastic synapses to recover an entire episode from a single feature. Units in this layer are disinhibited by a theta rhythm from an external source (septum or Papez circuit). The third layer (CA1 cells) implements a hetero-associative net with the previous layer, able to recover a sequence of episodes from the first one. During an encoding phase, simulating high-acetylcholine levels, the network is trained with Hebbian (synchronizing) and anti-Hebbian (desynchronizing) rules. During retrieval (low-acetylcholine), the network can correctly recover sequences from an initial cue using gamma oscillations nested inside the theta rhythm. Moreover, in high noise, the network isolated from the environment simulates a mind-wandering condition, randomly replicating previous sequences. Interestingly, in a state simulating sleep, with increased noise and reduced synapses, the network can "dream" by creatively combining sequences, exploiting features shared by different episodes. Finally, an irrational behavior (erroneous superimposition of features in various episodes, like "delusion") occurs after pathological-like reduction in fast inhibitory synapses. The model can represent a straightforward and innovative tool to help mechanistically understand the theta-gamma code in different mental states.

Multiscale neuro-inspired models for interpretation of EEG signals in patients with epilepsy

OBJECTIVE

The aim is to gain insight into the pathophysiological mechanisms underlying interictal epileptiform discharges observed in electroencephalographic (EEG) and stereo-EEG (SEEG, depth electrodes) recordings performed during pre-surgical evaluation of patients with drug-resistant epilepsy.

METHODS

We developed novel neuro-inspired computational models of the human cerebral cortex at three different levels of description: i) microscale (detailed neuron models), ii) mesoscale (neuronal mass models) and iii) macroscale (whole brain models). Although conceptually different, micro- and mesoscale models share some similar features, such as the typology of neurons (pyramidal cells and three types of interneurons), their spatial arrangement in cortical layers, and their synaptic connectivity (excitatory and inhibitory). The whole brain model consists of a large-scale network of interconnected neuronal masses, with connectivity based on the human connectome.

RESULTS

For these three levels of description, the fine-tuning of free parameters and the quantitative comparison with real data allowed us to reproduce interictal epileptiform discharges with a high degree of fidelity and to formulate hypotheses about the cell- and network-related mechanisms underlying the generation of fast ripples and SEEG-recorded epileptic spikes and spike-waves.

CONCLUSIONS

The proposed models provide valuable insights into the pathophysiological mechanisms underlying the generation of epileptic events. The knowledge gained from these models effectively complements the clinical analysis of SEEG data collected during the evaluation of patients with epilepsy.

SIGNIFICANCE

These models are likely to play a key role in the mechanistic interpretation of epileptiform activity.

Technological Modalities in the Assessment and Treatment of Disorders of Consciousness

Over the last 10 years, there have been rapid advances made in technologies that can be utilized in the diagnosis and treatment of patients with a disorder of consciousness (DoC). This article provides a comprehensive review of these modalities including the evidence supporting their potential use in DoC. This review specifically addresses diagnostic, non-invasive therapeutic, and invasive therapeutic technological modalities except for neuroimaging, which is discussed in another article. While technologic advances appear promising for both assessment and treatment of patients with a DoC, high-quality evidence supporting widespread clinical adoption remains limited.

Bi-sinusoidal light stimulation reveals an enhanced response power and reduced phase coherence at the visual cortex in migraine

Introduction

Migraine is associated with enhanced visual sensitivity during and outside attacks. Processing of visual information is a highly non-linear process involving complex interactions across (sub)cortical networks. In this exploratory study, we combined electroencephalography with bi-sinusoidal light stimulation to assess non-linear features of visual processing in participants with migraine.

Methods

Twenty participants with migraine (10 with aura, 10 without aura) and ten non-headache controls were measured (outside attacks). Participants received bi-sinusoidal 13 + 23 Hz red light visual stimulation. Electroencephalography spectral power and multi-spectral phase coherence were compared between groups at the driving stimulation frequencies together with multiples and combinations of these frequencies (harmonic and intermodulation frequencies) caused by non-linearities.

Results

Only at the driving frequency of 13 Hz higher spectral power was found in migraine with aura participants compared with those with migraine without aura and controls. Differences in phase coherence were present for 2nd, 4th, and 5th-order non-linearities in those with migraine (migraine with and without aura) compared with controls. Bi-sinusoidal light stimulation revealed evident non-linearities in the brain's electroencephalography response up to the 5th order with reduced phase coherence for higher order interactions in interictal participants with migraine.

Discussion

Insight into interictal non-linear visual processing may help understand brain dynamics underlying migraine attack susceptibility. Future research is needed to determine the clinical value of the results.

Identifying good practices for detecting inter-regional linear functional connectivity from EEG

Aggregating voxel-level statistical dependencies between multivariate time series is an important intermediate step when characterising functional connectivity (FC) between larger brain regions. However, there are numerous ways in which voxel-level data can be aggregated into inter-regional FC, and the advantages of each of these approaches are currently unclear. In this study we generate ground-truth data and compare the performances of various pipelines that estimate directed and undirected linear phase-to-phase FC between regions. We test the ability of several existing and novel FC analysis pipelines to identify the true regions within which connectivity was simulated. We test various inverse modelling algorithms, strategies to aggregate time series within regions, and connectivity metrics. Furthermore, we investigate the influence of the number of interactions, the signal-to-noise ratio, the noise mix, the interaction time delay, and the number of active sources per region on the ability of detecting phase-to-phase FC. Throughout all simulated scenarios, lowest performance is obtained with pipelines involving the absolute value of coherency. Further, the combination of dynamic imaging of coherent sources (DICS) beamforming with directed FC metrics that aggregate information across multiple frequencies leads to unsatisfactory results. Pipelines that show promising results with our simulated pseudo-EEG data involve the following steps: (1) Source projection using the linearly-constrained minimum variance (LCMV) beamformer. (2) Principal component analysis (PCA) using the same fixed number of components within every region. (3) Calculation of the multivariate interaction measure (MIM) for every region pair to assess undirected phase-to-phase FC, or calculation of time-reversed Granger Causality (TRGC) to assess directed phase-to-phase FC. We formulate recommendations based on these results that may increase the validity of future experimental connectivity studies. We further introduce the free ROIconnect plugin for the EEGLAB toolbox that includes the recommended methods and pipelines that are presented here. We show an exemplary application of the best performing pipeline to the analysis of EEG data recorded during motor imagery.

Computational modelling in disorders of consciousness: Closing the gap towards personalised models for restoring consciousness

Disorders of consciousness are complex conditions characterised by persistent loss of responsiveness due to brain injury. They present diagnostic challenges and limited options for treatment, and highlight the urgent need for a more thorough understanding of how human consciousness arises from coordinated neural activity. The increasing availability of multimodal neuroimaging data has given rise to a wide range of clinically- and scientifically-motivated modelling efforts, seeking to improve data-driven stratification of patients, to identify causal mechanisms for patient pathophysiology and loss of consciousness more broadly, and to develop simulations as a means of testing in silico potential treatment avenues to restore consciousness. As a dedicated Working Group of clinicians and neuroscientists of the international Curing Coma Campaign, here we provide our framework and vision to understand the diverse statistical and generative computational modelling approaches that are being employed in this fast-growing field. We identify the gaps that exist between the current state-of-the-art in statistical and biophysical computational modelling in human neuroscience, and the aspirational goal of a mature field of modelling disorders of consciousness; which might drive improved treatments and outcomes in the clinic. Finally, we make several recommendations for how the field as a whole can work together to address these challenges.

TMS-evoked responses are driven by recurrent large-scale network dynamics

A compelling way to disentangle the complexity of the brain is to measure the effects of spatially and temporally synchronized systematic perturbations. In humans, this can be non-invasively achieved by combining transcranial magnetic stimulation (TMS) and electroencephalography (EEG). Spatiotemporally complex and long-lasting TMS-EEG evoked potential (TEP) waveforms are believed to result from recurrent, re-entrant activity that propagates broadly across multiple cortical and subcortical regions, dispersing from and later re-converging on, the primary stimulation site. However, if we loosely understand the TEP of a TMS-stimulated region as the impulse response function of a noisy underdamped harmonic oscillator, then multiple later activity components (waveform peaks) should be expected even for an isolated network node in the complete absence of recurrent inputs. Thus emerges a critically important question for basic and clinical research on human brain dynamics: what parts of the TEP are due to purely local dynamics, what parts are due to reverberant, re-entrant network activity, and how can we distinguish between the two? To disentangle this, we used source-localized TMS-EEG analyses and whole-brain connectome-based computational modelling. Results indicated that recurrent network feedback begins to drive TEP responses from 100 ms post-stimulation, with earlier TEP components being attributable to local reverberatory activity within the stimulated region. Subject-specific estimation of neurophysiological parameters additionally indicated an important role for inhibitory GABAergic neural populations in scaling cortical excitability levels, as reflected in TEP waveform characteristics. The novel discoveries and new software technologies introduced here should be of broad utility in basic and clinical neuroscience research.

Effect of channel density, inverse solutions and connectivity measures on EEG resting-state networks reconstruction: A simulation study

Along with the study of brain activity evoked by external stimuli, the past two decades witnessed an increased interest in characterizing the spontaneous brain activity occurring during resting conditions. The identification of connectivity patterns in this so-called "resting-state" has been the subject of a great number of electrophysiology-based studies, using the Electro/Magneto-Encephalography (EEG/MEG) source connectivity method. However, no consensus has been reached yet regarding a unified (if possible) analysis pipeline, and several involved parameters and methods require cautious tuning. This is particularly challenging when different analytical choices induce significant discrepancies in results and drawn conclusions, thereby hindering the reproducibility of neuroimaging research. Hence, our objective in this study was to shed light on the effect of analytical variability on outcome consistency by evaluating the implications of parameters involved in the EEG source connectivity analysis on the accuracy of resting-state networks (RSNs) reconstruction. We simulated, using neural mass models, EEG data corresponding to two RSNs, namely the default mode network (DMN) and dorsal attentional network (DAN). We investigated the impact of five channel densities (19, 32, 64, 128, 256), three inverse solutions (weighted minimum norm estimate (wMNE), exact low-resolution brain electromagnetic tomography (eLORETA), and linearly constrained minimum variance (LCMV) beamforming) and four functional connectivity measures (phase-locking value (PLV), phase-lag index (PLI), and amplitude envelope correlation (AEC) with and without source leakage correction), on the correspondence between reconstructed and reference networks. We showed that, with different analytical choices related to the number of electrodes, source reconstruction algorithm, and functional connectivity measure, high variability is present in the results. More specifically, our results show that a higher number of EEG channels significantly increased the accuracy of the reconstructed networks. Additionally, our results showed significant variability in the performance of the tested inverse solutions and connectivity measures. Such methodological variability and absence of analysis standardization represent a critical issue for neuroimaging studies that should be prioritized. We believe that this work could be useful for the field of electrophysiology connectomics, by increasing awareness regarding the challenge of variability in methodological approaches and its implications on reported results.

Not with a “zap” but with a “beep”: measuring the origins of perinatal experience

When does the mind begin? Infant psychology is mysterious in part because we cannot remember our first months of life, nor can we directly communicate with infants. Even more speculative is the possibility of mental life prior to birth. The question of when consciousness, or subjective experience, begins in human development thus remains incompletely answered, though boundaries can be set using current knowledge from developmental neurobiology and recent investigations of the perinatal brain. Here, we offer our perspective on how the development of a sensory perturbational complexity index (sPCI) based on auditory ("beep-and-zip"), visual ("flash-and-zip"), or even olfactory ("sniff-and-zip") cortical perturbations in place of electromagnetic perturbations ("zap-and-zip") might be used to address this question. First, we discuss recent studies of perinatal cognition and consciousness using techniques such as functional magnetic resonance imaging (fMRI), electroencephalography (EEG), and, in particular, magnetoencephalography (MEG). While newborn infants are the archetypal subjects for studying early human development, researchers may also benefit from fetal studies, as the womb is, in many respects, a more controlled environment than the cradle. The earliest possible timepoint when subjective experience might begin is likely the establishment of thalamocortical connectivity at 26 weeks gestation, as the thalamocortical system is necessary for consciousness according to most theoretical frameworks. To infer at what age and in which behavioral states consciousness might emerge following the initiation of thalamocortical pathways, we advocate for the development of the sPCI and similar techniques, based on EEG, MEG, and fMRI, to estimate the perinatal brain's state of consciousness.

Analysis of task-related MEG functional brain networks using dynamic mode decomposition

Objective.Functional connectivity networks explain the different brain states during the diverse motor, cognitive, and sensory functions. Extracting connectivity network configurations and their temporal evolution is crucial for understanding brain function during diverse behavioral tasks.Approach.In this study, we introduce the use of dynamic mode decomposition (DMD) to extract the dynamics of brain networks. We compared DMD with principal component analysis (PCA) using real magnetoencephalography data during motor and memory tasks.Main results.The framework generates dominant connectivity brain networks and their time dynamics during simple tasks, such as button press and left-hand movement, as well as more complex tasks, such as picture naming and memory tasks. Our findings show that the proposed methodology with both the PCA-based and DMD-based approaches extracts similar dominant connectivity networks and their corresponding temporal dynamics.Significance.We believe that the proposed methodology with both the PCA and the DMD approaches has a very high potential for deciphering the spatiotemporal dynamics of electrophysiological brain network states during tasks.

First few seconds for flow: A comprehensive proposal of the neurobiology and neurodynamics of state onset

Flow is a cognitive state that manifests when there is complete attentional absorption while performing a task. Flow occurs when certain internal as well as external conditions are present, including intense concentration, a sense of control, feedback, and a balance between the challenge of the task and the relevant skillset. Phenomenologically, flow is accompanied by a loss of self-consciousness, seamless integration of action and awareness, and acute changes in time perception. Research has begun to uncover some of the neurophysiological correlates of flow, as well as some of the state's neuromodulatory processes. We comprehensively review this work and consider the neurodynamics of the onset of the state, considering large-scale brain networks, as well as dopaminergic, noradrenergic, and endocannabinoid systems. To accomplish this, we outline an evidence-based hypothetical situation, and consider the flow state in a broader context including other profound alterations in consciousness, such as the psychedelic state and the state of traumatic stress that can induce PTSD. We present a broad theoretical framework which may motivate future testable hypotheses.

Cortico-cortical and thalamo-cortical connectivity during non-REM and REM sleep: Insights from intracranial recordings in humans

OBJECTIVE

To study changes of thalamo-cortical and cortico-cortical connectivity during wakefulness, non-Rapid Eye Movement (non-REM) sleep, including N2 and N3 stages, and REM sleep, using stereoelectroencephalography (SEEG) recording in humans.

METHODS

We studied SEEG recordings of ten patients during wakefulness, non-REM sleep and REM sleep, in seven brain regions of interest including the thalamus. We calculated directed and undirected functional connectivity using a measure of non-linear correlation coefficient h².

RESULTS

The thalamus was more connected to other brain regions during N2 stage and REM sleep than during N3 stage during which cortex was more connected than the thalamus. We found two significant directed links: the first from the prefrontal region to the lateral parietal region in the delta band during N3 sleep and the second from the thalamus to the insula during REM sleep.

CONCLUSIONS

These results showed that cortico-cortical connectivity is more prominent in N3 stage than in N2 and REM sleep. During REM sleep we found significant thalamo-insular connectivity, with a driving role of the thalamus.

SIGNIFICANCE

We found a pattern of cortical connectivity during N3 sleep concordant with antero-posterior traveling slow waves. The thalamus seemed particularly involved as a hub of connectivity during REM sleep.

Rate and oscillatory switching dynamics of a multilayer visual microcircuit model

The neocortex is organized around layered microcircuits consisting of a variety of excitatory and inhibitory neuronal types which perform rate- and oscillation-based computations. Using modeling, we show that both superficial and deep layers of the primary mouse visual cortex implement two ultrasensitive and bistable switches built on mutual inhibitory connectivity motives between somatostatin, parvalbumin, and vasoactive intestinal polypeptide cells. The switches toggle pyramidal neurons between high and low firing rate states that are synchronized across layers through translaminar connectivity. Moreover, inhibited and disinhibited states are characterized by low- and high-frequency oscillations, respectively, with layer-specific differences in frequency and power which show asymmetric changes during state transitions. These findings are consistent with a number of experimental observations and embed firing rate together with oscillatory changes within a switch interpretation of the microcircuit.

Methods Used to Estimate EEG Source-Space Networks: A Comparative Simulation-Based Study

Along with the study of the brain activity evoked by external stimuli, an important advance in current neuroscience involves understanding the spontaneous brain activity that occurs during resting conditions. Interestingly, the identification of the connectivity patterns in "resting-state" has been the subject of a great number of electrophysiology-based studies. In this context, the Electroencephalography (EEG) source connectivity method enables estimating resting-state cortical networks from scalp-EEG recordings. However, there is still no consensus over a unified pipeline adapted in all cases (e.g., type of task, a priori on studied networks) and numerous methodological questions remain unanswered. In order to address this problem, we simulated, using neural mass models, EEG data corresponding to the default mode network (DMN), the most widely studied resting-state network, and tested the effect of different channel densities, two inverse solutions and two functional connectivity measures on the correspondence between the reconstructed networks and the reference networks. Results showed that increasing the number of electrodes enhances the accuracy of the network reconstruction, and that eLORETA/PLV led to better accuracy than other inverse solution/connectivity measure combinations in terms of the correlation between reconstructed and reference connectivity matrices. This work has a wide range of implications in the field of electrophysiology connectomics, and is a step towards a convergence and standardization of approaches in this emerging field.

Biophysical Model: A Promising Method in the Study of the Mechanism of Propofol: A Narrative Review

The physiological and neuroregulatory mechanism of propofol is largely based on very limited knowledge. It is one of the important puzzling issues in anesthesiology and is of great value in both scientific and clinical fields. It is acknowledged that neural networks which are comprised of a number of neural circuits might be involved in the anesthetic mechanism. However, the mechanism of this hypothesis needs to be further elucidated. With the progress of artificial intelligence, it is more likely to solve this problem through using artificial neural networks to perform temporal waveform data analysis and to construct biophysical computational models. This review focuses on current knowledge regarding the anesthetic mechanism of propofol, an intravenous general anesthetic, by constructing biophysical computational models.

Whole-brain modelling identifies distinct but convergent paths to unconsciousness in anaesthesia and disorders of consciousness

The human brain entertains rich spatiotemporal dynamics, which are drastically reconfigured when consciousness is lost due to anaesthesia or disorders of consciousness (DOC). Here, we sought to identify the neurobiological mechanisms that explain how transient pharmacological intervention and chronic neuroanatomical injury can lead to common reconfigurations of neural activity. We developed and systematically perturbed a neurobiologically realistic model of whole-brain haemodynamic signals. By incorporating PET data about the cortical distribution of GABA receptors, our computational model reveals a key role of spatially-specific local inhibition for reproducing the functional MRI activity observed during anaesthesia with the GABA-ergic agent propofol. Additionally, incorporating diffusion MRI data obtained from DOC patients reveals that the dynamics that characterise loss of consciousness can also emerge from randomised neuroanatomical connectivity. Our results generalise between anaesthesia and DOC datasets, demonstrating how increased inhibition and connectome perturbation represent distinct neurobiological paths towards the characteristic activity of the unconscious brain.

Computational Modeling of Information Propagation during the Sleep–Waking Cycle

Simple Summary

During the deep phases of sleep we do not normally wake up by a thunder, but we nevertheless notice it when awake. The exact same sound gets to our ears and cortex through the thalamus and still, it triggers two very different responses. There is growing experimental evidence that these two states of the brain—sleep and wakefulness—distribute sensory information in different ways across the cortex. In particular, during sleep, neuronal responses remain local and do not spread out across distant synaptically connected regions. On the contrary, during wakefulness, stimuli are able to elicit a wider spatial response. We have used a computational model of coupled cortical columns to study how these two propagation modes arise. Moreover, the transition from sleep-like to waking-like dynamics occurs in agreement with the synaptic homeostasis hypothesis and only requires the increase of excitatory conductances. We have found that, in order to reproduce the aforementioned observations, this parameter change has to be selectively applied: synaptic conductances between distinct columns have to be potentiated over local ones.

Abstract

Non-threatening familiar sounds can go unnoticed during sleep despite the fact that they enter our brain by exciting the auditory nerves. Extracellular cortical recordings in the primary auditory cortex of rodents show that an increase in firing rate in response to pure tones during deep phases of sleep is comparable to those evoked during wakefulness. This result challenges the hypothesis that during sleep cortical responses are weakened through thalamic gating. An alternative explanation comes from the observation that the spatiotemporal spread of the evoked activity by transcranial magnetic stimulation in humans is reduced during non-rapid eye movement (NREM) sleep as compared to the wider propagation to other cortical regions during wakefulness. Thus, cortical responses during NREM sleep remain local and the stimulus only reaches nearby neuronal populations. We aim at understanding how this behavior emerges in the brain as it spontaneously shifts between NREM sleep and wakefulness. To do so, we have used a computational neural-mass model to reproduce the dynamics of the sensory auditory cortex and corresponding local field potentials in these two brain states. Following the synaptic homeostasis hypothesis, an increase in a single parameter, namely the excitatory conductance g¯AMPA, allows us to place the model from NREM sleep into wakefulness. In agreement with the experimental results, the endogenous dynamics during NREM sleep produces a comparable, even higher, response to excitatory inputs to the ones during wakefulness. We have extended the model to two bidirectionally connected cortical columns and have quantified the propagation of an excitatory input as a function of their coupling. We have found that the general increase in all conductances of the cortical excitatory synapses that drive the system from NREM sleep to wakefulness does not boost the effective connectivity between cortical columns. Instead, it is the inter-/intra-conductance ratio of cortical excitatory synapses that should raise to facilitate information propagation across the brain.

Canard solutions in neural mass models: consequences on critical regimes

Mathematical models at multiple temporal and spatial scales can unveil the fundamental mechanisms of critical transitions in brain activities. Neural mass models (NMMs) consider the average temporal dynamics of interconnected neuronal subpopulations without explicitly representing the underlying cellular activity. The mesoscopic level offered by the neural mass formulation has been used to model electroencephalographic (EEG) recordings and to investigate various cerebral mechanisms, such as the generation of physiological and pathological brain activities. In this work, we consider a NMM widely accepted in the context of epilepsy, which includes four interacting neuronal subpopulations with different synaptic kinetics. Due to the resulting three-time-scale structure, the model yields complex oscillations of relaxation and bursting types. By applying the principles of geometric singular perturbation theory, we unveil the existence of the canard solutions and detail how they organize the complex oscillations and excitability properties of the model. In particular, we show that boundaries between pathological epileptic discharges and physiological background activity are determined by the canard solutions. Finally we report the existence of canard-mediated small-amplitude frequency-specific oscillations in simulated local field potentials for decreased inhibition conditions. Interestingly, such oscillations are actually observed in intracerebral EEG signals recorded in epileptic patients during pre-ictal periods, close to seizure onsets.

Mean-Field Modeling of Brain-Scale Dynamics for the Evaluation of EEG Source-Space Networks

Understanding the dynamics of brain-scale functional networks at rest and during cognitive tasks is the subject of intense research efforts to unveil fundamental principles of brain functions. To estimate these large-scale brain networks, the emergent method called "electroencephalography (EEG) source connectivity" has generated increasing interest in the network neuroscience community, due to its ability to identify cortical brain networks with satisfactory spatio-temporal resolution, while reducing mixing and volume conduction effects. However, no consensus has been reached yet regarding a unified EEG source connectivity pipeline, and several methodological issues have to be carefully accounted to avoid pitfalls. Thus, a validation toolbox that provides flexible "ground truth" models is needed for an objective methods/parameters evaluation and, thereby an optimization of the EEG source connectivity pipeline. In this paper, we show how a recently developed large-scale model of brain-scale activity, named COALIA, can provide to some extent such ground truth by providing realistic simulations of source-level and scalp-level activity. Using a bottom-up approach, the model bridges cortical micro-circuitry and large-scale network dynamics. Here, we provide an example of the potential use of COALIA to analyze, in the context of epileptiform activity, the effect of three key factors involved in the "EEG source connectivity" pipeline: (i) EEG sensors density, (ii) algorithm used to solve the inverse problem, and (iii) functional connectivity measure. Results showed that a high electrode density (at least 64 channels) is required to accurately estimate cortical networks. Regarding the inverse solution/connectivity measure combination, the best performance at high electrode density was obtained using the weighted minimum norm estimate (wMNE) combined with the weighted phase lag index (wPLI). Although those results are specific to the considered aforementioned context (epileptiform activity), we believe that this model-based approach can be successfully applied to other experimental questions/contexts. We aim at presenting a proof-of-concept of the interest of COALIA in the network neuroscience field, and its potential use in optimizing the EEG source-space network estimation pipeline.

Mechanisms Underlying Disorders of Consciousness: Bridging Gaps to Move Toward an Integrated Translational Science

AIM

In order to successfully detect, classify, prognosticate, and develop targeted therapies for patients with disorders of consciousness (DOC), it is crucial to improve our mechanistic understanding of how severe brain injuries result in these disorders.

METHODS

To address this need, the Curing Coma Campaign convened a Mechanisms Sub-Group of the Coma Science Work Group (CSWG), aiming to identify the most pressing knowledge gaps and the most promising approaches to bridge them.

RESULTS

We identified a key conceptual gap in the need to differentiate the neural mechanisms of consciousness per se, from those underpinning connectedness to the environment and behavioral responsiveness. Further, we characterised three fundamental gaps in DOC research: (1) a lack of mechanistic integration between structural brain damage and abnormal brain function in DOC; (2) a lack of translational bridges between micro- and macro-scale neural phenomena; and (3) an incomplete exploration of possible synergies between data-driven and theory-driven approaches.

CONCLUSION

In this white paper, we discuss research priorities that would enable us to begin to close these knowledge gaps. We propose that a fundamental step towards this goal will be to combine translational, multi-scale, and multimodal data, with new biomarkers, theory-driven approaches, and computational models, to produce an integrated account of neural mechanisms in DOC. Importantly, we envision that reciprocal interaction between domains will establish a "virtuous cycle," leading towards a critical vantage point of integrated knowledge that will enable the advancement of the scientific understanding of DOC and consequently, an improvement of clinical practice.

Computational Models in Electroencephalography

Computational models lie at the intersection of basic neuroscience and healthcare applications because they allow researchers to test hypotheses in silico and predict the outcome of experiments and interactions that are very hard to test in reality. Yet, what is meant by “computational model” is understood in many different ways by researchers in different fields of neuroscience and psychology, hindering communication and collaboration. In this review, we point out the state of the art of computational modeling in Electroencephalography (EEG) and outline how these models can be used to integrate findings from electrophysiology, network-level models, and behavior. On the one hand, computational models serve to investigate the mechanisms that generate brain activity, for example measured with EEG, such as the transient emergence of oscillations at different frequency bands and/or with different spatial topographies. On the other hand, computational models serve to design experiments and test hypotheses in silico. The final purpose of computational models of EEG is to obtain a comprehensive understanding of the mechanisms that underlie the EEG signal. This is crucial for an accurate interpretation of EEG measurements that may ultimately serve in the development of novel clinical applications.

In silico model reveals the key role of GABA in KCNT1-epilepsy in infancy with migrating focal seizures

OBJECTIVE

This study was undertaken to investigate how gain of function (GOF) of slack channel due to a KCNT1 pathogenic variant induces abnormal neuronal cortical network activity and generates specific electroencephalographic (EEG) patterns of epilepsy in infancy with migrating focal seizures.

METHODS

We used detailed microscopic computational models of neurons to explore the impact of GOF of slack channel (explicitly coded) on each subtype of neurons and on a cortical micronetwork. Then, we adapted a thalamocortical macroscopic model considering results obtained in detailed models and immature properties related to epileptic brain in infancy. Finally, we compared simulated EEGs resulting from the macroscopic model with interictal and ictal patterns of affected individuals using our previously reported EEG markers.

RESULTS

The pathogenic variants of KCNT1 strongly decreased the firing rate properties of γ-aminobutyric acidergic (GABAergic) interneurons and, to a lesser extent, those of pyramidal cells. This change led to hyperexcitability with increased synchronization in a cortical micronetwork. At the macroscopic scale, introducing slack GOF effect resulted in epilepsy of infancy with migrating focal seizures (EIMFS) EEG interictal patterns. Increased excitation-to-inhibition ratio triggered seizure, but we had to add dynamic depolarizing GABA between somatostatin-positive interneurons and pyramidal cells to obtain migrating seizure. The simulated migrating seizures were close to EIMFS seizures, with similar values regarding the delay between the different ictal activities (one of the specific EEG markers of migrating focal seizures due to KCNT1 pathogenic variants).

SIGNIFICANCE

This study illustrates the interest of biomathematical models to explore pathophysiological mechanisms bridging the gap between the functional effect of gene pathogenic variants and specific EEG phenotype. Such models can be complementary to in vitro cellular and animal models. This multiscale approach provides an in silico framework that can be further used to identify candidate innovative therapies.

A Connectome-Based, Corticothalamic Model of State- and Stimulation-Dependent Modulation of Rhythmic Neural Activity and Connectivity

Rhythmic activity in the brain fluctuates with behaviour and cognitive state, through a combination of coexisting and interacting frequencies. At large spatial scales such as those studied in human M/EEG, measured oscillatory dynamics are believed to arise primarily from a combination of cortical (intracolumnar) and corticothalamic rhythmogenic mechanisms. Whilst considerable progress has been made in characterizing these two types of neural circuit separately, relatively little work has been done that attempts to unify them into a single consistent picture. This is the aim of the present paper. We present and examine a whole-brain, connectome-based neural mass model with detailed long-range cortico-cortical connectivity and strong, recurrent corticothalamic circuitry. This system reproduces a variety of known features of human M/EEG recordings, including spectral peaks at canonical frequencies, and functional connectivity structure that is shaped by the underlying anatomical connectivity. Importantly, our model is able to capture state- (e.g., idling/active) dependent fluctuations in oscillatory activity and the coexistence of multiple oscillatory phenomena, as well as frequency-specific modulation of functional connectivity. We find that increasing the level of sensory drive to the thalamus triggers a suppression of the dominant low frequency rhythms generated by corticothalamic loops, and subsequent disinhibition of higher frequency endogenous rhythmic behaviour of intracolumnar microcircuits. These combine to yield simultaneous decreases in lower frequency and increases in higher frequency components of the M/EEG power spectrum during states of high sensory or cognitive drive. Building on this, we also explored the effect of pulsatile brain stimulation on ongoing oscillatory activity, and evaluated the impact of coexistent frequencies and state-dependent fluctuations on the response of cortical networks. Our results provide new insight into the role played by cortical and corticothalamic circuits in shaping intrinsic brain rhythms, and suggest new directions for brain stimulation therapies aimed at state-and frequency-specific control of oscillatory brain activity.

Neural mass modeling of slow-fast dynamics of seizure initiation and abortion

Epilepsy is a dynamic and complex neurological disease affecting about 1% of the worldwide population, among which 30% of the patients are drug-resistant. Epilepsy is characterized by recurrent episodes of paroxysmal neural discharges (the so-called seizures), which manifest themselves through a large-amplitude rhythmic activity observed in depth-EEG recordings, in particular in local field potentials (LFPs). The signature characterizing the transition to seizures involves complex oscillatory patterns, which could serve as a marker to prevent seizure initiation by triggering appropriate therapeutic neurostimulation methods. To investigate such protocols, neurophysiological lumped-parameter models at the mesoscopic scale, namely neural mass models, are powerful tools that not only mimic the LFP signals but also give insights on the neural mechanisms related to different stages of seizures. Here, we analyze the multiple time-scale dynamics of a neural mass model and explain the underlying structure of the complex oscillations observed before seizure initiation. We investigate population-specific effects of the stimulation and the dependence of stimulation parameters on synaptic timescales. In particular, we show that intermediate stimulation frequencies (>20 Hz) can abort seizures if the timescale difference is pronounced. Those results have the potential in the design of therapeutic brain stimulation protocols based on the neurophysiological properties of tissue.

Revealing Spectrum Features of Stochastic Neuron Spike Trains

Power spectra of spike trains reveal important properties of neuronal behavior. They exhibit several peaks, whose shape and position depend on applied stimuli and intrinsic biophysical properties, such as input current density and channel noise. The position of the spectral peaks in the frequency domain is not straightforwardly predictable from statistical averages of the interspike intervals, especially when stochastic behavior prevails. In this work, we provide a model for the neuronal power spectrum, obtained from Discrete Fourier Transform and expressed as a series of expected value of sinusoidal terms. The first term of the series allows us to estimate the frequencies of the spectral peaks to a maximum error of a few Hz, and to interpret why they are not harmonics of the first peak frequency. Thus, the simple expression of the proposed power spectral density (PSD) model makes it a powerful interpretative tool of PSD shape, and also useful for neurophysiological studies aimed at extracting information on neuronal behavior from spike train spectra.