Indirect Effects of Halorhodopsin Activation: Potassium Redistribution, Nonspecific Inhibition, and Spreading Depolarization.

The movement of ions in and out of neurons can exert significant effects on neighboring cells. Here we report several experimentally important consequences of activation of the optogenetic chloride pump, halorhodopsin. We recorded extracellular K+ concentration ([K+]extra) in neocortical brain slices prepared from young adult mice (both sexes) which express halorhodopsin in pyramidal cells. Strong halorhodopsin activation induced a pronounced drop in [K+]extra that persisted for the duration of illumination. Pharmacological blockade of K+ channels reduced the amplitude of this drop, indicating that it represents K+ redistribution into cells during the period of hyperpolarization. Halorhodopsin thus drives the inward movement of both Cl- directly, and K+ secondarily. When the illumination period ended, a rebound surge in extracellular [K+] developed over tens of seconds, partly reflecting the previous inward redistribution of K+, but additionally driven by clearance of Cl- coupled to K+ by the potassium-chloride cotransporter, KCC2. The drop in [K+]extra during light activation leads to a small (2-3 mV) hyperpolarization also of other cells that do not express halorhodopsin. Its activation therefore has both direct and indirect inhibitory effects. Finally, we show that persistent strong activation of halorhodopsin causes cortical spreading depolarizations (CSDs), both in vitro and in vivo This novel means of triggering CSDs is unusual, in that the events can arise during the actual period of illumination, when neurons are being hyperpolarized and [K+]extra is low. We suggest that this fundamentally different experimental model of CSDs will open up new avenues of research to explain how they occur naturally.SIGNIFICANCE STATEMENT Halorhodopsin is a light-activated electrogenic chloride pump, which has been widely used to inhibit neurons optogenetically. Here, we demonstrate three previously unrecognized consequences of its use: (1) intense activation leads to secondary movement of K+ ions into the cells; (2) the resultant drop in extracellular [K+] reduces excitability also in other, nonexpressing cells; and (3) intense persistent halorhodopsin activation can trigger cortical spreading depolarization (CSD). Halorhodopsin-induced CSDs can occur when neurons are hyperpolarized and extracellular [K+] is low. This contrasts with the most widely used experimental models that trigger CSDs with high [K+]. Both models, however, are consistent with the hypothesis that CSDs arise following net inward ionic movement into the principal neuron population.

Medial septal GABAergic neurons reduce seizure duration upon optogenetic closed-loop stimulation.

Seizures can emerge from multiple or large foci in temporal lobe epilepsy, complicating focally targeted strategies such as surgical resection or the modulation of the activity of specific hippocampal neuronal populations through genetic or optogenetic techniques. Here, we evaluate a strategy in which optogenetic activation of medial septal GABAergic neurons, which provide extensive projections throughout the hippocampus, is used to control seizures. We utilized the chronic intrahippocampal kainate mouse model of temporal lobe epilepsy, which results in spontaneous seizures and as is often the case in human patients, presents with hippocampal sclerosis. Medial septal GABAergic neuron populations were immunohistochemically labelled and were not reduced in epileptic conditions. Genetic labelling with mRuby of medial septal GABAergic neuron synaptic puncta and imaging across the rostral to caudal extent of the hippocampus, also indicated an unchanged number of putative synapses in epilepsy. Furthermore, optogenetic stimulation of medial septal GABAergic neurons consistently modulated oscillations across multiple hippocampal locations in control and epileptic conditions. Finally, wireless optogenetic stimulation of medial septal GABAergic neurons, upon electrographic detection of spontaneous hippocampal seizures, resulted in reduced seizure durations. We propose medial septal GABAergic neurons as a novel target for optogenetic control of seizures in temporal lobe epilepsy.

Excitatory GABAergic signalling is associated with benzodiazepine resistance in status epilepticus.

