Modification of brain conductivity in human focal epilepsy: A model-based estimation from stereoelectroencephalography.

OBJECTIVE: We have developed a novel method for estimating brain tissue electrical conductivity using low-intensity pulse stereoelectroencephalography (SEEG) stimulation coupled with biophysical modeling. We evaluated the hypothesis that brain conductivity is correlated with the degree of epileptogenicity in patients with drug-resistant focal epilepsy. METHODS: We used bipolar low-intensity biphasic pulse stimulation (.2 mA) followed by a postprocessing pipeline for estimating brain conductivity. This processing is based on biophysical modeling of the electrical potential induced in brain tissue between the stimulated contacts in response to pulse stimulation. We estimated the degree of epileptogenicity using a semi-automatic method quantifying the dynamic of fast discharge at seizure onset: the epileptogenicity index (EI). We also investigated how the location of stimulation within specific anatomical brain regions or within lesional tissue impacts brain conductivity. RESULTS: We performed 1034 stimulations of 511 bipolar channels in 16 patients. We found that brain conductivity was lower in the epileptogenic zone (EZ; unpaired median difference = .064, p < .001) and inversely correlated with the epileptogenic index value (p < .001, Spearman rho = -.32). Conductivity values were also influenced by anatomical site, location within lesion, and delay between SEEG electrode implantation and stimulation, and had significant interpatient variability. Mixed model multivariate analysis showed that conductivity is significantly associated with EI (F = 13.45, p < .001), anatomical regions (F = 5.586, p < .001), delay since implantation (F = 14.71, p = .003), and age at SEEG (F = 6.591, p = .027), but not with the type of lesion (F = .372, p = .773) or the delay since last seizure (F = 1.592, p = .235). SIGNIFICANCE: We provide a novel model-based method for estimating brain conductivity from SEEG low-intensity pulse stimulations. The brain tissue conductivity is lower in EZ as compared to non-EZ. Conductivity also varies significantly across anatomical brain regions. Involved pathophysiological processes may include changes in the extracellular space (especially volume or tortuosity) in epileptic tissue.

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.

Spatio-temporal dynamics of large-scale electrophysiological networks during cognitive action control in healthy controls and Parkinson's disease patients.

Among the cognitive symptoms that are associated with Parkinson's disease (PD), alterations in cognitive action control (CAC) are commonly reported in patients. CAC enables the suppression of an automatic action, in favor of a goal-directed one. The implementation of CAC is time-resolved and arguably associated with dynamic changes in functional brain networks. However, the electrophysiological functional networks involved, their dynamic changes, and how these changes are affected by PD, still remain unknown. In this study, to address this gap of knowledge, 10 PD patients and 10 healthy controls (HC) underwent a Simon task while high-density electroencephalography (HD-EEG) was recorded. Source-level dynamic connectivity matrices were estimated using the phase-locking value in the beta (12-25 Hz) and gamma (30-45 Hz) frequency bands. Temporal independent component analyses were used as a dimension reduction tool to isolate the task-related brain network states. Typical microstate metrics were quantified to investigate the presence of these states at the subject-level. Our results first confirmed that PD patients experienced difficulties in inhibiting automatic responses during the task. At the group-level, we found three functional network states in the beta band that involved fronto-temporal, temporo-cingulate and fronto-frontal connections with typical CAC-related prefrontal and cingulate nodes (e.g., inferior frontal cortex). The presence of these networks did not differ between PD patients and HC when analyzing microstates metrics, and no robust correlations with behavior were found. In the gamma band, five networks were found, including one fronto-temporal network that was identical to the one found in the beta band. These networks also included CAC-related nodes previously identified in different neuroimaging modalities. Similarly to the beta networks, no subject-level differences were found between PD patients and HC. Interestingly, in both frequency bands, the dominant network at the subject-level was never the one that was the most durably modulated by the task. Altogether, this study identified the dynamic functional brain networks observed during CAC, but did not highlight PD-related changes in these networks that might explain behavioral changes. Although other new methods might be needed to investigate the presence of task-related networks at the subject-level, this study still highlights that task-based dynamic functional connectivity is a promising approach in understanding the cognitive dysfunctions observed in PD and beyond.

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.

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.

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.

