Quantitative examination of stimulus-response relations in cortical networks in vitro

Activity-dependent recruitment of inhibition and excitation in the awake mammalian cortex during electrical stimulation

Electrical stimulation is an effective tool for mapping and altering brain connectivity, with applications ranging from treating pharmacology-resistant neurological disorders to providing sensory feedback for neural prostheses. Paramount to the success of these applications is the ability to manipulate electrical currents to precisely control evoked neural activity patterns. However, little is known about stimulation-evoked responses in inhibitory neurons nor how stimulation-evoked activity patterns depend on ongoing neural activity. In this study, we used 2-photon imaging and cell-type specific labeling to measure single-cell responses of excitatory and inhibitory neurons to electrical stimuli in the visual cortex of awake mice. Our data revealed strong interactions between electrical stimulation and pre-stimulus activity of single neurons in awake animals and distinct recruitment and response patterns for excitatory and inhibitory neurons. This work demonstrates the importance of cell-type-specific labeling of neurons in future studies.

Responses of Cortical Neurons to Intracortical Microstimulation in Awake Primates

Intracortical microstimulation (ICMS) is commonly used in many experimental and clinical paradigms; however, its effects on the activation of neurons are still not completely understood. To document the responses of cortical neurons in awake nonhuman primates to stimulation, we recorded single-unit activity while delivering single-pulse stimulation via Utah arrays implanted in primary motor cortex (M1) of three macaque monkeys. Stimuli between 5 and 50 μA delivered to single channels reliably evoked spikes in neurons recorded throughout the array with delays of up to 12 ms. ICMS pulses also induced a period of inhibition lasting up to 150 ms that typically followed the initial excitatory response. Higher current amplitudes led to a greater probability of evoking a spike and extended the duration of inhibition. The likelihood of evoking a spike in a neuron was dependent on the spontaneous firing rate as well as the delay between its most recent spike time and stimulus onset. Tonic repetitive stimulation between 2 and 20 Hz often modulated both the probability of evoking spikes and the duration of inhibition; high-frequency stimulation was more likely to change both responses. On a trial-by-trial basis, whether a stimulus evoked a spike did not affect the subsequent inhibitory response; however, their changes over time were often positively or negatively correlated. Our results document the complex dynamics of cortical neural responses to electrical stimulation that need to be considered when using ICMS for scientific and clinical applications.

Nanomaterial-based microelectrode arrays for in vitro bidirectional brain-computer interfaces: a review

A bidirectional in vitro brain-computer interface (BCI) directly connects isolated brain cells with the surrounding environment, reads neural signals and inputs modulatory instructions. As a noninvasive BCI, it has clear advantages in understanding and exploiting advanced brain function due to the simplified structure and high controllability of ex vivo neural networks. However, the core of ex vivo BCIs, microelectrode arrays (MEAs), urgently need improvements in the strength of signal detection, precision of neural modulation and biocompatibility. Notably, nanomaterial-based MEAs cater to all the requirements by converging the multilevel neural signals and simultaneously applying stimuli at an excellent spatiotemporal resolution, as well as supporting long-term cultivation of neurons. This is enabled by the advantageous electrochemical characteristics of nanomaterials, such as their active atomic reactivity and outstanding charge conduction efficiency, improving the performance of MEAs. Here, we review the fabrication of nanomaterial-based MEAs applied to bidirectional in vitro BCIs from an interdisciplinary perspective. We also consider the decoding and coding of neural activity through the interface and highlight the various usages of MEAs coupled with the dissociated neural cultures to benefit future developments of BCIs.

Comparative microelectrode array data of the functional development of hPSC-derived and rat neuronal networks

We present a dataset of microelectrode array (MEA) recordings from human pluripotent stem cell (hPSC)-derived and rat embryonic cortical neurons during their in vitro maturation. The data were prepared to assess extracellularly recorded spontaneous activity and to compare the functional development of these neuronal networks. In addition to recordings of spontaneous activity, we provide pharmacological responses of hPSC-derived and rat cortical cultures at their mature stage. Together with the recorded electrode raw data, we share the analysis code to form a comprehensive dataset including spike times, spike waveforms, burst activity and network synchronization metrics calculated with two different connectivity estimators. Moreover, we provide the analysis code that produced the key scientific findings published previously with this dataset. This large dataset enables investigation of the functional aspects of maturing cortical neuronal networks and provides substantial parameters to assess the differences and similarities between hPSC-derived and rat cortical networks in vitro. This publicly available dataset will be beneficial, especially for experimental and computational neuroscientists.

Measurement(s)	extracellular electrophysiological activity
Technology Type(s)	microelectrode array
Factor Type(s)	days in vitro • dosage
Sample Characteristic - Organism	hPSC-derived neuronal cultures • rat neuronal cultures
Sample Characteristic - Environment	functional development • pharmacological response

Modulation of dynamics in a pre-existing hippocampal network by neural stem cells on a microelectrode array

