Layer- and Cell-Specific Recruitment Dynamics during Epileptic Seizures In Vivo

Preferential superficial cortical layer activation during seizure propagation

OBJECTIVE

Focal cortical seizures travel long distances from the onset zone, but the long-distance propagation pathways are uncertain. In vitro and in vivo imaging techniques have investigated the local spread of seizures but did not elucidate long-distance spread. Furthermore, classical studies in slices suggested seizure spread locally along deep cortical layers, whereas more recent in vivo imaging studies posit a role for superficial cortical layers in local spread.

METHODS

We imaged seizure-activated neurons using activity reporter mice and measured local field potentials (LFPs) using microelectrode arrays to map cortical seizure propagation in awake mice.

RESULTS

Frontal lobe onset seizures activate more neurons in superficial layers 2-3 than deep layers 5-6 throughout the cortex. LFP recordings demonstrate that seizures spread faster through the superficial than deep layers over long cortical distances of 3.5 mm. We also show that monosynaptically connected long-distance neurons are in the seizure circuit.

SIGNIFICANCE

We propose that long-distance cortical seizure spread occurs preferentially via synaptically connected superficial cortical neurons.

NeoCoMM: A neocortical neuroinspired computational model for the reconstruction and simulation of epileptiform events

BACKGROUND

Understanding the pathophysiological dynamics that underline Interictal Epileptiform Events (IEEs) such as epileptic spikes, spike-and-waves or High-Frequency Oscillations (HFOs) is of major importance in the context of neocortical refractory epilepsy, as it paves the way for the development of novel therapies. Typically, these events are detected in Local Field Potential (LFP) recordings obtained through depth electrodes during pre-surgical investigations. Although essential, the underlying pathophysiological mechanisms for the generation of these epileptic neuromarkers remain unclear. The aim of this paper is to propose a novel neurophysiologically relevant reconstruction of the neocortical microcircuitry in the context of epilepsy. This reconstruction intends to facilitate the analysis of a comprehensive set of parameters encompassing physiological, morphological, and biophysical aspects that directly impact the generation and recording of different IEEs.

METHOD

a novel microscale computational model of an epileptic neocortical column was introduced. This model incorporates the intricate multilayered structure of the cortex and allows for the simulation of realistic interictal epileptic signals. The proposed model was validated through comparisons with real IEEs recorded using intracranial stereo-electroencephalography (SEEG) signals from both humans and animals. Using the model, the user can recreate epileptiform patterns observed in different species (human, rodent, and mouse) and study the intracellular activity associated with these patterns.

RESULTS

Our model allowed us to unravel the relationship between glutamatergic and GABAergic synaptic transmission of the epileptic neural network and the type of generated IEE. Moreover, sensitivity analyses allowed for the exploration of the pathophysiological parameters responsible for the transitions between these events. Finally, the presented modeling framework also provides an Electrode Tissue Model (ETI) that adds realism to the simulated signals and offers the possibility of studying their sensitivity to the electrode characteristics.

CONCLUSION

The model (NeoCoMM) presented in this work can be of great use in different applications since it offers an in silico framework for sensitivity analysis and hypothesis testing. It can also be used as a starting point for more complex studies.

Epilepsy insights revealed by intravital functional optical imaging

Despite an abundance of pharmacologic and surgical epilepsy treatments, there remain millions of patients suffering from poorly controlled seizures. One approach to closing this treatment gap may be found through a deeper mechanistic understanding of the network alterations that underly this aberrant activity. Functional optical imaging in vertebrate models provides powerful advantages to this end, enabling the spatiotemporal acquisition of individual neuron activity patterns across multiple seizures. This coupled with the advent of genetically encoded indicators, be them for specific ions, neurotransmitters or voltage, grants researchers unparalleled access to the intact nervous system. Here, we will review how in vivo functional optical imaging in various vertebrate seizure models has advanced our knowledge of seizure dynamics, principally seizure initiation, propagation and termination.

Differential cortical layer engagement during seizure initiation and spread in humans