Status epilepticus is defined as a state of unrelenting seizure activity. Generalized convulsive status epilepticus is associated with a rapidly rising mortality rate, and thus constitutes a medical emergency. Benzodiazepines, which act as positive modulators of chloride (Cl-) permeable GABAA receptors, are indicated as first-line treatment, but this is ineffective in many cases. We found that 48% of children presenting with status epilepticus were unresponsive to benzodiazepine treatment, and critically, that the duration of status epilepticus at the time of treatment is an important predictor of non-responsiveness. We therefore investigated the cellular mechanisms that underlie acquired benzodiazepine resistance, using rodent organotypic and acute brain slices. Removing Mg2+ ions leads to an evolving pattern of epileptiform activity, and eventually to a persistent state of repetitive discharges that strongly resembles clinical EEG recordings of status epilepticus. We found that diazepam loses its antiseizure efficacy and conversely exacerbates epileptiform activity during this stage of status epilepticus-like activity. Interestingly, a low concentration of the barbiturate phenobarbital had a similar exacerbating effect on status epilepticus-like activity, while a high concentration of phenobarbital was effective at reducing or preventing epileptiform discharges. We then show that the persistent status epilepticus-like activity is associated with a reduction in GABAA receptor conductance and Cl- extrusion capability. We explored the effect on intraneuronal Cl- using both gramicidin, perforated-patch clamp recordings and Cl- imaging. This showed that during status epilepticus-like activity, reduced Cl- extrusion capacity was further exacerbated by activity-dependent Cl- loading, resulting in a persistently high intraneuronal Cl-. Consistent with these results, we found that optogenetic stimulation of GABAergic interneurons in the status epilepticus-like state, actually enhanced epileptiform activity in a GABAAR dependent manner. Together our findings describe a novel potential mechanism underlying benzodiazepine-resistant status epilepticus, with relevance to how this life-threatening condition should be managed in the clinic.

Region-specific differences and areal interactions underlying transitions in epileptiform activity.

KEY POINTS: Local neocortical and hippocampal territories show different and sterotypical patterns of acutely evolving, epileptiform activity. Neocortical and entorhinal networks show tonic-clonic-like events, but the main hippocampal territories do not, unless it is relayed from the other areas. Transitions in the pattern of locally recorded epileptiform activity can be indicative of a shift in the source of pathological activity, and may spread through both synaptic and non-synaptic means. Hippocampal epileptiform activity is promoted by 4-aminopyridine and inhibited by GABAB receptor agonists, and appears far more sensitive to these drugs than neocortical activity. These signature features of local epileptiform activity can provide useful insight into the primary source of ictal activity, aiding both experimental and clinical investigation. ABSTRACT: Understanding the nature of epileptic state transitions remains a major goal for epilepsy research. Simple in vitro models offer unique experimental opportunities that we exploit to show that such transitions can arise from shifts in the ictal source of the activity. These transitions reflect the fact that cortical territories differ both in the type of epileptiform activity they can sustain and in their susceptibility to drug manipulation. In the zero-Mg2+ model, the earliest epileptiform activity is restricted to neocortical and entorhinal networks. Hippocampal bursting only starts much later, and triggers a marked transition in neo-/entorhinal cortical activity. Thereafter, the hippocampal activity acts as a pacemaker, entraining the other territories to their discharge pattern. This entrainment persists following transection of the major axonal pathways between hippocampus and cortex, indicating that it can be mediated through a non-synaptic route. Neuronal discharges are associated with large rises in extracellular [K+ ], but we show that these are very localized, and therefore are not the means of entraining distant cortical areas. We conclude instead that the entrainment occurs through weak field effects distant from the pacemaker, but which are highly effective at recruiting other brain territories that are already hyperexcitable. The hippocampal epileptiform activity appears unusually susceptible to drugs that impact on K+ conductances. These findings demonstrate that the local circuitry gives rise to stereotypical epileptic activity patterns, but these are also influenced by both synaptic and non-synaptic long-range effects. Our results have important implications for our understanding of epileptic propagation and anti-epileptic drug action.

Feedforward inhibition ahead of ictal wavefronts is provided by both parvalbumin- and somatostatin-expressing interneurons.