Subthalamic nucleus local field potentials recordings reveal subtle effects of promised reward during conflict resolution in Parkinson's disease.

Cognitive action control depends on cortical-subcortical circuits, involving notably the subthalamic nucleus (STN), as evidenced by local field potentials recordings (LFPs) studies. The STN consistently shows an increase in theta oscillations power during conflict resolution. Some studies have shown that cognitive action control in Parkinson's disease (PD) could be influenced by the occurrence of monetary reward. In this study, we investigated whether incentive motivation could modulate STN activity, and notably STN theta activity, during response conflict resolution. To achieve this objective, we recorded STN LFPs during a motivated Simon task in PD patients who had undergone deep brain stimulation surgery. Behavioral results revealed that promised rewards increased the difficulty in resolving conflict situations, thus replicating previous findings. Signal analyses locked on the imperative stimulus onset revealed the typical pattern of increased theta power in a conflict situation. However, this conflict-related modulation of theta power was not influenced by the size of the reward cued. We nonetheless identified a significant effect of the reward size on local functional organization (indexed by inter-trial phase clustering) of theta oscillations, with higher organization associated with high rewards while resolving conflict. When focusing on the period following the onset of the reward cue, we unveiled a stronger beta power decrease in higher reward conditions. However, these LFPs results were not correlated to behavioral results. Our study suggests that the STN is involved in how reward information can influence computations during conflict resolution. However, considering recent studies as well as the present results, we suspect that these effects are subtle.

Quantitative analysis and EEG markers of KCNT1 epilepsy of infancy with migrating focal seizures.

OBJECTIVE: We aimed to characterize epilepsy of infancy with migrating focal seizures (EIMFS), a rare, severe early onset developmental epilepsy related to KCNT1 mutation, and to define specific electroencephalography (EEG) markers using EEG quantitative analysis. The ultimate goal would be to improve early diagnosis and to better understand seizure onset and propagation of EIMFS as compared to other early onset developmental epilepsy. METHODS: EEG of 7 EIMFS patients with KCNT1 mutations (115 seizures) and 17 patients with other early onset epilepsies (30 seizures) was included in this study. After detection of seizure onset and termination, spatiotemporal characteristics were quantified. Seizure propagation dynamics were analyzed using chronograms and phase coherence. RESULTS: In patients with EIMFS, seizures started and were localized predominantly in temporal and occipital areas, and evolved with a stable frequency (4-10 Hz). Inter- and intrahemispheric migrations were present in 60% of EIMFS seizures with high intraindividual reproducibility of temporospatial dynamics. Interhemispheric migrating seizures spread in 71% from temporal or occipital channels to the homologous contralateral ones, whereas intrahemispheric seizures involved mainly frontotemporal, temporal, and occipital channels. Causality links were present between ictal activities detected under different channels during migrating seizures. Finally, time delay index (based on delays between the different ictal onsets) and phase correlation index (based on coherence of ictal activities) allowed discrimination of EIMFS and non-EIMFS seizures with a specificity of 91.2% and a sensitivity of 84.4%. SIGNIFICANCE: We showed that the migrating pattern in EIMFS is not a random process, as suggested previously, and that it is a particular propagation pattern that follows the classical propagation pathways. It is notable that this study reveals specific EEG markers (time delay and phase correlation) accessible to visual evaluation, which will improve EIMFS diagnosis.

Human exposure to power frequency magnetic fields up to 7.6 mT: An integrated EEG/fMRI study.

