Management of synchronized network activity by highly active neurons

Increasing evidence supports the idea that spontaneous brain activity may have an important functional role. Cultured neuronal networks provide a suitable model system to search for the mechanisms by which neuronal spontaneous activity is maintained and regulated. This activity is marked by synchronized bursting events (SBEs)--short time windows (hundreds of milliseconds) of rapid neuronal firing separated by long quiescent periods (seconds). However, there exists a special subset of rapidly firing neurons whose activity also persists between SBEs. It has been proposed that these highly active (HA) neurons play an important role in the management (i.e. establishment, maintenance and regulation) of the synchronized network activity. Here, we studied the dynamical properties and the functional role of HA neurons in homogeneous and engineered networks, during early network development, upon recovery from chemical inhibition and in response to electrical stimulations. We found that their sequences of inter-spike intervals (ISI) exhibit long time correlations and a unimodal distribution. During the network's development and under intense inhibition, the observed activity follows a transition period during which mostly HA neurons are active. Studying networks with engineered geometry, we found that HA neurons are precursors (the first to fire) of the spontaneous SBEs and are more responsive to electrical stimulations.     

Changing excitation and inhibition in simulated neural networks: effects on induced bursting behavior

 The development of synchronous bursting in neuronal ensembles represents an important change in network behavior. To determine the influences on development of such synchronous bursting behavior we study the dynamics of small networks of sparsely connected excitatory and inhibitory neurons using numerical simulations. The synchronized bursting activities in networks evoked by background spikes are investigated. Specifically, patterns of bursting activity are examined when the balance between excitation and inhibition on neuronal inputs is varied and the fraction of inhibitory neurons in the network is changed. For quantitative comparison of bursting activities in networks, measures of the degree of synchrony are used. We demonstrate how changes in the strength of excitation on inputs of neurons can be compensated by changes in the strength of inhibition without changing the degree of synchrony in the network. The effects of changing several network parameters on the network activity are analyzed and discussed. These changes may underlie the transition of network activity from normal to potentially pathologic (e.g., epileptic) states. Received: 21 May 2002 / Accepted in revised form: 3 December 2002 / Published online: 7 March 2003 Correspondence to: P. Kudela (e-mail: pkudela@jhmi.edu) Acknowledgements. This research was supported by NIH grant NS 38958.     

Firing synchronization and temporal order in noisy neuronal networks

Noise-induced complete synchronization and frequency synchronization in coupled spiking and bursting neurons are studied firstly. The effects of noise and coupling are discussed. It is found that bursting neurons are easier to achieve firing synchronization than spiking ones, which means that bursting activities are more important for information transfer in neuronal networks. Secondly, the effects of noise on firing synchronization in a noisy map neuronal network are presented. Noise-induced synchronization and temporal order are investigated by means of the firing rate function and the order index. Firing synchronization and temporal order of excitatory neurons can be greatly enhanced by subthreshold stimuli with resonance frequency. Finally, it is concluded that random perturbations play an important role in firing activities and temporal order in neuronal networks.     

Bursting Reverberation as a Multiscale Neuronal Network Process Driven by Synaptic Depression-Facilitation

Neuronal networks can generate complex patterns of activity that depend on membrane properties of individual neurons as well as on functional synapses. To decipher the impact of synaptic properties and connectivity on neuronal network behavior, we investigate the responses of neuronal ensembles from small (5–30 cells in a restricted sphere) and large (acute hippocampal slice) networks to single electrical stimulation: in both cases, a single stimulus generated a synchronous long-lasting bursting activity. While an initial spike triggered a reverberating network activity that lasted 2–5 seconds for small networks, we found here that it lasted only up to 300 milliseconds in slices. To explain this phenomena present at different scales, we generalize the depression-facilitation model and extracted the network time constants. The model predicts that the reverberation time has a bell shaped relation with the synaptic density, revealing that the bursting time cannot exceed a maximum value. Furthermore, before reaching its maximum, the reverberation time increases sub-linearly with the synaptic density of the network. We conclude that synaptic dynamics and connectivity shape the mean burst duration, a property present at various scales of the networks. Thus bursting reverberation is a property of sufficiently connected neural networks, and can be generated by collective depression and facilitation of underlying functional synapses.     

