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.

Spontaneous fast-ultradian dynamics of polymorphic interictal events in drug-resistant focal epilepsy.

OBJECTIVE: We studied the rate dynamics of interictal events occurring over fast-ultradian time scales, as commonly examined in clinics to guide surgical planning in epilepsy. METHODS: Stereo-electroencephalography (SEEG) traces of 35 patients with good surgical outcome (Engel I) were analyzed. For this we developed a general data mining method aimed at clustering the plethora of transient waveform shapes including interictal epileptiform discharges (IEDs) and assessed the temporal fluctuations in the capability of mapping the epileptogenic zone (EZ) of each type of event. RESULTS: We found that the fast-ultradian dynamics of the IED rate may effectively impair the precision of EZ identification, and appear to occur spontaneously, that is, not triggered by or exclusively associated with a particular cognitive task, wakefulness, sleep, seizure occurrence, post-ictal state, or antiepileptic drug withdrawal. Propagation of IEDs from the EZ to the propagation zone (PZ) could explain the observed fast-ultradian fluctuations in a reduced fraction of the analyzed patients, suggesting that other factors like the excitability of the epileptogenic tissue could play a more relevant role. A novel link was found between the fast-ultradian dynamics of the overall rate of polymorphic events and the rate of specific IEDs subtypes. We exploited this feature to estimate in each patient the 5 min interictal epoch for near-optimal EZ and resected-zone (RZ) localization. This approach produces at the population level a better EZ/RZ classification when compared to both (1) the whole time series available in each patient (p = .084 for EZ, p < .001 for RZ, Wilcoxon signed-rank test) and (2) 5 min epochs sampled randomly from the interictal recordings of each patient (p < .05 for EZ, p < .001 for RZ, 105 random samplings). SIGNIFICANCE: Our results highlight the relevance of the fast-ultradian IED dynamics in mapping the EZ, and show how this dynamics can be estimated prospectively to inform surgical planning in epilepsy.

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.

Modelling acute and lasting effects of tDCS on epileptic activity.

Transcranial Direct brain stimulation (tDCS) is commonly used in order to modulate cortical networks activity during physiological processes through the application of weak electrical fields with scalp electrodes. Cathodal stimulation has been shown to decrease brain excitability in the context of epilepsy, with variable success. However, the cellular mechanisms responsible for the acute and the long-lasting effect of tDCS remain elusive. Using a novel approach of computational modeling that combines detailed but functionally integrated neurons we built a physiologically-based thalamocortical column. This model comprises 10,000 individual neurons made of pyramidal cells, and 3 types of gamma-aminobutyric acid (GABA) -ergic cells (VIP, PV, and SST) respecting the anatomy, layers, projection, connectivity and neurites orientation. Simulating realistic electric fields in term of intensity, main results showed that 1) tDCS effects are best explained by modulation of the presynaptic probability of release 2) tDCS affects the dynamic of cortical network only if a sufficient number of neurons are modulated 3)VIP GABAergic interneurons of the superficial layer of the cortex are especially affected by tDCS 4) Long lasting effect depends on glutamatergic synaptic plasticity.

Cortico-striatal synchronization in human focal seizures.

Although a number of experimental and clinical studies have pointed out participation or an even more prominent role of basal ganglia in focal seizures, the mode of interaction between cortical and striatal signals remains unclear. In the present study, we took stereoelectroencephalographic (SEEG) recordings in drug-resistant epilepsy patients, to qualitatively and quantitatively analyse the ictal striatum activity as well as its synchronization with cerebral cortex. Eleven patients who underwent SEEG evaluation were prospectively included if they fulfilled two inclusion criteria: (i) at least one orthogonal intracerebral electrode contact explored the basal ganglia, in either their putaminal or caudate part; and (ii) at least two SEEG seizures were recorded. Cortical and subcortical regions of interest were defined and different periods of interest were analysed. SEEG was visually inspected and h2 non-linear correlation analysis performed to study functional connectivity between cortical region of interest and striatum. Six correlation indices were calculated. Two main patterns of striatal activation were recorded: the most frequent was characterized by an early alpha/beta activity that started within the first 5 s after seizure onset, sometimes concomitant with it. The second one was characterized by late, slower, theta/delta activity. A significant difference in h2 correlation indices was observed during the preictal and seizure onset period compared to background for global striatal index, mesio-temporal/striatal index, latero-temporal/striatal index, insular/striatal index, prefrontal/striatal index. In addition, a significant difference in h2 correlation indices was observed during the seizure termination period compared to all the other periods of interest for the six indices calculated. These results indicate that cortico-striatal synchronization can arise from the start of focal seizures. Depending on the ictal frequency pattern, desynchronization can occur later, but a late and terminal hypersynchronization progressively takes over. These changes in synchronization level between cortical and striatal activity might be part of an endogenous mechanism controlling the duration of abnormal oscillations within the striato-thalamo-cortical loop and thereby their termination. Pathophysiology of basal ganglia in focal seizures appears to be much more interlinked with the cortex than expected. Beyond the stereotypical features they could imprint to seizure semiology, their role in strengthening mechanisms underlying cessation of ictal propagation should inspire new rationales for deep brain stimulation in patients with intractable focal epilepsies.