Despite decades of research, we still do not understand how spontaneous human seizures start and spread - especially at the level of neuronal microcircuits. In this study, we used laminar arrays of micro-electrodes to simultaneously record the local field potentials and multi-unit neural activities across the six layers of the neocortex during focal seizures in humans. We found that, within the ictal onset zone, the discharges generated during a seizure consisted of current sinks and sources only within the infra-granular and granular layers. Outside of the seizure onset zone, ictal discharges reflected current flow in the supra-granular layers. Interestingly, these patterns of current flow evolved during the course of the seizure - especially outside the seizure onset zone where superficial sinks and sources extended into the deeper layers. Based on these observations, a framework describing cortical-cortical dynamics of seizures is proposed with implications for seizure localization, surgical targeting, and neuromodulation techniques to block the generation and propagation of seizures.

Early and widespread engagement of the cerebellum during hippocampal epileptiform activity Format: Brief Communication

Despite research illustrating the cerebellum may be a critical circuit element in the epilepsies, remarkably little is known about cerebellar engagement during seizures. We therefore implemented a novel method for repeated imaging of the cerebellum in awake, chronically epileptic animals. We found widespread changes in cerebellar calcium signals during behavioral seizures and during hippocampal seizures that remained electrographic only, arguing against cerebellar modulation simply reflecting motor components. Moreover, even brief interictal spikes produced widespread alterations in cerebellar activity. Changes were noted in the anterior and posterior cerebellum, along the midline, and both ipsilaterally and contralaterally to the seizure focus. Remarkably, changes in the cerebellum also occurred prior to any noticeable change in the hippocampal electrographic recordings, suggesting a special relationship between the cerebellum and hippocampal epileptiform activity. Together these results underscore the importance of the cerebellum in epilepsy, warranting a more consistent consideration of the cerebellum when evaluating epilepsy patients.

One Ring to Bind Them: The Annulus of GABAergic Inhibitory Restraint Fades at Seizure Emergence

Extracellular Glutamate and GABA Transients at the Transition From Interictal Spiking to Seizures Shimoda Y, Leite M, Graham RT, Marvin JS, Hasseman J, Kolb I, Looger LL, Magloire V, Kullmann DM. Brain . 2023: awad336. doi:10.1093/brain/awad336 Focal epilepsy is associated with intermittent brief population discharges (interictal spikes), which resemble sentinel spikes that often occur at the onset of seizures. Why interictal spikes self-terminate whilst seizures persist and propagate is incompletely understood. We used fluorescent glutamate and GABA sensors in an awake rodent model of neocortical seizures to resolve the spatiotemporal evolution of both neurotransmitters in the extracellular space. Interictal spikes were accompanied by brief glutamate transients which were maximal at the initiation site and rapidly propagated centrifugally. GABA transients lasted longer than glutamate transients and were maximal ∼1.5 mm from the focus where they propagated centripetally. Prior to seizure initiation GABA transients were attenuated, whilst glutamate transients increased, consistent with a progressive failure of local inhibitory restraint. As seizures increased in frequency, there was a gradual increase in the spatial extent of spike-associated glutamate transients associated with interictal spikes. Neurotransmitter imaging thus reveals a progressive collapse of an annulus of feed-forward GABA release, allowing seizures to escape from local inhibitory restraint.

A debate on the neuronal origin of focal seizures

A critical question regarding how focal seizures start is whether we can identify particular cell classes that drive the pathological process. This was the topic for debate at the recent International Conference for Technology and Analysis of Seizures (ICTALS) meeting (July 2022, Bern, CH) that we summarize here. The debate has been fueled in recent times by the introduction of powerful new ways to manipulate subpopulations of cells in relative isolation, mostly using optogenetics. The motivation for resolving the debate is to identify novel targets for therapeutic interventions through a deeper understanding of the etiology of seizures.

Seizure Event Detection Using Intravital Two-Photon Calcium Imaging Data

Significance

Genetic cellular calcium imaging has emerged as a powerful tool to investigate how different types of neurons interact at the microcircuit level to produce seizure activity, with newfound potential to understand epilepsy. Although many methods exist to measure seizure-related activity in traditional electrophysiology, few yet exist for calcium imaging.

Aim

To demonstrate an automated algorithmic framework to detect seizure-related events using calcium imaging - including the detection of pre-ictal spike events, propagation of the seizure wavefront, and terminal spreading waves for both population-level activity and that of individual cells.

Approach

