Propagation of epileptiform activity can be independent of synaptic transmission, gap junctions, or diffusion and is consistent with electrical field transmission

Does the brain behave like a (complex) network? I. Dynamics

Graph theory is now becoming a standard tool in system-level neuroscience. However, endowing observed brain anatomy and dynamics with a complex network structure does not entail that the brain actually works as a network. Asking whether the brain behaves as a network means asking whether network properties count. From the viewpoint of neurophysiology and, possibly, of brain physics, the most substantial issues a network structure may be instrumental in addressing relate to the influence of network properties on brain dynamics and to whether these properties ultimately explain some aspects of brain function. Here, we address the dynamical implications of complex network, examining which aspects and scales of brain activity may be understood to genuinely behave as a network. To do so, we first define the meaning of networkness, and analyse some of its implications. We then examine ways in which brain anatomy and dynamics can be endowed with a network structure and discuss possible ways in which network structure may be shown to represent a genuine organisational principle of brain activity, rather than just a convenient description of its anatomy and dynamics.

Identifying sources of human interictal discharges with travelling wave and white matter propagation

Interictal epileptiform discharges have been shown to propagate from focal epileptogenic sources as travelling waves or through more rapid white matter conduction. We hypothesize that both modes of propagation are necessary to explain interictal discharge timing delays. We propose a method that, for the first time, incorporates both propagation modes to identify unique potential sources of interictal activity. We retrospectively analysed 38 focal epilepsy patients who underwent intracranial EEG recordings and diffusion-weighted imaging for epilepsy surgery evaluation. Interictal discharges were detected and localized to the most likely source based on relative delays in time of arrival across electrodes, incorporating travelling waves and white matter propagation. We assessed the influence of white matter propagation on distance of spread, timing and clinical interpretation of interictal activity. To evaluate accuracy, we compared our source localization results to earliest spiking regions to predict seizure outcomes. White matter propagation helps to explain the timing delays observed in interictal discharge sequences, underlying rapid and distant propagation. Sources identified based on differences in time of receipt of interictal discharges are often distinct from the leading electrode location. Receipt of activity propagating rapidly via white matter can occur earlier than more local activity propagating via slower cortical travelling waves. In our cohort, our source localization approach was more accurate in predicting seizure outcomes than the leading electrode location. Inclusion of white matter in addition to travelling wave propagation in our model of discharge spread did not improve overall accuracy but allowed for identification of unique and at times distant potential sources of activity, particularly in patients with persistent postoperative seizures. Since distant white matter propagation can occur more rapidly than local travelling wave propagation, combined modes of propagation within an interictal discharge sequence can decouple the commonly assumed relationship between spike timing and distance from the source. Our findings thus highlight the clinical importance of recognizing the presence of dual modes of propagation during interictal discharges, as this may be a cause of clinical mislocalization.

Spike propagation mapping reveals effective connectivity and predicts surgical outcome in epilepsy

Neurosurgical intervention is the best available treatment for selected patients with drug resistant epilepsy. For these patients, surgical planning requires biomarkers that delineate the epileptogenic zone, the brain area that is indispensable for the generation of seizures. Interictal spikes recorded with electrophysiological techniques are considered key biomarkers of epilepsy. Yet, they lack specificity, mostly because they propagate across brain areas forming networks. Understanding the relationship between interictal spike propagation and functional connections among the involved brain areas may help develop novel biomarkers that can delineate the epileptogenic zone with high precision. Here, we reveal the relationship between spike propagation and effective connectivity among onset and areas of spread and assess the prognostic value of resecting these areas. We analysed intracranial EEG data from 43 children with drug resistant epilepsy who underwent invasive monitoring for neurosurgical planning. Using electric source imaging, we mapped spike propagation in the source domain and identified three zones: onset, early-spread and late-spread. For each zone, we calculated the overlap and distance from surgical resection. We then estimated a virtual sensor for each zone and the direction of information flow among them via Granger causality. Finally, we compared the prognostic value of resecting these zones, the clinically-defined seizure onset zone and the spike onset on intracranial EEG channels by estimating their overlap with resection. We observed a spike propagation in source space for 37 patients with a median duration of 95 ms (interquartile range: 34-206), a spatial displacement of 14 cm (7.5-22 cm) and a velocity of 0.5 m/s (0.3-0.8 m/s). In patients with good surgical outcome (25 patients, Engel I), the onset had higher overlap with resection [96% (40-100%)] than early-spread [86% (34-100%), P = 0.01] and late-spread [59% (12-100%), P = 0.002], and it was also closer to resection than late-spread [5 mm versus 9 mm, P = 0.007]. We found an information flow from onset to early-spread in 66% of patients with good outcomes, and from early-spread to onset in 50% of patients with poor outcome. Finally, resection of spike onset, but not area of spike spread or the seizure onset zone, predicted outcome with positive predictive value of 79% and negative predictive value of 56% (P = 0.04). Spatiotemporal mapping of spike propagation reveals information flow from onset to areas of spread in epilepsy brain. Surgical resection of the spike onset disrupts the epileptogenic network and may render patients with drug resistant epilepsy seizure-free without having to wait for a seizure to occur during intracranial monitoring.

The electrical properties of isolated microtubules

This study examines the electrical properties of isolated brain microtubules (MTs), which are long hollow cylinders assembled from αβ-tubulin dimers that form cytoskeletal structures engaged in several functions. MTs are implicated in sensory functions in cilia and flagella and cellular activities that range from cell motility, vesicular traffic, and neuronal processes to cell division in the centrosomes and centrioles. We determined the electrical properties of the MTs with the loose patch clamp technique in either the presence or absence of the MT stabilizer Paclitaxel. We observed electrical oscillations at different holding potentials that responded accordingly in amplitude and polarity. At zero mV in symmetrical ionic conditions, a single MT radiated an electrical power of 10-17 W. The spectral analysis of the time records disclosed a single fundamental peak at 39 Hz in the Paclitaxel-stabilized MTs. However, a richer oscillatory response and two mean conductances were observed in the non-Paclitaxel MTs. The findings evidence that the brain MTs are electrical oscillators that behave as "ionic-based" transistors to generate, propagate, and amplify electrical signals.

Visual cortical LFP in relation to the hippocampal theta rhythm in track running rats

Theta oscillations in the primary visual cortex (VC) have been observed during running tasks, but the mechanism behind their generation is not well understood. Some studies have suggested that theta in the VC is locally generated, while others have proposed that it is volume conducted from the hippocampus. The present study aimed to investigate the relationship between hippocampal and VC LFP dynamics. Analysis of power spectral density revealed that LFP in the VC was similar to that in the hippocampus, but with lower overall magnitude. As running velocity increased, both the power and frequency of theta and its harmonics increased in the VC, similarly to what is observed in the hippocampus. Current source density analysis triggered to theta did not identify distinct current sources and sinks in the VC, supporting the idea that theta in the VC is conducted from the adjacent hippocampus. Phase coupling between theta, its harmonics, and gamma is a notable feature in the hippocampus, particularly in the lacunosum moleculare. While some evidence of coupling between theta and its harmonics in the VC was found, bicoherence estimates did not reveal significant phase coupling between theta and gamma. Similar results were seen in the cross-region bicoherence analysis, where theta showed strong coupling with its harmonics with increasing velocity. Thus, theta oscillations observed in the VC during running tasks are likely due to volume conduction from the hippocampus.

Interictal discharges in the human brain are travelling waves arising from an epileptogenic source