Suppressing bursting synchronization in a modular neuronal network with synaptic plasticity

Excessive synchronization of neurons in cerebral cortex is believed to play a crucial role in the emergence of neuropsychological disorders such as Parkinson’s disease, epilepsy and essential tremor. This study, by constructing a modular neuronal network with modified Oja’s learning rule, explores how to eliminate the pathological synchronized rhythm of interacted busting neurons numerically. When all neurons in the modular neuronal network are strongly synchronous within a specific range of coupling strength, the result reveals that synaptic plasticity with large learning rate can suppress bursting synchronization effectively. For the relative small learning rate not capable of suppressing synchronization, the technique of nonlinear delayed feedback control including differential feedback control and direct feedback control is further proposed to reduce the synchronized bursting state of coupled neurons. It is demonstrated that the two kinds of nonlinear feedback control can eliminate bursting synchronization significantly when the control parameters of feedback strength and feedback delay are appropriately tuned. For the former control technique, the control domain of effective synchronization suppression is similar to a semi-elliptical domain in the simulated parameter space of feedback strength and feedback delay, while for the latter one, the effective control domain is similar to a fan-shaped domain in the simulated parameter space.     

Innate synchronous oscillations in freely-organized small neuronal circuits

Background

Information processing in neuronal networks relies on the network''s ability to generate temporal patterns of action potentials. Although the nature of neuronal network activity has been intensively investigated in the past several decades at the individual neuron level, the underlying principles of the collective network activity, such as the synchronization and coordination between neurons, are largely unknown. Here we focus on isolated neuronal clusters in culture and address the following simple, yet fundamental questions: What is the minimal number of cells needed to exhibit collective dynamics? What are the internal temporal characteristics of such dynamics and how do the temporal features of network activity alternate upon crossover from minimal networks to large networks?

Methodology/Principal Findings

We used network engineering techniques to induce self-organization of cultured networks into neuronal clusters of different sizes. We found that small clusters made of as few as 40 cells already exhibit spontaneous collective events characterized by innate synchronous network oscillations in the range of 25 to 100 Hz. The oscillation frequency of each network appeared to be independent of cluster size. The duration and rate of the network events scale with cluster size but converge to that of large uniform networks. Finally, the investigation of two coupled clusters revealed clear activity propagation with master/slave asymmetry.

Conclusions/Significance

The nature of the activity patterns observed in small networks, namely the consistent emergence of similar activity across networks of different size and morphology, suggests that neuronal clusters self-regulate their activity to sustain network bursts with internal oscillatory features. We therefore suggest that clusters of as few as tens of cells can serve as a minimal but sufficient functional network, capable of sustaining oscillatory activity. Interestingly, the frequencies of these oscillations are similar those observed in vivo.     

An integrate-and-fire model for synchronized bursting in a network of cultured cortical neurons

It has been suggested that spontaneous synchronous neuronal activity is an essential step in the formation of functional networks in the central nervous system. The key features of this type of activity consist of bursts of action potentials with associated spikes of elevated cytoplasmic calcium. These features are also observed in networks of rat cortical neurons that have been formed in culture. Experimental studies of these cultured networks have led to several hypotheses for the mechanisms underlying the observed synchronized oscillations. In this paper, bursting integrate-and-fire type mathematical models for regular spiking (RS) and intrinsic bursting (IB) neurons are introduced and incorporated through a small-world connection scheme into a two-dimensional excitatory network similar to those in the cultured network. This computer model exhibits spontaneous synchronous activity through mechanisms similar to those hypothesized for the cultured experimental networks. Traces of the membrane potential and cytoplasmic calcium from the model closely match those obtained from experiments. We also consider the impact on network behavior of the IB neurons, the geometry and the small world connection scheme. Action Editor: David Golomb     