KEY POINTS: There is a rapid interneuronal response to focal activity in cortex, which restrains laterally propagating activity, including spreading epileptiform activity. The interneuronal response involves intense activation of both parvalbumin- and somatostatin-expressing interneurons. Interneuronal bursting is time-locked to glutamatergic barrages in the pre-ictal period. Ca2+ imaging using conditional expression of GCaMP6f provides an accurate readout of the evolving firing patterns in both types of interneuron. The activation profiles of the two interneuronal classes are temporally offset, with the parvalbumin population being activated first, and typically, at higher rates. ABSTRACT: Previous work has described powerful restraints on laterally spreading activity in cortical networks, arising from a rapid feedforward interneuronal response to focal activity. This response is particularly prominent ahead of an ictal wavefront. Parvalbumin-positive interneurons are considered to be critically involved in this feedforward inhibition, but it is not known what role, if any, is provided by somatostatin-expressing interneurons, which target the distal dendrites of pyramidal cells. We used a combination of electrophysiology and cell class-specific Ca2+ imaging in mouse brain slices bathed in 0 Mg2+ medium to characterize the activity profiles of pyramidal cells and parvalbumin- and somatostatin-expressing interneurons during epileptiform activation. The GCaMP6f signal strongly correlates with the level of activity for both interneuronal classes. Both interneuronal classes participate in the feedfoward inhibition. This contrasts starkly with the pattern of pyramidal recruitment, which is greatly delayed. During these barrages, both sets of interneurons show intense bursting, at rates up to 300Hz, which is time-locked to the glutamatergic barrages. The activity of parvalbumin-expressing interneurons appears to peak early in the pre-ictal period, and can display depolarizing block during the ictal event. In contrast, somatostatin-expressing interneuronal activity peaks significantly later, and firing persists throughout the ictal events. Interictal events appear to be very similar to the pre-ictal period, albeit with slightly lower firing rates. Thus, the inhibitory restraint arises from a coordinated pattern of activity in the two main classes of cortical interneurons.

Pyramidal cell activity levels affect the polarity of activity-induced gene transcription changes in interneurons.

Changes in gene expression are an important mechanism by which activity levels are regulated in the nervous system. It is not known, however, how network activity influences gene expression in interneurons; since they themselves provide negative feedback in the form of synaptic inhibition, there exists a potential conflict between their cellular homeostatic tendencies and those of the network. We present a means of examining this issue, utilizing simple in vitro models showing different patterns of intense network activity. We found that the degree of concurrent pyramidal activation changed the polarity of the induced gene transcription. When pyramidal cells were quiescent, interneuronal activation led to an upregulation of glutamate decarboxylase 1 ( GAD1) and parvalbumin ( Pvalb) gene transcriptions, mediated by activation of the Ras/extracellular signal-related kinase mitogen-activated protein kinase (Ras/ERK MAPK) pathway. In contrast, coactivation of pyramidal cells led to an ionotropic glutamate receptor N-methyl-d-aspartate 2B-dependent decrease in transcription. Our results demonstrate a hitherto unrecognized complexity in how activity-dependent gene expression changes are manifest in cortical networks. NEW & NOTEWORTHY We demonstrate a novel feedback mechanism in cortical networks, by which glutamatergic drive, mediated through the Ras/ERK MAPK pathway, regulates gene transcription in interneurons. Using a unique feature of certain in vitro epilepsy models, we show that without this glutamatergic feedback, intense activation of interneurons causes parvalbumin and glutamate decarboxylase 1 mRNA expression to increase. If, on the other hand, pyramidal cells are coactivated with interneurons, this leads to a downregulation of these genes.

The contribution of raised intraneuronal chloride to epileptic network activity.