Objective.Neural stem cells (NSCs) are continuously produced throughout life in the hippocampus, which is a vital structure for learning and memory. NSCs in the brain incorporate into the functional hippocampal circuits and contribute to processing information. However, little is known about the mechanisms of NSCs' activity in a pre-existing neuronal network. Here, we investigate the role of NSCs in the neuronal activity of a pre-existing hippocampalin vitronetwork grown on microelectrode arrays.Approach.We assessed the change in internal dynamics of the network by additional NSCs based on spontaneous activity. We also evaluated the networks' ability to discriminate between different input patterns by measuring evoked activity in response to external inputs.Main results.Analysis of spontaneous activity revealed that additional NSCs prolonged network bursts with longer intervals, generated a lower number of initiating patterns, and decreased synchronization among neurons. Moreover, the network with NSCs showed higher synchronicity in close connections among neurons responding to external inputs and a larger difference in spike counts and cross-correlations during evoked response between two different inputs. Taken together, our results suggested that NSCs alter the internal dynamics of the pre-existing hippocampal network and produce more specific responses to external inputs, thus enhancing the ability of the network to differentiate two different inputs.Significance.We demonstrated that NSCs improve the ability to distinguish external inputs by modulating the internal dynamics of a pre-existing network in a hippocampal culture. Our results provide novel insights into the relationship between NSCs and learning and memory.

Characterization of neuronal viability and network activity under microfluidic flow

BACKGROUND

Microfluidics technology has the potential to allow precise control of the temporal and spatial aspects of solute concentration, making it highly relevant for the study of volume transmission mechanisms in neural tissue. However, full utilization of this technology depends on understanding how microfluidic flow at the rates needed for rapid solution exchange affects neuronal viability and network activity.

NEW METHOD

We designed a tape-based pressurized microfluidic flow system that is simple to fabricate and can be attached to commercial microelectrode arrays. The device is multi-layered, allowing the inclusion of a porous polycarbonate membrane to isolate neuronal cultures from shear forces while maintaining diffusive exchange of solutes. We used this system to investigate how flow affected survival and spiking patterns of cultured hippocampal neurons.

RESULTS

Viability and network activity of the cultures were reduced in proportion to flow rate. However, shear reduction measures did not improve survival or spiking activity; media conditioning in conjunction with culture age proved to be the critical factors for network stability. Diffusion simulations indicate that dilution of a small molecule accounts for the deleterious effects of flow on neuronal cultures.

COMPARISON WITH EXISTING METHODS

This work establishes the experimental conditions for real time measurement of network activity during rapid solution exchange, using multi-layered chambers with reversible bonding that allow for reuse of microelectrode arrays.

CONCLUSIONS

With correct media conditioning, the microfluidic flow system allows drug delivery on a subsecond timescale without disruption of network activity or viability, enabling in vitro reproduction of volume transmission mechanisms.

Unique Features of Network Bursts Emerge From the Complex Interplay of Excitatory and Inhibitory Receptors in Rat Neocortical Networks

Spontaneous network activity plays a fundamental role in the formation of functional networks during early development. The landmark of this activity is the recurrent emergence of intensive time-limited network bursts (NBs) rapidly spreading across the entire dissociated culture in vitro. The main excitatory mediators of NBs are glutamatergic alpha-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptors (AMPARs) and N-Methyl-D-aspartic-acid receptors (NMDARs) that express fast and slow ion channel kinetics, respectively. The fast inhibition of the activity is mediated through gamma-aminobutyric acid type A receptors (GABA_ARs). Although the AMPAR, NMDAR and GABA_AR kinetics have been biophysically characterized in detail at the monosynaptic level in a variety of brain areas, the unique features of NBs emerging from the kinetics and the complex interplay of these receptors are not well understood. The goal of this study is to analyze the contribution of fast GABA_ARs on AMPAR- and NMDAR- mediated spontaneous NB activity in dissociated neonatal rat cortical cultures at 3 weeks in vitro. The networks were probed by both acute and gradual application of each excitatory receptor antagonist and combinations of acute excitatory and inhibitory receptor antagonists. At the same time, the extracellular network-wide activity was recorded with microelectrode arrays (MEAs). We analyzed the characteristic NB measures extracted from NB rate profiles and the distributions of interspike intervals, interburst intervals, and electrode recruitment time as well as the similarity of spatio-temporal patterns of network activity under different receptor antagonists. We show that NBs were rapidly initiated and recruited as well as diversely propagated by AMPARs and temporally and spatially maintained by NMDARs. GABA_ARs reduced the spiking frequency in AMPAR-mediated networks and dampened the termination of NBs in NMDAR-mediated networks as well as slowed down the recruitment of activity in all networks. Finally, we show characteristic super bursts composed of slow NBs with highly repetitive spatio-temporal patterns in gradually AMPAR blocked networks. To the best of our knowledge, this study is the first to unravel in detail how the three main mediators of synaptic transmission uniquely shape the NB characteristics, such as the initiation, maintenance, recruitment and termination of NBs in cortical cell cultures in vitro.

Active High-Density Electrode Arrays: Technology and Applications in Neuronal Cell Cultures

Active high-density electrode arrays realized with complementary metal-oxide-semiconductor (CMOS) technology provide electrophysiological recordings from several thousands of closely spaced microelectrodes. This has drastically advanced the spatiotemporal recording resolution of conventional multielectrode arrays (MEAs). Thus, today's electrophysiology in neuronal cultures can exploit label-free electrical readouts from a large number of single neurons within the same network. This provides advanced capabilities to investigate the properties of self-assembling neuronal networks, to advance studies on neurotoxicity and neurodevelopmental alterations associated with human brain diseases, and to develop cell culture models for testing drug- or cell-based strategies for therapies.Here, after introducing the reader to this neurotechnology, we summarize the results of different recent studies demonstrating the potential of active high-density electrode arrays for experimental applications. We also discuss ongoing and possible future research directions that might allow for moving these platforms forward for screening applications.