While seizure activity may be electrographically widespread, increasing evidence has suggested that ictal discharges may in fact represent travelling waves propagated from a focal seizure source. Interictal epileptiform discharges (IEDs) are an electrographic manifestation of excessive hypersynchronization of cortical activity that occur between seizures and are considered a marker of potentially epileptogenic tissue. The precise relationship between brain regions demonstrating IEDs and those involved in seizure onset, however, remains poorly understood. Here, we hypothesize that IEDs likewise reflect the receipt of travelling waves propagated from the same regions which give rise to seizures. Forty patients from our institution who underwent invasive monitoring for epilepsy, proceeded to surgery and had at least one year of follow-up were included in our study. Interictal epileptiform discharges were detected using custom software, validated by a clinical epileptologist. We show that IEDs reach electrodes in sequences with a consistent temporal ordering, and this ordering matches the timing of receipt of ictal discharges, suggesting that both types of discharges spread as travelling waves. We use a novel approach for localization of ictal discharges, in which time differences of discharge receipt at nearby electrodes are used to compute source location; similar algorithms have been used in acoustics and geophysics. We find that interictal discharges co-localize with ictal discharges. Moreover, interictal discharges tend to localize to the resection territory in patients with good surgical outcome and outside of the resection territory in patients with poor outcome. The seizure source may originate at, and also travel to, spatially distinct IED foci. Our data provide evidence that interictal discharges may represent travelling waves of pathological activity that are similar to their ictal counterparts, and that both ictal and interictal discharges emerge from common epileptogenic brain regions. Our findings have important clinical implications, as they suggest that seizure source localizations may be derived from interictal discharges, which are much more frequent than seizures.

Subthreshold Oscillating Waves in Neural Tissue Propagate by Volume Conduction and Generate Interference

Subthreshold neural oscillations have been observed in several brain regions and can influence the timing of neural spikes. However, the spatial extent and function of these spontaneous oscillations remain unclear. To study the mechanisms underlying these oscillations, we use optogenetic stimulation to generate oscillating waves in the longitudinal hippocampal slice expressing optopatch proteins. We found that optogenetic stimulation can generate two types of neural activity: suprathreshold neural spikes and subthreshold oscillating waves. Both waves could propagate bidirectionally at similar speeds and go through a transection of the tissue. The propagating speed is independent of the oscillating frequency but increases with increasing amplitudes of the waves. The endogenous electric fields generated by oscillating waves are about 0.6 mV/mm along the dendrites and about 0.3 mV/mm along the cell layer. We also observed that these oscillating waves could interfere with each other. Optical stimulation applied simultaneously at each slice end generated a larger wave in the middle of the tissue (constructive interference) or destructive interference with laser signals in opposite phase. However, the suprathreshold neural spikes were annihilated when they collided. Finally, the waves were not affected by the NMDA blocker (APV) and still propagated in the presence of tetrodotoxin (TTX) but at a significantly lower amplitude. The role of these subthreshold waves in neural function is unknown, but the results show that at low amplitude, the subthreshold propagating waves lack a refractory period allowing a novel analog form of preprocessing of neural activity by interference independent of synaptic transmission.

Theta waves, neural spikes and seizures can propagate by ephaptic coupling in vivo

Electric field coupling has been shown to be responsible for non-synaptic neural activity propagation in hippocampal slices and cortical slices. Epileptiform and slow-wave sleep activity can propagate by electric field coupling without using synaptic connections at speeds of ~0.1 m/s in vitro. However, the characteristics of the events that can propagate using electric field coupling through a volume conductor in vivo have not been studied. Thus, we tested the hypothesis that various types of neural signals such as interictal spikes, theta waves and seizures could propagate in vivo across a transection in the hippocampus. We induced epileptiform activity in 4 rats under anesthesia by injecting 4-aminopyridine in the temporal region of the hippocampus, four recording electrodes were inserted along the longitudinal axis of the hippocampus. A transection was made between the electrodes to study the propagation of the neural activity. Although 54% of the interictal spikes could propagate through the cut, only those spikes with a high amplitude and short duration had a high probability to do so. 70% of seizure events could propagate through the cut but parameters distinguishing between propagating and non-propagating seizure events could not be identified. Theta activity was also observed to propagate at a mean speed of 0.16 ± 0.12 m/s in the characteristic range of propagation using electric field coupling through the transection. The electric field volume conduction mechanism was confirmed by showing that propagation was blocked by placing a dielectric layer within the cut. The speed of propagation was not affected by the transection thereby providing further evidence that various types of neural signals including activity in the theta range can propagate by electric field coupling in-vivo.

Repeated hippocampal seizures lead to brain-wide reorganization of circuits and seizure propagation pathways

Repeated seizure activity can lead to long-term changes in seizure dynamics and behavior. However, resulting changes in brain-wide dynamics remain poorly understood. This is due partly to technical challenges in precise seizure control and in vivo whole-brain mapping of circuit dynamics. Here, we developed an optogenetic kindling model through repeated stimulation of ventral hippocampal CaMKII neurons in adult rats. We then combined fMRI with electrophysiology to track brain-wide circuit dynamics resulting from non-afterdischarge (AD)-generating stimulations and individual convulsive seizures. Kindling induced widespread increases in non-AD-generating stimulation response and ipsilateral functional connectivity and elevated anxiety. Individual seizures in kindled animals showed more significant increases in brain-wide activity and bilateral functional connectivity. Onset time quantification provided evidence for kindled seizure propagation from the ipsilateral to the contralateral hemisphere. Furthermore, a core of slow-migrating hippocampal activity was identified in both non-kindled and kindled seizures, revealing a novel mechanism of seizure sustainment and propagation.

Mesoscopic Mapping of Ictal Neurovascular Coupling in Awake Behaving Mice Using Optical Spectroscopy and Genetically Encoded Calcium Indicators

Unambiguously identifying an epileptic focus with high spatial resolution is a challenge, especially when no anatomic abnormality can be detected. Neurovascular coupling (NVC)-based brain mapping techniques are often applied in the clinic despite a poor understanding of ictal NVC mechanisms, derived primarily from recordings in anesthetized animals with limited spatial sampling of the ictal core. In this study, we used simultaneous wide-field mesoscopic imaging of GCamp6f and intrinsic optical signals (IOS) to record the neuronal and hemodynamic changes during acute ictal events in awake, behaving mice. Similar signals in isoflurane-anesthetized mice were compared to highlight the unique characteristics of the awake condition. In awake animals, seizures were more focal at the onset but more likely to propagate to the contralateral hemisphere. The HbT signal, derived from an increase in cerebral blood volume (CBV), was more intense in awake mice. As a result, the “epileptic dip” in hemoglobin oxygenation became inconsistent and unreliable as a mapping signal. Our data indicate that CBV-based imaging techniques should be more accurate than blood oxygen level dependent (BOLD)-based imaging techniques for seizure mapping in awake behaving animals.

Neural recruitment by ephaptic coupling in epilepsy

OBJECTIVE

One of the challenges in treating patients with drug-resistant epilepsy is that the mechanisms of seizures are unknown. Most current interventions are based on the assumption that epileptic activity recruits neurons and progresses by synaptic transmission. However, several experimental studies have shown that neural activity in rodent hippocampi can propagate independently of synaptic transmission. Recent studies suggest these waves are self-propagating by electric field (ephaptic) coupling. In this study, we tested the hypothesis that neural recruitment during seizures can occur by electric field coupling.

METHODS

4-Aminopyridine was used in both in vivo and in vitro preparation to trigger seizures or epileptiform activity. A transection was made in the in vivo hippocampus and in vitro hippocampal and cortical slices to study whether the induced seizure activity can recruit neurons across the gap. A computational model was built to test whether ephaptic coupling alone can account for neural recruitment across the transection. The model prediction was further validated by in vitro experiments.

RESULTS

Experimental results show that electric fields generated by seizure-like activity in the hippocampus both in vitro and in vivo can recruit neurons locally and through a transection of the tissue. The computational model suggests that the neural recruitment across the transection is mediated by electric field coupling. With in vitro experiments, we show that a dielectric material can block the recruitment of epileptiform activity across a transection, and that the electric fields measured within the gap are similar to those predicted by model simulations. Furthermore, this nonsynaptic neural recruitment is also observed in cortical slices, suggesting that this effect is robust in brain tissue.

SIGNIFICANCE

These results indicate that ephaptic coupling, a nonsynaptic mechanism, can underlie neural recruitment by a small electric field generated by seizure activity and could explain the low success rate of surgical transections in epilepsy patients.