Altered inhibitory function is an important facet of epileptic pathology. A key concept is that GABAergic activity can become excitatory if intraneuronal chloride rises. However, it has proved difficult to separate the role of raised chloride from other contributory factors in complex network phenomena, such as epileptic pathology. Therefore, we asked what patterns of activity are associated with chloride dysregulation by making novel use of Halorhodopsin to load clusters of mouse pyramidal cells artificially with Cl(-). Brief (1-10 s) activation of Halorhodopsin caused substantial positive shifts in the GABAergic reversal potential that were proportional to the charge transfer during the illumination and in adult neocortical pyramidal neurons decayed with a time constant of τ = 8.0 ± 2.8s. At the network level, these positive shifts in EGABA produced a transient rise in network excitability, with many distinctive features of epileptic foci, including high-frequency oscillations with evidence of out-of-phase firing (Ibarz et al., 2010). We show how such firing patterns can arise from quite small shifts in the mean intracellular Cl(-) level, within heterogeneous neuronal populations. Notably, however, chloride loading by itself did not trigger full ictal events, even with additional electrical stimulation to the underlying white matter. In contrast, when performed in combination with low, subepileptic levels of 4-aminopyridine, Halorhodopsin activation rapidly induced full ictal activity. These results suggest that chloride loading has at most an adjunctive role in ictogenesis. Our simulations also show how chloride loading can affect the jitter of action potential timing associated with imminent recruitment to an ictal event (Netoff and Schiff, 2002).

RNA Sequencing Demonstrates Ex Vivo Neocortical Transcriptomic Changes Induced by Epileptiform Activity in Male and Female Mice.

Seizures are generally associated with epilepsy but may also be a symptom of many other neurological conditions. A hallmark of a seizure is the intensity of the local neuronal activation, which can drive large-scale gene transcription changes. Such changes in the transcriptional profile likely alter neuronal function, thereby contributing to the pathological process. Therefore, there is a strong clinical imperative to characterize how gene expression is changed by seizure activity. To this end, we developed a simplified ex vivo technique for studying seizure-induced transcriptional changes. We compared the RNA sequencing profile in mouse neocortical tissue with up to 3 h of epileptiform activity induced by 4-aminopyridine (4AP) relative to control brain slices not exposed to the drug. We identified over 100 genes with significantly altered expression after 4AP treatment, including multiple genes involved in MAPK, TNF, and neuroinflammatory signaling pathways, all of which have been linked to epilepsy previously. Notably, the patterns in male and female brain slices were almost identical. Various immediate early genes were among those showing the largest upregulation. The set of down-regulated genes included ones that might be expected either to increase or to decrease neuronal excitability. In summary, we found the seizure-induced transcriptional profile complex, but the changes aligned well with an analysis of published epilepsy-associated genes. We discuss how simple models may provide new angles for investigating seizure-induced transcriptional changes.

Synergistic Positive Feedback Mechanisms Underlying Seizure Initiation.

Investigations into seizure initiation, in recent years, have focused almost entirely upon alterations of interneuronal function, chloride homeostasis, and extracellular potassium levels. In contrast, little attention has been directed toward a possible role of dendritic plateau potentials in the actual ictogenic transition, despite a substantial literature dating back 40 years regarding its importance generally in epilepsy. Here, we argue that an increase in dendritic excitability, coordinated across the population of pyramidal cells, is a key stage in ictogenesis.

Xenon LFP Analysis Platform Is a Novel Graphical User Interface for Analysis of Local Field Potential From Large-Scale MEA Recordings.

High-density multi-electrode array (HD-MEA) has enabled neuronal measurements at high spatial resolution to record local field potentials (LFP), extracellular action potentials, and network-wide extracellular recording on an extended spatial scale. While we have advanced recording systems with over 4,000 electrodes capable of recording data at over 20 kHz, it still presents computational challenges to handle, process, extract, and view information from these large recordings. We have created a computational method, and an open-source toolkit built in Python, rendered on a web browser using Plotly's Dash for extracting and viewing the data and creating interactive visualization. In addition to extracting and viewing entire or small chunks of data sampled at lower or higher frequencies, respectively, it provides a framework to collect user inputs, analyze channel groups, generate raster plots, view quick summary measures for LFP activity, detect and isolate noise channels, and generate plots and visualization in both time and frequency domain. Incorporated into our Graphical User Interface (GUI), we also created a novel seizure detection method, which can be used to detect the onset of seizures in all or a selected group of channels and provide the following measures of seizures: distance, duration, and propagation across the region of interest. We demonstrate the utility of this toolkit, using datasets collected from an HD-MEA device comprising of 4,096 recording electrodes. For the current analysis, we demonstrate the toolkit and methods with a low sampling frequency dataset (300 Hz) and a group of approximately 400 channels. Using this toolkit, we present novel data demonstrating increased seizure propagation speed from brain slices of Scn1aHet mice compared to littermate controls. While there have been advances in HD-MEA recording systems with high spatial and temporal resolution, limited tools are available for researchers to view and process these big datasets. We now provide a user-friendly toolkit to analyze LFP activity obtained from large-scale MEA recordings with translatable applications to EEG recordings and demonstrate the utility of this new graphic user interface with novel biological findings.