Inhomogeneities in Network Structure and Excitability Govern Initiation and Propagation of Spontaneous Burst Activity

The mesoscale architecture of neuronal networks strongly influences the initiation of spontaneous activity and its pathways of propagation. Spontaneous activity has been studied extensively in networks of cultured cortical neurons that generate complex yet reproducible patterns of synchronous bursting events that resemble the activity dynamics in developing neuronal networks in vivo. Synchronous bursts are mostly thought to be triggered at burst initiation sites due to build-up of noise or by highly active neurons, or to reflect reverberating activity that circulates within larger networks, although neither of these has been observed directly. Inferring such collective dynamics in neuronal populations from electrophysiological recordings crucially depends on the spatial resolution and sampling ratio relative to the size of the networks assessed. Using large-scale microelectrode arrays with 1024 electrodes at 0.3 mm pitch that covered the full extent of in vitro networks on about 1 cm², we investigated where bursts of spontaneous activity arise and how their propagation patterns relate to the regions of origin, the network's structure, and to the overall distribution of activity. A set of alternating burst initiation zones (BIZ) dominated the initiation of distinct bursting events and triggered specific propagation patterns. Moreover, BIZs were typically located in areas with moderate activity levels, i.e., at transitions between hot and cold spots. The activity-dependent alternation between these zones suggests that the local networks forming the dominating BIZ enter a transient depressed state after several cycles (similar to Eytan et al., 2003), allowing other BIZs to take over temporarily. We propose that inhomogeneities in the network structure define such BIZs and that the depletion of local synaptic resources limit repetitive burst initiation.

The effect of pharmacological inhibition of Serine Proteases on neuronal networks in vitro

Neurons are embedded in an extracellular matrix (ECM), which functions both as a scaffold and as a regulator of neuronal function. The ECM is in turn dynamically altered through the action of serine proteases, which break down its constituents. This pathway has been implicated in the regulation of synaptic plasticity and of neuronal intrinsic excitability. In this study, we determined the short-term effects of interfering with proteolytic processes in the ECM, with a newly developed serine protease inhibitor. We monitored the spontaneous electrophysiological activity of in vitro primary rat cortical cultures, using microelectrode arrays. While pharmacological inhibition at a low dosage had no significant effect, at elevated concentrations it altered significantly network synchronization and functional connectivity but left unaltered single-cell electrical properties. These results suggest that serine protease inhibition affects synaptic properties, likely through its actions on the ECM.

Long-Term Activity Dynamics of Single Neurons and Networks

The firing rate of neuronal spiking in vitro and in vivo significantly varies over extended timescales, characterized by long-memory processes and complex statistics, and appears in spontaneous as well as evoked activity upon repeated stimulus presentation. These variations in response features and their statistics, in face of repeated instances of a given physical input, are ubiquitous in all levels of brain-behavior organization. They are expressed in single neuron and network response variability but even appear in variations of subjective percepts or psychophysical choices and have been described as stemming from history-dependent, stochastic, or rate-determined processes.But what are the sources underlying these temporally rich variations in firing rate? Are they determined by interactions of the nervous system as a whole, or do isolated, single neurons or neuronal networks already express these fluctuations independent of higher levels? These questions motivated the application of a method that allows for controlled and specific long-term activation of a single neuron or neuronal network, isolated from higher levels of cortical organization.This chapter highlights the research done in cultured cortical networks to study (1) the inherent non-stationarity of neuronal network activity, (2) single neuron response fluctuations and underlying processes, and (3) the interface layer between network and single cell, the non-stationary efficacy of the ensemble of synapses impinging onto the observed neuron.

State-dependent representation of stimulus-evoked activity in high-density recordings of neural cultures

Neuronal responses to external stimuli vary from trial to trial partly because they depend on continuous spontaneous variations of the state of neural circuits, reflected in variations of ongoing activity prior to stimulus presentation. Understanding how post-stimulus responses relate to the pre-stimulus spontaneous activity is thus important to understand how state dependence affects information processing and neural coding, and how state variations can be discounted to better decode single-trial neural responses. Here we exploited high-resolution CMOS electrode arrays to record simultaneously from thousands of electrodes in in-vitro cultures stimulated at specific sites. We used information-theoretic analyses to study how ongoing activity affects the information that neuronal responses carry about the location of the stimuli. We found that responses exhibited state dependence on the time between the last spontaneous burst and the stimulus presentation and that the dependence could be described with a linear model. Importantly, we found that a small number of selected neurons carry most of the stimulus information and contribute to the state-dependent information gain. This suggests that a major value of large-scale recording is that it individuates the small subset of neurons that carry most information and that benefit the most from knowledge of its state dependence.

Quantitative evaluation of extrinsic factors influencing electrical excitability in neuronal networks: Voltage Threshold Measurement Method (VTMM)