Rapidly spreading seizures arise from large-scale functional brain networks in focal epilepsy

It remains unclear whether epileptogenic networks in focal epilepsy develop on physiological networks. This work aimed to explore the association between the rapid spread of ictal fast activity (IFA), a proposed biomarker for epileptogenic networks, and the functional connectivity or networks of healthy subjects. We reviewed 45 patients with focal epilepsy who underwent electrocorticographic (ECoG) recordings to identify the patients showing the rapid spread of IFA. IFA power was quantified as normalized beta-gamma band power. Using published resting-state functional magnetic resonance imaging databases, we estimated resting-state functional connectivity of healthy subjects (RSFC-HS) and resting-state networks of healthy subjects (RSNs-HS) at the locations corresponding to the patients' electrodes. We predicted the IFA power of each electrode based on RSFC-HS between electrode locations (RSFC-HS-based prediction) using a recently developed method, termed activity flow mapping. RSNs-HS were identified using seed-based and atlas-based methods. We compared IFA power with RSFC-HS-based prediction or RSNs-HS using non-parametric correlation coefficients. RSFC and seed-based RSNs of each patient (RSFC-PT and seed-based RSNs-PT) were also estimated using interictal ECoG data and compared with IFA power in the same way as RSFC-HS and seed-based RSNs-HS. Spatial autocorrelation-preserving randomization tests were performed for significance testing. Nine patients met the inclusion criteria. None of the patients had reflex seizures. Six patients showed pathological evidence of a structural etiology. In total, we analyzed 49 seizures (2-13 seizures per patient). We observed significant correlations between IFA power and RSFC-HS-based prediction, seed-based RSNs-HS, or atlas-based RSNs-HS in 28 (57.1%), 21 (42.9%), and 28 (57.1%) seizures, respectively. Thirty-two (65.3%) seizures showed a significant correlation with either seed-based or atlas-based RSNs-HS, but this ratio varied across patients: 27 (93.1%) of 29 seizures in six patients correlated with either of them. Among atlas-based RSNs-HS, correlated RSNs-HS with IFA power included the default mode, control, dorsal attention, somatomotor, and temporal-parietal networks. We could not obtain RSFC-PT and RSNs-PT in one patient due to frequent interictal epileptiform discharges. In the remaining eight patients, most of the seizures showed significant correlations between IFA power and RSFC-PT-based prediction or seed-based RSNs-PT. Our study provides evidence that the rapid spread of IFA in focal epilepsy can arise from physiological RSNs. This finding suggests an overlap between epileptogenic and functional networks, which may explain why functional networks in patients with focal epilepsy frequently disrupt.

Travelling waves reveal a dynamic seizure source in human focal epilepsy

Treatment of patients with drug-resistant focal epilepsy relies upon accurate seizure localization. Ictal activity captured by intracranial EEG has traditionally been interpreted to suggest that the underlying cortex is actively involved in seizures. Here, we hypothesize that such activity instead reflects propagated activity from a relatively focal seizure source, even during later time points when ictal activity is more widespread. We used the time differences observed between ictal discharges in adjacent electrodes to estimate the location of the hypothesized focal source and demonstrated that the seizure source, localized in this manner, closely matches the clinically and neurophysiologically determined brain region giving rise to seizures. Moreover, we determined this focal source to be a dynamic entity that moves and evolves over the time course of a seizure. Our results offer an interpretation of ictal activity observed by intracranial EEG that challenges the traditional conceptualization of the seizure source.

Modulation of intercolumnar synchronization by endogenous electric fields in cerebral cortex

Neurons synaptically interacting in a conductive medium generate extracellular endogenous electric fields (EFs) that reciprocally affect membrane potential. Exogenous EFs modulate neuronal activity, and their clinical applications are being profusely explored. However, whether endogenous EFs contribute to network synchronization remains unclear. We analyzed spontaneously generated slow-wave activity in the cerebral cortex network in vitro, which allowed us to distinguish synaptic from nonsynaptic mechanisms of activity propagation and synchronization. Slow oscillations generated EFs that propagated independently of synaptic transmission. We demonstrate that cortical oscillations modulate spontaneous rhythmic activity of neighboring synaptically disconnected cortical columns if layers are aligned. We provide experimental evidence that these EF-mediated effects are compatible with electric dipoles. With a model of interacting dipoles, we reproduce the experimental measurements and predict that endogenous EF-mediated synchronizing effects should be relevant in the brain. Thus, experiments and models suggest that electric-dipole interactions contribute to synchronization of neighboring cortical columns.

The potential antiepileptogenic effect of neuronal Cx36 gap junction channel blockage

Epilepsy is one of the most prevalent neurological disorders and can result in neuronal injury and degeneration. Consequently, research into new antiepileptic drugs capable of providing protection against neuronal injury and degeneration is extremely important. Neuronal Cx36 gap junction channels have been found to play an important role in epilepsy; thus, pharmacological interference using Cx36 gap junction channel blockers may be a promising strategy for disrupting the synchronization of neurons during seizure activity and protecting neurons. Based on these promising findings, several in vivo and in vitro studies are ongoing and the first encouraging results have been published. The results bring hope that neurons can be protected from injury and degeneration in patients with epilepsy, which is currently impossible.

Brain Waves: Emergence of Localized, Persistent, Weakly Evanescent Cortical Loops

An inhomogeneous anisotropic physical model of the brain cortex is presented that predicts the emergence of nonevanescent (weakly damped) wave-like modes propagating in the thin cortex layers transverse to both the mean neural fiber direction and the cortex spatial gradient. Although the amplitude of these modes stays below the typically observed axon spiking potential, the lifetime of these modes may significantly exceed the spiking potential inverse decay constant. Full-brain numerical simulations based on parameters extracted from diffusion and structural MRI confirm the existence and extended duration of these wave modes. Contrary to the commonly agreed paradigm that the neural fibers determine the pathways for signal propagation in the brain, the signal propagation because of the cortex wave modes in the highly folded areas will exhibit no apparent correlation with the fiber directions. Nonlinear coupling of those linear weakly evanescent wave modes then provides a universal mechanism for the emergence of synchronized brain wave field activity. The resonant and nonresonant terms of nonlinear coupling between multiple modes produce both synchronous spiking-like high-frequency wave activity as well as low-frequency wave rhythms. Numerical simulation of forced multiple-mode dynamics shows that, as forcing increases, there is a transition from damped to oscillatory regime that can then transition quickly to a nonoscillatory state when a critical excitation threshold is reached. The resonant nonlinear coupling results in the emergence of low-frequency rhythms with frequencies that are several orders of magnitude below the linear frequencies of modes taking part in the coupling. The localization and persistence of these weakly evanescent cortical wave modes have significant implications in particular for neuroimaging methods that detect electromagnetic physiological activity, such as EEG and magnetoencephalography, and for the understanding of brain activity in general, including mechanisms of memory.

Realistic modeling of mesoscopic ephaptic coupling in the human brain

Several decades of research suggest that weak electric fields may influence neural processing, including those induced by neuronal activity and proposed as a substrate for a potential new cellular communication system, i.e., ephaptic transmission. Here we aim to model mesoscopic ephaptic activity in the human brain and explore its trajectory during aging by characterizing the electric field generated by cortical dipoles using realistic finite element modeling. Extrapolating from electrophysiological measurements, we first observe that modeled endogenous field magnitudes are comparable to those in measurements of weak but functionally relevant self-generated fields and to those produced by noninvasive transcranial brain stimulation, and therefore possibly able to modulate neuronal activity. Then, to evaluate the role of these fields in the human cortex in large MRI databases, we adapt an interaction approximation that considers the relative orientation of neuron and field to estimate the membrane potential perturbation in pyramidal cells. We use this approximation to define a simplified metric (EMOD1) that weights dipole coupling as a function of distance and relative orientation between emitter and receiver and evaluate it in a sample of 401 realistic human brain models from healthy subjects aged 16-83. Results reveal that ephaptic coupling, in the simplified mesoscopic modeling approach used here, significantly decreases with age, with higher involvement of sensorimotor regions and medial brain structures. This study suggests that by providing the means for fast and direct interaction between neurons, ephaptic modulation may contribute to the complexity of human function for cognition and behavior, and its modification across the lifespan and in response to pathology.