We developed an algorithm for precise recruitment detection of population and individual cells during seizure-associated events, which broadly leverages averaged population activity and high-magnitude slope features to detect single-cell pre-ictal spike and seizure recruitment. We applied this method to data recorded using awake in vivo two-photon calcium imaging during pentylenetetrazol induced seizures in mice.

Results

We demonstrate that our detected recruitment times are concordant with visually identified labels provided by an expert reviewer and are sufficiently accurate to model the spatiotemporal progression of seizure-associated traveling waves.

Conclusions

Our algorithm enables accurate cell recruitment detection and will serve as a useful tool for researchers investigating seizure dynamics using calcium imaging.

Microscale Neuronal Activity Collectively Drives Chaotic and Inflexible Dynamics at the Macroscale in Seizures

Neuronal activity propagates through the network during seizures, engaging brain dynamics at multiple scales. Such propagating events can be described through the avalanches framework, which can relate spatiotemporal activity at the microscale with global network properties. Interestingly, propagating avalanches in healthy networks are indicative of critical dynamics, where the network is organized to a phase transition, which optimizes certain computational properties. Some have hypothesized that the pathologic brain dynamics of epileptic seizures are an emergent property of microscale neuronal networks collectively driving the brain away from criticality. Demonstrating this would provide a unifying mechanism linking microscale spatiotemporal activity with emergent brain dysfunction during seizures. Here, we investigated the effect of drug-induced seizures on critical avalanche dynamics, using in vivo whole-brain two-photon imaging of GCaMP6s larval zebrafish (males and females) at single neuron resolution. We demonstrate that single neuron activity across the whole brain exhibits a loss of critical statistics during seizures, suggesting that microscale activity collectively drives macroscale dynamics away from criticality. We also construct spiking network models at the scale of the larval zebrafish brain, to demonstrate that only densely connected networks can drive brain-wide seizure dynamics away from criticality. Importantly, such dense networks also disrupt the optimal computational capacities of critical networks, leading to chaotic dynamics, impaired network response properties and sticky states, thus helping to explain functional impairments during seizures. This study bridges the gap between microscale neuronal activity and emergent macroscale dynamics and cognitive dysfunction during seizures.SIGNIFICANCE STATEMENT Epileptic seizures are debilitating and impair normal brain function. It is unclear how the coordinated behavior of neurons collectively impairs brain function during seizures. To investigate this we perform fluorescence microscopy in larval zebrafish, which allows for the recording of whole-brain activity at single-neuron resolution. Using techniques from physics, we show that neuronal activity during seizures drives the brain away from criticality, a regime that enables both high and low activity states, into an inflexible regime that drives high activity states. Importantly, this change is caused by more connections in the network, which we show disrupts the ability of the brain to respond appropriately to its environment. Therefore, we identify key neuronal network mechanisms driving seizures and concurrent cognitive dysfunction.

SCN1A channelopathies: Navigating from genotype to neural circuit dysfunction

The SCN1A gene is strongly associated with epilepsy and plays a central role for supporting cortical excitation-inhibition balance through the expression of Na_V1.1 within inhibitory interneurons. The phenotype of SCN1A disorders has been conceptualized as driven primarily by impaired interneuron function that predisposes to disinhibition and cortical hyperexcitability. However, recent studies have identified SCN1A gain-of-function variants associated with epilepsy, and the presence of cellular and synaptic changes in mouse models that point toward homeostatic adaptations and complex network remodeling. These findings highlight the need to understand microcircuit-scale dysfunction in SCN1A disorders to contextualize genetic and cellular disease mechanisms. Targeting the restoration of microcircuit properties may be a fruitful strategy for the development of novel therapies.

Dynamic Fluorescence Imaging in Experimental Epilepsy

The ability to develop effective new treatments for epilepsy may depend on improved understanding of seizure pathophysiology, about which many questions remain. Dynamic fluorescence imaging of activity at single-neuron resolution with fluorescent indicators in experimental model systems in vivo has revolutionized basic neuroscience and has the potential to do so for epilepsy research as well. Here, we review salient issues as they pertain to experimental imaging in basic epilepsy research, including commonly used imaging technologies, data processing and analysis, interpretation of results, and selected examples of how imaging-based approaches have revealed new insight into mechanisms of seizures and epilepsy.

Signal processing and computational modeling for interpretation of SEEG-recorded interictal epileptiform discharges in epileptogenic and non-epileptogenic zones