The electrical excitability of neural networks is influenced by different environmental factors. Effective and simple methods are required to objectively and quantitatively evaluate the influence of such factors, including variations in temperature and pharmaceutical dosage. The aim of this paper was to introduce ‘the voltage threshold measurement method’, which is a new method using microelectrode arrays that can quantitatively evaluate the influence of different factors on the electrical excitability of neural networks. We sought to verify the feasibility and efficacy of the method by studying the effects of acetylcholine, ethanol, and temperature on hippocampal neuronal networks and hippocampal brain slices. First, we determined the voltage of the stimulation pulse signal that elicited action potentials in the two types of neural networks under normal conditions. Second, we obtained the voltage thresholds for the two types of neural networks under different concentrations of acetylcholine, ethanol, and different temperatures. Finally, we obtained the relationship between voltage threshold and the three influential factors. Our results indicated that the normal voltage thresholds of the hippocampal neuronal network and hippocampal slice preparation were 56 and 31 mV, respectively. The voltage thresholds of the two types of neural networks were inversely proportional to acetylcholine concentration, and had an exponential dependency on ethanol concentration. The curves of the voltage threshold and the temperature of the medium for the two types of neural networks were U-shaped. The hippocampal neuronal network and hippocampal slice preparations lost their excitability when the temperature of the medium decreased below 34 and 33°C or increased above 42 and 43°C, respectively. These results demonstrate that the voltage threshold measurement method is effective and simple for examining the performance/excitability of neuronal networks.

Stimulation triggers endogenous activity patterns in cultured cortical networks

Cultures of dissociated cortical neurons represent a powerful trade-off between more realistic experimental models and abstract modeling approaches, allowing to investigate mechanisms of synchronized activity generation. These networks spontaneously alternate periods of high activity (i.e. network bursts) with periods of quiescence in a dynamic state which recalls the fluctuation of in vivo UP and DOWN states. Network bursts can also be elicited by external stimulation and their spatial propagation patterns tracked by means of multi-channel micro-electrode arrays. In this study, we used rat cortical cultures coupled to micro-electrode arrays to investigate the similarity between spontaneous and evoked activity patterns. We performed experiments by applying electrical stimulation to different network locations and demonstrated that the rank orders of electrodes during evoked and spontaneous events are remarkably similar independently from the stimulation source. We linked this result to the capability of stimulation to evoke firing in highly active and “leader” sites of the network, reliably and rapidly recruited within both spontaneous and evoked bursts. Our study provides the first evidence that spontaneous and evoked activity similarity is reliably observed also in dissociated cortical networks.

Investigating the impact of electrical stimulation temporal distribution on cortical network responses

Background

The brain is continuously targeted by a wealth of stimuli with complex spatio-temporal patterns and has presumably evolved in order to cope with those inputs in an optimal way. Previous studies investigating the response capabilities of either single neurons or intact sensory systems to external stimulation demonstrated that stimuli temporal distribution is an important, if often overlooked, parameter.

Results

In this study we investigated how cortical networks plated over micro-electrode arrays respond to different stimulation sequences in which inter-pulse intervals followed a 1/f^β distribution, for different values of β ranging from 0 to ∞. Cross-correlation analysis revealed that network activity preferentially synchronizes with external input sequences featuring β closer to 1 and, in any case, never for regular (i.e. fixed-frequency) stimulation sequences. We then tested the interplay between different average stimulation frequencies (based on the intrinsic firing/bursting frequency of the network) for two selected values of β, i.e. 1 (scale free) and ∞ (regular). In general, we observed no preference for stimulation frequencies matching the endogenous rhythms of the network. Moreover, we found that in case of regular stimulation the capability of the network to follow the stimulation sequence was negatively correlated to the absolute stimulation frequency, whereas using scale-free stimulation cross-correlation between input and output sequences was independent from average input frequency.

Conclusions

Our results point out that the preference for a scale-free distribution of the stimuli is observed also at network level and should be taken into account in designing more efficient protocols for neuromodulation purposes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12868-017-0366-z) contains supplementary material, which is available to authorized users.

Mesoscale Architecture Shapes Initiation and Richness of Spontaneous Network Activity

Spontaneous activity in the absence of external input, including propagating waves of activity, is a robust feature of neuronal networks in vivo and in vitro The neurophysiological and anatomical requirements for initiation and persistence of such activity, however, are poorly understood, as is their role in the function of neuronal networks. Computational network studies indicate that clustered connectivity may foster the generation, maintenance, and richness of spontaneous activity. Since this mesoscale architecture cannot be systematically modified in intact tissue, testing these predictions is impracticable in vivo Here, we investigate how the mesoscale structure shapes spontaneous activity in generic networks of rat cortical neurons in vitro In these networks, neurons spontaneously arrange into local clusters with high neurite density and form fasciculating long-range axons. We modified this structure by modulation of protein kinase C, an enzyme regulating neurite growth and cell migration. Inhibition of protein kinase C reduced neuronal aggregation and fasciculation of axons, i.e., promoted uniform architecture. Conversely, activation of protein kinase C promoted aggregation of neurons into clusters, local connectivity, and bundling of long-range axons. Supporting predictions from theory, clustered networks were more spontaneously active and generated diverse activity patterns. Neurons within clusters received stronger synaptic inputs and displayed increased membrane potential fluctuations. Intensified clustering promoted the initiation of synchronous bursting events but entailed incomplete network recruitment. Moderately clustered networks appear optimal for initiation and propagation of diverse patterns of activity. Our findings support a crucial role of the mesoscale architectures in the regulation of spontaneous activity dynamics.SIGNIFICANCE STATEMENT Computational studies predict richer and persisting spatiotemporal patterns of spontaneous activity in neuronal networks with neuron clustering. To test this, we created networks of varying architecture in vitro Supporting these predictions, the generation and spatiotemporal patterns of propagation were most variable in networks with intermediate clustering and lowest in uniform networks. Grid-like clustering, on the other hand, facilitated spontaneous activity but led to degenerating patterns of propagation. Neurons outside clusters had weaker synaptic input than neurons within clusters, in which increased membrane potential fluctuations facilitated the initiation of synchronized spike activity. Our results thus show that the intermediate level organization of neuronal networks strongly influences the dynamics of their activity.