Universal theory of brain waves: from linear loops to nonlinear synchronized spiking and collective brain rhythms

An inhomogeneous anisotropic physical model of the brain cortex is presented that predicts the emergence of non-evanescent (weakly damped) wave-like modes propagating in the thin cortex layers transverse to both the mean neural fiber direction and to the cortex spatial gradient. Although the amplitude of these modes stays below the typically observed axon spiking potential, the lifetime of these modes may significantly exceed the spiking potential inverse decay constant. Full brain numerical simulations based on parameters extracted from diffusion and structural MRI confirm the existence and extended duration of these wave modes. Contrary to the standard paradigm that the neural fibers determine the pathways for signal propagation in the brain, the signal propagation due to the cortex wave modes in highly folded areas will exhibit no apparent correlation with the fiber directions. The results are consistent with numerous recent experimental animal and human brain studies demonstrating the existence of electrostatic field activity in the form of traveling waves (including studies where neuronal connections were severed) and with wave loop induced peaks observed in EEG spectra. In addition, we demonstrate that the resonant and non-resonant terms of the nonlinear coupling between multiple modes produce both synchronous spiking-like high frequency wave activity as well as low frequency wave rhythms as a result of their unique dispersion properties. Numerical simulation of forced multiple mode dynamics shows that as forcing increases there is a transition from damped to oscillatory regime that subsequently decays away as over-excitation is reached. The resonant nonlinear coupling results in the emergence of low frequency rhythms with frequencies that are several orders of magnitude below the linear frequencies of modes taking part in the coupling. The localization and persistence of these cortical wave modes, and this new mechanism for understanding the nature of spiking behavior, have significant implications in particular for neuroimaging methods that detect electromagnetic physiological activity, such as EEG and MEG, and in general for the understanding of brain activity, including mechanisms of memory.

Ca2+ Signals in Astrocytes Facilitate Spread of Epileptiform Activity

Epileptic seizures are associated with increased astrocytic Ca²⁺ signaling, but the fine spatiotemporal kinetics of the ictal astrocyte–neuron interplay remains elusive. By using 2-photon imaging of awake head-fixed mice with chronic hippocampal windows we demonstrate that astrocytic Ca²⁺ signals precede neuronal Ca²⁺ elevations during the initial bout of kainate-induced seizures. On average, astrocytic Ca²⁺ elevations preceded neuronal activity in CA1 by about 8 s. In subsequent bouts of epileptic seizures, astrocytes and neurons were activated simultaneously. The initial astrocytic Ca²⁺ elevation was abolished in mice lacking the type 2 inositol-1,4,5-trisphosphate-receptor (Itpr2^−/−). Furthermore, we found that Itpr2^−/− mice exhibited 60% less epileptiform activity compared with wild-type mice when assessed by telemetric EEG monitoring. In both genotypes we also demonstrate that spreading depression waves may play a part in seizure termination. Our findings imply a role for astrocytic Ca²⁺ signals in ictogenesis.

Investigating the role of gap junctions in seizure wave propagation

The effect of gap junctions as well as the biological mechanisms behind seizure wave propagation is not completely understood. In this work, we use a simple neural field model to study the possible influence of gap junctions specifically on cortical wave propagation that has been observed in vivo preceding seizure termination. We consider a voltage-based neural field model consisting of an excitatory and an inhibitory population as well as both chemical and gap junction-like synapses. We are able to approximate important properties of cortical wave propagation previously observed in vivo before seizure termination. This model adds support to existing evidence from models and clinical data suggesting a key role of gap junctions in seizure wave propagation. In particular, we found that in this model gap junction-like connectivity determines the propagation of one-bump or two-bump traveling wave solutions with features consistent with the clinical data. For sufficiently increased gap junction connectivity, wave solutions cease to exist. Moreover, gap junction connectivity needs to be sufficiently low or moderate to permit the existence of linearly stable solutions of interest.

Optogenetically-Induced Population Discharge Threshold as a Sensitive Measure of Network Excitability

Network excitability is governed by synaptic efficacy, intrinsic excitability, and the circuitry in which these factors are expressed. The complex interplay between these factors determines how circuits function and, at the extreme, their susceptibility to seizure. We have developed a sensitive, quantitative estimate of network excitability in freely behaving mice using a novel optogenetic intensity-response procedure. Synchronous activation of deep sublayer CA1 pyramidal cells produces abnormal network-wide epileptiform population discharges (PDs) that are nearly indistinguishable from spontaneously-occurring interictal spikes (IISs). By systematically varying light intensity, and therefore the magnitude of the optogenetically-mediated current, we generated intensity-response curves using the probability of PD as the dependent variable. Manipulations known to increase excitability, such as sub-convulsive doses (20 mg/kg) of the chemoconvulsant pentylenetetrazol (PTZ), produced a leftward shift in the curve compared to baseline. The anti-epileptic drug levetiracetam (LEV; 40 mk/kg), in combination with PTZ, produced a rightward shift. Optogenetically-induced PD threshold (oPDT) baselines were stable over time, suggesting the metric is appropriate for within-subject experimental designs with multiple pharmacological manipulations.

From cells to sensations: A window to the physics of mind

Principles of methods for studying particles and fields that cannot be sensed by third-person observers by routine methods can be used to understand the physics of first-person properties of mind. Accordingly, whenever a system exhibits disparate features at multiple levels, unique combination of constraints offered by them direct us towards a solution that will be the first principle of that system. Using this method, it was possible to arrive at a third-person observable solution-point of brain-mind interface. Examination of this location identified a set of unique features that can allow an associatively learned (cue) stimulus to spark hallucinations that form units of first-person internal (inner) sensations reminiscent of stimuli from the associatively learned second item in timescales of milliseconds. It allows us to peep into a virtual space of mind where different modifications and integrations of units of internal sensations generate their different net conformations ranging from perception to an inner sense of hidden relationships that form a hypothesis. Since sparking of inner sensations of the late arriving (when far away) or non-arriving (when hidden) features of items started providing survival advantage, the focus of evolution might have been to optimize this property. Hence, the circuity that generates it can be considered as the primary circuitry of the system. The solution provides several testable predictions. By taking readers through the process of deriving the solution and by explaining how it interconnects disparate findings, it is hoped that the factors determining the physics of mind will become evident.

A hypothesis regarding how sleep can calibrate neuronal excitability in the central nervous system and thereby offer stability, sensitivity and the best possible cognitive function

The function of sleep in mammal and other vertebrates is one of the great mysteries of biology. Many hypotheses have been proposed, but few of these have made even the slightest attempt to explain the essence of sleep - the uncompromising need for reversible unconsciousness. During sleep, epiphenomena - often of a somatic character - occur, but these cannot explain the core function of sleep. One answer could be hidden in the observations made for long periods of time of the function of the central nervous system (CNS). The CNS is faced with conflicting requirements on stability and excitability. A high level of excitability is desirable, and is also a prerequisite for sensitivity and quick reaction times; however, it can also lead to instability and the risk of feedback, with life-threatening epileptic seizures. Activity-dependent negative feedback in neuronal excitability improves stability in the short term, but not to the degree that is required. A hypothesis is presented here demonstrating how calibration of individual neurons - an activity which occurs only during sleep - can establish the balanced and highest possible excitability while also preserving stability in the CNS. One example of a possible mechanism is the observation of slow oscillations in EEGs made on birds and mammals during slow wave sleep. Calibration to a genetically determined level of excitability could take place in individual neurons during the slow oscillation. This is only possible offline, which explains the need for sleep. The hypothesis can explain phenomena such as the need for unconsciousness during sleep, with the disconnection of sensory stimuli, slow EEG oscillations, the relationship of sleep and epilepsy, age, the effects of sleep on neuronal firing rate and the effects of sleep deprivation and sleep homeostasis. This is with regard primarily to mammals, including humans, but also all other vertebrates.

