Interaction of electrically evoked responses in networks of dissociated cortical neurons

Prediction in cultured cortical neural networks

Theory suggest that networks of neurons may predict their input. Prediction may underlie most aspects of information processing and is believed to be involved in motor and cognitive control and decision-making. Retinal cells have been shown to be capable of predicting visual stimuli, and there is some evidence for prediction of input in the visual cortex and hippocampus. However, there is no proof that the ability to predict is a generic feature of neural networks. We investigated whether random in vitro neuronal networks can predict stimulation, and how prediction is related to short- and long-term memory. To answer these questions, we applied two different stimulation modalities. Focal electrical stimulation has been shown to induce long-term memory traces, whereas global optogenetic stimulation did not. We used mutual information to quantify how much activity recorded from these networks reduces the uncertainty of upcoming stimuli (prediction) or recent past stimuli (short-term memory). Cortical neural networks did predict future stimuli, with the majority of all predictive information provided by the immediate network response to the stimulus. Interestingly, prediction strongly depended on short-term memory of recent sensory inputs during focal as well as global stimulation. However, prediction required less short-term memory during focal stimulation. Furthermore, the dependency on short-term memory decreased during 20 h of focal stimulation, when long-term connectivity changes were induced. These changes are fundamental for long-term memory formation, suggesting that besides short-term memory the formation of long-term memory traces may play a role in efficient prediction.

Plasticity and Adaptation in Neuromorphic Biohybrid Systems

Neuromorphic systems take inspiration from the principles of biological information processing to form hardware platforms that enable the large-scale implementation of neural networks. The recent years have seen both advances in the theoretical aspects of spiking neural networks for their use in classification and control tasks and a progress in electrophysiological methods that is pushing the frontiers of intelligent neural interfacing and signal processing technologies. At the forefront of these new technologies, artificial and biological neural networks are tightly coupled, offering a novel "biohybrid" experimental framework for engineers and neurophysiologists. Indeed, biohybrid systems can constitute a new class of neuroprostheses opening important perspectives in the treatment of neurological disorders. Moreover, the use of biologically plausible learning rules allows forming an overall fault-tolerant system of co-developing subsystems. To identify opportunities and challenges in neuromorphic biohybrid systems, we discuss the field from the perspectives of neurobiology, computational neuroscience, and neuromorphic engineering.

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.

Spike Code Flow in Cultured Neuronal Networks

We observed spike trains produced by one-shot electrical stimulation with 8 × 8 multielectrodes in cultured neuronal networks. Each electrode accepted spikes from several neurons. We extracted the short codes from spike trains and obtained a code spectrum with a nominal time accuracy of 1%. We then constructed code flow maps as movies of the electrode array to observe the code flow of “1101” and “1011,” which are typical pseudorandom sequence such as that we often encountered in a literature and our experiments. They seemed to flow from one electrode to the neighboring one and maintained their shape to some extent. To quantify the flow, we calculated the “maximum cross-correlations” among neighboring electrodes, to find the direction of maximum flow of the codes with lengths less than 8. Normalized maximum cross-correlations were almost constant irrespective of code. Furthermore, if the spike trains were shuffled in interval orders or in electrodes, they became significantly small. Thus, the analysis suggested that local codes of approximately constant shape propagated and conveyed information across the network. Hence, the codes can serve as visible and trackable marks of propagating spike waves as well as evaluating information flow in the neuronal network.

Simulation of Code Spectrum and Code Flow of Cultured Neuronal Networks

It has been shown that, in cultured neuronal networks on a multielectrode, pseudorandom-like sequences (codes) are detected, and they flow with some spatial decay constant. Each cultured neuronal network is characterized by a specific spectrum curve. That is, we may consider the spectrum curve as a “signature” of its associated neuronal network that is dependent on the characteristics of neurons and network configuration, including the weight distribution. In the present study, we used an integrate-and-fire model of neurons with intrinsic and instantaneous fluctuations of characteristics for performing a simulation of a code spectrum from multielectrodes on a 2D mesh neural network. We showed that it is possible to estimate the characteristics of neurons such as the distribution of number of neurons around each electrode and their refractory periods. Although this process is a reverse problem and theoretically the solutions are not sufficiently guaranteed, the parameters seem to be consistent with those of neurons. That is, the proposed neural network model may adequately reflect the behavior of a cultured neuronal network. Furthermore, such prospect is discussed that code analysis will provide a base of communication within a neural network that will also create a base of natural intelligence.

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.