Paper-based patterned 3D neural cultures as a tool to study network activity on multielectrode arrays

High-throughput platform targeting activity patterns of 3D neural cultures with arbitrary topology, by combining network-wide intracellular and local extracellular signals.

Spectral Entropy Based Neuronal Network Synchronization Analysis Based on Microelectrode Array Measurements

Synchrony and asynchrony are essential aspects of the functioning of interconnected neuronal cells and networks. New information on neuronal synchronization can be expected to aid in understanding these systems. Synchronization provides insight in the functional connectivity and the spatial distribution of the information processing in the networks. Synchronization is generally studied with time domain analysis of neuronal events, or using direct frequency spectrum analysis, e.g., in specific frequency bands. However, these methods have their pitfalls. Thus, we have previously proposed a method to analyze temporal changes in the complexity of the frequency of signals originating from different network regions. The method is based on the correlation of time varying spectral entropies (SEs). SE assesses the regularity, or complexity, of a time series by quantifying the uniformity of the frequency spectrum distribution. It has been previously employed, e.g., in electroencephalogram analysis. Here, we revisit our correlated spectral entropy method (CorSE), providing evidence of its justification, usability, and benefits. Here, CorSE is assessed with simulations and in vitro microelectrode array (MEA) data. CorSE is first demonstrated with a specifically tailored toy simulation to illustrate how it can identify synchronized populations. To provide a form of validation, the method was tested with simulated data from integrate-and-fire model based computational neuronal networks. To demonstrate the analysis of real data, CorSE was applied on in vitro MEA data measured from rat cortical cell cultures, and the results were compared with three known event based synchronization measures. Finally, we show the usability by tracking the development of networks in dissociated mouse cortical cell cultures. The results show that temporal correlations in frequency spectrum distributions reflect the network relations of neuronal populations. In the simulated data, CorSE unraveled the synchronizations. With the real in vitro MEA data, CorSE produced biologically plausible results. Since CorSE analyses continuous data, it is not affected by possibly poor spike or other event detection quality. We conclude that CorSE can reveal neuronal network synchronization based on in vitro MEA field potential measurements. CorSE is expected to be equally applicable also in the analysis of corresponding in vivo and ex vivo data analysis.

Extracellularly Recorded Somatic and Neuritic Signal Shapes and Classification Algorithms for High-Density Microelectrode Array Electrophysiology

High-density microelectrode arrays (HDMEA) have been recently introduced to study principles of neural function at high spatial resolution. However, the exact nature of the experimentally observed extracellular action potentials (EAPs) is still incompletely understood. The soma, axon and dendrites of a neuron can all exhibit regenerative action potentials that could be sensed with HDMEA electrodes. Here, we investigate the contribution of distinct neuronal sources of activity in HDMEA recordings from low-density neuronal cultures. We recorded EAPs with HDMEAs having 11,011 electrodes and then fixed and immunostained the cultures with β3-tubulin for high-resolution fluorescence imaging. Immunofluorescence images overlaid with the activity maps showed EAPs both at neuronal somata and distal neurites. Neuritic EAPs had mostly narrow triphasic shapes, consisting of a positive, a pronounced negative peak and a second positive peak. EAPs near somata had wide monophasic or biphasic shapes with a main negative peak, and following optional positive peak. We show that about 86% of EAP recordings consist of somatic spikes, while the remaining 14% represent neuritic spikes. Furthermore, the adaptation of the waveform shape during bursts of these neuritic spikes suggested that they originate from axons, rather than from dendrites. Our study improves the understanding of HDMEA signals and can aid in the identification of the source of EAPs.

Autonomous Optimization of Targeted Stimulation of Neuronal Networks