Role of inhibitory control in modulating focal seizure spread

Focal seizure propagation is classically thought to be spatially contiguous. However, distribution of seizures through a large-scale epileptic network has been theorized. Here, we used a multielectrode array, wide field calcium imaging, and two-photon calcium imaging to study focal seizure propagation pathways in an acute rodent neocortical 4-aminopyridine model. Although ictal neuronal bursts did not propagate beyond a 2-3-mm region, they were associated with hemisphere-wide field potential fluctuations and parvalbumin-positive interneuron activity outside the seizure focus. While bicuculline surface application enhanced contiguous seizure propagation, focal bicuculline microinjection at sites distant to the 4-aminopyridine focus resulted in epileptic network formation with maximal activity at the two foci. Our study suggests that both classical and epileptic network propagation can arise from localized inhibition defects, and that the network appearance can arise in the context of normal brain structure without requirement for pathological connectivity changes between sites.

Ionic and synaptic mechanisms of seizure generation and epileptogenesis

The biophysical mechanisms underlying epileptogenesis and the generation of seizures remain to be better understood. Among many factors triggering epileptogenesis are traumatic brain injury breaking normal synaptic homeostasis and genetic mutations disrupting ionic concentration homeostasis. Impairments in these mechanisms, as seen in various brain diseases, may push the brain network to a pathological state characterized by increased susceptibility to unprovoked seizures. Here, we review recent computational studies exploring the roles of ionic concentration dynamics in the generation, maintenance, and termination of seizures. We further discuss how ionic and synaptic homeostatic mechanisms may give rise to conditions which prime brain networks to exhibit recurrent spontaneous seizures and epilepsy.

Self-propagating, non-synaptic epileptiform activity recruits neurons by endogenous electric fields

It is well documented that synapses play a significant role in the transmission of information between neurons. However, in the absence of synaptic transmission, neural activity has been observed to continue to propagate. Previous studies have shown that propagation of epileptiform activity takes place in the absence of synaptic transmission and gap junctions and is outside the range of ionic diffusion and axonal conduction. Computer simulations indicate that electric field coupling could be responsible for the propagation of neural activity under pathological conditions such as epilepsy. Electric fields can modulate neuronal membrane voltage, but there is no experimental evidence suggesting that electric field coupling can mediate self-regenerating propagation of neural activity. Here we examine the role of electric field coupling by eliminating all forms of neural communications except electric field coupling with a cut through the neural tissue. We show that 4-AP induced activity generates an electric field capable of recruiting neurons on the distal side of the cut. Experiments also show that applied electric fields with amplitudes similar to endogenous values can induce propagating waves. Finally, we show that canceling the electrical field at a given point can block spontaneous propagation. The results from these in vitro electrophysiology experiments suggest that electric field coupling is a critical mechanism for non-synaptic neural propagation and therefore could contribute to the propagation of epileptic activity in the brain.

Slow periodic activity in the longitudinal hippocampal slice can self-propagate non-synaptically by a mechanism consistent with ephaptic coupling

KEY POINTS

Slow periodic activity can propagate with speeds around 0.1 m s-1 and be modulated by weak electric fields. Slow periodic activity in the longitudinal hippocampal slice can propagate without chemical synaptic transmission or gap junctions, but can generate electric fields which in turn activate neighbouring cells. Applying local extracellular electric fields with amplitude in the range of endogenous fields is sufficient to modulate or block the propagation of this activity both in the in silico and in the in vitro models. Results support the hypothesis that endogenous electric fields, previously thought to be too small to trigger neural activity, play a significant role in the self-propagation of slow periodic activity in the hippocampus. Experiments indicate that a neural network can give rise to sustained self-propagating waves by ephaptic coupling, suggesting a novel propagation mechanism for neural activity under normal physiological conditions.

ABSTRACT

Slow oscillations are a standard feature observed in the cortex and the hippocampus during slow wave sleep. Slow oscillations are characterized by low-frequency periodic activity (<1 Hz) and are thought to be related to memory consolidation. These waves are assumed to be a reflection of the underlying neural activity, but it is not known if they can, by themselves, be self-sustained and propagate. Previous studies have shown that slow periodic activity can be reproduced in the in vitro preparation to mimic in vivo slow oscillations. Slow periodic activity can propagate with speeds around 0.1 m s-1 and be modulated by weak electric fields. In the present study, we show that slow periodic activity in the longitudinal hippocampal slice is a self-regenerating wave which can propagate with and without chemical or electrical synaptic transmission at the same speeds. We also show that applying local extracellular electric fields can modulate or even block the propagation of this wave in both in silico and in vitro models. Our results support the notion that ephaptic coupling plays a significant role in the propagation of the slow hippocampal periodic activity. Moreover, these results indicate that a neural network can give rise to sustained self-propagating waves by ephaptic coupling, suggesting a novel propagation mechanism for neural activity under normal physiological conditions.

Recording and Modulation of Epileptiform Activity in Rodent Brain Slices Coupled to Microelectrode Arrays

Temporal lobe epilepsy (TLE) is the most common partial complex epileptic syndrome and the least responsive to medications. Deep brain stimulation (DBS) is a promising approach when pharmacological treatment fails or neurosurgery is not recommended. Acute brain slices coupled to microelectrode arrays (MEAs) represent a valuable tool to study neuronal network interactions and their modulation by electrical stimulation. As compared to conventional extracellular recording techniques, they provide the added advantages of a greater number of observation points and a known inter-electrode distance, which allow studying the propagation path and speed of electrophysiological signals. However, tissue oxygenation may be greatly impaired during MEA recording, requiring a high perfusion rate, which comes at the cost of decreased signal-to-noise ratio and higher oscillations in the experimental temperature. Electrical stimulation further stresses the brain tissue, making it difficult to pursue prolonged recording/stimulation epochs. Moreover, electrical modulation of brain slice activity needs to target specific structures/pathways within the brain slice, requiring that electrode mapping be easily and quickly performed live during the experiment. Here, we illustrate how to perform the recording and electrical modulation of 4-aminopyridine (4AP)-induced epileptiform activity in rodent brain slices using planar MEAs. We show that the brain tissue obtained from mice outperforms rat brain tissue and is thus better suited for MEA experiments. This protocol guarantees the generation and maintenance of a stable epileptiform pattern that faithfully reproduces the electrophysiological features observed with conventional field potential recording, persists for several hours, and outlasts sustained electrical stimulation for prolonged epochs. Tissue viability throughout the experiment is achieved thanks to the use of a small-volume custom recording chamber allowing for laminar flow and quick solution exchange even at low (1 mL/min) perfusion rates. Quick MEA mapping for real-time monitoring and selection of stimulating electrodes is performed by a custom graphic user interface (GUI).

Slow moving neural source in the epileptic hippocampus can mimic progression of human seizures

Fast and slow neural waves have been observed to propagate in the human brain during seizures. Yet the nature of these waves is difficult to study in a surgical setting. Here, we report an observation of two different traveling waves propagating in the in-vitro epileptic hippocampus at speeds similar to those in the human brain. A fast traveling spike and a slow moving wave were recorded simultaneously with a genetically encoded voltage sensitive fluorescent protein (VSFP Butterfly 1.2) and a high speed camera. The results of this study indicate that the fast traveling spike is NMDA-sensitive but the slow moving wave is not. Image analysis and model simulation demonstrate that the slow moving wave is moving slowly, generating the fast traveling spike and is, therefore, a moving source of the epileptiform activity. This slow moving wave is associated with a propagating neural calcium wave detected with calcium dye (OGB-1) but is independent of NMDA receptors, not related to ATP release, and much faster than those previously recorded potassium waves. Computer modeling suggests that the slow moving wave can propagate by the ephaptic effect like epileptiform activity. These findings provide an alternative explanation for slow propagation seizure wavefronts associated with fast propagating spikes.

Brain extracellular space, hyaluronan, and the prevention of epileptic seizures