Exploring the Correlation Between M/EEG Source-Space and fMRI Networks at Rest.

Magneto/electro-encephalography (M/EEG) source connectivity is an emerging approach to estimate brain networks with high temporal and spatial resolutions. Here, we aim to evaluate the effect of functional connectivity (FC) methods on the correlation between M/EEG source-space and fMRI networks at rest. Two main FC families are tested: (i) FC methods that do not remove zero-lag connectivity including Phase Locking Value (PLV) and Amplitude Envelope Correlation (AEC) and (ii) FC methods that remove zero-lag connections such as Phase Lag Index (PLI) and two orthogonalisation approaches combined with PLV (PLVCol, PLVPas) and AEC (AECCol, AECPas). Methods are evaluated on resting state M/EEG signals recorded from healthy participants at rest (N = 74). Networks obtained by each FC method are compared with fMRI networks (obtained from the Human Connectome Project). Results show low correlations for all FC methods, however PLV and AEC networks are significantly correlated with fMRI networks (ρ = 0.12, p = 1.93 × 10-8 and ρ = 0.06, p = 0.007, respectively), while other methods are not. These observations are consistent for all M/EEG frequency bands and for different FC matrices threshold. Our main message is to be careful in selecting FC methods when comparing or combining M/EEG with fMRI. We consider that more comparative studies based on simulation and real data and at different levels (node, module or sub networks) are still needed in order to improve our understanding on the relationships between M/EEG source-space networks and fMRI networks at rest.

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.

Brain regions and epileptogenicity influence epileptic interictal spike production and propagation during NREM sleep in comparison with wakefulness.

OBJECTIVE: Non-rapid eye movement (NREM) sleep is known to be a brain state associated with an activation of interictal epileptic activity. The goal of this work was to quantify topographic changes occurring during NREM sleep in comparison with wakefulness. METHOD: We studied intracerebral recordings of 20 patients who underwent stereo-electroencephalography (SEEG) during presurgical evaluation for pharmacoresistant focal epilepsy. We measured the number of interictal spikes (IS) and quantified the co-occurrence of IS between brain regions during 1 hour of NREM sleep and 1 hour of wakefulness. Co-occurrence is a method to estimate IS networks based on a temporal concordance between IS of different brain regions. Each studied region was labeled as "seizure-onset zone" (SOZ), "propagation zone" (PZ), or "not involved region" (NIR). RESULTS: During NREM sleep, the number of interictal spikes significantly increased in all regions (mean of 68%). This increase was higher in medial temporal regions than in other regions, whether involved in the SOZ. Spike co-occurrence increased significantly in all regions during NREM sleep in comparison with wakefulness but was greater in neocortical regions. Spike co-occurrence in medial temporal regions was not higher than in other regions, suggesting that the increase of the number of spikes in this region was in great part a local effect. SIGNIFICANCE: This study demonstrated that medial temporal regions show a greater propensity to spike production or propagation during NREM sleep compared to other brain regions, even when the medial temporal lobe is not involved in the SOZ.

Defining epileptogenic networks: Contribution of SEEG and signal analysis.