Driven by clinical needs and progress in neurotechnology, targeted interaction with neuronal networks is of increasing importance. Yet, the dynamics of interaction between intrinsic ongoing activity in neuronal networks and their response to stimulation is unknown. Nonetheless, electrical stimulation of the brain is increasingly explored as a therapeutic strategy and as a means to artificially inject information into neural circuits. Strategies using regular or event-triggered fixed stimuli discount the influence of ongoing neuronal activity on the stimulation outcome and are therefore not optimal to induce specific responses reliably. Yet, without suitable mechanistic models, it is hardly possible to optimize such interactions, in particular when desired response features are network-dependent and are initially unknown. In this proof-of-principle study, we present an experimental paradigm using reinforcement-learning (RL) to optimize stimulus settings autonomously and evaluate the learned control strategy using phenomenological models. We asked how to (1) capture the interaction of ongoing network activity, electrical stimulation and evoked responses in a quantifiable ‘state’ to formulate a well-posed control problem, (2) find the optimal state for stimulation, and (3) evaluate the quality of the solution found. Electrical stimulation of generic neuronal networks grown from rat cortical tissue in vitro evoked bursts of action potentials (responses). We show that the dynamic interplay of their magnitudes and the probability to be intercepted by spontaneous events defines a trade-off scenario with a network-specific unique optimal latency maximizing stimulus efficacy. An RL controller was set to find this optimum autonomously. Across networks, stimulation efficacy increased in 90% of the sessions after learning and learned latencies strongly agreed with those predicted from open-loop experiments. Our results show that autonomous techniques can exploit quantitative relationships underlying activity-response interaction in biological neuronal networks to choose optimal actions. Simple phenomenological models can be useful to validate the quality of the resulting controllers.

Electrical stimulation of the brain is increasingly used to alleviate the symptoms of a range of neurological disorders and as a means to artificially inject information into neural circuits in neuroprosthetic applications. Machine learning has been proposed to find optimal stimulation settings autonomously. However, this approach is impeded by the complexity of the interaction between the stimulus and the activity of the network, which makes it difficult to test how good the result actually is. We used phenomenological models of the interaction between stimulus and spontaneous activity in a neuronal network to design a testable machine learning challenge and evaluate the quality of the solution found by the algorithm. In this task, the learning algorithm had to find a solution that balances competing interdependencies of ongoing neuronal activity with opposing effects on the efficacy of stimulation. We show that machine learning can successfully solve this task and that the solutions found are close to the optimal settings to maximize the efficacy of stimulation. Since the paradigm involves several typical problems found in other settings, such concepts could help to formalize machine learning problems in more complex biological networks and to test the quality of their performance.

A comparison of computational methods for detecting bursts in neuronal spike trains and their application to human stem cell-derived neuronal networks

We provide an unbiased quantitative assessment of eight existing methods for identifying bursts in neuronal spike trains. We reveal limitations in a number of commonly used burst detection techniques and provide recommendations for the best practice for accurate identification of bursts using existing techniques. An analysis of the ontogeny of bursting activity in a novel data set of recordings from human induced pluripotent stem cell-derived neuronal networks, using the highest-performing burst detectors from our study, is also presented.

Accurate identification of bursting activity is an essential element in the characterization of neuronal network activity. Despite this, no one technique for identifying bursts in spike trains has been widely adopted. Instead, many methods have been developed for the analysis of bursting activity, often on an ad hoc basis. Here we provide an unbiased assessment of the effectiveness of eight of these methods at detecting bursts in a range of spike trains. We suggest a list of features that an ideal burst detection technique should possess and use synthetic data to assess each method in regard to these properties. We further employ each of the methods to reanalyze microelectrode array (MEA) recordings from mouse retinal ganglion cells and examine their coherence with bursts detected by a human observer. We show that several common burst detection techniques perform poorly at analyzing spike trains with a variety of properties. We identify four promising burst detection techniques, which are then applied to MEA recordings of networks of human induced pluripotent stem cell-derived neurons and used to describe the ontogeny of bursting activity in these networks over several months of development. We conclude that no current method can provide “perfect” burst detection results across a range of spike trains; however, two burst detection techniques, the MaxInterval and logISI methods, outperform compared with others. We provide recommendations for the robust analysis of bursting activity in experimental recordings using current techniques.

Quantification and automatized adaptive detection of in vivo and in vitro neuronal bursts based on signal complexity

In this paper, we propose employing entropy values to quantify action potential bursts in electrophysiological measurements from the brain and neuronal cultures. Conventionally in the electrophysiological signal analysis, bursts are quantified by means of conventional measures such as their durations, and number of spikes in bursts. Here our main aim is to device metrics for burst quantification to provide for enhanced burst characterization. Entropy is a widely employed measure to quantify regularity/complexity of time series. Specifically, we investigate the applicability and differences of spectral entropy and sample entropy in the quantification of bursts in in vivo rat hippocampal measurements and in in vitro dissociated rat cortical cell culture measurement done with microelectrode arrays. For the task, an automatized and adaptive burst detection method is also utilized. Whereas the employed metrics are known from other applications, they are rarely employed in the assessment of burst in electrophysiological field potential measurements. Our results show that the proposed metrics are potential for the task at hand.

Universality, complexity and the praxis of biology: Two case studies

The phenomenon of biology provides a prime example for a naturally occurring complex system. The approach to this complexity reflects the tension between a reductionist, reverse-engineering stance, and more abstract, systemic ones. Both of us are reductionists, but our observations challenge reductionism, at least the naive version of it. Here we describe the challenge, focusing on two universal characteristics of biological complexity: two-way microscopic-macroscopic degeneracy, and lack of time scale separation within and between levels of organization. These two features and their consequences for the praxis of experimental biology, reflect inherent difficulties in separating the dynamics of any given level of organization from the coupled dynamics of all other levels, including the environment within which the system is embedded. Where these difficulties are not deeply acknowledged, the impacts of fallacies that are inherent to naive reductionism are significant. In an era where technology enables experimental high-resolution access to numerous observables, the challenge faced by the mature reductionist-identification of relevant microscopic variables-becomes more demanding than ever. The demonstrations provided here are taken from two very different biological realizations: populations of microorganisms and populations of neurons, thus making the lesson potentially general.