Mutant mice deficient in hyaluronan (HA) have an epileptic phenotype. HA is one of the major constituents of the brain extracellular matrix. HA has a remarkable hydration capacity, and a lack of HA causes reduced extracellular space (ECS) volume in the brain. Reducing ECS volume can initiate or exacerbate epileptiform activity in many in vitro models of epilepsy. There is both in vitro and in vivo evidence of a positive feedback loop between reduced ECS volume and synchronous neuronal activity. Reduced ECS volume promotes epileptiform activity primarily via enhanced ephaptic interactions and increased extracellular potassium concentration; however, the epileptiform activity in many models, including the brain slices from HA synthase-3 knockout mice, may still require glutamate-mediated synaptic activity. In brain slice epilepsy models, hyperosmotic solution can effectively shrink cells and thus increase ECS volume and block epileptiform activity. However, in vivo, the intravenous administration of hyperosmotic solution shrinks both brain cells and brain ECS volume. Instead, manipulations that increase the synthesis of high-molecular-weight HA or decrease its breakdown may be used in the future to increase brain ECS volume and prevent seizures in patients with epilepsy. The prevention of epileptogenesis is also a future target of HA manipulation. Head trauma, ischemic stroke, and other brain insults that initiate epileptogenesis are known to be associated with an early decrease in high-molecular-weight HA, and preventing that decrease in HA may prevent the epileptogenesis.

Methodological standards for in vitro models of epilepsy and epileptic seizures. A TASK1-WG4 report of the AES/ILAE Translational Task Force of the ILAE

In vitro preparations are a powerful tool to explore the mechanisms and processes underlying epileptogenesis and ictogenesis. In this review, we critically review the numerous in vitro methodologies utilized in epilepsy research. We provide support for the inclusion of detailed descriptions of techniques, including often ignored parameters with unpredictable yet significant effects on study reproducibility and outcomes. In addition, we explore how recent developments in brain slice preparation relate to their use as models of epileptic activity.

Utilization of independent component analysis for accurate pathological ripple detection in intracranial EEG recordings recorded extra- and intra-operatively

OBJECTIVE

To develop and validate a detector that identifies ripple (80-200 Hz) events in intracranial EEG (iEEG) recordings in a referential montage and utilizes independent component analysis (ICA) to eliminate or reduce high-frequency artifact contamination. Also, investigate the correspondence of detected ripples and the seizure onset zone (SOZ).

METHODS

iEEG recordings from 16 patients were first band-pass filtered (80-600 Hz) and Infomax ICA was next applied to derive the first independent component (IC1). IC1 was subsequently pruned, and an artifact index was derived to reduce the identification of high-frequency events introduced by the reference electrode signal. A Hilbert detector identified ripple events in the processed iEEG recordings using amplitude and duration criteria. The identified ripple events were further classified and characterized as true or false ripple on spikes, or ripples on oscillations by utilizing a topographical analysis to their time-frequency plot, and confirmed by visual inspection.

RESULTS

The signal to noise ratio was improved by pruning IC1. The precision of the detector for ripple events was 91.27 ± 4.3%, and the sensitivity of the detector was 79.4 ± 3.0% (N = 16 patients, 5842 ripple events). The sensitivity and precision of the detector was equivalent in iEEG recordings obtained during sleep or intra-operatively. Across all the patients, true ripple on spike rates and also the rates of false ripple on spikes, that were generated due to filter ringing, classified the seizure onset zone (SOZ) with an area under the receiver operating curve (AUROC) of >76%. The magnitude and spectral content of true ripple on spikes generated in the SOZ was distinct as compared with the ripples generated in the NSOZ (p < .001).

CONCLUSIONS

Utilizing ICA to analyze iEEG recordings in referential montage provides many benefits to the study of high-frequency oscillations. The ripple rates and properties defined using this approach may accurately delineate the seizure onset zone.

SIGNIFICANCE

Strategies to improve the spatial resolution of intracranial EEG and reduce artifact can help improve the clinical utility of HFO biomarkers.

Reducing Extracellular Ca2+ Induces Adenosine Release via Equilibrative Nucleoside Transporters to Provide Negative Feedback Control of Activity in the Hippocampus

Neural circuit activity increases the release of the purine neuromodulator adenosine into the extracellular space leading to A₁ receptor activation and negative feedback via membrane hyperpolarization and inhibition of transmitter release. Adenosine can be released by a number of different mechanisms that include Ca²⁺ dependent processes such as the exocytosis of ATP. During sustained pathological network activity, ischemia and hypoxia the extracellular concentration of calcium ions (Ca²⁺) markedly falls, inhibiting exocytosis and potentially reducing adenosine release. However it has been observed that reducing extracellular Ca²⁺ can induce paradoxical neural activity and can also increase adenosine release. Here we have investigated adenosine signaling and release mechanisms that occur when extracellular Ca²⁺ is removed. Using electrophysiology and microelectrode biosensor measurements we have found that adenosine is directly released into the extracellular space by the removal of extracellular Ca²⁺ and controls the induced neural activity via A₁ receptor-mediated membrane potential hyperpolarization. Following Ca²⁺ removal, adenosine is released via equilibrative nucleoside transporters (ENTs), which when blocked leads to hyper-excitation. We propose that sustained action potential firing following Ca²⁺ removal leads to hydrolysis of ATP and a build-up of intracellular adenosine which then effluxes into the extracellular space via ENTs.

Induced mitochondrial membrane potential for modeling solitonic conduction of electrotonic signals

A cable model that includes polarization-induced capacitive current is derived for modeling the solitonic conduction of electrotonic potentials in neuronal branchlets with microstructure containing endoplasmic membranes. A solution of the nonlinear cable equation modified for fissured intracellular medium with a source term representing charge ‘soakage’ is used to show how intracellular capacitive effects of bound electrical charges within mitochondrial membranes can influence electrotonic signals expressed as solitary waves. The elastic collision resulting from a head-on collision of two solitary waves results in localized and non-dispersing electrical solitons created by the nonlinearity of the source term. It has been shown that solitons in neurons with mitochondrial membrane and quasi-electrostatic interactions of charges held by the microstructure (i.e., charge ‘soakage’) have a slower velocity of propagation compared with solitons in neurons with microstructure, but without endoplasmic membranes. When the equilibrium potential is a small deviation from rest, the nonohmic conductance acts as a leaky channel and the solitons are small compared when the equilibrium potential is large and the outer mitochondrial membrane acts as an amplifier, boosting the amplitude of the endogenously generated solitons. These findings demonstrate a functional role of quasi-electrostatic interactions of bound electrical charges held by microstructure for sustaining solitons with robust self-regulation in their amplitude through changes in the mitochondrial membrane equilibrium potential. The implication of our results indicate that a phenomenological description of ionic current can be successfully modeled with displacement current in Maxwell’s equations as a conduction process involving quasi-electrostatic interactions without the inclusion of diffusive current. This is the first study in which solitonic conduction of electrotonic potentials are generated by polarization-induced capacitive current in microstructure and nonohmic mitochondrial membrane current.

Focal cortical seizures start as standing waves and propagate respecting homotopic connectivity

Focal epilepsy involves excessive cortical activity that propagates both locally and distally. Does this propagation follow the same routes as normal cortical activity? We pharmacologically induced focal seizures in primary visual cortex (V1) of awake mice, and compared their propagation to the retinotopic organization of V1 and higher visual areas. We used simultaneous local field potential recordings and widefield imaging of a genetically encoded calcium indicator to measure prolonged seizures (ictal events) and brief interictal events. Both types of event are orders of magnitude larger than normal visual responses, and both start as standing waves: synchronous elevated activity in the V1 focus and in homotopic locations in higher areas, i.e. locations with matching retinotopic preference. Following this common beginning, however, seizures persist and propagate both locally and into homotopic distal regions, and eventually invade all of visual cortex and beyond. We conclude that seizure initiation resembles the initiation of interictal events, and seizure propagation respects the connectivity underlying normal visual processing.