Epileptogenic networks are defined by the brain regions involved in the production and propagation of epileptic activities. In this review we describe the historical, methodologic, and conceptual bases of this model in the analysis of electrophysiologic intracerebral recordings. In the context of epilepsy surgery, the determination of cerebral regions producing seizures (i.e., the "epileptogenic zone") is a crucial objective. In contrast with a traditional focal vision of focal drug-resistant epilepsies, the concept of epileptogenic networks has been progressively introduced as a model better able to describe the complexity of seizure dynamics and realistically describe the distribution of epileptogenic anomalies in the brain. The concept of epileptogenic networks is historically linked to the development of the stereoelectroencephalography (SEEG) method and subsequent introduction of means of quantifying the recorded signals. Seizures, and preictal and interictal discharges produce clear patterns on SEEG. These patterns can be analyzed utilizing signal analysis methods that quantify high-frequency oscillations or changes in functional connectivity. Dramatic changes in SEEG brain connectivity can be described during seizure genesis and propagation within cortical and subcortical regions, associated with the production of different patterns of seizure semiology. The interictal state is characterized by networks generating abnormal activities (interictal spikes) and also by modified functional properties. The introduction of novel approaches to large-scale modeling of these networks offers new methods in the goal of better predicting the effects of epilepsy surgery. The epileptogenic network concept is a key factor in identifying the anatomic distribution of the epileptogenic process, which is particularly important in the context of epilepsy surgery.

Update on the mechanisms and roles of high-frequency oscillations in seizures and epileptic disorders.

High-frequency oscillations (HFOs) are a type of brain activity that is recorded from brain regions capable of generating seizures. Because of the close association of HFOs with epileptogenic tissue and ictogenesis, understanding their cellular and network mechanisms could provide valuable information about the organization of epileptogenic networks and how seizures emerge from the abnormal activity of these networks. In this review, we summarize the most recent advances in the field of HFOs and provide a critical evaluation of new observations within the context of already established knowledge. Recent improvements in recording technology and the introduction of optogenetics into epilepsy research have intensified experimental work on HFOs. Using advanced computer models, new cellular substrates of epileptic HFOs were identified and the role of specific neuronal subtypes in HFO genesis was determined. Traditionally, the pathogenesis of HFOs was explored mainly in patients with temporal lobe epilepsy and in animal models mimicking this condition. HFOs have also been reported to occur in other epileptic disorders and models such as neocortical epilepsy, genetically determined epilepsies, and infantile spasms, which further support the significance of HFOs in the pathophysiology of epilepsy. It is increasingly recognized that HFOs are generated by multiple mechanisms at both the cellular and network levels. Future studies on HFOs combining novel high-resolution in vivo imaging techniques and precise control of neuronal behavior using optogenetics or chemogenetics will provide evidence about the causal role of HFOs in seizures and epileptogenesis. Detailed understanding of the pathophysiology of HFOs will propel better HFO classification and increase their information yield for clinical and diagnostic purposes.

Identification of Interictal Epileptic Networks from Dense-EEG.

Epilepsy is a network disease. The epileptic network usually involves spatially distributed brain regions. In this context, noninvasive M/EEG source connectivity is an emerging technique to identify functional brain networks at cortical level from noninvasive recordings. In this paper, we analyze the effect of the two key factors involved in EEG source connectivity processing: (i) the algorithm used in the solution of the EEG inverse problem and (ii) the method used in the estimation of the functional connectivity. We evaluate four inverse solutions algorithms (dSPM, wMNE, sLORETA and cMEM) and four connectivity measures (r 2, h 2, PLV, and MI) on data simulated from a combined biophysical/physiological model to generate realistic interictal epileptic spikes reflected in scalp EEG. We use a new network-based similarity index to compare between the network identified by each of the inverse/connectivity combination and the original network generated in the model. The method will be also applied on real data recorded from one epileptic patient who underwent a full presurgical evaluation for drug-resistant focal epilepsy. In simulated data, results revealed that the selection of the inverse/connectivity combination has a significant impact on the identified networks. Results suggested that nonlinear methods (nonlinear correlation coefficient, phase synchronization and mutual information) for measuring the connectivity are more efficient than the linear one (the cross correlation coefficient). The wMNE inverse solution showed higher performance than dSPM, cMEM and sLORETA. In real data, the combination (wMNE/PLV) led to a very good matching between the interictal epileptic network identified from noninvasive EEG recordings and the network obtained from connectivity analysis of intracerebral EEG recordings. These results suggest that source connectivity method, when appropriately configured, is able to extract highly relevant diagnostic information about networks involved in interictal epileptic spikes from non-invasive dense-EEG data.