We assessed the effects of power-line frequency (60 Hz in North America) magnetic fields (MF) in humans using simultaneous electroencephalography (EEG) and functional magnetic resonance imaging (fMRI). Twenty-five participants were enrolled in a pseudo-double-blind experiment involving "real" or "sham" exposure to sinusoidal 60 Hz MF exposures delivered using the gradient coil of an MRI scanner following two conditions: (i) 10 s exposures at 3 mT (10 repetitions); (ii) 2 s exposures at 7.6 mT (100 repetitions). Occipital EEG spectral power was computed in the alpha range (8-12 Hz, reportedly the most sensitive to MF exposure in the literature) with/without exposure. Brain functional activation was studied using fMRI blood oxygen level-dependent (BOLD, inversely correlated with EEG alpha power) maps. No significant effects were detected on occipital EEG alpha power during or post-exposure for any exposure condition. Consistent with EEG results, no effects were observed on fMRI BOLD maps in any brain region. Our results suggest that acute exposure (2-10 s) to 60 Hz MF from 3 to 7.6 mT (30,000 to 76,000 times higher than average public exposure levels for 60 Hz MF) does not induce detectable changes in EEG or BOLD signals. Combined with previous findings in which effects were observed on the BOLD signal after 1 h exposure to 3 mT, 60 Hz MF, this suggests that MF exposure in the low mT range (<10 mT) might require prolonged durations of exposure to induce detectable effects. Bioelectromagnetics. 38:425-435, 2017. © 2017 Wiley Periodicals, Inc.

Thresholds and mechanisms of human magnetophosphene perception induced by low frequency sinusoidal magnetic fields.

BACKGROUND: Virtually everyone is exposed to power-frequency MF (50/60 Hz), inducing in our body electric fields and currents, potentially modulating brain function. MF-induced electric fields within the central nervous system can generate flickering visual perceptions (magnetophosphenes), which form the basis of international MF exposure guidelines and recommendations protecting workers and the general public. However, magnetophosphene perception thresholds were estimated 40 years ago in a small, unreplicated study with significant uncertainties and leaving open the question of the involved interaction site. METHODS: We used a stimulation modality termed transcranial alternating magnetic stimulation (tAMS), delivering in situ sinusoidal electric fields comparable to transcranial alternating current stimulation (tACS). Magnetophosphene perception was quantified in 81 volunteers exposed to MF (eye or occipital exposure) between 0 and 50 mT at frequencies of 20, 50, 60 and 100 Hz. RESULTS: Reliable magnetophosphene perception was induced with tAMS without any scalp sensation, a major advantage as compared to tACS. Frequency-dependent thresholds were quantified using binary logistic regressions hence allowing to establish condition dependent probabilities of perception. Results support an interaction between induced current density and retinal rod cells. CONCLUSION: Beyond fundamental and immediate implications for international safety guidelines, and for identifying the interaction site underlying phosphene perception (ubiquitous in tACS experiments), our results support exploring the potential of tAMS for the differential diagnosis of retinal disorders and neuromodulation therapy.

Possible mechanisms of synaptic plasticity modulation by extremely low-frequency magnetic fields.

Understanding the biological mechanisms by which extremely low-frequency (ELF, < 300 Hz) magnetic fields (MFs) interact with human brain activity is an active field of research. Such knowledge is required by international agencies providing guidelines for general public and workers exposure to ELF MFs (such as ICNIRP, the International Commission on Non-Ionizing Radiation Protection). The identification of these interaction mechanisms is extremely challenging, since the effects of ELF MF exposure need to be monitored and understood at very different spatial (from micrometers to centimeters) and temporal (from milliseconds to minutes) scales. One possibility to overcome these issues is to develop biophysical models, based on the systems of mathematical equations describing the electric or metabolic activity of the brain tissue. Biophysical models of the brain activity offer the possibility to simulate how the brain tissue interacts with ELF MFs, in order to gain new insights into experimental data, and to test novel hypotheses regarding interaction mechanisms. This paper presents novel hypotheses regarding the effects of power line (60 Hz in North America) MFs on human brain activity, with arguments from biophysical models. We suggest a hypothetic chain of events that could bridge MF exposure with detectable effects on human neurophysiology. We also suggest novel directions of research in order to reach a convergence of biophysical models of brain activity and corresponding experimental data to identify interaction mechanisms.

Quasi-static approximation error of electric field analysis for transcranial current stimulation.

