- 1Surgical Neurology Branch, NINDS, National Institutes of Health, Baltimore, MD, United States
- 2Department of Neurosurgery, University of Maryland School of Medicine, Baltimore, MD, United States
- 3Department of Neurosurgery, Brigham and Women’s Hospital, Boston, MA, United States
- 4Department of Neurology, Yale University, New Haven, CT, United States
- 5Department of Neurosurgery, Yale University, New Haven, CT, United States
Epilepsy is a common neurological disorder, affecting over 65 million people worldwide. Unfortunately, despite resective surgery, over 30
Introduction
Epilepsy affects over 65 million people worldwide (Ngugi et al., 2010). When patients fail medical management with anti-seizure medications (ASMs), they must turn to surgical evaluation. However, despite resective surgery and/or neuromodulatory intervention, approximately 50
Epilepsy models suggest that seizures are characterized by excessive synchronous neuronal firing. Primary epilepsies (i.e., genetic epilepsies) are modeled by knock-out mice (ex. Kcnq3, Lgi1, Mecp2, etc) and secondary epilepsies are modeled using induction protocols (Marshall et al., 2021). These typically involve kindling via electrical stimulation or chemoconvulsant agents (ex. Pilocarpine, Pentylenetetrazole (PTZ), Kainate) that are introduced on repeated occasions, resulting in eventual spontaneous seizures (Curia et al., 2008; Lévesque and Avoli, 2013; Dhir, 2012; Goddard, 1967). Electrophysiology and seizure activity are typically recorded using invasive cortical/subcortical electrodes (Marshall et al., 2021).
In contrast to animal models, in vitro epilepsy models typically consist of mammalian brain slices, derived from whole brain or hippocampus, for acute or chronic (organotypic slices) study (Raimondo et al., 2017). These slices are obtained either from animal models of epilepsy or naive-animals. Electrophysiology from slices is recorded using a number of options, including grease gas chambers, ion-selective microelectrodes, functional microscopy and optogenetics, and most commonly, patch-clamp (Raimondo et al., 2017; Neher and Sakmann, 1976; Sakmann and Neher, 1984). Patch-clamp is an electrophysiological technique developed in the 1970s that allows for study of single-neuron electrical behavior using a micropipette in tight contact with the cell membrane (Segev et al., 2016). Patch-clamp recordings afford many advantages, including the ability to record from neurons in the context of their native, preserved circuitry and also capturing changes in activity in real-time. However, by virtue of how patch-clamp recordings are obtained, activity across multi-unit neuronal populations over long periods of time are not feasible. Furthermore, the quality of recordings from patch-clamp can vary greatly depending on cell-specific conditions, presenting a further challenge (Saleem et al., 2023). As such, in vitro slice models and patch-clamp recordings do not allow for effective study of network activity or changes.
Over the past few years, MEA technology has been increasingly used for in vitro investigations to better understand neuronal dynamics, balancing neuron-glia populations and a well-controlled environment to model neurologic disease (Mossink et al., 2021; Cerina et al., 2023). MEAs encompass microelectrodes at the base of tissue culture wells, on which neuronal populations are plated. These microelectrodes allow for capture of high-frequency neuronal spiking activity by recording of extracellular field potential from multiple electrodes, which studies have shown strongly resembles intracellular waveforms (Ashida et al., 2012; Funabiki et al., 2011; Kuokkanen et al., 2018). Additionally, MEA microelectrodes allow for delivery of electrical stimulation, allowing for intentional probing of in vitro neuronal networks. As discussed above, traditional in vitro epilepsy models use patch-clamp recordings from brain slices. Although patch-clamp directly captures intracellular neuronal action potentials, it fails to capture activity from populations of neurons. MEAs offer a unique advantage as recordings from extracellular field potential allow for useful estimates of population activity. In recent years, MEAs have been used to model neurological disease, such as Alzheimer’s disease (AD), Glioblastoma Multiforme (GBM), Parkinson’s disease (PD), Amyotrophic Lateral Sclerosis (ALS), and epilepsy (Li et al., 2020; Krishna et al., 2023; Wainger et al., 2014; Woodard et al., 2014; Tidball and Parent, 2016). Considering increasing evidence suggesting pathomechanisms surrounding epileptogenesis reflect distributed seizure networks, MEAs offer a robust tool in which such changes can be studied and modeled.
In a cohort of patients with drug-resistant epilepsy (DRE), we found evidence of increased functional connectivity (FC) between regions of seizure onset and spread. Based on these findings, it is plausible that years of repeated seizure insults may result in changes in underlying neuronal connectivity, that result in interictal FC changes. As such, to better understand network changes that may take place in epilepsy patients, we describe an in vitro epilepsy model, motivated by our findings in DRE patients, where we capture seizure-like hyperexcitability over time and monitor network changes longitudinally. We use 4-aminopyridine (4-AP), a known pro-convulsant agent, to capture hyperexcitable activity and model seizures. We use electrophysiology recorded from MEAs to characterize network changes after 2 days of chronic 4-AP treatment. We posited that we would be able to effectively capture network changes in our seizure model, serving as a proof-of-concept that our in vitro MEA seizure model could be used to study seizure network development over time. We find that 4-AP successfully captures hyperexcitable activity that persists at baseline. Furthermore, we find after 2 days of chronic treatment, network connectivity is significantly increased in 4-AP treated wells, compared to control, when considering functional connectivity (FC) and a novel technique using dimensionality reduction. These results support the idea that 4-AP can be used to effectively model hyperexcitable activity on MEAs and allow for in vitro study of network changes, that may provide insight into the pathogenesis of seizure networks.
Materials and methods
Functional connectivity in human epilepsy patients
We retrospectively identified three patients (1 female, 39.0
For each interictal epoch, we considered contacts recording from regions of seizure onset (SOZ), spread (SP), and uninvolved controls (Figure 1A). We were interested in comparing FC between SOZ-SP and SOZ-control. For each contact, we preprocessed signals by applying a local detrending procedure to remove slow fluctuations from the time series. We then used a regression-based approach to remove line noise at 60 and 120 Hz (Kutsy, 1999). We used a low-pass type 1 finite impulse response (FIR) filter (order = 180), to remove higher order line harmonics as well as high-frequency noise. After preprocessing signals, we applied a band-pass FIR filter (order = 390) designed using the Parks-McClellan algorithm in the high-gamma band (70–90 Hz) (Parks and McClellan, 1972). We were interested in the high-gamma band considering its relevance in seizure onset regions (Ren et al., 2015). We computed phase-locking coherence (PLC) based on raw real Hilbert transformed time series, comparing estimates across SOZ-SP and SOZ-control pairs (Lachaux et al., 1999). We considered PLC as our FC metric as it measures the degree of synchronization between icEEG signals, with high values reflecting greater synchronization, and this connectivity, and low values reflecting weaker (Figure 1B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1441345/fnetp-04-1441345-HTML/image_m/fnetp-04-1441345-g001.jpg)
Figure 1. Functional connectivity in human epilepsy patients. (A) Here, we show the schematic for how FC was computed in our cohort of DRE patients. Specifically, we considerd icEEG obtained during neuromonitoring for seizure localization. We selected for contacts recording from regions identified as SOZ, SP, or control, equidistant to the SOZ. We then computed PLC and compared across all SOZ-SP and SOZ-control contact pairs across all patients (left). (B) Here, we depict extreme cases of what PLC measures. Specifically, PLC captures how in-phase, or synchronized two signals are. High PLC values (top) mean signals are synchronized, and therefore suggests that the underlying recording sites may be functionally connected. On the contrary, low PLC values recording sites are likely unrelated to one another. (C) We average PLC computed for all SOZ-SP and SOZ-control contact pairs. We show that within a given contact pair, SOZ-SP contacts are more functionally connected than their respective SOZ-control pair (left). When we compare PLC across all patients, we find that SOZ-SP contact pairs are truly significantly more synchronized, and hence likely more strongly connected, compared to SOZ-control pairs (0.29 versus 0.25, t(2) = 6.72, p = 0.0215, paired t-test, right). Abbreviations: FC, functional connectivity; DRE, drug-resistant epilepsy; SOZ, seizure-onset zone; SP, seizure spread; icEEG, intracranial electrocorticography; PLC, phase-locking coherence.
Neuron-glia culture preparation
We established a mixed neuron-glia rat cortical cell culture on MEAs according to previously published protocols (Ksendzovsky et al., 2022; Mortazavi et al., 2022). Cortices were dissected from newborn P1 rat pups in a modified Puck’s dissociation medium D1 (5 mM HEPES, 16.5 mM glucose, 22 mM Sucrose, 137 mM NaCl, 0.32 mM
Induction of hyperexcitable activity using 4-AP
There are numerous techniques to capture hyperexcitable activity in dissociated cultures on MEAs (Grainger et al., 2018). One approach is introduction of proconvulsant agents (ex. Picrotoxin, Gabazine, PTZ, Bicuculline, Tutin, Tranexamic Acid, Endosulfan, 4-Aminopyridine (4-AP), SNC80, NMDA, Linopirdine, Strychnine HCL, Amoxapine, Pilocarpine HCL, Thioridazine HCL, Domoic Acid, Tetrodotoxin (TTX)) (Bradley et al., 2018). Another is by manipulation of ion concentrations (ex. low
Experiments with 4-AP began after neuron-glia culture maturation and differentiation, around day-in-vitro (DIV) 20, as suggested by previous literature (Latchoumane et al., 2018) and direct observation of stabilized mean firing rate (MFR) across at least 3 days. Experimental protocol consisted of obtaining 5-min baseline MEA recordings prior to 4-AP treatment. After baseline MEA recordings, concentrated 4-AP solution was applied to treatment wells, achieving a 500
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1441345/fnetp-04-1441345-HTML/image_m/fnetp-04-1441345-g002.jpg)
Figure 2. Experimental workflow. Here, we show the experimental workflow, delineating 4-AP treatment protocol. Treatment begins Day 0, after MFR has stabilized across neuron-glia cultures (see Methods). MEA recordings are collected prior to each day’s 4-AP wash-in using the Maestro Pro MEA system (Axion Biosystems, Atlanta, GA). 6-well MEA plates are used, with three wells serving as untreated controls and remaining used for 4-AP treatment. Protocol was repeated across multiple biological replicates, with each 6-well MEA plate reflecting a single biological replicate. During a single day’s 4-AP treatment, treatment wells (red) were spiked with concentrated stock solution of 4-AP. Simultaneously, control wells (grey) were treated with equivalent volume of Neurobasal medium. Treatment period lasts for 30-min, after which all wells (control and 4-AP-treated) undergo full wash-out and media change with Neurobasal medium. Treatment is repeated for 2 days, at approximately the same time each day. As such, each MEA recording is effectively collected 24-h apart, allowing for quantifying chronic changes. Abbreviations: 4-AP, 4-aminopyridine; MFR, mean firing rate; MEA, microelectrode arrays.
MEA recording acquisition and preprocessing
In vitro electrophysiology was collected using the Maestro Pro MEA system (Axion Biosystems, Atlanta, GA) and AxIS Navigator software v3.5.2 (Axion BioSystems, Atlanta, GA). We used six-well MEA plates, due to their high microelectrode density (64 microelectrodes). Microelectrodes are made of polymer poly (3,4-ethylenedioxythiophene) (PEDOT) and arranged in an 8 × 8 grid, with a 50
To extract spiking activity, we processed all electrophysiology recordings offline using custom MATLAB (MathWorks, Natick, MA) scripts. Preprocessing steps involved downsampling signals from 12.5 kHz to 2.5 kHz and band-pass filtering signals between 0.1 and 300 Hz (Butterworth, order = 20). These steps were done to eliminate the impacts of low-frequency drift, DC artifacts, and filter for multi-unit activity (MUA) to capture high-frequency neuronal spiking (Stark and Abeles, 2007; Waldert et al., 2013; Ahmadi et al., 2021). Subsequently, signals were filtered using a series of notch filters (Butterworth, order = 4) to eliminate 60 Hz line noise and its first two harmonics (i.e., 120, 180 Hz), with a narrow bandwidth within 1 Hz of the target frequency. All signals were subsequently z-scored.
Functional connectivity
We were interested in computing changes that estimate synchronized neural activity, or FC, across all multi-unit populations. To this end, we considered degree of correlation, a measure of statistical interdependence, between spike trains (Cohen et al., 2002; Selinger et al., 2004; Chiappalone et al., 2006; Chiappalone et al., 2011; Eisenman et al., 2015). For each spike train recorded from respective microelectrodes, we vectorized spike counts by binning spikes in non-overlapping 50-ms windows. FC was then quantified by computing Pearson’s correlations (Fisher z-transformed) across all unique pairwise binned spike trains in each well.
Of note, considering our choice of Pearson’s correlation (Fisher z-transformed) across pairs of binned spike trains, FC estimates would be biased when including periods of highly-synchronous firing, such as that observed during population bursts. This is a well understood confounder, in which correlation as a measure of spike train synchrony is positively biased in periods of higher firing rate (Cutts and Eglen, 2014). To assess whether inclusion of population burst periods in FC analysis would truly bias estimates in our data, we compared impact of including and excluding population burst periods on FC estimates. We observed in representative MEA recordings that FC was significantly greater when including burst periods, as anticipated (
Mapping network connectivity in low-dimension space
The aim of the presented method is to model seizure activity and be able to characterize network changes. To this end, we compute FC, as described above, to estimate connectivity between pairwise multiunit populations. However, network dynamics can be complex and therefore may not be appropriately captured by simply comparing difference of means of computed FC estimates. For example, if a treatment group exhibits a small number of pairwise connections that are preferentially strengthened, this difference may not be captured by comparing difference of means between treatment and control groups as connectivity distributions may be skewed. As such, it is important we devise a method that is able to capture relevant network changes.
In recent years, dimensionality reduction has emerged as a powerful tool for revealing patterns in complex neural data (Langdon et al., 2023; Cunningham and Yu, 2014). Here, we use the high-dimensional connectivity information we compute and represent this in a low-dimension embedding, using principal components analysis (PCA), allowing us to characterize how network connectivity evolves across a single well, over time.
We performed dimensionality reduction separately for each well considering each individual well has unique network connections that may evolve differently over time. Conducting dimensionality reduction across different wells represented in the same feature matrix would confuse the dimensionality reduction procedure into treating each pairwise connectivity estimate as a feature representing the same relative information across different wells, when in reality, pairwise connectivity estimates only serve as meaningful features in the context of individual wells.
We constructed a single m
where (i, j)
Statistical analyses
All data analysis was conducted using custom MATLAB scripts (MathWorks, Natick, MA) and Prism (GraphPad, San Diego, CA). Groups were compared using two-tailed t-tests. Parametric statistical tests were used considering data were normally distributed (p
Results
Evidence of seizure networks in epilepsy patients
To investigate seizure networks in epilepsy, we computed FC between regions involved in seizure onset and primary spread, and compared this to FC between seizure onset and equidistant, uninvolved control regions. We hypothesized that areas connected by seizure spread are more strongly connected compared to those that are not, lending evidence to the network theory of epilepsy. In our patient cohort, we find that regions involved in seizure spread are in fact more strongly synchronized in the high-gamma band (70–90 Hz), compared to control regions (0.29 versus 0.25, t (2) = 6.72, p = 0.0215, paired t-test, Figure 1C). From this, it is clear that on average, SOZ and SP regions are more synchronized than SOZ and control regions, suggesting that regions of seizure activity are more strongly functionally connected, and may comprise a pathologic seizure network that strengthens over time, serving as the primary motivation for our in vitro seizure model.
4-AP captures hyperexcitability in vitro
The presented model is predicated on effectively capturing hyperexcitable conditions in vitro to model high firing conditions observed during seizures. To assess whether 4-AP-induced hyperexcitability was reflective of this, we considered spike rasters and changes in neuronal firing activity. When considering representative spike rasters from control (Figure 3A) and 4-AP-treated (Figure 3B) and wells, apparent differences in firing activity are observed, with the 4-AP-treated well exhibiting increased spiking density and periods of synchronous firing (bursts). Changes in firing activity were quantified using mean firing rate (MFR) normalized to pre-treatment baseline (i.e., Day 0, see Figure 2). We observed that after just one treatment with 4-AP, the next day’s MFR was significantly elevated (1.90
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1441345/fnetp-04-1441345-HTML/image_m/fnetp-04-1441345-g003.jpg)
Figure 3. 4-aminopyridine successfully captures hyperexcitable activity in vitro. To assess if 4-AP could capture hyperexcitability in our neuron-glia cultures on MEAs, we considered spike rasters and MFR before and after chronic 4-AP treatment. (A) Here, we show spike rasters from a representative control well before 4-AP treatment protocol (left) and after (right). We observe no discernable changes in spiking density. (B) Here, we consider spike rasters before (left) and after (right) 4-AP treatment. Compared to control, we note significantly higher spiking density and more synchronous firing periods, consistent with expectations. (C) We show changes in MFR throughout 4-AP treatment protocol, normalizing to pre-treatment MFRs, across multiple technical (9) and biological (3) replicates. MFR is significantly moreso increased in 4-AP-treated wells compared to control after just 1 day of treatment (1.90
Increased functional connectivity after chronic 4-AP treatment
To show the presented model effectively captures network changes, we compared functional connectivity (FC) (Pearson’s correlation, Fisher z-transformed
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1441345/fnetp-04-1441345-HTML/image_m/fnetp-04-1441345-g004.jpg)
Figure 4. Increased baseline functional connectivity after chronic 4-aminopyridine treatment. To show our model can be used to study network changes, we considered FC changes before and after chronic 4-AP treatment. We quantified FC considering pairwise Pearson’s correlation (Fisher z-transformed
Network connectivity occupies distinct subspace after 4-AP treatment
We reduced high-dimensional functional connectivity information to a low-dimension embedding to allow for study of how a individual well’s network connectivity evolves over time (see Methods). Using this method, we are able to plot how a well’s network connectivity is evolving over time and compare network connectivity changes across different wells. We do so across all control and 4-AP-treated wells. Specifically, we observe that network connectivity patterns in control wells tend to cluster around pre-treatment observations, suggesting that networks are relatively stable, and unchanged over 2 days of recordings (Figure 5A). However, when we consider 4-AP-treated wells, we observe that after one and 2-days of 4-AP treatment, network connectivity patterns deviate from what is observed pre-treatment, suggesting that 4-AP causes network changes over time (Figure 5A). We quantify these network changes by computing distance of points to Day 0 recordings, on the low-dimension space. We find that after 2 days of chronic 4-AP treatment, network connectivity patterns occupy a unique subspace distinct from pre-treatment, compared to control (22.43
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1441345/fnetp-04-1441345-HTML/image_m/fnetp-04-1441345-g005.jpg)
Figure 5. Temporal evolution of network connectivity mapped in low-dimension space. As network connectivity may be complex and change dynamically, we developed a method to track how a well’s network connectivity evolves over time by reducing high-dimensional connectivity information to a low-dimension embedding (see Methods). (A) Here, we show how network connectivity evolves across 3 days in representative control (blue) and 4-AP-treated (red) wells. Timepoint represented is delineated next to each point. We observe that after 2 days of chronic 4-AP treatment, baseline network connectivity occupies a distinct subspace, far from pre-treatment network connectivity. In contrast, we observe that network connectivity in a representative control well occupies a smaller subspace, closer to what is observed Day 0. (B) To quantify change in network connectivity over time, we compute the Euclidean distance of each point to Day 0. We find that after 2 days, compared to controls, 4-AP-treated wells show significant changes in network connectivity compared to pre-treatment (22.43
Discussion
A significant challenge in studying neurologic disease is the availability of representative models. Although it is most optimal to study disease pathology at the patient level, oftentimes it is not possible to truly understand the pathomechanisms of a disease without a fundamental understanding of underlying mechanistic changes. In vivo and in vitro models offer the ability to study principal changes in the context of representative models of disease, greatly contributing to our understanding of them.
Epilepsy is a common neurological disorder, however, very little is understood of the mechanisms surrounding epileptogenesis. Until recently, epilepsy has largely been seen as a focal disorder, with seizure models reflecting this. However, over recent years, we have come to understand epilepsy as a network disorder, recruiting various cortical and subcortical regions into a pathologic seizure network (Kramer and Cash, 2012; Scharfman et al., 2018; Piper et al., 2022; Bröhl et al., 2023). Numerous intracranial studies considering network connectivity patterns have lent support to this idea, some even suggesting that disruption of seizure networks is likely critical in achieving seizure-freedom (Andrews et al., 2019; Neal et al., 2020; Jehi, 2017; Rijal et al., 2023; Taylor et al., 2018; Josephson et al., 2013).
In fact, here, we show evidence of high interictal FC between areas linked by seizure spread. These findings lend direct support for the network theory of epilepsy, suggesting that seizure activity travels along more strongly connected neuronal populations. This corroborates previous findings in support of the network theory of epilepsy and further suggests that longitudinal network changes are likely involved in seizure network formation, resulting in baseline changes that can be observed interictally. However, icEEG studies considering seizure networks only provide snapshots in time, and preclude the ability to study longitudinal network development. As such, it is imperative we have models that allow us to investigate network changes over time in the context of epileptogenesis. However, to date, there are no seizure models that consider network changes, with most simply focusing on studying molecular changes underlying hyperexcitability or focal alterations. As such, in this study, we created a novel in vitro seizure model in order to model changes in neuronal network connectivity associated with epileptogenesis, motivated by our human epilepsy findings.
As previously described, here we show 4-AP chronically increases neuronal firing, akin to an in vitro correlate of kindling, where residual hyperexcitability may lead to compensatory and pathologic changes within epileptic networks (Heuzeroth et al., 2019; Leite et al., 2005; Lillis et al., 2015). In addition, we demonstrate that 4-AP induces FC changes, likely secondary to induced hyperexcitability, which may replicate the changes that occur in human epileptic networks, providing a bridge between experimental models and clinical epilepsy. We further show changes in network activity using a dimensionality reduction-based method to characterize network connectivity changes in MEA wells.
Because network dynamics are complex, FC may fail to appropriately capture changes that could be informative of strengthening or weakening of specific pairwise connections. Additionally, it is challenging to capture complex network changes using a single measurement. As such, we reduced high-dimensional pairwise FC information into a low-dimension embedding, from which network connectivity information for an individual well can be simply visualized and changes across time can be quantified. Using this, we show changes in wells’ network connectivity over time and find that these results parallel changes observed when considering FC alone, suggesting that this method retains information conveyed by traditional FC measures. Ultimately, we believe using dimensionality reduction to reduce network connectivity information allows for more accurate capture of how a MEA well’s connectivity patterns change over time.
While exisiting in vitro models capture aspects of epileptic activity, to the best of our knowledge, none allow for the detailed study of chronic network changes of seizure-like states as effectively as our proposed method. Our protocol holds many advantages, including the unique ability to study network changes and circuit reorganization afforded by the composition of dissociated neuron-glia cultures. The inclusion of both neurons and their supporting glial cells not only better represents the in vivo environment, but also provides a more adaptable experimental setup with enhanced survival and cellular development (Kaech and Banker, 2007). This flexibility is crucial for investigating network dynamics in the absence of preexisting connections, allowing for precise manipulations and observations of how epileptic networks evolve and respond to interventions. For example, recent studies have suggested that specific regulatory pathways are altered in epileptic tissue (Ksendzovsky et al., 2022; Hammer et al., 2019; Pfisterer et al., 2020; Sumadewi et al., 2023; Meng et al., 2015). Using our model, we can define the molecular basis of network changes by introducing small molecular inhibitors that target these pathways and investigate subsequent network changes.
Nevertheless, this model’s reliance on initially naïve networks may also be seen as a limitation, given that seizures in humans affect preexisting neural networks. To address this, we begin our 4-AP treatments post-culture maturation, aiming to more accurately represent changes observed in human cortical/subcortical networks. An additional limitation of our in vitro model is the inherent inability to accurately replicate epileptiform activity, as is possible in more physiologically accurate animal models of epilepsy (i.e., electrical kindling, Kainic Acid). Our model uses 4-AP, a potassium-channel blocker that impairs neuronal repolarization and hyperpolarization, this increases overall neuronal excitability in the network allowing for modeling of seizure-like conditions. However, it is important to note that potassium-channel blockade does not underly representative seizure activity seen in chronic models of epilepsy i.e., Kainic Acid animal models of epilepsy, which demonstrate spontaneous recurrent seizures. Thus, it is important to restate that our model captures hyperexcitable neuronal firing activity which is observed as a result of seizure-like activity, and not direct seizure activity (Wong, 2011). As such, in vitro models preclude the ability to characterize true ictal activity, and rather limit us to correlates such as neuronal firing activity (i.e., firing rate, bursting rate). These are well-acknowledged limitations of in vitro epilepsy models, as such, conclusions must be considered in this context (Oblasov et al., 2024). However, despite this, our model still allows for study of what may happen to neuronal networks under hyperexcitable conditions, such as those captured by seizure activity. A key hallmark of epilepsy is repeated seizure insults which induce hyperexcitable conditions in cortical/subcortical regions. Our 4-AP model was created to simulate chronic changes that occur in the context of repeated seizure events, reflective of what happens in epilepsy in vivo. Daily 4-AP treatments not only replicate kindling activity, but also provide a mechanism by which to simulate the temporal progression of changes associated with epileptogenic stimuli. While many theories of epileptogenesis exist, hyperexcitable onslaughts leading to maladaptive network changes in a kindled fashion is certainly among them (Bromfield et al., 2006; Marques et al., 2022). As such, it is certainly plausible to use our 4-AP model to understand how changes in network activity advance throughout epileptogenesis. In the future, this in vitro model and its network findings can be corroborated by in vivo animal studies using high-density MEAs implanted directly on the cortical surface. Recent groups have shown that Utah Electrode Arrays (UEAs) (Blackrock Microsystems, Salt Lake City, UT) can be implanted directly on the cortical surface of adult rats, allowing for both multi-unit and single-unit recordings (Black et al., 2018; Nolta et al., 2015). Using this, similar to our in vitro model, high-resolution neuronal network connectivity patterns and changes can be characterized in the context of a more representative animal epilepsy model.
Additionally, our method’s integration of MEA recordings enable investigation of multi-unit populations, offering insights that are unattainable with traditional patch-clamp techniques used in slice models, which do not capture population-level dynamics. The use of microelectrodes in our model allows for capture of extracellular field potential at the level of multi-unit populations, enhancing our ability to characterize network dynamics, synaptic transmission, and plasticity in the in vitro setting (Hofmann and Bading, 2006).
Studying mechanistic network-level changes is crucial to advancing our pathophysiologic understanding of epilepsy and guiding therapeutic interventions. Here, we show that there likely exist interictal functional seizure networks in a cohort of DRE patients. We then show how these seizure networks may be modeled using in vitro neuron-glia cultures plated on MEAs. Given recent studies suggesting that disruption of seizure networks is important in achieving seizure control, our findings and in vitro model open a range of possible future interventions. For example, considering the persistence of seizure networks interictally, icEEG may be used to characterize critical nodes in a seizure network to guide subsequent surgical intervention (ex. resection, neuromodulation). An additional therapeutic application from our model is the optimization of neuromodulation for targeting seizure networks. Neuromodulation is an effective treatment option that uses brief pulse stimuli (BPS) to reduce seizure activity prior to its full-blown onset. Although offering a new and improved therapeutic approach to seizure control, neuromodulation still fails to achieve seizure control in a subset of patients (Boddeti et al., 2022). Recent studies suggest that this treatment failure may be due to underlying network connectivity patterns, which may inform those who may benefit from neuromodulation and those who may not. To this end, characterizing a patient’s interictal seizure network connectivity patterns using icEEG may inform those who may benefit from neuromodulation versus those who may not. Additionally, recent studies suggest that neuromodulatory devices exerts their therapeutic effects by modulating network activity, as opposed to focal activity. As such, our presented in vitro MEA model offers a unique platform in which this hypothesis can be tested, as MEA technology allows for applying direct electrical stimulation. Furthermore, this allows for testing various neuromodulation paradigms and their effects in 1) modulating network activity and 2) controlling seizure-like activity. Ultimately, using MEA technology to better understand network activity in seizure models allows for optimization of preexisting treatments such as neuromodulation and also discovery of new ones, by allowing for study of specific pathway inhibitors and pharmacologic interventions on network connectivity changes.
Conclusion
Here, we present a novel in vitro model to study seizure networks using neuron-glia populations cultured on MEAs. We show that our 4-AP model serves as a robust in vitro tool for modeling connectivity changes associated with epileptogenesis. This model captures not only acute hyperexcitablity, a hallmark of epileptic networks, but also chronic network adaptations over time, as evidenced by increased FC of underlying neuronal networks. Furthermore, we introduce a novel method for characterizing aggregate network connectivity and how it changes over time. We believe that our model and methods effectively capture network dynamics akin to those in epilepsy patients. Considering the recent idealogical shift and strong evidence for epilepsy as a network disorder, ultimately, we believe that our model can be used to not only better understand epileptogenesis, but also begin to develop new therapeutics, and adapt preexisting ones (i.e., neuromodulation) to target pathologic seizure networks.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving humans were approved by the Yale Institutional Review Board. The studies were conducted in accordance with the local legislation and institutional requirements. The ethics committee/institutional review board waived the requirement. The ethics committee/institutional review board waived the requirement of written informed consent for participation from the participants or the participants’ legal guardians/next of kin because this study used retrospective data that is deidentified. The animal study was approved by University of Maryland School of Medicine Institutional Animal Care and Use Committee. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
UB: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Resources, Software, Validation, Visualization, Writing - original draft, Writing - review and editing. JL: Writing - original draft, Writing - review and editing. DM: Data curation, Writing - review and editing. MA: Writing - review and editing. MB: Resources, Supervision, Writing - review and editing. HZ: Data curation, Writing - review and editing. DS: Data curation, Writing - review and editing. KZ: Conceptualization, Methodology, Resources, Supervision, Writing - review and editing. AK: Conceptualization, Funding acquisition, Methodology, Resources, Supervision, Writing - review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. Research support was provided by the NIH Medical Research Scholars Program, a public-private partnership supported jointly by the NIH and contributions to the Foundation for the NIH from private donors.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Ahmadi, N., Constandinou, T. G., and Bouganis, C.-S. (2021). Inferring entire spiking activity from local field potentials. Sci. Rep. 11, 19045. doi:10.1038/s41598-021-98021-9
Andrews, J. P., Gummadavelli, A., Farooque, P., Bonito, J., Arencibia, C., Blumenfeld, H., et al. (2019). Association of seizure spread with surgical failure in epilepsy. JAMA Neurol. 76, 462–469. issn: 2168-6149. doi:10.1001/jamaneurol.2018.4316
Ashida, G., Funabiki, K., Kuokkanen, P. T., Kempter, R., and Carr, C. E. (2012). Signal-to-noise ratio in the membrane potential of the owl’s auditory coincidence detectors. J. Neurophysiology 108, 2837–2845. issn: 1522-1598. doi:10.1152/jn.00366.2012
Black, B. J., Kanneganti, A., Joshi-Imre, A., Rihani, R., Chakraborty, B., Abbott, J., et al. (2018). Chronic recording and electrochemical performance of Utah microelectrode arrays implanted in rat motor cortex. J. neurophysiology 120, 2083–2090. issn: 1522-1598. doi:10.1152/jn.00181.2018
Boddeti, U., McAfee, D., Khan, A., Bachani, M., and Ksendzovsky, A. (2022). Responsive neurostimulation for seizure control: current status and future directions. Biomedicines 10, 2677–9059. doi:10.3390/biomedicines10112677
Bradley, J. A., Luithardt, H. H., Metea, M. R., and Strock, C. J. (2018). In vitro screening for seizure liability using microelectrode array technology. Toxicol. Sci. 163, 240–253. issn: 1096-6080. doi:10.1093/toxsci/kfy029
Bröhl, T., Rings, T., Pukropski, J., von Wrede, R., and Lehnertz, K. (2023). The time-evolving epileptic brain network: concepts, definitions, accomplishments, perspectives. Front. Netw. Physiology 3, 1338864. issn: 26740109. doi:10.3389/fnetp.2023.1338864
Bromfield, E. B., Cavazos, J. E., and Sirven, J. I. (2006). Basic mechanisms underlying seizures and epilepsy.
Cerina, M., Piastra, M. C., and Frega, M. (2023). The potential of in vitro neuronal networks cultured on micro electrode arrays for biomedical research. Prog. Biomed. Eng. 5, 032002. issn: 2516-1091. doi:10.1088/2516-1091/acce12
Chiang, C. C., Wei, X., Ananthakrishnan, A. K., Shivacharan, R. S., Gonzalez-Reyes, L. E., Zhang, M., et al. (2018). Slow moving neural source in the epileptic hippocampus can mimic progression of human seizures. Sci. Rep. 8. issn: 2045-2322. doi:10.1038/s41598-018-19925-7
Chiappalone, M., Bove, M., Vato, A., Tedesco, M., and Martinoia, S. (2006). Dissociated cortical networks show spontaneously correlated activity patterns during in vitro development. Brain Res. 1093, 41–53. issn: 0006-8993. doi:10.1016/j.brainres.2006.03.049
Chiappalone, M., Vato, A., Berdondini, L., Koudelka-Hep, M., and Martinoia, S. (2011). Network dynamics and synchronous activity in cultured cortical neurons. Int. J. Neural Syst. 17, 87–103. issn: 01290657. doi:10.1142/S0129065707000968
Choquet, D., and Korn, H. (1992). Mechanism of 4-aminopyridine action on voltage-gated potassium channels in lymphocytes. J. General Physiology 99, 217–240. issn: 0022-1295. doi:10.1085/jgp.99.2.217
Cohen, J., Cohen, P., West, S., and Aiken, L. (2002). Applied multiple regression/correlation analysis for the behavioral sciences. Mahwah, NJ: Routledge. isbn: 9780203774441.
Cramer, C. L., Stagnitto, M. L., Knowles, M. A., and Palmer, G. C. (1994). Kainic acid and 4-aminopyridine seizure models in mice: evaluation of efficacy of anti-epileptic agents and calcium antagonists. Life Sci. 54, PL271–PL275. issn: 0024-3205. doi:10.1016/0024-3205(94)00845-0
Cunningham, J. P., and Yu, B. M. (2014). Dimensionality reduction for large-scale neural recordings. Nat. Neurosci. 17, 1500–1509. issn: 1546-1726. doi:10.1038/nn.3776
Curia, G., Longo, D., Biagini, G., Jones, R. S., and Avoli, M. (2008). The pilocarpine model of temporal lobe epilepsy. J. Neurosci. Methods 172, 143–157. issn: 0165-0270. doi:10.1016/j.jneumeth.2008.04.019
Curran-Everett, D. (2000). Multiple comparisons: philosophies and illustrations. Am. J. Physiology - Regul. Integr. Comp. Physiology 279, R1–R8. issn: 03636119 (1 48-1. doi:10.1152/ajpregu.2000.279.1.R1
Cutts, C. S., and Eglen, X. S. J. (2014). Detecting pairwise correlations in spike trains: an objective comparison of methods and application to the study of retinal waves. J. Neurosci. 34, 14288–14303. issn: 15292401. doi:10.1523/JNEUROSCI.2767-14.2014
Dhir, A. (2012). Pentylenetetrazol (PTZ) kindling model of epilepsy. Curr. Protoc. Neurosci. 9, Unit9.37. issn: 1934-8576 (SUPPL.58 2012). doi:10.1002/0471142301.ns0937s58
Eisenman, L. N., Emnett, C. M., Mohan, J., Zorumski, C. F., and Mennerick, S. (2015). Quantification of bursting and synchrony in cultured hippocampal neurons. J. Neurophysiology 114, 1059–1071. issn: 1522-1598. doi:10.1152/jn.00079.2015
Funabiki, K., Ashida, G., and Konishi, M. (2011). Computation of interaural time difference in the owl’s coincidence detector neurons. J. Neurosci. 31, 15245–15256. issn: 1529-2401. doi:10.1523/JNEUROSCI.2127-11.2011
Goddard, G. V. (1967). Development of epileptic seizures through brain stimulation at low intensity. Nature 214, 1020–1021. issn: 0028-0836. doi:10.1038/2141020a0
Gonzalez-Sulser, A., Wang, J., Motamedi, G. K., Avoli, M., Vicini, S., and Dzakpasu, R. (2011). The 4-aminopyridine in vitro epilepsy model analyzed with a perforated multi-electrode array. Neuropharmacology 60, 1142–1153. doi:10.1016/j.neuropharm.2010.10.007
Grainger, A. I., King, M. C., Nagel, D. A., Parri, H. R., Coleman, M. D., and Hill, E. J. (2018). In vitro models for seizure-liability testing using induced pluripotent stem cells. Front. Neurosci. 12, 590. issn: 1662453X. doi:10.3389/fnins.2018.00590
Hammer, M. F., Sprissler, R., Bina, R. W., Lau, B., Johnstone, L., Walter, C. M., et al. (2019). Altered expression of signaling pathways regulating neuronal excitability in hippocampal tissue of temporal lobe epilepsy patients with low and high seizure frequency. Epilepsy Res. 155, 106145–106844. issn. doi:10.1016/j.eplepsyres.2019.05.013
Heuzeroth, H., Wawra, M., Fidzinski, P., Dag, R., and Holtkamp, M. (2019). The 4-aminopyridine model of acute seizures in vitro elucidates efficacy of new antiepileptic drugs. Front. Neurosci. 13, 677. doi:10.3389/fnins.2019.00677
Hofmann, F., and Bading, H. (2006). Long term recordings with microelectrode arrays: studies of transcription-dependent neuronal plasticity and axonal regeneration. J. Physiology-Paris 99, 125–132. issn: 0928-4257. doi:10.1016/j.jphysparis.2005.12.005
Jehi, L. (2017). Outcomes of epilepsy surgery for epileptic networks. Epilepsy Curr. 17, 160–162. issn: 15357511. doi:10.5698/1535-7511.17.3.160
Josephson, C. B., Dykeman, J., Fiest, K. M., Liu, X., Sadler, R. M., Jette, N., et al. (2013). Systematic review and meta-analysis of standard vs selective temporal lobe epilepsy surgery. Neurology 80, 1669–1676. issn: 1526-632X. doi:10.1212/WNL.0b013e3182904f82
Kaech, S., and Banker, G. (2007). Culturing hippocampal neurons. Nat. Protoc. 1, 2406–2415. issn: 1750-2799. doi:10.1038/nprot.2006.356
Kramer, M. A., and Cash, S. S. (2012). Epilepsy as a disorder of cortical network organization. Neurosci. 18, 360–372. issn: 10738584. doi:10.1177/1073858411422754
Krishna, S., Choudhury, A., Keough, M. B., Seo, K., Ni, L., Kakaizada, S., et al. (2023). Glioblastoma remodelling of human neural circuits decreases survival. Nature 617, 599–607. issn: 1476-4687. doi:10.1038/s41586-023-06036-1
Ksendzovsky, A., Bachani, M., Altshuler, M., Walbridge, S., Mortazavi, A., Moyer, M., et al. (2022). Chronic neuronal activation leads to elevated lactate dehydrogenase A through the AMP-activated protein kinase/hypoxia-inducible factor-1α hypoxia pathway. Brain Commun. 5, fcac298. issn: 26321297. doi:10.1093/braincomms/fcac298
Kuokkanen, P. T., Ashida, G., Kraemer, A., McColgan, T., Funabiki, K., Wagner, H., et al. (2018). Sensory Processing: contribution of action potentials to the extracellular field potential in the nucleus laminaris of barn owl. J. Neurophysiology 119, 1422–1436. issn: 15221598. doi:10.1152/jn.00175.2017
Kutsy, R. L. (1999). Focal extratemporal epilepsy: clinical features, EEG patterns, and surgical approach. J. neurological Sci. 166, 1–15. issn: 0022-510X. doi:10.1016/s0022-510x(99)00107-0
Lachaux, J.-P., Rodriguez, E., Martinerie, J., and Varela, F. J. (1999). Measuring phase synchrony in brain signals. Hum. Brain Mapp. 8, 194–208. doi:10.1002/(sici)1097-0193(1999)8:4<194::aid-hbm4>3.0.co;2-c
Langdon, C., Genkin, M., and Engel, T. A. (2023). A unifying perspective on neural manifolds and circuits for cognition. Nat. Rev. Neurosci. 24, 363–377. issn: 1471-0048. doi:10.1038/s41583-023-00693-x
Latchoumane, C. F. V., Jackson, L., Sendi, M. S. E., Tehrani, K. F., Mortensen, L. J., Stice, S. L., et al. (2018). Chronic electrical stimulation promotes the excitability and plasticity of ESC-derived neurons following glutamate-induced inhibition in vitro. Sci. Rep. 8, 10957. issn: 20452322. doi:10.1038/s41598-018-29069-3
Leite, J. P., Neder, L., Arisi, G. M., Carlotti, C. G., Assirati, J. A., and Moreira, J. E. (2005). Plasticity, synaptic strength, and epilepsy: what can we learn from ultrastructural data? Epilepsia 46 (Suppl. 5), 134–141. issn: 0013-9580. doi:10.1111/j.1528-1167.2005.01021.x
Lévesque, M., and Avoli, M. (2013). The kainic acid model of temporal lobe epilepsy. Neurosci. Biobehav. Rev. 37, 2887–2899. issn: 1873-7528 (10 Pt. doi:10.1016/j.neubiorev.2013.10.011
Li, S., Hayden, E. Y., Garcia, V. J., Fuchs, D. T., Sheyn, J., Daley, D. A., et al. (2020). Activated bone marrow-derived macrophages eradicate alzheimer's-related Aβ42 oligomers and protect synapses. Front. Immunol. 11. issn: 1664-3224. doi:10.3389/fimmu.2020.00049
Lillis, K. P., Wang, Z., Mail, M., Zhao, G. Q., Berdichevsky, Y., Bacskai, B., et al. (2015). Evolution of network synchronization during early epileptogenesis parallels synaptic circuit alterations. J. Neurosci. 35, 9920–9934. issn: 0270-6474. doi:10.1523/JNEUROSCI.4007-14.2015
Luhmann, H. J., Dzhala, V. I., and Ben-Ari, Y. (2000). Generation and propagation of 4-AP-induced epileptiform activity in neonatal intact limbic structures in vitro. Eur. J. Neurosci. 12, 2757–2768. doi:10.1046/j.1460-9568.2000.00156.x
Marques, F. V. B. S., Kawamura, L. R. d. S. M., Quintanilha, M. V. T., Mota, I. F. L., Vasconcelos, A. S., da Silva, J. F. F., et al. (2022). Pharmacological perspectives and mechanisms involved in epileptogenesis. Beni-Suef Univ. J. Basic Appl. Sci. 11, 97–10. issn: 23148543. doi:10.1186/s43088-022-00278-0
Marshall, G. F., Gonzalez-Sulser, A., and Abbott, C. M. (2021). Modelling epilepsy in the mouse: challenges and solutions. Dis. Models and Mech. 14. issn: dmm047449-8411. doi:10.1242/dmm.047449
Matsuda, Y., Yoshida, S., Fujimura, K., and Nakamura, M. (1986). Depression of spike adaptation and afterhyperpolarization by 4-aminopyridine in hippocampal neurons. Neurosci. Lett. 65, 316–320. issn: 0304-3940. doi:10.1016/0304-3940(86)90281-8
Meng, F., You, Y., Liu, Z., Liu, J., Ding, H., and Xu, R. (2015). Neuronal calcium signaling pathways are associated with the development of epilepsy. Mol. Med. Rep. 11, 196–202. doi:10.3892/mmr.2014.2756
Mortazavi, A., Fayed, I., Bachani, M., Dowdy, T., Jahanipour, J., Khan, A., et al. (2022). IDH-mutated gliomas promote epileptogenesis through d-2-hydroxyglutarate-dependent mTOR hyperactivation. Neuro-Oncology 24, 1423–1435. issn: 1522-8517 (9 Sept. doi:10.1093/neuonc/noac003
Mossink, B., Verboven, A. H. A., van Hugte, E. J. H., Klein Gunnewiek, T. M., Parodi, G., Linda, K., et al. (2021). Human neuronal networks on micro-electrode arrays are a highly robust tool to study disease-specific genotype-phenotype correlations in vitro. Stem Cell Rep. 16, 2182–2196. issn: 22136711 (9 Sept. doi:10.1016/j.stemcr.2021.07.001
Neal, E. G., Maciver, S., Schoenberg, M. R., and Vale, F. L. (2020). Surgical disconnection of epilepsy network correlates with improved outcomes. Seizure 76, 56–63. issn: 1532-2688. doi:10.1016/j.seizure.2020.01.018
Neher, E., and Sakmann, B. (1976). Single-channel currents recorded from membrane of denervated frog muscle fibres. Nature 260, 799–802. issn: 0028-0836 (5554. doi:10.1038/260799a0
Ngugi, A. K., Bottomley, C., Kleinschmidt, I., Sander, J. W., and Newton, C. R. (2010). Estimation of the burden of active and life-time epilepsy: a meta-analytic approach. Epilepsia 51, 883–890. issn: 1528-1167. doi:10.1111/j.1528-1167.2009.02481.x
Nolta, N. F., Christensen, M. B., Crane, P. D., Skousen, J. L., and Tresco, P. A. (2015). BBB leakage, astrogliosis, and tissue loss correlate with silicon microelectrode array recording performance. Biomaterials 53, 753–762. issn: 1878-5905. doi:10.1016/j.biomaterials.2015.02.081
Oblasov, I., Idzhilova, O., Balaban, P., and Nikitin, E. (2024). Cell culture models for epilepsy research and treatment. Explor. Med. 5, 65–75. issn: 2692-3106. doi:10.37349/emed.2024.00206
Parks, T. W., and McClellan, J. H. (1972). Chebyshev approximation for nonrecursive digital filters with linear phase. IEEE Trans. Circuit Theory 19, 189–194. issn: 00189324 (2. doi:10.1109/tct.1972.1083419
Pea, F., and Tapia, R. (2000). Seizures and neurodegeneration induced by 4-aminopyridine in rat hippocampus in vivo: role of glutamate- and GABA-mediated neurotransmission and of ion channels. Neuroscience 101, 547–561. issn: 03064522. doi:10.1016/s0306-4522(00)00400-0
Pfisterer, U., Petukhov, V., Demharter, S., Meichsner, J., Thompson, J. J., Batiuk, M. Y., et al. (2020). Identification of epilepsy-associated neuronal subtypes and gene expression underlying epileptogenesis. Nat. Commun. 11. issn: 2041-1723 (1 Oct. 2020). doi:10.1038/s41467-020-18752-7
Piper, R. J., Richardson, R. M., Worrell, G., Carmichael, D. W., Baldeweg, T., Litt, B., et al. (2022). Towards network-guided neuromodulation for epilepsy. Brain 145, 3347–3362. issn: 14602156. doi:10.1093/brain/awac234
Raimondo, J. V., Heinemann, U., de Curtis, M., Goodkin, H. P., Dulla, C. G., Janigro, D., et al. (2017). 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. Epilepsia 58 (Suppl. 4), 40–52. issn: 1528-1167. doi:10.1111/epi.13901
Ren, L., Kucewicz, M. T., Cimbalnik, J., Matsumoto, J. Y., Brinkmann, B. H., Hu, W., et al. (2015). Gamma oscillations precede interictal epileptiform spikes in the seizure onset zone. Neurology 84, 602–608. issn: 1526632X. doi:10.1212/WNL.0000000000001234
Rijal, S., Corona, L., Perry, M. S., Tamilia, E., Madsen, J. R., Stone, S. S. D., et al. (2023). Functional connectivity discriminates epileptogenic states and predicts surgical outcome in children with drug resistant epilepsy. Sci. Rep. 13. doi:10.1038/s41598-023-36551-0
Sakmann, B., and Neher, E. (1984). Patch clamp techniques for studying ionic channels in excitable membranes. Annu. Rev. Physiology 46, 455–472. issn: 0066-4278. doi:10.1146/annurev.ph.46.030184.002323
Saleem, A., Santos, A. C., Aquilino, M. S., Sivitilli, A. A., Attisano, L., and Carlen, P. L. (2023). Modelling hyperexcitability in human cerebral cortical organoids: oxygen/glucose deprivation most effective stimulant. Heliyon 9, e14999. issn: 24058440. doi:10.1016/j.heliyon.2023.e14999
Scharfman, H. E., Kanner, A. M., Friedman, A., Blümcke, I., Crocker, C. E., Cendes, F., et al. (2018). Epilepsy as a network disorder: what can we learn from other network disorders such as dementia and schizophrenia, and what are the implications for translational research? Epilepsy and Behav. 78, 302–312. issn: 15255069. doi:10.1016/j.yebeh.2017.09.016
Segev, A., Garcia-Oscos, F., and Kourrich, S. (2016). Whole-cell patch-clamp recordings in brain slices. J. Vis. Exp., 54024. issn: 1940087X (112 June 2016). doi:10.3791/54024
Selinger, J. V., Pancrazio, J. J., and Gross, G. W. (2004). Measuring synchronization in neuronal networks for biosensor applications. Biosens. Bioelectron. 19, 675–683. issn: 0956-5663. doi:10.1016/s0956-5663(03)00267-7
Stark, E., and Abeles, M. (2007). Predicting movement from multiunit activity. J. Neurosci. 27, 8387–8394. issn: 0270-6474. doi:10.1523/JNEUROSCI.1321-07.2007
Sumadewi, K. T., Harkitasari, S., and Tjandra, D. C. (2023). Biomolecular mechanisms of epileptic seizures and epilepsy: a review. Acta Epileptol. 5, 28–22. issn: 25244434. doi:10.1186/s42494-023-00137-0
Taylor, P. N., Sinha, N., Wang, Y., Vos, S. B., de Tisi, J., Miserocchi, A., et al. (2018). The impact of epilepsy surgery on the structural connectome and its relation to outcome. NeuroImage Clin. 18, 202–214. issn: 2213-1582. doi:10.1016/j.nicl.2018.01.028
Tidball, A. M., and Parent, J. M. (2016). Concise review: exciting cells: modeling genetic epilepsies with patient-derived induced pluripotent stem cells. Stem Cells 34, 27–33. issn: 1549-4918. doi:10.1002/stem.2203
Wainger, B. J., Kiskinis, E., Mellin, C., Wiskow, O., Han, S. S. W., Sandoe, J., et al. (2014). Intrinsic membrane hyperexcitability of amyotrophic lateral sclerosis patient-derived motor neurons. Cell Rep. 7 (1), 1–11. issn: 22111247. doi:10.1016/j.celrep.2014.03.019
Waldert, S., Lemon, R. N., and Kraskov, A. (2013). Influence of spiking activity on cortical local field potentials. J. Physiology 591, 5291–5303. issn: 00223751. doi:10.1113/jphysiol.2013.258228
Wong, M. (2011). Epilepsy in a dish: an in vitro model of epileptogenesis. Epilepsy Curr. 11, 153–154. issn: 15357597 (5 Sept. doi:10.5698/1535-7511-11.5.153
Woodard, C. M., Campos, B. A., Kuo, S. H., Nirenberg, M. J., Nestor, M. W., Zimmer, M., et al. (2014). iPSC-derived dopamine neurons reveal differences between monozygotic twins discordant for Parkinson’s disease. Cell Rep. 9, 1173–1182. issn: 2211-1247. doi:10.1016/j.celrep.2014.10.023
Yamaguchi, S.-I., and Rogawski, M. A. (1992). Effects of anticonvulsant drugs on 4-aminopyridine-induced seizures in mice. Epilepsy Res. 11, 9–16. doi:10.1016/0920-1211(92)90016-m
Yokoi, R., Shigemoto-Kuroda, T., Matsuda, N., Odawara, A., and Suzuki, I. (2022). Electrophysiological responses to seizurogenic compounds dependent on E/I balance in human iPSC-derived cortical neural networks. J. Pharmacol. Sci. 148, 267–278. issn: 1347-8613. doi:10.1016/j.jphs.2021.12.006
Keywords: 4-aminopyridine, epilepsy, epilepsy model, functional connectivity, microelectrode arrays, network changes, neuronal networks, seizure networks
Citation: Boddeti U, Langbein J, McAfee D, Altshuler M, Bachani M, Zaveri HP, Spencer D, Zaghloul KA and Ksendzovsky A (2024) Modeling seizure networks in neuron-glia cultures using microelectrode arrays. Front. Netw. Physiol. 4:1441345. doi: 10.3389/fnetp.2024.1441345
Received: 30 May 2024; Accepted: 20 August 2024;
Published: 03 September 2024.
Edited by:
Peter A. Tass, Stanford University, United StatesReviewed by:
Amrit Kashyap, Charité University Medicine Berlin, GermanyDaniel Ehrens, Stanford University, United States
Copyright © 2024 Boddeti, Langbein, McAfee, Altshuler, Bachani, Zaveri, Spencer, Zaghloul and Ksendzovsky. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Alexander Ksendzovsky, aksendzovsky@som.umaryland.edu