Functional clustering in hippocampal cultures: relating network structure and dynamics

In this work we investigate the relationship between gross anatomic structural network properties, neuronal dynamics and the resultant functional structure in dissociated rat hippocampal cultures. Specifically, we studied cultures as they developed under two conditions: the first supporting glial cell growth (high glial group), and the second one inhibiting it (low glial group). We then compared structural network properties and the spatio-temporal activity patterns of the neurons. Differences in dynamics between the two groups could be linked to the impact of the glial network on the neuronal network as the cultures developed. We also implemented a recently developed algorithm called the functional clustering algorithm (FCA) to obtain the resulting functional network structure. We show that this new algorithm is useful for capturing changes in functional network structure as the networks evolve over time. The FCA detects changes in functional structure that are consistent with expected dynamical differences due to the impact of the glial network. Cultures in the high glial group show an increase in global synchronization as the cultures age, while those in the low glial group remain locally synchronized. We additionally use the FCA to quantify the amount of synchronization present in the cultures and show that the total level of synchronization in the high glial group is stronger than in the low glial group. These results indicate an interdependence between the glial and neuronal networks present in dissociated cultures.     

Network bursts in cortical cultures are best simulated using pacemaker neurons and adaptive synapses

One of the most specific and exhibited features in the electrical activity of dissociated cultured neural networks (NNs) is the phenomenon of synchronized bursts, whose profiles vary widely in shape, width and firing rate. On the way to understanding the organization and behavior of biological NNs, we reproduced those features with random connectivity network models with 5,000 neurons. While the common approach to induce bursting behavior in neuronal network models is noise injection, there is experimental evidence suggesting the existence of pacemaker-like neurons. In our simulations noise did evoke bursts, but with an unrealistically gentle rising slope. We show that a small subset of ‘pacemaker’ neurons can trigger bursts with a more realistic profile. We found that adding pacemaker-like neurons as well as adaptive synapses yield burst features (shape, width, and height of the main phase) in the same ranges as obtained experimentally. Finally, we demonstrate how changes in network connectivity, transmission delays, and excitatory fraction influence network burst features quantitatively.     

Manifestation of function-follow-form in cultured neuronal networks

We expose hidden function-follow-form schemata in the recorded activity of cultured neuronal networks by comparing the activity with simulation results of a new modeling approach. Cultured networks grown from an arbitrary mixture of neuron and glia cells in the absence of external stimulations and chemical cues spontaneously form networks of different sizes (from 50 to several millions of neurons) that exhibit non-arbitrary complex spatio-temporal patterns of activity. The latter is marked by formation of a sequence of synchronized bursting events (SBEs)--short time windows (approximately 200 ms) of rapid neuron firing, separated by longer time intervals (seconds) of sporadic neuron firing. The new dynamical synapse and soma (DSS) model, used here, has been successful in generating sequences of SBEs with the same statistical scaling properties (over six time decades) as those of the small networks. Large networks generate statistically distinct sub-groups of SBEs, each with its own characteristic pattern of neuronal firing ('fingerprint'). This special function (activity) motif has been proposed to emanate from a structural (form) motif--self-organization of the large networks into a fabric of overlapping sub-networks of about 1 mm in size. Here we test this function-follow-form idea by investigating the influence of the connectivity architecture of a model network (form) on the structure of its spontaneous activity (function). We show that a repertoire of possible activity states similar to the observed ones can be generated by networks with proper underlying architecture. For example, networks composed of two overlapping sub-networks exhibit distinct types of SBEs, each with its own characteristic pattern of neuron activity that starts at a specific sub-network. We further show that it is possible to regulate the temporal appearance of the different sub-groups of SBEs by an additional non-synaptic current fed into the soma of the modeled neurons. The ability to regulate the relative temporal ordering of different SBEs might endow the networks with higher plasticity and complexity. These findings call for additional mechanisms yet to be discovered. Recent experimental observations indicate that glia cells coupled to neuronal soma might generate such non-synaptic regulating currents.     