Objective.Numerical modeling of electric fields induced by transcranial alternating current stimulation (tACS) is currently a part of the standard procedure to predict and understand neural response. Quasi-static approximation (QSA) for electric field calculations is generally applied to reduce the computational cost. Here, we aimed to analyze and quantify the validity of the approximation over a broad frequency range.Approach.We performed electromagnetic modeling studies using an anatomical head model and considered approximations assuming either a purely ohmic medium (i.e. static formulation) or a lossy dielectric medium (QS formulation). The results were compared with the solution of Maxwell's equations in the cases of harmonic and pulsed signals. Finally, we analyzed the effect of electrode positioning on these errors.Main results.Our findings demonstrate that the QSA is valid and produces a relative error below 1% up to 1.43 MHz. The largest error is introduced in the static case, where the error is over 1% across the entire considered spectrum and as high as 20% in the brain at 10 Hz. We also highlight the special importance of considering the capacitive effect of tissues for pulsed waveforms, which prevents signal distortion induced by the purely ohmic approximation. At the neuron level, the results point a difference of sense electric field as high as 22% at focusing point, impacting pyramidal cells firing times.Significance.QSA remains valid in the frequency range currently used for tACS. However, neglecting permittivity (static formulation) introduces significant error for both harmonic and non-harmonic signals. It points out that reliable low frequency dielectric data are needed for accurate transcranial current stimulation numerical modeling.

Biophysical modeling of the electric field magnitude and distribution induced by electrical stimulation with intracerebral electrodes.

Intracranial electrodes are used clinically for diagnostic or therapeutic purposes, notably in drug-refractory epilepsy (DRE) among others. Visualization and quantification of the energy delivered through such electrodes is key to understanding how the resulting electric fields modulate neuronal excitability, i.e. the ratio between excitation and inhibition. Quantifying the electric field induced by electrical stimulation in a patient-specific manner is challenging, because these electric fields depend on a number of factors: electrode trajectory with respect to folded brain anatomy, biophysical (electrical conductivity / permittivity) properties of brain tissue and stimulation parameters such as electrode contacts position and intensity. Here, we aimed to evaluate various biophysical models for characterizing the electric fields induced by electrical stimulation in DRE patients undergoing stereoelectroencephalography (SEEG) recordings in the context of pre-surgical evaluation. This stimulation was performed with multiple-contact intracranial electrodes used in routine clinical practice. We introduced realistic 3D models of electrode geometry and trajectory in the neocortex. For the electrodes, we compared point (0D) and line (1D) sources approximations. For brain tissue, we considered three configurations of increasing complexity: a 6-layer spherical model, a toy model with a sulcus representation, replicating results from previous approaches; and went beyond the state-of-the-art by using a realistic head model geometry. Electrode geometry influenced the electric field distribution at close distances (∼3 mm) from the electrode axis. For larger distances, the volume conductor geometry and electrical conductivity dominated electric field distribution. These results are the first step towards accurate and computationally tractable patient-specific models of electric fields induced by neuromodulation and neurostimulation procedures.

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.

A personalizable autonomous neural mass model of epileptic seizures.

Work in the last two decades has shown that neural mass models (NMM) can realistically reproduce and explain epileptic seizure transitions as recorded by electrophysiological methods (EEG, SEEG). In previous work, advances were achieved by increasing excitation and heuristically varying network inhibitory coupling parameters in the models. Based on these early studies, we provide a laminar NMM capable of realistically reproducing the electrical activity recorded by SEEG in the epileptogenic zone during interictal to ictal states. With the exception of the external noise input into the pyramidal cell population, the model dynamics are autonomous. By setting the system at a point close to bifurcation, seizure-like transitions are generated, including pre-ictal spikes, low voltage fast activity, and ictal rhythmic activity. A novel element in the model is a physiologically motivated algorithm for chloride dynamics: the gain of GABAergic post-synaptic potentials is modulated by the pathological accumulation of chloride in pyramidal cells due to high inhibitory input and/or dysfunctional chloride transport. In addition, in order to simulate SEEG signals for comparison with real seizure recordings, the NMM is embedded first in a layered model of the neocortex and then in a realistic physical model. We compare modeling results with data from four epilepsy patient cases. By including key pathophysiological mechanisms, the proposed framework captures succinctly the electrophysiological phenomenology observed in ictal states, paving the way for robust personalization methods based on NMMs.

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.

Design of optimal multi-site brain stimulation protocols via neuro-inspired epilepsy models for abatement of interictal discharges.