Shaping Neuronal Network Activity by Presynaptic Mechanisms

Neuronal microcircuits generate oscillatory activity, which has been linked to basic functions such as sleep, learning and sensorimotor gating. Although synaptic release processes are well known for their ability to shape the interaction between neurons in microcircuits, most computational models do not simulate the synaptic transmission process directly and hence cannot explain how changes in synaptic parameters alter neuronal network activity. In this paper, we present a novel neuronal network model that incorporates presynaptic release mechanisms, such as vesicle pool dynamics and calcium-dependent release probability, to model the spontaneous activity of neuronal networks. The model, which is based on modified leaky integrate-and-fire neurons, generates spontaneous network activity patterns, which are similar to experimental data and robust under changes in the model's primary gain parameters such as excitatory postsynaptic potential and connectivity ratio. Furthermore, it reliably recreates experimental findings and provides mechanistic explanations for data obtained from microelectrode array recordings, such as network burst termination and the effects of pharmacological and genetic manipulations. The model demonstrates how elevated asynchronous release, but not spontaneous release, synchronizes neuronal network activity and reveals that asynchronous release enhances utilization of the recycling vesicle pool to induce the network effect. The model further predicts a positive correlation between vesicle priming at the single-neuron level and burst frequency at the network level; this prediction is supported by experimental findings. Thus, the model is utilized to reveal how synaptic release processes at the neuronal level govern activity patterns and synchronization at the network level.

The activity of neuronal networks underlies basic neural functions such as sleep, learning and sensorimotor gating. Computational models of neuronal networks have been developed to capture the complexity of the network activity and predict how neuronal networks generate spontaneous activity. However, most computational models do not simulate the intricate synaptic release process that governs the interaction between neurons and has been shown to significantly impact neuronal network activity and animal behavior, learning and memory. Our paper demonstrates the importance of simulating the elaborate synaptic release process to understand how neuronal networks generate spontaneous activity and respond to manipulations of the release process. The model provides mechanistic explanations and predictions for experimental pharmacological and genetic manipulations. Thus, the model presents a novel computational platform to understand how mechanistic changes in the synaptic release process modulate network oscillatory activity that might impact basic neural functions.

Interaction of electrically evoked activity with intrinsic dynamics of cultured cortical networks with and without functional fast GABAergic synaptic transmission

The modulation of neuronal activity by means of electrical stimulation is a successful therapeutic approach for patients suffering from a variety of central nervous system disorders. Prototypic networks formed by cultured cortical neurons represent an important model system to gain general insights in the input–output relationships of neuronal tissue. These networks undergo a multitude of developmental changes during their maturation, such as the excitatory–inhibitory shift of the neurotransmitter GABA. Very few studies have addressed how the output properties to a given stimulus change with ongoing development. Here, we investigate input–output relationships of cultured cortical networks by probing cultures with and without functional GABA_Aergic synaptic transmission with a set of stimulation paradigms at various stages of maturation. On the cellular level, low stimulation rates (<15 Hz) led to reliable neuronal responses; higher rates were increasingly ineffective. Similarly, on the network level, lowest stimulation rates (<0.1 Hz) lead to maximal output rates at all ages, indicating a network wide refractory period after each stimulus. In cultures aged 3 weeks and older, a gradual recovery of the network excitability within tens of milliseconds was in contrast to an abrupt recovery after about 5 s in cultures with absent GABA_Aergic synaptic transmission. In these GABA deficient cultures evoked responses were prolonged and had multiple discharges. Furthermore, the network excitability changed periodically, with a very slow spontaneous change of the overall network activity in the minute range, which was not observed in cultures with absent GABA_Aergic synaptic transmission. The electrically evoked activity of cultured cortical networks, therefore, is governed by at least two potentially interacting mechanisms: A refractory period in the order of a few seconds and a very slow GABA dependent oscillation of the network excitability.

Slow dynamics in features of synchronized neural network responses

In this report trial-to-trial variations in the synchronized responses of neural networks are explored over time scales of minutes, in ex-vivo large scale cortical networks. We show that sub-second measures of the individual synchronous response, namely-its latency and decay duration, are related to minutes-scale network response dynamics. Network responsiveness is reflected as residency in, or shifting amongst, areas of the latency-decay plane. The different sensitivities of latency and decay durations to synaptic blockers imply that these two measures reflect aspects of inhibitory and excitatory activities. Taken together, the data suggest that trial-to-trial variations in the synchronized responses of neural networks might be related to effective excitation-inhibition ratio being a dynamic variable over time scales of minutes.

Synaptic dynamics contribute to long-term single neuron response fluctuations