OBJECTIVE

In partial epilepsies, interictal epileptiform discharges (IEDs) are paroxysmal events observed in epileptogenic and non-epileptogenic zones. IEDs' generation and recurrence are subject to different hypotheses: they appear through glutamatergic and GABAergic processes; they may trigger seizures or prevent seizure propagation. This paper focuses on a specific class of IEDs, spike-waves (SWs), characterized by a short-duration spike followed by a longer duration wave, both of the same polarity. Signal analysis and neurophysiological mathematical models are used to interpret puzzling IED generation.

APPROACH

Interictal activity was recorded by intracranial stereo-electroencephalography (SEEG) electrodes in five different patients. SEEG experts identified the epileptic and non-epileptic zones in which IEDs were detected. After quantifying spatial and temporal features of the detected IEDs, the most significant features for classifying epileptic and non-epileptic zones were determined. A neurophysiologically-plausible mathematical model was then introduced to simulate the IEDs and understand the underlying differences observed in epileptic and non-epileptic zone IEDs.

MAIN RESULTS

Two classes of SWs were identified according to subtle differences in morphology and timing of the spike and wave component. Results showed that type-1 SWs were generated in epileptogenic regions also involved at seizure onset, while type-2 SWs were produced in the propagation or non-involved areas. The modeling study indicated that synaptic kinetics, cortical organization, and network interactions determined the morphology of the simulated SEEG signals. Modeling results suggested that the IED morphologies were linked to the degree of preserved inhibition.

SIGNIFICANCE

This work contributes to the understanding of different mechanisms generating IEDs in epileptic networks. The combination of signal analysis and computational models provides an efficient framework for exploring IEDs in partial epilepsies and classifying epileptogenic and non-epileptogenic zones.

Temporal and morphological characteristics of high-frequency oscillations in an acute in vivo model of epilepsy

Approximately 30% of patients with epilepsy do not respond to anti-epileptogenic drugs. Surgical removal of the epileptogenic zone (EZ), the brain regions where the seizures originate and spread, can be a possible therapy for these patients, but localizing the EZ is challenging due to a variety of clinical factors. High-frequency oscillations (HFOs) in intracranial electroencephalography (EEG) are a promising biomarker of the EZ, but it is currently unknown whether HFO rates and HFO morphology modulate as pathological brain networks evolve in a way that gives rise to seizures. To address this question, we assessed the temporal evolution of the duration of HFO events, amplitude of HFO events, and rates of HFOs per minute. HFO events were quantified using the 4AP in vivo rodent model of epilepsy, inducing seizures in two different brain areas. We found that the duration and amplitude of HFO events were significantly increased for the cortex model when compared to the hippocampus model. Additionally, the duration and amplitude increased significantly between baseline and pre-ictal HFOs in both models. On the other hand, the two models did not display a consistent increasing or decreasing trend in amplitude, duration or rate when comparing ictal and postictal intervals. Clinical Relevance- We assessed the amplitude, duration, and rate of HFOs in two acute in vivo rodent models of epilepsy. The significant modulation of HFO morphology from baseline to pre-ictal periods suggests that these features may be a robust biomarker for pathological tissue involved in epileptogenesis. Moreover, the differences in HFO morphology observed between cortex and hippocampus animal models possibly indicate that different structural network characteristics of the EZ cause this modulation. In all, we found that HFO features modulate significantly with the onset of seizures, further highlighting the need to consider of HFO morphology in EZ-localizing studies.

Multimodal Characterization of Seizures in Zebrafish Larvae

Epilepsy accounts for a significant proportion of the world's disease burden. Indeed, many research efforts are produced both to investigate the basic mechanism ruling its genesis and to find more effective therapies. In this framework, the use of zebrafish larvae, owing to their peculiar features, offers a great opportunity. Here, we employ transgenic zebrafish larvae expressing GCaMP6s in all neurons to characterize functional alterations occurring during seizures induced by pentylenetetrazole. Using a custom two-photon light-sheet microscope, we perform fast volumetric functional imaging of the entire larval brain, investigating how different brain regions contribute to seizure onset and propagation. Moreover, employing a custom behavioral tracking system, we outline the progressive alteration of larval swim kinematics, resulting from different grades of seizures. Collectively, our results show that the epileptic larval brain undergoes transitions between diverse neuronal activity regimes. Moreover, we observe that different brain regions are progressively recruited into the generation of seizures of diverse severity. We demonstrate that midbrain regions exhibit highest susceptibility to the convulsant effects and that, during periods preceding abrupt hypersynchronous paroxysmal activity, they show a consistent increase in functional connectivity. These aspects, coupled with the hub-like role that these regions exert, represent important cues in their identification as epileptogenic hubs.