Paradoxical phase response of gamma rhythms facilitates their entrainment in heterogeneous networks

The synchronization of different γ-rhythms arising in different brain areas has been implicated in various cognitive functions. Here, we focus on the effect of the ubiquitous neuronal heterogeneity on the synchronization of ING (interneuronal network gamma) and PING (pyramidal-interneuronal network gamma) rhythms. The synchronization properties of rhythms depends on the response of their collective phase to external input. We therefore determine the macroscopic phase-response curve for finite-amplitude perturbations (fmPRC) of ING- and PING-rhythms in all-to-all coupled networks comprised of linear (IF) or quadratic (QIF) integrate-and-fire neurons. For the QIF networks we complement the direct simulations with the adjoint method to determine the infinitesimal macroscopic PRC (imPRC) within the exact mean-field theory. We show that the intrinsic neuronal heterogeneity can qualitatively modify the fmPRC and the imPRC. Both PRCs can be biphasic and change sign (type II), even though the phase-response curve for the individual neurons is strictly non-negative (type I). Thus, for ING rhythms, say, external inhibition to the inhibitory cells can, in fact, advance the collective oscillation of the network, even though the same inhibition would lead to a delay when applied to uncoupled neurons. This paradoxical advance arises when the external inhibition modifies the internal dynamics of the network by reducing the number of spikes of inhibitory neurons; the advance resulting from this disinhibition outweighs the immediate delay caused by the external inhibition. These results explain how intrinsic heterogeneity allows ING- and PING-rhythms to become synchronized with a periodic forcing or another rhythm for a wider range in the mismatch of their frequencies. Our results identify a potential function of neuronal heterogeneity in the synchronization of coupled γ-rhythms, which may play a role in neural information transfer via communication through coherence.     

Oscillations and synchrony in a cortical neural network

In this paper, the oscillations and synchronization status of two different network connectivity patterns based on Izhikevich model are studied. One of the connectivity patterns is a randomly connected neuronal network, the other one is a small-world neuronal network. This Izhikevich model is a simple model which can not only reproduce the rich behaviors of biological neurons but also has only two equations and one nonlinear term. Detailed investigations reveal that by varying some key parameters, such as the connection weights of neurons, the external current injection, the noise of intensity and the neuron number, this neuronal network will exhibit various collective behaviors in randomly coupled neuronal network. In addition, we show that by changing the number of nearest neighbor and connection probability in small-world topology can also affect the collective dynamics of neuronal activity. These results may be instructive in understanding the collective dynamics of mammalian cortex.     

Cortical network modeling: analytical methods for firing rates and some properties of networks of LIF neurons.

The circuitry of cortical networks involves interacting populations of excitatory (E) and inhibitory (I) neurons whose relationships are now known to a large extent. Inputs to E- and I-cells may have their origins in remote or local cortical areas. We consider a rudimentary model involving E- and I-cells. One of our goals is to test an analytic approach to finding firing rates in neural networks without using a diffusion approximation and to this end we consider in detail networks of excitatory neurons with leaky integrate-and-fire (LIF) dynamics. A simple measure of synchronization, denoted by S(q), where q is between 0 and 100 is introduced. Fully connected E-networks have a large tendency to become dominated by synchronously firing groups of cells, except when inputs are relatively weak. We observed random or asynchronous firing in such networks with diverse sets of parameter values. When such firing patterns were found, the analytical approach was often able to accurately predict average neuronal firing rates. We also considered several properties of E-E networks, distinguishing several kinds of firing pattern. Included were those with silences before or after periods of intense activity or with periodic synchronization. We investigated the occurrence of synchronized firing with respect to changes in the internal excitatory postsynaptic potential (EPSP) magnitude in a network of 100 neurons with fixed values of the remaining parameters. When the internal EPSP size was less than a certain value, synchronization was absent. The amount of synchronization then increased slowly as the EPSP amplitude increased until at a particular EPSP size the amount of synchronization abruptly increased, with S(5) attaining the maximum value of 100%. We also found network frequency transfer characteristics for various network sizes and found a linear dependence of firing frequency over wide ranges of the external afferent frequency, with non-linear effects at lower input frequencies. The theory may also be applied to sparsely connected networks, whose firing behaviour was found to change abruptly as the probability of a connection passed through a critical value. The analytical method was also found to be useful for a feed-forward excitatory network and a network of excitatory and inhibitory neurons.     