Simultaneous Intracranial EEG-fMRI Shows Inter-Modality Correlation in Time-Resolved Connectivity Within Normal Areas but Not Within Epileptic Regions.

For the first time in research in humans, we used simultaneous icEEG-fMRI to examine the link between connectivity in haemodynamic signals during the resting-state (rs) and connectivity derived from electrophysiological activity in terms of the inter-modal connectivity correlation (IMCC). We quantified IMCC in nine patients with drug-resistant epilepsy (i) within brain networks in 'healthy' non-involved cortical zones (NIZ) and (ii) within brain networks involved in generating seizures and interictal spikes (IZ1) or solely spikes (IZ2). Functional connectivity (h 2 ) estimates for 10 min of resting-state data were obtained between each pair of electrodes within each clinical zone for both icEEG and fMRI. A sliding window approach allowed us to quantify the variability over time of h 2 (vh 2) as an indicator of connectivity dynamics. We observe significant positive IMCC for h 2 and vh 2, for multiple bands in the NIZ only, with the strongest effect in the lower icEEG frequencies. Similarly, intra-modal h 2 and vh 2 were found to be differently modified as a function of different epileptic processes: compared to NIZ, [Formula: see text] was higher in IZ1, but lower in IZ2, while [Formula: see text] showed the inverse pattern. This corroborates previous observations of inter-modal connectivity discrepancies in pathological cortices, while providing the first direct invasive and simultaneous comparison in humans. We also studied time-resolved FC variability multimodally for the first time, finding that IZ1 shows both elevated internal [Formula: see text] and less rich dynamical variability, suggesting that its chronic role in epileptogenesis may be linked to greater homogeneity in self-sustaining pathological oscillatory states.

Epilepsy in young Tsc1(+/-) mice exhibits age-dependent expression that mimics that of human tuberous sclerosis complex.

OBJECTIVE: To describe the epileptic phenotype of Tsc1(+/-) mice pups in comparison with age-related seizures in human tuberous sclerosis complex (TSC). METHODS: Tsc1(+/-) and control mice underwent intracranial electroencephalography (EEG) recording at postnatal ages (P)8 to P33, with linear silicon probe implanted in the somatosensory cortex of one or both hemispheres for 8-24 h. Ictal events were classified visually by independent analyzers; distinct EEG patterns were related to age and analyzed to quantify field potential characteristics and signal dynamics between hemispheres. We collected retrospectively 20 infants with prenatally diagnosed TSC and EEG before seizure onset, and analyzed the electroclinical course of epilepsy, taking into account a first-line treatment by vigabatrin. RESULTS: Spontaneous seizures were disclosed in 55% of Tsc1(+/-) mice at P9-18. Three ictal patterns were identified: from P9 to P12 "spike clusters" consisted of recurring large spikes without clinical correlate; "spasm-like" discharges dominated from P13 to P16 consisting of high amplitude large field potential superimposed with or followed by fast activity repeated every 2-10 s for at least 20 s, accompanied by rhythmic limb contractions; from P14 to P18 a "tonic-clonic like" pattern comprised rhythmic spikes of increasing amplitude with tonic-clonic movements. Early onset "spike clusters" were mainly unilateral, whereas "spasm-like" and "tonic-clonic like" patterns were bilateral. Interhemispheric propagation was significantly faster for "tonic-clonic like" than for "spasm-like" events. In infants diagnosed prenatally with TSC, clusters of sharp waves or spikes preceded the first seizure, and vigabatrin prevented the development of seizures. Patients treated after seizure onset developed spasms or focal seizures that were pharmacoresistant in 66.7% of cases. SIGNIFICANCE: Tsc1(+/-) mice pups exhibit an age-dependent seizure pattern sequence mimicking early human TSC epilepsy features. Spike clusters before seizure onset in TSC should be considered as a first stage of epilepsy reinforcing the concept of preventive antiepileptic therapy.

Localization of Epileptogenic Zone on Pre-surgical Intracranial EEG Recordings: Toward a Validation of Quantitative Signal Analysis Approaches.