Objective. Electrical brain stimulation is recognized as a promising therapeutic approach for treating brain disorders such as epilepsy. However, the use of this technique is still largely empirical, since stimulation parameters and targets are chosen using a trial-and-error approach. Therefore, there is a pressing need to design optimal, rationale-based multi-site brain stimulation protocols to control epileptiform activity.Approach. Here, we developed biologically-inspired models of brain activity receiving stimulation at two levels of description (single- and multi-population epileptogenic networks). First, we used bifurcation analysis to determine optimal parameters able to abort epileptiform patterns. Second, we present a graph-theory based method to classify network populations in an epileptogenic network based on their contribution to seizure generation and propagation. Main results. The best therapeutic effects (i.e. reduction of epileptiform discharges duration and occurrence rate) were obtained by the specific targeting of populations with the highest eigenvector centrality values. The timing of stimulation was also found to be critical in seizure abortion impact.Significance. Overall, our results provide a proof-of-concept that using network neuroscience combined with physiology-based computational models of brain activity can provide an effective method for the rational design of brain stimulation protocols in epilepsy.

Decoding the circuitry of consciousness: From local microcircuits to brain-scale networks.

Identifying the physiological processes underlying the emergence and maintenance of consciousness is one of the most fundamental problems of neuroscience, with implications ranging from fundamental neuroscience to the treatment of patients with disorders of consciousness (DOCs). One major challenge is to understand how cortical circuits at drastically different spatial scales, from local networks to brain-scale networks, operate in concert to enable consciousness, and how those processes are impaired in DOC patients. In this review, we attempt to relate available neurophysiological and clinical data with existing theoretical models of consciousness, while linking the micro- and macrocircuit levels. First, we address the relationships between awareness and wakefulness on the one hand, and cortico-cortical and thalamo-cortical connectivity on the other hand. Second, we discuss the role of three main types of GABAergic interneurons in specific circuits responsible for the dynamical reorganization of functional networks. Third, we explore advances in the functional role of nested oscillations for neural synchronization and communication, emphasizing the importance of the balance between local (high-frequency) and distant (low-frequency) activity for efficient information processing. The clinical implications of these theoretical considerations are presented. We propose that such cellular-scale mechanisms could extend current theories of consciousness.

Probing the circuits of conscious perception with magnetophosphenes.

OBJECTIVE: We aimed at characterizing, in non-invasive human brain recordings, the large-scale, coordinated activation of distant brain regions thought to occur during conscious perception. This process is termed ignition in the Global Workspace Theory, and integration in Integrated Information Theory, which are two of the major theories of consciousness. APPROACH: Here, we provide evidence for this process in humans by combining a magnetically-induced phosphene perception task with electroencephalography. Functional cortical networks were identified and characterized using graph theory to quantify the impact of conscious perception on local (segregation) and distant (integration) processing. MAIN RESULTS: Conscious phosphene perception activated frequency-specific networks, each associated with a specific spatial scale of information processing. Integration increased within an alpha-band functional network, while changes in segregation occurred in the beta band. SIGNIFICANCE: These results bring novel evidence for the functional role of distinct brain oscillations and confirm the key role of integration processes for conscious perception in humans.

Brain Network Dynamics Correlate with Personality Traits.

Identifying the neural substrates underlying the personality traits is a topic of great interest. On the other hand, it is now established that the brain is a dynamic networked system that can be studied by using functional connectivity techniques. However, much of the current understanding of personality-related differences in functional connectivity has been obtained through the stationary analysis, which does not capture the complex dynamical properties of brain networks. In this study, we aimed at evaluating the feasibility of using dynamic network measures to predict personality traits. Using the electro-encephalography (EEG)/magneto-encephalography (MEG) source connectivity method combined with a sliding window approach, dynamic functional brain networks were reconstructed from two datasets: (1) resting-state EEG data acquired from 56 subjects; (2) resting-state MEG data provided from the Human Connectome Project. Then, several dynamic functional connectivity metrics were evaluated. Similar observations were obtained by the two modalities (EEG and MEG) according to the neuroticism, which showed a negative correlation with the dynamic variability of resting-state brain networks. In particular, a significant relationship between this personality trait and the dynamic variability of the temporal lobe regions was observed. Results also revealed that extraversion and openness are positively correlated with the dynamics of the brain networks. These findings highlight the importance of tracking the dynamics of functional brain networks to improve our understanding about the neural substrates of personality.