Understanding the temporal evolution of neuronal connectivity in cultured networks using statistical analysis

BACKGROUND

Micro-Electrode Array (MEA) technology allows researchers to perform long-term non-invasive neuronal recordings in-vitro while actively interacting with the cultured neurons. Despite numerous studies carried out using MEAs, many functional, chemical and structural mechanisms of how dissociated cortical neurons develop and respond to external stimuli are not yet well understood because of the lack of quantitative studies that assess how their development can be affected by chronic external stimulation.

METHODS

To investigate network changes, we analyzed a large MEA data set composed of neuron spikes recorded from cultures of dissociated rat cortical neurons plated on MEA dishes with 59 recording electrodes each. Neural network activity was recorded during the first five weeks of each culture's in-vitro development. Stimulation sessions were delivered to each of the 59 electrodes. The False Discovery Rate technique was used to quantify the temporal evolution of dissociated cortical neurons. Our analysis focused on network responses that occurred within selected time window durations, namely 50 ms, 100 ms and 150 ms after stimulus onset.

RESULTS

Our results show an evolution in dissociated cortical neuronal network activity over time, that reflects the network synaptic evolution. Furthermore, we tested the sensitivity of our technique to different observation time windows and found that varying the time windows, allows us to capture different dynamics of the observed responses. In addition, when selecting a 150 ms observation time window, our findings indicate that cultures dissociated from the same brain tissue display trends in their temporal evolution that are more similar than those obtained from different brains.

CONCLUSION

Our results emphasize that the FDR technique can be implemented without the need to make any particular assumptions about the data a priori. The proposed technique was able to capture the well-known dissociated cortical neuron networks' temporal evolution, that has been previously observed in in-vivo and in intact brain tissue studies. Furthermore, our findings suggest that the time window that is used to capture the stimulus-evoked network responses is a critical parameter to analyze the electrical behavioral and temporal evolution of dissociated cortical neurons.

Detection of M-sequences from spike sequence in neuronal networks

In circuit theory, it is well known that a linear feedback shift register (LFSR) circuit generates pseudorandom bit sequences (PRBS), including an M-sequence with the maximum period of length. In this study, we tried to detect M-sequences known as a pseudorandom sequence generated by the LFSR circuit from time series patterns of stimulated action potentials. Stimulated action potentials were recorded from dissociated cultures of hippocampal neurons grown on a multielectrode array. We could find several M-sequences from a 3-stage LFSR circuit (M3). These results show the possibility of assembling LFSR circuits or its equivalent ones in a neuronal network. However, since the M3 pattern was composed of only four spike intervals, the possibility of an accidental detection was not zero. Then, we detected M-sequences from random spike sequences which were not generated from an LFSR circuit and compare the result with the number of M-sequences from the originally observed raster data. As a result, a significant difference was confirmed: a greater number of “0–1” reversed the 3-stage M-sequences occurred than would have accidentally be detected. This result suggests that some LFSR equivalent circuits are assembled in neuronal networks.

Spatial-temporal dynamics of chaotic behavior in cultured hippocampal networks

Using multiple nonlinear techniques, we revealed the existence of chaos in the spontaneous activity of neuronal networks in vitro. The spatial-temporal dynamics of these networks indicated that emergent transition between chaotic behavior and superburst occurred periodically in low-frequency oscillations. An analysis of network-wide activity indicated that chaos was synchronized among different sites. Moreover, we found that the degree of chaos increased as the number of active sites in the network increased during long-term development (over three months in vitro). The chaotic behavior of the dissociated networks had similar spatial-temporal characteristics (rapid transition, periodicity, and synchronization) as the intact brain; however, the degree of chaos depended on the number of active sites at the mesoscopic level. This work could provide insight into neural coding and neurocybernetics.