In patients diagnosed with pharmaco-resistant epilepsy, cerebral areas responsible for seizure generation can be defined by performing implantation of intracranial electrodes. The identification of the epileptogenic zone (EZ) is based on visual inspection of the intracranial electroencephalogram (IEEG) performed by highly qualified neurophysiologists. New computer-based quantitative EEG analyses have been developed in collaboration with the signal analysis community to expedite EZ detection. The aim of the present report is to compare different signal analysis approaches developed in four different European laboratories working in close collaboration with four European Epilepsy Centers. Computer-based signal analysis methods were retrospectively applied to IEEG recordings performed in four patients undergoing pre-surgical exploration of pharmaco-resistant epilepsy. The four methods elaborated by the different teams to identify the EZ are based either on frequency analysis, on nonlinear signal analysis, on connectivity measures or on statistical parametric mapping of epileptogenicity indices. All methods converge on the identification of EZ in patients that present with fast activity at seizure onset. When traditional visual inspection was not successful in detecting EZ on IEEG, the different signal analysis methods produced highly discordant results. Quantitative analysis of IEEG recordings complement clinical evaluation by contributing to the study of epileptogenic networks during seizures. We demonstrate that the degree of sensitivity of different computer-based methods to detect the EZ in respect to visual EEG inspection depends on the specific seizure pattern.

Transcranial direct-current stimulation modulates synaptic mechanisms involved in associative learning in behaving rabbits.

Transcranial direct-current stimulation (tDCS) is a noninvasive brain stimulation technique that has been successfully applied for modulation of cortical excitability. tDCS is capable of inducing changes in neuronal membrane potentials in a polarity-dependent manner. When tDCS is of sufficient length, synaptically driven after-effects are induced. The mechanisms underlying these after-effects are largely unknown, and there is a compelling need for animal models to test the immediate effects and after-effects induced by tDCS in different cortical areas and evaluate the implications in complex cerebral processes. Here we show in behaving rabbits that tDCS applied over the somatosensory cortex modulates cortical processes consequent to localized stimulation of the whisker pad or of the corresponding area of the ventroposterior medial (VPM) thalamic nucleus. With longer stimulation periods, poststimulation effects were observed in the somatosensory cortex only after cathodal tDCS. Consistent with the polarity-specific effects, the acquisition of classical eyeblink conditioning was potentiated or depressed by the simultaneous application of anodal or cathodal tDCS, respectively, when stimulation of the whisker pad was used as conditioned stimulus, suggesting that tDCS modulates the sensory perception process necessary for associative learning. We also studied the putative mechanisms underlying immediate effects and after-effects of tDCS observed in the somatosensory cortex. Results when pairs of pulses applied to the thalamic VPM nucleus (mediating sensory input) during anodal and cathodal tDCS suggest that tDCS modifies thalamocortical synapses at presynaptic sites. Finally, we show that blocking the activation of adenosine A1 receptors prevents the long-term depression (LTD) evoked in the somatosensory cortex after cathodal tDCS.

Tuning Microelectrodes' Impedance to Improve Fast Ripples Recording.

Epilepsy is a chronic neurological disorder characterized by recurrent seizures resulting from abnormal neuronal hyperexcitability. In the case of pharmacoresistant epilepsy requiring resection surgery, the identification of the Epileptogenic Zone (EZ) is critical. Fast Ripples (FRs; 200-600 Hz) are one of the promising biomarkers that can aid in EZ delineation. However, recording FRs requires physically small electrodes. These microelectrodes suffer from high impedance, which significantly impacts FRs' observability and detection. In this study, we investigated the potential of a conductive polymer coating to enhance FR observability. We employed biophysical modeling to compare two types of microelectrodes: Gold (Au) and Au coated with the conductive polymer poly(3,4-ethylenedioxythiophene)-poly(styrene sulfonate) (Au/PEDOT:PSS). These electrodes were then implanted into the CA1 hippocampal neural network of epileptic mice to record FRs during epileptogenesis. The results showed that the polymer-coated electrodes had a two-order lower impedance as well as a higher transfer function amplitude and cut-off frequency. Consequently, FRs recorded with the PEDOT:PSS-coated microelectrode yielded significantly higher signal energy compared to the uncoated one. The PEDOT:PSS coating improved the observability of the recorded FRs and thus their detection. This work paves the way for the development of signal-specific microelectrode designs that allow for better targeting of pathological biomarkers.

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.