Activity-dependent clustering of functional synaptic inputs on developing hippocampal dendrites

During brain development, before sensory systems become functional, neuronal networks spontaneously generate repetitive bursts of neuronal activity, which are typically synchronized across many neurons. Such activity patterns have been described on the level of networks and cells, but the fine-structure of inputs received by an individual neuron during spontaneous network activity has not been studied. Here, we used calcium imaging to record activity at many synapses of hippocampal pyramidal neurons simultaneously to establish the activity patterns in the majority of synapses of an entire cell. Analysis of the spatiotemporal patterns of synaptic activity revealed a fine-scale connectivity rule: neighboring synapses (<16?μm intersynapse distance) are more likely to be coactive than synapses that are farther away from each other. Blocking spiking activity or NMDA receptor activation revealed that the clustering of synaptic inputs required neuronal activity, demonstrating a role of developmentally expressed spontaneous activity for connecting neurons with subcellular precision.     

Measuring synchronization in neuronal networks for biosensor applications

Cultures of neurons can be grown on microelectrode arrays (MEAs), so that their spike and burst activity can be monitored. These activity patterns are quite sensitive to changes in the environment, such as chemical exposure, and hence the cultures can be used as biosensors. One key issue in analyzing the data from neuronal networks is how to quantify the level of synchronization among different units, which represent different neurons in the network. In this paper, we propose a synchronization metric, based on the statistical distribution of unit-to-unit correlation coefficients. We show that this synchronization metric changes significantly when the networks are exposed to bicuculline, strychnine, or 2,3-dioxo-6-nitro-l,2,3,4-tetrahydrobenzoquinoxaline-7-sulphonamide (NBQX). For that reason, this metric can be used to characterize pharmacologically induced changes in a network, either for research or for biosensor applications.     

Learning-Induced Synchronization and Plasticity of a Developing Neural Network

Learning-induced synchronization of a neural network at various developing stages is studied by computer simulations using a pulse-coupled neural network model in which the neuronal activity is simulated by a one-dimensional map. Two types of Hebbian plasticity rules are investigated and their differences are compared. For both models, our simulations show a logarithmic increase in the synchronous firing frequency of the network with the culturing time of the neural network. This result is consistent with recent experimental observations. To investigate how to control the synchronization behavior of a neural network after learning, we compare the occurrence of synchronization for four networks with different designed patterns under the influence of an external signal. The effect of such a signal on the network activity highly depends on the number of connections between neurons. We discuss the synaptic plasticity and enhancement effects for a random network after learning at various developing stages.     

Shaping bursting by electrical coupling and noise