Focal cortical seizures result from local and widespread propagation of excitatory activity. Here the authors employ widefield calcium imaging in mouse visual areas to demonstrate that these seizures start as local synchronous activation and then propagate along the connectivity that underlies normal sensory processing.

Gamma Oscillations in the Rat Ventral Striatum Originate in the Piriform Cortex

Local field potentials (LFPs) recorded from the human and rodent ventral striatum (vStr) exhibit prominent, behaviorally relevant gamma-band oscillations. These oscillations are related to local spiking activity and transiently synchronize with anatomically related areas, suggesting a possible role in organizing vStr activity. However, the origin of vStr gamma is unknown. We recorded vStr gamma oscillations across a 1.4 mm² grid spanned by 64 recording electrodes as male rats rested and foraged for rewards, revealing a highly consistent power gradient originating in the adjacent piriform cortex. Phase differences across the vStr were consistently small (<15°) and current source density analysis further confirmed the absence of local sink-source pairs in the vStr. Reversible occlusions of the ipsilateral (but not contralateral) nostril, known to abolish gamma oscillations in the piriform cortex, strongly reduced vStr gamma power and the occurrence of transient gamma-band events. These results imply that local circuitry is not a major contributor to gamma oscillations in the vStr LFP and that piriform cortex is an important driver of gamma-band oscillations in the vStr and associated limbic areas.SIGNIFICANCE STATEMENT The ventral striatum (vStr) is an area of anatomical convergence in circuits underlying motivated behavior, but it remains unclear how its inputs from different sources interact. A major proposal about how neural circuits may switch dynamically between convergent inputs is through temporal organization reflected in local field potential (LFP) oscillations. Our results show that, in the rat, the mechanisms controlling gamma-band oscillations in the vStr LFP are primarily located in the in the adjacent piriform cortex rather than in the vStr itself, providing a novel interpretation of previous rodent work on gamma oscillations in the vStr and related circuits and an important consideration for future work seeking to use oscillations in these areas as biomarkers for behavioral and neurological disorders.

Changes of Ionic Concentrations During Seizure Transitions - A Modeling Study

Traditionally, it is considered that neuronal synchronization in epilepsy is caused by a chain reaction of synaptic excitation. However, it has been shown that synchronous epileptiform activity may also arise without synaptic transmission. In order to investigate the respective roles of synaptic interactions and nonsynaptic mechanisms in seizure transitions, we developed a computational model of hippocampal cells, involving the extracellular space, realistic dynamics of [Formula: see text], [Formula: see text], [Formula: see text] and [Formula: see text] ions, glial uptake and extracellular diffusion mechanisms. We show that the network behavior with fixed ionic concentrations may be quite different from the neurons' behavior when more detailed modeling of ionic dynamics is included. In particular, we show that in the extended model strong discharge of inhibitory interneurons may result in long lasting accumulation of extracellular [Formula: see text], which sustains the depolarization of the principal cells and causes their pathological discharges. This effect is not present in a reduced, purely synaptic network. These results point to the importance of nonsynaptic mechanisms in the transition to seizure.

Prediction of the Seizure Suppression Effect by Electrical Stimulation via a Computational Modeling Approach

In this paper, we identified factors that can affect seizure suppression via electrical stimulation by an integrative study based on experimental and computational approach. Preferentially, we analyzed the characteristics of seizure-like events (SLEs) using our previous in vitro experimental data. The results were analyzed in two groups classified according to the size of the effective region, in which the SLE was able to be completely suppressed by local stimulation. However, no significant differences were found between these two groups in terms of signal features or propagation characteristics (i.e., propagation delays, frequency spectrum, and phase synchrony). Thus, we further investigated important factors using a computational model that was capable of evaluating specific influences on effective region size. In the proposed model, signal transmission between neurons was based on two different mechanisms: synaptic transmission and the electrical field effect. We were able to induce SLEs having similar characteristics with differentially weighted adjustments for the two transmission methods in various noise environments. Although the SLEs had similar characteristics, their suppression effects differed. First of all, the suppression effect occurred only locally where directly received the stimulation effect in the high noise environment, but it occurred in the entire network in the low noise environment. Interestingly, in the same noise environment, the suppression effect was different depending on SLE propagation mechanism; only a local suppression effect was observed when the influence of the electrical field transmission was very weak, whereas a global effect was observed with a stronger electrical field effect. These results indicate that neuronal activities synchronized by a strong electrical field effect respond more sensitively to partial changes in the entire network. In addition, the proposed model was able to predict that stimulation of a seizure focus region is more effective for suppression. In conclusion, we confirmed the possibility of a computational model as a simulation tool to analyze the efficacy of deep brain stimulation (DBS) and investigated the key factors that determine the size of an effective region in seizure suppression via electrical stimulation.

Electrical Stimulation of the Human Cerebral Cortex by Extracranial Muscle Activity: Effect Quantification With Intracranial EEG and FEM Simulations

OBJECTIVE

Electric fields (EF) of approx. 0.2 V/m have been shown to be sufficiently strong to both modulate neuronal activity in the cerebral cortex and have measurable effects on cognitive performance. We hypothesized that the EF caused by the electrical activity of extracranial muscles during natural chewing may reach similar strength in the cerebral cortex and hence might act as an endogenous modality of brain stimulation. Here, we present first steps toward validating this hypothesis.

METHODS

Using a realistic volume conductor head model of an epilepsy patient having undergone intracranial electrode placement and utilizing simultaneous intracranial and extracranial electrical recordings during chewing, we derive predictions about the chewing-related cortical EF strength to be expected in healthy individuals.

RESULTS

We find that in the region of the temporal poles, the expected EF strength may reach amplitudes in the order of 0.1-1 V/m.

CONCLUSION

The cortical EF caused by natural chewing could be large enough to modulate ongoing neural activity in the cerebral cortex and influence cognitive performance.

SIGNIFICANCE

The present study lends first support for the assumption that extracranial muscle activity might represent an endogenous source of electrical brain stimulation. This offers a new potential explanation for the puzzling effects of gum chewing on cognition, which have been repeatedly reported in the literature.

Simulation of developing human neuronal cell networks

Background

Microelectrode array (MEA) is a widely used technique to study for example the functional properties of neuronal networks derived from human embryonic stem cells (hESC-NN). With hESC-NN, we can investigate the earliest developmental stages of neuronal network formation in the human brain.

Methods

In this paper, we propose an in silico model of maturating hESC-NNs based on a phenomenological model called INEX. We focus on simulations of the development of bursts in hESC-NNs, which are the main feature of neuronal activation patterns. The model was developed with data from developing hESC-NN recordings on MEAs which showed increase in the neuronal activity during the investigated six measurement time points in the experimental and simulated data.

Results

Our simulations suggest that the maturation process of hESC-NN, resulting in the formation of bursts, can be explained by the development of synapses. Moreover, spike and burst rate both decreased at the last measurement time point suggesting a pruning of synapses as the weak ones are removed.

Conclusions

To conclude, our model reflects the assumption that the interaction between excitatory and inhibitory neurons during the maturation of a neuronal network and the spontaneous emergence of bursts are due to increased connectivity caused by the forming of new synapses.

Propagating Neural Source Revealed by Doppler Shift of Population Spiking Frequency

Electrical activity in the brain during normal and abnormal function is associated with propagating waves of various speeds and directions. It is unclear how both fast and slow traveling waves with sometime opposite directions can coexist in the same neural tissue. By recording population spikes simultaneously throughout the unfolded rodent hippocampus with a penetrating microelectrode array, we have shown that fast and slow waves are causally related, so a slowly moving neural source generates fast-propagating waves at ∼0.12 m/s. The source of the fast population spikes is limited in space and moving at ∼0.016 m/s based on both direct and Doppler measurements among 36 different spiking trains among eight different hippocampi. The fact that the source is itself moving can account for the surprising direction reversal of the wave. Therefore, these results indicate that a small neural focus can move and that this phenomenon could explain the apparent wave reflection at tissue edges or multiple foci observed at different locations in neural tissue.

SIGNIFICANCE STATEMENT