Firing rate variability at the single neuron level is characterized by long-memory processes and complex statistics over a wide range of time scales (from milliseconds up to several hours). Here, we focus on the contribution of non-stationary efficacy of the ensemble of synapses–activated in response to a given stimulus–on single neuron response variability. We present and validate a method tailored for controlled and specific long-term activation of a single cortical neuron in vitro via synaptic or antidromic stimulation, enabling a clear separation between two determinants of neuronal response variability: membrane excitability dynamics vs. synaptic dynamics. Applying this method we show that, within the range of physiological activation frequencies, the synaptic ensemble of a given neuron is a key contributor to the neuronal response variability, long-memory processes and complex statistics observed over extended time scales. Synaptic transmission dynamics impact on response variability in stimulation rates that are substantially lower compared to stimulation rates that drive excitability resources to fluctuate. Implications to network embedded neurons are discussed.

Controlling neural network responsiveness: tradeoffs and constraints

In recent years much effort is invested in means to control neural population responses at the whole brain level, within the context of developing advanced medical applications. The tradeoffs and constraints involved, however, remain elusive due to obvious complications entailed by studying whole brain dynamics. Here, we present effective control of response features (probability and latency) of cortical networks in vitro over many hours, and offer this approach as an experimental toy for studying controllability of neural networks in the wider context. Exercising this approach we show that enforcement of stable high activity rates by means of closed loop control may enhance alteration of underlying global input-output relations and activity dependent dispersion of neuronal pair-wise correlations across the network.

Optimized temporally deconvolved Ca²⁺ imaging allows identification of spatiotemporal activity patterns of CA1 hippocampal ensembles

Hippocampal activity is characterized by the coordinated firing of a subset of neurons. Such neuronal ensembles can either be driven by external stimuli to form new memory traces or be reactivated by intrinsic mechanisms to reactivate and consolidate old memories. Hippocampal network oscillations orchestrate this coherent activity. One key question is how the topology, i.e. the functional connectivity of neuronal networks supports their desired function. Recently, this has been addressed by characterizing the intrinsic properties for the highly recurrently connected CA3 region using organotypic slice cultures and Ca(2+) imaging. In the present study, we aimed to determine the properties of CA1 hippocampal ensembles at high temporal and multiple single cell resolution. Thus, we performed Ca(2+) imaging using the chemical fluorescent Ca(2+) indicator Oregon Green BAPTA 1-AM. To achieve most physiological conditions, we used acute hippocampal slices that were recorded in a so-called interface chamber. To faithfully reconstruct firing patterns of multiple neurons in the field of view, we optimized deconvolution-based detection of action potential associated Ca(2+) events. Our approach outperformed currently available detection algorithms by its sensitivity and robustness. In combination with advanced network analysis, we found that acute hippocampal slices contain a median of 11 CA1 neuronal ensembles with a median size of 4 neurons. This apparently low number of neurons is likely due to the confocal imaging acquisition and therefore yields a lower limit. The distribution of ensemble sizes was compatible with a scale-free topology, as far as can be judged from data with small cell numbers. Interestingly, cells were more tightly clustered in large ensembles than in smaller groups. Together, our data show that spatiotemporal activity patterns of hippocampal neuronal ensembles can be reliably detected with deconvolution-based imaging techniques in mouse hippocampal slices. The here presented techniques are fully applicable to similar studies of distributed optical measurements of neuronal activity (in vivo), where signal-to-noise ratio is critical.

Parameters for burst detection

Bursts of action potentials within neurons and throughout networks are believed to serve roles in how neurons handle and store information, both in vivo and in vitro. Accurate detection of burst occurrences and durations are therefore crucial for many studies. A number of algorithms have been proposed to do so, but a standard method has not been adopted. This is due, in part, to many algorithms requiring the adjustment of multiple ad-hoc parameters and further post-hoc criteria in order to produce satisfactory results. Here, we broadly catalog existing approaches and present a new approach requiring the selection of only a single parameter: the number of spikes N comprising the smallest burst to consider. A burst was identified if N spikes occurred in less than T ms, where the threshold T was automatically determined from observing a probability distribution of inter-spike-intervals. Performance was compared vs. different classes of detectors on data gathered from in vitro neuronal networks grown over microelectrode arrays. Our approach offered a number of useful features including: a simple implementation, no need for ad-hoc or post-hoc criteria, and precise assignment of burst boundary time points. Unlike existing approaches, detection was not biased toward larger bursts, allowing identification and analysis of a greater range of neuronal and network dynamics.

DOC2B and Munc13-1 differentially regulate neuronal network activity

Alterations in the levels of synaptic proteins affect synaptic transmission and synaptic plasticity. However, the precise effects on neuronal network activity are still enigmatic. Here, we utilized microelectrode array (MEA) to elucidate how manipulation of the presynaptic release process affects the activity of neuronal networks. By combining pharmacological tools and genetic manipulation of synaptic proteins, we show that overexpression of DOC2B and Munc13-1, proteins known to promote vesicular maturation and release, elicits opposite effects on the activity of the neuronal network. Although both cause an increase in the overall number of spikes, the distribution of spikes is different. While DOC2B enhances, Munc13-1 reduces the firing rate within bursts of spikes throughout the network; however, Munc13-1 increases the rate of network bursts. DOC2B's effects were mimicked by Strontium that elevates asynchronous release but not by a DOC2B mutant that enhances spontaneous release rate. This suggests for the first time that increased asynchronous release on the single-neuron level promotes bursting activity in the network level. This innovative study demonstrates the complementary role of the network level in explaining the physiological relevance of the cellular activity of presynaptic proteins and the transformation of synaptic release manipulation from the neuron to the network level.