Gap-junctional coupling is an important way of communication between neurons and other excitable cells. Strong electrical coupling synchronizes activity across cell ensembles. Surprisingly, in the presence of noise synchronous oscillations generated by an electrically coupled network may differ qualitatively from the oscillations produced by uncoupled individual cells forming the network. A prominent example of such behavior is the synchronized bursting in islets of Langerhans formed by pancreatic β-cells, which in isolation are known to exhibit irregular spiking (Sherman and Rinzel, Biophys J 54:411–425, 1988; Sherman and Rinzel, Biophys J 59:547–559, 1991). At the heart of this intriguing phenomenon lies denoising, a remarkable ability of electrical coupling to diminish the effects of noise acting on individual cells. In this paper, building on an earlier analysis of denoising in networks of integrate-and-fire neurons (Medvedev, Neural Comput 21 (11):3057–3078, 2009) and our recent study of spontaneous activity in a closely related model of the Locus Coeruleus network (Medvedev and Zhuravytska, The geometry of spontaneous spiking in neuronal networks, submitted, 2012), we derive quantitative estimates characterizing denoising in electrically coupled networks of conductance-based models of square wave bursting cells. Our analysis reveals the interplay of the intrinsic properties of the individual cells and network topology and their respective contributions to this important effect. In particular, we show that networks on graphs with large algebraic connectivity (Fiedler, Czech Math J 23(98):298–305, 1973) or small total effective resistance (Bollobas, Modern graph theory, Graduate Texts in Mathematics, vol. 184, Springer, New York, 1998) are better equipped for implementing denoising. As a by-product of the analysis of denoising, we analytically estimate the rate with which trajectories converge to the synchronization subspace and the stability of the latter to random perturbations. These estimates reveal the role of the network topology in synchronization. The analysis is complemented by numerical simulations of electrically coupled conductance-based networks. Taken together, these results explain the mechanisms underlying synchronization and denoising in an important class of biological models.     

Robust emergence of small-world structure in networks of spiking neurons

Spontaneous activity in biological neural networks shows patterns of dynamic synchronization. We propose that these patterns support the formation␣of a small-world structure—network connectivity␣optimal for distributed information processing. We␣present numerical simulations with connected Hindmarsh–Rose neurons in which, starting from random connection distributions, small-world networks evolve as a result of applying an adaptive rewiring rule. The rule connects pairs of neurons that tend fire in synchrony, and disconnects ones that fail to synchronize. Repeated application of the rule leads to small-world structures. This mechanism is robustly observed for bursting and irregular firing regimes.     

Connecting the dots in ethology: applying network theory to understand neural and animal collectives

A major goal shared by neuroscience and collective behavior is to understand how dynamic interactions between individual elements give rise to behaviors in populations of neurons and animals, respectively. This goal has recently become within reach, thanks to techniques providing access to the connectivity and activity of neuronal ensembles as well as to behaviors among animal collectives. The next challenge using these datasets is to unravel network mechanisms generating population behaviors. This is aided by network theory, a field that studies structure–function relationships in interconnected systems. Here we review studies that have taken a network view on modern datasets to provide unique insights into individual and collective animal behaviors. Specifically, we focus on how analyzing signal propagation, controllability, symmetry, and geometry of networks can tame the complexity of collective system dynamics. These studies illustrate the potential of network theory to accelerate our understanding of behavior across ethological scales.     

Uncovering the synchronization dynamics from correlated neuronal activity quantifies assembly formation

 Synchronous network excitation is believed to play an outstanding role in neuronal information processing. Due to the stochastic nature of the contributing neurons, however, those synchronized states are difficult to detect in electrode recordings. We present a framework and a model for the identification of such network states and of their dynamics in a specific experimental situation. Our approach operationalizes the notion of neuronal groups forming assemblies via synchronization based on experimentally obtained spike trains. The dynamics of such groups is reflected in the sequence of synchronized states, which we describe as a renewal dynamics. We furthermore introduce a rate function which is dependent on the internal network phase that quantifies the activity of neurons contributing to the observed spike train. This constitutes a hidden state model which is formally equivalent to a hidden Markov model, and all its parameters can be accurately determined from the experimental time series using the Baum-Welch algorithm. We apply our method to recordings from the cat visual cortex which exhibit oscillations and synchronizations. The parameters obtained for the hidden state model uncover characteristic properties of the system including synchronization, oscillation, switching, background activity and correlations. In applications involving multielectrode recordings, the extracted models quantify the extent of assembly formation and can be used for a temporally precise localization of system states underlying a specific spike train. Received: 30 March 1993/Accepted in revised form: 16 April 1994