The use of novel techniques with an unfolded hippocampus and penetrating microelectrode array to record and analyze neural activity has revealed the existence of a source of neural signals that propagates throughout the hippocampus. The source itself is electrically silent, but its location can be inferred by building isochrone maps of population spikes that the source generates. The movement of the source can also be tracked by observing the Doppler frequency shift of these spikes. These results have general implications for how neural signals are generated and propagated in the hippocampus; moreover, they have important implications for the understanding of seizure generation and foci localization.

Computational modeling of epileptiform activities in medial temporal lobe epilepsy combined with in vitro experiments

In this paper, we propose a comprehensive computational model that is able to reproduce three epileptiform activities. The model targets a hippocampal formation that is known to be an important lesion in medial temporal lobe epilepsy. It consists of four sub-networks consisting of excitatory and inhibitory neurons and well-known signal pathways, with consideration of propagation delay. The three epileptiform activities involve fast and slow interictal discharge and ictal discharge, and those activities can be induced in vitro by application of 4-Aminopyridine in entorhinal cortex combined hippocampal slices. We model the three epileptiform activities upon previously reported biological mechanisms and verify the simulation results by comparing them with in vitro experimental data obtained using a microelectrode array. We use the results of Granger causality analysis of recorded data to set input gains of signal pathways in the model, so that the compatibility between the computational and experimental models can be improved. The proposed model can be expanded to evaluate the suppression effect of epileptiform activities due to new treatment methods.

Neuronal coupling by endogenous electric fields: cable theory and applications to coincidence detector neurons in the auditory brain stem

The ongoing activity of neurons generates a spatially and time-varying field of extracellular voltage (Ve). This Ve field reflects population-level neural activity, but does it modulate neural dynamics and the function of neural circuits? We provide a cable theory framework to study how a bundle of model neurons generates Ve and how this Ve feeds back and influences membrane potential (Vm). We find that these "ephaptic interactions" are small but not negligible. The model neural population can generate Ve with millivolt-scale amplitude, and this Ve perturbs the Vm of "nearby" cables and effectively increases their electrotonic length. After using passive cable theory to systematically study ephaptic coupling, we explore a test case: the medial superior olive (MSO) in the auditory brain stem. The MSO is a possible locus of ephaptic interactions: sounds evoke large (millivolt scale)Vein vivo in this nucleus. The Ve response is thought to be generated by MSO neurons that perform a known neuronal computation with submillisecond temporal precision (coincidence detection to encode sound source location). Using a biophysically based model of MSO neurons, we find millivolt-scale ephaptic interactions consistent with the passive cable theory results. These subtle membrane potential perturbations induce changes in spike initiation threshold, spike time synchrony, and time difference sensitivity. These results suggest that ephaptic coupling may influence MSO function.

Can Neural Activity Propagate by Endogenous Electrical Field?

It is widely accepted that synaptic transmissions and gap junctions are the major governing mechanisms for signal traveling in the neural system. Yet, a group of neural waves, either physiological or pathological, share the same speed of ∼0.1 m/s without synaptic transmission or gap junctions, and this speed is not consistent with axonal conduction or ionic diffusion. The only explanation left is an electrical field effect. We tested the hypothesis that endogenous electric fields are sufficient to explain the propagation with in silico and in vitro experiments. Simulation results show that field effects alone can indeed mediate propagation across layers of neurons with speeds of 0.12 ± 0.09 m/s with pathological kinetics, and 0.11 ± 0.03 m/s with physiologic kinetics, both generating weak field amplitudes of ∼2-6 mV/mm. Further, the model predicted that propagation speed values are inversely proportional to the cell-to-cell distances, but do not significantly change with extracellular resistivity, membrane capacitance, or membrane resistance. In vitro recordings in mice hippocampi produced similar speeds (0.10 ± 0.03 m/s) and field amplitudes (2.5-5 mV/mm), and by applying a blocking field, the propagation speed was greatly reduced. Finally, osmolarity experiments confirmed the model's prediction that cell-to-cell distance inversely affects propagation speed. Together, these results show that despite their weak amplitude, electric fields can be solely responsible for spike propagation at ∼0.1 m/s. This phenomenon could be important to explain the slow propagation of epileptic activity and other normal propagations at similar speeds.

Low- and high-gamma oscillations deviate in opposite directions from zero-phase synchrony in the limbic corticostriatal loop

The loop structure of cortico-striatal anatomy in principle enables both descending (cortico-striatal) and ascending (striato-cortical) influences, but the factors that regulate the flow of information in these loops are not known. We report that low- and high-gamma oscillations (∼50 and ∼80 Hz, respectively) in the local field potential of freely moving rats are highly synchronous between the infralimbic region of the medial prefrontal cortex (mPFC) and the ventral striatum (vStr). Strikingly, high-gamma oscillations in mPFC preceded those in vStr, whereas low-gamma oscillations in mPFC lagged those in vStr, with short (∼1 ms) time lags. These systematic deviations from zero-phase synchrony were consistent across measures based on amplitude cross-correlation and phase slopes and were robustly maintained between behavioral states and different individual subjects. Furthermore, low- and high-gamma oscillations were associated with distinct ensemble spiking patterns in vStr, even when controlling for overt behavioral differences and slow changes in neural activity. These results imply that neural activity in vStr and mPFC is tightly coupled at the gamma timescale and raise the intriguing possibility that frequency-specific deviations from this coupling may signal transient leader-follower switches.

Slow Spatial Recruitment of Neocortex during Secondarily Generalized Seizures and Its Relation to Surgical Outcome

Understanding the spatiotemporal dynamics of brain activity is crucial for inferring the underlying synaptic and nonsynaptic mechanisms of brain dysfunction. Focal seizures with secondary generalization are traditionally considered to begin in a limited spatial region and spread to connected areas, which can include both pathological and normal brain tissue. The mechanisms underlying this spread are important to our understanding of seizures and to improve therapies for surgical intervention. Here we study the properties of seizure recruitment-how electrical brain activity transitions to large voltage fluctuations characteristic of spike-and-wave seizures. We do so using invasive subdural electrode arrays from a population of 16 patients with pharmacoresistant epilepsy. We find an average delay of ∼30 s for a broad area of cortex (8 × 8 cm) to be recruited into the seizure, at an estimated speed of ∼4 mm/s. The spatiotemporal characteristics of recruitment reveal two categories of patients: one in which seizure recruitment of neighboring cortical regions follows a spatially organized pattern consistent from seizure to seizure, and a second group without consistent spatial organization of activity during recruitment. The consistent, organized recruitment correlates with a more regular, compared with small-world, connectivity pattern in simulation and successful surgical treatment of epilepsy. We propose that an improved understanding of how the seizure recruits brain regions into large amplitude voltage fluctuations provides novel information to improve surgical treatment of epilepsy and highlights the slow spread of massive local activity across a vast extent of cortex during seizure.

Seizure reduction through interneuron-mediated entrainment using low frequency optical stimulation

Low frequency electrical stimulation (LFS) can reduce neural excitability and suppress seizures in animals and patients with epilepsy. However the therapeutic outcome could benefit from the determination of the cell types involved in seizure suppression. We used optogenetic techniques to investigate the role of interneurons in LFS (1Hz) in the epileptogenic hippocampus. Optical low frequency stimulation (oLFS) was first used to activate the cation channel channelrhodopsin-2 (ChR2) in the Thy1-ChR2 transgenic mouse that expresses ChR2 in both excitatory and inhibitory neurons. We found that oLFS could effectively reduce epileptiform activity in the hippocampus through the activation of GAD-expressing hippocampal interneurons. This was confirmed using the VGAT-ChR2 transgenic mouse, allowing for selective optical activation of only GABA interneurons. Activating hippocampal interneurons through oLFS was found to cause entrainment of neural activity similar to electrical stimulation, but through a GABAA-mediated mechanism. These results confirm the robustness of the LFS paradigm and indicate that GABA interneurons play an unexpected role of shaping inter-ictal activity to decrease neural excitability in the hippocampus.