Input rate encoding and gain control in dendrites of neocortical pyramidal neurons

Elucidating how neurons encode network activity is essential to understanding how the brain processes information. Neocortical pyramidal cells receive excitatory input onto spines distributed along dendritic branches. Local dendritic branch nonlinearities can boost the response to spatially clustered and synchronous input, but how this translates into the integration of patterns of ongoing activity remains unclear. To examine dendritic integration under naturalistic stimulus regimes, we use two-photon glutamate uncaging to repeatedly activate multiple dendritic spines at random intervals. In the proximal dendrites of two populations of layer 5 pyramidal neurons in the mouse motor cortex, spatially restricted synchrony is not a prerequisite for dendritic boosting. Branches encode afferent inputs with distinct rate sensitivities depending upon cell and branch type. Thus, inputs distributed along a dendritic branch can recruit supralinear boosting and the window of this nonlinearity may provide a mechanism by which dendrites can preferentially amplify slow-frequency network oscillations.

State transitions through inhibitory interneurons in a cortical network model

Inhibitory interneurons shape the spiking characteristics and computational properties of cortical networks. Interneuron subtypes can precisely regulate cortical function but the roles of interneuron subtypes for promoting different regimes of cortical activity remains unclear. Therefore, we investigated the impact of fast spiking and non-fast spiking interneuron subtypes on cortical activity using a network model with connectivity and synaptic properties constrained by experimental data. We found that network properties were more sensitive to modulation of the fast spiking population, with reductions of fast spiking excitability generating strong spike correlations and network oscillations. Paradoxically, reduced fast spiking excitability produced a reduction of global excitation-inhibition balance and features of an inhibition stabilised network, in which firing rates were driven by the activity of excitatory neurons within the network. Further analysis revealed that the synaptic interactions and biophysical features associated with fast spiking interneurons, in particular their rapid intrinsic response properties and short synaptic latency, enabled this state transition by enhancing gain within the excitatory population. Therefore, fast spiking interneurons may be uniquely positioned to control the strength of recurrent excitatory connectivity and the transition to an inhibition stabilised regime. Overall, our results suggest that interneuron subtypes can exert selective control over excitatory gain allowing for differential modulation of global network state.

Inhibitory interneurons comprise a significant proportion of all cortical neurons and play a crucial role in sustaining normal neural activity in the brain. Although it is well established that there exist distinct subtypes of interneurons, the impact of different interneuron subtypes upon cortical function remains unclear. In this work, we explore the role of interneuron subtypes for modulating neural activity using a network model containing two of the most common interneuron subtypes. We find that one interneuron subtype, known as fast spiking interneurons, preferentially control the strength of activity between excitatory neurons to regulate changes in network state. These findings suggest that interneuron subtypes may selectively modulate cortical activity to promote different computational capabilities.

Two-photon calcium imaging of seizures in awake, head-fixed mice

Epilepsy is a severe neurological disorder defined by spontaneous seizures. Current treatment options fail in a large proportion of patients, while questions as to the basic mechanisms of seizure initiation and propagation remain. Advances in imaging of seizures in experimental model systems could lead to a better understanding of mechanisms of seizures and epilepsy. Recent studies have used two-photon calcium imaging (2 P imaging) in awake, behaving mice in head-fixed preparations to image seizures in vivo at high speed and cellular-level resolution to identify key seizure-related cell classes. Here, we discuss such advances and present 2 P imaging data of excitatory neurons and defined subsets of cerebral cortex GABAergic inhibitory interneurons during naturalistic seizures in a mouse model of Dravet syndrome (Scn1a+/- mice) along with other behavioral measures. Results demonstrate differential recruitment of discrete interneuron subclasses, which could inform mechanisms of seizure generation and propagation in Dravet syndrome and other epilepsies.

Intrinsic Disruption of the M1 Cortical Network in a Mouse Model of Parkinson's Disease