Divergent paths to seizure-like events.

Much debate exists about how the brain transitions into an epileptic seizure. One source of confusion is that there are likely to be critical differences between experimental seizure models. To address this, we have compared the evolving activity patterns in two widely used in vitro models of epileptic discharges. Brain slices from young adult mice were prepared in the same way and bathed either in 0 Mg2+ or 100 µmol/L 4AP artificial cerebrospinal fluid. We have found that while local field potential recordings of epileptiform discharges in the two models appear broadly similar, patch-clamp analysis reveals an important difference in the relative degree of glutamatergic involvement. 4AP affects parvalbumin-expressing interneurons more than other cortical populations, destabilizing their resting state and inducing spontaneous bursting behavior. Consequently, the most prominent pattern of transient discharge ("interictal event") in this model is almost purely GABAergic, although the transition to seizure-like events (SLEs) involves pyramidal recruitment. In contrast, interictal discharges in 0 Mg2+ are only maintained by a very large glutamatergic component that also involves transient discharges of the interneurons. Seizure-like events in 0 Mg2+ have significantly higher power in the high gamma frequency band (60-120Hz) than these events do in 4AP, and are greatly delayed in onset by diazepam, unlike 4AP events. We, therefore, conclude that the 0 Mg2+ and 4AP models display fundamentally different levels of glutamatergic drive, demonstrating how ostensibly similar pathological discharges can arise from different sources. We contend that similar interpretative issues will also be relevant to clinical practice.

Graphical user interface for simultaneous profiling of activity patterns in multiple neuronal subclasses.

We provide notes on how to use a graphical user interface (GUI), implemented with MATLAB, for aligning imaging datasets of biological tissue. The original use was for matching two imaging data sets, where one set was taken of the living preparation and another was taken post-fixation and following immunohistochemical processing. This technique is described in detail in an accompanying paper (Parrish et al., [1], where we also include information about the experimental procedures, and examples of the usage of the GUI.

Simultaneous profiling of activity patterns in multiple neuronal subclasses.

BACKGROUND: Neuronal networks typically comprise heterogeneous populations of neurons. A core objective when seeking to understand such networks, therefore, is to identify what roles these different neuronal classes play. Acquiring single cell electrophysiology data for multiple cell classes can prove to be a large and daunting task. Alternatively, Ca2+ network imaging provides activity profiles of large numbers of neurons simultaneously, but without distinguishing between cell classes. NEW METHOD: We therefore developed a strategy for combining cellular electrophysiology, Ca2+ network imaging, and immunohistochemistry to provide activity profiles for multiple cell classes at once. This involves cross-referencing easily identifiable landmarks between imaging of the live and fixed tissue, and then using custom MATLAB functions to realign the two imaging data sets, to correct for distortions of the tissue introduced by the fixation or immunohistochemical processing. RESULTS: We illustrate the methodology for analyses of activity profiles during epileptiform events recorded in mouse brain slices. We further demonstrate the activity profile of a population of parvalbumin-positive interneurons prior, during, and following a seizure-like event. COMPARISON WITH EXISTING METHODS: Current approaches to Ca2+ network imaging analyses are severely limited in their ability to subclassify neurons, and often rely on transgenic approaches to identify cell classes. In contrast, our methodology is a generic, affordable, and flexible technique to characterize neuronal behaviour with respect to classification based on morphological and neurochemical identity. CONCLUSIONS: We present a new approach for analysing Ca2+ network imaging datasets, and use this to explore the parvalbumin-positive interneuron activity during epileptiform events.