BACKGROUND

Parkinson's disease (PD) disrupts motor performance by affecting the basal ganglia system. Yet, despite the critical position of the primary motor cortex in linking basal ganglia computations with motor performance, its contribution to motor disability in PD is largely unknown. The objective of this study was to characterize the role of the primary motor cortex in PD-related motor disability.

METHODS

Two-photon calcium imaging and optogenetic stimulation of primary motor cortex neurons was done during performance of a dexterous reach-to-grasp motor task in control and 6-hydroxydopamine-induced PD mice.

RESULTS

Experimental PD disrupted performance of the reach-to-grasp motor task and especially initiation of the task, which was partially restored by optogenetic activation of the primary motor cortex. Two-photon calcium imaging during task performance revealed experimental-PD affected the primary motor cortex in a cell-type-specific manner. It suppressed activation of output layer 5 pyramidal tract neurons, with greater effects on freeze versus nonfreeze trials. In contrast, it did not attenuate the initial movement-related activation response of layer 2/3 pyramidal neurons while diminishing the late inhibitory phase of their response. At the network level, experimental PD disrupted movement-related population dynamics of the layer 5 pyramidal tract network while almost not affecting the dynamics of the layer 2/3 neuronal population. It also disrupted short- and long-term robustness and stability of the pyramidal tract subnetwork, with reduced intertrial temporal accuracy and diminished reproducibility of motor parameter encoding and temporal recruitment of the output pyramidal tract neurons over repeated daily sessions.

CONCLUSIONS

Experimental PD disrupts both external driving and intrinsic properties of the primary motor cortex. Motor disability in experimental PD results primarily from the inability to generate robust and stable output motor sequences in the parkinsonian primary motor cortex output layer 5 pyramidal tract subnetwork. © 2021 International Parkinson and Movement Disorder Society.

Multimodal in vivo recording using transparent graphene microelectrodes illuminates spatiotemporal seizure dynamics at the microscale

Neurological disorders such as epilepsy arise from disrupted brain networks. Our capacity to treat these disorders is limited by our inability to map these networks at sufficient temporal and spatial scales to target interventions. Current best techniques either sample broad areas at low temporal resolution (e.g. calcium imaging) or record from discrete regions at high temporal resolution (e.g. electrophysiology). This limitation hampers our ability to understand and intervene in aberrations of network dynamics. Here we present a technique to map the onset and spatiotemporal spread of acute epileptic seizures in vivo by simultaneously recording high bandwidth microelectrocorticography and calcium fluorescence using transparent graphene microelectrode arrays. We integrate dynamic data features from both modalities using non-negative matrix factorization to identify sequential spatiotemporal patterns of seizure onset and evolution, revealing how the temporal progression of ictal electrophysiology is linked to the spatial evolution of the recruited seizure core. This integrated analysis of multimodal data reveals otherwise hidden state transitions in the spatial and temporal progression of acute seizures. The techniques demonstrated here may enable future targeted therapeutic interventions and novel spatially embedded models of local circuit dynamics during seizure onset and evolution.

Network Fragility for Seizure Genesis in an Acute in vivo Model of Epilepsy

Epilepsy affects over 50 million people worldwide and 30% of patients' seizures are medically refractory. The process of localizing and removing the epileptogenic zone is error-prone and ill-posed in part because we do not understand how epilepsy manifests. It has recently been proposed that the epileptic cortex is fragile in the sense that seizures manifest through small perturbations in the synaptic connections that render the entire cortical network unstable. If the fragility of the cortical network could be computed over a period in which seizure genesis occurs, then it might elucidate network mechanisms correlated to the epileptogenic zone. In this study, we used local field potentials (LFP) from neocortex by implementing an acute model of epilepsy in mice. These recordings were used to develop a dynamical network model that quantifies the fragility of the nodes from LFP epochs of baseline activity, preictal and ictal states. Fragility was quantified by the generation of a linear time-varying model to which we then applied a perturbation to determine the sensitivity of nodes in the network. Spatiotemporal fragility maps showed clear quantifiable changes in the epileptogenic network's properties throughout different states of seizure genesis. We quantified this difference over a baseline, preictal and ictal periods to show that network fragility is modulated in the manifestation of epilepsy.

2P or not 2P: The Question of Seizure Initiation

[Box: see text].